- Research article
- Open Access
Dominant bacterial phyla in caves and their predicted functional roles in C and N cycle
BMC Microbiology volume 17, Article number: 90 (2017)
Bacteria present in cave often survive by modifying their metabolic pathway or other mechanism. Understanding these adopted bacteria and their survival strategy inside the cave is an important aspect of microbial ecology. Present study focuses on the bacterial community and geochemistry in five caves of Mizoram, Northeast India. The objective of this study was to explore the taxonomic composition and presumed functional diversity of cave sediment metagenomes using paired end Illumina sequencing using V3 region of 16S rRNA gene and bioinformatics pipeline.
Actinobacteria, Proteobacteria, Verrucomicrobia and Acidobacteria were the major phyla in all the five cave sediment samples. Among the five caves the highest diversity is found in Lamsialpuk with a Shannon index 12.5 and the lowest in Bukpuk (Shannon index 8.22). In addition, imputed metagenomic approach was used to predict the functional role of microbial community in biogeochemical cycling in the cave environments. Functional module showed high representation of genes involved in Amino Acid Metabolism in (20.9%) and Carbohydrate Metabolism (20.4%) in the KEGG pathways. Genes responsible for carbon degradation, carbon fixation, methane metabolism, nitrification, nitrate reduction and ammonia assimilation were also predicted in the present study.
The cave sediments of the biodiversity hotspot region possessing a oligotrophic environment harbours high phylogenetic diversity dominated by Actinobacteria and Proteobacteria. Among the geochemical factors, ferric oxide was correlated with increased microbial diversity. In-silico analysis detected genes involved in carbon, nitrogen, methane metabolism and complex metabolic pathways responsible for the survival of the bacterial community in nutrient limited cave environments. Present study with Paired end Illumina sequencing along with bioinformatics analysis revealed the essential ecological role of the cave bacterial communities. These results will be useful in documenting the biospeleology of this region and systematic understanding of bacterial communities in natural sediment environments as well.
Bacteria constitute the major portion of the cave biodiversity and plays a key role in maintaining cave ecosystem . Limited nutrient and energy sources create an oligotrophic environment inside the caves, wherein the primary production is carried out by autotrophic bacteria which inturn supports the growth of several chemo-organotrophic microbes . Bacteria present under this oligotrophic environment often survive by modifying their metabolic pathway or other mechanism . Understanding these adopted bacteria and their survival strategy inside the oligotrophic environment is an important aspect of microbial ecology.
Geomicrobial Investigations in nutrient limited caves are sparse and most of them have been carried out using culture based techniques. Such approach can only detect a minute portion of the total community. Such limitation is solved by the introduction of next generation sequencing (NGS) and expands our knowledge on uncultured microbes . Although the cost of amplicon sequencing (16S rDNA) used for the bacterial community composition studies has rapidly decreased, the functional study using the Shotgun approach or Geochip still remains expensive and thus, is restricted for selected studies . An indirect approach is to compare the uncultured bacterial sequences with closely related and well studied microbes to predict the functional role in the ecosystem. This is also useful to understand the unknown energy source required for metabolism [6, 7]. A computational approach, PICRUSt (phylogenetic investigation of communities by reconstruction of unobserved states) based on the relationship between phylogeny and function was developed to predict functional diversity using 16S rDNA data and a reference database and has been used to study in diverse environments .
Cave microorganisms contain a wide range of bacterial groups influenced by the geology, soil or sediment and other factors . Geochemistry parameter often drives the diversity and bacterial community composition inside the caves . Present study focuses on the bacterial community and geochemistry in five unexplored and unknown caves of Mizoram, Northeast India falling under the less- known biodiversity hotspot zone of the eastern Himalayan belt. The objective of this study was to explore the taxonomic composition and to understand how the bacterial communities respond to the cave oligotropic environments. This study was based on the hypothesis that the undisturbed and nutrient- limited cave habitats will host specific bacterial species and the cave geochemical parameters might favour species diversity and richness.
Sample collection and community DNA extraction
Cave sediment samples were collected from different sites of the caves – Bukpuk (CBP V3), Lamsialpuk (CLP V3) and Reiekpuk (CRP V3) followed by sieved and preserved at 4 °C (Fig. 1). The geochemical and molecular data of the sediment sample Lamsialpuk (CLP V3) and Khuangcherapuk (CKP V3) were collected from our previous study [11, 12]. All sites were not subjected to any human disturbances, except CLPV3 . The elevation, pH and other geochemical parameters of the caves are given in Table 1. The pH of the sediment samples was analysed using pH meter (Eutech, pH 510, USA). Major oxides and trace elements were measured using X-ray Fluorescence (XRF) (Bruker AXS, S4 Pioneer, Germany) at IIT Rookie, India.
DNA was extracted from the cave sediment samples using the Fast DNA spin kit (MP Biomedical, Solon, OH, USA) and the V3 hypervariable region of the 16S rRNA gene was amplified using 10 pmol/μl of each forward 341F (5′-CCTACGGGAGGCAGCAG-3′) and reverse 518R (5′-ATTACCGCGGCTGCTGG-3′) primer. PCR Master Mix will contain 2 μL each primers, 0.5 μL of 40 mM dNTP (NEB, USA), 5 μL of 5X Phusion HF reaction buffer (NEB, USA), 0.2 μL of 2 U/μL F-540Special Phusion HS DNA Polymerase (NEB, USA), 5 ng input DNA and water to make up the total volume to 25 μL. The PCR conditions were 98 °C for 30 s followed by 30 cycles of 98 °C for 10 s; 72 °C for 30 s and a final extension at 72 °C for 5 s followed by 4 °C hold.
Pre-processing and sequence analysis
Paired end Illumina sequencing (2 × 150 bp) was carried out at Scigenome Lab, Cochin, India. Raw sequence data for the two cave sediment samples, Farpuk (CFPV3) and Khuangcherapuk (CKPV3), were derived from our previous study [11, 12]. Raw fastq sequences were processed using the QIIME software package v.1.8.0 [13, 14]. Poor quality (quality score < 25) and smaller reads (read length < 100 bp) were filtered out using the split_libraries command. Pre-processed sequence reads were clustered to operational taxonomic units (OTU’s) using UCLUST method with similarity threshold of 97%  and were taxonomically classified using Greengenes database. Relative abundance of the bacterial phyla was calculated using QIIME. Statistical analysis was performed after rarefying the OTU table to 50,000 sequences per sample. Alpha and beta diversity plots were also generated using QIIME. Beta diversity between five bacterial cave communities was measured using unweighted UniFrac approach . Pearson correlations between soil characteristics and bacterial major phylum were estimated using PASW Statistics 18 (SPSS Inc., Chicago, IL, USA). Additionally, we performed imputed metagenomic analysis by the genome prediction software PICRUSt . The input used here was normalized OTU prepared by closed reference based approach. OTU’s were assigned at 97% similarity and were mapped to the Greengenes ver.13.5 database for functional prediction.
Multivariate principal component analysis (PCA) of 20 physicochemical parameters i.e., pH, Na2O, MgO, Al2O3, SiO2, P2O5, SO3, K2O, CaO, Fe2O3, Cr2O3, MnO, NiO, CuO, ZnO, Rb2O, SrO, ZrO2, BaO and Cl was carried out to determine which environmental variables best explained the observed community patterns using the PAST v3.02 software .
Geochemical characteristics of the cave sediment samples
The pH of the five cave sediment were recorded in the range of 6.7–7.5. The highest pH was recorded at CLPV3 (7.5) followed by CFPV3 (7.3) and CBPV3 (7.2), whereas the lowest pH was recorded at CKPV3 (6.7). The concentration of the oxides such as Na2O, MgO, Al2O3, SiO2, P2O5, SO3, K2O, CaO, CuO, ZnO, Fe2O3, Cr2O3, MnO, NiO, CuO, ZnO, Rb2O, SrO, ZrO2, BaO and Cl varied among the samples (Table 1). Soil samples from both CKPV3 and CRPV3 had similar, but relatively lower pH compared to the other three cave samples. Similarly, CLPV3 and CFPV3 were also geochemically similar with high concentration of Na2O (Additional file 1: Figure S1). CBPV3 showed the highest concentration of P2O5, SO3, CaO, MnO, CuO, ZnO and SrO, whereas the lowest concentration of Al2O3, SiO2, Na2O and Fe2O3 compared to other caves. Interestingly, the elevation of CBPV3 was lower than the other four caves under study. A principal component analysis (PCA) of the physicochemical parameters showed that the five caves were separated into four geochemically distinct habitats. The first two principal components explained 88.06% of the total variance. The sample CKPV3 and CRPV3 were found geochemically similar and were grouped together in the 2-dimensional PCA plot. The key influencing parameters for the geochemical diversity were Na2O and P2O5, while Cl and SO3 were the other influencing parameters in component 1 and component 2, respectively.
Analysis of bacterial community composition
The high throughput sequencing effort yielded a total of 54,90,239 paired end reads with an average of 9,15,040 paired end reads per sample. After assembly and quality assessment of the reads, a total of 54,88,530 high quality reads were obtained. A total of 48 phyla (AC1, Acidobacteria, Actinobacteria, AD3, Armatimonadetes, Bacteroidetes, BHI80–139, BRC1, Caldithrix, Chlorobi, Chloroflexi, Cyanobacteria, Deferribacteres, Elusimicrobia, FCPU426, Fibrobacteres, Firmicutes, Fusobacteria, GAL15, Gemmatimonadetes, GN02, GN04, MVP-21, NC10, Nitrospirae, NKB19, OD1, OP1, OP11, OP3, OP8, OP9, Planctomycetes, Proteobacteria, SBR1093, SC4, Spirochaetes, SR1, Synergistetes, Tenericutes, Thermi, TM6, TM7, Verrucomicrobia, WPS-2, WS2, WS3 and ZB3) were detected from different cave sediments (Fig. 2). The total bacterial community analysis showed that the phylum Actinobacteria was the most dominant contributing up to 65.1%, followed by Proteobacteria (24.8%), Acidobacteria (4.2%) and Firmicutes (3.6%) and the top ten phyla present in individual cave is shown in Fig. 3.
In the present study, the identified class under this phylum were Actinobacteria, Acidimicrobiia, Thermoleophilia, Rubrobacteria, MB-A2–108, Coriobacteriia, Nitriliruptoria,
OPB41 and KIST-JJY01. High abundance of dominant family (>0.01%) under Actinobacteria were Nocardiaceae, Streptomycetaceae, Micrococcaceae, Frankiaceae, Gaiellaceae, Pseudonocardiaceae, Streptomycetaceae, EB1017, Mycobacteriaceae, Actinosynnemataceae, Corynebacteriaceae, Rubrobacteraceae, Nocardioidaceae, Micromonosporacea, Geodermatophilaceae, Sporichthyaceae, Actinosynnemataceae, Nakamurellaceae, Pseudonocardiaceae, Cryptosporangiaceae, Kineosporiaceae and Ruaniaceae. Other dominant genus under Actinomycetes was Mycobacterium, Corynebacterium, Rubrobacter, Actinoplanes, Saccharothrix and Pseudonocardia.
Within the Proteobacteria, most phylotypes were classified under the class Alphaproteobacteria and Gammaproteobacteria. Other identified class were Betaproteobacteria, Deltaproteobacteria, Gammaproteobacteria, TA18, Epsilonproteobacteria and Zetaproteobacteria. Abundant genera (≥0.01%) under this phylum were Rhodoplanes, Kaistobacter, Sphingomonas, Bradyrhizobium, Alteromonas, Acidiphilium and Halomona. Under the class Alphaproteobacteria, two families (Hyphomicrobiaceae and Sphingomonadaceae) and three abundant genus (Rhodoplanes, Kaistobacter and Sphingomonas) were identified. Other detected genera, present in low abundance, under this class were Candidatus entotheonella, Plesiocystis, Desulfococcus, Nannocystis, Anaeromyxobacter, Sorangium, Haliangium, Geobacter, Cystobacter and Syntrophus. The dominant genera under the class Gammaproteobacteria were Alteromonas and Halomonas. Other genera (<0.01%) were Acinetobacter, Alcanivorax, Aquicella, Cronobacter, Dickeya, Dokdonella, Enhydrobacter, Enterobacter, Enterovibrio, Erwinia, Fulvimonas, Glaciecola, Hafnia, Halorhodospira, Idiomarina, Klebsiella, Legionella, Luteibacter, Luteimonas, Marinobacter, Marinobacterium, Oceanospirillum, Providencia, Psychrobacter, Pseudoxanthomonas, Pseudomonas, Pseudoalteromonas, Rheinheimera, Rhodanobacter, Salinisphaera, Salinivibrio, Serpens, Serratia, Shewanella, Stenotrophomonas, Steroidobacter and Thermomonas. The class Epsilonproteobacteria was present in low abundance which consisted three families (Helicobacteraceae, Campylobacteraceae and Helicobacteraceae) and two genus (Arcobacter and Sulfurimonas). However, no genus was identified under the class Zetaproteobacteria and TA18.
Acidobacteria was the third dominant phyla with eight families and 10 identified genera. Dominant families under this phylum were Solibacteraceae, Koribacteraceae and Acidobacteriaceae. Assigned genera under the family Acidobacteriaceae were Acidobacterium, Edaphobacter, Terriglobus, Acidicapsa and Acidopila.
Diversity estimates of the cave bacterial community
Based on the Shannon index, a high bacterial diversity was observed in CLPV3 (12.50) and low in CBPV3 (8.22) (Table 2). The principal coordinate analysis plot of the UniFrac distance matrix distinguish CBPV3 from rest of the samples suggesting the presence of different composition of the bacterial communities, whereas other four cave samples had similar community composition (Fig. 4).
Function prediction using PICRUSt
Analysis revealed five functional modules (i.e. metabolism, genetic information processing, environmental information processing, cellular process and organismal systems) where metabolism was the most represented, accounting for about 60% of the entire data set. A deeper analysis of the terms encompassed by the metabolism functional module showed high representation of Amino Acid Metabolism (24%), Carbohydrate Metabolism (23%), Energy Metabolism (10%), Lipid Metabolism (10%), Metabolism of Cofactors and Vitamins (8%), Metabolism of Terpenoids and Polyketides (6%), Nucleotide Metabolism (6%), Metabolism of Other Amino Acids (4%), Enzyme Families (3.0%), Glycan Biosynthesis and Metabolism (3%) and Biosynthesis of Other Secondary Metabolites (3%) (Additional file 1: Figure S2).
With the carbon metabolism, three reactions were involved (carbon degradation, carbon fixation, and methane metabolism). The identified carbon degradation enzymes included genes encoding alpha-amylase, glucoamylase, neopullulanase, and pullulanase (involved in starch degradation); alpha-glucosidase, endoglucanase; beta-glucanase, beta-glucosidase (involved in cellulose degradation); arabinofuranosidase, xylanase, and mannanase (involved in hemicellulose degradation); Chitinase; beta-hexosaminidase; alpha-mannosidase and beta-mannosidase acetyl-glucosaminidase, polygalacturonase (involved in chitin, pectin degradation) and other carbohydrate degradation enzymes. A rare fraction of the predicted metagenomes sequence was classified as 4-hydroxybutyryl-CoA dehydratase. The predicted carbohydrate degrading enzymes were shown in Additional file 1: Table S1.
Gene’s codes for the enzymes methenyltetrahydrofolate cyclohydrolase is also detected in our study. Predicted genes and enzymes show the prevalence of methane cycle in the caves (Additional file 1: Table S2). Analysis also revealed nitrogen cycling genes involved in nitrification, nitrate reduction and ammonia assimilation. Genes codes for the enzyme involved nitronate monooxygenase; nitrile hydratase; nitrate reductase; nitrilase; nitric oxide dioxygenase; nitric oxide reductase; nitric-oxide synthase; nitrite reductase; nitric-oxide reductase; nitrogenase; nitric nitrogen fixation protein; nitroreductase/dihydropteridine reductase; nitrous-oxide reductase; nitroreductase; nitrate reductase; nitrogenase; nitric oxide reductase; nitrogenise (Additional file 1: Table S3).
Association between bacterial communities with geochemical parameters
A correlation analysis was performed to study the association between the most abundant phyla identified (AD3, Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Gemmatimonadetes, Proteobacteria, TM7 and WPS-2) and the geochemical parameters. Analysis revealed that Al2O3 was positively correlated with Chloroflexi (r = 0.627, p = 0.060); and MnO was negatively correlated with Acidobacteria (r = −0.790, p = −0.060). No other relationship between geochemical parameters and the relative abundance of the major phyla was significant different among sampling sites. Within the candidate phyla, MgO was correlated with the relative abundance of the AD3 (r = 0.978, p = 0.001), TM7 (r = 0.974, p = −0.001); and WPS-2 (r = 0.938, p = −0.006) (Additional file 1: Table S4). Furthermore, the content of Fe2O3 showed highest positive correlation with the Shannon diversity index (r = 0.926, p = 0.001), followed by Al2O3, NiO and negative correlation with SO3 and MnO (Additional file 1: Table S5).
Speleological studies with NGS approaches are now becoming an important approach for analyzing the concealed microbial diversity in belowground ecosystems . Adaptation of the microorganism in cave ecosystem mostly involves interaction with the minerals, mobilizing inorganic phosphate, oxidizing methane and hydrogen, and deriving energy by hydrolyzing macromolecules derived from other cave microbial communities . High competition for resources in nutrient limited environment helps in natural selection leading to innovation and diversification of bacterial communities . Present study documents the bacterial community composition along with the geochemical analysis of the bacterial community from five different cave sediments in Mizoram, a state of northeast India, situated in Indo-Burma biodiversity hotspot zone.
Analysis of bacterial community composition
All the cave samples were dominated by the phylum Actinobacteria as seen by our previous study using V4 hypervariable region of 16SrRNA .The three most abundant bacterial phyla detected in this study were Actinobacteria a common cave inhabitant has been isolated in rock walls and bioflim of various caves . Isolation of rare and novel Actinobacteria from unexplored environment is an important area of research . Members of the dominant family Nocardiaceae have previously been reported in cave ecosystem, are oligotrophic, and can metabolize various substrates such as toluene, herbicides, naphthalene and PCBs [23–27]. The genus Streptomyces was the second highest genus falling under the family Streptomycetaceae. Members of this group can metabolize different compounds including alcohols, sugars, amino acids and aromatic compounds and capable of synthesizing clinically useful antibiotics .
Proteobacteria was dominated by Alphaproteobacterial species and Gammaproteobacteria. Some species under this subphylum can survive under extreme environment by using ABC (ATP-Binding Cassettes) and TRAP: (Tri-partite ATP-independent periplasmic transporters) mechanism . The genus Rhodoplanes under the subphylum Alphaproteobacteria accounts for 0.15% of the total bacterial community and possesses Photo- and chemo-organ heterotrophic growth . They can also produce hopanoids and carotenoids [31, 32]. Another identified genus under Alphaproteobacteria- Sphingomona (0.133%), a group commonly found in nutrient-limited subsurface environments can metabolize a large number of different aromatic compounds .
The most abundant genus under Gammaproteobacteria was Alteromonas, a gram negative heterotrophic bacteria capable of degrading aromatic carbon rings introduced through oil spill . Another dominant genus under this subphylum was Halomonas known to resist extreme conditions and also involve in sandstone formations . Among the Betaproteobacteria, the most abundant genera were Thiobacillus, Burkholderia and Delftia, but they were present in less number (<0.002%). Thiobacillus can obtain energy by oxidizingo sulfur and ferrous iron compounds . Most of the members under the genus Burkholderia were diazotrophs and degrades a variety of xenobiotic compounds .
The unique characteristics of the genus Bdellovibrio are that they can enter into the periplasmic space of other bacteria and feed on the biopolymers and thereby used as biocontrol purposes . The abundant genus under Acidobacteria was Candidatus Solibacter, an aerobic, chemoorganotrophic bacteria having a large number of anion: cation symporters which helps them to survive in nutrient limited condition . Other abundant genus Candidatus koribacter was primarily considered as heterotroph .
Metabolic prediction using PICRUSt
The cave environment is a diverse habitat harbouring organisms from all hierarchies starting from prokaryotes to higher eukaryotes . Phylogenetic analysis using 16S SSU-rDNAs data were also applied to assume the metabolic role of the identified bacterial in cave ecosystem by aligning the sequence information to its next nearest culturable representatives [41, 42]. More recently, PICRUSt software package was developed used to infer the potential functional role of bacterial communities in the cave sediment samples using 16S SSU-rDNAs data .
Microbial communities are well known key players of biogeochemical cycles and mainly contribute to the global biogeochemical cycling of carbon and nitrogen . Present study detected enzyme 4-hydroxybutyryl-CoA dehydratase involved in the CO2-fixation of Archea and fermentation in bacteria which supports the hypothesis that autotrophic archaea contribute to carbon assimilation in cave and other environments [43–45]. Analysis also detected methenyltetrahydrofolate cyclohydrolase which is involved in reverse methanogenesis prevalent in anaerobic methanotrophic archaea [46, 47]. The presence of genes encoding proteins for the phosphate recycling mechanism, such as phosphonate transpoters (PhnB, PhnG, PhnH, PhnI, PhnJ, and PhnM) in the cave samples suggest that they form carbonphosphorus lyase complex which is involved in methane production from methyl phosphonate .
Role of bacteria in nitrogen cycle have been well studied in soil and aquatic habitats, but information on cave sediment is limited. Some reports are available where microbes can accrue energy as well as nutrients in oligotrophic environments through nitrogen cycling processes. Most of the genes involved in nitrogen cycle were detected in the present study. Presence of the genes codes for hydroxylamine oxidase indicates the presence of a key ammonia oxidizing bacteria (AOB) . Presence of AOB and sulfur-oxidizing bacteria were also reported in chemolithotrophic Cave  and thus lithochemotrophy might be a survival strategy of the bacterial communities present in the cave sediments. Identified genus, Nitrospira and Nitrosospira were reported to perform autotrophic nitrification which is an indication of CO2-fixation-coupled ammonia oxidation process in the studied cave ecosystems .
Association between bacterial communities with geochemical parameters
Bacterial community structure is greatly influenced by the mineral substrates present in an environment . Present study observed the positive relationship between Fe2O3, Al2O3 and NiO with the Shannon diversity index. Fe (II) is produced on the subsurface under anoxic conditions by dissimilatory iron (III) reducing bacteria (DIRB) coupled with biotic/abiotic weathering of minerals. Reduced metals inside the cave serve as a source of electron donor for bacterial growth [52, 53]. Only certain organisms can survive in the presence of oligotrophic forces and a high-metal environment, and the natural selection favours adaptations in microbial communities to sustain in these environments.
Present study used Illumina sequencing to examine the taxonomical diversity of bacterial communities present in cave sediment samples, which were collected from Mizoram, an Indo-Burma Biodiversity Hotspot. These oligotrophic cave harbours a high phylogenetic diversity, including organisms from all hierarchies as well as a higher proportion of unclassified sequences indicating the possibility of novel species. The cave sediments were dominated by Actinobacteria and Proteobacteria. Fe2O3 content was correlated with increased microbial diversity in these cave environments. Bioinformatics analysis detected genes involved in various metabolic pathways which are essential for the survival of the community in nutrient limited cave environments. Further research by cultivating the uncultured communities or whole genome sequencing is needed to illustrate the actual survival strategies in the cave environments.
ammonia oxidizing bacteria
human microbiome project
next generation sequencing
principal component analysis
phylogenetic investigation of communities by reconstruction of unobserved states
tri-partite ATP-independent periplasmic transporters
Ortiz M, Neilson JW, Nelson WM, Legatzki A, Byrne A, Yu Y, Wing RA, Soderlund CA, Pryor BM, Pierson LS, Maier RM. Profiling bacterial diversity and taxonomic composition on Speleothem surfaces in Kartchner caverns, AZ. Microb Ecol. 2013;65:371–83.
Canaveras JC, Sanchez-Moral S, Soler V, Saiz-Jimenez C. Microorganisms and microbially induced fabrics in cave walls. Geomicrobiol J. 2001;18:223–40.
Ivanova V, Tomova I, Kamburov A, Tomova A, Vasileva-Tonkova E, Kambourova M. High phylogenetic diversity of bacteria in the area of prehistoric paintings in Magura Cave,Bulgaria. J Caves Karst Stud. 2013;75(3):218–28.
De Mandal S, Zothansanga PAK, Bisht SS, Kumar NS. MiSeq HV4 16S rRNA gene analysis of bacterial community composition among the cave sediments of Indo-Burma biodiversity hotspot. Environ Sci Pollut Res. 2016;23:12216–26.
Gilbert JA, Dupont CL. Microbial Metagenomics: beyond the genome. Annu Rev Mar Sci. 2011;3:347–71.
Barns SM, Fundyga RE, Jeffries MW, Pace NR. Remarkable archaeal diversity detected in a Yellowstone National Park hot spring environment. Proc Natl Acad Sci U S A. 1994;91(5):1609–13.
Hugenholtz P, Goebel BM, Pace NR. Impact of culture-independent studies on the emerging phylogenetic view of bacterial diversity. J Bacteriol. 1998;180:4765–74.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Thurber RL, Knight R, Beiko RG. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.
Adetutu EM, Thorpe K, Bourne S, Cao X, Shahsavari E, Kirby G, Ball AS. Phylogenetic diversity of fungal communities in areas accessible and not accessible to tourists in Naracoorte caves. Mycologia. 2011;103:959–68.
Barton HA, Taylor NM, Kreate MP, Springer AC, Oehrle SA, Bertog JL. The impact of host rock geochemistry on bacterial community structure in oligotrophic cave environments. Int J Speleol. 2007;36:93–104.
De Mandal S, Panda AK, Lalnunmawii E, Bisht SS, Kumar NS. Illumina-based analysis of bacterial community in Khuangcherapuk cave of Mizoram, Northeast India. Genomics data. 2015a;5:13–4.
De Mandal S, Sanga Z, Kumar NS. Metagenome sequencing reveals Rhodococcus dominance in Farpuk cave, Mizoram, India, an eastern Himalayan biodiversity hot spot region. Genome announcements. 2015b;3(3):e00610–5.
Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010a;26:266–7.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña 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. 2010b;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35.
Hammer O, Harper DAT, Ryan PD. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica. 2001;4:1-9.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB. Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005;437(7057):376–80.
Barton HA, Jurado V. What's up down there? Microbial diversity in caves. Microbe-American Society for Microbiology. 2007;2(3):132–8.
Schluter D. Ecological causes of adaptive radiation. Am Nat. 1996;148:S40.
Chelius MK, Moore JC. Molecular phylogenetic analysis of archaea and bacteria in wind cave, South Dakota. Geomicrobiol J. 2004;21(2):123–34.
Wu Y, Tan L, Liu W, Wang B, Wang J, Cai Y, Lin X. Profiling bacterial diversity in a limestone cave of the western loess plateau of China. Front Microbiol. 2015;6:244.
Vander Geize R, Dijkhuizen L. Harnessing the catabolic diversity of rhodococci for environmental and biotechnological applications. Curr Opin Microbiol. 2004;7(3):255–61.
Burkowski A. Corynebacteria – Genomics and Molecular Biology. Norfolk: Caister Academic Press; 2008.
McLeod MP, Warren RL, Hsiao WW, Araki N, Myhre M, Fernandes C, Miyazawa D, Wong W, Lillquist AL, Wang D, Dosanjh M. The complete genome of Rhodococcus sp. RHA1 provides insights into a catabolic powerhouse. PNAS. 2006;103(42):15582–7.
Groth I, Saiz-Jimenez C. Actinomycetes in hypogean environments. Geomicrobiol J. 1999;16:1–8.
Yoon JH, Cho YG, Lee ST, Suzuki KI, Nakase T, Park YH. Nocardioides nitrophenolicus sp. nov., a p-nitrophenol-degrading bacterium. Int J Systematic Bacteriol. 1999;49:675–80.
Madigan M, Martinko J. Brock biology of microorganisms. 11th ed; 2005. p. 149–52.
Kumbhare SV, Dhotre DP, Dhar SK, Jani K, Apte DA, Shouche YS, Sharma A. Insights into diversity and imputed metabolic potential of bacterial communities in the continental shelf of Agatti Island. PLoS One. 2015;10(6):e0129864.
Lakshmi KVNS, Sasikala C, Ramana CV. Rhodoplanes pokkaliisoli sp. nov., a phototrophic alphaproteobacterium isolated from a waterlogged brackish paddy soil. Int J Syst Evol Microbiol. 2009;59:2153–7.
Lodha TD, Srinivas A, Sasikala C, Ramana CV. Hopanoid inventory of Rhodoplanes spp. Arch Microbiol. 2015;197(6):861–7.
Takaichi S, Sasikala C, Ramana Ch V, Okamura K, Hiraishi A. Carotenoids in Rhodoplanes species: variation of compositions and substrate specificity of predicted carotenogenesis enzymes. Curr Microbiol. 2012;65:150–5.
Balkwill DL, Drake GR, Reeves RH, Fredrickson JK, White DC, Ringelberg DB, Chandler DP, Romine MF, Kennedy DW, Spadoni CM. Taxonomic study of aromatic-degrading bacteria from deep-terrestrial-subsurface sediments and description of Sphingomonas aromaticivorans sp. nov., Sphingomonas subterranea sp. nov., and Sphingomonas stygia sp. nov. Int J Syst Bacteriol. 1997;47(1):191–201.
Jin HM, Kim JM, Lee HJ, Madsen EL, Jeon CO. Alteromonas as a key agent of polycyclic aromatic hydrocarbon biodegradation in crude oil-contaminated coastal sediment. Environ Sci Technol. 2012 Jun 27;46(14):7731–40.
Dong Y, Kumar CG, Chia N, Kim PJ, Miller PA, Price ND, Cann IK, Flynn TM, Sanford RA, Krapac IG, Locke RA. Halomonas sulfidaeris dominated microbial community inhabits a 1.8 km deep subsurface Cambrian sandstone reservoir. Environ Microbiol. 2014;16(6):1695–708.
Frankel RB, Bazylinski DA. Biologically induced mineralization by bacteria. Rev Mineral Geochem. 2003;54(1):95–114.
Rusch A, Islam S, Savalia P, Amend JP. Burkholderia insulsa sp. nov., a facultatively chemolithotrophic bacterium isolated from an arsenic-rich shallow marine hydrothermal system. Int J Syst Evol Microbiol. 2015;65(1):189–94.
Yair S, Yaacov D, Susan K, Jurkevitch E. Small eats big: ecology and diversity of Bdellovibrio and like organisms, and their dynamics in predator-prey interactions. Agronomie. 2003;23:433–9.
Rowe OF, Sánchez España J, Hallberg KB, Johnson DB. Microbial communities and geochemical dynamics in an extremely acidic, metal rich stream at an abandoned sulfide mine (Huelva, Spain) underpinned by two functional primary production systems. Environ Microbiol. 9(7):1761–71.
Sait M, Davis KE, Janssen PH. Effect of pH on isolation and distribution of members of subdivision 1 of the phylum Acidobacteria occurring in soil. Appl Environ Microbiol. 2006;72(3):1852–7.
Pace NR. A molecular view of microbial diversity and the biosphere. Science. 1997;276:734–40.
Barton HA, Northup DE. Geomicrobiology in cave environments: past, current and future perspectives. J Cave Karst Stud. 2007;69(1):163–78.
Berg IA, Kockelkorn D, Buckel W, Fuchs G. A 3-hydroxypropionate/4-hydroxybutyrate autotrophic carbon dioxide assimilation pathway in archaea. Science. 2007;318:1782–6.
Zhang LM, Offre PR, He JZ, Verhamme DT, Nicol GW, Prosser JI. Autotrophic ammonia oxidation by soil thaumarchaea. Proc Natl Acad Sci U S A. 2010;107:17240–5.
Ortiz M, Legatzki A, Neilson JW, Fryslie B, Nelson WM, Wing RA, Soderlund CA, Pryor BM, Maier RM. Making a living while starving in the dark: metagenomic insights into the energy dynamics of a carbonate cave. ISME J. 2014;8:478–91.
Williams TJ, Cavicchioli R. Marine metaproteomics: deciphering the microbial metabolic food web. Trends Microbiol. 2014;22(5):248–60.
Wang DZ, Xie ZX, Zhang SF. Marine metaproteomics: current status and future directions. J Proteome. 2014;97:27–35.
Chen Y, Wu L, Boden R, Hillebrand A, Kumaresan D, Moussard H, et al. Life without light: microbial diversity and evidence of sulfur- and ammonium-based chemolithotrophy in Movile cave. ISME J. 2009;3:1093–104.
Kim BK, Jung MY, Yu DS, Park SJ, Oh TK, Rhee SK, Kim JF. Genome sequence of an ammonia-oxidizing soil archaeon, ‘Candidatus Nitrosoarchaeum koreensis’ MY1. J Bacteriol. 2011;193:5539–40.
Sarbu SM, Kane TC, Kinkle BK. A chemoautotrophically based cave ecosystem. Science. 1996;272:1953–5.
Carson JK, Campbell L, Rooney D, Clipson N, Gleeson DB. Minerals in soil select distinct bacterial communities in their microhabitats. FEMS Microbiol Ecol. 2009;67(3):381–8.
Cunningham KI, Northup DE, Pollastro RM, Wright WG, LaRock EJ. Bacteria, fungi and biokarst in Lechuguilla cave, Carlsbad caverns National Park, New Mexico. Environ Geol. 1995;25:2–8.
Northup DE, Dahm CN, Melim LA, Spilde MN, Crossey LJ, Lavoie KH, Mallory LM, Boston PJ, Cunningham KI, Barns SM. Evidence for geomicrobiological interactions in Guadalupe caves. J Cave Karst Stud. 2000;62:80–90.
This research was supported by a grant from the Bioinformatics Infrastructure Facility sponsored by Department of Biotechnology, Govt. of India, New Delhi. The grant has enabled us to establish computational facility to work on Metagenomics pipeline, and the Advanced State Biotech Hub grant has helped in sampling and generate NGS data. The DeLCON facility has enabled for study design, analysis and interpretation of data.
Availability of data and materials
All datasets are presented in the main text and the Additional file 1. Paired end Illumina sequence data from this study were submitted to the NCBI and the BioProject accession numbers are PRJNA335996; PRJNA335990 and PRJNA335989.
SDM performed the sample collection and drafted the manuscript. SDM and RC carried out bioinformatics and statistical analyses. NSK and RC supervised the Illumina sequencing and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
No specific permissions were required for this study, as no human subjects were involved and also did not include any endangered or protected species.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of the genes codes for enzymes involved in carbohydrate degradation identified using PICRUSt. Table S2. List of the homologs of methanogenesis-associated genes that were identified from the five cave sediments using PICRUSt. Table S3. List of the genes coding for enzymes involved in nitrogen cycle identified using PICRUSt. Table S4. Pearson correlation (PC) between physiochemical factors with the dominant bacterial phyla. Table S5. Pearson correlation (PC) between physiochemical factors with the bacterial diversity. Figure S1. Bioplot generated for the Principal Component Analysis (PCA) of 20 geochemical variables. Cave samples are shown as colored symbols and physicochemical variables are represented by green lines. Figure S2. Relative abundance of the functional genes present in the cave samples. (DOCX 122 kb)