Molecular alterations during larval development of Haemonchus contortus in vitro are under tight post-transcriptional control☆
Graphical abstract
Introduction
The free-living nematode Caenorhabditis elegans is one of the best-characterised metazoan organisms (Corsi et al., 2015). Integrated analyses of molecular data sets for this nematode have made it possible to gain deep, global insights into areas such as the developmental and reproductive biology, physiology, biochemistry, genetics and neurobiology of metazoans (Piano et al., 2006, van Assche et al., 2015). By comparison, little is known about these areas for parasitic nematodes. Since C. elegans is relatively closely related to the socioeconomically important parasitic nematodes of the Order Strongylida (Blaxter et al., 1998), the information and data sets available for C. elegans can significantly assist molecular investigations of strongylid nematodes (Gilleard, 2004, Howe et al., 2016a, Howe et al., 2016b). There have been major advances in the genomics and transcriptomics of worms such as Haemonchus contortus (barber’s pole worm), Necator americanus (hookworm) and Oesophagostomum dentatum (nodule worm) (see Laing et al., 2013, Schwarz et al., 2013, Tang et al., 2014, Tyagi et al., 2015), allowing novel insights into the systems biology of strongylid nematodes.
Haemonchus contortus is one of the most important parasitic nematodes of livestock animals worldwide (cf. Gasser and von Samson-Himmelstjerna, 2016). This parasite is transmitted orally from contaminated pasture to the host through a relatively complex, 3 week life cycle (Veglia, 1915): eggs are excreted in host faeces; the L1s develop inside the eggs and then hatch (within 1 day) and moult to the L2 and L3 stages within 1 week; the infective L3s are then ingested by the host, exsheath and, following a histotropic phase, develop through the L4 stage to dioecious adults which both feed on blood from vessels in the stomach wall.
Although some studies have explored changes or differences in mRNAs or proteins in H. contortus (see Campbell et al., 2008, Cantacessi et al., 2010, Laing et al., 2013, Schwarz et al., 2013, Wang et al., 2016), no investigation has yet integrated or compared transcriptomic and proteomic profiles during developmental transitions in this nematode or explored regulatory molecular processes. For instance, while small RNAs, such as microRNAs (miRNAs) (Kim, 2005, Grosshans and Filipowicz, 2008), play key roles in regulating genes via mRNA silencing or chromatin modification (Kim, 2005, Grosshans and Filipowicz, 2008), limited information is available for H. contortus (see Winter et al., 2012, Gillan et al., 2017, Gu et al., 2017), particularly in relation to this parasite’s developmental biology. Here, we explored mRNA, miRNA and protein alterations in H. contortus during the developmental transition from ensheathed L3 to xL3 and then to the L4 stage in vitro, utilising a combination of transcriptomic, proteomic and bioinformatic tools. The goal was to identify key molecular changes associated with this crucial phase of this nematode’s development using a data integration approach.
Section snippets
Procurement of parasite and storage
Haemonchus contortus (Haecon-5 strain) was maintained in experimental sheep (cf. Schwarz et al., 2013), in accordance with the institutional animal ethics guidelines (permit no. 1613878; The University of Melbourne, Australia). In brief, helminth-free Merino sheep (8 weeks of age) were inoculated with 7,000 L3s of H. contortus. Four weeks p.i., faecal samples were collected each day; L3s were produced by incubating faeces at 27 °C for 1 week (cf. Schwarz et al., 2013), purified through two
Transcriptome, miRNAome and proteome of L3s, xL3s and L4s of H. contortus
Three datasets (i.e. transcriptome, miRNAome and proteome) were produced for each L3, xL3 and L4 of H. contortus. For the transcriptome, a total of 9,754 mRNA transcripts were identified amongst >700 million pair-end reads for all three developmental stages (Supplementary Table S1). The numbers of mRNA transcripts detected in the L3, xL3 and L4 stages were 9,492, 9,581 and 9,487, respectively (Supplementary Table S1). Small RNA libraries yielded >400 million reads, with 88 miRNAs identified (
Discussion
Using transcriptomic, miRNA and proteomic datasets, we explored molecular alterations in the developmental switch from L3 to L4 of H. contortus in vitro. Key changes were seen in mRNA, miRNA and protein profiles, and plausible relationships were established between the transcription (mRNA and/or miRNA) and protein expression.
Qualitative and quantitative analyses performed in this study identified substantial molecular changes in the developmental switch of H. contortus. The molecular changes
Acknowledgements
Funding from the National Health and Medical Research Council (NHMRC) of Australia, the Australian Research Council, Melbourne Water Corporation and The University of Melbourne (BIP) is gratefully acknowledged (R.B.G.). Support from the Melbourne Bioinformatics Platform is gratefully acknowledged. P.K.K. holds an NHMRC Early Career Researcher Fellowship (ECRF), and N.D.Y. holds an NHMRC Career Development Fellowship (CDF1).
References (77)
MicroRNAs: target recognition and regulatory functions
Cell
(2009)- et al.
Transforming growth factor-beta and insulin-like signalling pathways in parasitic helminths
Int. J. Parasitol.
(2002) - et al.
Caenorhabditis elegans as a model for parasitic nematodes
Int. J. Parasitol.
(1998) - et al.
Gender-enriched transcripts in Haemonchus contortus–predicted functions and genetic interactions based on comparative analyses with Caenorhabditis elegans
Int. J. Parasitol.
(2008) The dauer hypothesis and the evolution of parasitism: 20 years on and still going strong
Int. J. Parasitol.
(2014)- et al.
Caenorhabditis elegans: how good a model for veterinary parasites?
Vet. Parasitol.
(2001) Small RNAs: classification, biogenesis, and function
Mol. Cell
(2005)- et al.
Systematic analysis of tissue-restricted miRISCs reveals a broad role for microRNAs in suppressing basal activity of the C. elegans pathogen response
Mol. Cell
(2012) - et al.
Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets
Cell
(2005) - et al.
Prospects for exploring molecular developmental processes in Haemonchus contortus
Int. J. Parasitol.
(2006)
A conserved endocrine mechanism controls the formation of dauer and infective larvae in nematodes
Curr. Biol.
Caenorhabditis elegans as a model for parasitic nematodes: a focus on the cuticle
Parasitol. Today
Low cost whole-organism screening of compounds for anthelmintic activity
Int. J. Parasitol.
Cracking the nodule worm code advances knowledge of parasite biology and biotechnology to tackle major diseases of livestock
Biotechnol. Adv.
Integrating -omics: systems biology as explored through C. elegans research
J. Mol. Biol.
How can we understand the genomic basis of nematode parasitism?
Trends Parasitol.
Identification of differentially expressed proteins between free-living and activated third-stage larvae of Haemonchus contortus
Vet. Parasitol.
Conserved miRNAs are candidate post-transcriptional regulators of developmental arrest in free-living and parasitic nematodes
Genome Biol. Evol.
miARma-Seq: a comprehensive tool for miRNA, mRNA and circRNA analysis
Sci. Rep.
Chemosensation in C. elegans
WormBook
Co-option of the hormone-signalling module dafachronic acid-DAF-12 in nematode evolution
Nature
A molecular evolutionary framework for the phylum Nematoda
Nature
Trimmomatic: a flexible trimmer for Illumina sequence data
Bioinformatics
Small-molecule pheromones and hormones controlling nematode development
Nat. Chem. Biol.
Differences in transcription between free-living and CO2-activated third-stage larvae of Haemonchus contortus
BMC Genomics
Zinc fingers in Caenorhabditis elegans: finding families and probing pathways
Science
A transparent window into biology: a primer on Caenorhabditis elegans
WormBook
Accounting for experimental noise reveals that mRNA levels, amplified by post-transcriptional processes, largely determine steady-state protein levels in yeast
PLoS Genet.
Regulation of the mammalian heat shock factor 1
FEBS J.
Transcriptional changes in the hookworm, Ancylostoma caninum, during the transition from a free-living to a parasitic larva
PLoS Negl. Trop. Dis.
A genome resequencing-based genetic map reveals the recombination landscape of an outbred parasitic nematode in the presence of polyploidy and polyandry
Genome Biol. Evol.
Gene-specific correlation of RNA and protein levels in human cells and tissues
Mol. Syst. Biol.
Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?
Nat. Rev. Genet.
Comparing Pearson, Spearman and Hoeffding's D measure for gene expression association analysis
J. Bioinform. Comput. Biol.
Haemonchus contortus and Haemonchosis – Past, Present and Future Trends
Increased expression of a microRNA correlates with anthelmintic resistance in parasitic nematodes
Front. Cell. Infect. Microbiol.
The use of Caenorhabditis elegans in parasitic nematode research
Parasitology
Full-length transcriptome assembly from RNA-Seq data without a reference genome
Nat. Biotechnol.
Cited by (31)
An informatic workflow for the enhanced annotation of excretory/secretory proteins of Haemonchus contortus
2023, Computational and Structural Biotechnology JournalUbiquitination pathway model for the barber's pole worm – Haemonchus contortus
2022, International Journal for ParasitologyHelminth lipidomics: Technical aspects and future prospects
2021, Current Research in Parasitology and Vector-Borne DiseasesCitation Excerpt :This haematophagous (‘barberʼs pole’) worm infects hundreds of millions of small ruminants (e.g. sheep and goats) worldwide, causing worldwide economic losses estimated at $3.3 billion per annum (Evans & Chapple, 2002; Gasser & von Samson-Himmelstjerna, 2016). Genomic (Laing et al., 2013; Schwarz et al., 2013; Palevich et al., 2019; Doyle et al., 2018, 2020), transcriptomic (Salle et al., 2018; Ma et al., 2018b, 2019; Stasiuk et al., 2019) and proteomic (Wang et al., 2019a, b, 2020b) investigations of H. contortus have not only detected a wide range of molecules (i.e. genes and proteins) involved in lipid metabolism and transportation, but also revealed marked alterations in transcription or expression of selected lipid-associated genes and proteins during the development, which could relate to dynamic changes in the lipidome of H. contortus. Some studies have taken the first step to study the lipidomes of H. contortus and its stages, organ systems and tissues.
N-glycome and N-glycoproteome of a hematophagous parasitic nematode Haemonchus
2021, Computational and Structural Biotechnology JournalToward integrative ‘omics of the barber's pole worm and related parasitic nematodes
2020, Infection, Genetics and EvolutionCitation Excerpt :In addition, detailed transcriptomic analysis has been widely employed in the areas of developmental biology, drug discovery, drug resistance and vaccine development (Ma et al., 2018; Preston et al., 2017; Sallé et al., 2018; Stasiuk et al., 2019). For instance, mRNA profiling has been conducted for L3s, exsheathed L3s and in vitro-cultured L4s of H. contortus (Ma et al., 2018), representing a comprehensive transcriptome dataset for the developmental transition from free-living to parasitic stage of this parasitic nematode. Specifically, molecular responses to compound (i.e., deguelin) exposure, anthelmintic (i.e., albendazole, mebendazole, thiabendazole and oxfendazole) treatment and vaccination (i.e., Barbervax®) have also been explored in H. contortus at the transcriptional level (Preston et al., 2017; Sallé et al., 2018; Stasiuk et al., 2019), providing some insights into the drug effects/metabolism and valuable information for vaccine improvement.
The developmental phosphoproteome of Haemonchus contortus
2020, Journal of Proteomics
- ☆
Note: Nucleic acid sequence data from this study have been deposited in the National Center for Biotechnology Information (NCBI) sequence read archive (SRA) – http://www.ncbi.nlm.nih.gov/sra/ under accession number SRP136037.