Functional analysis of putative operons in Brugia malayi

Long-read RNA sequencing of human and animal filarial parasites improves gene models and discovers operons

Filarial parasitic nematodes (Filarioidea) cause substantial disease burden to humans and animals around the world. Recently there has been a coordinated global effort to generate, annotate, and curate genomic data from nematode species of medical and veterinary importance. This has resulted in two chromosome-level assemblies (Brugia malayi and Onchocerca volvulus) and 11 additional draft genomes from Filarioidea. These reference assemblies facilitate comparative genomics to explore basic helminth biology and prioritize new drug and vaccine targets. While the continual improvement of genome contiguity and completeness advances these goals, experimental functional annotation of genes is often hindered by poor gene models. Short-read RNA sequencing data and expressed sequence tags, in cooperation with ab initio prediction algorithms, are employed for gene prediction, but these can result in missing clade-specific genes, fragmented models, imperfect mapping of gene ends, and lack of isoform resolution. Long-read RNA sequencing can overcome these drawbacks and greatly improve gene model quality. Here, we present Iso-Seq data for B. malayi and Dirofilaria immitis, etiological agents of lymphatic filariasis and canine heartworm disease, respectively. These data cover approximately half of the known coding genomes and substantially improve gene models by extending untranslated regions, cataloging novel splice junctions from novel isoforms, and correcting mispredicted junctions. Furthermore, we validated computationally predicted operons, manually curated new operons, and merged fragmented gene models. We carried out analyses of poly(A) tails in both species, leading to the identification of non-canonical poly(A) signals. Finally, we prioritized and assessed known and putative anthelmintic targets, correcting or validating gene models for molecular cloning and target-based anthelmintic screening efforts. Overall, these data significantly improve the catalog of gene models for two important parasites, and they demonstrate how long-read RNA sequencing should be prioritized for ongoing improvement of parasitic nematode genome assemblies.

Filarial parasitic nematodes are vector-borne parasites that infect humans and animals. Brugia malayi and Dirofilaria immitis are transmitted by mosquitoes and cause human lymphatic filariasis and canine heartworm disease, respectively. Recent years have seen a dramatic increase in genomic and transcriptomic data sets and the concomitant increase in innovative strategies for drug target identification, validation, and screening. However, while the completeness of genome assemblies of filarial parasitic nematodes has seen steady improvements, the reliability of gene models has not kept pace, hindering cloning efforts. Long-read RNA sequencing technologies are uniquely able to improve gene models, but have not been widely used for the causative agents of neglected tropical diseases. Here, we report the improvement of gene models in both B. malayi and D. immitis by long-read RNA sequencing. We identified novel operons, deprecated false positive operons, identified dozens of novel genes, and described the parameters of polyadenylation. We also focused on putative anthelmintic targets, identifying novel isoforms and correcting gene models. These data substantially increase the trustworthiness of gene models in these two species and demonstrate how long-read sequencing approaches should be prioritized in the continued improvement of genome assemblies and their gene annotations.

CRISPR/Cas9 Mutagenesis and Expression of Dominant Mutant Transgenes as Functional Genomic Approaches in Parasitic Nematodes

DNA transformation of parasitic nematodes enables novel approaches to validating predictions from genomic and transcriptomic studies of these important pathogens. Notably, proof of principle for CRISPR/Cas9 mutagenesis has been achieved in Strongyloides spp., allowing identification of molecules essential to the functions of sensory neurons that mediate behaviors comprising host finding, invasion, and location of predilection sites by parasitic nematodes. Likewise, CRISPR/Cas9 knockout of the developmental regulatory transcription factor Ss-daf-16 has validated its function in regulating morphogenesis of infective third-stage larvae in Strongyloides stercoralis. While encouraging, these studies underscore challenges that remain in achieving straightforward validation of essential intervention targets in parasitic nematodes. Chief among these is the likelihood that knockout of multifunctional regulators like Ss-DAF-16 or its downstream mediator, the nuclear receptor Ss-DAF-12, will produce phenotypes so complex as to defy interpretation and will render affected worms incapable of infecting their hosts, thus preventing establishment of stable mutant lines. Approaches to overcoming these impediments could involve refinements to current CRISPR/Cas9 methods in Strongyloides including regulatable Cas9 expression from integrated transgenes and CRISPR/Cas9 editing to ablate specific functional motifs in regulatory molecules without complete knockout. Another approach would express transgenes encoding regulatory molecules of interest with mutations designed to similarly ablate or degrade specific functional motifs such as the ligand binding domain of Ss-DAF-12 while preserving core functions such as DNA binding. Such mutant transgenes would be expected to exert a dominant interfering effect on their endogenous counterparts. Published reports validate the utility of such dominant-negative approaches in Strongyloides.

Rendering the Intractable More Tractable: Tools from Caenorhabditis elegans Ripe for Import into Parasitic Nematodes

Recent and rapid advances in genetic and molecular tools have brought spectacular tractability to Caenorhabditis elegans, a model that was initially prized because of its simple design and ease of imaging. C. elegans has long been a powerful model in biomedical research, and tools such as RNAi and the CRISPR/Cas9 system allow facile knockdown of genes and genome editing, respectively. These developments have created an additional opportunity to tackle one of the most debilitating burdens on global health and food security: parasitic nematodes. I review how development of nonparasitic nematodes as genetic models informs efforts to import tools into parasitic nematodes. Current tools in three commonly studied parasites (Strongyloides spp., Brugia malayi, and Ascaris suum) are described, as are tools from C. elegans that are ripe for adaptation and the benefits and barriers to doing so. These tools will enable dissection of a huge array of questions that have been all but completely impenetrable to date, allowing investigation into host-parasite and parasite-vector interactions, and the genetic basis of parasitism.

Operons are a conserved feature of nematode genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and "spliced leader" (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Biolistic transformation of Brugia malayi

Biolistics has become a versatile tool for direct gene transfer to various cell and tissue types. Following its successful use on the parasitic nematode Ascaris suum, we developed and evaluated biolistics in the transfection of the model filarial parasite Brugia malayi. Biolistics was proven to be an efficient strategy for transfection of all life stages of the parasite and paved the way for studies on elements essential for promoter function and gene regulation of filarial parasites. Here we present a biolistics protocol for the transfection of B. malayi based on the Biolistics PDS 1000/He system and gold microcarriers.

Functional genomics approaches in parasitic helminths

As research on parasitic helminths is moving into the post-genomic era, an enormous effort is directed towards deciphering gene function and to achieve gene annotation. The sequences that are available in public databases undoubtedly hold information that can be utilized for new interventions and control but the exploitation of these resources has until recently remained difficult. Only now, with the emergence of methods to genetically manipulate and transform parasitic worms will it be possible to gain a comprehensive understanding of the molecular mechanisms involved in nutrition, metabolism, developmental switches/maturation and interaction with the host immune system. This review focuses on functional genomics approaches in parasitic helminths that are currently used, to highlight potential applications of these technologies in the areas of cell biology, systems biology and immunobiology of parasitic helminths.

Nucleic acid transfection and transgenesis in parasitic nematodes

Transgenesis is an essential tool for assessing gene function in any organism, and it is especially crucial for parasitic nematodes given the dwindling armamentarium of effective anthelmintics and the consequent need to validate essential molecular targets for new drugs and vaccines. Two of the major routes of gene delivery evaluated to date in parasitic nematodes, bombardment with DNA-coated microparticles and intragonadal microinjection of DNA constructs, draw upon experience with the free-living nematode Caenorhabditis elegans. Bombardment has been used to transiently transfect Ascaris suum, Brugia malayi and Litomosoides sigmodontis with both RNA and DNA. Microinjection has been used to achieve heritable transgenesis in Strongyloides stercoralis, S. ratti and Parastrongyloides trichosuri and for additional transient expression studies in B. malayi. A third route of gene delivery revisits a classic method involving DNA transfer facilitated by calcium-mediated permeabilization of recipient cells in developing B. malayi larvae and results in transgene inheritance through host and vector passage. Assembly of microinjected transgenes into multi-copy episomal arrays likely results in their transcriptional silencing in some parasitic nematodes. Methods such as transposon-mediated transgenesis that favour low-copy number chromosomal integration may remedy this impediment to establishing stable transgenic lines. In the future, stable transgenesis in parasitic nematodes could enable loss-of-function approaches by insertional mutagenesis, in situ expression of inhibitory double-stranded RNA or boosting RNAi susceptibility through heterologous expression of dsRNA processing and transport proteins.

Incorporation of a horizontally transferred gene into an operon during cnidarian evolution

Genome sequencing has revealed examples of horizontally transferred genes, but we still know little about how such genes are incorporated into their host genomes. We have previously reported the identification of a gene (flp) that appears to have entered the Hydra genome through horizontal transfer. Here we provide additional evidence in support of our original hypothesis that the transfer was from a unicellular organism, and we show that the transfer occurred in an ancestor of two medusozoan cnidarian species. In addition we show that the gene is part of a bicistronic operon in the Hydra genome. These findings identify a new animal phylum in which trans-spliced leader addition has led to the formation of operons, and define the requirements for evolution of an operon in Hydra. The identification of operons in Hydra also provides a tool that can be exploited in the construction of transgenic Hydra strains.

Annotation of two large contiguous regions from the Haemonchus contortus genome using RNA-seq and comparative analysis with Caenorhabditis elegans

The genomes of numerous parasitic nematodes are currently being sequenced, but their complexity and size, together with high levels of intra-specific sequence variation and a lack of reference genomes, makes their assembly and annotation a challenging task. Haemonchus contortus is an economically significant parasite of livestock that is widely used for basic research as well as for vaccine development and drug discovery. It is one of many medically and economically important parasites within the strongylid nematode group. This group of parasites has the closest phylogenetic relationship with the model organism Caenorhabditis elegans, making comparative analysis a potentially powerful tool for genome annotation and functional studies. To investigate this hypothesis, we sequenced two contiguous fragments from the H. contortus genome and undertook detailed annotation and comparative analysis with C. elegans. The adult H. contortus transcriptome was sequenced using an Illumina platform and RNA-seq was used to annotate a 409 kb overlapping BAC tiling path relating to the X chromosome and a 181 kb BAC insert relating to chromosome I. In total, 40 genes and 12 putative transposable elements were identified. 97.5% of the annotated genes had detectable homologues in C. elegans of which 60% had putative orthologues, significantly higher than previous analyses based on EST analysis. Gene density appears to be less in H. contortus than in C. elegans, with annotated H. contortus genes being an average of two-to-three times larger than their putative C. elegans orthologues due to a greater intron number and size. Synteny appears high but gene order is generally poorly conserved, although areas of conserved microsynteny are apparent. C. elegans operons appear to be partially conserved in H. contortus. Our findings suggest that a combination of RNA-seq and comparative analysis with C. elegans is a powerful approach for the annotation and analysis of strongylid nematode genomes.

Inventory and analysis of ATP-binding cassette (ABC) systems inBrugia malayi

ABC systems are one of the largest described protein superfamilies. These systems have a domain organization that may contain 1 or more transmembrane domains (ABC_TM1F) and 1 or 2 ATP-binding domains (ABC_2). The functions (e.g., import, export and DNA repair) of these proteins distinguish the 3 classes of ABC systems. Mining and PCR-based cloning were used to identify 33 putative ABC systems from theBrugia malayigenome. There were 31 class 2 genes, commonly called ABC transporters, and 2 class 3 genes. The ABC transporters were divided into subfamilies. Three belonged to subfamily A, 16 to subfamily B, 5 to subfamily C, 1 to subfamily E and 3 to subfamilies F and G, respectively. None were placed in subfamilies D and H. Similar to other ABC systems, the ABC_2 domain ofB. malayigenes was conserved and contained the Walker A and B motifs, the signature sequence/linker region and the switch region with the conserved histidine. The ABC_TM1F domain was less conserved. The relative abundance of ABC systems was quantified using real-time reverse transcription PCR and was significantly higher in female adults ofB. malayithan in males and microfilaria, particularly those in subfamilies B and C, which are associated with drug resistance.