RIEMS: a software pipeline for sensitive and comprehensive taxonomic classification of reads from metagenomics datasets

Background

Fuelled by the advent and subsequent development of next generation sequencing technologies, metagenomics became a powerful tool for the analysis of microbial communities both scientifically and diagnostically. The biggest challenge is the extraction of relevant information from the huge sequence datasets generated for metagenomics studies. Although a plethora of tools are available, data analysis is still a bottleneck.

Results

To overcome the bottleneck of data analysis, we developed an automated computational workflow called RIEMS – Reliable Information Extraction from Metagenomic Sequence datasets. RIEMS assigns every individual read sequence within a dataset taxonomically by cascading different sequence analyses with decreasing stringency of the assignments using various software applications. After completion of the analyses, the results are summarised in a clearly structured result protocol organised taxonomically. The high accuracy and performance of RIEMS analyses were proven in comparison with other tools for metagenomics data analysis using simulated sequencing read datasets.

Conclusions

RIEMS has the potential to fill the gap that still exists with regard to data analysis for metagenomics studies. The usefulness and power of RIEMS for the analysis of genuine sequencing datasets was demonstrated with an early version of RIEMS in 2011 when it was used to detect the orthobunyavirus sequences leading to the discovery of Schmallenberg virus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0503-6) contains supplementary material, which is available to authorized users.     

The impact of genetic diversity in protozoa on molecular diagnostics

Detection of intestinal parasitic protists, commonly referred to as 'intestinal protozoa,' by PCR is increasingly used not only for identification or confirmation but also as a first-line diagnostic tool. Apart from the ability to sample correctly and extract parasite DNA directly from faeces, primer and probe specificity and sensitivity affect predictive values and hence the utility of diagnostic assays. Molecular characterization of intestinal protists is necessary to design primers and probes because this is the basic material for current and future improved diagnostic PCRs for either detecting all genetic variants or specifically differentiating among such variants. As an example, this paper highlights the existence of interspecific and intraspecific genetic diversity among intestinal, unicellular parasites and its implications for nucleic acid-based diagnostic assays.     

Integrating metagenomic and amplicon databases to resolve the phylogenetic and ecological diversity of the Chlamydiae

In the era of metagenomics and amplicon sequencing, comprehensive analyses of available sequence data remain a challenge. Here we describe an approach exploiting metagenomic and amplicon data sets from public databases to elucidate phylogenetic diversity of defined microbial taxa. We investigated the phylum Chlamydiae whose known members are obligate intracellular bacteria that represent important pathogens of humans and animals, as well as symbionts of protists. Despite their medical relevance, our knowledge about chlamydial diversity is still scarce. Most of the nine known families are represented by only a few isolates, while previous clone library-based surveys suggested the existence of yet uncharacterized members of this phylum. Here we identified more than 22 000 high quality, non-redundant chlamydial 16S rRNA gene sequences in diverse databases, as well as 1900 putative chlamydial protein-encoding genes. Even when applying the most conservative approach, clustering of chlamydial 16S rRNA gene sequences into operational taxonomic units revealed an unexpectedly high species, genus and family-level diversity within the Chlamydiae, including 181 putative families. These in silico findings were verified experimentally in one Antarctic sample, which contained a high diversity of novel Chlamydiae. In our analysis, the Rhabdochlamydiaceae, whose known members infect arthropods, represents the most diverse and species-rich chlamydial family, followed by the protist-associated Parachlamydiaceae, and a putative new family (PCF8) with unknown host specificity. Available information on the origin of metagenomic samples indicated that marine environments contain the majority of the newly discovered chlamydial lineages, highlighting this environment as an important chlamydial reservoir.     

Metaxa: a software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea,bacteria, eukaryotes,mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets

The ribosomal small subunit (SSU) rRNA gene has emerged as an important genetic marker for taxonomic identification in environmental sequencing datasets. In addition to being present in the nucleus of eukaryotes and the core genome of prokaryotes, the gene is also found in the mitochondria of eukaryotes and in the chloroplasts of photosynthetic eukaryotes. These three sets of genes are conceptually paralogous and should in most situations not be aligned and analyzed jointly. To identify the origin of SSU sequences in complex sequence datasets has hitherto been a time-consuming and largely manual undertaking. However, the present study introduces Metaxa (), an automated software tool to extract full-length and partial SSU sequences from larger sequence datasets and assign them to an archaeal, bacterial, nuclear eukaryote, mitochondrial, or chloroplast origin. Using data from reference databases and from full-length organelle and organism genomes, we show that Metaxa detects and scores SSU sequences for origin with very low proportions of false positives and negatives. We believe that this tool will be useful in microbial and evolutionary ecology as well as in metagenomics.     

Intestinal microsporidiosis in India: A two year study

Intestinal parasitic pathogens in HIV/AIDS patients include Cryptosporidium sp, Cystoisospora sp, microsporidia and less commonly other parasites. The two most common microsporidia causing intestinal infection are Enterocytozoon bieneusi and Encephalitozoon intestinalis. Most of the Indian studies for intestinal parasitic infections in HIV/AIDS patients have not included microsporidia, due to difficult staining and identification of the parasite. The aim of the present study was to find the prevalence of intestinal microsporidiosis and their species identification along with correlation of CD4 count with parasite positivity and diarrhoea in HIV positive individuals. Stool samples of 363 individuals including 125 HIV seropositive patients with diarrhoea, 158 HIV seropositive patients without diarrhoea, 55 HIV seronegative patients with diarrhoea and 25 healthy controls were obtained from various out-patient departments and in-patients admitted to a tertiary care hospital from August 2008 to October 2009. The stool samples were subjected to examination by wet mount, modified acid fast stain for coccidian parasites and multiplex nested PCR for microsporidia. The overall prevalence of all intestinal parasites among HIV patients in our study was 26.5%. The prevalence of intestinal parasitic pathogens in HIV positive patients with diarrhoea was 43.2%. Microsporidia were the most common parasites detected (14%) in all patients, while in HIV infected patients 15.9% patients had microsporidia infection. The most common species causing intestinal microsporidiosis in our study was E. intestinalis (10.5%). In HIV seropositive individuals with diarrhoea, E. intestinalis was 20.8% and E. bieneusi 8.0% while in HIV-seropositive individuals without diarrhoea, E. intestinalis was 3.8% and E. bieneusi 1.9%. E. intestinalis was present in 10.9% of HIV negative individuals with diarrhoea in whom E. bieneusi was not found. There was a significant association between CD4 count ≤ 200/μl and intestinal parasite positivity. Thus, it can be concluded that intestinal microsporidiosis is under reported but an important disease in India. The predominant species in our study is E. intestinalis , in contrast to other parts of the world where E. bieneusi is more common.     

Molecular phylogenetic evidence that the phylum Haplosporidia has an alveolate ancestry

The phylogenetic position of the phylum Haplosporidia among other protists was investigated with the complete 16S-like rRNA gene sequences from two species in the phylum: Haplosporidium nelsoni, a parasite of oysters, and Minchinia teredinis, a parasite of shipworms. Because the lack of obvious morphological homologies with other protists hampered decisions regarding taxonomic composition for sequence alignment and phylogenetic analysis, the complete sequences for these two haplosporidians were directed as search queries to the blast/ncbi.nlm.nih.gov electronic mail server. The results of this heuristic similarity search provided a basis for constructing a preliminary higher-taxonomic-level analysis comparing the haplosporidians with species from the slime molds, fungi, algae, amoebae, ciliates, dinoflagellates, and apicomplexans. Maximum parsimony yielded equivocal results, whereas transversionally weighted parsimony suggested an affinity with the alveolates (i.e., the ciliates, dinoflagellates, and apicomplexans). Multiple alignment of the two haplosporidian sequences against 17 taxa in a secondary analysis focusing on the alveolates and subsequent parsimony analysis placed the phylum Haplosporidia as a monophyletic group within the Alveolata and as a taxon of equal rank with the other three alveolate phyla. The precise placement within the Alveolata was sensitive to weighting.      

DNA from pre-erythrocytic stage malaria parasites is detectable by PCR in the faeces and blood of hosts

Following the bite of an infective mosquito, malaria parasites first invade the liver where they develop and replicate for a number of days before being released into the bloodstream where they invade red blood cells and cause disease. The biology of the liver stages of malaria parasites is relatively poorly understood due to the inaccessibility of the parasites to sampling during this phase of their life cycle. Here we report the detection in blood and faecal samples of malaria parasite DNA throughout their development in the livers of mice and before the parasites begin their growth in the blood circulation. It is shown that parasite DNA derived from pre-erythrocytic stage parasites reaches the faeces via the bile. We then show that different primate malaria species can be detected by PCR in blood and faecal samples from naturally infected captive macaque monkeys. These results demonstrate that pre-erythrocytic parasites can be detected and quantified in experimentally infected animals. Furthermore, these results have important implications for both molecular epidemiology and phylogenetics of malaria parasites. In the former case, individuals who are malaria parasite negative by microscopy, but PCR positive for parasite DNA in their blood, are considered to be “sub-microscopic” blood stage parasite carriers. We now propose that PCR positivity is not necessarily an indicator of the presence of blood stage parasites, as the DNA could derive from pre-erythrocytic parasites. Similarly, in the case of molecular phylogenetics based on DNA sequences alone, we argue that DNA amplified from blood or faeces does not necessarily come from a parasite species that infects the red blood cells of that particular host.     

Detection and discovery of crustacean parasites in blue crabs (Callinectes sapidus) by using 18S rRNA gene-targeted denaturing high-performance liquid chromatography

Recently, we described a novel denaturing high-performance liquid chromatography (DHPLC) approach useful for initial detection and identification of crustacean parasites. Because this approach utilizes general primers targeted to conserved regions of the 18S rRNA gene, a priori genetic sequence information on eukaryotic parasites is not required. This distinction provides a significant advantage over specifically targeted PCR assays that do not allow for the detection of unknown or unsuspected parasites. However, initial field evaluations of the DHPLC assay suggested that because of PCR-biased amplification of dominant host genes it was not possible to detect relatively rare parasite genes in infected crab tissue. Here, we describe the use of a peptide nucleic acid (PNA) PCR hybridization blocking probe in association with DHPLC (PNA-PCR DHPLC) to overcome inherent PCR bias associated with amplification of rare target genes by use of generic primers. This approach was utilized to detect infection of blue crabs (Callinectes sapidus) by the parasitic dinoflagellate Hematodinium sp. Evaluation of 76 crabs caught in Wassaw Sound, GA, indicated a 97% correspondence between detection of the parasite by use of a specific PCR diagnostic assay and that by use of PNA-PCR DHPLC. During these studies, we discovered one crab with an association with a previously undescribed protist symbiont. Phylogenetic analysis of the amplified symbiont 18S rRNA gene indicated that it is most closely related to the free-living kinetoplastid parasite Procryptobia sorokini. To our knowledge, this is the first report of this parasite group in a decapod crab and of this organism exhibiting a presumably parasitic life history.     

Loss of forest cover and host functional diversity increases prevalence of avian malaria parasites in the Atlantic Forest

Host phylogenetic relatedness and ecological similarity are thought to contribute to parasite community assembly and infection rates. However, recent landscape level anthropogenic changes may disrupt host-parasite systems by impacting functional and phylogenetic diversity of host communities. We examined whether changes in host functional and phylogenetic diversity, forest cover, and minimum temperature influence the prevalence, diversity, and distributions of avian haemosporidian parasites (genera Haemoproteus and Plasmodium) across 18 avian communities in the Atlantic Forest. To explore spatial patterns in avian haemosporidian prevalence and taxonomic and phylogenetic diversity, we surveyed 2241 individuals belonging to 233 avian species across a deforestation gradient. Mean prevalence and parasite diversity varied considerably across avian communities and parasites responded differently to host attributes and anthropogenic changes. Avian malaria prevalence (termed herein as an infection caused by Plasmodium parasites) was higher in deforested sites, and both Plasmodium prevalence and taxonomic diversity were negatively related to host functional diversity. Increased diversity of avian hosts increased local taxonomic diversity of Plasmodium lineages but decreased phylogenetic diversity of this parasite genus. Temperature and host phylogenetic diversity did not influence prevalence and diversity of haemosporidian parasites. Variation in the diversity of avian host traits that promote parasite encounter and vector exposure (host functional diversity) partially explained the variation in avian malaria prevalence and diversity. Recent anthropogenic landscape transformation (reduced proportion of native forest cover) had a major influence on avian malaria occurrence across the Atlantic Forest. This suggests that, for Plasmodium, host phylogenetic diversity was not a biotic filter to parasite transmission as prevalence was largely explained by host ecological attributes and recent anthropogenic factors. Our results demonstrate that, similar to human malaria and other vector-transmitted pathogens, prevalence of avian malaria parasites will likely increase with deforestation.     

Hypervariable regions (V1–V9) of the dinoflagellate 18S rRNA using a large dataset for marker considerations

DNA sequencing methods have been used for the molecular taxonomic discrimination of dinoflagellate protists, particularly using partial 18S rRNA sequences. This study evaluated the taxonomic discrimination power of rRNA gene hypervariable regions (V1 to V9) in dinoflagellates from a large dataset. These included 77 dinoflagellate species (9 orders, 17 families, 40 genera). The complete 18S rRNA sequences of the dinoflagellates ranged from 1,787 to 1,813?bp in length, and consisted of eight V regions with a total combined length of 678 to 699?bp. Regions longer than 100?bp were recoded for V2, V4, and V8 regions; high nucleotide divergences were detected in V1, V2, and V4 regions. Statistic tests showed that the divergences of individual V regions were significantly different (t-test, P?<?0.05) compared with the complete 18S rRNA. The V2 region showed the highest score (83.5%) for PI sites. Moreover, intra-genus DNA similarities of the V2 were considerably low (<93%). Neighbor-joining analyses showed that phylogenetic resolution in the V2–V4 region was 1.32-fold higher than that of the complete 18S rRNA. These results demonstrate that V2 has the highest taxonomic resolving power within the 18S rRNA gene of dinoflagellates, suggesting the V2 and adjacent regions (e.g., V1 to V4) may be the best for marker considerations.     

Circadian variation in shedding of the oocysts of Isospora turdi (Apicomplexa) in blackbirds (Turdusmerula): An adaptative trait against desiccation and ultraviolet radiation

Many parasite species spend part of their life cycle in the external environment waiting for a new host. Emergence of parasites often occurs once a day, which may help to minimise mortality in an inhospitable environment and increase transition rates. Many intestinal parasites in birds are released in faeces only in the late afternoon. However, the adaptative significance of this pattern is unclear. One hypothesis is that a particular time of emergence may prevent parasite desiccation and therefore increase the parasite’s life expectancy in the external environment. We tested this hypothesis experimentally using the blackbird (Turdus merula) infected with Isospora turdi (Protozoa: Apicomplexa). We found that short exposure of faeces to natural sunlight has a dramatic effect on oocyst survival. This appears to be due to the effect of warmth and ultraviolet (UV) radiation with UVB waves being more damaging than UVA. Oocysts contained in faeces shed in water are protected from the effect of sunlight. Together, these results suggest that the release of oocysts in the late afternoon is an adaptative trait to avoid desiccation and UV radiation, thus reducing mortality of the oocysts in the external environment.     

Caution on the assessment of intestinal parasitic load in studying parasite-mediated sexual selection: The case of Blackbirds coccidiosis

The parasite-mediated sexual selection (PMSS) theory has led to an increasing number of experimental studies essentially focussed on blood parasites. Currently, more research is being carried out on intestinal parasites in relationship to this theory. Before testing the theory with gastrointestinal parasites, it is important: (i) to determine an optimal research methodology to obtain an accurate assessment of parasite burden and (ii) to have information about life-history traits of the parasite to interpret data appropriately. In this study, we present data on oocyst output of Isosporaturdi in the faeces of blackbirds (Turdus merula) that illustrate the importance of developing methods that are relevant for a particular model system, instead of relying on existing methods that may work in other systems. Our results show that: (i) a single droplet of faeces will accurately indicate the parasitic load in blackbirds, (ii) oocyst shedding varies greatly within and between days, (iii) the course of infection is characterised with two successive peaks of oocyst shedding, (iv) infection lasts approximately 1 month and (iv) there is no effect of sex, size of infective dose or re-infection on the course of infection. We discuss the practical implications of these results in determining the reliability of, and in avoiding erroneous conclusions about, PMSS when using intestinal parasites as models. In particular, we emphasise that numerous measurements must be performed on different days. In addition, faeces must be collected at a particular time of day because there is a strong temporal variation in oocyst shedding. Finally, a standardised methodology that has been developed for a particular host-parasite assemblage does not necessarily work with other biological models. Together, our results should encourage future investigators to identify an accurate methodology for assessing parasitic load as a first step before testing hypotheses associated with the PMSS theory.     

Cryptosporidium: Genomic and biochemical features

Recent progress in understanding the unique biochemistry of the two closely related human enteric pathogens Cryptosporidium parvum and Cryptosporidium hominis has been stimulated by the elucidation of the complete genome sequences for both pathogens. Much of the work that has occurred since that time has been focused on understanding the metabolic pathways encoded by the genome in hopes of providing increased understanding of the parasite biology, and in the identification of novel targets for pharmacological interventions. However, despite identifying the genes encoding enzymes that participate in many of the major metabolic pathways, only a hand full of proteins have actually been the subjects of detailed scrutiny. Thus, much of the biochemistry of these parasites remains a true mystery.     

Molecular Identification of Nanoplanktonic Protists Based on Small Subunit Ribosomal RNA Gene Sequences for Ecological Studies1

Nanoplanktonic protists are comprised of a diverse assemblage of species which are responsible for a variety of trophic processes in marine and freshwater ecosystems. Current methods for identifying small protists by electron microscopy do not readily permit both identification and enumeration of nanoplanktonic protists in field samples. Thus, one major goal in the application of molecular approaches in protistan ecology has been the detection and quantification of individual species in natural water samples. Sequences of small subunit ribosomal RNA (SSU rRNA) genes have proven to be useful towards achieving this goal. Comparison of sequences from clone libraries of protistan SSU rRNA genes amplified from natural assemblages of protists by the polymerase chain reaction (PCR) can be used to examine protistan diversity. Furthermore, oligonucleotide probes complementary to short sequence regions unique to species of small protists can be designed by comparative analysis of rRNA gene sequences. These probes may be used to either detect the RNA of particular species of protists in total nucleic acid extracts immobilized on membranes, or the presence of target species in water samples via in situ hybridization of whole cells. Oligonucleotide probes may also serve as primers for the selective amplification of target sequences from total population DNA by PCR. Thus, molecular sequence information is becoming increasingly useful for identifying and enumerating protists, and for studying their spatial and temporal distribution in nature. Knowledge of protistan species composition, abundance and variability in an environment can ultimately be used to relate community structure to various aspects of community function and biogeochemical activity.     

Phylogenetic position of the genus Gonocerca Manter, 1925 (Trematoda,Hemiuroidea), based on partial sequences of 28S rRNA gene and a reconsideration of taxonomic status of Gonocercinae Skrjabin et Guschanskaja, 1955

Adult trematodes of the genus Gonocerca Manter, 1925, are parasites of marine fishes. Identification of the phylogenetic positions and a revision of the taxonomic status of the subfamily Gonocercinae Skrjabin et Guschanskaja, 1955 (Derogenidae) are the main purposes of this research article. Four Gonocerca species were used in the study, including the type-species G. phycidis Manter, 1925. Molecular phylogenetic analysis, based on partial sequences of 28S rRNA gene, revealed that Gonocerca spp. are phylogenetically distant from other hemiuroid trematodes, including Derogenes varicus (Müller, 1784), representative of the type-genus of the family Derogenidae. The taxonomic rank of Gonocercinae should be raised to the family level. The generic composition of the family Gonocercidae Skrjabin et Guschanskaja, 1955 stat. nov., requires further clarification as the molecular data do not support the inclusion of the genus Hemipera Nicoll, 1913, in this family.     

Intranuclear verrucomicrobial symbionts and evidence of lateral gene transfer to the host protist in the termite gut

In 1944, Harold Kirby described microorganisms living within nuclei of the protists Trichonympha in guts of termites; however, their taxonomic assignment remains to be accomplished. Here, we identified intranuclear symbionts of Trichonympha agilis in the gut of the termite Reticulitermes speratus. We isolated single nuclei of T. agilis, performed whole-genome amplification, and obtained bacterial 16S rRNA genes by PCR. Unexpectedly, however, all of the analyzed clones were from pseudogenes of 16S rRNA with large deletions and numerous sequence variations even within a single-nucleus sample. Authentic 16S rRNA gene sequences were finally recovered by digesting the nuclear DNA; these pseudogenes were present on the host Trichonympha genome. The authentic sequences represented two distinct bacterial species belonging to the phylum Verrucomicrobia, and the pseudogenes have originated from each of the two species. Fluorescence in situ hybridization confirmed that both species are specifically localized, and occasionally co-localized, within nuclei of T. agilis. Transmission electron microscopy revealed that they are distorted cocci with characteristic electron-dense and lucent regions, which resemble the intranuclear symbionts illustrated by Kirby. For these symbionts, we propose a novel genus and species, ‘Candidatus Nucleococcus trichonymphae'' and ‘Candidatus Nucleococcus kirbyi''. These formed a termite-specific cluster with database sequences, other members of which were also detected within nuclei of various gut protists, including both parabasalids and oxymonads. We suggest that this group is widely distributed as intranuclear symbionts of diverse protists in termite guts and that they might have affected the evolution of the host genome through lateral gene transfer.     

Assessment of Helminth Biodiversity in Wild Rats Using 18S rDNA Based Metagenomics

Parasite diversity has important implications in several research fields including ecology, evolutionary biology and epidemiology. Wide-ranging analysis has been restricted because of the difficult, highly specialised and time-consuming processes involved in parasite identification. In this study, we assessed parasite diversity in wild rats using 18S rDNA-based metagenomics. 18S rDNA PCR products were sequenced using an Illumina MiSeq sequencer and the analysis of the sequences using the QIIME software successfully classified them into several parasite groups. The comparison of the results with those obtained using standard methods including microscopic observation of helminth parasites in the rat intestines and PCR amplification/sequencing of 18S rDNA from isolated single worms suggests that this new technique is reliable and useful to investigate parasite diversity.     

GLOSSary: the GLobal Ocean 16S subunit web accessible resource

Background

Environmental metagenomics is a challenging approach that is exponentially spreading in the scientific community to investigate taxonomic diversity and possible functions of the biological components. The massive amount of sequence data produced, often endowed with rich environmental metadata, needs suitable computational tools to fully explore the embedded information. Bioinformatics plays a key role in providing methodologies to manage, process and mine molecular data, integrated with environmental metagenomics collections. One such relevant example is represented by the Tara Ocean Project.

Results

We considered the Tara 16S miTAGs released by the consortium, representing raw sequences from a shotgun metagenomics approach with similarities to 16S rRNA genes. We generated assembled 16S rDNA sequences, which were classified according to their lengths, the possible presence of chimeric reads, the putative taxonomic affiliation. The dataset was included in GLOSSary (the GLobal Ocean 16S Subunit web accessible resource), a bioinformatics platform to organize environmental metagenomics data. The aims of this work were: i) to present alternative computational approaches to manage challenging metagenomics data; ii) to set up user friendly web-based platforms to allow the integration of environmental metagenomics sequences and of the associated metadata; iii) to implement an appropriate bioinformatics platform supporting the analysis of 16S rDNA sequences exploiting reference datasets, such as the SILVA database. We organized the data in a next-generation NoSQL “schema-less” database, allowing flexible organization of large amounts of data and supporting native geospatial queries. A web interface was developed to permit an interactive exploration and a visual geographical localization of the data, either raw miTAG reads or 16S contigs, from our processing pipeline. Information on unassembled sequences is also available. The taxonomic affiliations of contigs and miTAGs, and the spatial distribution of the sampling sites and their associated sequence libraries, as they are contained in the Tara metadata, can be explored by a query interface, which allows both textual and visual investigations. In addition, all the sequence data were made available for a dedicated BLAST-based web application alongside the SILVA collection.

Conclusions

GLOSSary provides an expandable bioinformatics environment, able to support the scientific community in current and forthcoming environmental metagenomics analyses.