Immune genes undergo more adaptive evolution than non-immune system genes in Daphnia pulex

<p>Abstract</p> <p>Background</p> <p>Understanding which parts of the genome have been most influenced by adaptive evolution remains an unsolved puzzle. Some evidence suggests that selection has the greatest impact on regions of the genome that interact with other evolving genomes, including loci that are involved in host-parasite co-evolutionary processes. In this study, we used a population genetic approach to test this hypothesis by comparing DNA sequences of 30 putative immune system genes in the crustacean <it>Daphnia pulex</it> with 24 non-immune system genes.</p> <p>Results</p> <p>In support of the hypothesis, results from a multilocus extension of the McDonald-Kreitman (MK) test indicate that immune system genes as a class have experienced more adaptive evolution than non-immune system genes. However, not all immune system genes show evidence of adaptive evolution. Additionally, we apply single locus MK tests and calculate population genetic parameters at all loci in order to characterize the mode of selection (directional versus balancing) in the genes that show the greatest deviation from neutral evolution.</p> <p>Conclusions</p> <p>Our data are consistent with the hypothesis that immune system genes undergo more adaptive evolution than non-immune system genes, possibly as a result of host-parasite arms races. The results of these analyses highlight several candidate loci undergoing adaptive evolution that could be targeted in future studies.</p

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared to Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction

Novel insights into the insect trancriptome response to a natural DNA virus

ArticleCopyright © 2015 McTaggart et al.; licensee BioMed Central.Background Little is known about invertebrate responses to DNA viruses. Here, we infect a commercially important pest moth species Plodia interpunctella with its naturally infecting DNA virus. We sequenced, assembled and annotated the complete transcriptome of the moth, and a partial transcriptome of the virus. We then tested for differential gene expression between moths that were exposed to the virus and controls. Results We found 51 genes that were differentially expressed in moths exposed to a DNA baculovirus compared to controls. Gene set enrichment analysis revealed that cuticle proteins were significantly overrepresented in this group of genes. Interestingly, 6 of the 7 differentially expressed cuticle proteins were downregulated, suggesting that baculoviruses are able to manipulate its host’s response. In fact, an additional 29 of the 51 genes were also downregulated in exposed compared with control animals, including a gram-negative binding protein. In contrast, genes involved in transposable element movement were upregulated after infection. Conclusions We present the first experiment to measure genome-wide gene expression in an insect after infection with a natural DNA virus. Our results indicate that cuticle proteins might be key genes underpinning the response to DNA viruses. Furthermore, the large proportion of genes that were downregulated after viral exposure suggests that this virus is actively manipulating the insect immune response. Finally, it appears that transposable element activity might increase during viral invasion. Combined, these results provide much needed host candidate genes that respond to DNA viral invaders.NERC Biomolecular Analysis Facility (NBAF

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the “universal” genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3’UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids

The genome sequence of the Violet Carpenter Bee, Xylocopa violacea (Linnaeus, 1785): A hymenopteran species undergoing range expansion

We present a reference genome assembly from an individual male Violet Carpenter Bee (Xylocopa violacea, Linnaeus 1758). The assembly is 1.02 gigabases in span. 48% of the assembly is scaffolded into 17 pseudo-chromosomal units. The mitochondrial genome has also been assembled and is 21.8 kilobases in length. The genome is highly repetitive, likely representing a highly heterochromatic architecture expected of bees from the genus Xylocopa. We also use an evidence-based methodology to annotate 10,152 high confidence coding genes. This genome was sequenced as part of the pilot project of the European Reference Genome Atlas (ERGA) and represents an important addition to the genomic resources available for Hymenoptera

Transcriptome profiling during a natural host-parasite interaction

BACKGROUND: Infection outcome in some coevolving host-pathogens is characterised by host-pathogen genetic interactions, where particular host genotypes are susceptible only to a subset of pathogen genotypes. To identify candidate genes responsible for the infection status of the host, we exposed a Daphnia magna host genotype to two bacterial strains of Pasteuria ramosa, one of which results in infection, while the other does not. At three time points (four, eight and 12 h) post pathogen exposure, we sequenced the complete transcriptome of the hosts using RNA-Seq (Illumina). RESULTS: We observed a rapid and transient response to pathogen treatment. Specifically, at the four-hour time point, eight genes were differentially expressed. At the eight-hour time point, a single gene was differentially expressed in the resistant combination only, and no genes were differentially expressed at the 12-h time point. CONCLUSIONS: We found that pathogen-associated transcriptional activity is greatest soon after exposure. Genome-wide resistant combinations were more likely to show upregulation of genes, while susceptible combinations were more likely to be downregulated, relative to controls. Our results also provide several novel candidate genes that may play a pivotal role in determining infection outcomes. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-015-1838-0) contains supplementary material, which is available to authorized users

Rates of Recombination in the Ribosomal DNA of Apomictically Propagated Daphnia obtusa Lines

Ribosomal (r)DNA undergoes concerted evolution, the mechanisms of which are unequal crossing over and gene conversion. Despite the fundamental importance of these mechanisms to the evolution of rDNA, their rates have been estimated only in a few model species. We estimated recombination rate in rDNA by quantifying the relative frequency of intraindividual length variants in an expansion segment of the 18S rRNA gene of the cladoceran crustacean, Daphnia obtusa, in four apomictically propagated lines. We also used quantitative PCR to estimate rDNA copy number. The apomictic lines were sampled every 5 generations for 90 generations, and we considered each significant change in the frequency distribution of length variants between time intervals to be the result of a recombination event. Using this method, we calculated the recombination rate for this region to be 0.02–0.06 events/generation on the basis of three different estimates of rDNA copy number. In addition, we observed substantial changes in rDNA copy number within and between lines. Estimates of haploid copy number varied from 53 to 233, with a mean of 150. We also measured the relative frequency of length variants in 30 lines at generations 5, 50, and 90. Although length variant frequencies changed significantly within and between lines, the overall average frequency of each length variant did not change significantly between the three generations sampled, suggesting that there is little or no bias in the direction of change due to recombination

Nucleotide polymorphism and within-gene recombination in Daphnia magna and D. pulex, two cyclical parthenogens

Theory predicts that partially asexual organisms may make the "best of both worlds": for the most part, they avoid the costs of sexual reproduction, while still benefiting from an enhanced efficiency of selection compared to obligately asexual organisms. There is, however, little empirical data on partially asexual organisms to test this prediction. Here we examine patterns of nucleotide diversity at eight nuclear loci in continentwide samples of two species of cyclically parthenogenetic Daphnia to assess the effect of partial asexual reproduction on effective population size and amount of recombination. Both species have high nucleotide diversities and show abundant evidence for recombination, yielding large estimates of effective population sizes (300,000–600,000). This suggests that selection will act efficiently even on mutations with small selection coefficients. Divergence between the two species is less than one-tenth of previous estimates, which were derived using a mitochondrial molecular clock. As the two species investigated are among the most distantly related species of the genus, this suggests that the genus Daphnia may be considerably younger than previously thought. Daphnia has recently received increased attention because it is being developed as a model organism for ecological and evolutionary genomics. Our results confirm the attractiveness of Daphnia as a model organism, because the high nucleotide diversity and low linkage disequilibrium suggest that fine-scale mapping of genes affecting phenotypes through association studies should be feasible