BioMed Research International

BioMed Research International / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 1869748 | https://doi.org/10.1155/2021/1869748

Punam Yadav, Shikha Sharma, Tribikram Bhattarai, Lakshmaiah Sreerama, Gandham S. Prasad, Girish Sahni, Jyoti Maharjan, "Whole-Genome Sequence Data Analysis of Anoxybacillus kamchatkensis NASTPD13 Isolated from Hot Spring of Myagdi, Nepal", BioMed Research International, vol. 2021, Article ID 1869748, 11 pages, 2021. https://doi.org/10.1155/2021/1869748

Whole-Genome Sequence Data Analysis of Anoxybacillus kamchatkensis NASTPD13 Isolated from Hot Spring of Myagdi, Nepal

Academic Editor: Fengjie Sun
Received26 Apr 2020
Revised18 Mar 2021
Accepted17 Jun 2021
Published28 Jun 2021

Abstract

Anoxybacillus kamchatkensis NASTPD13 isolated from Paudwar hot spring of Myagdi, Nepal, upon morphological and biochemical analysis revealed to be Gram-positive, straight or slightly curved, rod-shaped, spore-forming, catalase, and oxidase-positive facultative anaerobes. It grows over a wide range of pH (5.0-11) and temperature (37-75°C), which showed growth in different reduced carbon sources such as starch raffinose, glucose, fructose, inositol, trehalose, sorbitol, mellobiose, and mannitol in aerobic conditions. Furthermore, the partial sequence obtained upon sequencing showed 99% sequence similarity in 16S rRNA gene sequence with A. kamchatkensis JW/VK-KG4 and was suggested to be Anoxybacillus kamchatkensis. Moreover, whole-genome analysis of NASTPD13 revealed 2,866,796 bp genome with a G+C content of 41.6%. Analysis of the genome revealed the presence of 102 RNA genes, which includes sequences coding for 19 rRNA and 79 tRNA genes. While the 16S rRNA gene sequence of strain NASTPD13 showed high similarity (>99%) to those of A. kamchatkensis JW/VK-KG4, RAST analysis of NASTPD13 genome suggested that A. kamchatkensis G10 is actually the closest neighbor in terms of sequence similarity. The genome annotation by RAST revealed various genes encoding glycoside hydrolases supporting that it can utilize several reduced carbon sources as observed and these genes could be important for carbohydrate-related industries. Xylanase pathway, particularly the genomic region encoding key enzymes for xylan depolymerization and xylose metabolism, further confirmed the presence of the complete gene in xylan metabolism. In addition, the complete xylose utilization gene locus analysis of NASTPD13 genome revealed all including D-xylose transport ATP-binding protein XylG and XylF, the xylose isomerase encoding gene XylA, and the gene XylB coding for a xylulokinase supported the fact that the isolate contains a complete set of genes related to xylan degradation, pentose transport, and metabolism. The results of the present study suggest that the isolated A. kamchatkensis NASTPD13 containing xylanase-producing genes could be useful in lignocellulosic biomass-utilizing industries where pentose polymers could also be utilized along with the hexose polymers.

1. Introduction

Anoxybacillus was first described as a separate genus by Pikuta et al. [1] in contrast to Bacillus and Geobacillus based on their phenotypic properties, 16S rDNA sequence, and DNA–DNA hybridization experiments. They are defined as rod-shaped, present in pairs or short chains with terminal endospores, Gram-positive, aerotolerant or facultative anaerobe, thermophilic, and alkalophilic or alkalitolerant and grew at a temperature of 37-62°C [2]. Mostly, Anoxybacillus spp. have been reported from hot springs but it has also been isolated from different environments such as animal manure, dairy products, and guts of animals ([3]).

The only genome of A. flavithermus WK1 (PRJNA59135) [4] is completely sequenced, while the draft genomes of other species, including A. ayderensis AB04T (PRJNA258494), A. thermarum AF/04T (PRJNA260786), A. gonensis G2T (PRJNA264351), A. tepidamans PS2 (PRJNA214279), A. kamchatkensisG10 ((PRJNA170961), and 5 strains of A. flavithermus, 25 (PRJNA258119), Kn10 (PRJDB1085), TNO-09.006 (PRJNA169174), and E13T (PRJNA213809), and many more has been reported [5].

Several studies has reported regarding the importance of Anoxybacillus sp., and mostly, the enzymes produced by bacilli that can tolerate the harsh industrial conditions such as alkaline pH and high temperature[2]. Many researches are seeking an alternative for renewable energy generation [6]. Certain features of Anoxybacillus spp. such as able to grow fast, tolerating extreme conditions, and producing various thermozymes are getting attention of many researchers as it can resolve the issues related to biomass hydrolysis [7]. In addition, Anoxybacillus sp. genomes are relatively small in size and can be used as a microbial cell factory for minimizing the difficulty during genome engineering and also serve as a host for certain thermostable enzyme applications [8].

Anoxybacillus kamchatkensis have first been reported from Geyser valley located in Kamchatka peninsula [9]. In the current work, we describe the features of Anoxybacillus kamchatkensis NASTPD13 isolated from Paudwar hot spring of Myagdi, Nepal, [10] and present its annotated draft genome. Additionally, we provide a comparative analysis of the GHs of strain NASTPD13 with four other sequenced Anoxybacillus spp. The current study describes the genomic level of thermophilic bacteria from the hot spring of Nepal and explores its cellular and molecular features in detail.

2. Materials and Methods

2.1. Isolation Details of NASTPD13

1 g of biomat sample collected from Paudwar hot spring [10] was suspended in 10 ml of nutrient broth (Himedia) and incubated at 60°C water bath, The diluted cultures (100 μl) were spread on nutrient agar (HiMedia) plate and incubated at 60°C for 24 to 48 hours [11]. Single colony was picked, and the pure cultures were stored as glycerol stocks for further study.

2.2. Growth Conditions

A. kamchatkensis NASTPD13 was cultured on nutrient agar at 60°C for 18 h. Single colony was transferred into the nutrient broth and incubated at 60°C at 200 rpm for 18 h. Cells were harvested by centrifugation at 10,000 × g for 5 min using a Microfuge® 16 centrifuge (Beckman Coulter, Brea, CA, USA).

2.3. Biochemical Characterization

Apart from some manual biochemical characteristics of NASTPD13 explained in [12], Biolog MicroPlates technique was used for detail biochemical characteristics of NASTPD13. In Biolog MicroPlates technique, each well is coated with a single carbon source and after incubation of the culture shows the detail biochemical characteristics of individual bacteria in a single test. Biolog GN microplates (Biolog, Hayward, CA, USA) were used for NASTPD13. As per the manufacturer protocol, pure colonies of NASTPD13 from nutrient agar were suspended in GN/GP fluid to a specified density. 150 μl of bacterial suspension was pipetted into each well, and the plates were incubated at 60°C up to 24 h. After incubation, plates were read in Biolog MicroStation™ system.

2.4. Antibiotic Susceptibility Testing

The antibiotic susceptibility test of NASTPD13 was performed by standard Kirby Bauer’s disc diffusion tested as per CSLI guidelines [13]. The antibiotics tested were vancomycin, erythromycin, gentamycin, chloramphenicol, penicillin G, cefixime, rifampicin, kanamycin, azithromycin, novobiocin, lomefloxacin, and bacitracin. The antibiotic discs (6 mm diameter) were purchased from HiMedia Laboratory Ltd., Mumbai (India).

2.5. Transmission Electron Microscope (TEM)

NASTPD13 grown in minimum salt medium [10] containing 0.5% of beechwood xylan and for control in nutrient broth were individually harvested. Cells were centrifuged at 12,000 × g for 10 min. After centrifugation, the pellets were washed twice with PBS and then diluted with PBS to a cell density of 108 CFU/ml. Washed cells were further stained with 0.1% (wt/vol) sodium phosphotungstate (PTA; Sigma, 32 USA) for 1 minute. Washed cells were negatively stained with 0.1% (wt/vol) sodium phosphotungstate (PTA; Sigma, USA) in water for 1 min on a carbon-coated copper grid (300 mesh; PolyScience cat# 71150). Grids were observed under JEOL JEM 2100, 200 kV transmission electron microscope (TEM) at a resolution of 0.1 to 0.2 m, and the TEM images were captured [14].

2.6. The Phylogenetic Analysis

NASTPD13 chromosomal DNA was extracted using the DNA Purification Kit (Promega Inc., Madison, WI, USA), and the 16S rRNA gene was amplified using universal primers 27f (5-AGAGTTTGATCCTGGCTCAG-3) and 1492r (5-GGTTACCTTGTTACGACTT-3). The amplified product was purified, using a QIA quick PCR purification kit (Qiagen Inc., San Diego, CA, USA), and sequenced in an ABI Prism 3700 automatic DNA sequencer by the use of a BigDye Terminator Cycle Sequencing kit (Applied Biosystems, Inc., Palo Alto, CA, USA). The sequence obtained was analyzed by the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/). The neighbor-joining phylogenetic tree was constructed with some species of the genus Anoxybacillus based on the 16S rDNA sequences by MEGA 6.0 [15].

2.7. Whole-Genome Sequencing, Assembly, and Annotation

Whole-genome sequence analysis of Anoxybacillus kamchatkensis NASTPD13 was done at CSIR-Institute of Microbial Technology, Chandigarh, India. Genomic DNA was extracted from the isolates using ZR Fungal/Bacterial DNA MiniPrep Kit (Zymo Research Corporation, Orange, CA, USA) as per the manufacturer’s guidelines. Contaminations, quality, and concentration of genomic DNA were checked by using NanoDrop (Thermo Fisher Scientific, MA, USA), and quantification was done by Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA). Library preparation was done by using the Nextera XT sample preparation kit (Illumina, Inc., San Diego, CA, USA) with dual indexing adaptors as per the instruction manual. Sequencing of libraries was done on in-house Illumina MiSeq sequencing platform (Illumina, Inc., San Diego, CA, USA), in a paired end run. Adapter trimming was done by MiSeq Control software (MCS). Quality of obtained reads was checked by using FastQC. Obtained reads were assembled using CLC Genomics Workbench v7 (CLC Bio-Qiagen, Aarhus, Denmark) with minimum contig length as 500 bp. 16S rRNA sequences were extracted from the assembled genomes using RNAmmer 1.2 server [16] and were characterized by Ez-taxon server [17], and tRNA was calculated by tRNAScan-SE ([18]) and CRISPRFinder [19] for finding CRISPR in the genome. The GHs were identified and verified using the dbCAN CAZy ([20]). Average nucleotide identity (ANI) was calculated by JSpecies [21], and Digital DNA-DNA hybridization (dDDH) was calculated by GGDC 2.0, DSMZ [22]. Assembled genomes were then annotated using Rapid Annotation using Subsystem Technology (RAST) [23] and BASys (a web server for automated bacterial genome annotation) [24]. BLAST Ring Image Generator (BRIG) software was used for comparison and visualization of genomes [25].

2.8. Comparative Genome Analysis

The comparative genomics helps to visualize the results and discovers correlations and trends in large datasets which makes understanding and interpretation of the data easier and generate figures for communicating results [26]. Venn Painter [26] and BLAST Ring Image Generator (BRIG) [25] were used for the comparative genome analysis of A. kamchatkensis NASTPD13.

3. Results and Discussions

3.1. Organism Information

Anoxybacillus kamchatkensis NASTPD13 was isolated from Paudwar hot spring of Myagdi district, Nepal [27]. The colonies of NASTPD13 were 1–2 mm in diameter, cream-colored, and regular in shape with round edges. The strain was Gram-positive, rod-shaped containing spherical endospore, and catalase and oxidase positive and has ability to reduce nitrate to nitrite [12]. NASPD13 is a facultative anaerobe and moderately thermophilic that can grow at a wide temperature range of 37–70°C (optimum 55-60°C) and pH 5.0–11.0 (optimum pH 6.5–7.5) [12]. TEM image showed that cells were 0.7–0.4 μm in size (Figure 1). NASTPD13 was able to utilize carbon sources including starch, gelatin, D-glucose, D-raffinose, D-sucrose, D-xylose, D-fructose, L-arabinose, maltose, and D-mannose, whereas there was no growth in the presence of vancomycin, erythromycin, gentamycin, chloramphenicol, penicillin G, cefixime, rifampicin, kanamycin, azithromycin, novobiocin, lomefloxacin, and bacitracin [10]. Phylogenetic analysis of 16S rRNA of NASTPD13 showed clusters together with Anoxybacillus kamchatkensis JW/VK-KG4 [9] (Figure 2).

3.2. TEM Imaging

The TEM image of NASTPD13 grown in nutrient broth (Figure 1(a)) showed rod-shaped (diameter 0.5 μm) bacilli with no flagella, and the cell was surrounded by electron-transparent extracellular zone, and the cytoplasm was homogeneous with a terminal round shape endospore. In the presence of xylan, the cell diameter was 0.2 μm with flagella (Figure 1(b)). An observable change of the electron-transparent material and a disorganized and distorted structure of cytoplasm were seen. Bacteria change their morphology depending on the carbon substrate used ([28]); phenotypic changes seen in NASTPD13 is due to the xylan present in the medium, and utilization of xylan explains the presence of xylan-degrading enzymes in NASTPD13. Various studies on the flagellar system of bacteria have reported that some bacteria have the unique property of switching on and off their flagellar system with some unknown mechanism for their survival in an unfavorable condition. Similar results have also been reported by [2] suggesting that Anoxybacillus spp. SK3-4 and DT3-1 were motile in certain media whereas nonmotile in certain growth media. Interestingly, this study also reports the TEM image of Anoxybacillus kamchatkensis NASTPD13 showing lots of flagella.

3.3. Whole-Genome Sequencing and Species Identification

The genome size of A. kamchatkensis NASTPD13 was 2,866,796 bp with GC 41.6% which assembled into 134 contigs with coverage of 184x. Genome features and assembly statistics of NASTPD13 are listed in Table 1. Species identification and phylogenetic relatedness are best studied by phylogenomic markers, average nucleotide identity (ANI) [29], and digital DNA-DNA hybridization (dDDH) [30]. Table 2 summarizes the pairwise ANI values of Anoxybacillus sp. A. kamchatkensis NASTPD13 with 14 other species. The whole-genome ANI value between Anoxybacillus sp. strain NASTPD13 and close relatives including A. kamchatkensis G10 [31], A. gonenesis G2 [32], A. mongoliensis MB4 [33], and A. thermarum AF/04 (T) [34] as 99.4%, 98.9%, 94.3%, and 93.6%, respectively. The dDDH of strain NASTPD13 neighbor strains was 95.1%, 85.5%, 73.9%, and 67.3%, respectively. A. kamchatkensis NASTPD13 showed the highest ANI of 99.4% and dDDH of 95.1% with A. kamchatkensis G10. The whole-genome shotgun project was deposited at DDBJ/EMBL/GenBank under the accession number NQLB000000000. The NCBI BioProject accession number is PRJNA397640, and the biosample accession number is SAMN07482658.


AttributeNASTPD13A. gonensis G2A. mongoliensis MB4A. kamchatkensis G10A. thermarum AF/04

Genome size (bp)2,866,7962,803,6682,807,5162,858,6572,736,908
G+C content (%)41.641.741.741.342.0
N50103,7592,803,66878,095130,03637,018
No. of contigs13419365159
Coding sequence29922808279829103004
Pseudogenes7556129107
Genes (RNA)10210611369117
rRna19243299
tRna7978775674
ncRNA44441
CRISPR arrays12112
Accession numberNQLB000000000JRZG01000000MRZM00000000ALJT00000000JXTH00000000


NASTPD13AB04AF/04BAA-2555BC01DSM15939DSM27374G2G10K1P3H1BPS2MB4SK3-4WK1

NASTPD1310094.2993.672.894.274.987.098.999.485.572.772.994.393.788.2
AB0494.310094.473.197.375.287.494.294.186.472.772.997.297.488.0
AF/0493.794.610073.993.974.887.193.693.586.672.773.093.994.188.0
BAA-255571.671.871.910071.873.971.271.671.671.988.275.271.871.871.8
BC0194.297.498.973.010074.987.293.994.086.272.572.896.596.887.7
DSM1593974.875.074.875.874.810074.274.674.974.475.976.975.074.975.1
DSM2737487.187.487.172.387.374.010087.287.085.972.372.688.287.488.4
G298.894.393.772.894.274.587.310098.785.972.672.894.494.088.4
G1099.494.293.572.894.274.887.098.710085.572.872.994.393.788.1
K185.586.486.373.286.374.686.085.885.610072.372.686.086.486.4
P3H1B72.472.372.598.172.275.472.272.272.572.310078.272.272.472.4
PS272.772.672.877.872.776.672.472.772.772.678.310072.772.672.7
MB494.297.293.672.996.575.088.294.194.186.072.772.910096.688.0
SK3-493.697.594.073.296.874.687.493.893.586.472.672.996.710088.1
WK188.387.987.873.087.775.388.488.288.286.373.072.988.08.1100

The reference protein sequence is denoted as 100%. AB04: A. ayderensis AB04 [5]; AF/04: A. thermarum AF/04 [42]; BAA-2555: Anoxybacillus geothermalis [43]; BCO1: Anoxybacillus sp. BCO1 [44]; DSM15939: Anoxybacillus amylolyticus DSM15939 [42]; DSM27374: Anoxybacillus suryakundensis strain DSM 27374 [45]; G2: A. gonensis G2 [32]; G10: A. kamchatkensis G10 [31]; K1: Anoxybacillus pushchinoensis K1 [46]; P3H1B: Anoxybacillus sp. P3H1B [47]; PS2: A. tepidamans PS2 [48]; MB4: Anoxybacillus mongoliensis MB4 [33]; SK3-4: Anoxybacillus sp. SK3-4 [49]; WK1: A. flavithermus WK1 [4].
3.4. Annotation of NASTPD13

RAST predicted 3169 genes among which 2992 were coding sequences (CDS) in the NASTPD13 genome. The subsystem distribution revealed a total of 436 subsystems among which amino acid and derivative subsystem featured with largest number 337 CDS. Other major subsystems annotated were carbohydrate 330, protein metabolism 232, cofactors, vitamins, prosthetic groups and pigments 232, RNA metabolism 139, sulfur metabolism 23, and motility or chemotaxis 84 (Figure 3). BASys annotated 3160 genes that were encoding various proteins and displayed in the form of circular DNA for easy representation of genomic data (Figure 4). The highest amino acid residue content was predicted for leucine followed by valine, alanine, glutamic acid, and isoleucine (Figure 5).

3.5. Comparative Analysis of NASTPD13

NASTPD13 genome and other four Anoxybacillus spp. shared 1469 genes as shown in the Venn diagram (Figure 6). NASTPD13 and G10 shared 3237 CDS exclusively. In addition, 421, 60, and 38 genes are shared between NASTPD13 and G2, AF/04, and WK1, respectively, showing that NASTPD13 apart from G10 are more closely related to G2 than the other strains.

Further, comparing the genome of G10, G2, MB4, and AF/04 (T) with NASTPD13 as a reference strain is sequenced in this study. The output image showed similarity between NASTPD13 and the other comparative four genomes as a set of concentric ring where BLAST matches are seen in different sliding colors defining the percentage of similarity among them. Even though small differences between the genome are not apparent in the BRIG concentric image, but it can be clearly seen in Figure 7 that Anoxybacillus kamchatkensis G10 genome is again the most similar to NASTPD13 among the strains analyzed.. The BRIG image also displays the bacteriophage genes in the NASTPD13 genome ranging in size between 33.7 kb and 64.8 kb. Phage 1 located at 1789826-1845786 bp, consisting of 229 ORFs, had the highest homology to phage Thermu_OH2_0127842 (28) of 33.7 kb; Phage 2 located at 2480822-2529216, consisting of 319 ORFs, had highest homology to Phage Bacter_Lily_NC_028841 (12) of 48.3 kb; and Phage 3 located at 2800662-2865493, consisting of 396 ORFs, had highest homology to Phage Entero_Lambda_NC_001416 (35) of 64.8 kb (Figure 7). Phage contributes to the evolution of their host bacteria by horizontal transfer of genes [35].

3.6. Analysis of the GHs in A. kamchatkensis NASTPD13 and Other Anoxybacillus Genomes

The photosynthetically fixed carbon in the plant cell wall is rich in organic carbon compound that can be converted to renewal biofuel. The degradation of plant polymers for various biotechnological uses requires combined action of several enzymes [36]. Glycoside hydrolase are potentially involved in the degradation of plant biomass [37]. The dbCAN CAZy database detected 21 genes in NASTPD13 genome encoding GH enzyme that belongs to GH families: 13, 23, 31, 32, 36, 42, 52, 65, and 109.The GHs are grouped based on their catalytic ability (Table 3). In NASTPD13, nine GH enzymes were active 275 on α-chain polysaccharides whereas five were specific for β-linked polysaccharides (i.e., cellulose and xylan). Along with GH families that hydrolyze the β-1,4-xyloside linkages of xylan, the CAZy database showed that NASTPD13 contains many other genes encoding for lignocellulose-degrading enzymes. The result showed 107 different CAZyme genes, encompassing all six CAZy families: 21 glycoside hydrolases (GHs), 21 glycosyltransferases (GTs), 15 carbohydrate esterases (CEs), 1 polysaccharide lyases (PLs), 4 auxiliary activities (AAs), and 23 carbohydrate-binding modules (CBMs). CBM (carbohydrate binding modules) families are capable of binding to xylan. The finding suggests that presences of various glycohydrolase genes in NASTPD13 hydrolyze the glycosidic bonds present in complex sugars.


GHsEnzymesNASTPD13G10MB4G2AF/04

131,4-Alpha-glucan branching enzyme100.00100.0099.7699.5496.14
13Oligo-1,6-glucosidase100.00100.0096.7699.5496.14
13Trehalose-6-phosphate hydrolase100.00100.0098.3899.1098.74
13Alpha-glucosidase100.0099.8299.1099.8298.92
13Cytoplasmic alpha-amylase100.00100.0093.96
13Pullulanase100.00100.0095.7499.3293.19
13Maltodextrin glucosidase100.0099.4598.39
23Murein hydrolase100.00100.0097.0493.10
31Alpha-glucosidase100.0099.8299.1099.8298.92
32Sucrose-6-phosphate hydrolase100.00100.0097.9997.9991.92
36Alpha-galactosidase100.00100.0055.13
42Beta-galactosidase100.00100.0099.69
52Beta-xylosidase100.00100.0089.0899.43
65Maltose phosphorylase100.00100.0096.4199.87
109Oxidoreductase100.00100.0090.9999.70

NASTPD13: A. kamchatkensis NASTPD13; G10: A. kamchatkensis G10; MB4: A. mongoliensis MB4; G2: A. gonensis G2; AF/04: A. thermarum AF/04.
3.7. Xylose Metabolism by A. kamchatkensis NASTPD13

Hemicellulosic sugars, especially D-xylose, are abundant in lignocellulosic biomass that can be recovered from by various physical and chemical pretreatments of the biomass, and also, D-xylose yields better from the biomass than D-glucose from cellulose. Xylose is one of the major fermentable sugars in nature after glucose; therefore, fermentation of xylose is necessary for producing biofuels, such as ethanol, from lignocellulosic biomass [38].

Hemicellulosic sugars has also been reported as a feedstock for production of ethanol and other chemicals but D-xylose is not so readily utilized as D-glucose due to the complex biochemical pathways for pentose and hexose metabolisms [39]. Various pathways have been employed by prokaryotes and eukaryotes for pentose assimilation [40].

Thermophilic bacteria capable of fermenting both glucose and xylose are on high industrial demand (Yejun [41]). Thermophiles like Anoxybacillus sp. are relevant to the industry, because they produce thermozymes that match the harsh industrial process [2]. The observed ability of A. kamchatkensis NASTPD13 to utilize xylose in our previous study [10] developed our keen interest in the investigation of the xylose metabolism pathway and performed the whole-genome sequence of A. kamchatkensis NASTPD13. NASTPD13 genome revealed the presence of genes related to xylose metabolism such as xylose-binding protein xylG (NAST_1706), xylF (NAST_1707), xylH (NAST_1705), ABC transporter system-associated protein (NAST_622, NAST_623, NAST_624), transcription regulator (NAST_627), xylose isomerase (NAST_619) encoding gene xylA, and gene xylB coding for a xylulokinase (NAST_618) (Figure 8). XylA deals with isomerization of xylose to xylulose, while XylB convert xylulose to xylulose-5-phosphate. NASTPD13 xylan ABC transporter system and β-xylosidase continue to further step of xylose metabolism.

To date, in the family Bacillaceae, genome sequencing of Geobacillus (>80 projects) and Bacillus (>1500 projects) have been undertaken and have been registered in the NCBI BioProject database, whereas genomic studies on Anoxybacillus are rather limited (16 registered projects) [5]. These findings on A. kamchatkensis NASTPD13 genome add knowledge on the industrially important Anoxybacillus genomes. This study provides additional information on the xylose metabolism in Anoxybacillus kamchatkensis that can further help in the detailed study of the species.

4. Conclusion

Anoxybacillus kamchatkensis NASTPD13, isolated from Paudwar hot spring of Nepal, is Gram-positive with terminal endospore rod-shaped bacilli [10]. The TEM image of NASTPD13 grown in the xylan-containing medium showed major changes in their cell morphology along with the interesting property of flagella production in a certain growth medium such as xylan. The genome size of NASTPD13 is 2.8 Mb. Comparison of the NASTPD13 genome with the closely related strain of Anoxybacillus spp. against various databases showed the closed relationship with Anoxybacillus kamchatkensis G10. Previous study explained in [10] related to NASTPD13 xylanase and combining the data from this study related to presence of genes related to metabolism of Xylan showed the presence of xylanase and β xylosidase gene that are essential for complete degradation of xylan to pentose sugar xylose. The gene related to utilization of xylose (xylG, xylF, and xylH) is also present in Anoxybacillus kamchatkensis NASTPD13. The presence of all the genes related to the metabolism of xylan in Anoxybacillus kamchatkensis NASTPD13 genome emphasizes its industrial importance along with promising sources of lignocellulase in second-generation biofuel production.

Data Availability

Anoxybacillus kamchatkensis NASTPD13 whole-genome sequence was deposited in DDBJ/EMBL/GenBank under the accession number NQLB000000000. The NCBI BioProject accession number is PRJNA397640, and the biosample accession number is SAMN07482658.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

PY isolated the strain NASTPD13 and drafted the manuscript. SS carried out the sequencing part. PY and JM were involved in data analysis. PY, JM, TB, LS, GSP and GS conceived and designed the experiment. All authors reviewed and approved the manuscript.

Acknowledgments

This work was supported by CSIR-TWAS Sandwich Postgraduate Fellowship Program (CSIR-TWAS PG, 2014). We thank Institute of CSIR-Institute of Microbial Technology, Microbial Type Culture Collection (MTCC) Chandigarh for all the technical support.

References

  1. E. Pikuta, A. Lysenko, N. Chuvilskaya et al., “Anoxybacillus pushchinensis gen. nov., sp. nov., a novel anaerobic, alkaliphilic, moderately thermophilic bacterium from manure, and description of Anoxybacillus flavitherms comb. nov,” International Journal of Systematic and Evolutionary Microbiology, vol. 50, no. 6, pp. 2109–2117, 2000. View at: Publisher Site | Google Scholar
  2. K. M. Goh, H. M. Gan, K. G. Chan et al., “Analysis of Anoxybacillus genomes from the aspects of lifestyle adaptations, prophage diversity, and carbohydrate metabolism,” PLoS One, vol. 9, no. 3, article e90549, 2014. View at: Publisher Site | Google Scholar
  3. X.-Q. Zhang, Z.-L. Zhang, N. Wu, X.-F. Zhu, and M. Wu, “Anoxybacillus vitaminiphilus sp. nov., a strictly aerobic and moderately thermophilic bacterium isolated from a hot spring,” International Journal of Systematic and Evolutionary Microbiology, vol. 63, no. Part_11, pp. 4064–4071, 2013. View at: Publisher Site | Google Scholar
  4. J. H. Saw, B. W. Mountain, L. Feng et al., “Encapsulated in silica: genome, proteome and physiology of the thermophilic bacterium Anoxybacillus flavithermus WK1,” Genome Biology, vol. 9, no. 11, p. R161, 2008. View at: Publisher Site | Google Scholar
  5. A. O. Belduz, S. Canakci, K. G. Chan et al., “Genome sequence of Anoxybacillus ayderensis AB04T isolated from the Ayder hot spring in Turkey,” Standards in Genomic Sciences, vol. 10, no. 1, pp. 1–12, 2015. View at: Publisher Site | Google Scholar
  6. E. S. Shuba and D. Kifle, “Microalgae to biofuels: 'promising' alternative and renewable energy, review,” Renewable and Sustainable Energy Reviews, vol. 81, pp. 743–755, 2018. View at: Publisher Site | Google Scholar
  7. G.-Q. Chen and X.-R. Jiang, “Next generation industrial biotechnology based on extremophilic bacteria,” Current Opinion in Biotechnology, vol. 50, pp. 94–100, 2018. View at: Publisher Site | Google Scholar
  8. K. M. Goh, U. M. Kahar, Y. Y. Chai et al., “Recent discoveries and applications of Anoxybacillus,” Applied Microbiology and Biotechnology, vol. 97, no. 4, pp. 1475–1488, 2013. View at: Publisher Site | Google Scholar
  9. V. V. Kevbrin, K. Zengler, A. M. Lysenko, and J. Wiegel, “Anoxybacillus kamchatkensis sp. nov., a novel thermophilic facultative aerobic bacterium with a broad pH optimum from the Geyser valley, Kamchatka,” Extremophiles: life under extreme conditions, vol. 9, no. 5, pp. 391–398, 2005. View at: Publisher Site | Google Scholar
  10. P. Yadav, J. Maharjan, S. Korpole et al., “Production, purification, and characterization of thermostable alkaline xylanase from Anoxybacillus kamchatkensis NASTPD13,” Frontiers in Bioengineering and Biotechnology, vol. 6, p. 65, 2018. View at: Publisher Site | Google Scholar
  11. H. N. Akmar, I. Asma, B. Venugopal, L. Y. Latha, and S. Sasidharan, “Identification of appropriate sample and culture method for isolation of new thermophilic bacteria from hot spring,” African Journal of Microbiology Research, vol. 5, no. 3, pp. 217–221, 2011. View at: Google Scholar
  12. P. Yadav, S. Korpole, G. S. Prasad et al., “Morphological, enzymatic screening and phylogenetic analysis of thermophilic bacilli isolated from five hot springs of Myagdi, Nepal,” Journal of Applied Biology and Biotechnology, vol. 6, no. 3, pp. 1–8, 2018. View at: Publisher Site | Google Scholar
  13. CLSI, “Clinical and Laboratory Standards Institute (CLSI),” in Performance Standards for Antimicrobial Susceptibility Testing, Clinical and Laboratory Standards Institute, Wayne, 28th edition, 2018. View at: Google Scholar
  14. P. Baindara, V. Chaudhry, G. Mittal et al., “Characterization of the antimicrobial peptide penisin, a class Ia novel lantibiotic from Paenibacillus sp. strain A3,” Antimicrobial Agents and Chemotherapy, vol. 60, no. 1, pp. 580–591, 2016. View at: Publisher Site | Google Scholar
  15. R. D. Kamble and A. R. Jadhav, “Isolation, purification, and characterization of xylanase produced by a new species of Bacillus in solid state fermentation,” International Journal of Microbiology, vol. 2012, Article ID 683193, 8 pages, 2012. View at: Publisher Site | Google Scholar
  16. K. Lagesen, P. Hallin, E. A. Rødland, H.-H. Stærfeldt, T. Rognes, and D. W. Ussery, “RNAmmer: consistent and rapid annotation of ribosomal RNA genes,” Nucleic Acids Research, vol. 35, no. 9, pp. 3100–3108, 2007. View at: Publisher Site | Google Scholar
  17. O.-S. Kim, Y.-J. Cho, K. Lee et al., “Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species,” International Journal of Systematic and Evolutionary Microbiology, vol. 62, no. Part_3, pp. 716–721, 2012. View at: Publisher Site | Google Scholar
  18. T. M. Lowe and S. R. Eddy, “tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence,” Nucleic Acids Research, vol. 25, no. 5, pp. 955–964, 1997. View at: Publisher Site | Google Scholar
  19. I. Grissa, G. Vergnaud, and C. Pourcel, “CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats,” Nucleic Acids Research, vol. 35, no. Web Server, pp. W52–W57, 2007. View at: Publisher Site | Google Scholar
  20. B. L. Cantarel, P. M. Coutinho, C. Rancurel, T. Bernard, V. Lombard, and B. Henrissat, “The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics,” Nucleic Acids Research, vol. 37, no. Database, pp. D233–D238, 2009. View at: Publisher Site | Google Scholar
  21. M. Richter, R. Rosselló-Móra, F. Oliver Glöckner, and J. Peplies, “JSpeciesWS: a web server for prokaryotic species circumscription based on pairwise genome comparison,” Bioinformatics, vol. 32, no. 6, pp. 929–931, 2016. View at: Publisher Site | Google Scholar
  22. J. P. Meier-Kolthoff, A. F. Auch, H. P. Klenk, and M. Goker, “Genome sequence-based species delimitation with confidence intervals and improved distance functions,” BMC Bioinformatics, vol. 14, no. 1, pp. 1–14, 2013. View at: Publisher Site | Google Scholar
  23. R. K. Aziz, D. Bartels, A. A. Best et al., “The RAST server: rapid annotations using subsystems technology,” BMC Genomics, vol. 9, no. 1, p. 75, 2008. View at: Publisher Site | Google Scholar
  24. G. H. Van Domselaar, P. Stothard, S. Shrivastava et al., “BASys: a web server for automated bacterial genome annotation,” Nucleic Acids Research, vol. 33, no. Web Server, pp. W455–W459, 2005. View at: Publisher Site | Google Scholar
  25. N. F. Alikhan, N. K. Petty, N. L. Ben Zakour, and S. A. Beatson, “BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons,” BMC Genomics, vol. 12, no. 1, p. 402, 2011. View at: Publisher Site | Google Scholar
  26. G. Lin, J. Chai, S. Yuan et al., “VennPainter: a tool for the comparison and identification of candidate genes based on Venn diagrams,” PLoS One, vol. 11, no. 4, article e0154315, 2016. View at: Publisher Site | Google Scholar
  27. P. Yadav, B. Rimal, J. Maharjan, S. Shrestha, L. Sreerama, and T. Bhattarai, “Physico-chemical analysis of hot springs of Myagdi District, Nepal for the balneotherapy and drinking purposes,” International Journal of Advanced Biotechnology and Research, vol. 8, no. 4, pp. 699–709, 2017. View at: Google Scholar
  28. M. V. Sizova, J. A. Izquierdo, N. S. Panikov, and L. R. Lynd, “Cellulose- and xylan-degrading thermophilic anaerobic bacteria from biocompost,” Applied and Environmental Microbiology, vol. 77, no. 7, pp. 2282–2291, 2011. View at: Publisher Site | Google Scholar
  29. N. Han, Y. Qiang, and W. Zhang, “ANItools web: a web tool for fast genome comparison within multiple bacterial strains,” Database: The Journal of Biological Databases and Curation, vol. 2016, article baw084, 2016. View at: Publisher Site | Google Scholar
  30. A. F. Auch, M. von Jan, H.-P. Klenk, and M. Göker, “Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison,” Standards in Genomic Sciences, vol. 2, no. 1, pp. 117–134, 2010. View at: Publisher Site | Google Scholar
  31. S. J. Lee, Y. J. Lee, N. Ryu et al., “Draft genome sequence of the thermophilic bacterium Anoxybacillus kamchatkensis G10,” Journal of Bacteriology, vol. 194, no. 23, pp. 6684-6685, 2012. View at: Publisher Site | Google Scholar
  32. A. O. Belduz, S. Dulger, and Z. Demirbag, “Anoxybacillus gonensis sp. nov., a moderately thermophilic, xylose-utilizing, endospore-forming bacterium,” International Journal of Systematic and Evolutionary Microbiology, vol. 53, no. 5, pp. 1315–1320, 2003. View at: Publisher Site | Google Scholar
  33. P. Mittal, R. Saxena, and V. K. Sharma, “Draft genome sequence of Anoxybacillus mongoliensis strain MB4, a sulfur-utilizing aerobic thermophile isolated from a hot spring in Tattapani, Central India,” Genome Announcements, vol. 5, no. 9, pp. e01709–e01716, 2017. View at: Publisher Site | Google Scholar
  34. A. Poli, I. Romano, P. Cordella, P. Orlando, B. Nicolaus, and C. Ceschi Berrini, “Anoxybacillus thermarum sp. nov., a novel thermophilic bacterium isolated from thermal mud in Euganean hot springs, Abano Terme, Italy,” Extremophiles: life under extreme conditions, vol. 13, no. 6, pp. 867–874, 2009. View at: Publisher Site | Google Scholar
  35. I. B. Naser, M. M. Hoque, M. A. Nahid, T. M. Tareq, M. K. Rocky, and S. M. Faruque, “Analysis of the CRISPR-Cas system in bacteriophages active on epidemic strains of _Vibrio cholerae_ in Bangladesh,” Scientific Reports, vol. 7, no. 1, p. 14880, 2017. View at: Publisher Site | Google Scholar
  36. P. Longoni, M. Rodolfi, L. Pantaleoni et al., “Functional analysis of the degradation of cellulosic substrates by a Chaetomium globosum endophytic isolate,” Applied and Environmental Microbiology, vol. 78, no. 10, pp. 3693–3705, 2012. View at: Publisher Site | Google Scholar
  37. S. E. Blumer-Schuette, D. L. Lewis, and R. M. Kelly, “Phylogenetic, microbiological, and glycoside hydrolase diversities within the extremely thermophilic, plant biomass-degrading genus Caldicellulosiruptor,” Applied and Environmental Microbiology, vol. 76, no. 24, pp. 8084–8092, 2010. View at: Publisher Site | Google Scholar
  38. D. N. Moysés, V. C. B. Reis, J. R. M. de Almeida, L. M. P. de Moraes, and F. A. G. Torres, “Xylose fermentation by Saccharomyces cerevisiae: challenges and prospects,” International Journal of Molecular Sciences, vol. 17, no. 3, pp. 207–207, 2016. View at: Publisher Site | Google Scholar
  39. T. W. Jeffries, “Utilization of xylose by bacteria, yeasts, and fungi,” in Pentoses and Lignin. Advances in Biochemical Engineering/Biotechnology, Springer, Berlin, Heidelberg, 1983. View at: Google Scholar
  40. X. Xiong, X. Wang, and S. Chen, “Engineering of a xylose metabolic pathway in Rhodococcus strains,” Applied and environmental microbiology, vol. 78, no. 16, pp. 5483–5491, 2012. View at: Publisher Site | Google Scholar
  41. Y. Han, V. Agarwal, D. Dodd et al., “Biochemical and structural insights into xylan utilization by the thermophilic bacterium Caldanaerobius polysaccharolyticus,” The Journal of Biological Chemistry, vol. 42, no. 42, pp. 34946–34960, 2012. View at: Google Scholar
  42. A. Poli, E. Esposito, L. Lama et al., “_Anoxybacillus amylolyticus_ sp. nov., a thermophilic amylase producing bacterium isolated from Mount Rittmann (Antarctica),” Systematic and Applied Microbiology, vol. 29, no. 4, pp. 300–307, 2006. View at: Publisher Site | Google Scholar
  43. S. Filippidou, M. Jaussi, T. Junier et al., “Genome sequence of Anoxybacillus geothermalis strain GSsed3, a novel thermophilic endospore-forming species,” Genome Announcement, vol. 3, no. 3, 2015. View at: Publisher Site | Google Scholar
  44. B. K. Patel, “Draft genome sequence of Anoxybacillus strain BCO1, isolated from a thermophilic microbial mat colonizing the outflow of a bore well of the Great Artesian Basin of Australia,” Genome Announcement, vol. 3, no. 1, 2015. View at: Publisher Site | Google Scholar
  45. K. Deep, A. Poddar, and S. K. Das, “Anoxybacillus suryakundensis sp. nov, a moderately thermophilic, alkalitolerant bacterium isolated from hot spring at Jharkhand, India,” PloS one, vol. 8, no. 12, p. e85493, 2013. View at: Publisher Site | Google Scholar
  46. E. Pikuta, D. Cleland, and J. Tang, “Aerobic growth of Anoxybacillus pushchinoensis K1: emended descriptions of A. pushchinoensis and the genus Anoxybacillus,” International journal of systematic and evolutionary microbiology, vol. 53, no. 5, 2003. View at: Google Scholar
  47. J. Adelskov and B. K. Patel, Draft genome sequence of Anoxybacillus strain P3H1B isolated from an Australian gas field environment, National Center for Biotechnology Information), 2015.
  48. C. A. Cordeiro, M. L. Martins, A. B. Luciano, and R. F. Silva, “Production and properties of xylanase from thermophilic Bacillus sp,” Brazilian Archives of Biology and Technology, vol. 45, no. 4, pp. 413–418, 2002. View at: Publisher Site | Google Scholar
  49. U. M. Kahar, K.-G. Chan, M. M. Salleh, S. M. Hii, and K. M. Goh, “A high molecular-mass Anoxybacillus sp. SK3-4 amylopullulanase: characterization and its relationship in carbohydrate utilization,” International Journal of Molecular Sciences, vol. 14, no. 6, pp. 11302–11318, 2013. View at: Publisher Site | Google Scholar

Copyright © 2021 Punam Yadav et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views99
Downloads88
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.