microbial samples from extreme habitats ( extremophiles ) represent a huge reservoir of genetic diversity and a potent source for industrially important enzymes   but the culturability of microorganisms from extreme habitats is generally very poor . for long time
, full length sequences generated by sanger sequencing of 16s rrna clone libraries were considered to be the gold standard for phylogenetic analysis but soon it was realized that this procedure is not only expensive but also have many limitations . to overcome the limitations and to gain high taxonomic resolution of bacterial population in extreme habitats , culture - independent 16s
both 454 ( roche ) and illumina platforms are now largely used to study metagenome and microbial diversity ,   but the most recent and advanced illumina miseq / hiseq sequencing platform provides a very distinguish and high quality view of microbial composition than other sequencing technologies .
hot springs are unique sites for extremophilic microorganisms and are of great interest for many years because enzymes obtained from them have been proved to be extremely valuable as biocatalysts for industrial and biotechnological purposes .
moreover , many unknown microbial species and genes have been revealed in culture - independent microbial diversity assessment of hot springs .
the peruvian amazon rainforest is one of the most biologically diverse areas on earth and a rich source of several novel microbial species .
peruvian amazon is endowed with few hot springs but till date none of them was explored in detail to generate high resolution microbial profile .
so , the main goal of this study is to generate a high resolution microbial phylogenetic profile of an amazonian hot spring . in this study water samples
were collected from four random sampling points of ac and mixed in an equal ratio as previously stated by chan et al .   during metagenomic study of a malaysian hot spring .
metagenomic dna was extracted using powerwater dna isolation kits ( mo bio laboratories ) following manufacturer 's protocol .
twenty five nanograms of nanodrop quantified dna was used for amplifying the v3v4 region of 16s rrna with specific primers which also have a 
tag sequence that are complementary to illumina sequence adapter and index primers from the nextera xt index kit v2 .
this round of pcr generates single amplicons of ~  530  bp . in the next round of pcr ( indexing pcr )
illumina sequencing adapters and dual indexing barcodes are added using limited cycle pcr to give a final product of ~  610  bp .
the libraries were cleaned using highprep pcr ( magbio , cat # ac-60050 ) magnetic beads , qubit quantified and validated for quality by running an aliquot on high sensitivity bioanalyzer chip ( agilent ) .
finally , the cleaned libraries were sequenced in illumina hiseq platform at genotypic technology private limited , bangalore , india .
qiime pipelines   was used for 16s rna detection , clustering and otu picking followed by biom file generation and statistical analysis .
sra files were deposited to ncbi database under accession number srx1809286 . in the present study functional analysis of 16s amplicons
after quality filtration and adapter trimming of raw reads , clean sequences were clustered into 5434 operational taxonomic units ( otus ) using a 97% similarity cut off .
good 's coverage estimator revealed that >  99% of the species were estimated , while high values of chao1 richness estimator ( 9762 ) and shannon diversity index ( 4.16 ) indicated that microbial communities in ac are highly rich and diverse .
all the resulting fragments were then classified into 58 phyla , 165 classes , 300 orders , 520 families and 954 genera .
the top 5 represented phyla were proteobacteria ( 88.06% ) , thermi ( 6.43% ) , firmicutes ( 3.41% ) , aquificae ( 1.10% ) and chloroflexi ( 0.41% ) ( fig .
1b ) ; however , an unidentified bacterial phylum ( 0.10% ) was also found among top 10 bacterial phyla .
an otu based phylogenetic tree displayed the genetic diversity among ac microbial community ( fig .
gammaproteobacteria ( 86.1% ) was found to be the highest represented class in ac microbial community followed by deinococci ( 6.43% ) , bacilli ( 3.22% ) , betaproteobacteria ( 1.28% ) and aquificae ( 1.10% ) .
a krona chart was constructed to illustrate the distribution pattern of phyla proteobacteria in ac hot spring ( fig .
further affiliation revealed that in ac microbial community the most abundant genus was acinetobacter ( 71.09% ) of moraxellaceae family in which most of the species were unidentified ( 48.29% ) ; remaining major genera included pseudomonas ( 8.59% ) of pseudomonadaceae family , thermus ( 5.99% ) of thermaceae family , enhydrobacter ( 3.38% ) of moraxellaceae family and brevibacillus ( 1.52% ) of paenibacillaceae family .
the dominance of proteobacteria , firmicutes and choloroflexi phyla in hot spring was reported earlier in few studies .
however , in contrast with our data , 16s rrna based microbial diversity analysis from little hot creek hot springs ( temperature 78.782.5  c and ph  6.756.97 ) , california , similar to ac hot spring showed the dominance of the phyla thermodesulfobacteria , deinococcus 
aquificae and thermotogae are two best known hyperermophilic bacterial phyla that show high heat tolerance . in this
study an unknown genus ( 1.09% ) of aquificae dominates the sample , the other two low abundant genera detected in this study were hydrogenobacter ( 0.01% ) and hydrogenivirga ( 0.00% ) ; while in thermotogae the genus fervidobacterium ( 0.09% ) was found to be the most represented one .
undoubtedly , microbial diversity in hot springs is greatly affected by the ph of the water .
water samples from an acidic hot spring ( ph  3.54 ) and a circumneutral hot spring ( ph  7.27.4 ) were analyzed previously   and the results revealed that thermotogae and gammaproteobacteria dominated the acidic hot spring , while thermodesulfobacteria , gammaproteobacteria and betaproteobacteria monopolized the circumneutral hot spring . in alkaline hot springs ,
thermus 
 , hydrogenobacter 
 , caldicellulosiruptor , dictyoglomus and fervidobacterium 
  were reported to be the most represented genera .
ac hotspring is slightly acidic to circumneutral in nature ( ph  5.07.0 ) and in agreement with previous reports monopolization of gammaproteobacteria was also found in this study .
regarding archaeal community structure euryarchaeota and crenarchaeota were the two major represented archaeal phyla detected in ac but with low read count . in a recent study thaumarchaeota and crenarchaeota
were found in both acidic and alkaline thermal spring environments , while euryarchaeota was only found in the acidic environment .
further affiliation of archaeal community in this study revealed that different unidentified families of the order pgrfc26 under crenarchaeota and few methanogenic members of the family methanosaetaceae , methanobacteriaceae and methanosarcinaceae under euryarchaeota were also present in ac hot spring . among which the complete genome sequence of a strain of methanobacterium thermoautotrophicum , one of the best methane producing thermophilic archaeon are available in the database . in this study ,
picrust , a modern tool designed to infer metagenomic information from 16s amplicon sequencing data , was used .
picrust uses the otu table of assigned taxa and their relative distribution to generate the relative abundance of functional categories based on sequenced genomes .
( 12.2% ) and  unknown(8.9% ) function categories , the highest represented category at second tier was  amino acid transport and metabolism 
( 8.5% ) , followed by  transcription  ( 7.1% ) and  energy production and conversion 
 genetic information processing  ( 16.2% ) ,  unclassified  ( 15.8% ) and  environmental information processing  ( 12.6% ) ( fig .
genetic information processing  was also reported previously in a malayasian hot spring microbial community .
 aguas calientes ( ac ) is an isolated geothermal spring located deep into the amazon rainforest ( 72112 s , 750054 w ) of peru . 
 this geothermal spring is slightly acidic ( ph  5.07.0 ) in nature , with temperatures varying from 45 to 90  c and continually fed by plant litter , resulting in a relatively high degree of total organic content ( toc ) . 
 pooled water sample was analyzed at 16s rrna v3v4 hypervariable region by amplicon metagenome sequencing on illumina hiseq platform . 
 a total of 2,976,534 paired ends reads were generated which were assigned into 5434 numbers of otus . 
 all the resulting 16s rrna fragments were then classified into 58 bacterial phyla and 2 archaeal phyla . 
 proteobacteria ( 88.06% ) was found to be the highest represented phyla followed by thermi ( 6.43% ) , firmicutes ( 3.41% ) and aquificae ( 1.10% ) , respectively . 
 crenarchaeota and euryarchaeota were the only 2 archaeal phyla detected in this study with low abundance . 
 metagenomic sequences were deposited to sra database which is available at ncbi with accession number srx1809286 . 
 functional categorization of the assigned otus was performed using picrust tool . in cog analysis 
  amino acid transport and metabolism  
 ( 8.5% ) was found to be the highest represented category whereas among predicted kegg pathways  metabolism  ( 50.6% ) was the most abundant . 
 this is the first report of a high resolution microbial phylogenetic profile of an amazonian hot spring . 
Direct link to deposited data
Experimental design, materials and methods
Nucleotide sequence accession number
Competing interest