Nonlinear expression patterns and multiple shifts in gene network interactions underlie robust phenotypic change in Drosophila melanogaster selected for night sleep duration.

All but the simplest phenotypes are believed to result from interactions between two or more genes forming complex networks of gene regulation. Sleep is a complex trait known to depend on the system of feedback loops of the circadian clock, and on many other genes; however, the main components regulating the phenotype and how they interact remain an unsolved puzzle. Genomic and transcriptomic data may well provide part of the answer, but a full account requires a suitable quantitative framework. Here we conducted an artificial selection experiment for sleep duration with RNA-seq data acquired each generation. The phenotypic results are robust across replicates and previous experiments, and the transcription data provides a high-resolution, time-course data set for the evolution of sleep-related gene expression. In addition to a Hierarchical Generalized Linear Model analysis of differential expression that accounts for experimental replicates we develop a flexible Gaussian Process model that estimates interactions between genes. 145 gene pairs are found to have interactions that are different from controls. Our method appears to be not only more specific than standard correlation metrics but also more sensitive, finding correlations not significant by other methods. Statistical predictions were compared to experimental data from public databases on gene interactions. Mutations of candidate genes implicated by our results affected night sleep, and gene expression profiles largely met predicted gene-gene interactions.

Individual variation in susceptibility or exposure to SARS-CoV-2 lowers the herd immunity threshold.

Individual variation in susceptibility and exposure is subject to selection by natural infection, accelerating the acquisition of immunity, and reducing herd immunity thresholds and epidemic final sizes. This is a manifestation of a wider population phenomenon known as "frailty variation". Despite theoretical understanding, public health policies continue to be guided by mathematical models that leave out considerable variation and as a result inflate projected disease burdens and overestimate the impact of interventions. Here we focus on trajectories of the coronavirus disease (COVID-19) pandemic in England and Scotland until November 2021. We fit models to series of daily deaths and infer relevant epidemiological parameters, including coefficients of variation and effects of non-pharmaceutical interventions which we find in agreement with independent empirical estimates based on contact surveys. Our estimates are robust to whether the analysed data series encompass one or two pandemic waves and enable projections compatible with subsequent dynamics. We conclude that vaccination programmes may have contributed modestly to the acquisition of herd immunity in populations with high levels of pre-existing naturally acquired immunity, while being crucial to protect vulnerable individuals from severe outcomes as the virus becomes endemic.

Individual variation in susceptibility or exposure to SARS-CoV-2 lowers the herd immunity threshold.

Individual variation in susceptibility and exposure is subject to selection by natural infection, accelerating the acquisition of immunity, and reducing herd immunity thresholds and epidemic final sizes. This is a manifestation of a wider population phenomenon known as "frailty variation". Despite theoretical understanding, public health policies continue to be guided by mathematical models that leave out considerable variation and as a result inflate projected disease burdens and overestimate the impact of interventions. Here we focus on trajectories of the coronavirus disease (COVID-19) pandemic in England and Scotland until November 2021. We fit models to series of daily deaths and infer relevant epidemiological parameters, including coefficients of variation and effects of non-pharmaceutical interventions which we find in agreement with independent empirical estimates based on contact surveys. Our estimates are robust to whether the analysed data series encompass one or two pandemic waves and enable projections compatible with subsequent dynamics. We conclude that vaccination programmes may have contributed modestly to the acquisition of herd immunity in populations with high levels of pre-existing naturally acquired immunity, while being critical to protect vulnerable individuals from severe outcomes as the virus becomes endemic.

Natural selection on sleep duration in Drosophila melanogaster.

Sleep is ubiquitous across animal species, but why it persists is not well understood. Here we observe natural selection act on Drosophila sleep by relaxing bi-directional artificial selection for extreme sleep duration for 62 generations. When artificial selection was suspended, sleep increased in populations previously selected for short sleep. Likewise, sleep decreased in populations previously selected for long sleep when artificial selection was relaxed. We measured the corresponding changes in the allele frequencies of genomic variants responding to artificial selection. The allele frequencies of these variants reversed course in response to relaxed selection, and for short sleepers, the changes exceeded allele frequency changes that would be expected under random genetic drift. These observations suggest that the variants are causal polymorphisms for sleep duration responding to natural selection pressure. These polymorphisms may therefore pinpoint the most important regions of the genome maintaining variation in sleep duration.

Variation in Wolbachia effects on Aedes mosquitoes as a determinant of invasiveness and vectorial capacity.

Wolbachia has been introduced into Aedes aegypti mosquitoes to control the spread of arboviruses, such as dengue, chikungunya and Zika. Studies showed that certain Wolbachia strains (such as wMel) reduce replication of dengue viruses in the laboratory, prompting the release of mosquitoes carrying the bacterium into the field, where vectorial capacity can be realistically assessed in relation to native non-carriers. Here we apply a new analysis to two published datasets, and show that wMel increases the mean and the variance in Ae. aegypti susceptibility to dengue infection when introgressed into Brazil and Vietnam genetic backgrounds. In the absence of other processes, higher mean susceptibility should lead to enhanced viral transmission. The increase in variance, however, widens the basis for selection imposed by unexplored natural forces, retaining the potential for reducing transmission overall.

Model-based inference from multiple dose, time course data reveals Wolbachia effects on infection profiles of type 1 dengue virus in Aedes aegypti.

Infection is a complex and dynamic process involving a population of invading microbes, the host and its responses, aimed at controlling the situation. Depending on the purpose and level of organization, infection at the organism level can be described by a process as simple as a coin toss, or as complex as a multi-factorial dynamic model; the former, for instance, may be adequate as a component of a population model, while the latter is necessary for a thorough description of the process beginning with a challenge with an infectious inoculum up to establishment or elimination of the pathogen. Experimental readouts in the laboratory are often static, snapshots of the process, assayed under some convenient experimental condition, and therefore cannot comprehensively describe the system. Different from the discrete treatment of infection in population models, or the descriptive summarized accounts of typical lab experiments, in this manuscript, infection is treated as a dynamic process dependent on the initial conditions of the infectious challenge, viral growth, and the host response along time. Here, experimental data is generated for multiple doses of type 1 dengue virus, and pathogen levels are recorded at different points in time for two populations of mosquitoes: either carrying endosymbiont bacteria Wolbachia or not. A dynamic microbe/host-response mathematical model is used to describe pathogen growth in the face of a host response like the immune system, and to infer model parameters for the two populations of insects, revealing a slight-but potentially important-protection conferred by the symbiont.

A theoretical framework to identify invariant thresholds in infectious disease epidemiology.

Setting global strategies and targets for disease prevention and control often involves mathematical models. Model structure is typically subject to intense scrutiny, such as confrontation with empirical data and alternative formulations, while a less frequently challenged aspect is the widely adopted reduction of parameters to their average values. Focusing on endemic diseases, we use a general transmission model to explain how mean field approximations decrease the estimated R0 from prevalence data, while threshold phenomena - such as the epidemic and reinfection thresholds - remain invariant. This results in an underestimation of the effort required to control disease, which may be particularly severe when the approximation inappropriately places transmission estimates below important thresholds. These concepts are widely applicable across endemic pathogen systems.

Unveiling time in dose-response models to infer host susceptibility to pathogens.

The biological effects of interventions to control infectious diseases typically depend on the intensity of pathogen challenge. As much as the levels of natural pathogen circulation vary over time and geographical location, the development of invariant efficacy measures is of major importance, even if only indirectly inferrable. Here a method is introduced to assess host susceptibility to pathogens, and applied to a detailed dataset generated by challenging groups of insect hosts (Drosophila melanogaster) with a range of pathogen (Drosophila C Virus) doses and recording survival over time. The experiment was replicated for flies carrying the Wolbachia symbiont, which is known to reduce host susceptibility to viral infections. The entire dataset is fitted by a novel quantitative framework that significantly extends classical methods for microbial risk assessment and provides accurate distributions of symbiont-induced protection. More generally, our data-driven modeling procedure provides novel insights for study design and analyses to assess interventions.