ORIGINAL RESEARCH article
Metatranscriptome Analysis Deciphers Multifunctional Genes and Enzymes Linked With the Degradation of Aromatic Compounds and Pesticides in the Wheat Rhizosphere
- 1ICAR-National Bureau of Agriculturally Important Microorganisms, Maunath Bhanjan, India
- 2Department of Bio-Medical Engineering and Bio-Informatics, Chhattisgarh Swami Vivekanand Technical University, Bhilai, India
- 3ERA Chair of Green Chemistry, Department of Chemistry and Biotechnology, School of Science, Tallinn University of Technology, Tallinn, Estonia
Agricultural soils are becoming contaminated with synthetic chemicals like polyaromatic compounds, petroleum hydrocarbons, polychlorinated biphenyls (PCBs), phenols, herbicides, insecticides and fungicides due to excessive dependency of crop production systems on the chemical inputs. Microbial degradation of organic pollutants in the agricultural soils is a continuous process due to the metabolic multifunctionalities and enzymatic capabilities of the soil associated communities. The plant rhizosphere with its complex microbial inhabitants and their multiple functions, is amongst the most live and dynamic component of agricultural soils. We analyzed the metatranscriptome data of 20 wheat rhizosphere samples to decipher the taxonomic microbial communities and their multifunctionalities linked with the degradation of organic soil contaminants. The analysis revealed a total of 21 different metabolic pathways for the degradation of aromatic compounds and 06 for the xenobiotics degradation. Taxonomic annotation of wheat rhizosphere revealed bacteria, especially the Proteobacteria, actinobacteria, firmicutes, bacteroidetes, and cyanobacteria, which are shown to be linked with the degradation of aromatic compounds as the dominant communities. Abundance of the transcripts related to the degradation of aromatic amin compounds, carbazoles, benzoates, naphthalene, ketoadipate pathway, phenols, biphenyls and xenobiotics indicated abundant degradation capabilities in the soils. The results highlighted a potentially dominant role of crop rhizosphere associated microbial communities in the remediation of contaminant aromatic compounds.
Toxic organic compounds such as polychlorinated biphenyls (PCBs), chlorophenols, polycyclic aromatic hydrocarbons (PAHs), petroleum hydrocarbons, dyes and pesticides are becoming major contaminants of agricultural soils (Diez, 2010; Valentìn et al., 2013). Polyaromatic, chlorinated and nitrobenzene compounds, dioxins, PCBs and the pesticides represent principal xenobiotics (Sinha et al., 2009; Agrawal and Dixit, 2015) that cause biospheric pollution (Rabinovich et al., 2004). PAHs and PCBs are resistant organic contaminants against biodegradation by native microflora (Chikere, 2013; Bisht et al., 2015). Crop roots and the rhizosphere are directly under the influence of hardly biodegradable organic pollutants and contaminants, especially the PAHs and pesticides under any contaminated soil environment (Chen et al., 2013; Lawal, 2017). This may create serious and threatening concerns for the healthy soil and crop environment and therefore, the biodegradation and/or rhizoremediation of the chemical wastes and organic pollutants are now attracting central focus worldwide (Chikere, 2013; Lu et al., 2014; Kronenberg et al., 2017).
The chemicals when used continuously, increasingly and excessively in the cropping system become contaminants raising the problems of soil and rhizosphere pollution (Mwangi et al., 2010; Chauhan and Singh, 2015). Abundance of organic contaminants and pesticides in the food chain and water is reported to cause carcinogenesis, neurotoxicity and reproduction and cell development disorders (Burrows et al., 2002; Damalas and Eleftherohorinos, 2011; Prüss-Ustün et al., 2011; Myers et al., 2016). Therefore, the process of biodegradation or biologically catalyzed transformation of organic contaminants is of prime importance (Karigar and Rao, 2011; Joutey et al., 2013; Kehinde and Isaac, 2016). Although abiotic pathways are known for the breakdown of organic pollutants and pesticides in the soils and water, microbe-mediated degradation is the primary mechanism of contaminant remediation and detoxification (Singh and Walker, 2006; Akbar and Sultan, 2016; Javaid et al., 2016). The role of microbial communities is of immense importance due to their multifunctional degradation capabilities (Das and Chandran, 2011; Dechesne et al., 2014). Microbial (bacterial and fungal) degradation of pesticide, xenobiotic compounds and organic chemicals has been widely reported (McGuinness and Dowling, 2009; Porto et al., 2011) since long, although the direct, indirect or intermediary linking of the metabolic pathways connected with the removal of contaminants is less understood (Ladino-Orjuela et al., 2016).
Rhizosphere represents the plant-root interface which is populated by distinctive microbial communities having a functional influence on the plants (McNear, 2013). In the rhizosphere, potential communities accelerate biodegradation processes and help plants enhance co-metabolism to degrade organic contaminants by (i) facilitating selective enrichment of biodegrader microorganisms for degradation of xenobiotics in the root free soils (Nichols et al., 1997) (ii) enhancing metabolism of microbial growth by secreting natural substrate depending on the quantity of xenobiotics (Haby and Crowley, 1996; Alexander, 1999) and (iii) enriching natural compounds that can provoke co-metabolism of xenobiotics in specific microbes exhibiting genes or plasmids with degradation functions (Gupta et al., 2016). In any way, crop rhizosphere resists or tolerates to the level of organic contaminants, xenobiotics and pesticides, majorly with the help of microbial functions linked with the metabolic degradation pathways that partially or totally detoxify or minimize quality and quantity of contaminants in the soils. Various studies mentioned microbial biodegradation of different PCBs and xenobiotics through plant roots (Pilon-Smits, 2005; Balcom et al., 2016) along with genes involved in degradation (Wittich, 1998; Arai et al., 2000; Duffner et al., 2000; Esteve-Núñez et al., 2001; Furukawa et al., 2004; Galvão et al., 2005) and transcriptional regulators responsible for degradation of aromatic compounds (Uchiyama and Miyazaki, 2013).
Metagenomics and metatranscriptomics approaches that allowed exploration of taxonomic communities and help linking of associated functions without cultivation of microorganisms (Schloss and Handelsman, 2003) specifically facilitated overall examination of microbial community-linked metabolic pathways related to specific metabolic functions. The catabolic genes related to the degradation of xenobiotics can be characteristically annotated and linked to the identified taxonomic group/genus of the microbial communities (Chikere, 2013; Sharma and Sharma, 2018). We have annotated and deciphered twenty wheat rhizosphere metatranscriptomes available from EBI Metagenomics database for digging out transcripts/genes related to direct degradation or intermediary degradation functions of aromatic compounds and xenobiotics and reported prospective rhizoremediation capacity linked with the rhizobiome of the crop plants.
Materials and Methods
Twenty wheat rhizosphere metatranscriptome samples (ERS025421, ERS238568, ERS238569, ERS238570, ERS238571, ERS238572, ERS238573, ERS238574, ERS238575, ERS238576, ERS238577, ERS238578, ERS238579, ERS238580, ERS238581, ERS238582, ERS238583, ERS238584, ERS238585, and ERS238586) were downloaded from the EBI Metagenomics database1. Raw sequence reads were taken for each dataset and processed further for the annotation and analysis.
Quality Check and Metatranscriptome Sample Processing
Raw sequence reads of each sample were processed separately and were uploaded to the MG-RAST (Metagenome Rapid Annotation using Subsystem Technology) server2 (Meyer et al., 2008). Further processing was carried out through the standard pipeline keeping default parameters. For identification and annotation of proteins and other sequences, sequence similarity searches were done against different databases associated with MG-RAST pipeline. Genbank (Benson et al., 2013, 2015), protein databases M5NR (Wilke et al., 2012), SEED (Chen et al., 2013), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2008) were used for annotations at various levels. The data analysis was performed using High Performance Supercomputing Facility of ICAR-National Bureau of Agriculturally Important Microorganisms, India.
Taxonomic and Functional Annotation
For identification of taxonomic abundance, all the datasets were individually aligned against M5NR database (Wilke et al., 2012) and Genbank (Benson et al., 2013, 2015) through the MG-RAST pipeline (Meyer et al., 2008). Parameters taken include maximum E-value: 1 × 10−5, minimum percentage identity: 60%, and minimum alignment length: 15. For the purpose of functional annotation, datasets were aligned against SEED subsystems taking into account similar parameters (Overbeek et al., 2014; Imchen et al., 2017).
Metabolic Pathway Analysis
For identification and functional characterization of diverse metabolic pathways related to the metabolism of aromatic compounds and xenobiotics biodegradation and metabolism, SEED Subsystems was applied from the MG-RAST pipeline (parameters were minimum alignment length 15 and E-value cutoff 1e-5). For obtaining entire metabolic pathways information, global gene expression was annotated and further visualized with KEGGmapper (an in-built tool on the basis of the KEGG pathway mapping system). Level 2 and Level 3 SEED subsystem of MG-RAST was applied for annotation. Level 3 of SEED subsystem corresponds to KEGG pathway3. Furthermore, the Enzyme Commission (EC) number was identified for the sequences of corresponding enzymes. Pathways, transcripts and enzymes for the biological processes were identified in each dataset and corresponding abundance of transcripts/enzymes was calculated.
Availability of Data and Associated Information
The metatranscriptome dataset are publically available from the EBI Metagenomics database4.
Results and Discussion
Organic chemical contaminants, toxic wastes and synthetic pesticides are reported to damage the quality of agricultural soils and disturb the intrinsic balance between biological and ecological functional services (Valentìn et al., 2013). Microbial communities with their efficient metabolic pathways and enzymatic capabilities help breakdown of contaminating organic chemicals and pesticides through biotransformation and biodegradation (Parte et al., 2017). The functions of microorganisms pertaining to biotransformation and biodegradation help improve soil quality under the conditions of xenobiotic contamination (Ortiz-Hernández et al., 2013). The metatranscriptome analysis reports such widely abundant microbial communities, which are known for their role of degradation and remediation of organic pollutants and prominently abundant transcripts, which were directly related to the enzymes that have potential to degrade various organic chemicals in 20 wheat rhizosphere metatranscriptomes.
Assessment of Taxonomic Abundance
A total of 688.8 Mb cDNA datasets varying in size from 2.8 to 102.3 Mb and reads from 6,832 to 249,64720 of 20 wheat rhizosphere metatranscriptomes were analyzed. Bacterial communities dominated taxonomic profiling in all the samples followed by eukaryota. Within the identified bacterial phylum, Proteobacteria and Actinobacteria remained the most dominant phyla. Firmicutes, Bacteroidetes, Streptophyta, Cyanobacteria, Acidobacteria, Verrucomicrobia, Planctomycetes, Ascomycota and unclassified microbes represented major phyla in wheat rhizosphere (Figure 1). Microbial degradation is reported to be the major process by which aromatic compounds are biodegraded/biotransformed into less toxic, persistent and complex metabolites (Haritash and Kaushik, 2009). The abundance of Proteobacteria, Actinobacteria, Acidobacteria, Bacteroidetes, Planctomycetes, and Firmicutes have earlier been reported from the soils (Janssen, 2006; Shade et al., 2012) and rhizosphere of many plants like Arabidopsis spp. (Lundberg et al., 2012; Bulgarelli et al., 2013; Schlaeppi et al., 2014), oak (Uroz et al., 2010), potato (Inceoglu et al., 2011), cucumber (Tian and Gao, 2014), soybean (Mendes et al., 2014), rice (Edwards et al., 2015), wheat (Ai et al., 2015; Donn et al., 2015), maize (Garcia-Salamanca et al., 2013; Peiffer et al., 2013) and others (Pascual et al., 2016).
Proteobacteria, being extremely diverse phylum comprising chemolithoautotrophs, heterotrophs and mixotrophs (Chao, 1984; Bent and Forney, 2008) have been extensively reported from different ecological zones in which they manage and regulate biogeochemical cycles (Ivarsoon and Holm, 2008). These communities are common inhabitant of soils to help functions for carbon, nitrogen and sulfur cycling (Spain et al., 2009; Castañeda and Barbosa, 2017; Yang et al., 2017). Their dominant abundance of the Proteobacteria as the largest phylum in the rhizosphere may possibly be attributed to their rapid growth rate (Johnston-Monje et al., 2016) and may have functional significance for the wheat in maintaining nutrient needs of the crop. The phyla Proteobacteria remained dominated by alpha-, beta- and gamma-proteobacterial communities that occupy diazotrophic status and possess symbiotic ability, root association and endophytic colonization with the plants (Pisa et al., 2011). Alpha-proteobacteria including the genera Sphingomonas and Azospirillum and gamma-proteobacteria, Xanthomonas remained dominant communities for the degradation of PAHs (Vinas et al., 2005; Simarro et al., 2013; Bisht et al., 2015). The study thus, indicated the dominance of diverse microbial communities, which were known for improved rhizoremediation against chemical pollutants in the soils contaminated with PAH, and other compounds (Kuiper et al., 2004).
The other dominant community of the wheat rhizosphere was the Actinobacteria, which are involved in the carbon cycle and secondary metabolite production (Jenkins et al., 2010; Castañeda and Barbosa, 2017). The phyla include different species with significant ecological importance and roles in decomposition, humus formation, nitrogen fixation and antibiotic production (Bond et al., 2004; Barka et al., 2016; Kielak et al., 2016). These communities are known biodegraders of PAHs (Leticia, 2006; Ghosal et al., 2016), pesticides (De Schrijver and De Mot, 1999; Fuentes et al., 2010), heavy metals (Alvarez et al., 2017), and xenobiotics (Polti et al., 2014). Firmicutes that also represented dominant community in the wheat rhizosphere inhabited varied range of environment as anaerobes (e.g., Clostridia) and aerobes (e.g., Bacillus). The phylum Bacteroidetes present in the wheat rhizosphere comprised of aerobes, anaerobes and facultative anaerobes (Janssen, 2006) that hold multifunctional potential in the soils (Wang et al., 2016). We also identified sequences linked with the cyanobacteria, the photosynthetic autotrophs commonly inhabiting agricultural soils (Oliveira et al., 2017). Cyanobacteria have also been predominantly implicated in the biodegradation of pesticides (Subramanian et al., 1994; Ibrahim et al., 2014), textile dyes (Dellamatrice et al., 2017) and aromatic compounds (Semple et al., 1999).
Aromatic Metabolism Pathway
Analysis of wheat rhizosphere metatranscriptomes for the pathways and transcripts related to the metabolism and degradation of aromatic compounds led to the identification of multiple enzymes associated with 21 different pathways (Table 1). Pathways identified for the metabolism of aromatic compounds were aromatic amin catabolism, Benzoate transport and degradation cluster, Gentisare degradation, Phenylacetyl-CoA catabolic pathway (core), Carbazol degradation cluster, Anaerobic benzoate metabolism, 4-Hydroxyphenylacetic acid catabolic pathway, Catechol branch of beta-ketoadipate pathway, Central meta-cleavage pathway of aromatic compound degradation, Homogentisate pathway of aromatic compound degradation, N-heterocyclic aromatic compound degradation, Protocatechuate branch of beta-ketoadipate pathway, Salicylate and gentisate catabolism, Benzoate catabolism, Benzoate degradation, Biphenyl Degradation, Naphtalene and anthracene degradation, Phenol hydroxylase and Phenylpropanoid compounds, Quinate and n-Phenylalkanoic acid degradation (Figure 2). Overall, 118 transcripts of 47 different enzymes linked with 21 pathways of aromatic compound metabolism were identified. The enzyme 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate hydrolase is commonly present in different pathways of aromatic metabolism namely carbazole degradation cluster, central meta-cleavage pathway of aromatic compound and biphenyl degradation. Members of this enzyme family unusually catalyze hydrolysis of C-C bonding in the biphenyl degradation pathway (Horsman et al., 2006). There were abundant (8) transcripts linked with the enzyme isoquinoline 1-oxidoreductase alpha subunit (EC 184.108.40.206) in the dataset for N-heterocyclic aromatic compound degradation pathway. The enzyme is specific for the degradation of N-heterocyclic compounds like isoquinoline and quinazoline and was reported from alphaproteobacteria Brevundimonas diminuta.5
TABLE 1. Identified aromatic degradation pathways and genes/enzymes involved in the metabolism of aromatic compounds in wheat rhizosphere metatranscriptome samples.
FIGURE 2. Abundance of pathways related to the metabolism of aromatic compounds in wheat metatranscriptomes. Name of the samples in which genes/enzymes were identified are written over the abundance of each pathway.
Various other pathways for aromatic compound degradation were also identified in the dataset. We identified transcripts for central meta-cleavage pathway, using which many organisms degrade high quantity of phenol (Mahiudddin et al., 2012). Besides, transcripts linked with homogentisate pathway and N-heterocyclic aromatic compound degradation was also identified. 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate hydrolase was identified for the central meta-cleavage pathway of aromatic compound degradation. Similarly, transcript for the isoquinoline 1-oxidoreductase alpha subunit (EC 220.127.116.11) was detected in the dataset for N-heterocyclic aromatic compound degradation (Table 1). For naphtalene and anthracene degradation, the transcripts related to the enzyme 2-hydroxychromene-2-carboxylate isomerase was also detected. The enzyme 2-hydroxychromene-2-carboxylate isomerase is known to degrade naphthalane sulfonates and anthracene and has been reported from bacterial species Pseudomonas testosteroni A3 and Sphingomonas yanoikuyae B1 (Kuhm et al., 1993; Kim et al., 2006).
Aromatic Amin Catabolism
The transcripts related to two prominent enzymes namely 4-hydroxyphenylacetate 3-monooxygenase (EC 18.104.22.168) and monoamine oxidase (22.214.171.124) were detected in wheat rhizosphere metatranscriptomes (Table 1). The enzyme 4-hydroxyphenylacetate 3-monooxygenase is involved in the breakdown of the aromatic ester 4-hydroxyphenylacetate, a ubiquitous compound that acts as the intermediate of aromatic amino acids, lignin and styrene (Mohamed et al., 2002). Members of this enzyme were reported from Escherichia coli (Xun and Sandvik, 2000). Monoamine oxidase family of flavoproteins catalyzes oxidative degradation of primary and secondary amines including polyamines, amino acids and tyramine (Díaz et al., 2001; Gaweska and Fitzpatrick, 2011). The cursing ecological role of polyclic and heterocyclic aromatic amines can be understood from the fact that this class of chemicals are carcinogens and environmental toxicants (Turesky and Le Marchand, 2011). Aromatic amines are also the component of polymers, dyes and other industrial products. Pesticide-derived aromatic amines commonly enter into the soil ecosystem following chemical control of insects and weeds in corn, soybean and other crops (Koutros et al., 2009; Luciana et al., 2015). Therefore, the presence of transcripts linked with the enzymes for the catabolism of aromatic amines signifies their role in the probable rhizoremediation of wheat rhizosphere.
Benzoate Degradation and Catabolism
Benzoate is a model compound to study anaerobic catabolism of aromatic compounds (Barragán et al., 2012). Benzoate is converted to benzoyl-Co-A following aerobic and anaerobic degradation pathways (Shinoda et al., 2004). Bacterial degradation of benzoate in the environment reflects biological strategy of microbial communities to increase cell fitness against these contaminants (Valderrama et al., 2012). Amongst the aromatic compounds and xenobiotic degradation linked transcripts found in this analysis, those related to the enzymes for degradation of benzoate were the most abundant in the wheat rhizosphere metatranscriptomes (Figure 2). Different enzyme families like decarboxylase, ligase, oxygenase, hydratase, kinase, isomerase, reductase, inner membrane translocater and alpha and beta subunits of dioxygenase were identified for the benzoate metabolism pathways including benzoate transport, degradation cluster and catabolism (Figure 2 and Table 1). Linked to the anaerobic benzoate degradation, we identified 15 sequences related to different enzymes including Acetyl-CoA acetyltransferase (EC 126.96.36.199). The enzyme was abundantly detected in 9 rhizosphere samples along with Benzoate-CoA ligase (EC 188.8.131.52) and Glutaryl-CoA dehydrogenase (EC 184.108.40.206), both of which were traced in two samples (Figure 2). Acetyl-CoA acetyltransferase enzyme participates in the benzoate degradation via CoA ligation and is also involved in the metabolic pathways including those of aromatic amino acid metabolism6. Aerobic degradation of benzoate is also reported in different microbial communities including gram-negative proteobacterium Azoarcus evansii and gram-positive Bacillus stearothermophilus (Zaar et al., 2001).
Seven different enzymes for benzoate transport and degradation cluster were detected in the rhizosphere samples including benzoate inner membrane translocator and its precursor (Figure 3). Benzoate-CoA ligase, which was identified in two samples (ERS025421 and ERS238578) is the initial enzyme of benzoate metabolism in anaerobic organism Geobacter metallireducens (Wischgoll et al., 2005). Furthermore, the enzymes/transcripts for benzoate 1,2-dioxygenase beta subunit (EC 220.127.116.11) and muconolactone isomerase (EC 18.104.22.168) were also detected in the samples ERS238578 and ERS238586, respectively (Figure 3). Benzoate 1,2-dioxygenase catalyzes aromatic ring hydroxylation (Yamaguchi and Fujisawa, 1980) while the enzyme muconolactone isomerase involves in the catabolism of catechol to succinate (Katti et al., 1989). The alpha and beta subunits of benzoate 1,2-dioxygenase (EC 22.214.171.124 and EC 126.96.36.199), benzoate dioxygenase, ferredoxin reductase component and ortho-halobenzoate 1,2-dioxygenase alpha-ISP protein OhbB were also traced in the metatranscriptome samples. All these enzymes contribute to the degradation of benzoates in the environment significantly and therefore, their abundance in the wheat rhizosphere indicates their rhizoremedial significance.
FIGURE 3. Abundance of transcripts related to benzoate degradation identified through KEGG analysis of the wheat rhizosphere metatranscriptome datasets. [a: 4-oxalocrotonate decarboxylase (EC 188.8.131.52); b: Benzoate transport, inner-membrane translocator; c: Benzoate transport, inner-membrane translocator precursor; d: Benzoate-CoA ligase (EC 184.108.40.206); e: Benzoyl-CoA oxygenase component B; f: Predicted 2-keto-4-pentenoate hydratase/2-oxohepta-3-ene-1,7-dioic acid hydratase; g: Shikimate kinase I (EC 220.127.116.11); h: 2-hydroxycyclohexanecarboxyl-CoA dehydrogenase (EC 1.1.1.-); i: 3-hydroxybutyryl-CoA dehydrogenase (EC 18.104.22.168); j: Acetyl-CoA acetyltransferase (EC 22.214.171.124); k: Benzoate-CoA ligase (EC 126.96.36.199); l: Glutaryl-CoA dehydrogenase (EC 188.8.131.52); m: Benzoate 1,2-dioxygenase beta subunit (EC 184.108.40.206); n: Muconolactone isomerase (EC 220.127.116.11); o: Benzoate 1,2-dioxygenase alpha subunit (EC 18.104.22.168); p: Benzoate 1,2-dioxygenase beta subunit (EC 22.214.171.124); q: benzoate dioxygenase, ferredoxin reductase component; r: Ortho-halobenzoate 1,2-dioxygenase alpha-ISP protein OhbB; s: 2,3-dihydroxybenzoate decarboxylase (EC:126.96.36.199); t: benzoate-CoA ligase (EC:188.8.131.52); u: benzoate/toluate 1,2-dioxygenase subunit beta (EC:184.108.40.206 1.14.12.-); v: benzoate/toluate 1,2-dioxygenase electron transfer component; w: muconolactone D-isomerase (EC:220.127.116.11); x: 2-pyrone-4,6-dicarboxylate lactonase (EC:18.104.22.168); y: 4-oxalmesaconate hydratase (EC:22.214.171.124); z: 3-carboxy-cis,cis-muconate cycloisomerase (EC:126.96.36.199); aa: 4-carboxymuconolactone decarboxylase (EC:188.8.131.52)].
Two transcripts for the enzyme 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate hydrolase (HOPD hydrolase) related to the carbazole degradation were identified in the wheat metatranscriptome (Figure 2). Degradation of carbazole, the N-heterocyclic aromatic compound, members of which occupy components of petroleum products is reported by dioxygenation, a process majorly implicated by Proteobacterial and Actinobacterial phyla (Salam et al., 2017). Many bacterial species possessing these genes degrade PCBs (Masai et al., 1997). Owing to these reasons, carbazole degradation is considered as the remediation strategies for rhizosphere soils using bacterial species (Yang et al., 2009).
Salicylate, Gentisare and Quinate Pathway
Salicylate and gentisate are recognized as the intermediates in the naphthalene metabolism (Balachandran et al., 2012). We report high abundance of the transcripts (05) linked with the gentisare degradation, followed by those of salicylate and gentisate catabolism (07) and quinate degradation (01). We have also detected high abundance of transcripts (14) linked with the enzymes of homogentisate pathway, a central catabolic pathway for the degradation of phenylalanine, tyrosine and hydroxyphenylacetic acid in microorganisms including Pseudomonas putida (Arias-Barrau et al., 2004). The transcripts related to the enzymes fumarylacetoacetase, fumarylacetoacetate hydrolase family protein, maleylacetoacetate isomerase and 3-dehydroquinate dehydratase were detected in multiple metatranscriptomes for salicylate, gentisare and quinate metabolism. Apart from these, the enzymes homogentisate 1,2-dioxygenase and 4-hydroxyphenylpyruvate dioxygenase linked with the homogentisate pathway of aromatic compound degradation were also detected (Table 1). The enzymes fumarylacetoacetase hydrolase and maleylacetoacetate isomerase catalyze hydrolytic breakdown of fumarylacetoacetate and metabolic degradation of phenylalanine and tyrosine (Polekhina et al., 2001). Homogentisate, after its formation from phenylalanine and tyrosine through hydroxylation is further catabolized to fumerate and acetoacetate via three enzymes homogentisate dioxygenase (HmgA), fumarylacetoacetate hydrolase (HmgB) and maleylacetoacetate isomerase (HmgC). Interestingly, the transcripts related to all these three enzymes were detected in the dataset, reflecting an evidence of utilization of homogentisate compounds by the native wheat rhizosphere microbial communities. In the soils, salicylate, quinate and gentisare are supposed to be derived from the organic matter residues of plant sources (Livermore et al., 2014). In the terrestrial habitat where the plant leaf litter and other decaying organic debris are consistently available as carbon source, the abundance of such transcripts linked with the potential enzymes having capability to degrade phenylalanine and tyrosine probably indicates more conserved carbon pathways of the native microbial communities.
n-Phenylalkanoic Acid Degradation
The esters of poly-phenylalkanoates (PHAs) are novel microbial biodegradable plastics usually formed from 3-hydroxy-n-phenylalkanoic acids (Garcia et al., 1999). Under nutrient limiting conditions but with the availability of excess carbon sources, PHAs are accumulated in the microorganisms either in the culture conditions or in the natural habitats like soil and water. The transcripts related to the enzymes 3-hydroxyacyl-CoA dehydrogenase (EC 184.108.40.206), 3-ketoacyl-CoA thiolase (EC 220.127.116.11), Enoyl-CoA hydratase (EC 18.104.22.168) and Long-chain-fatty-acid–CoA ligase (EC 22.214.171.124) that degrade PHAs are most abundantly (26) present in the wheat rhizopshere metatranscriptome samples (Table 1). PHAs are known to be produced intracellular by bacteria and archea as storage carbon and energy source (Koller et al., 2011). It can be inferred that after the decay of the bacterial cells with PHAs as storage compounds, living microbial communities may find access of these compounds under unfavorable conditions of nutrient limitations. The dominance of the transcripts related to PHA degrading enzymes in the wheat rhizosphere metatranscriptomes is supposed to enable associated microbial communities maintain the rhizospheric ecosystem through carbon recycling (Ray and Kalia, 2017).
Phenylacetyl-CoA Catabolic Pathway
Phenylacetyl-CoA, a component of various substrates like phenylalanine, lignin-related phenylpropane and phenylalkanoic acid is amongst the major environmental contaminants (Teufel et al., 2010) and intermediate product of microbial degradation (Díaz et al., 2001). This is also an effector molecule of the TetR family regulator isolated from the extreme thermophilic bacterium Thermus thermophilus HB8 and is supposed to be involved in phenylacetic acid degradation (Sakamoto et al., 2011). Various genes/enzymes such as Phenylacetate-CoA oxygenase (PaaG subunit), Phenylacetate-CoA oxygenase (PaaH subunit), Phenylacetate-CoA oxygenase/reductase (PaaK subunit), Phenylacetic acid degradation protein PaaD (thioesterase), Trancripts related to phenylacetic acid degradation protein PaaE, PaaN, ketothiolase and ring-opening aldehyde dehydrogenase (EC 126.96.36.199) were identified for phenylacetyl-CoA catabolic pathway in the metatranscriptome samples (Table 1). In the dataset, transcripts linked with the enzymes for 4-hydroxyphenylacetic acid catabolic pathway (2-hydroxyhepta-2,4-diene-1,7-dioate isomerase, 4-hydroxyphenylacetate 3-monooxygenase and 5-carboxymethyl-2-oxo-hex-3-ene-1,7-dioate decarboxylase were also detected (Table 1). Since phenylacetic acid is among the common intermediate in the pathway of structurally similar aromatic compounds, the presence of the transcripts linked with its degradation reflects the remediation capability of wheat rhizosphere microbial communities. Abundance of these transcripts is indication of microbial community-linked functions involved in phenylacetyl-CoA metabolism.
The beta-ketoadipate pathway is a chromosomally encoded convergent pathway involved in degradation of aromatic compounds and is in general present in soil bacteria and fungi (Harwood and Parales, 1996). Transcripts for two different branches of this pathway, i.e., catechol branch of the beta-ketoadipate pathway (muconolactone isomerase, EC 188.8.131.52) and protocatechuate branch of the beta-ketoadipate pathway (3-carboxy-cis,cis-muconate cycloisomerase, EC 184.108.40.206) were identified in the rhizosphere metatranscriptome dataset (Table 1). The β-ketoadipate pathway (β-KAP) is involved in the conversion of hazardous aromatic pollutants into simple metabolites like tricarboxylic acid. Studies have also shown that this pathway for conversion of aromatic lignocellulosic waste into bio-oils can be used for producing biodiesel (Wells and Ragauskas, 2012).
Phenylpropanoid and Biphenyl Pathway
Phenols and related phenolic compounds are major organic pollutants causing severe damage to the environment (Diez, 2010). Phenols also adversely affect microbial taxa richness, which can influence the degradation of other pollutants (Barrios-Martinez et al., 2006; Cordova-Rosa et al., 2009). Phenol hydroxylase is reported from many bacterial species to catalyze the catabolism of phenol (Silva et al., 2013). A transcript related to the enzyme phenol hydroxylase, P3 oxygenase component DmpN was identified in the metatranscriptome reflecting that wheat rhizosphere possesses microbial communities for rhizoremedial significance (Nordlund et al., 1993; Kirchner et al., 2003). Also transcripts for the enzymes linked with the phenylpropanoid compound degradation were also identified. Chlorogenate esterase, henylpropionate dioxygenase and related ring-hydroxylating dioxygenases, large terminal subunit, transcriptional regulator for ferulate or vanillate catabolism and vanillate O-demethylase oxygenase subunit were traced in the metatranscriptome. Apart from this, the transcript for related to the enzyme 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoate hydrolase responsible for biphenyl degradation was also identified in the wheat rhizosphere. Besides industrial source of phenolic compounds, a major source of phenolics in the agricultural soils are the plants and native organisms. Presence of transcripts for phenolics degrading enzymes suggests that the microbial communities in the wheat rhizosphere may potentially degrade these compounds as their source of carbon and energy. PCBs are persistent pollutants and their aerobic degradation by soil microorganisms involves biphenyl and benzoate pathways (Leewis et al., 2016).
Pathways for Xenobiotics Degradation
Apart from the identification of transcripts/enzymes in the pathways involved in the metabolism of aromatic compounds in wheat rhizosphere metatranscriptomes, we have also explored the transcripts and pathways linked with the xenobiotic degradation and metabolism. Overall, 37 transcripts of 16 different enzymes were identified from 6 pathways linked with the degradation of xenobiotic compounds (Table 2). Connected with the degradation of xenobiotics, six pathways and their associated enzymes were identified in 9 samples (ERR263035, ERR031114, ERR263026, ERR263028, ERR263029, ERR263033, ERR263031, ERR263036, and ERR263034) with varied levels of abundance (Table 2). These pathways were linked with the degradation of chlorocyclohexane, chlorobenzene, benzoate, aminobenzoate, nitrotoluene, styrene and caprolactam. Among both of these metabolic processes, i.e., metabolism of aromatic compounds and xenobiotic degradation, benzoate degradation pathways was commonly identified.
TABLE 2. Identified active pathways, genes/enzymes and their abundance in xenobiotics degradation in wheat rhizosphere metatranscriptome samples.
Analysis of the datasets further revealed that transcripts linked with the benzoate degradation genes/enzymes were high in abundance (Table 2). Enzymes like 2,3 dihydrobenzoate decarboxylate (DHDB) that reversibly catalyzes decarboxylation of 2,3-dihydroxybenzoate into catechol in Clostridium hydroxybenzoicum (He and Wiegel, 1996), benzoate-CoA ligase (badA) from Magnetospirillum sp. strain TS-6 that initially catalyzes benzoate degradation (Kawaguchi et al., 2006) and benB-xylZ-beta subunit of benzoate/toluene 1,2-dioxygenase that catalyze 1,2-dihydroxylation of different benzoate variants (Haddad et al., 2001) were detected (Figure 3).
We have also detected electron transfer gene component of the enzyme benzoate 1,2-dioxygenase encoded by benC as has been reported by Neidle et al. (1991). Similarly, catC; muconolactone D-isomerase which was traced in the metatranscriptome dataset was also identified from P. putida.7 Interestingly, the dataset showed transcript matching with the enzymes like ligI; 2-pyrone-4,6-dicarboxylate lactonase (Figure 3) that catalyzes the degradation of protocatechuates, methylgallate and syringate in Pseudomonas species8, ligJ; 4-oxalmesaconate hydratase from Sphingomonas sp. SKY69 and pcaB-3-carboxy-cis,cis-muconate cycloisomerase and pcaC-4-carboxymuconolactone decarboxylase, muconate lactonizing enzymes in microorganisms responsible in dehalogenation of muconate derivatives of haloaromatic xenobiotic compounds (Kajander et al., 2002) (Table 2). These results indicated the presence of diverse microbial communities and functions linked with the degradation capabilities of aromatic compounds and xenobiotics, which may otherwise have negative impact on the crop plants. Our analysis further supplements earlier observations (Kuiper et al., 2004; Thomas and Cébron, 2016) that suggested the role of rhizosphere microbial inhabitants in rhizoremediation for plant and soil benefits. Rhizosphere of crops being highly rich in diversity of microbial communities are supposed to be extremely metabolically active and multifunctional (Li et al., 2014). This study shows rich diversity of microbial phyla at both the taxonomic and functional level. The significance of the abundance of functional genes/transcripts in the wheat metatranscriptome samples may possibly be marked as an adaptation strategy of soil microbial communities in the rhizosphere in the presence of toxic aromatic compounds including pesticides (Xu et al., 2013).
While analyzing the datasets for identification of pathways related to xenobiotics degradation and metabolism, transcripts for the enzymes involved in the degradation of chlorocyclohexane, chlorobenzene, nitrotoluene, styrene and caprolactam were also identified (Table 2). Chlorobenzene causes narcosis, tremors, restlessness, and muscle spasms in animals after inhalation while in humans, its chronic or long term exposure influences the central nervous system (CNS)10. Styrene, on exposure, causes eye irritation, gastrointestinal effects, headache, fatigue, weakness, and depression etc. in humans, depending on the duration of exposure11. Similarly, caprolactam also influences the human health in a negative way on exposure12. Looking into the harmful ecological perspectives of the xenobiotics, this study has practical significance as it identifies diversely abundant taxonomic phyla in the wheat rhizosphere and relates identified transcripts to the enzymes linked with the metabolic pathways known for the organic contaminant metabolism or degradation in the soils.
Overall, the dominance of function-linked microbial communities in the rhizosphere of wheat further signifies their role in rhizoremediation of organic pollutants around the plant roots, an area which needs high attention under the conditions of increasing anthropogenic interactions belowground. The study also indicates that rhizosphere microbial communities possess capabilities to utilize the materials like PHAs produced and accumulated by microorganisms themselves under nutrient limiting conditions and deprived state of carbon and energy. The communities also degrade diverse phenylpropanoids introduced to the soils as the debris of leaf litter from plants for the source of carbon. Thus, the metatranscriptome analysis reflected both the pollutant remediation role for the microbial communities associated with such functions as well as the ecological role of maintaining the carbon cycle within the rhizosphere.
Identification of genes and enzymes involved in the metabolism and degradation of organic compounds and xenobiotics can serve as an efficient strategy for bioremediation approaches. Taxonomic exploration of the metatranscriptomes of wheat rhizosphere led to the identification of diverse taxonomic phyla having proven functional role in pollutant degradation. Genes, enzymes and pathways involved in degradation of different organic pollutants like aromatic amines, benzoates, carbazoles, styrene, phenols and biphenols were identified. Detection of particular genes and enzymes involved in pesticide and xenobiotics degradation and their further utilization for bioremediation of the environmental niches, particularly those contaminated with aromatic hydrocarbons is an important aspect. For the biological degradation of aromatic compounds, xenobiotics and pesticides, it is critical to understand the molecular mechanisms engaged in enzymatic catalysis, which will make feasible designing of novel substitutes along with proficient tools for the treatment of pesticide residues and bioremediation of contaminated locations. The study reflected the dominance of taxonomic phyla abundance of transcripts and enzymes linked with the ecological functions of degradation and metabolism. It further reflects that because of the dominance of potential microbial functionalities that link with the metabolic degradation processes for detoxifying or minimizing the impacts of the contaminants in the soils, crop rhizosphere becomes able to tolerate and/or resist harmful organic chemicals. The information on presence of the dominant microbial communities in the rhizosphere obtained through non-cultivable approach may open avenues for the exhaustive isolation and characterization of potential species with the soil bioremediation functions.
RP and DS designed the study, performed the analysis and wrote the manuscript. VG and MV edited the manuscript and made modifications.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
RP is thankful to the Science and Engineering Research Board (SERB) for financial support in terms of SERB National Post-Doctoral Fellowship (Fellowship Reference No. PDF/2016/000714). Funding support and help in data analysis from the Network Project on Agricultural Bioinformatics at Center of Agricultural Bioinformatics (CABin), ICAR-IASRI, New Delhi, India is gratefully acknowledged.
- ^ https://www.ebi.ac.uk/metagenomics/
- ^ http://metagenomics.anl.gov/
- ^ ftp://ftp.metagenomics.anl.gov/data/manual/mg-rast-tech-report-v3_r4.pdf
- ^ https://www.ebi.ac.uk/metagenomics/
- ^ http://www.uniprot.org/uniprot/Q51697
- ^ https://biocyc.org/META/NEW-IMAGE?type=PATHWAY&object=PWY-283
- ^ https://www.wikigenes.org/e/gene/e/1045925.html
- ^ https://biocyc.org/META/NEW-IMAGE?type=REACTION&object=RXN-2462
- ^ https://biocyc.org/gene?orgid=META&id=MONOMER-16015
- ^ https://www.epa.gov/sites/production/files/2016-09/documents/chlorobenzene.pdf
- ^ https://www.epa.gov/sites/production/files/2016-09/documents/styrene.pdf
- ^ https://www.epa.gov/sites/production/files/2016-09/documents/caprolactam.pdf
Ai, C., Liang, G., Sun, J., Wang, X., He, W., Zhou, W., et al. (2015). Reduced dependence if rhizospheric microbiome on plant-derived carbon in 32-year long-term inorganic and organic fertilized soils. Soil Biol. Biochem. 80, 70–78. doi: 10.1016/j.soilbio.2014.09.028
Alvarez, A., Saez, J. M., Costa, J. S. D., Colin, V. L., Fuentes, M. S., Cuozzo, A. S., et al. (2017). Actinobacteria: current research and perspectives for bioremediation of pesticides and heavy metals. Chemosphere 166, 41–62. doi: 10.1016/j.chemosphere.2016.09.070
Arai, H., Ohishi, T., Chang, M. Y., and Kudo, T. (2000). Arrangement and regulation of the genes for meta-pathway enzymes required for degradation of phenol in Comamonas testosteroni TA441. Microbiology 146, 1707–1715. doi: 10.1099/00221287-146-7-1707
Arias-Barrau, E., Olivera, E. R., Luengo, J. M., Fernández, C., Galán, B., García, J. L., et al. (2004). The homogentisate pathway: a central catabolic pathway involved in the degradation of l-phenylalanine, l-tyrosine, and 3-hydroxyphenylacetate in Pseudomonas putida. J. Bacteriol. 186, 5062–5077. doi: 10.1128/JB.186.15.5062-5077.2004
Balachandran, C., Duraipandiyan, V., Balakrishna, K., and Ignacimuthu, S. (2012). Petroleum and polycyclic aromatic hydrocarbons (PAHs) degradation and naphthalene metabolism in Streptomyces sp. (ERI-CPDA-1) isolated from oil contaminated soil. Bioresour. Technol. 112, 83–90. doi: 10.1016/j.biortech.2012.02.059
Balcom, I. N., Driscoll, H., Vincent, J., and Leduc, M. (2016). Metagenomic analysis of an ecological wastewater treatment plant’s microbial communities and their potential to metabolize pharmaceuticals. F1000Res. 5:1881. doi: 10.12688/f1000research.9157.1
Barka, E. A., Vatsa, P., Sanchez, L., Gaveau-Vaillant, N., Jacquard, C., Klenk, H. P., et al. (2016). Taxonomy, physiology, and natural products of Actinobacteria. Microbiol. Mol. Biol. Rev. 80, 1–43. doi: 10.1128/MMBR.00019-15
Barragán, M. J. L., Carmona, M., Zamarro, M. T., Thiele, B., Boll, M., Fuchs, G., et al. (2012). Plastics derived from biological sources: present, and future: a technical, and environmental review. Chem. Rev. 112, 2082–2099. doi: 10.1021/cr200162d
Barrios-Martinez, A., Barbot, E., Marrot, B., Moulin, P., and Roche, N. (2006). Degradation of synthetic phenol-containing wastewaters by MBR. J. Membr. Sci. 281, 288–296. doi: 10.1016/j.memsci.2006.03.048
Bisht, S., Pandey, P., Bhargava, B., Sharma, S., Kumar, V., and Sharma, K. D. (2015). Bioremediation of polyaromatic hydrocarbons (PAHs) using rhizosphere technology. Braz. J. Microbiol. 46, 7–21. doi: 10.1590/S1517-838246120131354
Bond, J. G., Marina, C. F., and Williams, T. (2004). The naturally derived insecticide spinosad is highly toxic to Aedes and Anopheles mosquito larvae. Med. Vet. Entomol. 18, 50–56. doi: 10.1111/j.0269-283X.2004.0480.x
Bulgarelli, D., Schlaeppi, K., Spaepen, S., Ver Loren van Themaat, E., and Schulze-Lefert, P. (2013). Structure and functions of the bacterial microbiota of plants. Annu. Rev. Plant Biol. 64, 807–838. doi: 10.1146/annurev-arplant-050312-120106
Burrows, H. D., Canle, L. M., Santaballa, J. A., and Steenken, S. (2002). Reaction pathways and mechanisms of photodegradation of pesticides. J. Photochem. Photobiol. B Biol. 67, 71–108. doi: 10.1016/S1011-1344(02)00277-4
Castañeda, L. E., and Barbosa, O. (2017). Metagenomic analysis exploring taxonomic and functional diversity of soil microbial communities in Chilean vineyards and surrounding native forests. PeerJ 5:e3098. doi: 10.7717/peerj.3098
Cordova-Rosa, S. M., Dams, R. I., Cordova-Rosa, E. V., Radetski, M. R., Corrêa, A. X. R., and Radetski, C. M. (2009). Remediation of phenol-contaminated soil by a bacterial consortium and Acinetobacter calcoaceticus isolated from an industrial wastewater treatment plant. J. Hazard. Mater. 164, 61–66. doi: 10.1016/j.jhazmat.2008.07.120
Dechesne, A., Badawi, N., Aamand, J., and Smets, B. F. (2014). Fine scale spatial variability of microbial pesticide degradation in soil: scales, controlling factors, and implications. Front. Microbiol. 5:667. doi: 10.3389/fmicb.2014.00667
Dellamatrice, D. M., Silva-Stenico, M. E., de Moraes, L. A. B., Fiore, M. F., and Monteiro, R. T. R. (2017). Degradation of textile dyes by cyanobacteria. Braz. J. Microbiol. 48, 25–31. doi: 10.1016/j.bjm.2016.09.012
Donn, S., Kirkegaard, J. A., Perera, G., Richardson, A. E., and Watt, R. (2015). Evolution of bacterial communities in the wheat crop rhizosphere. Environ. Microbiol. 17, 610–621. doi: 10.1111/1462-2920.12452
Duffner, F. M., Kirchner, U., Bauer, M. P., and Müller, R. (2000). Phenol/cresol degradation by the thermophilic Bacillus thermoglucosidasius A7: cloning and sequence analysis of five genes involved in the pathway. Gene 256, 215–221. doi: 10.1016/S0378-1119(00)00352-8
Edwards, J., Johnson, C., Santos-Medellín, C., Lurie, E., Podishetty, N. K., Bhatnagar, S., et al. (2015). Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. U.S.A. 112, E911–E920. doi: 10.1073/pnas.1414592112
Fuentes, M. S., Benimeli, C. S., Cuozzo, S. A., and Amoroso, M. J. (2010). Isolation of pesticide-degrading actinomycetes from a contaminated site: Bacterial growth, removal and dechlorination of organochlorine pesticides. Int. Biodeter. Biodegrad. 64, 434–441. doi: 10.1016/j.ibiod.2010.05.001
Garcia, B., Olivera, E. R., Miñambres, B., Fernández-Valverde, M., Cañedo, L. M., Prieto, M. A., et al. (1999). Novel biodegradable aromatic plastics from a bacterial source. J. Biol. Chem. 274, 29228–29241. doi: 10.1074/jbc.274.41.29228
Garcia-Salamanca, A., Molina-Henares, M. A., van Dillewijn, P., Solano, J., Pizarro-Tobias, P., Roca, A., et al. (2013). Bacterial diversity in the rhizosphere of maize and the surrounding carbonate-rich bulk soil. Microb. Biotechnol. 6, 36–44. doi: 10.1111/j.1751-7915.2012.00358.x
Ghosal, D., Ghosh, S., Dutta, T. K., and Ahn, Y. (2016). Current state of knowledge in microbial degradation of polycyclic aromatic hydrocarbons (PAHs): a review. Front. Microbiol. 7:1369. doi: 10.3389/fmicb.2016.01369
Haby, P. A., and Crowley, D. E. (1996). Biodegradation of 3-chlorobenzoate as affected by rhizodeposition and selected carbon substrates. J. Environ. Qual. 25, 304–310. doi: 10.2134/jeq1996.00472425002500020014x
Haddad, S., Eby, D. M., and Neidle, E. L. (2001). Cloning and expression of the benzoate dioxygenase genes from Rhodococcus sp. Strain 19070. Appl. Environ. Microbiol. 67, 2507–2514. doi: 10.1128/AEM.67.6.2507-2514.2001
He, Z., and Wiegel, J. (1996). Purification and characterization of an oxygen-sensitive, reversible 3,4-dihydroxybenzoate decarboxylase from Clostridium hydroxybenzoicum. J. Bacteriol. 178, 3539–3543.
Horsman, G. P., Ke, J., Dai, S., Seah, S. Y. K., Bolin, J. T., and Eltis, L. D. (2006). Kinetic and structural insight into the mechanism of BphD, a C-C bond hydrolase from the biphenyl degradation pathway. Biochemistry 45, 11071–11086. doi: 10.1021/bi0611098
Ibrahim, W. M., Karam, M. A., El-Shahat, R. M., and Adway, A. A. (2014). Biodegradation and utilization of organophosphorus pesticide malathion by cyanobacteria. Biomed Res. Int. 2014:392682. doi: 10.1155/2014/392682
Imchen, M., Kumavath, R., Barh, D., Azevedo, V., Ghosh, P., Viana, M., et al. (2017). Searching for signatures across microbial communities: metagenomic analysis of soil samples from mangrove and other ecosystems. Sci Rep. 7:8859. doi: 10.1038/s41598-017-09254-6
Inceoglu,Ö, Al-Soud, W. A., Salles, J. F., Semenov, A. V., and van Elsas, J. D. (2011). Comparative analysis of bacterial communities in a potato field as determined by pyrosequencing. PLoS One 6:e23321. doi: 10.1371/journal.pone.0023321
Ivarsoon, M., and Holm, N. G. (2008). “Microbial colonization of various habitable niches during alteration of ocean crust,” in Links between Geological Processes, Microbial Activities and Evolution of Life, Vol. 4, (Dordrecht: Springer), 69–111.
Jenkins, S. N., Waite, I. S., Blackburn, A., Husbund, R., Rushton, S. P., Manning, D. C., et al. (2010). Actinobacterial community dynamics in long term managed grasslands. Antonie Van Leeuwenhoek 95, 319–334. doi: 10.1007/s10482-009-9317-8
Johnston-Monje, D., Lundberg, D. S., Lazarovits, G., Reis, V. M., and Raizada, M. N. (2016). Bacterial populations in juvenile maize rhizospheres originate from both seed and soil. Plant Soil 405, 337–355. doi: 10.1007/s11104-016-2826-0
Joutey, N. T., Bahafid, W., Sayel, H., and Ghachtouli, N. E. (2013). “Biodegradation: involved microorganisms and genetically engineered microorganisms,” in Biodegradation - Life of Science, ed. R. Chamy (Rijeka: InTech)
Kajander, T., Merckel, M. C., Thompson, A., Deacon, A. M., Mazur, P., Kozarich, J. W., et al. (2002). The structure of Neurospora crassa 3-carboxy-cis,cis-muconate lactonizing enzyme, a beta propeller cycloisomerase. Structure 10, 483–492. doi: 10.1016/S0969-2126(02)00744-X
Kawaguchi, K., Shinoda, Y., Yurimoto, H., Sakai, Y., and Kato, N. (2006). Purification and characterization of benzoate-CoA ligase from Magnetospirillum sp. strain TS-6 capable of aerobic and anaerobic degradation of aromatic compounds. FEMS Microbiol. Lett. 257, 208–213. doi: 10.1111/j.1574-6968.2006.00165.x
Kehinde, F. O., and Isaac, S. A. (2016). Effectiveness of augmented consortia of Bacillus coagulans, Citrobacter koseri and Serratia ficaria in the degradation of diesel polluted soil supplemented with pig dung. Afr. J. Microbiol. Res. 10, 1637–1644. doi: 10.5897/AJMR2016.8249
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A., and Kuramae, E. E. (2016). The ecology of acidobacteria: moving beyond genes and genomes. Front. Microbiol. 7:744. doi: 10.3389/fmicb.2016.00744
Kim, E., Zylstra, G. J., Freeman, J. P., Heinze, T. M., Deck, J., and Cerniglia, C. E. (2006). Evidence for the role of 2-hydroxychromene-2-carboxylate isomerase in the degradation of anthracene by Sphingomonas yanoikuyae B1. FEMS Microbiol. Lett. 153, 479–484. doi: 10.1111/j.1574-6968.1997.tb12613.x
Kirchner, U., Westphal, A. H., Müller, R., and van Berkel, W. J. M. (2003). Phenol hydroxylase from Bacillus thermoglucosidasius A7, a two-protein component monooxygenase with a dual role for FAD. J. Biol. Chem. 278, 47545–47553. doi: 10.1074/jbc.M307397200
Koller, M., Gasser, I., Schmid, F., and Berg, G. (2011). Linking ecology with economy: insights into polyhydroxyalkanoate-producing microorganisms. Eng. Life Sci. 11, 222–237. doi: 10.1002/elsc.201000190
Koutros, S., Lynch, C. F., Ma, X., Lee, W. J., Hoppin, J. A., Christensen, C. H., et al. (2009). Aromatic amine pesticide use and human cancer risk: results from the U.S. Agricultural Health Study. Int. J. Cancer 124, 1206–1212. doi: 10.1002/ijc.24020
Kronenberg, M., Trably, E., Bernet, N., and Patureau, D. (2017). Biodegradation of polycyclic aromatic hydrocarbons: using microbial bioelectrochemical systems to overcome an impasse. Environ. Pollut. 231(Pt 1), 509–523. doi: 10.1016/j.envpol.2017.08.048
Kuiper, I., Lagendijk, E. L., Bloemberg, G. V., and Lugtenberg, B. J. (2004). Rhizoremediation: a beneficial plant-microbe interaction. Mol. Plant Microbe Interact. 17, 6–15. doi: 10.1094/MPMI.2004.17.1.6
Ladino-Orjuela, G., Gomes, E., da Silva, R., Salt, C., and Parsons, J. R. (2016). Metabolic pathways for degradation of aromatic hydrocarbons by bacteria. Rev. Environ. Contam. Toxicol. 237, 105–121. doi: 10.1007/978-3-319-23573-8_5
Leewis, M.-C., Uhlik, O., and Leigh, M. B. (2016). Synergistic processing of biphenyl and benzoate: carbon flow through the bacterial community in polychlorinated-biphenyl-contaminated soil. Sci. Rep. 6:22145. doi: 10.1038/srep22145
Lu, L., Huggins, T., Jin, S., Zuo, Y., and Ren, Z. J. (2014). Microbial metabolism and community structure in response to bioelectrochemically enhanced remediation of petroleum hydrocarbon-contaminated soil. Environ. Sci. Technol. 48, 4021–4029. doi: 10.1021/es4057906
Luciana, P., Mondal, P. K., and Alves, M. (2015). “Aromatic amines sources, environmental impact and remediation,” in Pollutants in Buildings, Water and Living Organisms, eds E. Lichtfouse, J. Schwarzbauer, and D. Robert (Berlin: Springer).
Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–90. doi: 10.1038/nature11237
Masai, E., Sugiyama, K., Iwashita, N., Shimizu, S., Hauschild, J. E., Hatta, T., et al. (1997). The bphDEF meta-cleavage pathway genes involved in biphenyl/polychlorinated biphenyl degradation are located on a linear plasmid and separated from the initial bphACB genes in Rhodococcus sp. strain RHA1. Gene 187, 141–149. doi: 10.1016/S0378-1119(96)00748-2
Mendes, L. W., Kuramae, E. E., Navarrete, A. A., van Veen, J. A., and Tsai, S. M. (2014). Taxonomical and functional microbial community selection in soybean rhizosphere. ISME 8, 1577–1587. doi: 10.1038/ismej.2014.17
Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386
Mwangi, K., Boga, H. I., Muigai, A. W., Kiiyukia, C., and Tsanuo, M. K. (2010). Degradation of dichlorodiphenyltrichloroethane (DDT) by bacterial isolates from cultivated and uncultivated soil. Afr. J. Microbiol. Res. 4, 185–196.
Myers, J. P., Antoniou, M. N., Blumberg, B., Carroll, L., Colborn, T., Everett, L. G., et al. (2016). Concerns over use of glyphosate-based herbicides and risks associated with exposures: a consensus statement. Environ. Health 15:19. doi: 10.1186/s12940-016-0117-0
Neidle, E. L., Hartnett, C., Ornston, L. N., Bairoch, A., Rekik, M., and Harayama, S. (1991). Nucleotide sequences of the Acinetobacter calcoaceticus benABC genes for benzoate 1,2-dioxygenase reveal evolutionary relationships among multicomponent oxygenases. J. Bacteriol. 173, 5385–5395. doi: 10.1128/jb.173.17.5385-5395.1991
Nichols, T. D., Wolf, D. C., Roders, H. B., Beyrouty, C. A., and Renolds, C. M. (1997). Rhizosphere microbial populations in contaminated soils. Water Air Soil Pollut. 95, 165–178. doi: 10.1007/BF02406163
Nordlund, I., Powlowski, J., Hagstrom, A., and Shingler, V. (1993). Conservation of regulatory and structural genes for a multi-component phenol hydroxylase within phenol-catabolizing bacteria that utilize a meta-cleavage pathway. J. Gen. Microbiol. 139, 2695–2703. doi: 10.1099/00221287-139-11-2695
Ortiz-Hernández, M. L., Sánchez-Salinas, E., Dantán-González, E., and Castrejón-Godínez, M. L. (2013). “Pesticide biodegradation: mechanisms, genetics and strategies to enhance the process,” in Biodegradation - Life of Science, ed. R. Chamy (Rijeka: InTech).
Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, D206–D214. doi: 10.1093/nar/gkt1226
Pascual, J., Blanco, S., García-López, M., García-Salamanca, A., Bursakov, S. A., Genilloud, O., et al. (2016). Assessing bacterial diversity in the rhizosphere of Thymus zygis growing in the sierra nevada national park (Spain) through culture-dependent and independent approaches. PLoS One 11:e0146558. doi: 10.1371/journal.pone.0146558
Peiffer, J., Spor, A., Koren, O., Jin, Z., Tringe, S. G., Dangl, J. L., et al. (2013). Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc. Natl. Acad. Sci. U.S.A. 110, 6548–6553. doi: 10.1073/pnas.1302837110
Pisa, G., Magnani, G. S., Weber, H., Souza, E. M., Faoro, H., Monteiro, R. A., et al. (2011). Diversity of 16S rRNA genes from bacteria of sugarcane rhizosphere soil. Braz. J. Med. Biol. Res. 44, 1215–1221. doi: 10.1590/S0100-879X2011007500148
Polekhina, G., Board, P. G., Blackburn, A. C., and Parker, M. W. (2001). Crystal structure of maleylacetoacetate isomerase/glutathione transferase zeta reveals the molecular basis for its remarkable catalytic promiscuity. Biochemistry 40, 1567–1576. doi: 10.1021/bi002249z
Polti, M. A., Aparicio, J. D., Benimeli, C. S., and Amoroso, M. J. (2014). “Role of actinobacteria in bioremediation,” in Microbial Biodegradation and Bioremediation, ed. S. Das (New York, NY: Elsevier), 269–286. doi: 10.1016/B978-0-12-800021-2.00011-X
Porto, A. L. M., Melgar, G. Z., Kasemodel, M. C., and Nitschke, M. (2011). “Biodegradation of pesticides, pesticides in the modern world,” in Pesticides Use and Management, ed. M. Stoytcheva (Rijeka: InTech).
Prüss-Ustün, A., Vickers, C., Haefliger, P., and Bertollini, R. (2011). Knowns and unknowns on burden of disease due to chemicals: a systematic review. Environ. Health. 10:9. doi: 10.1186/1476-069X-10-9
Rabinovich, M. L., Bolobova, A. V., and Vasil’chenko, L. G. (2004). Fungal decomposition of natural aromatic structures and xenobiotics: a review. Appl. Biochem. Microbiol. 40, 1–17. doi: 10.1023/B:ABIM.0000010343.73266.08
Ray, S., and Kalia, V. C. (2017). “Biological significance of degradation of polyhydroxyalkanoates,” in Microbial Applications: Bioremediation and Energy, Vol. 1 eds V. P. Kalia and P. Kumar (New York, NY: Springer) 125–139. doi: 10.1007/978-3-319-52666-9_5
Sakamoto, K., Agari, Y., Kuramitsu, S., and Shinkai, A. (2011). Phenylacetyl coenzyme A is an effector molecule of the TetR family transcriptional repressor PaaR from Thermus thermophilus HB8. J. Bacteriol. 193, 4388–4395. doi: 10.1128/JB.05203-11
Schlaeppi, K., Dombrowski, N., Oter, R. G., Ver Loren van Themaat, E., and Schulze-Lefert, P. (2014). Quantitative divergence of the bacterial root microbiota in Arabidopsis thaliana relatives. Proc. Natl. Acad. Sci. U.S.A. 111, 585–592. doi: 10.1073/pnas.1321597111
Shade, A., Hogan, C. S., Klimowicz, A. K., Linske, M., Mcmanus, P. S., and Handelsman, J. (2012). Culturing captures members of the soil rare biosphere. Environ. Microbiol. 14, 2247–2252. doi: 10.1111/j.1462-2920.2012.02817.x
Sharma, R., and Sharma, P. K. (2018). Metatranscriptome sequencing and analysis of agriculture soil provided significant insights about the microbial community structure and function. Ecol. Gen. Genomics 6, 9–15. doi: 10.1016/j.egg.2017.10.001
Shinoda, Y., Sakai, Y., Uenishi, H., Uchihashi, Y., Hiraishi, A., Yukawa, H., et al. (2004). Aerobic and anaerobic toluene degradation by a newly isolated denitrifying bacterium, Thauera sp. strain DNT-1. Appl. Environ. Microbiol. 70, 1385–1392. doi: 10.1128/AEM.70.3.1385-1392.2004
Silva, C. C., Hayden, H., Sawbridge, T., Mele, P., De Paula, S. O., Silva, L. C. F., et al. (2013). Identification of genes and pathways related to phenol degradation in metagenomic libraries from petroleum refinery wastewater. PLoS One 8:e61811. doi: 10.1371/journal.pone.0061811
Simarro, R., González, N., Bautista, L. F., and Molina, M. C. (2013). Assessment of the efficiency of in situ bioremediation techniques in a creosote polluted soil: change in bacterial community. J. Hazard. Mater. 262, 158–167. doi: 10.1016/j.jhazmat.2013.08.025
Sinha, S., Chattopadhyay, P., Pan, I., Chatterjee, S., Chanda, P., Bandyopadhyay, D., et al. (2009). Microbial transformation of xenobiotics for environmental bioremediation. Afr. J. Biotechnol. 8, 6016–6027. doi: 10.5897/AJB09.740
Subramanian, G., Sekar, S., and Sampoornam, S. (1994). Biodegradation and utilization of organophosphorus pesticides by cyanobacteria. Int. Biodeter. Biodegrad. 33, 129–143. doi: 10.1016/0964-8305(94)90032-9
Teufel, R., Mascaraque, V., Ismail, W., Voss, M., Perera, J., Eisenreich, W., et al. (2010). Bacterial phenylalanine and phenylacetate catabolic pathway revealed. Proc. Natl. Acad. Sci. U.S.A. 107, 14390–14395. doi: 10.1073/pnas.1005399107
Thomas, F., and Cébron, A. (2016). Short-Term rhizosphere effect on available carbon sources, phenanthrene degradation, and active microbiome in an aged-contaminated industrial soil. Front. Microbiol. 7:92. doi: 10.3389/fmicb.2016.00092
Tian, Y., and Gao, L. (2014). Bacterial diversity in the rhizosphere of cucumbers grown in soils covering a wide range of cucumber cropping histories and environmental conditions. Microb. Ecol. 68, 794–806. doi: 10.1007/s00248-014-0461-y
Turesky, R. J., and Le Marchand, L. (2011). Metabolism and biomarkers of heterocyclic aromatic amines in molecular epidemiology studies: Lessons learned from aromatic amines. Chem. Res. Toxicol. 24, 1169–1214. doi: 10.1021/tx200135s
Uroz, S., Buée, M., Murat, C., Frey-Klett, P., and Martin, F. (2010). Pyrosequencing reveals a contrasted bacterial diversity between oak rhizosphere and surrounding soil. Environ. Microbiol. Rep. 2, 281–288. doi: 10.1111/j.1758-2229.2009.00117.x
Valderrama, J. A., Durante-Rodríguez, G., Blázquez, B., García, J. L., Carmona, M., and Díaz, E. (2012). Bacterial degradation of benzoate: cross-regulation between aerobic and anaerobic pathways. J. Biol. Chem. 287, 10494–10508. doi: 10.1074/jbc.M111.309005
Valentìn, L., Nousiainen, A., and Mikkonen, A. (2013). “Fungal-mediated degradation of emerging pollutants in sewage sludge,” in Emerging Organic Contaminants in Sludges: Analysis, Fate and Biological Treatment Handbook of Environmental Chemistry, Vol. 24, eds T. Vicent, G. Caminal, and E. Eljarrat (Berlin: Springer-Verlag), 1–30.
Vinas, M., Sebate, J., Espuny, M. J., and Solanas, A. M. (2005). Bacterial community dynamics and polycyclic aromatic hydrocarbon degradation during bioremediation of heavily creosote-contaminated soil. Appl. Environ. Microbiol. 71, 7008–7018. doi: 10.1128/AEM.71.11.7008-7018.2005
Wang, N. F., Zhang, T., Yang, X., Wang, S., Yu, Y., Dong, L. L., et al. (2016). Diversity and composition of bacterial community in soils and lake sediments from an arctic lake area. Front. Microbiol. 7:1170. doi: 10.3389/fmicb.2016.01170
Wilke, A., Harrison, T., Wilkening, J., Field, D., Glass, E. M., Kyrpides, N., et al. (2012). The M5nr: a novel non-redundant database containing protein sequences and annotations from multiple sources and associated tools. BMC Bioinformatics 13:141. doi: 10.1186/1471-2105-13-141
Wischgoll, S., Heintz, D., Peters, F., Erxleben, A., Sarnighausen, E., Reski, R., et al. (2005). Gene clusters involved in anaerobic benzoate degradation of Geobacter metallireducens. Mol. Microbiol. 58, 1238–1252. doi: 10.1111/j.1365-2958.2005.04909.x
Xu, B., Xu, W., Yang, F., Li, J., Yang, Y., Tang, X., et al. (2013). Metagenomic analysis of the Pygmy loris fecal microbiome reveals unique functional capacity related to metabolism of aromatic compounds. PLoS One 8:e56565. doi: 10.1371/journal.pone.0056565
Xun, L., and Sandvik, E. R. (2000). Characterization of 4-hydroxyphenylacetate 3-hydroxylase (HpaB) of Escherichia coli as a reduced flavin adenine dinucleotide-utilizing monooxygenase. Appl. Environ. Microbiol. 66, 481–486. doi: 10.1128/AEM.66.2.481-486.2000
Yang, M., Li, W., Guo, X., Qu, Z., Zhu, X., and Wang, X. (2009). Isolation and identification of a carbazole degradation gene cluster from Sphingomonas sp. JS1. World J. Microbiol. Biotechnol. 25, 1625–1631. doi: 10.1007/s11274-009-0055-x
Yang, Y., Wang, N., Guo, X., Zhang, Y., and Ye, B. (2017). Comparative analysis of bacterial community structure in the rhizosphere of maize by high-throughput pyrosequencing. PLoS One 12:e0178425. doi: 10.1371/journal.pone.0178425
Keywords: metatranscriptome, xenobiotics, pesticides, polyaromatic hydrocarbons, polychlorinated biphenyls, petroleum hydrocarbons, microbial communities
Citation: Singh DP, Prabha R, Gupta VK and Verma MK (2018) Metatranscriptome Analysis Deciphers Multifunctional Genes and Enzymes Linked With the Degradation of Aromatic Compounds and Pesticides in the Wheat Rhizosphere. Front. Microbiol. 9:1331. doi: 10.3389/fmicb.2018.01331
Received: 13 April 2018; Accepted: 31 May 2018;
Published: 03 July 2018.
Edited by:Mariusz Cyco, Medical University of Silesia, Poland
Reviewed by:Jay Prakash Verma, Banaras Hindu University, India
Pushpender Kumar Sharma, Sri Guru Granth Sahib World University, India
Copyright © 2018 Singh, Prabha, Gupta and Verma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors are combined first authors.