Deep small RNA sequencing from the nematode Ascaris reveals conservation, functional diversification, and novel developmental profiles

End resection and telomere healing of DNA double-strand breaks during nematode programmed DNA elimination

Most DNA double-strand breaks (DSBs) are harmful to genome integrity. However, some forms of DSBs are essential to biological processes, such as meiotic recombination and V(D)J recombination. DSBs are also required for programmed DNA elimination (PDE) in ciliates and nematodes. In nematodes, the DSBs are healed with telomere addition. While telomere addition sites have been well characterized, little is known regarding the DSBs that fragment nematode chromosomes. Here, we used embryos from the human and pig parasitic nematode Ascaris to characterize the DSBs. Using END-seq, we demonstrate that DSBs are introduced before mitosis, followed by extensive end resection. The resection profile is unique for each break site, and the resection generates 3'-overhangs before the addition of neotelomeres. Interestingly, telomere healing occurs much more frequently on retained DSB ends than on eliminated ends. This biased repair of the DSB ends may be due to the sequestration of the eliminated DNA into micronuclei, preventing neotelomere formation at their ends. Additional DNA breaks occur within the eliminated DNA in both Ascaris and Parascaris, ensuring chromosomal breakage and providing a fail-safe mechanism for PDE. Overall, our data indicate that telomere healing of DSBs is specific to the break sites responsible for nematode PDE.

Diversity of Picorna-Like Viruses in the Teltow Canal, Berlin, Germany

The viromes of freshwater bodies are underexplored. The Picornavirales order, with 371 acknowledged species, is one of the most expansive and diverse groups of eukaryotic RNA viruses. In this study, we add 513 picorna-like viruses to the assemblage of more than 2000 unassigned picorna-like viruses. Our set of the aquatic Picornavirales virome of the Teltow Canal in Berlin, Germany, consists of 239 complete and 274 partial genomes. This urban freshwater body is characterized by the predominance of marna-like viruses (30.8%) and dicistro-like viruses (19.1%), whereas picornaviruses, iflaviruses, solinvi-like viruses, polycipi-like viruses, and nora-like viruses are considerably less prevalent. Caliciviruses and secoviruses were absent in our sample. Although presenting characteristic domains of Picornavirales, more than 100 viruses (20.8%) could not be assigned to any of the 9 Picornavirales families. Thirty-three viruses of the Marnaviridae-mostly locarna-like viruses-exhibit a monocistronic genome layout. Besides a wealth of novel virus sequences, viruses with peculiar features are reported. Among these is a clade of untypeable marna-like viruses with dicistronic genomes, but with the capsid protein-encoding open reading frame located at the 5' part of their RNA. A virus with a similar genome layout but clustering with dicistroviruses was also observed. We further detected monocistronic viruses with a polymerase gene related to aparaviruses. The detection of Aichi virus and five novel posa-like viruses indicates a slight burden in municipal wastewater.

End resection and telomere healing of DNA double-strand breaks during nematode programmed DNA elimination

Most DNA double-strand breaks (DSBs) are harmful to genome integrity. However, some forms of DSBs are essential to biological processes, such as meiotic recombination and V(D)J recombination. DSBs are also required for programmed DNA elimination (PDE) in ciliates and nematodes. In nematodes, the DSBs are healed with telomere addition. While telomere addition sites have been well-characterized, little is known regarding the DSBs that fragment nematode chromosomes. Here, we used embryos from the nematode Ascaris to study the timing of PDE breaks and examine the DSBs and their end processing. Using END-seq, we characterize the DSB ends and demonstrate that DNA breaks are introduced before mitosis, followed by extensive end resection. The resection profile is unique for each break site, and the resection generates 3' overhangs before the addition of telomeres. Interestingly, telomere healing occurs much more frequently on retained DSB ends than on eliminated ends. This biased repair of the DSB ends in Ascaris may be due to the sequestration of the eliminated DNA into micronuclei, preventing their ends from telomere healing. Additional DNA breaks occur within the eliminated DNA in both Ascaris and Parascaris, ensuring chromosomal breakage and providing a fail-safe mechanism for nematode PDE.

Getting around the roundworms: Identifying knowledge gaps and research priorities for the ascarids

The ascarids are a large group of parasitic nematodes that infect a wide range of animal species. In humans, they cause neglected diseases of poverty; many animal parasites also cause zoonotic infections in people. Control measures include hygiene and anthelmintic treatments, but they are not always appropriate or effective and this creates a continuing need to search for better ways to reduce the human, welfare and economic costs of these infections. To this end, Le Studium Institute of Advanced Studies organized a two-day conference to identify major gaps in our understanding of ascarid parasites with a view to setting research priorities that would allow for improved control. The participants identified several key areas for future focus, comprising of advances in genomic analysis and the use of model organisms, especially Caenorhabditis elegans, a more thorough appreciation of the complexity of host-parasite (and parasite-parasite) communications, a search for novel anthelmintic drugs and the development of effective vaccines. The participants agreed to try and maintain informal links in the future that could form the basis for collaborative projects, and to co-operate to organize future meetings and workshops to promote ascarid research.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Programmed DNA elimination in the parasitic nematode Ascaris

In most organisms, the whole genome is maintained throughout the life span. However, exceptions occur in some species where the genome is reduced during development through a process known as programmed DNA elimination (PDE). In the human and pig parasite Ascaris, PDE occurs during the 4 to 16 cell stages of embryogenesis, when germline chromosomes are fragmented and specific DNA sequences are reproducibly lost in all somatic cells. PDE was identified in Ascaris over 120 years ago, but little was known about its molecular details until recently. Genome sequencing revealed that approximately 1,000 germline-expressed genes are eliminated in Ascaris, suggesting PDE is a gene silencing mechanism. All germline chromosome ends are removed and remodeled during PDE. In addition, PDE increases the number of chromosomes in the somatic genome by splitting many germline chromosomes. Comparative genomics indicates that these germline chromosomes arose from fusion events. PDE separates these chromosomes at the fusion sites. These observations indicate that PDE plays a role in chromosome karyotype and evolution. Furthermore, comparative analysis of PDE in other parasitic and free-living nematodes illustrates conserved features of PDE, suggesting it has important biological significance. We summarize what is known about PDE in Ascaris and its relatives. We also discuss other potential functions, mechanisms, and the evolution of PDE in these parasites of humans and animals of veterinary importance.

Intracellular Compartmentalization: A Key Determinant of MicroRNA Functions

Being an integral part of the eukaryotic transcriptome, miRNAs are regarded as vital regulators of diverse developmental and physiological processes. Clearly, miRNA activity is kept in check by various regulatory mechanisms that control their biogenesis and decay pathways. With the increasing technical depth of RNA profiling technologies, novel insights have unravelled the spatial diversity exhibited by miRNAs inside a cell. Compartmentalization of miRNAs adds complexity to the regulatory circuits of miRNA expression, thereby providing superior control over the miRNA function. This review provides a bird's eye view of miRNAs expressed in different subcellular locations, thus affecting the gene regulatory pathways therein. Occurrence of miRNAs in diverse intracellular locales also reveals various unconventional roles played by miRNAs in different cellular organelles and expands the scope of miRNA functions beyond their traditionally known repressive activities.

microRNAs: Critical Players during Helminth Infections

microRNAs (miRNAs) are a group of small non-coding RNAs that regulate gene expression post-transcriptionally through their interaction with the 3' untranslated regions (3' UTR) of target mRNAs, affecting their stability and/or translation. Therefore, miRNAs regulate biological processes such as signal transduction, cell death, autophagy, metabolism, development, cellular proliferation, and differentiation. Dysregulated expression of microRNAs is associated with infectious diseases, where miRNAs modulate important aspects of the parasite-host interaction. Helminths are parasitic worms that cause various neglected tropical diseases affecting millions worldwide. These parasites have sophisticated mechanisms that give them a surprising immunomodulatory capacity favoring parasite persistence and establishment of infection. In this review, we analyze miRNAs in infections caused by helminths, emphasizing their role in immune regulation and its implication in diagnosis, prognosis, and the development of therapeutic strategies.

Transcriptome profiling of male and female Ascaris lumbricoides reproductive tissues

BACKGROUND

Ascaris lumbricoides causes human ascariasis, the most prevalent helminth disease, infecting approximately 1 billion individuals globally. In 2019 the global disease burden was estimated to be 754,000 DALYs and resulted in 2090 deaths. In the absence of a vaccination strategy, treatment of ascariasis has relied on anthelminthic chemotherapy, but drug resistance is a concern. The propensity for reinfection is also a major challenge to disease control; female worms lay up to 200,000 eggs daily, which contaminate surrounding environments and remain viable for years, resulting in high transmission rates. Understanding the molecular mechanisms of reproductive processes, including control of egg production, spermatogenesis, oogenesis and embryogenesis, will drive the development of new drugs and/or vaccine targets for future ascariasis control.

METHODS

Transcriptome profiles of discrete reproductive and somatic tissue samples were generated from adult male and female worms using Illumina HiSeq with 2 × 150 bp paired-end sequencing. Male tissues included: testis germinal zone, testis part of vas deferens, seminal vesicle and somatic tissue. Female tissues included: ovary germinal zone, ovary part of the oviduct, uterus and somatic tissue. Differentially expressed genes (DEGs) were identified from the fragments per kilobases per million reads (FPKM) profiles. Hierarchical analysis was performed to identify tissue-specific genes. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were employed to identify significant terms and pathways for the DEGs.

RESULTS

DEGs involved in protein phosphorylation and adhesion molecules were indicated to play a crucial role in spermatogenesis and fertilization, respectively. Those genes associated with the G-protein-coupled receptor (GPCR) signaling pathway and small GTPase-mediated signal transduction pathway play an essential role in cytoskeleton organization during oogenesis. Additionally, DEGs associated with the SMA genes and TGF-β signaling pathway are crucial in adult female embryogenesis. Some genes associated with particular biological processes and pathways that were identified in this study have been linked to defects in germline development, embryogenesis and reproductive behavior. In the enriched KEGG pathway analysis, Hippo signaling, oxytocin signaling and tight junction pathways were identified to play a role in Ascaris male and female reproductive systems.

CONCLUSIONS

This study has provided comprehensive transcriptome profiles of discrete A. lumbricoides reproductive tissue samples, revealing the molecular basis of these functionally important tissues. The data generated from this study will provide fundamental knowledge on the reproductive biology of Ascaris and will inform future target identification for anti-ascariasis drugs and/or vaccines.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

Identification of key interactions of benzimidazole resistance-associated amino acid mutations in Ascaris β-tubulins by molecular docking simulations

Ascaris species are soil-transmitted helminths that infect humans and livestock mainly in low and middle-income countries. Benzimidazole (BZ) class drugs have predominated for many years in the treatment of Ascaris infections, but persistent use of BZs has already led to widespread resistance in other nematodes, and treatment failure is emerging for Ascaris. Benzimidazoles act by binding to β-tubulin proteins and destabilising microtubules. Three mutations in the β-tubulin protein family are associated with BZ resistance. Seven shared β-tubulin isotypes were identified in Ascaris lumbricoides and A. suum genomes. Benzimidazoles were predicted to bind to all β-tubulin isotypes using in silico docking, demonstrating that the selectivity of BZs to interact with one or two β-tubulin isotypes is likely the result of isotype expression levels affecting the frequency of interaction. Ascaris β-tubulin isotype A clusters with helminth β-tubulins previously shown to interact with BZ. Molecular dynamics simulations using β-tubulin isotype A highlighted the key role of amino acid E198 in BZ-β-tubulin interactions. Simulations indicated that mutations at amino acids E198A and F200Y alter binding of BZ, whereas there was no obvious effect of the F167Y mutation. In conclusion, the key interactions vital for BZ binding with β-tubulins have been identified and show how mutations can lead to resistance in nematodes.

Changes in the expression of miR-103a and miR-21: a functional diagnosis of toxocariasis in rats

Introduction. Toxocariasis is a zoonotic parasitic disease caused by migrating nematode worms, Toxocara species larvae, within tissues. MicroRNAs (miRNAs) are small RNA molecules that regulate gene expression at a post-transcriptional level.Hypothesis/Gap Statement. miRNA-based diagnostic biomarkers for toxocariasis are emerging, but there is limited information about the role of many miRNAs and a more detailed diagnostic evaluation of miRNA expression patterns is needed to understand their immunobiological function.Aim. We investigated the expression levels of circulating miRNA 21 and miRNA 103a as potential biomarkers for the prediction and diagnosis of toxocariasis in Wistar rats infected with Toxocara canis.Methodology. Thirty Wistar rats were inoculated orally with 2500 T. canis embryonated eggs via gavage. Serum samples were collected from infected animals and were tested against T. canis antigens for 60 days post-infection. The plasma samples were isolated for quantitative real-time PCR (qPCR) assays and qPCR was used to assess transcription levels of miRNA 21 and miRNA 103a.Results. The prevalence of anti-Toxocara IgG was detected in 7/30 (23.3 %) infected rats. Molecular analysis of miRNAs 21 and 103a showed that expression levels of miRNAs in both groups of Toxocara-positive and negative samples were the same without significant association. The ratio of housekeeping gene expression (U6) to gene expression of miRNAs 21 and 103a indicated the rate of change (1/1.38 ≈ 0.75 and 1/0.751 ≈ 1.3, respectively).Conclusion. Our study revealed that miRNAs 21 and 103a might play fundamental roles as biomarkers and diagnostic tools for toxocariasis. However, the changes in expression of these miRNAs were not adequate to be used as biomarkers in diagnosis.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Small RNA pathways in the nematode Ascaris in the absence of piRNAs

Small RNA pathways play key and diverse regulatory roles in C. elegans, but our understanding of their conservation and contributions in other nematodes is limited. We analyzed small RNA pathways in the divergent parasitic nematode Ascaris. Ascaris has ten Argonautes with five worm-specific Argonautes (WAGOs) that associate with secondary 5’-triphosphate 22-24G-RNAs. These small RNAs target repetitive sequences or mature mRNAs and are similar to the C. elegans mutator, nuclear, and CSR-1 small RNA pathways. Even in the absence of a piRNA pathway, Ascaris CSR-1 may still function to “license” as well as fine-tune or repress gene expression. Ascaris ALG-4 and its associated 26G-RNAs target and likely repress specific mRNAs during testis meiosis. Ascaris WAGO small RNAs demonstrate target plasticity changing their targets between repeats and mRNAs during development. We provide a unique and comprehensive view of mRNA and small RNA expression throughout spermatogenesis. Overall, our study illustrates the conservation, divergence, dynamics, and flexibility of small RNA pathways in nematodes.

The parasitic nematode Ascaris lacks piRNAs. Here the authors compare Argonaute proteins and small RNAs from C. elegans and Ascaris, expanding our understanding of the conservation, divergence, and flexibility of Argonautes and small RNA pathways in nematodes.

Delete and survive: strategies of programmed genetic material elimination in eukaryotes

Genome stability is a crucial feature of eukaryotic organisms because its alteration drastically affects the normal development and survival of cells and the organism as a whole. Nevertheless, some organisms can selectively eliminate part of their genomes from certain cell types during specific stages of ontogenesis. This review aims to describe the phenomenon of programmed DNA elimination, which includes chromatin diminution (together with programmed genome rearrangement or DNA rearrangements), B and sex chromosome elimination, paternal genome elimination, parasitically induced genome elimination, and genome elimination in animal and plant hybrids. During programmed DNA elimination, individual chromosomal fragments, whole chromosomes, and even entire parental genomes can be selectively removed. Programmed DNA elimination occurs independently in different organisms, ranging from ciliate protozoa to mammals. Depending on the sequences destined for exclusion, programmed DNA elimination may serve as a radical mechanism of dosage compensation and inactivation of unnecessary or dangerous genetic entities. In hybrids, genome elimination results from competition between parental genomes. Despite the different consequences of DNA elimination, all genetic material destined for elimination must be first recognised, epigenetically marked, separated, and then removed and degraded.

MicroRNAs in helminth parasites: a systematic review

BACKGROUND

MicroRNAs (miRNAs) are about 22-nucleotide, small, non-coding RNAs that control gene expression post-transcriptionally. Helminth parasites usually express a unique repertoire of genes, including miRNAs, across different developmental stages with subtle regulatory mechanisms.

OBJECTIVE

There is a necessity to investigate the involvement of miRNAs in the development of parasites, host-parasite interaction, immune evasion and their abilities to govern infection in hosts. miRNAs present in helminth parasites have been summarized in the current systematic review (SR).

METHODS

Electronic databases, including PubMed, Scopus, ProQuest, Embase, and Google Scholar search engine, were searched to identify helminth miRNA studies published from February 1993 till December 2019. Only the published articles in English were included in the study.

RESULTS

A total of 1769 articles were preliminarily recorded. Following the strict inclusion and exclusion criteria, 105 studies were included in this SR. Most of these studies focused on the identification of miRNAs in helminth parasites and/or probing of differentially expressed host miRNA profiles in specific relevant tissues, while 12 studies aimed to detect parasite-derived miRNAs in host circulating system and 15 studies characterized extracellular vesicles (EV)-derived miRNAs secreted by parasites.

CONCLUSION

In the current SR, information regarding all miRNAs expressed in helminth parasites has been comprehensively provided and the utility of helminth parasites-derived miRNAs in diagnosis and control of parasitic infections has been discussed. Furthermore, functional studies on helminth-derived miRNAs have also been presented.

Programmed DNA elimination: silencing genes and repetitive sequences in somatic cells

In a multicellular organism, the genomes of all cells are in general the same. Programmed DNA elimination is a notable exception to this genome constancy rule. DNA elimination removes genes and repetitive elements in the germline genome to form a reduced somatic genome in various organisms. The process of DNA elimination within an organism is highly accurate and reproducible; it typically occurs during early embryogenesis, coincident with germline-soma differentiation. DNA elimination provides a mechanism to silence selected genes and repeats in somatic cells. Recent studies in nematodes suggest that DNA elimination removes all chromosome ends, resolves sex chromosome fusions, and may also promote the birth of novel genes. Programmed DNA elimination processes are diverse among species, suggesting DNA elimination likely has evolved multiple times in different taxa. The growing list of organisms that undergo DNA elimination indicates that DNA elimination may be more widespread than previously appreciated. These various organisms will serve as complementary and comparative models to study the function, mechanism, and evolution of programmed DNA elimination in metazoans.

Characterization of the β-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance

Background

The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding β-tubulin proteins, the principal targets of BZ drugs.

Methodology and principal findings

First, the β-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different β-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations.

Significance

This study demonstrated the presence of at least seven β-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the β-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of β-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.

Benzimidazole (BZ) drugs remain the standard of treatment in large-scale deworming programs that aim to control the morbidity caused by intestinal worms. As these deworming programs are expanding world-wide, there is an increasing risk of worms becoming resistant to BZ drugs, highlighting the necessity for tools to detect gene mutations associated with drug resistance. However, the development of such tools is impeded by a lack of insights into the genes that are coding for β-tubulin proteins, which are the principal targets of BZ drugs. The aim of this study was to comprehensively characterize these genes in the worm species Ascaris lumbricoides and Ascaris suum. The findings highlight that these species have at least seven β-tubulin genes. Only two genes are highly expressed throughout the different life stages of the worm, and hence are more likely to be involved in the development of BZ resistance. No mutations that have previously been associated with BZ resistance in other intestinal worms were found. This study provides a baseline towards more efficient and accurate monitoring of drug resistance in large-scale deworming programs.

Homology modelling, molecular dynamics simulation and docking evaluation of β-tubulin of Schistosoma mansoni

Schistosomiasis is one of the neglected diseases causing considerable morbidity and mortality throughout the world. Microtubules with its main component, tubulin play a vital role in helminthes including schistosomes. Benzimidazoles represent potential drug candidates by binding β-tubulin. The study aimed to generate a homology model for the β-tubulin of S. mansoni using the crystal structure of O visaries (Sheep) β-tubulin (PDB ID: 3N2G D) as a template, then different β-tubulin models were generated and two previously reported benzimidazole derivatives (NBTP-F and NBTP-OH) were docked to the generated models, the binding results indicated that both S. mansoni, S. haematobium were susceptible to the two NBTP derivatives. Additionally, three mutated versions of S. mansoni β-tubulin wild-type were generated and the mutation (F185Y) seems to slightly enhance the ligand binding. Dynamics simulation experiments showed S. haematobium β-tubulin is highly susceptible to the tested compounds; similar to S. mansoni, moreover, mutated models of S. mansoni β-tubulin altered its NBTPs susceptibility. Moreover, additional seven new benzimidazole derivatives were synthesized and tested by molecular docking on the generated model binding site of S. mansoni β-tubulin and were found to have good interaction inside the pocket.

First Evidence of Function for Schistosoma japonicumriok-1 and RIOK-1

Protein kinases are known as key molecules that regulate many biological processes in animals. The right open reading frame protein kinase (riok) genes are known to be essential regulators in model organisms such as the free-living nematode Caenorhabditis elegans. However, very little is known about their function in parasitic trematodes (flukes). In the present study, we characterized the riok-1 gene (Sj-riok-1) and the inferred protein (Sj-RIOK-1) in the parasitic blood fluke, Schistosoma japonicum. We gained a first insight into function of this gene/protein through double-stranded RNA interference (RNAi) and chemical inhibition. RNAi significantly reduced Sj-riok-1 transcription in both female and male worms compared with untreated control worms, and subtle morphological alterations were detected in the ovaries of female worms. Chemical knockdown of Sj-RIOK-1 with toyocamycin (a specific RIOK-1 inhibitor/probe) caused a substantial reduction in worm viability and a major accumulation of mature oocytes in the seminal receptacle (female worms), and of spermatozoa in the sperm vesicle (male worms). These phenotypic alterations indicate that the function of Sj-riok-1 is linked to developmental and/or reproductive processes in S. japonicum.

Parasite-derived microRNAs as a diagnostic biomarker: potential roles, characteristics, and limitations

MicroRNAs (miRNAs), a subclass of small regulatory RNAs that present from ancient unicellular protozoans to parasitic helminths and parasitic arthropods. MiRNAs' mode of action has attracted wide attention as a result of their unique functional importance. MiRNAs play a role in diverse physiological and pathological processes ranging from organ development, immune function to apoptosis and cancer at the post-transcription gene expression. Thus, miRNAs are known to be targets for clinical treatment and therapy. The discovery of the high stability of circulating miRNA in various types of host body fluids, such as whole blood, serum, plasma, saliva, and urine has increased great interest among researchers in the potential of circulating miRNA as a prognosis/diagnosis of infectious. Some circulating miRNAs biomarkers advanced to clinical applications related to human diseases. However, this idea starts to come only in the fields of infectious disease. The goal of this review is to enhance the current understanding of these molecules and their applicability in the field of medicine. A detailed review of the available literature consulting tools performed in online repositories such as NCBI, PubMed, Medline, ScienceDirect, and UpToDate. This review summarizes an overview of preclinical studies using circulating miRNAs biomarkers against infectious diseases affecting humans. The use of miRNA as a safe and potential tool is encouraging news, considering that until now, guidelines for the use of miRNA in clinical practice are still lacking.

Genomics of the Parasitic Nematode Ascaris and Its Relatives

Nematodes of the genus Ascaris are important parasites of humans and swine, and the phylogenetically related genera (Parascaris, Toxocara, and Baylisascaris) infect mammals of veterinary interest. Over the last decade, considerable genomic resources have been established for Ascaris, including complete germline and somatic genomes, comprehensive mRNA and small RNA transcriptomes, as well as genome-wide histone and chromatin data. These datasets provide a major resource for studies on the basic biology of these parasites and the host-parasite relationship. Ascaris and its relatives undergo programmed DNA elimination, a highly regulated process where chromosomes are fragmented and portions of the genome are lost in embryonic cells destined to adopt a somatic fate, whereas the genome remains intact in germ cells. Unlike many model organisms, Ascaris transcription drives early development beginning prior to pronuclear fusion. Studies on Ascaris demonstrated a complex small RNA network even in the absence of a piRNA pathway. Comparative genomics of these ascarids has provided perspectives on nematode sex chromosome evolution, programmed DNA elimination, and host-parasite coevolution. The genomic resources enable comparison of proteins across diverse species, revealing many new potential drug targets that could be used to control these parasitic nematodes.

Multiple clades of Husavirus in South America revealed by next generation sequencing

Husavirus (HuV) is an unclassified virus of the order Picornavirales that has already been identified worldwide in various locations. The genetic, epidemiological, and pathogenic characteristics are, however, little understood. In children with acute gastroenteritis, this study used next-generation sequencing to recognize unknown sources of viruses. In particular, 251 fecal samples obtained from individuals were sequenced in southern, northeastern, and northern Brazil. all samples were also analyzed using culture methods and parasitological tests to classify other enteric pathogens such as bacteria, parasites, and viruses. 1.9% of the samples tested positive for HuV, for a total of 5 positive children, with a mean age of 2 year, with three males and two females. Detailed molecular characterization of full genomes showed that Brazilian HuVs’ nucleotide divergence is less than 11%. The genetic gap between Brazilian sequences and the closest HuV reported previously, on the other hand, is 18%. The study showed that Brazilian sequences are closely related to the HuV defined in Viet Nam in 2013, further characterization based on phylogenetics. At least two divergent clades of HuV in South America were also seen in the phylogenetic study.

Expression profiling of Echinococcus multilocularis miRNAs throughout metacestode development in vitro

The neglected zoonotic disease alveolar echinococcosis (AE) is caused by the metacestode stage of the tapeworm parasite Echinococcus multilocularis. MicroRNAs (miRNAs) are small non-coding RNAs with a major role in regulating gene expression in key biological processes. We analyzed the expression profile of E. multilocularis miRNAs throughout metacestode development in vitro, determined the spatial expression of miR-71 in metacestodes cultured in vitro and predicted miRNA targets. Small cDNA libraries from different samples of E. multilocularis were sequenced. We confirmed the expression of 37 miRNAs in E. multilocularis being some of them absent in the host, such as miR-71. We found a few miRNAs highly expressed in all life cycle stages and conditions analyzed, whereas most miRNAs showed very low expression. The most expressed miRNAs were miR-71, miR-9, let-7, miR-10, miR-4989 and miR-1. The high expression of these miRNAs was conserved in other tapeworms, suggesting essential roles in development, survival, or host-parasite interaction. We found highly regulated miRNAs during the different transitions or cultured conditions analyzed, which might suggest a role in the regulation of developmental timing, host-parasite interaction, and/or in maintaining the unique developmental features of each developmental stage or condition. We determined that miR-71 is expressed in germinative cells and in other cell types of the germinal layer in E. multilocularis metacestodes cultured in vitro. MiRNA target prediction of the most highly expressed miRNAs and in silico functional analysis suggested conserved and essential roles for these miRNAs in parasite biology. We found relevant targets potentially involved in development, cell growth and death, lifespan regulation, transcription, signal transduction and cell motility. The evolutionary conservation and expression analyses of E. multilocularis miRNAs throughout metacestode development along with the in silico functional analyses of their predicted targets might help to identify selective therapeutic targets for treatment and control of AE.

Alveolar echinococcosis (AE) is a zoonotic disease caused by the metacestode stage of the helminth parasite Echinococcus multilocularis. Current treatment requires surgery and/or prolonged drug therapy. Thus, novel strategies for the treatment of AE are needed. MicroRNAs (miRNAs), a class of small ~22-nucleotide (nt) non-coding RNAs with a major role in regulating gene expression, have been suggested as potential therapeutic targets for treatment and control of helminth parasite infections. In this work, we analyzed the expression profile of E. multilocularis miRNAs throughout metacestode development in vitro. We predicted functional roles for highly expressed miRNAs and found that they could be involved in essential roles for survival and development in the host. We determined that E. multilocularis miR-71, a highly expressed miRNA that is absent in the human host, is expressed in germinative cells and in other cell types of the germinal layer in E. multilocularis metacestodes cultured in vitro. Germinative cells are a relevant cell type to target for anti-echinococcosis drug development. MiRNAs that are absent in the human host, involved in essential functions, highly expressed and/or expressed in germinative cells in E. multilocularis metacestodes may represent selective therapeutic targets for treatment and control of AE.

Helminth-derived cystatins: the immunomodulatory properties of an Ascaris lumbricoides cystatin

Helminth infections such as ascariasis elicit a type 2 immune response resembling that involved in allergic inflammation, but differing to allergy, they are also accompanied with strong immunomodulation. This has stimulated an increasing number of investigations, not only to better understand the mechanisms of allergy and helminth immunity but to find parasite-derived anti-inflammatory products that could improve the current treatments of chronic non-communicable inflammatory diseases such as asthma. A great number of helminth-derived immunomodulators have been discovered and some of them extensively analysed, showing their potential use as anti-inflammatory drugs in clinical settings. Since Ascaris lumbricoides is one of the most successful parasites, several groups have focused on the immunomodulatory properties of this helminth. As a result, several excretory/secretory components and purified molecules have been analysed, revealing interesting anti-inflammatory activities potentially useful as therapeutic tools. One of these molecules is A. lumbricoides cystatin, whose genomic, cellular, molecular, and immunomodulatory properties are described in this review.

Genome Analysis of Programmed DNA Elimination in Parasitic Nematodes

Two distinct groups of parasitic nematodes use programmed DNA elimination to silence germline-expressed genes in the somatic cells (ascarids) or for sex determination (Strongyloides spp.). In the ascarids, DNA is lost only in pre-somatic cells during early embryogenesis, leading to a reduced somatic genome compared to the intact germ cell genome. Comparative genome analysis has provided information on the retained vs. eliminated sequences, DNA breaks, a full chromosome view on DNA elimination, and the evolutionary conservation of DNA elimination among ascarids. These studies have revealed novel insights into the functions and mechanisms of DNA elimination and provided a reference for in-depth molecular analysis of DNA elimination. Here, I describe the genomics methods we used to study programmed DNA elimination, focusing on the parasitic nematode Ascaris.

Expression of Mir-21 and Mir-103a in Toxocara canis: Potential for Diagnosis of Human Toxocariasis

Background:

Toxocariasis is one of the most neglected zoonotic diseases, predominantly caused by Toxocara canis. We aimed to evaluate the expression of microRNAs 21 and 103a in seropositive individuals for human toxocariasis as diagnostic biomarkers.

Methods:

This study was conducted on 324 individuals for ELISA test on toxocariasis in Tehran and Karaj, Iran 2019. Then positive samples for anti-Toxocara IgG were obtained to quantitative Real-time PCR (qRT-PCR) assays to investigate the transcriptional profiles of miRNAs predicted to be involved in developmental and reproductive processes. qPCR was employed to assess levels of transcription for miRNAs of 103a and 21 in plasma samples.

Results:

After the experiments, the results were evaluated by REST software, Livak formula and quantitative t-test. The analyzes performed on human samples showed that in the case group compared to the control group, only in Tc-miR-21 gene, a 0.3-fold increase in expression was obtained with REST software (Fold change ≤ 1.5, P>0.05), which was statistically significant by t-test (P<0.05).

Conclusion:

To our knowledge, this is the first study to evaluate miR-21 and miR-103a in toxocariasis, which shed light on the fundamental role of it as a biomarker and diagnostic tool. However, due to the changes in expression of these miRNAs were not vast to be used as biomarkers in diagnosis. Despite of that the changes in the expression of these miRNAs were not vast but they could serve as novel promising biomarkers for diagnosis of toxocariasis.

Identification of oligo-adenylated small RNAs in the parasite Entamoeba and a potential role for small RNA control

Background

The RNA interference (RNAi) pathway is a gene regulation mechanism that utilizes small RNA (sRNA) and Argonaute (Ago) proteins to silence target genes. Our previous work identified a functional RNAi pathway in the protozoan parasite Entamoeba histolytica, including abundant 27 nt antisense sRNA populations which associate with EhAgo2–2 protein. However, there is lack of understanding about the sRNAs that are bound to two other EhAgos (EhAgo2–1 and 2–3), and the mechanism of sRNA regulation itself is unclear in this parasite. Therefore, identification of the entire pool of sRNA species and their sub-populations that associate with each individual EhAgo protein would be a major step forward.

Results

In the present study, we sequenced sRNA libraries from both total RNAs and EhAgo bound RNAs. We identified a new population of 31 nt sRNAs that results from the addition of a non-templated 3–4 adenosine nucleotides at the 3′-end of the 27 nt sRNAs, indicating a non-templated RNA-tailing event in the parasite. The relative abundance of these two sRNA populations is linked to the efficacy of gene silencing for the target gene when parasites are transfected with an RNAi-trigger construct, indicating that non-templated sRNA-tailing likely play a role in sRNA regulation in this parasite. We found that both sRNA populations (27 nt and 31 nt) are present in the related parasite Entamoeba invadens, and are unchanged during the development. In sequencing the sRNAs associating with the three EhAgo proteins, we observed that despite distinct cellular localization, all three EhAgo sRNA libraries contain 27 nt sRNAs with 5′-polyphosphate (5′-polyP) structure and share a largely overlapping sRNA repertoire. In addition, our data showed that a fraction of 31 nt sRNAs associate with EhAgo2–2 but not with its mutant protein (C-terminal deletion), nor other two EhAgos, indicating a specific EhAgo site may be required for sRNA modification process in the parasite.

Conclusion

We identified a new population of sRNA with non-templated oligo-adenylation modification, which is the first such observation amongst single celled protozoan parasites. Our sRNA sequencing libraries provide the first comprehensive sRNA dataset for all three Entamoeba Ago proteins, which can serve as a useful database for the amoeba community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07275-6.

Improving Transgenesis Efficiency and CRISPR-Associated Tools Through Codon Optimization and Native Intron Addition in Pristionchus Nematodes

A lack of appropriate molecular tools is one obstacle that prevents in-depth mechanistic studies in many organisms. Transgenesis, clustered regularly interspaced short palindromic repeats (CRISPR)-associated engineering, and related tools are fundamental in the modern life sciences, but their applications are still limited to a few model organisms. In the phylum Nematoda, transgenesis can only be performed in a handful of species other than Caenorhabditis elegans, and additionally, other species suffer from significantly lower transgenesis efficiencies. We hypothesized that this may in part be due to incompatibilities of transgenes in the recipient organisms. Therefore, we investigated the genomic features of 10 nematode species from three of the major clades representing all different lifestyles. We found that these species show drastically different codon usage bias and intron composition. With these findings, we used the species Pristionchus pacificus as a proof of concept for codon optimization and native intron addition. Indeed, we were able to significantly improve transgenesis efficiency, a principle that may be usable in other nematode species. In addition, with the improved transgenes, we developed a fluorescent co-injection marker in P. pacificus for the detection of CRISPR-edited individuals, which helps considerably to reduce associated time and costs.

Molecular evidence of hybridization between pig and human Ascaris indicates an interbred species complex infecting humans

Human ascariasis is a major neglected tropical disease caused by the nematode Ascaris lumbricoides. We report a 296 megabase (Mb) reference-quality genome comprised of 17,902 protein-coding genes derived from a single, representative Ascaris worm. An additional 68 worms were collected from 60 human hosts in Kenyan villages where pig husbandry is rare. Notably, the majority of these worms (63/68) possessed mitochondrial genomes that clustered closer to the pig parasite Ascaris suum than to A. lumbricoides. Comparative phylogenomic analyses identified over 11 million nuclear-encoded SNPs but just two distinct genetic types that had recombined across the genomes analyzed. The nuclear genomes had extensive heterozygosity, and all samples existed as genetic mosaics with either A. suum-like or A. lumbricoides-like inheritance patterns supporting a highly interbred Ascaris species genetic complex. As no barriers appear to exist for anthroponotic transmission of these 'hybrid' worms, a one-health approach to control the spread of human ascariasis will be necessary.

Genomes: Programmed DNA Elimination in a Parasitic Nematode

Programmed DNA elimination occurs in many eukaryotes. A new study provides a comprehensive view of programmed DNA elimination in a parasitic nematode, defining what sequences are eliminated from which chromosomal locations and presenting a new road map to investigate its molecular mechanism and evolution.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.

Increased viral read counts and metagenomic full genome characterization of porcine astrovirus 4 and Posavirus 1 in sows in a swine farm with unexplained neonatal piglet diarrhea

Neonatal diarrhea in piglets may cause major losses in affected pig herds. The present study used random high-throughput RNA sequencing (metagenomic next generation sequencing, mNGS) to investigate the virome of sows from a farm with persistent neonatal piglet diarrhea in comparison to two control farms without diarrhea problems. A variety of known swine gastrointestinal viruses was detected in the control farms as well as in the problem farm (Mamastrovirus, Enterovirus, Picobirnavirus, Posavirus 1, Kobuvirus, Proprismacovirus). A substantial increase in normalized viral read counts was observed in the affected farm compared to the control farms. The increase was attributable to a single viral species in each of the sampled sows (porcine astrovirus 4 and Posavirus 1). The complete genomes of a porcine astrovirus 4 and two co-infecting Posavirus 1 were de novo assembled and characterized. The 6734 nt single-stranded RNA genome of porcine astrovirus 4 (PoAstV-4) strain Belgium/2019 contains three overlapping open reading frames (nonstructural protein 1ab, nonstructural protein 1a, capsid protein). Posavirus 1 strains Belgium/01/2019 and Belgium/02/2019 have a 9814 nt single-stranded positive-sense RNA genome encoding a single open reading frame (polyprotein precursor) containing the five expected Picornavirales-conserved protein domains. The study highlights the potential of mNGS workflows to study unexplained neonatal diarrhea in piglets and contributes to the scarce availability of both PoAstV-4 and Posavirus-1 whole genome sequences from Western Europe.

Alternative splicing of coq-2 controls the levels of rhodoquinone in animals

Parasitic helminths use two benzoquinones as electron carriers in the electron transport chain. In normoxia, they use ubiquinone (UQ), but in anaerobic conditions inside the host, they require rhodoquinone (RQ) and greatly increase RQ levels. We previously showed the switch from UQ to RQ synthesis is driven by a change of substrates by the polyprenyltransferase COQ-2 (Del Borrello et al., 2019; Roberts Buceta et al., 2019); however, the mechanism of substrate selection is not known. Here, we show helminths synthesize two coq-2 splice forms, coq-2a and coq-2e, and the coq-2e-specific exon is only found in species that synthesize RQ. We show that in Caenorhabditis elegans COQ-2e is required for efficient RQ synthesis and survival in cyanide. Importantly, parasites switch from COQ-2a to COQ-2e as they transit into anaerobic environments. We conclude helminths switch from UQ to RQ synthesis principally via changes in the alternative splicing of coq-2.

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.

Toxocara "omics" and the promises it holds for medicine and veterinary medicine

Toxocariasis is one of the most neglected worldwide zoonoses that is caused by larval nematode parasites of the genus Toxocara, Toxocara canis, and to a lesser extent, Toxocara cati, whose migration mechanism is still largely unknown. Fortunately, some advanced tools have been employed, such as genomics, transcriptomics, and proteomics, to better understand the molecular biology and regulatory mechanisms of Toxocara. Using genomics and transcriptomics, we can identify a large number of genes that participate in the development of Toxocara and the interaction of parasites and their hosts and can predict the functions of unknown genes by comparing them with other relevant species. Using proteomics, we can identify somatic proteins and excretory and secretory (ES) proteins that perform specific biological functions in tissue degradation, pathogen invasion, immune evasion or modulation. These "omics" techniques also can contribute enormously to the development of new drugs, vaccines and diagnostic tools for toxocariasis. In a word, by utilizing "omics", we can better understand the Toxocara and toxocariasis. In this review, we summarized the representative achievements in Toxocara and the interaction between Toxocara spp. and their hosts based on expressed sequence tags (ESTs), microarray gene expression, next-generation sequencing (NGS) technologies and liquid chromatography-tandem mass spectrometry (LC-MS/MS), hoping to better understand the molecular biology of Toxocara, and contribute to new progress in the application areas of new drugs, vaccines and diagnostic tool for toxocariasis in the future.

Elucidating the molecular and developmental biology of parasitic nematodes: Moving to a multiomics paradigm

In the past two decades, significant progress has been made in the sequencing, assembly, annotation and analyses of genomes and transcriptomes of parasitic worms of socioeconomic importance. This progress has somewhat improved our knowledge and understanding of these pathogens at the molecular level. However, compared with the free-living nematode Caenorhabditis elegans, the areas of functional genomics, transcriptomics, proteomics and metabolomics of parasitic nematodes are still in their infancy, and there are major gaps in our knowledge and understanding of the molecular biology of parasitic nematodes. The information on signalling molecules, molecular pathways and microRNAs (miRNAs) that are known to be involved in developmental processes in C. elegans and the availability of some molecular resources (draft genomes, transcriptomes and some proteomes) for selected parasitic nematodes provide a basis to start exploring the developmental biology of parasitic nematodes. Indeed, some studies have identified molecules and pathways that might associate with developmental processes in related, parasitic nematodes, such as Haemonchus contortus (barber's pole worm). However, detailed information is often scant and 'omics resources are limited, preventing a proper integration of 'omic data sets and comprehensive analyses. Moreover, little is known about the functional roles of pheromones, hormones, signalling pathways and post-transcriptional/post-translational regulations in the development of key parasitic nematodes throughout their entire life cycles. Although C. elegans is an excellent model to assist molecular studies of parasitic nematodes, its use is limited when it comes to explorations of processes that are specific to parasitism within host animals. A deep understanding of parasitic nematodes, such as H. contortus, requires substantially enhanced resources and the use of integrative 'omics approaches for analyses. The improved genome and well-established in vitro larval culture system for H. contortus provide unprecedented opportunities for comprehensive studies of the transcriptomes (mRNA and miRNA), proteomes (somatic, excretory/secretory and phosphorylated proteins) and lipidomes (e.g., polar and neutral lipids) of this nematode. Such resources should enable in-depth explorations of its developmental biology at a level, not previously possible. The main aims of this review are (i) to provide a background on the development of nematodes, with a particular emphasis on the molecular aspects involved in the dauer formation and exit in C. elegans; (ii) to critically appraise the current state of knowledge of the developmental biology of parasitic nematodes and identify key knowledge gaps; (iii) to cover salient aspects of H. contortus, with a focus on the recent advances in genomics, transcriptomics, proteomics and lipidomics as well as in vitro culturing systems; (iv) to review recent advances in our knowledge and understanding of the molecular and developmental biology of H. contortus using an integrative multiomics approach, and discuss the implications of this approach for detailed explorations of signalling molecules, molecular processes and pathways likely associated with nematode development, adaptation and parasitism, and for the identification of novel intervention targets against these pathogens. Clearly, the multiomics approach established recently is readily applicable to exploring a wide range of interesting and socioeconomically significant parasitic worms (including also trematodes and cestodes) at the molecular level, and to elucidate host-parasite interactions and disease processes.

Small RNAs in parasitic nematodes - forms and functions

Small RNAs are important regulators of gene expression. They were first identified in Caenorhabditis elegans, but it is now apparent that the main small RNA silencing pathways are functionally conserved across diverse organisms. Availability of genome data for an increasing number of parasitic nematodes has enabled bioinformatic identification of small RNA sequences. Expression of these in different lifecycle stages is revealed by small RNA sequencing and microarray analysis. In this review we describe what is known of the three main small RNA classes in parasitic nematodes – microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs) and small interfering RNAs (siRNAs) – and their proposed functions. miRNAs regulate development in C. elegans and the temporal expression of parasitic nematode miRNAs suggest modulation of target gene levels as parasites develop within the host. miRNAs are also present in extracellular vesicles released by nematodes in vitro, and in plasma from infected hosts, suggesting potential regulation of host gene expression. Roles of piRNAs and siRNAs in suppressing target genes, including transposable elements, are also reviewed. Recent successes in RNAi-mediated gene silencing, and application of small RNA inhibitors and mimics will continue to advance understanding of small RNA functions within the parasite and at the host–parasite interface.

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 (CXCX₄C) 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.

Haemonchus contortus is a gastrointestinal parasitic nematode that causes major economic losses in small ruminants. Here, we investigated the structure and function of a TGF-β type II receptor homologue (Hc-TGFBR2) and its role in regulating H. contortus development. The results showed that the Hc-tgfbr2 gene was transcribed in all developmental stages of H. contortus, with the highest level in L3s and male adults; the encoded protein Hc-TGFBR2 was expressed in the intestine and gonads of adult stages of this nematode. The transcriptional abundance of Hc-tgfbr2 decreased significantly following knockdown by RNA interference in xL3s of H. contortus, which also caused a marked reduction in the number of xL3s developing to L4s in vitro. These findings reveal that the TGF-β type II receptor (Hc-TGFBR2) associates with development of H. contortus, particularly in its transition from the free-living to the parasitic stage.

Profiling microRNAs through development of the parasitic nematode Haemonchus identifies nematode-specific miRNAs that suppress larval development

Parasitic nematodes transition between dramatically different free-living and parasitic stages, with correctly timed development and migration crucial to successful completion of their lifecycle. However little is known of the mechanisms controlling these transitions. microRNAs (miRNAs) negatively regulate gene expression post-transcriptionally and regulate development of diverse organisms. Here we used microarrays to determine the expression profile of miRNAs through development and in gut tissue of the pathogenic nematode Haemonchus contortus. Two miRNAs, mir-228 and mir-235, were enriched in infective L3 larvae, an arrested stage analogous to Caenorhabditis elegans dauer larvae. We hypothesized that these miRNAs may suppress development and maintain arrest. Consistent with this, inhibitors of these miRNAs promoted H. contortus development from L3 to L4 stage, while genetic deletion of C. elegans homologous miRNAs reduced dauer arrest. Epistasis studies with C. elegans daf-2 mutants showed that mir-228 and mir-235 synergise with FOXO transcription factor DAF-16 in the insulin signaling pathway. Target prediction suggests that these miRNAs suppress metabolic and transcription factor activity required for development. Our results provide novel insight into the expression and functions of specific miRNAs in regulating nematode development and identify miRNAs and their target genes as potential therapeutic targets to limit parasite survival within the host.

Advances in Entamoeba histolytica Biology Through Transcriptomic Analysis

A large number of transcriptome-level studies in Entamoeba histolytica, the protozoan parasite that causes amoebiasis, have investigated gene expression patterns to help understand the pathology and biology of the organism. They have compared virulent and avirulent strains in lab culture and after tissue invasion, cells grown under different stress conditions, response to anti-amoebic drug treatments, and gene expression changes during the process of encystation. These studies have revealed interesting molecules/pathways that will help increase our mechanistic understanding of differentially expressed genes during growth perturbations and tissue invasion. Some of the important insights obtained from transcriptome studies include the observations that regulation of carbohydrate metabolism may be an important determinant for tissue invasion, while the novel up-regulated genes during encystation include phospholipase D, and meiotic genes, suggesting the possibility of meiosis during the process. Classification of genes according to expression levels showed that amongst the highly transcribed genes in cultured E. histolytica trophozoites were some virulence factors, raising the question of the role of these factors in normal parasite growth. Promoter motifs associated with differential gene expression and regulation were identified. Some of these motifs associated with high gene expression were located downstream of start codon, and were required for efficient transcription. The listing of E. histolytica genes according to transcript expression levels will help us determine the scale of post-transcriptional regulation, and the possible roles of predicted promoter motifs. The small RNA transcriptome is a valuable resource for detailed structural and functional analysis of these molecules and their regulatory roles. These studies provide new drug targets and enhance our understanding of gene regulation in E. histolytica.

Todos Santos small RNA symposium

Worm biologists from the United States, Canada, and the United Kingdom gathered at the Colorado State University Todos Santos Center in Baja California Sur, Mexico, April 3-5, 2019 for the Todos Santos Small RNA Symposium. Meeting participants, many of whom were still recovering from the bomb cyclone that struck a large swath of North America just days earlier, were greeted by the warmth and sunshine that is nearly ubiquitous in the sleepy seaside town of Todos Santos. With only 24 speakers, the meeting had the sort of laid-back vibe you might expect amongst the palm trees and ocean breeze of the Pacific coast of Mexico. The meeting started with tracing the laboratory lineages of participants. Not surprisingly, the most common parental lineages represented at the meeting were Dr. Craig Mello, Dr. Gary Ruvkun, and Dr. Victor Ambros, whom, together with Dr. Andy Fire and Dr. David Baulcombe, pioneered the small RNA field. In sad irony, on the closing day of the meeting, participants were met with the news of Dr. Sydney Brenner's passing. By establishing the worm, Caenorhabditis elegans, as a model system Dr. Brenner paved the way for much of the research discussed here.

Omics Driven Understanding of the Intestines of Parasitic Nematodes

The biological and molecular complexity of nematodes has impeded research on development of new therapies for treatment and control. We have focused on the versatility of the nematode intestine as a target for new therapies. To that end, it is desirable to establish a broad and deep understanding of the molecular architecture underlying intestinal cell functions at the pan-Nematoda level. Multiomics data were generated to uncover the evolutionary principles underlying both conserved and adaptable features of the nematode intestine. Whole genomes were used to reveal the functional potential of the nematodes, tissue-specific transcriptomes provided a deep assessment of genes that are expressed in the adult nematode intestine, and comparison of selected core species was used to determine a first approximation of the pan-Nematoda intestinal transcriptome. Differentially expressed transcripts were also identified among intestinal regions, with the largest number expressed at significantly higher levels in the anterior region, identifying this region as the most functionally unique compared to middle and posterior regions. Profiling intestinal miRNAs targeting these genes identified the conserved intestinal miRNAs. Proteomics of intestinal cell compartments assigned proteins to several different intestinal cell compartments (intestinal tissue, the integral and peripheral intestinal membranes, and the intestinal lumen). Finally, advanced bioinformatic approaches were used to predict intestinal cell functional categories of seminal importance to parasite survival, which can now be experimentally tested and validated. The data provide the most comprehensive compilation of constitutively and differentially expressed genes, predicted gene regulators, and proteins of the nematode intestine. The information provides knowledge that is essential to understand molecular features of nematode intestinal cells and functions of fundamental importance to the intestine of many, if not all, parasitic nematodes.

Phylogenetic analysis of novel posaviruses detected in feces of Japanese pigs with posaviruses and posa-like viruses of vertebrates and invertebrates

Posaviruses and posa-like viruses are unclassified viruses with sequence similarity to viruses of the order Picornavirales. They have been reported in various vertebrates and invertebrates. We identified 11 posavirus-like sequences in porcine feces and performed phylogenic analysis. Previously reported Japanese posaviruses and those identified in this study clustered with posavirus 1, 4, and 7 and husavirus 1, while five viruses branched into three independent lineages, tentatively named posavirus 10, 11, and 12. Interestingly, posaviruses, except for posavirus 8 and 9, husaviruses, and rasaviruses, formed a cluster consisting of viruses only from pigs, humans, and rats, while posavirus 8 and 9, fisavirus, and basaviruses clustered with posa-like viruses from invertebrates.

Function and Evolution of Nematode RNAi Pathways

Selfish genetic elements, like transposable elements or viruses, are a threat to genomic stability. A variety of processes, including small RNA-based RNA interference (RNAi)-like pathways, has evolved to counteract these elements. Amongst these, endogenous small interfering RNA and Piwi-interacting RNA (piRNA) pathways were implicated in silencing selfish genetic elements in a variety of organisms. Nematodes have several incredibly specialized, rapidly evolving endogenous RNAi-like pathways serving such purposes. Here, we review recent research regarding the RNAi-like pathways of Caenorhabditis elegans as well as those of other nematodes, to provide an evolutionary perspective. We argue that multiple nematode RNAi-like pathways share piRNA-like properties and together form a broad nematode toolkit that allows for silencing of foreign genetic elements.

Diverse RNA interference strategies in early-branching metazoans

BACKGROUND

Micro RNAs (miRNAs) and piwi interacting RNAs (piRNAs), along with the more ancient eukaryotic endogenous small interfering RNAs (endo-siRNAs) constitute the principal components of the RNA interference (RNAi) repertoire of most animals. RNAi in non-bilaterians - sponges, ctenophores, placozoans and cnidarians - appears to be more diverse than that of bilaterians, and includes structurally variable miRNAs in sponges, an enormous number of piRNAs in cnidarians and the absence of miRNAs in ctenophores and placozoans.

RESULTS

Here we identify thousands of endo-siRNAs and piRNAs from the sponge Amphimedon queenslandica, the ctenophore Mnemiopsis leidyi and the cnidarian Nematostella vectensis using a computational approach that clusters mapped small RNA sequences and annotates each cluster based on the read length and relative abundance of the constituent reads. This approach was validated on 11 small RNA libraries in Drosophila melanogaster, demonstrating the successful annotation of RNAi-associated loci with properties consistent with previous reports. In the non-bilaterians we uncover seven new miRNAs from Amphimedon and four from Nematostella as well as sub-populations of candidate cis-natural antisense transcript (cis-NAT) endo-siRNAs. We confirmed the absence of miRNAs in Mnemiopsis but detected an abundance of endo-siRNAs in this ctenophore. Analysis of putative piRNA structure suggests that conserved localised secondary structures in primary transcripts may be important for the production of mature piRNAs in Amphimedon and Nematostella, as is also the case for endo-siRNAs.

CONCLUSION

Together, these findings suggest that the last common ancestor of extant animals did not have the entrained RNAi system that typifies bilaterians. Instead it appears that bilaterians, cnidarians, ctenophores and sponges express unique repertoires and combinations of miRNAs, piRNAs and endo-siRNAs.

Complex virome in feces from Amerindian children in isolated Amazonian villages

The number of viruses circulating in small isolated human populations may be reduced by viral extinctions and rare introductions. Here we used viral metagenomics to characterize the eukaryotic virome in feces from healthy children from a large urban center and from three Amerindian villages with minimal outside contact. Numerous human enteric viruses, mainly from the Picornaviridae and Caliciviridae families, were sequenced from each of the sites. Multiple children from the same villages shed closely related viruses reflecting frequent transmission clusters. Feces of isolated villagers also contained multiple viral genomes of unknown cellular origin from the Picornavirales order and CRESS-DNA group and higher levels of nematode and protozoan DNA. Despite cultural and geographic isolation, the diversity of enteric human viruses was therefore not reduced in these Amazonian villages. Frequent viral introductions and/or increased susceptibility to enteric infections may account for the complex fecal virome of Amerindian children in isolated villages.

A TGF-β type I receptor-like molecule with a key functional role in Haemonchus contortus development

Here we investigated the gene of a transforming growth factor (TGF)-β type I receptor-like molecule in Haemonchus contortus, a highly pathogenic and economically important parasitic nematode of small ruminants. Designated Hc-tgfbr1, this gene is transcribed in all developmental stages of H. contortus, and the encoded protein has glycine-serine rich and kinase domains characteristic of a TGF-β family type I receptor. Expression of a GFP reporter driven by the putative Hc-tgfbr1 promoter localised to two intestinal rings, the anterior-most intestinal ring (int ring I) and the posterior-most intestinal ring (int ring IX) in Caenorhabditis elegans in vivo. Heterologous genetic complementation using a plasmid construct containing Hc-tgfbr1 genomic DNA failed to rescue the function of Ce-daf-1 (a known TGF-β type I receptor gene) in a daf-1-deficient mutant strain of C. elegans. In addition, a TGF-β type I receptor inhibitor, galunisertib, and double-stranded RNA interference (RNAi) were employed to assess the function of Hc-tgfbr1 in the transition from exsheathed L3 (xL3) to the L4 of H. contortus in vitro, revealing that both galunisertib and Hc-tgfbr1-specific double-stranded RNA could retard L4 development. Taken together, these results provide evidence that Hc-tgfbr1 is involved in developmental processes in H. contortus in the transition from the free-living to the parasitic stage.

A Single Mechanism of Biogenesis, Initiated and Directed by PIWI Proteins, Explains piRNA Production in Most Animals

In animals, PIWI-interacting RNAs (piRNAs) guide PIWI proteins to silence transposons and regulate gene expression. The mechanisms for making piRNAs have been proposed to differ among cell types, tissues, and animals. Our data instead suggest a single model that explains piRNA production in most animals. piRNAs initiate piRNA production by guiding PIWI proteins to slice precursor transcripts. Next, PIWI proteins direct the stepwise fragmentation of the sliced precursor transcripts, yielding tail-to-head strings of phased precursor piRNAs (pre-piRNAs). Our analyses detect evidence for this piRNA biogenesis strategy across an evolutionarily broad range of animals, including humans. Thus, PIWI proteins initiate and sustain piRNA biogenesis by the same mechanism in species whose last common ancestor predates the branching of most animal lineages. The unified model places PIWI-clade Argonautes at the center of piRNA biology and suggests that the ancestral animal-the Urmetazoan-used PIWI proteins both to generate piRNA guides and to execute piRNA function.