2.1 Microorganism enrichment, isolation and idenfication
The biogas slurry sample was collected from Shanghai Liming Resources Recycling Co., LTD, Shanghai city, PR China. Firstly, 1 mL slurry sample was incubated in 50 mL liquid medium on a rotary shaker with 200 rpm at 30℃. The liquid enrichment medium contained 50.0 mmol/L Na2S2O3, 1.0 g/L K2HPO4, 0.25 g/L KH2PO4, 0.2 g/L MgCl2, 2.0 g/L NH4Cl, and 10.0 g/L glucose. The pH was adjusted to 6.8. Every 48 h, 1 mL microbial suspension was subcultured into fresh enrichment medium and the conditions were maintained consistent. After 4 cycles of enrichment, serial dilution was adopted to isolate the strain on solid medium (liquid enrichment medium supplemented with 1% agar powder). The isolation process was repeated 3 times to ensure the purity of isolates and the single colony was inoculated into fresh liquid medium and the consumption of sodium thiosulfate was tested. The one isolate with the best growth and thiosulfate removal performance was selected and used for further investigation.
The genomic DNA of selected isolate was extracted with AxyPrep™ Bacterial Genomic DNA Miniprep Kit (AP-MN-BT-GDNA-50, AXYGEN) following the manufacture’s instruction. The 16S rDNA sequence was amplified by polymerase chain reaction (PCR) with KOD-plus-neo enzyme (KOD-401, TOYOBO) and primer pair of 27F (5’-AGAGTTTTGATCCTGGCTCAG-3’) and 1492R (5’-TACGGTTACCTTGTTACGACTT-3’). The PCR product was purified and sequenced by Sangon Biotech Co., Ltd., (Shanghai, China). The sequence acquired was compared with that of other microorganisms in GenBank by BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
2.2 Optimization of cultivation conditions
To improve the biodesulfurization ability, the effects of cultivation conditions (temperature and shaking speed) were investigated. The seed culture was prepared with LB medium (tryptone 10.0 g/L, yeast extract 5.0 g/L, NaCl 10.0 g/L) under standard conditions (pH=6.8, temperature=30 ℃ and shaking speed=200 rpm), and inoculated into growth medium which was the same as liquid enrichment medium. Different temperatures (20 ℃, 25 ℃, 30 ℃, 35 ℃ and 42 ℃) and shaking speeds (100 rpm, 150 rpm, 200 rpm and 250 rpm) were applied, while other parameters were consistent with standard conditions. The cultivation process lasted for 28 h. Samples were taken at different time points during culture and centrifuged at 12000 g for 10 min. The cell mass was resuspended in sterilized water and the optical density at 600 nm was measured by a UV-Vis spectrometer (DR2800, HACH, USA) to monitor cell growth. The supernatant was collected to analyze the consumption of thiosulfate, generation of tetrationate and sulfane sulfur. All treatments were performed in triplicate.
2.3 Analytical techniques
The concentration of thiosulfate was measured by using high performance liquid chromatography (HPLC) system (SCL-10Avp, Shimadzu, Japan) equipped with an electrical conductivity detector (CDD-10Avp, Shimadzu, Japan) and an anion exchange column (IC-A3, Shimadzu, Japan) (Mu et al. 2021). The mixture solution of 8 mmol/L p-Hydroxybenzoic acid, 3.2 mmol/L Bis-Tris and 50 mmol/L boric acid was used as the mobile phase, running at a flow rate of 1.0 mL/min. The column temperature was maintained at 40 ℃.
The concentration of tetrathionate was also measured by the same HPLC system and IC-A3 column, while with an UV detector (CDD-10Avp, Shimadzu, Japan) (Bak et al. 1993). The compounds were eluted with a elution of 40% mobile phases A and 60% mobile phases B at a flow rate of 1.0 mL/min. Solvent A consisted of deionized water and 50 mmol/L NaCl (ACS regent, purity ≥99%), and solvent B consisted of methanol (HPLC grade) and acetonitrile (HPLC grade) in a ratio of 1: 5. The column temperature was also maintained at 40 ℃ and the signal was monitored at 230nm.
The concentration of sulfane sulfur was determined by fluorescence reaction with SSP4 (PB10, DOJINDO) (Ikeda et al. 2017)(Olson et al. 2018). The total volume of reaction system was 100 µL, containing 80 µL DMSO, 10 µL 0.1mmol/L SSP4 solution, and 10 µL sample. After incubation at room temperature for 10 min, the fluorescence intensity was measured at 515 nm (λex=482nm) by a fluorophotometer (Synergy H1, BioTek, US). Potassium polysulfide (12665, Sigma) was used to construct the standard curve and the sulfane sulfur concentration was expressed as potassium polysulfide equivalent.
2.4 Genome sequencing
The whole-genome sequencing was performed by Grandomics Biosciences Co., Ltd, Wuhan, PR China. The genomic DNA was prepared by the CTAB method and followed by purification with QIAGEN Genomic kit (13343, QIAGEN). The NEBNext Ultra II End Repair/dA-tailing Kit (E7546, NEB) and adapter in the ligation library construction kit (SQK-LSK109, Oxford Nanopore Technologies) were used for constructing sequencing library, and then the libraries were launched on a Nanopore PromethION sequencer instrument (Oxford Nanopore Technologies, UK). The sequenced data were assembled and aligned by NextDenovo (https://github.com/Nextomics/NextDenovo) / NextPolish (https://github.com/Nextomics/NextPolish). Genes were predicted from the assembled scaffolds using the Prodigal software with default parameters. The Gene ontology (GO) annotations of predicted genes were determined by the Blast2GO program, and the corresponding function annotation was completed by blasting genes against Kyoto Encyclopedia of Genes and Genomes (KEGG) databases (http://www.genome.jp/kegg/). The genome has been deposited in the NCBI database with the accession number SAMN21246304.
2.5 Analysis of gene expression level
For gene expression analysis, samples were taken every 2 h starting from 4 h. Cells were collected by centrifugation at 12,000 g for 10 min at 4 ℃. Total RNA was extracted from collected cells with TaKaRa MiniBEST Universal RNA Extraction Kit (RR037A, TaKaRa) according to the manufacturer’s instruction. The concentrations of RNA samples were determined with the spectrophotometer (NanoDrop 2000C, ThermoFisher Scientific, USA). Only samples with the 260/280 value between 1.9~2.1 and 260/230 value greater than 2.0 were used for further analysis. The integrity of RNA samples was assessed by agarose gel electrophoresis. The cDNA synthesis was performed at 37℃ for 15 min using PrimeScript™ RT reagent Kit (RR047A, TaKaRa) following the manufacturer’s instructions. Prior to reverse transcription, the possible DNA contamination was removed by gDNA Eraser at 42℃ for 2 min.
Based on the whole genome sequencing and annotation, genes related to sulfur metabolism were selected including glpE, tsdA, TST1, TST2, cysJ1, cysJ2, cysI, SQR, cysM, cysK, SOD1, SOD2, iscS1, and iscS2. The expression level of those related genes were analyzed by RT-qPCR. The primers sequences and annealing temperature of the target genes and 16S rRNA as reference gene were listed in Table 1. The reliability of each primers pair was checked using genomic DNA as template. The 20 µL RT-qPCR reaction mixtures contained 10 µL TB Green® Premix Ex Taq™ (Tli RNaseH Plus) (RR420A, TaKaRa), 0.4 µL ROX Reference Dye II, 0.4 µL each of 10 µM forward and reverse primers, 1.0 µL cDNA template (the amount of template in mixture is less than 100 ng), and 7.8 µL ddH2O. Controls with no template were also performed for each primer pair. The RT-qPCR amplified reactions were conducted using a LightCycler 96 real-time fluorescent quantitative PCR instrument (QuantStudio™ 7 Flex, ThermoFisher Scientific, USA) with the following program: denaturing at 95 ℃ for 30 s (holding stage), followed by 40 cycles (95 ℃ for 5 s and 60 ℃ for 34 s; PCR stage) and a final extension step (95 ℃ for 15 s and 60 ℃ for 1 min; melting curve stage). Each gene was analyzed in triplicate and the relative expression level calculated by the 2−ΔΔCt method. The values of standard conditions (no thiosulfate addition) at 4 h were used as baseline.
Table 1
Primers and annealing temperature of the target genes
Genes | Primers | Sequences (5’-3’) | Annealing temperature (℃) |
glpE | forward | CGTGAGCAAGGTGCGGTGGTCGT | 60.0 |
reverse | CTCTGGCTGGAGTTGCCGTGGTA |
TST1 | forward | GGGCGCTTCATCTGGTTGCTCG | 55.0 |
reverse | CCAGTCGGCTTTGCAGGTATTCG |
TST2 | forward | GGACATGAAGCTGCTGATCGACGCC | 55.0 |
reverse | AGTTGGTCCGCAGGCAGGAAACG |
cysJ1 | forward | TGGGTTATGCGGATGGGTTCG | 55.0 |
reverse | AACGGGCTGCGGTCATTGTCG |
cysJ2 | forward | GTTCTACCTGCTGTTCGCCCTCACCG | 55.0 |
reverse | TTGGCCCAGATCGCGTCGTAGTCG |
cysI | forward | CCGGTGAGCTGAGCGAAGAAGAAT | 60.0 |
reverse | TGGCGGGTACTGATGTGGGCGTA |
sqr | forward | CCAGCCGCCGATGCCAATCAAGT | 55.0 |
reverse | CGCTACCAGGCGATGGGAATAGTTGA |
cysM | forward | GCCCGAGCAATACCTGCCGAAGA | 60.0 |
reverse | AGCATGGCTGCAACCGCACCAC |
cysK | forward | CCAAGATCGAAGGGCGTAACCCA | 60.0 |
reverse | GCATGGTCAGGAGCAGCTTGTAG |
SOD1 | forward | TGCACCAGCGCCCCTCGTGCGA | 55.0 |
reverse | CCTTGAAGTCCTCTGTCTTGAGC |
SOD2 | forward | CCTGGAGTATCACCACGACAAGC | 60.0 |
reverse | TGTTGAAGATGCCACCCGAAGA |
iscS1 | forward | GCAGCCGGGCCACCTGATCACCA | 60.0 |
reverse | TTGACCGCCATCAGCGACACCAG |
iscS2 | forward | AGACCAAGGGCAAGCACATCATCAC | 60.0 |
reverse | TCGAGCCGACTTCGTTGTTCACG |
tsdA | forward | CGGTGGGCAATGGCATGAACT | 55.0 |
reverse | TGGATACGCTCGGCGAAGGTG |
16S rRNA (Liang et al. 2021) | forward | CCTACGGGAGGCAGCAG | 58.0 |
reverse | ATTACCGCGGCTGCTGG |