Taxonomic and Functional Annotation of Termite Degraded Butea monosperma (Lam.) Kuntze (Flame of the Forest)
Ashwani Kumar1, *, Preeti Vyas1, Muneer A. Malla2, Anamika Dubey1
1 Department of Botany, Metagenomics and Secretomics Research Laboratory, Dr. Harisingh Gour University (A Central University), Sagar-470003, (M.P.), India
2 Department of Zoology, Dr. Harisingh Gour University (A Central University), Sagar-470003, (M.P.), India
Butea monosperma is an economically and medicinally important plant that grows all over India, however, the plant is highly susceptible to termite attack. The present study unravelled the bacterial community composition and their functional attributions from the termite degraded Butea.
Total genomic DNA from termite degraded Butea monosperma samples was extracted and subjected to sequencing on Illumina's Miseq. The raw and unassembled reads obtained from high-throughput sequencing were used for taxonomic and functional profiling using different online and stand-alone softwares. Moreover, to ascertain the effect of different geographical locations and environmental factors, comparative analysis was performed using four other publically available metagenomes.
The higher abundance of Actinobacteria (21.27%), Proteobacteria (14.18%), Firmicutes (10.46%), and Bacteroidetes (4.11%) was found at the phylum level. The genus level was dominated by Bacillus (4.33%), Gemmatimonas (3.13%), Mycobacterium (1.82%), Acidimicrobium (1.69%), Thermoleophilum (1.23%), Nocardioides (1.44%), Terrimonas and Acidithermus (1.09%) and Clostridium (1.05%). Functional annotation of the termite degraded B. monosperma metagenome revealed a high abundance of ammonia oxidizers, sulfate reducers, dehalogenators, nitrate reducers, sulfide oxidizers, xylan degraders, nitrogen fixers and chitin degraders.
The present study highlights the significance of the inherent microbiome of the degraded Butea shaping the microbial communities for effective degradation of biomass and different environmental toxicants. The unknown bacterial communities present in the sample can serve as enzyme sources for lignocelluloses degradation for biofuel production.
Keywords: Butea monosperma, Termites, DNA sequencing, Lignocellulose, Biofuel, Environmental toxicants.
open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: (https://creativecommons.org/licenses/by/4.0/legalcode). This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
* Address correspondence to this author at the Department of Botany, Metagenomics and Secretomics Research Laboratory, Dr. Harisingh Gour University (A Central University), Sagar-470003,(M.P.), India;
Tel: +917697432012; Fax: 7582264236; E-mail: firstname.lastname@example.org
Butea monosperma (Lam.) Kuntze is native to the tropical and sub-tropical parts of the Indian Subcontinent and Southeast Asia, and is distributed across India, Bangladesh, Nepal, Sri Lanka, Myanmar, Thailand, Laos, Cambodia, Vietnam, Malaysia, and western Indonesia. Butea is commonly known as ‘Flame of the forest’ and belongs to the family Fabaceae [1Patil MV, Pawar S, Patil DA. Ethnobotany of Butea monosperma (Lam.) Kuntze in North Maharashtra, India. Medicine (Baltimore) 2006; 5: 323-5.]. Locally, this plant is known as Palash, Palas, Mutthuga, Bijasneha, Khakara, Dhak Chichra, Bastard Teak, Bengal Kino, and Nourouc [2Chopra RN. Indigenous Drugs of India Their Medical and Economic Aspects 20 British Indian Street 1933.]. The crude extracts of various parts of this plant has been reported to have antifungal, anti-inflammatory, antibacterial and hypoglycemic activities. The plant has also been found to possess astringent, tonic, diuretics, and aphrodisiac properties. The widespread uses of Butea in traditional systems of medicine have resulted in their extensive chemical analysis for their bio-active components [3Michel T, Destandau E, Pecher V, et al. Two-step Centrifugal Partition Chromatography (CPC) fractionation of Butea monosperma (Lam.) biomarkers. Separ Purif Tech 2011; 80: 32-7.[http://dx.doi.org/10.1016/j.seppur.2011.04.002] ]. The plant holds a significant place because of its medicinal and other miscellaneous uses of economic importance [4Bhalerao S. A comprehensive review : Butea monosperma (LAM.) Kuntze. Int J Curr Res 2016; 4(2): 400007.]. Morphologically, the plant exhibits distinct features with large sized flowers that are arranged in a rigid racemes of 15 cm length. 3 flowers together form the tumid nodes of the dark olive-green velvety rachis. Flowers are used as vegetables by the tribal's and are also used for the preparation of dye for colouring garments and for making skin antiseptic ointments [4Bhalerao S. A comprehensive review : Butea monosperma (LAM.) Kuntze. Int J Curr Res 2016; 4(2): 400007.]. The above mentioned reasons make Butea an important and much sacred plant species.
Termites are ubiquitous and major constituents of the forest ecosystem in the tropical and subtropical areas [5Shanbhag RR, Sundararaj R. Host range pest status and distribution of wood destroying termites of India. J Trop Asian Entomol 2012; 2: 12-27.[http://dx.doi.org/10.1007/s00040-017-0572-9] ]. Apart from Antarctica, termites are found in almost all the continents [6Buczkowski G, Bertelsmeier C. Invasive termites in a changing climate: A global perspective. Ecol Evol 2017; 7(3): 974-85.[http://dx.doi.org/10.1002/ece3.2674] [PMID: 28168033] ] and are well known for their capacity to damage, and destroy wood and wooden products of human homes, building materials, forests, agriculture crops and other commercial products [7Peralta RCG, Menezes EB, Carvalho AG, Aguiar-Menezes E de L. Wood consumption rates of forest species by subterranean termites (isoptera) under field conditions 2004.[http://dx.doi.org/10.1590/S0100-67622004000200015] ]. Termites, besides being small insects, have a successful evolutionary history mainly because of their sociality and capability to efficiently digest lignocellulose [8Bourguignon T, Lo N, Cameron SL, et al. The evolutionary history of termites as inferred from 66 mitochondrial genomes. Mol Biol Evol 2015; 32(2): 406-21.[http://dx.doi.org/10.1093/molbev/msu308] [PMID: 25389205] ]. Presently, these tiny creatures contribute more than 95% of the insect biomass in the tropical soils [9Breznak JA, Brune A. Role of microorganisms in the digestion of lignocellulose by termites. Annu Rev Entomol 2003; 39(1): 453-87.[http://dx.doi.org/10.1146/annurev.en.39.010194.002321] ]. About 3,000 species of termites have been described so far, however, comparatively only few species pose threat to wooden structures. The principal component of their diet is wood, which is also the dominant structural element of the building construction industry. Termites have the ability to degrade 65-99% and 5-83% of the wood-cellulose and hemicelluloses and lignin within 24 hours under natural conditions [10Brune A. Microbiology: Woodworker’s digest. Nature 2007; 450(7169): 487-8.[http://dx.doi.org/10.1038/450487a] [PMID: 18033283] , 11Chandrasekharaiah M, Thulasi A, Bagath M, et al. Molecular cloning, expression and characterization of a novel feruloyl esterase enzyme from the symbionts of termite(Coptotermes formosanus) gut. BMB Rep 2011; 44(1): 52-7.[http://dx.doi.org/10.5483/BMBRep.2011.44.1.52] [PMID: 21266107] ]. In fact, termites possess the inimitable characteristics in order to handle lignin-carbohydrate complexes efficiently. A variety of lignocellulolytic enzymes do occur in the termite gut [12Chaffron S, von Mering C. Termites in the woodwork. Genome Biol 2007; 8(11): 229.[http://dx.doi.org/10.1186/gb-2007-8-11-229] [PMID: 18036268] , 13Zhang D, Lax AR, Bland JM, Allen AB. Characterization of a new endogenous endo-β-1,4-glucanase of Formosan subterranean termite (Coptotermes formosanus). Insect Biochem Mol Biol 2011; 41(4): 211-8.[http://dx.doi.org/10.1016/j.ibmb.2010.12.006] [PMID: 21195179] ]. Termites, in spite of their small body size, harbor an abundant and astonishingly diverse gut microbiome, which represents one of the most fascinating examples of symbiosis both within the microbes, as well as between animal and microbes [14Hongoh Y, Ohkuma M, Kudo T. Molecular analysis of bacterial microbiota in the gut of the termite Reticulitermes speratus (Isoptera; Rhinotermitidae). FEMS Microbiol Ecol 2003; 44(2): 231-42.[http://dx.doi.org/10.1016/S0168-6496(03)00026-6] [PMID: 19719 640] ]. It is widely accepted that lignocellulose digestion in termites is closely related to the host as well as the specific flora of symbiotic microbes [15KUDO T. Termite-Microbe symbiotic system and its efficient degradation of lignocellulose. Biosci Biotechnol Biochem 2009; 73: 2561-7.[http://dx.doi.org/10.1271/bbb.90304] [PMID: 19966490] , 16Scharf ME, Tartar A. Termite digestomes as sources for novel lignocellulases. Biofuels Bioprod Biorefin 2008; 540-52.[http://dx.doi.org/10.1002/bbb.107] ]. Various evidences do suggest that the lignin structure was modified in termites and that the modification process begins in the foregut and continues in the mid gut [17Geib SM, Filley TR, Hatcher PG, et al. Lignin degradation in wood-feeding insects. Proc Natl Acad Sci USA 2008; 105(35): 12932-7.[http://dx.doi.org/10.1073/pnas.0805257105] [PMID: 18725643] , 18Warnecke F, Luginbühl P, Ivanova N, et al. Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature 2007; 450(7169): 560-5.[http://dx.doi.org/10.1038/nature06269] [PMID: 18033299] ]. From the perspectives of both basic research as well as biotechnology, there arises a substantial interest in accomplishing a clear perceptive of the diversity of biological mechanisms that are employed during lignocellulose degradation [19Singh JK, Vyas P, Dubey A, et al. Assessment of different pretreatment technologies for efficient bioconversion of lignocellulose to ethanol. Frontiers in bioscience (Scholar edition) 2018; 10: 350-71.[http://dx.doi.org/10.2741/s521] , 20Tiwari G, Sharma A, Kumar A, Sharma S. Assessment of microwave-assisted alkali pretreatment for the production of sugars from banana fruit peel waste. Biofuels 2018; 1-8.[http://dx.doi.org/10.1080/17597269.2018.1442665] ]. Around the globe, termites comprise one of the extremely successful groups of organisms that are capable of degrading wood and are therefore vital for expanding their roles both in carbon turnover and as possible source of biocatalysts for the efforts that are aimed at converting the wood biomass into biofuels. The main aim of the present study is to unravel the structural composition and functional attributions of different microbial communities from the termite degraded B. monosperma and to study the variation in microbial community structure with respect to changing geographical locations and environmental factors. For this, we used the Next Generation Sequencing technology targeting the hypervariable region of the microbial 16S rRNA genes to unravel the bacterial community residing in termite degraded Butea sample. The complete work flow of this work is presented in Fig. (1).
2. MATERIALS AND METHODS
2.1. Sample Collection and Preparation
A total of nine samples from termites degraded Butea monosperma tree were randomly sampled from Jhira Ghati, Forest Area, Sagar, M.P. India (Latitude 23°12.070’N and longitude 079°02.357E’) as shown in Fig. (2). All the collected samples were properly labelled and packed in plastic bags, and were immediately transported to the laboratory and stored at 4°C till further use.
2.2. DNA Extraction and 16S rRNA Amplicon Sequencing
In order to profile the full-scale taxonomic diversity, the total bacterial genomic DNA was extracted from 250 mg/sample using the Mo Bio Power soil DNA isolation kit (MoBio Laboratories, Inc. CA, USA) as per the manufacturer’s instructions. The purity and integrity of the extracted genomic DNA were assessed through agarose gel electrophoresis (1%) and the concentration was quantified via NanoDrop 2000 spectrophotometer (Thermo Scientific) by measuring the absorbance at 260/280 nm. Total genomic DNA extracted from three replicates was pooled to one sample in order to cover the maximum diversity. The DNA samples were stored at -20°C till sequencing. Illumina’s MiSeq system was used for the 16S rRNA amplicon sequencing of extracted total genomic DNA as this instrument is a dominant platform for amplicon sequencing in microbial ecology studies due to its great flexibility, high-throughput, fast-turnaround time, longer sequence reads and higher accuracy.
2.3. Taxonomic Profiling of Termite 16S rRNA Gene Using MG-RAST
In this study, we have analysed the 16S rRNA gene sequences using Metagenomic Rapid Annotation Using Subsystem Technology (MG-RAST) server which is a comprehensive metagenome analysis server and offers a number of advantages like open access, publicly available metagenomes as well as comparative analysis of different metagenomes. The raw and unassembled sequences generated from high-throughput sequencing via Illumina were uploaded to MG-RAST server v 4.0.3 (http://metagenomics.anl.gov/) for taxonomic annotation and profiling under MG-RAST ID number mgm4738027.3 [21Wilke A, Glass EM, Bischof J, Braithwaite D, Gerlach W, Harrison T. MG-RAST Manual for version 36, revision 3 2017; 130]. The low quality regions were trimmed using SolexaQA [22Tang W, Wilkening J, Desai N, Gerlach W, Wilke A, Meyer F. A scalable data analysis platform for metagenomics. Big Data, 2013 IEEE Int Conf IEEE 2013; 21-6.[http://dx.doi.org/10.1109/BigData.2013.6691723] , 23Cox MP, Peterson DA, Biggs PJ, Solexa QA. At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinformatics. BioMed Central 2010; 11: 485.[http://dx.doi.org/10.1186/1471-2105-11-485] [PMID: 20875133] ] following downstream analysis by using RDP as the reference database with the parameters set as default.
2.4. Comparative Taxonomic and Functional Diversity for Selected samples by QIIME
QIIME is an open-source software pipeline built by making use of PyCogent toolkit [24Knight R, Maxwell P, Birmingham A, et al. PyCogent: A toolkit for making sense from sequence. Genome Biol 2007; 8(8): R171.[http://dx.doi.org/10.1186/gb-2007-8-8-r171] [PMID: 17708774] ] and is accessible at http://www.qiime.org/ or http://github.com/qiime/. For comparative taxonomic profiling and annotation of the termite metagenome with other termite metagenomes, the OTU table of the four publically available 16S rRNA amplicon metagenomes was downloaded from MG-RAST server v 4.0.3 (http://metagenomics.anl.gov/) [21Wilke A, Glass EM, Bischof J, Braithwaite D, Gerlach W, Harrison T. MG-RAST Manual for version 36, revision 3 2017; 130] and analysis was carried out by QIIME release 1.9.1 using default parameters by Caporaso et al. [25Caporaso JG, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010; 7(5): 335-6.[http://dx.doi.org/10.1038/nmeth.f.303] [PMID: 20383131] ]. The details of the four different geographically located termite metagenomes used for the comparison purposes are listed in Table 1. To ascertain functional profile of the termite degraded Butea monosperma, the taxonomic abundance dataset was uploaded on the METAGEN assist server [26Arndt D, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, et al. METAGENassist: A comprehensive web server for comparative metagenomics.2012; 40: W88-95.] (http://www. metagen assist. ca/METAGENassist/) which performs taxonomy to phenotypic mapping.
3.1. Illumina Sequencing
There are various approaches to study the microbiota of the termite degraded Butea such as culture-independent and culture-depended. The culture dependent molecular approaches can trace very less amount of bacterial diversity [27Malla MA, Dubey A, Yadav S, Kumar A, Hashem A, Abd Allah EF. Understanding and designing the strategies for the microbe-mediated remediation of environmental contaminants using omics approaches. Front Microbiol 2018; 9: 1132.[http://dx.doi.org/10.3389/fmicb.2018.01132] [PMID: 29915565] , 28Ndoye B, Rasolofo EA, LaPointe G, Roy D. A review of the molecular approaches to investigate the diversity and activity of cheese microbiota. Dairy Sci Technol 2011; 495-524.[http://dx.doi.org/10.1007/s13594-011-0031-8] ] whilst the culture independent methods still have some limitations such as they are time consuming, costly, low coverage and insufficient to reflect the true bacterial diversity, however recent advances in molecular approach of high throughput next-generation sequencing technology have provided large scale analysis with unprecedented depths and coverage [29Malla MA, Dubey A, Kumar A, Yadav S, Hashem A, Abd Allah EF. Exploring the human microbiome: The potential future role of next-generation sequencing in disease diagnosis and treatment. Front Immunol 2019; 9: 2868.[http://dx.doi.org/10.3389/fimmu.2018.02868] [PMID: 30666248] , 30Hutchison CA. DNA sequencing: Bench to bedside and beyond. Nucleic Acids Res. Oxford University Press; 2007;35:6227–37.]. Nowadays, the use of NGS technology has made it possible to study the descriptive metagenomics. Metagenomics study can reveal the community structure and variation of the microbiome and the relative microbial abundance is estimated based on different physiological and environmental conditions [31Mandal RS, Saha S, Das S. Metagenomic surveys of gut microbiota. Genomics, Proteomics Bioinf; 2015;13:148–58.[http://dx.doi.org/10.1016/j.gpb.2015.02.005] , 32Yang H, Huang X, Fang S, Xin W, Huang L, Chen C. Uncovering the composition of microbial community structure and metagenomics among three gut locations in pigs with distinct fatness. Sci Rep 2016; 6: 27427.[http://dx.doi.org/10.1038/srep27427] [PMID: 27255518] ]. Ultimately, these studies can help to reflect connections between the identity of a microbe or a community and their respective functions in the environment [33Chistoserdova L. Is metagenomics resolving identification of functions in microbial communities? Microb Biotechnol 2014; 7(1): 1-4.[http://dx.doi.org/10.1111/1751-7915.12077] [PMID: 23945370] ].
The total genomic DNA extracted from the microbial communities of the termite degraded Butea was subjected to Illumina sequencing for the bacterial 16S rRNA gene region (v3 and v4). Tagmentation (fragmentation and adapter tagging), PCR amplification and clean-up of the genomic DNA were performed as per Illumina Nextera library preparation guidelines (Illumina, USA). 4nM libraries were prepared, pooled and sequenced on the Illumina Miseq desktop sequencer (2×300 run). The reads generated were then subjected to quality check and were subsequently analyzed using the RDP database classifier with parameters set at default values [34Meyer F, Paarmann D, D’Souza M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 2008; 9: 386.[http://dx.doi.org/10.1186/1471-2105-9-386] [PMID: 18803844] ].
3.2. Taxonomic Composition of the Bacterial Metagenomes from Termite Sample
A total of 168,310 sequences, with 86,275,309 base pairs and an average length of 513 bps were obtained following the sequencing through Illumina sequencer of the termite degraded Butea sample (TO) (Table 2). The taxonomic profiling of the DNA samples extracted from the termite degraded Butea was performed against RDP database with the default parameters and the phylotypes were analyzed upto the genus level. The MG-RAST server used to assess the quality of this termite metagenome confirmed that 96.35% of the sequences belong to the bacteria, with 3.11% being Eukaroyata, followed by Archaea 0.04% and the unclassified being 0.55% (Fig. 3). The phylum level composition of the bacterial communities of the termite degraded Butea sample was composed of 9 phyla with one being unidentified. At the phylum level, Actinobacteria was predominant with a relative abundance of (21.27%) followed by Proteobacteria (14.18%), Firmicutes (10.46%), Bacteroidetes (4.11%) and Gemmatimonadetes (3.23%). The rest of the phyla showed a relatively minimum abundance with Chloroflexi showing a percentage abundance of (1.57%) followed by Verrucomicrobia (0.92%), and Acidobacteria (0.51%) (Fig. 4a). Level 3rd (class) distribution within these phyla can be observed in Fig. 4b. At the genus level, Bacillus were predominant, with the relative abundance of 4.33% followed by Gemmatimonas (3.13%), Mycobacterium (1.82%), Verticillium (1.70%), Acidimicrobium (1.69%), Thermoleophilum (1.23%), Nocardioides (1.44%), Terrimonas and Acidithermus (1.09) and Clostridium (1.05%) (Fig. 4c).
3.3. Rarefaction Analysis
Rarefaction curve is an index used to estimate the sampling depth. The adequacy of the sampling was estimated by means of rarefaction curve. The sampling efforts and the species richness of the termite degraded Butea sample were estimated by means of the rarefaction analysis. The current results depict that the rarefaction curves for this sample reached saturation, which corresponded to the true bacterial diversity within the sample (Fig. 5).
3.4. Comparative Taxonomic Profiling of Termite Degraded Butea with Other Publically Available Termite Metagenome Samples
In order to ascertain the variation in microbial diversity with respect to different geographical location and environmental conditions, the termite degraded Butea sample (TO) was then compared with four other publically available termite degraded plants metagenomes. The four different termite degraded plant metagenomes used for comparison purposes were (TA), (TFL), (TC) and (TFA) from Australia, la Coubre France, Soroa Cuba and Aledo lA France respectively, with each belonging to different geographical location and having different weather conditions. To our surprise, a varying trend was observed at the phylum level with Actinobacteria dominating in the present study (TO) while Bacteroidetes dominating in TA and Firmicutes in case of TC. Notably, Proteobacteria is the second most abundant taxa in three samples (TO, TA, TC) whereas the same (Proteobacteria) is the most dominant in case of TFL and TFA (Fig. 6). Comparative taxonomic profiles of the bacterial metagenomes from termite degraded Butea (TO) and other termite samples (TA, TC, TFL and TFA) at the phylum level using MG-RAST based-RDP database indicate a significant variation in the relative abundance for bacterial composition between different termite samples (Fig. 7a). The PCA plot generated illustrated the clustering of the metagenomes as individual samples, thereby indicting a significant variation in microbial composition within these termite metagenomes. The PCA plot of the 5 termite metagenomes using the taxonomic classification at phylum level showed that the termite metagenome (TO) did not cluster with the termite samples (TA, TC, TFL and TFA), which more or less cluster with each other. The present metagenomics study of the termite had significantly higher Actinobacterial but a much lesser Bacteroidetal reads compared to TA (Fig. 7b). Nonetheless, the termite metagenome from
Table 1 Selected 16S rRNA amplicons/metagenome that were used in the present study for comparative taxonomic profiling and annotation.
Table 2 Details of the sequence analysis of Termite degraded Butea sample (TO).
Schematic representation of whole process.
Location of sampling site using Google earth image of Jhira Ghati, Forest Area, Sagar, M.P. India.
Sequence distribution of the Termite degraded B. monosperma 16S rRNA amplicons at domain level.
Taxonomic profiles of the bacterial communities of Termite degraded B. monosperma (TO) (a) phylum level (b) class level and (c) Genus level.
Rarefaction curves showing the saturation level of the Termite degraded B. monosperma. The vertical axis corresponds to the species count and the horizontal to number of reads.
Comparative taxonomic profiles of bacterial communities of (TO), TA), TFL), TFA) and Cuba (TC) based on the MG-RAST RDP database.
3.5. Taxonomic to Phenotypic Mapping of the Termite Degrade B. Monosperma Metagenome
The extensive in silico analysis of the termite degraded B. monosperma metagenomic data by METAGENassist, a web based server, provides taxonomic phenotypic mappings of the datasets. The functional annotation was carried out by taking into consideration the metabolic compositions (Fig. 8). The functional mapping of the sample illustrated that termite gut can be viewed as a micro-ecological habitat, thereby serving as an environment for a number of bio-geochemical cycles. The representation of the functional abundance in the termite degraded B. monosperma showed a significance of the processes such as ammonia oxidation, chitin degradation, xylan degradation, nitrate reduction, nitrogen fixation, dehalogenation and sulfate reduction. The occurrence of these processes gives the indications of the operation of the sulphur and nitrogen cycles with in this micro-environment. In the present sample, the relative abundance of the different groups is as ammonia oxidizers (30.5%), sulfate reducers (23.5%), dehalogenators (20.1%), Nitrate reducers (15.6%), sulfide oxidizers (14.8%), xylan degraders (14.8%), nitrogen fixers (8.3%), and chitin degraders (7.4%), with a lower percentage abundance of sulphur metabolizing, lignin degraders and pollutant degraders.
Metagenomic studies provide a deeper knowledge of the complex microbial consortia inhabiting any environment [35Dubey A, Malla MA, Khan F, et al. Soil microbiome: A key player for conservation of soil health under changing climate. Biodivers Conserv. Springer 2019; 1-25.]. Assessment of the bacterial diversity using amplicon sequencing relies on the hyper-variable regions of 16S rRNA. The aim of the present study was to ascertain taxonomic composition of the termite degraded Butea using 16S rRNA amplicon sequencing. In order to obtain a meaningful data from termite sequencing, the data obtained from the present study was compared with the termite sequences that were available publically, so as to get an insight of the termite microbiome.
Fig. (7) (a) Heat map indicating a significant variation in the relative abundance of different bacterial communities of (TO), (TA), (TFL), (TFA), (TC) (b) PCA plot of (TO), (TA), (TFL), (TFA) and (TC).
Functional profiles of the bacterial metagenomes from Termite degraded B. monosperma (TO) from Sagar, (M. P.) India on the basis of metabolism.
The study of normal gut microbiota is essential as species composition and proportion changes with respect to varying environmental conditions. The extensive use of high throughput next generation sequencing helps to unravel the impact of external factors. Till now, there is very less knowledge regarding the biodiversity analysis of termite degraded Butea using HT-NGS. In our study, we used the 16S rDNA amplicon study of the termite degraded Butea microbiota and compared our results with the publically available termite degraded plants to characterize the bacterial diversity in termite gut and to investigate the diversity and species abundance variation with changing environmental conditions and geographical locations. As the gut microbiota can be influenced by a number of factors including diet, environment and locations, our data presented in this study should be considered as location-specific.
The most of the predominant genera found in our study are common residents of termite gut [36Rossmassler K, Dietrich C, Thompson C, et al. Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites. Microbiome 2015; 3: 56.[http://dx.doi.org/10.1186/s40168-015-0118-1] [PMID: 26607965] ]. Taxonomic profiling and categorization of the bacterial communities predicted that termite gut environment harbors a diverse array of bacterial communities that are unique to those of other environments [37Makonde HM, Mwirichia R, Osiemo Z, Boga HI, Klenk HP. 454 Pyrosequencing-based assessment of bacterial diversity and community structure in termite guts, mounds and surrounding soils. Springerplus 2015; 4: 471.[http://dx.doi.org/10.1186/s40064-015-1262-6] [PMID: 26355944] ]. In the present study, Actinobacteria dominated (21.27) out of the classified sequences (≥62%) followed by Protobacteria, Firmicutes, Gemmatimonadetes, Bacteroidetes, Chloroflexi, Verrucomicrobia, Planctomycetes and Cyanobacteria. Similarly, a different trend was observed in the bacterial diversity of other samples (Fig. 6). The variation in bacterial diversity may be associated with the change in environmental conditions (like pH, temperature, oxygen content, altitude and feeding habit) because the colonization of the microbes relies on the environmental factors [38Kumar A, Chanderman A, Makolomakwa M, Perumal K, Singh S. Microbial production of phytases for combating environmental phosphate pollution and other diverse applications. Crit Rev Environ Sci Technol 2016; 46(6): 556-91.[http://dx.doi.org/10.1080/10643389.2015.1131562] , 39Engel P, Moran NA. The gut microbiota of insects - diversity in structure and function. FEMS Microbiol Rev 2013; 37(5): 699-735.[http://dx.doi.org/10.1111/1574-6976.12025] [PMID: 23692388] ]. Protobacteria are usually characterized by their dominance in early stages of composting of plant materials and hence are responsible for cellulose degradation [40Singh NB, Kumar A, Rai S. Potential production of bioenergy from biomass in an Indian perspective. Renew Sustain Energy Rev 2014; 65-78.[http://dx.doi.org/10.1016/j.rser.2014.07.110] , 41Yu H, Zeng G, Huang H, et al. Microbial community succession and lignocellulose degradation during agricultural waste composting. Biodegradation 2007; 18(6): 793-802.[http://dx.doi.org/10.1007/s10532-007-9108-8] [PMID: 17308882] ]. Thus, the colonization of Protobacteria can be endorsed with their ability to hydrolyze the hemi-cellulose. The abundance of bacterial phyla Actinobacteria (TO), Bacteroidetes (TA), Spirochaetes (TFL), Proteobacteria (TFA) and Firmicutes (TC) has shown their contribution to symbiotic relation of the termite species, acetogenesis (acetate production) and nitrogen acquisition, and cellulose degradation [38Kumar A, Chanderman A, Makolomakwa M, Perumal K, Singh S. Microbial production of phytases for combating environmental phosphate pollution and other diverse applications. Crit Rev Environ Sci Technol 2016; 46(6): 556-91.[http://dx.doi.org/10.1080/10643389.2015.1131562] , 42Fall S, Hamelin J, Ndiaye F, et al. Differences between bacterial communities in the gut of a soil-feeding termite (Cubitermes niokoloensis) and its mounds. Appl Environ Microbiol 2007; 73(16): 5199-208.[http://dx.doi.org/10.1128/AEM.02616-06] [PMID: 17574999] ]. Termite microbiome incorporates abundant and morphologically distinct communities of Spirochetes, which own metabolic properties such as N2 fixation and H2/CO2-acetogenesis, thereby contributing to the nitrogen, carbon as well as energy requirements for the termite hosts. The presence of Spirochaetes, although relative lesser in abundance in the termite degraded Butea metagenome signifies that the plant is a good source of carbon and nitrogen and also reflects diverse degradation and bioconversion mechanisms and pathways. The presence of Acidobacteria (although relatively much lesser in abundance) in our sample (TO) and its absence in other samples hints towards the acidic nature of the environment. In the same way, their PCA is also differentiated by similar spread pattern. The STAMP generated PCA plot illustrated the clustering of the termite degraded Butea and other termite metagenomes as unique samples. Analysis revealed that the samples entirely had different bacterial compositions which may be due to varied dietary habitats as well as different geographical locations and environmental conditions (pH, carbon, nitrogen, oxygen, and temperature) [43Raghunandan K, Mchunu S, Kumar A, et al. Biodegradation of glycerol using bacterial isolates from soil under aerobic conditions J Environ Sci Heal - Part A 2014; 49: 85-92.[http://dx.doi.org/10.1080/10934529.2013.824733.] ]. This seems to be a foreseeable result as other studies hold the same view [44Lauber CL, Strickland MS, Bradford MA, Fierer N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol Biochem 2008; 40(9): 2407-15.[http://dx.doi.org/10.1016/j.soilbio.2008.05.021] ]. As environmental conditions can directly alter the physiological status of the endemic microbial communities [45Haruta S, Kanno N. Survivability of Microbes in Natural Environments and Their Ecological Impacts. Microbes Environ Environ 2015; 30(2): 123-5.[http://dx.doi.org/10.1264/jsme2.ME3002rh] , 46Abreu NA, Taga ME. Decoding molecular interactions in microbial communities. FEMS Microbiol Rev 2016; 40(5): 648-63.[http://dx.doi.org/10.1093/femsre/fuw019] [PMID: 27417261] ]. Here, we conclude that the shift within the microbial diversity to a large extent depends on environmental conditions.
After sequencing on Illumina platform and comparative analysis of termite degraded Butea collected from forest area of Sagar, M.P. India (TO), the results obtained indicate that this study is ideal for the taxonomic profiling of different bacterial communities responsible for the degradation of Butea. Taxonomic profiling of the termite metagenome showed that Actinobacteria was dominant followed by Proteobacteria, Firmicutes, Gemmatimonadetes, Bacteroidetes. Preponderance of the Protobacteria suggests that termite is the chief degrader of lignocellulose on the Butea. This study is the first of its kind to explore the taxonomic diversity of the termite responsible for the degradation of the Butea. The microbial community is influenced by various environmental factors (moisture and increase in temperature) within the degraded pile of Butea. As the location and environmental factors change, so does the bacterial diversity. A number of the sequences from the sample are only classified upto the domain level or either the phylum level, signifying the presence of these uncultivable or novel bacterial species. The isolation of the bacteria from the termite degraded Butea monosperma with diverse functions such as nitrogen cycling and dinitrogen fixation, production of hydrolytic enzymes, and pollutant degradation may prove to be of great importance to industries and scientific communities that are focused towards the degradation of ligno-cellulosic biomass for biofuel production. To develop approaches in order to cultivate these unknown and novel bacteria sets a big stage for utilizing these novel organisms for various biotechnological applications.
AUTHOR'S CONTRIBUTION STATEMENT
AK executed the experimental work, designed and wrote the manuscript. MM and AD analysed and arranged the data, PV gave some suggestions.
ETHICAL APPROVAL AND CONSENT TO PARTICIPATE
HUMAN AND ANIMAL RIGHTS
No animals/humans were used for studies that are the basis of this research.
CONSENT FOR PUBLICATION
AVAILABILITY OF DATA AND MATERIALS
UGC-BSR Start-up grant by University Grant Commission, New Delhi, India.
CONFLICT OF INTEREST
The authors declare no conflict of interest, financial or otherwise.
AK would like to acknowledge UGC-BSR Start-up grant by University Grant Commission, New Delhi, India for providing financial support for this study. AD would like to acknowledge the Department of Science and Technology, New Delhi, India for providing financial support through the DST Inspire Ph.D. Fellowship (IF160797).
Patil MV, Pawar S, Patil DA. Ethnobotany of Butea monosperma (Lam.) Kuntze in North Maharashtra, India. Medicine (Baltimore) 2006; 5: 323-5.
Chopra RN. Indigenous Drugs of India Their Medical and Economic Aspects 20 British Indian Street 1933.
Chandrasekharaiah M, Thulasi A, Bagath M, et al. Molecular cloning, expression and characterization of a novel feruloyl esterase enzyme from the symbionts of termite(Coptotermes formosanus) gut. BMB Rep 2011; 44(1): 52-7.[http://dx.doi.org/10.5483/BMBRep.2011.44.1.52] [PMID: 21266107]
Zhang D, Lax AR, Bland JM, Allen AB. Characterization of a new endogenous endo-β-1,4-glucanase of Formosan subterranean termite (Coptotermes formosanus). Insect Biochem Mol Biol 2011; 41(4): 211-8.[http://dx.doi.org/10.1016/j.ibmb.2010.12.006] [PMID: 21195179]
Singh JK, Vyas P, Dubey A, et al. Assessment of different pretreatment technologies for efficient bioconversion of lignocellulose to ethanol. Frontiers in bioscience (Scholar edition) 2018; 10: 350-71.[http://dx.doi.org/10.2741/s521]
Arndt D, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, et al. METAGENassist: A comprehensive web server for comparative metagenomics.2012; 40: W88-95.
Malla MA, Dubey A, Yadav S, Kumar A, Hashem A, Abd Allah EF. Understanding and designing the strategies for the microbe-mediated remediation of environmental contaminants using omics approaches. Front Microbiol 2018; 9: 1132.[http://dx.doi.org/10.3389/fmicb.2018.01132] [PMID: 29915565]
Ndoye B, Rasolofo EA, LaPointe G, Roy D. A review of the molecular approaches to investigate the diversity and activity of cheese microbiota. Dairy Sci Technol 2011; 495-524.[http://dx.doi.org/10.1007/s13594-011-0031-8]
Malla MA, Dubey A, Kumar A, Yadav S, Hashem A, Abd Allah EF. Exploring the human microbiome: The potential future role of next-generation sequencing in disease diagnosis and treatment. Front Immunol 2019; 9: 2868.[http://dx.doi.org/10.3389/fimmu.2018.02868] [PMID: 30666248]
Hutchison CA. DNA sequencing: Bench to bedside and beyond. Nucleic Acids Res. Oxford University Press; 2007;35:6227–37.
Yang H, Huang X, Fang S, Xin W, Huang L, Chen C. Uncovering the composition of microbial community structure and metagenomics among three gut locations in pigs with distinct fatness. Sci Rep 2016; 6: 27427.[http://dx.doi.org/10.1038/srep27427] [PMID: 27255518]
Meyer F, Paarmann D, D’Souza M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 2008; 9: 386.[http://dx.doi.org/10.1186/1471-2105-9-386] [PMID: 18803844]
Dubey A, Malla MA, Khan F, et al. Soil microbiome: A key player for conservation of soil health under changing climate. Biodivers Conserv. Springer 2019; 1-25.
Makonde HM, Mwirichia R, Osiemo Z, Boga HI, Klenk HP. 454 Pyrosequencing-based assessment of bacterial diversity and community structure in termite guts, mounds and surrounding soils. Springerplus 2015; 4: 471.[http://dx.doi.org/10.1186/s40064-015-1262-6] [PMID: 26355944]
Kumar A, Chanderman A, Makolomakwa M, Perumal K, Singh S. Microbial production of phytases for combating environmental phosphate pollution and other diverse applications. Crit Rev Environ Sci Technol 2016; 46(6): 556-91.[http://dx.doi.org/10.1080/10643389.2015.1131562]
Fall S, Hamelin J, Ndiaye F, et al. Differences between bacterial communities in the gut of a soil-feeding termite (Cubitermes niokoloensis) and its mounds. Appl Environ Microbiol 2007; 73(16): 5199-208.[http://dx.doi.org/10.1128/AEM.02616-06] [PMID: 17574999]
Lauber CL, Strickland MS, Bradford MA, Fierer N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol Biochem 2008; 40(9): 2407-15.[http://dx.doi.org/10.1016/j.soilbio.2008.05.021]
Waidele L, Korb J, Voolstra CR, Künzel S, Dedeine F, Staubach F. Differential ecological specificity of protist and bacterial microbiomes across a set of termite species. Front Microbiol 2017; 8: 1-13.