- Research note
- Open Access
A constrained SSU-rRNA phylogeny reveals the unsequenced diversity of photosynthetic Cyanobacteria (Oxyphotobacteria)
BMC Research Notesvolume 11, Article number: 435 (2018)
Cyanobacteria are an ancient phylum of prokaryotes that contain the class Oxyphotobacteria. This group has been extensively studied by phylogenomics notably because it is widely accepted that Cyanobacteria were responsible for the spread of photosynthesis to the eukaryotic domain. The aim of this study was to evaluate the fraction of the oxyphotobacterial diversity for which sequenced genomes are available for genomic studies. For this, we built a phylogenomic-constrained SSU rRNA (16S) tree to pinpoint unexploited clusters of Oxyphotobacteria that should be targeted for future genome sequencing, so as to improve our understanding of Oxyphotobacteria evolution.
We show that only a little fraction of the oxyphotobacterial diversity has been sequenced so far. Indeed 31 rRNA clusters of the 60 composing the photosynthetic Cyanobacteria have a fraction of sequenced genomes < 1%. This fraction remains low (min = 1%, median = 11.1%, IQR = 7.3%) within the remaining “sequenced” clusters that already contain some representative genomes. The “unsequenced” clusters are scattered across the whole Oxyphotobacteria tree, at the exception of very basal clades. Yet, these clades still feature some (sub)clusters without any representative genome. This last result is especially important, as these basal clades are prime candidate for plastid emergence.
Oxyphotobacteria are a class of the phylum Cyanobacteria, one of the most important groups of prokaryotes . This morphologically diversified class is the only group of Bacteria that can carry out oxygenic photosynthesis. It is commonly assumed that this bioenergetic process appeared in Oxyphotobacteria (or their ancestor) and caused the Great Oxidation Event (GOE) 2.4 billion years ago [2,3,4,5]. The atmospheric increase in free oxygen had a critical impact on early Earth evolution, including on Oxyphotobacteria themselves, with the appearance of new morphological features (e.g., akinetes) . Oxyphotobacteria were also responsible for the spread of photosynthesis to eukaryotic lineages, through an initial event of primary endosymbiosis, followed by higher-order endosymbiosis . In consequence of their great scientific interest, numerous phylogenies of Oxyphotobacteria have been published during the last decade (e.g., [5, 8,9,10,11,12,13,14,15,16,17,18,19]), using ever more genomes as it appears in public repositories. The diversity within Oxyphotobacteria is very high, both morphologically and genetically, and their GC-content ranges from 35 to 71 . Therefore, it is important to assess how well the existing genomes represent the whole diversity of the class and detect the lineages for which representative genomes are not yet available. Indeed, only 659 genomes of Oxyphotobacteria have been sequenced, whereas 6997 non-redundant rRNA sequences can be found in the SILVA database (as of October 2017). In this short communication, we built a SSU rRNA (16S) phylogeny to investigate the distribution of representative genomes across oxyphotobacterial diversity, so as to guide the selection of taxa for future sequencing projects.
We took advantage of 7603 different SSU rRNA (16S) sequences [6997 from SILVA 128 Ref NR , 466 from genome assemblies publicly available on the NCBI  and 140 from the BCCM/ULC public collection of Cyanobacteria (http://bccm.belspo.be/about-us/bccm-ulc)], leading to a dataset of 2302 unique Operational Taxonomic Units (OTUs) after dereplication at an identity threshold of 97.5%. Because multiple sequence alignment (MSA) quality is essential in single-gene analyses [23,24,25], we used an in-house hand-curated MSA as a starting point to carefully align all the OTUs. This MSA represented a high fraction of known Oxyphotobacteria diversity, as it completely included the broadly sampled dataset of Schirrmeister et al. . To alleviate the stochastic error plaguing single-gene analyses , we constrained our SSU rRNA (16S) phylogeny using a more reliable multi-gene phylogeny computed on a subset of the OTUs at hand. To this end, we used the phylogeny of Oxyphotobacteria developed in Cornet et al. , the largest ever-built in terms of conserved sites (> 170,000 unambiguously aligned amino-acid positions from 675 genes). This strategy allowed us to easily compare the genome diversity available for phylogenomic studies to the more completely sampled rRNA diversity of Oxyphotobacteria.
As shown in Fig. 1, 31 rRNA clusters have a fraction of sequenced genomes < 1% and are thus nearly completely ignored in phylogenomic studies. These “unsequenced” clusters are scattered across the whole oxyphotobacterial tree, at the exception of very basal clades (G, F, E) and the Oscillatoriales clade (A), which have higher fractions of representative genomes. The group C (C1, C2, C3) appears paraphyletic in our tree, and is especially rich in “unsequenced” clusters (14 out of 31), including the largest one (212 OTUs). In all the remaining “‘sequenced” clusters that contain representative genomes, the fraction of sequenced genomes remains low (min = 1%, median = 11.1%, IQR = 7.3%). It is thus clear that genomes have not been evenly sampled across the tree. To help researchers to select organisms for future sequencing projects, our final dataset and tree are made available (see below). In these files, SSU rRNA (16S) genes can be readily identified by their accession number, whether for SILVA/ULC strains or for rRNA genes predicted from genome assemblies.
The emergence of eukaryotic plastids from Oxyphotobacteria has been the subject of different hypotheses during the last decade, the most frequent being an early origin around clade E, composed of Acaryochloris and Thermosynechococcus species [17, 18]. By focusing on the base of our oxyphotobacterial phylogeny, Fig. 2 shows that in spite of the presence of a number of representative genomes, very basal clades (G, F, E) still feature some clusters without any genome (highlighted in yellow in Fig. 2). Therefore, these clusters should be targeted in priority for genome sequencing if one wants to resolve the early diversification of Oxyphotobacteria and the origin of plastids by assembling phylogenomic datasets more faithful to the true diversity.
Cyanobacterial (651 Oxyphotobacteria and 8 Melainabacteria) genome assemblies were downloaded from the NCBI FTP server on October 1st, 2017. SSU rRNA (16S) genes were predicted for all assemblies with RNAmmer v1.2 . Predicted rRNA genes were taxonomically classified with SINA v1.2.11 , using the release 128 of SILVA database composed of 1,922,213 SSU rRNA (16S) reference sequences . Among the 651 Oxyphotobacteria assemblies, 193 were devoid of Oxyphotobacteria rRNA genes (176 without any prediction and 17 with only non-cyanobacterial predicted genes). All Oxyphotobacteria rRNA genes (7001 sequences) of SILVA 128 Ref NR (non-redundant) were downloaded on October 1st, 2017. Four SSU rRNA sequences with a SILVA-provided incongruent lineage (e.g., Proteobacteria for Prochlorothrix hollandica PCC 9006) were deleted from the dataset, whereas the 466 sequences predicted with RNAmmer from genome assemblies were added, along with 140 new sequences from strains of the BBMC/ULC Cyanobacteria collection (http://bccm.belspo.be/about-us/bccm-ulc), hence leading to a total of 7603 SSU rRNA (16S) sequences (including 8 Melainabacteria).
The dataset was dereplicated with CD-HIT (cd-hit-est) v4.6 , using an identity threshold of 97.5%, an alignment coverage control of 0.0 for the longest sequence and 0.9 for the shortest sequence. The value of 97.5% for the identity threshold was chosen by following the OTU definition of Taton et al. . It probably underestimates the true number of different species, since Edgar  has recently shown that the bacterial species delineation threshold is closer to 99% for full-length sequences. However, for the purpose of this study aimed at pinpointing large patches of unsequenced diversity, a coarse-grained clustering was sufficient. The cluster file produced by CD-HIT was retrieved and 13 representative sequences initially selected by CD-HIT were replaced by genomic sequences from the same CD-HIT clusters, so as to be able to enforce the phylogenomic constraints from Cornet et al.  when inferring the rRNA tree. This process led to a dataset of 2302 OTUs.
The 2302 sequences were aligned using the BLASTN-based software “two-scalp” (D. Baurain; https://metacpan.org/release/Bio-MUST-Apps-TwoScalp), starting from a hand-curated MSA (1252 sites post-BMGE; 1.52% missing character states) of 1128 SSU rRNA (16S) sequences (L. Cornet; unpublished). The 1128 sequences had been selected in order to be as close as possible to the dataset of Schirrmeister et al. , which represented all the Oxyphotobacteria diversity known at that time. Unambiguously aligned sites were selected with BMGE v1.12 , using moderately severe settings (entropy cut-off 0.5, gap cut-off 0.2) and yielding a high-quality MSA of 1233 conserved sites with only 1.83% missing character states. In comparison, building a MSA from scratch on the same dataset with MAFFT v7.273 , using its automatic mode, led to only 1013 conserved sites (2.52% missing character states) after applying the same protocol to select the unambiguously aligned sites.
The tree was inferred with RAxML v8.1.17  under a GTR + Γ4 model with a 100× rapid bootstrap analysis, using a constraining tree of 75 different Oxyphotobacteria. The constraining tree was a phylogenomic tree based on a concatenation of 675 genes (> 170,000 unambiguously aligned amino-acid positions ). However, four organisms (Leptolyngbya sp. GCF_000733415.1, Oscillatoriales cyanobacterium GCF_000309945.1, Leptolyngbya sp. OTC1/1, Phormidium priestleyi ANT.PROGRESS2.5) had to be pruned before inference due to the lack of genuine rRNA genes in their assemblies. The rRNA tree was then rooted on the Melainabacteria clan (sister-group of Oxyphotobacteria, composed of eight sequences), manually collapsed using FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree) and further arranged in InkScape v0.92 . For each collapsed subtree, a “sequenced genome fraction” statistic was computed by taking the ratio between the number of OTUs having a publicly available genome in their CD-HIT cluster to the total number of OTUs in the subtree.
Among the 659 available Oxyphotobacteria genomes, 193 genomes were either completely devoid of rRNA genes (or only featured taxonomically incongruent rRNA genes) and were thus unusable for phylogeny (see “Methods” for details). Obviously, the results of this study might have been different if we had used these genomes to evaluate the sequenced genome fraction. The absence of rRNA genes is explained by the fact that the majority of these genomes (143 out of 193) were assembled with a metagenomic pipeline, according to NCBI metadata. This phenomenon was noticed in Cornet et al.  and is due to the frequent loss of rRNA sequences during metagenomic assembly . Given that neither the isolation source nor the assembly pipeline was easy to determine, we decided not to consider these 193 genomes in our analysis by precautionary principle. Hence, a significant number of these metagenomes were found in the genome section of RefSeq, even though it is not allowed. This decision is moreover consolidated by the observation that only a small number (about 20) of SSU rRNA (16S) sequences could be recovered after extensive searches on the NCBI using the “strain” names of these 193 genomes.
Operational Taxonomic Units
multiple sequence alignment
Soo RM, Hemp J, Parks DH, Fischer WW, Hugenholtz P. On the origins of oxygenic photosynthesis and aerobic respiration in Cyanobacteria. Science. 2017;355:1436–40.
Bekker A, Holland HD, Wang P-L, Rumble D, Stein HJ, Hannah JL, et al. Dating the rise of atmospheric oxygen. Nature. 2004;427:117–20.
Knoll AH. The geological consequences of evolution. Geobiology. 2003;1:3–14.
Kopp RE, Kirschvink JL, Hilburn IA, Nash CZ. The Paleoproterozoic snowball Earth: a climate disaster triggered by the evolution of oxygenic photosynthesis. Proc Natl Acad Sci USA. 2005;102:11131–6.
de Alda JAGO, Esteban R, Diago ML, Houmard J. The plastid ancestor originated among one of the major cyanobacterial lineages. Nat Commun. 2014;5:4937.
Gonzalez-Esquer CR, Smarda J, Rippka R, Axen SD, Guglielmi G, Gugger M, et al. Cyanobacterial ultrastructure in light of genomic sequence data. Photosynth Res. 2016;129:147–57.
Archibald JM. The puzzle of plastid evolution. Curr Biol. 2009;19:R81–8.
Rodríguez-Ezpeleta N, Brinkmann H, Burey SC, Roure B, Burger G, Löffelhardt W, et al. Monophyly of primary photosynthetic eukaryotes: green plants, red algae, and glaucophytes. Curr Biol. 2005;15:1325–30.
Blank CE, Sánchez-Baracaldo P. Timing of morphological and ecological innovations in the cyanobacteria—a key to understanding the rise in atmospheric oxygen. Geobiology. 2010;8:1–23.
Criscuolo A, Gribaldo S. Large-scale phylogenomic analyses indicate a deep origin of primary plastids within Cyanobacteria. Mol Biol Evol. 2011;28:3019–32.
Dagan T, Roettger M, Stucken K, Landan G, Koch R, Major P, et al. Genomes of Stigonematalean cyanobacteria (subsection V) and the evolution of oxygenic photosynthesis from prokaryotes to plastids. Genome Biol Evol. 2013;5:31–44.
Shih PM, Wu D, Latifi A, Axen SD, Fewer DP, Talla E, et al. Improving the coverage of the cyanobacterial phylum using diversity-driven genome sequencing. Proc Natl Acad Sci. 2013;110:1053–8.
Blank CE. Origin and early evolution of photosynthetic eukaryotes in freshwater environments: reinterpreting proterozoic paleobiology and biogeochemical processes in light of trait evolution. J Phycol. 2013;49:1040–55.
Li B, Lopes JS, Foster PG, Embley TM, Cox CJ. Compositional biases among synonymous substitutions cause conflict between gene and protein trees for plastid origins. Mol Biol Evol. 2014;31:1697–709.
Sánchez-Baracaldo P. Origin of marine planktonic cyanobacteria. Sci Rep. 2015;5. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4665016/. Accessed 7 Nov 2017.
Uyeda JC, Harmon LJ, Blank CE. A comprehensive study of cyanobacterial morphological and ecological evolutionary dynamics through deep geologic time. PLoS ONE. 2016;11:e0162539.
Ponce-Toledo RI, Deschamps P, López-García P, Zivanovic Y, Benzerara K, Moreira D. An early-branching freshwater cyanobacterium at the origin of plastids. Curr Biol. 2017;27:386–91.
Sánchez-Baracaldo P, Raven JA, Pisani D, Knoll AH. Early photosynthetic eukaryotes inhabited low-salinity habitats. Proc Natl Acad Sci. 2017;114:E7737–45.
Walter JM, Coutinho FH, Dutilh BE, Swings J, Thompson FL, Thompson CC. Ecogenomics and taxonomy of Cyanobacteria phylum. Front Microbiol. 2017;8:2132. https://doi.org/10.3389/fmicb.2017.02132.
Herdman M, Janvier M, Waterbury JB, Rippka R, Stanier RY, Mandel M. Deoxyribonucleic acid base composition of cyanobacteria. Microbiology. 1979;111:63–71.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, et al. The NCBI biosystems database. Nucleic Acids Res. 2010;38:D492–6.
Lunter G, Rocco A, Mimouni N, Heger A, Caldeira A, Hein J. Uncertainty in homology inferences: assessing and improving genomic sequence alignment. Genome Res. 2008;18:298–309.
Wong KM, Suchard MA, Huelsenbeck JP. Alignment uncertainty and genomic analysis. Science. 2008;319:473–6.
Dessimoz C, Gil M. Phylogenetic assessment of alignments reveals neglected tree signal in gaps. Genome Biol. 2010;11:R37.
Schirrmeister BE, Anisimova M, Antonelli A, Bagheri HC. Evolution of cyanobacterial morphotypes. Commun Integr Biol. 2011;4:424–7.
Delsuc F, Brinkmann H, Philippe H. Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet. 2005;6:nrg1603.
Cornet L, Bertrand AR, Hanikenne M, Javaux EJ, Wilmotte A, Baurain D. Metagenomic assembly of new (sub)arctic Cyanobacteria and their associated microbiome from non-axenic cultures. bioRxiv. 2018. https://doi.org/10.1101/287730.
Lagesen K, Hallin P, Rødland EA, Stærfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.
Pruesse E, Peplies J, Glöckner FO. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012;28:1823–9.
Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.
Taton A, Grubisic S, Brambilla E, Wit RD, Wilmotte A. Cyanobacterial diversity in natural and artificial microbial mats of Lake Fryxell (McMurdo Dry Valleys, Antarctica): a morphological and molecular approach. Appl Environ Microbiol. 2003;69:5157–69.
Edgar RC. Accuracy of taxonomy prediction for 16S rRNA and fungal ITS sequences. PeerJ. 2018;6:e4652.
Criscuolo A, Gribaldo S. BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol Biol. 2010;10:210.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.
Bah T. Inkscape: guide to a vector drawing program (Digital Short Cut). New York: Pearson Education; 2009.
Cornet L, Meunier L, Vlierberghe MV, Leonard RR, Durieu B, Lara Y, et al. Consensus assessment of the contamination level of publicly available cyanobacterial genomes. PLOS One. 2018. https://doi.org/10.1101/301788.
LC and DB designed the experiments, AW generated new SSU rRNA (16S) sequences for 140 strains from the BCCM/ULC public collection of Cyanobacteria, LC performed the computational analyses and drew the figures, EJJ and AW helped with the biogeological interpretations of the results, LC and DB wrote the manuscript, with the assistance of AW and EJJ. All authors read and approved the final manuscript.
We thank Yannick Lara (ULiège) for the cultivation and sequencing of some the BCCM/ULC strains. Rosa Gago is acknowledged for her help when drawing the figures. Finally, we are grateful to the editor and the anonymous reviewer for their insightful comments.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset (hand-curated MSA, final raw MSA, final MSA post-processed by BMGE, phylogenomic constraining tree, and final SSU rRNA ML tree) supporting the conclusions of this article is available as a compressed archive in the figshare repository through the following link (https://doi.org/10.6084/m9.figshare.6225641).
Consent for publication
Ethics approval and consent to participate
This work was supported by operating funds from FRS-FNRS (National Fund for Scientific Research of Belgium) (AW), the European Research Council Stg ELITE FP7/308074 (EJJ), the BELSPO project CCAMBIO (SD/BA/03A) (AW), the BELSPO Interuniversity Attraction Pole Planet TOPERS (EJJ and LC). LC was a FRIA fellow of the FRS-FNRS. AW is a Research Associate of the FRS-FNRS. Computational resources were provided through two grants to DB (University of Liège “Crédit de démarrage 2012” SFRD-12/04; FRS-FNRS “Crédit de recherche 2014” CDR J.0080.15).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.