Non-coding RNA in the gut of the blood-feeding parasitic worm, Haemonchus contortus.

The intestine of Haemonchus contortus is an essential tissue that has been indicated to be a major target for the prevention of haemonchosis caused by this parasitic nematode of small ruminants. Biological peculiarities of the intestine warrant in-depth exploitation, which can be leveraged for future disease control efforts. Here, we determined the intestinal ncRNA (lncRNA, circRNA and miRNA) atlas using whole-transcriptome sequencing and bioinformatics approaches. In total, 4846 novel lncRNA, 982 circRNA, 96 miRNA (65 known and 31 novel) and 8821 mRNA were identified from the H. contortus intestine. The features of lncRNA, circRNA and miRNA were fully characterized. Comparison of miRNA from the intestines and extracellular vesicles supported the speculation that the miRNA from the latter were of intestinal origin in H. contortus. Further function analysis suggests that the cis-lncRNA targeted genes were involved in protein binding, intracellular anatomical structure, organelle and cellular process, whereas the circRNA parental genes were mainly enriched in molecular function categories, such as ribonucleotide binding, nucleotide binding, ATP binding and carbohydrate derivative binding. The miRNA target genes were related to the cellular process, cellular response to stimulus, cellular protein modification process and signal transduction. Moreover, competing endogenous RNA network analysis revealed that the majority of lncRNA, circRNA and mRNA only have one or two binding sites with specific miRNA. Lastly, randomly selected circRNA, lncRNA and miRNA were verified successfully using RT-PCR. Collectively, these data provide the most comprehensive compilation of intestinal transcripts and their functions, and it will be helpful to decipher the biological and molecular complexity of the intestine and lay the foundation for further functional research.

A serine/threonine-specific protein kinase of Haemonchus contortus with a role in the development.

In the free-living nematode Caenorhabditis elegans, the serine/threonine-specific protein kinase, AKT, is known to play a key role in dauer formation, life-span, and stress-resistance through the insulin-like signaling pathway. Although the structure and function of AKT-coding genes of C. elegans are understood, this is not the case for homologous genes in parasitic nematodes. In the present study, we explored a C. elegans akt-1 gene homolog in the parasitic nematode Haemonchus contortus, investigated its transcript isoforms (Hc-akt-1a and Hc-akt-1b), and studied expression and function using both homologous and heterologous functional genomic tools. In C. elegans, we showed that the predicted promoter of Hc-akt-1 drives substantial expression in ASJ neurons of the N2 (wild-type) strain. In H. contortus (Haecon-5 stain), RNAi (soaking) led to a significantly decreased transcript abundance for both Hc-akt-1a and Hc-akt-1b, and reduced larval development in larval stages in vitro. Chemical inhibition was also shown to block larval development. Taken together, the evidence from this study points to a key functional role for Hc-akt-1 in H. contortus.

NPM1 mutation reprograms leukemic transcription network via reshaping TAD topology.

C-terminal mutation of Nucleophosmin 1 (NPM1C+) was thought to be a primary driving event in acute myeloid leukemia (AML) that reprograms leukemic-associated transcription programs to transform hematopoietic stem and progenitor cells (HSPCs). However, molecular mechanisms underlying NPM1C+-driven leukemogenesis remain elusive. Here, we report that NPM1C+ activates signature HOX genes and reprograms cell cycle regulators by altering CTCF-driven topologically associated domains (TADs). Hematopoietic-specific NPM1C+ knock-in alters TAD topology leading to disrupted regulation of the cell cycle as well as aberrant chromatin accessibility and homeotic gene expression, which results in myeloid differentiation block. Restoration of NPM1 within the nucleus re-establishes differentiation programs by reorganizing TADs critical for myeloid TFs and cell cycle regulators that switch the oncogenic MIZ1/MYC regulatory axis in favor of interacting with coactivator NPM1/p300, and prevents NPM1C+-driven leukemogenesis. In sum, our data reveal that NPM1C+ reshapes CTCF-defined TAD topology to reprogram signature leukemic transcription programs required for cell cycle progression and leukemic transformation.

Comparative analysis on transcriptomics of ivermectin resistant and susceptible strains of Haemonchus contortus.

BACKGROUND: Ivermectin (IVM) is one of the most important and widely used anthelmintics in veterinary medicine. However, its efficacy is increasingly compromised by widespread resistance, and the exact mechanism of IVM resistance remains unclear for most parasitic nematodes, including Haemonchus contortus, a blood-sucking parasitic nematode of small ruminants. METHODS: In this study, an H. contortus IVM-resistant strain from Zhaosu, Xinjiang, China, was isolated and assessed by the control test, faecal egg count reduction test (FECRT) and the larval development assay (LDA). Subsequently, comparative analyses on the transcriptomics of IVM-susceptible and IVM-resistant adult worms of this parasite were carried out using RNA sequencing (RNA-seq) and bioinformatics. RESULTS: In total, 543 (416 known, 127 novel) and 359 (309 known, 50 novel) differentially expressed genes (DEGs) were identified in male and female adult worms of the resistant strain compared with those of the susceptible strain, respectively. In addition to several previously known candidate genes which were supposed to be associated with IVM resistance and whose functions were involved in receptor activity, transport, and detoxification, we found some new potential target genes, including those related to lipid metabolism, structural constituent of cuticle, and important pathways such as antigen processing and presentation, lysosome, autophagy, apoptosis, and NOD1-like receptor signalling pathways. Finally, the results of quantitative real-time polymerase chain reaction confirmed that the transcriptional profiles of selected DEGs (male: 8 genes, female: 10 genes) were consistent with those obtained by the RNA-seq. CONCLUSIONS: Our results indicate that IVM has multiple effects, including both neuromuscular and non-neuromuscular targets, and provide valuable information for further studies on the IVM resistance mechanism in H. contortus.

A transcription factor DAF-5 functions in Haemonchus contortus development.

BACKGROUND: Abnormal dauer formation gene (daf-5), located downstream of the DAF-7 signalling pathway, mainly functions in dauer formation and reproductive processes in the free-living nematode Caenorhabditis elegans. Although the structure and function of daf-5 have been clarified in C. elegans, they still remain totally unknown in Haemonchus contortus, a socio-economically important parasitic nematode of gastric ruminants. METHODS: A homologue of daf-5, Hc-daf-5, and its inferred product (Hc-DAF-5) in H. contortus were identified and characterized in this study. Then the transcriptional profiles of Hc-daf-5 and the anatomical expression of Hc-DAF-5 in H. contortus were studied using an integrated molecular approach. RNA interference (RNAi) was performed to explore its function in transition from the exsheathed third-stage larvae (xL3s) to the fourth-stage larvae (L4s) in vitro. Finally, the interaction between Hc-DAF-5 and Hc-DAF-3 (a co-Smad) was detected by bimolecular fluorescence complementation (BiFc) in vitro. RESULTS: It was shown that Hc-DAF-5 was a member of the Sno/Ski superfamily. Hc-daf-5 was transcribed in all developmental stages of H. contortus, with significant upregulation in L3s. Native Hc-DAF-5 was localized in the reproductive organs, cuticle, and intestine via immunohistochemistry. RNAi revealed that specific small interfering RNAs (siRNAs) could retard xL3 development. In addition, the interaction between Hc-DAF-5 and Hc-DAF-3 indicated that the SDS box of Hc-DAF-5 was dispensable for the binding of Hc-DAF-5 to Hc-DAF-3, and the MH2 domain was the binding region between Hc-DAF-3 and Hc-DAF-5. CONCLUSIONS: In summary, these findings show that Hc-daf-5 functions in the developmental processes of H. contortus, and this study is the first attempt to characterize the daf-5 gene in parasitic nematodes.

Genome-Wide Identification of CircRNAs of Infective Larvae and Adult Worms of Parasitic Nematode, Haemonchus contortus.

CircRNAs, a novel class of ncRNA family, are endogenous transcriptional products involved in various biological and physiological processes in plants and animals. However, almost no information is available for circRNAs of parasitic helminths. In the present study, the circRNAs repertoire was comprehensively explored in Haemonchus contortus, a blood-sucking parasitic nematode of ruminants. In total, 20073 circRNAs were identified and annotated from three key developmental stages/genders of H. contortus including the free-living infective third-stage larvae (L3, 18883), parasitic adult female (Af, 3491), and male worms (Am, 2550) via deep-sequencing technology and bioinformatic analysis. Among these identified circRNAs, 71% were derived from exonic regions of protein-coding genes. The number of circRNAs transcribed from the X chromosome (4704) was higher than that from Chromosome I-V (3143, 3273, 3041, 3030, 2882). The amount of highly expressed circRNAs in third-stage larvae was significantly more abundant than that in adult stage. 15948 and 16847 circRNAs were differentially expressed between Af and L3s and between Am and L3, respectively. Among them, 13409 circRNAs existed in both comparisons. Furthermore, 1119 circRNAs were differentially expressed between Af_and_Am. GO enrichment analysis indicated that source genes of circRNAs differentially expressed between Am and L3 as well as between Af and L3 were significantly enriched in many biological processes, primarily including signaling, signal transduction and cell communication terms. KEGG analysis revealed that parental genes of differentially expressed circRNAs were mainly related to metabolism (pyruvate metabolism, glycerophospholipid metabolism, and carbon metabolism), MAPK signaling pathway, and phosphatidylinositol signaling system. Moreover, many circRNAs contained one or more miRNA potential binding sites, suggesting that they could regulate gene expression at the post-transcriptional level. Furthermore, the correctness of head-to-tail back splicing site and alternative circularization events were verified by Sanger sequencing using both divergent and convergent primers. Finally, the reliability of RNA-Seq data and the resistance of circRNAs to RNase R digestion were confirmed by quantitative RT-PCR. Taken together, our findings provide a foundation for elucidating the regulatory mechanisms of circRNAs in H. contortus, which will advance the understanding of circRNAs in parasitic nematodes.

Identification and characterization of an R-Smad homologue (Hco-DAF-8) from Haemonchus contortus.

BACKGROUND: Smad proteins are essential cellular mediators within the transforming growth factor-ß (TGF-ß) superfamily. They directly transmit incoming signals from the cell surface receptors to the nucleus. In spite of their functional importance, almost nothing is known about Smad proteins in parasitic nematodes including Haemonchus contortus, an important blood-sucking nematode of small ruminants. METHODS: Based on genomic and transcriptome data for H. contortus and using bioinformatics methods, a Smad homologue (called Hco-daf-8) was inferred from H. contortus and the structural characteristics of this gene and its encoded protein Hco-DAF-8 established. Using real-time PCR and immunofluorescence assays, temporal transcriptional and spatial expression profiles of Hco-daf-8 were studied. Gene rescue in Caenorhabditis elegans was then applied to assess the function of Hco-daf-8 and a specific inhibitor of human Smad3 (called SIS3) was employed to evaluate the roles of Hco-DAF-8 in H. contortus development. RESULTS: The features of Hco-DAF-8 (502 amino acids), including conserved R-Smad domains and residues of the L3-loop that determine pathway specificity, are consistent with a TGF-ß type I receptor-activated R-Smad. The Hco-daf-8 gene was transcribed in all developmental stages of H. contortus studied, with a higher level of transcription in the fourth-stage larval (L4) females and the highest level in adult males. Hco-DAF-8 was expressed in the platymyarian muscular cells, intestine and reproductive system of adult stages. Gene rescue experiments showed that Hco-daf-8 was able to partially rescue gene function in a daf-8 deficient mutant strain of C. elegans, leading to a resumption of normal development. In H. contortus, SIS3 was shown to affect H. contortus development from the exsheathed third-stage larvae (L3s) to L4s in vitro. CONCLUSIONS: These findings suggest that Hco-DAF-8, encoded by the gene Hco-daf-8, is an important cellular mediator of H. contortus development via the TGF-ß signalling pathway. They provide a basis for future explorations of Hco-DAF-8 and associated pathways in H. contortus and other important parasitic nematodes.

A daf-7-related TGF-ß ligand (Hc-tgh-2) shows important regulations on the development of Haemonchus contortus.

BACKGROUND: In most multicellular organisms, the transforming growth factor-ß (TGF-ß) signalling pathway is involved in regulating the growth and stem cell differentiation. Previous studies have demonstrated the importance of three key molecules in this pathway in the parasitic nematode Haemonchus contortus, including one TGF-ß type I receptor (Hc-tgfbr1), one TGF-ß type II receptor (Hc-tgfbr2), and one co-Smad (Hc-daf-3), which regulated the developmental transition from the free-living to the parasitic stages of this parasite. However, almost nothing is known about the function of the TGF-ß ligand (Hc-tgh-2) of H. contortus. METHODS: Here, the temporal transcription profiles of Hc-tgh-2 at eight different developmental stages and spatial expression patterns of Hc-TGH-2 in adult female and male worms of H. contortus have been examined by real-time PCR and immunohistochemistry, respectively. In addition, RNA interference (RNAi) by soaking was employed to assess the importance of Hc-tgh-2 in the development from exsheathed third-stage larvae (xL3s) to fourth-stage larvae (L4s) in H. contortus. RESULTS: Hc-tgh-2 was continuously transcribed in all eight developmental stages of H. contortus studied with the highest level in the infective third-stage larvae (iL3) and Hc-TGH-2 was located in the muscle of the body wall, intestine, ovary of adult females and testes of adult males. Silencing Hc-tgh-2 by the specific double-stranded RNA (dsRNA), decreased the transcript level of Hc-tgh-2 and resulted in fewer xL3s developing to L4s in vitro. CONCLUSIONS: These results suggested that the TGF-ß ligand, Hc-TGH-2, could play important roles in the developmental transition from the free-living (L3s) to the parasitic stage (L4s). Furthermore, it may also take part in the processes such as digestion, absorption, host immune response and reproductive development in H. contortus adults.

Development and evaluation of a loop-mediated isothermal amplification (LAMP) assay for the detection of the E198A SNP in the isotype-1 ß-tubulin gene of Haemonchus contortus populations in China.

Haemonchus contortus is one of the most important gastrointestinal nematodes of small ruminants around the world, seriously hampering the healthy development of the sheep industry. The control of this parasite mainly depends on anthelmintics, however, drug resistance of H. contortus has become a serious problems worldwide. Previous studies demonstrated that the E198A (GAA to GCA), a single nucleotide polymorphism (SNP) in the isotype-1 ß-tubulin gene is associated with benzimidazole resistance in H. contortus. However, only PCR-RFLP and ARMS-PCR methods have been previously used for the detection of the E198A mutation. In the present study, a loop-mediated isothermal amplification (LAMP) assay was established for rapid detection of the E198A SNP in H. contortus. The results showed that optimization of LAMP reaction reagents and conditions could achieve this. The resulting amplicons were visualized by adding hydroxynaphthol blue dye (HNB) prior to amplification. The color of LAMP products amplified without DNA or from DNA from worms with the E198A homozygous susceptible genotype was still violet, but the products with DNA from worms with the E198A heterozygous genotype or the E198A resistant homozygous genotype changed to sky blue. The specificity of this method was further verified by sequencing, which confirmed the successful LAMP detection of the E198A mutation with high specificity. In conclusion, the developed LAMP method has high specificity and good reproducibility for screening the E198A SNP of isotype-1 ß-tubulin gene of H. contortus of field samples without using sophisticated equipment, providing useful technique for the rapid detection and thus prevention and control of benzimidazole resistant H. contortus infections.

A TGF-ß type II receptor that associates with developmental transition in Haemonchus contortus in vitro.

BACKGROUND: The TGF-ß signalling pathway plays a key role in regulating dauer formation in the free-living nematode Caenorhabditis elegans, and previous work has shown that TGF-ß receptors are involved in parasitic nematodes. Here, we explored the structure and function of a TGF-ß type II receptor homologue in the TGF-ß signalling pathway in Haemonchus contortus, a highly pathogenic, haematophagous parasitic nematode. METHODOLOGY/PRINCIPAL FINDINGS: Amino acid sequence and phylogenetic analyses revealed that the protein, called Hc-TGFBR2 (encoded by the gene Hc-tgfbr2), is a member of TGF-ß type II receptor family and contains conserved functional domains, both in the extracellular region containing cysteine residues that form a characteristic feature (CXCX4C) of TGF-ß type II receptor and in the intracellular regions containing a serine/threonine kinase domain. The Hc-tgfbr2 gene was transcribed in all key developmental stages of H. contortus, with particularly high levels in the infective third-stage larvae (L3s) and male adults. Immunohistochemical results revealed that Hc-TGFBR2 was expressed in the intestine, ovary and eggs within the uterus of female adults, and also in the testes of male adults of H. contortus. Double-stranded RNA interference (RNAi) in this nematode by soaking induced a marked decrease in transcription of Hc-tgfbr2 and in development from the exsheathed L3 to the fourth-stage larva (L4) in vitro. CONCLUSIONS/SIGNIFICANCE: These results indicate that Hc-TGFBR2 plays an important role in governing developmental processes in H. contortus via the TGF-ß signalling pathway, particularly in the transition from the free-living to the parasitic stages.

Characterization of the Complete Mitochondrial Genome of Ostertagia trifurcata of Small Ruminants and its Phylogenetic Associations for the Trichostrongyloidea Superfamily.

The complete mitochondrial (mt) genome of Ostertagia trifurcata, a parasitic nematode of small ruminants, has been sequenced and its phylogenetic relationship with selected members from the superfamily Trichostrongyloidea was investigated on the basis of deduced datasets of mt amino acid sequences. The entire mt genome of Ostertagia trifurcata is circular and 14,151 bp in length. It consists of a total of 36 genes comprising 12 genes coding for proteins (PCGs), 2 genes for ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes and 2 non-coding regions, since all genes are transcribed in the same direction. The phylogenetic analysis based on the concatenated datasets of predicted amino acid sequences of the 12 protein coding genes supported monophylies of the Haemonchidae, Dictyocaulidae and Molineidae families, but rejected monophylies of the Trichostrongylidae family. The complete characterization and provision of the mtDNA sequence of Ostertagia trifurcata provides novel genetic markers for molecular epidemiological investigations, systematics, diagnostics and population genetics of Ostertagia trifurcata and its correspondents.

Characterization of the Complete Mitochondrial Genome of a Whipworm Trichuris skrjabini (Nematoda: Trichuridae).

The complete mitochondrial (mt) genome of Trichuris skrjabini has been determined in the current study and subsequently compared with closely related species by phylogenetic analysis based on concatenated datasets of mt amino acid sequences. The whole mt genome of T. skrjabini is circular and 14,011 bp in length. It consists of a total of 37 genes including 13 protein coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNAs) genes, and two non-coding regions. The gene arrangement and contents were consistent with other members of the Trichuridae family including Trichuris suis, Trichuris trichiura, Trichuris ovis, and Trichuris discolor. Phylogenetic analysis based on concatenated datasets of amino acids of the 12 PCGs predicted the distinctiveness of Trichuris skrjabini as compared to other members of the Trichuridae family. Overall, our study supports the hypothesis that T. skrjabini is a distinct species. The provision of molecular data of whole mt genome of T. skrjabini delivers novel genetic markers for future studies of diagnostics, systematics, population genetics, and molecular epidemiology of T. skrjabini.