WormBase 2024: status and transitioning to Alliance infrastructure

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

Recombination in bdelloid rotifer genomes: asexuality, transfer and stress

Bdelloid rotifers constitute a class of microscopic animals living in freshwater habitats worldwide. Several strange features of bdelloids have drawn attention: their ability to tolerate desiccation and other stresses, a lack of reported males across the clade despite centuries of study, and unusually high numbers of horizontally acquired, non-metazoan genes. Genome sequencing is transforming our understanding of their lifestyle and its consequences, while in turn providing wider insights about recombination and genome organisation in animals. Many questions remain, not least how to reconcile apparent genomic signatures of sex with the continued absence of reported males, why bdelloids have so many horizontally acquired genes, and how their remarkable ability to survive stress interacts with recombination and other genomic processes.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

Expression of Hairpin-Enriched Mitochondrial DNA in Two Hairworm Species (Nematomorpha)

Nematomorpha (hairworms) is a phylum of parasitic ecdysozoans, best known for infecting arthropods and guiding their hosts toward water, where the parasite can complete its life cycle. Over 350 species of nematomorphs have been described, yet molecular data for the group remain scarce. The few available mitochondrial genomes of nematomorphs are enriched with long inverted repeats, which are embedded in the coding sequences of their genes-a remarkably unusual feature exclusive to this phylum. Here, we obtain and annotate the repeats in the mitochondrial genome of another nematomorph species-Parachordodes pustulosus. Using genomic and transcriptomic libraries, we investigate the impact of inverted repeats on the read coverage of the mitochondrial genome. Pronounced drops in the read coverage coincide with regions containing long inverted repeats, denoting the 'blind spots' of short-fragment sequencing libraries. Phylogenetic inference with the novel data reveals multiple disagreements between the traditional system of Nematomorpha and molecular data, rendering several genera paraphyletic, including Parachordodes.

Population differentiation and epidemic tracking of Bursaphelenchus xylophilus in China based on chromosome-level assembly and whole-genome sequencing data

BACKGROUND

Bursaphelenchus xylophilus, the pinewood nematode, kills millions of pine trees worldwide every year, and causes enormous economic and ecological losses. Despite extensive research on population variation, there is little understanding of the population-wide variation spectrum in China.

RESULTS

We sequenced an inbred B. xylophilus strain using Pacbio+Illumina+Bionano+Hi-C and generated a chromosome-level assembly (AH1) with six chromosomes of 77.1 Mb (chromosome N50: 12 Mb). The AH1 assembly shows very high continuity and completeness, and contains novel genes with potentially important functions compared with previous assemblies. Subsequently, we sequenced 181 strains from China and the USA and found ~7.8 million single nucleotide polymorphisms (SNPs). Analysis shows that the B. xylophilus population in China can be divided into geographically bounded subpopulations with severe cross-infection and potential migrations. In addition, distribution of B. xylophilus is dominated by temperature zones while geographically associated SNPs are mainly located on adaptation related GPCR gene families, suggesting the nematode has been evolving to adapt to different temperatures. A machine-learning based epidemic tracking method has been established to predict their geographical origins, which can be applied to any other species.

CONCLUSION

Our study provides the community with the first high-quality chromosome-level assembly which includes a comprehensive catalogue of genetic variations. It provides insights into population structure and effective tracking method for this invasive species, which facilitates future studies to address a variety of applied, genomic and evolutionary questions in B. xylophilus as well as related species.

Approaches for CRISPR/Cas9 Genome Editing in C. elegans

The clustered, regularly interspaced, short, palindromic repeat (CRISPR)-associated (CAS) nuclease Cas9 has been used in many organisms to generate specific mutations and transgene insertions. Here we describe our most up-to-date protocols using the S. pyogenes Cas9 in C. elegans that provides a convenient and effective approach for making heritable changes to the worm genome. We present several considerations when deciding which strategy best suits the needs of the experiment.

Genome-wide transcriptome analysis of the early developmental stages of Echinococcus granulosus protoscoleces reveals extensive alternative splicing events in the spliceosome pathway

BACKGROUND

The complex life cycle of Echinococcus granulosus involves numerous environmental conditions within different intermediate and definitive hosts. This requires adaptation at different levels of transcript regulation. Alternative splicing (AS) and the related cellular functions as one of the major fields of post-genomics has been poorly studied in tapeworms. In the present study, we investigated AS events and their potential biological effects in E. granulosus.

METHODS

Whole transcriptome sequencing data of four groups of protoscoleces were prepared for RNA-seq library construction. Fresh protoscoleces were either used as non-induced controls (NT group) or incubated for 15 min with pepsin (PEP group) and cultivated in a biphasic medium for 12 and 24 h (12 and 24 h groups). The frequency and different types of AS events were identified using rMATS software. Functional annotations and gene ontology of differential AS (DAS) genes were performed using Blast2GO software. AS events were experimentally validated by PCR on the protoscolex cDNAs using specific primers for each gene.

RESULTS

At least one AS event was found in 38.1% of the genes (3904 out of 10,245) in the protoscoleces during early strobilar development. The genes were associated primarily with cellular and metabolic processes and binding and catalytic activity. KEGG pathway analysis of DAS events revealed a number of genes belonging to different components of the spliceosome complex. These genes tended to belong to common SR proteins, U1-related factors, U2-related factors, complex A-specific factors and other splicing-related proteins.

CONCLUSIONS

The high number of AS events in the transcriptome regulatory mechanisms indicates the essential rapid molecular processes required by the parasite for adaptation in different environments.

High-efficiency CRISPR gene editing in C. elegans using Cas9 integrated into the genome

Gene editing in C. elegans using plasmid-based CRISPR reagents requires microinjection of many animals to produce a single edit. Germline silencing of plasmid-borne Cas9 is a major cause of inefficient editing. Here, we present a set of C. elegans strains that constitutively express Cas9 in the germline from an integrated transgene. These strains markedly improve the success rate for plasmid-based CRISPR edits. For simple, short homology arm GFP insertions, 50-100% of injected animals typically produce edited progeny, depending on the target locus. Template-guided editing from an extrachromosomal array is maintained over multiple generations. We have built strains with the Cas9 transgene on multiple chromosomes. Additionally, each Cas9 locus also contains a heatshock-driven Cre recombinase for selectable marker removal and a bright fluorescence marker for easy outcrossing. These integrated Cas9 strains greatly reduce the workload for producing individual genome edits.

Hookworm infection in central China: morphological and molecular diagnosis

BACKGROUND

Necator americanus is one of the major etiological agents of human ancylostomiasis. Historically, the epidemiology of ancylostomiasis in Henan Province of central China and the molecular characteristics of N. americanus have been poorly understood.

METHODS

In this study, we report a case of ancylostomiasis in Zhengzhou city of Henan Province. We also review the epidemiology of ancylostomiasis in Henan Province from 1949 to 2020. In addition, the complete mitochondrial (mt) genome of one clinical isolate is fully characterized using Illumina sequencing. All available mt genomes of hookworms in GenBank were included to reconstruct the phylogeny using both maximum likelihood (ML) and Bayesian inference (BI) methods.

RESULTS

A total of three worms were collected from the patient. These worms were identified as N. americanus based on morphological characteristics as well as confirmed by genotyping with the barcoding gene cox1. Although ancylostomiasis cases have dropped substantially in recent years, hookworm infection is still a public health problem in underdeveloped areas and remote rural areas in Henan Province. The mt genome features of the N. americanus contained 12 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a major non-coding region. The nad1 gene showed high sequence variability among isolates, which is worth considering for future genetic studies of N. americanus. Phylogenetic analyses support the monophyly of hookworm isolates from different hosts and distinct geographical locations.

CONCLUSIONS

The mt genome of N. americanus presented here will serve as a useful data set for studying population genetics and phylogenetic relationships of hookworms. Positive measures for preventing and controlling ancylostomiasis are required by both health services and individuals in Henan Province.

C. elegans genome-wide analysis reveals DNA repair pathways that act cooperatively to preserve genome integrity upon ionizing radiation

Ionizing radiation (IR) is widely used in cancer therapy and accidental or environmental exposure is a major concern. However, little is known about the genome-wide effects IR exerts on germ cells and the relative contribution of DNA repair pathways for mending IR-induced lesions. Here, using C. elegans as a model system and using primary sequencing data from our recent high-level overview of the mutagenic consequences of 11 genotoxic agents, we investigate in detail the genome-wide mutagenic consequences of exposing wild-type and 43 DNA repair and damage response defective C. elegans strains to a Caesium (Cs-137) source, emitting γ-rays. Cs-137 radiation induced single nucleotide variants (SNVs) at a rate of ~1 base substitution per 3 Gy, affecting all nucleotides equally. In nucleotide excision repair mutants, this frequency increased 2-fold concurrently with increased dinucleotide substitutions. As observed for DNA damage induced by bulky DNA adducts, small deletions were increased in translesion polymerase mutants, while base changes decreased. Structural variants (SVs) were augmented with dose, but did not arise with significantly higher frequency in any DNA repair mutants tested. Moreover, 6% of all mutations occurred in clusters, but clustering was not significantly altered in any DNA repair mutant background. Our data is relevant for better understanding how DNA repair pathways modulate IR-induced lesions.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

High-quality reference genome of Fasciola gigantica: Insights into the genomic signatures of transposon-mediated evolution and specific parasitic adaption in tropical regions

Fasciola gigantica and Fasciola hepatica are causative pathogens of fascioliasis, with the widest latitudinal, longitudinal, and altitudinal distribution; however, among parasites, they have the largest sequenced genomes, hindering genomic research. In the present study, we used various sequencing and assembly technologies to generate a new high-quality Fasciola gigantica reference genome. We improved the integration of gene structure prediction, and identified two independent transposable element expansion events contributing to (1) the speciation between Fasciola and Fasciolopsis during the Cretaceous-Paleogene boundary mass extinction, and (2) the habitat switch to the liver during the Paleocene-Eocene Thermal Maximum, accompanied by gene length increment. Long interspersed element (LINE) duplication contributed to the second transposon-mediated alteration, showing an obvious trend of insertion into gene regions, regardless of strong purifying effect. Gene ontology analysis of genes with long LINE insertions identified membrane-associated and vesicle secretion process proteins, further implicating the functional alteration of the gene network. We identified 852 predicted excretory/secretory proteins and 3300 protein-protein interactions between Fasciola gigantica and its host. Among them, copper/zinc superoxide dismutase genes, with specific gene copy number variations, might play a central role in the phase I detoxification process. Analysis of 559 single-copy orthologs suggested that Fasciola gigantica and Fasciola hepatica diverged at 11.8 Ma near the Middle and Late Miocene Epoch boundary. We identified 98 rapidly evolving gene families, including actin and aquaporin, which might explain the large body size and the parasitic adaptive character resulting in these liver flukes becoming epidemic in tropical and subtropical regions.

Linked surveillance and genetic data uncovers programmatically relevant geographic scale of Guinea worm transmission in Chad

Background

Guinea worm (Dracunculus medinensis) was detected in Chad in 2010 after a supposed ten-year absence, posing a challenge to the global eradication effort. Initiation of a village-based surveillance system in 2012 revealed a substantial number of dogs infected with Guinea worm, raising questions about paratenic hosts and cross-species transmission.

Methodology/principal findings

We coupled genomic and surveillance case data from 2012-2018 to investigate the modes of transmission between dog and human hosts and the geographic connectivity of worms. Eighty-six variants across four genes in the mitochondrial genome identified 41 genetically distinct worm genotypes. Spatiotemporal modeling revealed worms with the same genotype (‘genetically identical’) were within a median range of 18.6 kilometers of each other, but largely within approximately 50 kilometers. Genetically identical worms varied in their degree of spatial clustering, suggesting there may be different factors that favor or constrain transmission. Each worm was surrounded by five to ten genetically distinct worms within a 50 kilometer radius. As expected, we observed a change in the genetic similarity distribution between pairs of worms using variants across the complete mitochondrial genome in an independent population.

Conclusions/significance

In the largest study linking genetic and surveillance data to date of Guinea worm cases in Chad, we show genetic identity and modeling can facilitate the understanding of local transmission. The co-occurrence of genetically non-identical worms in quantitatively identified transmission ranges highlights the necessity for genomic tools to link cases. The improved discrimination between pairs of worms from variants identified across the complete mitochondrial genome suggests that expanding the number of genomic markers could link cases at a finer scale. These results suggest that scaling up genomic surveillance for Guinea worm may provide additional value for programmatic decision-making critical for monitoring cases and intervention efficacy to achieve elimination.

The global eradication effort for Guinea worm disease has dramatically decreased the global burden of the disease and enabled 187 countries to be certified by the World Health Organization to be free of endemic transmission. Despite this progress, several countries continue to have endemic transmission. In Chad, a long absence of reported cases was interrupted with the identification of new Guinea worm cases, prompting a substantial scale up of surveillance and intervention efforts. Here, we study the value of increasing genomic surveillance as a tool for programmatic evaluation of surveillance and intervention efforts in Chad. Linking surveillance and genomic samples, parsimonious spatial models help reveal a consistent geographic clustering of similar genetic sequences across Chad. We also demonstrate that expanding the sequencing can offer better resolution for distinguishing Guinea worm samples. In this retrospective study, we found evidence that scaling up genomic surveillance can be an important monitoring and evaluation tool for the eradication program in Chad.

Genome-Wide Characterization and Identification of Long Non-Coding RNAs during the Molting Process of a Spider Mite, Panonychus citri

Molting is essential for arthropods to grow. As one of the important arthropod pests in agriculture, key spider mite species (Tetranychus and Panonychus) can normally molt three times from the larva to adult stage within a week. This physiological strategy results in the short lifecycle of spider mites and difficulties in their control in the field. Long non-coding RNAs (lncRNAs) regulate transcriptional editing, cellular function, and biological processes. Thus, analysis of the lncRNAs in the spider mite molting process may provide new insights into their roles in the molting mechanism. For this purpose, we used high-throughput RNA-seq to examine the expression dynamics of lncRNAs and mRNAs in the molting process of different development stages in Panonychus citri. We identified 9199 lncRNAs from 18 transcriptomes. Analysis of the lncRNAs suggested that they were shorter and had fewer exons and transcripts than mRNAs. Among these, 356 lncRNAs were differentially expressed during three molting processes: late larva to early protonymph, late protonymph to early deutonymph, and late deutonymph to early adult. A time series profile analysis of differentially expressed lncRNAs showed that 77 lncRNAs were clustered into two dynamic expression profiles (Pattern a and Pattern c), implying that lncRNAs were involved in the molting process of spider mites. Furthermore, the lncRNA-mRNA co-expression networks showed that several differentially expressed hub lncRNAs were predicted to be functionally associated with typical molting-related proteins, such as cuticle protein and chitin biosynthesis. These data reveal the potential regulatory function of lncRNAs in the molting process and provide datasets for further analysis of lncRNAs and mRNAs in spider mites.

DYF-4 regulates patched-related/DAF-6-mediated sensory compartment formation in C. elegans

Coordination of neurite extension with surrounding glia development is critical for neuronal function, but the underlying molecular mechanisms remain poorly understood. Through a genome-wide mutagenesis screen in C. elegans, we identified dyf-4 and daf-6 as two mutants sharing similar defects in dendrite extension. DAF-6 encodes a glia-specific patched-related membrane protein that plays vital roles in glial morphogenesis. We cloned dyf-4 and found that DYF-4 encodes a glia-secreted protein. Further investigations revealed that DYF-4 interacts with DAF-6 and functions in a same pathway as DAF-6 to regulate sensory compartment formation. Furthermore, we demonstrated that reported glial suppressors of daf-6 could also restore dendrite elongation and ciliogenesis in both dyf-4 and daf-6 mutants. Collectively, our data reveal that DYF-4 is a regulator for DAF-6 which promotes the proper formation of the glial channel and indirectly affects neurite extension and ciliogenesis.

In C. elegans sensory organ, the ciliated neuronal endings are wrapped in a luminal channel formed by glial cells, forming a specialized sensory compartment critical for sensory activity. Coordination of glial channel formation, dendritic tip anchoring and ciliogenesis are critical for the formation of a functional sensory compartment. DAF-6, a patched-related glial membrane protein, was reported to play an important role in glial channel morphogenesis, but the molecular function and regulatory mechanism of DAF-6 remain poorly understood. Here, we found that DYF-4, a glia-secreted protein, interacts and colocalizes with DAF-6, and functions in a same pathway as DAF-6 to regulate sensory compartment formation. We propose that DYF-4 is a novel regulator for DAF-6 to control sensory compartment formation.

Linking genome size variation to population phenotypic variation within the rotifer, Brachionus asplanchnoidis

Eukaryotic organisms usually contain much more genomic DNA than expected from their biological complexity. In explaining this pattern, selection-based hypotheses suggest that genome size evolves through selection acting on correlated life history traits, implicitly assuming the existence of phenotypic effects of (extra) genomic DNA that are independent of its information content. Here, we present conclusive evidence of such phenotypic effects within a well-mixed natural population that shows heritable variation in genome size. We found that genome size is positively correlated with body size, egg size, and embryonic development time in a population of the monogonont rotifer Brachionus asplanchnoidis. The effect on embryonic development time was mediated partly by an indirect effect (via egg size), and a direct effect, the latter indicating an increased replication cost of the larger amounts of DNA during mitosis. Our results suggest that selection-based change of genome size can operate in this population, provided it is strong enough to overcome drift or mutational change of genome size.

Protection of the C. elegans germ cell genome depends on diverse DNA repair pathways during normal proliferation

Maintaining genome integrity is particularly important in germ cells to ensure faithful transmission of genetic information across generations. Here we systematically describe germ cell mutagenesis in wild-type and 61 DNA repair mutants cultivated over multiple generations. ~44% of the DNA repair mutants analysed showed a >2-fold increased mutagenesis with a broad spectrum of mutational outcomes. Nucleotide excision repair deficiency led to higher base substitution rates, whereas polh-1(Polη) and rev-3(Polζ) translesion synthesis polymerase mutants resulted in 50-400 bp deletions. Signatures associated with defective homologous recombination fall into two classes: 1) brc-1/BRCA1 and rad-51/RAD51 paralog mutants showed increased mutations across all mutation classes, 2) mus-81/MUS81 and slx-1/SLX1 nuclease, and him-6/BLM, helq-1/HELQ or rtel-1/RTEL1 helicase mutants primarily accumulated structural variants. Repetitive and G-quadruplex sequence-containing loci were more frequently mutated in specific DNA repair backgrounds. Tandem duplications embedded in inverted repeats were observed in helq-1 helicase mutants, and a unique pattern of 'translocations' involving homeologous sequences occurred in rip-1 recombination mutants. atm-1/ATM checkpoint mutants harboured structural variants specifically enriched in subtelomeric regions. Interestingly, locally clustered mutagenesis was only observed for combined brc-1 and cep-1/p53 deficiency. Our study provides a global view of how different DNA repair pathways contribute to prevent germ cell mutagenesis.

Parallel genetics of regulatory sequences using scalable genome editing in vivo

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

Identifying and validating the presence of Guanine-Quadruplexes (G4) within the blood fluke parasite Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that currently affects over 250 million individuals worldwide. In the absence of an immunoprophylactic vaccine and the recognition that mono-chemotherapeutic control of schistosomiasis by praziquantel has limitations, new strategies for managing disease burden are urgently needed. A better understanding of schistosome biology could identify previously undocumented areas suitable for the development of novel interventions. Here, for the first time, we detail the presence of G-quadruplexes (G4) and putative quadruplex forming sequences (PQS) within the Schistosoma mansoni genome. We find that G4 are present in both intragenic and intergenic regions of the seven autosomes as well as the sex-defining allosome pair. Amongst intragenic regions, G4 are particularly enriched in 3´ UTR regions. Gene Ontology (GO) term analysis evidenced significant G4 enrichment in the wnt signalling pathway (p<0.05) and PQS oligonucleotides synthetically derived from wnt-related genes resolve into parallel and anti-parallel G4 motifs as elucidated by circular dichroism (CD) spectroscopy. Finally, utilising a single chain anti-G4 antibody called BG4, we confirm the in situ presence of G4 within both adult female and male worm nuclei. These results collectively suggest that G4-targeted compounds could be tested as novel anthelmintic agents and highlights the possibility that G4-stabilizing molecules could be progressed as candidates for the treatment of schistosomiasis.

The Genetics of Human Schistosomiasis Infection Intensity and Liver Disease: A Review

Schistosomiasis remains the fourth most prevalent parasitic disease affecting over 200 million people worldwide. Control efforts have focussed on the disruption of the life cycle targeting the parasite, vector and human host. Parasite burdens are highly skewed, and the majority of eggs are shed into the environment by a minority of the infected population. Most morbidity results from hepatic fibrosis leading to portal hypertension and is not well-correlated with worm burden. Genetics as well as environmental factors may play a role in these skewed distributions and understanding the genetic risk factors for intensity of infection and morbidity may help improve control measures. In this review, we focus on how genetic factors may influence parasite load, hepatic fibrosis and portal hypertension. We found 28 studies on the genetics of human infection and 20 studies on the genetics of pathology in humans. S. mansoni and S. haematobium infection intensity have been showed to be controlled by a major quantitative trait locus SM1, on chromosome 5q31-q33 containing several genes involved in the Th2 immune response, and three other loci of smaller effect on chromosomes 1, 6, and 7. The most common pathology associated with schistosomiasis is hepatic and portal vein fibroses and the SM2 quantitative trait locus on chromosome six has been linked to intensity of fibrosis. Although there has been an emphasis on Th2 cytokines in candidate gene studies, we found that four of the five QTL regions contain Th17 pathway genes that have been included in schistosomiasis studies: IL17B and IL12B in SM1, IL17A and IL17F in 6p21-q2, IL6R in 1p21-q23 and IL22RA2 in SM2. The Th17 pathway is known to be involved in response to schistosome infection and hepatic fibrosis but variants in this pathway have not been tested for any effect on the regulation of these phenotypes. These should be priorities for future studies.

A first look at the 'repeatome' of Benedenia humboldti, a major pathogen in yellowtail aquaculture: Repetitive element characterization, nuclear rRNA operon assembly, and microsatellite discovery

The monogenean Benedenia humboldti is a pathogen of the yellowtail Seriola lalandi in the South-Eastern Pacific ocean. Using low-coverage short Illumina 150bp pair-end reads sequencing, this study examines, for the first time, the 'repeatome' (= repetitive genomic elements), including the 45S ribosomal RNA DNA operon and microsatellites, in B. humboldti. Repetitive elements comprised a large fraction of the nuclear genome and a considerable proportion of them could not be assigned to known repeat element families. Taking into account only annotated repetitive elements, the most frequent belonged to the 45S ribosomal RNA operon or were classified as satellite DNA and Class I - Long Interspersed Nuclear Elements (LINEs) which were considerably more abundant than Class I - LTR elements. The ribosomal RNA gene operon in B. humboldti is comprised of, in the following order, a 5' ETS (length = 233 bp), ssrDNA (2082 bp), ITS1 (346 bp), 5.8S rDNA (150 bp), ITS2 (572 bp), lsrDNA (3887 bp), and a 3' ETS (1097 bp). A total of 15 SSRs were identified. These newly developed genomic resources will contribute to the better understanding of meta-population connectivity in this species, cryptic species in the genus, and will advance pest management strategies.

Resolution of polycistronic RNA by SL2 trans-splicing is a widely conserved nematode trait

Spliced leader trans-splicing is essential for the processing and translation of polycistronic RNAs generated by eukaryotic operons. In C. elegans, a specialized spliced leader, SL2, provides the 5' end for uncapped pre-mRNAs derived from polycistronic RNAs. Studies of other nematodes suggested that SL2-type trans-splicing is a relatively recent innovation, confined to Rhabditina, the clade containing C. elegans and its close relatives. Here we conduct a survey of transcriptome-wide spliced leader trans-splicing in Trichinella spiralis, a distant relative of C. elegans with a particularly diverse repertoire of 15 spliced leaders. By systematically comparing the genomic context of trans-splicing events for each spliced leader, we identified a subset of T. spiralis spliced leaders that are specifically used to process polycistronic RNAs-the first examples of SL2-type spliced leaders outside of Rhabditina. These T. spiralis spliced leader RNAs possess a perfectly conserved stem-loop motif previously shown to be essential for SL2-type trans-splicing in C. elegans We show that genes trans-spliced to these SL2-type spliced leaders are organized in operonic fashion, with short intercistronic distances. A subset of T. spiralis operons show conservation of synteny with C. elegans operons. Our work substantially revises our understanding of nematode spliced leader trans-splicing, showing that SL2 trans-splicing is a major mechanism for nematode polycistronic RNA processing, which may have evolved prior to the radiation of the Nematoda. This work has important implications for the improvement of genome annotation pipelines in nematodes and other eukaryotes with operonic gene organization.

Population genomic evidence that human and animal infections in Africa come from the same populations of Dracunculus medinensis

BACKGROUND

Guinea worm-Dracunculus medinensis-was historically one of the major parasites of humans and has been known since antiquity. Now, Guinea worm is on the brink of eradication, as efforts to interrupt transmission have reduced the annual burden of disease from millions of infections per year in the 1980s to only 54 human cases reported globally in 2019. Despite the enormous success of eradication efforts to date, one complication has arisen. Over the last few years, hundreds of dogs have been found infected with this previously apparently anthroponotic parasite, almost all in Chad. Moreover, the relative numbers of infections in humans and dogs suggests that dogs are currently the principal reservoir on infection and key to maintaining transmission in that country.

PRINCIPAL FINDINGS

In an effort to shed light on this peculiar epidemiology of Guinea worm in Chad, we have sequenced and compared the genomes of worms from dog, human and other animal infections. Confirming previous work with other molecular markers, we show that all of these worms are D. medinensis, and that the same population of worms are causing both infections, can confirm the suspected transmission between host species and detect signs of a population bottleneck due to the eradication efforts. The diversity of worms in Chad appears to exclude the possibility that there were no, or very few, worms present in the country during a 10-year absence of reported cases.

CONCLUSIONS

This work reinforces the importance of adequate surveillance of both human and dog populations in the Guinea worm eradication campaign and suggests that control programs aiming to interrupt disease transmission should stay aware of the possible emergence of unusual epidemiology as pathogens approach elimination.

Complete representation of a tapeworm genome reveals chromosomes capped by centromeres, necessitating a dual role in segregation and protection

BACKGROUND

Chromosome-level assemblies are indispensable for accurate gene prediction, synteny assessment, and understanding higher-order genome architecture. Reference and draft genomes of key helminth species have been published, but little is yet known about the biology of their chromosomes. Here, we present the complete genome of the tapeworm Hymenolepis microstoma, providing a reference quality, end-to-end assembly that represents the first fully assembled genome of a spiralian/lophotrochozoan, revealing new insights into chromosome evolution.

RESULTS

Long-read sequencing and optical mapping data were added to previous short-read data enabling complete re-assembly into six chromosomes, consistent with karyology. Small genome size (169 Mb) and lack of haploid variation (1 SNP/3.2 Mb) contributed to exceptionally high contiguity with only 85 gaps remaining in regions of low complexity sequence. Resolution of repeat regions reveals novel gene expansions, micro-exon genes, and spliced leader trans-splicing, and illuminates the landscape of transposable elements, explaining observed length differences in sister chromatids. Syntenic comparison with other parasitic flatworms shows conserved ancestral linkage groups indicating that the H. microstoma karyotype evolved through fusion events. Strikingly, the assembly reveals that the chromosomes terminate in centromeric arrays, indicating that these motifs play a role not only in segregation, but also in protecting the linear integrity and full lengths of chromosomes.

CONCLUSIONS

Despite strong conservation of canonical telomeres, our results show that they can be substituted by more complex, species-specific sequences, as represented by centromeres. The assembly provides a robust platform for investigations that require complete genome representation.

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.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

Schistosome TRP channels: An appraisal

Ion channels underlie electrical excitability in cells and are essential for a variety of functions, most notably neuromuscular and sensory activity. They are also validated targets for a preponderance of approved anthelmintic compounds. Transient receptor potential (TRP) channels constitute an ion channel superfamily whose members play important roles in sensory signaling, regulation of ion homeostasis, organellar trafficking, and other key cellular and organismal activities. Unlike most other ion channels, TRP channels are often polymodal, gated by a variety of mechanisms. Furthermore, TRP channels fall into several classes or subtypes based on sequence and structure. Until recently, there had been very little investigation of the properties and functions of TRP channels from parasitic helminths, including schistosomes, but that situation has changed in the past few years. Indeed, it is now clear that at least some schistosome TRP channels exhibit unusual pharmacological properties, and, intriguingly, both a mammalian and a schistosome TRP channel are activated by praziquantel, the current antischistosomal drug of choice. With the latest release of the Schistosoma mansoni genome database, several changes in predicted TRP channel sequences appeared, some of which were significant. This review updates and reassesses the TRP channel repertoire in S. mansoni, examines recent findings regarding these potential therapeutic targets, and provides guideposts for some of the physiological functions that may be mediated by these channels in schistosomes.

Regional sequence expansion or collapse in heterozygous genome assemblies

High levels of heterozygosity present a unique genome assembly challenge and can adversely impact downstream analyses, yet is common in sequencing datasets obtained from non-model organisms. Here we show that by re-assembling a heterozygous dataset with variant parameters and different assembly algorithms, we are able to generate assemblies whose protein annotations are statistically enriched for specific gene ontology categories. While total assembly length was not significantly affected by assembly methodologies tested, the assemblies generated varied widely in fragmentation level and we show local assembly collapse or expansion underlying the enrichment or depletion of specific protein functional groups. We show that these statistically significant deviations in gene ontology groups can occur in seemingly high-quality assemblies, and result from difficult-to-detect local sequence expansion or contractions. Given the unpredictable interplay between assembly algorithm, parameter, and biological sequence data heterozygosity, we highlight the need for better measures of assembly quality than N50 value, including methods for assessing local expansion and collapse.

In the genomic era, genomes must be reconstructed from fragments using computational methods, or assemblers. How do we know that a new genome assembly is correct? This is important because errors in assembly can lead to downstream problems in gene predictions and these inaccurate results can contaminate databases, affecting later comparative studies. A particular challenge occurs when a diploid organism inherits two highly divergent genome copies from its parents. While it is widely appreciated that this type of data is difficult for assemblers to handle properly, here we show that the process is prone to more errors than previously appreciated. Specifically, we document examples of regional expansion and collapse, affecting downstream gene prediction accuracy, but without changing the overall genome assembly size or other metrics of accuracy. Our results suggest that assembly evaluation methods should be altered to identify whether regional expansions and collapses are present in the genome assembly.

Cloning and functional complementation of ten Schistosoma mansoni phosphodiesterases expressed in the mammalian host stages

Only a single drug against schistosomiasis is currently available and new drug development is urgently required but very few drug targets have been validated and characterised. However, regulatory systems including cyclic nucleotide metabolism are emerging as primary candidates for drug discovery. Here, we report the cloning of ten cyclic nucleotide phosphodiesterase (PDE) genes of S. mansoni, out of a total of 11 identified in its genome. We classify these PDEs by homology to human PDEs. Male worms displayed higher expression levels for all PDEs, in mature and juvenile worms, and schistosomula. Several functional complementation approaches were used to characterise these genes. We constructed a Trypanosoma brucei cell line in which expression of a cAMP-degrading PDE complements the deletion of TbrPDEB1/B2. Inhibitor screens of these cells expressing only either SmPDE4A, TbrPDEB1 or TbrPDEB2, identified highly potent inhibitors of the S. mansoni enzyme that elevated the cellular cAMP concentration. We further expressed most of the cloned SmPDEs in two pde1Δ/pde2Δ strains of Saccharomyces cerevisiae and some also in a specialised strain of Schizosacharomyces pombe. Five PDEs, SmPDE1, SmPDE4A, SmPDE8, SmPDE9A and SmPDE11 successfully complemented the S. cerevisiae strains, and SmPDE7var also complemented to a lesser degree, in liquid culture. SmPDE4A, SmPDE8 and SmPDE11 were further assessed in S. pombe for hydrolysis of cAMP and cGMP; SmPDE11 displayed considerable preferrence for cGMP over cAMP. These results and tools enable the pursuit of a rigorous drug discovery program based on inhibitors of S. mansoni PDEs.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Comprehensive analysis of the secreted proteome of adult Necator americanus hookworms

The human hookworm Necator americanus infects more than 400 million people worldwide, contributing substantially to the poverty in these regions. Adult stage N. americanus live in the small intestine of the human host where they inject excretory/secretory (ES) products into the mucosa. ES products have been characterized at the proteome level for a number of animal hookworm species, but until now, the difficulty in obtaining sufficient live N. americanus has been an obstacle in characterizing the secretome of this important human pathogen. Herein we describe the ES proteome of N. americanus and utilize this information along with RNA Seq data to conduct the first proteogenomic analysis of a parasitic helminth, significantly improving the available genome and thereby generating a robust description of the parasite secretome. The genome annotation resulted in a revised prediction of 3,425 fewer genes than initially reported, accompanied by a significant increase in the number of exons and introns, total gene length and the percentage of the genome covered by genes. Almost 200 ES proteins were identified by LC-MS/MS with SCP/TAPS proteins, ‘hypothetical’ proteins and proteases among the most abundant families. These proteins were compared to commonly used model species of human parasitic infections, including Ancylostoma caninum, Nippostrongylus brasiliensis and Heligmosomoides polygyrus. SCP/TAPS proteins are immunogenic in nematode infections, so we expressed four of those identified in this study in recombinant form and showed that they are all recognized to varying degrees by serum antibodies from hookworm-infected subjects from a disease-endemic area of Brazil. Our findings provide valuable information on important families of proteins with both known and unknown functions that could be instrumental in host-parasite interactions, including protein families that might be key for parasite survival in the onslaught of robust immune responses, as well as vaccine and diagnostic targets.

Hookworms infect hundreds of millions of people in tropical regions of the world. Adult worms reside in the small bowel where they feed on blood, causing iron-deficiency anemia when present in large numbers and contributing substantially to the poverty in these regions. Hookworms inject excretory/secretory (ES) products into the gut tissue when they feed, and while the protein constituents of ES products have been characterized for a number of animal hookworm species, difficulty in obtaining sufficient live human hookworms has thus far precluded characterization of the secreted proteome. Herein we describe the ES proteins of the major human hookworm, Necator americanus, and utilize this information to significantly improve the available genome sequence. Almost 200 ES proteins were identified and compared to the secreted proteomes of other parasitic roundworms to provide a molecular snapshot of the host-parasite interface. We produced recombinant forms of some of the identified proteins and showed that they are all recognized to varying degrees by antibodies from hookworm-infected subjects. Our work sheds light on important families of proteins that might be key for parasite survival in the human host, and presents a dataset that can now be mined in the search for vaccine, drug and diagnostic targets.

Screening for CRISPR/Cas9-induced mutations using a co-injection marker in the nematode Pristionchus pacificus

CRISPR/Cas9 genome-editing methods are used to reveal functions of genes and molecular mechanisms underlying biological processes in many species, including nematodes. In evolutionary biology, the nematode Pristionchus pacificus is a satellite model and has been used to understand interesting phenomena such as phenotypic plasticity and self-recognition. In P. pacificus, CRISPR/Cas9-mediated mutations are induced by microinjecting a guide RNA (gRNA) and Cas9 protein into the gonads. However, mutant screening is laborious and time-consuming due to the absence of visual markers. In this study, we established a Co-CRISPR strategy by using a dominant roller marker in P. pacificus. We found that heterozygous mutations in Ppa-prl-1 induced the roller phenotype, which can be used as an injection marker. After the co-injection of Ppa-prl-1 gRNA, target gRNA, and the Cas9 protein, roller progeny and their siblings were examined using the heteroduplex mobility assay and DNA sequencing. We found that some of the roller and non-roller siblings had mutations at the target site. We used varying Cas9 concentrations and found that a higher concentration of Cas9 did not increase genome-editing events. The Co-CRISPR strategy promotes the screening for genome-editing events and will facilitate the development of new genome-editing methods in P. pacificus.

Chromosome-Level Assembly of the Caenorhabditis remanei Genome Reveals Conserved Patterns of Nematode Genome Organization

The nematode Caenorhabditis elegans is one of the key model systems in biology, including possessing the first fully assembled animal genome. Whereas C. elegans is a self-reproducing hermaphrodite with fairly limited within-population variation, its relative C. remanei is an outcrossing species with much more extensive genetic variation, making it an ideal parallel model system for evolutionary genetic investigations. Here, we greatly improve on previous assemblies by generating a chromosome-level assembly of the entire C. remanei genome (124.8 Mb of total size) using long-read sequencing and chromatin conformation capture data. Like other fully assembled genomes in the genus, we find that the C. remanei genome displays a high degree of synteny with C. elegans despite multiple within-chromosome rearrangements. Both genomes have high gene density in central regions of chromosomes relative to chromosome ends and the opposite pattern for the accumulation of repetitive elements. C. elegans and C. remanei also show similar patterns of interchromosome interactions, with the central regions of chromosomes appearing to interact with one another more than the distal ends. The new C. remanei genome presented here greatly augments the use of the Caenorhabditis as a platform for comparative genomics and serves as a basis for molecular population genetics within this highly diverse species.

Comparative Transcriptomics Reveals Clues for Differences in Pathogenicity between Hysterothylacium aduncum, Anisakis simplex sensu stricto and Anisakis pegreffii

Ascaridoid nematodes are widespread in marine fishes. Despite their major socioeconomic importance, mechanisms associated to the fish-borne zoonotic disease anisakiasis are still obscure. RNA-Seq and de-novo assembly were herein applied to RNA extracted from larvae and dissected pharynx of Hysterothylacium aduncum (HA), a non-pathogenic nematode. Assembled transcripts in HA were annotated and compared to the transcriptomes of the zoonotic species Anisakis simplex sensu stricto (AS) and Anisakis pegreffii (AP). Approximately 60,000,000 single-end reads were generated for HA, AS and AP. Transcripts in HA encoded for 30,254 putative peptides while AS and AP encoded for 20,574 and 20,840 putative peptides, respectively. Differential gene expression analyses yielded 471, 612 and 526 transcripts up regulated in the pharynx of HA, AS and AP. The transcriptomes of larvae and pharynx of HA were enriched in transcripts encoding collagen, peptidases, ribosomal proteins and in heat-shock motifs. Transcripts encoding proteolytic enzymes, anesthetics, inhibitors of primary hemostasis and virulence factors, anticoagulants and immunomodulatory peptides were up-regulated in AS and AP pharynx. This study represents the first transcriptomic characterization of a marine parasitic nematode commonly recovered in fish and probably of negligible concern for public health.

Point-of-care diagnostic (POCD) method for detecting Bursaphelenchus xylophilus in pinewood using recombinase polymerase amplification (RPA) with the portable optical isothermal device (POID)

The pinewood nematode (PWN), Bursaphelenchus xylophilus, is a causative agent of pine wilt disease (PWD). To date, although several molecular diagnostic methods have been developed, rapid on-site diagnostic tools for detecting PWN in pinewood are limited. In this study, a point of care diagnostic (POCD) method for detecting PWN in pinewood using recombinase polymerase amplification (RPA) assay was developed. This method comprises quick gDNA extraction buffer (DAP buffer) for the direct extraction of gDNA of PWN from pinewood and a battery-mounted portable optical isothermal device (POID) for the detection of PWD in the field. The RPA assay can distinguish between the PWN and its conspecies which exist in pinewood and can complete diagnostic procedures within 25 min in the field. Moreover, the RPA assay can detect PWN in old wood samples in both natural and stored conditions. The POCD-RPA assay to detect PWN will be useful for epidemiological investigations in the field as well as for quarantine processes in the wood trade.

Hybrid de novo whole-genome assembly and annotation of the model tapeworm Hymenolepis diminuta

Despite the use of Hymenolepis diminuta as a model organism in experimental parasitology, a full genome description has not yet been published. Here we present a hybrid de novo genome assembly based on complementary sequencing technologies and methods. The combination of Illumina paired-end, Illumina mate-pair and Oxford Nanopore Technology reads greatly improved the assembly of the H. diminuta genome. Our results indicate that the hybrid sequencing approach is the method of choice for obtaining high-quality data. The final genome assembly is 177 Mbp with contig N50 size of 75 kbp and a scaffold N50 size of 2.3 Mbp. We obtained one of the most complete cestode genome assemblies and annotated 15,169 potential protein-coding genes. The obtained data may help explain cestode gene function and better clarify the evolution of its gene families, and thus the adaptive features evolved during millennia of co-evolution with their hosts.

Evidence for densovirus integrations into tapeworm genomes

BACKGROUND

Tapeworms lack a canonical piRNA-pathway, raising the question of how they can silence existing mobile genetic elements (MGE). Investigation towards the underlying mechanisms requires information on tapeworm transposons which is, however, presently scarce.

METHODS

The presence of densovirus-related sequences in tapeworm genomes was studied by bioinformatic approaches. Available RNA-Seq datasets were mapped against the Echinococcus multilocularis genome to calculate expression levels of densovirus-related genes. Transcription of densovirus loci was further analyzed by sequencing and RT-qPCR.

RESULTS

We herein provide evidence for the presence of densovirus-related elements in a variety of tapeworm genomes. In the high-quality genome of E. multilocularis we identified more than 20 individual densovirus integration loci which contain the information for non-structural and structural virus proteins. The majority of densovirus loci are present as head-to-tail concatemers in isolated repeat containing regions of the genome. In some cases, unique densovirus loci have integrated close to histone gene clusters. We show that some of the densovirus loci of E. multilocularis are actively transcribed, whereas the majority are transcriptionally silent. RT-qPCR data further indicate that densovirus expression mainly occurs in the E. multilocularis stem cell population, which probably forms the germline of this organism. Sequences similar to the non-structural densovirus genes present in E. multilocularis were also identified in the genomes of E. canadensis, E. granulosus, Hydatigera taeniaeformis, Hymenolepis diminuta, Hymenolepis microstoma, Hymenolepis nana, Taenia asiatica, Taenia multiceps, Taenia saginata and Taenia solium.

CONCLUSIONS

Our data indicate that densovirus integration has occurred in many tapeworm species. This is the first report on widespread integration of DNA viruses into cestode genomes. Since only few densovirus integration sites were transcriptionally active in E. multilocularis, our data are relevant for future studies into gene silencing mechanisms in tapeworms. Furthermore, they indicate that densovirus-based vectors might be suitable tools for genetic manipulation of cestodes.

A Chemosensory GPCR as a Potential Target to Control the Root-Knot Nematode Meloidogyne incognita Parasitism in Plants

Root-knot nematodes (RKN), from the Meloidogyne genus, have a worldwide distribution and cause severe economic damage to many life-sustaining crops. Because of their lack of specificity and danger to the environment, most chemical nematicides have been banned from use. Thus, there is a great need for new and safe compounds to control RKN. Such research involves identifying beforehand the nematode proteins essential to the invasion. Since G protein-coupled receptors GPCRs are the target of a large number of drugs, we have focused our research on the identification of putative nematode GPCRs such as those capable of controlling the movement of the parasite towards (or within) its host. A datamining procedure applied to the genome of Meloidogyne incognita allowed us to identify a GPCR, belonging to the neuropeptide GPCR family that can serve as a target to carry out a virtual screening campaign. We reconstructed a 3D model of this receptor by homology modeling and validated it through extensive molecular dynamics simulations. This model was used for large scale molecular dockings which produced a filtered limited set of putative antagonists for this GPCR. Preliminary experiments using these selected molecules allowed the identification of an active compound, namely C260-2124, from the ChemDiv provider, which can serve as a starting point for further investigations.

The genome of the marine monogonont rotifer Brachionus plicatilis: Genome-wide expression profiles of 28 cytochrome P450 genes in response to chlorpyrifos and 2-ethyl-phenanthrene

Brachionus spp. (Rotifera: Monogononta) are globally distributed in aquatic environments and play important roles in the aquatic ecosystem. The marine monogonont rotifer Brachionus plicatilis is considered a suitable model organism for ecology, evolution, and ecotoxicology. In this study, we assembled and characterized the B. plicatilis genome. The total length of the assembled genome was 106.9 Mb and the number of final scaffolds was 716 with an N50 value of 1.15 Mb and a GC content of 26.75%. A total of 20,154 genes were annotated after manual curation. To demonstrate the use of whole genome data, we targeted one of the main detoxifying enzyme of phase I detoxification system and identified in a total of 28 cytochrome P450 s (CYPs). Based on the phylogenetic analysis using the maximum likelihood, 28 B. plicatilis-CYPs were apparently separated into five different clans, namely, 2, 3, 4, mitochondrial (MT), and 46 clans. To better understand the CYPs-mediated xenobiotic detoxification, we measured the mRNA expression levels of 28 B. plicatilis CYPs in response to chlorpyrifos and 2-ethyl-phenanthrene. Most B. plicatilis CYPs were significantly modulated (P < 0.05) in response to chlorpyrifos and 2-ethyl-phenanthrene. In addition, xenobiotic-sensing nuclear receptor (XNR) response element sequences were identified in the 5 kb upstream of promoter regions of 28 CYPs from the genome of B. plicatilis, indicating that these XNR can be associated with detoxification of xenobiotics. Overall, the assembled B. plicatilis genome presented here will be a useful resource for a better understanding the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification in this species.

High-quality Schistosoma haematobium genome achieved by single-molecule and long-range sequencing

BACKGROUND

Schistosoma haematobium causes urogenital schistosomiasis, a neglected tropical disease affecting >100 million people worldwide. Chronic infection with this parasitic trematode can lead to urogenital conditions including female genital schistosomiasis and bladder cancer. At the molecular level, little is known about this blood fluke and the pathogenesis of the disease that it causes. To support molecular studies of this carcinogenic worm, we reported a draft genome for S. haematobium in 2012. Although a useful resource, its utility has been somewhat limited by its fragmentation.

FINDINGS

Here, we systematically enhanced the draft genome of S. haematobium using a single-molecule and long-range DNA-sequencing approach. We achieved a major improvement in the accuracy and contiguity of the genome assembly, making it superior or comparable to assemblies for other schistosome species. We transferred curated gene models to this assembly and, using enhanced gene annotation pipelines, inferred a gene set with as many or more complete gene models as those of other well-studied schistosomes. Using conserved, single-copy orthologs, we assessed the phylogenetic position of S. haematobium in relation to other parasitic flatworms for which draft genomes were available.

CONCLUSIONS

We report a substantially enhanced genomic resource that represents a solid foundation for molecular research on S. haematobium and is poised to better underpin population and functional genomic investigations and to accelerate the search for new disease interventions.

Mitochondrial Genomes of Two Thaparocleidus Species (Platyhelminthes: Monogenea) Reveal the First rRNA Gene Rearrangement among the Neodermata

Phylogenetic framework for the closely related Ancylodiscoidinae and Ancyrocephalinae subfamilies remains contentious. As this issue was never studied using a large molecular marker, we sequenced the first two Ancylodiscoidinae mitogenomes: Thaparocleidus asoti and Thaparocleidus varicus. Both mitogenomes had two non-coding regions (NCRs) that contained a number of repetitive hairpin-forming elements (RHE). Due to these, the mitogenome of T. asoti (16,074 bp) is the longest among the Monogenea; especially large is its major NCR, with 3500 bp, approximately 1500 bp of which could not be sequenced (thus, the total mitogenome size is ≈ 17,600 bp). Although RHEs have been identified in other monopisthocotyleans, they appear to be independently derived in different taxa. The presence of RHEs may have contributed to the high gene order rearrangement rate observed in the two mitogenomes, including the first report of a transposition of rRNA genes within the Neodermata. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrinae embedded within the Ancyrocephalinae (paraphyly), whereas Ancylodiscoidinae formed a sister-group with them. This was also supported by the gene order analysis. 28S rDNA dataset produced polyphyletic Dactylogyridae and Ancyrocephalinae. The phylogeny of the two subfamilies shall have to be further evaluated with more data.

The Helitron family classification using SVM based on Fourier transform features applied on an unbalanced dataset

Helitrons are mobile sequences which belong to the class 2 of eukaryotic transposons. Their specificity resides in their mechanism of transposition: the rolling circle mechanism. They play an important role in remodeling proteomes due to their ability to modify existing genes and introducing new ones. A major difficulty in identifying and classifying Helitron families comes from the complex structure, the unspecified length, and the unbalanced appearance number of each Helitron type. The Helitron's recognition is still not solved in literature. The purpose of this paper is to characterize and classify Helitron types using spectral features and support vector machine (SVM) classification technique. Thus, the helitronic DNA is transformed into a numerical form using the FCGS₂ coding technique. Then, a set of spectral features is extracted from the smoothed Fourier transform applied on the FCGS₂ signals. Based on the spectral signature and the classification's confusion matrix, we demonstrated that some specific classes which do not show similarities, such as HelitronY2 and NDNAX3, are easily discriminated with important accuracy rates exceeding 90%. However, some Helitron types have great similarities such as the following: Helitron1, HelitronY1, HelitronY1A, and HelitronY4. Our system is also able to predict them with promising values reaching 70%. Graphical abstract The Helitron recognizer based on features extracted from smoothed Fourier transform.

An improved genome assembly of the fluke Schistosoma japonicum

BACKGROUND

Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China and the Philippines. A single draft genome was available for S. japonicum, yet this assembly is very fragmented and only covers 90% of the genome, which make it difficult to be applied as a reference in functional genome analysis and genes discovery.

FINDINGS

In this study, we present a high-quality assembly of the fluke S. japonicum genome by combining 20 G (~53X) long single molecule real time sequencing reads with 80 G (~ 213X) Illumina paired-end reads. This improved genome assembly is approximately 370.5 Mb, with contig and scaffold N50 length of 871.9 kb and 1.09 Mb, representing 142.4-fold and 6.2-fold improvement over the released WGS-based assembly, respectively. Additionally, our assembly captured 85.2% complete and 4.6% partial eukaryotic Benchmarking Universal Single-Copy Orthologs. Repetitive elements account for 46.80% of the genome, and 10,089 of the protein-coding genes were predicted from the improved genome, of which 96.5% have been functionally annotated. Lastly, using the improved assembly, we identified 20 significantly expanded gene families in S. japonicum, and those genes were primarily enriched in functions of proteolysis and protein glycosylation.

CONCLUSIONS

Using the combination of PacBio and Illumina Sequencing technologies, we provided an improved high-quality genome of S. japonicum. This improved genome assembly, as well as the annotation, will be useful for the comparative genomics of the flukes and more importantly facilitate the molecular studies of this important parasite in the future.

Small, but surprisingly repetitive genomes: transposon expansion and not polyploidy has driven a doubling in genome size in a metazoan species complex

BACKGROUND

The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110-422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content).

RESULTS

Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers.

CONCLUSIONS

Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy.

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.