- 1Department of Solar Materials, Helmholtz-Centre for Environmental Research (UFZ), Leipzig, Germany
- 2Department of Environmental Microbiology, Helmholtz-Centre for Environmental Research (UFZ), Leipzig, Germany
Cyanobacteria are highly promising microorganisms in forthcoming biotechnologies. Besides the systematic development of molecular tools for genetic engineering, the design of chassis strains and novel reactor concepts are in focus. The latter includes capillary biofilm reactors (CBR), which offer a high surface area-to-volume ratio and very high cell densities. In this context, Tolypothrix sp. PCC 7712 was found to be highly suited for this reactor system due to maximal surface coverage, extraordinarily strong biofilm attachment, and high biomass formation. Here, we provide the genome sequence of Tolypothrix sp. PCC 7712 to potentially allow targeted strain engineering. Surprisingly, it was almost identical to an available incomplete genome draft of Tolypothrix sp. PCC 7601. Thus, we completely sequenced this strain as well and compared it in detail to strain PCC 7712. Comparative genome analysis revealed 257 and 80 unique protein-coding sequences for strains PCC 7601 and PCC 7712, respectively. Clustering genomes based on average nucleotide identity (ANI) and 16S rRNA homology showed 99.98% similarity and only minor distance, respectively, between the two strains in contrast to 21 other cyanobacterial genomes. Despite these high similarities, both strains differ in the ability to fix atmospheric nitrogen and show specific sequence variations, which are discussed in the paper.
Introduction
Cyanobacteria are the only prokaryotes performing oxygenic photosynthesis; i.e., oxygen is released as a side product of light-driven water oxidation. The obtained electrons are used to drive an autotrophic metabolism based on CO2 fixation. Furthermore, multiple species are capable of fixing dinitrogen gas (N2; Tsygankov, 2007; Bharti et al., 2017). These features form the basis for a sustainable biotech-workhorse being independent of organic carbon and reduced nitrogen compounds, which usually add significantly to the ecological footprint of biotech processes. Although numerous proof-of-concept studies show the feasibility of using cyanobacteria as solar cell factories for producing commodity products (Angermayr et al., 2015; Betterle and Melis, 2019; Hoschek et al., 2019; Xie and Lindblad, 2022), only a few examples exist where cyanobacterial biocatalysts are applied at an economic scale. All these processes are based on biomass, which needs to be harvested to extract, e.g., pigments or lipids (Jones and Mayfield, 2012; Garlapati et al., 2019; Nagappan et al., 2020). Persisting challenges for applying cyanobacteria in production processes are low productivity and reaction stability, low cell densities due to light limitation, and insufficient light energy for efficient product conversion (Posten, 2009).
Recently, a novel reactor concept has been introduced, enabling long-term and high cell-density cultivation of cyanobacteria and potentially allowing for continuous production processes. In this regard, the unicellular model strain Synechocystis sp. PCC 6803 was grown as a biofilm in a capillary biofilm reactor (CBR; Hoschek et al., 2019). Biofilms are surface-attached microbial communities, supported and protected by a self-produced extracellular matrix containing mainly polysaccharides and other biopolymers like proteins, DNA, or lipids. They are widespread in nature with cyanobacteria playing a key role as primary producers in complex biofilms also termed microbial mats. In a biotechnological context, biofilms can be regarded as a robust biocatalyst, naturally immobilized to a given surface, enabling continuous bioprocessing (Halan et al., 2012).
In addition to Synechocystis, other cyanobacteria were screened for their utilization in CBRs (Bozan et al., 2022). In this survey, the filamentous, diazotrophic strain Tolypothrix sp. PCC 7712 (also known as Gloeotrichia sp.), first isolated from a soil sample collected in New York, United States, was identified as a top-performing organism. It was superior to all other strains investigated, e.g. in biofilm biomass formation, a low biofilm detachment rate, and high surface coverage. As processes that utilize cyanobacteria as solar cell factories aim to maximize biomass and maintain the cells in the reactor system, the above-mentioned parameters are important properties of an effective and productive catalytic biofilm. Nevertheless, to enable targeted engineering, including genetic modification of this promising cyanobacterium and to develop it further to become an established chassis strain, the genome sequence of Tolypothrix sp. PCC 7712 is required.
In order to establish this organism as a photo-biotech workhorse, we set out to sequence and analyze the genome of Tolypothrix sp. PCC 7712. Surprisingly, its genome sequence showed a high similarity to the available genome sequence of Tolypothrix sp. PCC 7601, also known as Fremyella diplosiphon, which was originally isolated from a freshwater sample at another location. Here, we present a comparative genome analysis of these two strains, which present distinct deviations in gene composition and arrangement causing substantial physiological differences between both strains. As for Tolypothrix sp. PCC 7601 only a permanent draft genome was available, we also provided a completed genome sequence for this strain.
Materials and methods
Cultivation and maintenance of strains
Tolypothrix sp. PCC 7712 and Tolypothrix sp. PCC 7601 were obtained from Pasteur Culture Collection of Cyanobacteria (PCC). Both strains were maintained on an agar-solidified BG11 medium (Rippka et al., 1979) in growth chambers (INFORS) at 25 μE m−2 s−1 illumination at 30°C. For the experiments, they were transferred to liquid media, either to standard BG11 or BG11-0 (nitrate omitted) in 250 ml flasks with 20 ml of culture volume. They were incubated at constant illumination of 25 μE m−2 s−1 at 30°C without shaking.
For determining the ability of chromatic light adaptation, bacterial pre-cultures growing in BG11-0 or BG11 media were transferred to fresh media after 3 weeks, covered with gray, red, and blue foil, and cultivated for another 3 weeks under 200 μE cm−2 s−1 light-emitting diode (LED) illumination. Whole-cell absorption spectra were analyzed using a Cary 300 UV–Vis spectrophotometer (Agilent Technologies, Santa Clara, United States).
DNA isolation and quality assessment
Genomic DNA (gDNA) was isolated from Tolypothrix sp. PCC 7601 and PCC 7712 cells using established extraction protocols (Wilson, 2001) with few modifications. Briefly, 1 ml cell culture grown for 2 weeks was centrifuged at 13,000 g for 5 min, and the pellet was resuspended in 467 μl TE buffer. After addition of 100 μl Lysozyme (10 mg/ml), the pellet was resuspended by incubating it in a thermomixer (37°C, 15 min, 500 rpm), followed by the addition of 30 μl 10% SDS. The sample was incubated at the same conditions for another 15 min. Subsequently, 10 μl proteinase K (20 mg/ml) and 6 μl RNase (10 mg/ml) were added and the final solution was incubated for 1 h at 37°C while 400 rpm shaking. Then, pre-heated 80 μl 10% cetyltrimethylammonium bromide (CTAB) in 5 M NaCl solution was added together with 100 μl of pre-heated 5 M NaCl solution and incubated for 20 min at 65°C in a water bath. An equal volume of commercially obtained phenol/chloroform/isoamyl alcohol (25:24:1) was added to the mixture before it was centrifuged at 13,000 g for 10 min. The upper phase was transferred to a new tube followed by the addition of an equal volume of chloroform/isoamyl alcohol (24:1). After centrifugation (13,000 g, 10 min), the upper phase was transferred to a new tube, and 0.7 volume of isopropanol was added to the mixture. DNA was precipitated in isopropanol solution with centrifugation (13,000 g, 15 min). The pellet was washed with 70% ethanol. After removing the ethanol from the mixture by centrifugation, it was left open at 37°C for 1 h in order to remove residual ethanol from the gDNA pellet. The final pellet was resuspended with ddH2O and stored at 4°C. The quality and quantity of gDNA were checked via Nanodrop OneC Spectrophotometer (Thermo Fisher Scientific, Waltham, United States) at 260 and 280 nm.
Whole-genome sequencing, assembly, and annotation
For Illumina sequencing, isolated gDNA was fragmented and a sequencing library was prepared using the NEBNext® Ultra™ II DNA Library Prep Kit for Illumina® (New England Biolabs) according to the manufacturer’s instructions. Sequencing was performed on an Illumina MiSeq platform with a MiSeq Reagent Kit v3 (600-cycle; Illumina). Adapter sequences from Illumina raw sequencing data (PRJNA625426 for PCC7712 and PRJNA625641 for PCC7601) were trimmed using BBDuk of the bbmap suite v38.33.1 In addition, whole genome sequencing was performed by PCR-free Nanopore sequencing using an R9.4.1 flow cell on a MinIon MK1B device (Oxford Nanopore Technologies) controlled by MinKnow software release 19.12.5. gDNA isolates were prepared for sequencing using an SQK-LSK 109 Ligation Sequencing Kit in combination with an EXP-NBD104 Native Barcoding Expansion Kit according to the manufacturer’s instructions with the following exceptions. The incubation times for the end-repair step were increased to 30 min at room temperature and 30 min at 65°C. The time for the ligation step was increased to 60 min. Raw Nanopore sequence data were base called and demultiplexed using guppy version 3.6.0 and the provided high accuracy model. Adapter sequences were trimmed using Porechop version 0.2.4.2 Genomes were assembled using (a) unicycler v0.4.8 in hybrid mode using Illumina and Nanopore reads (Wick et al., 2017) and (b) Flye v2.8 using Nanopore reads only (Kolmogorov et al., 2020). The resulting assemblies were polished by medaka v1.0.3 using the respective Nanopore reads3 and by four rounds of pilon v1.22 (Walker et al., 2014) using the respective Illumina reads. Final assembly quality was checked with CheckM v1.011 using the lineage-specific workflow (Parks et al., 2015) and quast v5.0.2 (Gurevich et al., 2013). After completing assemblies, final sequences were submitted to Genbank (PRJNA625426 for PCC7712 and PRJNA625641 for PCC7601) and annotated via PGAP pipeline.4 Plasmid types were identified with MOB-Recon (Robertson and Nash, 2018).
Genome comparison
Both genomes were compared using Diffseq (Aggeli et al., 2018) to get information about genome variations between the two strains. MAUVE alignment (Darling et al., 2004) was applied to align chromosomes and plasmids to rearrange their initial locations for further Diffseq analyses. The genomes and generation of the annotation list were visualized using Geneious R10.0.5 (Kearse et al., 2012). To identify unique proteins BLAST RBH implemented to the Galaxy server5 was applied by comparing the encoding nucleotide sequences. After obtaining a “first match” list, respective sequences were extracted as a fasta file and BLAST RBH was used again against this extracted list for each genome. This process was repeated three times in total to avoid errors caused by multi-copy genes.
Pairwise comparison of different genome sets
Twenty-one other cyanobacterial species listed in Supplementary Table S3, which are either well-known in the biotechnology field or compose other Tolypothrix species were selected from the NCBI genome database to compare them to the newly sequenced two genomes used in this study. After accessing the selected genomes, dRep (Olm et al., 2017) MASH ANI clustering was applied in the Galaxy Server https://usegalaxy.eu/. The primary clustering ANI threshold was set to 90% and the secondary clustering ANI threshold was set to 99%.
Acetylene reduction assay
The acetylene reduction assay was applied as described (Stewart et al., 1968; Yoon and Golden, 2001) with some modifications to quantify nitrogenase activity in vivo. The respective cultures were adapted to nitrate-omitted media (BG11-0) for 1 week. Cells were harvested and transferred to 15 ml fresh BG11-0 medium in 20 ml GC vials; cell density was adjusted based on Chla content to 1.5 μM Chla which was measured according to the procedure described in a previous study (Zavřel et al., 2015). Briefly, a cell pellet obtained from 1 ml of cell culture was exposed to 1 ml of 100% methanol solution for 20 min in the dark followed by measuring the absorbance of the supernatant at 470, 665, and 720 nm wavelength via a visible range spectrophotometer (Libra S11, Biochrom, Cambridge, United Kingdom). After adjusting cell density via Chla measurement, the 5 ml headspace was filled with acetylene (0.5%), oxygen (20.9%), and nitrogen (78.6%) gas mixture. Vials were left for incubation at 50 μE m−2 s−1 for 24 h before measuring ethylene production via gas chromatography (TRACE 1310; Thermo Scientific, Waltham, USA). The device was equipped with a 30 m long TracePLOT TG-BOND Q+ column with 0.32 mm inner diameter and 10 μm film thickness (Thermo Scientific, Waltham, USA). The temperature of the flame ionization detector and the oven were adjusted to 300 and 60°C, respectively. A volume of 10 μl was injected via a Thermo TriPlus RSH autosampler. The flow rate of the carrier gas (nitrogen) was set to 10 ml min−1, with a total running time of 2 min. The calibration curve was set by injecting 0.5, 1, and 2 μl of ethylene gas.
Results
Characterization and classification of the whole genome sequence of Tolypothrix sp. PCC 7712
The genome sequence of Tolypothrix sp. PCC 7712 was obtained by using a hybrid sequencing approach based on two different techniques, namely Illumina and Nanopore. The full genome sequence of Tolypothrix sp. PCC 7712 consisted of one large contig covering 9 Mbp and 15 smaller contigs of 0.03–0.2 Mbp in length (Table 1). The overall GC content was 40.7%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-t001.jpg)
Table 1. Summary of the Tolypothrix sp. PCC 7712 genome as obtained from Illumina and Nanopore sequencing.
The largest contig (CP063785) of around 9 Mbp represents the chromosome, whereas several of the other contigs were predicted to be mobilizable and thus could represent plasmids. For the latter, we used the plasmid prediction tool MOB-Recon (Robertson and Nash, 2018). It should be noted that still several of these contigs might represent parts of the chromosome, especially when considering typically essential elements such as tRNAs. While the chromosome contig CP063785 contains 78 tRNA genes, 24 additional tRNA genes are found on contig CP063788, which possibly represents a plasmid. Nevertheless, plasmids that harbor tRNA arrays have also been reported for other bacterial species (Tran et al., 2016). Altogether, a set of 7,133 genes were identified on the largest contig and 978 genes, including 24 tRNA genes were found on the smaller contigs (Table 1).
Random sequences from the obtained PCC 7712 genome were manually analyzed using BLASTN. Remarkably, they appeared to be identical to an available genome sequence, namely that of Tolypothrix sp. PCC 7601 even though both strains were isolated independently and at different geographical locations. To exclude possible cross-contaminations, both strains were reordered from the PCC, gDNA was isolated and the sequencing approach was performed again for both strains (resulting in the read and assembly statistics given in Supplementary Table S1). The average nucleotide identity (ANI) of both obtained genome sequences was determined using the method FastANI (Jain et al., 2018). Indeed, both strains showed 99.98% ANI, which makes the two strains closely related but also indicates particular sequence alterations. As expected from high ANI, dRep (Olm et al., 2017) MASH ANI clustering with 21 selected cyanobacterial genomes showed only minor distance between Tolypothrix sp. PCC 7601 and PCC 7712, which is also reflected by a phylogenetic tree based on 16S rRNA comparison (Figure 1). The most related strains based on ANI clustering and the 16S rRNA based phylogenetic tree were Tolypothrix tenuis PCC 7101 followed by Tolypothrix sp. PCC 7910.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-g001.jpg)
Figure 1. (A) The tree shows genome sequence similarity of 23 cyanobacterial strains including PCC 7601 and PCC 7712 as calculated by dRep MASH ANI clustering. The dotted line shows the primary ANI clustering threshold, which was set to 90%. (B) The phylogenetic tree indicates the distances between these strains based on 16S rRNA homology.
Despite these high similarities of Tolypothrix sp. PCC 7712 and PCC 7601, one strain might have genes that the other does not have, or the order of genes might differ significantly. Therefore, gene order and genome synteny were examined using MAUVE to reveal potential differences. However, also in this respect both genomes showed high similarity in their genome arrangement (Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-g002.jpg)
Figure 2. (A) Chromosome alignment together with the direction of sequences (arrows) and (B) DNA plotters of the two chromosomes of Tolypothrix sp. PCC 7712 and PCC 7601 showing GC and GC skew in the middle followed by RNAs (blue and red lines), and the outer line shows the sequence itself with its complementary strand (green).
Comparative analysis of the genomes of Tolypothrix sp. PCC 7712 and PCC 7601 at global scale
Tolypothrix sp. PCC 7601 is the closest relative of strain PCC 7712. So far only an incomplete genome sequence of PCC 7601 containing 74 assembly gaps was available as a permanent draft deposited at GenBank (GCA_002368275.1). We, therefore, performed a hybrid-sequencing approach and achieved a complete circular chromosome of PCC 7601, as well as for strain PCC 7712 without any assembly gaps (Supplementary Figure S1), allowing for a detailed comparison of both strains using different bioinformatics tools.
To enable a detailed genome comparison of PCC 7712 and PCC 7601, including gene arrangement and composition, and identification of differences at gene as well as at nucleotide level, both strains’ genomes were annotated using the NCBI prokaryotic genome annotation pipeline (PGAP). This resulted in the annotation given in Supplementary Table S2. The annotation revealed that strain PCC 7712 lacks one tRNA, which was identified as a tRNA-Glu. This difference appeared to originate from the insertion of an IS701 family transposase (Figure 3A). Nevertheless, the number of encoded tRNAs already indicated that there is tRNA redundancy in Tolypothrix genomes, which also includes multiple tRNA-Glu with the same codon usage.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-g003.jpg)
Figure 3. Comparative analysis of the genomes of Tolypothrix sp. PCC 7712 and PCC 7601 at global scale. (A) Insertion of a transposase belonging to the IS701 family located in the gene encoding the tRNA-Glu in the PCC 7712 chromosome. (B) Comparison of contigs PCC 7601_CP064884 and PCC 7712_CP063789, which differ by a large insertion sequence containing a number of protein-coding regions like five hypothetical proteins (HP), a helix-turn-helix domain-containing protein, an ASCH domain-containing protein, a DDE-type integrase, and an ExeA family protein.
Moreover, a whole-genome comparison was carried out using Diffseq (Aggeli et al., 2018), which revealed differences between both genomes in gene composition as well as at the single nucleotide level. Both strains harbor several unique genes (Table 2). For instance, the possible plasmid CP063789 of strain PCC 7712 contains an insertion that significantly alters the gene composition compared to strain PCC 7601 (Figure 3B). It contains several unique genes such as for an activating signal cointegrator 1 homology (ASCH) domain-containing protein, a helix-turn-helix domain-containing protein, an AAA family ATPase, a helicase, a DUF4186 family protein, or an ExeA family protein and several hypothetical proteins. The region of insertion showed a high query coverage to proteins also found in closely related cyanobacterial strains of the genera Anabaena, Nostoc or Calothrix. ExeA family protein was reported to be involved in type II secretion systems (Vanderlinde et al., 2014), which might play a role in plasmid maintenance (Zhang et al., 2020). Further open reading frames (ORF) were uniquely identified either in the genome of strain PCC 7712 or PCC 7601 (see Table 2). Altogether, 80 unique CDS annotations were identified in strain PCC 7712 and 257 in strain PCC 7601. However, 199 CDS out of 257 (77%) in the genome of PCC 7601 and 41 out of 80 (51%) in the genome of PCC 7712 were identified as coding for hypothetical proteins (for a full list of unique ORFs see Supplementary Excel File S1).
Comparative analysis of both genomes at the single nucleotide level
Moreover, numerous single nucleotide polymorphisms (SNPs) were detected, mostly in the locations of mobile elements such as IS5, IS1634, and ISKra4 family transposases. Table 3 summarizes selected examples of the SNPs found in the genome of PCC 7712. However, most of the SNPs have been found in the non-coding regions of the genome rather than in the coding regions. These could become interesting if they affect cis- or trans- genetic elements, promotors or ribosomal binding sites. Some were also located within protein-coding regions, e.g., affecting an ATP/GTP-binding protein, HAMP domain-containing histidine kinase, and SDR family NAD(P)-dependent oxidoreductase all having a non-synonymous SNP, resulting in amino acid sequence alterations. In contrast, a gene for a tetratricopeptide repeat protein had synonymous SNP, not affecting the corresponding amino acid sequence.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-t003.jpg)
Table 3. Examples of SNPs and insertions found in the genome of PCC 7712 compared to PCC 7601. For a complete list please refer to Supplementary Excel File S1.
Tolypothrix strains PCC 7712 and PCC 7601 are closely related but differ in nitrogen fixation ability
Microscopy analyses showed that both PCC 7712 and PCC 7601 have similar morphology (Figure 4A). Both are filamentous and the filaments are generally entangled in some regions, which explains their hairy flocs appearance in suspended batch cultures (Figures 4A,B). Furthermore, pigmentation and the ability of chromatic adaptation are fairly similar in both strains (Figures 4B,C), which is also indicated by the sequence identity in the respective proteins involved such as phytochrome superfamily photoreceptors RcaE and DpxA (Bordowitz and Montgomery, 2008; Wiltbank and Kehoe, 2016). Tolypothrix sp. PCC 7601 is a model organism for investigating the mechanism behind complementary chromatic acclimation (CCA; Grossman et al., 2001), while this phenomenon has not been described for PCC 7712 so far, but is obviously also existing in this strain.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1042437/fmicb-13-1042437-HTML-r1/image_m/fmicb-13-1042437-g004.jpg)
Figure 4. (A) CLSM images of PCC 7601 and PCC 7712 cultures growing in BG11-0 or BG11 media; white circles indicate heterocyst cells, which are less fluorescent as they lack pigmentation. (B) Cultures growing in different light spectra in shake flasks in either BG11-0 or BG11. (C) Chromatic light adaptation. Whole-cell absorption spectra of cultures covered with different colored foils—red, blue, and gray.
A major physiological difference between the two strains is the ability of PCC 7712 to fix dinitrogen gas and thus being able to grow in a medium lacking an assimilable nitrogen source like nitrate (BG11-0). Under nitrate-limiting conditions only strain PCC 7712 formed terminal heterocysts (Figure 4A; white circles), which are specialized cells harboring the nitrogenase needed for N2 fixation. In contrast, strain PCC 7601 showed stressed filaments in BG11-0 lacking clear terminal heterocysts. Furthermore, no cell growth was observed which is consistent with previous reports (de Alda et al., 2004).
We performed an acetylene reduction assay to quantify the nitrogenase activity under aerobic conditions to confirm these observations. As expected, strain PCC 7601 did not show any activity, while a fairly high nitrogenase activity was detected in strain PCC 7712 converting up to 26.1 μmol H2 mgChla−1 h−1 under here used assay conditions. Strikingly, this major difference was not observable on the genomic level based on current state-of-the-art knowledge about nitrogen fixation in Tolypothrix. Even though there were several unique proteins, insertions, and SNPs observed in both genomes, none of the differences were identified in the regions currently known to be related to nitrogen fixation or heterocyst formation genes (Tables 2, 3).
Discussion
Tolypothrix sp. PCC 7712 is a novel candidate strain for cyanobacterial photo-biotechnology research due to its high biofilm biomass formation, low detachment rate, and high surface coverage in a recently introduced CBR system (Bozan et al., 2022). On the way to establish this organism as a potential chassis strain, the genome of Tolypothrix sp. PCC 7712 was sequenced and analyzed. We used the ANI approach, which, in contrast to the 16S rRNA-based classification, considers the total genome, and is regarded as a much more reliable tool for deciphering the degrees of bacterial relation (Jain et al., 2018). Even though strain PCC 7712 originated from a different geographical region it showed 99.9871% sequence identity with strain PCC 7601. Therefore, one could assume that both strains would have a similar phenotype despite minor differences in their genomes. This, for instance, is the case for complementary chromatic acclimation (CCA). Via CCA, cyanobacteria are able to change the arrangement of proteins and pigments in their phycobilisomes as well as in the accessory light-harvesting complexes connected to their photosystems when exposed to an altered light quality such as red or green light. While strain PCC 7601 has been widely used as a model strain for understanding the mechanism behind CCA (Kehoe and Gutu, 2006; Gutu and Kehoe, 2012; Montgomery, 2016), this phenomenon was so far not described for strain PCC 7712.
The most striking physiological difference was the ability of PCC7712 to fix dinitrogen gas via terminal heterocysts, while PCC 7601 was relying on dissolved inorganic nitrogen compounds like nitrate. This came as a surprise, as the genome of PCC 7601 harbors all essential genes for nitrogen fixation and heterocyst formation with 100% sequence identity to PCC 7712. Besides differences listed in Tables 2, 3 and Supplementary Excel File S1, there are no other differences on nucleotide level, which allow a conclusion on the reason for this observation.
Interestingly, a spontaneous revertant, Tolypothrix sp. PCC 7601/1 was reported previously, which is able to form heterocysts and fix dinitrogen (de Alda et al., 2004). Unfortunately, there are no sequencing data available for this strain but likely SNPs play a role in this regard. Moreover, it has been shown previously that the phycobilisome degradation protein NblA1 might play a role in the heterocyst differentiation of Tolypothrix sp. PCC 7601/1 but the underlying mechanism is not clearly identified yet (de Alda et al., 2004). This protein was found highly abundant under nitrogen limiting conditions in PCC 7601/1; however, it has been shown that it is also present in nitrogen replete conditions. Therefore, it was hypothesized that differential expression of the nblA1 gene could be involved in the complex heterocyst differentiation in filamentous cyanobacteria (de Alda et al., 2004). However, when comparing both PCC 7712 and PCC 7601, we could not determine any differences neither in the nblA1 gene nor in respective up- and downstream sequences. Moreover, there are several reports on small non-coding RNAs that are important during heterocyst differentiation in Nostoc sp. PCC 7120 (Brenes-Álvarez et al., 2020, 2022). Yet, a whole-genome comparison of PCC 7712 and PCC 7601 did not show any differences for small non-coding RNAs at the genome level. Nevertheless, it might be worth investigating the transcriptomes of PCC 7712 and PCC 7601 under different conditions to reveal non-coding RNAs and to unravel further details about the complex heterocyst formation process together with nitrogen fixation. Furthermore, regulatory effects based on epigenetic DNA modifications should also be considered as we identified a DNA cytosine methyltransferase (HGR01_40180) as a unique protein-coding sequence in the genome of PCC 7712 (Table 2). The presence of such a protein likely affects the DNA methylation pattern, which in turn could also affect the expression of specific genes similar to previous reports for Synechocystis (Gärtner et al., 2019). In this context, DNA methylation and its impact on nitrogen fixation were investigated in Trichodesmium erythraeum IMS101 (Walworth et al., 2017). However, in the latter cytosine methylation was not involved in the regulation of the nitrogen fixation (nif) genes directly. Nevertheless, the authors proposed a possible indirect effect on nitrogen fixation by influencing the expression of other nitrogen assimilatory genes. However, this is rather speculative and hence all these aspects need to be further investigated in PCC 7712 and PCC 7601. Meanwhile, there are also numerous hypothetical proteins defined uniquely in the PCC 7712 genome. Therefore, further comprehensive in silico analyses and physiological experiments are needed to understand the mechanism behind the difference in nitrogen fixation ability.
Conclusion
Here, we evaluated the genome of Tolypothrix sp. PCC 7712 which was previously introduced as a novel biofilm-forming strain with high application potential for photobiotechnological applications. Thereby, we discovered its high similarity to the genome of Tolypothrix sp. PCC 7601 (often referred to as Fremyella diplosiphon). Despite these significant similarities at the genome level, PCC 7712 was able to reduce N2, while PCC 7601 neither shows nitrogenase activity, nor heterocyst formation. This is surprising, as the genome encodes for all necessary proteins currently known to be involved in nitrogen fixation. Thereby, these strains may become interesting models for research focusing on understanding the process of nitrogen fixation in cyanobacteria. It also shows, that a classification based on genome comparison only does not necessarily end up with close relatives also exhibiting the same physiological behavior.
Data availability statement
The datasets presented in this study are available at the Genbank database under the following accession numbers: PRJNA625426 for PCC7712 and PRJNA625641 for PCC7601.
Author contributions
MB: data curation, experiment, and writing—original draft preparation. MB, DP, RK, UR, SK, and KB: methodology and writing—review and editing. KB and SK: conceptualization, methodology, formal analysis, supervision, and writing—review and editing. All authors contributed to the article and approved the submitted version.
Acknowledgments
We acknowledge the use of the facilities of the Centre for Biocatalysis (MiKat) at the Helmholtz Centre for Environmental Research, which is supported by European Regional Development Funds (EFRE, Europe funds Saxony), and the Helmholtz Association. Cloud computing facilities used for base calling of Nanopore reads and assembly were provided by the BMBF-funded de.NBI Cloud within the German Network for Bioinformatics Infrastructure (de.NBI; 031A537B, 031A533A, 031A538A, 031A533B, 031A535A, 031A537C, 031A534A, and 031A532B).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1042437/full#supplementary-material
Footnotes
1. ^https://sourceforge.net/projects/bbmap/
2. ^https://github.com/rrwick/Porechop
3. ^https://github.com/nanoporetech/medaka
References
Aggeli, D., Karas, V. O., Sinnott-Armstrong, N. A., Varghese, V., Shafer, R. W., Greenleaf, W. J., et al. (2018). Diff-seq: a high throughput sequencing-based mismatch detection assay for DNA variant enrichment and discovery. Nucleic Acids Res. 46:e42. doi: 10.1093/nar/gky022
Angermayr, S. A., Gorchs Rovira, A., and Hellingwerf, K. J. (2015). Metabolic engineering of cyanobacteria for the synthesis of commodity products. Trends Biotechnol. 33, 352–361. doi: 10.1016/j.tibtech.2015.03.009
Betterle, N., and Melis, A. (2019). Photosynthetic generation of heterologous terpenoids in cyanobacteria. Biotechnol. Bioeng. 116, 2041–2051. doi: 10.1002/bit.26988
Bharti, A., Velmourougane, K., and Prasanna, R. (2017). Phototrophic biofilms: diversity, ecology and applications. J. Appl. Phycol. 29, 2729–2744. doi: 10.1007/s10811-017-1172-9
Bordowitz, J. R., and Montgomery, B. L. (2008). Photoregulation of cellular morphology during complementary chromatic adaptation requires sensor-kinase-class protein Rca E in Fremyella diplosiphon. J. Bacteriol. 190, 4069–4074. doi: 10.1128/JB.00018-08
Bozan, M., Schmid, A., and Bühler, K. (2022). Evaluation of self-sustaining cyanobacterial biofilms for technical applications. Biofilms 4:100073. doi: 10.1016/j.bioflm.2022.100073
Brenes-Álvarez, M., Minguet, M., Vioque, A., and Muro-Pastor, A. M. (2020). NsiR1, a small RNA with multiple copies, modulates heterocyst differentiation in the cyanobacterium Nostoc sp. PCC 7120. Environ. Microbiol. 22, 3325–3338. doi: 10.1111/1462-2920.15103
Brenes-Álvarez, M., Vioque, A., and Muro-Pastor, A. M. (2022). The heterocyst-specific small RNA NsiR1 regulates the commitment to differentiation in Nostoc. Microbiol. Spectr. 10:e0227421. doi: 10.1128/spectrum.02274-21
Darling, A. C. E., Mau, B., Blattner, F. R., and Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704
de Alda, J. A. G. O., Lichtlé, C., Thomas, J.-C., and Houmard, J. (2004). Immunolocalization of NblA, a protein involved in phycobilisome turnover, during heterocyst differentiation in cyanobacteria. Microbiology 150, 1377–1384. doi: 10.1099/mic.0.26992-0
Garlapati, D., Chandrasekaran, M., Devanesan, A., Mathimani, T., and Pugazhendhi, A. (2019). Role of cyanobacteria in agricultural and industrial sectors: an outlook on economically important byproducts. Appl. Microbiol. Biotechnol. 103, 4709–4721. doi: 10.1007/s00253-019-09811-1
Gärtner, K., Klähn, S., Watanabe, S., Mikkat, S., Scholz, I., Hess, W. R., et al. (2019). Cytosine N4-methylation via M.Ssp6803II is involved in the regulation of transcription, fine- tuning of DNA replication and DNA repair in the cyanobacterium Synechocystis sp. PCC 6803. Front. Microbiol. 10:1233. doi: 10.3389/fmicb.2019.01233
Grossman, A. R., Bhaya, D., and He, Q. (2001). Tracking the light environment by cyanobacteria and the dynamic nature of light harvesting. J. Biol. Chem. 276, 11449–11452. doi: 10.1074/jbc.R100003200
Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086
Gutu, A., and Kehoe, D. M. (2012). Emerging perspectives on the mechanisms, regulation, and distribution of light color acclimation in cyanobacteria. Mol. Plant 5, 1–13. doi: 10.1093/mp/ssr054
Halan, B., Buehler, K., and Schmid, A. (2012). Biofilms as living catalysts in continuous chemical syntheses. Trends Biotechnol. 30, 453–465. doi: 10.1016/j.tibtech.2012.05.003
Hoschek, A., Heuschkel, I., Schmid, A., Bühler, B., Karande, R., and Bühler, K. (2019). Mixed-species biofilms for high-cell-density application of Synechocystis sp. PCC 6803 in capillary reactors for continuous cyclohexane oxidation to cyclohexanol. Bioresour. Technol. 282, 171–178. doi: 10.1016/j.biortech.2019.02.093
Jain, C., Rodriguez-R, L. M., Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9:5114. doi: 10.1038/s41467-018-07641-9
Jones, C. S., and Mayfield, S. P. (2012). Algae biofuels: versatility for the future of bioenergy. Curr. Opin. Biotechnol. 23, 346–351. doi: 10.1016/j.copbio.2011.10.013
Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199
Kehoe, D. M., and Gutu, A. (2006). Responding to color: the regulation of complementary chromatic adaptation. Annu. Rev. Plant Biol. 57, 127–150. doi: 10.1146/annurev.arplant.57.032905.105215
Kolmogorov, M., Bickhart, D. M., Behsaz, B., Gurevich, A., Rayko, M., Shin, S. B., et al. (2020). metaFlye: scalable long-read metagenome assembly using repeat graphs. Nat. Methods 17, 1103–1110. doi: 10.1038/s41592-020-00971-x
Montgomery, B. L. (2016). Mechanisms and fitness implications of photomorphogenesis during chromatic acclimation in cyanobacteria. J. Exp. Bot. 67, 4079–4090. doi: 10.1093/jxb/erw206
Nagappan, S., Bhosale, R., Duc Nguyen, D., Pugazhendhi, A., Tsai, P.-C., Chang, S. W., et al. (2020). Nitrogen-fixing cyanobacteria as a potential resource for efficient biodiesel production. Fuel 279:118440. doi: 10.1016/j.fuel.2020.118440
Olm, M. R., Brown, C. T., Brooks, B., and Banfield, J. F. (2017). dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 11, 2864–2868. doi: 10.1038/ismej.2017.126
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Posten, C. (2009). Design principles of photo-bioreactors for cultivation of microalgae. Eng. Life Sci. 9, 165–177. doi: 10.1002/elsc.200900003
Rippka, R., Deruelles, J., Waterbury, J. B., Herdman, M., and Stanier, R. Y. (1979). Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology 111, 1–61. doi: 10.1099/00221287-111-1-1
Robertson, J., and Nash, J. H. E. (2018). MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb. Genomics 4:e000206. doi: 10.1099/mgen.0.000206
Stewart, W. D. P., Fitzgerald, G. P., and Burris, R. H. (1968). Acetylene reduction by nitrogen-fixing blue-green algae. Arch. Mikrobiol. 62, 336–348. doi: 10.1007/BF00425639
Tran, T. T. T., Belahbib, H., Bonnefoy, V., and Talla, E. (2016). A comprehensive tRNA genomic survey unravels the evolutionary history of tRNA arrays in prokaryotes. Genome Biol. Evol. 8, 282–295. doi: 10.1093/gbe/evv254
Tsygankov, A. A. (2007). Nitrogen-fixing cyanobacteria: a review. Appl. Biochem. Microbiol. 43, 250–259. doi: 10.1134/S0003683807030040
Vanderlinde, E. M., Zhong, S., Li, G., Martynowski, D., Grochulski, P., and Howard, S. P. (2014). Assembly of the type two secretion system in Aeromonas hydrophila involves direct interaction between the periplasmic domains of the assembly factor ExeB and the secretin ExeD. PLoS One 9:e102038. doi: 10.1371/journal.pone.0102038
Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9:e112963. doi: 10.1371/journal.pone.0112963
Walworth, N. G., Hutchins, D. A., Dolzhenko, E., Lee, M. D., Fu, F., Smith, A. D., et al. (2017). Biogeographic conservation of the cytosine epigenome in the globally important marine, nitrogen-fixing cyanobacterium Trichodesmium. Environ. Microbiol. 19, 4700–4713. doi: 10.1111/1462-2920.13934
Wick, R. R., Judd, L. M., Gorrie, C. L., and Holt, K. E. (2017). Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 13:e1005595. doi: 10.1371/journal.pcbi.1005595
Wilson, K. (2001). Preparation of genomic DNA from bacteria. Curr. Protoc. Mol. Biol. 56, 2–4. doi: 10.1002/0471142727.mb0204s56
Wiltbank, L. B., and Kehoe, D. M. (2016). Two cyanobacterial photoreceptors regulate photosynthetic light harvesting by sensing teal, green, yellow, and red light. MBio 7, e02130–e02115. doi: 10.1128/mBio.02130-15
Xie, H., and Lindblad, P. (2022). Expressing 2-keto acid pathway enzymes significantly increases photosynthetic isobutanol production. Microb. Cell Factories 21:17. doi: 10.1186/s12934-022-01738-z
Yoon, H.-S., and Golden, J. W. (2001). PatS and products of nitrogen fixation control heterocyst pattern. J. Bacteriol. 183, 2605–2613. doi: 10.1128/JB.183.8.2605-2613.2001
Zavřel, T., Sinetova, M., and Červený, J. (2015). Measurement of chlorophyll a and carotenoids concentration in cyanobacteria. Bio-Protocol 5, 1–5. doi: 10.21769/BioProtoc.1467
Keywords: next-generation sequencing, comparative genomics, cyanobacteria, Tolypothrix, Fremyella diplosiphon
Citation: Bozan M, Popp D, Kallies R, da Rocha UN, Klähn S and Bühler K (2022) Whole-genome sequence of the filamentous diazotrophic cyanobacterium Tolypothrix sp. PCC 7712 and its comparison with non-diazotrophic Tolypothrix sp. PCC 7601. Front. Microbiol. 13:1042437. doi: 10.3389/fmicb.2022.1042437
Edited by:
Xuefeng Lu, Qingdao Institute of Bioenergy and Bioprocess Technology (CAS), ChinaReviewed by:
Satoru Watanabe, Tokyo University of Agriculture, JapanShengwei Hou, Southern University of Science and Technology, China
Copyright © 2022 Bozan, Popp, Kallies, da Rocha, Klähn and Bühler. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stephan Klähn, stephan.klaehn@ufz.de
†Present address: Denny Popp, Institute of Human Genetics, University of Leipzig Medical Center, Leipzig, Germany
‡These authors have contributed equally to this work and share last authorship