Indoleacrylic Acid Produced by Commensal Peptostreptococcus Species Suppresses Inflammation

Host factors in the intestine help select for bacteria that promote health. Certain commensals can utilize mucins as an energy source, thus promoting their colonization. However, health conditions such as inflammatory bowel disease (IBD) are associated with a reduced mucus layer, potentially leading to dysbiosis associated with this disease. We characterize the capability of commensal species to cleave and transport mucin-associated monosaccharides and identify several Clostridiales members that utilize intestinal mucins. One such mucin utilizer, Peptostreptococcus russellii, reduces susceptibility to epithelial injury in mice. Several Peptostreptococcus species contain a gene cluster enabling production of the tryptophan metabolite indoleacrylic acid (IA), which promotes intestinal epithelial barrier function and mitigates inflammatory responses. Furthermore, metagenomic analysis of human stool samples reveals that the genetic capability of microbes to utilize mucins and metabolize tryptophan is diminished in IBD patients. Our data suggest that stimulating IA production could promote anti-inflammatory responses and have therapeutic benefits.

Mining novel effector proteins from the esophageal gland cells of Meloidogyne incognita

Meloidogyne incognita is one of the most economically damaging plant pathogens in agriculture and horticulture. Identifying and characterizing the effector proteins which M. incognita secretes into its host plants during infection is an important step toward finding new ways to manage this pest. In this study, we have identified the cDNAs for 18 putative effectors (i.e., proteins that have the potential to facilitate M. incognita parasitism of host plants). These putative effectors are secretory proteins that do not contain transmembrane domains and whose genes are specifically expressed in the secretory gland cells of the nematode, indicating that they are likely secreted from the nematode through its stylet. We have determined that, in the plant cells, these putative effectors are likely to localize to the cytoplasm. Furthermore, the transcripts of many of these novel effectors are specifically upregulated during different stages of the nematode's life cycle, indicating that they function at specific stages during M. incognita parasitism. The predicted proteins showed little to no homology to known proteins from free-living nematode species, suggesting that they evolved recently to support the parasitic lifestyle. On the other hand, several of the effectors are part of gene families within the M. incognita genome as well as that of M. hapla, which points to an important role that these putative effectors are playing in both parasites. With the discovery of these putative effectors, we have increased our knowledge of the effector repertoire utilized by root-knot nematodes to infect, feed on, and reproduce on their host plants. Future studies investigating the roles that these proteins play in planta will help mitigate the effects of this damaging pest.

The prevalence of species and strains in the human microbiome: a resource for experimental efforts

Experimental efforts to characterize the human microbiota often use bacterial strains that were chosen for historical rather than biological reasons. Here, we report an analysis of 380 whole-genome shotgun samples from 100 subjects from the NIH Human Microbiome Project. By mapping their reads to 1,751 reference genome sequences and analyzing the resulting relative strain abundance in each sample we present metrics and visualizations that can help identify strains of interest for experimentalists. We also show that approximately 14 strains of 10 species account for 80% of the mapped reads from a typical stool sample, indicating that the function of a community may not be irreducibly complex. Some of these strains account for >20% of the sequence reads in a subset of samples but are absent in others, a dichotomy that could underlie biological differences among subjects. These data should serve as an important strain selection resource for the community of researchers who take experimental approaches to studying the human microbiota.

Identification and characterization of alternative splicing in parasitic nematode transcriptomes

Background

Alternative splicing (AS) of mRNA is a vital mechanism for enhancing genomic complexity in eukaryotes. Spliced isoforms of the same gene can have diverse molecular and biological functions and are often differentially expressed across various tissues, times, and conditions. Thus, AS has important implications in the study of parasitic nematodes with complex life cycles. Transcriptomic datasets are available from many species, but data must be revisited with splice-aware assembly protocols to facilitate the study of AS in helminthes.

Methods

We sequenced cDNA from the model worm Caenorhabditis elegans using 454/Roche technology for use as an experimental dataset. Reads were assembled with Newbler software, invoking the cDNA option. Several combinations of parameters were tested and assembled transcripts were verified by comparison with previously reported C. elegans genes and transcript isoforms and with Illumina RNAseq data.

Results

Thoughtful adjustment of program parameters increased the percentage of assembled transcripts that matched known C. elegans sequences, decreased mis-assembly rates (i.e., cis- and trans-chimeras), and improved the coverage of the geneset. The optimized protocol was used to update de novo transcriptome assemblies from nine parasitic nematode species, including important pathogens of humans and domestic animals. Our assemblies indicated AS rates in the range of 20-30%, typically with 2-3 transcripts per AS locus, depending on the species. Transcript isoforms from the nine species were translated and searched for similarity to known proteins and functional domains. Some 21 InterPro domains, including several involved in nucleotide and chromatin binding, were statistically correlated with AS genetic loci. In most cases, the Roche/454 data explored in this study are the only sequences available from the species in question; however, the recently published genome of the human hookworm Necator americanus provided an additional opportunity to validate our results.

Conclusions

Our optimized assembly parameters facilitated the first survey of AS among parasitic nematodes. The nine transcriptome assemblies, their protein translations, and basic annotations are available from Nematode.net as a resource for the research community. These should be useful for studies of specific genes and gene families of interest as well as for curating draft genome assemblies as they become available.

Sepsis from the gut: the enteric habitat of bacteria that cause late-onset neonatal bloodstream infections

BACKGROUND

Late-onset sepsis is a major problem in neonatology, but the habitat of the pathogens before bloodstream invasion occurs is not well established.

METHODS

We examined prospectively collected stools from premature infants with sepsis to find pathogens that subsequently invaded their bloodstreams, and sought the same organisms in stools of infants without sepsis. Culture-based techniques were used to isolate stool bacteria that provisionally matched the bloodstream organisms, which were then genome sequenced to confirm or refute commonality.

RESULTS

Of 11 children with late-onset neonatal bloodstream infections, 7 produced at least 1 stool that contained group B Streptococcus (GBS), Serratia marcescens, or Escherichia coli before their sepsis episode with provisionally matching organisms. Of 96 overlap comparison subjects without sepsis temporally associated with these cases, 4 were colonized with provisionally matching GBS or S. marcescens. Of 175 comparisons of stools from randomly selected infants without sepsis, 1 contained a GBS (this infant had also served as an overlap comparison subject and both specimens contained provisionally matching GBS). Genome sequencing confirmed common origin of provisionally matching fecal and blood isolates. The invasive E. coli were present in all presepticemic stools since birth, but gut colonization with GBS and S. marcescens occurred closer to time of bloodstream infection.

CONCLUSIONS

The neonatal gut harbors sepsis-causing pathogens, but such organisms are not inevitable members of the normal microbiota. Surveillance microbiology, decolonization, and augmented hygiene might prevent dissemination of invasive bacteria between and within premature infants.

Discovery of anthelmintic drug targets and drugs using chokepoints in nematode metabolic pathways

Parasitic roundworm infections plague more than 2 billion people (1/3 of humanity) and cause drastic losses in crops and livestock. New anthelmintic drugs are urgently needed as new drug resistance and environmental concerns arise. A “chokepoint reaction” is defined as a reaction that either consumes a unique substrate or produces a unique product. A chokepoint analysis provides a systematic method of identifying novel potential drug targets. Chokepoint enzymes were identified in the genomes of 10 nematode species, and the intersection and union of all chokepoint enzymes were found. By studying and experimentally testing available compounds known to target proteins orthologous to nematode chokepoint proteins in public databases, this study uncovers features of chokepoints that make them successful drug targets. Chemogenomic screening was performed on drug-like compounds from public drug databases to find existing compounds that target homologs of nematode chokepoints. The compounds were prioritized based on chemical properties frequently found in successful drugs and were experimentally tested using Caenorhabditis elegans. Several drugs that are already known anthelmintic drugs and novel candidate targets were identified. Seven of the compounds were tested in Caenorhabditis elegans and three yielded a detrimental phenotype. One of these three drug-like compounds, Perhexiline, also yielded a deleterious effect in Haemonchus contortus and Onchocerca lienalis, two nematodes with divergent forms of parasitism. Perhexiline, known to affect the fatty acid oxidation pathway in mammals, caused a reduction in oxygen consumption rates in C. elegans and genome-wide gene expression profiles provided an additional confirmation of its mode of action. Computational modeling of Perhexiline and its target provided structural insights regarding its binding mode and specificity. Our lists of prioritized drug targets and drug-like compounds have potential to expedite the discovery of new anthelmintic drugs with broad-spectrum efficacy.

The World Health Organization estimates that 2.9 million people are infected with parasitic roundworms, causing high-morbidity and mortality rates, developmental delays in children, and low productivity of affected individuals. The agricultural industry experiences drastic losses in crop and livestock due to parasitic worm infections. Therefore, there is an urgent need to identify new targets and drugs to fight parasitic nematode infection. This study identified metabolic chokepoint compounds that were either produced or consumed by a single reaction and elucidated the chokepoint enzyme that drives the reaction. If the enzyme that catalyzes that reaction is blocked, a toxic build-up of a compound or lack of compound necessary for subsequent reaction will occur, potentially causing adverse effects to the parasite organism. Compounds that target some of the chokepoint enzymes were tested in C. elegans and several compounds showed efficacy. One drug-like compound, Perhexiline, showed efficacy in two different parasitic worms and yielded expected physiological effects, indicating that this drug-like compound may have efficacy on a pan-phylum level through the predicted mode of action. The methodology to find and prioritize metabolic chokepoint targets and prioritize compounds could be applied to other parasites.

mBLAST: Keeping up with the sequencing explosion for (meta)genome analysis

Recent advances in next-generation sequencing technologies require alignment algorithms and software that can keep pace with the heightened data production. Standard algorithms, especially protein similarity searches, represent significant bottlenecks in analysis pipelines. For metagenomic approaches in particular, it is now often necessary to search hundreds of millions of sequence reads against large databases. Here we describe mBLAST, an accelerated search algorithm for translated and/or protein alignments to large datasets based on the Basic Local Alignment Search Tool (BLAST) and retaining the high sensitivity of BLAST. The mBLAST algorithms achieve substantial speed up over the National Center for Biotechnology Information (NCBI) programs BLASTX, TBLASTX and BLASTP for large datasets, allowing analysis within reasonable timeframes on standard computer architectures. In this article, the impact of mBLAST is demonstrated with sequences originating from the microbiota of healthy humans from the Human Microbiome Project. mBLAST is designed as a plug-in replacement for BLAST for any study that involves short-read sequences and includes high-throughput analysis. The mBLAST software is freely available to academic users at www.multicorewareinc.com.

Widespread colonization of the lung by Tropheryma whipplei in HIV infection

RATIONALE

Lung infections caused by opportunistic or virulent pathogens are a principal cause of morbidity and mortality in HIV infection. It is unknown whether HIV infection leads to changes in basal lung microflora, which may contribute to chronic pulmonary complications that increasingly are being recognized in individuals infected with HIV.

OBJECTIVES

To determine whether the immunodeficiency associated with HIV infection resulted in alteration of the lung microbiota.

METHODS

We used 16S ribosomal RNA targeted pyrosequencing and shotgun metagenomic sequencing to analyze bacterial gene sequences in bronchoalveolar lavage (BAL) and mouths of 82 HIV-positive and 77 HIV-negative subjects.

MEASUREMENTS AND MAIN RESULTS

Sequences representing Tropheryma whipplei, the etiologic agent of Whipple's disease, were significantly more frequent in BAL of HIV-positive compared with HIV-negative individuals. T. whipplei dominated the community (>50% of sequence reads) in 11 HIV-positive subjects, but only 1 HIV-negative individual (13.4 versus 1.3%; P = 0.0018). In 30 HIV-positive individuals sampled longitudinally, antiretroviral therapy resulted in a significantly reduced relative abundance of T. whipplei in the lung. Shotgun metagenomic sequencing was performed on eight BAL samples dominated by T. whipplei 16S ribosomal RNA. Whole genome assembly of pooled reads showed that uncultured lung-derived T. whipplei had similar gene content to two isolates obtained from subjects with Whipple's disease.

CONCLUSIONS

Asymptomatic subjects with HIV infection have unexpected colonization of the lung by T. whipplei, which is reduced by effective antiretroviral therapy and merits further study for a potential pathogenic role in chronic pulmonary complications of HIV infection.

A nuclear single-nucleotide polymorphism (SNP) potentially useful for the separation of Rhodnius prolixus from members of the Rhodnius robustus cryptic species complex (Hemiptera: Reduviidae)

The design and application of rational strategies that rely on accurate species identification are pivotal for effective vector control. When morphological identification of the target vector species is impractical, the use of molecular markers is required. Here we describe a non-coding, single-copy nuclear DNA fragment that contains a single-nucleotide polymorphism (SNP) with the potential to distinguish the important domestic Chagas disease vector, Rhodnius prolixus, from members of the four sylvatic Rhodnius robustus cryptic species complex. A total of 96 primer pairs obtained from whole genome shotgun sequencing of the R. prolixus genome (12,626 random reads) were tested on 43 R. prolixus and R. robustus s.l. samples. One of the seven amplicons selected (AmpG) presented a SNP, potentially diagnostic for R. prolixus, on the 280th site. The diagnostic nature of this SNP was then confirmed based on the analysis of 154 R. prolixus and R. robustus s.l. samples representing the widest possible geographic coverage. The results of a 60% majority-rule Bayesian consensus tree and a median-joining network constructed based on the genetic variability observed reveal the paraphyletic nature of the R. robustus species complex, with respect to R. prolixus. The AmpG region is located in the fourth intron of the Transmembrane protein 165 gene, which seems to be in the R. prolixus X chromosome. Other possible chromosomal locations of the AmpG region in the R. prolixus genome are also presented and discussed.

Transcriptomic and proteomic analyses of a Wolbachia-free filarial parasite provide evidence of trans-kingdom horizontal gene transfer

Most filarial parasites in the subfamilies Onchocercinae and Dirofilariinae depend on Wolbachia endobacteria to successfully carry out their life cycle. Recently published data indicate that the few Wolbachia-free species in these subfamilies were infected in the distant past and have subsequently shed their endosymbionts. We used an integrated transcriptomic and proteomic analysis of Onchocerca flexuosa to explore the molecular mechanisms that allow worms of this species to survive without a bacterial partner. Roche/454 sequencing of the adult transcriptome produced 16,814 isogroup and 47,252 singleton sequences that are estimated to represent approximately 41% of the complete gene set. Sequences similar to 97 Wolbachia genes were identified from the transcriptome, some of which appear on the same transcripts as sequences similar to nematode genes. Computationally predicted peptides, including those with similarity to Wolbachia proteins, were classified at the domain and pathway levels in order to assess the metabolic capabilities of O. flexuosa and compare against the Wolbachia-dependent model filaria, Brugia malayi. Transcript data further facilitated a shotgun proteomic analysis of O. flexuosa adult worm lysate, resulting in the identification of 1,803 proteins. Three of the peptides detected by mass spectroscopy map to two ABC transport-related proteins from Wolbachia. Antibodies raised to one of the Wolbachia-like peptides labeled a single 38 kDa band on Western blots of O. flexuosa lysate and stained specific worm tissues by immunohistology. Future studies will be required to determine the exact functions of Wolbachia-like peptides and proteins in O. flexuosa and to assess their roles in worm biology.

Exploring metazoan evolution through dynamic and holistic changes in protein families and domains

Background

Proteins convey the majority of biochemical and cellular activities in organisms. Over the course of evolution, proteins undergo normal sequence mutations as well as large scale mutations involving domain duplication and/or domain shuffling. These events result in the generation of new proteins and protein families. Processes that affect proteome evolution drive species diversity and adaptation. Herein, change over the course of metazoan evolution, as defined by birth/death and duplication/deletion events within protein families and domains, was examined using the proteomes of 9 metazoan and two outgroup species.

Results

In studying members of the three major metazoan groups, the vertebrates, arthropods, and nematodes, we found that the number of protein families increased at the majority of lineages over the course of metazoan evolution where the magnitude of these increases was greatest at the lineages leading to mammals. In contrast, the number of protein domains decreased at most lineages and at all terminal lineages. This resulted in a weak correlation between protein family birth and domain birth; however, the correlation between domain birth and domain member duplication was quite strong. These data suggest that domain birth and protein family birth occur via different mechanisms, and that domain shuffling plays a role in the formation of protein families. The ratio of protein family birth to protein domain birth (domain shuffling index) suggests that shuffling had a more demonstrable effect on protein families in nematodes and arthropods than in vertebrates. Through the contrast of high and low domain shuffling indices at the lineages of Trichinella spiralis and Gallus gallus, we propose a link between protein redundancy and evolutionary changes controlled by domain shuffling; however, the speed of adaptation among the different lineages was relatively invariant. Evaluating the functions of protein families that appeared or disappeared at the last common ancestors (LCAs) of the three metazoan clades supports a correlation with organism adaptation. Furthermore, bursts of new protein families and domains in the LCAs of metazoans and vertebrates are consistent with whole genome duplications.

Conclusion

Metazoan speciation and adaptation were explored by birth/death and duplication/deletion events among protein families and domains. Our results provide insights into protein evolution and its bearing on metazoan evolution.

A case study for large-scale human microbiome analysis using JCVI's metagenomics reports (METAREP)

As metagenomic studies continue to increase in their number, sequence volume and complexity, the scalability of biological analysis frameworks has become a rate-limiting factor to meaningful data interpretation. To address this issue, we have developed JCVI Metagenomics Reports (METAREP) as an open source tool to query, browse, and compare extremely large volumes of metagenomic annotations. Here we present improvements to this software including the implementation of a dynamic weighting of taxonomic and functional annotation, support for distributed searches, advanced clustering routines, and integration of additional annotation input formats. The utility of these improvements to data interpretation are demonstrated through the application of multiple comparative analysis strategies to shotgun metagenomic data produced by the National Institutes of Health Roadmap for Biomedical Research Human Microbiome Project (HMP) (http://nihroadmap.nih.gov). Specifically, the scalability of the dynamic weighting feature is evaluated and established by its application to the analysis of over 400 million weighted gene annotations derived from 14 billion short reads as predicted by the HMP Unified Metabolic Analysis Network (HUMAnN) pipeline. Further, the capacity of METAREP to facilitate the identification and simultaneous comparison of taxonomic and functional annotations including biological pathway and individual enzyme abundances from hundreds of community samples is demonstrated by providing scenarios that describe how these data can be mined to answer biological questions related to the human microbiome. These strategies provide users with a reference of how to conduct similar large-scale metagenomic analyses using METAREP with their own sequence data, while in this study they reveal insights into the nature and extent of variation in taxonomic and functional profiles across body habitats and individuals. Over one thousand HMP WGS datasets and the latest open source code are available at http://www.jcvi.org/hmp-metarep.

Metabolic reconstruction for metagenomic data and its application to the human microbiome

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/humann. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies.

Effects of doxycycline on gene expression in Wolbachia and Brugia malayi adult female worms in vivo

Background

Most filarial nematodes contain Wolbachia symbionts. The purpose of this study was to examine the effects of doxycycline on gene expression in Wolbachia and adult female Brugia malayi.

Methods

Brugia malayi infected gerbils were treated with doxycycline for 6-weeks. This treatment largely cleared Wolbachia and arrested worm reproduction. RNA recovered from treated and control female worms was labeled by random priming and hybridized to the Version 2- filarial microarray to obtain expression profiles.

Results and discussion

Results showed significant changes in expression for 200 Wolbachia (29% of Wolbachia genes with expression signals in untreated worms) and 546 B. malayi array elements after treatment. These elements correspond to known genes and also to novel genes with unknown biological functions. Most differentially expressed Wolbachia genes were down-regulated after treatment (98.5%). In contrast, doxycycline had a mixed effect on B. malayi gene expression with many more genes being significantly up-regulated after treatment (85% of differentially expressed genes). Genes and processes involved in reproduction (gender-regulated genes, collagen, amino acid metabolism, ribosomal processes, and cytoskeleton) were down-regulated after doxycycline while up-regulated genes and pathways suggest adaptations for survival in response to stress (energy metabolism, electron transport, anti-oxidants, nutrient transport, bacterial signaling pathways, and immune evasion).

Conclusions

Doxycycline reduced Wolbachia and significantly decreased bacterial gene expression. Wolbachia ribosomes are believed to be the primary biological target for doxycycline in filarial worms. B. malayi genes essential for reproduction, growth and development were also down-regulated; these changes are consistent with doxycycline effects on embryo development and reproduction. On the other hand, many B. malayi genes involved in energy production, electron-transport, metabolism, anti-oxidants, and others with unknown functions had increased expression signals after doxycycline treatment. These results suggest that female worms are able to compensate in part for the loss of Wolbachia so that they can survive, albeit without reproductive capacity. This study of doxycycline induced changes in gene expression has provided new clues regarding the symbiotic relationship between Wolbachia and B. malayi.

Nematode.net update 2011: addition of data sets and tools featuring next-generation sequencing data

Nematode.net (http://nematode.net) has been a publicly available resource for studying nematodes for over a decade. In the past 3 years, we reorganized Nematode.net to provide more user-friendly navigation through the site, a necessity due to the explosion of data from next-generation sequencing platforms. Organism-centric portals containing dynamically generated data are available for over 56 different nematode species. Next-generation data has been added to the various data-mining portals hosted, including NemaBLAST and NemaBrowse. The NemaPath metabolic pathway viewer builds associations using KOs, rather than ECs to provide more accurate and fine-grained descriptions of proteins. Two new features for data analysis and comparative genomics have been added to the site. NemaSNP enables the user to perform population genetics studies in various nematode populations using next-generation sequencing data. HelmCoP (Helminth Control and Prevention) as an independent component of Nematode.net provides an integrated resource for storage, annotation and comparative genomics of helminth genomes to aid in learning more about nematode genomes, as well as drug, pesticide, vaccine and drug target discovery. With this update, Nematode.net will continue to realize its original goal to disseminate diverse bioinformatic data sets and provide analysis tools to the broad scientific community in a useful and user-friendly manner.

Targeting protein-protein interactions for parasite control

Finding new drug targets for pathogenic infections would be of great utility for humanity, as there is a large need to develop new drugs to fight infections due to the developing resistance and side effects of current treatments. Current drug targets for pathogen infections involve only a single protein. However, proteins rarely act in isolation, and the majority of biological processes occur via interactions with other proteins, so protein-protein interactions (PPIs) offer a realm of unexplored potential drug targets and are thought to be the next-generation of drug targets. Parasitic worms were chosen for this study because they have deleterious effects on human health, livestock, and plants, costing society billions of dollars annually and many sequenced genomes are available. In this study, we present a computational approach that utilizes whole genomes of 6 parasitic and 1 free-living worm species and 2 hosts. The species were placed in orthologous groups, then binned in species-specific orthologous groups. Proteins that are essential and conserved among species that span a phyla are of greatest value, as they provide foundations for developing broad-control strategies. Two PPI databases were used to find PPIs within the species specific bins. PPIs with unique helminth proteins and helminth proteins with unique features relative to the host, such as indels, were prioritized as drug targets. The PPIs were scored based on RNAi phenotype and homology to the PDB (Protein DataBank). EST data for the various life stages, GO annotation, and druggability were also taken into consideration. Several PPIs emerged from this study as potential drug targets. A few interactions were supported by co-localization of expression in M. incognita (plant parasite) and B. malayi (H. sapiens parasite), which have extremely different modes of parasitism. As more genomes of pathogens are sequenced and PPI databases expanded, this methodology will become increasingly applicable.

Gender-associated genes in filarial nematodes are important for reproduction and potential intervention targets

Background

A better understanding of reproductive processes in parasitic nematodes may lead to development of new anthelmintics and control strategies for combating disabling and disfiguring neglected tropical diseases such as lymphatic filariasis and onchocerciasis. Transcriptomatic analysis has provided important new insights into mechanisms of reproduction and development in other invertebrates. We have performed the first genome-wide analysis of gender-associated (GA) gene expression in a filarial nematode to improve understanding of key reproductive processes in these parasites.

Methodology/Principal Findings

The Version 2 Filarial Microarray with 18,104 elements representing ∼85% of the filarial genome was used to identify GA gene transcripts in adult Brugia malayi worms. Approximately 19% of 14,293 genes were identified as GA genes. Many GA genes have potential Caenorhabditis elegans homologues annotated as germline-, oogenesis-, spermatogenesis-, and early embryogenesis- enriched. The potential C. elegans homologues of the filarial GA genes have a higher frequency of severe RNAi phenotypes (such as lethal and sterility) than other C. elegans genes. Molecular functions and biological processes associated with GA genes were gender-segregated. Peptidase, ligase, transferase, regulator activity for kinase and transcription, and rRNA and lipid binding were associated with female GA genes. In contrast, catalytic activity from kinase, ATP, and carbohydrate binding were associated with male GA genes. Cell cycle, transcription, translation, and biological regulation were increased in females, whereas metabolic processes of phosphate and carbohydrate metabolism, energy generation, and cell communication were increased in males. Significantly enriched pathways in females were associated with cell growth and protein synthesis, whereas metabolic pathways such as pentose phosphate and energy production pathways were enriched in males. There were also striking gender differences in environmental information processing and cell communication pathways. Many proteins encoded by GA genes are secreted by Brugia malayi, and these encode immunomodulatory molecules such as antioxidants and host cytokine mimics. Expression of many GA genes has been recently reported to be suppressed by tetracycline, which blocks reproduction in female Brugia malayi. Our localization of GA transcripts in filarial reproductive organs supports the hypothesis that these genes encode proteins involved in reproduction.

Conclusions/Significance

Genome-wide expression profiling coupled with a robust bioinformatics analysis has greatly expanded our understanding of the molecular biology of reproduction in filarial nematodes. This study has highlighted key molecules and pathways associated with reproductive and other biological processes and identified numerous potential candidates for rational drug design to target reproductive processes.

Lymphatic filariasis is a neglected tropical disease that is caused by thread-like parasitic worms that live and reproduce in lymphatic vessels of the human host. There are no vaccines to prevent filariasis, and available drugs are not effective against all stages of the parasite. In addition, recent reports suggest that the filarial nematodes may be developing resistance to key medications. Therefore, there is an urgent need to identify new drug targets in filarial worms. The purpose of this study was to perform a genome-wide analysis of gender-associated gene transcription to improve understanding of key reproductive processes in filarial nematodes. Our results indicate that thousands of genes are differentially expressed in male and female adult worms. Many of those genes are involved in specific reproductive processes such as embryogenesis and spermatogenesis. In addition, expression of some of those genes is suppressed by tetracycline, a drug that leads to sterilization of adult female worms in many filarial species. Thus, gender-associated genes represent priority targets for design of vaccines and drugs that interfere with reproduction of filarial nematodes. Additional work with this type of integrated systems biology approach should lead to important new tools for controlling filarial diseases.

Novel venom gene discovery in the platypus

BACKGROUND

To date, few peptides in the complex mixture of platypus venom have been identified and sequenced, in part due to the limited amounts of platypus venom available to study. We have constructed and sequenced a cDNA library from an active platypus venom gland to identify the remaining components.

RESULTS

We identified 83 novel putative platypus venom genes from 13 toxin families, which are homologous to known toxins from a wide range of vertebrates (fish, reptiles, insectivores) and invertebrates (spiders, sea anemones, starfish). A number of these are expressed in tissues other than the venom gland, and at least three of these families (those with homology to toxins from distant invertebrates) may play non-toxin roles. Thus, further functional testing is required to confirm venom activity. However, the presence of similar putative toxins in such widely divergent species provides further evidence for the hypothesis that there are certain protein families that are selected preferentially during evolution to become venom peptides. We have also used homology with known proteins to speculate on the contributions of each venom component to the symptoms of platypus envenomation.

CONCLUSIONS

This study represents a step towards fully characterizing the first mammal venom transcriptome. We have found similarities between putative platypus toxins and those of a number of unrelated species, providing insight into the evolution of mammalian venom.

A practical, bioinformatic workflow system for large data sets generated by next-generation sequencing

Transcriptomics (at the level of single cells, tissues and/or whole organisms) underpins many fields of biomedical science, from understanding the basic cellular function in model organisms, to the elucidation of the biological events that govern the development and progression of human diseases, and the exploration of the mechanisms of survival, drug-resistance and virulence of pathogens. Next-generation sequencing (NGS) technologies are contributing to a massive expansion of transcriptomics in all fields and are reducing the cost, time and performance barriers presented by conventional approaches. However, bioinformatic tools for the analysis of the sequence data sets produced by these technologies can be daunting to researchers with limited or no expertise in bioinformatics. Here, we constructed a semi-automated, bioinformatic workflow system, and critically evaluated it for the analysis and annotation of large-scale sequence data sets generated by NGS. We demonstrated its utility for the exploration of differences in the transcriptomes among various stages and both sexes of an economically important parasitic worm (Oesophagostomum dentatum) as well as the prediction and prioritization of essential molecules (including GTPases, protein kinases and phosphatases) as novel drug target candidates. This workflow system provides a practical tool for the assembly, annotation and analysis of NGS data sets, also to researchers with a limited bioinformatic expertise. The custom-written Perl, Python and Unix shell computer scripts used can be readily modified or adapted to suit many different applications. This system is now utilized routinely for the analysis of data sets from pathogens of major socio-economic importance and can, in principle, be applied to transcriptomics data sets from any organism.

Characterizing Ancylostoma caninum transcriptome and exploring nematode parasitic adaptation

BACKGROUND

Hookworm infection is one of the most important neglected diseases in developing countries, with approximately 1 billion people infected worldwide. To better understand hookworm biology and nematode parasitism, the present study generated a near complete transcriptome of the canine hookworm Ancylostoma caninum to a very high coverage using high throughput technology, and compared it to those of the free-living nematode Caenorhabditis elegans and the parasite Brugia malayi.

RESULTS

The generated transcripts from four developmental stages, infective L3, serum stimulated L3, adult male and adult female, covered 93% of the A. caninum transcriptome. The broad diversity among nematode transcriptomes was confirmed, and an impact of parasitic adaptation on transcriptome diversity was inferred. Intra-population analysis showed that A. caninum has higher coding sequence diversity than humans. Examining the developmental expression profiles of A. caninum revealed major transitions in gene expression from larval stages to adult. Adult males expressed the highest number of selectively expressed genes, but adult female expressed the highest number of selective parasitism-related genes. Genes related to parasitism adaptation and A. caninum specific genes exhibited more expression selectivity while those conserved in nematodes tend to be consistently expressed. Parasitism related genes were expressed more selectively in adult male and female worms. The comprehensive analysis of digital expression profiles along with transcriptome comparisons enabled identification of a set of parasitism genes encoding secretory proteins in animal parasitic nematode.

CONCLUSIONS

This study validated the usage of deep sequencing for gene expression profiling. Parasitic adaptation of the canine hookworm is related to its diversity and developmental dynamics. This comprehensive comparative genomic and expression study substantially improves our understanding of the basic biology and parasitism of hookworms and, is expected, in the long run, to accelerate research toward development of vaccines and novel anthelmintics.

Transcripts analysis of the entomopathogenic nematode Steinernema carpocapsae induced in vitro with insect haemolymph

Steinernema carpocapsae is an insect parasitic nematode widely used in pest control programs. The efficacy of this nematode in controlling insects has been found to be related to the pathogenicity of the infective stage. In order to study the parasitic mechanisms exhibited by this parasite, a cDNA library of the induced S. carpocapsae parasitic phase was generated. A total of 2500 clones were sequenced and 2180 high-quality ESTs were obtained from this library. Cluster analysis generated a total of 1592 unique sequences including 1393 singletons. About 63% of the unique sequences had significant hits (e</=1e-05) to the non-redundant protein database. The remaining sequences most likely represent putative novel protein coding genes. Comparative analysis identified 377 homologs in C. elegans, 431 in C. briggsae and 75 in other nematodes. Classification of the predicted proteins revealed involvement in diverse cellular, metabolic and extracellular functions. One hundred and nineteen clusters were predicted to encode putative secreted proteins such as proteases, proteases inhibitors, lectins, saposin-like proteins, acetylcholinesterase, anti-oxidants, and heat-shock proteins, which can possibly have host interactions. This dataset provides a basis for genomic studies towards a better understanding of the events that occur in the parasitic process of this entomopathogenic nematode, including invasion of the insect haemocoelium, adaptations to insect innate immunity and stress responses, and production of virulence factors. The identification of key genes in the parasitic process provides useful tools for the improvement of S. carpocapsae as a biological agent.

Molecular determinants archetypical to the phylum Nematoda

Background

Nematoda diverged from other animals between 600–1,200 million years ago and has become one of the most diverse animal phyla on earth. Most nematodes are free-living animals, but many are parasites of plants and animals including humans, posing major ecological and economical challenges around the world.

Results

We investigated phylum-specific molecular characteristics in Nematoda by exploring over 214,000 polypeptides from 32 nematode species including 27 parasites. Over 50,000 nematode protein families were identified based on primary sequence, including ~10% with members from at least three different species. Nearly 1,600 of the multi-species families did not share homology to Pfam domains, including a total of 758 restricted to Nematoda. Majority of the 462 families that were conserved among both free-living and parasitic species contained members from multiple nematode clades, yet ~90% of the 296 parasite-specific families originated only from a single clade. Features of these protein families were revealed through extrapolation of essential functions from observed RNAi phenotypes in C. elegans, bioinformatics-based functional annotations, identification of distant homology based on protein folds, and prediction of expression at accessible nematode surfaces. In addition, we identified a group of nematode-restricted sequence features in energy-generating electron transfer complexes as potential targets for new chemicals with minimal or no toxicity to the host.

Conclusion

This study identified and characterized the molecular determinants that help in defining the phylum Nematoda, and therefore improved our understanding of nematode protein evolution and provided novel insights for the development of next generation parasite control strategies.

The transcriptomes of the cattle parasitic nematode Ostertagia ostartagi

Ostertagia ostertagi is a gastrointestinal parasitic nematode that affects cattle and leads to a loss of production. In this study, we present the first large-scale genomic survey of O. ostertagi by the analysis of expressed transcripts from three stages of the parasite: third-stage larvae, fourth-stage larvae and adult worms. Using an in silico approach, 2284 genes were identified from over 7000 expressed sequence tags and abundant transcripts were analyzed and characterized by their functional profile. Of the 2284 genes, 66% had similarity to other known or predicted genes while the rest were novel and potentially represent genes specific to the species and/or stages. Furthermore, a subset of the novel proteins were structurally annotated and assigned putative function based on orthologs in Caenorhabditis elegans and corresponding RNA interference phenotypes. Hence, over 70% of the genes were annotated using protein sequences, domains and pathway databases. Differentially expressed transcripts from the two larval stages and their functional profiles were also studied leading to a more detailed understanding of the parasite's life-cycle. The identified transcripts are a valuable resource for genomic studies of O. ostertagi and can facilitate the design of control strategies and vaccine programs.

Systematic analysis of insertions and deletions specific to nematode proteins and their proposed functional and evolutionary relevance

Background

Amino acid insertions and deletions in proteins are considered relatively rare events, and their associations with the evolution and adaptation of organisms are not yet understood. In this study, we undertook a systematic analysis of over 214,000 polypeptides from 32 nematode species and identified insertions and deletions unique to nematode proteins in more than 1000 families and provided indirect evidence that these alterations are linked to the evolution and adaptation of nematodes.

Results

Amino acid alterations in sequences of nematodes were identified by comparison with homologous sequences from a wide range of eukaryotic (metzoan) organisms. This comparison revealed that the proteins inferred from transcriptomic datasets for nematodes contained more deletions than insertions, and that the deletions tended to be larger in length than insertions, indicating a decreased size of the transcriptome of nematodes compared with other organisms. The present findings showed that this reduction is more pronounced in parasitic nematodes compared with the free-living nematodes of the genus Caenorhabditis. Consistent with a requirement for conservation in proteins involved in the processing of genetic information, fewer insertions and deletions were detected in such proteins. On the other hand, more insertions and deletions were recorded for proteins inferred to be involved in the endocrine and immune systems, suggesting a link with adaptation. Similarly, proteins involved in multiple cellular pathways tended to display more deletions and insertions than those involved in a single pathway. The number of insertions and deletions shared by a range of plant parasitic nematodes were higher for proteins involved in lipid metabolism and electron transport compared with other nematodes, suggesting an association between metabolic adaptation and parasitism in plant hosts. We also identified three sizable deletions from proteins found to be specific to and shared by parasitic nematodes, which, given their uniqueness, might serve as target candidates for drug design.

Conclusion

This study illustrates the significance of using comparative genomics approaches to identify molecular elements unique to parasitic nematodes, which have adapted to a particular host organism and mode of existence during evolution. While the focus of this study was on nematodes, the approach has applicability to a wide range of other groups of organisms.

NemaPath: online exploration of KEGG-based metabolic pathways for nematodes

Background

Nematode.net is a web-accessible resource for investigating gene sequences from parasitic and free-living nematode genomes. Beyond the well-characterized model nematode C. elegans, over 500,000 expressed sequence tags (ESTs) and nearly 600,000 genome survey sequences (GSSs) have been generated from 36 nematode species as part of the Parasitic Nematode Genomics Program undertaken by the Genome Center at Washington University School of Medicine. However, these sequencing data are not present in most publicly available protein databases, which only include sequences in Swiss-Prot. Swiss-Prot, in turn, relies on GenBank/Embl/DDJP for predicted proteins from complete genomes or full-length proteins.

Description

Here we present the NemaPath pathway server, a web-based pathway-level visualization tool for navigating putative metabolic pathways for over 30 nematode species, including 27 parasites. The NemaPath approach consists of two parts: 1) a backend tool to align and evaluate nematode genomic sequences (curated EST contigs) against the annotated Kyoto Encyclopedia of Genes and Genomes (KEGG) protein database; 2) a web viewing application that displays annotated KEGG pathway maps based on desired confidence levels of primary sequence similarity as defined by a user. NemaPath also provides cross-referenced access to nematode genome information provided by other tools available on Nematode.net, including: detailed NemaGene EST cluster information; putative translations; GBrowse EST cluster views; links from nematode data to external databases for corresponding synonymous C. elegans counterparts, subject matches in KEGG's gene database, and also KEGG Ontology (KO) identification.

Conclusion

The NemaPath server hosts metabolic pathway mappings for 30 nematode species and is available on the World Wide Web at . The nematode source sequences used for the metabolic pathway mappings are available via FTP , as provided by the Genome Center at Washington University School of Medicine.

Nematode.net update 2008: improvements enabling more efficient data mining and comparative nematode genomics

Nematode.net (http://nematode.net) is a publicly available resource dedicated to the study of parasitic nematodes. In 2000, the Genome Center at Washington University (GC) joined a consortium including the Nematode Genomics group in Edinburgh, and the Pathogen Sequencing Unit of the Sanger Institute to generate expressed sequence tags (ESTs) as an inexpensive and efficient solution for gene discovery in parasitic nematodes. As of 2008 the GC, sampling key parasites of humans, animals and plants, has generated over 500 000 ESTs and 1.2 million genome survey sequences from more than 30 non-Caenorhabditis elegans nematodes. Nematode.net was implemented to offer user-friendly access to data produced by this project. In addition to sequence data, the site hosts: assembled NemaGene clusters in GBrowse views characterizing composition and protein homology, functional Gene Ontology annotations presented via the AmiGO browser, KEGG-based graphical display of NemaGene clusters mapped to metabolic pathways, codon usage tables, NemFam protein families which represent conserved nematode-restricted coding sequences not found in public protein databases, a web-based WU-BLAST search tool that allows complex querying and other assorted resources. The primary aim of Nematode.net is the dissemination of this diverse collection of information to the broader scientific community in a way that is useful, consistent, centralized and enduring.

Intestinal transcriptomes of nematodes: comparison of the parasites Ascaris suum and Haemonchus contortus with the free-living Caenorhabditis elegans

Background

The nematode intestine is a major organ responsible for nutrient digestion and absorption; it is also involved in many other processes, such as reproduction, innate immunity, stress responses, and aging. The importance of the intestine as a target for the control of parasitic nematodes has been demonstrated. However, the lack of detailed knowledge on the molecular and cellular functions of the intestine and the level of its conservation across nematodes has impeded breakthroughs in this application.

Methods and Findings

As part of an extensive effort to investigate various transcribed genomes from Ascaris suum and Haemonchus contortus, we generated a large collection of intestinal sequences from parasitic nematodes by identifying 3,121 A. suum and 1,755 H. contortus genes expressed in the adult intestine through the generation of expressed sequence tags. Cross-species comparisons to the intestine of the free-living C. elegans revealed substantial diversification in the adult intestinal transcriptomes among these species, suggesting lineage- or species-specific adaptations during nematode evolution. In contrast, significant conservation of the intestinal gene repertories was also evident, despite the evolutionary distance of ∼350 million years separating them. A group of 241 intestinal protein families (IntFam-241), each containing members from all three species, was identified based on sequence similarities. These conserved proteins accounted for ∼20% of the sampled intestinal transcriptomes from the three nematodes and are proposed to represent conserved core functions in the nematode intestine. Functional characterizations of the IntFam-241 suggested important roles in molecular functions such as protein kinases and proteases, and biological pathways of carbohydrate metabolism, energy metabolism, and translation. Conservation in the core protein families was further explored by extrapolating observable RNA interference phenotypes in C. elegans to their parasitic counterparts.

Conclusions

Our study has provided novel insights into the nematode intestine and lays foundations for further comparative studies on biology, parasitism, and evolution within the phylum Nematoda.

Biological properties of the nematode intestine warrant in-depth investigation, the results of which can be utilized in the control of parasitic nematodes that infect humans, livestock, and plants. Both the importance of intestinal antigens from Haemonchus contortus in immunity and the damage to H. contortus intestine by anthelmintic fenbendazole have highlighted the versatility of the intestine as an emerging target. However, biological information regarding fundamental intestinal cell functions and mechanisms is currently limited. Conserved intestinal genes across nematode pathogens could offer molecular targets for broad parasite control. Furthermore, qualitative and quantitative comparisons on intestinal gene expression among species and lineages can identify basic adaptations relative to a critical selective force, the nutrient acquisition. This study begins to identify intestinal cell characteristics that are conserved across representatives of two clades of nematodes (V and III) and further clarifies diversities that likely reflect species- or lineage-specific adaptations. Results consistent with functional data on digestive enzymes from H. contortus and RNAi in Caenorhabditis elegans, as examples, support the potential for the comparative genomics approach to produce practical applications. This study provides a platform on which extensive investigation of intestinal genes and a more comprehensive understanding of the Nematoda can be gained.

The canine hookworm genome: analysis and classification of Ancylostoma caninum survey sequences

Hookworms infect nearly a billion people. The Ancylostoma caninum hookworm of canids is a model for studying human infections and information from its genome coupled with functional genomics and proteomics can accelerate progress towards hookworm control. As a step towards a full-scale A. caninum genome project, we generated 104,000 genome survey sequences (GSSs) and determined the genome size of the canine hookworm. GSSs assembled into 57.6 Mb of unique sequence from a genome that we estimate by flow cytometry of isolated nuclei to be 347 +/- 1.2 Mb, substantially larger than other Rhabditina species. Gene finding identified 5538 genes in the GSS assembly, for a total of 9113 non-redundant A. caninum genes when EST sequences are also considered. Functional classifications of many of the 70% of genes with homology to genes in other species are provided based on gene ontology and KEGG associations and secreted and membrane-bound proteins are also identified.