Daphnia magna egg piRNA cluster expression profiles change as mothers age

Background: PiRNAs prevent transposable elements wreaking havoc on the germline genome. Changes in piRNA expression over the lifetime of an individual may impact on ageing through continued suppression, or release, of transposable element expression. We identified piRNA producing clusters in the genome of Daphnia magna by a combination of bioinformatic methods, and then contrasted their expression between parthenogenetically produced eggs representing maternally-deposited germline piRNAs of young (having their 1st clutch) and old (having their 5th clutch) mothers. Results from eggs were compared to cluster expression in three generations of adults. Results: As for other arthropods, D. magna encodes long uni-directionally transcribed non-coding RNAs consisting of fragmented transposable elements which account for most piRNAs expressed. Egg tissues showed extensive differences between clutches from young mothers and those from old mothers, with 578 and 686 piRNA clusters upregulated, respectively. Most log fold-change differences for significant clusters were modest, however. When considering only highly expressed clusters, there was a bias towards 1st clutch eggs at 41 upregulated versus eight clusters in the eggs from older mothers. F0 generation differences between young and old mothers were fewer than eggs, as 179 clusters were up-regulated in young versus 170 old mothers. This dropped to 31 versus 22 piRNA clusters when comparing adults in the F1 generation, and no differences were detected in the F3 generation. Inter-generational losses of differential piRNA cluster were similar to that observed for D. magna micro-RNA expression. Conclusions: Little overlap in differentially expressed clusters was found between adults containing mixed somatic and germline (ovary) tissues and germ-line representing eggs. A cluster encompassing a Tudor domain containing gene important in the piRNA pathway was upregulated in the eggs from old mothers. We hypothesise that regulation of this gene could form part of a feedback loop that reduces piRNA pathway activity explaining the reduced number of highly-expressed clusters in eggs from old mothers

Exploring population history and gall induction in cynipid gall wasps using genomics and transcriptomics

Cynipid gall wasps have fascinating biology that has piqued the interest of naturalists throughout history. They induce morphologically complex, sometimes spectacular, gall structures on plants in which the larval stages develop. Gall wasps have therefore evolved an intimate association with their hosts - both metabolically, and in terms of their population histories. Gall wasps must both interact physiologically with their hosts to induce galls, and track their host plants through space and time. My thesis centres on two uses of genomic data in understanding the biology of the oak apple gall wasp Biorhiza pallida. I provide a comprehensive investigation into patterns of oak and gall wasp gene expression associated with gall induction, and a population genomic reconstruction of the population history of this species across the Western Palaearctic. While advances in sequencing technology and reduced costs have made these aims possible, analysis of the massive resulting datasets generated creates new challenges. Firstly, in reconstructing the population history of B. pallida, I describe the use of shotgun sequencing and an informatic pipeline to generate alignments of several thousand loci for three B. pallida individuals sampled from putative glacial refugia across the Western Palaearctic in Iberia, the Balkans and Iran. This dataset was analysed using a new maximum likelihood method capable of estimating population splitting and admixture among refugia across very large numbers of loci. The results showed an ancient divide between Iberia and the other two refugia, followed by very recent admixture between easternmost and westernmost regions. This suggests that gall wasps have migrated westwards along the North African coast as well as through mainland Europe. Second, I compare the gene expression profiles of gall wasp and oak tissues sampled from each of three stages of gall development, leading to new insights into potential mechanisms of gall wasp-oak interaction. A highly expressed gall wasp protein was identified that is hypothesised to stimulate somatic embryogenesis-like development of the gall through interaction with oak tissue glycoproteins. Highly expressed oak genes include those coding for nodulin-like proteins similar to those involved in legume nodule formation. Finally, analysis of the gall wasp genome has revealed potential, but as yet unconfirmed, horizontal gene transfer events into gall wasp genomes. Genes discovered in three gall wasp genomes and expressed in three transcriptomes encode plant cell wall degrading enzymes. They are not of hymenopteran origin, and are most homologous to genes of plant pathogenic bacteria. These genes could be involved in several aspects of gall wasp biology, including feeding and developmental manipulation of host plant tissue

DNA methylation differs extensively between strains of the same geographical origin and changes with age in Daphnia magna

Background Patterns of methylation influence lifespan, but methylation and lifespan may also depend on diet, or differ between genotypes. Prior to this study, interactions between diet and genotype have not been explored together to determine their influence on methylation. The invertebrate Daphnia magna is an excellent choice for testing the epigenetic response to the environment: parthenogenetic offspring are identical to their siblings (making for powerful genetic comparisons), they are relatively short lived and have well-characterised inter-strain life-history trait differences. We performed a survival analysis in response to caloric restriction and then undertook a 47-replicate experiment testing the DNA methylation response to ageing and caloric restriction of two strains of D. magna. Results Methylated cytosines (CpGs) were most prevalent in exons two to five of gene bodies. One strain exhibited a significantly increased lifespan in response to caloric restriction, but there was no effect of food-level CpG methylation status. Inter-strain differences dominated the methylation experiment with over 15,000 differently methylated CpGs. One gene, Me31b, was hypermethylated extensively in one strain and is a key regulator of embryonic expression. Sixty-one CpGs were differentially methylated between young and old individuals, including multiple CpGs within the histone H3 gene, which were hypermethylated in old individuals. Across all age-related CpGs, we identified a set that are highly correlated with chronological age. Conclusions Methylated cytosines are concentrated in early exons of gene sequences indicative of a directed, non-random, process despite the low overall DNA methylation percentage in this species. We identify no effect of caloric restriction on DNA methylation, contrary to our previous results, and established impacts of caloric restriction on phenotype and gene expression. We propose our approach here is more robust in invertebrates given genome-wide CpG distributions. For both strain and ageing, a single gene emerges as differentially methylated that for each factor could have widespread phenotypic effects. Our data showed the potential for an epigenetic clock at a subset of age positions, which is exciting but requires confirmation

Low coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp

Population divergence and gene flow are key processes in evolution and ecology. Model‐based analysis of genome‐wide datasets allows discrimination between alternative scenarios for these processes even in non‐model taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum (bSFS), the second on the Pairwise Sequentially Markovian Coalescent (PSMC)) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia, and are separated into northern and southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of post divergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca . 196kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca . 57kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome‐level data for small population samples of tiny, non‐model organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.Output Status: Forthcoming/Available Onlin

An Africa-wide genomic evolution of insecticide resistance in the malaria vector Anopheles funestus involves selective sweeps, copy number variations, gene conversion and transposons

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improve control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweep, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes

From inquilines to gall inducers:Genomic signature of a life-style transition in synergus gall wasps

Gall wasps (Hymenoptera: Cynipidae) induce complex galls on oaks, roses and other plants, but the mechanism of gall induction is still unknown. Here we take a comparative genomic approach to revealing the genetic basis of gall induction. We focus on Synergus itoensis, a species that induces galls inside oak acorns. Previous studies suggested that this species evolved the ability to initiate gall formation recently, as it is deeply nested within the genus Synergus, whose members are mostly inquilines that develop inside the galls of other species. We compared the genome of S. itoensis to that of three related Synergus inquilines to identify genomic changes associated with the origin of gall induction. We used a novel Bayesian selection analysis, which accounts for branch-specific and gene-specific selection effects, to search for signatures of selection in 7,600 single-copy orthologous genes shared by the four Synergus species. We found that the terminal branch leading to S. itoensis had more genes with a significantly elevated dN/dS ratio (Positive Signature Genes, PSGs) than the other terminal branches in the tree; the S. itoensis branch also had more genes with a significantly decreased dN/dS ratio. Gene set enrichment analysis showed that the PSG set of S. itoensis, unlike those of the inquiline species, is enriched in several biological process Gene Ontology terms, the most prominent of which is “Ovarian Follicle Cell Development”. Our results indicate that the origin of gall induction is associated with distinct genomic changes, and provide a good starting point for further characterization of the genes involved

Phylogenomic analysis of protein-coding genes resolves complex gall wasp relationships

International audienceAbstract Gall wasps (Hymenoptera: Cynipidae) comprise 13 distinct tribes whose interrelationships remain incompletely understood. Recent analyses of ultra‐conserved elements (UCEs) represent the first attempt at resolving these relationships using phylogenomics. Here, we present the first analysis based on protein‐coding sequences from genome and transcriptome assemblies. Unlike UCEs, these data allow more sophisticated substitution models, which can potentially resolve issues with long‐branch attraction. We include data for 37 cynipoid species, including two tribes missing in the UCE analysis: Aylacini (s. str.) and Qwaqwaiini. Our results confirm the UCE result that Cynipidae are not monophyletic. Specifically, the Paraulacini and Diplolepidini + Pediaspidini fall outside a core clade (Cynipidae s. str.), which is more closely related to the insect‐parasitic Figitidae, and this result is robust to the exclusion of long‐branch taxa that could mislead the analysis. Given this, we here divide the Cynipidae into three families: the Paraulacidae stat. prom., Diplolepididae stat. prom. and Cynipidae (s. str.). Our results suggest that the Eschatocerini are the sister group of the remaining Cynipidae (s. str.). Within the Cynipidae (s. str.), the Aylacini (s. str.) are more closely related to oak gall wasps (Cynipini) and some of their inquilines (Ceroptresini) than to other herb gallers (Aulacideini and Phanacidini), and the Qwaqwaiini likely form a clade together with Synergini (s. str.) and Rhoophilini. Several alternative scenarios for the evolution of cynipid life histories are compatible with the relationships suggested by our analysis, but all are complex and require multiple shifts among parasitoids, inquilines and gall inducers