Survival Rate and Transcriptional Response upon Infection with the Generalist Parasite Beauveria bassiana in a World-Wide Sample of Drosophila melanogaster

Pathogen infection alters the gene expression landscape of transposable elements in Drosophila melanogaster

Transposable elements make up substantial proportions of eukaryotic genomes and many are thought to be remnants of ancient viral infections. Current research has begun to highlight the role transposable elements can play in the immune system response to infections. However, most of our knowledge about transposable element expression during infection is limited by the specific host and pathogen factors from each study, making it difficult to compare studies and develop broader patterns regarding the role of transposable elements during infection. Here, we use the tools and resources available in the model, Drosophila melanogaster, to analyze multiple gene expression datasets of flies subject to bacterial, fungal, and viral infections. We analyzed differences in pathogen species, host genotype, host tissue, and sex to understand how these factors impact transposable element expression during infection. Our results highlight both shared and unique transposable element expression patterns between pathogens and suggest a larger effect of pathogen factors over host factors for influencing transposable element expression.

A single amino acid polymorphism in natural Metchnikowin alleles of Drosophila results in systemic immunity and life history tradeoffs

Antimicrobial peptides (AMPs) are at the interface of interactions between hosts and microbes and are therefore expected to be rapidly evolving in a coevolutionary arms race with pathogens. In contrast, previous work demonstrated that insect AMPs tend to evolve more slowly than the genome average. Metchikowin (Mtk) is a Drosophila AMP that has a single amino acid residue that segregates as either proline (P) or arginine (R) in populations of four different species, some of which diverged more than 10 million years ago. These results suggest that there is a distinct functional importance to each allele. The most likely hypotheses are driven by two main questions: does each allele have a different efficacy against different specific pathogens (specificity hypothesis)? Or, is one allele a more potent antimicrobial, but with a host fitness cost (autoimmune hypothesis)? To assess their functional differences, we created D. melanogaster lines with the P allele, R allele, or Mtk null mutation using CRISPR/Cas9 genome editing and performed a series of life history and infection assays to assess them. In males, testing of systemic immune responses to a repertoire of bacteria and fungi demonstrated that the R allele performs as well or better than the P and null alleles with most infections. Females show some results that contrast with males, with Mtk alleles either not contributing to survival or with the P allele outperforming the R allele. In addition, measurements of life history traits demonstrate that the R allele is more costly in the absence of infection for both sexes. These results are consistent with both the specificity hypothesis (either allele can perform better against certain pathogens depending on context), and the autoimmune hypothesis (the R allele is generally the more potent antimicrobial in males, and carries a fitness cost). These results provide strong in vivo evidence that differential fitness with or without infection and sex-based functional differences in alleles may be adaptive mechanisms of maintaining immune gene polymorphisms in contrast with expectations of rapid evolution. Therefore, a complex interplay of forces including pathogen species and host sex may lead to balancing selection for immune genotypes. Strikingly, this selection may act on even a single amino acid polymorphism in an AMP.

Exploring cross-protective effects between cold and immune stress in Drosophila melanogaster

It is well established that environmental and biotic stressors like temperature and pathogens/parasites are essential for the life of small ectotherms. There are complex interactions between cold stress and pathogen infection in insects. Possible cross-protective mechanisms occur between both stressors, suggesting broad connectivity in insect stress responses. In this study, the functional significance of these interactions was tested, as well as the potential role of newly uncovered candidate genes, turandot. This was done using an array of factorial experiments exposing Drosophila melanogaster flies to a combination of different cold stress regimes (acute or chronic) and infections with the parasitic fungus Beauveria bassiana. Following these crossed treatments, phenotypic and molecular responses were assessed by measuring 1) induced cold tolerance, 2) immune resistance to parasitic fungus, and 3) activation of turandot genes. We found various responses in the phenotypic outcomes according to the various treatment combinations with higher susceptibility to infection following cold stress, but also significantly higher acute cold survival in flies that were infected. Regarding molecular responses, we found overexpression of turandot genes in response to most treatments, suggesting reactivity to both cold and infection. Moreover, maximum peak expressions were distinctly observed in the combined treatments (infection plus cold), indicating a marked synergistic effect of the stressors on turandot gene expression patterns. These results reflect the great complexity of cross-tolerance reactions between infection and abiotic stress, but could also shed light on the mechanisms underlying the activation of these responses.

Time-series gene expression patterns and their characteristics of Beauveria bassiana in the process of infecting pest insects

Beauveria bassiana has been widely used as an important biological control fungus for agricultural and forest pests, and clarifying the interaction mechanism between B. bassiana and its host will help to better exert the efficacy of the mycoinsecticide. Here, we proposed a novel pattern analysis (PA) method for analyzing time-series data and applied it to a transcriptomic data set of B. bassiana infecting Galleria mellonella. We screened out 14 patterns including 868 genes, which had some characteristics that were not inferior to differentially expressed genes (DEGs). Compared with the previous analysis of this data set, we had three novel discoveries during B. bassiana infection, including overall downregulation of gene expression, the more critical first 24 h, and enrichment of regulatory functions of downregulated genes. Our new PA method promises to be an important complement to DEGs analysis for time-series transcriptomic data, and our findings enrich our knowledge of molecular mechanisms of fungal-host interactions.

The molecular architecture of Drosophila melanogaster defense against Beauveria bassiana explored through evolve and resequence and quantitative trait locus mapping

Little is known about the genetic architecture of antifungal immunity in natural populations. Using two population genetic approaches, quantitative trait locus (QTL) mapping and evolve and resequence (E&R), we explored D. melanogaster immune defense against infection with the fungus Beauveria bassiana. The immune defense was highly variable both in the recombinant inbred lines from the Drosophila Synthetic Population Resource used for our QTL mapping and in the synthetic outbred populations used in our E&R study. Survivorship of infection improved dramatically over just 10 generations in the E&R study, and continued to increase for an additional nine generations, revealing a trade-off with uninfected longevity. Populations selected for increased defense against B. bassiana evolved cross resistance to a second, distinct B. bassiana strain but not to bacterial pathogens. The QTL mapping study revealed that sexual dimorphism in defense depends on host genotype, and the E&R study indicated that sexual dimorphism also depends on the specific pathogen to which the host is exposed. Both the QTL mapping and E&R experiments generated lists of potentially causal candidate genes, although these lists were nonoverlapping.

Spatial and temporal transcriptomic analyses reveal locust initiation of immune responses to Metarhizium acridum at the pre-penetration stage

Insect hemocyte and fat body tissues play critical functional roles in insect immunity. Little, however, is known concerning the dynamic responses of these tissues to fungal infection. Here, we report on a time course of locust hemocyte and fat body transcriptomic responses to infection by the acridid specific fungal pathogen, Metarhizium acridum. Fat body responses were more pronounced at all infection stages as compared to hemocytes. Immune and other related genes were induced far earlier than previously considered including at pre-penetration stages. Differential expression in hemocyte and fat body tissues persisted throughout the course of infection up until host death. Our data indicate selective pressure on the host to recognize the infection as early as possible in order to limit its spread. Overall, fat body and hemocyte tissues launch a robust multi-tiered response to combat the fungal pathogen, with our data providing potential host targets for exploitation in pest control.

Achilles-Mediated and Sex-Specific Regulation of Circadian mRNA Rhythms in Drosophila

The circadian clock is an evolutionarily conserved mechanism that generates the rhythmic expression of downstream genes. The core circadian clock drives the expression of clock-controlled genes, which in turn play critical roles in carrying out many rhythmic physiological processes. Nevertheless, the molecular mechanisms by which clock output genes orchestrate rhythmic signals from the brain to peripheral tissues are largely unknown. Here we explored the role of one rhythmic gene, Achilles, in regulating the rhythmic transcriptome in the fly head. Achilles is a clock-controlled gene in Drosophila that encodes a putative RNA-binding protein. Achilles expression is found in neurons throughout the fly brain using fluorescence in situ hybridization (FISH), and legacy data suggest it is not expressed in core clock neurons. Together, these observations argue against a role for Achilles in regulating the core clock. To assess its impact on circadian mRNA rhythms, we performed RNA sequencing (RNAseq) to compare the rhythmic transcriptomes of control flies and those with diminished Achilles expression in all neurons. Consistent with previous studies, we observe dramatic upregulation of immune response genes upon knock-down of Achilles. Furthermore, many circadian mRNAs lose their rhythmicity in Achilles knock-down flies, suggesting that a subset of the rhythmic transcriptome is regulated either directly or indirectly by Achilles. These Achilles-mediated rhythms are observed in genes involved in immune function and in neuronal signaling, including Prosap, Nemy and Jhl-21. A comparison of RNAseq data from control flies reveals that only 42.7% of clock-controlled genes in the fly brain are rhythmic in both males and females. As mRNA rhythms of core clock genes are largely invariant between the sexes, this observation suggests that sex-specific mechanisms are an important, and heretofore under-appreciated, regulator of the rhythmic transcriptome.

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.

Sources of variance in immunological traits: evidence of congruent latitudinal trends across species

Among-population differences in immunological traits allow assessment of both evolutionary and plastic changes in organisms' resistance to pathogens. Such knowledge also provides information necessary to predict responses of such traits to environmental changes. Studies on latitudinal trends in insect immunity have so far yielded contradictory results, suggesting that multispecies approaches with highly standardised experimental conditions are needed. Here, we studied among-population differences of two parameters reflecting constitutive immunity-phenoloxidase (PO) and lytic activity, using common-garden design on three distantly related moth species represented by populations ranging from northern Finland to Georgia (Caucasus). The larvae were reared at different temperatures and on different host plants under a crossed factors experimental design. Haemolymph samples for measurement of immune status were taken from the larvae strictly synchronously. Clear among-population differences could be shown only for PO activity in one species (elevated activity in the northern populations). There was some indication that the cases of total absence of lytic activity were more common in southern populations. The effects of temperature, host and sex on the immunological traits studied remained highly species specific. Some evidence was found that lytic activity may be involved in mediating trade-offs between immunity and larval growth performance. In contrast, PO activity rarely covaried with fitness-related traits, and neither were the values of PO and lytic activity correlated with each other. The relatively inconsistent nature of the detected patterns suggests that studies on geographic differences in immunological traits should involve multiple species, and rely on several immunological indices if general trends are a point of interest.

Leger RJ. The genetic basis for variation in resistance to infection in the Drosophila melanogaster genetic reference panel

Individuals vary extensively in the way they respond to disease but the genetic basis of this variation is not fully understood. We found substantial individual variation in resistance and tolerance to the fungal pathogen Metarhizium anisopliae Ma549 using the Drosophila melanogaster Genetic Reference Panel (DGRP). In addition, we found that host defense to Ma549 was correlated with defense to the bacterium Pseudomonas aeruginosa Pa14, and several previously published DGRP phenotypes including oxidative stress sensitivity, starvation stress resistance, hemolymph glucose levels, and sleep indices. We identified polymorphisms associated with differences between lines in both their mean survival times and microenvironmental plasticity, suggesting that lines differ in their ability to adapt to variable pathogen exposures. The majority of polymorphisms increasing resistance to Ma549 were sex biased, located in non-coding regions, had moderately large effect and were rare, suggesting that there is a general cost to defense. Nevertheless, host defense was not negatively correlated with overall longevity and fecundity. In contrast to Ma549, minor alleles were concentrated in the most Pa14-susceptible as well as the most Pa14-resistant lines. A pathway based analysis revealed a network of Pa14 and Ma549-resistance genes that are functionally connected through processes that encompass phagocytosis and engulfment, cell mobility, intermediary metabolism, protein phosphorylation, axon guidance, response to DNA damage, and drug metabolism. Functional testing with insertional mutagenesis lines indicates that 12/13 candidate genes tested influence susceptibility to Ma549. Many candidate genes have homologs identified in studies of human disease, suggesting that genes affecting variation in susceptibility are conserved across species.

Canalization of gene expression is a major signature of regulatory cold adaptation in temperate Drosophila melanogaster

BACKGROUND

Transcriptome analysis may provide means to investigate the underlying genetic causes of shared and divergent phenotypes in different populations and help to identify potential targets of adaptive evolution. Applying RNA sequencing to whole male Drosophila melanogaster from the ancestral tropical African environment and a very recently colonized cold-temperate European environment at both standard laboratory conditions and following a cold shock, we seek to uncover the transcriptional basis of cold adaptation.

RESULTS

In both the ancestral and the derived populations, the predominant characteristic of the cold shock response is the swift and massive upregulation of heat shock proteins and other chaperones. Although we find ~25 % of the genome to be differentially expressed following a cold shock, only relatively few genes (n = 16) are up- or down-regulated in a population-specific way. Intriguingly, 14 of these 16 genes show a greater degree of differential expression in the African population. Likewise, there is an excess of genes with particularly strong cold-induced changes in expression in Africa on a genome-wide scale.

CONCLUSIONS

The analysis of the transcriptional cold shock response most prominently reveals an upregulation of components of a general stress response, which is conserved over many taxa and triggered by a plethora of stressors. Despite the overall response being fairly similar in both populations, there is a definite excess of genes with a strong cold-induced fold-change in Africa. This is consistent with a detrimental deregulation or an overshooting stress response. Thus, the canalization of European gene expression might be responsible for the increased cold tolerance of European flies.

Balancing selection on immunity genes: review of the current literature and new analysis in Drosophila melanogaster

Balancing selection has been widely assumed to be an important evolutionary force, yet even today little is known about its abundance and its impact on the patterns of genetic diversity. Several studies have shown examples of balancing selection in humans, plants or parasites, and many genes under balancing selection are involved in immunity. It has been proposed that host-parasite coevolution is one of the main forces driving immune genes to evolve under balancing selection. In this paper, we review the literature on balancing selection on immunity genes in several organisms, including Drosophila. Furthermore, we performed a genome scan for balancing selection in an African population of Drosophila melanogaster using coalescent simulations of a demographic model with and without selection. We find very few genes under balancing selection and only one novel candidate gene related to immunity. Finally, we discuss the possible causes of the low number of genes under balancing selection.