the thermal springs are high - temperature aquatic ecosystems and are broadly distributed throughout the world . microbial life in hot springs is representing a stable and broad scientific interest in the past decades . the diversity of unculturable bacteria demonstrates a vast gene pool for biotechnological exploration and creates a major face for microbiologists to know their phylogenetic correlation and ecological implication . besides that metagenomic investigation using high throughput sequencing or library construction shows an enormous consequence for demonstrating microbial composition and functionality in extreme environments ,   and for discovering novel genes , , . despite the possibility of the presence of novel microbes with high economic and industrial prospective , there are no detailed reports available on the microbial diversity of most of hot spring of indian subcontinent . therefore , the present research was designed to study community and functional microbial diversity with metagenomic approach . india is one of the most tectonically active areas in the world and has many hot springs . the lasundra hot spring ( n 2255'30 " and e 7312'33 " ) is one of the less explored sulfuric hot water spring in kheda district , gujarat state , india . the temperature of the hot spring reservoir water is reported ranging from 42  c to 52  c . a total discharge of water from the spring was about 18,000  l / h , which is used by visitors to take a bath for therapeutic belief . the maximum temperature of the central reservoir was 52  c and the nearby reservoirs have a temperature varying from 42  c to 49  c . the base of the reservoir contains stones and algal growth along with soil sediments . in the present investigation , the total community dna was isolated from water samples of the spring using a hipura water dna isolation kit ( mb547 ; hi - media , mumbai , india ) . shotgun metagenomic sequencing was performed using the 318 chip and 300-bp chemistry ion torrent pgm platform as per the manufacturer 's instructions . analysis and annotation of output data were carried out by mg - rast v 3.5 online server using default parameter . the output file contained a total of 606,867 sequences comprised of 98,567,305  bps with 46% g  +  c content . the community analysis revealed predominance of domain bacteria ( 99.21% ) , followed by eukaryota ( 0.43% ) and archaea ( 0.11% ) . a total of 33 prokaryotic phyla including 29 bacterial and 4 archaeal phyla along with 20 eukaryotic phyla were detected . abundant prokaryotic phyla were firmicutes ( 95.5% ) followed by proteobacteria ( 2.0% ) , actinobacteria ( 0.8% ) , bacteroidetes ( 0.1% ) , cyanobacteria ( 0.1% ) and euryarchaeota ( 0.09% ) . at the family level , bacillaceae ( 90.1% ) followed by paenibacillaceae ( 1.3% ) , clostridiaceae ( 0.8% ) , listeriaceae ( 0.5% ) and staphylococcaceae ( 0.5% ) were major dominating bacterial families . most leading genera were bacillus ( 86.7% ) , geobacillus ( 2.4% ) , paenibacillus ( 1.0% ) , clostridium ( 0.7% ) and listeria ( 0.5% ) ( fig .  1 ) . a total of 388,662 sequences have passed quality control and out of them 76.9% produced 183,408 predicted protein coding regions . from this 183,408 predicted protein features , 66.7% of features assigned an annotation using m5nr database whereas , 33.3% of features have no significant similarities to the protein database . a total of 104,110 features were assigned to functional categories with cog approach and interestingly 45.4% was metabolism connected and 19.6% falls in poorly characterized group that is higher than the earlier metagenome study of saline desert   indicating the possibilities of gaining novel gene from the metagenome . similarly , in subsystem based annotation approach 14.0% was carbohydrates , 7.0% was protein metabolism and 3.0% genes were belonged to stress responses i.e. oxidative stress , periplasmic stress , osmotic stress , heat shock , cold shock , acid stress , and detoxification ( fig .  2 ) . furthermore , the detection of genes encoded for secondary metabolism and metabolism of aromatic compounds indicates that the hot spring is prosperous with commercially valuable enzymes ( fig .  2 ) . lipase producing bacillus licheniformis and bacillus subtilis were reported from this hot spring earlier , which support our study . the detection of the photosynthetic bacteria , heterotrophs and autotrophs in the hot spring metagenome suggested the nutritive interaction among the microorganism . to the best of our knowledge , this is the first study that describes the complete profiling of microbial diversity from the lasundra hot spring using the next generation sequencing approach . metagenomic sequence analysis may provide significant clues in understanding the taxonomic structure and functional potential of lasundra hot spring with the possibility of getting novel genes and microbes for research and commercial application .
this is the first report on the metagenomic approach for unveiling the microbial diversity of lasundra hot spring , gujarat state , india . 
 high - throughput sequencing of community dna was performed on an ion torrent pgm platform . 
 metagenome consisted of 606,867 sequences represent 98,567,305  bps size with an average length of 162  bps and 46% g  +  c content . 
 metagenome sequence information is available at ebi under ebi metagenomic database with accession no . 
 erp009313 . 
 mg - rast assisted community analysis revealed that 99.21% sequences were bacterial origin , 0.43% was fit to eukaryotes and 0.11% belongs to archaea . 
 a total of 29 bacterial , 20 eukaryotic and 4 archaeal phyla were detected . 
 abundant genera were bacillus ( 86.7% ) , geobacillus ( 2.4% ) , paenibacillus ( 1.0% ) , clostridium ( 0.7% ) and listeria ( 0.5% ) , that represent 91.52% in metagenome . 
 in functional analysis , cluster of orthologous group ( cog ) based annotation revealed that 45.4% was metabolism connected and 19.6% falls in poorly characterized group . 
 subsystem based annotation approach suggests that the 14.0% was carbohydrates , 7.0% was protein metabolism and 3.0% genes for various stress responses together with the versatile presence of commercially useful traits .