PLoS Computational Biology最新文献

Differences in the apparent 1/f component of neural power spectra require correction depending on the underlying neural mechanisms, which remain incompletely understood. Past studies suggest that neuronal spiking produces broadband signals and shapes the spectral trend of invasive macroscopic recordings, but it is unclear to what extent action potentials (APs) influence scalp EEG. Here, we combined biophysical simulations with statistical modelling to examine the amplitude and spectral content of scalp potentials generated by the electric fields from spiking activity. In physiological parameter regimes, we found that APs contribute negligibly to the EEG spectral trend. Consistent with this, comparing our biophysical simulations with previously published data from pharmacologically paralyzed subjects suggested that the EEG spectral trend can be explained by a combination of synaptic timescales and electromyogram contamination. We also modelled rhythmic EEG generation, finding that APs can generate detectable narrowband power between approximately 60 and 1000 Hz, reaching frequencies much faster than would be possible from synaptic currents. Finally, we show that different spectral detrending strategies are required for AP generated oscillations compared to synaptically generated oscillations, suggesting that existing detrending methods for EEG spectra need to be modified for high frequency signals.

Bacteriophage (or 'phage' - viruses that infect and kill bacteria) are increasingly considered as a therapeutic alternative to treat antibiotic-resistant bacterial infections. However, bacteria can evolve resistance to phage, presenting a significant challenge to the near- and long-term success of phage therapeutics. Application of mixtures of multiple phages (i.e., 'cocktails') has been proposed to limit the emergence of phage-resistant bacterial mutants that could lead to therapeutic failure. Here, we combine theory and computational models of in vivo phage therapy to study the efficacy of a phage cocktail, composed of two complementary phages motivated by the example of Pseudomonas aeruginosa facing two phages that exploit different surface receptors, LUZ19v and PAK_P1. As confirmed in a Luria-Delbrück fluctuation test, this motivating example serves as a model for instances where bacteria are extremely unlikely to develop simultaneous resistance mutations against both phages. We then quantify therapeutic outcomes given single- or double-phage treatment models, as a function of phage traits and host immune strength. Building upon prior work showing monophage therapy efficacy in immunocompetent hosts, here we show that phage cocktails comprised of phage targeting independent bacterial receptors can improve treatment outcome in immunocompromised hosts and reduce the chance that pathogens simultaneously evolve resistance against phage combinations. The finding of phage cocktail efficacy is qualitatively robust to differences in virus-bacteria interactions and host immune dynamics. Altogether, the combined use of theory and computational analysis highlights the influence of viral life history traits and receptor complementarity when designing and deploying phage cocktails in immunocompetent and immunocompromised hosts.

Animal behavior spans many timescales, from short, seconds-scale actions to daily rhythms over many hours to life-long changes during aging. To access longer timescales of behavior, we continuously recorded individual Drosophila melanogaster at 100 frames per second for up to 7 days at a time in featureless arenas on sucrose-agarose media. We use the deep learning framework SLEAP to produce a full-body postural dataset for 47 individuals resulting in nearly 2 billion pose instances. We identify stereotyped behaviors such as grooming, proboscis extension, and locomotion and use the resulting ethograms to explore how the flies' behavior varies across time of day and days in the experiment. We find distinct daily patterns in all stereotyped behaviors, adding specific information about trends in different grooming modalities, proboscis extension duration, and locomotion speed to what is known about the D. melanogaster circadian cycle. Using our holistic measurements of behavior, we find that the hour after dawn is a unique time point in the flies' daily pattern of behavior, and that the behavioral composition of this hour tracks well with other indicators of health such as locomotion speed and the fraction of time spend moving vs. resting. The method, data, and analysis presented here give us a new and clearer picture of D. melanogaster behavior across timescales, revealing novel features that hint at unexplored underlying biological mechanisms.

Given the high global seroprevalence of SARS-CoV-2, understanding the risk of reinfection has become increasingly important. Models developed to track trends in reinfection risk should be robust against possible biases arising from imperfect data observation processes. We performed simulation-based validation of an existing catalytic model designed to detect changes in the risk of reinfection by SARS-CoV-2. The catalytic model assumes the risk of reinfection is proportional to observed infections. Validation involved using simulated primary infections, consistent with the number of observed infections in South Africa. To assess the performance of the catalytic model, we simulated reinfection datasets that incorporated different processes that may bias inference, including imperfect observation and mortality. A Bayesian approach was used to fit the model to simulated data, assuming a negative binomial distribution around the expected number of reinfections, and model projections were compared to the simulated data using different magnitudes of change in reinfection risk. We assessed the model's ability to accurately detect changes in reinfection risk when included in the simulations, as well as the occurrence of false positives when reinfection risk remained constant. The model parameters converged in most scenarios leading to model outputs aligning with anticipated outcomes. The model successfully detected changes in the risk of reinfection when such a change was introduced to the data. Low observation probabilities (10%) of both primary- and reinfections resulted in low numbers of observed cases from the simulated data and poor convergence. The model's performance was assessed on simulated data representative of the South African SARS-CoV-2 epidemic, reflecting its timing of waves and outbreak magnitude. Model performance under similar scenarios may be different in settings with smaller epidemics (and therefore smaller numbers of reinfections). Ensuring model parameter convergence is essential to avoid false-positive detection of shifts in reinfection risk. While the model is robust in most scenarios of imperfect observation and mortality, further simulation-based validation for regions experiencing smaller outbreaks is recommended. Caution must be exercised in directly extrapolating results across different epidemiological contexts without additional validation efforts.

Population-based serological surveys are a key tool in epidemiology to characterize the level of population immunity and reconstruct the past circulation of pathogens. A variety of serocatalytic models have been developed to estimate the force of infection (FOI) (i.e., the rate at which susceptible individuals become infected) from age-stratified seroprevalence data. However, few tool currently exists to easily implement, combine, and compare these models. Here, we introduce an R package, Rsero, that implements a series of serocatalytic models and estimates the FOI from age-stratified seroprevalence data using Bayesian methods. The package also contains a series of features to perform model comparison and visualise model fit. We introduce new serocatalytic models of successive outbreaks and extend existing models of seroreversion to any transmission model. The different features of the package are illustrated with simulated and real-life data. We show we can identify the correct epidemiological scenario and recover model parameters in different epidemiological settings. We also show how the package can support serosurvey study design in a variety of epidemic situations. This package provides a standard framework to epidemiologists and modellers to study the dynamics of past pathogen circulation from cross-sectional serological survey data.

Vaccines which can transmit from vaccinated to unvaccinated animals may be especially useful for increasing immunity in hard to reach populations or in populations where achieving high coverage is logistically infeasible. However, gauging the public health utility for future use of such transmissible vaccines and assessing their risk-benefit tradeoff, given their potential for unintended evolution, hinges on accurate estimates of their indirect protective effect. Here, we establish the conditions under which a two-stage randomized field trial can characterize the protective effects of a transmissible vaccine relative to a traditional vaccine. We contrast the sample sizes required to adequately power these trials when the vaccine is weakly and strongly transmissible. We also identify how required sample sizes change based on the characteristics of host ecology such as the overdispersion of the contact structure of the population, as well as the efficacy of the vaccine and timing of vaccination. Our results indicate the range of scenarios where two-stage randomized field trial designs are feasible and appropriate to capture the protective effects of transmissible vaccines. Our estimates identify the protective benefit of using transmissible vaccines compared to traditional vaccines, and thus can be used to weigh against evolutionary risks.

We develop a Bayesian approach, BayesIso, to identify differentially expressed isoforms from RNA-seq data. The approach features a novel joint model of the sample variability and the deferential state of isoforms. Specifically, the within-sample variability and the between-sample variability of each isoform are modeled by a Poisson-Lognormal model and a Gamma-Gamma model, respectively. Using a Bayesian framework, the differential state of each isoform and the model parameters are jointly estimated by a Markov Chain Monte Carlo (MCMC) method. Extensive studies using simulation and real data demonstrate that BayesIso can effectively detect isoforms of less differentially expressed and differential transcripts for genes with multiple isoforms. We applied the approach to breast cancer RNA-seq data and uncovered a unique set of isoforms that form key pathways associated with breast cancer recurrence. First, PI3K/AKT/mTOR signaling and PTEN signaling pathways are identified as being involved in breast cancer development. Further integrated with protein-protein interaction data, pathways of Jak-STAT, mTOR, MAPK and Wnt signaling are revealed in association with breast cancer recurrence. Finally, several pathways are activated in the early recurrence of breast cancer. In tumors that occur early, members of pathways of cellular metabolism and cell cycle (such as CD36 and TOP2A) are upregulated, while immune response genes such as NFATC1 are downregulated.

Personalized cancer drug treatment is emerging as a frontier issue in modern medical research. Considering the genomic differences among cancer patients, determining the most effective drug treatment plan is a complex and crucial task. In response to these challenges, this study introduces the Adaptive Sparse Graph Contrastive Learning Network (ASGCL), an innovative approach to unraveling latent interactions in the complex context of cancer cell lines and drugs. The core of ASGCL is the GraphMorpher module, an innovative component that enhances the input graph structure via strategic node attribute masking and topological pruning. By contrasting the augmented graph with the original input, the model delineates distinct positive and negative sample sets at both node and graph levels. This dual-level contrastive approach significantly amplifies the model's discriminatory prowess in identifying nuanced drug responses. Leveraging a synergistic combination of supervised and contrastive loss, ASGCL accomplishes end-to-end learning of feature representations, substantially outperforming existing methodologies. Comprehensive ablation studies underscore the efficacy of each component, corroborating the model's robustness. Experimental evaluations further illuminate ASGCL's proficiency in predicting drug responses, offering a potent tool for guiding clinical decision-making in cancer therapy.

Persistent homology applied to the activity of grid cells in the Medial Entorhinal Cortex suggests that this activity lies on a toroidal manifold. By analyzing real data and a simple model, we show that neural oscillations play a key role in the appearance of this toroidal topology. To quantitatively monitor how changes in spike trains influence the topology of the data, we first define a robust measure for the degree of toroidality of a dataset. Using this measure, we find that small perturbations ( ~ 100 ms) of spike times have little influence on both the toroidality and the hexagonality of the ratemaps. Jittering spikes by  ~ 100-500 ms, however, destroys the toroidal topology, while still having little impact on grid scores. These critical jittering time scales fall in the range of the periods of oscillations between the theta and eta bands. We thus hypothesized that these oscillatory modulations of neuronal spiking play a key role in the appearance and robustness of toroidal topology and the hexagonal spatial selectivity is not sufficient. We confirmed this hypothesis using a simple model for the activity of grid cells, consisting of an ensemble of independent rate-modulated Poisson processes. When these rates were modulated by oscillations, the network behaved similarly to the real data in exhibiting toroidal topology, even when the position of the fields were perturbed. In the absence of oscillations, this similarity was substantially lower. Furthermore, we find that the experimentally recorded spike trains indeed exhibit temporal modulations at the eta and theta bands, and that the ratio of the power in the eta band to that of the theta band, [Formula: see text], correlates with the critical jittering time at which the toroidal topology disappears.