PLoS Computational Biology最新文献

The rapid advance of large-scale atlas-level single cell RNA sequences and single-cell chromatin accessibility data provide extraordinary avenues to broad and deep insight into complex biological mechanism. Leveraging the datasets and transfering labels from scRNA-seq to scATAC-seq will empower the exploration of single-cell omics data. However, the current label transfer methods have limited performance, largely due to the lower capable of preserving fine-grained cell populations and intrinsic or extrinsic heterogeneity between datasets. Here, we present a robust deep transfer model based graph convolutional network, scTGCN, which achieves versatile performance in preserving biological variation, while achieving integration hundreds of thousands cells in minutes with low memory consumption. We show that scTGCN is powerful to the integration of mouse atlas data and multimodal data generated from APSA-seq and CITE-seq. Thus, scTGCN shows high label transfer accuracy and effectively knowledge transfer across different modalities.

This paper describes and validates an algorithm to solve optimal control problems for agent-based models (ABMs). For a given ABM and a given optimal control problem, the algorithm derives a surrogate model, typically lower-dimensional, in the form of a system of ordinary differential equations (ODEs), solves the control problem for the surrogate model, and then transfers the solution back to the original ABM. It applies to quite general ABMs and offers several options for the ODE structure, depending on what information about the ABM is to be used. There is a broad range of applications for such an algorithm, since ABMs are used widely in the life sciences, such as ecology, epidemiology, and biomedicine and healthcare, areas where optimal control is an important purpose for modeling, such as for medical digital twin technology.

Chronic pain is a wide-spread condition that is debilitating and expensive to manage, costing the United States alone around $600 billion in 2010. In a common symptom of chronic pain called allodynia, non-painful stimuli produce painful responses with highly variable presentations across individuals. While the specific mechanisms remain unclear, allodynia is hypothesized to be caused by the dysregulation of excitatory-inhibitory (E-I) balance in pain-processing neural circuitry in the dorsal horn of the spinal cord. In this work, we analyze biophysically-motivated subcircuit structures that represent common motifs in neural circuits in laminae I-II of the dorsal horn. These circuits are hypothesized to be part of the neural pathways that mediate two different types of allodynia: static and dynamic. We use neural firing rate models to describe the activity of populations of excitatory and inhibitory interneurons within each subcircuit. By accounting for experimentally-observed responses under healthy conditions, we specify model parameters defining populations of subcircuits that yield typical behavior under normal conditions. Then, we implement a sensitivity analysis approach to identify the mechanisms most likely to cause allodynia-producing dysregulation of the subcircuit's E-I signaling. We find that disruption of E-I balance generally occurs either due to downregulation of inhibitory signaling so that excitatory neurons are "released" from inhibitory control, or due to upregulation of excitatory neuron responses so that excitatory neurons "escape" their inhibitory control. Which of these mechanisms is most likely to occur, the subcircuit components involved in the mechanism, and the proportion of subcircuits exhibiting the mechanism can vary depending on the subcircuit structure. These results suggest specific hypotheses about diverse mechanisms that may be most likely responsible for allodynia, thus offering predictions for the high interindividual variability observed in allodynia and identifying targets for further experimental studies on the underlying mechanisms of this chronic pain symptom.

Astrocytes critically shape whole-brain structure and function by forming extensive gap junctional networks that intimately and actively interact with neurons. Despite their importance, existing computational models of whole-brain activity ignore the roles of astrocytes while primarily focusing on neurons. Addressing this oversight, we introduce a biophysical neural mass network model, designed to capture the dynamic interplay between astrocytes and neurons via glutamatergic and GABAergic transmission pathways. This network model proposes that neural dynamics are constrained by a two-layered structural network interconnecting both astrocytic and neuronal populations, allowing us to investigate astrocytes' modulatory influences on whole-brain activity and emerging functional connectivity patterns. By developing a simulation methodology, informed by bifurcation and multilayer network theories, we demonstrate that the dialogue between astrocytic and neuronal networks manifests over fast-slow fluctuation mechanisms as well as through phase-amplitude connectivity processes. The findings from our research represent a significant leap forward in the modeling of glial-neuronal collaboration, promising deeper insights into their collaborative roles across health and disease states.

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm's behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl.

Intracellular liquid-liquid phase separation (LLPS) of proteins and nucleic acids is a fundamental mechanism by which cells compartmentalize their components and perform essential biological functions. Molecular simulations play a crucial role in providing microscopic insights into the physicochemical processes driving this phenomenon. In this study, we systematically compare six state-of-the-art sequence-dependent residue-resolution models to evaluate their performance in reproducing the phase behaviour and material properties of condensates formed by seven variants of the low-complexity domain (LCD) of the hnRNPA1 protein (A1-LCD)-a protein implicated in the pathological liquid-to-solid transition of stress granules. Specifically, we assess the HPS, HPS-cation-π, HPS-Urry, CALVADOS2, Mpipi, and Mpipi-Recharged models in their predictions of the condensate saturation concentration, critical solution temperature, and condensate viscosity of the A1-LCD variants. Our analyses demonstrate that, among the tested models, Mpipi, Mpipi-Recharged, and CALVADOS2 provide accurate descriptions of the critical solution temperatures and saturation concentrations for the multiple A1-LCD variants tested. Regarding the prediction of material properties for condensates of A1-LCD and its variants, Mpipi-Recharged stands out as the most reliable model. Overall, this study benchmarks a range of residue-resolution coarse-grained models for the study of the thermodynamic stability and material properties of condensates and establishes a direct link between their performance and the ranking of intermolecular interactions these models consider.

Respiratory diseases represent a significant healthcare burden, as evidenced by the devastating impact of COVID-19. Biophysical models offer the possibility to anticipate system behavior and provide insights into physiological functions, advancements which are comparatively and notably nascent when it comes to pulmonary mechanics research. In this context, an Inverse Finite Element Analysis (IFEA) pipeline is developed to construct the first continuously ventilated three-dimensional structurally representative pulmonary model informed by both organ- and tissue-level breathing experiments from a cadaveric human lung. Here we construct a generalizable computational framework directly validated by pressure, volume, and strain measurements using a novel inflating apparatus interfaced with adapted, lung-specific, digital image correlation techniques. The parenchyma, pleura, and airways are represented with a poroelastic formulation to simulate pressure flows within the lung lobes, calibrating the model's material properties with the global pressure-volume response and local tissue deformations strains. The optimization yielded the following shear moduli: parenchyma (2.8 kPa), airways (0.2 kPa), and pleura (1.7 Pa). The proposed complex multi-material model with multi-experimental inputs was successfully developed using human lung data, and reproduced the shape of the inflating pressure-volume curve and strain distribution values associated with pulmonary deformation. This advancement marks a significant step towards creating a generalizable human lung model for broad applications across animal models, such as porcine, mouse, and rat lungs to reproduce pathological states and improve performance investigations regarding medical therapeutics and intervention.

Computational models are complex scientific constructs that have become essential for us to better understand the world. Many models are valuable for peers within and beyond disciplinary boundaries. However, there are no widely agreed-upon standards for sharing models. This paper suggests 10 simple rules for you to both (i) ensure you share models in a way that is at least "good enough," and (ii) enable others to lead the change towards better model-sharing practices.

Complex deep learning models trained on very large datasets have become key enabling tools for current research in natural language processing and computer vision. By providing pre-trained models that can be fine-tuned for specific applications, they enable researchers to create accurate models with minimal effort and computational resources. Large scale genomics deep learning models come in two flavors: the first are large language models of DNA sequences trained in a self-supervised fashion, similar to the corresponding natural language models; the second are supervised learning models that leverage large scale genomics datasets from ENCODE and other sources. We argue that these models are the equivalent of foundation models in natural language processing in their utility, as they encode within them chromatin state in its different aspects, providing useful representations that allow quick deployment of accurate models of gene regulation. We demonstrate this premise by leveraging the recently created Sei model to develop simple, interpretable models of intron retention, and demonstrate their advantage over models based on the DNA language model DNABERT-2. Our work also demonstrates the impact of chromatin state on the regulation of intron retention. Using representations learned by Sei, our model is able to discover the involvement of transcription factors and chromatin marks in regulating intron retention, providing better accuracy than a recently published custom model developed for this purpose.

Transfer learning aims to integrate useful information from multi-source datasets to improve the learning performance of target data. This can be effectively applied in genomics when we learn the gene associations in a target tissue, and data from other tissues can be integrated. However, heavy-tail distribution and outliers are common in genomics data, which poses challenges to the effectiveness of current transfer learning approaches. In this paper, we study the transfer learning problem under high-dimensional linear models with t-distributed error (Trans-PtLR), which aims to improve the estimation and prediction of target data by borrowing information from useful source data and offering robustness to accommodate complex data with heavy tails and outliers. In the oracle case with known transferable source datasets, a transfer learning algorithm based on penalized maximum likelihood and expectation-maximization algorithm is established. To avoid including non-informative sources, we propose to select the transferable sources based on cross-validation. Extensive simulation experiments as well as an application demonstrate that Trans-PtLR demonstrates robustness and better performance of estimation and prediction when heavy-tail and outliers exist compared to transfer learning for linear regression model with normal error distribution. Data integration, Variable selection, T distribution, Expectation maximization algorithm, Genotype-Tissue Expression, Cross validation.