PLoS Computational Biology最新文献

Over the past two decades, extensive studies, particularly in cancer analysis through large datasets like The Cancer Genome Atlas (TCGA), have aimed at improving patient therapies and precision medicine. However, limited overlap and inconsistencies among gene signatures across different cohorts pose challenges. The dynamic nature of the transcriptome, encompassing diverse RNA species and functional complexities at gene and isoform levels, introduces intricacies, and current gene signatures face reproducibility issues due to the unique transcriptomic landscape of each patient. In this context, discrepancies arising from diverse sequencing technologies, data analysis algorithms, and software tools further hinder consistency. While careful experimental design, analytical strategies, and standardized protocols could enhance reproducibility, future prospects lie in multiomics data integration, machine learning techniques, open science practices, and collaborative efforts. Standardized metrics, quality control measures, and advancements in single-cell RNA-seq will contribute to unbiased gene signature identification. In this perspective article, we outline some thoughts and insights addressing challenges, standardized practices, and advanced methodologies enhancing the reliability of gene signatures in disease transcriptomic research.

In natural behaviors, multiple neural signals simultaneously drive activation across overlapping brain networks. Due to limitations in the amount of data that can be acquired in common experimental designs, the determination of these interactions is commonly inferred via modeling approaches, which reduce overfitting by finding appropriate regularizing hyperparameters. However, it is unclear whether these hyperparameters can also be related to any aspect of the underlying biological phenomena and help interpret them. We applied a state-of-the-art regularization procedure-automatic locality determination-to interacting neural activations in the mouse posterior cortex associated with movements of the body and eyes. As expected, regularization significantly improved the determination and interpretability of the response interactions. However, regularizing hyperparameters also changed considerably, and seemingly unpredictably, from animal to animal. We found that these variations were not random; rather, they correlated with the variability in visually evoked responses and with the variability in the state of arousal of the animals measured by pupillometry-both pieces of information that were not included in the modeling framework. These observations could be generalized to another commonly used-but potentially less informative-regularization method, ridge regression. Our findings demonstrate that optimal model hyperparameters can be discovery tools that are informative of factors not a priori included in the model's design.

Genome-wide association studies (GWASs) have been widely applied in the neuroimaging field to discover genetic variants associated with brain-related traits. So far, almost all GWASs conducted in neuroimaging genetics are performed on univariate quantitative features summarized from brain images. On the other hand, powerful deep learning technologies have dramatically improved our ability to classify images. In this study, we proposed and implemented a novel machine learning strategy for systematically identifying genetic variants that lead to detectable nuances on Magnetic Resonance Images (MRI). For a specific single nucleotide polymorphism (SNP), if MRI images labeled by genotypes of this SNP can be reliably distinguished using machine learning, we then hypothesized that this SNP is likely to be associated with brain anatomy or function which is manifested in MRI brain images. We applied this strategy to a catalog of MRI image and genotype data collected by the Alzheimer's Disease Neuroimaging Initiative (ADNI) consortium. From the results, we identified novel variants that show strong association to brain phenotypes.

Owing to its specialised methodology, palaeoecology is often regarded as a separate field from ecology, even though it is essential for understanding long-term ecological processes that have shaped the ecosystems that ecologists study and manage. Despite advances in ecological modelling, sample dating, and proxy-based reconstructions facilitating direct comparison of palaeoecological data with neo-ecological data, most of the scientific knowledge derived from palaeoecological studies remains siloed. We surveyed a group of palaeo-researchers with experience in crossing the divide between palaeoecology and neo-ecology, to develop Ten Simple Rules for publishing your palaeoecological research in non-palaeo journals. Our 10 rules are divided into the preparation phase, writing phase, and finalising phase when the article is submitted to the target journal. These rules provide a suite of strategies, including improved networking early in the process, building effective collaborations, transmitting results more efficiently and cross-disciplinary, and integrating concepts and methodologies that appeal to ecologists and a wider readership. Adhering to these Ten Simple Rules can ensure palaeoecologists' findings are more accessible and impactful among ecologists and the wider scientific community. Although this article primarily shows examples of how palaeoecological studies were published in journals for a broader audience, the rules apply to anyone who aims to publish outside specialised journals.

From a game of darts to neurorehabilitation, the ability to explore and fine tune our movements is critical for success. Past work has shown that exploratory motor behaviour in response to reinforcement (reward) feedback is closely linked with the basal ganglia, while movement corrections in response to error feedback is commonly attributed to the cerebellum. While our past work has shown these processes are dissociable during adaptation, it is unknown how they uniquely impact exploratory behaviour. Moreover, converging neuroanatomical evidence shows direct and indirect connections between the basal ganglia and cerebellum, suggesting that there is an interaction between reinforcement-based and error-based neural processes. Here we examine the unique roles and interaction between reinforcement-based and error-based processes on sensorimotor exploration in a neurotypical population. We also recruited individuals with Parkinson's disease to gain mechanistic insight into the role of the basal ganglia and associated reinforcement pathways in sensorimotor exploration. Across three reaching experiments, participants were given either reinforcement feedback, error feedback, or simultaneously both reinforcement & error feedback during a sensorimotor task that encouraged exploration. Our reaching results, a re-analysis of a previous gait experiment, and our model suggests that in isolation, reinforcement-based and error-based processes respectively boost and suppress exploration. When acting in concert, we found that reinforcement-based and error-based processes interact by mutually opposing one another. Finally, we found that those with Parkinson's disease had decreased exploration when receiving reinforcement feedback, supporting the notion that compromised reinforcement-based processes reduces the ability to explore new motor actions. Understanding the unique and interacting roles of reinforcement-based and error-based processes may help to inform neurorehabilitation paradigms where it is important to discover new and successful motor actions.

Exquisite binding specificity is essential for many protein functions but is difficult to engineer. Many biotechnological or biomedical applications require the discrimination of very similar ligands, which poses the challenge of designing protein sequences with highly specific binding profiles. Experimental methods for generating specific binders rely on in vitro selection, which is limited in terms of library size and control over specificity profiles. Additional control was recently demonstrated through high-throughput sequencing and downstream computational analysis. Here we follow such an approach to demonstrate the design of specific antibodies beyond those probed experimentally. We do so in a context where very similar epitopes need to be discriminated, and where these epitopes cannot be experimentally dissociated from other epitopes present in the selection. Our approach involves the identification of different binding modes, each associated with a particular ligand against which the antibodies are either selected or not. Using data from phage display experiments, we show that the model successfully disentangles these modes, even when they are associated with chemically very similar ligands. Additionally, we demonstrate and validate experimentally the computational design of antibodies with customized specificity profiles, either with specific high affinity for a particular target ligand, or with cross-specificity for multiple target ligands. Overall, our results showcase the potential of leveraging a biophysical model learned from selections against multiple ligands to design proteins with tailored specificity, with applications to protein engineering extending beyond the design of antibodies.

A tight regulation of morphogen production is key for morphogen gradient formation and thereby for reproducible and organised organ development. Although many genetic interactions involved in the establishment of morphogen production domains are known, the biophysical mechanisms of morphogen source formation are poorly understood. Here we addressed this by focusing on the morphogen Sonic hedgehog (Shh) in the vertebrate neural tube. Shh is produced by the adjacently located notochord and by the floor plate of the neural tube. Using a data-constrained computational screen, we identified different possible mechanisms by which floor plate formation can occur, only one of which is consistent with experimental data. In this mechanism, the floor plate is established rapidly in response to Shh from the notochord and the dynamics of regulatory interactions within the neural tube. In this process, uniform activators and Shh-dependent repressors are key for establishing the floor plate size. Subsequently, the floor plate becomes insensitive to Shh and increases in size due to tissue growth, leading to scaling of the floor plate with neural tube size. In turn, this results in scaling of the Shh amplitude with tissue growth. Thus, this mechanism ensures a separation of time scales in floor plate formation, so that the floor plate domain becomes growth-dependent after an initial rapid establishment phase. Our study raises the possibility that the time scale separation between specification and growth might be a common strategy for scaling the morphogen gradient amplitude in growing organs. The model that we developed provides a new opportunity for quantitative studies of morphogen source formation in growing tissues.

Conventional approaches to enhance movement coordination, such as providing instructions and visual feedback, are often inadequate in complex motor tasks with multiple degrees of freedom (DoFs). To effectively address coordination deficits in such complex motor systems, it becomes imperative to develop interventions grounded in a model of human motor learning; however, modeling such learning processes is challenging due to the large DoFs. In this paper, we present a computational motor learning model that leverages the concept of motor synergies to extract low-dimensional learning representations in the high-dimensional motor space and the internal model theory of motor control to capture both fast and slow motor learning processes. We establish the model's convergence properties and validate it using data from a target capture game played by human participants. We study the influence of model parameters on several motor learning trade-offs such as speed-accuracy, exploration-exploitation, satisficing, and flexibility-performance, and show that the human motor learning system tunes these parameters to optimize learning and various output performance metrics.

The fate choices of stem cells between self-renewal and differentiation are often tightly regulated by juxtacrine (cell-cell contact) signalling. Here, we assess how the interplay between cell division, cell fate choices, and juxtacrine signalling can affect the macroscopic ordering of cell types in self-renewing epithelial sheets, by studying a simple spatial cell fate model with cells being arranged on a 2D lattice. We show in this model that if cells commit to their fate directly upon cell division, macroscopic patches of cells of the same type emerge, if at least a small proportion of divisions are symmetric, except if signalling interactions are laterally inhibiting. In contrast, if cells are first 'licensed' to differentiate, yet retaining the possibility to return to their naive state, macroscopic order only emerges if the signalling strength exceeds a critical threshold: if then the signalling interactions are laterally inducing, macroscopic patches emerge as well. Lateral inhibition, on the other hand, can in that case generate periodic patterns of alternating cell types (checkerboard pattern), yet only if the proportion of symmetric divisions is sufficiently low. These results can be understood theoretically by an analogy to phase transitions in spin systems known from statistical physics.

Behavioral speech tasks have been widely used to understand the mechanisms of speech motor control in typical speakers as well as in various clinical populations. However, determining which neural functions differ between typical speakers and clinical populations based on behavioral data alone is difficult because multiple mechanisms may lead to the same behavioral differences. For example, individuals with cerebellar ataxia (CA) produce atypically large compensatory responses to pitch perturbations in their auditory feedback, compared to typical speakers, but this pattern could have many explanations. Here, computational modeling techniques were used to address this challenge. Bayesian inference was used to fit a state feedback control (SFC) model of voice fundamental frequency (fo) control to the behavioral pitch perturbation responses of speakers with CA and typical speakers. This fitting process resulted in estimates of posterior likelihood distributions for five model parameters (sensory feedback delays, absolute and relative levels of auditory and somatosensory feedback noise, and controller gain), which were compared between the two groups. Results suggest that the speakers with CA may proportionally weight auditory and somatosensory feedback differently from typical speakers. Specifically, the CA group showed a greater relative sensitivity to auditory feedback than the control group. There were also large group differences in the controller gain parameter, suggesting increased motor output responses to target errors in the CA group. These modeling results generate hypotheses about how CA may affect the speech motor system, which could help guide future empirical investigations in CA. This study also demonstrates the overall proof-of-principle of using this Bayesian inference approach to understand behavioral speech data in terms of interpretable parameters of speech motor control models.