Divergent copies of the large inverted repeat in the chloroplast genomes of ulvophycean green algae

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.

Chloroplast-nuclear genome interaction drives asymmetric inverted repeats in chloroplast genome of Huperzia

Huperzia plants are important natural medicine species that have attracted increased amounts of attention due to the presence of huperzine A (Hup A), which efficiently, selectively, and reversibly inhibits acetylcholinesterase activity. However, due to the complexity of the genome structure and lack of genetic information, the phylogenetic and evolutionary relationships of Huperzia species are still unclear. In this study, seven chloroplast genome structures of Huperzia were verified and comparatively analyzed, and the phylogenetic relationships between Huperzia and other related ferns were evaluated. The results showed that a 58 kb insertion in the middle of the IRb region in the seven chloroplast genomes led to large differences in sequence size, GC content and tRNA number between IRa and IRb. All seven chloroplast genomes exhibited asymmetric structural conformations in the IR regions. We used sequence alignment and screening of the organelle and nuclear genomes to identify the source of the gene. The inserted sequence had no homologous region with the mitochondrial genome, and the same sequence as another segment was found in the nuclear genome. Then, the chloroplast genomes of Huperzia were comparatively analyzed. Phylogenetic tree showed that Phlegmariurus and Huperzia were closely related and belonged to the same branch. H. javanica and Huperzia serrata are closely related. From the perspective of phylogenetic tree, their phylogenetic relationship is closely related to geographical distribution. In conclusion, this study provides insights into the interaction mechanism between chloroplast genome and nuclear genome and the phylogeny of Huperzia. In addition, by comparing the chloroplast genomes of Huperzia, we analyzed its genetic diversity and provided important molecular information for the protection of these endangered species.

The New Genus Caulinema Revealed New Insights into the Generic Relationship of the Order Ulotrichales (Ulvophyceae, Chlorophyta)

Traditionally, the order Ulotrichales comprised green algae of an unbranched, uniseriate, filamentous morphology. However, since the establishment of ultrastructural features, the circumscription of this order has dramatically changed. Some genera and species have been excluded from this order and others with different morphologies (sarcinoid, branched filaments or even parenchymatous taxa) have been included. Phylogenetic analyses have confirmed the monophyly of this order, but its differentiation from the Ulvales and Acrosiphoniales remains difficult because of the lack of synapomorphies at every level (morphology, molecular signatures). To demonstrate the difficulties of placement into genera and orders, we investigated two sarcinoid taxa with the absence of zoospore formation. SSU and ITS rDNA tree topology and the ITS-2/CBC approach revealed that both strains SAG 2661 and CCAP 312/1 belong to Ulosarcina terrestrica and the newly erected genus Caulinema, respectively. The species conception using this approach was evaluated by sequencing the plastid-coding gene tufA, a commonly used barcode marker for green algae. All three molecular markers resulted in similar topologies at the generic and species levels, which is consistent with the ITS-2/CBC approach and tufA for barcoding. The reevaluation of the ultrastructural features revealed that the presence of organic scales on the surfaces of motile cells is characteristic for the order Ulotrichales and can be used for separation from the closely related orders. As a consequence of our study, we propose the new genus Caulinema for strain CCAP 312/1.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Global exact optimisations for chloroplast structural haplotype scaffolding

BACKGROUND

Scaffolding is an intermediate stage of fragment assembly. It consists in orienting and ordering the contigs obtained by the assembly of the sequencing reads. In the general case, the problem has been largely studied with the use of distances data between the contigs. Here we focus on a dedicated scaffolding for the chloroplast genomes. As these genomes are small, circular and with few specific repeats, numerous approaches have been proposed to assemble them. However, their specificities have not been sufficiently exploited.

RESULTS

We give a new formulation for the scaffolding in the case of chloroplast genomes as a discrete optimisation problem, that we prove the decision version to be [Formula: see text]-Complete. We take advantage of the knowledge of chloroplast genomes and succeed in expressing the relationships between a few specific genomic repeats in mathematical constraints. Our approach is independent of the distances and adopts a genomic regions view, with the priority on scaffolding the repeats first. In this way, we encode the structural haplotype issue in order to retrieve several genome forms that coexist in the same chloroplast cell. To solve exactly the optimisation problem, we develop an integer linear program that we implement in Python3 package khloraascaf. We test it on synthetic data to investigate its performance behaviour and its robustness against several chosen difficulties.

CONCLUSIONS

We succeed to model biological knowledge on genomic structures to scaffold chloroplast genomes. Our results suggest that modelling genomic regions is sufficient for scaffolding repeats and is suitable for finding several solutions corresponding to several genome forms.

Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)

BACKGROUND

Neltuma pallida is a tree that grows in arid soils in northwestern Peru. As a predominant species of the Equatorial Dry Forest ecoregion, it holds significant economic and ecological value for both people and environment. Despite this, the species is severely threatened and there is a lack of genetic and genomic research, hindering the proposal of evidence-based conservation strategies.

RESULTS

In this work, we conducted the assembly, annotation, analysis and comparison of the chloroplast genome of a N. pallida specimen with those of related species. The assembled chloroplast genome has a length of 162,381 bp with a typical quadripartite structure (LSC-IRA-SSC-IRB). The calculated GC content was 35.97%. However, this is variable between regions, with a higher GC content observed in the IRs. A total of 132 genes were annotated, of which 19 were duplicates and 22 contained at least one intron in their sequence. A substantial number of repetitive sequences of different types were identified in the assembled genome, predominantly tandem repeats (> 300). In particular, 142 microsatellites (SSR) markers were identified. The phylogenetic reconstruction showed that N. pallida grouped with the other Neltuma species and with Prosopis cineraria. The analysis of sequence divergence between the chloroplast genome sequences of N. pallida, N. juliflora, P. farcta and Strombocarpa tamarugo revealed a high degree of similarity.

CONCLUSIONS

The N. pallida chloroplast genome was found to be similar to those of closely related species. With a size of 162,831 bp, it had the classical chloroplast quadripartite structure and GC content of 35.97%. Most of the 132 identified genes were protein-coding genes. Additionally, over 800 repetitive sequences were identified, including 142 SSR markers. In the phylogenetic analysis, N. pallida grouped with other Neltuma spp. and P. cineraria. Furthermore, N. pallida chloroplast was highly conserved when compared with genomes of closely related species. These findings can be of great potential for further diversity studies and genetic improvement of N. pallida.

Hydrocytium expands the phylogenetic, morphological, and genomic diversity of the poorly known green algal order Chaetopeltidales

PREMISE

Chaetopeltidales is a small, understudied order of the green algal class Chlorophyceae, that is slowly expanding with the occasional discoveries of novel algae. Here we demonstrate that hitherto unrecognized chaetopeltidaleans also exist among previously described but neglected and misclassified species.

METHODS

Strain SAG 40.91 of Characium acuminatum, shown by previous preliminary evidence to have affinities with the orders Oedogoniales, Chaetophorales, and Chaetopeltidales (together constituting the OCC clade), was investigated with light and electron microscopy to characterize its morphology and ultrastructure. Sequence assemblies of the organellar and nuclear genomes were obtained and utilized in bioinformatic and phylogenetic analyses to address the phylogenetic position of the alga and its salient genomic features.

RESULTS

The characterization of strain SAG 40.91 and a critical literature review led us to reinstate the forgotten genus Hydrocytium A.Braun 1855, with SAG 40.91 representing its type species, Hydrocytium acuminatum. Independent molecular markers converged on placing H. acuminatum as a deeply diverged lineage of the order Chaetopeltidales, formalized as the new family Hydrocytiaceae. Both chloroplast and mitochondrial genomes shared characteristics with other members of Chaetopeltidales and were bloated by repetitive sequences. Notably, the mitochondrial cox2a gene was transferred into the nuclear genome in the H. acuminatum lineage, independently of the same event in Volvocales. The nuclear genome data from H. acuminatum and from another chaetopeltidalean that was reported by others revealed endogenized viral sequences corresponding to novel members of the phylum Nucleocytoviricota.

CONCLUSIONS

The resurrected genus Hydrocytium expands the known diversity of chaetopeltidalean algae and provides the first glimpse into their virosphere.

IRplus: An Augmented Tool to Detect Inverted Repeats in Plastid Genomes

High-throughput sequencing methods have increased the accessibility of plastid genomes, which are crucial for clarifying phylogenetic relationships. Current large sequencing efforts require software tools for routine display of their distinctive quadripartite structure, which is denoted by four junction sites. By concentrating on these junctions and their close vicinity, IRscope has emerged as the standard tool for detection of this structure and creating simplified comparative graphical maps of plastid genomes. Here, we provide an augmented version (IRplus) that encompasses a novel set of functions such as integrated error detection, flexible color schemes, and an upgraded method to detect inverted repeats in genomic sequences. Spanning across the plant tree of life, IRplus allows the quick visualization of various sets of plastid genomes and features, next to smooth interoperability with other widely used annotation file formats and platforms. The IRplus can be accessed at https://irscope.shinyapps.io/IRplus/, and source codes are freely available at https://github.com/AmiryousefiLab/IRplus.

Seven newly sequenced chloroplast genomes from the order Watanabeales (Trebouxiophyceae, Chlorophyta): Phylogenetic and comparative analysis

The little-known order Watanabeales currently includes 10 genera with Chlorella-like species that reproduce by unequal-sized autospores and are predominantly solitary or terrestrial. The taxonomic scheme of Watanabeales has only been primarily inferred by short and less informative rDNA phylogenetic analysis. In the present study, seven newly sequenced genomes and one reported chloroplast genome representing the existing major branches of Watanabeales were harvested to phylogenetically reconstruct this order and to further understand its evolution. All chloroplast genomes of Watanabeales ranging from 133 to 274 kb were circular mapping and lacked a quadripartite structure. The chloroplast genome size, GC content, number of introns, and length of intergenic region proportion of the Watanabeales showed consistent trends, with Calidiella yingdensis D201 and Kalinella pachyderma 2601 having the lowest and highest values, respectively, echoing the positive correlation between organismal size and genome size. Phylogenetic analysis of Watanabeales based on 76 protein-coding genes coupled with the establishment of various complex analytical methods determined the unique robust taxonomic scheme which was incongruent with rDNA. Comparative genomic analyses revealed that the chloroplast genomes of Watanabeales accounted for numerous complex rearrangements and inversions which indicated high cryptic diversity. Substitution rate estimation indicated that the chloroplast genomes of Watanabeales were under purifying selection and similar evolutionary pressure and supported the view that genus Symbiochloris should be excluded from Watanabeales. Our results enrich the chloroplast genome resources of Watanabeales, clarify the phylogenetic status of species within this order, and provide more reference information for subsequent taxonomic and phylogenetic study.

Characterization of the chloroplast genome of Symbiochlorum hainanensis (Ulvophyceae, Chlorophyta) and its phylogenetic analysis

Symbiochlorum hainandiae S.Q. Gong & Z.Y. Li, 2018 is a unicellular green alga belonging to Ulvophyceae, Chlorophyta, and plays important roles in coral reef ecosystem. In this study, high-throughput sequencing technology is used to sequence and assemble the chloroplast genome of S. hainandiae. The complete chloroplast genome of S. hainandiae was 158, 96 bp with the GC content of 32.86%. A total of 126 genes were identified, including 98 protein-coding genes, 26 tRNA, and 2 rRNA genes. The inverted repeat region was lost in the complete chloroplast genome of S. hainandiae. The phylogenetic analysis supports that S. hainandiae is a new sister lineage to the genus Ignatius within the class Ulvophyceae.

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.

Challenging the Importance of Plastid Genome Structure Conservation: New Insights From Euglenophytes

Plastids, similar to mitochondria, are organelles of endosymbiotic origin, which retained their vestigial genomes (ptDNA). Their unique architecture, commonly referred to as the quadripartite (four-part) structure, is considered to be strictly conserved; however, the bulk of our knowledge on their variability and evolutionary transformations comes from studies of the primary plastids of green algae and land plants. To broaden our perspective, we obtained seven new ptDNA sequences from freshwater species of photosynthetic euglenids-a group that obtained secondary plastids, known to have dynamically evolving genome structure, via endosymbiosis with a green alga. Our analyses have demonstrated that the evolutionary history of euglenid plastid genome structure is exceptionally convoluted, with a patchy distribution of inverted ribosomal operon (rDNA) repeats, as well as several independent acquisitions of tandemly repeated rDNA copies. Moreover, we have shown that inverted repeats in euglenid ptDNA do not share their genome-stabilizing property documented in chlorophytes. We hypothesize that the degeneration of the quadripartite structure of euglenid plastid genomes is connected to the group II intron expansion. These findings challenge the current global paradigms of plastid genome architecture evolution and underscore the often-underestimated divergence between the functionality of shared traits in primary and complex plastid organelles.

Structural Variation of Plastomes Provides Key Insight Into the Deep Phylogeny of Ferns

Structural variation of plastid genomes (plastomes), particularly large inversions and gene losses, can provide key evidence for the deep phylogeny of plants. In this study, we investigated the structural variation of fern plastomes in a phylogenetic context. A total of 127 plastomes representing all 50 recognized families and 11 orders of ferns were sampled, making it the most comprehensive plastomic analysis of fern lineages to date. The samples included 42 novel plastomes of 15 families with a focus on Hymenophyllales and Gleicheniales. We reconstructed a well-supported phylogeny of all extant fern families, detected significant structural synapomorphies, including 9 large inversions, 7 invert repeat region (IR) boundary shifts, 10 protein-coding gene losses, 7 tRNA gene losses or anticodon changes, and 19 codon indels (insertions or deletions) across the deep phylogeny of ferns, particularly on the backbone nodes. The newly identified inversion V5, together with the newly inferred expansion of the IR boundary R5, can be identified as a synapomorphy of a clade composed of Dipteridaceae, Matoniaceae, Schizaeales, and the core leptosporangiates, while a unique inversion V4, together with an expansion of the IR boundary R4, was verified as a synapomorphy of Gleicheniaceae. This structural evidence is in support of our phylogenetic inference, thus providing key insight into the paraphyly of Gleicheniales. The inversions of V5 and V7 together filled the crucial gap regarding how the "reversed" gene orientation in the IR region characterized by most extant ferns (Schizaeales and the core leptosporangiates) evolved from the inferred ancestral type as retained in Equisetales and Osmundales. The tRNA genes trnR-ACG and trnM-CAU were assumed to be relicts of the early-divergent fern lineages but intact in most Polypodiales, particularly in eupolypods; and the loss of the tRNA genes trnR-CCG, trnV-UAC, and trnR-UCU in fern plastomes was much more prevalent than previously thought. We also identified several codon indels in protein-coding genes within the core leptosporangiates, which may be identified as synapomorphies of specific families or higher ranks. This study provides an empirical case of integrating structural and sequence information of plastomes to resolve deep phylogeny of plants.

Phylotranscriptomic insights into a Mesoproterozoic-Neoproterozoic origin and early radiation of green seaweeds (Ulvophyceae)

The Ulvophyceae, a major group of green algae, is of particular evolutionary interest because of its remarkable morphological and ecological diversity. Its phylogenetic relationships and diversification timeline, however, are still not fully resolved. In this study, using an extensive nuclear gene dataset, we apply coalescent- and concatenation-based approaches to reconstruct the phylogeny of the Ulvophyceae and to explore the sources of conflict in previous phylogenomic studies. The Ulvophyceae is recovered as a paraphyletic group, with the Bryopsidales being a sister group to the Chlorophyceae, and the remaining taxa forming a clade (Ulvophyceae sensu stricto). Molecular clock analyses with different calibration strategies emphasize the large impact of fossil calibrations, and indicate a Meso-Neoproterozoic origin of the Ulvophyceae (sensu stricto), earlier than previous estimates. The results imply that ulvophyceans may have had a profound influence on oceanic redox structures and global biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition.

“Ulvophyceae is a remarkably morphologically and ecologically diverse clade of green algae. Here, the authors reconstruct the Ulvophyceae phylogeny, showing that these algae originated earlier than expected and may have influenced biogeochemical cycles at the Mesoproterozoic-Neoproterozoic transition.”

Comparison of Auxenochlorella protothecoides and Chlorella spp. Chloroplast Genomes: Evidence for Endosymbiosis and Horizontal Virus-like Gene Transfer

Resequencing of the chloroplast genome (cpDNA) of Auxenochlorella protothecoides UTEX 25 was completed (GenBank Accession no. KC631634.1), revealing a genome size of 84,576 base pairs and 30.8% GC content, consistent with features reported for the previously sequenced A. protothecoides 0710, (GenBank Accession no. KC843975). The A. protothecoides UTEX 25 cpDNA encoded 78 predicted open reading frames, 32 tRNAs, and 4 rRNAs, making it smaller and more compact than the cpDNA genome of C. variabilis (124,579 bp) and C. vulgaris (150,613 bp). By comparison, the compact genome size of A. protothecoides was attributable primarily to a lower intergenic sequence content. The cpDNA coding regions of all known Chlorella species were found to be organized in conserved colinear blocks, with some rearrangements. The Auxenochlorella and Chlorella species genome structure and composition were similar, and of particular interest were genes influencing photosynthetic efficiency, i.e., chlorophyll synthesis and photosystem subunit I and II genes, consistent with other biofuel species of interest. Phylogenetic analysis revealed that Prototheca cutis is the closest known A. protothecoides relative, followed by members of the genus Chlorella. The cpDNA of A. protothecoides encodes 37 genes that are highly homologous to representative cyanobacteria species, including rrn16, rrn23, and psbA, corroborating a well-recognized symbiosis. Several putative coding regions were identified that shared high nucleotide sequence identity with virus-like sequences, suggestive of horizontal gene transfer. Despite these predictions, no corresponding transcripts were obtained by RT-PCR amplification, indicating they are unlikely to be expressed in the extant lineage.

Standard Candles for Dating Microbial Lineages

Molecular clock analyses are challenging for microbial phylogenies, due to a lack of fossil calibrations that can reliably provide absolute time constraints. An alternative source of temporal constraints for microbial groups is provided by the inheritance of proteins that are specific for the utilization of eukaryote-derived substrates, which have often been dispersed across the Tree of Life via horizontal gene transfer. In particular, animal, algal, and plant-derived substrates are often produced by groups with more precisely known divergence times, providing an older-bound on their availability within microbial environments. Therefore, these ages can serve as "standard candles" for dating microbial groups across the Tree of Life, expanding the reach of informative molecular clock investigations. Here, we formally develop the concept of substrate standard candles and describe how they can be propagated and applied using both microbial species trees and individual gene family phylogenies. We also provide detailed evaluations of several candidate standard candles and discuss their suitability in light of their often complex evolutionary and metabolic histories.

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.

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.

Mitochondrial and Plastid Genomes of the Monoraphid Diatom Schizostauron trachyderma

We provide for the first time the complete plastid and mitochondrial genomes of a monoraphid diatom: Schizostauron trachyderma. The mitogenome is 41,957 bp in size and displays two group II introns in the cox1 gene. The 187,029 bp plastid genome features the typical quadripartite architecture of diatom genomes. It contains a group II intron in the petB gene that overlaps the large single-copy and the inverted repeat region. There is also a group IB4 intron encoding a putative LAGLIDADG homing endonuclease in the rnl gene. The multigene phylogenies conducted provide more evidence of the proximity between S. trachyderma and fistula-bearing species of biraphid diatoms.

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

The widespread IS200/IS605 transposon family encodes diverse programmable RNA-guided endonucleases

IscB proteins are putative nucleases encoded in a distinct family of IS200/IS605 transposons and are likely ancestors of the RNA-guided endonuclease Cas9, but the functions of IscB and its interactions with any RNA remain uncharacterized. Using evolutionary analysis, RNA sequencing, and biochemical experiments, we reconstructed the evolution of CRISPR-Cas9 systems from IS200/IS605 transposons. We found that IscB uses a single noncoding RNA for RNA-guided cleavage of double-stranded DNA and can be harnessed for genome editing in human cells. We also demonstrate the RNA-guided nuclease activity of TnpB, another IS200/IS605 transposon-encoded protein and the likely ancestor of Cas12 endonucleases. This work reveals a widespread class of transposon-encoded RNA-guided nucleases, which we name OMEGA (obligate mobile element–guided activity), with strong potential for developing as biotechnologies.

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.

airpg: automatically accessing the inverted repeats of archived plastid genomes

BACKGROUND

In most flowering plants, the plastid genome exhibits a quadripartite genome structure, comprising a large and a small single copy as well as two inverted repeat regions. Thousands of plastid genomes have been sequenced and submitted to public sequence repositories in recent years. The quality of sequence annotations in many of these submissions is known to be problematic, especially regarding annotations that specify the length and location of the inverted repeats: such annotations are either missing or portray the length or location of the repeats incorrectly. However, many biological investigations employ publicly available plastid genomes at face value and implicitly assume the correctness of their sequence annotations.

RESULTS

We introduce airpg, a Python package that automatically assesses the frequency of incomplete or incorrect annotations of the inverted repeats among publicly available plastid genomes. Specifically, the tool automatically retrieves plastid genomes from NCBI Nucleotide under variable search parameters, surveys them for length and location specifications of inverted repeats, and confirms any inverted repeat annotations through self-comparisons of the genome sequences. The package also includes functionality for automatic identification and removal of duplicate genome records and accounts for taxa that genuinely lack inverted repeats. A survey of the presence of inverted repeat annotations among all plastid genomes of flowering plants submitted to NCBI Nucleotide until the end of 2020 using airpg, followed by a statistical analysis of potential associations with record metadata, highlights that release year and publication status of the genome records have a significant effect on the frequency of complete and equal-length inverted repeat annotations.

CONCLUSION

The number of plastid genomes on NCBI Nucleotide has increased dramatically in recent years, and many more genomes will likely be submitted over the next decade. airpg enables researchers to automatically access and evaluate the inverted repeats of these plastid genomes as well as their sequence annotations and, thus, contributes to increasing the reliability of publicly available plastid genomes. The software is freely available via the Python package index at http://pypi.python.org/pypi/airpg .

Phylogenomic analysis restructures the ulvophyceae

Here, we present new transcriptome sequencing data from seven species of Dasycladales (Ulvophyceae) and a phylogenomic analysis of the Chlorophyta with a particular focus on Ulvophyceae. We have focused on a broad selection of green algal groups and carefully selected genes suitable for reconstructing deep eukaryote evolutionary histories. Increasing the taxon sampling of Dasycladales restructures the Ulvophyceae by identifying Dasycladales as closely related to Scotinosphaerales and Oltmannsiellopsidales. Contrary to previous studies, we do not find support for a close relationship between Dasycladales and a group with Cladophorales and Trentepohliales. Instead, the latter group is sister to the remainder of the Ulvophyceae. Furthermore, our analyses show high and consistent statistical support for a sister relationship between Bryopsidales and Chlorophyceae in trees generated with both homogeneous and heterogeneous (heterotachy) evolutionary models. Our study provides a new framework for interpreting the evolutionary history of Ulvophyceae and the evolution of cellular morphologies.

Molecular Phylogeny of Unicellular Marine Coccoid Green Algae Revealed New Insights into the Systematics of the Ulvophyceae (Chlorophyta)

Most marine coccoid and sarcinoid green algal species have traditionally been placed within genera dominated by species from freshwater or soil habitats. For example, the genera Chlorocystis and Halochlorococcum contain exclusively marine species; however, their familial and ordinal affinities are unclear. They are characterized by a vegetative cell with lobated or reticulated chloroplast, formation of quadriflagellated zoospores and living epi- or endophytically within benthic macroalgae. They were integrated into the family Chlorochytriaceae which embraces all coccoid green algae with epi- or endophytic life phases. Later, they were excluded from the family of Chlorococcales based on studies of their life histories in culture, and transferred to their newly described order, Chlorocystidales of the Ulvophyceae. Both genera form a "Codiolum"-stage that serves as the unicellular sporophyte in their life cycles. Phylogenetic analyses of SSU and ITS rDNA sequences confirmed that these coccoid taxa belong to the Chlorocystidales, together with the sarcinoid genus Desmochloris. The biflagellated coccoid strains were members of the genus Sykidion, which represented its own order, Sykidiales, among the Ulvophyceae. Considering these results and the usage of the ITS-2/CBC approach revealed three species of Desmochloris, six of Chlorocystis, and three of Sykidion. Three new species and several new combinations were proposed.

Genome sequencing and de novo assembly of the giant unicellular alga Acetabularia acetabulum using droplet MDA

The macroscopic single-celled green alga Acetabularia acetabulum has been a model system in cell biology for more than a century. However, no genomic information is available from this species. Since the alga has a long life cycle, is difficult to grow in dense cultures, and has an estimated diploid genome size of almost 2 Gb, obtaining sufficient genomic material for genome sequencing is challenging. Here, we have attempted to overcome these challenges by amplifying genomic DNA using multiple displacement amplification (MDA) combined with microfluidics technology to distribute the amplification reactions across thousands of microscopic droplets. By amplifying and sequencing DNA from five single cells we were able to recover an estimated ~ 7–11% of the total genome, providing the first draft of the A. acetabulum genome. We highlight challenges associated with genome recovery and assembly of MDA data due to biases arising during genome amplification, and hope that our study can serve as a reference for future attempts on sequencing the genome from non-model eukaryotes.

Large Phylogenomic Data sets Reveal Deep Relationships and Trait Evolution in Chlorophyte Green Algae

The chlorophyte green algae (Chlorophyta) are species-rich ancient groups ubiquitous in various habitats with high cytological diversity, ranging from microscopic to macroscopic organisms. However, the deep phylogeny within core Chlorophyta remains unresolved, in part due to the relatively sparse taxon and gene sampling in previous studies. Here we contribute new transcriptomic data and reconstruct phylogenetic relationships of core Chlorophyta based on four large data sets up to 2,698 genes of 70 species, representing 80% of extant orders. The impacts of outgroup choice, missing data, bootstrap-support cutoffs, and model misspecification in phylogenetic inference of core Chlorophyta are examined. The species tree topologies of core Chlorophyta from different analyses are highly congruent, with strong supports at many relationships (e.g., the Bryopsidales and the Scotinosphaerales-Dasycladales clade). The monophyly of Chlorophyceae and of Trebouxiophyceae as well as the uncertain placement of Chlorodendrophyceae and Pedinophyceae corroborate results from previous studies. The reconstruction of ancestral scenarios illustrates the evolution of the freshwater-sea and microscopic–macroscopic transition in the Ulvophyceae, and the transformation of unicellular→colonial→multicellular in the chlorophyte green algae. In addition, we provided new evidence that serine is encoded by both canonical codons and noncanonical TAG code in Scotinosphaerales, and stop-to-sense codon reassignment in the Ulvophyceae has originated independently at least three times. Our robust phylogenetic framework of core Chlorophyta unveils the evolutionary history of phycoplast, cyto-morphology, and noncanonical genetic codes in chlorophyte green algae.

Haslea silbo, A Novel Cosmopolitan Species of Blue Diatoms

Specimens of a new species of blue diatoms from the genus Haslea Simonsen were discovered in geographically distant sampling sites, first in the Canary Archipelago, then North Carolina, Gulf of Naples, the Croatian South Adriatic Sea, and Turkish coast of the Eastern Mediterranean Sea. An exhaustive characterization of these specimens, using a combined morphological and genomic approach led to the conclusion that they belong to a single new to science cosmopolitan species, Haslea silbo sp. nov. A preliminary characterization of its blue pigment shows similarities to marennine produced by Haslea ostrearia, as evidenced by UV-visible spectrophotometry and Raman spectrometry. Life cycle stages including auxosporulation were also observed, providing data on the cardinal points of this species. For the two most geographically distant populations (North Carolina and East Mediterranean), complete mitochondrial and plastid genomes were sequenced. The mitogenomes of both strains share a rare atp6 pseudogene, but the number, nature, and positions of the group II introns inside its cox1 gene differ between the two populations. There are also two pairs of genes fused in single ORFs. The plastid genomes are characterized by large regions of recombination with plasmid DNA, which are in both cases located between the ycf35 and psbA genes, but whose content differs between the strains. The two sequenced strains hosts three plasmids coding for putative serine recombinase protein whose sequences are compared, and four out of six of these plasmids were highly conserved.

A gene-rich and compact chloroplast genome of the green alga Nephroselmis pyriformis (N.Carter) Ettl 1982 from the shores of Mersin (Eastern Mediterranean Sea)

We report the complete chloroplast genome of the MED1 strain of Nephroselmis pyriformis from the Eastern Mediterranean Sea. At 111,026 bp, this genome is smaller and more compact than those of Nephroselmis olivacea and Nephroselmis astigmatica, and in contrast to the latter taxa, its inverted repeat contains no complete protein-coding genes. It encodes 3 rRNAs, 33 tRNAs and 94 proteins. Maximum likelihood analysis of a concatenated set of chloroplast genes from green algae belonging to deep-diverging lineages positioned the three Nephroselmis species in a strongly supported clade in which N. pyriformis is sister to N. astigmatica.

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

Cephaleuros is often known as an algal pathogen with 19 taxonomically valid species, some of which are responsible for red rust and algal spot diseases in vascular plants. No chloroplast genomes have yet been reported in this genus, and the limited genetic information is an obstacle to understanding the evolution of this genus. In this study, we sequenced six new Trentepohliales chloroplast genomes, including four Cephaleuros and two Trentepohlia. The chloroplast genomes of Trentepohliales are large compared to most green algae, ranging from 216 to 408 kbp. They encode between 93 and 98 genes and have a GC content of 26-36%. All new chloroplast genomes were circular-mapping and lacked a quadripartite structure, in contrast to the previously sequenced Trentepohlia odorata, which does have an inverted repeat. The duplicated trnD -GTC, petD, and atpA genes in C. karstenii may be remnants of the IR region and shed light on its reduction. Chloroplast genes of Trentepohliales show elevated rates of evolution, strong rearrangement dynamics and several genes display an alternative genetic code with reassignment of the UGA/UAG codon presumably coding for arginine. Our results present the first whole chloroplast genome of the genus Cephaleuros and enrich the chloroplast genome resources of Trentepohliales.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

Comparative Plastid Genomics of Cryptomonas Species Reveals Fine-Scale Genomic Responses to Loss of Photosynthesis

Loss of photosynthesis is a recurring theme in eukaryotic evolution. In organisms that have lost the ability to photosynthesize, nonphotosynthetic plastids are retained because they play essential roles in processes other than photosynthesis. The unicellular algal genus Cryptomonas contains both photosynthetic and nonphotosynthetic members, the latter having lost the ability to photosynthesize on at least three separate occasions. To elucidate the evolutionary processes underlying the loss of photosynthesis, we sequenced the plastid genomes of two nonphotosynthetic strains, Cryptomonas sp. CCAC1634B and SAG977-2f, as well as the genome of the phototroph Cryptomonas curvata CCAP979/52. These three genome sequences were compared with the previously sequenced plastid genome of the nonphotosynthetic species Cryptomonas paramecium CCAP977/2a as well as photosynthetic members of the Cryptomonadales, including C. curvata FBCC300012D. Intraspecies comparison between the two C. curvata strains showed that although their genome structures are stable, the substitution rates of their genes are relatively high. Although most photosynthesis-related genes, such as the psa and psb gene families, were found to have disappeared from the nonphotosynthetic strains, at least ten pseudogenes are retained in SAG977-2f. Although gene order is roughly shared among the plastid genomes of photosynthetic Cryptomonadales, genome rearrangements are seen more frequently in the smaller genomes of the nonphotosynthetic strains. Intriguingly, the light-independent protochlorophyllide reductase comprising chlB, L, and N is retained in nonphotosynthetic SAG977-2f and CCAC1634B. On the other hand, whereas CCAP977/2a retains ribulose-1,5-bisphosphate carboxylase/oxygenase-related genes, including rbcL, rbcS, and cbbX, the plastid genomes of the other two nonphotosynthetic strains have lost the ribulose-1,5-bisphosphate carboxylase/oxygenase protein-coding genes.

Gene Loss and Evolution of the Plastome

Chloroplasts are unique organelles within the plant cells and are responsible for sustaining life forms on the earth due to their ability to conduct photosynthesis. Multiple functional genes within the chloroplast are responsible for a variety of metabolic processes that occur in the chloroplast. Considering its fundamental role in sustaining life on the earth, it is important to identify the level of diversity present in the chloroplast genome, what genes and genomic content have been lost, what genes have been transferred to the nuclear genome, duplication events, and the overall origin and evolution of the chloroplast genome. Our analysis of 2511 chloroplast genomes indicated that the genome size and number of coding DNA sequences (CDS) in the chloroplasts genome of algae are higher relative to other lineages. Approximately 10.31% of the examined species have lost the inverted repeats (IR) in the chloroplast genome that span across all the lineages. Genome-wide analyses revealed the loss of the Rbcl gene in parasitic and heterotrophic plants occurred approximately 56 Ma ago. PsaM, Psb30, ChlB, ChlL, ChlN, and Rpl21 were found to be characteristic signature genes of the chloroplast genome of algae, bryophytes, pteridophytes, and gymnosperms; however, none of these genes were found in the angiosperm or magnoliid lineage which appeared to have lost them approximately 203-156 Ma ago. A variety of chloroplast-encoded genes were lost across different species lineages throughout the evolutionary process. The Rpl20 gene, however, was found to be the most stable and intact gene in the chloroplast genome and was not lost in any of the analyzed species, suggesting that it is a signature gene of the plastome. Our evolutionary analysis indicated that chloroplast genomes evolved from multiple common ancestors ~1293 Ma ago and have undergone vivid recombination events across different taxonomic lineages.

Chloroplast Dual Divergent Promoter Plasmid for Heterologous Protein Expression in Tetraselmis suecica (Chlorophyceae, Chlorodendrales)

The eukaryotic green microalga Tetraselmis suecica is commonly used for aquaculture purposes because of its high stress tolerance and ease of culture in a wide spectrum of environments; they are therefore suitable candidates for biotechnology applications. To date, no data are available regarding chloroplast transformation vectors based on specific endogenous promoters and homologous targeting regions. We report on the identification of Tetraselmis suecica genes encoding the ribulose bisphosphate carboxylase/oxygenase large subunit protein, the photosystem II D1 protein and the ATP synthase CF1-beta subunit protein together with their untranslated regions (5'UTR, 3'UTR). The full-length ORFs of the putative genes with their regulatory sequences were obtained. We were also able to identify the downstream 3' end of the large subunit ribosomal RNA gene (23S) along with the 5S RNA end-to-end with the psbA gene on the complementary strand. The intergenic region between these genes appears to be a good target site for the integration of target proteins. Moreover, we identified a back-to-back promoter region among the rbcL and atpB genes. To assess the bidirectionality activities of both promoters, a dual reporter vector was constructed for Tetraselmis suecica transformation containing the cat and TurboGFP genes driven by the 5'rbcL/5'atpB divergent promoter. The vector included the 23S-5S and psbA nucleotide sequences as flanking regions. These flanking regions provided suitable insertion sites within the chloroplast genome for cassette integration via homologous recombination. Simultaneous expression of the chloramphenicol-resistant conferring gene and the gene coding for TurboGFP driven by 5'rbcL/5'atpB showed a potent natural bidirectional promoter as a reliable genetic tool.

PACVr: plastome assembly coverage visualization in R

BACKGROUND

Plastid genomes typically display a circular, quadripartite structure with two inverted repeat regions, which challenges automatic assembly procedures. The correct assembly of plastid genomes is a prerequisite for the validity of subsequent analyses on genome structure and evolution. The average coverage depth of a genome assembly is often used as an indicator of assembly quality. Visualizing coverage depth across a draft genome is a critical step, which allows users to inspect the quality of the assembly and, where applicable, identify regions of reduced assembly confidence. Despite the interplay between genome structure and assembly quality, no contemporary, user-friendly software tool can visualize the coverage depth of a plastid genome assembly while taking its quadripartite genome structure into account. A software tool is needed that fills this void.

RESULTS

We introduce 'PACVr', an R package that visualizes the coverage depth of a plastid genome assembly in relation to the circular, quadripartite structure of the genome as well as the individual plastome genes. By using a variable window approach, the tool allows visualizations on different calculation scales. It also confirms sequence equality of, as well as visualizes gene synteny between, the inverted repeat regions of the input genome. As a tool for plastid genomics, PACVr provides the functionality to identify regions of coverage depth above or below user-defined threshold values and helps to identify non-identical IR regions. To allow easy integration into bioinformatic workflows, PACVr can be invoked from a Unix shell, facilitating its use in automated quality control. We illustrate the application of PACVr on four empirical datasets and compare visualizations generated by PACVr with those of alternative software tools.

CONCLUSIONS

PACVr provides a user-friendly tool to visualize (a) the coverage depth of a plastid genome assembly on a circular, quadripartite plastome map and in relation to individual plastome genes, and (b) gene synteny across the inverted repeat regions. It contributes to optimizing plastid genome assemblies and increasing the reliability of publicly available plastome sequences. The software, example datasets, technical documentation, and a tutorial are available with the package at https://cran.r-project.org/package=PACVr.

The Complete Chloroplast Genome of Two Important Annual Clover Species, Trifolium alexandrinum and T. resupinatum: Genome Structure, Comparative Analyses and Phylogenetic Relationships with Relatives in Leguminosae

Trifolium L., which belongs to the IR lacking clade (IRLC), is one of the largest genera in the Leguminosae and contains several economically important fodder species. Here, we present whole chloroplast (cp) genome sequencing and annotation of two important annual grasses, Trifolium alexandrinum (Egyptian clover) and T. resupinatum (Persian clover). Abundant single nucleotide polymorphisms (SNPs) and insertions/deletions (In/Dels) were discovered between those two species. Global alignment of T. alexandrinum and T. resupinatum to a further thirteen Trifolium species revealed a large amount of rearrangement and repetitive events in these fifteen species. As hypothetical cp open reading frame (ORF) and RNA polymerase subunits, ycf1 and rpoC2 in the cp genomes both contain vast repetitive sequences and observed high Pi values (0.7008, 0.3982) between T. alexandrinum and T. resupinatum. Thus they could be considered as the candidate genes for phylogenetic analysis of Trifolium species. In addition, the divergence time of those IR lacking Trifolium species ranged from 84.8505 Mya to 4.7720 Mya. This study will provide insight into the evolution of Trifolium species.

Neoproterozoic origin and multiple transitions to macroscopic growth in green seaweeds

The Neoproterozoic Era records the transition from a largely bacterial to a predominantly eukaryotic phototrophic world, creating the foundation for the complex benthic ecosystems that have sustained Metazoa from the Ediacaran Period onward. This study focuses on the evolutionary origins of green seaweeds, which play an important ecological role in the benthos of modern sunlit oceans and likely played a crucial part in the evolution of early animals by structuring benthic habitats and providing novel niches. By applying a phylogenomic approach, we resolve deep relationships of the core Chlorophyta (Ulvophyceae or green seaweeds, and freshwater or terrestrial Chlorophyceae and Trebouxiophyceae) and unveil a rapid radiation of Chlorophyceae and the principal lineages of the Ulvophyceae late in the Neoproterozoic Era. Our time-calibrated tree points to an origin and early diversification of green seaweeds in the late Tonian and Cryogenian periods, an interval marked by two global glaciations with strong consequent changes in the amount of available marine benthic habitat. We hypothesize that unicellular and simple multicellular ancestors of green seaweeds survived these extreme climate events in isolated refugia, and diversified in benthic environments that became increasingly available as ice retreated. An increased supply of nutrients and biotic interactions, such as grazing pressure, likely triggered the independent evolution of macroscopic growth via different strategies, including true multicellularity, and multiple types of giant-celled forms.

The chloroplast genome of the lichen-symbiont microalga Trebouxia sp. Tr9 (Trebouxiophyceae, Chlorophyta) shows short inverted repeats with a single gene and loss of the rps4 gene, which is encoded by the nucleus

The Trebouxiophyceae is the class of Chlorophyta algae from which the highest number of chloroplast genome (cpDNA) sequences has been obtained. Several species in this class participate in symbioses with fungi to form lichens. However, no cpDNA has been obtained from any Trebouxia lichen-symbiont microalgae, which are present in approximately half of all lichens. Here, we report the sequence of the completely assembled cpDNA from Trebouxia sp. TR9 and a comparative study with other Trebouxio-phyceae. The organization of the chloroplast genome of Trebouxia sp. TR9 has certain features that are unusual in the Trebouxiophyceae and other green algae. The most remarkable characteristics are the presence of long intergenic spacers, a quadripartite structure with short inverted repeated sequences (IRs), and the loss of the rps4 gene. The presence of long intergenic spacers accounts for a larger cpDNA size in comparison to other closely related Trebouxiophyceae. The IRs, which were thought to be lost in the Trebouxiales, are distinct from most of cpDNAs since they lack the rRNA operon and uniquely includes the rbcL gene. The functional transfer of the rps4 gene to the nuclear genome has been confirmed by sequencing and examination of the gene architecture, which includes three spliceosomal introns as well as the verification of the presence of the corresponding transcript. This is the first documented transfer of the rps4 gene from the chloroplast to the nucleus among Viridiplantae. Additionally, a fairly well-resolved phylogenetic reconstruction, including Trebouxia sp. TR9 along with other Trebouxiophyceae, was obtained based on a set of conserved chloroplast genes.

Highly Reduced Plastid Genomes of the Non-photosynthetic Dictyochophyceans Pteridomonas spp. (Ochrophyta, SAR) Are Retained for tRNA-Glu-Based Organellar Heme Biosynthesis

Organisms that have lost their photosynthetic capabilities are present in a variety of eukaryotic lineages, such as plants and disparate algal groups. Most of such non-photosynthetic eukaryotes still carry plastids, as these organelles retain essential biological functions. Most non-photosynthetic plastids possess genomes with varied protein-coding contents. Such remnant plastids are known to be present in the non-photosynthetic, bacteriovorous alga Pteridomonas danica (Dictyochophyceae, Ochrophyta), which, regardless of its obligatory heterotrophic lifestyle, has been reported to retain the typically plastid-encoded gene for ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) large subunit (rbcL). The presence of rbcL without photosynthetic activity suggests that investigating the function of plastids in Pteridomonas spp. would likely bring unique insights into understanding the reductive evolution of plastids, their genomes, and plastid functions retained after the loss of photosynthesis. In this study, we demonstrate that two newly established strains of the non-photosynthetic genus Pteridomonas possess highly reduced plastid genomes lacking rbcL gene, in contrast to the previous report. Interestingly, we discovered that all plastid-encoded proteins in Pteridomonas spp. are involved only in housekeeping processes (e.g., transcription, translation and protein degradation), indicating that all metabolite synthesis pathways in their plastids are supported fully by nuclear genome-encoded proteins. Moreover, through an in-depth survey of the available transcriptomic data of another strain of the genus, we detected no candidate sequences for nuclear-encoded, plastid-directed Fe-S cluster assembly pathway proteins, suggesting complete loss of this pathway in the organelle, despite its widespread conservation in non-photosynthetic plastids. Instead, the transcriptome contains plastid-targeted components of heme biosynthesis, glycolysis, and pentose phosphate pathways. The retention of the plastid genomes in Pteridomonas spp. is not explained by the Suf-mediated constraint against loss of plastid genomes, previously proposed for Alveolates, as they lack Suf genes. Bearing all these findings in mind, we propose the hypothesis that plastid DNA is retained in Pteridomonas spp. for the purpose of providing glutamyl-tRNA, encoded by trnE gene, as a substrate for the heme biosynthesis pathway.

Complete mitogenome of the streptophyte green alga Coleochaete scutata (Coleochaetophyceae)

We have sequenced the mitogenome of Coleochaete scutata strain SAG 110.80M. This mitogenome is the largest among the streptophyte green algae examined to date. At 242,024 bp, it is 4.3-fold larger than the mitogenome of Chaetosphaeridium globosum, the only other mitogenome available for the Coleochaetophyceae. This size difference is mainly explained by differences in the abundance of introns and in the length of intergenic regions containing vestiges of coding sequences thought to be of foreign origin. With 31 group I and 26 group II introns, the C. scutata mitogenome is the most intron-rich organelle genome known among streptophyte algae.

Promising prospects of nanopore sequencing for algal hologenomics and structural variation discovery

BACKGROUND

The MinION Access Program (MAP, 2014-2016) allowed selected users to test the prospects of long nanopore reads for diverse organisms and applications through the rapid development of improving chemistries. In 2014, faced with a fragmented Illumina assembly for the chloroplast genome of the green algal holobiont Caulerpa ashmeadii, we applied to the MAP to test the prospects of nanopore reads to investigate such intricacies, as well as further explore the hologenome of this species with native and hybrid approaches.

RESULTS

The chloroplast genome could only be resolved as a circular molecule in nanopore assemblies, which also revealed structural variants (i.e. chloroplast polymorphism or heteroplasmy). Signal and Illumina polishing of nanopore-assembled organelle genomes (chloroplast and mitochondrion) reflected the importance of coverage on final quality and current limitations. In hybrid assembly, our modest nanopore data sets showed encouraging results to improve assembly length, contiguity, repeat content, and binning of the larger nuclear and bacterial genomes. Profiling of the holobiont with nanopore or Illumina data unveiled a dominant Rhodospirillaceae (Alphaproteobacteria) species among six putative endosymbionts. While very fragmented, the cumulative hybrid assembly length of C. ashmeadii's nuclear genome reached 24.4 Mbp, including 2.1 Mbp in repeat, ranging closely with GenomeScope's estimate (> 26.3 Mbp, including 4.8 Mbp in repeat).

CONCLUSION

Our findings relying on a very modest number of nanopore R9 reads as compared to current output with newer chemistries demonstrate the promising prospects of the technology for the assembly and profiling of an algal hologenome and resolution of structural variation. The discovery of polymorphic 'chlorotypes' in C. ashmeadii, most likely mediated by homing endonucleases and/or retrohoming by reverse transcriptases, represents the first report of chloroplast heteroplasmy in the siphonous green algae. Improving contiguity of C. ashmeadii's nuclear and bacterial genomes will require deeper nanopore sequencing to greatly increase the coverage of these larger genomic compartments.

Complete chloroplast genome of the tiny marine diatom Nanofrustulum shiloi (Bacillariophyta) from the Adriatic Sea

We report the chloroplast genome sequence of Nanofrustulum shiloi, a tiny araphid pennate diatom collected from the Adriatic Sea. The 160,994-bp N. shiloi genome displays a quadripartite structure and its gene repertoire resembles those of other diatom chloroplast genomes. Besides the genes located in the inverted repeat, psbY is duplicated. A gene-poor region in the large single-copy region contains multiple ORFs sharing sequence similarities with plasmids and chloroplast ORFs found in other diatom species. The genome features a single intron, a group II intron in petB. Phylogenomic analysis identified N. shiloi at a basal position within the araphid 2 clade.

De Novo Assembly Discovered Novel Structures in Genome of Plastids and Revealed Divergent Inverted Repeats in Mammillaria (Cactaceae, Caryophyllales)

The complete sequence of chloroplast genome (cpDNA) has been documented for single large columnar species of Cactaceae, lacking inverted repeats (IRs). We sequenced cpDNA for seven species of the short-globose cacti of Mammillaria and de novo assembly revealed three novel structures in land plants. These structures have a large single copy (LSC) that is 2.5 to 10 times larger than the small single copy (SSC), and two IRs that contain strong differences in length and gene composition. Structure 1 is distinguished by short IRs of <1 kb composed by rpl23-trnI-CAU-ycf2; with a total length of 110,189 bp and 113 genes. In structure 2, each IR is approximately 7.2 kb and is composed of 11 genes and one Intergenic Spacer-(psbK-trnQ)-trnQ-UUG-rps16-trnK-UUU-matK-trnK-UUU-psbA-trnH-GUG-rpl2-rpl23-trnI-CAU-ycf2; with a total size of 116,175 bp and 120 genes. Structure 3 has divergent IRs of approximately 14.1 kb, where IRA is composed of 20 genes: psbA-trnH-GUG-rpl23-trnI-CAU-ycf2-ndhB-rps7-rps12-trnV-GAC-rrn16-ycf68-trnI-GAU-trnA-AGC-rrn23-rrn4.5-rrn5-trnR-ACG-trnN-GUU-ndhF-rpl32; and IRB is identical to the IRA, but lacks rpl23. This structure has 131 genes and, by pseudogenization, it is shown to have the shortest cpDNA, of just 107,343 bp. Our findings show that Mammillaria bears an unusual structural diversity of cpDNA, which supports the elucidation of the evolutionary processes involved in cacti lineages.

Tracing the Evolution of the Plastome and Mitogenome in the Chloropicophyceae Uncovered Convergent tRNA Gene Losses and a Variant Plastid Genetic Code

The tiny green algae belonging to the Chloropicophyceae play a key role in marine phytoplankton communities; this newly erected class of prasinophytes comprises two genera (Chloropicon and Chloroparvula) containing each several species. We sequenced the plastomes and mitogenomes of eight Chloropicon and five Chloroparvula species to better delineate the phylogenetic affinities of these taxa and to infer the suite of changes that their organelle genomes sustained during evolution. The relationships resolved in organelle-based phylogenomic trees were essentially congruent with previously reported rRNA trees, and similar evolutionary trends but distinct dynamics were identified for the plastome and mitogenome. Although the plastome sustained considerable changes in gene content and order at the time the two genera split, subsequently it remained stable and maintained a very small size. The mitogenome, however, was remodeled more gradually and showed more fluctuation in size, mainly as a result of expansions/contractions of intergenic regions. Remarkably, the plastome and mitogenome lost a common set of three tRNA genes, with the trnI(cau) and trnL(uaa) losses being accompanied with important variations in codon usage. Unexpectedly, despite the disappearance of trnI(cau) from the plastome in the Chloroparvula lineage, AUA codons (the codons recognized by this gene product) were detected in certain plastid genes. By comparing the sequences of plastid protein-coding genes from chloropicophycean and phylogenetically diverse chlorophyte algae with those of the corresponding predicted proteins, we discovered that the AUA codon was reassigned from isoleucine to methionine in Chloroparvula. This noncanonical genetic code has not previously been uncovered in plastids.

The chloroplast and mitochondrial genomes of the green algae Pediastrum duplex isolated from Central Georgia (USA)

A Pediastrum duplex (Chlorophyta) strain was isolated from a freshwater system in Milledgeville, Georgia and its chloroplast and mitochondrial genomes sequenced. The chloroplast genome was 199,241 bp with 136 genes and the mitochondrial 40,756 bp with 40 genes, both were circular. Comparison of the ‘Milledgeville’ plastome to other P. duplex isolates revealed a nearly identical sequence identity to archived genes and genomic fragments from the strain UTEX1364 which was isolated from Lake Machovo in 1962. These sequences provide chloroplast and mitochondrial genomes from a wild P. duplex isolate and provide more organelle genomes for a genus with cryptic phylogenetic relationships.

Dynamic evolution of mitochondrial genomes in Trebouxiophyceae, including the first completely assembled mtDNA from a lichen-symbiont microalga (Trebouxia sp. TR9)

Trebouxiophyceae (Chlorophyta) is a species-rich class of green algae with a remarkable morphological and ecological diversity. Currently, there are a few completely sequenced mitochondrial genomes (mtDNA) from diverse Trebouxiophyceae but none from lichen symbionts. Here, we report the mitochondrial genome sequence of Trebouxia sp. TR9 as the first complete mtDNA sequence available for a lichen-symbiont microalga. A comparative study of the mitochondrial genome of Trebouxia sp. TR9 with other chlorophytes showed important organizational changes, even between closely related taxa. The most remarkable change is the enlargement of the genome in certain Trebouxiophyceae, which is principally due to larger intergenic spacers and seems to be related to a high number of large tandem repeats. Another noticeable change is the presence of a relatively large number of group II introns interrupting a variety of tRNA genes in a single group of Trebouxiophyceae, which includes Trebouxiales and Prasiolales. In addition, a fairly well-resolved phylogeny of Trebouxiophyceae, along with other Chlorophyta lineages, was obtained based on a set of seven well-conserved mitochondrial genes.

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.