Recognition of nonself is necessary to activate Drosophila's immune response against an insect parasite

BACKGROUND

Innate immune responses can be activated by pathogen-associated molecular patterns (PAMPs), danger signals released by damaged tissues, or the absence of self-molecules that inhibit immunity. As PAMPs are typically conserved across broad groups of pathogens but absent from the host, it is unclear whether they allow hosts to recognize parasites that are phylogenetically similar to themselves, such as parasitoid wasps infecting insects.

RESULTS

Parasitoids must penetrate the cuticle of Drosophila larvae to inject their eggs. In line with previous results, we found that the danger signal of wounding triggers the differentiation of specialized immune cells called lamellocytes. However, using oil droplets to mimic infection by a parasitoid wasp egg, we found that this does not activate the melanization response. This aspect of the immune response also requires exposure to parasite molecules. The unidentified factor enhances the transcriptional response in hemocytes and induces a specific response in the fat body.

CONCLUSIONS

We conclude that a combination of danger signals and the recognition of nonself molecules is required to activate Drosophila's immune response against parasitic insects.

The evolution of constitutively active humoral immune defenses in Drosophila populations under high parasite pressure

Both constitutive and inducible immune mechanisms are employed by hosts for defense against infection. Constitutive immunity allows for a faster response, but it comes with an associated cost that is always present. This trade-off between speed and fitness costs leads to the theoretical prediction that constitutive immunity will be favored where parasite exposure is frequent. We selected populations of Drosophila melanogaster under high parasite pressure from the parasitoid wasp Leptopilina boulardi. With RNA sequencing, we found the evolution of resistance in these populations was associated with them developing constitutively active humoral immunity, mediated by the larval fat body. Furthermore, these evolved populations were also able to induce gene expression in response to infection to a greater level, which indicates an overall more activated humoral immune response to parasitization. The anti-parasitoid immune response also relies on the JAK/STAT signaling pathway being activated in muscles following infection, and this induced response was only seen in populations that had evolved under high parasite pressure. We found that the cytokine Upd3, which induces this JAK/STAT response, is being expressed by immature lamellocytes. Furthermore, these immune cells became constitutively present when populations evolved resistance, potentially explaining why they gained the ability to activate JAK/STAT signaling. Thus, under intense parasitism, populations evolved resistance by increasing both constitutive and induced immune defenses, and there is likely an interplay between these two forms of immunity.

Activation of immune defences against parasitoid wasps does not underlie the cost of infection

Parasites reduce the fitness of their hosts, and different causes of this damage have fundamentally different consequences for the evolution of immune defences. Damage to the host may result from the parasite directly harming its host, often due to the production of virulence factors that manipulate host physiology. Alternatively, the host may be harmed by the activation of its own immune defences, as these can be energetically demanding or cause self-harm. A well-studied model of the cost of infection is Drosophila melanogaster and its common natural enemy, parasitoid wasps. Infected Drosophila larvae rely on humoral and cellular immune mechanisms to form a capsule around the parasitoid egg and kill it. Infection results in a developmental delay and reduced adult body size. To disentangle the effects of virulence factors and immune defences on these costs, we artificially activated anti-parasitoid immune defences in the absence of virulence factors. Despite immune activation triggering extensive differentiation and proliferation of immune cells together with hyperglycaemia, it did not result in a developmental delay or reduced body size. We conclude that the costs of infection do not result from these aspects of the immune response and may instead result from the parasite directly damaging the host.

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

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

Wolbachia protects Drosophila melanogaster against two naturally occurring and virulent viral pathogens

Wolbachia is a common endosymbiont that can protect insects against viral pathogens. However, whether the antiviral effects of Wolbachia have a significant effect on fitness remains unclear. We have investigated the interaction between Drosophila melanogaster, Wolbachia and two viruses that we recently isolated from wild flies, La Jolla virus (LJV; Iflaviridae) and Newfield virus (NFV; Permutotetraviridae). Flies infected with these viruses have increased mortality rates, and NFV partially sterilizes females. These effects on fitness were reduced in Wolbachia-infected flies, and this was associated with reduced viral titres. However, Wolbachia alone also reduces survival, and under our experimental conditions these costs of the symbiont can outweigh the benefits of antiviral protection. In contrast, protection against the sterilizing effect of NFV leads to a net benefit of Wolbachia infection after exposure to the virus. These results support the hypothesis that Wolbachia is an important defense against the natural pathogens of D. melanogaster. Furthermore, by reducing the cost of Wolbachia infection, the antiviral effects of Wolbachia may aid its invasion into populations and help explain why it is so common in nature.

Wolbachia endosymbiotic bacteria alter the gut microbiome in the fly Drosophila nigrosparsa

Wolbachia are known to cause reproductive manipulations and in some arthropod species, Wolbachia were reported to cause changes in gut microbiome. However, the effects of Wolbachia bacteria on the microbiomes of their hosts, including Drosophila flies, have not been fully accessed. Here, we checked the bacterial microbiome in guts of Wolbachia-uninfected and of Wolbachia-infected Drosophila nigrosparsa, both separated into a bleach-only (embryos bleached) and a gnotobiotic (embryos bleached and inoculated with bacteria) treatment. We observed a clear separation between the Wolbachia-infected and the Wolbachia-uninfected samples, and the infected samples had higher variation in alpha diversity than the uninfected ones. There were reductions in the abundances of Proteobacteria (Pseudomonadota), especially Acetobacter, in the infected samples of both treatments. These findings highlight that Wolbachia change the gut microbiome in D. nigrosparsa as well as that the interactions between Wolbachia and bacteria like Acetobacter need to be investigated.

Trans-regulatory changes underpin the evolution of the Drosophila immune response

When an animal is infected, the expression of a large suite of genes is changed, resulting in an immune response that can defend the host. Despite much evidence that the sequence of proteins in the immune system can evolve rapidly, the evolution of gene expression is comparatively poorly understood. We therefore investigated the transcriptional response to parasitoid wasp infection in Drosophila simulans and D. sechellia. Although these species are closely related, there has been a large scale divergence in the expression of immune-responsive genes in their two main immune tissues, the fat body and hemocytes. Many genes, including those encoding molecules that directly kill pathogens, have cis regulatory changes, frequently resulting in large differences in their expression in the two species. However, these changes in cis regulation overwhelmingly affected gene expression in immune-challenged and uninfected animals alike. Divergence in the response to infection was controlled in trans. We argue that altering trans-regulatory factors, such as signalling pathways or immune modulators, may allow natural selection to alter the expression of large numbers of immune-responsive genes in a coordinated fashion.

A fundamental question in biology is the nature of the genetic changes underlying evolutionary change, and immune systems provide an ideal system to examine this as they tend to evolve fast as animals adapt to an ever-changing array of parasites and pathogens. Comparing two species of the fruit fly Drosophila, we found that the transcriptional response to infection evolves extremely fast. However, changes in cis (where the genetic change is on the same DNA molecule as the gene in question) and trans (where the genetic change can be elsewhere in the genome) are playing different roles. Changes in cis frequently caused large differences in immune gene expression between species, but these differences were seen regardless of whether the animal was infected. In contrast, changes in trans were responsible for altering how gene expression changes in response to infection. Immune responses are complex and multifaceted, requiring the expression of many genes to be altered in a tightly regulated manner when the animal is infected. Natural selection acting on trans regulatory factors may allow the expression of many downstream genes to be altered in a coordinated fashion.

Wolbachia reduces virus infection in a natural population of Drosophila

Wolbachia is a maternally transmitted bacterial symbiont that is estimated to infect approximately half of arthropod species. In the laboratory it can increase the resistance of insects to viral infection, but its effect on viruses in nature is unknown. Here we report that in a natural population of Drosophila melanogaster, individuals that are infected with Wolbachia are less likely to be infected by viruses. By characterising the virome by metagenomic sequencing and then testing individual flies for infection, we found the protective effect of Wolbachia was virus-specific, with the prevalence of infection being up to 15% greater in Wolbachia-free flies. The antiviral effects of Wolbachia may contribute to its extraordinary ecological success, and in nature the symbiont may be an important component of the antiviral defences of insects.

DROP: Molecular voucher database for identification of Drosophila parasitoids

Molecular identification is increasingly used to speed up biodiversity surveys and laboratory experiments. However, many groups of organisms cannot be reliably identified using standard databases such as GenBank or BOLD due to lack of sequenced voucher specimens identified by experts. Sometimes a large number of sequences are available, but with too many errors to allow identification. Here, we address this problem for parasitoids of Drosophila by introducing a curated open-access molecular reference database, DROP (Drosophila parasitoids). Identifying Drosophila parasitoids is challenging and poses a major impediment to realize the full potential of this model system in studies ranging from molecular mechanisms to food webs, and in biological control of Drosophila suzukii. In DROP, genetic data are linked to voucher specimens and, where possible, the voucher specimens are identified by taxonomists and vetted through direct comparison with primary type material. To initiate DROP, we curated 154 laboratory strains, 856 vouchers, 554 DNA sequences, 16 genomes, 14 transcriptomes, and six proteomes drawn from a total of 183 operational taxonomic units (OTUs): 114 described Drosophila parasitoid species and 69 provisional species. We found species richness of Drosophila parasitoids to be heavily underestimated and provide an updated taxonomic catalogue for the community. DROP offers accurate molecular identification and improves cross-referencing between individual studies that we hope will catalyse research on this diverse and fascinating model system. Our effort should also serve as an example for researchers facing similar molecular identification problems in other groups of organisms.

Life and Death of Selfish Genes: Comparative Genomics Reveals the Dynamic Evolution of Cytoplasmic Incompatibility

Cytoplasmic incompatibility is a selfish reproductive manipulation induced by the endosymbiont Wolbachia in arthropods. In males Wolbachia modifies sperm, leading to embryonic mortality in crosses with Wolbachia-free females. In females, Wolbachia rescues the cross and allows development to proceed normally. This provides a reproductive advantage to infected females, allowing the maternally transmitted symbiont to spread rapidly through host populations. We identified homologs of the genes underlying this phenotype, cifA and cifB, in 52 of 71 new and published Wolbachia genome sequences. They are strongly associated with cytoplasmic incompatibility. There are up to seven copies of the genes in each genome, and phylogenetic analysis shows that Wolbachia frequently acquires new copies due to pervasive horizontal transfer between strains. In many cases, the genes have subsequently acquired loss-of-function mutations to become pseudogenes. As predicted by theory, this tends to occur first in cifB, whose sole function is to modify sperm, and then in cifA, which is required to rescue the cross in females. Although cif genes recombine, recombination is largely restricted to closely related homologs. This is predicted under a model of coevolution between sperm modification and embryonic rescue, where recombination between distantly related pairs of genes would create a self-incompatible strain. Together, these patterns of gene gain, loss, and recombination support evolutionary models of cytoplasmic incompatibility.

Constitutive activation of cellular immunity underlies the evolution of resistance to infection in Drosophila

Organisms rely on inducible and constitutive immune defences to combat infection. Constitutive immunity enables a rapid response to infection but may carry a cost for uninfected individuals, leading to the prediction that it will be favoured when infection rates are high. When we exposed populations of Drosophila melanogaster to intense parasitism by the parasitoid wasp Leptopilina boulardi, they evolved resistance by developing a more reactive cellular immune response. Using single-cell RNA sequencing, we found that immune-inducible genes had become constitutively upregulated. This was the result of resistant larvae differentiating precursors of specialized immune cells called lamellocytes that were previously only produced after infection. Therefore, populations evolved resistance by genetically hard-wiring the first steps of an induced immune response to become constitutive.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.

Age, tissue, genotype and virus infection regulate Wolbachia levels in Drosophila

The bacterial symbiont Wolbachia can protect insects against viral pathogens, and the varying levels of antiviral protection are correlated with the endosymbiont load within the insects. To understand why Wolbachia strains differ in their antiviral effects, we investigated the factors controlling Wolbachia density in five closely related strains in their natural Drosophila hosts. We found that Wolbachia density varied greatly across different tissues and between flies of different ages, and these effects depended on the host-symbiont association. Some endosymbionts maintained largely stable densities as flies aged while others increased, and these effects in turn depended on the tissue being examined. Measuring Wolbachia rRNA levels in response to viral infection, we found that viral infection itself also altered Wolbachia levels, with Flock House virus causing substantial reductions in symbiont loads late in the infection. This effect, however, was virus-specific as Drosophila C virus had little impact on Wolbachia in all of the five host systems. Because viruses have strong tissue tropisms and antiviral protection is thought to be cell-autonomous, these effects are likely to affect the virus-blocking phenomenon. However, we were unable to find any evidence of a correlation between Wolbachia and viral titres within the same tissues. We conclude that Wolbachia levels within flies are regulated in a complex host-symbiont-virus-dependent manner and this trinity is likely to influence the antiviral effects of Wolbachia.

Wolbachia affect behavior and possibly reproductive compatibility but not thermoresistance, fecundity, and morphology in a novel transinfected host, Drosophila nigrosparsa

Wolbachia, intracellular endosymbionts, are estimated to infect about half of all arthropod species. These bacteria manipulate their hosts in various ways for their maximum benefits. The rising global temperature may accelerate species migration, and thus, horizontal transfer of Wolbachia may occur across species previously not in contact. We transinfected and then cured the alpine fly Drosophila nigrosparsa with Wolbachia strain wMel to study its effects on this species. We found low Wolbachia titer, possibly cytoplasmic incompatibility, and an increase in locomotion of both infected larvae and adults compared with cured ones. However, no change in fecundity, no impact on heat and cold tolerance, and no change in wing morphology were observed. Although Wolbachia increased locomotor activities in this species, we conclude that D. nigrosparsa may not benefit from the infection. Still, D. nigrosparsa can serve as a host for Wolbachia because vertical transmission is possible but may not be as high as in the native host of wMel, Drosophila melanogaster.

Virus evolution in Wolbachia-infected Drosophila

Wolbachia, a common vertically transmitted symbiont, can protect insects against viral infection and prevent mosquitoes from transmitting viral pathogens. For this reason, Wolbachia-infected mosquitoes are being released to prevent the transmission of dengue and other arboviruses. An important question for the long-term success of these programmes is whether viruses can evolve to escape the antiviral effects of Wolbachia. We have found that Wolbachia altered the outcome of competition between strains of the DCV virus in Drosophila. However, Wolbachia still effectively blocked the virus genotypes that were favoured in the presence of the symbiont. We conclude that Wolbachia did cause an evolutionary response in viruses, but this has little or no impact on the effectiveness of virus blocking.

Independent effects on cellular and humoral immune responses underlie genotype-by-genotype interactions between Drosophila and parasitoids

It is common to find abundant genetic variation in host resistance and parasite infectivity within populations, with the outcome of infection frequently depending on genotype-specific interactions. Underlying these effects are complex immune defenses that are under the control of both host and parasite genes. We have found extensive variation in Drosophila melanogaster's immune response against the parasitoid wasp Leptopilina boulardi. Some aspects of the immune response, such as phenoloxidase activity, are predominantly affected by the host genotype. Some, such as upregulation of the complement-like protein Tep1, are controlled by the parasite genotype. Others, like the differentiation of immune cells called lamellocytes, depend on the specific combination of host and parasite genotypes. These observations illustrate how the outcome of infection depends on independent genetic effects on different aspects of host immunity. As parasite-killing results from the concerted action of different components of the immune response, these observations provide a physiological mechanism to generate phenomena like epistasis and genotype-interactions that underlie models of coevolution.

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.

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.

Adaptive introgression underlies polymorphic seasonal camouflage in snowshoe hares

Snowshoe hares (Lepus americanus) maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation controlling seasonal expression of the Agouti gene underlies this adaptive winter camouflage polymorphism. Genetic variation at Agouti clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.

Vertically transmitted rhabdoviruses are found across three insect families and have dynamic interactions with their hosts

A small number of free-living viruses have been found to be obligately vertically transmitted, but it remains uncertain how widespread vertically transmitted viruses are and how quickly they can spread through host populations. Recent metagenomic studies have found several insects to be infected with sigma viruses (Rhabdoviridae). Here, we report that sigma viruses that infect Mediterranean fruit flies (Ceratitis capitata), Drosophila immigrans, and speckled wood butterflies (Pararge aegeria) are all vertically transmitted. We find patterns of vertical transmission that are consistent with those seen in Drosophila sigma viruses, with high rates of maternal transmission, and lower rates of paternal transmission. This mode of transmission allows them to spread rapidly in populations, and using viral sequence data we found the viruses in D. immigrans and C. capitata had both recently swept through host populations. The viruses were common in nature, with mean prevalences of 12% in C. capitata, 38% in D. immigrans and 74% in P. aegeria. We conclude that vertically transmitted rhabdoviruses may be widespread in a broad range of insect taxa, and that these viruses can have dynamic interactions with their hosts.

Pan-arthropod analysis reveals somatic piRNAs as an ancestral defence against transposable elements

In animals, small RNA molecules termed PIWI-interacting RNAs (piRNAs) silence transposable elements (TEs), protecting the germline from genomic instability and mutation. piRNAs have been detected in the soma in a few animals, but these are believed to be specific adaptations of individual species. Here, we report that somatic piRNAs were probably present in the ancestral arthropod more than 500 million years ago. Analysis of 20 species across the arthropod phylum suggests that somatic piRNAs targeting TEs and messenger RNAs are common among arthropods. The presence of an RNA-dependent RNA polymerase in chelicerates (horseshoe crabs, spiders and scorpions) suggests that arthropods originally used a plant-like RNA interference mechanism to silence TEs. Our results call into question the view that the ancestral role of the piRNA pathway was to protect the germline and demonstrate that small RNA silencing pathways have been repurposed for both somatic and germline functions throughout arthropod evolution.

Addicted? Reduced host resistance in populations with defensive symbionts

Heritable symbionts that protect their hosts from pathogens have been described in a wide range of insect species. By reducing the incidence or severity of infection, these symbionts have the potential to reduce the strength of selection on genes in the insect genome that increase resistance. Therefore, the presence of such symbionts may slow down the evolution of resistance. Here we investigated this idea by exposing Drosophila melanogaster populations to infection with the pathogenic Drosophila C virus (DCV) in the presence or absence of Wolbachia, a heritable symbiont of arthropods that confers protection against viruses. After nine generations of selection, we found that resistance to DCV had increased in all populations. However, in the presence of Wolbachia the resistant allele of pastrel-a gene that has a major effect on resistance to DCV-was at a lower frequency than in the symbiont-free populations. This finding suggests that defensive symbionts have the potential to hamper the evolution of insect resistance genes, potentially leading to a state of evolutionary addiction where the genetically susceptible insect host mostly relies on its symbiont to fight pathogens.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L)

Background

Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex F_ST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings.

Results

Intersex standardized variance (F_ST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with F_ST values ≥ 0.100 and/or with increased expected heterozygosity (H_exp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex F_ST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex F_ST genes proximal to the M locus, whereas SenAae had inter-sex F_ST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H_exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation.

Conclusions

The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4348-4) contains supplementary material, which is available to authorized users.

Symbiont strain is the main determinant of variation in Wolbachia-mediated protection against viruses across Drosophila species

Wolbachia is a common heritable bacterial symbiont in insects. Its evolutionary success lies in the diverse phenotypic effects it has on its hosts coupled to its propensity to move between host species over evolutionary timescales. In a survey of natural host-symbiont associations in a range of Drosophila species, we found that 10 of 16 Wolbachia strains protected their hosts against viral infection. By moving Wolbachia strains between host species, we found that the symbiont genome had a much greater influence on the level of antiviral protection than the host genome. The reason for this was that the level of protection depended on the density of the symbiont in host tissues, and Wolbachia rather than the host-controlled density. The finding that virus resistance and symbiont density are largely under the control of symbiont genes in this system has important implications both for the evolution of these traits and for public health programmes using Wolbachia to prevent mosquitoes from transmitting disease.

Population genomics reveals that an anthropophilic population of Aedes aegypti mosquitoes in West Africa recently gave rise to American and Asian populations of this major disease vector

BACKGROUND

The mosquito Aedes aegypti is the main vector of dengue, Zika, chikungunya and yellow fever viruses. This major disease vector is thought to have arisen when the African subspecies Ae. aegypti formosus evolved from being zoophilic and living in forest habitats into a form that specialises on humans and resides near human population centres. The resulting domestic subspecies, Ae. aegypti aegypti, is found throughout the tropics and largely blood-feeds on humans.

RESULTS

To understand this transition, we have sequenced the exomes of mosquitoes collected from five populations from around the world. We found that Ae. aegypti specimens from an urban population in Senegal in West Africa were more closely related to populations in Mexico and Sri Lanka than they were to a nearby forest population. We estimate that the populations in Senegal and Mexico split just a few hundred years ago, and we found no evidence of Ae. aegypti aegypti mosquitoes migrating back to Africa from elsewhere in the tropics. The out-of-Africa migration was accompanied by a dramatic reduction in effective population size, resulting in a loss of genetic diversity and rare genetic variants.

CONCLUSIONS

We conclude that a domestic population of Ae. aegypti in Senegal and domestic populations on other continents are more closely related to each other than to other African populations. This suggests that an ancestral population of Ae. aegypti evolved to become a human specialist in Africa, giving rise to the subspecies Ae. aegypti aegypti. The descendants of this population are still found in West Africa today, and the rest of the world was colonised when mosquitoes from this population migrated out of Africa. This is the first report of an African population of Ae. aegypti aegypti mosquitoes that is closely related to Asian and American populations. As the two subspecies differ in their ability to vector disease, their existence side by side in West Africa may have important implications for disease transmission.

Latitudinal clines in gene expression and cis-regulatory element variation in Drosophila melanogaster

BACKGROUND

Organisms can rapidly adapt to their environment when colonizing a new habitat, and this could occur by changing protein sequences or by altering patterns of gene expression. The importance of gene expression in driving local adaptation is increasingly being appreciated, and cis-regulatory elements (CREs), which control and modify the expression of the nearby genes, are predicted to play an important role. Here we investigate genetic variation in gene expression in immune-challenged Drosophila melanogaster from temperate and tropical or sub-tropical populations in Australia and United States.

RESULTS

We find parallel latitudinal changes in gene expression, with genes involved in immunity, insecticide resistance, reproduction, and the response to the environment being especially likely to differ between latitudes. By measuring allele-specific gene expression (ASE), we show that cis-regulatory variation also shows parallel latitudinal differences between the two continents and contributes to the latitudinal differences in gene expression.

CONCLUSIONS

Both Australia and United States were relatively recently colonized by D. melanogaster, and it was recently shown that introductions of both African and European flies occurred, with African genotypes contributing disproportionately to tropical populations. Therefore, both the demographic history of the populations and local adaptation may be causing the patterns that we see.

The genetic architecture of resistance to virus infection in Drosophila

Variation in susceptibility to infection has a substantial genetic component in natural populations, and it has been argued that selection by pathogens may result in it having a simpler genetic architecture than many other quantitative traits. This is important as models of host-pathogen co-evolution typically assume resistance is controlled by a small number of genes. Using the Drosophila melanogaster multiparent advanced intercross, we investigated the genetic architecture of resistance to two naturally occurring viruses, the sigma virus and DCV (Drosophila C virus). We found extensive genetic variation in resistance to both viruses. For DCV resistance, this variation is largely caused by two major-effect loci. Sigma virus resistance involves more genes - we mapped five loci, and together these explained less than half the genetic variance. Nonetheless, several of these had a large effect on resistance. Models of co-evolution typically assume strong epistatic interactions between polymorphisms controlling resistance, but we were only able to detect one locus that altered the effect of the main effect loci we had mapped. Most of the loci we mapped were probably at an intermediate frequency in natural populations. Overall, our results are consistent with major-effect genes commonly affecting susceptibility to infectious diseases, with DCV resistance being a near-Mendelian trait.

Genome Sequencing of the Phytoseiid Predatory Mite Metaseiulus occidentalis Reveals Completely Atomized Hox Genes and Superdynamic Intron Evolution

Metaseiulus occidentalis is an eyeless phytoseiid predatory mite employed for the biological control of agricultural pests including spider mites. Despite appearances, these predator and prey mites are separated by some 400 Myr of evolution and radically different lifestyles. We present a 152-Mb draft assembly of the M. occidentalis genome: Larger than that of its favored prey, Tetranychus urticae, but considerably smaller than those of many other chelicerates, enabling an extremely contiguous and complete assembly to be built-the best arachnid to date. Aided by transcriptome data, genome annotation cataloged 18,338 protein-coding genes and identified large numbers of Helitron transposable elements. Comparisons with other arthropods revealed a particularly dynamic and turbulent genomic evolutionary history. Its genes exhibit elevated molecular evolution, with strikingly high numbers of intron gains and losses, in stark contrast to the deer tick Ixodes scapularis Uniquely among examined arthropods, this predatory mite's Hox genes are completely atomized, dispersed across the genome, and it encodes five copies of the normally single-copy RNA processing Dicer-2 gene. Examining gene families linked to characteristic biological traits of this tiny predator provides initial insights into processes of sex determination, development, immune defense, and how it detects, disables, and digests its prey. As the first reference genome for the Phytoseiidae, and for any species with the rare sex determination system of parahaploidy, the genome of the western orchard predatory mite improves genomic sampling of chelicerates and provides invaluable new resources for functional genomic analyses of this family of agriculturally important mites.

Wolbachia Blocks Viral Genome Replication Early in Infection without a Transcriptional Response by the Endosymbiont or Host Small RNA Pathways

The intracellular endosymbiotic bacterium Wolbachia can protect insects against viral infection, and is being introduced into mosquito populations in the wild to block the transmission of arboviruses that infect humans and are a major public health concern. To investigate the mechanisms underlying this antiviral protection, we have developed a new model system combining Wolbachia-infected Drosophila melanogaster cell culture with the model mosquito-borne Semliki Forest virus (SFV; Togaviridae, Alphavirus). Wolbachia provides strong antiviral protection rapidly after infection, suggesting that an early stage post-infection is being blocked. Wolbachia does appear to have major effects on events distinct from entry, assembly or exit as it inhibits the replication of an SFV replicon transfected into the cells. Furthermore, it causes a far greater reduction in the expression of proteins from the 3´ open reading frame than the 5´ non-structural protein open reading frame, indicating that it is blocking the replication of viral RNA. Further to this separation of the replicase proteins and viral RNA in transreplication assays shows that uncoupling of viral RNA and replicase proteins does not overcome Wolbachia’s antiviral activity. This further suggests that replicative processes are disrupted, such as translation or replication, by Wolbachia infection. This may occur by Wolbachia mounting an active antiviral response, but the virus did not cause any transcriptional response by the bacterium, suggesting that this is not the case. Host microRNAs (miRNAs) have been implicated in protection, but again we found that host cell miRNA expression was unaffected by the bacterium and neither do our findings suggest any involvement of the antiviral siRNA pathway. We conclude that Wolbachia may directly interfere with early events in virus replication such as translation of incoming viral RNA or RNA transcription, and this likely involves an intrinsic (as opposed to an induced) mechanism.

The intracellular endosymbiotic bacterium Wolbachia can protect insects against viral infection. However, the mechanisms underlying this antiviral activity are poorly understood. We have developed a new model system combining Wolbachia-infected Drosophila melanogaster cell culture and the model mosquito-borne virus, Semliki Forest virus. Wolbachia confers strong antiviral activity against SFV. Our study indicates that viral replication appears to be inhibited at a very early stage, such as initial translation or replication. Results indicate that Wolbachia does not mount a transcriptional response to SFV infection and that host small RNA pathways are not involved in Wolbachia mediated antiviral activity in our system. We conclude that Wolbachia may directly interfere with early events in virus replication such as translation of incoming viral RNA or RNA transcription, and this likely involves an intrinsic (as opposed to an induced) mechanism.

A gene associated with social immunity in the burying beetle Nicrophorus vespilloides

Some group-living species exhibit social immunity, where the immune response of one individual can protect others in the group from infection. In burying beetles, this is part of parental care. Larvae feed on vertebrate carcasses which their parents smear with exudates that inhibit microbial growth. We have sequenced the transcriptome of the burying beetle Nicrophorus vespilloides and identified six genes that encode lysozymes—a type of antimicrobial enzyme that has previously been implicated in social immunity in burying beetles. When females start breeding and producing antimicrobial anal exudates, we found that the expression of one of these genes was increased by approximately 1000 times to become one of the most abundant transcripts in the transcriptome. Females varied considerably in the antimicrobial properties of their anal exudates, and this was strongly correlated with the expression of this lysozyme. We conclude that we have likely identified a gene encoding a key effector molecule in social immunity and that it was recruited during evolution from a function in personal immunity.

A Polymorphism in the Processing Body Component Ge-1 Controls Resistance to a Naturally Occurring Rhabdovirus in Drosophila

Hosts encounter an ever-changing array of pathogens, so there is continual selection for novel ways to resist infection. A powerful way to understand how hosts evolve resistance is to identify the genes that cause variation in susceptibility to infection. Using high-resolution genetic mapping we have identified a naturally occurring polymorphism in a gene called Ge-1 that makes Drosophila melanogaster highly resistant to its natural pathogen Drosophila melanogaster sigma virus (DMelSV). By modifying the sequence of the gene in transgenic flies, we identified a 26 amino acid deletion in the serine-rich linker region of Ge-1 that is causing the resistance. Knocking down the expression of the susceptible allele leads to a decrease in viral titre in infected flies, indicating that Ge-1 is an existing restriction factor whose antiviral effects have been increased by the deletion. Ge-1 plays a central role in RNA degradation and the formation of processing bodies (P bodies). A key effector in antiviral immunity, the RNAi induced silencing complex (RISC), localises to P bodies, but we found that Ge-1-based resistance is not dependent on the small interfering RNA (siRNA) pathway. However, we found that Decapping protein 1 (DCP1) protects flies against sigma virus. This protein interacts with Ge-1 and commits mRNA for degradation by removing the 5’ cap, suggesting that resistance may rely on this RNA degradation pathway. The serine-rich linker domain of Ge-1 has experienced strong selection during the evolution of Drosophila, suggesting that this gene may be under long-term selection by viruses. These findings demonstrate that studying naturally occurring polymorphisms that increase resistance to infections enables us to identify novel forms of antiviral defence, and support a pattern of major effect polymorphisms controlling resistance to viruses in Drosophila.

Hosts and their pathogens are engaged in a never-ending arms race, and hosts must continually evolve new defences to protect themselves from infection. In the fruit fly Drosophila melanogaster we show that virus resistance can evolve through a single mutation. In flies that are highly resistant to a naturally occurring virus called sigma virus we identified a deletion in the protein-coding region of a gene called Ge-1. We experimentally confirmed that this was the cause of resistance by deleting this region in transgenic flies. Furthermore, we show that even the susceptible allele of Ge-1 helps protect flies against the virus, suggesting that this mutation has made an existing antiviral defence more effective. Ge-1 plays a central role in RNA degradation in regions of the cytoplasm called P bodies, and our results suggest that this pathway has been recruited during evolution to protect D. melanogaster against sigma virus. The protein domain that contains the deletion has experienced strong selection during its evolution, suggesting that it may be involved in an ongoing arms race with viruses.

The evolution, diversity, and host associations of rhabdoviruses

Metagenomic studies are leading to the discovery of a hidden diversity of RNA viruses. These new viruses are poorly characterized and new approaches are needed predict the host species these viruses pose a risk to. The rhabdoviruses are a diverse family of RNA viruses that includes important pathogens of humans, animals, and plants. We have discovered thirty-two new rhabdoviruses through a combination of our own RNA sequencing of insects and searching public sequence databases. Combining these with previously known sequences we reconstructed the phylogeny of 195 rhabdovirus sequences, and produced the most in depth analysis of the family to date. In most cases we know nothing about the biology of the viruses beyond the host they were identified from, but our dataset provides a powerful phylogenetic approach to predict which are vector-borne viruses and which are specific to vertebrates or arthropods. By reconstructing ancestral and present host states we found that switches between major groups of hosts have occurred rarely during rhabdovirus evolution. This allowed us to propose seventy-six new likely vector-borne vertebrate viruses among viruses identified from vertebrates or biting insects. Based on currently available data, our analysis suggests it is likely there was a single origin of the known plant viruses and arthropod-borne vertebrate viruses, while vertebrate- and arthropod-specific viruses arose at least twice. There are also few transitions between aquatic and terrestrial ecosystems. Viruses also cluster together at a finer scale, with closely related viruses tending to be found in closely related hosts. Our data therefore suggest that throughout their evolution, rhabdoviruses have occasionally jumped between distantly related host species before spreading through related hosts in the same environment. This approach offers a way to predict the most probable biology and key traits of newly discovered viruses.

Comparative Genomics Reveals the Origins and Diversity of Arthropod Immune Systems

Insects are an important model for the study of innate immune systems, but remarkably little is known about the immune system of other arthropod groups despite their importance as disease vectors, pests, and components of biological diversity. Using comparative genomics, we have characterized the immune system of all the major groups of arthropods beyond insects for the first time--studying five chelicerates, a myriapod, and a crustacean. We found clear traces of an ancient origin of innate immunity, with some arthropods having Toll-like receptors and C3-complement factors that are more closely related in sequence or structure to vertebrates than other arthropods. Across the arthropods some components of the immune system, such as the Toll signaling pathway, are highly conserved. However, there is also remarkable diversity. The chelicerates apparently lack the Imd signaling pathway and beta-1,3 glucan binding proteins--a key class of pathogen recognition receptors. Many genes have large copy number variation across species, and this may sometimes be accompanied by changes in function. For example, we find that peptidoglycan recognition proteins have frequently lost their catalytic activity and switch between secreted and intracellular forms. We also find that there has been widespread and extensive duplication of the cellular immune receptor Dscam (Down syndrome cell adhesion molecule), which may be an alternative way to generate the high diversity produced by alternative splicing in insects. In the antiviral short interfering RNAi pathway Argonaute 2 evolves rapidly and is frequently duplicated, with a highly variable copy number. Our results provide a detailed analysis of the immune systems of several important groups of animals for the first time and lay the foundations for functional work on these groups.

Should Symbionts Be Nice or Selfish? Antiviral Effects of Wolbachia Are Costly but Reproductive Parasitism Is Not

Symbionts can have mutualistic effects that increase their host's fitness and/or parasitic effects that reduce it. Which of these strategies evolves depends in part on the balance of their costs and benefits to the symbiont. We have examined these questions in Wolbachia, a vertically transmitted endosymbiont of insects that can provide protection against viral infection and/or parasitically manipulate its hosts' reproduction. Across multiple symbiont strains we find that the parasitic phenotype of cytoplasmic incompatibility and antiviral protection are uncorrelated. Strong antiviral protection is associated with substantial reductions in other fitness-related traits, whereas no such trade-off was detected for cytoplasmic incompatibility. The reason for this difference is likely that antiviral protection requires high symbiont densities but cytoplasmic incompatibility does not. These results are important for the use of Wolbachia to block dengue virus transmission by mosquitoes, as natural selection to reduce these costs may lead to reduced symbiont density and the loss of antiviral protection.

Exome and transcriptome sequencing of Aedes aegypti identifies a locus that confers resistance to Brugia malayi and alters the immune response

Many mosquito species are naturally polymorphic for their abilities to transmit parasites, a feature which is of great interest for controlling vector-borne disease. Aedes aegypti, the primary vector of dengue and yellow fever and a laboratory model for studying lymphatic filariasis, is genetically variable for its capacity to harbor the filarial nematode Brugia malayi. The genome of Ae. aegypti is large and repetitive, making genome resequencing difficult and expensive. We designed exome captures to target protein-coding regions of the genome, and used association mapping in a wild Kenyan population to identify a single, dominant, sex-linked locus underlying resistance. This falls in a region of the genome where a resistance locus was previously mapped in a line established in 1936, suggesting that this polymorphism has been maintained in the wild for the at least 80 years. We then crossed resistant and susceptible mosquitoes to place both alleles of the gene into a common genetic background, and used RNA-seq to measure the effect of this locus on gene expression. We found evidence for Toll, IMD, and JAK-STAT pathway activity in response to early stages of B. malayi infection when the parasites are beginning to die in the resistant genotype. We also found that resistant mosquitoes express anti-microbial peptides at the time of parasite-killing, and that this expression is suppressed in susceptible mosquitoes. Together, we have found that a single resistance locus leads to a higher immune response in resistant mosquitoes, and we identify genes in this region that may be responsible for this trait.

The causes and consequences of changes in virulence following pathogen host shifts

Emerging infectious diseases are often the result of a host shift, where the pathogen originates from a different host species. Virulence--the harm a pathogen does to its host-can be extremely high following a host shift (for example Ebola, HIV, and SARs), while other host shifts may go undetected as they cause few symptoms in the new host. Here we examine how virulence varies across host species by carrying out a large cross infection experiment using 48 species of Drosophilidae and an RNA virus. Host shifts resulted in dramatic variation in virulence, with benign infections in some species and rapid death in others. The change in virulence was highly predictable from the host phylogeny, with hosts clustering together in distinct clades displaying high or low virulence. High levels of virulence are associated with high viral loads, and this may determine the transmission rate of the virus.

The intracellular bacterium Wolbachia uses parasitoid wasps as phoretic vectors for efficient horizontal transmission

Facultative bacterial endosymbionts are associated with many arthropods and are primarily transmitted vertically from mother to offspring. However, phylogenetic affiliations suggest that horizontal transmission must also occur. Such horizontal transfer can have important biological and agricultural consequences when endosymbionts increase host fitness. So far horizontal transmission is considered rare and has been difficult to document. Here, we use fluorescence in situ hybridization (FISH) and multi locus sequence typing (MLST) to reveal a potentially common pathway of horizontal transmission of endosymbionts via parasitoids of insects. We illustrate that the mouthparts and ovipositors of an aphelinid parasitoid become contaminated with Wolbachia when this wasp feeds on or probes Wolbachia-infected Bemisia tabaci AsiaII7, and non-lethal probing of uninfected B. tabaci AsiaII7 nymphs by parasitoids carrying Wolbachia resulted in newly and stably infected B. tabaci matrilines. After they were exposed to infected whitefly, the parasitoids were able to transmit Wolbachia efficiently for the following 48 h. Whitefly infected with Wolbachia by parasitoids had increased survival and reduced development times. Overall, our study provides evidence for the horizontal transmission of Wolbachia between insect hosts by parasitic wasps, and the enhanced survival and reproductive abilities of insect hosts may adversely affect biological control programs.

Environmental and genetic factors determine whether the mosquito Aedes aegypti lays eggs without a blood meal

Some mosquito strains or species are able to lay eggs without taking a blood meal, a trait named autogeny. This may allow populations to persist through times or places where vertebrate hosts are scarce. Autogenous egg production is highly dependent on the environment in some species, but the ideal conditions for its expression in Aedes aegypti mosquitoes are unknown. We found that 3.2% of females in a population of Ae. aegypti from Kenya were autogenous. Autogeny was strongly influenced by temperature, with many more eggs laid at 28°C compared with 22°C. Good nutrition in larval stages and feeding on higher concentrations of sugar solution during the adult stage both result in more autogenous eggs being produced. The trait also has a genetic basis, as not all Ae. aegypti genotypes can lay autogenously. We conclude that Ae. aegypti requires a favorable environment and a suitable genotype to be able to lay eggs without a blood meal.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

Vector competence of Aedes aegypti mosquitoes for filarial nematodes is affected by age and nutrient limitation

Mosquitoes are one of the most important vectors of human disease. The ability of mosquitoes to transmit disease is dependent on the age structure of the population, as mosquitoes must survive long enough for the parasites to complete their development and infect another human. Age could have additional effects due to mortality rates and vector competence changing as mosquitoes senesce, but these are comparatively poorly understood. We have investigated these factors using the mosquito Aedes aegypti and the filarial nematode Brugia malayi. Rather than observing any effects of immune senescence, we found that older mosquitoes were more resistant, but this only occurred if they had previously been maintained on a nutrient-poor diet of fructose. Constant blood feeding reversed this decline in vector competence, meaning that the number of parasites remained relatively unchanged as mosquitoes aged. Old females that had been maintained on fructose also experienced a sharp spike in mortality after an infected blood meal ("refeeding syndrome") and few survived long enough for the parasite to develop. Again, this effect was prevented by frequent blood meals. Our results indicate that old mosquitoes may be inefficient vectors due to low vector competence and high mortality, but that frequent blood meals can prevent these effects of age.

The evolution and genetics of virus host shifts

Emerging viral diseases are often the product of a host shift, where a pathogen jumps from its original host into a novel species. Phylogenetic studies show that host shifts are a frequent event in the evolution of most pathogens, but why pathogens successfully jump between some host species but not others is only just becoming clear. The susceptibility of potential new hosts can vary enormously, with close relatives of the natural host typically being the most susceptible. Often, pathogens must adapt to successfully infect a novel host, for example by evolving to use different cell surface receptors, to escape the immune response, or to ensure they are transmitted by the new host. In viruses there are often limited molecular solutions to achieve this, and the same sequence changes are often seen each time a virus infects a particular host. These changes may come at a cost to other aspects of the pathogen's fitness, and this may sometimes prevent host shifts from occurring. Here we examine how these evolutionary factors affect patterns of host shifts and disease emergence.

Estimates of allele-specific expression in Drosophila with a single genome sequence and RNA-seq data

MOTIVATION

Genetic variation in cis-regulatory elements is an important cause of variation in gene expression. Cis-regulatory variation can be detected by using high-throughput RNA sequencing (RNA-seq) to identify differences in the expression of the two alleles of a gene. This requires that reads from the two alleles are equally likely to map to a reference genome(s), and that single-nucleotide polymorphisms (SNPs) are accurately called, so that reads derived from the different alleles can be identified. Both of these prerequisites can be achieved by sequencing the genomes of the parents of the individual being studied, but this is often prohibitively costly.

RESULTS

In Drosophila, we demonstrate that biases during read mapping can be avoided by mapping reads to two alternative genomes that incorporate SNPs called from the RNA-seq data. The SNPs can be reliably called from the RNA-seq data itself, provided any variants not found in high-quality SNP databases are filtered out. Finally, we suggest a way of measuring allele-specific expression (ASE) by crossing the line of interest to a reference line with a high-quality genome sequence. Combined with our bioinformatic methods, this approach minimizes mapping biases, allows poor-quality data to be identified and removed and aides in the biological interpretation of the data as the parent of origin of each allele is known. In conclusion, our results suggest that accurate estimates of ASE do not require the parental genomes of the individual being studied to be sequenced.

AVAILABILITY AND IMPLEMENTATION

Scripts used to perform our analysis are available at https://github.com/d-quinn/bio_quinn2013.

Flies on the move: an inherited virus mirrors Drosophila melanogaster's elusive ecology and demography

Vertically transmitted parasites rely on their host's reproduction for their transmission, leading to the evolutionary histories of both parties being intimately entwined. Parasites can thus serve as a population genetic magnifying glass for their host's demographic history. Here, we study the fruitfly Drosophila melanogaster's vertically transmitted sigma virus DMelSV. The virus has a high mutation rate and low effective population size, allowing us to reconstruct at a fine scale how the combined forces of the movement of flies and selection on the virus have shaped its migration patterns. We found that the virus is likely to have spread to Europe from Africa, mirroring the colonization route of Drosophila. The North American DMelSV population appears to be the result of a recent single immigration from Europe, invading together with its host in the late 19th century. Across Europe, DMelSV migration rates are low and populations are highly genetically structured, likely reflecting limited fly movement. Despite being intolerant of extreme cold, viral diversity suggests that fly populations can persist in harsh continental climates and that recolonization from the warmer south plays a minor role. In conclusion, studying DMelSV can provide insights into the poorly understood ecology of D. melanogaster, one of the best-studied organisms in biology.

Wolbachia variants induce differential protection to viruses in Drosophila melanogaster: a phenotypic and phylogenomic analysis

Wolbachia are intracellular bacterial symbionts that are able to protect various insect hosts from viral infections. This tripartite interaction was initially described in Drosophila melanogaster carrying wMel, its natural Wolbachia strain. wMel has been shown to be genetically polymorphic and there has been a recent change in variant frequencies in natural populations. We have compared the antiviral protection conferred by different wMel variants, their titres and influence on host longevity, in a genetically identical D. melanogaster host. The phenotypes cluster the variants into two groups--wMelCS-like and wMel-like. wMelCS-like variants give stronger protection against Drosophila C virus and Flock House virus, reach higher titres and often shorten the host lifespan. We have sequenced and assembled the genomes of these Wolbachia, and shown that the two phenotypic groups are two monophyletic groups. We have also analysed a virulent and over-replicating variant, wMelPop, which protects D. melanogaster even better than the closely related wMelCS. We have found that a ~21 kb region of the genome, encoding eight genes, is amplified seven times in wMelPop and may be the cause of its phenotypes. Our results indicate that the more protective wMelCS-like variants, which sometimes have a cost, were replaced by the less protective but more benign wMel-like variants. This has resulted in a recent reduction in virus resistance in D. melanogaster in natural populations worldwide. Our work helps to understand the natural variation in wMel and its evolutionary dynamics, and inform the use of Wolbachia in arthropod-borne disease control.

Previous exposure to an RNA virus does not protect against subsequent infection in Drosophila melanogaster

BACKGROUND

Immune priming has been shown to occur in a wide array of invertebrate taxa, with individuals exposed to a pathogen showing increased protection upon subsequent exposure. However, the mechanisms underlying immune priming are poorly understood. The antiviral RNAi response in Drosophila melanogaster is an ideal candidate for providing a specific and acquired response to subsequent infection. We exposed D. melanogaster to two challenges of a virus known to produce an antiviral RNAi response, to examine whether any protective effects of prior exposure on survival were observed.

RESULTS

In this experiment we found no evidence that prior exposure to Drosophila C Virus (DCV) protects flies from a subsequent lethal challenge, with almost identical levels of mortality in flies previously exposed to DCV or a control.

CONCLUSIONS

Our results confirm the finding that 'acquired' immune responses are not ubiquitous across all invertebrate-pathogen interactions. We discuss why we may have observed no effect in this study, with focus on the mechanistic basis of the RNAi pathway.

Identification of Wolbachia strains in mosquito disease vectors

Wolbachia bacteria are common endosymbionts of insects, and some strains are known to protect their hosts against RNA viruses and other parasites. This has led to the suggestion that releasing Wolbachia-infected mosquitoes could prevent the transmission of arboviruses and other human parasites. We have identified Wolbachia in Kenyan populations of the yellow fever vector Aedes bromeliae and its relative Aedes metallicus, and in Mansonia uniformis and Mansonia africana, which are vectors of lymphatic filariasis. These Wolbachia strains cluster together on the bacterial phylogeny, and belong to bacterial clades that have recombined with other unrelated strains. These new Wolbachia strains may be affecting disease transmission rates of infected mosquito species, and could be transferred into other mosquito vectors as part of control programs.

The dynamics of reciprocal selective sweeps of host resistance and a parasite counter-adaptation in Drosophila

Host-parasite coevolution can result in consecutive selective sweeps of host resistance alleles and parasite counter-adaptations. To illustrate the dynamics of this important but little studied form of coevolution, we have modeled an ongoing arms race between Drosophila melanogaster and the vertically transmitted sigma virus, using parameters we estimated in the field. We integrate these results with previous work showing that the spread of a resistance allele of the ref(2)P gene in the host was followed by the spread of a virus genotype, which overcomes this resistance. In line with these observations, our model predicts that there can be rapid selective sweeps in both the host and parasite, which can drive large changes in the prevalence of infection. The virus will tend to be ahead in the arms race, as incomplete dominance slows down host adaptation and selection for host resistance is weaker than selection for parasites to overcome resistance--the "life-dinner" principle. This asymmetry in the adaptation rates results in a partial sweep of the host resistance allele, as it loses its advantage part way through the selective sweep. This well-understood natural system illustrates how the outcome of host-parasite coevolution is determined by different population genetic parameters in the field.

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

Variation in susceptibility to infectious disease often has a substantial genetic component in animal and plant populations. We have used genome-wide association studies (GWAS) in Drosophila melanogaster to identify the genetic basis of variation in susceptibility to viral infection. We found that there is substantially more genetic variation in susceptibility to two viruses that naturally infect D. melanogaster (DCV and DMelSV) than to two viruses isolated from other insects (FHV and DAffSV). Furthermore, this increased variation is caused by a small number of common polymorphisms that have a major effect on resistance and can individually explain up to 47% of the heritability in disease susceptibility. For two of these polymorphisms, it has previously been shown that they have been driven to a high frequency by natural selection. An advantage of GWAS in Drosophila is that the results can be confirmed experimentally. We verified that a gene called pastrel--which was previously not known to have an antiviral function--is associated with DCV-resistance by knocking down its expression by RNAi. Our data suggest that selection for resistance to infectious disease can increase genetic variation by increasing the frequency of major-effect alleles, and this has resulted in a simple genetic basis to variation in virus resistance.