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.

Genome Evolution of Plant-Parasitic Nematodes

Plant parasitism has evolved independently on at least four separate occasions in the phylum Nematoda. The application of next-generation sequencing (NGS) to plant-parasitic nematodes has allowed a wide range of genome- or transcriptome-level comparisons, and these have identified genome adaptations that enable parasitism of plants. Current genome data suggest that horizontal gene transfer, gene family expansions, evolution of new genes that mediate interactions with the host, and parasitism-specific gene regulation are important adaptations that allow nematodes to parasitize plants. Sequencing of a larger number of nematode genomes, including plant parasites that show different modes of parasitism or that have evolved in currently unsampled clades, and using free-living taxa as comparators would allow more detailed analysis and a better understanding of the organization of key genes within the genomes. This would facilitate a more complete understanding of the way in which parasitism has shaped the genomes of plant-parasitic nematodes.

Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Mapping-by-sequencing has become a standard method to map and identify phenotype-causing mutations in model species. Here, we show that a fragmented draft assembly is sufficient to perform mapping-by-sequencing in nonmodel species. We generated a draft assembly and annotation of the genome of the free-living nematode Oscheius tipulae, a distant relative of the model Caenorhabditis elegans We used this draft to identify the likely causative mutations at the O. tipulae cov-3 locus, which affect vulval development. The cov-3 locus encodes the O. tipulae ortholog of C. elegans mig-13, and we further show that Cel-mig-13 mutants also have an unsuspected vulval-development phenotype. In a virtuous circle, we were able to use the linkage information collected during mutant mapping to improve the genome assembly. These results showcase the promise of genome-enabled forward genetics in nonmodel species.

Reliable CRISPR/Cas9 Genome Engineering in Caenorhabditis elegans Using a Single Efficient sgRNA and an Easily Recognizable Phenotype

CRISPR/Cas9 genome engineering strategies allow the directed modification of the Caenorhabditis elegans genome to introduce point mutations, generate knock-out mutants, and insert coding sequences for epitope or fluorescent tags. Three practical aspects, however, complicate such experiments. First, the efficiency and specificity of single-guide RNAs (sgRNA) cannot be reliably predicted. Second, the detection of animals carrying genome edits can be challenging in the absence of clearly visible or selectable phenotypes. Third, the sgRNA target site must be inactivated after editing to avoid further double-strand break events. We describe here a strategy that addresses these complications by transplanting the protospacer of a highly efficient sgRNA into a gene of interest to render it amenable to genome engineering. This sgRNA targeting the dpy-10 gene generates genome edits at comparatively high frequency. We demonstrate that the transplanted protospacer is cleaved at the same time as the dpy-10 gene. Our strategy generates scarless genome edits because it no longer requires the introduction of mutations in endogenous sgRNA target sites. Modified progeny can be easily identified in the F1 generation, which drastically reduces the number of animals to be tested by PCR or phenotypic analysis. Using this strategy, we reliably generated precise deletion mutants, transcriptional reporters, and translational fusions with epitope tags and fluorescent reporter genes. In particular, we report here the first use of the new red fluorescent protein mScarlet in a multicellular organism. wrmScarlet, a C. elegans-optimized version, dramatically surpassed TagRFP-T by showing an eightfold increase in fluorescence in a direct comparison.

Histone deacetylase inhibition modulates histone acetylation at gene promoter regions and affects genome-wide gene transcription in Schistosoma mansoni

BACKGROUND

Schistosomiasis is a parasitic disease infecting hundreds of millions of people worldwide. Treatment depends on a single drug, praziquantel, which kills the Schistosoma spp. parasite only at the adult stage. HDAC inhibitors (HDACi) such as Trichostatin A (TSA) induce parasite mortality in vitro (schistosomula and adult worms), however the downstream effects of histone hyperacetylation on the parasite are not known.

METHODOLOGY/PRINCIPAL FINDINGS

TSA treatment of adult worms in vitro increased histone acetylation at H3K9ac and H3K14ac, which are transcription activation marks, not affecting the unrelated transcription repression mark H3K27me3. We investigated the effect of TSA HDACi on schistosomula gene expression at three different time points, finding a marked genome-wide change in the transcriptome profile. Gene transcription activity was correlated with changes on the chromatin acetylation mark at gene promoter regions. Moreover, combining expression data with ChIP-Seq public data for schistosomula, we found that differentially expressed genes having the H3K4me3 mark at their promoter region in general showed transcription activation upon HDACi treatment, compared with those without the mark, which showed transcription down-regulation. Affected genes are enriched for DNA replication processes, most of them being up-regulated. Twenty out of 22 genes encoding proteins involved in reducing reactive oxygen species accumulation were down-regulated. Dozens of genes encoding proteins with histone reader motifs were changed, including SmEED from the PRC2 complex. We targeted SmEZH2 methyltransferase PRC2 component with a new EZH2 inhibitor (GSK343) and showed a synergistic effect with TSA, significantly increasing schistosomula mortality.

CONCLUSIONS/SIGNIFICANCE

Genome-wide gene expression analyses have identified important pathways and cellular functions that were affected and may explain the schistosomicidal effect of TSA HDACi. The change in expression of dozens of histone reader genes involved in regulation of the epigenetic program in S. mansoni can be used as a starting point to look for possible novel schistosomicidal targets.

Comparative analyses within Gyrodactylus (Platyhelminthes: Monogenea) mitochondrial genomes and conserved polymerase chain reaction primers for gyrodactylid mitochondrial DNA

In this study, we describe the complete mitochondrial genomes of Gyrodactylus brachymystacis and Gyrodactylus parvae infecting rainbow trout (Oncorhynchus mykiss) and the invasive topmouth gudgeon (Pseudorasbora parva), respectively. The two circular genomes have a common genome organization found in other Gyrodactylus species. Comparative analyses of mitochondrial genomes from six Gyrodactylus species were carried out to determine base composition, codon usage, transfer RNA and ribosomal RNA genes, major non-coding regions, and nucleotide diversity within the genus. We also provide the first universal models of the secondary structures of rrnS and rrnL for this group thereby promoting utilization of these genetic markers. Universal primers provided herein can be used to obtain more mitochondrial information for pathogen identification and may reveal different levels of molecular phylogenetic inferences for this lineage.

Lessons from the genomes and transcriptomes of filarial nematodes

Human filarial infections are a leading cause of morbidity in the developing world. While a small arsenal of drugs exists to treat these infections, there remains a tremendous need for the development of additional interventions. Recent genome sequences and transcriptome analyses of filarial nematodes have provided novel biological insight and allowed for the prediction of novel drug targets as well as potential vaccine candidates. In this review, we discuss the currently available data, insights gained into the metabolism of these organisms, and how the filaria field can move forward by leveraging these data.

Genomes of Fasciola hepatica from the Americas Reveal Colonization with Neorickettsia Endobacteria Related to the Agents of Potomac Horse and Human Sennetsu Fevers

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis' gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans.

PENELOPE-LIKE RETROTRANSPOSONS IN THE GENOME OF ASIAN BLOOD FLUKE SCHISTOSOMA JAPONICUM (TREMATODA: SCHISTOSOMATIDAE)

Mobile genetic elements (in particular, retrotransposons) constitute a considerable part of eukaryotic genomes. These elements are one of the main sources of variability and instability. Although mobile genetic elements had been discovered more than 60 years ago, certain classes of these elements, such as Penelope-like elements (PLE) remain insufficiently understood. In this work, the results of in silico research into PLE of the S. japonicum genome (one of the most important human parasites from the epidemiological viewpoint) are presented. A new estimate of PLE in the S. japonicum genome was made. It was shown that their high heterogeneity presumably reflects several transposition events. Two groups of canonical PLE copies differing by the amino sequence in the domains RT and EN were found to be present in the Schistosoma genome. The data on the representation, structural features, and possible functional load and evolution of Penelope-like retroelements of the parasite are discussed.

Genome-Wide Analyses of Individual Strongyloides stercoralis (Nematoda: Rhabditoidea) Provide Insights into Population Structure and Reproductive Life Cycles

The helminth Strongyloides stercoralis, which is transmitted through soil, infects 30–100 million people worldwide. S. stercoralis reproduces sexually outside the host as well as asexually within the host, which causes a life-long infection. To understand the population structure and transmission patterns of this parasite, we re-sequenced the genomes of 33 individual S. stercoralis nematodes collected in Myanmar (prevalent region) and Japan (non-prevalent region). We utilised a method combining whole genome amplification and next-generation sequencing techniques to detect 298,202 variant positions (0.6% of the genome) compared with the reference genome. Phylogenetic analyses of SNP data revealed an unambiguous geographical separation and sub-populations that correlated with the host geographical origin, particularly for the Myanmar samples. The relatively higher heterozygosity in the genomes of the Japanese samples can possibly be explained by the independent evolution of two haplotypes of diploid genomes through asexual reproduction during the auto-infection cycle, suggesting that analysing heterozygosity is useful and necessary to infer infection history and geographical prevalence.

Strongyloides stercoralis, one of the most neglected helminths causes strongyloidiasis mainly in tropical and subtropical regions worldwide. The parasite’s complex lifecycle includes sexual and asexual reproduction outside and inside the host, respectively. The parasite can also asexually complete a life cycle within the host's body, which is called autoinfection causing life-long infection. In order to investigate the population structure and transmission patterns of this parasite we sequenced individual nematodes isolated from human faeces in Japan and Myanmar, where the parasite is present at low and high frequencies, respectively. Whole genome sequencing of small parasites is generally difficult because the amount of DNA is limiting. However, we overcame this problem by combining whole genome amplification with next-generation sequencing. Sequence comparisons revealed 0.6% of the genome is variable among samples, and the variants showed clear separation by the location of their origin. We found that heterozygosity within the genomes was higher in Japan, which is likely explained by the predominance of asexual reproduction through auto-infection, suggesting that analyses of heterozygosity are required to better understand the history of a population.

Fasciola hepatica demonstrates high levels of genetic diversity, a lack of population structure and high gene flow: possible implications for drug resistance

•

Self-fertilisation does occur but is rare in field populations of Fasciola hepatica.

•

Some hosts harboured genotypically identical parasites (clones).

•

The presence of clones is consistent with clonal expansion and clumped transmission.

•

84% of 1579 F. hepatica had unique genotypes, indicating high genetic diversity.

•

We found high gene flow, no population structure and low self-fertilisation rate.

Fasciola hepatica, the liver fluke, is a trematode parasite of considerable economic importance to the livestock industry and is a re-emerging zoonosis that poses a risk to human health in F. hepatica-endemic areas worldwide. Drug resistance is a substantial threat to the current and future control of F. hepatica, yet little is known about how the biology of the parasite influences the development and spread of resistance. Given that F. hepatica can self-fertilise and therefore inbreed, there is the potential for greater population differentiation and an increased likelihood of recessive alleles, such as drug resistance genes, coming together. This could be compounded by clonal expansion within the snail intermediate host and aggregation of parasites of the same genotype on pasture. Alternatively, widespread movement of animals that typically occurs in the UK could promote high levels of gene flow and prevent population differentiation. We identified clonal parasites with identical multilocus genotypes in 61% of hosts. Despite this, 84% of 1579 adult parasites had unique multilocus genotypes, which supports high levels of genotypic diversity within F. hepatica populations. Our analyses indicate a selfing rate no greater than 2%, suggesting that this diversity is in part due to the propensity for F. hepatica to cross-fertilise. Finally, although we identified high genetic diversity within a given host, there was little evidence for differentiation between populations from different hosts, indicating a single panmictic population. This implies that, once those emerge, anthelmintic resistance genes have the potential to spread rapidly through liver fluke populations.