Characteristics and function of sulfur dioxygenase in Echiuran worm Urechis unicinctus

Identification and validation of the reference genes in the echiuran worm Urechis unicinctus based on transcriptome data

BACKGROUND

Real-time quantitative PCR (RT-qPCR) is a crucial and widely used method for gene expression analysis. Selecting suitable reference genes is extremely important for the accuracy of RT-qPCR results. Commonly used reference genes are not always stable in various organisms or under different environmental conditions. With the increasing application of high-throughput sequencing, transcriptome analysis has become an effective method for identifying novel stable reference genes.

RESULTS

In this study, we identified candidate reference genes based on transcriptome data covering embryos and larvae of early development, normal adult tissues, and the hindgut under sulfide stress using the coefficient of variation (CV) method in the echiuran Urechis unicinctus, resulting in 6834 (15.82%), 7110 (16.85%) and 13880 (35.87%) candidate reference genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the candidate reference genes were significantly enriched in cellular metabolic process, protein metabolic process and ribosome in early development and normal adult tissues as well as in cellular localization and endocytosis in the hindgut under sulfide stress. Subsequently, ten genes including five new candidate reference genes and five commonly used reference genes, were validated by RT-qPCR. The expression stability of the ten genes was analyzed using four methods (geNorm, NormFinder, BestKeeper, and ∆Ct). The comprehensive results indicated that the new candidate reference genes were more stable than most commonly used reference genes. The commonly used ACTB was the most unstable gene. The candidate reference genes STX12, EHMT1, and LYAG were the most stable genes in early development, normal adult tissues, and hindgut under sulfide stress, respectively. The log₂(TPM) of the transcriptome data was significantly negatively correlated with the Ct values of RT-qPCR (Ct =  - 0.5405 log₂(TPM) + 34.51), which made it possible to estimate the Ct value before RT-qPCR using transcriptome data.

CONCLUSION

Our study is the first to select reference genes for RT-qPCR from transcriptome data in Echiura and provides important information for future gene expression studies in U. unicinctus.

Microbial carbon, sulfur, iron, and nitrogen cycling linked to the potential remediation of a meromictic acidic pit lake

Cueva de la Mora is a permanently stratified acidic pit lake and a model system for extreme acid mine drainage (AMD) studies. Using a combination of amplicon sequencing, metagenomics and metatranscriptomics we performed a taxonomically resolved analysis of microbial contributions to carbon, sulfur, iron, and nitrogen cycling. We found that active green alga Coccomyxa onubensis dominated the upper layer and chemocline. The chemocline had activity for iron(II) oxidation carried out by populations of Ca. Acidulodesulfobacterium, Ferrovum, Leptospirillium, and Armatimonadetes. Predicted activity for iron(III) reduction was only detected in the deep layer affiliated with Proteobacteria. Activity for dissimilatory nitrogen cycling including nitrogen fixation and nitrate reduction was primarily predicted in the chemocline. Heterotrophic archaeal populations with predicted activity for sulfide oxidation related to uncultured Thermoplasmatales dominated in the deep layer. Abundant sulfate-reducing Desulfomonile and Ca. Acidulodesulfobacterium populations were active in the chemocline. In the deep layer, uncultured populations from the bacterial phyla Actinobacteria, Chloroflexi, and Nitrospirae contributed to both sulfate reduction and sulfide oxidation. Based on this information we evaluated the potential for sulfide mineral precipitation in the deep layer as a tool for remediation. We argue that sulfide precipitation is not limited by microbial genetic potential but rather by the quantity and quality of organic carbon reaching the deep layer as well as by oxygen additions to the groundwater enabling sulfur oxidation. Addition of organic carbon and elemental sulfur should stimulate sulfate reduction and limit reoxidation of sulfide minerals.

Genome-Wide Analyses of Heat Shock Protein Superfamily Provide New Insights on Adaptation to Sulfide-Rich Environments in Urechis unicinctus (Annelida, Echiura)

The intertidal zone is a transitional area of the land-sea continuum, in which physical and chemical properties vary during the tidal cycle and highly toxic sulfides are rich in sediments due to the dynamic regimes. As a typical species thriving in this habitat, Urechis unicinctus presents strong sulfide tolerance and is expected to be a model species for sulfide stress research. Heat shock proteins (HSPs) consist of a large group of highly conserved molecular chaperones, which play important roles in stress responses. In this study, we systematically analyzed the composition and expression of HSPs in U. unicinctus. A total of eighty-six HSP genes from seven families were identified, in which two families, including sHSP and HSP70, showed moderate expansion, and this variation may be related to the benthic habitat of the intertidal zone. Furthermore, expression analysis revealed that almost all the HSP genes in U. unicinctus were significantly induced under sulfide stress, suggesting that they may be involved in sulfide stress response. Weighted gene co-expression network analysis (WGCNA) showed that 12 HSPs, including 5 sHSP and 4 HSP70 family genes, were highly correlated with the sulfide stress response which was distributed in steelblue and green modules. Our data indicate that HSPs, especially sHSP and HSP70 families, may play significant roles in response to sulfide stress in U. unicinctus. This systematic analysis provides valuable information for further understanding of the function of the HSP gene family for sulfide adaptation in U. unicinctus and contributes a better understanding of the species adaptation strategies of marine benthos in the intertidal zone.

The critical roles of mitochondrial alternative chains in juvenile ark shells (Anadara broughtonii) exposed to acute hypoxia with or without sulfide

Reduced oxygen levels and increased sulfide concentrations have become a concern for marine animals. This study examines respiratory and energetic adaption to acute (0-96 h) hypoxia (0.5 mg/L dissolved oxygen) with or without sulfide (0.2 mM, 1 mM, 3 mM) in the hypoxia-resistant and sulfide-tolerant ark shell, Anadara broughtonii. The different states of aerobic respiration, energy-balance, and activity of the mitochondrial sulfide oxidation chain (MSOC) under these conditions were evaluated. The results indicated that the anaerobic pathway was activated by hypoxia at 24 h without sulfide, but was activated in the presence of sulfide at only 2 h. Exposure to sulfide resulted in significant accumulation of ATP, probably due to the activated MSOC and lowered metabolism via suppression of Na⁺-K⁺ ATPase activity and protein synthesis. During hypoxia, both enzyme activity and mRNA levels of alternative oxidase (AOX) increased while the key enzymes in MSOC, sulfide: quinone oxidoreductase (SQR) and sulfur dioxygenase (SDO), were not altered. With additional sulfide, the enzyme activity and mRNA levels of AOX, SQR, SDO significantly increased. Classical aerobic respiration was significantly inhibited, and induction of alternative respiration was detected. The corresponding alternative electron transport chain (AETC) accepted the electrons originating from both the tricarboxylic acid cycle and MSOC during the challenge, indicating that the capacity of aerobic respiration and sulfide-oxidation under a reduction state might greatly depend on AETC. The synergistically induced alternative chains (AETC and MSOC) and anaerobic pathway suggested energy-balance between respiration and sulfide-oxidation, which might contribute to the endurance of ark shells to acute sulfide exposure.

Analysis of Raman Spectra by Using Deep Learning Methods in the Identification of Marine Pathogens

The need for efficient and accurate identification of pathogens in seafood and the environment has become increasingly urgent, given the current global pandemic. Traditional methods are not only time consuming but also lead to sample wastage. Here, we have proposed two new methods that involve Raman spectroscopy combined with a long short-term memory (LSTM) neural network and compared them with a method using a normal convolutional neural network (CNN). We used eight strains isolated from the marine organism Urechis unicinctus, including four kinds of pathogens. After the models were configured and trained, the LSTM methods that we proposed achieved average isolation-level accuracies exceeding 94%, not only meeting the requirement for identification but also indicating that the proposed methods were faster and more accurate than the normal CNN models. Finally, through a computational approach, we designed a loss function to explore the mechanism reflected by the Raman data, finding the Raman segments that most likely exhibited the characteristics of nucleic acids. These novel experimental results provide insights for developing additional deep learning methods to accurately analyze complex Raman data.

Defense responses of sulfur dioxygenase to sulfide stress in the razor clam Sinonovacula constricta

BACKGROUND

Sulfide is a well-known toxicant widely distributed in the culture environment. As a representative burrowing benthic bivalve, the razor clam Sinonovacula constricta is highly sulfide tolerant. Mitochondrial sulfide oxidation is an important way for sulfide detoxification, where sulfur dioxygenase (SDO) is the second key enzyme.

OBJECTIVE

To investigate the mechanism of sulfide tolerance in S. constricta, the molecular characterization of its SDO (designated as ScSDO) was studied.

METHODS

The cDNA sequence of ScSDO was cloned by RACE technique. The response of ScSDO in gills and livers of S. constricta was investigated during sulfide exposure (50, 150, and 300 μM sulfide) for 0, 3, 6, 12, 24, 48, 72, and 96 h by qRT-PCR. Moreover, the temporal expression of ScSDO protein in S. constricta gills after exposure to 150 μM sulfide was detected by Western blot. The subcellular location of ScSDO was identified by TargetP 1.1 prediction and Western Blot analysis.

RESULTS

The full-length cDNA of ScSDO was 2914 bp, encoding a protein of 304 amino acids. The deduced ScSDO protein was highly conserved, containing the signature HXHXDH motif of the metallo-β-lactamase superfamily and two metal-binding sites, of which metal-binding site I is known to be the catalytically active center. Subcellular localization confirmed that ScSDO was located only in the mitochondria. Responding to the sulfide exposure, distinct time-dependent increases in ScSDO expression were detected at both mRNA and protein levels. Moreover, the gills exhibited a higher ScSDO expression level than the livers.

CONCLUSIONS

All of our results suggest that ScSDO plays an important role in mitochondrial sulfide oxidation during sulfide stress, making S. constricta highly sulfide tolerant. In addition, as a respiratory tissue, the gills play a more critical role in sulfide detoxification.

New Insights Into Acidithiobacillus thiooxidans Sulfur Metabolism Through Coupled Gene Expression, Solution Chemistry, Microscopy, and Spectroscopy Analyses

Here, we experimentally expand understanding of the reactions and enzymes involved in Acidithiobacillus thiooxidans ATCC 19377 S0 andS 2 ⁢ O 3 2 - metabolism by developing models that integrate gene expression analyzed by RNA-Seq, solution sulfur speciation, electron microscopy and spectroscopy. The A. thiooxidansS 2 ⁢ O 3 2 - metabolism model involves the conversion ofS 2 ⁢ O 3 2 - to SO 4 2 - , S0 andS 4 ⁢ O 6 2 - , mediated by the sulfur oxidase complex (Sox), tetrathionate hydrolase (TetH), sulfide quinone reductase (Sqr), and heterodisulfate reductase (Hdr) proteins. These same proteins, with the addition of rhodanese (Rhd), were identified to convert S0 to SO 3 2 - ,S 2 ⁢ O 3 2 - and polythionates in the A. thiooxidans S0 metabolism model. Our combined results shed light onto the important role specifically of TetH inS 2 ⁢ O 3 2 - metabolism. Also, we show that activity of Hdr proteins rather than Sdo are likely associated with S0 oxidation. Finally, our data suggest that formation of intracellularS 2 ⁢ O 3 2 - is a critical step in S0 metabolism, and that recycling of internally generated SO 3 2 - occurs, through comproportionating reactions that result inS 2 ⁢ O 3 2 - . Electron microscopy and spectroscopy confirmed intracellular production and storage of S0 during growth on both S0 andS 2 ⁢ O 3 2 - substrates.

Transcriptome Analysis of Larval Segment Formation and Secondary Loss in the Echiuran Worm Urechis unicinctus

The larval segment formation and secondary loss in echiurans is a special phenomenon, which is considered to be one of the important characteristics in the evolutionary relationship between the Echiura and Annelida. To better understand the molecular mechanism of this phenomenon, we revealed the larval transcriptome profile of the echiuran worm Urechis unicinctus using RNA-Seq technology. Twelve cDNA libraries of U. unicinctus larvae, late-trochophore (LT), early-segmentation larva (ES), segmentation larva (SL), and worm-shaped larva (WL) were constructed. Totally 243,381 unigenes were assembled with an average length of 1125 bp and N50 of 1836 bp, and 149,488 unigenes (61.42%) were annotated. We obtained 70,517 differentially expressed genes (DEGs) by pairwise comparison of the larval transcriptome data at different developmental stages and clustered them into 20 gene expression profiles using STEM software. Based on the typical profiles during the larval segment formation and secondary loss, eight signaling pathways were enriched, and five of which, mTOR, PI3K-AKT, TGF-β, MAPK, and Dorso-ventral axis formation signaling pathway, were proposed for the first time to be involved in the segment formation. Furthermore, we identified 119 unigenes related to the segment formation of annelids, arthropods, and chordates, in which 101 genes were identified in Drosophila and annelids. The function of most segment polarity gene homologs (hedgehog, wingless, engrailed, etc.) was conserved in echiurans, annelids, and arthropods based on their expression profiles, while the gap and pair-rule gene homologs were not. Finally, we verified that strong positive signals of Hedgehog were indeed located on the boundary of larval segments using immunofluorescence. Data in this study provide molecular evidence for the understanding of larval segment development in echiurans and may serve as a blueprint for segmented ancestors in future research.

The response of sulfur dioxygenase to sulfide in the body wall of Urechis unincinctus

Background

In some sedimentary environments, such as coastal intertidal and subtidal mudflats, sulfide levels can reach millimolar concentrations (2–5 mM) and can be toxic to marine species. Interestingly, some organisms have evolved biochemical strategies to overcome and tolerate high sulfide conditions, such as the echiuran worm, Urechis unicinctus. Mitochondrial sulfide oxidation is important for detoxification, in which sulfur dioxygenase (SDO) plays an indispensable role. Meanwhile, the body wall of the surface of the worm is in direct contact with sulfide. In our study, we chose the body wall to explore the SDO response to sulfide.

Methods

Two sulfide treatment groups (50 µM and 150 µM) and a control group (natural seawater) were used. The worms, U. unicinctus, were collected from the intertidal flat of Yantai, China, and temporarily reared in aerated seawater for three days without feeding. Finally, sixty worms with similar length and mass were evenly assigned to the three groups. The worms were sampled at 0, 6, 24, 48 and 72 h after initiation of sulfide exposure. The body walls were excised, frozen in liquid nitrogen and stored at −80 °C for RNA and protein extraction. Real-time quantitative RT-PCR, enzyme-linked immunosorbent assay and specific activity detection were used to explore the SDO response to sulfide in the body wall.

Results

The body wall of U. unicinctus consists of a rugal epidermis, connective tissue, outer circular muscle and middle longitudinal muscle. SDO protein is mainly located in the epidermis. When exposed to 50 µM sulfide, SDO mRNA and protein contents almost remained stable, but SDO activity increased significantly after 6 h (P < 0.05). However, in the 150 µM sulfide treatment group, SDO mRNA and protein contents and activity all increased with sulfide exposure time; significant increases all began to occur at 48 h (P < 0.05).

Discussion

All the results indicated that SDO activity can be enhanced by sulfide in two regulation mechanisms: allosteric regulation, for low concentrations, and transcription regulation, which is activated with an increase in sulfide concentration.

Persulfide Dioxygenase From Acidithiobacillus caldus: Variable Roles of Cysteine Residues and Hydrogen Bond Networks of the Active Site

Persulfide dioxygenases (PDOs) are abundant in Bacteria and also crucial for H2S detoxification in mitochondria. One of the two pdo-genes of the acidophilic bacterium Acidithiobacillus caldus was expressed in Escherichia coli. The protein (AcPDO) had 0.77 ± 0.1 Fe/subunit and an average specific sulfite formation activity of 111.5 U/mg protein (Vmax) at 40°C and pH 7.5 with sulfur and GSH following Michaelis-Menten kinetics. KM for GSH and Kcat were 0.5 mM and 181 s-1, respectively. Glutathione persulfide (GSSH) as substrate gave a sigmoidal curve with a Vmax of 122.3 U/mg protein, a Kcat of 198 s-1 and a Hill coefficient of 2.3 ± 0.22 suggesting positive cooperativity. Gel permeation chromatography and non-denaturing gels showed mostly tetramers. The temperature optimum was 40-45°C, the melting point 63 ± 1.3°C in thermal unfolding experiments, whereas low activity was measurable up to 95°C. Site-directed mutagenesis showed that residues located in the predicted GSH/GSSH binding site and in the central hydrogen bond networks including the iron ligands are essential for activity. Among these, the R139A, D141A, and H171A variants were inactive concomitant to a decrease of their melting points by 3-8 K. Other variants were inactivated without significant melting point change. Two out of five cysteines are likewise essential, both of which lie presumably in close proximity at the surface of the protein (C87 and C224). MalPEG labeling experiments suggests that they form a disulfide bridge. The reducing agent Tris(2-carboxyethyl)phosphine was inhibitory besides N-ethylmaleimide and iodoacetamide suggesting an involvement of cysteines and the disulfide in catalysis and/or protein stabilization. Mass spectrometry revealed modification of C87, C137, and C224 by 305 mass units equivalent to GSH after incubation with GSSH and with GSH in case of the C87A and C224A variants. The results of this study suggest that disulfide formation between the two essential surface-exposed cysteines and Cys-S-glutathionylation serve as a protective mechanism against uncontrolled thiol oxidation and the associated loss of enzyme activity.

Genetic diversity and population structure of penis fish (Urechis unicinctus) based on mitochondrial and nuclear gene markers

Urechis unicinctus is distributed only in Bohai Gulf of China and Korean and Japanese coast. The wild populations of this species have sharply declined in China and Japan. We collected 105 samples from six localities of Bohai Gulf and Korea coast, and investigated genetic diversity and population structure with mitochondrial COI, 16S-rRNA and nuclear 28S-rRNA gene fragments. Genetic diversity of U. unicinctus based on COI sequences was still high (Hd: 0.9595, π: 0.0101), however, 28S-rRNA gene sequences showed low level of genetic diversity (Hd: 0.4084, π: 0.0007). Based on COI sequences, F_ST values between populations ranged from -0.00204 to 0.05210, and 99.12% genetic diversity was contributed by different individuals within population. Both phylogenetic trees and median-joining network did not show clear geographic cluster, haplotypes from different populations were mixed. Our results indicated low level of genetic divergence among different localities of U. unicinctus, and this species should be treated as a whole population among China, Japan and Korea coast during species conservation.

Discovery of a new subgroup of sulfur dioxygenases and characterization of sulfur dioxygenases in the sulfur metabolic network of Acidithiobacillus caldus

Acidithiobacillus caldus is a chemolithoautotrophic sulfur-oxidizing bacterium that is widely used for bioleaching processes. Acidithiobacillus spp. are suggested to contain sulfur dioxygenases (SDOs) that facilitate sulfur oxidation. In this study, two putative sdo genes (A5904_0421 and A5904_1112) were detected in the genome of A. caldus MTH-04 by BLASTP searching with the previously identified SDO (A5904_0790). We cloned and expressed these genes, and detected the SDO activity of recombinant protein A5904_0421 by a GSH-dependent in vitro assay. Phylogenetic analysis indicated that A5904_0421and its homologous SDOs, mainly found in autotrophic bacteria, were distantly related to known SDOs and were categorized as a new subgroup of SDOs. The potential functions of genes A5904_0421 (termed sdo1) and A5904_0790 (termed sdo2) were investigated by generating three knockout mutants (Δsdo1, Δsdo2 and Δsdo1&2), two sdo overexpression strains (OE-sdo1 and OE-sdo2) and two sdo complemented strains (Δsdo1/sdo1’ and Δsdo2/sdo2’) of A. caldus MTH-04. Deletion or overexpression of the sdo genes did not obviously affect growth of the bacteria on S⁰, indicating that the SDOs did not play an essential role in the oxidation of extracellular elemental sulfur in A. caldus. The deletion of sdo1 resulted in complete inhibition of growth on tetrathionate, slight inhibition of growth on thiosulfate and increased GSH-dependent sulfur oxidation activity on S⁰. Transcriptional analysis revealed a strong correlation between sdo1 and the tetrathionate intermediate pathway. The deletion of sdo2 promoted bacterial growth on tetrathionate and thiosulfate, and overexpression of sdo2 altered gene expression patterns of sulfide:quinone oxidoreductase and rhodanese. Taken together, the results suggest that sdo1 is essential for the survival of A. caldus when tetrathionate is used as the sole energy resource, and sdo2 may also play a role in sulfur metabolism.

The NF-κB pathway participates in the response to sulfide stress in Urechis unicinctus

The NF-κB pathway is known to be involved in regulating apoptosis, inflammation and immunity in organisms. In this study, we first identified full-length cDNA sequences of two key molecules in the NF-κB pathway, namely, NEMO and p65, and characterized their responses in the hindgut of Urechis unicinctus (Echiura, Urechidae) exposed to sulfide. The full-length of cDNA was 2491 bp for U. unicinctus NEMO (UuNEMO) and 1971 bp for U. unicinctus p65 (Uup65), and both polyclonal antibodies were prepared using UuNEMO or Uup65 expressed prokaryotically with the sequence of their whole open reading frame. Immunoprecipitation and Western blotting showed that the NF-κB pathway was activated in U. unicinctus exposed to sulfide, in which the content of UuNEMO ubiquitination and nuclear Uup65 increased significantly (p < 0.05) in hindgut tissue of U. unicinctus exposed to sulfide. Furthermore, the mRNA level of UuBcl-xL, a downstream anti-apoptosis gene of the NF-κB pathway, increased significantly (p < 0.05) from 48 h to 72 h and the mRNA level of UuBax, a Bcl-xL antagonist gene, decreased significantly (p < 0.05) at 48 h in the hindgut of U. unicinctus exposed to 50 μM sulfide. During the 150 μM sulfide exposure, the level of UuBcl-xL showed no obvious change, whereas the UuBax mRNA level increased significantly (p < 0.05) at 72 h post-exposure to 150 μM sulfide. We suggested that the activated NF-κB pathway up-regulates UuBcl-xL expression, and evokes an anti-apoptotic response to resist sulfide damage at 50 μM in U. unicinctus. Meanwhile, a Bax-mediated pro-apoptotic response occurs when U. unicinctus is exposed to 150 μM sulfide.

Sulfide exposure results in enhanced sqr transcription through upregulating the expression and activation of HSF1 in echiuran worm Urechis unicinctus

Sulfide is a natural, widely distributed, poisonous substance. Sulfide: quinone oxidoreductase (SQR) is responsible for the initial oxidation of sulfide in mitochondria. To study transcriptional regulation of sqr after sulfide exposure, a 2.6-kb sqr upstream sequence from echiuran worm Urechis unicinctus was cloned by genome walking. Bioinformatics analysis showed 3 heat shock elements (HSEs) in proximal promoter region of the sqr upstream sequence. Moreover, an Hsf1 cDNA in U. unicinctus (UuHsf1) was isolated with a full-length sequence of 2334 bp and its polyclonal antibody was prepared using U. unicinctus HSF1 (UuHSF1) expressed prokaryotically with whole sequence of its open reading frame (ORF). In vivo ChIP and in vitro EMSA assays revealed UuHSF1 could interact with the sqr proximal promoter region. Transient transfection and mutation assays indicated that UuHSF1 bound specifically to HSE (-155bp to -143bp) and enhanced the transcription of sqr. Furthermore, sulfide treatment experiments demonstrated that sulfide could increase the expression of HSF1 protein, and induce trimerization of the protein which binds to HSEs and then activate sqr transcription. Quantitative real-time PCR analysis revealed sqr mRNA level increased significantly after U. unicinctus was exposed to sulfide for 6h, which corresponded to content changes of both trimeric HSF1 and HSF1-HSE complex. We concluded that UuHSF1 is a transcription factor of sqr and sulfide could induce sqr transcription by upregulating the expression and activation of HSF1 in U. unicinctus exposed to sulfide.

Transcriptional response to sulfide in the Echiuran Worm Urechis unicinctus by digital gene expression analysis

Background

Urechis unicinctus, an echiuran worm inhabiting the U-shaped burrows in the coastal mud flats, is an important commercial and ecological invertebrate in Northeast Asian countries, which has potential applications in the study of animal evolution, coastal sediment improvement and marine drug development. Furthermore, the worm can tolerate and utilize well-known toxicant-sulfide. However, knowledge is limited on the molecular mechanism of U. unicinctus responding to sulfide due to deficiency of its genetic information.

Methods

In this study, we performed Illumina sequencing to obtain the first Urechis unicinctus transcriptome data. Sequenced reads were assembled and then annotated using blast searches against Nr, Nt, Swiss-Prot, KEGG and COG. The clean tags from four digital gene expression (DGE) libraries were mapped to the U. unicinctus transcriptome. DGE analysis and functional annotation were then performed to reveal its response to sulfide. The expressions of 12 candidate genes were validated using quantitative real-time PCR. The results of qRT-PCR were regressed against the DGE analysis, with a correlation coefficient and p-value reported for each of them.

Results

Here we first present a draft of U. unicinctus transcriptome using the Illumina HiSeq^TM 2000 platform and 52,093 unique sequences were assembled with the average length of 738 bp and N50 of 1131 bp. About 51.6 % of the transcriptome were functionally annotated based on the databases of Nr, Nt, Swiss-Prot, KEGG and COG. Then based on the transcriptome, the digital gene expression analysis was conducted to examine the transcriptional response to sulfide during 6, 24 and 48 h exposure, and finally 1705, 1181 and 1494 tag-mapped genes were identified as differentially expressed genes in the 6-h, 24-h and 48-h libraries, then were further subjected to pathway analyses.

Conclusions

In the DGE database of U. unicinctus, the alterations in certain known sulfide-related pathways indicate similar changes in response to sulfide. For more than 80 % of the identified pathway members, this is the first report on their association with sulfide stress, among which glycolysis pathway and PIDD involving pathways were unique and discussed in details, and were thought to play important roles in the sulfide tolerance of U. unicinctus. All the results are helpful to explain the mechanism of sulfide tolerance and detoxification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2094-z) contains supplementary material, which is available to authorized users.

Biosynthesis of a central intermediate in hydrogen sulfide metabolism by a novel human sulfurtransferase and its yeast ortholog

Human sulfide:quinone oxidoreductase (SQOR) catalyzes the conversion of H2S to thiosulfate, the first step in mammalian H2S metabolism. SQOR's inability to produce the glutathione persulfide (GSS(-)) substrate for sulfur dioxygenase (SDO) suggested that a thiosulfate:glutathione sulfurtransferase (TST) was required to provide the missing link between the SQOR and SDO reactions. Although TST could be purified from yeast, attempts to isolate the mammalian enzyme were not successful. We used bioinformatic approaches to identify genes likely to encode human TST (TSTD1) and its yeast ortholog (RDL1). Recombinant TSTD1 and RDL1 catalyze a predicted thiosulfate-dependent conversion of glutathione to GSS(-). Both enzymes contain a rhodanese homology domain and a single catalytically essential cysteine, which is converted to cysteine persulfide upon reaction with thiosulfate. GSS(-) is a potent inhibitor of TSTD1 and RDL1, as judged by initial rate accelerations and ≥25-fold lower Km values for glutathione observed in the presence of SDO. The combined action of GSS(-) and SDO is likely to regulate the biosynthesis of the reactive metabolite. SDO drives to completion p-toluenethiosulfonate:glutathione sulfurtransferase reactions catalyzed by TSTD1 and RDL1. The thermodynamic coupling of the irreversible SDO and reversible TST reactions provides a model for the physiologically relevant reaction with thiosulfate as the sulfane donor. The discovery of bacterial Rosetta Stone proteins that comprise fusions of SDO and TSTD1 provides phylogenetic evidence of the association of these enzymes. The presence of adjacent bacterial genes encoding SDO-TSTD1 fusion proteins and human-like SQORs suggests these prokaryotes and mammals exhibit strikingly similar pathways for H2S metabolism.

Whole-genome sequencing reveals novel insights into sulfur oxidation in the extremophile Acidithiobacillus thiooxidans

BACKGROUND

Acidithiobacillus thiooxidans (A. thiooxidans), a chemolithoautotrophic extremophile, is widely used in the industrial recovery of copper (bioleaching or biomining). The organism grows and survives by autotrophically utilizing energy derived from the oxidation of elemental sulfur and reduced inorganic sulfur compounds (RISCs). However, the lack of genetic manipulation systems has restricted our exploration of its physiology. With the development of high-throughput sequencing technology, the whole genome sequence analysis of A. thiooxidans has allowed preliminary models to be built for genes/enzymes involved in key energy pathways like sulfur oxidation.

RESULTS

The genome of A. thiooxidans A01 was sequenced and annotated. It contains key sulfur oxidation enzymes involved in the oxidation of elemental sulfur and RISCs, such as sulfur dioxygenase (SDO), sulfide quinone reductase (SQR), thiosulfate:quinone oxidoreductase (TQO), tetrathionate hydrolase (TetH), sulfur oxidizing protein (Sox) system and their associated electron transport components. Also, the sulfur oxygenase reductase (SOR) gene was detected in the draft genome sequence of A. thiooxidans A01, and multiple sequence alignment was performed to explore the function of groups of related protein sequences. In addition, another putative pathway was found in the cytoplasm of A. thiooxidans, which catalyzes sulfite to sulfate as the final product by phosphoadenosine phosphosulfate (PAPS) reductase and adenylylsulfate (APS) kinase. This differs from its closest relative Acidithiobacillus caldus, which is performed by sulfate adenylyltransferase (SAT). Furthermore, real-time quantitative PCR analysis showed that most of sulfur oxidation genes were more strongly expressed in the S0 medium than that in the Na2S2O3 medium at the mid-log phase.

CONCLUSION

Sulfur oxidation model of A. thiooxidans A01 has been constructed based on previous studies from other sulfur oxidizing strains and its genome sequence analyses, providing insights into our understanding of its physiology and further analysis of potential functions of key sulfur oxidation genes.

Identification and characterization of an ETHE1-like sulfur dioxygenase in extremely acidophilic Acidithiobacillus spp

Elemental sulfur (S(0)) oxidation in Acidithiobacillus spp. is an important process in metal sulfide bioleaching. However, the gene that encodes the sulfur dioxygenase (SDO) for S(0) oxidation has remained unclarified in Acidithiobacillus spp. By BLASTP with the eukaryotic mitochondrial sulfur dioxygenases (ETHE1s), the putative sdo genes (AFE_0269 and ACAL_0790) were recovered from the genomes of Acidithiobacillus ferrooxidans ATCC 23270 and Acidithiobacillus caldus MTH-04. The purified recombinant proteins of AFE_0269 and ACAL_0790 exhibited remarkable SDO activity at optimal mildly alkaline pH by using the GSH-dependent in vitro assay. Then, a sdo knockout mutant and a sdo overexpression strain of A. ferrooxidans ATCC 23270 were constructed and characterized. By overexpressing sdo in A. ferrooxidans ATCC 23270, a significantly increased transcriptional level of sdo (91-fold) and a 2.5-fold increase in SDO activity were observed when S(0) was used as sole energy source. The sdo knockout mutant of A. ferrooxidans displayed a slightly reduced growth capacity in S(0)-medium compared with the wild type but still maintained high S(0)-oxidizing activity, suggesting that there is at least one other S(0)-oxidizing enzyme besides SDO in A. ferrooxidans ATCC 23270 cells. In addition, no obvious changes in transcriptional levels of selected genes related to sulfur oxidation was observed in response to the sdo overexpression or knockout in A. ferrooxidans when cultivated in S(0)-medium. All the results might suggest that SDO is involved in sulfide detoxification rather than bioenergetic S(0) oxidation in chemolithotrophic bacteria.