Synechococcus sp. Strain PCC7002 Uses Sulfide:Quinone Oxidoreductase To Detoxify Exogenous Sulfide and To Convert Endogenous Sulfide to Cellular Sulfane Sulfur.

Cyanobacteria are a major force for primary production via oxygenic photosynthesis in the ocean. A marine cyanobacterium, PCC7002, is actively involved in sulfide metabolism. It uses SQR to detoxify exogenous sulfide, enabling it to survive better than its Δsqr mutant in sulfide-rich environments. PCC7002 also uses SQR to oxidize endogenously generated sulfide to S0, which is required for the proper expression of key genes involved in photosynthesis. Thus, SQR has at least two physiological functions in PCC7002. The observation provides a new perspective for the interplays of C and S cycles.

The toxic effect of sulfide is well known, inhibiting respiration by acting on cytochrome c oxidase in heterotrophic bacteria (12)(13)(14)(15) and photosynthesis by binding to metalloproteins of photosynthesis system II (PSII) (16)(17)(18)(19). Sulfide can reach high concentrations in specific habitats, such as hydrothermal vents and seeps (10) and coastal mudflats (11). Sulfide is mainly derived from sulfate as the terminal electron acceptor for organic mineralization by sulfur-reducing bacteria in OMZs (10).
The cyanobacteria play a vital role in global primary production as the most ancient and abundant phytoplankton in the ocean (39,40). Cyanobacteria possess both photosynthesis system I (PSI) and photosynthesis system II (PSII), and they usually perform oxygenic photosynthesis but can do anoxygenic photosynthesis (41,42). They are effective participants in marine carbon cycles (43). Cyanobacteria are found to inhabit multiple marine habitats, and Synechococcus spp. dominate the picocyanobacterial communities in eutrophic coastal and mesotrophic open ocean waters (44,45). In some habitats, they will encounter sulfide. The sulfur cycle is generally coupled with the carbon cycle (46)(47)(48); however, it is still unclear how sulfide affects Synechococcus spp. and how they deal with it.
Here, we report that Synechococcus sp. strain PCC7002 (PCC7002) used SQR to oxidize self-produced sulfide to polysulfide, maintaining a relatively high level of cellular sulfane sulfur, which offers growth advantages to the wild type over the Δsqr mutant. Further, PCC7002 also survived better in sulfide-rich environments via the detoxification role of SQR and PDO. The two enzymes collectively oxidize sulfide to sulfite and thiosulfate, resembling the newly reported pathway in heterotrophic bacteria (25).
The Δsqr mutant contained a sharply decreased level of cellular sulfane sulfur in comparison with the wild type and PCC7002 Δsqr::sqr (Fig. 1B). The results indicate that PCC7002 uses SQR to oxidize endogenously produced sulfide, allowing the bacterium to keep a relatively high level of cellular sulfane sulfur.
PCC7002 has a competitive advantage over the ⌬sqr mutant in coculture. To investigate the effect of SQR, we monitored the growth and photosynthesis of the wild type and the Δsqr mutant. The deletion of sqr had no effect on the growth of PCC7002 (Fig. S2). We then cocultured PCC7002 and PCC7002Δsqr to access their competitive capacity. The relative abundance of PCC7002 and PCC7002Δsqr in the coculturing system was detected by real-time quantitative reverse transcription-PCR (qRT-PCR) with specific primer pairs (Table S2). Different ratios of PCC7002 and PCC7002Δsqr were mixed, and the qRT-PCR results correlated with the ratios accordingly (Fig. S3), confirming the approach. PCC7002 and PCC7002Δsqr showed little difference in the coculture system during the initial 48 h; however, the advantage of the wild type became apparent after continuous transfer (Fig. 1C). We also cultured PCC7002 and PCC7002Δsqr separately using the same method, but results showed no difference between the wild type and the mutant (Fig. S4).
The effects of sqr inactivation on the photosynthesis process. The photosynthesis parameters of the wild type and the Δsqr mutant were tested. The Δsqr mutant had a higher oxygen evolution rate than that of the wild type ( Fig. 2A). The increased oxygen evolution rate upon the Δsqr mutant was consistent with an acceleration of H 2 O oxidation to make up for the loss of H 2 S oxidation. The maximal PSII quantum yield The effect of SQR on H 2 S accumulation, cellular sulfane sulfur, and competitiveness of PCC7002 with its mutant. (A) H 2 S accumulation during culture as detected in the gas phase with lead-acetate paper. 1, PCC7002; 2, PCC7002Δsqr; 3, PCC7002Δsqr::sqr. (B) Sulfane sulfur contents of PCC7002 (black bar), PCC7002Δsqr (red bar), and PCC7002Δsqr::sqr (blue bar)were detected by using SSP4 (sulfane sulfur probe 4). (C) Competition experiment of PCC7002 and PCC7002Δsqr. The abundance of PCC7002 (black bars)showed an apparent advantage over PCC7002Δsqr (red bars) after transfer 4 times. The mixing of two strains at equal cell numbers was defined as the first transfer. The time interval between each transfer was 48 h, and real-time quantitative PCR (qPCR) was done at 48 h. All data are averages from three samples with standard deviation (error bar). The experiment was repeated at least three times. ****, P Ͻ 0.0001 (paired t test).

Functions of SQR in Cyanobacteria
® was calculated as Fv/Fm, a measure of the conversion efficiency of the light energy of the PSII reaction center. Here, the Fv/Fm value showed little difference between the wild type and the mutant (Fig. 2B). The relative electron transport rate (rETR) increased in the mutant (Fig. 2C). Thus, photosynthesis efficiency was increased in the Δsqr mutant. However, the apparently increased efficiency in photosynthesis did not affect the growth (Fig. S2). We then checked the expression levels of major genes involved in photosynthesis in the wild type and the Δsqr mutant. PCC7002 had three psbA genes, encoding the isoforms of D1 protein, a key protein of the PSII core, and the isoforms allowed cyanobacteria to adapt to different light intensities. The expression levels of psbA1, psbA2, and psbA3 were all upregulated (1.5-to 2-fold) in the Δsqr mutant (Fig. S5). The rbcX, rbcL, and rbcS genes encoding RuBisCO and the tkt gene encoding transketolase of the Calvin-Benson-Bassham (CBB) cycle all showed apparent transcriptional increases (Ͼ2-fold) in the Δsqr mutant in comparison with the wild type (Fig. 2D). The mRNA level of petA encoding apocytochrome f precursor increased by 1.8-fold, and the transcripts of petC encoding a Rieske-FeS protein increased by 1.5-fold in the Δsqr mutant. However, the transcripts of the petB gene encoding cytochrome b 6 showed little change (Fig. S5). Overall, the qRT-PCR results indicated that the inactivation of sqr in PCC7002 affected many physiological processes, including PSII, the CBB cycle, and the photosynthetic electron transport chain ( Fig. 2D and Fig. S5).
SQR plays a detoxification role in PCC7002. We treated the wild type, the Δsqr mutant, and the complementation strain with different concentrations of NaHS for 6 h and then washed, diluted, and inoculated cells onto the Aϩ medium for incubation under 50 mol photons m Ϫ2 s Ϫ1 at 30°C. Without sulfide, all the three strains grew well The Δsqr mutant (red squares)showed a higher oxygen evolution rate than the wild type (black circles). PAR, photosynthetically active radiation. ***, P Ͻ 0.001 (two-way analysis of variance [ANOVA]). (B) The Fv/Fm value of the Δsqr mutant (red bar)showed a 16% increase over that of the wild type (black bar). (C) The relative electron transport rate (rETR) of the Δsqr mutant (red squares)increased over that of the wild type (black circles). ***, P Ͻ 0.001 (two-way ANOVA). (D) The rbcS, rbcL, rbcX, and tkt genes in PCC7002Δsqr (red bars)were all upregulated compared with the wild type (black bars)as detected by qRT-PCR. All data are averages from six samples with standard deviation (error bar). The experiment was repeated at least three times. ****, P Ͻ 0.0001 (paired t test). on the Aϩ agar plate. After treatment with 3 mM NaHS, the Δsqr mutant grew worse than the wild type and the complementation strain (Fig. 3A). The growth of the Δsqr mutant was fully inhibited while the wild type and the complementation strain were partly inhibited after treatment with 4 mM and 5 mM NaHS.
Whether PCC7002 oxidized exogenous sulfide was also tested. Cells at log phase of growth (optical density at 730 nm [OD 730 ] ϭ 1) were harvested and resuspended in 50 mM Tris buffer (pH 8.0) at an OD 730 of 8. Without induction, all the three strains (PCC7002, PCC7002Δsqr, and PCC7002Δsqr::sqr) showed no oxidation activity (data not shown). After induction with 2 mM NaHS for 6 h, cell suspensions of the wild type and the complementation strain showed apparent activity with about 500 M NaHS being fully oxidized in 8 h (Fig. 3B). However, the Δsqr mutant showed a similar oxidation rate as the buffer control (Fig. 3B). Thus, SQR plays a detoxification role in PCC7002 by oxidizing excessive sulfide.
Recombinant Escherichia coli with SQR and PDO from PCC7002 oxidized sulfide to sulfite and thiosulfate. PCC7002 also has a persulfide dioxygenase (PDO). We constructed recombinant strains E. coli(pBBR5-sqr), containing only sqr, and E. coli  The presence of sqr and pdo in cyanobacteria. The previously reported SQRs and PDOs were used to search for potential SQRs and PDOs from 130 sequenced cyanobacterial genomes (NCBI database, updated to 17 June 2019) (21). Sixty-seven SQRs from 59 cyanobacterial genomes and 101 PDOs from 85 cyanobacterial genomes were identified after checking via phylogenetic analysis; 54 genomes possessed both SQRs and PDOs (Fig. 5). All the cyanobacterial SQRs belonged to the type I and type II SQRs (Fig. S6).

DISCUSSION
Our data indicate two physiological roles of SQR in PCC7002 (Fig. 6). First, SQR oxidizes endogenous H 2 S to polysulfide, which makes PCC7002 more competitive in coculture of the wild type and the Δsqr mutant (Fig. 1). It is likely that polysulfide is   (50). Recently, it was shown that the signaling role of H 2 S is via sulfane sulfur, which is sensed by gene regulators that activate the genes involved in H 2 S oxidation in bacteria (51)(52)(53). OxyR also senses high levels of cellular sulfane sulfur to produce enzymes for its removal (54). Our evidence suggests that sulfane sulfur is likely involved in regulating photosynthetic genes ( Fig. 2D; see also Second, SQR plays a detoxification role. Sulfide may accumulate at high concentrations in OMZs (10). With SQR, PCC7002 survives better than the Δsqr mutant in the presence of high concentrations of sulfide (Fig. 3A). The detoxification also requires PDO that prevents the excessive accumulation of polysulfide, and the sulfide oxidation pathway is the same as in heterotrophic bacteria (25). Bioinformatics analysis shows that about 45% of the sequenced cyanobacterial genomes contain SQR (Fig. 5), while only 21% of the sequenced bacterial genomes contain SQR (21), indicating the wide distribution and importance of SQR in cyanobacteria. Thus, SQR may facilitate PCC7002 and possibly other cyanobacteria to adapt to changing environments and increase their growth fitness.
SQR may have a third role in cyanobacteria, allowing them to use H 2 S as the electron donor for photosynthesis. In the sulfur photosynthetic bacteria, SQR oxidizes H 2 S and passes the electrons to drive anaerobic photosynthesis (52). Some cyanobacteria also perform anoxygenic photosynthesis in sulfide-rich environments (18, 31, 41), partly FIG 6 SQR helps PCC7002 deal with high sulfide and improves its growth fitness. On the one hand, SQR oxidizes exogenous H 2 S to sulfane sulfur, which spontaneously reacts with GSH and generates GSSH. Then, PDO oxidizes GSSH to sulfite and generates a new GSH. At last, sulfite spontaneously reacts with sulfane sulfur to generate thiosulfate, which is transported out of the cell. Hence, SQR coupling with PDO plays a detoxification role by oxidizing H 2 S to sulfite and thiosulfate. On the other hand, SQR oxidizes endogenous H 2 S to sulfane sulfur, which may act as a signal molecule and participates in the regulation of photosynthesis, especially in the CBB cycle (the proposed signaling function of sulfane sulfur is indicated with a dashed black line). Meanwhile, H 2 S oxidation may donate to the photosynthetic electron transport (proposed electron transfer from SQR to plastoquinone [PQ] is indicated with a dashed gray line). As a result, SQR enables PCC7002 to increase its growth fitness as shown when culturing the wild type and the mutant in a mixture.

Functions of SQR in Cyanobacteria
® because sulfide is a better electron donor and partly because sulfide at high concentrations inhibits photosystem II (16). We observed that the Δsqr mutant had a higher oxygen evolution rate than that of the wild type ( Fig. 2A). The observation may be a result of increased expression of photosynthetic genes or the oxidation of endogenously produced sulfide by SQR. Further investigation is needed to clarify whether sulfide provides electrons for photosynthesis under aerobic conditions. The physiological roles of SQR reported here may guide the investigation of cyanobacteria in ecological niches with varied O 2 and sulfide contents, as many sequenced cyanobacterial genomes contain SQR (Fig. 5).
The source of endogenous H 2 S in PCC7002 is likely from the metabolism of cysteine or homocysteine, similar to H 2 S production by heterotrophic bacteria during normal growth (21), or from assimilatory sulfate reduction (55). The endogenous H 2 S has been shown to act as a signal molecule at low concentrations (12). Cellular sulfane sulfur, including polysulfide and GSSH, is believed to be the actual mediator of H 2 S signaling (56)(57)(58)(59); sulfane sulfur can also protect Escherichia coli cells from oxidative stress (35).
When SQR is inactivated in PCC7002, an interesting phenomenon is that the growth advantage of PCC7002 occurs only in the coculture process ( Fig. 1C and Fig. S2). This may be caused by the intracellular polysulfide coming from endogenous H 2 S oxidation by SQR. The mutant contains less sulfane sulfur (Fig. 1B) but has increased oxygen evolution rates and rETRs with elevated expression of photosynthesis genes (Fig. 2). The gene expression variation between the mutant and the wild type may also be a result of the lowered sulfane sulfur in the mutant, as sulfane sulfur is a known signaling molecule participating in many physiological processes in eukaryotes (58). And from our results, we infer that carbon and sulfur metabolisms can be linked via sulfane sulfur just like carbon and nitrogen metabolisms in cyanobacteria (60,61). As we know, sulfane sulfur can perform as a reductant to help relieve oxidative stress (35,37). Here, increased photosynthesis efficiency with lower sulfane sulfur content in the Δsqr mutant may cause oxidative stress (62), thus making it less competitive in coculture. However, it is unclear how polysulfide regulates key genes involved in photosynthesis processes, such as PSII, the CBB cycle, and photosynthesis electron transport.
In summary, SQR helps PCC7002 to convert endogenous H 2 S to polysulfide, which is used to maintain a relatively high level of cellular sulfane sulfur, and to detoxify exogenous H 2 S. The relatively high cellular sulfane sulfur offers the wild type a growth advantage over the Δsqr mutant in a mixed culture. The detoxification role of SQR may enable PCC7002 to survive better in OMZs, thus helping relieve deoxygenation by generating more O 2 . Thus, sulfide oxidation plays a vital role in the adaption process of PCC7002 and possibly other cyanobacteria. We present a new perspective for understanding the coupling of carbon and sulfur cycles, which are key components of the marine biogeochemical cycle.

MATERIALS AND METHODS
Strains and culture conditions. PCC7002 cultures were grown in conical flasks containing medium A, supplemented with 1 mg of NaNO 3 ml Ϫ1 (designated medium Aϩ) (63) under continuous illumination of 50 mol photons m Ϫ2 s Ϫ1 at 30°C. Glycerol (10 mM) was added as a supplement in the Aϩ medium to serve as the carbon and energy source. For mutant strains, the concentrations of appropriate antibiotics required by the mutant strains were as following: 100 g/ml kanamycin and 20 g/ml chloramphenicol. E. coli was cultured in LB medium at 37°C. The strains and plasmids used in this paper are all listed in Table S1 in the supplemental material.
Strain construction. The Δsqr mutant of PCC7002 was constructed by homologous recombination. Two gene segments immediately upstream and downstream of the sqr gene (CyanoBase: SYNPCC7002_G0075), about 1,000 bp long, were amplified by PCR from genomic DNA of the wild-type strain with the primer sets sqr-del-1/sqr-del-2 and sqr-del-5/sqr-del-6 (Table S2). A kanamycin resistance cartridge was excised from the plasmid pET30a with the primers sqr-del-3/sqr-del-4. All the above three fragments were joined via fusion PCR, with the kanamycin resistance cartridge in the middle. Then, the fusion fragment was cloned into pJET1.2 blunt vector (Thermo, Beijing, China). The sqr gene was inserted at neutral site 1 (CyanoBase: SYNPCC7002_A0933) of PCC7002 to form the complementation strain (64). Wild-type cells were transformed with these constructs, and transformants were selected using appropriate antibiotics. The primers used in this paper are all listed in Table S2.
The sqr gene and pdo gene (CyanoBase: SYNPCC7002_A2866) of PCC7002 were also cloned and expressed in E. coli BL21(DE3). The fragment of sqr was amplified from the PCC7002 genome using primers pbr-sqr-F and pbr-sqr-R (listed in Table S1) containing 20-bp extensions overlapping the adjacent fragment. The sqr fragment was cloned into the plasmid pBBR1mcs-5 (pBBR5) by using the TEDA assembly (65). The plasmid pBBR5-sqr-pdo containing both sqr and pdo was constructed using the same protocol with primers listed in Table S1. The two plasmids pBBR5-sqr and pBBR5-sqr-pdo were transformed to E. coli BL21(DE3).
Sulfide accumulation detection. The method used for the detection of sulfide accumulation was the same as described previously (65,66). PCC7002 cultures were transferred into 3 ml of Aϩ medium in a 15-ml glass tube and incubated with shaking for 96 h with lead acetate [Pb(Ac) 2 ] paper strips at the top of the tube, which would turn black in the presence of sulfide as the production of PbS black precipitates.
Toxicity analysis of sulfide. PCC7002 and its mutants at log phase at an OD 730 of 0.6 to 0.7 were treated with 3 mM, 4 mM, or 5 mM NaHS for 6 h in the sealed tubes. After incubation, cells were washed and resuspended in fresh Aϩ medium. The cells were diluted with Aϩ medium at an OD 730 of 0.05, and 10 l was placed on the Aϩ agar plate. The differences between PCC7002 and its mutants appeared after cultivation at 30°C under continuous illumination of 50 mol photons · m Ϫ2 · s Ϫ1 for 7 days.
Sulfide oxidation and product analysis. Sulfide oxidation by PCC7002 and its mutants was assayed. Cells were harvested after induction with 2 mM NaHS for 6 h at log phase. They were resuspended in 50 mM Tris buffer (pH 8.0) at an OD 730 of 8 after being washed twice with the same buffer. The reaction was performed in 30-ml serum bottles sealed with a rubber stopper to minimize the loss of sulfide. NaHS (0.5 mM) was added to the bottle to initiate the reaction.
Sulfide oxidation by recombinant E.coli was similar, but cells were harvested and resuspended to an OD 600 of 2. The concentrations of sulfide, polysulfide, sulfite, and thiosulfate were analyzed as described previously (25,66). Briefly, sulfide was detected by a colorimetric method; polysulfide, sulfite, and thiosulfate were derivatized with monobromobimane and then measured by high-performance liquid chromatography (HPLC) with a fluorescence detector.
Competition experiment. PCC7002 and PCC7002Δsqr at log phase with an OD 730 of 0.6 to 0.7 were transferred to fresh Aϩ medium to a final concentration with an OD 730 of 0.05 each as a mixture. The total genomic DNA of the mixed culture was extracted and used as the template for real-time quantitative PCR (qPCR) analysis after culturing for 48 h. The abundance of the sqr gene was used as the marker for the relative abundance of PCC7002, and the abundance of the kan gene was used for PCC7002Δsqr, in which the sqr gene was replaced by the kan gene. Then, the mixed culture was transferred to fresh Aϩ medium at an OD 730 of 0.05 for another 48 h and analysis. The mixing of two strains at equal cell numbers was defined as the first transfer. The mixture was transferred 4 times.
Endogenous sulfane sulfur analysis. Sulfane sulfur probe 4 (SSP4) was used to detect endogenous sulfane sulfur as previously described (35,49). Cells were harvested, washed, and resuspended in 50 mM Tris buffer (pH 8.0). Cells were disrupted by ultrasonication, and debris was removed by centrifugation at 12,500 ϫ g for 10 min. SSP4 (10 M) was added to the supernatant and incubated in the dark at 30°C for 15 min. The fluorescence resulting from the SSP4 reaction with sulfane sulfur was detected with excitation of 482 nm and emission of 515 nm by using a Synergy H1 microplate reader. The fluorescence values were expressed as the intensity per milligram of protein.
Oxygen evolution and chlorophyll fluorescence assay. Cells at log phase with an OD 730 of 0.6 to 0.7 were collected for oxygen evolution measurement via a Clark oxygen electrode (Chlorolab 2ϩ; Hansatech, Norfolk, United Kingdom). The samples were illuminated under continuous stirring (900 rpm) at 30°C with increasing light intensities (from 0 to 1,000 mol photons m Ϫ2 s Ϫ1 ). The rate of oxygen evolution was recorded continuously for 2 min at each light level. The oxygen evolution rate was then normalized to OD 730 .
Chlorophyll fluorescence was measured by a lightweight, hand-held fluorometer (AP-C 100; AquaPen, Drasov, Czech Republic). Cells at log phase with an OD 730 of 0.6 to 0.7 were collected and then adapted to darkness for half an hour. Fv/Fm and rETR were all automatically detected by AP-C 100.
RNA extraction and qRT-PCR analysis. Cells at log phase with an OD 73 of 0.6 to 0.7 were harvested by centrifugation at 10,000 ϫ g, 4°C, for 10 min. Total RNA was isolated by using the TaKaRa MiniBEST universal RNA extraction kit, and the concentration of RNA was verified by Qubit 4 (Thermo Fisher). The cDNA was acquired by using the Prime Script RT reagent kit with genomic DNA (gDNA) eraser (TaKaRa, Beijing, China). The SYBR Premix Ex Taq II kit (TaKaRa) was used for qRT-PCR, and the reactions were run in a Light Cycler 480 II sequence detection system (Roche, Shanghai, China). Primers for target genes are all shown in Table S1, and rnpA (SYNPCC7002_A0989) was used as the reference gene (67). The results were analyzed according to the threshold cycle (2 ϪΔΔCT ) method (68).
Bioinformatics analysis. One hundred thirty cyanobacterial genomes were downloaded from the NCBI database (update to 17 June 2019). The query sequences of SQR and PDO were based on our previous work (21). The SQR and PDO candidates were obtained by searching the database with the standalone BLASTP algorithm, using conventional criteria (E value of Յ1e Ϫ5 , coverage of Ն45%, and identity of Ն30%). The candidates were analyzed by using ClustalW for alignment and MEGA 7.0 for neighbor-joining tree building with the following parameters: pairwise deletion, p-distance distribution, and bootstrap analysis of 1,000 repeats (69). The candidates in the same clade as the seeds were selected. The published flavocytochrome c sulfide dehydrogenase (FCSD) sequences were used as the outgroup of SQR, and the GloB (GloB hydroxyacylglutathione hydrolase) sequences were used as the outgroup of PDO.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only.