Growth characteristics of A. ferrooxidans YQ-N3
In this study, sediment was collected from a river polluted by the AMD of an abandoned mine in Shanxi, China. The sediment samples were treated with sterile water and inoculated in 9K-Fe2+ liquid medium for enrichment culture. After cultivation, all three replicate cultures exhibited the same characteristics. The color of the medium gradually changed from light blue-green to turbid, after which it became clear red-brown after continued cultivation. After culturing on 9K-Fe2+ solid medium for 30 days via the dilution coating method, colonies as shown in Fig. 1-A appeared on the plate. The colonies were round, with a prominent center, neat edges, and yellow surrounding. After repeated solid-liquid alternating culture, the bacteria were collected by high-speed freezing and centrifugation, and DNA was extracted for gene sequencing. After gene amplification, there was a single clear target band, which was tentatively named A. ferrooxidans YQ-N3 according to the comparison between the sequencing results and BLAST. The collected bacteria were diluted several times with normal saline and then observed by Gram staining, thus confirming that the bacteria were gram-negative. Optical microscope observations showed that the bacteria were rod-shaped (both as single cells and aggregates) and were able to swim rapidly. SEM analyses indicated that the bacteria had a short rod shape with blunt rounded ends, as shown in Fig. 1-B, with a length of approximately 0.8–1.2 µm and a width of 0.2–0.5 µm.
Upon comparing the 16S rRNA sequences in the GenBank database, the 19 strains that were closest at the species level were selected, and the NJ (Neighbor-Joining) method was used to construct the phylogenetic tree of A. ferrooxidans YQ-N3 using the MEGA 6.0 software, as shown in Fig. 2. Phylogenetic tree analyses demonstrated that A. ferrooxidans YQ-N3 was distinct from Pseudomonas, Oceanicoccus, Luteimonas, and Xanthomonas, but appeared to be related to A. thiooxidans, A. ferridurans, and A. ferrivorans. Particularly, the 16S rDNA sequence of A. ferrooxidans YQ-N3 was more than 99.93% similar to that of A. ferrooxidans ATCC23270, and it was therefore concluded that the isolated strains were A. ferrooxidans.
After screening and obtaining pure A. ferrooxidans YQ-N3 cultures, the optimum growth conditions were explored. The results showed that A. ferrooxidans YQ-N3 could grow at a pH range of 1.6–2.4 and a temperature range of 20–35 ℃, which was consistent with previous studies. The strain achieved optimal growth at an initial pH of 1.8, a 30 ℃ temperature, and an inoculum volume of 10%. Under these conditions, A. ferrooxidans YQ-N3 could reach a density of up to 2.3 × 108 cells/mL after 16 h.
Genome overview of A. ferrooxidans YQ-N3
Upon sequencing the whole genome of A. ferrooxidans YQ-N3, short-length circular consensus sequencing (ccs) reads were filtered out from the raw data as a quality control measure. The filtered ccs reads were denovo assembled, and the assembled genome was corrected for errors. Once the assembly was completed, genome analysis and functional annotation were performed. Whole-genome analysis indicated that the genome size of A. ferrooxidans YQ-N3 was 3,217,720 bp, including one circular chromosome and five circular plasmids (plasmid A, plasmid B, plasmid C, plasmid D, and plasmid E). The sequence length of the circular chromosome was 3,043,496 bp and the GC content was 58.64%. To fully demonstrate the genome features of A. ferrooxidans YQ-N3, the genome circle was drawn using the Circos software (Fig. 3). The coding sequences (CDs) in the genome were predicted using the glimmer, GeneMarkS, and prodigal software. A total of 3,200 predicted CD genes, 6 rRNAs, 46 tRNAs, and 18 sRNAs were predicted. The sequence length of the circular Plasmid A was 79,659 bp and the GC content was 61.64%. The length of Plasmid B was 34,460 bp and the GC content was 60.54%. The sequence length of Plasmid C was 29,178 bp and the GC content was 62.67%. The sequence length of plasmid D was 23,017 bp and the GC content was 60.69%. The sequence length of plasmid E was 7,910 bp and the GC content was 52.40%. Among them, the Plasmid E sequence had not been annotated in the reference database. The complete genome sequences for the main chromosome and plasmids A-E of A. ferrooxidans YQ-N3 were submitted to the GenBank database under accession numbers CP084172.1, CP084173.1, CP084174.1, CP084175.1, CP084176.1, and CP084177.1, respectively.
Genes associated with iron and sulfur metabolism
In-depth characterization of its iron and sulfur metabolism system from a functional gene perspective would not only provide key insights into the physiology of these microorganisms but could establish a theoretical basis for the efficient application of these bacteria in industry. According to the KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation results of A. ferrooxidans YQ-N3, the genes related to iron and sulfur metabolism are summarized in Table 2. Our findings demonstrated that this strain possesses multiple genes related to Fe2+ oxidation metabolism, including cyc2, rus, petA, petB, coxA, coxB, and coxC. A. ferrooxidans can oxidize various reductive inorganic sulfides, including sulfur, sulfides, thiosulfates, tetrathionates, and sulfites. A. ferrooxidans YQ-N3 possesses multiple genes related to sulfur metabolism, including sqr, doxDA, moaD, cysD, hdrA2, hdrB2, hdrC2, and thiS.
Table 2
General features and genomic comparison between A. ferrooxidans YQ-N3 and selected representatives.
Strain
|
Geographic origin
|
BioProject
|
Genome size (Mb)
|
GC%
|
CDS
|
Plasmids
|
Scaffolds
|
A. ferrooxidans YQ-N3
|
Shanxi, China
|
PRJNA543567
|
3.21772
|
58.7316%
|
3132
|
5
|
6
|
A. ferrooxidans YNTRS-40
|
-
|
PRJNA543563
|
3.25704
|
58.4696%
|
3168
|
1
|
2
|
A. ferrooxidans ATCC23270
|
Bituminous coal
mine effluent
|
PRJNA543564
|
2.9824
|
58.8%
|
2927
|
-
|
1
|
A. ferrooxidans BYM
|
Baiyin, China
|
PRJNA543565
|
3.2555
|
58.4696%
|
3134
|
1
|
2
|
A. ferrooxidans NFP31
|
Volcanic ash deposits on Miyake-jima, Japan
|
PRJNA543566
|
3.24985
|
58.4724%
|
3193
|
1
|
2
|
A. ferrooxidans ATCC53993
|
-
|
PRJNA543568
|
2.88504
|
58.9%
|
2811
|
-
|
1
|
A. ferrooxidans CCM4253
|
Mine waters, Czech Republic
|
PRJNA5435690
|
3.19656
|
58.6%
|
3073
|
-
|
15
|
A. ferrooxidans DSM16786
|
Wudalianchi Heilongjiang, China
|
PRJNA543570
|
3.67579
|
58.4%
|
3636
|
-
|
49
|
A. ferrooxidans YQH-1
|
Wudalianchi volcano water, China
|
PRJNA543571
|
3.11122
|
58.60%
|
3012
|
-
|
96
|
A. ferrooxidans Hel18
|
Flue dust
|
PRJNA543572
|
3.10916
|
58.6%
|
3065
|
-
|
123
|
A. ferrooxidans PQ506
|
Santiago, Chile
|
PRJNA543573
|
3.37146
|
58.3%
|
3342
|
-
|
277
|
A. ferrooxidans PQ505
|
Santiago, Chile
|
PRJNA543574
|
3.51657
|
58.4%
|
3534
|
-
|
305
|
A. ferrooxidans CF3
|
Santiago, Chile
|
PRJNA543575
|
3.01139
|
58.7%
|
3057
|
-
|
310
|
A. ferrooxidans F221
|
-
|
PRJNA543576
|
3.00995
|
58.7%
|
3006
|
-
|
360
|
A. ferrooxidans BY0502
|
-
|
PRJNA543577
|
2.97667
|
56.8%
|
3026
|
-
|
295
|
A. ferrooxidans BY-3
|
Gansu, China
|
PRJNA543578
|
3.83234
|
57.8%
|
3777
|
-
|
194
|
A. ferrooxidans RVS1
|
Andacollo gold mining area, Argentina
|
PRJNA543579
|
2.82631
|
58.8%
|
2705
|
-
|
49
|
A. ferrooxidans DLC-5
|
-
|
PRJNA543580
|
4.18422
|
57.6%
|
-
|
-
|
2090
|
A. ferrooxidans COP1
|
-
|
PRJNA543581
|
3.008
|
58.8%
|
3855
|
-
|
1561
|
A. ferrooxidans S10
|
Santiago, Chile
|
PRJNA543582
|
2.953
|
58.8%
|
3998
|
-
|
1827
|
Note: “-” indicates no data. |
Table 3
List of selected genes identified in strain YQ-N3 strain via KEGG annotation, including genes for iron and sulfur metabolism.
Gene name
|
Gene length
|
KEGG gene ID
|
KEGG orthology description
|
fdxA
|
327bp
|
afr: AFE_0014
|
ferredoxin
|
feoA
|
288bp
|
afr:AFE_2523
|
ferrous iron transport protein A
|
feoB
|
2349bp
|
afr:AFE_2524
|
ferrous iron transport protein B
|
hemH
|
1014bp
|
afr:AFE_0179
|
protoporphyrin/coproporphyrin ferrochelatase
|
cyc2
|
1458bp
|
afr:AFE_3153
|
iron: rusticyanin reductase [EC:1.16.9.1]
|
resB
|
1845bp
|
afr:AFE_3112
|
cytochrome c biogenesis protein
|
hyaC
|
756bp
|
afr:AFE_2429
|
Ni/Fe-hydrogenase 1 B-type cytochrome subunit
|
coxA
|
1884bp
|
afe:Lferr_2747
|
cytochrome c oxidase subunit I
|
coxB
|
765bp
|
afe:Lferr_2748
|
cytochrome c oxidase subunit II
|
coxC
|
555bp
|
afr:AFE_3148
|
cytochrome c oxidase subunit III
|
porA
|
984bp
|
tti:THITH_06955
|
pyruvate ferredoxin oxidoreductase alpha subunit
|
porB
|
1164bp
|
tig:THII_3692
|
pyruvate ferredoxin oxidoreductase beta subunit
|
porC
|
603bp
|
tti:THITH_06960
|
pyruvate ferredoxin oxidoreductase gamma subunit
|
-
|
2340bp
|
afe:Lferr_1935
|
iron complex outermembrane receptor protein
|
fdxA
|
621bp
|
afr:AFE_1844
|
ferredoxin
|
fdx
|
306bp
|
afr:AFE_1541
|
ferredoxin, 2Fe-2S
|
rus
|
564bp
|
afr:AFE_3146
|
rustycanin
|
-
|
1035bp
|
afe:Lferr_1212
|
iron complex transport system substrate-binding protein
|
petA
|
621bp
|
afe:Lferr_2707
|
ubiquinol-cytochrome c reductase iron-sulfur subunit
|
petB
|
1209bp
|
afe:Lferr_2708
|
ubiquinol-cytochrome c reductase cytochrome b subunit
|
petC
|
729bp
|
afr:AFE_3111
|
ubiquinol-cytochrome c reductase cytochrome c1 subunit
|
erpA
|
372bp
|
afj:AFERRID_10140
|
iron-sulfur cluster insertion protein
|
moaD
|
243bp
|
afr:AFE_0975
|
sulfur-carrier protein
|
iscA
|
324bp
|
afr:AFE_0675
|
iron-sulfur cluster assembly protein
|
doxDA
|
1083bp
|
afe: Lferr_0045
|
thiosulfate dehydrogenase (quinone)
|
doxA
|
1083bp
|
afj:AFERRID_13680
|
thiosulfate dehydrogenase (quinone) small subunit
|
cysN
|
1353bp
|
afr:AFE_3125
|
sulfate adenylyltransferase subunit 1
|
cysD
|
939bp
|
afe:Lferr_2723
|
sulfate adenylyltransferase subunit 2
|
cysJ
|
1767bp
|
afr:AFE_3121
|
sulfite reductase (NADPH) flavoprotein alpha-component
|
cysI
|
1692bp
|
afr:AFE_3122
|
sulfite reductase (NADPH) hemoprotein beta-component
|
cysH
|
738bp
|
afr:AFE_3123
|
phosphoadenosine phosphosulfate reductase
|
sqr
|
1140bp
|
afr:AFE_2601
|
sulfide:quinone oxidoreductase
|
hdrA2
|
1056bp
|
afr:AFE_2553
|
heterodisulfide reductase subunit A2
|
hdrB2
|
912bp
|
afr:AFE_2550
|
heterodisulfide reductase subunit B2
|
hdrC2
|
720bp
|
afr:AFE_2551
|
heterodisulfide reductase subunit C2
|
thiS
|
201bp
|
afr:AFE_0642
|
sulfur carrier protein
|
Gene traits related to environmental adaptability
Based on protein sequence alignments, the coding genes predicted by A. ferrooxidans YQ-N3 were compared with the COG (Clusters of Orthologous Groups of proteins) and KEGG databases for functional annotation, and the corresponding functional annotations were obtained Fig. 4 shows the COG annotation classification statistics of A. ferrooxidans YQ-N3. Our findings indicated that A. ferrooxidans YQ-N3 has a total of 2,571 genes annotated in COG, accounting for approximately 80.34% of the total number of genes, and these genes are classified into 22 COG types. Among these genes, a total of 728 with unknown function were identified. The proportion of genes related to replication, recombination and repair (L), cell wall/membrane/envelope biogenesis (M), energy production and conversion (C), and inorganic ion transport and metabolism (P) was slightly higher than that of genes with other functions, accounting for 9.45%, 6.76%, 7.0%, and 6.2% of all annotated genes, respectively. Additionally, defense mechanisms accounted for 5.98% of all annotated genes, suggesting that this strain has a strong ability to self-repair and resist harsh environments (Zhang et al. 2019).
The results of A. ferrooxidans YQ-N3 genome annotation using the KEGG pathway database indicated that the gene functions of this bacterium could be mainly divided into six categories, including Metabolism, Cellular Processes, Human Diseases, Genetic Information Processing, Organismal Systems, and Environmental Information Processing, as shown in Fig. 5. Among these classifications, 1,353 functional genes related to Metabolism could be divided into 12 types, 81 functional genes related to Cellular Processes could be divided into 3 types, 124 genes related to Human Diseases could be divided into 10 types, 183 functional genes related to Genetic Information Processing could be divided into 6 categories, 56 genes related to Organic Systems could be divided into 8 categories, and 177 genes related to Environmental Information Processing could be divided into two categories.
Gene island prediction of the A. ferrooxidans YQ-N3 genome was performed using IslandViewer (Bertelli et al. 2017). The prophage prediction was performed using Phage_Finder (Fouts 2006). CRISPR-Cas (Clustered Regularly Interspersed Short Palindromic Repeats) prediction was performed using Minced (Bland et al. 2007). The predicted results indicated that the chromosome of A. ferrooxidans YQ-N3 contained 10 gene islands, whereas Plasmid A contained one gene island, and these gene islands contained a total of 181 CDs. Interestingly, our analyses indicated that a prophage genome had integrated into the chromosome of A. ferrooxidans YQ-N3. The sequence length of the prophage was 17,089 bp, and had 19 CDs, and a 58.64% GC content. The chromosome of A. ferrooxidans YQ-N3 was predicted to contain four CRISPRs with average repeat lengths of 27, 28, 23, and 27 bp, respectively. The average lengths of the spacer sequences were 27, 32, 43, and 45 bp, respectively.
Oxidation of Fe 2+ , S0, and pyrite by A. ferrooxidans YQ-N3
During Fe2+ oxidation by A. ferrooxidans YQ-N3, the color of the medium gradually shifted from clear light green to yellow-green at first, and after 16 h, the medium completely changed to reddish-brown. Then, a yellow precipitate gradually appeared in the reaction system. The blank control group did not exhibit any obvious color change. Figure 6 shows the trends of pH, ORP, Fe2+ content, and total iron content of the medium during Fe2+ oxidation. According to Fig. 1a, the pH in the experiment involving A. ferrooxidans YQ-N3 increased from 1.8 to 2.16 after 48 h, and then rapidly decreased to 1.7. ORP also increased throughout the experiment but this trend decelerated slightly after 48 h. In the blank control group, pH and ORP showed a slight upward trend throughout the experiment. This can be attributed to H+ consumption by oxidation of Fe2+ into Fe3+ (Offeddu et al. 2015). Fe3+, SO42−, and other cations (e.g, K+, Na+, NH4+) begin to react as the Fe3+ content increases, which results in the production of H+ and the minerals jarosite (KFe3(SO4)2(OH)6) and ammoniojarosite ((NH4)2Fe6(SO4)4(OH)12) (see Fig. 7), which is consistent with previous studies (Nazari et al. 2014). Therefore, the pH rapidly decreased and the ORP increase rate slowed down in the later stage of the experimental group (Li et al. 2021). Figure 1b shows that the Fe2+ content in the experiment involving A. ferrooxidans YQ-N3 continued to decrease rapidly, and the Fe2+ content reached 0 mg/L when the reaction reached 60 h. Total iron concentration decreased by approximately 2700 mg/L during the whole reaction period, whereas only 14.5% Fe2+ was oxidized and the total iron concentration in the blank control group remained unchanged for 96 h. This can be attributed to the fact that A. ferrooxidans YQ-N3 accelerates the oxidation of Fe2+, and the Fe3+ generated by the reaction leads to a decrease in the total iron content in the system. Fe2+ is naturally oxidized and generates free Fe3+ in the blank control group during the reaction, and therefore the total iron content remained unchanged in the system.
During oxidation of S0 by A. ferrooxidans YQ-N3, the medium gradually changed from clear and colorless to light yellow, the particle size of sulfur powder gradually became smaller, and its hydrophobicity gradually weakened. However, no obvious change was observed during the reaction of the blank control group. This phenomenon was mainly because A. ferrooxidans adsorbed on the surface of S when oxidizing it and secreted hydrophilic organic substances that covered the surface of S to enhance its hydrophilicity (Konishi et al. 1995; Knickerbocker 2000). During the oxidation of S0 by A. ferrooxidans YQ-N3, the pH value of the medium decreased from 2.2 to 1.74, whereas the SO42− concentration of the medium increased to 3896.66 mg/L after 45 d (Fig. 8). This can be attributed to the release of H+ through the oxidation of S0 to SO42−.
Figure 9 illustrates the trends of pH, ORP, SO42−, Fe2+, and Fet during pyrite oxidation by A. ferrooxidans YQ-N3. In the experimental group, the pH value decreased from 1.9 to 1.2 after 45 d, the ORP increased from 353 to 632 mV after 45 d, the SO42− concentration increased to 5433.8 mg/L after 45 d, the Fe2+ content decreased drastically from 43.6 mg/L to a negligible level after 23 d, and the Fet content increased to 3,853 mg/L after 45 d. Characterization by XRD revealed that pyrite oxidized by A. ferrooxidans YQ-N3 contains jarosite and FeOOH, in addition to quartz impurities (Fig. 10). In the blank control group, the pH value decreased slightly, the ORP decreased and then slightly increased, the SO42− concentration increased to 563 mg/L after 45 d, the Fe2+ content increased to 178.2 mg/L after 45 d, and the Fet concentration increased to 197.7 mg/L after 45 d. This could be because the H+ released by the oxidation of S0 into SO42− exceeded the levels consumed by Fe2+ oxidation. A. ferrooxidans YQ-N3 can accelerate the oxidation of Fe2+ to Fe3+ and pyrite oxidation can be carried out continuously and rapidly. The dissolved Fe mainly exists in the form of Fe3+ in this system, and therefore ORP continues to increase.