Next Article in Journal
The Endocannabinoid System in the Mediterranean Mussel Mytilus galloprovincialis: Possible Mediators of the Immune Activity?
Previous Article in Journal
IgM and IgG Immunoreactivity of SARS-CoV-2 Recombinant M Protein
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Representative Bacillus sp. AM1 from Gut Microbiota Harbor Versatile Molecular Pathways for Bisphenol A Biodegradation

1
Department of Microbiology, Faculty of Pharmacy, University of Granada, Campus of Cartuja, 18071 Granada, Spain
2
Instituto de Nutrición y Tecnología de los Alimentos, INYTA-Granada, 18100 Granada, Spain
3
Department of Biochemistry and Molecular Biology, Faculty of Pharmacy, University of Granada, Campus of Cartuja, 18071 Granada, Spain
4
IBS: Instituto de Investigación Biosanitaria ibs., 18012 Granada, Spain
*
Author to whom correspondence should be addressed.
Authors contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(9), 4952; https://doi.org/10.3390/ijms22094952
Submission received: 7 February 2021 / Revised: 24 February 2021 / Accepted: 24 February 2021 / Published: 7 May 2021
(This article belongs to the Section Molecular Microbiology)

Abstract

:
Human gut microbiota harbors numerous microbial species with molecular enzymatic potential that impact on the eubiosis/dysbiosis and health/disease balances. Microbiota species isolation and description of their specific molecular features remain largely unexplored. In the present study, we focused on the cultivation and selection of species able to tolerate or biodegrade the endocrine disruptor bisphenol A (BPA), a xenobiotic extensively found in food plastic containers. Chemical xenobiotic addition methods for the directed isolation, culturing, Whole Genome Sequencing (WGS), phylogenomic identification, and specific gene-encoding searches have been applied to isolate microorganisms, assess their BPA metabolization potential, and describe encoded catabolic pathways. BPA-tolerant strains were isolated from 30% of infant fecal microbial culture libraries analyzed. Most isolated strains were phylogenetically related to the operational taxonomic group Bacillus amyloliquefaciens spp. Importantly, WGS analysis of microbial representative strain, Bacillus sp. AM1 identified the four complete molecular pathways involved on BPA degradation indicating its versatility and high potential to degrade BPA. Pathways for Exopolysaccharide (EPS) and Polyhydroxyalkanates (PHA) biopolymer synthesis were also identified and phenotypically confirmed by transmission electronic microscopy (TEM). These microbial biopolymers could generally contribute to capture and/or deposit xenobiotics.

1. Introduction

The microorganisms that colonize the human gut are functionally diverse and highly variable [1]. Four phyla, Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria, dominate the human gut. Daily exposure to industrial xenobiotics found in dermal products, food, and beverage containers trigger alterations in gut microbiota composition [2]. Dysbiosis, an imbalance in the composition and metabolic capacity of our microbiota, has been related to the cumulative intake of a wide variety of xenobiotic substances that induce obesogenic phenotypes [3], and act as endocrine disruptors through microbiota interaction [4,5]. Bisphenol A (BPA) and its analogs (bisphenol-B, -S, -F and -E), used as building blocks of polycarbonate plastics, as well as parabens (methyl-, ethyl-, propyl- and phenyl-paraben), used as antimicrobial preservatives, have been associated with metabolic and endocrine/hormonal diseases [6,7]. Phthalates, such as di-(2-ethylhexyl) phthalate, are used as plasticizers and have been shown to exert obesogenic effects, affect the reproduction system, and modify key glycemic parameters [8,9]. Perchlorates, organophosphates [10], and other pesticides trigger symptoms of diabetes and obesity [11].
Numerous gut microorganisms belonging to the Firmicutes phylum possess a diverse and versatile array of enzymes. Species within the Bacillus genus are consistently identified in the human gut and are a source of enzymes with potential technological applications in the modern food sector [12]. The use of bacterial enzymes in the current food industry has gained increasing recognition in recent years as a result of the demand for new and better products [13]. In this sense, due to the increase in the human population and the growing need for products to satisfy the nutritional requirements of an increasingly demanding population, enzymes of bacterial origin are presented as among the most versatile alternatives in the current industrial context [14]. Commensal microorganisms isolated from human microbiota could fulfill the criteria of safety assessment and the status of Qualified Presumption of Safety (QPS) [15,16]. In this sense, most species within the Bacillus subtilis cluster are considered QPS [17] and are increasingly considered as feed additives [18].
Taxonomic classification of Bacillus species and the specific WGS description of their encoded biochemical pathways are needed to consider the food and feed safety of Bacillus strains [19]. Ecological and phenotypical heterogeneity within the Bacillus genus have rendered pure 16S rDNA-based taxonomy inefficient to discern between Bacillus species, even leading to misclassifications in the past [20]. Next-Generation Sequencing platforms and in silico tools have been developed to perform appropriate phylogenomic classifications and predict metabolic pathways involved in industrial and environmental processes, expanding the molecular knowledge of the enzymatic arsenal for BPA degradation in Bacillus species. Few microbial cultivation studies explored the potential BPA degrading capacity of Bacillus species for environmental bioremediation [21]. Recent approaches to assess the degradation pathways of BPA in Bacillus and other microorganisms have provided a new field of research to understand in the context of industrial and environmental biotechnology [22,23].
BPA is a dietary xenobiotic banned at specific levels. Because gut microbiota is a potential source for screening the capacities of tolerance and biodegradation capacities, particular microorganisms might provide endogenous bioremediation and could become detoxifying probiotics. Thus, the main objective of the present study was to describe, for the first time, BPA-tolerant microorganisms from the human gut and the specific pathways that allow them to behave as active metabolizers of this endocrine disruptor, possibly modulating its physiological impact. Culturing isolation strategies and WGS descriptions are used to find the specific BPA-degrading genes involved in the degradation pathways.

2. Results and Discussion

2.1. BPA-Tolerant Microorganisms Isolated from Human Microbiota: 16S RNA Partial Gene Identification and Phylogenetic Analysis

Our hypothesis was to seek human microbes capable of tolerating and metabolizing the endocrine disruptor BPA [4], the most important bisphenol in a large family of chemicals that are utilized to produce polycarbonates and epoxy resins and found in foods, drinks, and baby bottles. A comprehensive approach involved the culturing and isolation of gut microbial species in BPA-enriched complete or minimal media for further morphology, phylogenetical, biochemical, as well as WGS analyses was undertaken. It is an untested innovative strategy for a rational and functional search of potential next generation probiotics for biotechnological use in BPA detoxification. This study sets the basis for future linked omics and multidisciplinary studies aimed at the identification of new BPA-degrading microbes in the human gut and their association with the short and long-term daily intake of BPA.
Isolation and identification of BPA-tolerant strains from microbiota samples were successful from the different concentrations used [0.5; 10; 20 and 50 ppm] and demonstrated a clear predominance of Firmicutes phylum taxa. This is one of the most populated phyla in microbiota with wide compositional, clinical, and functional variability [24,25]. Out of these 20 isolates analyzed, 83.3% were related to the genus Bacillus, and 16.7% of the strains were distributed within Staphylococcus, Enterococcus, and Streptococcus genera (Table S1). However, the analysis of these minority genera was beyond the objective of the present work. The phylogenetic tree based on partial 16S RNA gene of Bacillus spp. isolated strains showed a main cluster related to B. amyloliquefaciens, B. velezensis and B. siamensis (Figure 1), being only C1, C2 and B5 isolates different to the rest of Bacillus spp. The repeated presence of closely related B. amyloliquefaciens isolates in approximately 30% of the analyzed fecal samples suggested differential metabolic detoxifying potential of these species in individuals under high dietary exposure level to BPA [7] or food endocrine microbiota disruptors (EMD) [4].
Moreover, preliminary data on parallel metagenomics analyses were in concordance with this percentage of Bacillus sp. presence (data not shown). We hypothesize that the long-term daily intake of xenobiotic compounds such as BPA or parabens present in food containers may select specific taxa groups in the intestinal microbiota, thus increasing or reducing microbial populations according to their enzymatic capacity to tolerate and metabolize BPA. Increased isolation of species Bacillus amyloliquefaciens in individuals might reveal a clinical potential impact through high tolerance to BPA [26].
Interestingly, the potential BPA metabolization by Bacillus spp. has solely been investigated in relation to environmental and bioremediation resources [22,23]. The present study constitutes the first report on BPA-tolerant human microbiota species that might became a new extensive biotechnological source of probiotics for plants, food, and clinical purposes. In this sense, several Bacillus strains are proposed for being anti-obesogenic probiotics [27,28]. Still, the major limitation for using these new potential probiotics is the lack of standardized criteria and protocols for target disorders. Next-generation probiotics for obesity-related diseases promise to improve the rational of clinical trial design and outcomes [29]. Additional microbial food safety and risk assessments would be needed for strains proposed for food and feed, including strain characterization and genome sequencing, screening for undesirable attributes and metabolites, and experimental evidence of safety by appropriately designed safety evaluation studies [18,30].

2.2. Identification of the Representative BPA-Tolerant Microbiota Strain and Morphological Features

All isolated and typified strains were able to tolerate between 0.5 and 50 ppm of BPA. Bacillus sp. AM1 (B1) was the BPA-tolerant selected strain representing the main phylogenetic cluster (Figure 1) for being further analyzed (phenotypic features and WGS). This selection was done based on the following enzymatic and phenotypic characteristics preliminary analyzed: the B1 strain was able to tolerate higher concentration of BPA; B1 showed special morphological and phenotypic features, such as biopolymers EPS and PHA production that can protect from the xenobiotic exposure (Figure 2), specific gene-encoding enzymes for both biopolymers were also found in the WGS of Bacillus sp. AM1 (Tables S2 and S3), and the phylogenetic identification as shown in Figure 1, B1 may represent the major close-related group. This cluster represented by the B1 strain and the grouping of B3, B4, C3, C4, and C5 strains was closely related to the operational B. amyloliquefaciens group [20], and became the most frequently identified related species from the microbiota samples analyzed.
Bacillus sp. AM1 is an anaerobic facultative rod with subterminal swollen spores. Pale yellow colonies forming oleaginous drops in MRS medium are typical (Figure 2a). Similarly, lipid, wax ester, or PHA producers showed similar morphologies and characteristics [31,32]. Transmission electron microscopy (TEM) revealed the presence of clear cytoplasmatic lipid accumulations similar to PHA granule accumulation. We also observed biopolymers type capsular EPS (Figure 2b).
The hydrolytic enzyme capacities of Bacillus sp. AM1 (Table 1) revealed similar results to those described for Bacillus amyloliquefaciens, B. velezensis, and B. siamensis, including inulinase activity. Bacillus genus is considered of interest for its enzyme productivity for many manufacturing processes [33,34,35]. The antibiotic sensitivity testing revealed that Bacillus sp. AM1 was resistant to ampicillin, ceftazidime, penicillin, and amoxicillin/clavulanic acid (Table 2). These resistances were previously described as intrinsic resistance showed commonly by Bacillus spp., not posing any safety concern [18].
Bacillus sp. AM1 was able to grow in BPA-enriched M9 medium (0.5 ppm) showing slower colony development compared to the control and M9 containing glucose as sole carbon source (0.5 ppm). Conversely, Bacillus sp. AM1 was not able to grow on M9 medium containing BPS (0.5 ppm), an analogue of BPA, indicating that BPS may be toxic or cannot be used as sole carbon source by this strain (Figure S1). Bacillus sp. AM1 was capable of degrading 84.68% of BPA after 72 h, using a starting BPA concentration of 25 ppb (Figure 2c). A similar pattern of biodegradation was shown by Bacillus sp. GZB [22].

2.3. Whole Genome Sequence of Bacillus sp. AM1 Analysis

WGS of Bacillus sp. AM1 (GenBank CP047644.1) showed a total length (bp) of 4,203,926 and a G+C content value of 45.98 % (Figure S2). The genome of strain Bacillus sp. AM1 shared 67.67–98.36% ANI values with the Bacillus subtilis cluster (Table S4). DNA–DNA hybridization (isDDH) relatedness between strain Bacillus sp. AM1 and the known Bacillus species ranged from 28.2 to 83%. Results from the ANI and isDDH analyses and the phylogenomic trees based on 16S RNA gene (Figure 3a), gyrB (Figure 3b) and rpoB (Figure 3c) suggested that strain Bacillus sp. AM1 belongs to the safe operational group B. amyloliquefaciens [20], with no pathogenic or virulence factors or toxin genes. Moreover, rMLST analysis showed match to 55 genes and, within them, to 6 genes unique and specific of Bacillus sp. AM1: rpmD, rplM, rplI, rplE, rplA, rpsG. We are aware that many species of Bacillus are difficult to be taxonomically classified under the universal rules; thus additional phylogenomic and phenotypic data will be obtained for a more precise classification. Moreover, EPS and PHA biosynthesis complete gene pathways (Tables S2 and S3) were also identified, confirming the features observed by TEM analyses (Figure 2b).

2.4. Description of Putative BPA Biodegradation Pathways by Microbiota BPA-Tolerant Isolate

BPA biodegradation occurs due to the contribution and coordination of many expressed protein-coding genes. Complete gene-encoded enzymes for the four reference pathways of BPA degradation were found in the Bacillus sp. AM1 genome, a major difference with Bacillus sp. GZB strain, isolated from recalcitrant environment, whose description seemed to lack 30% of these enzymes [22]. Description of reaction steps, specific enzymes, EC number, genes loci, and protein ID from each pathway are shown in Figure 4. This versatile gene arsenal for BPA metabolization enzymes should be confirmed with functional and chemical analyses, but underscores the huge potential of BPA biodegradation, possibly generating intermediate metabolites with more or less estrogenic-linked effects [22,36]. Our results suggest that the establishment of routes linked to the degradation of estrogenic xenobiotics, such as BPA, involves the participation of a large number of microbial enzymes whose joint activity entails the production of a multitude of intermediate less toxic byproducts which can even be found related to other metabolic routes or destined to the energy production through the Cycle of Tricarboxylic Acids [23]. Likewise, numerous studies related to the metabolism of BPA highlighted the central role of cytochrome P450 in bacteria [37,38].
Lobos et al. [39] and Sasaki et al. [40] showed the existence of microorganisms belonging to Sphingomonas and Bacillus genera capable of using BPA as a carbon source through its degradation in multiple stages, with the production of end byproducts of variable toxicity.

3. Materials and Methods

3.1. Tolerant-BPA bacterial Isolation from Human Gut Microbiota with Biotechnological Interest

Ten microbial isolates from fecal human microbiota collections of 0 to 1 year old infants (Isolates B- Project INFABIO) and 6–8-year-old children (Isolates C-Project OBEMIRISK) were obtained by a serial dilution method, exposition to different BPA concentrations [0.5, 10, 20, and 50 ppm] during 72 h and further spreading in BHI/MRS media incubated under aerobic and anaerobic conditions (anaerobic jars anaerocult®) at 37 °C.
The BPA-tolerant bacterial colonies with distinguishing features were isolated as pure culture for subsequent morphological, phenotypic and genotypic identifications: bacterial cell counts, gram staining, spore staining, capsule staining, catalase activity, oxidase, and motility tests.

3.2. Genomic DNA from Tolerant-BPA Bacterial Isolates, Taxonomy Identification and Phylogenetic Analysis

Genomic DNA was extracted using DNAeasy columns (Qiagen®, Germany) following the manufacturing instructions. The isolated DNA was quantified using Nanodrop (Thermo Scientific) and biophotometer (Eppendorf® D30). The quality of DNA was monitored through gel electrophoreses. 16S RNA gene sequencing of selected bacterial strains was done by Sanger method (IPBLN Service). Forward and reverse sequences were provided separately. Reverse sequence was converted to complementary sequence with Chromas Pro 2.0 software (Technelysium Pty Ltd., Tewantin, Australia). Sequences were inspected for maximum homology against GenBank using NCBI’s BLASTn program. The collection of 16S RNA partial gene sequences is done using the Ezbiocloud platform.

3.3. Morphological Features and Identification of BPA-Tolerant Bacterial Strains

The selection of a representative sample from fecal BPA-tolerant microorganisms for further phenotypic and WGS analysis was done according to specific phenotypic features and specific phylogenetic classification.
Transmission electron microscopy was used to examine possible structural changes in the selected strain. The samples were obtained after exposing the bacteria to the different BPA treatments. The techniques detailed below were carried out in the Laboratory for the Preparation of Biological Samples of the Microscopy Service of the Center for Scientific Instrumentation (CIC) of the University of Granada. Samples were fixed with 2.5% glutaraldehyde in sodium cacodylate buffer (pH7.4, 0.05M) for 24 h at 4 °C. After washing them three times in sodium cacodylate buffer (pH7.4, 0.1M) for 15 min, the post-fixation was carried out with 1% osmium tetroxide + 1% potassium ferricyanide in water for 2 h. It was washed again three times in water 15 min. And dehydration was realized using ascending ethanol gradients (50%, 70%, 90%, and 100%) for 15 min. Subsequently, the infiltration was carried out in a mixture of ethanol/resin (Epon 812) 2:1, 2:2 and 1:2 for 1 h and then it was the inclusion in pure resin overnight. After polymerization at 60 °C for 48 h, ultra-fine sections of 700 Angstrom were made (in a Leica Ultracut R ultramicrotome) deposited on copper grids and contrasted with uranyl acetate and lead citrate.
The enzymatic assay were also determinate for starch, carboxymethyl cellulose, gelatin, inulin, tween 80, DNase supplemented media were used to test several enzymes production according to complementary methodologies previously described [41,42].
The antimicrobial resistance patterns of Bacillus sp. AM1 to 17 clinically important antibiotics, included those antibiotics recommending on the EFSA guidelines for testing antimicrobial susceptibility of the Bacillus species [43], were conducted using the agar disc-diffusion method. The strains were characterized as sensitive or resistant based on the size of the inhibition zones around each disc, according to the National Committee for Clinical Laboratory Standards (CLSI) criteria [44].

3.4. BPA Tolerance Testing

BPA biodegradation capacity was tested in TSB liquid cultures of Bacillus sp. AM1 exposed to 25 ppb concentration of BPA at 30 °C during 72 h. Microbial preparations were measured through LC-MS/MS system for BPA quantification. Minimal medium 9 (M9) and M9 containing BPA (0.5 ppm) as sole carbon source were used for the isolation of BPA tolerant or biodegradation strain. M9 containing glucose (0.5 ppm) and M9 containing BPS (0.5 ppm) were also used to compare the growth ability. Chemicals, reagents, instrumentation and software for bisphenols determination were carried out by CIC services under validated procedures previously described by García-Córcoles et al. [45].

3.5. Whole Genome Sequence of BPA-Tolerant Bacterial Strain

Genomic DNA extraction for WGS: Bacillus sp. AM1 (B1) selected for its phenotypic characteristics and its potential in biotechnology was inoculated in 2 mL nutrient broth, followed by overnight incubation at 37 °C. The cells were harvested by centrifugation at 5000× g for 10 min. Genomic DNA was extracted using DNAeasy columns (Qiagen®, Germany) as previously described.
High-throughput sequencing: Shotgun library preparation of Bacillus sp. AM1 (B1) strain was performed using the Nextera DNA Sample Preparation kit (Illumina, USA). Clusters generation and the sequencing were performed using the TruSeq™ PE Cluster kit v2-cBot-HS and TruSeq SBS v3-HS (200 cycle) kit respectively (Ilumina, USA). The sequencing was performed on the HiScan SQ sequencer (Illumina, USA).
Genome assembly and annotation: Hybrid de novo assembly was performed with Illumina short reads (c.a. 440X coverage, 150bp pair end reads) and Oxford Nanopore long reads (c.a. 1800X coverage) using Unicycler1 with the “conservative” option and threads 10 parameter [46].

3.6. Whole Genome Sequence Analyses

3.6.1. Phylogenomic Analysis

Genome-based taxonomic classification and genome assembly were searched at Microbial Genomes Atlas (MiGA) online server (http://www.microbial-genomes.org accessed on 20 February 2021) that uses NCBI non-redundant prokaryotic genomes database, JSpeciesWS (http://jspecies.ribohost.com/jspeciesws accessed on 20 February 2021), and TrueBacTM ID system from the EzBioCloud server (https://www.ezbiocloud.net/contents/genome accessed on 20 February 2021). The species identification carried out in TrueBac-ID was based on the comparison of average nucleotide sequence identity (ANI) between the genomes of Bacillus sp. AM1 and the related type strains. For a reliable identification of the strains, the Type (Strain) Genome Server (TYGS) was employed. Prediction from identified Ribosomal Multilocus Sequence Typing (rMLST) alleles linked to genomes (53 rps genes) was used as a means of integrating microbial taxonomy and typing [47]. 16S RNA gene and the housekeeping genes gyrB and rpoB sequences were used for phylogenomic identification.
Prokaryotic genome comparisons were performed by ANI calculator using the OrthoANIu algorithm available on Ezbiocloud [48]. Digital or isDDH analyses were carried out by Genome-to-Genome Distance Calculator (GGDC) (formula 2) [49].

3.6.2. Molecular and Gene Prediction and Annotations

Draft Bacillus sp. AM1 genome was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP). The functional annotation of genes was also carried out using web-based Rapid Annotation Using Subsystem Technology (RAST) annotation server (http://rast.theseed.org/FIG/rast.cgi accessed on 20 February 2021). Protein coding genes were scanned for their organization into operons using the web server Operon-mapper. Functional annotation of genes in terms of Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology assignments and predictions of KEGG pathways were carried out through KEGG Automatic Annotation Server (KAAS) server (https://www.genome.jp/kegg/kaas/ accessed on 20 February 2021) using bi-directional best hit (BBH) method in GhostZ. Similarly, the Bacillus sp. AM1 genome were also scanned for cluster of orthologous groups (COGs) annotations using eggNOG-mapper v2 (http://eggnog-mapper.embl.de accessed on 20 February 2021). Automated Carbohydrate-active enzyme ANnotation web server, (dbCAN meta server (http://bcb.unl.edu/dbCAN2/blast.php accessed on 20 February 2021), was employed to annotate Carbohydrate-Active enZYmes (CAZy).

3.6.3. Identification of Microbial Enzymes Isolated from Bisphenol A Tolerant Microorganisms

Specific genes and amino acid sequences of the BPA-tolerant representative microorganism were searched by screening according to phenotypic and chemical characteristics such as EPS, PHA, and other enzymes related to the BPA degradation pathways: quinol-enzymes, laccase, through BLAST (NCBI), EZBIOCLOUD, BRENDA, UNIPROT, and EXPASY. A more exhaustive and combined analysis using the KEGG database (Kyoto Encyclopedia of Genes and Genomes) allowed corroborating the identity of the enzymes involved in metabolic pathways of BPA degradation previously describe by Li et al. [22]. EC numbers results essential for the screening in WGS Bacillus sp. AM1 protein encoding genes.

4. Conclusions

This study is the first report on the isolation of BPA-degrading microbes from the human gut and its complete description of BPA biodegradation pathways. The isolation of strains belonging to the Firmicutes phyla from human microbiota able to biodegrade BPA was successful through a BPA-enriched media and a directed cultivation approach. WGS analyses of representative Bacillus sp. AM1 allowed taxonomic and phylogenomic identifications, relevant phenotypic features confirmation, and safety assessment aspects. Bacillus sp. AM1 WGS revealed the presence of a complete arsenal of genes that encode enzymes for versatile BPA degradation; this might open avenues for the investigation of the use of these strains for environment, feed and food industries, and/or to explore the clinical modulation impact on endocrine pathogenesis modulation by BPA-degrader species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22094952/s1.

Author Contributions

Conceptualization, M.A.; Methodology, all authors; Microbiota Laboratory analysis, A.L.-M., and A.T.-S.; data curation, subcontracting P.Manfredi; Funding Acquisition, M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by grant number GP/EFSA/ENCO/380 2018/03/G04: OBEMIRISK: Knowledge platform for assessing the risk of Bisphenols on gut microbiota and its role in obesogenic phenotype: looking for biomarkers. This research was also funded by FEDER-Infrastructure: IE_2019-198. The APC was funded by EIN-2019-103082.

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to previous approvals for fecal microbiota collection libraries used.

Acknowledgments

Part of results is from Ana López-Moreno doctoral thesis, Biomedicine Doctorate Program of the University of Granada. A.L.-M. under the programme UGR-Plan Propio Intensificación de la Investigación.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Patents

IPR-823 Application in progress.

References

  1. Gibson, M.K.; Crofts, T.S.; Dantas, G. Antibiotics and the developing infant gut microbiota and resistome. Curr. Opin. Microbiol. 2015, 27, 51–56. [Google Scholar] [CrossRef] [Green Version]
  2. Pagliari, D.; Gambassi, G.; Piccirillo, C.A.; Cianci, R. The intricate link among gut “Immunological Niche,” Microbiota, and Xenobiotics in intestinal pathology. Mediators Inflamm. 2017, 8390595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Andújar, N.; Gálvez-Ontiveros, Y.; Zafra-Gómez, A.; Rodrigo, L.; Álvarez-Cubero, M.J.; Aguilera, M.; Monteagudo, C.; Rivas, A.A. Bisphenol A analogues in food and their hormonal and obesogenic effects: A Review. Nutrients 2019, 11, 2136. [Google Scholar] [CrossRef] [Green Version]
  4. Aguilera, M.; Gálvez-Ontiveros, Y.; Rivas, A. Endobolome, a new concept for determining the influence of microbiota disrupting chemicals (MDC) in relation to specific endocrine pathogenesis. Front. Microbiol. 2020, 11, 578007. [Google Scholar] [CrossRef] [PubMed]
  5. Gálvez-Ontiveros, Y.; Páez, S.; Monteagudo, C.; Rivas, A. Endocrine disruptors in food: Impact on gut microbiota and metabolic diseases. Nutrients 2020, 12, 1158. [Google Scholar] [CrossRef]
  6. Verbanck, M.; Canouil, M.; Leloire, A.; Dhennin, V.; Coumoul, X.; Yengo, L.; Froguel, P.; Poulain-Godefroy, O. Low-dose exposure to bisphenols A, F and S of human primary adipocyte impacts coding and non-coding RNA profiles. PLoS ONE 2017, 12, e0179583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Gálvez-Ontiveros, Y.; Moscoso-Ruiz, I.; Rodrigo, L.; Aguilera, M.; Rivas, A.; Zafra-Gómez, A. Presence of parabens and bisphenols in food commonly consumed in Spain. Foods 2021, 10, 92. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Dong, T.; Hu, W.; Wang, X.; Xu, B.; Lin, Z.; Hofer, T.; Stefanoff, P.; Chen, Y.; Wang, X.; et al. Association between exposure to a mixture of Phenols, Pesticides, and Phthalates and Obesity: Comparison of three statistical models. Environ. Int. 2019, 123, 325–336. [Google Scholar] [CrossRef]
  9. Milošević, N.; Milanović, M.; Sudji, J.; Bosić Živanović, D.; Stojanoski, S.; Vuković, B.; Milić, N.; Medić Stojanoska, M. Could Phthalates exposure contribute to the development of metabolic syndrome and liver disease in humans? Environ. Sci. Pollut. Res. Int. 2020, 27, 772–784. [Google Scholar] [CrossRef] [PubMed]
  10. Roman, P.; Cardona, D.; Sempere, L.; Carvajal, F. Microbiota and Organophosphates. Neurotoxicology 2019, 75, 200–208. [Google Scholar] [CrossRef] [PubMed]
  11. Yuan, X.; Pan, Z.; Jin, C.; Ni, Y.; Fu, Z.; Jin, Y. Gut Microbiota: An underestimated and unintended recipient for Pesticide-induced toxicity. Chemosphere 2019, 227, 425–434. [Google Scholar] [CrossRef] [PubMed]
  12. Woldemariam, K.; Wan, Z.; Yu, Q.; Li, H.; Wei, X.; Liu, Y.; Wang, J.; Sun, B. Prebiotic, probiotic, antimicrobial, and functional food applications of Bacillus amyloliquefaciens. J. Agric. Food Chem. 2020, 68, 14709–14727. [Google Scholar] [CrossRef] [PubMed]
  13. Casado, V.; Martín, D.; Torres, C.; Reglero, G. Phospholipases in food industry: A Review. Methods Mol. Biol. 2012, 861, 495–523. [Google Scholar] [CrossRef] [PubMed]
  14. Moran, K.; de Lange, C.F.M.; Ferket, P.; Fellner, V.; Wilcock, P.; van Heugten, E. Enzyme supplementation to improve the nutritional value of fibrous feed ingredients in swine diets fed in dry or liquid form. J. Anim. Sci. 2016, 94, 1031–1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Gómez-Gallego, C.; Pohl, S.; Salminen, S.; De Vos, W.M.; Kneifel, W. Akkermansia muciniphila: A novel functional microbe with probiotic properties. Benef. Microbes 2016, 7, 571–584. [Google Scholar] [CrossRef]
  16. Brodmann, T.; Endo, A.; Gueimonde, M.; Vinderola, G.; Kneifel, W.; de Vos, W.M.; Salminen, S.; Gómez-Gallego, C. Safety of novel microbes for human consumption: Practical examples of assessment in the European Union. Front. Microbiol. 2017, 8, 1725. [Google Scholar] [CrossRef] [PubMed]
  17. EFSA Panel on Biological Hazards (BIOHAZ); Koutsoumanis, K.; Allende, A.; Alvarez-Ordóñez, A.; Bolton, D.; Bover-Cid, S.; Chemaly, M.; Davies, R.; De Cesare, A.; Hilbert, F.; et al. Scientific opinion on the update of the list of QPS-recommended biological agents intentionally added to food or feed as notified to EFSA (2017–2019). EFSA J. 2020, 18, e05966. [Google Scholar] [CrossRef] [Green Version]
  18. EFSA Panel on Additives and Products or Substances used in Animal Feed (FEEDAP); Bampidis, V.; Azimonti, G.; de Bastos, M.L.; Christensen, H.; Dusemund, B.; Kouba, M.; Fašmon Durjava, M.; López-Alonso, M.; López Puente, S.; et al. Safety and efficacy of Bacillus subtilis PB6 (Bacillus velezensis ATCC PTA-6737) as a feed additive for chickens for fattening, chickens reared for laying, minor poultry species (except for laying purposes), ornamental, sporting and game birds. EFSA J. 2020, 18, e06280. [Google Scholar] [CrossRef]
  19. EFSA Panel on Additives and Products or Substances used in Animal Feed (FEEDAP); Rychen, G.; Aquilina, G.; Azimonti, G.; Bampidis, V.; de Bastos, M.L.; Bories, G.; Chesson, A.; Cocconcelli, P.S.; Flachowsky, G.; et al. Guidance on the characterisation of microorganisms used as feed additives or as production organisms. EFSA J. 2018, 16, e05206. [Google Scholar] [CrossRef]
  20. Fan, B.; Blom, J.; Klenk, H.-P.; Borriss, R. Bacillus amyloliquefaciens, Bacillus velezensis, and Bacillus siamensis Form an “Operational group B. amyloliquefaciens” within the B. subtilis species complex. Front. Microbiol. 2017, 8, 22. [Google Scholar] [CrossRef] [Green Version]
  21. Yamanaka, H.; Moriyoshi, K.; Ohmoto, T.; Ohe, T.; Sakai, K. Degradation of Bisphenol A by Bacillus pumilus isolated from Kimchi, a traditionally fermented food. Appl. Biochem. Biotechnol. 2007, 136, 39–51. [Google Scholar] [CrossRef] [PubMed]
  22. Li, G.; Zu, L.; Wong, P.-K.; Hui, X.; Lu, Y.; Xiong, J.; An, T. Biodegradation and detoxification of Bisphenol A with one newly-isolated strain Bacillus sp. GZB: Kinetics, mechanism and estrogenic transition. Bioresour. Technol. 2012, 114, 224–230. [Google Scholar] [CrossRef] [PubMed]
  23. Das, R.; Liang, Z.; Li, G.; Mai, B.; An, T. Genome sequence of a spore-laccase forming, BPA-degrading Bacillus sp. GZB isolated from an electronic-waste recycling site reveals insights into BPA degradation pathways. Arch. Microbiol. 2019, 201, 623–638. [Google Scholar] [CrossRef] [PubMed]
  24. Jandhyala, S.M.; Talukdar, R.; Subramanyam, C.; Vuyyuru, H.; Sasikala, M.; Nageshwar Reddy, D. Role of the normal gut microbiota. World J. Gastroenterol. 2015, 21, 8787–8803. [Google Scholar] [CrossRef] [PubMed]
  25. Indiani, C.; Rizzardi, K.F.; Castelo, P.M.; Ferraz, L.F.C.; Darrieux, M.; Parisotto, T.M. Childhood Obesity and Firmicutes/Bacteroidetes ratio in the gut microbiota: A systematic review. Child Obes. 2018, 14, 501–509. [Google Scholar] [CrossRef] [PubMed]
  26. Andra, S.S.; Austin, C.; Yang, J.; Patel, D.; Arora, M. Recent advances in simultaneous analysis of Bisphenol A and its conjugates in human matrices: Exposure biomarker perspectives. Sci. Total Environ. 2016, 572, 770–781. [Google Scholar] [CrossRef] [Green Version]
  27. Cao, G.T.; Dai, B.; Wang, K.L.; Yan, Y.; Xu, Y.L.; Wang, Y.X.; Yang, C.M. Bacillus licheniformis, a potential probiotic, inhibits Obesity by modulating colonic microflora in C57BL/6J mice model. J. Appl. Microbiol. 2019, 127, 880–888. [Google Scholar] [CrossRef]
  28. Wang, Y.; Wu, Y.; Wang, B.; Xu, H.; Mei, X.; Xu, X.; Zhang, X.; Ni, J.; Li, W. Bacillus amyloliquefaciens SC06 protects mice against High-Fat Diet-Induced Obesity and liver injury via regulating host metabolism and gut microbiota. Front. Microbiol. 2019, 10. [Google Scholar] [CrossRef]
  29. López-Moreno, A.; Suárez, A.; Avanzi, C.; Monteoliva-Sánchez, M.; Aguilera, M. Probiotic strains and intervention total doses for modulating Obesity-related microbiota dysbiosis: A systematic review and meta-analysis. Nutrients 2020, 12, 1912. [Google Scholar] [CrossRef]
  30. Pariza, M.W.; Gillies, K.O.; Kraak-Ripple, S.F.; Leyer, G.; Smith, A.B. Determining the safety of microbial cultures for consumption by humans and animals. Regul. Toxicol. Pharmacol. 2015, 73, 164–171. [Google Scholar] [CrossRef]
  31. Thomson, N.; Summers, D.; Sivaniah, E. Synthesis, properties and uses of bacterial storage lipid granules as naturally occurring nanoparticles. Soft Matter 2010, 6, 4045–4057. [Google Scholar] [CrossRef]
  32. Zhang, K.; Zhou, C.; Zhou, Z.; Li, X.; Li, Z.; Zhang, M.; Zhang, X.; Zhang, C.; Wei, T.; Zhang, S.; et al. Identification of lipid droplets in gut microbiota. bioRxiv 2020. [Google Scholar] [CrossRef]
  33. Devaraj, K.; Aathika, S.; Periyasamy, K.; Manickam Periyaraman, P.; Palaniyandi, S.; Subramanian, S. Production of thermostable multiple enzymes from Bacillus amyloliquefaciens KUB29. Nat. Prod. Res. 2019, 33, 1674–1677. [Google Scholar] [CrossRef]
  34. Deb, P.; Talukdar, S.A.; Mohsina, K.; Sarker, P.K.; Sayem, S.A. Production and partial characterization of extracellular amylase enzyme from Bacillus amyloliquefaciens P-001. Springerplus 2013, 2, 154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Latorre, J.D.; Hernandez-Velasco, X.; Wolfenden, R.E.; Vicente, J.L.; Wolfenden, A.D.; Menconi, A.; Bielke, L.R.; Hargis, B.M.; Tellez, G. Evaluation and selection of Bacillus species based on enzyme production, antimicrobial activity, and biofilm synthesis as direct-fed microbial candidates for poultry. Front. Vet. Sci. 2016, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ohore, O.E.; Zhang, S. Endocrine Disrupting Effects of Bisphenol A exposure and recent advances on its removal by water treatment systems: A Review. Sci. Afr. 2019, 5, e00135. [Google Scholar] [CrossRef]
  37. Sasaki, M.; Tsuchido, T.; Matsumura, Y. Molecular cloning and characterization of cytochrome P450 and ferredoxin genes involved in Bisphenol A degradation in Sphingomonas bisphenolicum strain AO1. J. Appl. Microbiol. 2008, 105, 1158–1169. [Google Scholar] [CrossRef]
  38. Nakamura, S.; Tezuka, Y.; Ushiyama, A.; Kawashima, C.; Kitagawara, Y.; Takahashi, K.; Ohta, S.; Mashino, T. Ipso Substitution of Bisphenol A catalyzed by microsomal cytochrome P450 and enhancement of estrogenic activity. Toxicol. Lett. 2011, 203, 92–95. [Google Scholar] [CrossRef] [PubMed]
  39. Lobos, J.H.; Leib, T.K.; Su, T.M. Biodegradation of Bisphenol A and other bisphenols by a gram-negative aerobic bacterium. Appl. Environ. Microbiol. 1992, 58, 1823–1831. [Google Scholar] [CrossRef] [Green Version]
  40. Sasaki, M.; Akahira, A.; Oshiman, K.; Tsuchido, T.; Matsumura, Y. Purification of Cytochrome P450 and ferredoxin, involved in Bisphenol A degradation, from Sphingomonas sp. strain AO1. Appl. Environ. Microbiol. 2005, 71, 8024–8030. [Google Scholar] [CrossRef] [Green Version]
  41. Menasria, T.; Aguilera, M.; Hocine, H.; Benammar, L.; Ayachi, A.; Si Bachir, A.; Dekak, A.; Monteoliva-Sánchez, M. Diversity and bioprospecting of extremely halophilic archaea isolated from Algerian arid and semi-arid wetland ecosystems for halophilic-active hydrolytic enzymes. Microbiol. Res. 2018, 207, 289–298. [Google Scholar] [CrossRef]
  42. Menasria, T.; Monteoliva-Sánchez, M.; Benammar, L.; Benhadj, M.; Ayachi, A.; Hacène, H.; Gonzalez-Paredes, A.; Aguilera, M. Culturable halophilic bacteria inhabiting algerian saline ecosystems: A source of promising features and potentialities. World J. Microbiol. Biotechnol. 2019, 35, 132. [Google Scholar] [CrossRef] [PubMed]
  43. EFSA Panel on additives and products or substances used in animal feed (FEEDAP) 2012; guidance on the assessment of bacterial susceptibility to antimicrobials of human and veterinary importance (question no EFSA-Q-2011-01108, adopted: 23 May 2012 by European Food Safety Authority). EFSA J. 2012, 10, 10. [CrossRef]
  44. Humphries, R.M.; Ambler, J.; Mitchell, S.L.; Castanheira, M.; Dingle, T.; Hindler, J.A.; Koeth, L.; Sei, K. CLSI Methods development and standardization working group best practices for evaluation of antimicrobial susceptibility tests. J. Clin. Microbiol. 2018, 56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. García-Córcoles, M.T.; Cipa, M.; Rodríguez-Gómez, R.; Rivas, A.; Olea-Serrano, F.; Vílchez, J.L.; Zafra-Gómez, A. Determination of bisphenols with estrogenic activity in plastic packaged baby food samples using solid-liquid extraction and clean-up with dispersive sorbents followed by gas chromatography tandem mass spectrometry analysis. Talanta 2018, 178, 441–448. [Google Scholar] [CrossRef]
  46. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 2017, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Jolley, K.A.; Bliss, C.M.; Bennett, J.S.; Bratcher, H.B.; Brehony, C.; Colles, F.M.; Wimalarathna, H.; Harrison, O.B.; Sheppard, S.K.; Cody, A.J.; et al. Ribosomal multilocus sequence typing: Universal characterization of bacteria from domain to strain. Microbiology (Reading) 2012, 158 (Pt 4), 1005–1015. [Google Scholar] [CrossRef]
  48. Yoon, S.-H.; Ha, S.-M.; Lim, J.; Kwon, S.; Chun, J. A Large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 2017, 110, 1281–1286. [Google Scholar] [CrossRef]
  49. Meier-Kolthoff, J.P.; Auch, A.F.; Klenk, H.-P.; Göker, M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 2013, 14, 60. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A phylogenomic tree based on partial 16S RNA gene (Table S1) of isolated strains of Bacillus sp. and closely related species. The tree was inferred by the maximum-likelihood method. The species and strain names are shown, and accession numbers are indicated in parentheses. Bootstrap values after 1000 resamplings are indicated at branch nodes. Bar: 0.0020 nucleotide substitution per site.
Figure 1. A phylogenomic tree based on partial 16S RNA gene (Table S1) of isolated strains of Bacillus sp. and closely related species. The tree was inferred by the maximum-likelihood method. The species and strain names are shown, and accession numbers are indicated in parentheses. Bootstrap values after 1000 resamplings are indicated at branch nodes. Bar: 0.0020 nucleotide substitution per site.
Ijms 22 04952 g001
Figure 2. (a) Colony morphology of Bacillus sp. in MRS medium; (b) Image of Bacillus sp. AM1 by TEM showing lipidic PHA granules, EPS, and Endocrine disruptor ED-BPA. 1µm scale; (c) BPA degradation curve by Bacillus sp. AM1. BPA biodegradation capacity was tested in TSB liquid cultures of Bacillus sp. AM1 exposed to 25 ppb concentration of BPA at 30ºC for 72h (LC-MS/MS system for quantification).
Figure 2. (a) Colony morphology of Bacillus sp. in MRS medium; (b) Image of Bacillus sp. AM1 by TEM showing lipidic PHA granules, EPS, and Endocrine disruptor ED-BPA. 1µm scale; (c) BPA degradation curve by Bacillus sp. AM1. BPA biodegradation capacity was tested in TSB liquid cultures of Bacillus sp. AM1 exposed to 25 ppb concentration of BPA at 30ºC for 72h (LC-MS/MS system for quantification).
Ijms 22 04952 g002
Figure 3. Phylogenomic trees based on gene sequences of Bacillus sp. AM1: (a) 16S RNA gene; (b) gyrA gene, similarity distances were computed using JTT’s model; (c) rpoB gene, similarity distances were computed using JTT’s model. The tree was inferred by the maximum-likelihood method. The species and strain names are shown, and accession numbers are indicated in parentheses. Bootstrap values shown after 1000 resamplings.
Figure 3. Phylogenomic trees based on gene sequences of Bacillus sp. AM1: (a) 16S RNA gene; (b) gyrA gene, similarity distances were computed using JTT’s model; (c) rpoB gene, similarity distances were computed using JTT’s model. The tree was inferred by the maximum-likelihood method. The species and strain names are shown, and accession numbers are indicated in parentheses. Bootstrap values shown after 1000 resamplings.
Ijms 22 04952 g003
Figure 4. Proposed Biodegradation Pathways of BPA based on the WGS analysis of Bacillus sp. AM1: Reaction steps, Enzymes, EC number, protein ID, and specific genes loci.
Figure 4. Proposed Biodegradation Pathways of BPA based on the WGS analysis of Bacillus sp. AM1: Reaction steps, Enzymes, EC number, protein ID, and specific genes loci.
Ijms 22 04952 g004
Table 1. Results derived from the enzymatic analysis of Bacillus sp. AM1.
Table 1. Results derived from the enzymatic analysis of Bacillus sp. AM1.
Enzyme TestMicroorganisms
Bacillus sp. AM1B. amyloliquefaciens (DSM7)B. velezensis (CR-502)B. siamensis (KCTC13613) B.cereus (B7)B.mycoides
(B12)
Starch +++++++++
Carboxymethylcelulose
Gelatine
Inuline+++
Tween 80 ++
DNase++++
Table 2. Antibiotic sensitivity testing.
Table 2. Antibiotic sensitivity testing.
AntibioticConcentration
(μg)
Resistance (R)
Sensitive (S)
Zone inhibition
(mm)
Ampicillin10R-
Amikacin30S14
Ceftazidime30R-
Ciprofloxacin5S40
Chloramphenicol30S42
Clindamycin2S22
Erythromycin15S25
Gentamicin10S22
Imipenem10S50
Penicillin10R-
Rifampicin30S29
Vancomycin30S22
Kanamycin30S19
Streptomycin10R-
Tetracycline30S17
Co-Trimoxazole25S39
Amoxicillin/
Clavulanic acid
20/10R-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

López-Moreno, A.; Torres-Sánchez, A.; Acuña, I.; Suárez, A.; Aguilera, M. Representative Bacillus sp. AM1 from Gut Microbiota Harbor Versatile Molecular Pathways for Bisphenol A Biodegradation. Int. J. Mol. Sci. 2021, 22, 4952. https://doi.org/10.3390/ijms22094952

AMA Style

López-Moreno A, Torres-Sánchez A, Acuña I, Suárez A, Aguilera M. Representative Bacillus sp. AM1 from Gut Microbiota Harbor Versatile Molecular Pathways for Bisphenol A Biodegradation. International Journal of Molecular Sciences. 2021; 22(9):4952. https://doi.org/10.3390/ijms22094952

Chicago/Turabian Style

López-Moreno, Ana, Alfonso Torres-Sánchez, Inmaculada Acuña, Antonio Suárez, and Margarita Aguilera. 2021. "Representative Bacillus sp. AM1 from Gut Microbiota Harbor Versatile Molecular Pathways for Bisphenol A Biodegradation" International Journal of Molecular Sciences 22, no. 9: 4952. https://doi.org/10.3390/ijms22094952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop