# Genomic Analysis Unravels Reduced Inorganic Sulfur Compound Oxidation of Heterotrophic AcidophilicAcidicaldus sp. Strain DX-1

Authors: Yuanyuan Liu; Hongying Yang; Xian Zhang; Yunhua Xiao; Xue Guo; Xueduan Liu
Journal: BioMed Research International
 (2016)

---

## Abstract

Although reduced inorganic sulfur compound (RISC) oxidation in many chemolithoautotrophic sulfur oxidizers has been investigated in recent years, there is little information about RISC oxidation in heterotrophic acidophiles. In this study,Acidicaldus sp. strain DX-1, a heterotrophic sulfur-oxidizing acidophile, was isolated. Its genome was sequenced and then used for comparative genomics. Furthermore, real-time quantitative PCR was performed to identify the expression of genes involved in the RISC oxidation. Gene encoding thiosulfate: quinone oxidoreductase was present inAcidicaldus sp. strain DX-1, while no candidate genes with significant similarity to tetrathionate hydrolase were found. Additionally, there were genes encoding heterodisulfide reductase complex, which was proposed to play a crucial role in oxidizing cytoplasmic sulfur. Like many heterotrophic sulfur oxidizers,Acidicaldus sp. strain DX-1 had no genes encoding enzymes essential for the direct oxidation of sulfite. An indirect oxidation of sulfite via adenosine-5′-phosphosulfate was proposed inAcidicaldusstrain DX-1. However, compared to other closely related bacteriaAcidiphilium cryptum andAcidiphilium multivorum, which harbored the genes encoding Sox system, almost all of these genes were not detected inAcidicaldus sp. strain DX-1. This study might provide some references for the future study of RISC oxidation in heterotrophic sulfur-oxidizing acidophiles.


---

## Body

## 1. Introduction

Due to the evolutionary importance and applied perspectives of extreme acidophiles (optimal growth at or below pH 3.0), researchers have paid much attention to these microorganisms [1–3]. Extreme acidophiles are widely used in the bioprocessing of minerals and bioremediation of acidic and metal-enriched water [4–6]. Many chemolithoautotrophic acidophiles, such asLeptospirillum spp. andSulfolobus metallicus, can accelerate the dissolution of sulfide minerals by oxidizing ferrous iron or reduced inorganic sulfur compound (RISC) [7]. In addition to obligate autotrophs, many heterotrophic and mixotrophic acidophiles also play key roles in iron and sulfur cycling, such asFerrimicrobium acidiphilum, Sulfobacillusspp.,Thermoplasma sp., andAcidiphiliumspp. [7–12].Species from theAcidicaldus genus are moderately thermophilic and obligate heterotrophic acidophiles within the Alphaproteobacteria class. The type speciesAcidicaldus organivorans Y008, which was isolated from a geothermal site in Yellowstone National Park, is a member of facultatively anaerobes and can grow by ferric iron respiration in the anaerobic environment [13]. The closest phylogenetic relatives ofAc. organivoransY008 are acidophilic heterotrophs, includingAcidiphilium[14],Acidisphaera[15], andAcetobacter[16]. Previous studies showed thatAcidicaldus organivorus (the synonym ofAc.organivorans) harbored some unique physiological characteristics, which differed from other acidophilic heterotrophs. Compared toAcidisphaera rubrifaciens, which grows at pH 4.5–5.0 and optimal temperature 30–35°C,Acidicaldus organivorus requires a lower pH (2.5–3.0) and higher temperature (optima 50–55°C).Ac. organivorus can oxidize elemental sulfur to sulfate, though it cannot perform autotrophic growth on sulfur in organic-free media. Additionally,Ac. organivorus can catalyze the dissimilatory reduction of ferric iron, whereas the type strain ofAs. rubrifaciens does not grow by ferric iron respiration in the absence of oxygen [13].Although knowledge in the last decade has greatly advanced our understanding on the acidophilic bacteria, it is scarce to attempt to study the physiology and genetics of the heterotrophic sulfur-oxidizing acidophiles. Compared to autotrophic sulfur-oxidizing bacteria, metabolisms of these heterotrophs are more versatile [17]. They can not only use organic matters as energy sources but also obtain energy derived from the oxidation of elemental sulfur or RISC. These traits can confer the heterotrophic sulfur-oxidizing bacteria dual contributions to microbial communities in bioleaching systems: (i) while organic matters inhibit the growth of some facultative autotrophic bacteria [18, 19], heterotrophs which can consume the organic matters in the microbial communities might relieve the inhibiting effect of organic matters on the autotrophs; (ii) RISC oxidation would accelerate dissolution of acid soluble metal sulfides. In this study, we chose a strain, isolated from Dexing Copper Mine (Jiangxi, China), to explore its physiologic and genetic properties. The genome sequence of this strain was acquired using the whole genome sequencing and comparative genomics was then performed, aiming to propose its putative RISC oxidation model.


## 2. Materials and Methods

### 2.1. Isolation and Cultivation of Bacteria

Samples from the surface of effusion pool were inoculated in 9K liquid media [20] supplemented with 5 g/L elemental sulfur, 4.5% (w/v) ferrous sulfate, and 0.02% yeast extract (Oxoid, UK). Strains were purified by repeated single colony (in triplicate) isolation on 9K plates solidified with 0.02% yeast extract and grown at 45°C in 9K liquid media, pH 3.0, with 0.02% (w/v) yeast extract and 10 mM glucose. All procedures were operated aseptically.


### 2.2. DNA Extraction, 16S rRNA Gene Amplification, and Phylogenetic Analysis

At later exponential phase (generally 4th day), strain DX-1 was harvested by centrifugation (12,000 g for 10 min at 4°C). Genomic DNA was extracted from the pelleted cells using TIANamp Bacteria DNA Kit (Tiangen Biotech, China) according to the manufacturer’s instruction and finally was resuspended in TE buffer. 16S rRNA gene sequence was PCR amplified from the genomic DNA extracts of strain DX-1. Amplification was performed in 50μL reaction mixtures containing 1 μL of DNA extracts, 1 μL each of 10 mM forward and reverse primers, 25 μL of universal Taq PCR Master Mix (Tiangen Biotech, China), and 22 μL of deionized water. The PCR conditions for amplification were as follows: 94°C for 4 min, then 32 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 45 s, followed by a final extension at 72°C for 10 min. PCR product of 16S rRNA gene was purified directly with the QIAquick PCR Purification Kit (Qiagen, Germany) and sequenced.16S rRNA gene sequence of strain DX-1 was aligned with other acidophilic Alphaproteobacteria sequences from the GenBank database using the CLUSTAL program [21]. This alignment was used to make a distance matrix. Phylogenetic trees were constructed using the neighbor joining method by Molecular Evolutionary Genetics Analysis 5.2 (MEGA, version 5.2) [22]. Bootstrap analysis was carried out on 1000 replicate input data sets.


### 2.3. Genome Sequencing, Assembly, and Annotation

Genomic library construction, sequencing, and assembly were performed according to previous methods [23]. Given the high GC content of genome ofAcidicaldusspp., more than 600 Mb pair-ends reads with a depth of over 200-fold coverage were obtained. After genome assembly, 375 scaffolds were acquired using SOAPdenovo package [24]. The completeness (95.34%) of strain DX-1 genome was estimated using the CheckM [25]. Coding sequences (CDSs) were then predicted with the ORF finders Glimmer [26]. All CDSs were annotated by comparison with the public available databases nonredundant NCBI [27] and KEGG [28] using the annotation software BLAST [29]. The unassigned CDSs were further annotated using the HMMPfam program [30]. And the hidden Markov models for the protein domains were obtained from the Pfam database [31]. The software programs tRNAscan-SE and RNAmmer were used for the identifications of tRNA and rRNA, respectively [32, 33].


### 2.4. Comparative Genomics

The genomes ofAp. cryptum JF-5 andAp. multivorum AIU301 were retrieved from the NCBI database. Orthologs betweenAcidicaldus strain DX-1,Ap.cryptum JF-5, andAp.multivorum AIU301 were detected via an all-versus-all reciprocal BLASTP search against their own protein sets, respectively [29]. The best sequence similarities were obtained using two cut-off values: E-value ≤ 1
e
-
05 and minimal coverage by local alignment = 70%.


### 2.5. RNA Extraction and Real-Time Quantitative PCR Analysis

To investigate the transcript profiles of genes associated with RISC oxidation,Acidicaldus sp. strain DX-1 was cultured at 45°C in 9K liquid media, pH 3.0 with 0.02% (w/v) yeast extract and 10 mM glucose. And the additional 1% (w/v) elemental sulfur (S0) was added in another group. Microbial cells were harvested at mid-exponential growth phase according to bacterial growth curve. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, USA), treated with RNase-free DNase I (Qiagen, Valencia, USA) and purified with RNeasy Kit (Qiagen, Valencia, USA). Subsequently, single-stranded cDNA was synthesized with ReverTra Ace qPCR RT Kit (Toyobo, Japan), according to the manufacturer’s protocol. The cDNA was stored at −80°C, which was used for further RT-qPCR analysis. All experiments under the same conditions were performed in triplicate. The gene expression level with fold change ≥ 1 were upregulated and otherwise downregulated.Specific primers targeting selected genes putatively involved in RISC oxidation were designed for real-time PCR analysis (Table1). The real-time PCR amplification was performed with iCycler iQ Real-Time PCR detection system (Bio-Rad Laboratories Inc., Hercules, USA). 25 μL reaction mixtures contained 12.5 μL SYBR Green Real-Time PCR Master Mix (Toyobo Co., Ltd., Osaka, Japan), 0.5 μL single-stranded cDNA, 1 μL each of 10 μM forward and reverse primers, and 10 μL deionized water. The specific amplification protocol was as follows: 95°C for 5 min, 40 cycles of 95°C for 20 s, 55°C for 15 s, 72°C for 15 s, and a final incubation of 72°C for 10 min. The expression of each gene was determined from triplicate reactions in a single real-time PCR amplification. To standardize the quantification of the selected target genes, 16S rRNA and glyceraldehyde 3-phosphate dehydrogenase (gapdh) were used as transcription controls to regulate the random and systematic errors. The expressions of selected genes in the medium with S0 and without S0 were calculated for further analysis.Table 1
Primers used for real-time quantitative PCR in this study.
Targeted gene
Primers
Product size [bp]
Melting temp. [°C]
Primer sequence
hdrA
FwdRev
262
55
GCACCGGCTTCACCCATTTCAGCTGCTCGCGGATCTCCAT
hdrB
FwdRev
184
55
ACTCGACCTGGTATGATTGCTGCTTGCCGATCCACTGGCTCTT
hdrC
FwdRev
161
55
CACGGCGTGGGAGGTCAATCCGGCGAGCCAGTCTTCGTAAT
hdrB2
FwdRev
201
55
ACAACGAACACGCGAAAGAGGGACGGACGTACATGCAGCCATA
rhd
FwdRev
238
55
CCCAACGAGGGCAAGGACGGAGAGGCGGAGCAGCCACCAG
doxD
FwdRev
103
55
CAACTGGATGGCGCACAAACGCATAGAGCAGCGGGAAA
sat
FwdRev
119
55
GACCTTCTCGAAAGCCAGAGCTTGCGGGCAAGCCAAACC
sqr
FwdRev
182
56
CGATTTCGGCGACTCGGGTGTGCTGGCGGAGTAGGAATTTCTCAT
sqr2
FwdRev
258
55
AAACCGGGTGCGCTTCGTTGCGGCTCCTTCGGATTGC


### 2.6. Nucleotide Sequence Accession Number

The draft genome sequence ofAcidicaldus strain DX-1 has been deposited at DDBJ/EMBL/GenBank under accession number JPYW00000000.


## 2.1. Isolation and Cultivation of Bacteria

Samples from the surface of effusion pool were inoculated in 9K liquid media [20] supplemented with 5 g/L elemental sulfur, 4.5% (w/v) ferrous sulfate, and 0.02% yeast extract (Oxoid, UK). Strains were purified by repeated single colony (in triplicate) isolation on 9K plates solidified with 0.02% yeast extract and grown at 45°C in 9K liquid media, pH 3.0, with 0.02% (w/v) yeast extract and 10 mM glucose. All procedures were operated aseptically.


## 2.2. DNA Extraction, 16S rRNA Gene Amplification, and Phylogenetic Analysis

At later exponential phase (generally 4th day), strain DX-1 was harvested by centrifugation (12,000 g for 10 min at 4°C). Genomic DNA was extracted from the pelleted cells using TIANamp Bacteria DNA Kit (Tiangen Biotech, China) according to the manufacturer’s instruction and finally was resuspended in TE buffer. 16S rRNA gene sequence was PCR amplified from the genomic DNA extracts of strain DX-1. Amplification was performed in 50μL reaction mixtures containing 1 μL of DNA extracts, 1 μL each of 10 mM forward and reverse primers, 25 μL of universal Taq PCR Master Mix (Tiangen Biotech, China), and 22 μL of deionized water. The PCR conditions for amplification were as follows: 94°C for 4 min, then 32 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 45 s, followed by a final extension at 72°C for 10 min. PCR product of 16S rRNA gene was purified directly with the QIAquick PCR Purification Kit (Qiagen, Germany) and sequenced.16S rRNA gene sequence of strain DX-1 was aligned with other acidophilic Alphaproteobacteria sequences from the GenBank database using the CLUSTAL program [21]. This alignment was used to make a distance matrix. Phylogenetic trees were constructed using the neighbor joining method by Molecular Evolutionary Genetics Analysis 5.2 (MEGA, version 5.2) [22]. Bootstrap analysis was carried out on 1000 replicate input data sets.


## 2.3. Genome Sequencing, Assembly, and Annotation

Genomic library construction, sequencing, and assembly were performed according to previous methods [23]. Given the high GC content of genome ofAcidicaldusspp., more than 600 Mb pair-ends reads with a depth of over 200-fold coverage were obtained. After genome assembly, 375 scaffolds were acquired using SOAPdenovo package [24]. The completeness (95.34%) of strain DX-1 genome was estimated using the CheckM [25]. Coding sequences (CDSs) were then predicted with the ORF finders Glimmer [26]. All CDSs were annotated by comparison with the public available databases nonredundant NCBI [27] and KEGG [28] using the annotation software BLAST [29]. The unassigned CDSs were further annotated using the HMMPfam program [30]. And the hidden Markov models for the protein domains were obtained from the Pfam database [31]. The software programs tRNAscan-SE and RNAmmer were used for the identifications of tRNA and rRNA, respectively [32, 33].


## 2.4. Comparative Genomics

The genomes ofAp. cryptum JF-5 andAp. multivorum AIU301 were retrieved from the NCBI database. Orthologs betweenAcidicaldus strain DX-1,Ap.cryptum JF-5, andAp.multivorum AIU301 were detected via an all-versus-all reciprocal BLASTP search against their own protein sets, respectively [29]. The best sequence similarities were obtained using two cut-off values: E-value ≤ 1
e
-
05 and minimal coverage by local alignment = 70%.


## 2.5. RNA Extraction and Real-Time Quantitative PCR Analysis

To investigate the transcript profiles of genes associated with RISC oxidation,Acidicaldus sp. strain DX-1 was cultured at 45°C in 9K liquid media, pH 3.0 with 0.02% (w/v) yeast extract and 10 mM glucose. And the additional 1% (w/v) elemental sulfur (S0) was added in another group. Microbial cells were harvested at mid-exponential growth phase according to bacterial growth curve. Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, USA), treated with RNase-free DNase I (Qiagen, Valencia, USA) and purified with RNeasy Kit (Qiagen, Valencia, USA). Subsequently, single-stranded cDNA was synthesized with ReverTra Ace qPCR RT Kit (Toyobo, Japan), according to the manufacturer’s protocol. The cDNA was stored at −80°C, which was used for further RT-qPCR analysis. All experiments under the same conditions were performed in triplicate. The gene expression level with fold change ≥ 1 were upregulated and otherwise downregulated.Specific primers targeting selected genes putatively involved in RISC oxidation were designed for real-time PCR analysis (Table1). The real-time PCR amplification was performed with iCycler iQ Real-Time PCR detection system (Bio-Rad Laboratories Inc., Hercules, USA). 25 μL reaction mixtures contained 12.5 μL SYBR Green Real-Time PCR Master Mix (Toyobo Co., Ltd., Osaka, Japan), 0.5 μL single-stranded cDNA, 1 μL each of 10 μM forward and reverse primers, and 10 μL deionized water. The specific amplification protocol was as follows: 95°C for 5 min, 40 cycles of 95°C for 20 s, 55°C for 15 s, 72°C for 15 s, and a final incubation of 72°C for 10 min. The expression of each gene was determined from triplicate reactions in a single real-time PCR amplification. To standardize the quantification of the selected target genes, 16S rRNA and glyceraldehyde 3-phosphate dehydrogenase (gapdh) were used as transcription controls to regulate the random and systematic errors. The expressions of selected genes in the medium with S0 and without S0 were calculated for further analysis.Table 1
Primers used for real-time quantitative PCR in this study.
Targeted gene
Primers
Product size [bp]
Melting temp. [°C]
Primer sequence
hdrA
FwdRev
262
55
GCACCGGCTTCACCCATTTCAGCTGCTCGCGGATCTCCAT
hdrB
FwdRev
184
55
ACTCGACCTGGTATGATTGCTGCTTGCCGATCCACTGGCTCTT
hdrC
FwdRev
161
55
CACGGCGTGGGAGGTCAATCCGGCGAGCCAGTCTTCGTAAT
hdrB2
FwdRev
201
55
ACAACGAACACGCGAAAGAGGGACGGACGTACATGCAGCCATA
rhd
FwdRev
238
55
CCCAACGAGGGCAAGGACGGAGAGGCGGAGCAGCCACCAG
doxD
FwdRev
103
55
CAACTGGATGGCGCACAAACGCATAGAGCAGCGGGAAA
sat
FwdRev
119
55
GACCTTCTCGAAAGCCAGAGCTTGCGGGCAAGCCAAACC
sqr
FwdRev
182
56
CGATTTCGGCGACTCGGGTGTGCTGGCGGAGTAGGAATTTCTCAT
sqr2
FwdRev
258
55
AAACCGGGTGCGCTTCGTTGCGGCTCCTTCGGATTGC


## 2.6. Nucleotide Sequence Accession Number

The draft genome sequence ofAcidicaldus strain DX-1 has been deposited at DDBJ/EMBL/GenBank under accession number JPYW00000000.


## 3. Results and Discussion

### 3.1. Isolation and Phylogenetic Analysis

Several strains (named DX-1, DX-5, DX-20, and DX-22) from Dexing Copper Mine sampling sites, which showed similar colony and cellular morphologies, were isolated on solid medium with 0.02% yeast extract. All of these strains grew as small, white, round to convex-shaped colonies. Considering that 16S rRNA sequence analysis showed 100% sequence identities, only one (strain DX-1) was selected for further analysis. The rooted phylogenetic tree indicated that 16S rRNA gene of strain DX-1 was assigned into a phyletic cluster withAcidicaldusspp. (Figure 1). Furthermore, 16S rRNA gene of strain DX-1 shared 99% sequence identity with that of type strainAcidicaldus organivoransY008 [13], indicating that strain DX-1 belongs toAcidicaldus species.Figure 1
Phylogenetic tree of 16S rRNA genes showing the relationship of strain DX-1 (inbold) to other acidophilic α-Proteobacteria. Phylogenetic tree was constructed based on an alignment of 1283 bp using the neighbor joining method. The scale bar represents 0.01 nucleotide substitutions per 100. And the database accession numbers of the gene sequences used are given in parentheses.Acidisphaera rubrifaciens HS-AP3 was used to root the tree.


### 3.2. Genome Assembly and Annotation

The draft genome sequence ofAcidicaldus sp. strain DX-1 contained 2,990,377 bp with GC content of 68.48% distributed in 376 scaffolds, which ranged from 1,000 bp to 67,370 bp. Given the high sequence coverage (200-fold), it was likely to identify most genes in the draft genome ofAcidicaldus sp. strain DX-1. Results showed that genome ofAcidicaldus sp. strain DX-1 harbored 3,259 predicted coding sequences (CDSs), which represented 89.83% of the genome, 1 rRNA operon, and 43 tRNA genes (Table 2). The number of CDSs in this strain was less than that inAcidiphilium multivorum AIU301 but more than that inAcidiphilium cryptum JF-5 andAcetobacter pasteurianus.Table 2
General features ofAcidicaldus sp. strain DX-1 in comparison with other heterotrophic acidophiles.
Organism
Acidicaldus sp. strain DX-1
Acidiphilium cryptum JF-5
Acidiphilium multivorumAIU301
Acetobacter pasteurianus
Isolated from
Acid effusion of mill tailing
Acidic Sediments
Pyritic acid mine drainage
Spontaneous cocoa bean heap fermentation
Temperature range
Moderate thermophilic (50~55°C)
Mesophilic (30~35°C)
Mesophilic (30–35°C)
Mesophilic (25–30°C)
pH range
2.5–3.0
3.0–3.5
3.0–3.5
5.4–6.3
Nutrition type
Heterotrophic
Heterotrophic
Heterotrophic
Heterotrophic
Oxygen requirement
Facultatively anaerobic
Facultatively anaerobic
Facultatively anaerobic
Obligately aerobic
Fe(III) reduction
Yes
Yes
Yes
No
Sulfur oxidation
Yes
No
No
No
Genome size (Mb)
2.99
3.39
3.75
2.91
GC content
68.48%
67.99%
67.58%
53.04%
Protein-coding genes
3259
3063
3448
2875
rRNA operon
1
2
2
5
Number of tRNAs
43
47
48
57


### 3.3. Hydrogen Sulfide Oxidation

Phototrophic bacteria such asAllochromatium vinosum and someChlorobium spp. can anaerobically oxidize hydrogen sulfide to sulfur via a flavocytochrome c-sulfide dehydrogenase [34]. In this process, cytochrome c
550 is used as electron acceptor for energy production. However, flavocytochrome c inA. vinosum is not required for phototrophic growth with hydrogen sulfide. Therefore, sulfide: quinone reductase (SQR) is considered to be essential for growth ofA. vinosum with hydrogen sulfide. In chemotrophic bacteria, such as someAcidithiobacillus spp., SQR is also responsible for the oxidation of hydrogen sulfide [35]. InAcidicaldus sp. strain DX-1, no candidate genes encoding flavocytochrome c-sulfide dehydrogenase were identified, while two genes encoding SQR were found (Table S1, see Supplementary Material available online at http://dx.doi.org/10.1155/2016/8137012). One putative SQR shared 64% amino acid similarity with that ofAt. ferrooxidans ATCC 53993, and the other shared 72% amino acid similarity with that ofBradyrhizobium oligotrophicum S58. Both of them had a small NADH binding domain within a larger FAD binding domain. These results suggested that SQR might be responsible for the oxidation of hydrogen sulfide inAcidicaldus sp. strain DX-1.


### 3.4. S4I Pathway

In some strains of the genusAcidithiobacillus, the membrane-bound thiosulfate: quinone oxidoreductase (TQO) and tetrathionate hydrolase (TetH) are responsible for the S4I pathway [36, 37]. InAt. ferrooxidans, TQO which is constituted of subunits DoxA and DoxD catalyzes the conversion of thiosulfate to tetrathionate, and TetH catalyzes the hydrolysis of tetrathionate to generate sulfur, sulfate and thiosulfate. The consecutive reactions catalyzed by TetH and TQO promote the RISC oxidation inAt. ferrooxidans[38]. One ortholog ofdoxDA encoding TQO was identified inAcidicaldus sp. strain DX-1, which shared 61% similarity with that inAt. ferrooxidans (Table S1). Further analysis revealed that the putative TQO had conserved DoxD domain (Pfam: PF04173) and DoxA domain (Pfam: PF07680). Although previous experiments indicated thatAc. organivoransY008 can grow on tetrathionate overlay medium [13], no candidate genes with significant similarity to TetH were found in the draft genome ofAcidicaldus strain DX-1. Future analyses will be necessary to determine whether TQO indeed catalyzes thiosulfate to tetrathionate and TetH gene is on the portion of the missing draft genome ofAcidicaldus sp. strain DX-1.


### 3.5. Sox Pathway

In the alphaproteobacteriumParacoccus pantotrophus, the periplasmic Sox system encoded by 15 genes [39] are comprised of four subunits, including SoxXA, SoxYZ, SoxB, and Sox(CD)2. Each Sox subunit directing its own function was elaborated in previous studies [40]. BothAcidiphilium cryptum andAcidiphilium multivorum have one gene cluster ofsoxXYZABCD to code for four periplasmic proteins. However, no candidate genes encoding Sox system were identified in the genome ofAcidicaldus sp. strain DX-1. Considering the Sox system was absent, S4I pathway was probably the sole way of further oxidation of thiosulfate in periplasm.


### 3.6. Sulfur Oxidation

The cytoplasmic sulfur oxygenase reductase (SOR), which mediates biological oxidation of inorganic sulfur, can catalyze the disproportionation reaction of sulfur to generate sulfite, thiosulfate, and hydrogen sulfide [23, 38, 40–42]. Sulfur dioxygenase (SDO) is proposed to be an important elemental sulfur oxidation enzyme in the genusAcidithiobacillus[43], and the relevant gene of SDO and its enzyme activity have been identified inAt. caldus[41]. Result showed that the genessor andsdo were not identified in heterotrophAcidicaldussp. strain DX-1. Another enzyme directing sulfur oxidation is the cytoplasmic heterodisulfide reductase complex [37], which catalyzes the oxidation of sulfane sulfate (RSSH) to produce sulfite. InAt. caldus, RSSH is the product of thio proteins (RSH) with a sulfur atom in the catalysis of thiosulfate by rhodanese (TST) [41]. Results showed that one ortholog of TST was found inAcidicaldus sp. strain DX-1. Additionally, compared to other heterodisulfide reductase complex reported in previous studies [38, 44, 45], only two copies ofhdrB gene, one cope of each HdrAC subunit, anddsrE gene were distributed in the draft genome ofAcidicaldus sp. strain DX-1 (Table S1). We cannot confirm whether other components were present based on the draft genome alone. The putative HdrA subunit was flavoprotein that contained FAD binding site and conserved residues (C-X-G-X-R-D-X6–8-C-S-X2-C-C) for binding of Fe-S cluster, which shared 72% sequence similarity with that ofAt. ferrooxidans. Additionally, the putative HdrB having one typical cysteine-rich regions was distributed in the Hdr gene cluster, and the HdrC subunit having the 4Fe-4S ferredoxin iron-sulfur binding domain shared 74% and 71% sequence similarities with that ofAt. ferrooxidans, respectively.


### 3.7. Sulfite Oxidation

Another important step in the biological RISC oxidation was sulfite oxidation, which was involved in two different pathways. In these models, sulfite was (i) directly oxidized to sulfate, which was catalyzed by a molybdenum-containing sulfite: acceptor oxidoreductase or (ii) indirectly oxidized by the intermediate adenosine-5′-phosphosulfate (APS) [40]. For indirect sulfite oxidation, sulfite was catalyzed by APS reductase to produce APS and then oxidized to generate sulfate by sulfate adenylyltransferase (SAT). Most studies support the fact that the indirect sulfite oxidation exists in bothAt. ferrooxidansandAt. caldus, which lack gene encoding APS reductase [36, 38]. InAcidicaldus sp. strain DX-1, genes coding for sulfite: acceptor oxidoreductase and APS reductase were not found. However, a putative gene encoding SAT, which has a conservative phosphoadenosine phosphosulfate (PAPS) reductase family domain, was identified inAcidicaldus sp. strain DX-1 (Table S1). Thus, indirect sulfite oxidation might also play an important role inAcidicaldussp. strain DX-1.


### 3.8. Terminal Oxidases in RISC Oxidation

In chemolithoautotrophic sulfur-oxidizing bacteria, the RISC oxidation closely links with electron transfer through terminal oxidases. InAt. ferrooxidans, electrons from RISC oxidation are transferred via the quinol pool (QH2) (i) either directly to terminal oxidasebd orbo
3 to produce a proton gradient or indirectly through abc
1 complex and cytochrome c or a high potential iron-sulfur protein (HiPIP) toaa
3 oxidase where the concentration of O2 is low or (ii) to NADH complex I to generate reducing power [38]. Genomic analysis showed thatAcidicaldus sp. strain DX-1 had one copy ofbd ubiquinol oxidase genes (cydAB), one gene cluster encodingbc
1 complex, and a gene cluster encoding theaa
3 oxidase (Table S2). Genes encoding 14 subunits of NADH complex I were also identified in its genome. However, no candidate genes encoding cytochromebo
3 ubiquinol oxidase were found. Therefore, a putative electron transfer chain inAcidicaldus sp. strain DX-1 was proposed: electrons from SQR, HDR, and TQO were transferred via the QH2 (i) either tobd oxidase orbc
1 complex to produce the proton gradient or (ii) to NADH complex I to generate reducing power. Additionally, electrons inbc
1 complex were probably transferred via cytochrome c to theaa
3 oxidase, where the concentration of O2 is low, to produce the proton gradient.


### 3.9. Expression of Key Genes Involved in RISC Oxidation

Genomic analysis provided evidence that most genes associated with RISC oxidation inAcidicaldussp. strain DX-1 were similar to that in other acidophilic sulfur oxidizers. Furthermore, the expressions of relevant genes in both organic matter and organic matter supplemented with sulfur (S0) medium were validated by real-time quantitative PCR (Figure 2). Results indicated that all selected genes were upregulated in organic matter + S0 culture, suggesting that these encoding proteins are likely involved in the RISC oxidation in this strain. Genes encoding TST and three subunits of HdrABC complex forAcidicaldus strain DX-1 were significantly upregulated, whilehdrB2 gene has relatively lower expression. The conserved family domains of each HdrABC subunit and change of gene expression pattern strongly suggested that the HdrABC complex was concluded to oxidize cytoplasmic sulfur inAcidicaldus sp. strain DX-1. In particular, genes encoding HdrABC complex were not found in the closely related bacteriaAp. cryptum JF-5 andAp. multivorumAIU301, which might be the reason why many strains ofAcidiphiliumspp. could not directly oxidize sulfur. Additionally, compared withsqr gene,sqr2 gene was significantly upregulated, indicating its importance to metabolize the periplasmic sulfide which might be generated from other unknown metabolic pathways.Figure 2
RT-qPCR results depicting the expression difference of targeted genes involved in RISC oxidation inAcidicaldus sp. strain DX-1 (in organic matter supplemented with sulfur (S0) medium or without sulfur). These genes used in this experiment (the encoding proteins and genomic loci are shown in parentheses) includehdrA (pyridine nucleotidedisulfide oxidoreductase; scaffold254: 1-941),hdrB (heterodisulfide reductase subunit B, homolog; scaffold254: 827-2329),hdrC (iron-sulfur cluster-binding protein; scaffold254: 2500-3246),hdrB2 (heterodisulfide reductase subunit B, homolog; scaffold9: 22130-22789),rhd (rhodanese, thiosulfate sulfurtransferase; scaffold153: 15026-15202),doxD (thiosulfate: quinone oxidoreductase; scaffold68: 10659-11696),sat (sulfate adenylyltransferase; scaffold217: 14172-14891),sqr (sulfide quinone reductase; scaffold18: 14834-16117), andsqr2 (sulfide quinone reductase; scaffold92: 4297-5430).


## 3.1. Isolation and Phylogenetic Analysis

Several strains (named DX-1, DX-5, DX-20, and DX-22) from Dexing Copper Mine sampling sites, which showed similar colony and cellular morphologies, were isolated on solid medium with 0.02% yeast extract. All of these strains grew as small, white, round to convex-shaped colonies. Considering that 16S rRNA sequence analysis showed 100% sequence identities, only one (strain DX-1) was selected for further analysis. The rooted phylogenetic tree indicated that 16S rRNA gene of strain DX-1 was assigned into a phyletic cluster withAcidicaldusspp. (Figure 1). Furthermore, 16S rRNA gene of strain DX-1 shared 99% sequence identity with that of type strainAcidicaldus organivoransY008 [13], indicating that strain DX-1 belongs toAcidicaldus species.Figure 1
Phylogenetic tree of 16S rRNA genes showing the relationship of strain DX-1 (inbold) to other acidophilic α-Proteobacteria. Phylogenetic tree was constructed based on an alignment of 1283 bp using the neighbor joining method. The scale bar represents 0.01 nucleotide substitutions per 100. And the database accession numbers of the gene sequences used are given in parentheses.Acidisphaera rubrifaciens HS-AP3 was used to root the tree.


## 3.2. Genome Assembly and Annotation

The draft genome sequence ofAcidicaldus sp. strain DX-1 contained 2,990,377 bp with GC content of 68.48% distributed in 376 scaffolds, which ranged from 1,000 bp to 67,370 bp. Given the high sequence coverage (200-fold), it was likely to identify most genes in the draft genome ofAcidicaldus sp. strain DX-1. Results showed that genome ofAcidicaldus sp. strain DX-1 harbored 3,259 predicted coding sequences (CDSs), which represented 89.83% of the genome, 1 rRNA operon, and 43 tRNA genes (Table 2). The number of CDSs in this strain was less than that inAcidiphilium multivorum AIU301 but more than that inAcidiphilium cryptum JF-5 andAcetobacter pasteurianus.Table 2
General features ofAcidicaldus sp. strain DX-1 in comparison with other heterotrophic acidophiles.
Organism
Acidicaldus sp. strain DX-1
Acidiphilium cryptum JF-5
Acidiphilium multivorumAIU301
Acetobacter pasteurianus
Isolated from
Acid effusion of mill tailing
Acidic Sediments
Pyritic acid mine drainage
Spontaneous cocoa bean heap fermentation
Temperature range
Moderate thermophilic (50~55°C)
Mesophilic (30~35°C)
Mesophilic (30–35°C)
Mesophilic (25–30°C)
pH range
2.5–3.0
3.0–3.5
3.0–3.5
5.4–6.3
Nutrition type
Heterotrophic
Heterotrophic
Heterotrophic
Heterotrophic
Oxygen requirement
Facultatively anaerobic
Facultatively anaerobic
Facultatively anaerobic
Obligately aerobic
Fe(III) reduction
Yes
Yes
Yes
No
Sulfur oxidation
Yes
No
No
No
Genome size (Mb)
2.99
3.39
3.75
2.91
GC content
68.48%
67.99%
67.58%
53.04%
Protein-coding genes
3259
3063
3448
2875
rRNA operon
1
2
2
5
Number of tRNAs
43
47
48
57


## 3.3. Hydrogen Sulfide Oxidation

Phototrophic bacteria such asAllochromatium vinosum and someChlorobium spp. can anaerobically oxidize hydrogen sulfide to sulfur via a flavocytochrome c-sulfide dehydrogenase [34]. In this process, cytochrome c
550 is used as electron acceptor for energy production. However, flavocytochrome c inA. vinosum is not required for phototrophic growth with hydrogen sulfide. Therefore, sulfide: quinone reductase (SQR) is considered to be essential for growth ofA. vinosum with hydrogen sulfide. In chemotrophic bacteria, such as someAcidithiobacillus spp., SQR is also responsible for the oxidation of hydrogen sulfide [35]. InAcidicaldus sp. strain DX-1, no candidate genes encoding flavocytochrome c-sulfide dehydrogenase were identified, while two genes encoding SQR were found (Table S1, see Supplementary Material available online at http://dx.doi.org/10.1155/2016/8137012). One putative SQR shared 64% amino acid similarity with that ofAt. ferrooxidans ATCC 53993, and the other shared 72% amino acid similarity with that ofBradyrhizobium oligotrophicum S58. Both of them had a small NADH binding domain within a larger FAD binding domain. These results suggested that SQR might be responsible for the oxidation of hydrogen sulfide inAcidicaldus sp. strain DX-1.


## 3.4. S4I Pathway

In some strains of the genusAcidithiobacillus, the membrane-bound thiosulfate: quinone oxidoreductase (TQO) and tetrathionate hydrolase (TetH) are responsible for the S4I pathway [36, 37]. InAt. ferrooxidans, TQO which is constituted of subunits DoxA and DoxD catalyzes the conversion of thiosulfate to tetrathionate, and TetH catalyzes the hydrolysis of tetrathionate to generate sulfur, sulfate and thiosulfate. The consecutive reactions catalyzed by TetH and TQO promote the RISC oxidation inAt. ferrooxidans[38]. One ortholog ofdoxDA encoding TQO was identified inAcidicaldus sp. strain DX-1, which shared 61% similarity with that inAt. ferrooxidans (Table S1). Further analysis revealed that the putative TQO had conserved DoxD domain (Pfam: PF04173) and DoxA domain (Pfam: PF07680). Although previous experiments indicated thatAc. organivoransY008 can grow on tetrathionate overlay medium [13], no candidate genes with significant similarity to TetH were found in the draft genome ofAcidicaldus strain DX-1. Future analyses will be necessary to determine whether TQO indeed catalyzes thiosulfate to tetrathionate and TetH gene is on the portion of the missing draft genome ofAcidicaldus sp. strain DX-1.


## 3.5. Sox Pathway

In the alphaproteobacteriumParacoccus pantotrophus, the periplasmic Sox system encoded by 15 genes [39] are comprised of four subunits, including SoxXA, SoxYZ, SoxB, and Sox(CD)2. Each Sox subunit directing its own function was elaborated in previous studies [40]. BothAcidiphilium cryptum andAcidiphilium multivorum have one gene cluster ofsoxXYZABCD to code for four periplasmic proteins. However, no candidate genes encoding Sox system were identified in the genome ofAcidicaldus sp. strain DX-1. Considering the Sox system was absent, S4I pathway was probably the sole way of further oxidation of thiosulfate in periplasm.


## 3.6. Sulfur Oxidation

The cytoplasmic sulfur oxygenase reductase (SOR), which mediates biological oxidation of inorganic sulfur, can catalyze the disproportionation reaction of sulfur to generate sulfite, thiosulfate, and hydrogen sulfide [23, 38, 40–42]. Sulfur dioxygenase (SDO) is proposed to be an important elemental sulfur oxidation enzyme in the genusAcidithiobacillus[43], and the relevant gene of SDO and its enzyme activity have been identified inAt. caldus[41]. Result showed that the genessor andsdo were not identified in heterotrophAcidicaldussp. strain DX-1. Another enzyme directing sulfur oxidation is the cytoplasmic heterodisulfide reductase complex [37], which catalyzes the oxidation of sulfane sulfate (RSSH) to produce sulfite. InAt. caldus, RSSH is the product of thio proteins (RSH) with a sulfur atom in the catalysis of thiosulfate by rhodanese (TST) [41]. Results showed that one ortholog of TST was found inAcidicaldus sp. strain DX-1. Additionally, compared to other heterodisulfide reductase complex reported in previous studies [38, 44, 45], only two copies ofhdrB gene, one cope of each HdrAC subunit, anddsrE gene were distributed in the draft genome ofAcidicaldus sp. strain DX-1 (Table S1). We cannot confirm whether other components were present based on the draft genome alone. The putative HdrA subunit was flavoprotein that contained FAD binding site and conserved residues (C-X-G-X-R-D-X6–8-C-S-X2-C-C) for binding of Fe-S cluster, which shared 72% sequence similarity with that ofAt. ferrooxidans. Additionally, the putative HdrB having one typical cysteine-rich regions was distributed in the Hdr gene cluster, and the HdrC subunit having the 4Fe-4S ferredoxin iron-sulfur binding domain shared 74% and 71% sequence similarities with that ofAt. ferrooxidans, respectively.


## 3.7. Sulfite Oxidation

Another important step in the biological RISC oxidation was sulfite oxidation, which was involved in two different pathways. In these models, sulfite was (i) directly oxidized to sulfate, which was catalyzed by a molybdenum-containing sulfite: acceptor oxidoreductase or (ii) indirectly oxidized by the intermediate adenosine-5′-phosphosulfate (APS) [40]. For indirect sulfite oxidation, sulfite was catalyzed by APS reductase to produce APS and then oxidized to generate sulfate by sulfate adenylyltransferase (SAT). Most studies support the fact that the indirect sulfite oxidation exists in bothAt. ferrooxidansandAt. caldus, which lack gene encoding APS reductase [36, 38]. InAcidicaldus sp. strain DX-1, genes coding for sulfite: acceptor oxidoreductase and APS reductase were not found. However, a putative gene encoding SAT, which has a conservative phosphoadenosine phosphosulfate (PAPS) reductase family domain, was identified inAcidicaldus sp. strain DX-1 (Table S1). Thus, indirect sulfite oxidation might also play an important role inAcidicaldussp. strain DX-1.


## 3.8. Terminal Oxidases in RISC Oxidation

In chemolithoautotrophic sulfur-oxidizing bacteria, the RISC oxidation closely links with electron transfer through terminal oxidases. InAt. ferrooxidans, electrons from RISC oxidation are transferred via the quinol pool (QH2) (i) either directly to terminal oxidasebd orbo
3 to produce a proton gradient or indirectly through abc
1 complex and cytochrome c or a high potential iron-sulfur protein (HiPIP) toaa
3 oxidase where the concentration of O2 is low or (ii) to NADH complex I to generate reducing power [38]. Genomic analysis showed thatAcidicaldus sp. strain DX-1 had one copy ofbd ubiquinol oxidase genes (cydAB), one gene cluster encodingbc
1 complex, and a gene cluster encoding theaa
3 oxidase (Table S2). Genes encoding 14 subunits of NADH complex I were also identified in its genome. However, no candidate genes encoding cytochromebo
3 ubiquinol oxidase were found. Therefore, a putative electron transfer chain inAcidicaldus sp. strain DX-1 was proposed: electrons from SQR, HDR, and TQO were transferred via the QH2 (i) either tobd oxidase orbc
1 complex to produce the proton gradient or (ii) to NADH complex I to generate reducing power. Additionally, electrons inbc
1 complex were probably transferred via cytochrome c to theaa
3 oxidase, where the concentration of O2 is low, to produce the proton gradient.


## 3.9. Expression of Key Genes Involved in RISC Oxidation

Genomic analysis provided evidence that most genes associated with RISC oxidation inAcidicaldussp. strain DX-1 were similar to that in other acidophilic sulfur oxidizers. Furthermore, the expressions of relevant genes in both organic matter and organic matter supplemented with sulfur (S0) medium were validated by real-time quantitative PCR (Figure 2). Results indicated that all selected genes were upregulated in organic matter + S0 culture, suggesting that these encoding proteins are likely involved in the RISC oxidation in this strain. Genes encoding TST and three subunits of HdrABC complex forAcidicaldus strain DX-1 were significantly upregulated, whilehdrB2 gene has relatively lower expression. The conserved family domains of each HdrABC subunit and change of gene expression pattern strongly suggested that the HdrABC complex was concluded to oxidize cytoplasmic sulfur inAcidicaldus sp. strain DX-1. In particular, genes encoding HdrABC complex were not found in the closely related bacteriaAp. cryptum JF-5 andAp. multivorumAIU301, which might be the reason why many strains ofAcidiphiliumspp. could not directly oxidize sulfur. Additionally, compared withsqr gene,sqr2 gene was significantly upregulated, indicating its importance to metabolize the periplasmic sulfide which might be generated from other unknown metabolic pathways.Figure 2
RT-qPCR results depicting the expression difference of targeted genes involved in RISC oxidation inAcidicaldus sp. strain DX-1 (in organic matter supplemented with sulfur (S0) medium or without sulfur). These genes used in this experiment (the encoding proteins and genomic loci are shown in parentheses) includehdrA (pyridine nucleotidedisulfide oxidoreductase; scaffold254: 1-941),hdrB (heterodisulfide reductase subunit B, homolog; scaffold254: 827-2329),hdrC (iron-sulfur cluster-binding protein; scaffold254: 2500-3246),hdrB2 (heterodisulfide reductase subunit B, homolog; scaffold9: 22130-22789),rhd (rhodanese, thiosulfate sulfurtransferase; scaffold153: 15026-15202),doxD (thiosulfate: quinone oxidoreductase; scaffold68: 10659-11696),sat (sulfate adenylyltransferase; scaffold217: 14172-14891),sqr (sulfide quinone reductase; scaffold18: 14834-16117), andsqr2 (sulfide quinone reductase; scaffold92: 4297-5430).


## 4. Conclusions

In acidophilic microbial communities, heterotrophic sulfur-oxidizing bacteria have versatile functions, which not only can utilize numerous organic matters as energy sources but can also acquire energy from the oxidation of RISC. In this study, we isolated a heterotrophic sulfur-oxidizing bacterium,Acidicaldus sp. strain DX-1, from the effusion of mill tailings in Dexing Copper Mine. Comparative genomics and RT-qPCR analysis revealed that several genes associated with RISC oxidation inAcidicaldussp. strain DX-1 were similar to other acidophilic sulfur-oxidizing bacteria. Based on the aforementioned analyses, a model for RISC oxidation inAcidicaldus sp. strain DX-1 was hypothesized (Figure 3), which might provide new insights and guides for the RISC oxidation of heterotrophs in the future.Figure 3
Genome-based model of RISC oxidation inAcidicaldus sp. strain DX-1. The figure was adapted from previous sulfur oxidation models [23, 38]. RISC oxidation involves various enzymes and a number of electron carriers. The solid lines denote the oxidation of RISC catalyzed by various enzymes, and the dashed lines denote the direction of electron transfer. SQR, sulfide quinone reductase; TQO, thiosulfate quinone oxidoreductase; TST, rhodanese; HdrABC, heterodisulfide reductase complex; SAT, sulfate adenylyltransferase;bd, terminal oxidasebd;bc
1, terminal oxidasebc
1 complex; NADH, NADH complex I;aa
3, cytochrome c oxidaseaa
3-type; QH2, quinol pool; and CycA, cytochromec.



---

*Source: 8137012-2016-04-28.xml*