3.1. Physiochemistry, Cell Morphology and Activity of Epipelon Biofilm
The physicochemical properties at the sampling site of epipelon biofilm were: ammonia 0.66 mg L
−1, nitrate 0.78 mg L
−1, total phosphorus (TP) 0.23 mg L
−1, total nitrogen (TN) 2.12 mg L
−1 and pH 7.78. Epipelon are an integration of a wide variety of microorganisms, including algae, bacteria and protozoa, that collaboratively live in or on fine-grained substrata or sand. However, these biofilms under ambient methane concentration were predominantly composed of cyanobacteria and phytoplankton (a). Moreover, after incubation under 120,000 ppm methane, bacterial species outnumbered cyanobacteria and phytoplankton as shown by PCM (b). SEM showed that the structure of epipelon biofilm was intertwined bundles, with voids of algae, and diatoms, while bacteria attached to these were evident in these structures. The spaces and tunnels found within epipelonic biofilms are likely to provide physical space for bacterial attachment and metals sorption sites and are clearly visible in the SEM (c,d) [
41].
Figure 1. Phase Contrast Microscopy (PCM) and Scanning Electron microscopy (SEM) were performed to differentiate microbial community change (a) and (b) community composition of epipelon observed under PCM before and after 12% methane treatment, respectively. (c) and (d) structure of periphytic biofilm observed under SEM before and after 12% methane treatment, respectively. (e) The average well color development (AWCD) of the Biolog EcoPlates at 590 nm for periphytic biofilm and diversity indices (E) for species diversity.
The capability of different types of periphytic biofilms in the bioremediation and biosorption of different contaminants from physiochemically contrasting lakes has been intensively reported [
42,
43]. However, no studies have yet been reported about enhanced heavy metals removal by epipelon under a high methane regime which could likely represent naturally occurring conditions in a wetland.
In Biolog® experiments, the Average Well Color Development (AWCD) is an indicator of the degree of activity of epipelon biofilm. The darker color of the wells represents the higher activity of the microorganisms as confirmed from the slope of the graph (e). The diversity indices showed higher community abundance, higher functional activity of epipelon, uniformity and species eveneness of the microbial community (e), e.g., after incubation of Biolog® for 7 days, the Shannon index was found to be 3.6. This result indicates higher microbial diversity and species evenness in epipelon when compared to the control.
3.2. Effect of Methane Oxidation on Heavy Metal Bioremediation
Epipelon was found to be highly active in removing Pb
2+ and Cd
2+ under 12% (120,000 ppm) methane, with 100% removal efficiency within 336 h (a), although only 20% of Cr
6+ was removed within same interval of time. It should also be noted that the low removal efficiency of Cr
6+ could not entirely reflect the operating conditions in natural environments because Cr
6+ is toxic in laboratory scale experiments without continuous recirculation, as previously suggested [
44]. Meanwhile, hexavalent chromium is very toxic due to its high water-solubility and strong oxidizing ability. Chromium (VI) posed a toxicity risk to periphytic biofilm via inhibition of its photosynthetic activity and carbon utilization function [
45], thus inducing destruction of the biofilm, which ultimately results in releasing Cr
6+ back into the environment, making it more toxic. For instance, hexavalent chromium removal efficiency decreased from 90% to 50% for
E. crassipes when compared to trivalent chromium [
46]. A free water surface wetland made from sediment was found to be only capable of 7% chromium removal [
47]. Furthermore, the isotope fractionation factor (ε) of Cr
6+ was 2.62 ± 0.20% after 509 days of continuous incubation with methane based on theoretical estimation using the Rayleigh distillation mode [
48] during slow bioremediation of Cr
6+ by methanotrophs. In one similar periphytic study, Cd(II) was more secluded in comparison to Cr(VI), while field validation of this observation remained an intriguing issue for further research [
49].
Figure 2. Methane oxidation and heavy metals bioremediation potential of epipelon (a) the ability of epipelonic biofilm to remove heavy metals at different time intervals coupled with methane oxidation (b) Methane oxidation at different time intervals under heavy metals effect.
It was further observed that under ambient methane concentration (6 ppm), epipelon completely removed amended Pb
2+ or Cd
2+ (both at 50 mg L
−1) within 480 hrs, but the incubation under high methane concentration (120,000 ppm) greatly enhanced the metal removal efficiencies, which took only 336 hrs to complete the process (a). This might be attributed to saturation effect of the active sites on the surface of the epipelonic biofilms. Comparatively, limited sensitivity of methane oxidation was found after adding Pb(II), which might have been due to soil water content as the primary key force for Pb(II) removal [
50].
The patterns of metal tolerance, as shown by methane oxidizers, are comparatively distinct to chloride salts of Cd(II), Cr(III), Hg(II), Cu(II)and Zn(II) [
51]. Another interesting aspect of the novel isolates of methanotrophs was their resistance to various types of pollutions (heavy metals, arsenate or organic) and were considered to be real ‘super-bugs’ due to their capability to withstand remarkably high concentration of various pollutants [
52]. The abundance of extracellular polymeric substances, plays a key role in the higher removal efficacies of periphytic biofilms [
21] as well as in-tank microbial activity [
53].
The potential implication of methane monooxygenases, for bioremediation of different contaminants including heavy metals, is accentuated by their broad substrate profiles [
54]. Hence, epipelonic biofilms have proven to be considerably effectual towards the removal of heavy metals within a shorter time period under high methane doses compared to previous studies, thus providing a reliable source for the removal of heavy metals from aquatic environments.
3.3. Potential of Epipelon for Methane Oxidation and 13C Assimilation
All treatments were found to exhibit strong methane oxidation activity and consumed more than 99% of amended methane except for 6 ppm 13CH4 treatment, which shows no significant methane oxidation, and 50 mg L−1 Cr6+ + 120,000 ppm 13CH4 only utilized 23.94% of available methane (b). No significant difference in methane oxidation rate was observed between 12CH4 and 13CH4 treatment in microcosms with high concentration.
The assimilation of methane derived carbon obtained by methane oxidation was increased in organic carbon in high dose (12%) methane treatments. The percentage of
13C atom was significantly increased from 1.08%
13C atom in the near-atmospheric dose
13C methane (6 ppm) treatment and in high-dose
12C methane treatments to 16.24%, 12.14% and 13.78%
13C atoms in 120,000 ppm
13CH
4, 50 mgL
−1 Cd
2+ + 120,000 ppm
13CH
4 and 50 mgL
−1 Pb
2+ + 120,000 ppm
13CH
4 treatments, respectively (
Figure S1). These results confirm the potential of epipelonic biofilms for the assimilation of methane in higher amounts, thus confirming their potential towards effective methane oxidation in aquatic environments.
3.4. Methanotrophic Community Abundance in Response to Heavy Metals
The abundance of the total microbiome and methanotrophs in epipelon was evaluated by real-time quantitative PCR of 16S rRNA and
pmoA genes, respectively. Heavy metals significantly inhibit the abundance of total microbial communities (d). Compared to control treatment (120,000 ppm
13CH
4), the number of 16S rRNA genes was 37.94% and 26.49% lower in microcosm amended with Cd
2 and Pb
2+, respectively (d). In previous studies, single cultures of methanotrophic bacteria were verified for bioremediation of heavy metals such as chromium [
17,
52,
55]. Although the favorable conditions were revealed for a single type of strain, it is not a cost-effective and sustainable solution for heavy metal removal on a large scale. In the present study, the presence of Cr
6+ in 50 mg L
−1 Cr
6+ + 120,000 ppm
13CH
4 treatment resulted in a lower abundance (6.97 less folds) of methanotrophs due to Cr
6+ toxicity as compared to 120,000 ppm
13CH
4 whereas the treatments with lead (1.20 less folds) and cadmium (1.57 less folds) did not show much effect on the methanotrophic community in biofilms.
Figure 3. Illumina sequencing was carried out at the whole epipelon microbial community level in microcosms targeting total methanotrophs from total 16S rRNA genes (a), and the relative abundances of type I (b), type II (c), The pmoA gene copy numbers were assessed using real-time quantitative PCR (d).
On the other hand, 120,000 ppm 13CH4 got 884, 12.13 and 9.46 fold higher MOB than 6 ppm 13CH4, 50 mg L−1 Cd2++ 120,000 ppm 13CH4, 50 mg L−1 Pb2++ 120,000 ppm 13CH4, respectively (d).
According to high-throughput sequencing of 16S rRNA genes, high methane treatments (120,000 ppm 13CH4, 50 mg L−1 Cd2+ + 120,000 ppm 13CH4 and 50 mg L−1 Pb2++ 120,000 ppm 13CH4 ) got 9.7, 7.35, 7.69 more folds of MOB genes, respectively, compared to near-atmospheric methane treatment (6 ppm 13CH4) (a). Chromium-augmented treatment showed only 1.41 more folds of MOB genes than near atmospheric methane treatment. This might be due to toxic nature of chromium that suppressed the abundance of MOB genes in the following treatments.
The abundance of type I and type II MOB in the whole microbial community was determined using high-throughput sequencing of 16S rRNA genes data. Type I methanotrophs copy numbers were increased by factors of 12.16, 8.39, 8.69 fold in the high methane treatments (120,000 ppm 13CH4, 50 mg L−1 Cd2+ + 120,000 ppm 13CH4 and 50 mg L−1 Pb2+ + 120,000 ppm 13CH4) compared to near-atmospheric methane treatments (6 ppm) (b). Type II MOB was also stimulated by 5.43, 4.02, 3.50 and 1.01 folds in 120,000 ppm 13CH4, 50 mg L−1 Cd2+ + 120,000 ppm 13CH4, 50 mg L−1 Pb2+ + 120,000 ppm 13CH4 and 50 mg L−1 Cr + 120,000 ppm 13CH4, respectively, as compared to 6 ppm 13CH4 (c).
The chromium-augmented treatment showed only a 1.61-fold increase for the high methane as compared to near atmospheric methane treatment (6 ppm). The sequence analysis of the microbial community from lake Biwa, Japan [
3,
56], stratified eutrophic Swiss lakes [
57] and Lacamas lake in Washington, USA [
58] showed Type I methanotroph dominances, similarly to our findings in this study.
3.5. SIP for Identification of Active Methanotrophs
The heavy fraction resolved by isopycnic centrifugation of the DNA was employed for qPCR of
pmoA genes and clearly indicates active cell propagation during the consumption of
13C-methane in all the high methane dose treatments, except for the Cr
6+ amended (10.05%) treatment and ambient
13CH
4 (6 ppm) treatment. A peak shift of relative
pmoA gene abundances (a) towards heavy fractions (buoyant density around 1.745 g mL
−1) was very evident in all high-dose
13CH
4 treatments, while the
12CH
4 pmoA gene abundances peaked only in the light fractions (buoyant density around 1.72 g mL
−1) [
29]. In
13C-labeled heavy DNA fractions (buoyant density around 1.745 g mL
−1), the number of
pmoA genes was significantly increased to 57.73%, 61.79% and 62.62% in 120,000 ppm
13CH
4, 50 mg L
−1 Cd
2+ + 120,000 ppm
13CH
4 and 50 mg L
−1 Pb
2+ + 120,000 ppm
13CH
4, respectively (b) as compared to the non-significant
13CH
4 (6 ppm) treatment.
Figure 4. The enrichment of 13Clabeled methanotrophs on the basis of qPCR of pmoA gene (a). The relative pmoA gene abundance was the proportion of MOB like gene copy numbers in each fraction to the total abundance through the gradient in a treatment (b). The red columns represent the percentage of MOB affiliated reads in total 16S rRNA genes in the heavy fractions (buoyant density 1.745 g mL−1) while the black column represents light fraction (buoyant density 1.72 g mL−1).
3.6. Microbial Community Response to Heavy Metals Bioremediation and Methane Oxidation
The structures and abundances of microbial communities were drastically changed due to the selection of distinct phylotypes in the different heavy metal treatments, along with methane oxidation (
Figure S2). The taxonomic analysis of heavy metal resistant genes revealed that the most common bacterial phyla in the treatments was
Proteobacteria followed by
Actinobacteria [
59] which is in consistent with our study. Microbial community composition and interaction was analyzed by high throughput sequencing to determine the impact of metals and methane oxidation on microbial community structure. The taxonomic classification of OTUs at different levels produced 28 phyla, 92 classes, 178 orders, 334 families and 751 genera. Approximately 28 different phyla were found to be present in sequencing of the epipelonic biofilm samples, among these Proteobacteria, Cyanobacteria, Saccharibacteria and Gemmatimonadetes were the most abundant phyla in all treatments (
Figure S3). The phylum Proteobacteria was found to be the most dominant in 120,000 ppm
13CH
4. In comparison to other treatments, it outnumbered 50 mg L
−1 Pb
2+ +120,000 ppm
13CH
4, 50 mg L
−1 Cd
2+ + 120,000 ppm
13CH
4 and 6 ppm
13CH
4 by 13.66%, 17.23% and 22%, respectively. Alpha and gamma Proteobacteria have recently been shown to be effectual for the bioremediation of different heavy metals such as cadmium in an iron reducing bacterial immobilization system [
60].
On the basis of Spearman correlation analysis, a co-occurrence network of different co-existing communities 120,000 ppm
13CH
4 comprised of 287 nodes, 5128 edges, average degree 1.87, highest graph density 0.182, highest modularity 0.193 and average network distance 2.58 was generated, indicating strong microbial interactions (). Ambient 6 ppm
13CH
4 and 50 mg L
−1 Cr
6+ + 120,000 ppm
13CH
4 treatments were dominated by alpha-Proteobacteria but the scenario was changed to gamma Proteobacteria in high methane dose treatments as represented by colored dots in co-occurrence network. Based on the level of interactions from all the calculated modules, we can arrange treatments in the following decreasing order of interaction: 120,000 ppm
13CH
4 > 50 mg L
−1 Cd
2+ + 120,000 ppm
13CH
4 > 50 mg L
−1 Pb
2+ + 120,000 ppm
13CH
4 > 6 ppm
13CH
4 > 50 mg L
−1 Cr
6+ + 120,000 ppm
13CH
4 (
Table S4). The modularity of all the networks differs from each other, demonstrating the verity that the modularity of different networks with different doses of methane or metals or contaminants changed [
61]. However, the strong correlation between nodes illustrates that the microbe–microbe interaction stayed strong after the intrusion of external sources or by the bioadsorption of heavy metals or methane oxidation.
Figure 5. The co-occurrence network analysis of bacterial communities in the order (a) 6 ppm 13CH4 (b) 120,000 ppm 13CH4 (c) 50 mg L−1 Pb2+ + 120,000 ppm 13CH4 (d) 50 mg L−1 Cd2+ + 120,000 ppm 13CH4 (e) 50 mg L−1 Cr6+ + 120,000 ppm 13CH4. The network represents strong (Spearman’s r > 0.8) and significant (p value < 0.01) correlation. Each node represents a unique sequence in the whole data set, colored nodes (excluding light gray) depicts the ten major groups while light gray represents all other groups.
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) analysis represented functional attributes of all treatments, 120,000 ppm
13CH
4, 50 mg L
−1 Cd
2+ + 120,000 ppm
13CH
4 and 50 mg L
−1 Pb
2+ + 120,000 ppm
13CH
4 was dominated by methanotrophy (
Figure S4) followed by chemohetrotrophy, aerobic chemohetrotrophy, and phototrophy. Pheatmap confirms the presence of genes engaged in numerous processes such as cellular processes, environmental information processing, metabolism, and genetic information processing [
62].
We further investigated the presence of dominant species at phylum, class, order, and family level by linear discriminant analysis effect size pipeline (LEfSe) analysis. Proteobacteria dominated among all treatments while the 50 mg L
−1 Cr
6+ + 120,000 ppm
13CH
4 community was dominated by Cyanobacteria and Acidobacteria (
Figure S5). The microbial community in treatments, 50 mg L
−1 Cr
6+ + 120,000 ppm
13CH
4 was predominately composed of Verrucomicrobia,
Methylomonas and Fimbriimonadales.
3.7. Methanotrophic Community Composition
Taxonomic analysis of 13C-labeled methanotrophic 16S rRNA genes further demonstrated distinct proportions of active type I and II methane oxidizers in all treatments (b). Particularly, Methylobacter a type I methanotroph was 51.81%, 48.09%, 49.59% abundant in 50 mg L−1 Pb2+ + 120,000 ppm 13CH4, 50 mg L−1 Cd2+ + 120,000 ppm 13CH4, 120,000 ppm 13CH4 treatments, respectively (a). On the contrary, Methylomicrobium was 80.04% abundant in 50 mg L−1 Cr6+ + 120,000 ppm 13CH4. These results were further confirmed by neighbor-joining tree-clustering around Methylobacter spp. (c).
Figure 6. Methanotrophic community compositions of 13C-labeled MOB based on phylogenetic tree analysis of 16S rRNA and pmoA genes. (a) Percentage of type Ia, Ib, II and IIa methanotrophs were computed from the 13C-labeled genes. (b) Percentages of type I and II MOB were determined from16S rRNA high throughput sequencing of the 13C-labeled genes (c) Neighbor-joining phylogenetic tree depicts the phylogenetic relationship of dominant OTU sequences with the known MOB sequences from GenBank. The percentages illustrate the relative sequence abundance of each OTU to the total MOB-affiliated 16S rRNA or pmoA genes. Bootstrap values higher than 60% are pointed at the branch nodes.
Thus, it was demonstrated that the abundance of Methylobacter sp. by the 13C-labeled pmoA genes. Type I outnumbered Type II with percentages of 90.65%, 90.91% and 99.80% of total methanotrophs in 50 mg L−1 Pb2+ + 120,000 ppm 13CH4, 50 mg L−1 Cd2+ + 120,000 ppm 13CH4 and 120,000 ppm 13CH4, respectively. Type II MOBs, Methylosinus and Methylocystis, were stimulated by high methane dose and metal contamination in treatments 50 mg L−1 Pb2+ + 120,000 ppm 13CH4 and 50 mg L−1 Cd2+ + 120,000 ppm 13CH4 by 9.14% and 8.8% as compared to 120,000 ppm 13CH4.
Methylobacter spp. and
Rhodoferax spp. were found to be the most dominant species in heavy metals contaminated landfill [
63]. Zolgharnein, et al. [
64] have evaluated the role of
Delfetia and
Methylobacter spp. in the removal of heavy metals such as lead, zinc, copper and cadmium, for the first time. They reported that in the absence of Cu(II) or Cu(I), methanobactin from
Methylosinus trichosporium OB3b will bind Ag(I), Co(II), Au(III), Cd(II), Hg(II), Fe(III), Ni(II), Mn(II), Pb(II), Zn(II), or U(VI) which represents the ability of
Methylosinus to bioremediate heavy metals [
55]. Similarly,
Methylococcus capsulatus has been found to detoxify lead and zinc at low pH [
65] and chromium [
17]. The type II methanotrophic strains identified in this study have been widely utilized for the bioremediation of heavy metals [
16,
66].
These results imply that many methanotrophs, including Methylobacter, Methylococcus, Methylosinus, Methylocystis and Methylmicrobium have an unacknowledged metabolic potential and multifaceted capabilities for bio-remediating heavy metals. Higher type I methanotrophic activities were confirmed in all treatments, which might be due to the higher impact of intrinsic biotic and abiotic characteristics of Xuan Wu Lake.