Parallel evolution of highly conserved plastid genome architecture in red seaweeds and seed plants

The complete chloroplast genome of Erodium cicutarium (Linnaeus) l' Héritier ex Aiton 1789 (Geraniaceae): genome characterization and phylogenetic consideration

Erodium cicutarium is an annual herbaceous plant valued for its applications in traditional medicine. However, the chloroplast genome of E. cicutarium has yet to be reported. In this study, we assembled chloroplast genomes of Erodium cicutarium using Illumina sequencing reads. The chloroplast genome was 114,652 bp long, harbored 111 complete genes, and its overall GC content was 39.1%. In Maximum Likelihood (ML) and Bayesian inference (BI) trees, the 13 Erodium species divided into three main clades, with E. cicutarium and E. carvifolium forming a monophyletic group, suggesting a close relationship between the two species. The E. cicutarium cp genome presented in this study lays a good foundation for the Erodium.

Of sea, rivers and symbiosis: Diversity, systematics, biogeography and evolution of the deeply diverging florideophycean order Hildenbrandiales (Rhodophyta)

The Hildenbrandiales, a typically saxicolous red algal order, is an early diverging florideophycean group with global significance in marine and freshwater ecosystems across diverse temperature zones. To comprehensively elucidate the diversity, phylogeny, biogeography, and evolution of this order, we conducted a thorough re-examination employing molecular data derived from nearly 700 specimens. Employing a species delimitation method, we identified Evolutionary Species Units (ESUs) within the Hildenbrandiales aiming to enhance our understanding of species diversity and generate the first time-calibrated tree and ancestral area reconstruction for this order. Mitochondrial cox1 and chloroplast rbcL markers were used to infer species boundaries, and subsequent phylogenetic reconstructions involved concatenated sequences of cox1, rbcL, and 18S rDNA. Time calibration of the resulting phylogenetic tree used a fossil record from a Triassic purportedly freshwater Hildenbrandia species and three secondary time points from the literature. Our species delimitation analysis revealed an astounding 97 distinct ESUs, quintupling the known diversity within this order. Our time-calibration analysis placed the origin of Hildenbrandiales (crown age) in the Ediacaran period, with freshwater species emerging as a monophyletic group during the later Permian to early Triassic. Phylogenetic reconstructions identified seven major clades, experiencing early diversification during the Silurian to Carboniferous period. Two major evolutionary events-colonization of freshwater habitats and obligate systemic symbiosis with a marine fungus-marked this order, leading to significant morphological alterations without a commensurate increase in species diversification. Despite the remarkable newly discovered diversity, the extant taxon diversity appears relatively constrained when viewed against an evolutionary timeline spanning over 800 million years. This limitation may stem from restricted geographic sampling or the prevalence of asexual reproduction. However, species richness estimation and rarefaction analyses suggest a substantially larger diversity yet to be uncovered-potentially four times greater. These findings drastically reshape our understanding of the deeply diverging florideophycean order Hildenbrandiales species diversity, and contribute valuable insights into this order's evolutionary history and ecological adaptations. Supported by phylogenetic, ecological and morphological evidence, we established the genus Riverina gen. nov. to accommodate freshwater species of Hildenbrandiales, which form a monophyletic clade in our analyses. This marks the first step toward refining the taxonomy of the Hildenbrandiales, an order demanding thorough revisions, notably with the creation of several genera to address the polyphyletic status of Hildenbrandia. However, the limited diagnostic features pose a challenge, necessitating a fresh approach to defining genera. A potential solution lies in embracing a molecular systematic perspective, which can offer precise delineations of taxonomic boundaries.

Genome-Wide Identification and Expression Analysis of the Sweet Cherry Whirly Gene Family

Sweet cherry (Prunus avium) is one of the economically valuable horticultural fruit trees and it is widely cultivated throughout the world. Whirly (WHY) genes are a unique gene family with few members and have important biological functions in plant growth, development, and response to abiotic stress. This study utilized whole-genome identification to conduct a comprehensive analysis of the WHY genes in sweet cherry and examined their transcription levels in different tissues and under abiotic stress to explore their functions. Two WHY genes were identified in the sweet cherry genome and named PavWHY1 and PavWHY2, respectively, based on their homology with those in Arabidopsis thaliana. Both genes have theoretical isoelectric points greater than seven and are hydrophilic proteins, suggesting that they may be localized in plastids. The two genes are evolutionarily classified into two categories, with large differences in gene structure, and highly similar protein tertiary structures, and both have conserved domains of WHY. PavWHY1 and PavWHY2 are collinear with AtWHY1 and AtWHY2, respectively. The promoter sequence contains cis-acting elements related to hormones and abiotic stress, which are differentially expressed during flower bud differentiation, fruit development, and cold accumulation. qRT-PCR showed that PavWHY1 and PavWHY2 were differentially expressed in flower and fruit development and responded to low temperature and exogenous ABA treatment. The recombinant plasmid pGreenII-0800-Luc with the promoters of these two genes can activate luciferase expression in tobacco. Protein interaction predictions indicate that these gene products may interact with other proteins. This study reveals the molecular features, evolutionary relationships, and expression patterns of sweet cherry WHY genes, and investigates the activities of their promoters, which lays the foundation for further exploration of their biological functions and provides new insights into the WHY gene family in Rosaceae.

From dusk till dawn: cell cycle progression in the red seaweed Gracilariopsis chorda (Rhodophyta)

The conserved eukaryotic functions of cell cycle genes have primarily been studied using animal/plant models and unicellular algae. Cell cycle progression and its regulatory components in red (Rhodophyta) seaweeds are poorly understood. We analyzed diurnal gene expression data to investigate the cell cycle in the red seaweed Gracilariopsis chorda. We identified cell cycle progression and transitions in G. chorda which are induced by interactions of key regulators such as E2F/DP, RBR, cyclin-dependent kinases, and cyclins from dusk to dawn. However, several typical CDK inhibitor proteins are absent in red seaweeds. Interestingly, the G1-S transition in G. chorda is controlled by delayed transcription of GINS subunit 3. We propose that the delayed S phase entry in this seaweed may have evolved to minimize DNA damage (e.g., due to UV radiation) during replication. Our results provide important insights into cell cycle-associated physiology and its molecular mechanisms in red seaweeds.

Taxonomic reinvestigation of the genus Tetradesmus (Scenedesmaceae; Sphaeropleales) based on morphological characteristics and chloroplast genomes

The genus Tetradesmus (Scenedesmaceae; Sphaeropleales) comprises one of the most abundant green algae in freshwater environments. It includes morphologically diverse species that exhibit bundle-like, plane-arranged coenobia, and unicells, because several different Scenedesmus-like groups were integrated into this genus based on phylogenetic analysis. Nevertheless, there is no clear information regarding the phylogenetic relationship of Tetradesmus species, determined using several marker genes, because of low phylogenetic support and insufficient molecular data. Currently, genome information is available from diverse taxa, which could provide high-resolution evolutionary relationships. In particular, phylogenetic studies using chloroplast genomes demonstrated the potential to establish high-resolution phylogenetic relationships. However, only three chloroplast genomes are available from the genus Tetradesmus. In this study, we newly generated 9 chloroplast genomes from Tetradesmus and constructed a high-resolution phylogeny using a concatenated alignment of 69 chloroplast protein sequences. We also report one novel species (T. lancea), one novel variety (T. obliquus var. spiraformis), and two novel formae (T. dissociatus f. oviformis, T. obliquus f. rectilineare) within the genus Tetradesmus based on morphological characteristics (e.g., cellular arrangements and coenobial types) and genomic features (e.g., different exon-intron structures in chloroplast genomes). Moreover, we taxonomically reinvestigated the genus Tetradesmus based on these results. Altogether, our study can provide a comprehensive understanding of the taxonomic approaches for investigating this genus.

Diurnal Rhythms in the Red Seaweed Gracilariopsis chorda are Characterized by Unique Regulatory Networks of Carbon Metabolism

Cellular and physiological cycles are driven by endogenous pacemakers, the diurnal and circadian rhythms. Key functions such as cell cycle progression and cellular metabolism are under rhythmic regulation, thereby maintaining physiological homeostasis. The photoreceptors phytochrome and cryptochrome, in response to light cues, are central input pathways for physiological cycles in most photosynthetic organisms. However, among Archaeplastida, red algae are the only taxa that lack phytochromes. Current knowledge about oscillatory rhythms is primarily derived from model species such as Arabidopsis thaliana and Chlamydomonas reinhardtii in the Viridiplantae, whereas little is known about these processes in other clades of the Archaeplastida, such as the red algae (Rhodophyta). We used genome-wide expression profiling of the red seaweed Gracilariopsis chorda and identified 3,098 rhythmic genes. Here, we characterized possible cryptochrome-based regulation and photosynthetic/cytosolic carbon metabolism in this species. We found a large family of cryptochrome genes in G. chorda that display rhythmic expression over the diurnal cycle and may compensate for the lack of phytochromes in this species. The input pathway gates regulatory networks of carbon metabolism which results in a compact and efficient energy metabolism during daylight hours. The system in G. chorda is distinct from energy metabolism in most plants, which activates in the dark. The green lineage, in particular, land plants, balance water loss and CO2 capture in terrestrial environments. In contrast, red seaweeds maintain a reduced set of photoreceptors and a compact cytosolic carbon metabolism to thrive in the harsh abiotic conditions typical of intertidal zones.

Red macroalgae in the genomic era

Rhodophyta (or red algae) are a diverse and species-rich group that forms one of three major lineages in the Archaeplastida, a eukaryotic supergroup whose plastids arose from a single primary endosymbiosis. Red algae are united by several features, such as relatively small intron-poor genomes and a lack of cytoskeletal structures associated with motility like flagella and centrioles, as well as a highly efficient photosynthetic capacity. Multicellular red algae (or macroalgae) are one of the earliest diverging eukaryotic lineages to have evolved complex multicellularity, yet despite their ecological, evolutionary, and commercial importance, they have remained a largely understudied group of organisms. Considering the increasing availability of red algal genome sequences, we present a broad overview of fundamental aspects of red macroalgal biology and posit on how this is expected to accelerate research in many domains of red algal biology in the coming years.

Characterization and organelle genome sequencing of Pyropia species from Myanmar

Pyropia is a genus comprising red algae of the Bangiaceae family that is commonly found in intertidal zones worldwide. However, understanding of Pyropia species that are prone to tropical regions remains limited despite recent breakthroughs in genomic research. Within the realm of Pyropia species thriving in tropical regions, P. vietnamensis stands out as a widely recognized species. In this study, we aimed to investigate Pyropia species in the southwest coast of Myanmar using physiological and molecular approaches, culture-based analyses, chloroplast rbcL and nuclear SSU gene sequencing, and whole chloroplast and mitochondrial genome sequencing. Physiological analysis showed that the Myanmar samples were more heat-tolerant than their Japanese counterparts, including those of subtropical origin. Additionally, molecular characterization revealed that the Myanmar samples were closely related to P. vietnamensis from India. This study is the first to sequence the chloroplast and mitochondrial genomes of Pyropia species from tropical regions. A unique deletion event was observed within a ribosomal RNA gene cluster in the chloroplast genome of the studied Pyropia species, which is a deviation from the usual characteristics of most Pyropia species. This study improves current understanding of the physiological and molecular characteristics of this comparatively understudied Pyropia species that grows in tropical regions.

Revised classification of the Cyanidiophyceae based on plastid genome data with descriptions of the Cavernulicolales ord. nov. and Galdieriales ord. nov. (Rhodophyta)

The Cyanidiophyceae, an extremophilic red algal class, is distributed worldwide in extreme environments. Species grow either in acidic hot environments or in dim light conditions (e.g., "cave Cyanidium"). The taxonomy and classification systems are currently based on morphological, eco-physiological, and molecular phylogenetic characters; however, previous phylogenetic results showed hidden diversity of the Cyanidiophyceae and suggested a revision of the classification system. To clarify phylogenetic relationships within this red algal class, we employ a phylogenomic approach based on 15 plastomes (10 new) and 15 mitogenomes (seven new). Our phylogenies show consistent relationships among four lineages (Galdieria, "cave Cyanidium", Cyanidium, and Cyanidioschyzon lineages). Each lineage is distinguished by organellar genome characteristics. The "cave Cyanidium" lineage is a distinct clade that diverged after the Galdieria clade but within a larger monophyletic clade that included the Cyanidium and Cyanidioschyzon lineages. Because the "cave Cyanidium" lineage is a mesophilic lineage that differs substantially from the other three thermoacidophilic lineages, we describe it as a new order (Cavernulicolales). Based on this evidence, we reclassified the Cyanidiophyceae into four orders: Cyanidiales, Cyanidioschyzonales, Cavernulicolales ord. nov., and Galdieriales ord. nov. The genetic distance among these four orders is comparable to, or greater than, the distances found between other red algal orders and subclasses. Three new genera (Cavernulicola, Gronococcus, Sciadococcus), five new species (Galdieria javensis, Galdieria phlegrea, Galdieria yellowstonensis, Gronococcus sybilensis, Sciadococcus taiwanensis), and a new nomenclatural combination (Cavernulicola chilensis) are proposed.

Algae obscura: The potential of rare species as model systems

Model organism research has provided invaluable knowledge about foundational biological principles. However, most of these studies have focused on species that are in high abundance, easy to cultivate in the lab, and represent only a small fraction of extant biodiversity. Here, we present three examples of rare algae with unusual features that we refer to as "algae obscura." The Cyanidiophyceae (Rhodophyta), Glaucophyta, and Paulinella (rhizarian) lineages have all transitioned out of obscurity to become models for fundamental evolutionary research. Insights have been gained into the prevalence and importance of eukaryotic horizontal gene transfer, early Earth microbial community dynamics, primary plastid endosymbiosis, and the origin of Archaeplastida. By reviewing the research that has come from the exploration of these organisms, we demonstrate that underappreciated algae have the potential to help us formulate, refine, and substantiate core hypotheses and that such organisms should be considered when establishing future model systems.

Phylogenomics and taxon-rich phylogenies of new and historical specimens shed light on the systematics of Hypnea (Cystocloniaceae, Rhodophyta)

Cystocloniacae is a highly diverse family of Rhodophyta, including species of ecological and economic importance, whose phylogeny remains largely unresolved. Species delimitation is unclear, particularly in the most speciose genus, Hypnea, and cryptic diversity has been revealed by recent molecular assessments, especially in the tropics. Here, we carried out the first phylogenomic investigation of Cystocloniaceae, focused on the genus Hypnea, inferred from chloroplast and mitochondrial genomes including taxa sampled from new and historical collections. In this work, molecular synapomorphies (gene losses, InDels and gene inversions) were identified to better characterize clades in our congruent organellar phylogenies. We also present taxon-rich phylogenies based on plastid and mitochondrial markers. Molecular and morphological comparisons of historic collections with contemporary specimens revealed the need for taxonomic updates in Hypnea, the synonymization of H. marchantae to a later heterotypic synonym of H. cervicornis and the description of three new species: H. davisiana sp. nov., H. djamilae sp. nov. and H. evaristoae sp. nov.

The origin and early evolution of plants

Plant (archaeplastid) evolution has transformed the biosphere, but we are only now beginning to learn how this took place through comparative genomics, phylogenetics, and the fossil record. This has illuminated the phylogeny of Archaeplastida, Viridiplantae, and Streptophyta, and has resolved the evolution of key characters, genes, and genomes - revealing that many key innovations evolved long before the clades with which they have been casually associated. Molecular clock analyses estimate that Streptophyta and Viridiplantae emerged in the late Mesoproterozoic to late Neoproterozoic, whereas Archaeplastida emerged in the late-mid Palaeoproterozoic. Together, these insights inform on the coevolution of plants and the Earth system that transformed ecology and global biogeochemical cycles, increased weathering, and precipitated snowball Earth events, during which they would have been key to oxygen production and net primary productivity (NPP).

Life cycle and reproduction dynamics of Bangiales in response to environmental stresses

Red algae of the order Bangiales are notable for exhibiting flexible promotion of sexual and asexual reproductive processes by environmental stresses. This flexibility indicates that a trade-off between vegetative growth and reproduction occurs in response to environmental stresses that influence the timing of phase transition within the life cycle. Despite their high phylogenetic divergence, both filamentous and foliose red alga in the order Bangiales exhibit a haploid-diploid life cycle, with a haploid leafy or filamentous gametophyte (thallus) and a diploid filamentous sporophyte (conchocelis). Unlike haploid-diploid life cycles in other orders, the gametophyte in Bangiales is generated independently of meiosis; the regulation of this generation transition is not fully understood. Based on transcriptome and gene expression analyses, the originally proposed biphasic model for alternation of generations in Bangiales was recently updated to include a third stage. Along with the haploid gametophyte and diploid sporophyte, the triphasic framework recognizes a diploid conchosporophyte-a conchosporangium generated on the conchocelis-phase and previously considered to be part of the sporophyte. In addition to this sexual life cycle, some Bangiales species have an asexual life cycle in which vegetative cells of the thallus develop into haploid asexual spores, which are then released from the thallus to produce clonal thalli. Here, we summarize the current knowledge of the triphasic life cycle and life cycle trade-off in Neopyropia yezoensis and 'Bangia' sp. as model organisms for the Bangiales.

Extremophilic red algae as models for understanding adaptation to hostile environments and the evolution of eukaryotic life on the early earth

Extremophiles have always garnered great interest because of their exotic lifestyles and ability to thrive at the physical limits of life. In hot springs environments, the Cyanidiophyceae red algae are the only photosynthetic eukaryotes able to live under extremely low pH (0-5) and relatively high temperature (35ºC to 63ºC). These extremophiles live as biofilms in the springs, inhabit acid soils near the hot springs, and form endolithic populations in the surrounding rocks. Cyanidiophyceae represent a remarkable source of knowledge about the evolution of extremophilic lifestyles and their genomes encode specialized enzymes that have applied uses. Here we review the evolutionary origin, taxonomy, genome biology, industrial applications, and use of Cyanidiophyceae as genetic models. Currently, Cyanidiophyceae comprise a single order (Cyanidiales), three families, four genera, and nine species, including the well-known Cyanidioschyzon merolae and Galdieria sulphuraria. These algae have small, gene-rich genomes that are analogous to those of prokaryotes they live and compete with. There are few spliceosomal introns and evidence exists for horizontal gene transfer as a driver of local adaptation to gain access to external fixed carbon and to extrude toxic metals. Cyanidiophyceae offer a variety of commercial opportunities such as phytoremediation to detoxify contaminated soils or waters and exploitation of their mixotrophic lifestyles to support the efficient production of bioproducts such as phycocyanin and floridosides. In terms of exobiology, Cyanidiophyceae are an ideal model system for understanding the evolutionary effects of foreign gene acquisition and the interactions between different organisms inhabiting the same harsh environment on the early Earth. Finally, we describe ongoing research with C. merolae genetics and summarize the unique insights they offer to the understanding of algal biology and evolution.

Long-read sequencing characterizes mitochondrial and plastid genome variants in Arabidopsis msh1 mutants

The abundant repeats in plant mitochondrial genomes can cause rapid genome rearrangements and are also a major obstacle in short-read sequencing studies. Nuclear-encoded proteins such as MSH1 are known to suppress the generation of repeat-associated mitochondrial genome variants, but our understanding of these mechanisms has been constrained by the limitations of short-read technologies. Here, we used highly accurate long-read sequencing (PacBio HiFi) to characterize mitochondrial and plastid genome variants in Arabidopsis thaliana msh1 mutant individuals. The HiFi reads provided a global view of recombination dynamics with detailed quantification of parental and crossover recombination products for both large and small repeats. We found that recombination breakpoints were distributed relatively evenly across the length of repeated sequences and detected widespread internal exchanges of sequence variants between pairs of imperfect repeats in the mitochondrial genome of msh1 mutants. Long-read assemblies of mitochondrial genomes from seven other A. thaliana wild-type accessions differed by repeat-mediated structural rearrangements similar to those observed in msh1 mutants, but they were all in a simple low-heteroplasmy state. The Arabidopsis plastid genome generally lacks small repeats and exhibited a very different pattern of variant accumulation in msh1 mutants compared with the mitochondrial genome. Our data illustrate the power of HiFi technology in studying repeat-mediated recombination in plant organellar genomes and improved the sequence resolution for recombinational processes suppressed by MSH1. Plant organellar genomes can undergo rapid rearrangements. Long-read sequencing provides a detailed and quantitative view of mitochondrial and plastid genome variants normally suppressed by MSH1, advancing our understanding of plant organellar genome dynamics.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

Organelle Genome Variation in the Red Algal Genus Ahnfeltia (Florideophyceae)

The agarophyte Ahnfeltia (Ahnfeltiales, Rhodophyta) is a globally widespread genus with 11 accepted species names. Two of the most widespread species in this genus, A. plicata and A. fastigiata, may have diverged genetically due to past geographic changes and subsequent geographic isolation. To investigate this genomic and genetic diversity, we generated new plastid (ptDNAs) and mitochondrial genomes (mtDNAs) of these Ahnfeltia species from four different regions (A. plicata - Chile and UK and A. fastigiata - Korea and Oregon). Two architecture variations were found in the Ahnfeltia genomes: in ptDNA of A. fastigiata Oregon, the hypothetical pseudogene region was translocated, likely due to recombination with palindromic repeats or a gene transfer from a red algal plasmid. In mtDNA of A. fastigiata Korea, the composition of the group II intronic ORFs was distinct from others suggesting different scenarios of gain and loss of group II intronic ORFs. These features resulted in genome size differences between the two species. Overall gene contents of organelle genomes of Ahnfeltia were conserved. Phylogenetic analysis using concatenated genes from ptDNAs and mtDNAs supported the monophyly of the Ahnfeltiophycidae. The most probable individual gene trees showed that the Ahnfeltia populations were genetically diversified. These trees, the cox1 haplotype network, and a dN/dS analysis all supported the theory that these Ahnfeltia populations have diversified genetically in accordance with geographic distribution.

Phylogenomics, divergence time estimation and trait evolution provide a new look into the Gracilariales (Rhodophyta)

The Gracilariales is a highly diverse, widely distributed order of red algae (Rhodophyta) that forms a well-supported clade. Aside from their ecological importance, species of Gracilariales provide important sources of agarans and possess bioactive compounds with medicinal and pharmaceutical use. Recent phylogenetic analyses from a small number of genes have greatly advanced our knowledge of evolutionary relationships in this clade, yet several key nodes were not especially well resolved. We assembled a phylogenomic data set containing 79 nuclear genes, 195 plastid genes, and 24 mitochondrial genes from species representing all three major Gracilariales lineages, including: Melanthalia, Gracilariopsis, and Gracilaria sensu lato. This data set leads to a fully-resolved phylogeny of Gracilariales, which is highly-consistent across genomic compartments. In agreement with previous findings, Melanthalia obtusata was sister to a clade including Gracilaria s.l. and Gracilariopsis, which were each resolved as well-supported clades. Our results also clarified the long-standing uncertainty about relationships in Gracilaria s.l., not resolved in single and multi-genes approaches. We further characterized the divergence time, organellar genome architecture, and morphological trait evolution in Gracilarales to better facilitate its taxonomic treatment. Gracilariopsis and Gracilaria s.l. are comparable taxonomic ranks, based on the overlapping time range of their divergence. The genomic structure of plastid and mitochondria is highly conserved within each clade but differs slightly among these clades in gene contents. For example, the plastid gene petP is lost in Gracilaria s.l. and the mitochondrial gene trnH is in different positions in the genome of Gracilariopsis and Gracilaria s.l. Our analyses of ancestral character evolution provide evidence that the main characters used to delimitate genera in Gracilariales, such as spermatangia type and features of the cystocarp's anatomy, overlap in subclades of Gracilaria s.l. We discuss the taxonomy of Gracilariales in light of these results and propose an objective and practical classification, which is in agreement with the criteria of monophyly, exclusive characters, predictability and nomenclatural stability.

Complete plastid genome of Cumathamnion serrulatum (Ceramiales, Rhodophyta)

We report the complete plastid genome of Cumathamnion serrulatum, also known as Delesseria serrulata. The plastid genome was 174,192 bp in size. Annotation showed there were 193 protein coding genes, three ribosomal RNAs, and 29 transfer RNAs. One intron was found, and the GC content was 27.2%. The maximum likelihood tree with the concatenated 177 plastid coding genes showed a strong monophyletic relationship to Membranoptera spp. within the Ceramiales.

Divergence times and plastid phylogenomics within the intron-rich order Erythropeltales (Compsopogonophyceae, Rhodophyta)

The advent of high-throughput sequencing (HTS) has allowed for the use of large numbers of coding regions to produce robust phylogenies. These phylogenies have been used to highlight relationships at ancient diversifications (subphyla, class) and highlight the evolution of plastid genome structure. The Erythropeltales are an order in the Compsopogonophyceae, a group with unusual plastid genomes but with low taxon sampling. We use HTS to produce near complete plastid genomes of all genera, and multiple species within some genera, to produce robust phylogenies to investigate character evolution, dating of divergence in the group, and plastid organization, including intron patterns. Our results produce a fully supported phylogeny of the genera in the Erythropeltales and suggest that morphologies (upright versus crustose) have evolved multiple times. Our dated phylogeny also indicates that the order is very old (~800 Ma), with diversification occurring after the ice ages of the Cryogenian period (750-635 Ma). Plastid gene order is congruent with phylogenetic relationships and suggests that genome architecture does not change often. Our data also highlight the abundance of introns in the plastid genomes of this order. We also produce a nearly complete plastid genome of Tsunamia transpacifica (Stylonematophyceae) to add to the taxon sampling of genomes of this class. The use of plastid genomes clearly produces robust phylogenetic relationships that can be used to infer evolutionary events, and increased taxon sampling, especially in less well-known red algal groups, will provide additional insights into their evolution.

Plastid genomes and phylogenomics of liverworts (Marchantiophyta): Conserved genome structure but highest relative plastid substitution rate in land plants

With some 7300 species of small nonvascular spore-producing plants, liverworts represent one of the major lineages of land plants. Although multi-locus molecular phylogenetic studies have elucidated relationships of liverworts at different taxonomic categories, the backbone phylogeny of liverworts is still to be fully resolved, especially for the placement of Ptilidiales and the relationships within Jungermanniales and Marchantiales. Here, we provided phylogenomic inferences of liverworts based on 42 newly sequenced and 24 published liverwort plastid genomes representing all but two orders of liverworts, and characterized the evolution of the plastome in liverworts. The structure of the plastid genome is overall conserved across the phylogeny of liverworts, with only two structural variants detected from simple thalloids, besides 18 out of 43 liverwort genera showing intron variations in their plastomes. Complex thalloid liverworts maintain the most plastid genes, and seem to undergo fewer gene deletions and pseudogenization events than other liverworts. Plastid phylogenetic inferences yielded mostly robustly supported relationships, and consistently resolved Ptilidiales as the sister to Porellales. The relative ratio of silent substitutions across the three genetic compartments (i.e., 1:15:10, for mitochondrial:plastid:nuclear) suggests that liverwort plastid genes have the potential to evolve faster than their nuclear counterparts, unlike in any other major land plant lineages where the mutation rate of nuclear genes overwhelm those of their plastid and mitochondrial counterparts.

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.

Comparative Genome Analysis Reveals Cyanidiococcus gen. nov., A New Extremophilic Red Algal Genus Sister to Cyanidioschyzon (Cyanidioschyzonaceae, Rhodophyta)

The taxonomic placement of strains belonging to the extremophilic red alga Galdieria maxima has been controversial due to the inconsistent phylogenetic position inferred from molecular phylogenetic analyses. Galdieria maxima nom. inval. was classified in this genus based on morphology and molecular data in the early work, but some subsequent molecular phylogenetic analyses have inferred strains of G. maxima to be closely related to the genus Cyanidioschyzon. To address this controversy, an isolated strain identified as G. maxima using the rbcL gene sequence as the genetic barcode was examined using a comprehensive analysis across morphological, physiological, and genomic traits. Herein are reported the chloroplast-, mitochondrion-, and chromosome-level nuclear genome assemblies. Comparative analysis of orthologous gene clusters and genome arrangements suggested that the genome structure of this strain was more similar to that of the generitype of Cyanidioschyzon, C. merolae than to the generitype of Galdieria, G. sulphuraria. While the ability to uptake various forms of organic carbon for growth is an important physiological trait of Galdieria, this strain was identified as an ecologically obligate photoautotroph (i.e., the inability to utilize the natural concentrations of organic carbons) and lacked various gene models predicted as sugar transporters. Based on the genomic, morphological, and physiological traits, we propose this strain to be a new genus and species, Cyanidiococcus yangmingshanensis. Re-evaluation of the 18S rRNA and rbcL gene sequences of the authentic strain of G. maxima, IPPAS-P507, with those of C. yangmingshanensis suggests that the rbcL sequences of "G. maxima" deposited in GenBank correspond to misidentified isolates.

OGDA: a comprehensive organelle genome database for algae

Algae are the oldest taxa on Earth, with an evolutionary relationship that spans prokaryotes (Cyanobacteria) and eukaryotes. A long evolutionary history has led to high algal diversity. Their organelle DNAs are characterized by uniparental inheritance and a compact genome structure compared with nuclear genomes; thus, they are efficient molecular tools for the analysis of gene structure, genome structure, organelle function and evolution. However, an integrated organelle genome database for algae, which could enable users to both examine and use relevant data, has not previously been developed. Therefore, to provide an organelle genome platform for algae, we have developed a user-friendly database named Organelle Genome Database for Algae (OGDA, http://ogda.ytu.edu.cn/). OGDA contains organelle genome data either retrieved from several public databases or sequenced in our laboratory (Laboratory of Genetics and Breeding of Marine Organism [MOGBL]), which are continuously updated. The first release of OGDA contains 1055 plastid genomes and 755 mitochondrial genomes. Additionally, a variety of applications have been integrated into this platform to analyze the structural characteristics, collinearity and phylogeny of organellar genomes for algae. This database represents a useful tool for users, enabling the rapid retrieval and analysis of information related to organellar genomes for biological discovery.

Comparison of the transcriptomes of different life history stages of the freshwater Rhodophyte Thorea hispida

Thorea hispida exclusively inhabits freshwater environments and is characterized by a triphasic life history. In this study, the organelle genomes and transcriptomes of different life history stages of T. hispida were examined using next generation sequencing. The chloroplast and mitochondrial genomes of the chantransia stage were 175,747 and 25,411 bp in length, respectively. The chantransia stage was highly similar to the gametophyte stage based on comparisons of organelle genomes and phylogenetic reconstruction. Transcriptomic comparisons of two stages found that ribosome-related genes were the most up-regulated in the gametophyte stage of T. hispida. Seven meiosis-specific genes, including SPO11 initiator of meiotic double-stranded breaks(spo11), meiotic nuclear divisions 1(mnd1), RAD51 recombinase(rad51), mutS homolog 4(msh4), mutS homolog 5(msh5), REC8 meiotic recombination protein(rec8), and DNA helicase Mer3(mer3), were differentially regulated between the two life history stages. The organelle genomes and transcriptomes from T. hispida provided in this study will be valuable for future studies of freshwater red algae.

Genetic and Biochemical Reconstitution of Bromoform Biosynthesis in Asparagopsis Lends Insights into Seaweed Reactive Oxygen Species Enzymology

Marine macroalgae, seaweeds, are exceptionally prolific producers of halogenated natural products. Biosynthesis of halogenated molecules in seaweeds is inextricably linked to reactive oxygen species (ROS) signaling as hydrogen peroxide serves as a substrate for haloperoxidase enzymes that participate in the construction these halogenated molecules. Here, using red macroalga Asparagopsis taxiformis, a prolific producer of the ozone depleting molecule bromoform, we provide the discovery and biochemical characterization of a ROS-producing NAD(P)H oxidase from seaweeds. This discovery was enabled by our sequencing of Asparagopsis genomes, in which we find the gene encoding the ROS-producing enzyme to be clustered with genes encoding bromoform-producing haloperoxidases. Biochemical reconstitution of haloperoxidase activities establishes that fatty acid biosynthesis can provide viable hydrocarbon substrates for bromoform production. The ROS production haloperoxidase enzymology that we describe here advances seaweed biology and biochemistry by providing the molecular basis for decades worth of physiological observations in ROS and halogenated natural product biosyntheses.

Plastid-targeted forms of restriction endonucleases enhance the plastid genome rearrangement rate and trigger the reorganization of its genomic architecture

Plant cells have acquired chloroplasts (plastids) with a unique genome (ptDNA), which developed during the evolution of endosymbiosis. The gene content and genome structure of ptDNAs in land plants are considerably stable, although those of algal ptDNAs are highly varied. Plant cells seem, therefore, to be intolerant of any structural or organizational changes in the ptDNA. Genome rearrangement functions as a driver of genomic evolutionary divergence. Here, we aimed to create various types of rearrangements in the ptDNA of Arabidopsis genomes using plastid-targeted forms of restriction endonucleases (pREs). Arabidopsis plants expressing each of the three specific pREs, i.e., pTaqI, pHinP1I, and pMseI, were generated; they showed the leaf variegation phenotypes associated with impaired chloroplast development. We confirmed that these pREs caused double-stranded breaks (DSB) at their recognition sites in ptDNAs. Genome-wide analysis of ptDNAs revealed that the transgenic lines exhibited a large number of rearrangements such as inversions and deletions/duplications, which were dominantly repaired by microhomology-mediated recombination and microhomology-mediated end-joining, and less by non-homologous end-joining. Notably, pHinP1I, which recognized a small number of sites in ptDNA, induced drastic structural changes, including regional copy number variations throughout ptDNAs. In contrast, the transient expression of either pTaqI or pMseI, whose recognition site numbers were relatively larger, resulted in small-scale changes at the whole genome level. These results indicated that DSB frequencies and their distribution are major determinants in shaping ptDNAs.

The draft mitochondrial genome of Magnolia biondii and mitochondrial phylogenomics of angiosperms

The mitochondrial genomes of flowering plants are well known for their large size, variable coding-gene set and fluid genome structure. The available mitochondrial genomes of the early angiosperms show extreme genetic diversity in genome size, structure, and sequences, such as rampant HGTs in Amborella mt genome, numerous repeated sequences in Nymphaea mt genome, and conserved gene evolution in Liriodendron mt genome. However, currently available early angiosperm mt genomes are still limited, hampering us from obtaining an overall picture of the mitogenomic evolution in angiosperms. Here we sequenced and assembled the draft mitochondrial genome of Magnolia biondii Pamp. from Magnoliaceae (magnoliids) using Oxford Nanopore sequencing technology. We recovered a single linear mitochondrial contig of 967,100 bp with an average read coverage of 122 × and a GC content of 46.6%. This draft mitochondrial genome contains a rich 64-gene set, similar to those of Liriodendron and Nymphaea, including 41 protein-coding genes, 20 tRNAs, and 3 rRNAs. Twenty cis-spliced and five trans-spliced introns break ten protein-coding genes in the Magnolia mt genome. Repeated sequences account for 27% of the draft genome, with 17 out of the 1,145 repeats showing recombination evidence. Although partially assembled, the approximately 1-Mb mt genome of Magnolia is still among the largest in angiosperms, which is possibly due to the expansion of repeated sequences, retention of ancestral mtDNAs, and the incorporation of nuclear genome sequences. Mitochondrial phylogenomic analysis of the concatenated datasets of 38 conserved protein-coding genes from 91 representatives of angiosperm species supports the sister relationship of magnoliids with monocots and eudicots, which is congruent with plastid evidence.

Relative expression analysis of light-harvesting genes in the freshwater alga Lympha mucosa (Batrachospermales, Rhodophyta)

Members of the freshwater red algal order Batrachospermales are often described as shade-adapted. Nevertheless, recent ecophysiological studies have demonstrated species-level differences in acclimation to a range of irradiances. Lympha mucosa occurs in open and shaded portions of temperate streams and is abundant during summer months, suggesting it tolerates high and low irradiances. Specimens of L. mucosa were collected from open (sun-acclimated) or shaded (shade-acclimated) sites and exposed to low (<20 μmol photons · m^-2 · s^-1 ) or high (220 μmol photon · m^-2 · s^-1 ) light for 72 h to examine mechanisms of photoacclimation at the transcriptional level. High-throughput sequence data were used to design specific primers for genes involved with light harvesting and these were quantified with qPCR. The greatest significant difference in transcript abundances was observed in the psaA gene (Photosystem I P700 apoprotein), and site-type had an effect on these responses. Shade-acclimated thalli were 22-fold down-regulated at high light, whereas sun-acclimated thalli were only 5-fold down-regulated. Another gene involved with Photosystem I (petF ferredoxin) was down-regulated at high light, but only individuals from the shaded site were significantly different (4-fold). In thalli from both sites, cpeA (Phycoerythrin alpha chain) was down-regulated at high light. Although not statistically significant, patterns consistent with previous physiological and transcriptomic studies were uncovered, namely the inverse response of transcriptional activity in genes that encode phycobiliproteins. In support of previous ecophysiological studies of freshwater red algae, these data indicate significant transcriptional changes involving Photosystem I and phycobiliprotein synthesis are required to tolerate and grow at various irradiances.

Analysis of an improved Cyanophora paradoxa genome assembly

Glaucophyta are members of the Archaeplastida, the founding group of photosynthetic eukaryotes that also includes red algae (Rhodophyta), green algae, and plants (Viridiplantae). Here we present a high-quality assembly, built using long-read sequences, of the ca. 100 Mb nuclear genome of the model glaucophyte Cyanophora paradoxa. We also conducted a quick-freeze deep-etch electron microscopy (QFDEEM) analysis of C. paradoxa cells to investigate glaucophyte morphology in comparison to other organisms. Using the genome data, we generated a resolved 115-taxon eukaryotic tree of life that includes a well-supported, monophyletic Archaeplastida. Analysis of muroplast peptidoglycan (PG) ultrastructure using QFDEEM shows that PG is most dense at the cleavage-furrow. Analysis of the chlamydial contribution to glaucophytes and other Archaeplastida shows that these foreign sequences likely played a key role in anaerobic glycolysis in primordial algae to alleviate ATP starvation under night-time hypoxia. The robust genome assembly of C. paradoxa significantly advances knowledge about this model species and provides a reference for exploring the panoply of traits associated with the anciently diverged glaucophyte lineage.

Expansion of phycobilisome linker gene families in mesophilic red algae

The common ancestor of red algae (Rhodophyta) has undergone massive genome reduction, whereby 25% of the gene inventory has been lost, followed by its split into the species-poor extremophilic Cyanidiophytina and the broadly distributed mesophilic red algae. Success of the mesophile radiation is surprising given their highly reduced gene inventory. To address this latter issue, we combine an improved genome assembly from the unicellular red alga Porphyridium purpureum with a diverse collection of other algal genomes to reconstruct ancient endosymbiotic gene transfers (EGTs) and gene duplications. We find EGTs associated with the core photosynthetic machinery that may have played important roles in plastid establishment. More significant are the extensive duplications and diversification of nuclear gene families encoding phycobilisome linker proteins that stabilize light-harvesting functions. We speculate that the origin of these complex families in mesophilic red algae may have contributed to their adaptation to a diversity of light environments.

Widely distributed red algae have experienced massive genome reduction during evolution. Here, using an improved genome assembly of Porphyridium purpureum, Lee et al. show the role of endosymbiotic gene transfer in plastid evolution and the correlation between phycobilisome linker diversification and the red algal radiation.

Horizontal and endosymbiotic gene transfer in early plastid evolution

Plastids evolved from a cyanobacterium that was engulfed by a heterotrophic eukaryotic host and became a stable organelle. Some of the resulting eukaryotic algae entered into a number of secondary endosymbioses with diverse eukaryotic hosts. These events had major consequences on the evolution and diversification of life on Earth. Although almost all plastid diversity derives from a single endosymbiotic event, the analysis of nuclear genomes of plastid-bearing lineages has revealed a mosaic origin of plastid-related genes. In addition to cyanobacterial genes, plastids recruited for their functioning eukaryotic proteins encoded by the host nucleus and also bacterial proteins of noncyanobacterial origin. Therefore, plastid proteins and plastid-localised metabolic pathways evolved by tinkering and using gene toolkits from different sources. This mixed heritage seems especially complex in secondary algae containing green plastids, the acquisition of which appears to have been facilitated by many previous acquisitions of red algal genes (the 'red carpet hypothesis').

Complete chloroplast genome of red seaweedHalymenia maculata(J. Agardh, 1885)

In this study, we first assembled and characterized the complete chloroplast genome of red seaweed Halymenia maculata. It is 190,431 bp in length, with an AT content of 70.56%. In total, 122 genes were identified, and they consisted of 87 protein-coding genes, 33 tRNA genes, and two rRNA genes. The chloroplast genome of H. maculata did not show an obvious quadripartite structure. A total of five microsatellites (SSRs) were identified in the genome using MISA. Phylogenetic analysis revealed that H. maculata was first clustered with Grateloupia filicina and Grateloupia taiwanensis in a monophyletic clade that provides useful data for the phylogeny and taxonomy of Rhodymeniophycidae.

Evolution: King-Size Plastid Genomes in a New Red Algal Clade

Plastids, the photosynthetic organelles of eukaryotes, exhibit remarkably stable genome architecture. However, a recent study of microscopic red algae has found new record-sized plastid genomes with unusual architectures. These species form a new branch in the tree of red algae.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

In Silico Analysis of Relationship between Proteins from Plastid Genome of Red Alga Palmaria sp. (Japan) and Angiotensin I Converting Enzyme Inhibitory Peptides

Plastid proteins are one of the main components in red algae. In order to clarify the angiotensin I converting enzyme (ACE) inhibitory peptides from red alga Palmaria sp. (Japan), we determined the plastid genome sequence. The genome possesses 205 protein coding genes, which were classified as genetic systems, ribosomal proteins, photosystems, adenosine triphosphate (ATP) synthesis, metabolism, transport, or unknown. After comparing ACE inhibitory peptides between protein sequences and a database, photosystems (177 ACE inhibitory peptides) were found to be the major source of ACE inhibitory peptides (total of 751). Photosystems consist of phycobilisomes, photosystem I, photosystem II, cytochrome complex, and a redox system. Among them, photosystem I (53) and II (51) were the major source of ACE inhibitory peptides. We found that the amino acid sequence of apcE (14) in phycobilisomes, psaA (18) and psaB (13) in photosystem I, and psbB (11) and psbC (10) in photosystem II covered a majority of bioactive peptide sequences. These results are useful for evaluating the bioactive peptides from red algae.

Evolutionary dynamics of the chromatophore genome in three photosynthetic Paulinella species

The thecate amoeba Paulinella is a valuable model for understanding plastid organellogenesis because this lineage has independently gained plastids (termed chromatophores) of alpha-cyanobacterial provenance. Plastid primary endosymbiosis in Paulinella occurred relatively recently (90–140 million years ago, Mya), whereas the origin of the canonical Archaeplastida plastid occurred >1,500 Mya. Therefore, these two events provide independent perspectives on plastid formation on vastly different timescales. Here we generated the complete chromatophore genome sequence from P. longichromatophora (979,356 bp, GC-content = 38.8%, 915 predicted genes) and P. micropora NZ27 (977,190 bp, GC-content = 39.9%, 911 predicted genes) and compared these data to that from existing chromatophore genomes. Our analysis suggests that when a basal split occurred among photosynthetic Paulinella species ca. 60 Mya, only 35% of the ancestral orthologous gene families from the cyanobacterial endosymbiont remained in chromatophore DNA. Following major gene losses during the early stages of endosymbiosis, this process slowed down significantly, resulting in a conserved gene content across extant taxa. Chromatophore genes faced relaxed selection when compared to homologs in free-living alpha-cyanobacteria, likely reflecting the homogeneous intracellular environment of the Paulinella host. Comparison of nucleotide substitution and insertion/deletion events among different P. micropora strains demonstrates that increases in AT-content and genome reduction are ongoing and dynamic processes in chromatophore evolution.

Flip-flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae, Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein-coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip-flop organization. The flip-flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70-gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data-rich approach (i.e., organelle genome data) from broader taxon sampling.

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

BACKGROUND

The Synurophyceae is one of most important photosynthetic stramenopile algal lineages in freshwater ecosystems. They are characterized by siliceous scales covering the cell or colony surface and possess plastids of red-algal secondary or tertiary endosymbiotic origin. Despite their ecological and evolutionary significance, the relationships amongst extant Synurophyceae are unclear, as is their relationship to most other stramenopiles.

RESULTS

Here we report a comparative analysis of plastid genomes sequenced from five representative synurophycean algae. Most of these plastid genomes are highly conserved with respect to genome structure and coding capacity, with the exception of gene re-arrangements and partial duplications at the boundary of the inverted repeat and single-copy regions. Several lineage-specific gene loss/gain events and intron insertions were detected (e.g., cemA, dnaB, syfB, and trnL).

CONCLUSIONS

Unexpectedly, the cemA gene of Synurophyceae shows a strong relationship with sequences from members of the green-algal lineage, suggesting the occurrence of a lateral gene transfer event. Using a molecular clock approach based on silica fossil record data, we infer the timing of genome re-arrangement and gene gain/loss events in the plastid genomes of Synurophyceae.

Plastid Genomes from Diverse Glaucophyte Genera Reveal a Largely Conserved Gene Content and Limited Architectural Diversity

Plastid genome (ptDNA) data of Glaucophyta have been limited for many years to the genus Cyanophora. Here, we sequenced the ptDNAs of Gloeochaete wittrockiana, Cyanoptyche gloeocystis, Glaucocystis incrassata, and Glaucocystis sp. BBH. The reported sequences are the first genome-scale plastid data available for these three poorly studied glaucophyte genera. Although the Glaucophyta plastids appear morphologically “ancestral,” they actually bear derived genomes not radically different from those of red algae or viridiplants. The glaucophyte plastid coding capacity is highly conserved (112 genes shared) and the architecture of the plastid chromosomes is relatively simple. Phylogenomic analyses recovered Glaucophyta as the earliest diverging Archaeplastida lineage, but the position of viridiplants as the first branching group was not rejected by the approximately unbiased test. Pairwise distances estimated from 19 different plastid genes revealed that the highest sequence divergence between glaucophyte genera is frequently higher than distances between species of different classes within red algae or viridiplants. Gene synteny and sequence similarity in the ptDNAs of the two Glaucocystis species analyzed is conserved. However, the ptDNA of Gla. incrassata contains a 7.9-kb insertion not detected in Glaucocystis sp. BBH. The insertion contains ten open reading frames that include four coding regions similar to bacterial serine recombinases (two open reading frames), DNA primases, and peptidoglycan aminohydrolases. These three enzymes, often encoded in bacterial plasmids and bacteriophage genomes, are known to participate in the mobilization and replication of DNA mobile elements. It is therefore plausible that the insertion in Gla. incrassata ptDNA is derived from a DNA mobile element.

Organellar genomics: a useful tool to study evolutionary relationships and molecular evolution in Gracilariaceae (Rhodophyta)

Gracilariaceae has a worldwide distribution including numerous economically important species. We applied high-throughput sequencing to obtain organellar genomes (mitochondria and chloroplast) from 10 species of Gracilariaceae and, combined with published genomes, to infer phylogenies and compare genome architecture among species representing main lineages. We obtained similar topologies between chloroplast and mitochondrial genomes phylogenies. However, the chloroplast phylogeny was better resolved with full support. In this phylogeny, Melanthalia intermedia is sister to a monophyletic clade including Gracilaria and Gracilariopsis, which were both resolved as monophyletic genera. Mitochondrial and chloroplast genomes were highly conserved in gene synteny, and variation mainly occurred in regions where insertions of plasmid-derived sequences (PDS) were found. In mitochondrial genomes, PDS insertions were observed in two regions where the transcription direction changes: between the genes cob and trnL, and trnA and trnN. In chloroplast genomes, PDS insertions were in different positions, but generally found between psdD and rrs genes. Gracilariaceae is a good model system to study the impact of PDS in genome evolution due to the frequent presence of these insertions in organellar genomes. Furthermore, the bacterial leuC/leuD operon was found in chloroplast genomes of Gracilaria tenuistipitata, G. chilensis, and M. intermedia, and in extrachromosomal plasmid of G. vermiculophylla. Phylogenetic trees show two different origins of leuC/leuD: genes found in chloroplast and plasmid were placed with proteobacteria, and genes encoded in the nucleus were close to Viridiplantae and cyanobacteria.

Mitochondrial and Plastid Genomes from Coralline Red Algae Provide Insights into the Incongruent Evolutionary Histories of Organelles

Mitochondria and plastids are generally uniparentally inherited and have a conserved gene content over hundreds of millions of years, which makes them potentially useful phylogenetic markers. Organelle single gene-based trees have long been the basis for elucidating interspecies relationships that inform taxonomy. More recently, high-throughput genome sequencing has enabled the construction of massive organelle genome databases from diverse eukaryotes, and these have been used to infer species relationships in deep evolutionary time. Here, we test the idea that despite their expected utility, conflicting phylogenetic signal may exist in mitochondrial and plastid genomes from the anciently diverged coralline red algae (Rhodophyta). We generated complete organelle genome data from five coralline red algae (Lithothamnion sp., Neogoniolithon spectabile, Renouxia sp., Rhodogorgon sp., and Synarthrophyton chejuensis) for comparative analysis with existing organelle genome data from two other species (Calliarthron tuberculosum and Sporolithon durum). We find strong evidence for incongruent phylogenetic signal from both organelle genomes that may be explained by incomplete lineage sorting that has maintained anciently derived gene copies or other molecular evolutionary processes such as hybridization or gene flow during the evolutionary history of coralline red algae.

The first plastid genome of a filamentous taxon 'Bangia' sp. OUCPT-01 in the Bangiales

Red algae are important primary photosynthetic organisms. The Bangiales comprise a morphologically diverse order of red algae. Until now, complete plastid genomes of the Bangiales were only mapped for foliose species. To date, no filamentous plastomes have been published. The aim of this study was to determine and analyze the complete plastid genome of the filamentous marine species ‘Bangia’ sp. OUCPT-01. It is a circular molecule, 196,913 bps in length with a guanine-cytosine (GC) content of 33.5%. It has a quadripartite structure with two single copy regions separated by two direct non-identical repeats. It has 205 protein-coding genes, 37 tRNAs, and 6 rRNAs. Therefore, it has a high coding capacity and is highly similar to other Bangiales species in terms of content and structure. In particular, it reveals that the genera in the Bangiales have highly conserved gene content and plastome synteny. This plastome and existing data provide insights into the phylogenetic relationships among the Bangiales genera of the Rhodophyta. According to its plastid- and mitochondrial genomes, ‘Bangia 2′ is a sister group to Porphyra. However, the position of Wildemania schizophylla in the Bangiales is still controversial. Our results show that the Bangiales divergence time was ~225 million years ago.

Complete mitochondrial genomes of six species of the freshwater red algal order Batrachospermales (Rhodophyta)

Only two mitochondrial (mt) genomes had been reported in members of the red algal order Batrachospermales, which are confined to freshwater habitats. Additional mt genomes of six representative members (Batrachospermum macrosporum, Kumanoa ambigua, K. mahlacensis, Paralemanea sp., Sheathia arcuata, and Sirodotia delicatula) were sequenced aiming to gain insights on the evolution of their mt genomes from a comparative analysis with other red algal groups. Mt genomes sequenced had the following characteristics: lengths ranging between 24,864 nt and 29,785 nt, 22 to 26 protein-coding genes, G + C contents of 21.3 to 30.7%, number of tRNA of 16 to 37, non-coding DNA from 3.8% to 14.8%. Comparative analysis revealed that mt genomes in Batrachospermales are highly conserved in terms of genome size and gene content and synteny. Phylogenetic analyses based on COI nucleotide data revealed high bootstrap support only for the genera usually recovered in the phylogenetic analyses but no support for supra-generic groups. The insertion of a group II intron carrying an ORF coding for the corresponding intron maturase interrupting the COI gene was observed in Paralamenea sp. and accounted for its larger genome in comparison to the other Batrachospermales mt genomes.

Plastid genome analysis of three Nemaliophycidae red algal species suggests environmental adaptation for iron limited habitats

The red algal subclass Nemaliophycidae includes both marine and freshwater taxa that contribute to more than half of the freshwater species in Rhodophyta. Given that these taxa inhabit diverse habitats, the Nemaliophycidae is a suitable model for studying environmental adaptation. For this purpose, we characterized plastid genomes of two freshwater species, Kumanoa americana (Batrachospermales) and Thorea hispida (Thoreales), and one marine species Palmaria palmata (Palmariales). Comparative genome analysis identified seven genes (ycf34, ycf35, ycf37, ycf46, ycf91, grx, and pbsA) that were different among marine and freshwater species. Among currently available red algal plastid genomes (127), four genes (pbsA, ycf34, ycf35, ycf37) were retained in most of the marine species. Among these, the pbsA gene, known for encoding heme oxygenase, had two additional copies (HMOX1 and HMOX2) that were newly discovered and were reported from previously red algal nuclear genomes. Each type of heme oxygenase had a different evolutionary history and special modifications (e.g., plastid targeting signal peptide). Based on this observation, we suggest that the plastid-encoded pbsA contributes to the iron controlling system in iron-deprived conditions. Thus, we highlight that this functional requirement may have prevented gene loss during the long evolutionary history of red algal plastid genomes.

Organization of plastid genomes in the freshwater red algal order Batrachospermales (Rhodophyta)

Little is known about genome organization in members of the order Batrachospermales, and the infra-ordinal relationship remains unresolved. Plastid (cp) genomes of seven members of the freshwater red algal order Batrachospermales were sequenced, with the following aims: (i) to describe the characteristics of cp genomes and compare these with other red algal groups; (ii) to infer the phylogenetic relationships among these members to better understand the infra-ordinal classification. Cp genomes of Batrachospermales are large, with several cases of gene loss, they are gene-dense (high gene content for the genome size and short intergenic regions) and have highly conserved gene order. Phylogenetic analyses based on concatenated nucleotide genome data roughly supports the current taxonomic system for the order. Comparative analyses confirm data for members of the class Florideophyceae that cp genomes in Batrachospermales is highly conserved, with little variation in gene composition. However, relevant new features were revealed in our study: genome sizes in members of Batrachospermales are close to the lowest values reported for Florideophyceae; differences in cp genome size within the order are large in comparison with other orders (Ceramiales, Gelidiales, Gracilariales, Hildenbrandiales, and Nemaliales); and members of Batrachospermales have the lowest number of protein-coding genes among the Florideophyceae. In terms of gene loss, apcF, which encodes the allophycocyanin beta subunit, is absent in all sequenced taxa of Batrachospermales. We reinforce that the interordinal relationships between the freshwater orders Batrachospermales and Thoreales within the Nemaliophycidae is not well resolved due to limited taxon sampling.