Six Newly Sequenced Chloroplast Genomes From Trentepohliales: The Inflated Genomes, Alternative Genetic Code and Dynamic Evolution

Six Newly Sequenced Chloroplast Genomes From Quadriflagellate Chlamydomonadales (Chlorophyceae): Phylogeny and Comparative Genome Analyses

Several quadriflagellate genera were revealed as the members of deep branches in Chlamydomonadales. However, the phylogenetic relationships among these quadriflagellate genera remained unresolved, and little information is known about the chloroplast genome structure for the chlamydomonadalean early-diverging lineages due to the limited data. In this study, we conducted phylogenetic and comparative genomic analyses with 6 newly sequenced quadriflagellate chlamydomonadalean chloroplast genomes. Four phylogenetic inferences based on different datasets recovered the robust topology, with Staurocarteria-Hafniomonas as the earliest-diverging lineage, followed by Corbierea within Chlamydomonadales and Spermatozopsis included in Sphaeropleales. The amino acid dataset combined with the site-heterogeneous model received the highest support for key nodes and may better fit the inferences of the deep relationships in Chlamydomonadales. Moreover, Sp. similis chloroplast genome is also more structurally similar to its close relatives than to other quadriflagellate chlamydomonadaleans. Phylogeny and genome structure features both indicated the taxonomic position of Spermatozopsis should be reconsidered. The loss of large inverted repeats (IRs) was first reported in chlamydomonadaleans (Co. pseudopalmata), and may occur at least 4 times in Chlamydomonadales. Comparative genome analyses demonstrated the highly divergent large IRs and a high level of rearrangements across the entire genome. IR expansions/contractions and inversions contribute to changes in gene content and gene order in this region. This study provides a foundation for future research on the phylogenetic relationships as well as chloroplast genomic features and evolution of the entire Chlamydomonadales.

Snow alga Sanguina aurantia as revealed through de novo genome assembly and annotation

To thrive on melting alpine and polar snow, some Chlorophytes produce an abundance of astaxanthin, causing red blooms, often dominated by genus Sanguina. The red cells have not been cultured, but we recently grew a green biciliate conspecific with Sanguina aurantia from a sample of watermelon snow. This culture provided source material for Oxford Nanopore Technology and Illumina sequencing. Our assembly pipeline exemplifies the value of a hybrid long- and short-read approach for the complexities of working with a culture grown from a field sample. Using bioinformatic tools, we separated assembled contigs into 2 genomic pools based on a difference in GC content (57.5 and 55.1%). We present the data as 2 assemblies of S. aurantia variants but explore other possibilities. High-throughput chromatin conformation capture analysis (Hi-C sequencing) was used to scaffold the assemblies into a 96-Mb genome designated as "A" and a 102-Mb genome designated as "B." Both assemblies are highly contiguous: genome A consists of 38 scaffolds with an N50 of 5.4 Mb, while genome B has 50 scaffolds with an N50 of 6.4 Mb. RNA sequencing was used to improve gene annotation.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

Comparative Analysis of Chloroplast Genomes in Cephaleuros and Its Related Genus (Trentepohlia): Insights into Adaptive Evolution

Cephaleuros species are well-known as plant pathogens that cause red rust or algae spot diseases in many economically cultivated plants that grow in shady and humid environments. Despite their prevalence, the adaptive evolution of these pathogens remains poorly understood. We sequenced and characterized three Cephaleuros (Cephaleuros lagerheimii, Cephaleuros diffusus, and Cephaleuros virescens) chloroplast genomes, and compared them with seven previously reported chloroplast genomes. The chloroplast sequences of C. lagerheimii, C. diffusus, and C. virescens were 480,613 bp, 383,846 bp, and 472,444 bp in length, respectively. These chloroplast genomes encoded 94 genes, including 27 tRNA genes, 3 rRNA genes, and 64 protein-coding genes. Comparative analysis uncovered that the variation in genome size was principally due to the length of intergenic spacer sequences, followed by introns. Furthermore, several highly variable regions (trnY-GTA, trnL-TAG, petA, psbT, trnD-GTC, trnL-TAA, ccsA, petG, psaA, psaB, rps11, rps2, and rps14) were identified. Codon bias analysis revealed that the codon usage pattern of Cephaleuros is predominantly shaped by natural selection. Additionally, six chloroplast protein-coding genes (atpF, chlN, psaA, psaB, psbA, and rbcL) were determined to be under positive selection, suggesting they may play a vital roles in the adaptation of Cephaleuros to low-light intensity habitats.

Novel insights into chloroplast genome evolution in the green macroalgal genus Ulva (Ulvophyceae, Chlorophyta)

To understand the evolutionary driving forces of chloroplast (or plastid) genomes (plastomes) in the green macroalgal genus Ulva (Ulvophyceae, Chlorophyta), in this study, we sequenced and constructed seven complete chloroplast genomes from five Ulva species, and conducted comparative genomic analysis of Ulva plastomes in Ulvophyceae. Ulva plastome evolution reflects the strong selection pressure driving the compactness of genome organization and the decrease of overall GC composition. The overall plastome sequences including canonical genes, introns, derived foreign sequences and non-coding regions show a synergetic decrease in GC content at varying degrees. Fast degeneration of plastome sequences including non-core genes (minD and trnR3), derived foreign sequences, and noncoding spacer regions was accompanied by the marked decrease of their GC composition. Plastome introns preferentially resided in conserved housekeeping genes with high GC content and long length, as might be related to high GC content of target site sequences recognized by intron-encoded proteins (IEPs), and to more target sites contained by long GC-rich genes. Many foreign DNA sequences integrated into different intergenic regions contain some homologous specific orfs with high similarity, indicating that they could have been derived from the same origin. The invasion of foreign sequences seems to be an important driving force for plastome rearrangement in these IR-lacking Ulva cpDNAs. Gene partitioning pattern has changed and distribution range of gene clusters has expanded after the loss of IR, indicating that genome rearrangement was more extensive and more frequent in Ulva plastomes, which was markedly different from that in IR-containing ulvophycean plastomes. These new insights greatly enhance our understanding of plastome evolution in ecologically important Ulva seaweeds.

Chloroplast Genome Traits Correlate With Organismal Complexity and Ecological Traits in Chlorophyta

A positive relationship between cell size and chloroplast genome size within chloroplast-bearing protists has been hypothesized in the past and shown in some case studies, but other factors influencing chloroplast genome size during the evolution of chlorophyte algae have been less studied. We study chloroplast genome size and GC content as a function of habitats and cell size of chlorophyte algae. The chloroplast genome size of green algae in freshwater, marine and terrestrial habitats was differed significantly, with terrestrial algae having larger chloroplast genome sizes in general. The most important contributor to these enlarged genomes in terrestrial species was the length of intergenic regions. There was no clear difference in the GC content of chloroplast genomes from the three habitats categories. Functional morphological categories also showed differences in chloroplast genome size, with filamentous algae having substantially larger genomes than other forms of algae, and foliose algae had lower GC content than other groups. Chloroplast genome size showed no significant differences among the classes Ulvophyceae, Trebouxiophyceae, and Chlorophyceae, but the GC content of Chlorophyceae chloroplast genomes was significantly lower than that of Ulvophyceae and Trebouxiophyceae. There was a certain positive relationship between chloroplast genome size and cell size for the Chlorophyta as a whole and within each of three major classes. Our data also confirmed previous reports that ancestral quadripartite architecture had been lost many times independently in Chlorophyta. Finally, the comparison of the phenotype of chlorophytes algae harboring plastids uncovered that most of the investigated Chlorophyta algae housed a single plastid per cell.