Survey of Global Genetic Diversity Within the Drosophila Immune System

Natural selection has driven the recurrent loss of an immunity gene that protects Drosophila against a major natural parasite

Polymorphisms in immunity genes can have large effects on susceptibility to infection. To understand the origins of this variation, we have investigated the genetic basis of resistance to the parasitoid wasp Leptopilina boulardi in Drosophila melanogaster. We found that increased expression of the gene lectin-24A after infection by parasitic wasps was associated with a faster cellular immune response and greatly increased rates of killing the parasite. lectin-24A encodes a protein that is strongly up-regulated in the fat body after infection and localizes to the surface of the parasite egg. In certain susceptible lines, a deletion upstream of the lectin-24A has largely abolished expression. Other mutations predicted to abolish the function of this gene have arisen recurrently in this gene, with multiple loss-of-expression alleles and premature stop codons segregating in natural populations. The frequency of these alleles varies greatly geographically, and in some southern African populations, natural selection has driven them near to fixation. We conclude that natural selection has favored the repeated loss of an important component of the immune system, suggesting that in some populations, a pleiotropic cost to lectin-24A expression outweighs the benefits of resistance.

Signatures of Adaptation, Constraints, and Potential Redundancy in the Canonical Immune Genes of a Key Pollinator

All organisms require an immune system to recognize, differentiate, and defend against pathogens. From an evolutionary perspective, immune systems evolve under strong selective pressures exerted by fast-evolving pathogens. However, the functional diversity of the immune system means that different immune components and their associated genes may evolve under varying forms of selection. Insect pollinators, which provide essential ecosystem services, are an important system in which to understand how selection has shaped immune gene evolution as their populations are experiencing declines with pathogens highlighted as a potential contributing factor. To improve our understanding of the genetic variation found in the immune genes of an essential pollinator, we performed whole-genome resequencing of wild-caught Bombus terrestris males. We first assessed nucleotide diversity and extended haplotype homozygosity for canonical immune genes finding the strongest signatures of positive selection acting on genes involved in pathogen recognition and antiviral defense, possibly driven by growing pathogen spread in wild populations. We also identified immune genes evolving under strong purifying selection, highlighting potential constraints on the bumblebee immune system. Lastly, we highlight the potential loss of function alleles present in the immune genes of wild-caught haploid males, suggesting that such genes are potentially less essential for development and survival and represent redundancy in the gene repertoire of the bumblebee immune system. Collectively, our analysis provides novel insights into the recent evolutionary history of the immune system of a key pollinator, highlighting targets of selection, constraints to adaptation, and potential redundancy.

Antimicrobial peptides do not directly contribute to aging in Drosophila, but improve lifespan by preventing dysbiosis

Antimicrobial peptides (AMPs) are innate immune effectors first studied for their role in host defence. Recent studies have implicated these peptides in the clearance of aberrant cells and in neurodegenerative syndromes. In Drosophila, many AMPs are produced downstream of Toll and Imd NF-κB pathways upon infection. Upon aging, AMPs are upregulated, drawing attention to these molecules as possible causes of age-associated inflammatory diseases. However, functional studies overexpressing or silencing these genes have been inconclusive. Using an isogenic set of AMP gene deletions, we investigated the net impact of AMPs on aging. Overall, we found no major effect of individual AMPs on lifespan, with the possible exception of Defensin. However, ΔAMP14 flies lacking seven AMP gene families displayed reduced lifespan. Increased bacterial load in the food of aged ΔAMP14 flies suggested that their lifespan reduction was due to microbiome dysbiosis, consistent with a previous study. Moreover, germ-free conditions extended the lifespan of ΔAMP14 flies. Overall, our results did not point to an overt role of individual AMPs in lifespan. Instead, we found that AMPs collectively impact lifespan by preventing dysbiosis during aging.

Nuclear access of DNlg3 c-terminal fragment and its function in regulating innate immune response genes

Neuroligins (NLGNs) are one of the autism susceptibility genes, however, the mechanism that how dysfunction of NLGNs leads to Autism remains unclear. More and more studies have shown that the transcriptome alteration may be one of the important factors to generate Autism. Therefore, we are very concerned about whether Neuroligins would affect transcriptional regulation, which may at last lead to Autism. As a single-transmembrane receptor, proteolytic cleavage is one of the most important posttranslational modifications of NLGN proteins. In this study, we demonstrated the existence of DNlg3 C-terminal fragment. Studies in the S2 cells and HEK293T cells showed the evidence for nuclear access of the DNlg3 C-terminal fragment. Then we identified the possible targets of DNlg3 C-terminal fragment after its nuclear access by RNA-seq. The bioinformatics analysis indicated the transcriptome alteration between dnlg3 null flies and wild type flies focused on genes for the innate immune responses. These results were consistent with the infection hypotheses for autism. Our study revealed the nuclear access ability of DNlg3 c-terminal fragment and its possible function in transcriptional regulation of the innate immune response genes. This work provides the new links between synaptic adhesion molecule NLGNs and immune activation, which may help us to get a deeper understanding on the relationship between NLGNs and Autism.

Evolving and assembling to pierce through: Evolutionary and structural aspects of antimicrobial peptides

The burgeoning menace of antimicrobial resistance across the globe has necessitated investigations into other chemotherapeutic strategies to combat infections. Antimicrobial peptides, or host defense peptides, are a set of promising therapeutic candidates in this regard. Most of them cause membrane permeabilization and are a key component of the innate immune response to pathogenic invasion. It has also been reported that peptide self-assembly is a driving factor governing the microbicidal activity of these peptide candidates. While efforts have been made to develop novel synthetic peptides against various microbes, many clinical trials of such peptides have failed due to toxicity and hemolytic activity to the host. A function-guided rational peptide engineering, based on evolutionary principles, physicochemical properties and activity determinants of AMP activity, is expected to help in targeting specific microbes. Furthermore, it is important to develop a unified understanding of the evolution of AMPs in order to fully appreciate their importance in host defense. This review seeks to explore the evolution of AMPs and the physicochemical determinants of AMP activity. The specific interactions driving AMP self-assembly have also been reviewed, emphasizing implications of this self-assembly on microbicidal and immunomodulatory activity.

Draft genome of the blister beetle, Epicauta chinensis

Existence of cantharidin (CTD) in blister beetles is a significant ecological adaptive mechanism that defends against predators and regulates courtship and mating behaviors. To better understand CTD biosynthetic information as well as its biology and pharmacology, we assembled a genome of 151.88 Mb for Epicauta chinensis using PacBio sequencing technology. Gene annotation yielded 249,238 repeats, 527 non-coding RNAs and 12,520 protein-coding genes. Compared to other 11 insects, expansions of gene families in E. chinensis for most core gene families likely associated with environmental adaptation, such as chemoreception, immunity, and detoxification. We further annotated P450s and immune-related genes, a total of 117 putative P450s comprising 7 CYP2, 67 CYP3, 36 CYP4, and 7 mitochondrial P450s and 281 immune-related genes were identified. Comparative analysis of the insect immune repertoires indicated presence of immune genes detected only from Coleopteran insects such as MD2-like. This suggested a lineage-specific gene evolution for Coleopteran insects. Based on the gene family evolution analysis, we identified two probable candidate genes including CYP4TT1 and phytanoyl-CoA dioxygenase for CTD biosynthesis. The high-quality reference genome of E. chinensis provides the genetic basis for further investigation of CTD biosynthesis and in-depth studies of the development and evolution of blister beetles.

Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila

Background

Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements.

Results

In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context.

Conclusions

We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02471-3.

Population genomics reveals variable patterns of immune gene evolution in monarch butterflies (Danaus plexippus)

Humoral and cellular immune responses provide animals with major defences against harmful pathogens. While it is often assumed that immune genes undergo rapid diversifying selection, this assumption has not been tested in many species. Moreover, it is likely that different classes of immune genes experience different levels of evolutionary constraint, resulting in varying selection patterns. We examined the evolutionary patterns for a set of 91 canonical immune genes of North American monarch butterflies (Danaus plexippus), using as an outgroup the closely related soldier butterfly (Danaus eresimus). As a comparison to these immune genes, we selected a set of control genes that were paired with each immune for approximate size and genomic location. As a whole, these immune genes had a significant but modest reduction in Tajima's D relative to paired-control genes, but otherwise did not show distinct patterns of population genetic variation or evolutionary rates. When further partitioning these immune genes into four functional classes (recognition, signalling, modulation, and effector), we found distinct differences among these groups. Relative to control genes, recognition genes exhibit increased nonsynonymous diversity and divergence, suggesting reduced constraints on evolution, and supporting the notion that coevolution with pathogens results in diversifying selection. In contrast, signalling genes showed an opposite pattern of reduced diversity and divergence, suggesting evolutionary constraints and conservation. Modulator and effector genes showed no statistical differences from controls. These results are consistent with patterns found in immune genes in fruit flies and Pieris butterflies, suggesting that consistent selective pressures on different classes of immune genes broadly govern the evolution of innate immunity among insects.

Dense time-course gene expression profiling of the Drosophila melanogaster innate immune response

BACKGROUND

Immune responses need to be initiated rapidly, and maintained as needed, to prevent establishment and growth of infections. At the same time, resources need to be balanced with other physiological processes. On the level of transcription, studies have shown that this balancing act is reflected in tight control of the initiation kinetics and shutdown dynamics of specific immune genes.

RESULTS

To investigate genome-wide expression dynamics and trade-offs after infection at a high temporal resolution, we performed an RNA-seq time course on D. melanogaster with 20 time points post Imd stimulation. A combination of methods, including spline fitting, cluster analysis, and Granger causality inference, allowed detailed dissection of expression profiles, lead-lag interactions, and functional annotation of genes through guilt-by-association. We identified Imd-responsive genes and co-expressed, less well characterized genes, with an immediate-early response and sustained up-regulation up to 5 days after stimulation. In contrast, stress response and Toll-responsive genes, among which were Bomanins, demonstrated early and transient responses. We further observed a strong trade-off with metabolic genes, which strikingly recovered to pre-infection levels before the immune response was fully resolved.

CONCLUSIONS

This high-dimensional dataset enabled the comprehensive study of immune response dynamics through the parallel application of multiple temporal data analysis methods. The well annotated data set should also serve as a useful resource for further investigation of the D. melanogaster innate immune response, and for the development of methods for analysis of a post-stress transcriptional response time-series at whole-genome scale.

A Population Genomic Investigation of Immune Cell Diversity and Phagocytic Capacity in a Butterfly

Insects rely on their innate immune system to successfully mediate complex interactions with their internal microbiota, as well as the microbes present in the environment. Given the variation in microbes across habitats, the challenges to respond to them are likely to result in local adaptations in the immune system. Here we focus upon phagocytosis, a mechanism by which pathogens and foreign particles are engulfed in order to be contained, killed, and processed. We investigated the phenotypic and genetic variation related to phagocytosis in two allopatric populations of the butterfly Pieris napi. Populations were found to differ in their hemocyte composition and overall phagocytic capability, driven by the increased phagocytic propensity of each cell type. Yet, genes annotated to phagocytosis showed no large genomic signal of divergence. However, a gene set enrichment analysis on significantly divergent genes identified loci involved in glutamine metabolism, which recently have been linked to immune cell differentiation in mammals. Together these results suggest that heritable variation in phagocytic capacity arises via a quantitative trait architecture with variation in genes affecting the activation and/or differentiation of phagocytic cells, suggesting them as potential candidate genes underlying these phenotypic differences.

Recurrent Collection of Drosophila melanogaster from Wild African Environments and Genomic Insights into Species History

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

Genomic landscape and genetic manipulation of the black soldier fly Hermetia illucens, a natural waste recycler

The black soldier fly (BSF), Hermetia illucens (Diptera: Stratiomyidae), is renowned for its bioconversion of organic waste into a sustainable source of animal feed. We report a high-quality genome of 1.1 Gb and a consensus set of 16,770 gene models for this beneficial species. Compared to those of other dipteran species, the BSF genome has undergone a substantial expansion in functional modules related to septic adaptation, including immune system factors, olfactory receptors, and cytochrome P450s. We further profiled midgut transcriptomes and associated microbiomes of BSF larvae fed with representative types of organic waste. We find that the pathways related to digestive system and fighting infection are commonly enriched and that Firmicutes bacteria dominate the microbial community in BSF across all diets. To extend its potential practical applications, we further developed an efficient CRISPR/Cas9-based gene editing approach and implemented this to yield flightless and enhanced feeding capacity phenotypes, both of which could expand BSF production capabilities. Our study provides valuable genomic and technical resources for optimizing BSF lines for industrialization.

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.

Mapping Selection within Drosophila melanogaster Embryo's Anatomy

We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.

Demographic History of the Human Commensal Drosophila melanogaster

The cohabitation of Drosophila melanogaster with humans is nearly ubiquitous. Though it has been well established that this fly species originated in sub-Saharan Africa, and only recently has spread globally, many details of its swift expansion remain unclear. Elucidating the demographic history of D. melanogaster provides a unique opportunity to investigate how human movement might have impacted patterns of genetic diversity in a commensal species, as well as providing neutral null models for studies aimed at identifying genomic signatures of local adaptation. Here, we use whole-genome data from five populations (Africa, North America, Europe, Central Asia, and the South Pacific) to carry out demographic inferences, with particular attention to the inclusion of migration and admixture. We demonstrate the importance of these parameters for model fitting and show that how previous estimates of divergence times are likely to be significantly underestimated as a result of not including them. Finally, we discuss how human movement along early shipping routes might have shaped the present-day population structure of D. melanogaster.

Draft genome of the cotton aphid Aphis gossypii

The cotton aphid Aphis gossypii Glover is a worldwide agricultural pest that feeds on cotton, melon, and other landscape plants, causing a high level of economic loss. In addition to the common characteristics shared with other aphids, Ap. gossypii has evolved multiple biotypes that present substantial differences in host adaption. These intriguing biological features are of interest from both a fundamental and applied perspective. However, the molecular studies of Ap. gossypii have been restrained by the lack of a reference genome. Furthermore, in order to establish a platform for the development of novel and sustainable control methods, it is necessary to generate genomic resources for Ap. gossypii. Here, we present a 294 Mb draft genome sequence of Ap. gossypii, which consists of 4,724 scaffolds with an N50 size of 438 kb. Compared to other aphid species with published genomes, Ap. gossypii presents the most compact genome size. A total of 14,694 protein-coding genes were predicted and annotated in the consensus gene set, 98.03% of CEGMA genes and 93.5% of BUSCO genes were captured respectively. Genome-wide selection analyses revealed that significantly evolving pathways in the genus Aphis are related to biological processes of detoxification, steroid biosynthesis, and ethylbenzene degradation. The acquisition of the genome of Ap. gossypii makes it possible to understand the molecular mechanism of intricate biological traits of this species, and will further facilitate the study of aphid evolution.

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.

Immune genes are hotspots of shared positive selection across birds and mammals

Consistent patterns of positive selection in functionally similar genes can suggest a common selective pressure across a group of species. We use alignments of orthologous protein-coding genes from 39 species of birds to estimate parameters related to positive selection for 11,000 genes conserved across birds. We show that functional pathways related to the immune system, recombination, lipid metabolism, and phototransduction are enriched for positively selected genes. By comparing our results with mammalian data, we find a significant enrichment for positively selected genes shared between taxa, and that these shared selected genes are enriched for viral immune pathways. Using pathogen-challenge transcriptome data, we show that genes up-regulated in response to pathogens are also enriched for positively selected genes. Together, our results suggest that pathogens, particularly viruses, consistently target the same genes across divergent clades, and that these genes are hotspots of host-pathogen conflict over deep evolutionary time.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

Dynamic miRNA-mRNA regulations are essential for maintaining Drosophila immune homeostasis during Micrococcus luteus infection

Pathogen bacteria infections can lead to dynamic changes of microRNA (miRNA) and mRNA expression profiles, which may control synergistically the outcome of immune responses. To reveal the role of dynamic miRNA-mRNA regulation in Drosophila innate immune responses, we have detailedly analyzed the paired miRNA and mRNA expression profiles at three time points during Drosophila adult males with Micrococcus luteus (M. luteus) infection using RNA- and small RNA-seq data. Our results demonstrate that differentially expressed miRNAs and mRNAs represent extensively dynamic changes over three time points during Drosophila with M. luteus infection. The pathway enrichment analysis indicates that differentially expressed genes are involved in diverse signaling pathways, including Toll and Imd as well as orther signaling pathways at three time points during Drosophila with M. luteus infection. Remarkably, the dynamic change of miRNA expression is delayed by compared to mRNA expression change over three time points, implying that the "time" parameter should be considered when the function of miRNA/mRNA is further studied. In particular, the dynamic miRNA-mRNA regulatory networks have shown that miRNAs may synergistically regulate gene expressions of different signaling pathways to promote or inhibit innate immune responses and maintain homeostasis in Drosophila, and some new regulators involved in Drosophila innate immune response have been identified. Our findings strongly suggest that miRNA regulation is a key mechanism involved in fine-tuning cooperatively gene expressions of diverse signaling pathways to maintain innate immune response and homeostasis in Drosophila. Taken together, the present study reveals a novel role of dynamic miRNA-mRNA regulation in immune response to bacteria infection, and provides a new insight into the underlying molecular regulatory mechanism of Drosophila innate immune responses.

Genomic signatures of local adaptation in the Drosophila immune response

As environments and pathogen landscapes shift, host defenses must evolve to remain effective. Due to this selection pressure, among-species comparisons of genetic sequence data often find immune genes to be among the fastest evolving genes across the genome. The full extent and nature of these immune adaptations, however, remain largely unexplored. In a recent study, we analyzed patterns of selection within distinct components of the Drosophila melanogaster immune pathway. While we found evidence of positive selection within some immune processes, immune genes were not universally characterized by signatures of strong selection. On the contrary, we even found that some immune functions show greater than expected constraint. Overall these results highlight 2 major factors that appear to play an outsize role in determining a gene's evolutionary rate: the type of pathogen the gene targets and the gene's position within the immune network. These results join a growing body of literature that highlight the complexity of immune adaptation. Rather than there being uniformly strong selection across all immune genes, a combination of pathogen-specificity and host genetic constraints appear to play key roles in determining each immune gene's individual evolutionary trajectory.

Costs and benefits of sublethal Drosophila C virus infection

Viruses are major evolutionary drivers of insect immune systems. Much of our knowledge of insect immune responses derives from experimental infections using the fruit fly Drosophila melanogaster. Most experiments, however, employ lethal pathogen doses through septic injury, frequently overwhelming host physiology. While this approach has revealed several immune mechanisms, it is less informative about the fitness costs hosts may experience during infection in the wild. Using both systemic and oral infection routes, we find that even apparently benign, sublethal infections with the horizontally transmitted Drosophila C virus (DCV) can cause significant physiological and behavioural morbidity that is relevant for host fitness. We describe DCV-induced effects on fly reproductive output, digestive health and locomotor activity, and we find that viral morbidity varies according to the concentration of pathogen inoculum, host genetic background and sex. Notably, sublethal DCV infection resulted in a significant increase in fly reproduction, but this effect depended on host genotype. We discuss the relevance of sublethal morbidity for Drosophila ecology and evolution, and more broadly, we remark on the implications of deleterious and beneficial infections for the evolution of insect immunity.