Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

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.

Whole-genome sequence of the oriental lung fluke Paragonimus westermani

Background

Foodborne infections caused by lung flukes of the genus Paragonimus are a significant and widespread public health problem in tropical areas. Approximately 50 Paragonimus species have been reported to infect animals and humans, but Paragonimus westermani is responsible for the bulk of human disease. Despite their medical and economic importance, no genome sequence for any Paragonimus species is available.

Results

We sequenced and assembled the genome of P. westermani, which is among the largest of the known pathogen genomes with an estimated size of 1.1 Gb. A 922.8 Mb genome assembly was generated from Illumina and Pacific Biosciences (PacBio) sequence data, covering 84% of the estimated genome size. The genome has a high proportion (45%) of repeat-derived DNA, particularly of the long interspersed element and long terminal repeat subtypes, and the expansion of these elements may explain some of the large size. We predicted 12,852 protein coding genes, showing a high level of conservation with related trematode species. The majority of proteins (80%) had homologs in the human liver fluke Opisthorchis viverrini, with an average sequence identity of 64.1%. Assembly of the P. westermani mitochondrial genome from long PacBio reads resulted in a single high-quality circularized 20.6 kb contig. The contig harbored a 6.9 kb region of non-coding repetitive DNA comprised of three distinct repeat units. Our results suggest that the region is highly polymorphic in P. westermani, possibly even within single worm isolates.

Conclusions

The generated assembly represents the first Paragonimus genome sequence and will facilitate future molecular studies of this important, but neglected, parasite group.

SCNBase: a genomics portal for the soybean cyst nematode (Heterodera glycines)

Soybean is an important worldwide crop, and farmers continue to experience significant yield loss due to the soybean cyst nematode (SCN), Heterodera glycines. This soil-borne roundworm parasite is rated the most important pathogen problem in soybean production. The infective nematodes enter into complex interactions with their host plant by inducing the development of specialized plant feeding cells that provide the parasites with nourishment. Addressing the SCN problem will require the development of genomic resources and a global collaboration of scientists to analyze and use these resources. SCNBase.org was designed as a collaborative hub for the SCN genome. All data and analyses are downloadable and can be analyzed with three integrated genomic tools: JBrowse, Feature Search and BLAST. At the time of this writing, a number of genomic and transcriptomic data sets are already available, with 43 JBrowse tracks and 21 category pages describing SCN genomic analyses on gene predictions, transcriptome and read alignments, effector-like genes, expansion and contraction of genomic repeats, orthology and synteny with related nematode species, Single Nucleotide Polymorphism (SNPs) from 15 SCN populations and novel splice sites. Standard functional gene annotations were supplemented with orthologous gene annotations using a comparison to nine related plant-parasitic nematodes, thereby enabling functional annotations for 85% of genes. These annotations led to a greater grasp on the SCN effectorome, which include over 3324 putative effector genes. By designing SCNBase as a hub, future research findings and genomic resources can easily be uploaded and made available for use by others with minimal needs for further curation. By providing these resources to nematode research community, scientists will be empowered to develop novel, more effective SCN management tools.

Genome aware CRISPR gRNA target prediction for parasitic nematodes

The pace of research towards a genetic model to understand the unique molecular biology of parasitic nematodes has increased recently. This research has developed a diverse suite of genetic tools for a variety of parasitic nematodes. CRISPR/Cas9 technology in particular offers much promise as a game changing tool for researchers studying parasitic nematodes. Unlike RNAi, which depends on diverse nematode effectors to silence gene expression, the effectors for CRISPR/Cas9 mutations are typically supplied by the experimenter, making gene editing via CRISPR/Cas9 ideal for testing on genetically intractable nematode systems. To facilitate the development of CRISPR/Cas9 technology for parasitic nematodes, I here describe a tool for identifying gRNA targets and diagnostic primers to a user supplied sequence. The software attempts to minimize non-specific targets by interrogating the genomes of parasitic nematodes. This software is freely available online and features an intuitive interface to help researchers design effective CRISPR experiments for parasitic nematodes.

Whole genome amplification and exome sequencing of archived schistosome miracidia

Adult schistosomes live in the blood vessels and cannot easily be sampled from humans, so archived miracidia larvae hatched from eggs expelled in feces or urine are commonly used for population genetic studies. Large collections of archived miracidia on FTA cards are now available through the Schistosomiasis Collection at the Natural History Museum (SCAN). Here we describe protocols for whole genome amplification of Schistosoma mansoni and Schistosome haematobium miracidia from these cards, as well as real time PCR quantification of amplified schistosome DNA. We used microgram quantities of DNA obtained for exome capture and sequencing of single miracidia, generating dense polymorphism data across the exome. These methods will facilitate the transition from population genetics, using limited numbers of markers to population genomics using genome-wide marker information, maximising the value of collections such as SCAN.

De novo genome and transcriptome analyses provide insights into the biology of the trematode human parasite Fasciolopsis buski

Many trematode parasites cause infection in humans and are thought to be a major public health problem. Their ecological diversity in different regions provides challenging questions on evolution of these organisms. In this report, we perform transcriptome analysis of the giant intestinal fluke, Fasciolopsis buski, using next generation sequencing technology. Short read sequences derived from polyA containing RNA of this organism were assembled into 30,677 unigenes that led to the annotation of 12,380 genes. Annotation of the assembled transcripts enabled insight into processes and pathways in the intestinal fluke, such as RNAi pathway and energy metabolism. The expressed kinome of the organism was characterized by identifying all protein kinases. A rough draft genome assembly for Fasciolopsis buski is also reported herewith with SRA accessions for crosschecking the findings in the analyzed transcriptome data. Transcriptome data also helped us to identify some of the expressed transposable elements. Though many Long Interspersed elements (LINEs) were identified, only two Short Interspersed Elements (SINEs) were visible. Overall transcriptome and draft genome analysis of F. buski helped us to characterize some of its important biological characteristics and provided enormous resources for development of a suitable diagnostic system and anti-parasitic therapeutic molecules.

The small GTPase RAB-35 defines a third pathway that is required for the recognition and degradation of apoptotic cells

In metazoans, apoptotic cells are swiftly engulfed by phagocytes and degraded inside phagosomes. Multiple small GTPases in the Rab family are known to function in phagosome maturation by regulating vesicle trafficking. We discovered rab-35 as a new gene important for apoptotic cell clearance from a genetic screen targeting putative Rab GTPases in Caenorhabditis elegans. We further identified TBC-10 as a putative GTPase-activating protein (GAP), and FLCN-1 and RME-4 as two putative Guanine Nucleotide Exchange Factors (GEFs), for RAB-35. We found that RAB-35 was required for the efficient incorporation of early endosomes to phagosomes and for the timely degradation of apoptotic cell corpses. More specifically, RAB-35 promotes two essential events that initiate phagosome maturation: the switch of phagosomal membrane phosphatidylinositol species from PtdIns(4,5)P2 to PtdIns(3)P, and the recruitment of the small GTPase RAB-5 to phagosomal surfaces. These functions of RAB-35 were previously unknown. Remarkably, although the phagocytic receptor CED-1 regulates these same events, RAB-35 and CED-1 appear to function independently. Upstream of degradation, RAB-35 also facilitates the recognition of apoptotic cells independently of the known CED-1 and CED-5 pathways. RAB-35 localizes to extending pseudopods and is further enriched on nascent phagosomes, consistent with its dual roles in regulating apoptotic cell-recognition and phagosome maturation. Epistasis analyses indicate that rab-35 acts in parallel to both of the canonical ced-1/6/7 and ced-2/5/10/12 clearance pathways. We propose that RAB-35 acts as a robustness factor, defining a novel pathway that aids these canonical pathways in both the recognition and degradation of apoptotic cells.

Streptothricin acetyl transferase 2 (Sat2): A dominant selection marker for Caenorhabditis elegans genome editing

Studies on Caenorhabditis elegans would benefit from the introduction of new selectable markers to allow more complex types of experiments to be conducted with this model animal. We established a new antibiotic selection marker for C. elegans transformation based on nourseothricin (NTC) and its resistance-encoding gene, streptothricin-acetyl transferase 2 (Sat2). NTC was able to efficiently prevent worm development at very low concentrations, and the worms expressing Sat2 were able to survive on the selection plates without any developmental defects. Using CRISPR/Cas9 and NTC selection, we were able to easily insert a 13-kb expression cassette into a defined locus in C. elegans. The structure and spectrum of NTC differs from other antibiotics like hygromycin B and geneticin, making it possible to use NTC alongside them. Indeed, we confirmed NTC-sat2 selection could work with the hygromycin B selection system simultaneously. Thus, the new NTC–Sat2 system can act as a useful dominant marker for gene transfer and genome editing in C. elegans.

Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups

Despite the substantial amount of genomic and transcriptomic data available for a wide range of eukaryotic organisms, most genomes are still in a draft state and can have inaccurate gene predictions. To gain a sound understanding of the biology of an organism, it is crucial that inferred protein sequences are accurately identified and annotated. However, this can be challenging to achieve, particularly for organisms such as parasitic worms (helminths), as most gene prediction approaches do not account for substantial phylogenetic divergence from model organisms, such as Caenorhabditis elegans and Drosophila melanogaster, whose genomes are well-curated. In this paper, we describe a bioinformatic strategy for the curation of gene families and subsequent annotation of encoded proteins. This strategy relies on pairwise gene curation between at least two closely related species using genomic and transcriptomic data sets, and is built on recent work on kinase complements of parasitic worms. Here, we discuss salient technical aspects of this strategy and its implications for the curation of protein families more generally.

G-protein-coupled receptor genes of Dirofilaria immitis

The diversity and uniqueness of nematode heterotrimeric G-protein-coupled receptors (GPCRs) provides impetus for identifying ligands that can be used as therapeutics for treating diseases caused by parasitic nematode infections. In human medicine, GPCRs have represented the largest group of 'drugable' targets exploited in the market today. In the filarial nematode Dirofilaria immitis, which causes heartworm disease, the macrocyclic lactones (ML) have been used as the sole preventatives for more than 25 years and now there is confirmed ML resistance in this parasite. A novel anthelmintic emodepside, with antifilarial activity, can act on a GPCR. In view of the ML resistance, there is an urgent need to identify new drug targets and GPCRs of D. immitis may be promising receptors. Knowledge of polymorphism within the GPCR superfamily is of interest. A total of 127 GPCR genes have been identified, so far, in the genome of D. immitis. Whole genome sequencing data from four ML susceptible and four ML loss of efficacy populations was used to identify 393 polymorphic loci in 35 D. immitis GPCR genes. Out of 57 SNPs in exonic regions, 36 of them caused a change in an amino acid, out of which 2 changed the predicted secondary structure of the protein. Knowledge about GPCR genes and their polymorphism is valuable information for drug design processes. Further studies need to be carried out to more fully understand the implications of each of the SNPs identified by this study.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

Profiling G protein-coupled receptors of Fasciola hepatica identifies orphan rhodopsins unique to phylum Platyhelminthes

G protein-coupled receptors (GPCRs) are established drug targets. Despite their considerable appeal as targets for next-generation anthelmintics, poor understanding of their diversity and function in parasitic helminths has thwarted progress towards GPCR-targeted anti-parasite drugs. This study facilitates GPCR research in the liver fluke, Fasciola hepatica, by generating the first profile of GPCRs from the F. hepatica genome. Our dataset describes 147 high confidence GPCRs, representing the largest cohort of GPCRs, and the largest set of in silico ligand-receptor predictions, yet reported in any parasitic helminth. All GPCRs fall within the established GRAFS nomenclature; comprising three glutamate, 135 rhodopsin, two adhesion, five frizzled, one smoothened, and one secretin GPCR. Stringent annotation pipelines identified 18 highly diverged rhodopsins in F. hepatica that maintained core rhodopsin signatures, but lacked significant similarity with non-flatworm sequences, providing a new sub-group of potential flukicide targets. These facilitated identification of a larger cohort of 76 related sequences from available flatworm genomes, representing new members of existing groups (PROF1/Srfb, Rho-L, Rho-R, Srfa, Srfc) of flatworm-specific rhodopsins. These receptors imply flatworm specific GPCR functions, and/or co-evolution with unique flatworm ligands, and could facilitate the development of exquisitely selective anthelmintics. Ligand binding domain sequence conservation relative to deorphanised rhodopsins enabled high confidence ligand-receptor matching of seventeen receptors activated by acetylcholine, neuropeptide F/Y, octopamine or serotonin. RNA-Seq analyses showed expression of 101 GPCRs across various developmental stages, with the majority expressed most highly in the pathogenic intra-mammalian juvenile parasites. These data identify a broad complement of GPCRs in F. hepatica, including rhodopsins likely to have key functions in neuromuscular control and sensory perception, as well as frizzled and adhesion/secretin families implicated, in other species, in growth, development and reproduction. This catalogue of liver fluke GPCRs provides a platform for new avenues into our understanding of flatworm biology and anthelmintic discovery.

Novel global effector mining from the transcriptome of early life stages of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN) Heterodera glycines is an obligate parasite that relies on the secretion of effector proteins to manipulate host cellular processes that favor the formation of a feeding site within host roots to ensure its survival. The sequence complexity and co-evolutionary forces acting upon these effectors remain unknown. Here we generated a de novo transcriptome assembly representing the early life stages of SCN in both a compatible and an incompatible host interaction to facilitate global effector mining efforts in the absence of an available annotated SCN genome. We then employed a dual effector prediction strategy coupling a newly developed nematode effector prediction tool, N-Preffector, with a traditional secreted protein prediction pipeline to uncover a suite of novel effector candidates. Our analysis distinguished between effectors that co-evolve with the host genotype and those conserved by the pathogen to maintain a core function in parasitism and demonstrated that alternative splicing is one mechanism used to diversify the effector pool. In addition, we confirmed the presence of viral and microbial inhabitants with molecular sequence information. This transcriptome represents the most comprehensive whole-nematode sequence currently available for SCN and can be used as a tool for annotation of expected genome assemblies.

P-glycoprotein-9 and macrocyclic lactone resistance status in selected strains of the ovine gastrointestinal nematode, Teladorsagia circumcincta

The Teladorsagia circumcincta P-glycoprotein-9 (Tci-pgp-9) gene has previously been implicated in multiple-anthelmintic resistance in this parasite. Here we further characterise genetic diversity in Tci-pgp-9 and its possible role in ivermectin (IVM) and multi-drug resistance using two UK field isolates of T. circumcincta, one susceptible to anthelmintics (MTci2) and the other resistant to most available anthelmintics including IVM (MTci5). A comparison of full-length Tci-pgp-9 cDNA transcripts from the MTci2 and MTci5 isolates (∼3.8 kb in both cases) indicated that they shared 95.6% and 99.5% identity at the nucleotide and amino acid levels, respectively. Nine non-synonymous SNPs were found in the MTci5 sequences relative to their MTci2 counterparts. Twelve genomic sequence variants of the first internucleotide binding domain of Tci-pgp-9 were identified and up to 10 of these were present in some individual worms, strongly supporting previous evidence that amplification of this gene has occurred in T. circumcincta. On average, fewer distinct sequence variants of Tci-pgp-9 were present in individual worms of the MTci5 isolate than in those of the MTci2 isolate. A further reduction in the number of sequence variants was observed in individuals derived from an IVM-treated sub-population of MTci5. These findings suggest that Tci-pgp-9 was under purifying selection in the face of IVM treatment in T. circumcincta, with some sequence variants being selected against.

•

Tci-pgp-9 cDNA sequences from the MTci2 and MTci5 isolates shared 95.6% identity.

•

Comparison of Tci-pgp-9 sequences from MTci5 and MTci2 showed 9 non-synonymous SNPs.

•

Multiple Tci-pgp-9-IBDA sequence variants occurred in the majority of individual larvae.

•

Tci-pgp-9 appears to be under purifying selection after anthelmintics are applied.

Rapid genome shrinkage in a self-fertile nematode reveals sperm competition proteins

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.

The Substrates of Nonsense-Mediated mRNA Decay in Caenorhabditis elegans

Nonsense-mediated mRNA decay (NMD) is a conserved pathway that strongly influences eukaryotic gene expression. Inactivating or inhibiting NMD affects the abundance of a substantial fraction of the transcriptome in numerous species. Transcripts whose abundance is altered in NMD-deficient cells may represent either direct substrates of NMD or indirect effects of inhibiting NMD. We present a genome-wide investigation of the direct substrates of NMD in Caenorhabditis elegans Our goals were (i) to identify mRNA substrates of NMD and (ii) to distinguish those mRNAs from others whose abundance is indirectly influenced by the absence of NMD. We previously demonstrated that Upf1p/SMG-2, the central effector of NMD in all studied eukaryotes, preferentially associates with mRNAs that contain premature translation termination codons. We used this preferential association to distinguish direct from indirect effects by coupling immunopurification of Upf1/SMG-2 with high-throughput mRNA sequencing of NMD-deficient mutants and NMD-proficient controls. We identify 680 substrates of NMD, 171 of which contain novel spliced forms that (i) include sequences of annotated introns and (ii) have not been previously documented in the C. elegans transcriptome. NMD degrades unproductively spliced mRNAs with sufficient efficiency in NMD-proficient strains that such mRNAs were not previously known. Two classes of genes are enriched among the identified NMD substrates: (i) mRNAs of expressed pseudogenes and (ii) mRNAs of gene families whose gene number has recently expanded in the C. elegans genome. Our results identify novel NMD substrates and provide a context for understanding NMD's role in normal gene expression and genome evolution.

Comparative Genomics of Gene Loss and Gain in Caenorhabditis and Other Nematodes

Nematodes, such as Caenorhabditis elegans, form one of the most species-rich animal phyla. By now more than 30 nematode genomes have been published allowing for comparative genomic analyses at various different time-scales. The majority of a nematode's gene repertoire is represented by either duplicated or so-called orphan genes of unknown origin. This indicates the importance of mechanisms that generate new genes during the course of evolution. While it is certain that nematodes have acquired genes by horizontal gene transfer from various donors, this process only explains a small portion of the nematode gene content. As evolutionary genomic analyses strongly support that most orphan genes are indeed protein-coding, future studies will have to decide, whether they are result from extreme divergence or evolved de novo from previously noncoding sequences. In this contribution, I summarize several studies investigating gene loss and gain in nematodes and discuss the strengths and weaknesses of individual approaches and datasets. These approaches can be used to ask nematode-specific questions such as associated with the evolution of parasitism or with switches in mating systems, but also can complement studies in other animal phyla like vertebrates and insects to broaden our general view on genome evolution.

A Multi-Omics Database for Parasitic Nematodes and Trematodes

Helminth.net ( www.helminth.net ) is a web-based resource that was launched in 2000 as simply " Nematode.net " to host and investigate gene sequences from nematode genomes. Over the years it has evolved to become the moniker for a collection of databases: Nematode.net and Trematode.net . These databases host information for 73 nematode (roundworms) and 17 trematode (flatworms) species and serve as backbone for a number of tools that allow users to query slices of the data for multifactorial combinations of species-omics properties. Recent focus has been on inclusion of gene and protein expression data, population genomics and cross-kingdom interactions (metagenomics datasets). This chapter describes the website, the available tools and some of the new features.

Using WormBase: A Genome Biology Resource for Caenorhabditis elegans and Related Nematodes

WormBase ( www.wormbase.org ) provides the nematode research community with a centralized database for information pertaining to nematode genes and genomes. As more nematode genome sequences are becoming available and as richer data sets are published, WormBase strives to maintain updated information, displays, and services to facilitate efficient access to and understanding of the knowledge generated by the published nematode genetics literature. This chapter aims to provide an explanation of how to use basic features of WormBase, new features, and some commonly used tools and data queries. Explanations of the curated data and step-by-step instructions of how to access the data via the WormBase website and available data mining tools are provided.

Tissue-specific transcriptome analyses provide new insights into GPCR signalling in adult Schistosoma mansoni

Schistosomes are blood-dwelling trematodes with global impact on human and animal health. Because medical treatment is currently based on a single drug, praziquantel, there is urgent need for the development of alternative control strategies. The Schistosoma mansoni genome project provides a platform to study and connect the genetic repertoire of schistosomes to specific biological functions essential for successful parasitism. G protein-coupled receptors (GPCRs) form the largest superfamily of transmembrane receptors throughout the Eumetazoan phyla, including platyhelminths. Due to their involvement in diverse biological processes, their pharmacological importance, and proven druggability, GPCRs are promising targets for new anthelmintics. However, to identify candidate receptors, a more detailed understanding of the roles of GPCR signalling in schistosome biology is essential. An updated phylogenetic analysis of the S. mansoni GPCR genome (GPCRome) is presented, facilitated by updated genome data that allowed a more precise annotation of GPCRs. Additionally, we review the current knowledge on GPCR signalling in this parasite and provide new insights into the potential roles of GPCRs in schistosome reproduction based on the findings of a recent tissue-specific transcriptomic study in paired and unpaired S. mansoni. According to the current analysis, GPCRs contribute to gonad-specific functions but also to nongonad, pairing-dependent processes. The latter may regulate gonad-unrelated functions during the multifaceted male-female interaction. Finally, we compare the schistosome GPCRome to that of another parasitic trematode, Fasciola, and discuss the importance of GPCRs to basic and applied research. Phylogenetic analyses display GPCR diversity in free-living and parasitic platyhelminths and suggest diverse functions in schistosomes. Although their roles need to be substantiated by functional studies in the future, the data support the selection of GPCR candidates for basic and applied studies, invigorating the exploitation of this important receptor class for drug discovery against schistosomes but also other trematodes.

Conservation of a microRNA cluster in parasitic nematodes and profiling of miRNAs in excretory-secretory products and microvesicles of Haemonchus contortus

microRNAs are small non-coding RNAs that are important regulators of gene expression in a range of animals, including nematodes. We have analysed a cluster of four miRNAs from the pathogenic nematode species Haemonchus contortus that are closely linked in the genome. We find that the cluster is conserved only in clade V parasitic nematodes and in some ascarids, but not in other clade III species nor in clade V free-living nematodes. Members of the cluster are present in parasite excretory-secretory products and can be detected in the abomasum and draining lymph nodes of infected sheep, indicating their release in vitro and in vivo. As observed for other parasitic nematodes, H. contortus adult worms release extracellular vesicles (EV). Small RNA libraries were prepared from vesicle-enriched and vesicle-depleted supernatants from both adult worms and L4 stage larvae. Comparison of the miRNA species in the different fractions indicated that specific miRNAs are packaged within vesicles, while others are more abundant in vesicle-depleted supernatant. Hierarchical clustering analysis indicated that the gut is the likely source of vesicle-associated miRNAs in the L4 stage, but not in the adult worm. These findings add to the growing body of work demonstrating that miRNAs released from parasitic helminths may play an important role in host-parasite interactions.

Different species of parasitic worms release microRNAs, small non-coding RNA molecules, some of which are known to interact with host genes to alter the immune response. We characterized a cluster of four microRNAs from Haemonchus contortus, an important parasitic nematode of livestock. The miRNA cluster appeared to be present only in nematode worms that inhabit the intestinal tract and members of the cluster could be detected in the excretory-secretory products of adult H. contortus. Some of these miRNAs were also detected from the site of infection within the sheep host and from the draining lymph node suggesting that they may be released in vivo as well as in vitro. Analysis of excretory-secretory products released from adult worms and L4 larval stages demonstrated that specific miRNAs were contained within EV. These findings support a role for secreted miRNAs in the host-parasite relationship.

Sequencing of the complete mitochondrial genome of a fish-parasitic flatworm Paratetraonchoides inermis (Platyhelminthes: Monogenea): tRNA gene arrangement reshuffling and implications for phylogeny

BACKGROUND

Paratetraonchoides inermis (Monogenea: Tetraonchoididae) is a flatworm parasitising the gills of uranoscopid fishes. Its morphological characteristics are ambiguous, and molecular data have never been used to study its phylogenetic relationships, which makes its taxonomic classification controversial. Also, several decades of unsuccessful attempts to resolve the relationships within the Monogenea present a strong indication that morphological datasets may not be robust enough to be used to infer evolutionary histories. As the use of molecular data is currently severely limited by their scarcity, we have sequenced and characterized the complete mitochondrial (mt) genome of P. inermis. To investigate its phylogenetic position, we performed phylogenetic analyses using Bayesian inference and maximum likelihood approaches using concatenated amino acid sequences of all 12 protein-coding genes on a dataset containing all available monogenean mt genomes.

RESULTS

The circular mt genome of P. inermis (14,654 bp) contains the standard 36 genes: 22 tRNAs, two rRNAs, 12 protein-encoding genes (PCGs; Atp8 is missing) and a major non-coding region (mNCR). All genes are transcribed from the same strand. The A + T content of the whole genome (82.6%), as well as its elements, is the highest reported among the monogeneans thus far. Three tRNA-like cloverleaf structures were found in mNCR. Several results of the phylogenomic analysis are in disagreement with previously proposed relationships: instead of being closely related to the Gyrodactylidea, Tetraonchidea exhibit a phylogenetic affinity with the Dactylogyridea + Capsalidea clade; and the order Capsalidea is neither basal within the subclass Monopisthocotylea, nor groups with the Gyrodactylidea, but instead forms a sister clade with the Dactylogyridea. The mt genome of P. inermis exhibits a unique gene order, with an extensive reorganization of tRNAs. Monogenea exhibit exceptional gene order plasticity within the Neodermata.

CONCLUSIONS

This study shows that gene order within monopisthocotylid mt genomes is evolving at uneven rates, which creates misleading evolutionary signals. Furthermore, our results indicate that all previous attempts to resolve the evolutionary history of the Monogenea may have produced at least partially erroneous relationships. This further corroborates the necessity to generate more molecular data for this group of parasitic animals.

Genome Architecture and Evolution of a Unichromosomal Asexual Nematode

Asexual reproduction in animals, though rare, is the main or exclusive mode of reproduction in some long-lived lineages. The longevity of asexual clades may be correlated with the maintenance of heterozygosity by mechanisms that rearrange genomes and reduce recombination. Asexual species thus provide an opportunity to gain insight into the relationship between molecular changes, genome architecture, and cellular processes. Here we report the genome sequence of the parthenogenetic nematode Diploscapter pachys with only one chromosome pair. We show that this unichromosomal architecture is shared by a long-lived clade of asexual nematodes closely related to the genetic model organism Caenorhabditis elegans. Analysis of the genome assembly reveals that the unitary chromosome arose through fusion of six ancestral chromosomes, with extensive rearrangement among neighboring regions. Typical nematode telomeres and telomeric protection-encoding genes are lacking. Most regions show significant heterozygosity; homozygosity is largely concentrated to one region and attributed to gene conversion. Cell-biological and molecular evidence is consistent with the absence of key features of meiosis I, including synapsis and recombination. We propose that D. pachys preserves heterozygosity and produces diploid embryos without fertilization through a truncated meiosis. As a prelude to functional studies, we demonstrate that D. pachys is amenable to experimental manipulation by RNA interference.

HSF-1 is a regulator of miRNA expression in Caenorhabditis elegans

The ability of an organism to sense and adapt to environmental stressors is essential for proteome maintenance and survival. The highly conserved heat shock response is a survival mechanism employed by all organisms, including the nematode Caenorhabditis elegans, upon exposure to environmental extremes. Transcriptional control of the metazoan heat shock response is mediated by the heat shock transcription factor HSF-1. In addition to regulating global stress-responsive genes to promote stress-resistance and survival, HSF-1 has recently been shown to regulate stress-independent functions in controlling development, metabolism, and longevity. However, the indirect role of HSF-1 in coordinating stress-dependent and -independent processes through post-transcriptional regulation is largely unknown. MicroRNAs (miRNAs) have emerged as a class of post-transcriptional regulators that control gene expression through translational repression or mRNA degradation. To determine the role of HSF-1 in regulating miRNA expression, we have performed high-throughput small RNA-sequencing in C. elegans grown in the presence and absence of hsf-1 RNAi followed by treatment with or without heat shock. This has allowed us to uncover the miRNAs regulated by HSF-1 via heat-dependent and -independent mechanisms. Integrated miRNA/mRNA target-prediction analyses suggest HSF-1 as a post-transcriptional regulator of development, metabolism, and longevity through regulating miRNA expression. This provides new insight into the possible mechanism by which HSF-1 controls these processes. We have also uncovered oxidative stress response factors and insulin-like signaling factors as a common link between processes affected by HSF-1-regulated miRNAs in stress-dependent and -independent mechanisms, respectively. This may provide a role for miRNAs in regulating cross-talk between various stress responses. Our work therefore uncovers an interesting potential role for HSF-1 in post-transcriptionally controlling gene expression in C. elegans, and suggests a mechanism for cross-talk between stress responses.