Skip to main content
  • Software review
  • Open access
  • Published:

MALINA: a web service for visual analytics of human gut microbiota whole-genome metagenomic reads

Abstract

MALINA is a web service for bioinformatic analysis of whole-genome metagenomic data obtained from human gut microbiota sequencing. As input data, it accepts metagenomic reads of various sequencing technologies, including long reads (such as Sanger and 454 sequencing) and next-generation (including SOLiD and Illumina). It is the first metagenomic web service that is capable of processing SOLiD color-space reads, to authors’ knowledge. The web service allows phylogenetic and functional profiling of metagenomic samples using coverage depth resulting from the alignment of the reads to the catalogue of reference sequences which are built into the pipeline and contain prevalent microbial genomes and genes of human gut microbiota. The obtained metagenomic composition vectors are processed by the statistical analysis and visualization module containing methods for clustering, dimension reduction and group comparison. Additionally, the MALINA database includes vectors of bacterial and functional composition for human gut microbiota samples from a large number of existing studies allowing their comparative analysis together with user samples, namely datasets from Russian Metagenome project, MetaHIT and Human Microbiome Project (downloaded fromhttp://hmpdacc.org). MALINA is made freely available on the web athttp://malina.metagenome.ru. The website is implemented in JavaScript (using Ext JS), Microsoft .NET Framework, MS SQL, Python, with all major browsers supported.

Background

Whole-genome sequencing of environmental samples is producing data at an increasing pace. With the advent of high-throughput next-generation sequencing (NGS) technologies, a deeper insight into phylogenetic and functional composition of metagenomes has become feasible. The research community has a need for robust data analysis tools that allow efficient description of composition, classification and clustering coupled with comprehensive visualization of results, while providing means for comparative analysis within the context of all accumulated metagenomic data for same type of environment. There is a number of existing web services (including CAMERA[1], IMG/M[2], MG-RAST[3], METAGENassist[4] and others) and stand-alone applications (including QIIME[5] and SmashCommunity[6]) that integrate data visualization and statistical analysis functionalities with databases of publicly available metagenomic data, allowing the user to compare his/her own samples with those of other researches. However, the number of pre-loaded human gut metagenomic samples in the repertoire of these tools is limited.

The human gut microbiome is one of the most extensively studied subjects in metagenomic research. It is of particular interest to scientists because of its significant role in host health status. Representative reference genomes for many taxa have been sequenced, and a catalogue of prevalent gut microbial genes has already been established[7]. MALINA exploits this accumulated knowledge in the form of reference sequence sets to provide a means for analyzing human gut whole-genome reads within the context of world public metagenomic datasets. The inclusion of a vast set of existing human gut metagenomic datasets allows the user to check which datasets are most similar to his/her own data, and, if present, to examine the metadata of those and pose hypotheses based on similarities. Features of MALINA and existing software allowing human gut whole-genome metagenomic reads analysis are compared in Table 1.

Table 1 Comparison of MALINA and existing software allowing analysis of human gut metagenomic reads

Implementation

The MALINA workflow is shown in Figure1. As input, MALINA accepts short nucleotide reads of length starting with 35 bp. Color-space (SOLiD), as well as long (such as Sanger, 454) reads are supported. To our knowledge, MALINA is the first metagenomic analysis web-service supporting SOLiD color-space reads. It is beneficial, considering the increasing volume of metagenomic data sequenced using this technology. Files with reads are uploaded by FTP. Through the web interface, the user creates groups of samples, with each sample including one or more read sets. The files for a given sample are associated with appropriate read sets and prepared for analysis.

Figure 1
figure 1

MALINA workflow. Input data and the main stages of analysis are illustrated.

The MALINA analysis pipeline characterizes metagenomic composition in two ways: phylogenetic and functional - by assessing relative abundance of microbial genera and genes, correspondingly. In case of genes, total metabolic potential of all microbes is described. The quantitative profiling is based on alignment of reads to reference genomes and gene catalogue. The genome catalogue contains more than 440 genomes of human gastrointestinal bacteria obtained from HMP, NCBI and relevant studies of human gut microbiome. The gene catalogue of prevalent human gut microbial genes discovered by MetaHIT project consists of 3.3 million genes. After the reads are aligned to reference set, the resulting position-wise coverage of each sequence is normalized by its length and total number of reads in read set. Summed over genera (for genomes) or functional groups (clusters of orthologous groups, COGs[8]) for genes, it yields relative abundance of phylogenetic and functional units. For functional profiling, COG annotation from MetaHIT gene catalogue is used. Each metagenomic read set is thus described by two feature vectors.

Feature vectors of the read sets selected by the user are subject to statistical analysis and visualization: boxplots of the most abundant genera/COGs, principal components analysis (PCA), clustering (partitioning around medoids [PAM] and hierarchical clustering)[9], multidimensional scaling (MDS)[10] and between-class analysis (BCA)[11] for the results of clustering. PCA plot shows 2D projection of feature vectors along the directions of maximum variance in the data, with arrows showing genera “drivers” that contribute most strongly to variation between samples. For COG groups, PCA plots are constructed separately for several functional classes: antibiotic resistance (COGs were collected from ARDB database[12]), transcription factors (COGs selected from total COG list according to description, i.e. “transcription regulator/factor/repressor”) and vitamin metabolism (COGs from KEGG[13] vitamin synthesis pathways). PAM clustering calculates the optimal number of clusters and assigns the samples to the clusters. BCA is a special case of PCA with respect to an instrumental variable (that is represented by cluster number here) producing plot that highlights differences between the clusters. The second implemented clustering algorithm, hierarchical clustering, produces dendrogram heatmap of abundance. Sample visualizations produced by MALINA are shown in Figure2. Moreover, statistical analysis includes detection of genera and gene categories discriminatory among the clusters by the Mann–Whitney test[14] and Random Forests algorithm[15] as well as taxa co-occurrence analysis based on the abundance values correlation.

Figure 2
figure 2

Examples of visualizations by MALINA. Metagenomic samples are denoted by text tags. (A) BCA plot of phylogenetic abundance, with the clusters visualized as ellipses and the samples connected with the centers of their respective clusters by lines. (B) PCA of COGs associated with vitamin metabolism, with drivers shown by red arrows. (C) Dendrogram heatmap of phylogenetic abundance based on hierarchical clustering algorithm (with samples as columns and genera as rows).

The plots can be downloaded as PDF files, and the relative abundance, clustering results and other output tables can be downloaded as tabulated text files. All features of MALINA are available without registration, via guest account. The user can register a dedicated account for free to provide privacy of uploaded data and results, as well as “analysis complete” notifications by e-mail.

An important functionality of MALINA is that besides the user’s own data, it is possible to co-analyze it with cohorts from large existing human gut metagenomic studies: 85 Illumina samples and 37 Sanger samples from MetaHIT study[7], as well as 139 Illumina samples from HMPDACC and 96 SOLiD samples from a new, previously unpublished Russian metagenomic study. Thus clustering functionality is of particular interest to researchers exploring human gut microbiota in relation to the concept of enterotypes[16] across a large number of samples. A stand-alone analysis of pre-loaded datasets is also available.

The user interface is implemented using Ext JS framework. Read alignment is performed using Bowtie[17]. In the interest of performance, MALINA does not filter reads using raw quality score, as the experience showed that filtration does not significantly increase the fraction of mapped reads. However, such preprocessing can be performed by the user manually. Coverage statistics are calculated using BEDtools[18]. Statistical analysis is implemented in R[19] using ade4[11], cluster[20], ecodist[21], fpc[22] and randomForest[23] packages. The pipeline steps are integrated using Oracle database, Microsoft .NET framework and Python.

Conclusions

MALINA allows an easy and intuitive way to infer metagenomic composition from reads and to analyze similarity of samples and organization into clusters within the global context of human gut metagenomic datasets. The features include statistical analysis methods like clustering and group comparison, as well as illustrative visualizations of phylogenetic and functional composition. The support for color-space SOLiD reads is a unique feature that makes MALINA a particularly valuable service to the growing community of researchers using SOLiD technology for metagenomic analysis.

The reference gene catalogue used in MALINA will be updated regularly as new version of MetaHIT data becomes available. In the future, it is planned that additional detailed metadata will be associated with the samples, allowing the user to check if newly sequenced samples are similar to certain groups distinguished by medical, ethno-geographic or dietary factors. The further development of the web service will include updates of the human gut samples database from Russian population as well as from other new studies. Support for diverse types of environment profiles besides human gut and additional methods for statistical analysis and visualization will be added.

Availability and requirements

Project name: MALINA

Project page:http://malina.metagenome.ru

Operating system: platform independent web site

Programming languages: Microsoft C# .NET, JavaScript, R, Python

Other requirements: None

License: FreeBSD

Any restrictions to use by non-academics: None

Funding

This work was supported by State Contracts 16.512.11.2111, 16.552.11.7034 and RFBR Grant 12-07-90008.

Abbreviations

BCA:

Between-Class Analysis

HMPDACC:

Human Microbiome Project Data Acquisition and Coordination Center

MDS:

Multidimensional scaling

PCA:

Principal Components Analysis.

References

  1. Seshadri R, Kravitz SA, Smarr L, Gilna P, Frazier M: CAMERA: a community resource for metagenomics. PLoS Biol. 2007, 5: e75-10.1371/journal.pbio.0050075.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Markowitz VM, Ivanova NN, Szeto E, Palaniappan K, Chu K, Dalevi D, Chen IA, Grechkin Y, Dubchak I, Anderson I, Lykidis A, Mavromatis K, Hugenholtz P, Kyrpides NC: IMG/M: a data management and analysis system for metagenomes. Nucleic Acids Res. 2008, 36: D534-D538.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, Wilkening J, Edwards RA: The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinforma. 2008, 9: 386-10.1186/1471-2105-9-386.

    Article  CAS  Google Scholar 

  4. Arndt D, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, Sinelnikov I, Budwill K, Nesbo CL, Wishart DS: METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012, 40: 88-95. 10.1093/nar/gkr734.

    Article  Google Scholar 

  5. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, Mcdonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Arumugam M, Harrington ED, Foerstner KU, Raes J, Bork P: SmashCommunity: a metagenomic annotation and analysis tool. Bioinformatics. 2010, 26: 2977-2978. 10.1093/bioinformatics/btq536.

    Article  CAS  PubMed  Google Scholar 

  7. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto JM, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464: 59-65. 10.1038/nature08821.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science. 1997, 278: 631-637. 10.1126/science.278.5338.631.

    Article  CAS  PubMed  Google Scholar 

  9. Kaufman L, Rousseeuw PJ: Finding groups in data: An introduction to Cluster Analysis. 1990, New York: Wiley

    Book  Google Scholar 

  10. Torgerson WS: Theory & methods of scaling. 1958, New York: Wiley

    Google Scholar 

  11. Dray S, Dufour AB: The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007, 22 (4): 1-20.

    Article  Google Scholar 

  12. ARDB:http://ardb.cbcb.umd.edu/.

  13. KEGG: Kyoto Encyclopedia of Genes and Genomes.http://www.genome.jp/kegg.

  14. Frank W: Individual comparisons by ranking methods. Biom Bull. 1945, 1 (6): 80-83. 10.2307/3001968.

    Article  Google Scholar 

  15. Breiman L: Random forests. Mach Learn. 2001, 45 (1): 5-32. 10.1023/A:1010933404324.

    Article  Google Scholar 

  16. Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S: Enterotypes of the human gut microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Quinlan A, Hall I: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842. 10.1093/bioinformatics/btq033.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. R Development Core Team: R: A Language and Environment for Statistical Computing. 2010, Vienna, Austria: R Foundation for Statistical Computing

    Google Scholar 

  20. Maechler M, Rousseeuw P, Struyf A, Hubert M, Hornik K: Cluster: cluster analysis basics and extensions.http://cran.r-project.org/web/packages/cluster/index.html.

  21. Goslee SC, Urban DL: The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw. 2007, 22 (7): 1-19.

    Article  Google Scholar 

  22. Hennig C: fpc: Flexible procedures for clustering.http://cran.r-project.org/web/packages/fpc/index.html.

  23. Liaw A, Wiener M: Classification and regression by randomForest. R News. 2002, 2 (3): 18-22.

    Google Scholar 

Download references

Acknowledgements

We thank Prof. Vadim Govorun who participated in study design and coordination. Russian metagenomic samples were sequenced by Genomic Center of Research Institute of Physico-Chemical Medicine.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander V Tyakht.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ESK, OVS, AKL and IYK extracted metagenomic DNA, prepared the libraries for sequencing and performed sequencing on SOLiD 4, yielding readsets used as part of web-service. DGA, AVT and MSB designed the pipeline, performed coding and processed sequenced data. IAA designed the database and developed Web-interface. ASP maintained the database, developed statistical analysis and visualization module. AVT, AVP and ASP wrote the manuscript text. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Tyakht, A.V., Popenko, A.S., Belenikin, M.S. et al. MALINA: a web service for visual analytics of human gut microbiota whole-genome metagenomic reads. Source Code Biol Med 7, 13 (2012). https://doi.org/10.1186/1751-0473-7-13

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1751-0473-7-13

Keywords