A genome sequence survey of the filarial nematode Brugia malayi: repeats, gene discovery, and comparative genomics

The diversity of ABC transporter genes across the Phylum Nematoda

It is estimated that one billion people globally are infected by parasitic nematodes, with children, pregnant women, and the elderly particularly susceptible to morbidity from infection. Control methods are limited to de-worming, which is hampered by rapid re-infection and the inevitable development of anthelmintic resistance. One family of proteins that has been implicated in nematode anthelmintic resistance are the ATP binding cassette (ABC) transporters. ABC transporters are characterized by a highly conserved ATP-binding domain and variable transmembrane regions. A growing number of studies have associated ABC transporters in anthelmintic resistance through a protective mechanism of drug efflux. Genetic deletion of P glycoprotein type ABC transporters in Caenorhabditis elegans demonstrated increased sensitivity to anthelmintics, while in the livestock parasite, Haemonchus contortus, anthelmintic use has been shown to increase the expression of ATP transporter genes. These studies as well as others, provide evidence for a potential role of ABC transporters in drug resistance in nematodes. In order to understand more about the family of ABC transporters, we used hidden Markov models to predict ABC transporter proteins from 108 species across the phylum Nematoda and use these data to analyze patterns of diversification and loss in diverse nematode species. We also examined temporal patterns of expression for the ABC transporter family within the filarial nematode Brugia malayi and identify cases of differential expression across diverse life-cycle stages. Taken together, our data provide a comprehensive overview of ABC transporters in diverse nematode species and identify examples of gene loss and diversification in nematodes based on lifestyle and taxonomy.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

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.

Genomic and proteomic approaches for Chagas’ disease: critical analysis of diagnostic methods

Trypanosoma cruzi is the etiologic agent of Chagas' disease, a chronic inflammatory condition that results in heart and digestive complications. The first draft of the parasite genome is now complete and it is expected that, along with the published genomic and proteomic analyses discussed herein, it will lead to the identification of crucial genes and proteins directly associated with disease. This article reviews the current research trends addressing host-parasite interaction, parasite genetic variability and diagnosis. These advances will certainly bring about major developments not only in our understanding of Trypanosoma cruzi biology, but also in the application of new technologies to disease prevention and control.

Inter and intra-specific diversity of parasites that cause lymphatic filariasis

Lymphatic filariasis is caused by three closely related nematode parasites: Wuchereria bancrofti, Brugia malayi and Brugia timori. These species have many ecological variants that differ in several aspects of their biology such as mosquito vector species, host range, periodicity, and morphology. Although the genome of B. malayi (the first genome sequenced from a parasitic nematode) has been available for more than five years, very little is known about genetic variability among the lymphatic dwelling filariae. The genetic diversity among these worms is not only interesting from a biological perspective, but it may have important practical implications for the Global Program to Eliminate Lymphatic Filariasis, as the parasites may respond differently to diagnostic tests and/or medical interventions. Therefore, better information on their genetic variability is urgently needed. With improved methods for nucleic acid extraction and recent advances in sequencing chemistry and instrumentation, this gap can be filled relatively inexpensively. Improved information on filarial genetic diversity may increase the chances of success for lymphatic filariasis elimination programs.

Diversity in parasitic nematode genomes: the microRNAs of Brugia pahangi and Haemonchus contortus are largely novel

Background

MicroRNAs (miRNAs) play key roles in regulating post-transcriptional gene expression and are essential for development in the free-living nematode Caenorhabditis elegans and in higher organisms. Whether microRNAs are involved in regulating developmental programs of parasitic nematodes is currently unknown. Here we describe the the miRNA repertoire of two important parasitic nematodes as an essential first step in addressing this question.

Results

The small RNAs from larval and adult stages of two parasitic species, Brugia pahangi and Haemonchus contortus, were identified using deep-sequencing and bioinformatic approaches. Comparative analysis to known miRNA sequences reveals that the majority of these miRNAs are novel. Some novel miRNAs are abundantly expressed and display developmental regulation, suggesting important functional roles. Despite the lack of conservation in the miRNA repertoire, genomic positioning of certain miRNAs within or close to specific coding genes is remarkably conserved across diverse species, indicating selection for these associations. Endogenous small-interfering RNAs and Piwi-interacting (pi)RNAs, which regulate gene and transposon expression, were also identified. piRNAs are expressed in adult stage H. contortus, supporting a conserved role in germline maintenance in some parasitic nematodes.

Conclusions

This in-depth comparative analysis of nematode miRNAs reveals the high level of divergence across species and identifies novel sequences potentially involved in development. Expression of novel miRNAs may reflect adaptations to different environments and lifestyles. Our findings provide a detailed foundation for further study of the evolution and function of miRNAs within nematodes and for identifying potential targets for intervention.

Cloning and bioinformatic identification of small RNAs in the filarial nematode, Brugia malayi

Characterization of small RNAs from the filarial nematode Brugia malayi is the initial step in understanding their role in gene silencing. Both RNA cloning and bioinformatics were used to identify 32 microRNAs (miRNAs) belonging to 24 families. One family, miR-36 only occurs in helminths including B. malayi. Several of the miRNAs are arranged in clusters and are coordinately expressed as determined by northern blot analysis. In addition, small RNAs were identified from Pao/Bleo retrotransposons and their associated repeat sequences indicating that B. malayi uses an RNAi mechanism to maintain genome integrity. Analysis of these data provides a first glimpse into how small RNA-mediated silencing pathways regulate the parasitic life cycle of B. malayi.

Genomics-based drug design targets the AT-rich malaria parasite: implications for antiparasite chemotherapy

Evaluation of: Woynarowski JM, Krugliak M, Ginsburg H: Pharmacogenomic analyses of targeting the AT-rich malaria parasite genome with AT-specific alkylating drugs. Mol. Biochem. Parasitol. 154(1), 70-81 (2007) [1] . The sequencing of the malaria genome sought to expose the parasite's ability to cause disease and identify new targets for antimalarial drugs and vaccines. In this study, the authors discovered how malaria genomic DNA, which is unusually rich in adenine and thymine nucleotides, is intrinsically a target for a selective class of compounds. AT-specific DNA-binding agents have previously been shown to have potent antimalarial activity in vitro. The authors used high-resolution bioinformatic tools to explore the genomic basis for this drug susceptibility, first at the level of individual DNA-binding sites, then expanding to the entire genomic context of each malaria chromosome. Their findings revealed a nonrandom distribution and organization of drug-binding sites that can be further exploited to target these AT sequences. Based on these findings, comparative bioinformatics analyses with other parasite genomes may lead to the identification of new target organisms for these AT-specific drugs and have wide implications for the treatment of human and animal parasitic diseases.

Genes encoding putative biogenic amine receptors in the parasitic nematode Brugia malayi

Filarial nematodes, such as Brugia malayi, cause major health problems worldwide. The lack of a vaccine against B. malayi, combined with ineffective chemotherapy against the adult has prompted the examination of biogenic amine receptors (BARs) as possible targets for drug discovery. We employed bioinformatics to identify genes encoding putative B. malayi BARs. Surprisingly, the B. malayi genome contains half of the genes predicted to encode BARs in the genomes of free-living nematodes such as Caenorhabditis elegans or C. briggsae; however, all of the predicted B. malayi receptors have clear orthologues in C. elegans. The B. malayi genes encode each of the major BAR subclasses, including three serotonin, two dopamine and two tyramine/octopamine receptors and the structure of orthologous BAR genes is conserved. We find that potential G-protein coupling and ligand-specificity of individual BARs may be predicted by phylogenetic comparisons. Our results provide a framework for how G-protein coupled receptors may be targeted for drug development in medically important parasitic nematodes.

Response of Armigeres subalbatus (Diptera: Culicidae) to intraperitoneally isolated Brugia spp. microfilariae

The relationship between mosquito and parasite involves a delicate balance that is influenced not only by the mosquito but also by parasite determinants. Using the biologically and morphologically similar parasites Brugia malayi and Brugia pahangi and the mosquito Armigeres subalbatus (Coquillett) (Diptera: Culicidae), it should be possible to dissect out the key elements involved in initiating or avoiding an immune response, known as melanotic encapsulation, because in this mosquito B. malayi microfilariae (mf) are melanized and destroyed, but B. pahangi mf develop normally into infective-stage larvae. Because of limitations in isolating sufficient mf from the circulation of an infected mammalian host, Brugia spp. mf that can be obtained in large numbers from the peritoneal cavity of an infected host were tested to ascertain the immune response of Ar. subalbatus to this source of mf. Results indicate that the immune response of Ar. subalbatus against intraperitoneal (i.p.) Brugia spp. mf mimics that which is observed when this mosquito is exposed to mf-infected animals, indicating that i.p. mf are similar to those mf that circulate naturally in the blood of the vertebrate host. Therefore, the i.p. mf should serve as an excellent source of material for genomic and proteomic studies designed to analyze the role of the parasite in influencing the immune response of the mosquito.

Penetration of the mosquito midgut is not required for Brugia pahangi microfilariae to avoid the melanotic encapsulation response of Armigeres subalbatus

Insect vectors of disease have the capacity to respond to, and prevent further development of, parasites and pathogens using a response known as melanotic encapsulation. The naturally-occurring Armigeres subalbatus-Brugia spp. system provides an excellent way to investigate melanotic encapsulation and immune recognition in a mosquito host, because Brugia malayi microfilariae (mf) acquired via a blood meal are rapidly melanized in the body cavity of Ar. subalbatus, but Brugia pahangi mf evade or suppress the immune response and develop normally into infective stage larvae. Previous studies have suggested that B. pahangi mf are changed in some manner in the process of exiting the mosquito gut, thereby facilitating escape from, or suppression of, the melanotic encapsulation response. By inoculating mosquitoes with parasites, thus circumventing the midgut, we show that approximately 88% of B. pahangi mf escape the melanotic encapsulation response while approximately 90% of inoculated B. malayi mf are melanized. Methods to isolate parasites for this procedure are described. These results mimic those observed in Ar. subalbatus against Brugia spp. mf that are ingested following blood feeding, and demonstrate that midgut penetration is not required for B. pahangi mf to avoid the melanotic encapsulation response of Ar. subalbatus.

Operon conservation and the evolution of trans-splicing in the phylum Nematoda

The nematode Caenorhabditis elegans is unique among model animals in that many of its genes are cotranscribed as polycistronic pre-mRNAs from operons. The mechanism by which these operonic transcripts are resolved into mature mRNAs includes trans-splicing to a family of SL2-like spliced leader exons. SL2-like spliced leaders are distinct from SL1, the major spliced leader in C. elegans and other nematode species. We surveyed five additional nematode species, representing three of the five major clades of the phylum Nematoda, for the presence of operons and the use of trans-spliced leaders in resolution of polycistronic pre-mRNAs. Conserved operons were found in Pristionchus pacificus, Nippostrongylus brasiliensis, Strongyloides ratti, Brugia malayi, and Ascaris suum. In nematodes closely related to the rhabditine C. elegans, a related family of SL2-like spliced leaders is used for operonic transcript resolution. However, in the tylenchine S. ratti operonic transcripts are resolved using a family of spliced leaders related to SL1. Non-operonic genes in S. ratti may also receive these SL1 variants. In the spirurine nematodes B. malayi and A. suum operonic transcripts are resolved using SL1. Mapping these phenotypes onto the robust molecular phylogeny for the Nematoda suggests that operons evolved before SL2-like spliced leaders, which are an evolutionary invention of the rhabditine lineage.

The genome of the nematode worm Caenorhabditis elegans was the first of any animal to be completely sequenced. One surprising finding in this worm's genome was that about one-fifth of its genes were organised as sets of from two to eight genes expressed from the same promoter, similar to bacterial “operons.” The pre-mRNAs made from these operons are processed by an intermolecular ligation process called SL trans-splicing. Other animal genomes, such as the human genome or that of the fruit fly contain neither operons nor SL trans-splicing. In this article, Guiliano and Blaxter have investigated whether this curious facet of genome organisation is peculiar to C. elegans and close relatives by examining the genomes of a wide range of parasitic and free-living nematodes. The authors find that both operons and trans-splicing are present across the nematodes, and that operons evolve as other genome features do. All of the species surveyed use trans-splicing to resolve their multigene pre-mRNAs into single-gene mRNAs, but the details differ significantly from the process in C. elegans. In particular, the short piece of RNA that is attached to the beginning of operon-derived mRNAs has changed independently in many nematode groups.

Testing the efficacy of RNA interference in Haemonchus contortus

RNA interference (RNAi) is widely used in Caenorhabditis elegans to identify gene function and has been adapted as a high throughput screening method to identify genes involved in essential processes. We have been examining whether RNAi could also be used on the strongylid parasitic nematode Haemonchus contortus to study gene function. Eleven genes were targeted in L1 and exsheathed L3 H. contortus larvae with RNAi methodologies which have been shown to be effective in C. elegans and parasitic nematodes-feeding, soaking and electroporation. Reverse transcriptase-PCR and, where possible, protein assays were carried out to examine decreases in mRNA and protein levels. RNAi soaking in dsRNA to beta-tubulin and sec-23, a gene involved in vesicle transport, resulted in specific decreases in mRNA levels in exsheathed L3 larvae. No signs of specific decreases in expression levels were observed for the other nine genes tested. Following electroporation of dsRNA in L1 stage larvae, significant decreases were observed for two out of four genes tested. These findings suggest that the RNAi pathway is functional in H. contortus and that, under certain conditions, it is possible to suppress gene expression by RNAi. However, it only works on a limited number of genes and in some cases the effect is small and difficult to reproduce. This indicates that the RNAi approaches established for C. elegans and other nematodes have limited efficacy in H. contortus. This may reflect differences between nematode species in dsRNA uptake and transport into cells and between cells.

Identification and analysis of genes expressed in the adult filarial parasitic nematode Dirofilaria immitis

The heartworm Dirofilaria immitis is a filarial parasitic nematode infecting dogs and other mammals worldwide causing fatal complications. Here, we present the first large-scale survey of the adult heartworm transcriptome by generation and analysis of 4005 expressed sequence tags, identifying about 1800 genes and expanding the available sequence information for the parasite significantly. Brugia malayi genomic data offered the most valuable information to interpret heartworm genes, with about 70% of D. immitis genes showing significant similarities to the assembly. Comparative genomic analyses revealed both genes common to metazoans or nematodes and genes specific to filarial parasites that may relate to parasitism. Characterization of abundant transcripts suggested important roles for genes involved in energy generation and antioxidant defense in adults. In particular, we proposed that adult heartworm likely adopted an anaerobic electron transfer-based energy generation system distinct from the aerobic pathway utilized by its mammalian host, making it a promising target in developing next generation macrofilaricides and other treatments. Our survey provided novel insights into the D. immitis transcriptome and laid a foundation for further comparative studies on biology, parasitism and evolution within the phylum Nematoda.

The biology and genomics of Strongyloides

Parasitic nematodes are widespread and important pathogens of humans and other animals. The parasitic nematodes Strongyloides have an unusual life cycle in which there is a facultative free-living generation in addition to the obligate parasitic generation. The genomes of many species of parasitic nematodes, including Strongyloides ratti and Strongyloides stercoralis, have been investigated, principally by expressed sequence tag (EST) analyses. These have discovered very many genes from these parasites but, in so doing, have also revealed how different these species are from each other and from other organisms. Understanding the role and function of these newly discovered genes is now the challenge, made more difficult by the parasitic lifestyle. The genomic information available for parasitic nematodes is allowing new approaches for the control of parasitic nematodes to be considered.

Comparative genomics of nematodes

Recent transcriptome and genome projects have dramatically expanded the biological data available across the phylum Nematoda. Here we summarize analyses of these sequences, which have revealed multiple unexpected results. Despite a uniform body plan, nematodes are more diverse at the molecular level than was previously recognized, with many species- and group-specific novel genes. In the genus Caenorhabditis, changes in chromosome arrangement, particularly local inversions, are also rapid, with breakpoints occurring at 50-fold the rate in vertebrates. Tylenchid plant parasitic nematode genomes contain several genes closely related to genes in bacteria, implicating horizontal gene transfer events in the origins of plant parasitism. Functional genomics techniques are also moving from Caenorhabditis elegans to application throughout the phylum. Soon, eight more draft nematode genome sequences will be available. This unique resource will underpin both molecular understanding of these most abundant metazoan organisms and aid in the examination of the dynamics of genome evolution in animals.

An expressed sequence tag analysis of the life-cycle of the parasitic nematode Strongyloides ratti

14,761 expressed sequence tags (ESTs) were generated, representing five stages during the parasitic and free-living phases of the life-cycle of the parasitic nematode Strongyloides ratti. These ESTs formed 4152 clusters, of which 97% contained 10 or fewer ESTs and 66% were singletons. These 4152 clusters are likely to represent approximately 20% of S. ratti's genes. The clusters' consensus sequences were used to assign each cluster to one of three databases: (i) Caenorhabditis elegans and C. briggsae sequences; (ii) other nematode sequences; (iii) non-nematode sequences. This approach has identified putative nematode-specific genes, that may be targets for developing approaches for parasitic nematode control. Approximately 25% of the clusters have no significant alignments and may therefore represent novel genes. The EST representation between the libraries was used to analyse stage-specific or -biased expression in silico. This showed that 81% of clusters are present in only one library and 12% are present in any two libraries, indicating substantial stage-specificity of gene expression. The 30-most abundantly expressed clusters were analysed in further detail. Many of these have significantly different parasitic- or free-living-specific or -biased expression. Many of the parasitic-specific genes are, as yet, uncharacterised: one of these represents 25% of all ESTs obtained from the parasitic stage.

Heterologous expression of the filarial nematode alt gene products reveals their potential to inhibit immune function

Background

Parasites exploit sophisticated strategies to evade host immunity that require both adaptation of existing genes and evolution of new gene families. We have addressed this question by testing the immunological function of novel genes from helminth parasites, in which conventional transgenesis is not yet possible. We investigated two such novel genes from Brugia malayi termed abundant larval transcript (alt), expression of which reaches ~5% of total transcript at the time parasites enter the human host.

Results

To test the hypothesis that ALT proteins modulate host immunity, we adopted an alternative transfection strategy to express these products in the protozoan parasite Leishmania mexicana. We then followed the course of infection in vitro in macrophages and in vivo in mice. Expression of ALT proteins, but not a truncated mutant, conferred greater infectivity of macrophages in vitro, reaching 3-fold higher parasite densities. alt-transfected parasites also caused accelerated disease in vivo, and fewer mice were able to clear infection of organisms expressing ALT. alt-transfected parasites were more resistant to IFN-γ-induced killing by macrophages. Expression profiling of macrophages infected with transgenic L. mexicana revealed consistently higher levels of GATA-3 and SOCS-1 transcripts, both associated with the Th2-type response observed in in vivo filarial infection.

Conclusion

Leishmania transfection is a tractable and informative approach to determining immunological functions of single genes from heterologous organisms. In the case of the filarial ALT proteins, our data suggest that they may participate in the Th2 bias observed in the response to parasite infection by modulating cytokine-induced signalling within immune system cells.

A transcriptomic analysis of the phylum Nematoda

The phylum Nematoda occupies a huge range of ecological niches, from free-living microbivores to human parasites. We analyzed the genomic biology of the phylum using 265,494 expressed-sequence tag sequences, corresponding to 93,645 putative genes, from 30 species, including 28 parasites. From 35% to 70% of each species' genes had significant similarity to proteins from the model nematode Caenorhabditis elegans. More than half of the putative genes were unique to the phylum, and 23% were unique to the species from which they were derived. We have not yet come close to exhausting the genomic diversity of the phylum. We identified more than 2,600 different known protein domains, some of which had differential abundances between major taxonomic groups of nematodes. We also defined 4,228 nematode-specific protein families from nematode-restricted genes: this class of genes probably underpins species- and higher-level taxonomic disparity. Nematode-specific families are particularly interesting as drug and vaccine targets.