PLoS Computational Biology最新文献

The SARS-CoV-2 pandemic has brought molecular biology and genomic sequencing into the public consciousness and lexicon. With an emphasis on rapid turnaround, genomic data informed both diagnostic and surveillance decisions for the current pandemic at a previously unheard-of scale. The surge in the submission of genomic data to publicly available databases proved essential as comparing different genome sequences offers a wealth of knowledge, including phylogenetic links, modes of transmission, rates of evolution, and the impact of mutations on infection and disease severity. However, the scale of the pandemic has meant that sequencing runs are rarely repeated due to limited sample material and/or the availability of sequencing resources, resulting in the upload of some imperfect runs to public repositories. As a result, it is crucial to investigate the data obtained from these imperfect runs to determine whether the results are reliable prior to depositing them in a public database. Numerous studies have identified a variety of sources of contamination in public next-generation sequencing (NGS) data as the number of NGS studies increases along with the diversity of sequencing technologies and procedures. For this study, we conducted an in silico experiment with known SARS-CoV-2 sequences produced from Oxford Nanopore Technologies sequencing to investigate the effect of contamination on lineage calls and single nucleotide variants (SNVs). A contamination threshold below which runs are expected to generate accurate lineage calls and maintain genome-relatedness and integrity was identified. Together, these findings provide a benchmark below which imperfect runs may be considered robust for reporting results to both stakeholders and public repositories and reduce the need for repeat or wasted runs.

Community organisation permeates both social and biological complex systems. To study its interplay with behaviour emergence, we model mobile structured populations with multiplayer interactions. We derive general analytical methods for evolutionary dynamics under high home fidelity when populations self-organise into networks of asymptotically isolated communities. In this limit, community organisation dominates over the network structure and emerging behaviour is independent of network topology. We obtain the rules of multiplayer cooperation in networks of communities for different types of social dilemmas. The success of cooperation is a result of the benefits shared among communal cooperators outperforming the benefits reaped by defectors in mixed communities. Under weak selection, cooperation can evolve and be stable for any size (Q) and number (M) of communities if the reward-to-cost ratio (V/K) of public goods is higher than a critical value. Community organisation is a solid mechanism for sustaining the evolution of cooperation under public goods dilemmas, particularly when populations are organised into a higher number of smaller communities. Contrary to public goods dilemmas relating to production, the multiplayer Hawk-Dove (HD) dilemma is a commons dilemma focusing on the fair consumption of preexisting resources. This game yields mixed results but tends to favour cooperation under larger communities, highlighting that the two types of social dilemmas might lead to solid differences in the behaviour adopted under community structure.

Human cytomegalovirus (HCMV) is the most common congenital infection. Several HCMV vaccines are in development, but none have yet been approved. An understanding of the kinetics of CMV replication and transmission may inform the rational design of vaccines to prevent this infection. The salivary glands (SG) are an important site of sustained CMV replication following primary infection and during viral reactivation from latency. As such, the strength of the immune response in the SG likely influences viral dissemination within and between hosts. To study the relationship between the immune response and viral replication in the SG, and viral dissemination from the SG to other tissues, mice were infected with low doses of murine CMV (MCMV). Following intra-SG inoculation, we characterized the viral and immunological dynamics in the SG, blood, and spleen, and identified organ-specific immune correlates of protection. Using these data, we constructed compartmental mathematical models of MCMV infection. Model fitting to data and analysis indicate the importance of cellular immune responses in different organs and point to a threshold of infection within the SG necessary for the establishment and spread of infection.

The central complex of insects contains cells, organised as a ring attractor, that encode head direction. The 'bump' of activity in the ring can be updated by idiothetic cues and external sensory information. Plasticity at the synapses between these cells and the ring neurons, that are responsible for bringing sensory information into the central complex, has been proposed to form a mapping between visual cues and the heading estimate which allows for more accurate tracking of the current heading, than if only idiothetic information were used. In Drosophila, ring neurons have well characterised non-linear receptive fields. In this work we produce synthetic versions of these visual receptive fields using a combination of excitatory inputs and mutual inhibition between ring neurons. We use these receptive fields to bring visual information into a spiking neural network model of the insect central complex based on the recently published Drosophila connectome. Previous modelling work has focused on how this circuit functions as a ring attractor using the same type of simple visual cues commonly used experimentally. While we initially test the model on these simple stimuli, we then go on to apply the model to complex natural scenes containing multiple conflicting cues. We show that this simple visual filtering provided by the ring neurons is sufficient to form a mapping between heading and visual features and maintain the heading estimate in the absence of angular velocity input. The network is successful at tracking heading even when presented with videos of natural scenes containing conflicting information from environmental changes and translation of the camera.

Reducing spillover of zoonotic pathogens is an appealing approach to preventing human disease and minimizing the risk of future epidemics and pandemics. Although the immediate human health benefit of reducing spillover is clear, over time, spillover reduction could lead to counterintuitive negative consequences for human health. Here, we use mathematical models and computer simulations to explore the conditions under which unanticipated consequences of spillover reduction can occur in systems where the severity of disease increases with age at infection. Our results demonstrate that, because the average age at infection increases as spillover is reduced, programs that reduce spillover can actually increase population-level disease burden if the clinical severity of infection increases sufficiently rapidly with age. If, however, immunity wanes over time and reinfection is possible, our results reveal that negative health impacts of spillover reduction become substantially less likely. When our model is parameterized using published data on Lassa virus in West Africa, it predicts that negative health outcomes are possible, but likely to be restricted to a small subset of populations where spillover is unusually intense. Together, our results suggest that adverse consequences of spillover reduction programs are unlikely but that the public health gains observed immediately after spillover reduction may fade over time as the age structure of immunity gradually re-equilibrates to a reduced force of infection.

We study the problem of mass screening of heterogeneous populations under limited testing budget. Mass screening is an essential tool that arises in various settings, e.g., the COVID-19 pandemic. The objective of mass screening is to classify the entire population as positive or negative for a disease as efficiently and accurately as possible. Under limited budget, testing facilities need to allocate a portion of the budget to target sub-populations (i.e., proactive screening) while reserving the remaining budget to screen for symptomatic cases (i.e., reactive screening). This paper addresses this decision problem by taking advantage of accessible population-level risk information to identify the optimal set of sub-populations for proactive/reactive screening. The framework also incorporates two widely used testing schemes: Individual and Dorfman group testing. By leveraging the special structure of the resulting bilinear optimization problem, we identify key structural properties, which in turn enable us to develop efficient solution schemes. Furthermore, we extend the model to accommodate customized testing schemes across different sub-populations and introduce a highly efficient heuristic solution algorithm for the generalized model. We conduct a comprehensive case study on COVID-19 in the US, utilizing geographically-based data. Numerical results demonstrate a significant improvement of up to 52% in total misclassifications compared to conventional screening strategies. In addition, our case study offers valuable managerial insights regarding the allocation of proactive/reactive measures and budget across diverse geographic regions.

Translocation in protein synthesis entails the efficient and accurate movement of the mRNA-[tRNA]2 substrate through the ribosome after peptide bond formation. An essential conformational change during this process is the swiveling of the small subunit head domain about two rRNA 'hinge' elements. Using iterative selection and molecular dynamics simulations, we derive alternate hinge elements capable of translocation in vitro and in vivo and describe their effects on the conformational trajectory of the EF-G-bound, translocating ribosome. In these alternate conformational pathways, we observe a diversity of swivel kinetics, hinge motions, three-dimensional head domain trajectories and tRNA dynamics. By finding alternate conformational pathways of translocation, we identify motions and intermediates that are essential or malleable in this process. These findings highlight the plasticity of protein synthesis and provide a more thorough understanding of the available sequence and conformational landscape of a central biological process.

Surprise is a key component of many learning experiences, and yet its precise computational role, and how it changes with age, remain debated. One major challenge is that surprise often occurs jointly with other variables, such as uncertainty and outcome probability. To assess how humans learn from surprising events, and whether aging affects this process, we studied choices while participants learned from bandits with either Gaussian or bi-modal outcome distributions, which decoupled outcome probability, uncertainty, and surprise. A total of 102 participants (51 older, aged 50-73; 51 younger, 19-30 years) chose between three bandits, one of which had a bimodal outcome distribution. Behavioral analyses showed that both age-groups learned the average of the bimodal bandit less well. A trial-by-trial analysis indicated that participants performed choice reversals immediately following large absolute prediction errors, consistent with heightened sensitivity to surprise. This effect was stronger in older adults. Computational models indicated that learning rates in younger as well as older adults were influenced by surprise, rather than uncertainty, but also suggested large interindividual variability in the process underlying learning in our task. Our work bridges between behavioral economics research that has focused on how outcomes with low probability affect choice in older adults, and reinforcement learning work that has investigated age differences in the effects of uncertainty and suggests that older adults overly adapt to surprising events, even when accounting for probability and uncertainty effects.

The experimental study and transplantation of pancreatic islets requires their isolation from the surrounding tissue, and therefore, from the vasculature. Under these conditions, avascular islets rely on the diffusion of peripheral oxygen and nutrients to comply with the requirements of islet cells while responding to changes in body glucose. As a complement to the experimental work, computational models have been widely used to estimate how avascular islets would be affected by the hypoxic conditions found both in culture and transplant sites. However, previous models have been based on simplified representations of pancreatic islets which has limited the reach of the simulations performed. Aiming to contribute with a more realistic model of avascular human islets, in this work we used architectures of human islets reconstructed from experimental data to simulate the availability of oxygen for α, β and δ-cells, emulating culture and transplant conditions at different glucose concentrations. The modeling approach proposed allowed us to quantitatively estimate how the loss of cells due to severe hypoxia would impact interactions between islet cells, ultimately segregating the islet into disconnected subnetworks. According to the simulations performed, islet encapsulation, by reducing the oxygen available within the islets, could severely compromise cell viability. Moreover, our model suggests that even without encapsulation, only microislets composed of less than 100 cells would remain viable in oxygenation conditions found in transplant sites. Overall, in this article we delineate a novel modeling methodology to simulate detailed avascular islets in experimental and transplant conditions with potential applications in the field of islet encapsulation.