Tarini Shankar Ghosh

Classification of metagenomic sequences: methods and challenges

Characterizing the taxonomic diversity of microbial communities is one of the primary objectives of metagenomic studies. Taxonomic analysis of microbial communities, a process referred to as binning, is challenging for the following reasons. Primarily, query sequences originating from the genomes of most microbes in an environmental sample lack taxonomically related sequences in existing reference databases. This absence of a taxonomic context makes binning a very challenging task. Limitations of current sequencing platforms, with respect to short read lengths and sequencing errors/artifacts, are also key factors that determine the overall binning efficiency. Furthermore, the sheer volume of metagenomic datasets also demands highly efficient algorithms that can operate within reasonable requirements of compute power. This review discusses the premise, methodologies, advantages, limitations and challenges of various methods available for binning of metagenomic datasets obtained using the shotgun sequencing approach. Various parameters as well as strategies used for evaluating binning efficiency are then reviewed.

Metagenome of the gut of a malnourished child

BackgroundMalnutrition, a major health problem, affects a significant proportion of preschool children in developing countries. The devastating consequences of malnutrition include diarrhoea, malabsorption, increased intestinal permeability, suboptimal immune response, etc. Nutritional interventions and dietary solutions have not been effective for treatment of malnutrition till date. Metagenomic procedures allow one to access the complex cross-talk between the gut and its microbial flora and understand how a different community composition affects various states of human health. In this study, a metagenomic approach was employed for analysing the differences between gut microbial communities obtained from a malnourished and an apparently healthy child.ResultsOur results indicate that the malnourished child gut has an abundance of enteric pathogens which are known to cause intestinal inflammation resulting in malabsorption of nutrients. We also identified a few functional sub-systems from these pathogens, which probably impact the overall metabolic capabilities of the malnourished child gut.ConclusionThe present study comprehensively characterizes the microbial community resident in the gut of a malnourished child. This study has attempted to extend the understanding of the basis of malnutrition beyond nutrition deprivation.

Computational analysis of the ESX-1 region of Mycobacterium tuberculosis: insights into the mechanism of type VII secretion system.

Type VII secretion system (T7SS) is a recent discovery in bacterial secretion systems. First identified in Mycobacterium tuberculosis, this secretion system has later been reported in organisms belonging to the Actinomycetales order and even to distant phyla like Firmicutes. The genome of M. tuberculosis H37Rv contains five gene clusters that have evolved through gene duplication events and include components of the T7SS secretion machinery. These clusters are called ESAT-6 secretion system (ESX) 1 through 5. Out of these, ESX-1 has been the most widely studied region because of its pathological importance. In spite of this, the overall mechanism of protein translocation through ESX-1 secretion machinery is not clearly understood. Specifically, the structural components contributing to the translocation through the mycomembrane have not been characterized yet. In this study, we have carried out a comprehensive in silico analysis of the genes known to be involved in ESX-1 secretion pathway and identified putative proteins having high probability to be associated with this particular pathway. Our study includes analysis of phylogenetic profiles, identification of domains, transmembrane helices, 3D folds, signal peptides and prediction of protein-protein associations. Based on our analysis, we could assign probable novel functions to a few of the ESX-1 components. Additionally, we have identified a few proteins with probable role in the initial activation and formation of mycomembrane translocon of ESX-1 secretion machinery. We also propose a probable working model of T7SS involving ESX-1 secretion pathway.

In Silico Analysis of Antibiotic Resistance Genes in the Gut Microflora of Individuals from Diverse Geographies and Age-Groups

The spread of antibiotic resistance, originating from the rampant and unrestrictive use of antibiotics in humans and livestock over the past few decades has emerged as a global health problem. This problem has been further compounded by recent reports implicating the gut microbial communities to act as reservoirs of antibiotic resistance. We have profiled the presence of probable antibiotic resistance genes in the gut flora of 275 individuals from eight different nationalities. For this purpose, available metagenomic data sets corresponding to 275 gut microbiomes were analyzed. Sequence similarity searches of the genomic fragments constituting each of these metagenomes were performed against genes conferring resistance to around 240 antibiotics. Potential antibiotic resistance genes conferring resistance against 53 different antibiotics were detected in the human gut microflora analysed in this study. In addition to several geography/country-specific patterns, four distinct clusters of gut microbiomes, referred to as ‘Resistotypes’, exhibiting similarities in their antibiotic resistance profiles, were identified. Groups of antibiotics having similarities in their resistance patterns within each of these clusters were also detected. Apart from this, mobile multi-drug resistance gene operons were detected in certain gut microbiomes. The study highlighted an alarmingly high abundance of antibiotic resistance genes in two infant gut microbiomes. The results obtained in the present study presents a holistic ‘big picture’ on the spectra of antibiotic resistance within our gut microbiota across different geographies. Such insights may help in implementation of new regulations and stringency on the existing ones.

Taxonomic and functional characteristics of microbial communities and their correlation with physicochemical properties of four geothermal springs in Odisha, India

This study describes microbial diversity in four tropical hot springs representing moderately thermophilic environments (temperature range: 40–58°C; pH: 7.2–7.4) with discrete geochemistry. Metagenome sequence data showed a dominance of Bacteria over Archaea; the most abundant phyla were Chloroflexi and Proteobacteria, although other phyla were also present, such as Acetothermia, Nitrospirae, Acidobacteria, Firmicutes, Deinococcus-Thermus, Bacteroidetes, Thermotogae, Euryarchaeota, Verrucomicrobia, Ignavibacteriae, Cyanobacteria, Actinobacteria, Planctomycetes, Spirochaetes, Armatimonadetes, Crenarchaeota, and Aquificae. The distribution of major genera and their statistical correlation analyses with the physicochemical parameters predicted that the temperature, aqueous concentrations of ions (such as sodium, chloride, sulfate, and bicarbonate), total hardness, dissolved solids and conductivity were the main environmental variables influencing microbial community composition and diversity. Despite the observed high taxonomic diversity, there were only little variations in the overall functional profiles of the microbial communities in the four springs. Genes involved in the metabolism of carbohydrates and carbon fixation were the most abundant functional class of genes present in these hot springs. The distribution of genes involved in carbon fixation predicted the presence of all the six known autotrophic pathways in the metagenomes. A high prevalence of genes involved in membrane transport, signal transduction, stress response, bacterial chemotaxis, and flagellar assembly were observed along with genes involved in the pathways of xenobiotic degradation and metabolism. The analysis of the metagenomic sequences affiliated to the candidate phylum Acetothermia from spring TB-3 provided new insight into the metabolism and physiology of yet-unknown members of this lineage of bacteria.

HabiSign: a novel approach for comparison of metagenomes and rapid identification of habitat-specific sequences

BackgroundOne of the primary goals of comparative metagenomic projects is to study the differences in the microbial communities residing in diverse environments. Besides providing valuable insights into the inherent structure of the microbial populations, these studies have potential applications in several important areas of medical research like disease diagnostics, detection of pathogenic contamination and identification of hitherto unknown pathogens. Here we present a novel and rapid, alignment-free method called HabiSign, which utilizes patterns of tetra-nucleotide usage in microbial genomes to bring out the differences in the composition of both diverse and related microbial communities.ResultsValidation results show that the metagenomic signatures obtained using the HabiSign method are able to accurately cluster metagenomes at biome, phenotypic and species levels, as compared to an average tetranucleotide frequency based approach and the recently published dinucleotide relative abundance based approach. More importantly, the method is able to identify subsets of sequences that are specific to a particular habitat. Apart from this, being alignment-free, the method can rapidly compare and group multiple metagenomic data sets in a short span of time.ConclusionsThe proposed method is expected to have immense applicability in diverse areas of metagenomic research ranging from disease diagnostics and pathogen detection to bio-prospecting. A web-server for the HabiSign algorithm is available at http://metagenomics.atc.tcs.com/HabiSign/.

ProViDE: A software tool for accurate estimation of viral diversity in metagenomic samples

Given the absence of universal marker genes in the viral kingdom, researchers typically use BLAST (with stringent E-values) for taxonomic classification of viral metagenomic sequences. Since majority of metagenomic sequences originate from hitherto unknown viral groups, using stringent e-values results in most sequences remaining unclassified. Furthermore, using less stringent e-values results in a high number of incorrect taxonomic assignments. The SOrt-ITEMS algorithm provides an approach to address the above issues. Based on alignment parameters, SOrt-ITEMS follows an elaborate work-flow for assigning reads originating from hitherto unknown archaeal/bacterial genomes. In SOrt-ITEMS, alignment parameter thresholds were generated by observing patterns of sequence divergence within and across various taxonomic groups belonging to bacterial and archaeal kingdoms. However, many taxonomic groups within the viral kingdom lack a typical Linnean-like taxonomic hierarchy. In this paper, we present ProViDE (Program for Viral Diversity Estimation), an algorithm that uses a customized set of alignment parameter thresholds, specifically suited for viral metagenomic sequences. These thresholds capture the pattern of sequence divergence and the non-uniform taxonomic hierarchy observed within/across various taxonomic groups of the viral kingdom. Validation results indicate that the percentage of ‘correct’ assignments by ProViDE is around 1.7 to 3 times higher than that by the widely used similarity based method MEGAN. The misclassification rate of ProViDE is around 3 to 19% (as compared to 5 to 42% by MEGAN) indicating significantly better assignment accuracy. ProViDE software and a supplementary file (containing supplementary figures and tables referred to in this article) is available for download from http://metagenomics.atc.tcs.com/binning/ProViDE/

C16S — A Hidden Markov Model based algorithm for taxonomic classification of 16S rRNA gene sequences

Recent advances in high throughput sequencing technologies and concurrent refinements in 16S rDNA isolation techniques have facilitated the rapid extraction and sequencing of 16S rDNA content of microbial communities. The taxonomic affiliation of these 16S rDNA fragments is subsequently obtained using either BLAST-based or word frequency based approaches. However, the classification accuracy of such methods is observed to be limited in typical metagenomic scenarios, wherein a majority of organisms are hitherto unknown. In this study, we present a 16S rDNA classification algorithm, called C16S, that uses genus-specific Hidden Markov Models for taxonomic classification of 16S rDNA sequences. Results obtained using C16S have been compared with the widely used RDP classifier. The performance of C16S algorithm was observed to be consistently higher than the RDP classifier. In some scenarios, this increase in accuracy is as high as 34%. A web-server for the C16S algorithm is available at http://metagenomics.atc.tcs.com/C16S/.

Eu-Detect: An algorithm for detecting eukaryotic sequences in metagenomic data sets

Physical partitioning techniques are routinely employed (during sample preparation stage) for segregating the prokaryotic and eukaryotic fractions of metagenomic samples. In spite of these efforts, several metagenomic studies focusing on bacterial and archaeal populations have reported the presence of contaminating eukaryotic sequences in metagenomic data sets. Contaminating sequences originate not only from genomes of micro-eukaryotic species but also from genomes of (higher) eukaryotic host cells. The latter scenario usually occurs in the case of host-associated metagenomes. Identification and removal of contaminating sequences is important, since these sequences not only impact estimates of microbial diversity but also affect the accuracy of several downstream analyses. Currently, the computational techniques used for identifying contaminating eukaryotic sequences, being alignment based, are slow, inefficient, and require huge computing resources. In this article, we present Eu-Detect, an alignment-free algorithm that can rapidly identify eukaryotic sequences contaminating metagenomic data sets. Validation results indicate that on a desktop with modest hardware specifications, the Eu-Detect algorithm is able to rapidly segregate DNA sequence fragments of prokaryotic and eukaryotic origin, with high sensitivity. A Web server for the Eu-Detect algorithm is available at http://metagenomics.atc.tcs.com/Eu-Detect/.

Xenobiotic Metabolism and Gut Microbiomes

Humans are exposed to numerous xenobiotics, a majority of which are in the form of pharmaceuticals. Apart from human enzymes, recent studies have indicated the role of the gut bacterial community (microbiome) in metabolizing xenobiotics. However, little is known about the contribution of the plethora of gut microbiome in xenobiotic metabolism. The present study reports the results of analyses on xenobiotic metabolizing enzymes in various human gut microbiomes. A total of 397 available gut metagenomes from individuals of varying age groups from 8 nationalities were analyzed. Based on the diversities and abundances of the xenobiotic metabolizing enzymes, various bacterial taxa were classified into three groups, namely, least versatile, intermediately versatile and highly versatile xenobiotic metabolizers. Most interestingly, specific relationships were observed between the overall drug consumption profile and the abundance and diversity of the xenobiotic metabolizing repertoire in various geographies. The obtained differential abundance patterns of xenobiotic metabolizing enzymes and bacterial genera harboring them, suggest their links to pharmacokinetic variations among individuals. Additional analyses of a few well studied classes of drug modifying enzymes (DMEs) also indicate geographic as well as age specific trends.