PLoS Computational Biology最新文献

Patient-derived tumor organoids (PDTOs) are novel cellular models that maintain the genetic, phenotypic and structural features of patient tumor tissue and are useful for studying tumorigenesis and drug response. When integrated with advanced 3D imaging and analysis techniques, PDTOs can be used to establish physiologically relevant high-throughput and high-content drug screening platforms that support the development of patient-specific treatment strategies. However, in order to effectively leverage high-throughput PDTO observations for clinical predictions, it is critical to establish a quantitative understanding of the basic properties and variability of organoid growth dynamics. In this work, we introduced an innovative workflow for analyzing and understanding PDTO growth dynamics, by integrating a high-throughput imaging deep learning platform with mathematical modeling, incorporating flexible growth laws and variable dormancy times. We applied the workflow to colon cancer organoids and demonstrated that organoid growth is well-described by the Gompertz model of growth. Our analysis showed significant intrapatient heterogeneity in PDTO growth dynamics, with the initial exponential growth rate of an organoid following a lognormal distribution within each dataset. The level of intrapatient heterogeneity varied between patients, as did organoid growth rates and dormancy times of single seeded cells. Our work contributes to an emerging understanding of the basic growth characteristics of PDTOs, and it highlights the heterogeneity in organoid growth both within and between patients. These results pave the way for further modeling efforts aimed at predicting treatment response dynamics and drug resistance timing.

Sounds are temporal stimuli decomposed into numerous elementary components by the auditory nervous system. For instance, a temporal to spectro-temporal transformation modelling the frequency decomposition performed by the cochlea is a widely adopted first processing step in today's computational models of auditory neural responses. Similarly, increments and decrements in sound intensity (i.e., of the raw waveform itself or of its spectral bands) constitute critical features of the neural code, with high behavioural significance. However, despite the growing attention of the scientific community on auditory OFF responses, their relationship with transient ON, sustained responses and adaptation remains unclear. In this context, we propose a new general model, based on a pair of linear filters, named AdapTrans, that captures both sustained and transient ON and OFF responses into a unifying and easy to expand framework. We demonstrate that filtering audio cochleagrams with AdapTrans permits to accurately render known properties of neural responses measured in different mammal species such as the dependence of OFF responses on the stimulus fall time and on the preceding sound duration. Furthermore, by integrating our framework into gold standard and state-of-the-art machine learning models that predict neural responses from audio stimuli, following a supervised training on a large compilation of electrophysiology datasets (ready-to-deploy PyTorch models and pre-processed datasets shared publicly), we show that AdapTrans systematically improves the prediction accuracy of estimated responses within different cortical areas of the rat and ferret auditory brain. Together, these results motivate the use of our framework for computational and systems neuroscientists willing to increase the plausibility and performances of their models of audition.

Budding yeast, Saccharomyces cerevisiae, is widely used as a model organism to study the genetics underlying eukaryotic cellular processes and growth critical to cancer development, such as cell division and cell cycle progression. The budding yeast cell cycle is also one of the best-studied dynamical systems owing to its thoroughly resolved genetics. However, the dynamics underlying the crucial cell cycle decision point called the START transition, at which the cell commits to a new round of DNA replication and cell division, are under-studied. The START machinery involves a central cyclin-dependent kinase; cyclins responsible for starting the transition, bud formation, and initiating DNA synthesis; and their transcriptional regulators. However, evidence has shown that the mechanism is more complicated than a simple irreversible transition switch. Activating a key transcription regulator SBF requires the phosphorylation of its inhibitor, Whi5, or an SBF/MBF monomeric component, Swi6, but not necessarily both. Also, the timing and mechanism of the inhibitor Whi5's nuclear export, while important, are not critical for the timing and execution of START. Therefore, there is a need for a consolidated model for the budding yeast START transition, reconciling regulatory and spatial dynamics. We built a detailed mathematical model (START-BYCC) for the START transition in the budding yeast cell cycle based on established molecular interactions and experimental phenotypes. START-BYCC recapitulates the underlying dynamics and correctly emulates key phenotypic traits of ~150 known START mutants, including regulation of size control, localization of inhibitor/transcription factor complexes, and the nutritional effects on size control. Such a detailed mechanistic understanding of the underlying dynamics gets us closer towards deconvoluting the aberrant cellular development in cancer.

Neural population responses in sensory systems are driven by external physical stimuli. This stimulus-response relationship is typically characterized by receptive fields, which have been estimated by neural system identification approaches. Such models usually require a large amount of training data, yet, the recording time for animal experiments is limited, giving rise to epistemic uncertainty for the learned neural transfer functions. While deep neural network models have demonstrated excellent power on neural prediction, they usually do not provide the uncertainty of the resulting neural representations and derived statistics, such as most exciting inputs (MEIs), from in silico experiments. Here, we present a Bayesian system identification approach to predict neural responses to visual stimuli, and explore whether explicitly modeling network weight variability can be beneficial for identifying neural response properties. To this end, we use variational inference to estimate the posterior distribution of each model weight given the training data. Tests with different neural datasets demonstrate that this method can achieve higher or comparable performance on neural prediction, with a much higher data efficiency compared to Monte Carlo dropout methods and traditional models using point estimates of the model parameters. At the same time, our variational method provides us with an effectively infinite ensemble, avoiding the idiosyncrasy of any single model, to generate MEIs. This allows us to estimate the uncertainty of stimulus-response function, which we have found to be negatively correlated with the predictive performance at model level and may serve to evaluate models. Furthermore, our approach enables us to identify response properties with credible intervals and to determine whether the inferred features are meaningful by performing statistical tests on MEIs. Finally, in silico experiments show that our model generates stimuli driving neuronal activity significantly better than traditional models in the limited-data regime.

Like other tropical and subtropical regions, influenza viruses can circulate year-round in Hong Kong. However, during the COVID-19 pandemic, there was a significant decrease in influenza activity. The objective of this study was to retrospectively forecast influenza activity during the year 2020 and assess the impact of COVID-19 public health social measures (PHSMs) on influenza activity and hospital admissions in Hong Kong. Using weekly surveillance data on influenza virus activity in Hong Kong from 2010 to 2019, we developed a statistical modeling framework to forecast influenza virus activity and associated hospital admissions. We conducted short-term forecasts (1-4 weeks ahead) and medium-term forecasts (1-13 weeks ahead) for the year 2020, assuming no PHSMs were implemented against COVID-19. We estimated the reduction in transmissibility, peak magnitude, attack rates, and influenza-associated hospitalization rate resulting from these PHSMs. For short-term forecasts, mean ambient ozone concentration and school holidays were found to contribute to better prediction performance, while absolute humidity and ozone concentration improved the accuracy of medium-term forecasts. We observed a maximum reduction of 44.6% (95% CI: 38.6% - 51.9%) in transmissibility, 75.5% (95% CI: 73.0% - 77.6%) in attack rate, 41.5% (95% CI: 13.9% - 55.7%) in peak magnitude, and 63.1% (95% CI: 59.3% - 66.3%) in cumulative influenza-associated hospitalizations during the winter-spring period of the 2019/2020 season in Hong Kong. The implementation of PHSMs to control COVID-19 had a substantial impact on influenza transmission and associated burden in Hong Kong. Incorporating information on factors influencing influenza transmission improved the accuracy of our predictions.

The pulsatile activity of gonadotropin-releasing hormone neurons (GnRH neurons) is a key factor in the regulation of reproductive hormones. This pulsatility is orchestrated by a network of neurons that release the neurotransmitters kisspeptin, neurokinin B, and dynorphin (KNDy neurons), and produce episodic bursts of activity driving the GnRH neurons. We show in this computational study that the features of coordinated KNDy neuron activity can be explained by a neural network in which connectivity among neurons is modular. That is, a network structure consisting of clusters of highly-connected neurons with sparse coupling among the clusters. This modular structure, with distinct parameters for intracluster and intercluster coupling, also yields predictions for the differential effects on synchronization of changes in the coupling strength within clusters versus between clusters.

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials.

Melanoma showcases a complex interplay of genetic alterations and intra- and inter-cellular morphological changes during metastatic transformation. While pivotal, the role of specific mutations in dictating these changes still needs to be fully elucidated. Telomerase promoter mutations (TERTp mutations) significantly influence melanoma's progression, invasiveness, and resistance to various emerging treatments, including chemical inhibitors, telomerase inhibitors, targeted therapy, and immunotherapies. We aim to understand the morphological and phenotypic implications of the two dominant monoallelic TERTp mutations, C228T and C250T, enriched in melanoma metastasis. We developed isogenic clonal cell lines containing the TERTp mutations and utilized dual-color expression reporters steered by the endogenous Telomerase promoter, giving us allelic resolution. This approach allowed us to monitor morpholomic variations induced by these mutations. TERTp mutation-bearing cells exhibited significant morpholome differences from their wild-type counterparts, with increased allele expression patterns, augmented wound-healing rates, and unique spatiotemporal dynamics. Notably, the C250T mutation exerted more pronounced changes in the morpholome than C228T, suggesting a differential role in metastatic potential. Our findings underscore the distinct influence of TERTp mutations on melanoma's cellular architecture and behavior. The C250T mutation may offer a unique morpholomic and systems-driven advantage for metastasis. These insights provide a foundational understanding of how a non-coding mutation in melanoma metastasis affects the system, manifesting in cellular morpholome.

Finding communities in gene co-expression networks is a common first step toward extracting biological insight from these complex datasets. Most community detection algorithms expect genes to be organized into assortative modules, that is, groups of genes that are more associated with each other than with genes in other groups. While it is reasonable to expect that these modules exist, using methods that assume they exist a priori is risky, as it guarantees that alternative organizations of gene interactions will be ignored. Here, we ask: can we find meaningful communities without imposing a modular organization on gene co-expression networks, and how modular are these communities? For this, we use a recently developed community detection method, the weighted degree corrected stochastic block model (SBM), that does not assume that assortative modules exist. Instead, the SBM attempts to efficiently use all information contained in the co-expression network to separate the genes into hierarchically organized blocks of genes. Using RNAseq gene expression data measured in two tissues derived from an outbred population of Drosophila melanogaster, we show that (a) the SBM is able to find ten times as many groups as competing methods, that (b) several of those gene groups are not modular, and that (c) the functional enrichment for non-modular groups is as strong as for modular communities. These results show that the transcriptome is structured in more complex ways than traditionally thought and that we should revisit the long-standing assumption that modularity is the main driver of the structuring of gene co-expression networks.

Collective alignment of cell populations is a commonly observed phenomena in biology. An important example are aligning fibroblasts in healthy or scar tissue. In this work we derive and simulate a mechanistic agent-based model of the collective behaviour of actively moving and interacting cells, with a focus on understanding collective alignment. The derivation strategy is based on energy minimisation. The model ingredients are motivated by data on the behaviour of different populations of aligning fibroblasts and include: Self-propulsion, overlap avoidance, deformability, cell-cell junctions and cytoskeletal forces. We find that there is an optimal ratio of self-propulsion speed and overlap avoidance that maximises collective alignment. Further we find that deformability aids alignment, and that cell-cell junctions by themselves hinder alignment. However, if cytoskeletal forces are transmitted via cell-cell junctions we observe strong collective alignment over large spatial scales.