“Thiodictyon syntrophicum” sp. nov. strain Cad16T is a photoautotrophic purple sulfur bacterium belonging to the family of Chromatiaceae in the class of Gammaproteobacteria. The type strain Cad16T was isolated from the chemocline of the alpine meromictic Lake Cadagno in Switzerland. Strain Cad16T represents a key species within this sulfur-driven bacterial ecosystem with respect to carbon fixation. The 7.74-Mbp genome of strain Cad16T has been sequenced and annotated. It encodes 6237 predicted protein sequences and 59 RNA sequences. Phylogenetic comparison based on 16S rRNA revealed that Thiodictyon elegans strain DSM 232T the most closely related species. Genes involved in sulfur oxidation, central carbon metabolism and transmembrane transport were found. Noteworthy, clusters of genes encoding the photosynthetic machinery and pigment biosynthesis are found on the 0.48 Mb plasmid pTs485. We provide a detailed insight into the Cad16T genome and analyze it in the context of the microbial ecosystem of Lake Cadagno.
PSB belonging to the family of Chromatiaceae are generally found at the interface of aerobic and sulfidic-anaerobic zones that are exposed to sunlight such as stagnant, hypertrophic water bodies, littoral zones and bacterial mats . The genus Thiodictyon was first described by Winogradsky in 1888  and comprises two type strains, Thiodictyon elegans strain DSM 232T and Thiodictyon bacillosum strain DSM 234T. “Thiodictyon syntrophicum” sp. nov. strain Cad16T is the proposed type strain of the species “Thiodictyon syntrophicum”  within the family of Chromatiaceae of the genus Thiodictyon . Cultures of strain Cad16T were isolated from the chemocline of the alpine meromictic Lake Cadagno (Ticino, Switzerland). This lake is characterized by high influx of sulfate, magnesium and calcium in the euxinic monimolimnion which favors the formation of a steep chemocline at 10 to 14 m depth [5, 6]. Within this zone a dense population (up to 107 cells per ml in summer) of mainly anaerobic phototrophic sulfur bacteria belonging to the PSB genera Chromatium,Lamprocystis,Thiodictyon,Thiocystis, and the GSB Chlorobium  is responsible for up to 40% of the total CO2 fixation measured in Lake Cadagno . Strain Cad16T has been shown to be highly active in CO2 fixation both in situ and in vitro . Furthermore, aggregation of strain Cad16T with SRBof the genus Desulfocapsa has been described . In this publication we describe the first complete genome of strain Cad16T providing details especially on CO2 fixation, sulfur metabolism and on CRISPRs. The sequencing of strain Cad16T is part of a larger sequencing project that includes the key species of the microbial community from the anoxic layers of Lake Cadagno.
Classification and features
Strain Cad16T is Gram-negative, the cells are oval-sphere shaped and 1.4–2.4 μm in diameter, non-motile, vacuolated and contain BChl a. Isolate Cad16T can grow as single cells, as well as in cell aggregates with up to 100 cells contained in EPS layer (Fig. 1). It was isolated from the chemocline of Lake Cadagno in a depth of 10–14 m where it grows in a non-obligate mutualistic association with sulfur-reducing bacteria of the genus Desulfocapsa . Based upon morphology and partial 16S rRNA sequence analysis, the strain Cad16T was classified as a member of the genus Thiodictyon within the family Chromatiaceae before . Figure 2 shows the phylogenetic placement of strain Cad16T (complete 16S rRNA sequence) in a 16S rRNA based maximum likelihood phylogenetic tree. The closest relatives of isolate Cad16T are T. bacillosumDSM 234T and T. elegansDSM 232T with 99% sequence identity (partial 16S rRNA sequences). A comparison of the strain Cad16T core genome with other whole genome sequenced PSB confirmed the phylogenetic placemant (Additional file 1: Figure S1).
Strain Cad16T was anaerobically grown in Pfennigs medium , containing per liter: 0.25 g KH2PO4, 0.34 g NH4Cl, 0.5 g MgSO4·7H2O, 0.25 g CaCl2·2H2O, 0.34 g KCl, 1.5 g NaHCO3, 0.5 ml trace element solution SL10, and 0.02 mg vitamin B12 with 2 mM acetate in 100 mL serum bottles with rubber stoppers. The medium was prepared in a 2 l bottle with a N2/CO2 (80%/20%) gas phase. The medium was then reduced with 0.3 g l− 1 Na2S·9H2O (1.10 mM final concentration) and adjusted to a pH of 7.2. Cultures were incubated at 20–23 °C under photoheterotrophic conditions with 6 h light/dark photoperiods with a 40-W tungsten bulb placed at a distance of 60 cm from the cultures (ca. 10 μE m− 2 s− 1).
Different electron donors and carbon substrates were tested under phototautotrophic conditions by Peduzzi et al. [3, 10]. Photolithoautotrophic growth was observed under anoxic conditions with hydrogen sulfide, thiosulfate and elemental sulfur as electron donors. Thereby, elemental sulfur is stored within the periplasma as intermediate oxidation product (Fig. 1). The carbon sources acetate, butyrate, ethanol, formate, fructose, fumarate, glucose, glycerol, lactate, malate, propanol, propionate, pyruvate and succinate were added at 5 mM concentration, respectively. Strain Cad16T was observed to assimilate only acetate, pyruvate and fructose in the presence of sulfide and bicarbonate. Strain Cad16T was additionally tested for chemolithoautrophic growth with bicarbonate under a headspace atmosphere containing 5% O2, 10% CO2 and 85% N2, in the dark. Growth was observed with 0.02% hydrogen sulfid and 0.07% thiosulfate, or with 0.07% sulfide only, respectively. The pigments responsible for the purple-red color of strain Cad16T were analysed spectrometrically in vivo by Peduzzi et al. . Local absorption maxima at 833 nm, 582 nm and 374 nm gave evidence for the presence of BChl a, and at 528 nm for the carotenoid okenone, respectively .
A further characterization of strain Cad16T can be found in Table 1.
A circular representation of the genome sequence and annotation according to the COG criteria is shown in Fig. 3.
Genome sequencing information
Genome project history
Sampling was done in August 2001 using a Friedinger-type bottle on Lake Cadagno. Subsequent isolation and cultivation of strain Cad16T was done in Pfennig’s medium I . gDNA was isolated in November 2014 and sequencing was performed in January 2015. Raw data was assembled in with the SMRTview assembly platform and annotated using the NCBI Prokaryotic Genome Annotation Pipeline. Completeness of the isolate Cad16T sequence was verified using the 31 single copy genes of the Amphora Net analysis platform .
The genome sequence was deposited in GenBank under BioProject PRJNA354524, with the accession numbers CP020370-CP020372. The key elements of the genome studied are listed in Table 2.
Growth conditions and genomic DNA preparation
Strain Cad16T was anaerobically grown in Pfennigs medium  Cells were collected by centrifugation for 15 min at 10,600 g. DNA was extracted using phenol/chloroform/isoamylalcohol solution (25:24:1, v/v, Sigma, Buchs, Switzerland) following the protocol provided by Pacific Biosciences  in combination with phase lock gels (VWR International). gDNA was purified using AMPure beads (Agencourt, Beckman Coulter Life Sciences, Indianapolis, USA) following the E2612 protocol form New England Biolabs . Purity of the DNA was tested using the Qbit UV/VIS absorption reader (Thermo Fisher Scientific, Rheinach, Switzerland).
Genome sequencing and assembly
The library construction and genome sequencing was done on the Pacific Biosciences RS II platform at the Functional Genomic Center Zurich, Zurich, Switzerland. A 10 kb SMRTbell library was constructed using the DNA Template Prep Kit 1.0 (Pacific Biosciences, Menlo Park, USA). SMRTbell template fragments over 10 kb length were used for creating a SMRT bell-Polymerase Complex with P6-C4 chemistry (Pacific Biosciences) according to the manufacturer instructions.
Four SMRT cells v3.0 (Pacific Biosciences) for PacBio RS II chemistry were used for sequencing. Separate sequencing quality reports for all four cells were created through the SMRT portal software.
The SMRT web portal was used for genome assembly with the RS_HGAP_Assembly.2 pipeline from the SMRT Analysis 2.3 server. The polished assembly consists of 153 scaffolds with a mean coverage of 175× and a N50 value of 6,849,178. Thereof, three scaffolds were distinctly longer (6.85, 0.50 and 0.43 Mb, respectively) and showed a coverage greater than 200×, whereas mean coverage dropped below a value of 50× for the remaining 150 scaffolds.
These three scaffolds showed self-similar ends in dot-plot graphs and could be circularized manually.
The genome was manually corrected for SNPs using MiSeq Illumina 300-bp paired-end reads from previous sequencing (unpublished data, N. Storelli, J.F. Pothier, M. Tonolla).
NCBI Prokaryotic Genome Annotation Pipeline (Annotation Software revision 4.1) NCBI Prokaryotic Genome Annotation Pipeline (Annotation Software revision 4.1) was used for gene calling and gene annotation. To identify CRISPR-Cas sequences the CRISPRfinder server was used . The Pfam-A v29 database was used to predict Pfam domains . Transmembrane domains were predicted with the webserver based TMHMM2 program  and signal peptides were predicted with SignalP 4.1 server .
The complete genome of strain Cad16T comprises one circular chromosome (6,837,296 bp) and two circular plasmids pTs485 (484,824 bp) and pTs417 (416,864 bp) (Table 3). The average GC content for the chromosome, and plasmids pTs485 and pTs417, is 66.28%, 65.59 and 65.97%, respectively. A total of 6601 coding sequences were predicted. Thereof, 6237 were predicted to encode proteins whereas six rRNA, 49 tRNA and four ncRNA sequences were predicted. A putative function is assigned for 3471 (46.57%) protein encoding genes (Table 4). The classification of genes into COGs functional categories is given in Table 5. The replicons pTs485 and pTs417 could be made circular, have their own origin of replication each, but do not contain any RNA or house-keeping genes. Therefore, to our understanding, both pTs485 and pTs417 fulfill the plasmid definition.
Extended insights from the genome sequence
PSB typically transform light energy into chemical energy with the membrane bound type 2 photochemical reaction center. The chromosome of strain Cad16T encodes the core antenna proteins LH1, subunits PufA and PufB (THSYN_31145 and THSYN_31140), and the regulatory protein PufQ (THSYN_31110) upstream to the reaction center genes composed of reaction RC subunits PufL, PufM, and PufC (THSYN_31125–31,135). Additional two copies of subunits LH2 alpha and beta (THSYN_31115 and THSYN_31120), respectively, are encoded further downstream, as well as pairwise in two other clusters (THSYN_30995/31005/31030/31040 and THSYN_31000/3100531010/31035/31045), similar as described for the PSB Allochromatium vinosumDSM 180T . The photosynthetic reaction center H subunit PuhA (THSYN_31405) and PucC (THSYN_31410) are clustered upstream with genes encoding RC-LH1 auxiliary proteins (THSYN_31390–31,400). Furthermore, a homologousHiPIP (THSYN_25970) is found in strain Cad16T. It may function as the main electron donor to the photosynthetic reaction center similar as in A. vinosum .
The absorption spectrum of strain Cad16T shows strong absorption peaks at 374 nm, 582 nm and 833 nm which are characteristic for BChl a . The genes for the complete enzymatic pathway from protoporphyrin to chlorophyllide, and further to BChl a (THSYN_31090–31,105, THSYN_31375, THSYN_31385, THSYN_31415–31,445, THSYN_31555, THSYN_32265–32,270), are clustered on pTs485. BChl a formation is thereby catalyzed by an anaerobic type of the Mg-protoporphyrin IX monomethyl ester oxidative cyclase (ChlE) (THSYN_31385) and a light independent protochlorophyllide reductase complex (ChlLNB) (THSYN_31420–31,430) in strain Cad16T.
Strain Cad16T produces okenone as its sole carotenoid  and Crt proteins involved in carotenoid biosynthesis are found on pTs485. The complete synthesis of this keto-carotenoid is mediated through two novel types of carotenoid ketolases, the C-4/4′ ketolase CruO (THSYN_31065) and the oxygen dependent CruS bifunctional desaturase (THSYN_31070) . The characteristic χ-ring of okenone is introduced through the key enzymes CrtY and CrtU (THSYN_31055 and THSYN_31050) [21, 22].
Remarkably, most of the proteins involved in photosynthesis are encoded on plasmid pTs485, forming a PGC (Fig. 3) . The highly modular character of the pufLM and pufC genes of α, β and γ-proteobacteria has been demonstrated previously [24, 25]. To our knowledge, this is the first description of a PGC being localized on a plasmid in a PSB species. Interestingly, the gene cluster is similarly organized as in the γ-proteobacterium Congregibacter litoralis strain KT71T and as in members from the α-proteobacteria families Rhodobacteraceae and Rhodospirillaceae, respectively.
For the photoautotrophic process of CO2 assimilation in PSB, electrons derived from the oxidation of reduced sulfur compounds, are transferred to electron carriers NAD(P)+ and ferredoxin through light energy. During photolithoautotrophic growth under anaerobic conditions, strain Cad16T uses electrons from the oxidation of sulfide, thiosulfate and elemental sulfur as reducing equivalents . Strain Cad16T can use thiosulfate as an electron source during phototrophic growth . No homologous genes for the thiosulfate oxidizing multi-enzyme complex SoxAX, could be found in the strain Cad16T genome. However, soxB (THSYN_26690) and clustered genes encoding SoxYZ (THSYN_09005–09010) that binds thiosulfate were identified in the genome. Remarkably, this gene combination is found in several genome sequenced Ectothiorhodospiraceae. In contrast to the PSB A. vinosumDSM 180T , no homologous sequence for the tetrathionate-forming thiosulfate dehydrogenase TsdA was found. However, a c4 cytochrome type TsdB homolog (THSYN_17090) was identified. Due to this unusual combination of genes involved in thiosulfate oxidation, further studies are needed to elucidate the thiosulfate oxidation pathways in strain Cad16T.
Initial sulfide and thiosulfate oxidation is immediately followed SGB formation in strain Cad16T (Fig. 1). In strain Cad16T the SGB structure is mediated through envelope SGP homologues to SgpA and SgpB (THSYN_20250 and THSYN_05960) from “Thioflavicoccus mobilis” and Thiocystis violascens, respectively. The sequence of SgpC (THSYN_11025) shows homology to Marichromatium species SgpC/CV3. Predicted signal peptides suggest export of for all three SGP proteins into the periplasm in Cad16T, as proposed for A. vinosumDSM 180T .
Moreover, the genome of strain Cad16T encodes the membrane-bound sulfide: quinone oxidoreductases SqrD (THSYN_04215) and SqrF (THSYN_09305). These are possibly involved in the oxidation of sulfide in the periplasm.
The mode of sulfur transport across the inner membrane is not known for PSBs . Organic persulfides such as glutathione or glutathione amide persulfide are proposed as possible candidates. In a next step, the rhodanese-like protein Rhd transfers the sulfur from the persulfide-carrier to the TusA protein in the cytoplasm. The further oxidation steps from sulfur to sulfite are typically mediated through the reverse acting dsr genes in PSB . The strain Cad16T genes in the dsr cluster (THSYN_22480, THSYN_22490–22,545) are arranged in a highly conserved organization similar to A. vinosumDSM 180T, only missing dsrS that is non-essential for sulfur oxidation . The DsrEFH complex mediates persulfate transfer from TusA onto DsrC. The persulfurated form of DsrC is then substrate for the cytoplasmic reverse-acting dissimilatory sulfite reductase DsrAB that catalyzes the formation of sulfite. Finally, DsrMKJOP complex reduces DsrC .
The genome harbors three additional sulfur relay proteins similar to DsrC (THSYN_09485, THSYN_18820 and THSYN_22565) that could function as TusA homologues. In A. vinosumDSM 180T DsrC is able to bind DNA upstream the dsr cluster .
In strain Cad16T, soeABC (THSYN_16370–16,380) encode the sulfur-iron molybdoprotein complex that further oxidizes sulfite to sulfate on the cytoplasmic site of the membrane . Alternatively, strain Cad16T oxidizes sulfite via APS by APS-reductase AprBA (THSYN_16395 and THSYN_16400) and ATP sulfurylase Sat (THSYN_16390), as in other PSB [33, 34]. Thereby, the membrane-bound QmoABHdrCB-complex  (THSYN_16425–6440) possibly functions as an electron acceptor for the AprAB reductase complex since no aprM homolog was found in the strain Cad16T sequence. For the extra-cytoplasmic export of the final oxidation product sulfate, a SulP sulfate permease (THSYN_14085) homolog to A. vinosumDSM 180T is encoded in the strain Cad16T sequence.
Hydrogen uptake and consumption has been shown to be linked to sulfur metabolism in Thiocapsa roseopersicina BBS [36, 37]. Thereby, electrons from hydrogen oxidation in the periplasm by the hyn-type hydrogenase HydSL could be transferred via the Isp membrane complex to the disulfide bound to DsrC. In A. vinosumDSM 180T, transcription of isp1 and isp2 encoding the Isp hydrogenase subunits is upregulated during growth on sulfide . The Isp complex is composed of two subunits, Isp1 and Isp2, that contain similar catalytic domains as DsrM and DsrK, respectively. Similarly, homologous Isp1 and Isp2 proteins (THSYN_28105 and THSYN_28100) may link sulfur to hydrogen metabolisms in strain Cad16T. In accordance, an increase in the sulfide concentration was observed while SGB were consumed by strain Cad16T during incubation in the dark (unpublished results, F. Danza).
Additionally, other [NiFe]-hydrogenases of the Hox and Hup type (THSYN_22655, THSYN_22660 and THSYN_28115) are found in the sequence that could mediate light-dependent H2 evolution as proposed for T. roseopersicina [39, 40].
The Cad16T genome also harbors cys genes (THSYN_05020–05035) that are probably involved in sulfate assimilation under sulfur-limiting conditions. Furthermore, the genome also encompasses genes encoding the CydDC (THSYN_18930 and THSYN_18935) ATP-driven cysteine transport proteins .
In PSB, CO2 fixation is essentially achieved through the reductive pentose phosphate also known as the CBB cycle. In accordance, the strain Cad16T genome harbors the complete CBB enzymatic pathway. On the chromosome, the dimeric RuBis-CO form II (THSYN_13250) clusters with RuBis-CO activation protein subunits CbbR, CbbQ and CbbO, (THSYN_13245, THSYN_13255 and THSYN_13285). Interestingly, small and large RuBis-CO subunits form I (THSYN_29475 and THSYN_29480) cluster together with carboxysome shell and auxiliary proteins on plasmid pTs417 (THSYN_29485–29,520 and THSYN_29530–29,535). The carboxysome may allow efficient photoassimilation across varying CO2 concentrations as proposed for A. vinosumDSM 180T . Previous studies showed different expression regulation for RuBis-CO type I and type II genes from Cad16T suggesting that only the type II is involved in the process of CO2 fixation . Interestingly, the plasmid pTs485 also harbors a RuBis-CO -like protein form III gene (THSYN_31160) upstream the PGC.
The missing sedoheptulose-1,7-bisphosphatase SBP is possibly bypassed by via the fructose-1,6-bisphosphatase (THSYN_25630). The genes gltA citrate synthase (THSYN_12620), fumA fumarate hydratase (THSYN_24360) and sucCD succinyl-CoA ligase (THSYN_00880 and THSYN_00885) that are essential for the TCA cycle, and isocitrate lyase (THSYN_16275) and malate synthase (THSYN_15655) that are essential for the glyoxylate cycle, respectively, are identified in the strain Cad16T sequence. Recently a proteomic study about the capacity of Cad16T to fix CO2 in the dark suggested the presence of a particular archael DC/HB cycle . However, nofurther genes coding for this DC/HB cycle were found. Also a complete set of genes coding for polyhydroxyalkanoic acid synthase PhaC (THSYN_06910) and poly-(3-hydroxybutyrate) depolymerase PhaE (THSYN_06905) are found in the strain Cad16T genome.
Strain Cad16T additionally encodes genes necessary for glycogen polymerisation. The glucose 1-phosphate adenylyltransferase GlgC (THSYN_00810), the glycogen synthase GlgA (THSYN_11615) and the 1,4-alpha-glucan branching enzyme GlgB (THSYN_00805) allow the synthesis of glycogen.
Interestingely, strain Cad16T also has the potential to produce the storage compound cyanophycin normally found in caynobacteria , since the two subunits of the enzyme cyanophycin synthetase (THSYN_26990 and THSYN_26995) are found.
Togheter, these finding provide genetic evidence for the high carbon fixation potential of strain Cad16T in the dark [8, 44].
Anaerobic Fe(II)-oxidation was described for other Thiodictyon strains [45, 46] and evidence of cryptic in situ iron cycling has been demonstrated recently . In accordance with these findings, we found cbb3 type terminal cytochrome C oxidases (THSYN_06760–08775) possibly involved in Fe(II) driven carbon fixation in strain Cad16T genome.
Strain Cad16T grows chemoautotrophically under microaerobic conditions (5% O2) with sulfide, thiosulfate, or sulfide only , as also observed in other PSB in vitro studies with Lamprocystis purpurea [10, 48], Thiocystis violacea and A. vinosum . In situ, strain Cad16T is possibly exposed to low concentration of oxygen produced by oxygenic microbiota at the mixolimnion-chemocline interface . Accordingly, we observe genes encoding sod-type superoxide dismutases (THSYN_20405 and THSYN_22720), as well as fnr and fur-type transcriptional regulators involved in peroxide stress response. In situ, strain Cad16T is possibly exposed to oxygen produced by oxygenic microbiota at the mixolimnion-chemocline interface .
Furthermore, with the genes encoding NifB (THSYN_03975), NifD (THSYN_08880), NifH (THSYN_08885), NifK (THSYN_08875), NifT (THSYN_08870) NifW, NifZ and NifM (THSYN_10720, THSYN_10725 and THSYN_10730), NifX (THSYN_21435) and NifL (THSYN_24590) strain Cad16T could possibly fix nitrogen. Genes encoding the multisubunit urease UreDEFG (THSYN_03745, 03750, 03760 and 03765) and the urea transporter UrtABCDE (THSYN_07940–07955, 03760, 07975) indicate the possible utilisation of urea.
Transmembrane transport proteins
Several membrane transport genes were found in the strain Cad16T genome, including protein secretion system Type II, genes encoding the TAT pathway and several TRAP transporter genes, as well as genes encoding Ton-Tol type and ABC-type transporter complexes. Additionally, a complete TSS4 pilus machinery is encoded in six clusters dispersed on the strain Cad16T chromosome. Notably, also structural components of TSS6 secretion system are found in two clusters on the chromosome (THSYN_11395–11,410) and on pTs485 (THSYN_32540-THSYN_32580). Two effector proteins of the VrgG family were identified. THSYN_15360 belongs to the vgr_GE type Rhs family proteins similar sequences found in β-proteobacterial family of the Burkholderiaceae whereas THSYN_32425 is conserved in γ-proteobacteria and contains a type IV Rhs element. Togheter, the secretion machinery allows strain Cad16T to interact within the highly populated chemocline with up to 107 bacterial cells per milliliter. The secretion and uptake mechanism may also play a key role in the cell-to-cell contact with Desulfocapsa thiozymogenes.
Buoyancy regulation and chemotaxis
Strain Cad16T can possibly regulate buyoncy by gas vesicles that are formed with the encoded structural gas vesicle proteins. Whereas GvpA proteins forms the vesicle core (THSYN_11790, THSYN_11825, THSYN_15290, THSYN_18705 and THSYN_31215), GvpFL (THSYN_11800 and THSYN_18685), GvpK (THSYN_11785) and GvpN (THSYN_11815 and THSYN_18695) further stabilize the structure. Proteins homologoues to the transcriptional regulatory factors GvrA (THSYN_11850) and GvrC (THSYN_11830) from the enterobacterium Serratia sp. ATCC 39006 are also found in Cad16T.
The diurnal and sesonal behavior of vacuolated Chromatiaceae has been described for different lakes [50, 51]. In strain Cad16T a diguanylate cyclase (THSYN_19835) is found upstream the circadian clock genes kaiCBB (THSYN_19820–19,830). These genes act togheter  and may synchronize optimal flotation within the chemocline.
Bacterial CRISPR-Cas systems provide a mechanism against bacteriophage infection and plasmid transformation . A CRISPR locus is composed out of a 300–500 bp leader sequence, spacer sequences (21–72 bp), complementary to foreign DNA, and direct repeats (DRs, 24–40 bp) flanking them [53,54,55]. Adjacent cas genes encode protein that are co-transcribed with the CRISPR locus and interfere with invading DNA guided by the specific spacers [56, 57].
Five CRISPR repeat regions (CRR1-CRR5) were identified in the genome of strain Cad16T, four being located on the chromosome and one on the plasmid pTs485 (Fig. 3). The number of DRs ranges from 19 (CRR4) to 146 (CRR2) as seen in Table 6.
Furthermore, three CRISPR-Cas loci were identified in the strain Cad16T sequence, containing cas3 genes that are characteristic for type I CRIPSR-Cas systems . A complete CRISPR-Cas loci (THSYN_08045–08070) is located 201 bp upstream of CRR2 and is assigned to subtype I-U, containing the signature protein (THSYN_08055) of the GSU0054 family (TIGR02165 and a cas3, THSYN_08070) with C-terminal HD domain (TIGR01596) . Another CRISPR array (THSYN_19240–19,290) is located 182 bp upstream of CRR3 and is classified as subtype I-C due to the cas8c gene and the lack of a cas6 sequence. Additionally, an incomplete CRISPR-Cas locus (CRR5) is identified on plasmid pTs485, encoding for Cas2, Cas1, (THSYN_19240–19,245, THSYN_19265, THSYN_19275, THSYN_19285 and THSYN_19,290).
We report on the first complete genome sequence of “Thiodictyon syntrophicum” sp. nov. strain Cad16T and the metabolic versatility of this environmentally relevant organism. The observed carbon fixation potential can be explained by the highly developed photosynthesis machinery that is coupled to the sulfur and carbon metabolism. Within the changing conditions in the chemocline, strain Cad16T is able to optimally use light, different organic and inorganic carbon compounds, reduced sulfur, nitrogen and oxygen. The two 0.4 Mb plasmids found in Cad16T are unique for known PSB species and we report structural similarity to sequences from α- and γ-proteobacterial phototrophs. The availability of the complete genome sequence of strain Cad16T will facilitate further studies that elucidate its role as key species of the chemocline and the tight association with the Desulfocapsa sp. and the interaction with different PSB and GSB species present in the anoxic part of Lake Cadagno. Due to the limited molecular data on other Thiodictyon strains and no reference strains available, no (digital) DNA-DNA hybridization experiments could be performed. However, the result from phylogenetic analyses on 16S rRNA sequence level, comparative genomic analyses as well a morphological and physiological differences (see above) indicate a novel species within the genus Thiodictyon.
“Thiodictyon syntrophicum” (syn.tro’phi.cum. Gr. pref. Syn, together with; Gr. adj. Trophikos, nursing, tending or feeding; N.L. neut. Adj. syntrophicum, syntrophic).
Gram-negative, cells are oval-round shaped and 1.4–2.4 μm in diameter, non-motile, vacuolated and contain BChl a and okeneone. Growth as single cells, as well as in aggregates with up to 100 cells in a EPS layer. Assimilation of elemental sulfur in intracellular sulfur globules. Grow photoautotrophically in Pfennig's minimal medium with a doubling time of 121 h at 20–23 °C, a pH of 6.8–7.2, at 1 mM sulfide and a photoperiode of 12 h dark/ 12 h light. Dense cultures show a milky purple-red and milky color. Carbon assimilation via Calvin cycle. Following carbon substrates were utilized at a concentration of 5 mM: acetate, fructose and pyruvate. No growth was observed with 5 mM butyrate, ethanol, formate, fumarate, glucose, glycerol, lactate, malate, propanol, propionate and succinate, respectively. Chemolitoautotrophic growth was observed with 5% Oxygen and 0.02% hydrogen sulfid and 0.07% thiosulfate, or with 0.07% sulfide only, respectively.
The type strain Cad16T (=JCM 15483T =KCTC5955T) was isolated from a sulfidic chemocline in the alpine Lake Cadagno in Switzerland. The genome size of the type strain is 6.84 Mb (chromosome), contains two plasmids, pTs485 (0.49 Mb) and pTs417 (0.42 Mb) and the G + C content of the genome is 66.22%. The 16S RNA gene sequence of strain Cad16T is deposited under the GenBank/EMBL/DDBJ accession number AJ511274. The complete genome sequence of the type strain Cad16T is deposited under the GenBank ID CP020370, CP020371 and CP020372. The type strain has been deposited both at the Japan Collection of Microorganisms (JCM 15483T) and at the Korean Collection for Type Cultures (KCTC 5955T).
clustered regularly interspaced short palindromic repeats
Peduzzi S, Storelli N, Welsh A, Peduzzi R, Hahn D, Perret X, et al. Candidatus “Thiodictyon syntrophicum”, sp. nov., a new purple sulfur bacterium isolated from the chemocline of Lake Cadagno forming aggregates and specific associations with Desulfocapsa sp. Syst Appl Microbiol. 2012;35:139–44.
Imhoff JF. Chromatiales Ord. Nov. In: Brenner DJ, Krieg NR, Staley JT, Garrity GM, Boone DR, De Vos P, et al., editors. Bergey’s manual® Syst. Bacteriol. Vol. two Proteobacteria part B Gammaproteobacteria. Boston: Springer US; 2005. p. 1–59. Available from: http://dx.doi.org/10.1007/0-387-28022-7_1
Del Don C, Hanselmann KW, Peduzzi R, Bachofen R. The meromictic alpine Lake Cadagno: orographical and biogeochemical description. Aquat Sci. 2001;63:70–90.
Tonolla M, Storelli N, Danza F, Ravasi D, Peduzzi S, Posth NR, et al. Lake Cadagno: Microbial Life in Crenogenic Meromixis. In: Gulati RD, Zadereev ES, Degermendzhi AG, editors. Ecol Meromictic Lakes. Cham: Springer International Publishing; 2017 [cited 2017 Jun 16]. p. 155–86. Available from: http://link.springer.com/10.1007/978-3-319-49143-1_7
Tonolla M, Peduzzi R, Hahn D. Long-term population dynamics of phototrophic sulfur Bacteria in the chemocline of Lake Cadagno, Switzerland. Appl Environ Microbiol. 2005;71:3544–50.
Camacho A, Erez J, Chicote A, Florín M, Squires MM, Lehmann C, et al. Microbial microstratification, inorganic carbon photoassimilation and dark carbon fixation at the chemocline of the meromictic Lake Cadagno (Switzerland) and its relevance to the food web. Aquat Sci. 2001;63:91–106.
Storelli N, Peduzzi S, Saad MM, Frigaard N-U, Perret X, Tonolla M. CO2 assimilation in the chemocline of Lake Cadagno is dominated by a few types of phototrophic purple sulfur bacteria. FEMS Microbiol Ecol. 2013;84:421–32.
Peduzzi S, Tonolla M, Hahn D. Isolation and characterization of aggregate-forming sulfate-reducing and purple sulfur bacteria from the chemocline of meromictic Lake Cadagno, Switzerland. FEMS Microbiol Ecol. 2003;45:29–37.
AmphoraNet: The webserver implementation of the AMPHORA2 metagenomic workflow suite [Internet]. [cited 2016 Jun 5]. Available from: http://www.sciencedirect.com/science/article/pii/S0378111913014091
SharedProtocol-Extracting-DNA-usinig-Phenol-Chloroform.pdf [Internet]. [cited 2017 Jan 5]. Available from: http://www.pacb.com/wp-content/uploads/2015/09/SharedProtocol-Extracting-DNA-usinig-Phenol-Chloroform.pdf
Agencourt AMPure XP Bead Clean-up - NEBNext Microbiome DNA Enrichment Kit (E2612) | NEB [Internet]. [cited 2017 Jan 5]. Available from: https://www.neb.com/protocols/2013/04/18/agencourt-ampure-xp-bead-clean-up-e2612
Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.
Verméglio A, Li J, Schoepp-Cothenet B, Pratt N, Knaff DB. The role of high-potential Iron protein and cytochrome c8 as alternative Electron donors to the reaction Center of Chromatium vinosum†. Biochemistry (Mosc). 2002;41:8868–75.
Bauer CE, Buggy JJ, Yang ZM, Marrs BL. The superoperonal organization of genes for pigment biosynthesis and reaction center proteins is a conserved feature in Rhodobacter capsulatus: analysis of overlapping bchB and puhA transcripts. Mol Gen Genet. 1991;228:433–44.
Nagashima KVP, Verméglio A, Fusada N, Nagashima S, Shimada K, Inoue K. Exchange and complementation of genes coding for photosynthetic reaction center Core subunits among purple Bacteria. J Mol Evol. 2014;79:52–62.
Frigaard N-U, Dahl C. Sulfur Metabolism in Phototrophic Sulfur Bacteria. In: Poole RK, editor. Adv. Microb. Physiol. [Internet]. Academic Press; 2008 [cited 2016 Mar 16]. p. 103–200. Available from: http://www.sciencedirect.com/science/article/pii/S0065291108000027
Pott AS, Dahl C. Sirohaem sulfite reductase and other proteins encoded by genes at the dsr locus of Chromatium vinosum are involved in the oxidation of intracellular sulfur. Microbiology. 1998;144:1881–94.
Grein F, Pereira IAC, Dahl C. Biochemical characterization of individual components of the Allochromatium vinosum DsrMKJOP transmembrane complex aids understanding of complex function in vivo. J Bacteriol. 2010;192:6369–77.
Dahl C, Franz B, Hensen D, Kesselheim A, Zigann R. Sulfite oxidation in the purple sulfur bacterium Allochromatium vinosum: identification of SoeABC as a major player and relevance of SoxYZ in the process. Microbiology. 2013;159:2626–38.
Hipp WM, Pott AS, Thum-Schmitz N, Faath I, Dahl C, Trüper HG. Towards the phylogeny of APS reductases and sirohaem sulfite reductases in sulfate-reducing and sulfur-oxidizing prokaryotes. Microbiology. 1997;143:2891–902.
Tengölics R, Mészáros L, Győri E, Doffkay Z, Kovács KL, Rákhely G. Connection between the membrane electron transport system and Hyn hydrogenase in the purple sulfur bacterium, Thiocapsa roseopersicina BBS. Biochim Biophys Acta BBA - Bioenerg. 1837;2014:1691–8.
Laurinavichene T. The effect of sulfur compounds on H2 evolution/consumption reactions, mediated by various hydrogenases, in the purple sulfur bacterium, Thiocapsa roseopersicina - Springer [Internet]. [cited 2016 Jun 20]. Available from: http://link.springer.com/article/10.1007%2Fs00203-007-0260-7
Weissgerber T, Sylvester M, Kröninger L, Dahl C. A comparative quantitative proteomic study identifies new proteins relevant for sulfur oxidation in the purple sulfur bacterium Allochromatium vinosum. Appl Environ Microbiol. 2014;80:2279–92.
Pittman MS, Corker H, Wu G, Binet MB, Moir AJG, Poole RK. Cysteine is exported from the Escherichia coli cytoplasm by CydDC, an ATP-binding cassette-type transporter required for cytochrome assembly. J Biol Chem. 2002;277:49841–9.
Storelli N, Saad MM, Frigaard N-U, Perret X, Tonolla M. Proteomic analysis of the purple sulfur bacterium Candidatus “Thiodictyon syntrophicum” strain Cad16T isolated from Lake Cadagno. EuPA Open Proteomics [Internet]. [cited 2013 Nov 28]; Available from: http://www.sciencedirect.com/science/article/pii/S2212968513000172
Ehrenreich A, Widdel F. Anaerobic oxidation of ferrous iron by purple bacteria, a new type of phototrophic metabolism. Appl Environ Microbiol. 1994;60:4517–26.
Berg J. Intensive cryptic microbial iron cycling in the low iron water column of the meromictic Lake Cadagno - Berg - 2016 - Environmental Microbiology - Wiley Online Library [Internet]. [cited 2016 Oct 23]. Available from: http://onlinelibrary.wiley.com/doi/10.1111/1462-2920.13587/full
Overmann J, Pfennig N. Continuous chemotrophic growth and respiration of Chromatiaceae species at low oxygen concentrations. Arch Microbiol. 1992;158:59–67.
Egli K, Wiggli M, Fritz M, Klug J, Gerss J, Bachofen R. Spatial and temporal dynamics of a plume of phototrophic microorganisms in a meromictic alpine lake using turbidity as a measure of cell density. Aquat Microb Ecol. 2004;35:105–13.
Al-Attar S, Westra ER, van der Oost J, SJJ B. Clustered regularly interspaced short palindromic repeats (CRISPRs): the hallmark of an ingenious antiviral defense mechanism in prokaryotes. Biol Chem. 2011;392:277–89.
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, Ashburner M, Axelrod N, Baldauf S, Ballard S, Boore J, Cochrane G, Cole J, Dawyndt P, De Vos P, dePamphilis C, Edwards R, Faruque N, Feldman R, Gilbert J, Gilna P, Glöckner FO, Goldstein P, Guralnick R, Haft D, Hancock D, Hermjakob H, Hertz-Fowler C, Hugenholtz P, Joint I, Kagan L, Kane M, Kennedy J, Kowalchuk G, Kottmann R, Kolker E, Kravitz S, Kyrpides N, Leebens-Mack J, Lewis JE, Li K, Lister AL, Lord P, Maltsev N, Markowitz V, Martiny J, Methe B, Mizrachi I, Moxon R, Nelson K, Parkhill J, Proctor L, White O, Sansone S-A, Spiers A, Stevens R, Swift P, Taylor C, Tateno Y, Tett A, Turner S, Ussery D, Vaughan B, Ward N, Whetzel T, Gil IS, Wilson G, Wipat A. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26(5):541–547.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–1797.
We kindly thank Nicole Liechti, Corinne Oechslin, Pierre H.H. Schneeberger and Christian Beuret. We also thank the assistance of Andrea Patrignani from the Functional Genomics Center Zurich for performing the PacBio RS II sequencing.
Authors and Affiliations
University of Geneva, Sciences III, Department of Botany and Plant Biology, Microbiology Unit, 1211, Geneva, Switzerland
Samuel M. Luedin, Francesco Danza, Nicola Storelli & Mauro Tonolla
University of Applied Sciences of Southern Switzerland (SUPSI), Department of Environment, Constructions and Design (DACD), Laboratory of Applied Microbiology (LMA), Via Mirasole 22A, 6500, Bellinzona, Switzerland
Samuel M. Luedin, Francesco Danza, Nicola Storelli & Mauro Tonolla
Federal Office for Civil Protection, Spiez Laboratory, Biology Division, Spiez, Switzerland
Samuel M. Luedin & Matthias Wittwer
Zurich University of Applied Sciences (ZHAW), Institute of Natural Resource Sciences, Environmental Genomics and System Biology Research Group, Wädenswil, Switzerland
Joël F. Pothier
University of Copenhagen, Department of Biology, Helsingør, Denmark
SML and FD drafted the manuscript, carried out cultivation, DNA extraction and purification and microscopy. SML, JFP and MW participated in the genome assembly, correction and annotation. SML, JFP, FD, NS, MW, NUF and MT discussed and analyzed the data and revised the manuscript. MT, JFP and MW conceived of and supervised the study. All authors read and approved the final manuscript.
Figure S1. Phylogenetic placement of “T. syntrophicum” strain Cad16T within the other 12 Chromatiaceae species with a publicly available whole genome sequences. Additionally, the closely related phylogenetic lineages Nitrosococcus, Rheinheimera and Arsukibacterium are also included. Strain Cad16T is most closely related to L. purpurea DSM 4197. The maximum likelihood tree was inferred from 100 concatenated single-copy orthologues sequences  and a total of 1000 bootstrap replicates were performed. Numbers at the nodes indicate the SH-aLRT support (%) and ultrafast bootstrap support (%). OrthoMCL , was used to define at set orthologues proteins between these 23 species. Hundred single-copy orthologues were randomly chosen and aligned with MUSCLE  . The best-fit phylogenetic model and subsequent consensus tree computation, based on maximum-likelihood and 1000 bootstrap iterations, was performed with the IQ-TREE software . Nodes with both, 100% SH-aLRT and ultrafast bootstrap support, are indicated with filled black circle symbols for convenience. (TIF 57220 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Luedin, S.M., Pothier, J.F., Danza, F. et al. Complete genome sequence of “Thiodictyon syntrophicum” sp. nov. strain Cad16T, a photolithoautotrophic purple sulfur bacterium isolated from the alpine meromictic Lake Cadagno.
Stand in Genomic Sci13, 14 (2018). https://doi.org/10.1186/s40793-018-0317-z