Genome-wide association studies reveal a simple genetic basis of resistance to naturally coevolving viruses in Drosophila melanogaster

Effects of spontaneous mutations on survival and reproduction of Drosophila serrata infected with Drosophila C virus

The impact of selection on host immune function genes has been widely documented. However, it remains essentially unknown how mutation influences the quantitative immune traits that selection acts on. Applying a classical mutation accumulation (MA) experimental design in Drosophila serrata, we found the mutational variation in susceptibility (median time of death, LT50) to Drosophila C virus (DCV) was of similar magnitude to that reported for intrinsic survival traits. Mean LT50 did not change as mutations accumulated, suggesting no directional bias in mutational effects. Maintenance of genetic variance in immune function is hypothesized to be influenced by pleiotropic effects on immunity and other traits that contribute to fitness. To investigate this, we assayed female reproductive output for a subset of MA lines with relatively long or short survival times under DCV infection. Longer survival time tended to be associated with lower reproductive output, suggesting that mutations affecting susceptibility to DCV had pleiotropic effects on investment in reproductive fitness. Further studies are needed to uncover the general patterns of mutational effect on immune responses and other fitness traits, and to determine how selection might typically act on new mutations via their direct and pleiotropic effects.

Evolutionary immunology to explore original antiviral strategies

Over the past 25 years, the field of evolutionary developmental biology (evo-devo) has used genomics and genetics to gain insight on the developmental mechanisms underlying the evolution of morphological diversity of animals. Evo-devo exploits the key insight that conserved toolkits of development (e.g. Hox genes) are used in animals to produce genetic novelties that provide adaptation to a new environment. Like development, immunity is forged by interactions with the environment, namely the microbial world. Yet, when it comes to the study of immune defence mechanisms in invertebrates, interest primarily focuses on evolutionarily conserved molecules also present in humans. Here, focusing on antiviral immunity, we argue that immune genes not conserved in humans represent an unexplored resource for the discovery of new antiviral strategies. We review recent findings on the cGAS-STING pathway and explain how cyclic dinucleotides produced by cGAS-like receptors may be used to investigate the portfolio of antiviral genes in a broad range of species. This will set the stage for evo-immuno approaches, exploiting the investment in antiviral defences made by metazoans over hundreds of millions of years of evolution. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Investigating the Evolution of Drosophila STING-Dependent Antiviral Innate Immunity by Multispecies Comparison of 2'3'-cGAMP Responses

Viruses represent a major threat to all animals, which defend themselves through induction of a large set of virus-stimulated genes that collectively control the infection. In vertebrates, these genes include interferons that play a critical role in the amplification of the response to infection. Virus- and interferon-stimulated genes include restriction factors targeting the different steps of the viral replication cycle, in addition to molecules associated with inflammation and adaptive immunity. Predictably, antiviral genes evolve dynamically in response to viral pressure. As a result, each animal has a unique arsenal of antiviral genes. Here, we exploit the capacity to experimentally activate the evolutionarily conserved stimulator of IFN genes (STING) signaling pathway by injection of the cyclic dinucleotide 2'3'-cyclic guanosine monophosphate-adenosine monophosphate into flies to define the repertoire of STING-regulated genes in 10 Drosophila species, spanning 40 million years of evolution. Our data reveal a set of conserved STING-regulated factors, including STING itself, a cGAS-like-receptor, the restriction factor pastel, and the antiviral protein Vago, but also 2 key components of the antiviral RNA interference pathway, Dicer-2, and Argonaute2. In addition, we identify unknown species- or lineage-specific genes that have not been previously associated with resistance to viruses. Our data provide insight into the core antiviral response in Drosophila flies and pave the way for the characterization of previously unknown antiviral effectors.

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

Parallel evolution and enhanced virulence upon in vivo passage of an RNA virus in Drosophila melanogaster

Virus evolution is strongly affected by antagonistic co-evolution of virus and host. Host immunity positively selects for viruses that evade the immune response, which in turn may drive counter-adaptations in host immune genes. We investigated how host immune pressure shapes virus populations, using the fruit fly Drosophila melanogaster and its natural pathogen Drosophila C virus (DCV), as a model. We performed an experimental evolution study in which DCV was serially passaged for ten generations in three fly genotypes differing in their antiviral RNAi response: wild-type flies and flies in which the endonuclease gene Dicer-2 was either overexpressed or inactivated. All evolved virus populations replicated more efficiently in vivo and were more virulent than the parental stock. The number of polymorphisms increased in all three host genotypes with passage number, which was most pronounced in Dicer-2 knockout flies. Mutational analysis showed strong parallel evolution, as mutations accumulated in a specific region of the VP3 capsid protein in every lineage in a host genotype-independent manner. The parental tyrosine at position ninety-five of VP3 was substituted with either one of five different amino acids in fourteen out of fifteen lineages. However, no consistent amino acid changes were observed in the viral RNAi suppressor gene 1A, nor elsewhere in the genome in any of the host backgrounds. Our study indicates that the RNAi response restricts the sequence space that can be explored by viral populations. Moreover, our study illustrates how evolution towards higher virulence can be a highly reproducible, yet unpredictable process.

The genetic basis of variation in immune defense against Lysinibacillus fusiformis infection in Drosophila melanogaster

The genetic causes of phenotypic variation often differ depending on the population examined, particularly if the populations were founded by relatively small numbers of genotypes. Similarly, the genetic causes of phenotypic variation among similar traits (resistance to different xenobiotic compounds or pathogens) may also be completely different or only partially overlapping. Differences in genetic causes for variation in the same trait among populations suggests context dependence for how selection acts on those traits. Similarities in the genetic causes of variation for different traits, on the other hand, suggests pleiotropy which would also influence how natural selection shapes variation in a trait. We characterized immune defense against a natural Drosophila pathogen, the Gram-positive bacterium Lysinibacillus fusiformis, in three different populations and found almost no overlap in the genetic architecture of variation in survival post infection. However, when comparing our results to a similar experiment with the fungal pathogen, B. bassiana, we found a convincing shared QTL peak for both pathogens. This peak contains the Bomanin cluster of Drosophila immune effectors. Loss of function mutants and RNAi knockdown experiments confirms a role of some of these genes in immune defense against both pathogens. This suggests that natural selection may act on the entire cluster of Bomanin genes (and the linked region under the QTL) or specific peptides for specific pathogens.

Investigating the outcomes of virus coinfection within and across host species

Interactions between coinfecting pathogens have the potential to alter the course of infection and can act as a source of phenotypic variation in susceptibility between hosts. This phenotypic variation may influence the evolution of host-pathogen interactions within host species and interfere with patterns in the outcomes of infection across host species. Here, we examine experimental coinfections of two Cripaviruses-Cricket Paralysis Virus (CrPV), and Drosophila C Virus (DCV)-across a panel of 25 Drosophila melanogaster inbred lines and 47 Drosophilidae host species. We find that interactions between these viruses alter viral loads across D. melanogaster genotypes, with a ~3 fold increase in the viral load of DCV and a ~2.5 fold decrease in CrPV in coinfection compared to single infection, but we find little evidence of a host genetic basis for these effects. Across host species, we find no evidence of systematic changes in susceptibility during coinfection, with no interaction between DCV and CrPV detected in the majority of host species. These results suggest that phenotypic variation in coinfection interactions within host species can occur independently of natural host genetic variation in susceptibility, and that patterns of susceptibility across host species to single infections can be robust to the added complexity of coinfection.

Dengue and chikungunya virus loads in the mosquito Aedes aegypti are determined by distinct genetic architectures

Aedes aegypti is the primary vector of the arboviruses dengue (DENV) and chikungunya (CHIKV). These viruses exhibit key differences in their vector interactions, the latter moving more quicky through the mosquito and triggering fewer standard antiviral pathways. As the global footprint of CHIKV continues to expand, we seek to better understand the mosquito's natural response to CHIKV-both to compare it to DENV:vector coevolutionary history and to identify potential targets in the mosquito for genetic modification. We used a modified full-sibling design to estimate the contribution of mosquito genetic variation to viral loads of both DENV and CHIKV. Heritabilities were significant, but higher for DENV (40%) than CHIKV (18%). Interestingly, there was no genetic correlation between DENV and CHIKV loads between siblings. These data suggest Ae. aegypti mosquitoes respond to the two viruses using distinct genetic mechanisms. We also examined genome-wide patterns of gene expression between High and Low CHIKV families representing the phenotypic extremes of viral load. Using RNAseq, we identified only two loci that consistently differentiated High and Low families: a long non-coding RNA that has been identified in mosquito screens post-infection and a distant member of a family of Salivary Gland Specific (SGS) genes. Interestingly, the latter gene is also associated with horizontal gene transfer between mosquitoes and the endosymbiotic bacterium Wolbachia. This work is the first to link the SGS gene to a mosquito phenotype. Understanding the molecular details of how this gene contributes to viral control in mosquitoes may, therefore, also shed light on its role in Wolbachia.

Effects of Wolbachia on Transposable Element Expression Vary Between Drosophila melanogaster Host Genotypes

Transposable elements (TEs) are repetitive DNA sequences capable of changing position in host genomes, thereby causing mutations. TE insertions typically have deleterious effects but they can also be beneficial. Increasing evidence of the contribution of TEs to adaptive evolution further raises interest in understanding what factors impact TE activity. Based on previous studies associating the bacterial endosymbiont Wolbachia with changes in the abundance of piRNAs, a mechanism for TE repression, and to transposition of specific TEs, we hypothesized that Wolbachia infection would interfere with TE activity. We tested this hypothesis by studying the expression of 14 TEs in a panel of 25 Drosophila melanogaster host genotypes, naturally infected with Wolbachia and annotated for TE insertions. The host genotypes differed significantly in Wolbachia titers inside individual flies, with broad-sense heritability around 20%, and in the number of TE insertions, which depended greatly on TE identity. By removing Wolbachia from the target host genotypes, we generated a panel of 25 pairs of Wolbachia-positive and Wolbachia-negative lines in which we quantified transcription levels for our target TEs. We found variation in TE expression that was dependent on Wolbachia status, TE identity, and host genotype. Comparing between pairs of Wolbachia-positive and Wolbachia-negative flies, we found that Wolbachia removal affected TE expression in 21.1% of the TE-genotype combinations tested, with up to 2.3 times differences in the median level of transcript. Our data show that Wolbachia can impact TE activity in host genomes, underscoring the importance this endosymbiont can have in the generation of genetic novelty in hosts.

Intraspecific genetic variation in host vigour, viral load and disease tolerance during Drosophila C virus infection

Genetic variation for resistance and disease tolerance has been described in a range of species. In Drosophila melanogaster, genetic variation in mortality following systemic Drosophila C virus (DCV) infection is driven by large-effect polymorphisms in the restriction factor pastrel (pst). However, it is unclear if pst contributes to disease tolerance. We investigated systemic DCV challenges spanning nine orders of magnitude, in males and females of 10 Drosophila Genetic Reference Panel lines carrying either a susceptible (S) or resistant (R) pst allele. We find among-line variation in fly survival, viral load and disease tolerance measured both as the ability to maintain survival (mortality tolerance) and reproduction (fecundity tolerance). We further uncover novel effects of pst on host vigour, as flies carrying the R allele exhibited higher survival and fecundity even in the absence of infection. Finally, we found significant genetic variation in the expression of the JAK-STAT ligand upd3 and the epigenetic regulator of JAK-STAT G9a. However, while G9a has been previously shown to mediate tolerance of DCV infection, we found no correlation between the expression of either upd3 or G9a on fly tolerance or resistance. Our work highlights the importance of both resistance and tolerance in viral defence.

Drosophila Ectoderm-expressed 4 modulates JAK/STAT pathway and protects flies against Drosophila C virus infection

Sterile alpha and HEAT/Armadillo motif-containing protein (SARM) is conserved in evolution and negatively regulates TRIF-dependent Toll signaling in mammals. The SARM protein from Litopenaeus vannamei and its Drosophila orthologue Ectoderm-expressed (Ect4) are also involved in immune defense against pathogen infection. However, the functional mechanism of the protective effect remains unclear. In this study, we show that Ect4 is essential for the viral load in flies after a Drosophila C virus (DCV) infection. Viral load is increased in Ect4 mutants resulting in higher mortality rates than wild-type. Overexpression of Ect4 leads to a suppression of virus replication and thus improves the survival rate of the animals. Ect4 is required for the viral induction of STAT-responsive genes, TotA and TotM. Furthermore, Ect4 interacts with Stat92E, affecting the tyrosine phosphorylation and nuclear translocation of Stat92E in S2 cells. Altogether, our study identifies the adaptor protein Ect4 of the Toll pathway contributes to resistance to viral infection and regulates JAK/STAT signaling pathway.

Discovery of malathion resistance QTL in Drosophila melanogaster using a bulked phenotyping approach

Drosophila melanogaster has proved an effective system with which to understand the evolutionary genetics and molecular mechanisms of insecticide resistance. Insecticide use has left signatures of selection in the fly genome, and both functional and quantitative genetic studies in the system have identified genes and variants associated with resistance. Here, we use D. melanogaster and leverage a bulk phenotyping and pooled sequencing "extreme quantitative trait loci" approach to genetically dissect variation in resistance to malathion, an organophosphate insecticide. We resolve 2 quantitative trait loci, one of which implicates allelic variation at the cytochrome P450 gene Cyp6g1, a strong candidate based on previous work. The second shows no overlap with hits from a previous genome-wide association study for malathion resistance, recapitulating other studies showing that different strategies for complex trait dissection in flies can yield apparently different architectures. Notably, we see no genetic signal at the Ace gene. Ace encodes the target of organophosphate insecticide inhibition, and genome-wide association studies have identified strong Ace-linked associations with resistance in flies. The absence of quantitative trait locus implicating Ace here is most likely because our mapping population does not segregate for several of the known functional polymorphisms impacting resistance at Ace, perhaps because our population is derived from flies collected prior to the widespread use of organophosphate insecticides. Our fundamental approach can be an efficient, powerful strategy to dissect genetic variation in resistance traits. Nonetheless, studies seeking to interrogate contemporary insecticide resistance variation may benefit from deriving mapping populations from more recently collected strains.

A novel transposable element-mediated mechanism causes antiviral resistance in Drosophila through truncating the Veneno protein

Hosts are continually selected to evolve new defenses against an ever-changing array of pathogens. To understand this process, we examined the genetic basis of resistance to the Drosophila A virus in Drosophila melanogaster. In a natural population, we identified a polymorphic transposable element (TE) insertion that was associated with an ∼19,000-fold reduction in viral titers, allowing flies to largely escape the harmful effects of infection by this virulent pathogen. The insertion occurs in the protein-coding sequence of the gene Veneno, which encodes a Tudor domain protein. By mutating Veneno with CRISPR-Cas9 in flies and expressing it in cultured cells, we show that the ancestral allele of the gene has no effect on viral replication. Instead, the TE insertion is a gain-of-function mutation that creates a gene encoding a novel resistance factor. Viral titers remained reduced when we deleted the TE sequence from the transcript, indicating that resistance results from the TE truncating the Veneno protein. This is a novel mechanism of virus resistance and a new way by which TEs can contribute to adaptation.

Powerful, efficient QTL mapping in Drosophila melanogaster using bulked phenotyping and pooled sequencing

Despite the value of recombinant inbred lines for the dissection of complex traits, large panels can be difficult to maintain, distribute, and phenotype. An attractive alternative to recombinant inbred lines for many traits leverages selecting phenotypically extreme individuals from a segregating population, and subjecting pools of selected and control individuals to sequencing. Under a bulked or extreme segregant analysis paradigm, genomic regions contributing to trait variation are revealed as frequency differences between pools. Here, we describe such an extreme quantitative trait locus, or extreme quantitative trait loci, mapping strategy that builds on an existing multiparental population, the Drosophila Synthetic Population Resource, and involves phenotyping and genotyping a population derived by mixing hundreds of Drosophila Synthetic Population Resource recombinant inbred lines. Simulations demonstrate that challenging, yet experimentally tractable extreme quantitative trait loci designs (≥4 replicates, ≥5,000 individuals/replicate, and selecting the 5-10% most extreme animals) yield at least the same power as traditional recombinant inbred line-based quantitative trait loci mapping and can localize variants with sub-centimorgan resolution. We empirically demonstrate the effectiveness of the approach using a 4-fold replicated extreme quantitative trait loci experiment that identifies 7 quantitative trait loci for caffeine resistance. Two mapped extreme quantitative trait loci factors replicate loci previously identified in recombinant inbred lines, 6/7 are associated with excellent candidate genes, and RNAi knock-downs support the involvement of 4 genes in the genetic control of trait variation. For many traits of interest to drosophilists, a bulked phenotyping/genotyping extreme quantitative trait loci design has considerable advantages.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells

Background

Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown.

Results

The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D.

Conclusion

Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08057-4.

Host genetic factors associated with the range limit of a European hantavirus

The natural host ranges of many viruses are restricted to very specific taxa. Little is known about the molecular barriers between species that lead to the establishment of this restriction or generally prevent virus emergence in new hosts. Here, we identify genomic polymorphisms in a natural rodent host associated with a strong genetic barrier to the transmission of European Tula orthohantavirus (TULV). We analysed the very abrupt spatial transition between two major phylogenetic clades in TULV across the comparatively much wider natural hybrid zone between evolutionary lineages of their reservoir host, the common vole (Microtus arvalis). Genomic scans of 79,225 single nucleotide polymorphisms (SNPs) in 323 TULV‐infected host individuals detected 30 SNPs that were consistently associated with the TULV clades CEN.S or EST.S in two replicate sampling transects. Focusing the analysis on 199 voles with evidence of genomic admixture at the individual level (0.1–0.9) supported statistical significance for all 30 loci. Host genomic variation at these SNPs explained up to 37.6% of clade‐specific TULV infections. Genes in the vicinity of associated SNPs include SAHH, ITCH and two members of the Syngr gene family, which are involved in functions related to immune response or membrane transport. This study demonstrates the relevance of natural hybrid zones as systems not only for studying processes of evolutionary divergence and speciation, but also for the detection of evolving genetic barriers for specialized parasites.

Between virus correlations in the outcome of infection across host species: Evidence of virus by host species interactions

Virus host shifts are a major source of outbreaks and emerging infectious diseases, and predicting the outcome of novel host and virus interactions remains a key challenge for virus research. The evolutionary relationships between host species can explain variation in transmission rates, virulence, and virus community composition between hosts, but it is unclear if correlations exist between related viruses in infection traits across novel hosts. Here, we measure correlations in viral load of four Cripavirus isolates across experimental infections of 45 Drosophilidae host species. We find positive correlations between every pair of viruses tested, suggesting that some host clades show broad susceptibility and could act as reservoirs and donors for certain types of viruses. Additionally, we find evidence of virus by host species interactions, highlighting the importance of both host and virus traits in determining the outcome of virus host shifts. Of the four viruses tested here, those that were more closely related tended to be more strongly correlated, providing tentative evidence that virus evolutionary relatedness may be a useful proxy for determining the likelihood of novel virus emergence, which warrants further research.

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10⁻⁵) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-β signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

A Two-Locus System with Strong Epistasis Underlies Rapid Parasite-Mediated Evolution of Host Resistance

Parasites are a major evolutionary force, driving adaptive responses in host populations. Although the link between phenotypic response to parasite-mediated natural selection and the underlying genetic architecture often remains obscure, this link is crucial for understanding the evolution of resistance and predicting associated allele frequency changes in the population. To close this gap, we monitored the response to selection during epidemics of a virulent bacterial pathogen, Pasteuria ramosa, in a natural host population of Daphnia magna. Across two epidemics, we observed a strong increase in the proportion of resistant phenotypes as the epidemics progressed. Field and laboratory experiments confirmed that this increase in resistance was caused by selection from the local parasite. Using a genome-wide association study, we built a genetic model in which two genomic regions with dominance and epistasis control resistance polymorphism in the host. We verified this model by selfing host genotypes with different resistance phenotypes and scoring their F1 for segregation of resistance and associated genetic markers. Such epistatic effects with strong fitness consequences in host–parasite coevolution are believed to be crucial in the Red Queen model for the evolution of genetic recombination.

The Antiviral Effects of the Symbiont Bacteria Wolbachia in Insects

Wolbachia is a maternally transmitted bacterium that lives inside arthropod cells. Historically, it was viewed primarily as a parasite that manipulates host reproduction, but more recently it was discovered that Wolbachia can also protect Drosophila species against infection by RNA viruses. Combined with Wolbachia's ability to invade insect populations due to reproductive manipulations, this provides a way to modify mosquito populations to prevent them transmitting viruses like dengue. In this review, we discuss the main advances in the field since Wolbachia's antiviral effect was discovered 12 years ago, identifying current research gaps and potential future developments. We discuss that the antiviral effect works against a broad range of RNA viruses and depends on the Wolbachia lineage. We describe what is known about the mechanisms behind viral protection, and that recent studies suggest two possible mechanisms: activation of host immunity or competition with virus for cellular resources. We also discuss how association with Wolbachia may influence the evolution of virus defense on the insect host genome. Finally, we investigate whether the antiviral effect occurs in wild insect populations and its ecological relevance as a major antiviral component in insects.

Dissecting genetic and sex-specific sources of host heterogeneity in pathogen shedding and spread

Host heterogeneity in disease transmission is widespread but precisely how different host traits drive this heterogeneity remains poorly understood. Part of the difficulty in linking individual variation to population-scale outcomes is that individual hosts can differ on multiple behavioral, physiological and immunological axes, which will together impact their transmission potential. Moreover, we lack well-characterized, empirical systems that enable the quantification of individual variation in key host traits, while also characterizing genetic or sex-based sources of such variation. Here we used Drosophila melanogaster and Drosophila C Virus as a host-pathogen model system to dissect the genetic and sex-specific sources of variation in multiple host traits that are central to pathogen transmission. Our findings show complex interactions between genetic background, sex, and female mating status accounting for a substantial proportion of variance in lifespan following infection, viral load, virus shedding, and viral load at death. Two notable findings include the interaction between genetic background and sex accounting for nearly 20% of the variance in viral load, and genetic background alone accounting for ~10% of the variance in viral shedding and in lifespan following infection. To understand how variation in these traits could generate heterogeneity in individual pathogen transmission potential, we combined measures of lifespan following infection, virus shedding, and previously published data on fly social aggregation. We found that the interaction between genetic background and sex explained ~12% of the variance in individual transmission potential. Our results highlight the importance of characterising the sources of variation in multiple host traits to understand the drivers of heterogeneity in disease transmission.

2'3'-cGAMP triggers a STING- and NF-κB-dependent broad antiviral response in Drosophila

We previously reported that an ortholog of STING regulates infection by picorna-like viruses in Drosophila In mammals, STING is activated by the cyclic dinucleotide 2'3'-cGAMP produced by cGAS, which acts as a receptor for cytosolic DNA. Here, we showed that injection of flies with 2'3'-cGAMP induced the expression of dSTING-regulated genes. Coinjection of 2'3'-cGAMP with a panel of RNA or DNA viruses resulted in substantially reduced viral replication. This 2'3'-cGAMP-mediated protection was still observed in flies with mutations in Atg7 and AGO2, genes that encode key components of the autophagy and small interfering RNA pathways, respectively. By contrast, this protection was abrogated in flies with mutations in the gene encoding the NF-κB transcription factor Relish. Transcriptomic analysis of 2'3'-cGAMP-injected flies revealed a complex response pattern in which genes were rapidly induced, induced after a delay, or induced in a sustained manner. Our results reveal that dSTING regulates an NF-κB-dependent antiviral program that predates the emergence of interferons in vertebrates.

Host-shift as the cause of emerging infectious diseases: Experimental approaches using Drosophila-virus interactions

Host shifts, when a cross-species transmission of a pathogen can lead to successful infections, are the main cause of emerging infectious diseases, such as COVID-19. A complex challenge faced by the scientific community is to address the factors that determine whether the cross-species transmissions will result in spillover or sustained onwards infections. Here we review recent literature and present a perspective on current approaches we are using to understand the mechanisms underlying host shifts. We highlight the usefulness of the interactions between Drosophila species and viruses as an ideal study model. Additionally, we discuss how cross-infection experiments - when pathogens from a natural reservoir are intentionally injected in novel host species- can test the effect cross-species transmissions may have on the fitness of virus and host, and how the host phylogeny may influence this response. We also discuss experiments evaluating how cooccurrence with other viruses or the presence of the endosymbiont bacteria Wolbachia may affect the performance of new viruses in a novel host. Finally, we discuss the need of surveys of virus diversity in natural populations using next-generation sequencing technologies. In the long term, these approaches can contribute to a better understanding of the basic biology of host shifts.

Partitiviruses Infecting Drosophila melanogaster and Aedes aegypti Exhibit Efficient Biparental Vertical Transmission

Partitiviruses are segmented, multipartite double-stranded RNA (dsRNA) viruses that until recently were only known to infect fungi, plants, and protozoans. Metagenomic surveys have revealed that partitivirus-like sequences are also commonly associated with arthropods. One arthropod-associated partitivirus, galbut virus, is common in wild populations of Drosophila melanogaster To begin to understand the processes that underlie this virus's high global prevalence, we established colonies of wild-caught infected flies. Infection remained at stably high levels over 3 years, with between 63 and 100% of individual flies infected. Galbut virus infects fly cells and replicates in tissues throughout infected adults, including reproductive tissues and the gut epithelium. We detected no evidence of horizontal transmission via ingestion, but vertical transmission from either infected females or infected males was ∼100% efficient. Vertical transmission of a related partitivirus, verdadero virus, that we discovered in a laboratory colony of Aedes aegypti mosquitoes was similarly efficient. This suggests that efficient biparental vertical transmission may be a feature of at least a subset of insect-infecting partitiviruses. To study the impact of galbut virus infection free from the confounding effect of other viruses, we generated an inbred line of flies with galbut virus as the only detectable virus infection. We were able to transmit infection experimentally via microinjection of homogenate from these galbut-only flies. This sets the stage for experiments to understand the biological impact and possible utility of partitiviruses infecting model organisms and disease vectors.IMPORTANCE Galbut virus is a recently discovered partitivirus that is extraordinarily common in wild populations of the model organism Drosophila melanogaster Like for most viruses discovered through metagenomics, most of the basic biological questions about this virus remain unanswered. We found that galbut virus, along with a closely related partitivirus found in Aedes aegypti mosquitoes, is transmitted from infected females or males to offspring with ∼100% efficiency and can be maintained in laboratory colonies over years. This efficient transmission mechanism likely underlies the successful spread of these viruses through insect populations. We created Drosophila lines that contained galbut virus as the only virus infection and showed that these flies can be used as a source for experimental infections. This provides insight into how arthropod-infecting partitiviruses may be maintained in nature and sets the stage for exploration of their biology and potential utility.

Genetic variation for resistance to the specific fly pathogen Entomophthora muscae

We found substantial variation in resistance to the fly-specific pathogen Entomophthora muscae 'Berkeley' (Entomophthoromycota), in 20 lines from the Drosophila melanogaster Genetic Reference Panel (DGRP). Resistance to E. muscae is positively (r = 0.55) correlated with resistance to the broad host range ascomycete entomopathogen Metarhizium anisopliae (Ma549), indicative of generalist (non-specific) defenses. Most of the lines showing above average resistance to Ma549 showed cross-resistance to E. muscae. However, lines that succumbed quickly to Ma549 exhibited the full range of resistance to E. muscae. This suggests fly populations differ in E. muscae-specific resistance mechanisms as well as generic defences effective against both Ma549 and E. muscae. We looked for trade-offs that could account for inter-line variation, but increases (decreases) in disease resistance to E. muscae are not consistently associated with increases (decreases) of resistance to oxidative stress, starvation stress and sleep indices. That these pathogens are dynamic agents of selection on hosts is reflected in this genetic variation for resistance in lines derived from wild populations.

Resource quality determines the evolution of resistance and its genetic basis

Parasites impose strong selection on their hosts, but the level of any evolved resistance may be constrained by the availability of resources. However, studies identifying the genomic basis of such resource-mediated selection are rare, particularly in nonmodel organisms. Here, we investigated the role of nutrition in the evolution of resistance to a DNA virus (PiGV), and any associated trade-offs in a lepidopteran pest species (Plodia interpunctella). Through selection experiments and whole-genome resequencing, we identify genetic markers of resistance that vary between the nutritional environments during selection. We do not find consistent evolution of resistance in the presence of virus but rather see substantial variation among replicate populations. Resistance in a low-nutrition environment is negatively correlated with growth rate, consistent with an established trade-off between immunity and development, but this relationship is highly context dependent. Whole-genome resequencing of the host shows that resistance mechanisms are likely to be highly polygenic and although the underlying genetic architecture may differ between high and low-nutrition environments, similar mechanisms are commonly used. As a whole, our results emphasize the importance of the resource environment on influencing the evolution of resistance.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Temporal specificity and heterogeneity of Drosophila immune cells

Immune cells provide defense against non-self and have recently been shown to also play key roles in diverse processes such as development, metabolism, and tumor progression. The heterogeneity of Drosophila immune cells (hemocytes) remains an open question. Using bulk RNA sequencing, we find that the hemocytes display distinct features in the embryo, a closed and rapidly developing system, compared to the larva, which is exposed to environmental and metabolic challenges. Through single-cell RNA sequencing, we identify fourteen hemocyte clusters present in unchallenged larvae and associated with distinct processes, e.g., proliferation, phagocytosis, metabolic homeostasis, and humoral response. Finally, we characterize the changes occurring in the hemocyte clusters upon wasp infestation, which triggers the differentiation of a novel hemocyte type, the lamellocyte. This first molecular atlas of hemocytes provides insights and paves the way to study the biology of the Drosophila immune cells in physiological and pathological conditions.

The Genetic Basis of Natural Variation in Drosophila melanogaster Immune Defense against Enterococcus faecalis

Dissecting the genetic basis of natural variation in disease response in hosts provides insights into the coevolutionary dynamics of host-pathogen interactions. Here, a genome-wide association study of Drosophila melanogaster survival after infection with the Gram-positive entomopathogenic bacterium Enterococcus faecalis is reported. There was considerable variation in defense against E. faecalis infection among inbred lines of the Drosophila Genetics Reference Panel. We identified single nucleotide polymorphisms associated with six genes with a significant (p < 10^-08, corresponding to a false discovery rate of 2.4%) association with survival, none of which were canonical immune genes. To validate the role of these genes in immune defense, their expression was knocked-down using RNAi and survival of infected hosts was followed, which confirmed a role for the genes krishah and S6k in immune defense. We further identified a putative role for the Bomanin gene BomBc1 (also known as IM23), in E. faecalis infection response. This study adds to the growing set of association studies for infection in Drosophila melanogaster and suggests that the genetic causes of variation in immune defense differ for different pathogens.

In Vivo Characterization of a Bank Vole-Derived Cowpox Virus Isolate in Natural Hosts and the Rat Model

Cowpox virus (CPXV) belongs to the genus Orthopoxvirus in the Poxviridae family and is endemic in western Eurasia. Based on seroprevalence studies in different voles from continental Europe and UK, voles are suspected to be the major reservoir host. Recently, a CPXV was isolated from a bank vole (Myodes glareolus) in Germany that showed a high genetic similarity to another isolate originating from a Cotton-top tamarin (Saguinus oedipus). Here we characterize this first bank vole-derived CPXV isolate in comparison to the related tamarin-derived isolate. Both isolates grouped genetically within the provisionally called CPXV-like 3 clade. Previous phylogenetic analysis indicated that CPXV is polyphyletic and CPXV-like 3 clade represents probably a different species if categorized by the rules used for other orthopoxviruses. Experimental infection studies with bank voles, common voles (Microtusarvalis) and Wistar rats showed very clear differences. The bank vole isolate was avirulent in both common voles and Wistar rats with seroconversion seen only in the rats. In contrast, inoculated bank voles exhibited viral shedding and seroconversion for both tested CPXV isolates. In addition, bank voles infected with the tamarin-derived isolate experienced a marked weight loss. Our findings allow for the conclusion that CPXV isolates might differ in their replication capacity in different vole species and rats depending on their original host. Moreover, the results indicate host-specific differences concerning CPXV-specific virulence. Further experiments are needed to identify individual virulence and host factors involved in the susceptibility and outcome of CPXV-infections in the different reservoir hosts.

Drosophila immunity against natural and nonnatural viral pathogens

The fruit fly Drosophila melanogaster is extensively used as a model species for molecular biology and genetics. It is also widely studied for its innate immune system to expand our understanding of immune host defenses against numerous pathogens. More precisely, studies using both natural and nonnatural Drosophila pathogens have provided a better perspective of viral infection strategies and immunity processes than any other invertebrate. This has made significant advances in identifying and characterizing the innate immune mechanisms by which hosts can combat viral pathogens. However, in-depth studies on antiviral immunity are still lacking due in part to the narrow research focus on the evolution and conservation of antiviral strategies to combat infections caused by both natural and nonnatural viruses. In this review, we will cover three major areas. First, we will describe the well-characterized antiviral immune mechanisms in Drosophila. Second, we will survey the specific pathways induced by natural viruses that have been studied in Drosophila. Finally, we will discuss the pathways activated by nonnatural viruses, drawing comparisons to natural viruses and giving an unprecedented insight into the virus community of Drosophila that is necessary to understand the evolutionary and immune context needed to develop Drosophila as a model for virus research.

Use of DNA Markers for Grape Phylloxera Population and Evolutionary Genetics: From RAPDs to SSRs and Beyond

Grape phylloxera (Daktulosphaira vitifoliae Fitch) is a major pest of cultivated grapevines (Vitis spp.), occurring in virtually all viticultural regions around the world. Different grape phylloxera strains can be found at varying levels on leaves and roots on both own-rooted plants and in plants grafted onto partially resistant rootstocks. Considering its relevance for the adequate management of the pest in infested vineyards, the analysis of its genetic diversity has received considerable attention from the scientific community in the last decades. Here, we review 25 years of DNA-based molecular markers applied to the analysis of the genetic structure and the reproductive mode of grape phylloxera in its native range and in different introduced regions. The use given to RAPD, AFLP, mtDNA sequencing and microsatellite (SSR) genetic markers for the analysis of grape phylloxera diversity is discussed, and an overview of the main findings obtained after their application to different populations collected in diverse regions all around the world is shown. Lastly, we explore how recent advancements in molecular biology and in modern high throughput genotyping technologies may be applied to better understand grape phylloxera natural diversity at a molecular level.

An assessment of the immune costs associated with meiotic drive elements in Drosophila

Most organisms are constantly adapting to pathogens and parasites that exploit their host for their own benefit. Less studied, but perhaps more ubiquitous, are intragenomic parasites or selfish genetic elements. These include transposable elements, selfish B chromosomes and meiotic drivers that promote their own replication without regard to fitness effects on hosts. Therefore, intragenomic parasites are also a constant evolutionary pressure on hosts. Gamete-killing meiotic drive elements are often associated with large chromosomal inversions that reduce recombination between the drive and wild-type chromosomes. This reduced recombination is thought to reduce the efficacy of selection on the drive chromosome and allow for the accumulation of deleterious mutations. We tested whether gamete-killing meiotic drive chromosomes were associated with reduced immune defence against two bacterial pathogens in three species of Drosophila. We found little evidence of reduced immune defence in lines with meiotic drive. One line carrying the Drosophila melanogaster autosomal Segregation Distorter did show reduced defence, but we were unable to attribute that reduced defence to either genotype or immune gene expression differences. Our results suggest that though gamete-killing meiotic drive chromosomes probably accumulate deleterious mutations, those mutations do not result in reduced capacity for immune defence.

Viral infection causes sex-specific changes in fruit fly social aggregation behaviour

Host behavioural changes following infection are common and could be important determinants of host behavioural competence to transmit pathogens. Identifying potential sources of variation in sickness behaviours is therefore central to our understanding of disease transmission. Here, we test how group social aggregation and individual locomotor activity vary between different genotypes of male and female fruit flies (Drosophila melanogaster) following septic infection with Drosophila C virus (DCV). We find genetic-based variation in both locomotor activity and social aggregation, but we did not detect an effect of DCV infection on fly activity or sleep patterns within the initial days following infection. However, DCV infection caused sex-specific effects on social aggregation, as male flies in most genetic backgrounds increased the distance to their nearest neighbour when infected. We discuss possible causes for these differences in the context of individual variation in immunity and their potential consequences for disease transmission.

Dissecting the genetic architecture of a stepwise infection process

How a host fights infection depends on an ordered sequence of steps, beginning with attempts to prevent a pathogen from establishing an infection, through to steps that mitigate a pathogen's control of host resources or minimize the damage caused during infection. Yet empirically characterizing the genetic basis of these steps remains challenging. Although each step is likely to have a unique genetic and environmental signature, and may therefore respond to selection in different ways, events that occur earlier in the infection process can mask or overwhelm the contributions of subsequent steps. In this study, we dissect the genetic architecture of a stepwise infection process using a quantitative trait locus (QTL) mapping approach. We control for variation at the first line of defence against a bacterial pathogen and expose downstream genetic variability related to the host's ability to mitigate the damage pathogens cause. In our model, the water-flea Daphnia magna, we found a single major effect QTL, explaining 64% of the variance, that is linked to the host's ability to completely block pathogen entry by preventing their attachment to the host oesophagus; this is consistent with the detection of this locus in previous studies. In susceptible hosts allowing attachment, however, a further 23 QTLs, explaining between 5% and 16% of the variance, were mapped to traits related to the expression of disease. The general lack of pleiotropy and epistasis for traits related to the different stages of the infection process, together with the wide distribution of QTLs across the genome, highlights the modular nature of a host's defence portfolio, and the potential for each different step to evolve independently. We discuss how isolating the genetic basis of individual steps can help to resolve discussion over the genetic architecture of host resistance.

The Genome of Drosophila innubila Reveals Lineage-Specific Patterns of Selection in Immune Genes

Pathogenic microbes can exert extraordinary evolutionary pressure on their hosts. They can spread rapidly and sicken or even kill their host to promote their own proliferation. Because of this strong selective pressure, immune genes are some of the fastest evolving genes across metazoans, as highlighted in mammals and insects. Drosophila melanogaster serves as a powerful model for studying host/pathogen evolution. While Drosophila melanogaster are frequently exposed to various pathogens, little is known about D. melanogaster's ecology, or if they are representative of other Drosophila species in terms of pathogen pressure. Here, we characterize the genome of Drosophila innubila, a mushroom-feeding species highly diverged from D. melanogaster and investigate the evolution of the immune system. We find substantial differences in the rates of evolution of immune pathways between D. innubila and D. melanogaster. Contrasting what was previously found for D. melanogaster, we find little evidence of rapid evolution of the antiviral RNAi genes and high rates of evolution in the Toll pathway. This suggests that, while immune genes tend to be rapidly evolving in most species, the specific genes that are fastest evolving may depend either on the pathogens faced by the host and/or divergence in the basic architecture of the host's immune system.

Expanding the canon: Non-classical mosquito genes at the interface of arboviral infection

Mosquito transmitted viruses cause significant morbidity and mortality in human populations. Despite the use of insecticides and other measures of vector control, arboviral diseases are on the rise. One potential solution for limiting disease transmission to humans is to render mosquitoes refractory to viral infection through genetic modification. Substantial research effort in Drosophila, Aedes and Anopheles has helped to define the major innate immune pathways, including Toll, IMD, Jak/Stat and RNAi, however we still have an incomplete picture of the mosquito antiviral response. Transcriptional profiles of virus-infected insects reveal a much wider range of pathways activated by the process of infection. Within these lists of genes are unexplored mosquito candidates of viral defense. Wolbachia species are endosymbiotic bacteria that naturally limit arboviral infection in mosquitoes. Our understanding of the Wolbachia-mediated viral blocking mechanism is poor, but it does not appear to operate via the classical immune pathways. Herein, we reviewed the transcriptomic response of mosquitoes to multiple viral species and put forth consensus gene types/families outside the immune canon whose expression responds to infection, including cytoskeleton and cellular trafficking, the heat shock response, cytochromes P450, cell proliferation, chitin and small RNAs. We then examine emerging evidence for their functional role in viral resistance in diverse insect and mammalian hosts and their potential role in Wolbachia-mediated viral blocking. These candidate gene families offer novel avenues for research into the nature of insect viral defense.

Host-pathogen coevolution increases genetic variation in susceptibility to infection

It is common to find considerable genetic variation in susceptibility to infection in natural populations. We have investigated whether natural selection increases this variation by testing whether host populations show more genetic variation in susceptibility to pathogens that they naturally encounter than novel pathogens. In a large cross-infection experiment involving four species of Drosophila and four host-specific viruses, we always found greater genetic variation in susceptibility to viruses that had coevolved with their host. We went on to examine the genetic architecture of resistance in one host species, finding that there are more major-effect genetic variants in coevolved host-pathogen interactions. We conclude that selection by pathogens has increased genetic variation in host susceptibility, and much of this effect is caused by the occurrence of major-effect resistance polymorphisms within populations.

Interaction of RNA viruses of the natural virome with the African malaria vector, Anopheles coluzzii

Mosquitoes are colonized by a little-studied natural virome. Like the bacterial microbiome, the virome also probably influences the biology and immunity of mosquito vector populations, but tractable experimental models are lacking. We recently discovered two novel viruses in the virome of wild Anopheles and in colonies of the malaria vector Anopheles coluzzii: Anopheles C virus and Anopheles cypovirus. Here, we describe biological interactions between these two viruses and An. coluzzii mosquitoes. Viral abundance varies reproducibly during mosquito development. DNA forms of these viruses were not detected, and thus viral persistence is likely based on vertical transmission of RNA genomes. At least Anopheles C virus is vertically transmitted by an intraembryonic route. Relative abundance of the two viruses is inversely correlated in individual mosquitoes. One possible mechanism for this could be interactions with host immunity, and functional genomic analysis indicated differential influence of at least the Toll and JAK/STAT immune signaling pathways upon the viruses. The nonrandom distributions and interactions with host immunity suggest that these and other members of the natural virome may constitute a source of unrecognized heterogeneity in mosquito vector populations.

Peroxisome-associated Sgroppino links fat metabolism with survival after RNA virus infection in Drosophila

The fruit fly Drosophila melanogaster is a valuable model organism for the discovery and characterization of innate immune pathways, but host responses to virus infection remain incompletely understood. Here, we describe a novel player in host defense, Sgroppino (Sgp). Genetic depletion of Sgroppino causes hypersensitivity of adult flies to infections with the RNA viruses Drosophila C virus, cricket paralysis virus, and Flock House virus. Canonical antiviral immune pathways are functional in Sgroppino mutants, suggesting that Sgroppino exerts its activity via an as yet uncharacterized process. We demonstrate that Sgroppino localizes to peroxisomes, organelles involved in lipid metabolism. In accordance, Sgroppino-deficient flies show a defect in lipid metabolism, reflected by higher triglyceride levels, higher body mass, and thicker abdominal fat tissue. In addition, knock-down of Pex3, an essential peroxisome biogenesis factor, increases sensitivity to virus infection. Together, our results establish a genetic link between the peroxisomal protein Sgroppino, fat metabolism, and resistance to virus infection.

Parallel adaptation of rabbit populations to myxoma virus

In the 1950s the myxoma virus was released into European rabbit populations in Australia and Europe, decimating populations and resulting in the rapid evolution of resistance. We investigated the genetic basis of resistance by comparing the exomes of rabbits collected before and after the pandemic. We found a strong pattern of parallel evolution, with selection on standing genetic variation favoring the same alleles in Australia, France, and the United Kingdom. Many of these changes occurred in immunity-related genes, supporting a polygenic basis of resistance. We experimentally validated the role of several genes in viral replication and showed that selection acting on an interferon protein has increased the protein's antiviral effect.

Genetic redundancy fuels polygenic adaptation in Drosophila

The genetic architecture of adaptive traits is of key importance to predict evolutionary responses. Most adaptive traits are polygenic-i.e., result from selection on a large number of genetic loci-but most molecularly characterized traits have a simple genetic basis. This discrepancy is best explained by the difficulty in detecting small allele frequency changes (AFCs) across many contributing loci. To resolve this, we use laboratory natural selection to detect signatures for selective sweeps and polygenic adaptation. We exposed 10 replicates of a Drosophila simulans population to a new temperature regime and uncovered a polygenic architecture of an adaptive trait with high genetic redundancy among beneficial alleles. We observed convergent responses for several phenotypes-e.g., fitness, metabolic rate, and fat content-and a strong polygenic response (99 selected alleles; mean s = 0.059). However, each of these selected alleles increased in frequency only in a subset of the evolving replicates. We discerned different evolutionary paradigms based on the heterogeneous genomic patterns among replicates. Redundancy and quantitative trait (QT) paradigms fitted the experimental data better than simulations assuming independent selective sweeps. Our results show that natural D. simulans populations harbor a vast reservoir of adaptive variation facilitating rapid evolutionary responses using multiple alternative genetic pathways converging at a new phenotypic optimum. This key property of beneficial alleles requires the modification of testing strategies in natural populations beyond the search for convergence on the molecular level.

Readapting to DCV Infection without Wolbachia: Frequency Changes of Drosophila Antiviral Alleles Can Replace Endosymbiont Protection

There is now ample evidence that endosymbionts can contribute to host adaptation to environmental challenges. However, how endosymbiont presence affects the adaptive trajectory and outcome of the host is yet largely unexplored. In Drosophila, Wolbachia confers protection to RNA virus infection, an effect that differs between Wolbachia strains and can be targeted by selection. Adaptation to RNA virus infections is mediated by both Wolbachia and the host, raising the question of whether adaptive genetic changes in the host vary with the presence/absence of the endosymbiont. Here, we address this question using a polymorphic D. melanogaster population previously adapted to DCV infection for 35 generations in the presence of Wolbachia, from which we removed the endosymbiont and followed survival over the subsequent 20 generations of infection. After an initial severe drop, survival frequencies upon DCV selection increased significantly, as seen before in the presence of Wolbachia. Whole-genome sequencing, revealed that the major genes involved in the first selection experiment, pastrel and Ubc-E2H, continued to be selected in Wolbachia-free D. melanogaster, with the frequencies of protective alleles being closer to fixation in the absence of Wolbachia. Our results suggest that heterogeneity in Wolbachia infection status may be sufficient to maintain polymorphisms even in the absence of costs.

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.

Bub1 Facilitates Virus Entry through Endocytosis in a Model of Drosophila Pathogenesis

In this work, we identify for the first time that the nuclear protein Bub1 (budding uninhibited by benzimidazoles 1), a highly conserved subunit of the kinetochore complex regulating chromosome congression, has a novel and important function on the cell membrane to facilitate the virus to enter host cells. Bub1 deficiency empowers the host to have the ability to resist viral infection in Drosophila and a human cell line. Bub1 is involved in the virus entry step through regulating endocytosis. The DCV capsid protein can recruit Bub1, and DCV infection can strengthen the interaction between Bub1 and a clathrin-dependent endocytosis component. The restricted entry of vesicular stomatitis virus (VSV) and Listeria monocytogenes in bub1-deficient flies and cell lines was also observed. Therefore, our data implicate a previously unknown function of Bub1 that can be hijacked by pathogens to facilitate their entry, and Bub1 may serve as a potential antiviral therapy target for limiting viral entry.

In order to establish productive infection and dissemination, viruses usually evolve a number of strategies to hijack and/or subvert the host defense systems. However, host factors utilized by the virus to facilitate infection remain poorly characterized. In this work, we found that Drosophila melanogaster deficient in budding uninhibited by benzimidazoles 1 (bub1), a highly conserved subunit of the kinetochore complex regulating chromosome congression (1), became resistant to Drosophila C virus (DCV) infection, evidenced in increased survival rates and reduced viral loads, compared to the wild-type control. Mechanistic analysis further showed that Bub1 also functioned in the cytoplasm and was essentially involved in clathrin-dependent endocytosis of DCV and other pathogens, thus limiting pathogen entry. DCV infection potentially had strengthened the interaction between Bub1 and the clathrin adaptor on the cell membrane. Furthermore, the conserved function of Bub1 was also verified in a mammalian cell line. Thus, our data demonstrated a previously unknown function of Bub1 that could be hijacked by pathogens to facilitate their infection and spread.

IMPORTANCE In this work, we identify for the first time that the nuclear protein Bub1 (budding uninhibited by benzimidazoles 1), a highly conserved subunit of the kinetochore complex regulating chromosome congression, has a novel and important function on the cell membrane to facilitate the virus to enter host cells. Bub1 deficiency empowers the host to have the ability to resist viral infection in Drosophila and a human cell line. Bub1 is involved in the virus entry step through regulating endocytosis. The DCV capsid protein can recruit Bub1, and DCV infection can strengthen the interaction between Bub1 and a clathrin-dependent endocytosis component. The restricted entry of vesicular stomatitis virus (VSV) and Listeria monocytogenes in bub1-deficient flies and cell lines was also observed. Therefore, our data implicate a previously unknown function of Bub1 that can be hijacked by pathogens to facilitate their entry, and Bub1 may serve as a potential antiviral therapy target for limiting viral entry.

Isolation of a natural DNA virus of Drosophila melanogaster, and characterisation of host resistance and immune responses

Drosophila melanogaster has played a key role in our understanding of invertebrate immunity. However, both functional and evolutionary studies of host-virus interaction in Drosophila have been limited by a dearth of native virus isolates. In particular, despite a long history of virus research, DNA viruses of D. melanogaster have only recently been described, and none have been available for experimental study. Here we report the isolation and comprehensive characterisation of Kallithea virus, a large double-stranded DNA virus, and the first DNA virus to have been reported from wild populations of D. melanogaster. We find that Kallithea virus infection is costly for adult flies, reaching high titres in both sexes and disproportionately reducing survival in males, and movement and late fecundity in females. Using the Drosophila Genetic Reference Panel, we quantify host genetic variance for virus-induced mortality and viral titre and identify candidate host genes that may underlie this variation, including Cdc42-interacting protein 4. Using full transcriptome sequencing of infected males and females, we examine the transcriptional response of flies to Kallithea virus infection and describe differential regulation of virus-responsive genes. This work establishes Kallithea virus as a new tractable model to study the natural interaction between D. melanogaster and DNA viruses, and we hope it will serve as a basis for future studies of immune responses to DNA viruses in insects.

Defense Mechanisms against Viral Infection in Drosophila: RNAi and Non-RNAi

RNAi is considered a major antiviral defense mechanism in insects, but its relative importance as compared to other antiviral pathways has not been evaluated comprehensively. Here, it is attempted to give an overview of the antiviral defense mechanisms in Drosophila that involve both RNAi and non-RNAi. While RNAi is considered important in most viral infections, many other pathways can exist that confer antiviral resistance. It is noted that very few direct recognition mechanisms of virus infections have been identified in Drosophila and that the activation of immune pathways may be accomplished indirectly through cell damage incurred by viral replication. In several cases, protection against viral infection can be obtained in RNAi mutants by non-RNAi mechanisms, confirming the variability of the RNAi defense mechanism according to the type of infection and the physiological status of the host. This analysis is aimed at more systematically investigating the relative contribution of RNAi in the antiviral response and more specifically, to ask whether RNAi efficiency is affected when other defense mechanisms predominate. While Drosophila can function as a useful model, this issue may be more critical for economically important insects that are either controlled (agricultural pests and vectors of diseases) or protected from parasite infection (beneficial insects as bees) by RNAi products.

A Genomic Reference Panel for Drosophila serrata

Here we describe a collection of re-sequenced inbred lines of Drosophila serrata, sampled from a natural population situated deep within the species endemic distribution in Brisbane, Australia. D. serrata is a member of the speciose montium group whose members inhabit much of south east Asia and has been well studied for aspects of climatic adaptation, sexual selection, sexual dimorphism, and mate recognition. We sequenced 110 lines that were inbred via 17-20 generations of full-sib mating at an average coverage of 23.5x with paired-end Illumina reads. 15,228,692 biallelic SNPs passed quality control after being called using the Joint Genotyper for Inbred Lines (JGIL). Inbreeding was highly effective and the average levels of residual heterozygosity (0.86%) were well below theoretical expectations. As expected, linkage disequilibrium decayed rapidly, with r² dropping below 0.1 within 100 base pairs. With the exception of four closely related pairs of lines which may have been due to technical errors, there was no statistical support for population substructure. Consistent with other endemic populations of other Drosophila species, preliminary population genetic analyses revealed high nucleotide diversity and, on average, negative Tajima’s D values. A preliminary GWAS was performed on a cuticular hydrocarbon trait, 2-Me-C₂₈ revealing 4 SNPs passing Bonferroni significance residing in or near genes. One gene Cht9 may be involved in the transport of CHCs from the site of production (oenocytes) to the cuticle. Our panel will facilitate broader population genomic and quantitative genetic studies of this species and serve as an important complement to existing D. melanogaster panels that can be used to test for the conservation of genetic architectures across the Drosophila genus.