Recent and recurrent selective sweeps of the antiviral RNAi gene Argonaute-2 in three species of Drosophila

Blessing or curse: how the epigenetic resolution of host-transposable element conflicts shapes their evolutionary dynamics

Transposable elements (TEs) are selfish genetic elements whose antagonistic interactions with hosts represent a common genetic conflict in eukaryotes. To resolve this conflict, hosts have widely adopted epigenetic silencing that deposits repressive marks at TEs. However, this mechanism is imperfect and fails to fully halt TE replication. Furthermore, TE epigenetic silencing can inadvertently spread repressive marks to adjacent functional sequences, a phenomenon considered a 'curse' of this conflict resolution. Here, we used forward simulations to explore how TE epigenetic silencing and its harmful side effects shape the evolutionary dynamics of TEs and their hosts. Our findings reveal that epigenetic silencing allows TEs and their hosts to stably coexist under a wide range of conditions, because the underlying molecular mechanisms give rise to copy-number dependency of the strength of TE silencing. Interestingly, contrary to intuitive expectations that TE epigenetic silencing should evolve to be as strong as possible, we found a selective benefit for modifier alleles that weaken TE silencing under biologically feasible conditions. These results reveal that the dual nature of TE epigenetic silencing, with both positive and negative effects, complicates its evolutionary trajectory and makes it challenging to determine whether TE epigenetic silencing is a 'blessing' or a 'curse'.

Bacterial origins of human cell-autonomous innate immune mechanisms

The cell-autonomous innate immune system enables animal cells to resist viral infection. This system comprises an array of sensors that, after detecting viral molecules, activate the expression of antiviral proteins and the interferon response. The repertoire of immune sensors and antiviral proteins has long been considered to be derived from extensive evolutionary innovation in vertebrates, but new data challenge this dogma. Recent studies show that central components of the cell-autonomous innate immune system have ancient evolutionary roots in prokaryotic genes that protect bacteria from phages. These include the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, Toll/IL-1 receptor (TIR) domain-containing pathogen receptors, the viperin family of antiviral proteins, SAMHD1-like nucleotide-depletion enzymes, gasdermin proteins and key components of the RNA interference pathway. This Perspective details current knowledge of the elements of antiviral immunity that are conserved from bacteria to humans, and presents possible evolutionary scenarios to explain the observed conservation.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

Insecticidal RNA interference, thinking beyond long dsRNA

Over 20 years ago double-stranded RNA (dsRNA) was described as the trigger of RNAi interference (RNAi)-based gene silencing. This paradigm has held since, especially for insect biopesticide technologies where dsRNAs, similar to those described in 1998, are used to inhibit gene expression. In the intervening years, investigation of RNAi pathways has revealed the small RNA effectors of RNAi are diverse and rapidly evolving. The rich biology of insect small RNAs suggests potential to use multiple RNAi modes for manipulating gene expression. By exploiting different RNAi pathways, the menu of options for pest control can be expanded and could lead to better tailored solutions. Fortunately, basic delivery strategies used for dsRNA such as direct application or transgenic expression will translate well between RNAs transiting different RNAi pathways. Importantly, further engineering of RNAi-based biopesticides may provide an opportunity to address dsRNA insensitivity seen in some pests. Characterization of RNAi pathways unique to target species will be indispensable to this end and may require thinking beyond long dsRNA. © 2020 Society of Chemical Industry.

Population Genomics on the Fly: Recent Advances in Drosophila

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

Telomeric TART elements target the piRNA machinery in Drosophila

Coevolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by reestablishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. The TART-A TE functions as an important component of Drosophila telomeres but has also reportedly inserted into the Drosophila melanogaster nuclear export factor gene nxf2. We find that, rather than inserting into nxf2, TART-A has actually captured a portion of nxf2 sequence. We show that TART-A produces abundant Piwi-interacting small RNAs (piRNAs), some of which are antisense to the nxf2 transcript, and that the TART-like region of nxf2 is evolving rapidly. Furthermore, in D. melanogaster, TART-A is present at higher copy numbers, and nxf2 shows reduced expression, compared to the closely related species Drosophila simulans. We propose that capturing nxf2 sequence allowed TART-A to target the nxf2 gene for piRNA-mediated repression and that these 2 elements are engaged in antagonistic coevolution despite the fact that TART-A is serving a critical role for its host genome.

Co-evolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by re-establishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. This study shows that a specialized Drosophila retrotransposon that functions as a telomere has captured a portion of a host piRNA gene which may allow it to evade silencing.

Inference of coevolutionary dynamics and parameters from host and parasite polymorphism data of repeated experiments

There is a long-standing interest in understanding host-parasite coevolutionary dynamics and associated fitness effects. Increasing amounts of genomic data for both interacting species offer a promising source to identify candidate loci and to infer the main parameters of the past coevolutionary history. However, so far no method exists to perform the latter. By coupling a gene-for-gene model with coalescent simulations, we first show that three types of biological costs, namely, resistance, infectivity and infection, define the allele frequencies at the internal equilibrium point of the coevolution model. These in return determine the strength of selective signatures at the coevolving host and parasite loci. We apply an Approximate Bayesian Computation (ABC) approach on simulated datasets to infer these costs by jointly integrating host and parasite polymorphism data at the coevolving loci. To control for the effect of genetic drift on coevolutionary dynamics, we assume that 10 or 30 repetitions are available from controlled experiments or several natural populations. We study two scenarios: 1) the cost of infection and population sizes (host and parasite) are unknown while costs of infectivity and resistance are known, and 2) all three costs are unknown while populations sizes are known. Using the ABC model choice procedure, we show that for both scenarios, we can distinguish with high accuracy pairs of coevolving host and parasite loci from pairs of neutrally evolving loci, though the statistical power decreases with higher cost of infection. The accuracy of parameter inference is high under both scenarios especially when using both host and parasite data because parasite polymorphism data do inform on costs applying to the host and vice-versa. As the false positive rate to detect pairs of genes under coevolution is small, we suggest that our method complements recently developed methods to identify host and parasite candidate loci for functional studies.

Arboviruses and the Challenge to Establish Systemic and Persistent Infections in Competent Mosquito Vectors: The Interaction With the RNAi Mechanism

Arboviruses are capable to establish long-term persistent infections in mosquitoes that do not affect significantly the physiology of the insect vectors. Arbovirus infections are controlled by the RNAi machinery via the production of viral siRNAs and the formation of RISC complexes targeting viral genomes and mRNAs. Engineered arboviruses that contain cellular gene sequences can therefore be transformed to "viral silencing vectors" for studies of gene function in reverse genetics approaches. More specifically, "ideal" viral silencing vectors must be competent to induce robust RNAi effects while other interactions with the host immune system should be kept at a minimum to reduce non-specific effects. Because of their inconspicuous nature, arboviruses may approach the "ideal" viral silencing vectors in insects and it is therefore worthwhile to study the mechanisms by which the interactions with the RNAi machinery occur. In this review, an analysis is presented of the antiviral RNAi response in mosquito vectors with respect to the major types of arboviruses (alphaviruses, flaviviruses, bunyaviruses, and others). With respect to antiviral defense, the exo-RNAi pathway constitutes the major mechanism while the contribution of both miRNAs and viral piRNAs remains a contentious issue. However, additional mechanisms exist in mosquitoes that are capable to enhance or restrict the efficiency of viral silencing vectors such as the amplification of RNAi effects by DNA forms, the existence of incorporated viral elements in the genome and the induction of a non-specific systemic response by Dicer-2. Of significance is the observation that no major "viral suppressors of RNAi" (VSRs) seem to be encoded by arboviral genomes, indicating that relatively tight control of the activity of the RNA-dependent RNA polymerase (RdRp) may be sufficient to maintain the persistent character of arbovirus infections. Major strategies for improvement of viral silencing vectors therefore are proposed to involve engineering of VSRs and modifying of the properties of the RdRp. Because of safety issues (pathogen status), however, arbovirus-based silencing vectors are not well suited for practical applications, such as RNAi-based mosquito control. In that case, related mosquito-specific viruses that also establish persistent infections and may cause similar RNAi responses may represent a valuable alternative solution.

Coevolution of hytrosaviruses and host immune responses

BACKGROUND

Hytrosaviruses (SGHVs; Hytrosaviridae family) are double-stranded DNA (dsDNA) viruses that cause salivary gland hypertrophy (SGH) syndrome in flies. Two structurally and functionally distinct SGHVs are recognized; Glossina pallidipes SGHV (GpSGHV) and Musca domestica SGHV (MdSGHV), that infect the hematophagous tsetse fly and the filth-feeding housefly, respectively. Genome sizes and gene contents of GpSGHV (~ 190 kb; 160-174 genes) and MdSGHV (~ 124 kb; 108 genes) may reflect an evolution with the SGHV-hosts resulting in differences in pathobiology. Whereas GpSGHV can switch from asymptomatic to symptomatic infections in response to certain unknown cues, MdSGHV solely infects symptomatically. Overt SGH characterizes the symptomatic infections of SGHVs, but whereas MdSGHV induces both nuclear and cellular hypertrophy (enlarged non-replicative cells), GpSGHV induces cellular hyperplasia (enlarged replicative cells). Compared to GpSGHV's specificity to Glossina species, MdSGHV infects other sympatric muscids. The MdSGHV-induced total shutdown of oogenesis inhibits its vertical transmission, while the GpSGHV's asymptomatic and symptomatic infections promote vertical and horizontal transmission, respectively. This paper reviews the coevolution of the SGHVs and their hosts (housefly and tsetse fly) based on phylogenetic relatedness of immune gene orthologs/paralogs and compares this with other virus-insect models.

RESULTS

Whereas MdSGHV is not vertically transmitted, GpSGHV is both vertically and horizontally transmitted, and the balance between the two transmission modes may significantly influence the pathogenesis of tsetse virus. The presence and absence of bacterial symbionts (Wigglesworthia and Sodalis) in tsetse and Wolbachia in the housefly, respectively, potentially contributes to the development of SGH symptoms. Unlike MdSGHV, GpSGHV contains not only host-derived proteins, but also appears to have evolutionarily recruited cellular genes from ancestral host(s) into its genome, which, although may be nonessential for viral replication, potentially contribute to the evasion of host's immune responses. Whereas MdSGHV has evolved strategies to counteract both the housefly's RNAi and apoptotic responses, the housefly has expanded its repertoire of immune effector, modulator and melanization genes compared to the tsetse fly.

CONCLUSIONS

The ecologies and life-histories of the housefly and tsetse fly may significantly influence coevolution of MdSGHV and GpSGHV with their hosts. Although there are still many unanswered questions regarding the pathogenesis of SGHVs, and the extent to which microbiota influence expression of overt SGH symptoms, SGHVs are attractive 'explorers' to elucidate the immune responses of their hosts, and the transmission modes of other large DNA viruses.

Interaction of entomopathogenic fungi with the host immune system

Entomopathogenic fungi can invade wide range of insect hosts in the natural world and have been used as environmentally friendly alternatives to chemical insecticides for pest control. Studies of host-pathogen interactions provide valuable insights into the coevolutionay arms race between fungal pathogens and their hosts. Entomopathogenic fungi have evolved a series of sophisticated strategies to counter insect immune defenses. In response to fungal infection, insect hosts rely on behavior avoidance, physical barrier and innate immune defenses in the fight against invading pathogens. The insect cuticle acts as the first physical barrier against pathogens. It is an inhospitable physiological environment that contains chemicals (e.g., antimicrobial peptides and reactive oxygen species), which inhibit fungal growth. In addition, innate immune responses, including cellular immunity and humoral immunity, play critical roles in preventing fungal infection. In this review, we outline the current state of our knowledge of insect defenses to fungal infection and discuss the strategies by which entomopathogenic fungi counter the host immune system. Increased knowledge regarding the molecular interactions between entomopathogenic fungi and the insect host could provide new strategies for pest management.

Natural Variation in Resistance to Virus Infection in Dipteran Insects

The power and ease of Drosophila genetics and the medical relevance of mosquito-transmitted viruses have made dipterans important model organisms in antiviral immunology. Studies of virus-host interactions at the molecular and population levels have illuminated determinants of resistance to virus infection. Here, we review the sources and nature of variation in antiviral immunity and virus susceptibility in model dipteran insects, specifically the fruit fly Drosophila melanogaster and vector mosquitoes of the genera Aedes and Culex. We first discuss antiviral immune mechanisms and describe the virus-specificity of these responses. In the following sections, we review genetic and microbiota-dependent variation in antiviral immunity. In the final sections, we explore less well-studied sources of variation, including abiotic factors, sexual dimorphism, infection history, and endogenous viral elements. We borrow from work on other pathogen types and non-dipteran species when it parallels or complements studies in dipterans. Understanding natural variation in virus-host interactions may lead to the identification of novel restriction factors and immune mechanisms and shed light on the molecular determinants of vector competence.

RNA-Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates

Conflict between organisms can lead to a reciprocal adaptation that manifests as an increased evolutionary rate in genes mediating the conflict. This adaptive signature has been observed in RNA-interference (RNAi) pathway genes involved in the suppression of viruses and transposable elements in Drosophila melanogaster, suggesting that a subset of Drosophila RNAi genes may be locked in an arms race with these parasites. However, it is not known whether rapid evolution of RNAi genes is a general phenomenon across invertebrates, or which RNAi genes generally evolve adaptively. Here we use population genomic data from eight invertebrate species to infer rates of adaptive sequence evolution, and to test for past and ongoing selective sweeps in RNAi genes. We assess rates of adaptive protein evolution across species using a formal meta-analytic framework to combine data across species and by implementing a multispecies generalized linear mixed model of mutation counts. Across species, we find that RNAi genes display a greater rate of adaptive protein substitution than other genes, and that this is primarily mediated by positive selection acting on the genes most likely to defend against viruses and transposable elements. In contrast, evidence for recent selective sweeps is broadly spread across functional classes of RNAi genes and differs substantially among species. Finally, we identify genes that exhibit elevated adaptive evolution across the analyzed insect species, perhaps due to concurrent parasite-mediated arms races.

Simulated climate change, epidemic size, and host evolution across host-parasite populations

Climate change is causing warmer and more variable temperatures as well as physical flux in natural populations, which will affect the ecology and evolution of infectious disease epidemics. Using replicate seminatural populations of a coevolving freshwater invertebrate-parasite system (host: Daphnia magna, parasite: Pasteuria ramosa), we quantified the effects of ambient temperature and population mixing (physical flux within populations) on epidemic size and population health. Each population was seeded with an identical suite of host genotypes and dose of parasite transmission spores. Biologically reasonable increases in environmental temperature caused larger epidemics, and population mixing reduced overall epidemic size. Mixing also had a detrimental effect on host populations independent of disease. Epidemics drove parasite-mediated selection, leading to a loss of host genetic diversity, and mixed populations experienced greater evolution due to genetic drift over the season. These findings further our understanding of how diversity loss will reduce the host populations' capacity to respond to changes in selection, therefore stymying adaptation to further environmental change.

Environmental variation causes different (co) evolutionary routes to the same adaptive destination across parasite populations

Epidemics are engines for host-parasite coevolution, where parasite adaptation to hosts drives reciprocal adaptation in host populations. A key challenge is to understand whether parasite adaptation and any underlying evolution and coevolution is repeatable across ecologically realistic populations that experience different environmental conditions, or if each population follows a completely unique evolutionary path. We established twenty replicate pond populations comprising an identical suite of genotypes of crustacean host, Daphnia magna, and inoculum of their parasite, Pasteuria ramosa. Using a time-shift experiment, we compared parasite infection traits before and after epidemics and linked patterns of parasite evolution with shifts in host genotype frequencies. Parasite adaptation to the sympatric suite of host genotypes came at a cost of poorer performance on foreign genotypes across populations and environments. However, this consistent pattern of parasite adaptation was driven by different types of frequency-dependent selection that was contingent on an ecologically relevant environmental treatment (whether or not there was physical mixing of water within ponds). In unmixed ponds, large epidemics drove rapid and strong host-parasite coevolution. In mixed ponds, epidemics were smaller and host evolution was driven mainly by the mixing treatment itself; here, host evolution and parasite evolution were clear, but coevolution was absent. Population mixing breaks an otherwise robust coevolutionary cycle. These findings advance our understanding of the repeatability of (co)evolution across noisy, ecologically realistic populations.

Genomic signatures of local adaptation in the Drosophila immune response

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

Repeated Duplication of Argonaute2 Is Associated with Strong Selection and Testis Specialization in Drosophila

Argonaute2 (Ago2) is a rapidly evolving nuclease in the Drosophila melanogaster RNA interference (RNAi) pathway that targets viruses and transposable elements in somatic tissues. Here we reconstruct the history of Ago2 duplications across the D. obscura group and use patterns of gene expression to infer new functional specialization. We show that some duplications are old, shared by the entire species group, and that losses may be common, including previously undetected losses in the lineage leading to D. pseudoobscura We find that while the original (syntenic) gene copy has generally retained the ancestral ubiquitous expression pattern, most of the novel Ago2 paralogs have independently specialized to testis-specific expression. Using population genetic analyses, we show that most testis-specific paralogs have significantly lower genetic diversity than the genome-wide average. This suggests recent positive selection in three different species, and model-based analyses provide strong evidence of recent hard selective sweeps in or near four of the six D. pseudoobscura Ago2 paralogs. We speculate that the repeated evolution of testis specificity in obscura group Ago2 genes, combined with their dynamic turnover and strong signatures of adaptive evolution, may be associated with highly derived roles in the suppression of transposable elements or meiotic drive. Our study highlights the lability of RNAi pathways, even within well-studied groups such as Drosophila, and suggests that strong selection may act quickly after duplication in RNAi pathways, potentially giving rise to new and unknown RNAi functions in nonmodel species.

Variation and Evolution in the Glutamine-Rich Repeat Region of Drosophila Argonaute-2

RNA interference pathways mediate biological processes through Argonaute-family proteins, which bind small RNAs as guides to silence complementary target nucleic acids . In insects and crustaceans Argonaute-2 silences viral nucleic acids, and therefore acts as a primary effector of innate antiviral immunity. Although the function of the major Argonaute-2 domains, which are conserved across most Argonaute-family proteins, are known, many invertebrate Argonaute-2 homologs contain a glutamine-rich repeat (GRR) region of unknown function at the N-terminus . Here we combine long-read amplicon sequencing of Drosophila Genetic Reference Panel (DGRP) lines with publicly available sequence data from many insect species to show that this region evolves extremely rapidly and is hyper-variable within species. We identify distinct GRR haplotype groups in Drosophila melanogaster, and suggest that one of these haplotype groups has recently risen to high frequency in a North American population. Finally, we use published data from genome-wide association studies of viral resistance in D. melanogaster to test whether GRR haplotypes are associated with survival after virus challenge. We find a marginally significant association with survival after challenge with Drosophila C Virus in the DGRP, but we were unable to replicate this finding using lines from the Drosophila Synthetic Population Resource panel.

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.

The diversity of insect antiviral immunity: insights from viruses

Insects represent over 70% of all animal species. Recent virome analyses reveal unprecedented genetic diversity of insect viruses, which appears to match that of their hosts. Thus, insect-virus interactions may provide information on a vast repertoire of antiviral immune mechanisms. Tapping into this diversity is challenging because of several constraints imposed by the uniqueness of each insect model. Nevertheless, it is clear that many conserved and divergent pathways participate in the control of viral infection in insects. Co-evolution between hosts and viruses favors the development of immune evasion mechanisms by the pathogen. Viral suppressors can offer unique perspective on host pathways and emphasize the importance of RNA interference, apoptosis, but also NF-κB pathways and translation control in insect antiviral immunity.

Early origin and adaptive evolution of the GW182 protein family, the key component of RNA silencing in animals

The GW182 proteins are a key component of the miRNA-dependent post-transcriptional silencing pathway in animals. They function as scaffold proteins to mediate the interaction of Argonaute (AGO)-containing complexes with cytoplasmic poly(A)-binding proteins (PABP) and PAN2-PAN3 and CCR4-NOT deadenylases. The AGO-GW182 complexes mediate silencing of the target mRNA through induction of translational repression and/or mRNA degradation. Although the GW182 proteins are a subject of extensive experimental research in the recent years, very little is known about their origin and evolution. Here, based on complex functional annotation and phylogenetic analyses, we reveal 448 members of the GW182 protein family from the earliest animals to humans. Our results indicate that a single-copy GW182/TNRC6C progenitor gene arose with the emergence of multicellularity and it multiplied in the last common ancestor of vertebrates in 2 rounds of whole genome duplication (WGD) resulting in 3 genes. Before the divergence of vertebrates, both the AGO- and CCR4-NOT-binding regions of GW182s showed significant acceleration in the accumulation of amino acid changes, suggesting functional adaptation toward higher specificity to the molecules of the silencing complex. We conclude that the silencing ability of the GW182 proteins improves with higher position in the taxonomic classification and increasing complexity of the organism. The first reconstruction of the molecular journey of GW182 proteins from the ancestral metazoan protein to the current mammalian configuration provides new insight into development of the miRNA-dependent post-transcriptional silencing pathway in animals.

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.

Recurrent Gene Duplication Diversifies Genome Defense Repertoire in Drosophila

Transposable elements (TEs) comprise large fractions of many eukaryotic genomes and imperil host genome integrity. The host genome combats these challenges by encoding proteins that silence TE activity. Both the introduction of new TEs via horizontal transfer and TE sequence evolution requires constant innovation of host-encoded TE silencing machinery to keep pace with TEs. One form of host innovation is the adaptation of existing, single-copy host genes. Indeed, host suppressors of TE replication often harbor signatures of positive selection. Such signatures are especially evident in genes encoding the piwi-interacting-RNA pathway of gene silencing, for example, the female germline-restricted TE silencer, HP1D/Rhino Host genomes can also innovate via gene duplication and divergence. However, the importance of gene family expansions, contractions, and gene turnover to host genome defense has been largely unexplored. Here, we functionally characterize Oxpecker, a young, tandem duplicate gene of HP1D/rhino We demonstrate that Oxpecker supports female fertility in Drosophila melanogaster and silences several TE families that are incompletely silenced by HP1D/Rhino in the female germline. We further show that, like Oxpecker, at least ten additional, structurally diverse, HP1D/rhino-derived daughter and "granddaughter" genes emerged during a short 15-million year period of Drosophila evolution. These young paralogs are transcribed primarily in germline tissues, where the genetic conflict between host genomes and TEs plays out. Our findings suggest that gene family expansion is an underappreciated yet potent evolutionary mechanism of genome defense diversification.

The Discovery, Distribution, and Evolution of Viruses Associated with Drosophila melanogaster

Drosophila melanogaster is a valuable invertebrate model for viral infection and antiviral immunity, and is a focus for studies of insect-virus coevolution. Here we use a metagenomic approach to identify more than 20 previously undetected RNA viruses and a DNA virus associated with wild D. melanogaster. These viruses not only include distant relatives of known insect pathogens but also novel groups of insect-infecting viruses. By sequencing virus-derived small RNAs, we show that the viruses represent active infections of Drosophila. We find that the RNA viruses differ in the number and properties of their small RNAs, and we detect both siRNAs and a novel miRNA from the DNA virus. Analysis of small RNAs also allows us to identify putative viral sequences that lack detectable sequence similarity to known viruses. By surveying >2,000 individually collected wild adult Drosophila we show that more than 30% of D. melanogaster carry a detectable virus, and more than 6% carry multiple viruses. However, despite a high prevalence of the Wolbachia endosymbiont--which is known to be protective against virus infections in Drosophila--we were unable to detect any relationship between the presence of Wolbachia and the presence of any virus. Using publicly available RNA-seq datasets, we show that the community of viruses in Drosophila laboratories is very different from that seen in the wild, but that some of the newly discovered viruses are nevertheless widespread in laboratory lines and are ubiquitous in cell culture. By sequencing viruses from individual wild-collected flies we show that some viruses are shared between D. melanogaster and D. simulans. Our results provide an essential evolutionary and ecological context for host-virus interaction in Drosophila, and the newly reported viral sequences will help develop D. melanogaster further as a model for molecular and evolutionary virus research.

Evolutionary genetics of insect innate immunity

Patterns of evolution in immune defense genes help to understand the evolutionary dynamics between hosts and pathogens. Multiple insect genomes have been sequenced, with many of them having annotated immune genes, which paves the way for a comparative genomic analysis of insect immunity. In this review, I summarize the current state of comparative and evolutionary genomics of insect innate immune defense. The focus is on the conserved and divergent components of immunity with an emphasis on gene family evolution and evolution at the sequence level; both population genetics and molecular evolution frameworks are considered.

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host-parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host-parasite interactions. In the QTL panel used here, Daphnia's resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host-parasite systems. Only the Pasteuria-Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium-Daphnia system remains unclear.

The genetics of host-virus coevolution in invertebrates

Although viral infection and antiviral defence are ubiquitous, genetic data are currently unavailable from the vast majority of animal phyla-potentially biasing our overall perspective of the coevolutionary process. Rapid adaptive evolution is seen in some insect antiviral genes, consistent with invertebrate-virus 'arms-race' coevolution, but equivalent signatures of selection are hard to detect in viruses. We find that, despite the large differences in vertebrate, invertebrate, and plant immune responses, comparison of viral evolution fails to identify any difference among these hosts in the impact of positive selection. The best evidence for invertebrate-virus coevolution is currently provided by large-effect polymorphisms for host resistance and/or viral evasion, as these often appear to have arisen and spread recently, and can be favoured by virus-mediated selection.

Recent insights into the evolution of innate viral sensing in animals

•

Viral sensor evolution may be constrained by highly conserved viral immune-elicitors.

•

Alternatively, viral sensor evolution may be driven by ‘arms race’ coevolution.

•

We find viral-sensing Toll-like receptors evolve more slowly than other TLRs.

•

In contrast, viral-sensing helicase related genes often evolve rapidly and adaptively.

The evolution of viral sensors is likely to be shaped by the constraint imposed through high conservation of viral Pathogen-Associated Molecular Patterns (PAMPs), and by the potential for ‘arms race’ coevolution with more rapidly evolving viral proteins. Here we review the recent progress made in understanding the evolutionary history of two types of viral sensor, RNA helicases and Toll-like receptors. We find differences both in their rates of evolution, and in the levels of positive selection they experience. We suggest that positive selection has been the primary driver of the rapid evolution of the RNA helicases, while selective constraint has been a stronger influence shaping the slow evolution of the Toll-like receptors.

Novel Drosophila viruses encode host-specific suppressors of RNAi

The ongoing conflict between viruses and their hosts can drive the co-evolution between host immune genes and viral suppressors of immunity. It has been suggested that an evolutionary ‘arms race’ may occur between rapidly evolving components of the antiviral RNAi pathway of Drosophila and viral genes that antagonize it. We have recently shown that viral protein 1 (VP1) of Drosophila melanogaster Nora virus (DmelNV) suppresses Argonaute-2 (AGO2)-mediated target RNA cleavage (slicer activity) to antagonize antiviral RNAi. Here we show that viral AGO2 antagonists of divergent Nora-like viruses can have host specific activities. We have identified novel Nora-like viruses in wild-caught populations of D. immigrans (DimmNV) and D. subobscura (DsubNV) that are 36% and 26% divergent from DmelNV at the amino acid level. We show that DimmNV and DsubNV VP1 are unable to suppress RNAi in D. melanogaster S2 cells, whereas DmelNV VP1 potently suppresses RNAi in this host species. Moreover, we show that the RNAi suppressor activity of DimmNV VP1 is restricted to its natural host species, D. immigrans. Specifically, we find that DimmNV VP1 interacts with D. immigrans AGO2, but not with D. melanogaster AGO2, and that it suppresses slicer activity in embryo lysates from D. immigrans, but not in lysates from D. melanogaster. This species-specific interaction is reflected in the ability of DimmNV VP1 to enhance RNA production by a recombinant Sindbis virus in a host-specific manner. Our results emphasize the importance of analyzing viral RNAi suppressor activity in the relevant host species. We suggest that rapid co-evolution between RNA viruses and their hosts may result in host species-specific activities of RNAi suppressor proteins, and therefore that viral RNAi suppressors could be host-specificity factors.

Viruses and their hosts can engage in an evolutionary arms race. Viruses may select for hosts with more effective immune responses, whereas the immune response of the host may select for viruses that evade the immune system. These viral counter-defenses may in turn drive adaptations in host immune genes. A potential outcome of this perpetual cycle is that the interaction between virus and host becomes more specific. In insects, the host antiviral RNAi machinery exerts strong evolutionary pressure that has led to the evolution of viral proteins that can antagonize the RNAi response. We have identified novel viruses that infect different fruit fly species and we show that the RNAi suppressor proteins of these viruses can be specific to their host. Furthermore, we show that these proteins can enhance virus replication in a host-specific manner. These results are in line with the hypothesis that virus-host co-evolution shapes the genomes of both virus and host. Moreover, our results suggest that RNAi suppressor proteins have the potential to determine host specificity of viruses.

A recent global selective sweep on the age-1 phosphatidylinositol 3-OH kinase regulator of the insulin-like signaling pathway within Caenorhabditis remanei

The discovery that genetic pathways can be manipulated to extend lifespan has revolutionized our understanding of aging, yet their function within natural populations remains poorly characterized. In particular, evolutionary theories of aging predict tradeoffs in resource investment toward somatic maintenance vs. reproductive output that should impose strong natural selection on genetic components that influence this balance. To explore such selective pressure at the molecular level, we examine population genetic variation in the insulin-like signaling pathway of the nematode Caenorhabditis remanei. We document a recent global selective sweep on the phosphoinositide-3-kinase pathway regulator, age-1, the first life-extension gene to have been identified. In particular, we find that age-1 has 5−20 times less genetic variation than any other insulin-like signaling pathway components and that evolutionary signatures of selection center on the age-1 locus within its genomic environment. These results demonstrate that critical components of aging-related pathways can be subject to shifting patterns of strong selection, as predicted by theory. This highly polymorphic outcrossing species offers high-resolution, population-level analyses of molecular variation as a complement to functional genetic studies within the self-reproducing C. elegans model system.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

Suppressors of RNAi from plant viruses are subject to episodic positive selection

Viral suppressors of RNAi (VSRs) are proteins that actively inhibit the antiviral RNA interference (RNAi) immune response, providing an immune evasion route for viruses. It has been hypothesized that VSRs are engaged in a molecular ‘arms race’ with RNAi pathway genes. Two lines of evidence support this. First, VSRs from plant viruses display high sequence diversity, and are frequently gained and lost over evolutionary time scales. Second, Drosophila antiviral RNAi genes show high rates of adaptive evolution. Here, we investigate whether VSRs diversify faster than other genes and, if so, whether this is a result of positive selection, as might be expected in an arms race. By analysis of 12 plant RNA viruses, we show that the relative rate of protein evolution is higher for VSRs than for other genes, but that this is not attributable to pervasive positive selection. We argue that, because evolutionary time scales are extremely different for viruses and eukaryotes, it is improbable that viral adaptation (as measured by the ratio of non-synonymous to synonymous change) will be dominated by one-to-one coevolution with eukaryotes. Instead, for plant virus VSRs, we find strong evidence of episodic selection—diversifying selection that acts on a subset of lineages—which might be attributable to frequent shifts between different host genotypes or species.

Detection of selective sweeps in cattle using genome-wide SNP data

Background

The domestication and subsequent selection by humans to create breeds and biological types of cattle undoubtedly altered the patterning of variation within their genomes. Strong selection to fix advantageous large-effect mutations underlying domesticability, breed characteristics or productivity created selective sweeps in which variation was lost in the chromosomal region flanking the selected allele. Selective sweeps have now been identified in the genomes of many animal species including humans, dogs, horses, and chickens. Here, we attempt to identify and characterise regions of the bovine genome that have been subjected to selective sweeps.

Results

Two datasets were used for the discovery and validation of selective sweeps via the fixation of alleles at a series of contiguous SNP loci. BovineSNP50 data were used to identify 28 putative sweep regions among 14 diverse cattle breeds. Affymetrix BOS 1 prescreening assay data for five breeds were used to identify 85 regions and validate 5 regions identified using the BovineSNP50 data. Many genes are located within these regions and the lack of sequence data for the analysed breeds precludes the nomination of selected genes or variants and limits the prediction of the selected phenotypes. However, phenotypes that we predict to have historically been under strong selection include horned-polled, coat colour, stature, ear morphology, and behaviour.

Conclusions

The bias towards common SNPs in the design of the BovineSNP50 assay led to the identification of recent selective sweeps associated with breed formation and common to only a small number of breeds rather than ancient events associated with domestication which could potentially be common to all European taurines. The limited SNP density, or marker resolution, of the BovineSNP50 assay significantly impacted the rate of false discovery of selective sweeps, however, we found sweeps in common between breeds which were confirmed using an ultra-high-density assay scored in a small number of animals from a subset of the breeds. No sweep regions were shared between indicine and taurine breeds reflecting their divergent selection histories and the very different environmental habitats to which these sub-species have adapted.

Polymorphisms at the innate immune receptor TLR2 are associated with Borrelia infection in a wild rodent population

The discovery of the key role of Toll-like receptors (TLRs) in initiating innate immune responses and modulating adaptive immunity has revolutionized our understanding of vertebrate defence against pathogens. Yet, despite their central role in pathogen recognition and defence initiation, there is little information on how variation in TLRs influences disease susceptibility in natural populations. Here, we assessed the extent of naturally occurring polymorphisms at TLR2 in wild bank voles (Myodes glareolus) and tested for associations between TLR2 variants and infection with Borrelia afzelii, a common tick-transmitted pathogen in rodents and one of the causative agents of human Lyme disease. Bank voles in our population had 15 different TLR2 haplotypes (10 different haplotypes at the amino acid level), which grouped in three well-separated clusters. In a large-scale capture-mark-recapture study, we show that voles carrying TLR2 haplotypes of one particular cluster (TLR2c2) were almost three times less likely to be Borrelia infected than animals carrying other haplotypes. Moreover, neutrality tests suggested that TLR2 has been under positive selection. This is, to our knowledge, the first demonstration of an association between TLR polymorphism and parasitism in wildlife, and a striking example that genetic variation at innate immune receptors can have a large impact on host resistance.

Recurrent and recent selective sweeps in the piRNA pathway

Uncontrolled transposable element (TE) insertions and excisions can cause chromosome breaks and mutations with dramatic deleterious effects. The PIWI interacting RNA (piRNA) pathway functions as an adaptive TE silencing system during germline development. Several essential piRNA pathway proteins appear to be rapidly evolving, suggesting that TEs and the silencing machinery may be engaged in a classical "evolutionary arms race." Using a variety of molecular evolutionary and population genetic approaches, we find that the piRNA pathway genes rhino, krimper, and aubergine show patterns suggestive of extensive recurrent positive selection across Drosophila species. We speculate that selection on these proteins reflects crucial roles in silencing unfamiliar elements during vertical and horizontal transmission of TEs into naïve populations and species, respectively.

Genetic adaptation associated with genome-doubling in autotetraploid Arabidopsis arenosa

Genome duplication, which results in polyploidy, is disruptive to fundamental biological processes. Genome duplications occur spontaneously in a range of taxa and problems such as sterility, aneuploidy, and gene expression aberrations are common in newly formed polyploids. In mammals, genome duplication is associated with cancer and spontaneous abortion of embryos. Nevertheless, stable polyploid species occur in both plants and animals. Understanding how natural selection enabled these species to overcome early challenges can provide important insights into the mechanisms by which core cellular functions can adapt to perturbations of the genomic environment. Arabidopsis arenosa includes stable tetraploid populations and is related to well-characterized diploids A. lyrata and A. thaliana. It thus provides a rare opportunity to leverage genomic tools to investigate the genetic basis of polyploid stabilization. We sequenced the genomes of twelve A. arenosa individuals and found signatures suggestive of recent and ongoing selective sweeps throughout the genome. Many of these are at genes implicated in genome maintenance functions, including chromosome cohesion and segregation, DNA repair, homologous recombination, transcriptional regulation, and chromatin structure. Numerous encoded proteins are predicted to interact with one another. For a critical meiosis gene, ASYNAPSIS1, we identified a non-synonymous mutation that is highly differentiated by cytotype, but present as a rare variant in diploid A. arenosa, indicating selection may have acted on standing variation already present in the diploid. Several genes we identified that are implicated in sister chromatid cohesion and segregation are homologous to genes identified in a yeast mutant screen as necessary for survival of polyploid cells, and also implicated in genome instability in human diseases including cancer. This points to commonalities across kingdoms and supports the hypothesis that selection has acted on genes controlling genome integrity in A. arenosa as an adaptive response to genome doubling.

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.

Rapid intraspecific evolution of miRNA and siRNA genes in the mosquito Aedes aegypti

RNA silencing, or RNA interference (RNAi) in metazoans mediates development, reduces viral infection and limits transposon mobility. RNA silencing involves 21-30 nucleotide RNAs classified into microRNA (miRNA), exogenous and endogenous small interfering RNAs (siRNA), and Piwi-interacting RNA (piRNA). Knock-out, silencing and mutagenesis of genes in the exogenous siRNA (exo-siRNA) regulatory network demonstrate the importance of this RNAi pathway in antiviral immunity in Drosophila and mosquitoes. In Drosophila, genes encoding components for processing exo-siRNAs are among the fastest evolving 3% of all genes, suggesting that infection with pathogenic RNA viruses may drive diversifying selection in their host. In contrast, paralogous miRNA pathway genes do not evolve more rapidly than the genome average. Silencing of exo-siRNA pathway genes in mosquitoes orally infected with arboviruses leads to increased viral replication, but little is known about the comparative patterns of molecular evolution among the exo-siRNA and miRNA pathways genes in mosquitoes. We generated nearly complete sequences of all exons of major miRNA and siRNA pathway genes dicer-1 and dicer-2, argonaute-1 and argonaute-2, and r3d1 and r2d2 in 104 Aedes aegypti mosquitoes collected from six distinct geographic populations and analyzed their genetic diversity. The ratio of replacement to silent amino acid substitutions was 1.4 fold higher in dicer-2 than in dicer-1, 27.4 fold higher in argonaute-2 than in argonaute-1 and similar in r2d2 and r3d1. Positive selection was supported in 32% of non-synonymous sites in dicer-1, in 47% of sites in dicer-2, in 30% of sites in argonaute-1, in all sites in argonaute-2, in 22% of sites in r3d1 and in 55% of sites in r2d2. Unlike Drosophila, in Ae. aegypti, both exo-siRNA and miRNA pathway genes appear to be undergoing rapid, positive, diversifying selection. Furthermore, refractoriness of mosquitoes to infection with dengue virus was significantly positively correlated for nucleotide diversity indices in dicer-2.

Long-term and short-term evolutionary impacts of transposable elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host-TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host-pathogen coevolution.

Rapid genetic change underpins antagonistic coevolution in a natural host-pathogen metapopulation

Antagonistic coevolution is a critical force driving the evolution of diversity, yet the selective processes underpinning reciprocal adaptive changes in nature are not well understood. Local adaptation studies demonstrate partner impacts on fitness and adaptive change, but do not directly expose genetic processes predicted by theory. Specifically, we have little knowledge of the relative importance of fluctuating selection vs. arms-race dynamics in maintaining polymorphism in natural systems where metapopulation processes predominate. We conducted cross-year epidemiological, infection and genetic studies of multiple wild host and pathogen populations in the Linum-Melampsora association. We observed asynchronous phenotypic fluctuations in resistance and infectivity among demes. Importantly, changes in allelic frequencies at pathogen infectivity loci, and in host recognition of these genetic variants, correlated with disease prevalence during natural epidemics. These data strongly support reciprocal coevolution maintaining balanced resistance and infectivity polymorphisms, and highlight the importance of characterising spatial and temporal dynamics in antagonistic interactions.

Co-evolution between transposable elements and their hosts: a major factor in genome size evolution?

Most models of genome size evolution emphasize changes in relative rates of and/or the efficacy of selection on insertions and deletions. However, transposable elements (TEs) are a major contributor to genome size evolution, and since they experience their own selective pressures for expansion, genome size changes may in part be driven by the dynamics of co-evolution between TEs and their hosts. Under this perspective, predictions about the conditions that allow for genome expansion may be altered. In this review, we outline the evidence for TE-host co-evolution, discuss the conditions under which these dynamics can change, and explore the possible contribution to the evolution of genome size. Aided partly by advances in our understanding of the mechanisms of TE silencing via small RNAs, there is growing evidence that the evolution of transposition rates can be important in driving genome expansion and contraction. Shifts in genome size and transposon abundance associated with interspecific hybridization and changes in mating system are consistent with an important role for transposition rate evolution, although other possible explanations persist. More understanding of the potential for the breakdown of host silencing mechanisms and/or the potential for TEs to evade host immune responses will improve our understanding of the importance of changes in TE activity in driving genome size evolution.

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.

Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance

Visceral leishmaniasis is a potentially fatal disease endemic to large parts of Asia and Africa, primarily caused by the protozoan parasite Leishmania donovani. Here, we report a high-quality reference genome sequence for a strain of L. donovani from Nepal, and use this sequence to study variation in a set of 16 related clinical lines, isolated from visceral leishmaniasis patients from the same region, which also differ in their response to in vitro drug susceptibility. We show that whole-genome sequence data reveals genetic structure within these lines not shown by multilocus typing, and suggests that drug resistance has emerged multiple times in this closely related set of lines. Sequence comparisons with other Leishmania species and analysis of single-nucleotide diversity within our sample showed evidence of selection acting in a range of surface- and transport-related genes, including genes associated with drug resistance. Against a background of relative genetic homogeneity, we found extensive variation in chromosome copy number between our lines. Other forms of structural variation were significantly associated with drug resistance, notably including gene dosage and the copy number of an experimentally verified circular episome present in all lines and described here for the first time. This study provides a basis for more powerful molecular profiling of visceral leishmaniasis, providing additional power to track the drug resistance and epidemiology of an important human pathogen.

Host phylogeny determines viral persistence and replication in novel hosts

Pathogens switching to new hosts can result in the emergence of new infectious diseases, and determining which species are likely to be sources of such host shifts is essential to understanding disease threats to both humans and wildlife. However, the factors that determine whether a pathogen can infect a novel host are poorly understood. We have examined the ability of three host-specific RNA-viruses (Drosophila sigma viruses from the family Rhabdoviridae) to persist and replicate in 51 different species of Drosophilidae. Using a novel analytical approach we found that the host phylogeny could explain most of the variation in viral replication and persistence between different host species. This effect is partly driven by viruses reaching a higher titre in those novel hosts most closely related to the original host. However, there is also a strong effect of host phylogeny that is independent of the distance from the original host, with viral titres being similar in groups of related hosts. Most of this effect could be explained by variation in general susceptibility to all three sigma viruses, as there is a strong phylogenetic correlation in the titres of the three viruses. These results suggest that the source of new emerging diseases may often be predictable from the host phylogeny, but that the effect may be more complex than simply causing most host shifts to occur between closely related hosts.

Emerging infectious diseases such as SARS, HIV and swine-origin influenza have all been recently acquired by humans from other species. Understanding the reasons why parasites jump between different host species is essential to allow us to predict future threats and understand the causes of disease emergence. Here we ask how host-relatedness might determine when host-shifts can occur in the most important group of emerging diseases—RNA viruses. We show that the relationship between host species is the primary factor in determining a virus's ability to persist and replicate in a novel host following exposure. This can be broken down into two components. Firstly, species closely related to the virus's natural host are more susceptible than distantly related species. Secondly, independent of the distance effect, groups of closely related host species have similar levels of susceptibility. This has important implications for our understanding of disease-emergence, and until now the only large-scale studies of viruses have been correlative rather than experimental. We also found groups of related species that are susceptible to these viruses but are distantly related to the natural hosts, which may explain why viruses sometimes jump between distantly related species.