Microbiome of highly polluted coal mine drainage from Onyeama, Nigeria, and its potential for sequestrating toxic heavy metals
Geochemistry and ecotoxicology of AMDAMD systems are an important source of metal/metalloid pollution to the receiving hydrosphere with devastating consequences on the biological drivers of affected ecosystems. Environmental menaces of AMD have not been exhaustively reported worldwide. Scanty information exists across Africa and many developing economies. The homogenised mixture of detached biofilm and AMD samples from a derelict coal mine at three sampling periods were assayed for geochemical delineation and analysed for pollution intensity against reference background geochemical values. The measured values of the physical properties and contents of selected HMs in drains from a coal mine in Nigeria were as presented in Supplementary Table A.1. Virtually all the measured parameters exceeded the permissible limits of WHO guidelines for potable water. The AMD water was acidic (pH = 3.1 ± 0.265), and contained characteristic anions that are common to AMD including dissolved sulphides (1.37 ± 0.233 mg l−1), sulphates (313.0 ± 15.9 mg l−1), carbonate (253.0 ± 22.4 mg l−1) and nitrate (86.6 ± 41.0 mg l−1) above the allowable limits of WHO. Although the acidic pH of AMD in the present study compares well with those associated with mines in Russia14, more extreme acidic pH values have been reported in other climes. Negative pH values of − 1.56 and − 3.6 were observed in AMD from Iberian Pyrite Belt20 and Richmond Mine at Iron Mountain, USA21, respectively. The values of physicochemical parameters associated with the AMD from Onyeama were similar to data reported for other mine wastewaters in Nigeria22 and elsewhere4. It is known that sulphide minerals, in presence of water and oxygen, oxidise to sulphate as observed in the elevated sulphate concentration (313 ± 15.9 mg l−1) in the present study. The low pH observed in the AMD is due to the formation of sulphuric acid from sulphate in presence of protons (H+). This consequently causes the leaching of metal/metalloid ions into the drains. The concentrations of dissolved organic matter in AMD tends to be relatively low ( Co > Pb > As > Ni > Cr > Fe (Table 1). Enrichment of five HMs was exceptionally high (Cd > Co > Pb > As > Ni), while Cr and Fe were very high and moderately enriched the AMD water, respectively. The astronomically high contamination and enrichment factors of the AMD signified the enrichment potentials the AMD portends on receiving surface waters. The AMD from the Onyeama coal mine has been reportedly impacting the water qualities of rivers within the location25. It is assumed that the extremely high concentrations of toxic metals/metalloids in the AMD dilutes out upon discharges into nearby rivers, contaminating the surface water and raising the bioavailable metals/metalloids beyond safe thresholds. Further reports of toxic metals/metalloids enrichment of surface waters via inflow of AMDs from other mines in Nigeria26 and other climes3,4,27 are worrisome and oblige mitigations.Table 1 Physico-chemistry, pollution and ecological impact determinants of heavy metals and metalloid contained in the AMD from coal mine.Full size tableThe HMs-enriched environments inadvertently exert ecotoxicity unto the drivers of the ecosystems. The level of HMs accumulation to the organic matter in the AMD, through geo-accumulation (Igeo) index of Fe (7.60 ± 0.779) to Cd (20.9 ± 0.075) (Supplementary Table A.2), was very severe and in a similar order to CF. It possibly implies organic matter in the AMD harbours the mobile toxic metal/metalloid concentrations and make them available to the food web28. Thus, biomagnification of the toxic metals/metalloids along the trophic level becomes palpable and a challenge to the biota of any surface water receiving the AMD and to public health21,28. Ecological risk assessments define and categorise the pollution status of ecosystems with the HMs contained in the AMD. Based on the potential ecological risk factor (Er), Cd exerted an extremely high-risk index (36.3 ± 1.96 × 106), and none of the metals/metalloids exercised less than 1000 risk index (Supplementary Table A.2). All the HMs/metalloid contained in the AMD posed very high ecological risks and could be categorised in the order of Cd > Co > Pb > As > Ni > Cr > Fe. The modified potential ecological risk factor (MEr), however, stipulated that five HMs posed a very high risk in the order: Cd > Co > Pb > As > Ni, whereas Cr and Fe were determined to be of considerate and low risks, respectively. The HMs exerted high risk to the AMD ecosystem as calculated by ecological risk quotient (RQ) in the order: Pb > Cd > As > Ni > Co > Fe > Cr. The ecological risk index of all the HMs as a whole was very high (375,000 ± 22,400) index as stipulated by the modified potential ecological risk index (Table 1). The prodigiously high ecological risks indexes of the HMs/metalloid in the AMD indicated grave danger the AMD would portend on the surface- and ground-waters.Microbial community structure of AMD from Onyeama coal mineA total of 26,160 and 40,403 valid (filtered) sequence reads were obtained for bacteria and eukarya, respectively, after a quality check of biofilm-water amplicon sequence data. The valid sequences were clustered into 2036 and 1002 operational taxonomic units (OTUs) of bacteria and eukarya domains of life, respectively, as presented in Table 2. Microbial community structures are sensitive descriptors of ecological stressors pivotal to understanding ecosystem functions29. The number of clustered high quality, non-chimeric sequences as OTUs based on UCLUST and CD-HIT against the sequence reads was depicted as asymptotic rarefaction curves (Supplementary Fig. A.1). The curves revealed that higher numbers of OTUs were delineated from valid sequence reads of 16S rRNA genes, unlike the lesser number of OTUs obtained from valid sequence reads of ITS2 region located between 5.8S and 28S rRNA genes of eukaryotes. The OTU richness observed in the rarefaction curves established coverage of the majority of species and was further validated with the richness and diversity estimations presented in Table 2. Despite the higher number of valid sequence reads obtained from the amplified ITS2 (40,403) than that of 16S rRNA genes (26,160), the observed OTUs were more in 16S rRNA genes (2036) than those of ITS2 (1002). More than 99.8% and about 98.5% of the microbial community in AMD from the Onyeama coal mine represented eukarya and bacteria OTUs, respectively, based on estimated Good’s library coverage. The coverage degree of the MiSeq sequencing corroborated the rarefaction curves. Furthermore, the estimated OTU richness (based on higher values obtained from ACE, Chao1 and JackKnife indexes) showed that bacterial phylotypes were richer than those of eukarya. Alpha diversity indexes (NPShannon, Shannon, and inverse Simpson) phylogenetic diversity index revealed that bacteria in the AMD were more diverse than eukarya OTUs.Table 2 Alpha diversity of microbiome evenness, richness and varieties of species in the sediments.Full size tableTaxonomy and phylogeny of microbial OTUs in AMD from coal mineThe taxonomic composition and relative abundances of the AMD microbiome, as shown in Fig. 1, revealed that the bacterial community spanned 10 phyla whose sequence reads were at least 1% (Fig. 1a). Whereas the eukarya domain of life (with sequence reads ≥ 1%) found in the AMD include Fungi, Plantae and Animalia kingdoms (Fig. 1b). Ascomycota, unclassified Fungi phylum (Fungi_p), Basidiomycota, and Mucoromycota represented Fungi kingdom, while Ciliophora and Arthropoda phyla were Animalia and Chlorophyta phylum epitomised Plantae kingdom. Association of the domain Eukarya (comprising Alveolates, Chlorophyta and Fungi as observed in this study) with AMD is reported to a lesser extent when compared with Bacteria30. The Fungi, largely represented by Ascomycota and Basidiomycota, are primarily found in sub-surface low-pH biofilms thriving in AMD31. While the Alveolates are suggested to have acted as primary/secondary consumers, the amoebae were secondary grazers in the AMD ecosystem29,32. Fungi taxa must have participated in carbon cycling as the main decomposers in the microbial community of the AMD. The taxonomic composition and relative abundance of phyla regarded as ‘Others’ (sequence reads 50%). Evolutionary analyses were conducted in MEGA6.Full size imageUrease-producing bacteria instigate insoluble metal-carbonate micro-precipitation through urease activity16. The growth-time courses and urease activities of the bacteria consortium in simulated AMD were presented as curves (Fig. 5). It was observed that the impact of high concentrations of HMs cocktails was not pronounced beyond the early 6 h post-inoculation, which was regarded as the lag phase. The bacteria consortium might have activated necessary genes needed to tolerate and sequester the metals/metalloids toxicity during the lag phase without cell multiplications. Afterwards, the bacteria consortium grew steadily with the production of urease, based on increasing measurement of urease activity, as incubation continued. At 30 h post-inoculation, 245.3 (± 23.7) U ml−1 activity of urease was observed in broth without a toxic metal cocktail. However, more urease activity (255 ± 7.6 U ml−1) by the bacteria consortium was observed in medium amended with low concentrations of metal cocktails unlike lesser activities of 235 (± 7.6) U ml−1 and 193.7 (± 10.7) U ml−1 associated with medium and high metal concentrations, respectively. As the growth remains stationary and pH further increased to > 8.2, urease activities were at least 253 U ml−1 in all the cultures. Although urease activities at acidic pH have been reported in acid-tolerant human pathogens19, the findings in this report were assumedly the first amongst bacterial strains from AMD-polluted environments. The urease activities at acidic pH compared favourably with activities at alkaline pH in previous studies7,16,42,44. Moreover, the pH of the culture system kept increasing, alleviating the acidity condition that initially prevailed in the AMD system.Figure 5Growth kinetics of bacterial consortium via viable counts extrapolated into optical density at 600 nm wavelength (a) and growth-dependent urease activity of bacterial consortium (b) in TGYM broth without heavy metals (HMs) cocktail, and with low, medium, and high concentrations of HMs cocktails. Low HMs concentrations cocktail comprised (per liter) Cd, 27.9 mg; Pb, 118.7 mg; Co, 16.2 mg; Ni, 16.2 mg; and As, 61.5 mg. While medium HMs concentration contained (per liter) Cd, 55.7 mg; Pb, 237.3 mg; Co, 32.4 mg; Ni, 32.3 mg; and As, 123.1 mg. High HMs concentration contained (per liter) Cd, 139.3 mg; Pb, 593.3 mg; Co, 81.1 mg; Ni, 80.7 mg; and As, 307.6 mg. The mean pH at the beginning of experiment was 3.5 and rose to 8.2–8.4 at 48 h post-inoculation. Growth kinetics at exponential growth phase are in the inserts of panel (a), where ‘Td’ represents doubling time and ‘K’ is the growth rate at exponential growth phase. Error bars represent standard error mean (SEM) of triplicate experiments. The culture conditions were as explained in the “Methods” Section (Growth kinetics and urease activity of bacteria consortium; Determination of bacterial growth-dependent HMs/metalloid sequestration in simulated and natural AMD).Full size imageInterestingly, urease activity was observed in low quantity at acidic pH, unlike higher activity when the pH inclined towards alkaline (Fig. 5). It is proposed that urea finds its way into Onyeama coal mine drains through runoff from agricultural soils fortified with urea fertilizers and animal manures, which are common agricultural practices in Nigeria. The products of urea hydrolysis might have equilibrated in water to form bicarbonate, ammonium and hydroxyl ions that serially increased the culture pH. Ultimately, the bicarbonate equilibrium might have shifted to form carbonate ions (HCO3− + H+ + 2NH4+ + 2OH− ↔ CO32− + NH4+ + 2H2O) that enhanced the metal-carbonate micro-precipitation (Me2+ + Cell → Cell-Me2+ + CO32− → Cell-MeCO3). The gradual increase in pH could have further indulged the formation of CO32− from HCO3−, leading to metal-CO3 precipitation around cells and in culture media. Bicarbonates enrichment with inherent ammonia production was thought to have provided additional acid neutralization of the AMD. The growth kinetics after the presumed lag phase in the early 6 h to late exponential phase at 18 h showed that a low concentration of HMs cocktails did not have an impact on the growth of the bacteria consortium. Consequently, the bacteria consortium exhibited excellent sequestration of multi-component toxic HMs in both the simulated toxic metal-rich AMD and the actual AMD obtained from the Onyeama coal mine (Table 3).Table 3 Growth associated sequestration and precipitation of heavy metals/metalloid cocktail and AMD from Onyeama coal mine.Full size tableThe bacteria consortium displayed more than 94% efficiency of Cd and Pb sequestration in natural AMD, while 100% efficiency was observed in all the simulated AMD treatments (Table 3). Low performance was found with Ni and As, but not less than 70% sequestration efficiency was observed in all treatments. Efficient sequestrations of toxic metals, up to 100% removal efficiency of most toxic metals, observed with the bacteria consortium were similar to findings in a previous study13. Mixed-bacterial cultures are known to be able to perform more complex tasks and survive in more unstable environments than a monoculture. Nevertheless, 89.3–98% removal efficiencies of Ni, Pb, Co, and Cd from solution have been reportedly achievable with urease-producing Sporosarcina koreensis45. Similarly, Bacillus sp. KK1 reportedly mitigated lead-contaminated mines tailings containing mobile Pb (1050 mg kg−1) to form insoluble precipitates of PbS and PbSiO334. Growth-dependent sequestration of HMs cocktails by the bacteria consortium was adduced to be via precipitation. The weight of the precipitates was evaluated to be proportional to concentrations of HMs cocktail present. The bacteria consortium was observed to drive the formation of as much as 15.6 (± 0.92) mg ml−1 precipitates (Table 3) that were assumed to be in form of HMs-carbonates in TGYM supplemented with high concentrations of HMs cocktail within 24 h post-inoculation. In natural AMD bio-stimulated with urea and seeded with bacteria consortium for 24 h, 10.5 (± 0.52) mg ml−1 HMs precipitates was observed unlike 8.57 (± 2.52) mg ml−1 precipitates obtained from natural AMD toxic metals sequestration without urea fortification. It appeared that the quantity of toxic metal precipitate was proportional to quantities of available toxic metals, which corresponded to the number of heterogeneous nucleation sites on the surface of the bacterial cells. Omoregie et al.42 reported a relatively similar quantum of precipitation as CaCO3 with species of ureolytic Firmicutes isolated from limestone caves. As such, there was no correlation between urease activity and quantum of toxic metal precipitation since there is a likelihood that other metabolic activities may be linked to urease activities. Nevertheless, the bioremediation strategies demonstrated in the present study exhibited excellent toxic metal sequestrations unlike insignificant (p > 0.05) natural attenuation process of the autochthonous community without augmentation with bacteria consortium and stimulation with nutrients (as presented in Table 3).In conclusion, AMD from the Onyeama coal mine is a point source of pollution to the surrounding environments because of its richness in anions and toxic metals/metalloids. It has a high potential of enriching the receiving hydrosphere with toxic metals/metalloids and exerts severe ecological risks (Er > 320) with Cd and Pb wielding a huge critical risk index (38.1 ± 2.18 × 106) on the biological elements of the ecosystems. The dominance of Proteobacteria (50.8%), Bacteroidetes (18.9%), Ascomycota (60.8%), and Ciliophora (12.6%) characterised the microbial community of the AMD, where unclassified OTUs occurred mostly among the species. Enrichment of the AMDs skewed the bacterial community as depicted in the alpha diversity indexes against that of coal AMD leading to the selection of bacteria consortium with an excellent potential of stemming the toxicants in the AMD. The bacteria consortium efficiently removed toxic metals/metalloids ( > 70%) through precipitation and simultaneously neutralised AMD acidity. The bacteria consortium exhibited appreciable urease activity ( > 190 U ml−1), through which the precipitation was assumed possible via the formation of metal/metalloid-carbonates. The bacteria consortium is suggested to be a sustainable biotechnological candidate in designing a bioremediation strategy for decommissioning AMD before discharge into the surrounding environment. More