Comparative plastid genomics of Synurophyceae: inverted repeat dynamics and gene content variation

Comparative structure and evolution of the organellar genomes of Padina usoehtunii (Dictyotales) with the brown algal crown radiation clade

BACKGROUND

Organellar genomes have become increasingly essential for studying genetic diversity, phylogenetics, and evolutionary histories of seaweeds. The order Dictyotales (Dictyotophycidae), a highly diverse lineage within the Phaeophyceae, is long-term characterized by a scarcity of organellar genome datasets compared to orders of the brown algal crown radiation (Fucophycidae).

RESULTS

We sequenced the organellar genomes of Padina usoehtunii, a representative of the order Dictyotales, to investigate the structural and evolutionary differences by comparing to five other major brown algal orders. Our results confirmed previously reported findings that the rate of structural rearrangements in chloroplast genomes is higher than that in mitochondria, whereas mitochondrial sequences exhibited a higher substitution rate compared to chloroplasts. Such evolutionary patterns contrast with land plants and green algae. The expansion and contraction of the inverted repeat (IR) region in the chloroplast correlated with the changes in the number of boundary genes. Specifically, the size of the IR region influenced the position of the boundary gene rpl21, with complete rpl21 genes found within the IR region in Dictyotales, Sphacelariales and Ectocarpales, while the rpl21 genes in Desmarestiales, Fucales, and Laminariales span both the IR and short single copy (SSC) regions. The absence of the rbcR gene in the Dictyotales may indicate an endosymbiotic transfer from the chloroplast to the nuclear genome. Inversion of the SSC region occurred at least twice in brown algae. Once in a lineage only represented by the Ectocarpales in the present study and once in a lineage only represented by the Fucales. Photosystem genes in the chloroplasts experienced the strongest signature of purifying selection, while ribosomal protein genes in both chloroplasts and mitochondria underwent a potential weak purifying selection.

CONCLUSIONS

Variations in chloroplast genome structure among different brown algal orders are evolutionarily linked to their phylogenetic positions in the Phaeophyceae tree. Chloroplast genomes harbor more structural rearrangements than the mitochondria, despite mitochondrial genes exhibiting faster mutation rates. The position and the change in the number of boundary genes likely shaped the IR regions in the chloroplast, and the produced structural variability is important mechanistically to create gene diversity in brown algal chloroplast.

Organellar genome comparisons of Sargassum polycystum and S. plagiophyllum (Fucales, Phaeophyceae) with other Sargassum species

BACKGROUND

Sargassum polycystum C. Agardh and Sargassum plagiophyllum C. Agardh are inhabitants of tropical coastal areas, their populations are negatively influenced by global warming and marine environment changes. The mitochondrial and chloroplast genomes of these species have not been sequenced.

RESULTS

The mitochondrial genomes of S. polycystum and S. plagiophyllum were 34,825 bp and 34,862 bp, respectively, and their corresponding chloroplast genomes were 124,493 bp and 124,536 bp, respectively. The mitochondrial and chloroplast genomes of these species share conserved synteny, sequence regions and gene number when compared with the organellar genomes of other Sargassum species. Based on sequence analysis of 35 protein-coding genes, we deduced that S. polycystum and S. plagiophyllum were closely related with S. ilicifolium; these species diverged approximately 0.3 million years ago (Ma; 0.1-0.53 Ma) during the Pleistocene period (0.01-2.59 Ma). Rates of synonymous and non-synonymous substitutions in the mitochondrial genome of the Sargassum genus were 3 times higher than those in the chloroplast genome. In the mitochondrial genome, rpl5, rpl31 and rps11 had the highest synonymous substitution rates. In the chloroplast genome, psaE, rpl14 and rpl27 had the highest synonymous substitution rates.

CONCLUSIONS

Phylogenetic analysis confirms the close relationship between the two sequenced species and S. ilicifolium. Both synonymous and non-synonymous substitution rates show significant divergence between the group of mitochondrial genomes versus the group of chloroplast genomes. The deciphering of complete mitochondrial and chloroplast genomes is significant as it advances our understanding of the evolutionary and phylogenetic relationships between species of brown seaweeds.

Evolutionary Dynamics and Lateral Gene Transfer in Raphidophyceae Plastid Genomes

The Raphidophyceae is an ecologically important eukaryotic lineage of primary producers and predators that inhabit marine and freshwater environments worldwide. These organisms are of great evolutionary interest because their plastids are the product of eukaryote-eukaryote endosymbiosis. To obtain deeper insight into the evolutionary history of raphidophycean plastids, we sequenced and analyzed the plastid genomes of three freshwater and three marine species. Our comparison of these genomes, together with the previously reported plastid genome of Heterosigma akashiwo, revealed unexpected variability in genome structure. Unlike the genomes of other analyzed species, the plastid genome of Gonyostomum semen was found to contain only a single rRNA operon, presumably due to the loss of genes from the inverted repeat (IR) region found in most plastid genomes. In contrast, the marine species Fibrocapsa japonica contains the largest IR region and overall plastid genome for any raphidophyte examined thus far, mainly due to the presence of four large gene-poor regions and foreign DNA. Two plastid genes, tyrC in F. japonica and He. akashiwo and serC in F. japonica, appear to have arisen via lateral gene transfer (LGT) from diatoms, and several raphidophyte open reading frames are demonstrably homologous to sequences in diatom plasmids and plastid genomes. A group II intron in the F. japonica psbB gene also appears to be derived by LGT. Our results provide important insights into the evolutionary history of raphidophyte plastid genomes via LGT from the plastids and plasmid DNAs of diatoms.

Olisthodiscus represents a new class of Ochrophyta

The phylogenetic diversity of Ochrophyta, a diverse and ecologically important radiation of algae, is still incompletely understood even at the level of the principal lineages. One taxon that has eluded simple classification is the marine flagellate genus Olisthodiscus. We investigated Olisthodiscus luteus K-0444 and documented its morphological and genetic differences from the NIES-15 strain, which we described as Olisthodiscus tomasii sp. nov. Phylogenetic analyses of combined 18S and 28S rRNA sequences confirmed that Olisthodiscus constitutes a separate, deep, ochrophyte lineage, but its position could not be resolved. To overcome this problem, we sequenced the plastid genome of O. luteus K-0444 and used the new data in multigene phylogenetic analyses, which suggested that Olisthodiscus is a sister lineage of the class Pinguiophyceae within a broader clade additionally including Chrysophyceae, Synchromophyceae, and Eustigmatophyceae. Surprisingly, the Olisthodiscus plastid genome contained three genes, ycf80, cysT, and cysW, inherited from the rhodophyte ancestor of the ochrophyte plastid yet lost from all other ochrophyte groups studied so far. Combined with nuclear genes for CysA and Sbp proteins, Olisthodiscus is the only known ochrophyte possessing a plastidial sulfate transporter SulT. In addition, the finding of a cemA gene in the Olisthodiscus plastid genome and an updated phylogenetic analysis ruled out the previously proposed hypothesis invoking horizontal cemA transfer from a green algal plastid into Synurales. Altogether, Olisthodiscus clearly represents a novel phylogenetically distinct ochrophyte lineage, which we have proposed as a new class, Olisthodiscophyceae.

Extreme Enlargement of the Inverted Repeat Region in the Plastid Genomes of Diatoms from the Genus Climaconeis

We sequenced the plastid genomes of three diatoms from the genus Climaconeis, including two strains formerly designated as Climaconeis scalaris. At 208,097 and 216,580 bp, the plastid genomes of the latter strains are the largest ever sequenced among diatoms and their increased size is explained by the massive expansion of the inverted repeat region. Important rearrangements of gene order were identified among the two populations of Climaconeis cf. scalaris. The other sequenced Climaconeis chloroplast genome is 1.5 times smaller compared with those of the Climaconeis cf. scalaris strains and it features an usual quadripartite structure. The extensive structural changes reported here for the genus Climaconeis are compared with those previously observed for other algae and plants displaying large plastid genomes.

High Sequence Divergence but Limited Architectural Rearrangements in Organelle Genomes of Cyanophora (Glaucophyta) Species

Cyanophora is the glaucophyte model taxon. Following the sequencing of the nuclear genome of C. paradoxa, studies based on single organelle and nuclear molecular markers revealed previously unrecognized species diversity within this glaucophyte genus. Here, we present the complete plastid (ptDNA) and mitochondrial (mtDNA) genomes of C. kugrensii, C. sudae, and C. biloba. The respective sizes and coding capacities of both ptDNAs and mtDNAs are conserved among Cyanophora species with only minor differences due to specific gene duplications. Organelle phylogenomic analyses consistently recover the species C. kugrensii and C. paradoxa as a clade and C. sudae and C. biloba as a separate group. The phylogenetic affiliations of the four Cyanophora species are consistent with architectural similarities shared at the organelle genomic level. Genetic distance estimations from both organelle sequences are also consistent with phylogenetic and architecture evidence. Comparative analyses confirm that the Cyanophora mitochondrial genes accumulate substitutions at 3-fold higher rates than plastid counterparts, suggesting that mtDNA markers are more appropriate to investigate glaucophyte diversity and evolutionary events that occur at a population level. The study of complete organelle genomes is becoming the standard for species delimitation and is particularly relevant to study cryptic diversity in microbial groups.

Substantial intraspecific genome size variation in golden-brown algae and its phenotypic consequences

BACKGROUND AND AIMS

While nuclear DNA content variation and its phenotypic consequences have been well described for animals, vascular plants and macroalgae, much less about this topic is known regarding unicellular algae and protists in general. The dearth of data is especially pronounced when it comes to intraspecific genome size variation. This study attempts to investigate the extent of intraspecific variability in genome size and its adaptive consequences in a microalgal species.

METHODS

Propidium iodide flow cytometry was used to estimate the absolute genome size of 131 strains (isolates) of the golden-brown alga Synura petersenii (Chrysophyceae, Stramenopiles), identified by identical internal transcribed spacer (ITS) rDNA barcodes. Cell size, growth rate and genomic GC content were further assessed on a sub-set of strains. Geographic location of 67 sampling sites across the Northern hemisphere was used to extract climatic database data and to evaluate the ecogeographical distribution of genome size diversity.

KEY RESULTS

Genome size ranged continuously from 0.97 to 2.02 pg of DNA across the investigated strains. The genome size was positively associated with cell size and negatively associated with growth rate. Bioclim variables were not correlated with genome size variation. No clear trends in the geographical distribution of strains of a particular genome size were detected, and strains of different genome size occasionally coexisted at the same locality. Genomic GC content was significantly associated only with genome size via a quadratic relationship.

CONCLUSIONS

Genome size variability in S. petersenii was probably triggered by an evolutionary mechanism operating via gradual changes in genome size accompanied by changes in genomic GC content, such as, for example, proliferation of transposable elements. The variation was reflected in cell size and relative growth rate, possibly with adaptive consequences.

Comparative Plastid Genomics of Non-Photosynthetic Chrysophytes: Genome Reduction and Compaction

Spumella-like heterotrophic chrysophytes are important eukaryotic microorganisms that feed on bacteria in aquatic and soil environments. They are characterized by their lack of pigmentation, naked cell surface, and extremely small size. Although Spumella-like chrysophytes have lost their photosynthetic ability, they still possess a leucoplast and retain a plastid genome. We have sequenced the plastid genomes of three non-photosynthetic chrysophytes, Spumella sp. Baeckdong012018B8, Pedospumella sp. Jangsampo120217C5 and Poteriospumella lacustris Yongseonkyo072317C3, and compared them to the previously sequenced plastid genome of "Spumella" sp. NIES-1846 and photosynthetic chrysophytes. We found the plastid genomes of Spumella-like flagellates to be generally conserved with respect to genome structure and housekeeping gene content. We nevertheless also observed lineage-specific gene rearrangements and duplication of partial gene fragments at the boundary of the inverted repeat and single copy regions. Most gene losses correspond to genes for proteins involved in photosynthesis and carbon fixation, except in the case of petF. The newly sequenced plastid genomes range from ~55.7 kbp to ~62.9 kbp in size and share a core set of 45 protein-coding genes, 3 rRNAs, and 32 to 34 tRNAs. Our results provide insight into the evolutionary history of organelle genomes via genome reduction and gene loss related to loss of photosynthesis in chrysophyte evolution.

Dictyochophyceae Plastid Genomes Reveal Unusual Variability in Their Organization

Dictyochophyceae (silicoflagellates) are unicellular freshwater and marine algae (Heterokontophyta, stramenopiles). Despite their abundance in global oceans and potential ecological significance, discovered in recent years, neither nuclear nor organellar genomes of representatives of this group were sequenced until now. Here, we present the first complete plastid genome sequences of Dictyochophyceae, obtained from four species: Dictyocha speculum, Rhizochromulina marina, Florenciella parvula and Pseudopedinella elastica. Despite their comparable size and genetic content, these four plastid genomes exhibit variability in their organization: plastid genomes of F. parvula and P. elastica possess conventional quadripartite structure with a pair of inverted repeats, R. marina instead possesses two direct repeats with the same orientation and D. speculum possesses no repeats at all. We also observed a number of unusual traits in the plastid genome of D. speculum, including expansion of the intergenic regions, presence of an intron in the otherwise non-intron-bearing psaA gene, and an additional copy of the large subunit of RuBisCO gene (rbcL), the last of which has never been observed in any plastid genome. We conclude that despite noticeable gene content similarities between the plastid genomes of Dictyochophyceae and their relatives (pelagophytes, diatoms), the number of distinctive features observed in this lineage strongly suggests that additional taxa require further investigation.

Draft genome assembly and transcriptome sequencing of the golden algae Hydrurus foetidus (Chrysophyceae)

Hydrurus foetidus is a freshwater chrysophyte alga. It thrives in cold rivers in polar and high alpine regions. It has several morphological traits reminiscent of single-celled eukaryotes, but can also form macroscopic thalli. Despite its ability to produce polyunsaturated fatty acids, its life under cold conditions and its variable morphology, very little is known about its genome and transcriptome. Here, we present an extensive set of next-generation sequencing data, including genomic short reads from Illumina sequencing and long reads from Nanopore sequencing, as well as full length cDNAs from PacBio IsoSeq sequencing and a small RNA dataset (smaller than 200 bp) sequenced with Illumina. The genome sequences were combined  to produce an assembly consisting of 5069 contigs, with a total assembly size of 171 Mb and a 77% BUSCO completeness. The new data generated here may contribute to a better understanding of the evolution and ecological roles of chrysophyte algae, as well as to resolve the branching patterns at a larger phylogenetic scale.

Draft genome assembly and transcriptome sequencing of the golden algae Hydrurus foetidus (Chrysophyceae)

Hydrurusfoetidus is a freshwater chrysophyte alga. It thrives in cold rivers in polar and high alpine regions. It has several morphological traits reminiscent of single-celled eukaryotes, but can also form macroscopic thalli. Despite its ability to produce polyunsaturated fatty acids, its life under cold conditions and its variable morphology, very little is known about its genome and transcriptome. Here, we present an extensive set of next-generation sequencing data, including genomic short reads from Illumina sequencing and long reads from Nanopore sequencing, as well as full length cDNAs from PacBio IsoSeq sequencing and a small RNA dataset (smaller than 200 bp) sequenced with Illumina. The genome sequences were combined  to produce an assembly consisting of 5069 contigs, with a total assembly size of 171 Mb and a 77% BUSCO completeness. The new data generated here may contribute to a better understanding of the evolution and ecological roles of chrysophyte algae, as well as to resolve the branching patterns at a larger phylogenetic scale.

Principles of plastid reductive evolution illuminated by nonphotosynthetic chrysophytes

The division of life into producers and consumers is blurred by evolution. For example, eukaryotic phototrophs can lose the capacity to photosynthesize, although they may retain vestigial plastids that perform other essential cellular functions. Chrysophyte algae have undergone a particularly large number of photosynthesis losses. Here, we present a plastid genome sequence from a nonphotosynthetic chrysophyte, "Spumella" sp. NIES-1846, and show that it has retained a nearly identical set of plastid-encoded functions as apicomplexan parasites. Our transcriptomic analysis of 12 different photosynthetic and nonphotosynthetic chrysophyte lineages reveals remarkable convergence in the functions of these nonphotosynthetic plastids, along with informative lineage-specific retentions and losses. At one extreme, Cornospumella fuschlensis retains many photosynthesis-associated proteins, although it appears to have lost the reductive pentose phosphate pathway and most plastid amino acid metabolism pathways. At the other extreme, Paraphysomonas lacks plastid-targeted proteins associated with gene expression and all metabolic pathways that require plastid-encoded partners, indicating a complete loss of plastid DNA in this genus. Intriguingly, some of the nucleus-encoded proteins that once functioned in the expression of the Paraphysomonas plastid genome have been retained. These proteins were likely to have been dual targeted to the plastid and mitochondria of the chrysophyte ancestor, and are uniquely targeted to the mitochondria in Paraphysomonas Our comparative analyses provide insights into the process of functional reduction in nonphotosynthetic plastids.