Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster

Improvement of starvation resistance via periodic fasting is genetically variable in Drosophila melanogaster

Organisms subjected to periodic nutrient limitation early in life exhibit improvements in aspects of survival, including resistance to some environmental stressors. Recent findings indicate that forms of periodic fasting such as intermittent fasting and time restricted feeding can improve starvation resistance. However, it remains unclear to what extent this survival improvement persists across different genetic backgrounds. In this study, we examine fasting-induced starvation resistance across a broad survey of wild-derived lineages and document genetic variation within this trait. We adopt a standard dietary intervention and show improvement to starvation resistance within a common laboratory lineage, replicating previous results. Next, we examine fasting-induced starvation resistance across isofemale lines collected across latitudes and in different seasons, and among inbred lines derived from flies collected on different continents. We discover genetic variation of fasting-induced starvation resistance, and show that fasting improved starvation resistance as often as it worsened starvation resistance. Fasted flies generally showed reduced fat concentration, and their starvation survival varied with sex, season of collection, and geographic origin. While specific lineages common to the laboratory can show a specific fasting-induced phenotype, we show that this result is not consistent across genetic backgrounds, reinforcing the idea that phenotypes observed in historic laboratory strains may not be conserved across a species.

Drosophila melanogaster pigmentation demonstrates adaptive phenotypic parallelism but genomic unpredictability over multiple timescales

Populations are capable of responding to environmental change over ecological timescales via adaptive tracking. However, the translation from patterns of allele frequency change to rapid adaptation of complex traits remains unresolved. We used abdominal pigmentation in Drosophila melanogaster as a model phenotype to address the nature, genetic architecture, and repeatability of rapid adaptation in the field. We show that D. melanogaster pigmentation evolves as a highly parallel and deterministic response to shared environmental gradients across latitude and season in natural North American populations. We then experimentally evolved replicate, genetically diverse fly populations in field mesocosms to remove any confounding effects of demography and/or cryptic structure that may drive patterns in wild populations; we show that pigmentation rapidly responds, in parallel, in fewer than ten generations. Thus, pigmentation evolves concordantly in response to spatial and temporal climatic gradients. We next examined whether phenotypic differentiation was associated with allele frequency change at loci with established links to genetic variance in pigmentation in natural populations. We found that across all spatial and temporal scales, phenotypic patterns were associated with variation at pigmentation-related loci, and the sets of genes we identified in each context were largely nonoverlapping. Therefore, our findings suggest that parallel phenotypic evolution is associated with an unpredictable genomic response, with distinct components of the polygenic architecture shifting across each environmental gradient to produce redundant adaptive patterns.

Significance Statement

Shifts in global climate conditions have heightened our need to understand the dynamics and pace of adaptation in natural populations. In order to anticipate the population-level response to rapidly changing environmental conditions, we need to understand whether trait evolution is predictable over short timescales, and whether the genetic basis of adaptation is shared or distinct across multiple timescales. Here, we explored parallelism in the adaptive response of a complex phenotype, D. melanogaster pigmentation, to shared conditions that varied over multiple spatiotemporal scales. Our results demonstrate that while phenotypic adaptation proceeds as a predictable response to environmental gradients, even over short timescales, the genetic basis of the adaptive response is variable and nuanced across spatial and temporal contexts.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

Rapid evolutionary change, constraints and the maintenance of polymorphism in natural populations of Drosophila melanogaster

Allele frequencies can shift rapidly within natural populations. Under certain conditions, repeated rapid allele frequency shifts can lead to the long-term maintenance of polymorphism. In recent years, studies of the model insect Drosophila melanogaster have suggested that this phenomenon is more common than previously believed and is often driven by some form of balancing selection, such as temporally fluctuating or sexually antagonistic selection. Here we discuss some of the general insights into rapid evolutionary change revealed by large-scale population genomic studies, as well as the functional and mechanistic causes of rapid adaptation uncovered by single-gene studies. As an example of the latter, we consider a regulatory polymorphism of the D. melanogaster fezzik gene. Polymorphism at this site has been maintained at intermediate frequency over an extended period of time. Regular observations from a single population over a period of 7 years revealed significant differences in the frequency of the derived allele and its variance across collections between the sexes. These patterns are highly unlikely to arise from genetic drift alone or from the action of sexually antagonistic or temporally fluctuating selection individually. Instead, the joint action of sexually antagonistic and temporally fluctuating selection can best explain the observed rapid and repeated allele frequency shifts. Temporal studies such as those reviewed here further our understanding of how rapid changes in selection can lead to the long-term maintenance of polymorphism as well as improve our knowledge of the forces driving and limiting adaptation in nature.

Activation of immune pathways in common bed bugs, Cimex lectularius, in response to bacterial immune challenges - a transcriptomics analysis

The common bed bug, Cimex lectularius, is an urban pest of global health significance, severely affecting the physical and mental health of humans. In contrast to most other blood-feeding arthropods, bed bugs are not major vectors of pathogens, but the underlying mechanisms for this phenomenon are largely unexplored. Here, we present the first transcriptomics study of bed bugs in response to immune challenges. To study transcriptional variations in bed bugs following ingestion of bacteria, we extracted and processed mRNA from body tissues of adult male bed bugs after ingestion of sterile blood or blood containing the Gram-positive (Gr+) bacterium Bacillus subtilis or the Gram-negative (Gr-) bacterium Escherichia coli. We analyzed mRNA from the bed bugs' midgut (the primary tissue involved in blood ingestion) and from the rest of their bodies (RoB; body minus head and midgut tissues). We show that the midgut exhibits a stronger immune response to ingestion of bacteria than the RoB, as indicated by the expression of genes encoding antimicrobial peptides (AMPs). Both the Toll and Imd signaling pathways, associated with immune responses, were highly activated by the ingestion of bacteria. Bacterial infection in bed bugs further provides evidence for metabolic reconfiguration and resource allocation in the bed bugs' midgut and RoB to promote production of AMPs. Our data suggest that infection with particular pathogens in bed bugs may be associated with altered metabolic pathways within the midgut and RoB that favors immune responses. We further show that multiple established cellular immune responses are preserved and are activated by the presence of specific pathogens. Our study provides a greater understanding of nuances in the immune responses of bed bugs towards pathogens that ultimately might contribute to novel bed bug control tactics.

Genetic basis of variation in thermal developmental plasticity for Drosophila melanogaster body pigmentation

Seasonal differences in insect pigmentation are attributed to the influence of ambient temperature on pigmentation development. This thermal plasticity is adaptive and heritable, and thereby capable of evolving. However, the specific genes contributing to the variation in plasticity that can drive its evolution remain largely unknown. To address this, we analysed pigmentation and pigmentation plasticity in Drosophila melanogaster. We measured two components of pigmentation in the thorax and abdomen: overall darkness and the proportion of length covered by darker pattern elements (a trident in the thorax and bands in the abdomen) in females from two developmental temperatures (17 or 28°C) and 191 genotypes. Using a GWAS approach to identify the genetic basis of variation in pigmentation and its response to temperature, we identified numerous dispersed QTLs, including some mapping to melanogenesis genes (yellow, ebony, and tan). Remarkably, we observed limited overlap between QTLs for variation within specific temperatures and those influencing thermal plasticity, as well as minimal overlap between plasticity QTLs across pigmentation components and across body parts. For most traits, consistent with selection favouring the retention of plasticity, we found that lower plasticity alleles were often at lower frequencies. The functional analysis of selected candidate QTLs and pigmentation genes largely confirmed their contributions to variation in pigmentation and/or pigmentation plasticity. Overall, our study reveals the existence and underlying basis of extensive and trait-specific genetic variation for pigmentation and pigmentation plasticity, offering a rich reservoir of raw material for natural selection to shape the evolution of these traits independently.

A cosmopolitan inversion facilitates seasonal adaptation in overwintering Drosophila

Fluctuations in the strength and direction of natural selection through time are a ubiquitous feature of life on Earth. One evolutionary outcome of such fluctuations is adaptive tracking, wherein populations rapidly adapt from standing genetic variation. In certain circumstances, adaptive tracking can lead to the long-term maintenance of functional polymorphism despite allele frequency change due to selection. Although adaptive tracking is likely a common process, we still have a limited understanding of aspects of its genetic architecture and its strength relative to other evolutionary forces such as drift. Drosophila melanogaster living in temperate regions evolve to track seasonal fluctuations and are an excellent system to tackle these gaps in knowledge. By sequencing orchard populations collected across multiple years, we characterized the genomic signal of seasonal demography and identified that the cosmopolitan inversion In(2L)t facilitates seasonal adaptive tracking and shows molecular footprints of selection. A meta-analysis of phenotypic studies shows that seasonal loci within In(2L)t are associated with behavior, life history, physiology, and morphological traits. We identify candidate loci and experimentally link them to phenotype. Our work contributes to our general understanding of fluctuating selection and highlights the evolutionary outcome and dynamics of contemporary selection on inversions.

A Simple Protocol for Isolating Hemolymph from Single Drosophila melanogaster Adult Flies

Drosophila melanogaster is an excellent model for dissecting innate immune signaling and functions. Humoral and cellular immune mechanisms in the fly take place in the hemolymph, where host defense components are secreted and act in response to microbial invaders. Studying hemolymph factors is critical for understanding the regulation of the host's antimicrobial immune system. Therefore, methods for extracting the fly hemolymph efficiently and in sufficient quantities are essential for isolating and characterizing immune proteins and peptides. Here, we describe a novel and simple hemolymph isolation protocol for single D. melanogaster male and female adults. This procedure substantially improves the already used technique and allows fly immunologists to explore innate immune hemolymph activity in D. melanogaster individuals.

Current stewardship practices in invasion biology limit the value and secondary use of genomic data

Invasive species threaten native biota, putting fragile ecosystems at risk and having a large-scale impact on primary industries. Growing trade networks and the popularity of personal travel make incursions a more frequent risk, one only compounded by global climate change. With increasing publication of whole-genome sequences lies an opportunity for cross-species assessment of invasive potential. However, the degree to which published sequences are accompanied by satisfactory spatiotemporal data is unclear. We assessed the metadata associated with 199 whole-genome assemblies of 89 invasive terrestrial invertebrate species and found that only 38% of these were derived from field-collected samples. Seventy-six assemblies (38%) reported an 'undescribed' sample origin and, while further examination of associated literature closed this gap to 23.6%, an absence of spatial data remained for 47 of the total assemblies. Of the 76 assemblies that were ultimately determined to be field-collected, associated metadata relevant for invasion studies was predominantly lacking: only 35% (27 assemblies) provided granular location data, and 33% (n = 25) lacked sufficient collection date information. Our results support recent calls for standardized metadata in genome sequencing data submissions, highlighting the impact of missing metadata on current research in invasion biology (and likely other fields). Notably, large-scale consortia tended to provide the most complete metadata submissions in our analysis-such cross-institutional collaborations can foster a culture of increased adherence to improved metadata submission standards and a standard of metadata stewardship that enables reuse of genomes in invasion science.

Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes

Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.

The genomic basis of copper tolerance in Drosophila is shaped by a complex interplay of regulatory and environmental factors

BACKGROUND

Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors.

RESULTS

In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response.

CONCLUSIONS

Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.

Genome-wide selection signatures reveal widespread synergistic effects of two different stressors in Drosophila melanogaster

Experimental evolution combined with whole-genome sequencing (evolve and resequence (E&R)) is a powerful approach to study the adaptive architecture of selected traits. Nevertheless, so far the focus has been on the selective response triggered by a single stressor. Building on the highly parallel selection response of founder populations with reduced variation, we evaluated how the presence of a second stressor affects the genomic selection response. After 20 generations of adaptation to laboratory conditions at either 18°C or 29°C, strong genome-wide selection signatures were observed. Only 38% of the selection signatures can be attributed to laboratory adaptation (no difference between temperature regimes). The remaining selection responses are either caused by temperature-specific effects, or reflect the joint effects of temperature and laboratory adaptation (same direction, but the magnitude differs between temperatures). The allele frequency changes resulting from the combined effects of temperature and laboratory adaptation were more extreme in the hot environment for 83% of the affected genomic regions-indicating widespread synergistic effects of the two stressors. We conclude that E&R with reduced genetic variation is a powerful approach to study genome-wide fitness consequences driven by the combined effects of multiple environmental factors.

Seasonal variation in wing size and shape of Drosophila melanogaster reveals rapid adaptation to environmental changes

Populations in seasonal fluctuating environments receive multiple environmental cues and must deal with this heterogenic environment to survive and reproduce. An enlarged literature shows that this situation can be resolved through rapid adaptation in Drosophila melanogaster populations. Long-term monitoring of a population in its natural habitat and quantitative measurement of its responses to seasonal environmental changes are important for understanding the adaptive response of D. melanogaster to temporal variable selection. Here, we use inbred lines of a D. melanogaster population collected at monthly intervals between May to October over a temporal scale spanning three consecutive years to understand the variation in wing size and wing shape over these timepoints. The wing size and shape of this population changed significantly between months and a seasonal cycle of this traits is repeated for three years. Our results suggest that the effects of environmental variables that generated variation in body size between populations such as latitudinal clines, are a selective pressure in a different manner in terms of seasonal variation. Temperature related variable have a significant nonlinear relation to this fluctuating pattern in size and shape, whereas precipitation and humidity have a sex-specific effect which is more significant in males.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Adaptive Evolution in TRIF Leads to Discordance between Human and Mouse Innate Immune Signaling

The TIR domain-containing adapter inducing IFN-β (TRIF) protein is an innate immune system protein that mediates the MyD88-independent toll-like receptor response pathway in mice and humans. Previously, we identified positive selection at seven distinct residues in mouse TRIF (mTRIF), as compared with human and other mammalian orthologs, thus predicting protein functional shift in mTRIF. We reconstructed TRIF for the most recent common ancestor of mouse and human, and mutated this at the seven sites to their extant mouse/human states. We overexpressed these TRIF mutants in immortalized human and mouse cell lines and monitored TRIF-dependent cytokine production and gene expression induction. We show that optimal TRIF function in human and mouse is dependent on the identity of the positively selected sites. These data provide us with molecular data relating observed differences in response between mouse and human MyD88-independent signaling in the innate immune system with protein functional change.

easyFulcrum: An R package to process and analyze ecological sampling data generated using the Fulcrum mobile application

Large-scale ecological sampling can be difficult and costly, especially for organisms that are too small to be easily identified in a natural environment by eye. Typically, these microscopic floral and fauna are sampled by collecting substrates from nature and then separating organisms from substrates in the laboratory. In many cases, diverse organisms can be identified to the species-level using molecular barcodes. To facilitate large-scale ecological sampling of microscopic organisms, we used a geographic data-collection platform for mobile devices called Fulcrum that streamlines the organization of geospatial sampling data, substrate photographs, and environmental data at natural sampling sites. These sampling data are then linked to organism isolation data from the laboratory. Here, we describe the easyFulcrum R package, which can be used to clean, process, and visualize ecological field sampling and isolation data exported from the Fulcrum mobile application. We developed this package for wild nematode sampling, but it can be used with other organisms. The advantages of using Fulcrum combined with easyFulcrum are (1) the elimination of transcription errors by replacing manual data entry and/or spreadsheets with a mobile application, (2) the ability to clean, process, and visualize sampling data using a standardized set of functions in the R software environment, and (3) the ability to join disparate data to each other, including environmental data from the field and the molecularly defined identities of individual specimens isolated from samples.

Horizontal gene transfer-mediated bacterial strain variation affects host fitness in Drosophila

Background

How microbes affect host fitness and environmental adaptation has become a fundamental research question in evolutionary biology. To better understand the role of microbial genomic variation for host fitness, we tested for associations of bacterial genomic variation and Drosophila melanogaster offspring number in a microbial Genome Wide Association Study (GWAS).

Results

We performed a microbial GWAS, leveraging strain variation in the genus Gluconobacter, a genus of bacteria that are commonly associated with Drosophila under natural conditions. We pinpoint the thiamine biosynthesis pathway (TBP) as contributing to differences in fitness conferred to the fly host. While an effect of thiamine on fly development has been described, we show that strain variation in TBP between bacterial isolates from wild-caught D. melanogaster contributes to variation in offspring production by the host. By tracing the evolutionary history of TBP genes in Gluconobacter, we find that TBP genes were most likely lost and reacquired by horizontal gene transfer (HGT).

Conclusion

Our study emphasizes the importance of strain variation and highlights that HGT can add to microbiome flexibility and potentially to host adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01124-y.

Frequent Drivers, Occasional Passengers: Signals of Symbiont-Driven Seasonal Adaptation and Hitchhiking in the Pea Aphid, Acyrthosiphon pisum

Insects harbor a variety of maternally inherited bacterial symbionts. As such, variation in symbiont presence/absence, in the combinations of harbored symbionts, and in the genotypes of harbored symbiont species provide heritable genetic variation of potential use in the insects' adaptive repertoires. Understanding the natural importance of symbionts is challenging but studying their dynamics over time can help to elucidate the potential for such symbiont-driven insect adaptation. Toward this end, we studied the seasonal dynamics of six maternally transferred bacterial symbiont species in the multivoltine pea aphid (Acyrthosiphon pisum). Our sampling focused on six alfalfa fields in southeastern Pennsylvania, and spanned 14 timepoints within the 2012 growing season, in addition to two overwintering periods. To test and generate hypotheses on the natural relevance of these non-essential symbionts, we examined whether symbiont dynamics correlated with any of ten measured environmental variables from the 2012 growing season, including some of known importance in the lab. We found that five symbionts changed prevalence across one or both overwintering periods, and that the same five species underwent such frequency shifts across the 2012 growing season. Intriguingly, the frequencies of these dynamic symbionts showed robust correlations with a subset of our measured environmental variables. Several of these trends supported the natural relevance of lab-discovered symbiont roles, including anti-pathogen defense. For a seventh symbiont-Hamiltonella defensa-studied previously across the same study periods, we tested whether a reported correlation between prevalence and temperature stemmed not from thermally varying host-level fitness effects, but from selection on co-infecting symbionts or on aphid-encoded alleles associated with this bacterium. In general, such "hitchhiking" effects were not evident during times with strongly correlated Hamiltonella and temperature shifts. However, we did identify at least one time period in which Hamiltonella spread was likely driven by selection on a co-infecting symbiont-Rickettsiella viridis. Recognizing the broader potential for such hitchhiking, we explored selection on co-infecting symbionts as a possible driver behind the dynamics of the remaining six species. Out of twelve examined instances of symbiont dynamics unfolding across 2-week periods or overwintering spans, we found eight in which the focal symbiont underwent parallel frequency shifts under single infection and one or more co-infection contexts. This supported the idea that phenotypic variation created by the presence/absence of individual symbionts is a direct target for selection, and that symbiont effects can be robust under co-habitation with other symbionts. Contrastingly, in two cases, we found that selection may target phenotypes emerging from symbiont co-infections, with specific species combinations driving overall trends for the focal dynamic symbionts, without correlated change under single infection. Finally, in three cases-including the one described above for Hamiltonella-our data suggested that incidental co-infection with a (dis)favored symbiont could lead to large frequency shifts for "passenger" symbionts, conferring no apparent cost or benefit. Such hitchhiking has rarely been studied in heritable symbiont systems. We propose that it is more common than appreciated, given the widespread nature of maternally inherited bacteria, and the frequency of multi-species symbiotic communities across insects.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

Genomic Responses to Climate Change: Making the Most of the Drosophila Model

It is pressing to understand how animal populations evolve in response to climate change. We argue that new sequencing technologies and the use of historical samples are opening unprecedented opportunities to investigate genome-wide responses to changing environments. However, there are important challenges in interpreting the emerging findings. First, it is essential to differentiate genetic adaptation from phenotypic plasticity. Second, it is extremely difficult to map genotype, phenotype, and fitness. Third, neutral demographic processes and natural selection affect genetic variation in similar ways. We argue that Drosophila melanogaster, a classical model organism with decades of climate adaptation research, is uniquely suited to overcome most of these challenges. In the near future, long-term time series genome-wide datasets of D. melanogaster natural populations will provide exciting opportunities to study adaptation to recent climate change and will lay the groundwork for related research in non-model systems.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

Anticipated effects of abiotic environmental change on intraspecific social interactions

Social interactions are ubiquitous across the animal kingdom. A variety of ecological and evolutionary processes are dependent on social interactions, such as movement, disease spread, information transmission, and density-dependent reproduction and survival. Social interactions, like any behaviour, are context dependent, varying with environmental conditions. Currently, environments are changing rapidly across multiple dimensions, becoming warmer and more variable, while habitats are increasingly fragmented and contaminated with pollutants. Social interactions are expected to change in response to these stressors and to continue to change into the future. However, a comprehensive understanding of the form and magnitude of the effects of these environmental changes on social interactions is currently lacking. Focusing on four major forms of rapid environmental change currently occurring, we review how these changing environmental gradients are expected to have immediate effects on social interactions such as communication, agonistic behaviours, and group formation, which will thereby induce changes in social organisation including mating systems, dominance hierarchies, and collective behaviour. Our review covers intraspecific variation in social interactions across environments, including studies in both the wild and in laboratory settings, and across a range of taxa. The expected responses of social behaviour to environmental change are diverse, but we identify several general themes. First, very dry, variable, fragmented, or polluted environments are likely to destabilise existing social systems. This occurs as these conditions limit the energy available for complex social interactions and affect dissimilar phenotypes differently. Second, a given environmental change can lead to opposite responses in social behaviour, and the direction of the response often hinges on the natural history of the organism in question. Third, our review highlights the fact that changes in environmental factors are not occurring in isolation: multiple factors are changing simultaneously, which may have antagonistic or synergistic effects, and more work should be done to understand these combined effects. We close by identifying methodological and analytical techniques that might help to study the response of social interactions to changing environments, highlight consistent patterns among taxa, and predict subsequent evolutionary change. We expect that the changes in social interactions that we document here will have consequences for individuals, groups, and for the ecology and evolution of populations, and therefore warrant a central place in the study of animal populations, particularly in an era of rapid environmental change.

Broad geographic sampling reveals the shared basis and environmental correlates of seasonal adaptation in Drosophila

To advance our understanding of adaptation to temporally varying selection pressures, we identified signatures of seasonal adaptation occurring in parallel among Drosophila melanogaster populations. Specifically, we estimated allele frequencies genome-wide from flies sampled early and late in the growing season from 20 widely dispersed populations. We identified parallel seasonal allele frequency shifts across North America and Europe, demonstrating that seasonal adaptation is a general phenomenon of temperate fly populations. Seasonally fluctuating polymorphisms are enriched in large chromosomal inversions, and we find a broad concordance between seasonal and spatial allele frequency change. The direction of allele frequency change at seasonally variable polymorphisms can be predicted by weather conditions in the weeks prior to sampling, linking the environment and the genomic response to selection. Our results suggest that fluctuating selection is an important evolutionary force affecting patterns of genetic variation in Drosophila.

Allelic polymorphism at foxo contributes to local adaptation in Drosophila melanogaster

The insulin/insulin-like growth factor signalling pathway has been hypothesized as a major determinant of life-history profiles that vary adaptively in natural populations. In Drosophila melanogaster, multiple components of this pathway vary predictably with latitude; this includes foxo, a conserved gene that regulates insulin signalling and has pleiotropic effects on a variety of fitness-associated traits. We hypothesized that allelic variation at foxo contributes to genetic variance for size-related traits that vary adaptively with latitude. We first examined patterns of variation among natural populations along a latitudinal transect in the eastern United States and show that thorax length, wing area, wing loading, and starvation tolerance exhibit significant latitudinal clines for both males and females but that development time does not vary predictably with latitude. We then generated recombinant outbred populations and show that naturally occurring allelic variation at foxo, which exhibits stronger clinality than expected, is associated with the same traits that vary with latitude in the natural populations. Our results suggest that allelic variation at foxo contributes to adaptive patterns of life-history variation in natural populations of this genetic model.

Does getting defensive get you anywhere?-Seasonal balancing selection, temperature, and parasitoids shape real-world, protective endosymbiont dynamics in the pea aphid

Facultative, heritable endosymbionts are found at intermediate prevalence within most insect species, playing frequent roles in their hosts' defence against environmental pressures. Focusing on Hamiltonella defensa, a common bacterial endosymbiont of aphids, we tested the hypothesis that such pressures impose seasonal balancing selection, shaping a widespread infection polymorphism. In our studied pea aphid (Acyrthosiphon pisum) population, Hamiltonella frequencies ranged from 23.2% to 68.1% across a six-month longitudinal survey. Rapid spikes and declines were often consistent across fields, and we estimated that selection coefficients for Hamiltonella-infected aphids changed sign within this field season. Prior laboratory research suggested antiparasitoid defence as the major Hamiltonella benefit, and costs under parasitoid absence. While a prior field study suggested these forces can sometimes act as counter-weights in a regime of seasonal balancing selection, our present survey showed no significant relationship between parasitoid wasps and Hamiltonella prevalence. Field cage experiments provided some explanation: parasitoids drove modest ~10% boosts to Hamiltonella frequencies that would be hard to detect under less controlled conditions. They also showed that Hamiltonella was not always costly under parasitoid exclusion, contradicting another prediction. Instead, our longitudinal survey - and two overwintering studies - showed temperature to be the strongest predictor of Hamiltonella prevalence. Matching some prior lab discoveries, this suggested that thermally sensitive costs and benefits, unrelated to parasitism, can shape Hamiltonella dynamics. These results add to a growing body of evidence for rapid, seasonal adaptation in multivoltine organisms, suggesting that such adaptation can be mediated through the diverse impacts of heritable bacterial endosymbionts.

Using seasonal genomic changes to understand historical adaptation to new environments: Parallel selection on stickleback in highly-variable estuaries

Parallel evolution is considered strong evidence for natural selection. However, few studies have investigated the process of parallel selection as it plays out in real time. The common approach is to study historical signatures of selection in populations already well adapted to different environments. Here, to document selection under natural conditions, we study six populations of threespine stickleback (Gasterosteus aculeatus) inhabiting bar-built estuaries that undergo seasonal cycles of environmental changes. Estuaries are periodically isolated from the ocean due to sandbar formation during dry summer months, with concurrent environmental shifts that resemble the long-term changes associated with postglacial colonization of freshwater habitats by marine populations. We used pooled whole-genome sequencing to track seasonal allele frequency changes in six of these populations and search for signatures of natural selection. We found consistent changes in allele frequency across estuaries, suggesting a potential role for parallel selection. Functional enrichment among candidate genes included transmembrane ion transport and calcium binding, which are important for osmoregulation and ion balance. The genomic changes that occur in threespine stickleback from bar-built estuaries could provide a glimpse into the early stages of adaptation that have occurred in many historical marine to freshwater transitions.

Competitive history shapes rapid evolution in a seasonal climate

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

Population Genomics on the Fly: Recent Advances in Drosophila

Drosophila melanogaster, a small dipteran of African origin, represents one of the best-studied model organisms. Early work in this system has uniquely shed light on the basic principles of genetics and resulted in a versatile collection of genetic tools that allow to uncover mechanistic links between genotype and phenotype. Moreover, given its worldwide distribution in diverse habitats and its moderate genome-size, Drosophila has proven very powerful for population genetics inference and was one of the first eukaryotes whose genome was fully sequenced. In this book chapter, we provide a brief historical overview of research in Drosophila and then focus on recent advances during the genomic era. After describing different types and sources of genomic data, we discuss mechanisms of neutral evolution including the demographic history of Drosophila and the effects of recombination and biased gene conversion. Then, we review recent advances in detecting genome-wide signals of selection, such as soft and hard selective sweeps. We further provide a brief introduction to background selection, selection of noncoding DNA and codon usage and focus on the role of structural variants, such as transposable elements and chromosomal inversions, during the adaptive process. Finally, we discuss how genomic data helps to dissect neutral and adaptive evolutionary mechanisms that shape genetic and phenotypic variation in natural populations along environmental gradients. In summary, this book chapter serves as a starting point to Drosophila population genomics and provides an introduction to the system and an overview to data sources, important population genetic concepts and recent advances in the field.

The Combined Analysis of Pleiotropy and Epistasis (CAPE)

Epistasis, or gene-gene interaction, contributes substantially to trait variation in organisms ranging from yeast to humans, and modeling epistasis directly is critical to understanding the genotype-phenotype map. However, inference of genetic interactions is challenging compared to inference of individual allele effects due to low statistical power. Furthermore, genetic interactions can appear inconsistent across different quantitative traits, presenting a challenge for the interpretation of detected interactions. Here we present a method called the Combined Analysis of Pleiotropy and Epistasis (CAPE) that combines information across multiple quantitative traits to infer directed epistatic interactions. By combining information across multiple traits, CAPE not only increases power to detect genetic interactions but also interprets these interactions across traits to identify a single interaction that is consistent across all observed data. This method generates informative, interpretable interaction networks that explain how variants interact with each other to influence groups of related traits. This method could potentially be used to link genetic variants to gene expression, physiological endophenotypes, and higher-level disease traits.

Unique genetic signatures of local adaptation over space and time for diapause, an ecologically relevant complex trait, in Drosophila melanogaster

Organisms living in seasonally variable environments utilize cues such as light and temperature to induce plastic responses, enabling them to exploit favorable seasons and avoid unfavorable ones. Local adapation can result in variation in seasonal responses, but the genetic basis and evolutionary history of this variation remains elusive. Many insects, including Drosophila melanogaster, are able to undergo an arrest of reproductive development (diapause) in response to unfavorable conditions. In D. melanogaster, the ability to diapause is more common in high latitude populations, where flies endure harsher winters, and in the spring, reflecting differential survivorship of overwintering populations. Using a novel hybrid swarm-based genome wide association study, we examined the genetic basis and evolutionary history of ovarian diapause. We exposed outbred females to different temperatures and day lengths, characterized ovarian development for over 2800 flies, and reconstructed their complete, phased genomes. We found that diapause, scored at two different developmental cutoffs, has modest heritability, and we identified hundreds of SNPs associated with each of the two phenotypes. Alleles associated with one of the diapause phenotypes tend to be more common at higher latitudes, but these alleles do not show predictable seasonal variation. The collective signal of many small-effect, clinally varying SNPs can plausibly explain latitudinal variation in diapause seen in North America. Alleles associated with diapause are segregating in Zambia, suggesting that variation in diapause relies on ancestral polymorphisms, and both pro- and anti-diapause alleles have experienced selection in North America. Finally, we utilized outdoor mesocosms to track diapause under natural conditions. We found that hybrid swarms reared outdoors evolved increased propensity for diapause in late fall, whereas indoor control populations experienced no such change. Our results indicate that diapause is a complex, quantitative trait with different evolutionary patterns across time and space.

In Vitro and In Vivo Models for Evaluating the Oral Toxicity of Nanomedicines

Toxicity studies for conventional oral drug formulations are standardized and well documented, as required by the guidelines of administrative agencies such as the US Food & Drug Administration (FDA), the European Medicines Agency (EMA) or European Medicines Evaluation Agency (EMEA), and the Japanese Pharmaceuticals and Medical Devices Agency (PMDA). Researchers tend to extrapolate these standardized protocols to evaluate nanoformulations (NFs) because standard nanotoxicity protocols are still lacking in nonclinical studies for testing orally delivered NFs. However, such strategies have generated many inconsistent results because they do not account for the specific physicochemical properties of nanomedicines. Due to their tiny size, accumulated surface charge and tension, sizeable surface-area-to-volume ratio, and high chemical/structural complexity, orally delivered NFs may generate severe topical toxicities to the gastrointestinal tract and metabolic organs, including the liver and kidney. Such toxicities involve immune responses that reflect different mechanisms than those triggered by conventional formulations. Herein, we briefly analyze the potential oral toxicity mechanisms of NFs and describe recently reported in vitro and in vivo models that attempt to address the specific oral toxicity of nanomedicines. We also discuss approaches that may be used to develop nontoxic NFs for oral drug delivery.

The Cyclically Seasonal Drosophila subobscura Inversion O7 Originated From Fragile Genomic Sites and Relocated Immunity and Metabolic Genes

Chromosome inversions are important contributors to standing genetic variation in Drosophila subobscura. Presently, the species is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. Yet not all low-latitude inversions are correlated with the ongoing warming trend. This is particularly unexpected in the case of O7 because it shows a regular seasonal cycle that peaks in summer and rose with a heatwave. The inconsistent behavior of O7 across components of the ambient temperature suggests that is causally more complex than simply due to temperature alone. In order to understand the dynamics of O7, high-quality genomic data are needed to determine both the breakpoints and the genetic content. To fill this gap, here we generated a PacBio long read-based chromosome-scale genome assembly, from a highly homozygous line made isogenic for an O3 + 4 + 7 chromosome. Then we isolated the complete continuous sequence of O7 by conserved synteny analysis with the available reference genome. Main findings include the following: (i) the assembled O7 inversion stretches 9.936 Mb, containing > 1,000 annotated genes; (ii) O7 had a complex origin, involving multiple breaks associated with non-B DNA-forming motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes; (iii) the O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene; and (iv) the O7 inversion relocated the major insulin signaling forkhead box subgroup O (foxo) gene in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1) and disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. Our findings suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Furthermore, they are relevant for fragility in genome rearrangement evolution and for current views on the contribution of breakage versus repair in shaping inversion-breakpoint junctions.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Humic-acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding

Background

Next generation sequencing (NGS) technologies are extensively used to dissect the molecular mechanisms of host-parasite interactions in human pathogens. However, ecological studies have yet to fully exploit the power of NGS as a rich source for formulating and testing new hypotheses.

Methods

We studied Eurasian perch (Perca fluviatilis) and its eye parasite (Trematoda, Diplostomidae) communities in 14 lakes that differed in humic content in order to explore host-parasite-environment interactions. We hypothesised that high humic content along with low pH would decrease the abundance of the intermediate hosts (gastropods), thus limiting the occurrence of diplostomid parasites in humic lakes. This hypothesis was initially invoked by whole eye RNA-seq data analysis and subsequently tested using PCR-based detection and a novel targeted metabarcoding approach.

Results

Whole eye transcriptome results revealed overexpression of immune-related genes and the presence of eye parasite sequences in RNA-seq data obtained from perch living in clear-water lakes. Both PCR-based and targeted-metabarcoding approach showed that perch from humic lakes were completely free from diplostomid parasites, while the prevalence of eye flukes in clear-water lakes that contain low amounts of humic substances was close to 100%, with the majority of NGS reads assigned to Tylodelphys clavata.

Conclusions

High intraspecific diversity of T. clavata indicates that massively parallel sequencing of naturally pooled samples represents an efficient and powerful strategy for shedding light on cryptic diversity of eye parasites. Our results demonstrate that perch populations in clear-water lakes experience contrasting eye parasite pressure compared to those from humic lakes, which is reflected by prevalent differences in the expression of immune-related genes in the eye. This study highlights the utility of NGS to discover novel host-parasite-environment interactions and provide unprecedented power to characterize the molecular diversity of cryptic parasites.

Prior exposure to long-day photoperiods alters immune responses and increases susceptibility to parasitic infection in stickleback

Seasonal disease and parasitic infection are common across organisms, including humans, and there is increasing evidence for intrinsic seasonal variation in immune systems. Changes are orchestrated through organisms' physiological clocks using cues such as day length. Ample research in diverse taxa has demonstrated multiple immune responses are modulated by photoperiod, but to date, there have been few experimental demonstrations that photoperiod cues alter susceptibility to infection. We investigated the interactions among photoperiod history, immunity and susceptibility in laboratory-bred three-spined stickleback (a long-day breeding fish) and its external, directly reproducing monogenean parasite Gyrodactylus gasterostei. We demonstrate that previous exposure to long-day photoperiods (PLD) increases susceptibility to infection relative to previous exposure to short days (PSD), and modifies the response to infection for the mucin gene muc2 and Treg cytokine foxp3a in skin tissues in an intermediate 12 L : 12 D photoperiod experimental trial. Expression of skin muc2 is reduced in PLD fish, and negatively associated with parasite abundance. We also observe inflammatory gene expression variation associated with natural inter-population variation in resistance, but find that photoperiod modulation of susceptibility is consistent across host populations. Thus, photoperiod modulation of the response to infection is important for host susceptibility, highlighting new mechanisms affecting seasonality of host-parasite interactions.

Common structuring principles of the Drosophila melanogaster microbiome on a continental scale and between host and substrate

The relative importance of host control, environmental effects and stochasticity in the assemblage of host-associated microbiomes is being debated. We analysed the microbiome among fly populations that were sampled across Europe by the European Drosophila Population Genomics Consortium (DrosEU). In order to better understand the structuring principles of the natural D. melanogaster microbiome, we combined environmental data on climate and food-substrate with dense genomic data on host populations and microbiome profiling. Food-substrate, temperature, and host population structure correlated with microbiome structure. Microbes, whose abundance was co-structured with host populations, also differed in abundance between flies and their substrate in an independent survey. This finding suggests common, host-related structuring principles of the microbiome on different spatial scales.

Speculations on the evolution of humoral adaptive immunity

The protection of a multicellular organism from infection, at both cell and humoral levels, has been a tremendous driver of gene selection and cellular response strategies. Here we focus on a critical event in the development of humoral immunity: The transition from principally innate responses to a system of adaptive cell selection, with all the attendant mechanical problems that must be solved in order for it to work effectively. Here we review recent advances, but our major goal is to highlight that the development of adaptive immunity resulted from the adoption, reuse and repurposing of an ancient, autonomous cellular program that combines and exploits three titratable cellular fate timers. We illustrate how this common cell machinery recurs and appears throughout biology, and has been essential for the evolution of complex organisms, at many levels of scale.

Solution of the Crow-Kimura model with a periodically changing (two-season) fitness function

Since the origin of life, both evolutionary dynamics and rhythms have played a key role in the functioning of living systems. The Crow-Kimura model of periodically changing fitness function has been solved exactly, using integral equation with time-ordered exponent. We also found a simple approximate solution for the two-season case. The evolutionary dynamics accompanied by the rhythms provide important insights into the properties of certain biological systems and processes.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.

Phenotypic plasticity, but not adaptive tracking, underlies seasonal variation in post-cold hardening freeze tolerance of Drosophila melanogaster

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short-lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post-cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm-caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post-cold hardening freeze tolerance. Comparing indoor and field-caught flies and their laboratory-reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field- and laboratory-based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post-cold hardening freeze tolerance results from environmental influences and not genetic changes.

Accurate Allele Frequencies from Ultra-low Coverage Pool-Seq Samples in Evolve-and-Resequence Experiments

Evolve-and-resequence (E+R) experiments leverage next-generation sequencing technology to track the allele frequency dynamics of populations as they evolve. While previous work has shown that adaptive alleles can be detected by comparing frequency trajectories from many replicate populations, this power comes at the expense of high-coverage (>100x) sequencing of many pooled samples, which can be cost-prohibitive. Here, we show that accurate estimates of allele frequencies can be achieved with very shallow sequencing depths (<5x) via inference of known founder haplotypes in small genomic windows. This technique can be used to efficiently estimate frequencies for any number of bi-allelic SNPs in populations of any model organism founded with sequenced homozygous strains. Using both experimentally-pooled and simulated samples of Drosophila melanogaster, we show that haplotype inference can improve allele frequency accuracy by orders of magnitude for up to 50 generations of recombination, and is robust to moderate levels of missing data, as well as different selection regimes. Finally, we show that a simple linear model generated from these simulations can predict the accuracy of haplotype-derived allele frequencies in other model organisms and experimental designs. To make these results broadly accessible for use in E+R experiments, we introduce HAF-pipe, an open-source software tool for calculating haplotype-derived allele frequencies from raw sequencing data. Ultimately, by reducing sequencing costs without sacrificing accuracy, our method facilitates E+R designs with higher replication and resolution, and thereby, increased power to detect adaptive alleles.

A clinal polymorphism in the insulin signaling transcription factor foxo contributes to life-history adaptation in Drosophila

A fundamental aim of adaptation genomics is to identify polymorphisms that underpin variation in fitness traits. In Drosophila melanogaster, latitudinal life-history clines exist on multiple continents and make an excellent system for dissecting the genetics of adaptation. We have previously identified numerous clinal single-nucleotide polymorphism in insulin/insulin-like growth factor signaling (IIS), a pathway known from mutant studies to affect life history. However, the effects of natural variants in this pathway remain poorly understood. Here we investigate how two clinal alternative alleles at foxo, a transcriptional effector of IIS, affect fitness components (viability, size, starvation resistance, fat content). We assessed this polymorphism from the North American cline by reconstituting outbred populations, fixed for either the low- or high-latitude allele, from inbred DGRP lines. Because diet and temperature modulate IIS, we phenotyped alleles across two temperatures (18°C, 25°C) and two diets differing in sugar source and content. Consistent with clinal expectations, the high-latitude allele conferred larger body size and reduced wing loading. Alleles also differed in starvation resistance and expression of insulin-like receptor, a transcriptional target of FOXO. Allelic reaction norms were mostly parallel, with few GxE interactions. Together, our results suggest that variation in IIS makes a major contribution to clinal life-history adaptation.

Balancing Selection Drives the Maintenance of Genetic Variation in Drosophila Antimicrobial Peptides

Genes involved in immune defense against pathogens provide some of the most well-known examples of both directional and balancing selection. Antimicrobial peptides (AMPs) are innate immune effector genes, playing a key role in pathogen clearance in many species, including Drosophila. Conflicting lines of evidence have suggested that AMPs may be under directional, balancing, or purifying selection. Here, we use both a linear model and control-gene-based approach to show that balancing selection is an important force shaping AMP diversity in Drosophila. In Drosophila melanogaster, this is most clearly observed in ancestral African populations. Furthermore, the signature of balancing selection is even more striking once background selection has been accounted for. Balancing selection also acts on AMPs in Drosophila mauritiana, an isolated island endemic separated from D. melanogaster by about 4 Myr of evolution. This suggests that balancing selection may be broadly acting to maintain adaptive diversity in Drosophila AMPs, as has been found in other taxa.

Rhodnius prolixus: Identification of missing components of the IMD immune signaling pathway and functional characterization of its role in eliminating bacteria

The innate immune system in insects is regulated by specific signalling pathways. Most immune related pathways were identified and characterized in holometabolous insects such as Drosophila melanogaster, and it was assumed they would be highly conserved in all insects. The hemimetabolous insect, Rhodnius prolixus, has served as a model to study basic insect physiology, but also is a major vector of the human parasite, Trypanosoma cruzi, that causes 10,000 deaths annually. The publication of the R. prolixus genome revealed that one of the main immune pathways, the Immune-deficiency pathway (IMD), was incomplete and probably non-functional, an observation shared with other hemimetabolous insects including the pea aphid (Acyrthosiphon pisum) and the bedbug (Cimex lectularius). It was proposed that the IMD pathway is inactive in R. prolixus as an adaptation to prevent eliminating beneficial symbiont gut bacteria. We used bioinformatic analyses based on reciprocal BLAST and HMM-profile searches to find orthologs for most of the "missing" elements of the IMD pathway and provide data that these are regulated in response to infection with Gram-negative bacteria. We used RNAi strategies to demonstrate the role of the IMD pathway in regulating the expression of specific antimicrobial peptides (AMPs) in the fat body of R. prolixus. The data indicate that the IMD pathway is present and active in R. prolixus, which opens up new avenues of research on R. prolixus-T. cruzi interactions.

Comparative analyses of sperm DNA methylomes among human, mouse and cattle provide insights into epigenomic evolution and complex traits

Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.

Rapid seasonal evolution in innate immunity of wild Drosophila melanogaster

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of Drosophila melanogaster to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens Enterococcus faecalis and Providencia rettgeri to measure post-infection survival and bacterial load of wild D. melanogaster populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in Thioester-containing protein 3 (Tep3) have different functional responses to infection and that epistatic interactions among seasonal Tep3 and Drosomycin-like 6 (Dro6) alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.