Transcriptome profiling during a natural host-parasite interaction

Host behaviour alteration by its parasite: from brain gene expression to functional test

Many parasites with complex life cycles modify their intermediate hosts' behaviour, presumably to increase transmission to their final host. The threespine stickleback (Gasterosteus aculeatus) is an intermediate host in the cestode Schistocephalus solidus life cycle, which ends in an avian host, and shows increased risky behaviours when infected. We studied brain gene expression profiles of sticklebacks infected with S. solidus to determine the proximal causes of these behavioural alterations. We show that infected fish have altered expression levels in genes involved in the inositol pathway. We thus tested the functional implication of this pathway and successfully rescued normal behaviours in infected sticklebacks using lithium exposure. We also show that exposed but uninfected fish have a distinct gene expression profile from both infected fish and control individuals, allowing us to separate gene activity related to parasite exposure from consequences of a successful infection. Finally, we find that selective serotonin reuptake inhibitor-treated sticklebacks and infected fish do not have similarly altered gene expression, despite their comparable behaviours, suggesting that the serotonin pathway is probably not the main driver of phenotypic changes in infected sticklebacks. Taken together, our results allow us to predict that if S. solidus directly manipulates its host, it could target the inositol pathway.

Transcriptome and Proteome Response of Rhipicephalus annulatus Tick Vector to Babesia bigemina Infection

A system biology approach was used to gain insight into tick biology and interactions between vector and pathogen. Rhipicephalus annulatus is one of the main vectors of Babesia bigemina which has a massive impact on animal health. It is vital to obtain more information about this relationship, to better understand tick and pathogen biology, pathogen transmission dynamics, and new potential control approaches. In ticks, salivary glands (SGs) play a key role during pathogen infection and transmission. RNA sequencing obtained from uninfected and B. bigemina infected SGs obtained from fed female ticks resulted in 6823 and 6475 unigenes, respectively. From these, 360 unigenes were found to be differentially expressed (p < 0.05). Reversed phase liquid chromatography-mass spectrometry identified a total of 3679 tick proteins. Among them 406 were differently represented in response to Babesia infection. The omics data obtained suggested that Babesia infection lead to a reduction in the levels of mRNA and proteins (n = 237 transcripts, n = 212 proteins) when compared to uninfected controls. Integrated transcriptomics and proteomics datasets suggested a key role for stress response and apoptosis pathways in response to infection. Thus, six genes coding for GP80, death-associated protein kinase (DAPK-1), bax inhibitor-1 related (BI-1), heat shock protein (HSP), heat shock transcription factor (PHSTF), and queuine trna-ribosyltransferase (QtRibosyl) were selected and RNA interference (RNAi) performed. Gene silencing was obtained for all genes except phstf. Knockdown of gp80, dapk-1, and bi-1 led to a significant increase in Babesia infection levels while hsp and QtRibosyl knockdown resulted in a non-significant decrease of infection levels when compared to the respective controls. Gene knockdown did not affect tick survival, but engorged female weight and egg production were affected in the gp80, dapk-1, and QtRibosyl-silenced groups in comparison to controls. These results advanced our understanding of tick-Babesia molecular interactions, and suggested new tick antigens as putative targets for vaccination to control tick infestations and pathogen infection/transmission.

Daphnia galeata responds to the exposure to an ichthyosporean gut parasite by down-regulation of immunity and lipid metabolism

BACKGROUND

Regulatory circuits of infection in the emerging experimental model system, water flea Daphnia and their microparasites, remain largely unknown. Here we provide the first molecular insights into the response of Daphnia galeata to its highly virulent and common parasite Caullerya mesnili, an ichthyosporean that infects the gut epithelium. We generated a transcriptomic dataset using RNAseq from parasite-exposed (vs. control) Daphnia, at two time points (4 and 48 h) after parasite exposure.

RESULTS

We found a down-regulation of metabolism and immunity-related genes, at 48 h (but not 4 h) after parasite exposure. These genes are involved in lipid metabolism and fatty acid biosynthesis, as well as microbe recognition (e.g. c-type lectins) and pathogen attack (e.g. gut chitin).

CONCLUSIONS

General metabolic suppression implies host energy shift from reproduction to survival, which is in agreement with the known drastic reduction in Daphnia fecundity after Caullerya infection. The down-regulation of gut chitin indicates a possible interaction between the peritrophic matrix and the evading host immune system. Our study provides the first description of host transcriptional responses in this very promising host-parasite experimental system.

A transcriptomic snapshot of early molecular communication between Pasteuria penetrans and Meloidogyne incognita

BACKGROUND

Southern root-knot nematode Meloidogyne incognita (Kofoid and White, 1919), Chitwood, 1949 is a key pest of agricultural crops. Pasteuria penetrans is a hyperparasitic bacterium capable of suppressing the nematode reproduction, and represents a typical coevolved pathogen-hyperparasite system. Attachment of Pasteuria endospores to the cuticle of second-stage nematode juveniles is the first and pivotal step in the bacterial infection. RNA-Seq was used to understand the early transcriptional response of the root-knot nematode at 8 h post Pasteuria endospore attachment.

RESULTS

A total of 52,485 transcripts were assembled from the high quality (HQ) reads, out of which 582 transcripts were found differentially expressed in the Pasteuria endospore encumbered J2 s, of which 229 were up-regulated and 353 were down-regulated. Pasteuria infection caused a suppression of the protein synthesis machinery of the nematode. Several of the differentially expressed transcripts were putatively involved in nematode innate immunity, signaling, stress responses, endospore attachment process and post-attachment behavioral modification of the juveniles. The expression profiles of fifteen selected transcripts were validated to be true by the qRT PCR. RNAi based silencing of transcripts coding for fructose bisphosphate aldolase and glucosyl transferase caused a reduction in endospore attachment as compared to the controls, whereas, silencing of aspartic protease and ubiquitin coding transcripts resulted in higher incidence of endospore attachment on the nematode cuticle.

CONCLUSIONS

Here we provide evidence of an early transcriptional response by the nematode upon infection by Pasteuria prior to root invasion. We found that adhesion of Pasteuria endospores to the cuticle induced a down-regulated protein response in the nematode. In addition, we show that fructose bisphosphate aldolase, glucosyl transferase, aspartic protease and ubiquitin coding transcripts are involved in modulating the endospore attachment on the nematode cuticle. Our results add new and significant information to the existing knowledge on early molecular interaction between M. incognita and P. penetrans.

Genomic evidence for population-specific responses to co-evolving parasites in a New Zealand freshwater snail

Reciprocal co-evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population- or genotype-specific evolutionary change. These host-parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co-evolving host-parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co-evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus-infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population-level co-evolution is leading to population-specific host-parasite interactions and evolutionary trajectories. We also used an F_ST -based approach to identify a set of loci that represent promising candidates for targets of parasite-mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population-specific responses to co-evolving infection in the P. antipodarum-Microphallus interaction and provide new insights into the genomic basis of co-evolutionary interactions in nature.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

The genetic basis of resistance and matching-allele interactions of a host-parasite system: The Daphnia magna-Pasteuria ramosa model

Negative frequency-dependent selection (NFDS) is an evolutionary mechanism suggested to govern host-parasite coevolution and the maintenance of genetic diversity at host resistance loci, such as the vertebrate MHC and R-genes in plants. Matching-allele interactions of hosts and parasites that prevent the emergence of host and parasite genotypes that are universally resistant and infective are a genetic mechanism predicted to underpin NFDS. The underlying genetics of matching-allele interactions are unknown even in host-parasite systems with empirical support for coevolution by NFDS, as is the case for the planktonic crustacean Daphnia magna and the bacterial pathogen Pasteuria ramosa. We fine-map one locus associated with D. magna resistance to P. ramosa and genetically characterize two haplotypes of the Pasteuria resistance (PR-) locus using de novo genome and transcriptome sequencing. Sequence comparison of PR-locus haplotypes finds dramatic structural polymorphisms between PR-locus haplotypes including a large portion of each haplotype being composed of non-homologous sequences resulting in haplotypes differing in size by 66 kb. The high divergence of PR-locus haplotypes suggest a history of multiple, diverse and repeated instances of structural mutation events and restricted recombination. Annotation of the haplotypes reveals striking differences in gene content. In particular, a group of glycosyltransferase genes that is present in the susceptible but absent in the resistant haplotype. Moreover, in natural populations, we find that the PR-locus polymorphism is associated with variation in resistance to different P. ramosa genotypes, pointing to the PR-locus polymorphism as being responsible for the matching-allele interactions that have been previously described for this system. Our results conclusively identify a genetic basis for the matching-allele interaction observed in a coevolving host-parasite system and provide a first insight into its molecular basis.

Transcriptome analysis of Meloidogyne incognita encumbered by Pasteuria penetrans endospores provides new insights into bacteria and nematode interaction

Root-knot nematodes (RKN,Meloidogynespp.) are one of the most harmful pests in agriculture.Pasteuria penetrans, an obligate hyperparasite of RKN, is an effective biological control agent. However, little is known about the molecular mechanisms of interaction betweenP. penetransand the second-stage juvenile (J2) ofMeloidogyne incognita. Here, we used transcriptome sequencing to characterise the differential expression profiles between control J2 ofM. incognitaand J2 encumbered byP. penetransendospores. A total of 445 genes were found to be differentially expressed, including 406 up-regulated and 39 down-regulated genes. Thirty-seven putative immune-related genes encoding collagens, cytochrome P450, lysozymes and other active proteins were identified. Genes involved in the ‘biosynthesis of unsaturated fatty acids’ pathway and several core sets of immune effectors were up-regulated, indicating conserved immune mechanisms among different nematodes. Down-regulation of cytochrome P450-related genes might suggest a specific defence response ofM. incognitaencumbered byP. penetransendospores.

Reproductive switching analysis of Daphnia similoides between sexual female and parthenogenetic female by transcriptome comparison

The water flea Daphnia are planktonic crustaceans commonly found in freshwater environment that can switch their reproduction mode from parthenogenesis to sexual reproduction to adapt to the external environment. As such, Daphnia are great model organisms to study the mechanism of reproductive switching, the underlying mechanism of reproduction and development in cladocerans and other animals. However, little is known about the Daphnia's reproductive behaviour at a molecular level. We constructed a genetic database of the genes expressed in a sexual female (SF) and a parthenogenetic female (PF) of D. similoides using Illumina HiSeq 2500. A total of 1,763 differentially expressed genes (865 up- and 898 down-regulated) were detected in SF. Of the top 30 up-regulated SF unigenes, the top 4 unigenes belonged to the Chitin_bind_4 family. In contrast, of the top down-regulated SF unigenes, the top 3 unigenes belonged to the Vitellogenin_N family. This is the first study to indicate genes that may have a crucial role in reproductive switching of D. similoides, which could be used as candidate genes for further functional studies. Thus, this study provides a rich resource for investigation and elucidation of reproductive switching in D. similoides.

Mechanisms Underlying Adaptation to Life in Hydrogen Sulfide-Rich Environments

Hydrogen sulfide (H2S) is a potent toxicant interfering with oxidative phosphorylation in mitochondria and creating extreme environmental conditions in aquatic ecosystems. The mechanistic basis of adaptation to perpetual exposure to H2S remains poorly understood. We investigated evolutionarily independent lineages of livebearing fishes that have colonized and adapted to springs rich in H2S and compared their genome-wide gene expression patterns with closely related lineages from adjacent, nonsulfidic streams. Significant differences in gene expression were uncovered between all sulfidic and nonsulfidic population pairs. Variation in the number of differentially expressed genes among population pairs corresponded to differences in divergence times and rates of gene flow, which is consistent with neutral drift driving a substantial portion of gene expression variation among populations. Accordingly, there was little evidence for convergent evolution shaping large-scale gene expression patterns among independent sulfide spring populations. Nonetheless, we identified a small number of genes that was consistently differentially expressed in the same direction in all sulfidic and nonsulfidic population pairs. Functional annotation of shared differentially expressed genes indicated upregulation of genes associated with enzymatic H2S detoxification and transport of oxidized sulfur species, oxidative phosphorylation, energy metabolism, and pathways involved in responses to oxidative stress. Overall, our results suggest that modification of processes associated with H2S detoxification and toxicity likely complement each other to mediate elevated H2S tolerance in sulfide spring fishes. Our analyses allow for the development of novel hypotheses about biochemical and physiological mechanisms of adaptation to extreme environments.