The exceptionally large chloroplast genome of the green alga Floydiella terrestris illuminates the evolutionary history of the Chlorophyceae

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.

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.

Complete Chloroplast Genomes of Pterodon emarginatus Vogel and Pterodon pubescens Benth: Comparative and Phylogenetic Analyses

Background

The species Pterodon emarginatus and P. pubescens, popularly known as white sucupira or faveira, are native to the Cerrado biome and have the potential for medicinal use and reforestation. They are sister species with evolutionary proximity.

Objective

Considering that the chloroplast genome exhibits a conserved structure and genes, the analysis of its sequences can contribute to the understanding of evolutionary, phylogenetic, and diversity issues.

Methods

The chloroplast genomes of P. emarginatus and P. pubescens were sequenced on the Illumina MiSeq platform. The genomes were assembled based on the de novo strategy. We performed the annotation of the genes and the repetitive regions of the genomes. The nucleotide diversity and phylogenetic relationships were analyzed using the gene sequences of these species and others of the Leguminosae family, whose genomes are available in databases.

Results

The complete chloroplast genome of P. emarginatus is 159,877 bp, and that of P. pubescens is 159,873 bp. The genomes of both species have circular and quadripartite structures. A total of 127 genes were predicted in both species, including 110 single-copy genes and 17 duplicated genes in the inverted regions. 141 microsatellite regions were identified in P. emarginatus and 140 in P. pubescens. The nucleotide diversity estimates of the gene regions in twenty-one species of the Leguminosae family were 0.062 in LSC, 0.086 in SSC, and 0.036 in IR. The phylogenetic analysis demonstrated the proximity between the genera Pterodon and Dipteryx, both from the clade Dipterygeae. Ten pairs of primers with potential for the development of molecular markers were designed.

Conclusion

The genetic information obtained on the chloroplast genomes of P. emarginatus and P. pubescens presented here reinforces the similarity and evolutionary proximity between these species, with a similarity percentage of 99.8%.

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.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

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.

Mitochondrial genomes in the iconic reindeer lichens: Architecture, variation, and synteny across multiple evolutionary scales

Variation in mitochondrial genome composition across intraspecific, interspecific, and higher taxonomic scales has been little studied in lichen obligate symbioses. Cladonia is one of the most diverse and ecologically important lichen genera, with over 500 species representing an array of unique morphologies and chemical profiles. Here, we assess mitochondrial genome diversity and variation in this flagship genus, with focused sampling of two clades of the "true" reindeer lichens, Cladonia subgenus Cladina, and additional genomes from nine outgroup taxa. We describe composition and architecture at the gene and the genome scale, examining patterns in organellar genome size in larger taxonomic groups in Ascomycota. Mitochondrial genomes of Cladonia, Pilophorus, and Stereocaulon were consistently larger than those of Lepraria and contained more introns, suggesting a selective pressure in asexual morphology in Lepraria driving it toward genomic simplification. Collectively, lichen mitochondrial genomes were larger than most other fungal life strategies, reaffirming the notion that coevolutionary streamlining does not correlate to genome size reductions. Genomes from Cladonia ravenelii and Stereocaulon pileatum exhibited ATP9 duplication, bearing paralogs that may still be functional. Homing endonuclease genes (HEGs), though scarce in Lepraria, were diverse and abundant in Cladonia, exhibiting variable evolutionary histories that were sometimes independent of the mitochondrial evolutionary history. Intraspecific HEG diversity was also high, with C. rangiferina especially bearing a range of HEGs with one unique to the species. This study reveals a rich history of events that have transformed mitochondrial genomes of Cladonia and related genera, allowing future study alongside a wealth of assembled genomes.

The Chloroplast Genome of the Lichen Photobiont Trebouxiophyceae sp. DW1 and Its Phylogenetic Implications

Lichens are symbiotic associations of algae and fungi. The genetic mechanism of the symbiosis of lichens and the influence of symbiosis on the size and composition of the genomes of symbiotic algae have always been intriguing scientific questions explored by lichenologists. However, there were limited data on lichen genomes. Therefore, we isolated and purified a lichen symbiotic alga to obtain a single strain (Trebouxiophyceae sp. DW1), and then obtained its chloroplast genome information by next-generation sequencing (NGS). The chloroplast genome is 129,447 bp in length, and the GC content is 35.2%. Repetitive sequences with the length of 30–35 bp account for 1.27% of the total chloroplast genome. The simple sequence repeats are all mononucleotide repeats. Codon usage analysis showed that the genome tended to use codon ending in A/U. By comparing the length of different regions of Trebouxiophyceae genomes, we found that the changes in the length of exons, introns, and intergenic sequences affect the size of genomes. Trebouxiophyceae had an unstable chloroplast genome structure, with IRs repeatedly losing during evolution. Phylogenetic analysis showed that Trebouxiophyceae is paraphyletic, and Trebouxiophyceae sp. DW1 is sister to the clade of Koliella longiseta and Pabia signiensis.

Dynamic plastid and mitochondrial genomes in Chaetopeltidales (Chlorophyceae) and characterization of a new chlorophyte taxon

PREMISE

Chaetopeltidales is a poorly characterized order in the Chlorophyceae, with only two plastid and no mitochondrial genomes published. Here we describe a new taxon in Chaetopeltidales, Gormaniella terricola gen. et sp. nov. and characterize both of its organellar genomes.

METHODS

Gormaniella terricola was inadvertently isolated from a surface-sterilized hornwort thallus. Light microscopy was used to characterize its vegetative morphology. Organellar genomes were assembled, annotated, and analyzed using a variety of software packages.

RESULTS

The mitochondrial genome (66,927 bp) represents the first complete mitochondrial genome published for Chaetopeltidales. The chloroplast genome, measuring 428,981 bp, is one of the largest plastid genomes published to date and shares this large size and an incredible number of short, dispersed repeats with the other sequenced chloroplast genomes in Chaetopeltidales. Despite these shared features, the chloroplast genomes of Chaetopeltidales appear to be highly rearranged when compared to one another, with numerous inversions, translocations, and duplications, suggesting a particularly dynamic chloroplast genome. Both the chloroplast and mitochondrial genomes of G. terricola contain a number of mobile group I and group II introns, which appear to have invaded separately. Three of the introns within the mitochondrial genome encode homing endonucleases that are phylogenetically nested within those found in fungi, rather than algae, suggesting a possible case of horizontal gene transfer.

CONCLUSIONS

These results help to shed light on a poorly understood group of algae and their unusual organellar genomes, raising additional questions about the unique patterns of genome evolution within Chaetopeltidales.

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.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

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.

Complete Chloroplast Genome Sequence of Erigeron breviscapus and Characterization of Chloroplast Regulatory Elements

Erigeron breviscapus is a famous medicinal plant. However, the limited chloroplast genome information of E. breviscapus, especially for the chloroplast DNA sequence resources, has hindered the study of E. breviscapus chloroplast genome transformation. Here, the complete chloroplast (cp) genome of E. breviscapus was reported. This genome was 152,164bp in length, included 37.2% GC content and was structurally arranged into two 24,699bp inverted repeats (IRs) and two single-copy areas. The sizes of the large single-copy region and the small single-copy region were 84,657 and 18,109bp, respectively. The E. breviscapus cp genome consisted of 127 coding genes, including 83 protein coding genes, 36 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. For those genes, 95 genes were single copy genes and 16 genes were duplicated in two inverted regions with seven tRNAs, four rRNAs, and five protein coding genes. Then, genomic DNA of E. breviscapus was used as a template, and the endogenous 5' and 3' flanking sequences of the trnI gene and trnA gene were selected as homologous recombinant fragments in vector construction and cloned through PCR. The endogenous 5' flanking sequences of the psbA gene and rrn16S gene, the endogenous 3' flanking sequences of the psbA gene, rbcL gene, and rps16 gene and one sequence element from the psbN-psbH chloroplast operon were cloned, and certain chloroplast regulatory elements were identified. Two homologous recombination fragments and all of these elements were constructed into the cloning vector pBluescript SK (+) to yield a series of chloroplast expression vectors, which harbored the reporter gene EGFP and the selectable marker aadA gene. After identification, the chloroplast expression vectors were transformed into Escherichia coli and the function of predicted regulatory elements was confirmed by a spectinomycin resistance test and fluorescence intensity measurement. The results indicated that aadA gene and EGFP gene were efficiently expressed under the regulation of predicted regulatory elements and the chloroplast expression vector had been successfully constructed, thereby providing a solid foundation for establishing subsequent E. breviscapus chloroplast transformation system and genetic improvement of E. breviscapus.

Identification of polycistronic transcriptional units and non-canonical introns in green algal chloroplasts based on long-read RNA sequencing data

BACKGROUND

Chloroplasts are important semi-autonomous organelles in plants and algae. Unlike higher plants, the chloroplast genomes of green algal linage have distinct features both in organization and expression. Despite the architecture of chloroplast genome having been extensively studied in higher plants and several model species of algae, little is known about the transcriptional features of green algal chloroplast-encoded genes.

RESULTS

Based on full-length cDNA (Iso-Seq) sequencing, we identified widely co-transcribed polycistronic transcriptional units (PTUs) in the green alga Caulerpa lentillifera. In addition to clusters of genes from the same pathway, we identified a series of PTUs of up to nine genes whose function in the plastid is not understood. The RNA data further allowed us to confirm widespread expression of fragmented genes and conserved open reading frames, which are both important features in green algal chloroplast genomes. In addition, a newly fragmented gene specific to C. lentillifera was discovered, which may represent a recent gene fragmentation event in the chloroplast genome. With the newly annotated exon-intron boundary information, gene structural annotation was greatly improved across the siphonous green algae lineages. Our data also revealed a type of non-canonical Group II introns, with a deviant secondary structure and intronic ORFs lacking known splicing or mobility domains. These widespread introns have conserved positions in their genes and are excised precisely despite lacking clear consensus intron boundaries.

CONCLUSION

Our study fills important knowledge gaps in chloroplast genome organization and transcription in green algae, and provides new insights into expression of polycistronic transcripts, freestanding ORFs and fragmented genes in algal chloroplast genomes. Moreover, we revealed an unusual type of Group II intron with distinct features and conserved positions in Bryopsidales. Our data represents interesting additions to knowledge of chloroplast intron structure and highlights clusters of uncharacterized genes that probably play important roles in plastids.

The Chloroplast Trans-Splicing RNA-Protein Supercomplex from the Green Alga Chlamydomonas reinhardtii

In eukaryotes, RNA trans-splicing is a significant RNA modification process for the end-to-end ligation of exons from separately transcribed primary transcripts to generate mature mRNA. So far, three different categories of RNA trans-splicing have been found in organisms within a diverse range. Here, we review trans-splicing of discontinuous group II introns, which occurs in chloroplasts and mitochondria of lower eukaryotes and plants. We discuss the origin of intronic sequences and the evolutionary relationship between chloroplast ribonucleoprotein complexes and the nuclear spliceosome. Finally, we focus on the ribonucleoprotein supercomplex involved in trans-splicing of chloroplast group II introns from the green alga Chlamydomonas reinhardtii. This complex has been well characterized genetically and biochemically, resulting in a detailed picture of the chloroplast ribonucleoprotein supercomplex. This information contributes substantially to our understanding of the function of RNA-processing machineries and might provide a blueprint for other splicing complexes involved in trans- as well as cis-splicing of organellar intron RNAs.

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.

Taxonomic scheme of the order Chaetophorales (Chlorophyceae, Chlorophyta) based on chloroplast genomes

Background

Order Chaetophorales currently includes six families, namely Schizomeridaceae, Aphanochaetaceae, Barrancaceae, Uronemataceae, Fritschiellaceae, and Chaetophoraceae. The phylogenetic relationships of Chaetophorales have been inferred primarily based on short and less informative rDNA sequences. This study aimed to phylogenetically reconstruct order Chaetophorales and determine the taxonomic scheme, and to further understand the evolution of order Chaetophorales.

Results

In the present study, seven complete and five fragmentary chloroplast genomes were harvested. Phylogenomic and comparative genomic analysis were performed to determine the taxonomic scheme within Chaetophorales. Consequently, Oedogoniales was found to be a sister to a clade linking Chaetophorales and Chaetopeltidales. Schizomeriaceae, and Aphanochaetaceae clustered into a well-resolved basal clade in Chaetophorales, inconsistent with the results of phylogenetic analysis based on rDNA sequences. Comparative genomic analyses revealed that the chloroplast genomes of Schizomeriaceae and Aphanochaetaceae were highly conserved and homologous, highlighting the closest relationship in this order. Germination types of zoospores precisely correlated with the phylogenetic relationships.

Conclusions

chloroplast genome structure analyses, synteny analyses, and zoospore germination analyses were concurrent with phylogenetic analyses based on the chloroplast genome, and all of them robustly determined the unique taxonomic scheme of Chaetophorales and the relationships of Oedogoniales, Chaetophorales, and Chaetopeltidales.

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.

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.

Characterization of the Chloroplast Genome of Trentepohlia odorata (Trentepohliales, Chlorophyta), and Discussion of its Taxonomy

Trentepohliales is an aerial order of Chlorophyta with approximately 80 species distributed mainly in tropical and subtropical regions. The taxonomy of this genus is quite difficult and presents a challenge for many phycologists. Although plentiful molecular data is available, most of the sequences are not identified at the species level. In the present study, we described a new specimen with detailed morphological data and identified it as Trentepohlia odorata. A phylogenetic analysis showed T. odorata as a novel lineage in Trentepohliales. T. odorata has the closest relationship with T. annulata, which is expected since sporangia of both species are without stalk cell and with dorsal pore. Species with such morphological characteristics may represent deep lineages in Trentepohliales. Although an increasing number of chloroplast genomes of Ulvophyceae have been reported in recent years, the whole plastome of Trentepohliales has not yet been reported. Thus, the chloroplast genome of Trentepohlia odorata was reported in the present study. The whole plastome was 399,372 bp in length, with 63 predicted protein-coding genes, 31 tRNAs, and 3 rRNAs. Additionally, we annotated 95 free-standing open reading frames, of which seven were annotated with plastid origins, 16 with eukaryotic genome origins, and 33 with bacterial genome origins. Four rpo genes (rpoA, rpoB, rpoC1, and rpoC2) were annotated within ORF clusters. These four genes were fragmented into several (partial) ORFs by in-frame stop codons. Additionally, we detected a frame shift mutation in the rpoB gene. The phylogenetic analysis supported that Trentepohliales clustered with Dasycladales and nested into the BDT clade (Bryopsidales, Dasycladales and Trentepohliales). Our results present the first whole chloroplast genome of a species of Trentepohliales and provided new data for understanding the evolution of the chloroplast genome in Ulvophyceae.

Large Diversity of Nonstandard Genes and Dynamic Evolution of Chloroplast Genomes in Siphonous Green Algae (Bryopsidales, Chlorophyta)

Chloroplast genomes have undergone tremendous alterations through the evolutionary history of the green algae (Chloroplastida). This study focuses on the evolution of chloroplast genomes in the siphonous green algae (order Bryopsidales). We present five new chloroplast genomes, which along with existing sequences, yield a data set representing all but one families of the order. Using comparative phylogenetic methods, we investigated the evolutionary dynamics of genomic features in the order. Our results show extensive variation in chloroplast genome architecture and intron content. Variation in genome size is accounted for by the amount of intergenic space and freestanding open reading frames that do not show significant homology to standard plastid genes. We show the diversity of these nonstandard genes based on their conserved protein domains, which are often associated with mobile functions (reverse transcriptase/intron maturase, integrases, phage- or plasmid-DNA primases, transposases, integrases, ligases). Investigation of the introns showed proliferation of group II introns in the early evolution of the order and their subsequent loss in the core Halimedineae, possibly through RT-mediated intron loss.

Exploring the Limits and Causes of Plastid Genome Expansion in Volvocine Green Algae

Plastid genomes are not normally celebrated for being large. But researchers are steadily uncovering algal lineages with big and, in rare cases, enormous plastid DNAs (ptDNAs), such as volvocine green algae. Plastome sequencing of five different volvocine species has revealed some of the largest, most repeat-dense plastomes on record, including that of Volvox carteri (∼525 kb). Volvocine algae have also been used as models for testing leading hypotheses on organelle genome evolution (e.g., the mutational hazard hypothesis), and it has been suggested that ptDNA inflation within this group might be a consequence of low mutation rates and/or the transition from a unicellular to multicellular existence. Here, we further our understanding of plastome size variation in the volvocine line by examining the ptDNA sequences of the colonial species Yamagishiella unicocca and Eudorina sp. NIES-3984 and the multicellular Volvox africanus, which are phylogenetically situated between species with known ptDNA sizes. Although V. africanus is closely related and similar in multicellular organization to V. carteri, its ptDNA was much less inflated than that of V. carteri. Synonymous- and noncoding-site nucleotide substitution rate analyses of these two Volvox ptDNAs suggest that there are drastically different plastid mutation rates operating in the coding versus intergenic regions, supporting the idea that error-prone DNA repair in repeat-rich intergenic spacers is contributing to genome expansion. Our results reinforce the idea that the volvocine line harbors extremes in plastome size but ultimately shed doubt on some of the previously proposed hypotheses for ptDNA inflation within the lineage.

Next-Generation Sequencing of Haematococcus lacustris Reveals an Extremely Large 1.35-Megabase Chloroplast Genome

Haematococcus lacustris is an industrially relevant microalga that is used for the production of the carotenoid astaxanthin. Here, we report the use of PacBio long-read sequencing to assemble the chloroplast genome of H. lacustris strain UTEX:2505. At 1.35 Mb, this is the largest assembled chloroplast of any plant or alga known to date.

Organellar phylogenomics inform systematics in the green algal family Hydrodictyaceae (Chlorophyceae) and provide clues to the complex evolutionary history of plastid genomes in the green algal tree of life

PREMISE OF THE STUDY

Phylogenomic analyses across the green algae are resolving relationships at the class, order, and family levels and highlighting dynamic patterns of evolution in organellar genomes. Here we present a within-family phylogenomic study to resolve genera and species relationships in the family Hydrodictyaceae (Chlorophyceae), for which poor resolution in previous phylogenetic studies, along with divergent morphological traits, have precluded taxonomic revisions.

METHODS

Complete plastome sequences and mitochondrial protein-coding gene sequences were acquired from representatives of the Hydrodictyaceae using next-generation sequencing methods. Plastomes were characterized, and gene order and content were compared with plastomes spanning the Sphaeropleales. Single-gene and concatenated-gene phylogenetic analyses of plastid and mitochondrial genes were performed.

KEY RESULTS

The Hydrodictyaceae contain the largest sphaeroplealean plastomes thus far fully sequenced. Conservation of plastome gene order within Hydrodictyaceae is striking compared with more dynamic patterns revealed across Sphaeropleales. Phylogenetic analyses resolve Hydrodictyon sister to a monophyletic Pediastrum, though the morphologically distinct P. angulosum and P. duplex continue to be polyphyletic. Analyses of plastid data supported the neochloridacean genus Chlorotetraëdron as sister to Hydrodictyaceae, while conflicting signal was found in the mitochondrial data.

CONCLUSIONS

A phylogenomic approach resolved within-family relationships not obtainable with previous phylogenetic analyses. Denser taxon sampling across Sphaeropleales is necessary to capture patterns in plastome evolution, and further taxa and studies are needed to fully resolve the sister lineage to Hydrodictyaceae and polyphyly of Pediastrum angulosum and P. duplex.

The Plastid Genome in Cladophorales Green Algae Is Encoded by Hairpin Chromosomes

Virtually all plastid (chloroplast) genomes are circular double-stranded DNA molecules, typically between 100 and 200 kb in size and encoding circa 80-250 genes. Exceptions to this universal plastid genome architecture are very few and include the dinoflagellates, where genes are located on DNA minicircles. Here we report on the highly deviant chloroplast genome of Cladophorales green algae, which is entirely fragmented into hairpin chromosomes. Short- and long-read high-throughput sequencing of DNA and RNA demonstrated that the chloroplast genes of Boodlea composita are encoded on 1- to 7-kb DNA contigs with an exceptionally high GC content, each containing a long inverted repeat with one or two protein-coding genes and conserved non-coding regions putatively involved in replication and/or expression. We propose that these contigs correspond to linear single-stranded DNA molecules that fold onto themselves to form hairpin chromosomes. The Boodlea chloroplast genes are highly divergent from their corresponding orthologs, and display an alternative genetic code. The origin of this highly deviant chloroplast genome most likely occurred before the emergence of the Cladophorales, and coincided with an elevated transfer of chloroplast genes to the nucleus. A chloroplast genome that is composed only of linear DNA molecules is unprecedented among eukaryotes, and highlights unexpected variation in plastid genome architecture.

Strategies for complete plastid genome sequencing

Plastid sequencing is an essential tool in the study of plant evolution. This high-copy organelle is one of the most technically accessible regions of the genome, and its sequence conservation makes it a valuable region for comparative genome evolution, phylogenetic analysis and population studies. Here, we discuss recent innovations and approaches for de novo plastid assembly that harness genomic tools. We focus on technical developments including low-cost sequence library preparation approaches for genome skimming, enrichment via hybrid baits and methylation-sensitive capture, sequence platforms with higher read outputs and longer read lengths, and automated tools for assembly. These developments allow for a much more streamlined assembly than via conventional short-range PCR. Although newer methods make complete plastid sequencing possible for any land plant or green alga, there are still challenges for producing finished plastomes particularly from herbarium material or from structurally divergent plastids such as those of parasitic plants.

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C-169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.

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

The chloroplast genomes of many algae and almost all land plants carry two identical copies of a large inverted repeat (IR) sequence that can pair for flip-flop recombination and undergo expansion/contraction. Although the IR has been lost multiple times during the evolution of the green algae, the underlying mechanisms are still largely unknown. A recent comparison of IR-lacking and IR-containing chloroplast genomes of chlorophytes from the Ulvophyceae (Ulotrichales) suggested that differential elimination of genes from the IR copies might lead to IR loss. To gain deeper insights into the evolutionary history of the chloroplast genome in the Ulvophyceae, we analyzed the genomes of Ignatius tetrasporus and Pseudocharacium americanum (Ignatiales, an order not previously sampled), Dangemannia microcystis (Oltmannsiellopsidales), Pseudoneochloris marina (Ulvales) and also Chamaetrichon capsulatum and Trichosarcina mucosa (Ulotrichales). Our comparison of these six chloroplast genomes with those previously reported for nine ulvophyceans revealed unsuspected variability. All newly examined genomes feature an IR, but remarkably, the copies of the IR present in the Ignatiales, Pseudoneochloris, and Chamaetrichon diverge in sequence, with the tRNA genes from the rRNA operon missing in one IR copy. The implications of this unprecedented finding for the mechanism of IR loss and flip-flop recombination are discussed.

Comparative DNA sequence analyses of Pyramimonas parkeae (Prasinophyceae) chloroplast genomes

Prasinophytes form a paraphyletic assemblage of early diverging green algae, which have the potential to reveal the traits of the last common ancestor of the main two green lineages: (i) chlorophyte algae and (ii) streptophyte algae. Understanding the genetic composition of prasinophyte algae is fundamental to understanding the diversification and evolutionary processes that may have occurred in both green lineages. In this study, we sequenced the chloroplast genome of Pyramimonas parkeae NIES254 and compared it with that of P. parkeae CCMP726, the only other fully sequenced P. parkeae chloroplast genome. The results revealed that P. parkeae chloroplast genomes are surprisingly variable. The chloroplast genome of NIES254 was larger than that of CCMP726 by 3,204 bp, the NIES254 large single copy was 288 bp longer, the small single copy was 5,088 bp longer, and the IR was 1,086 bp shorter than that of CCMP726. Similarity values of the two strains were almost zero in four large hot spot regions. Finally, the strains differed in copy number for three protein-coding genes: ycf20, psaC, and ndhE. Phylogenetic analyses using 16S and 18S rDNA and rbcL sequences resolved a clade consisting of these two P. parkeae strains and a clade consisting of these plus other Pyramimonas isolates. These results are consistent with past studies indicating that prasinophyte chloroplast genomes display a higher level of variation than is commonly found among land plants. Consequently, prasinophyte chloroplast genomes may be less useful for inferring the early history of Viridiplantae than has been the case for land plant diversification.

The Plastid Genome of Polytoma uvella Is the Largest Known among Colorless Algae and Plants and Reflects Contrasting Evolutionary Paths to Nonphotosynthetic Lifestyles

The loss of photosynthesis is frequently associated with parasitic or pathogenic lifestyles, but it also can occur in free-living, plastid-bearing lineages. A common consequence of becoming nonphotosynthetic is the reduction in size and gene content of the plastid genome. In exceptional circumstances, it can even result in the complete loss of the plastid DNA (ptDNA) and its associated gene expression system, as reported recently in several lineages, including the nonphotosynthetic green algal genus Polytomella Closely related to Polytomella is the polyphyletic genus Polytoma, the members of which lost photosynthesis independently of Polytomella Species from both genera are free-living organisms that contain nonphotosynthetic plastids, but unlike Polytomella, Polytoma members have retained a genome in their colorless plastid. Here, we present the plastid genome of Polytoma uvella: to our knowledge, the first report of ptDNA from a nonphotosynthetic chlamydomonadalean alga. The P. uvella ptDNA contains 25 protein-coding genes, most of which are related to gene expression and none are connected to photosynthesis. However, despite its reduced coding capacity, the P. uvella ptDNA is inflated with short repeats and is tens of kilobases larger than the ptDNAs of its closest known photosynthetic relatives, Chlamydomonas leiostraca and Chlamydomonas applanata In fact, at approximately 230 kb, the ptDNA of P. uvella represents the largest plastid genome currently reported from a nonphotosynthetic alga or plant. Overall, the P. uvella and Polytomella plastid genomes reveal two very different evolutionary paths following the loss of photosynthesis: expansion and complete deletion, respectively. We hypothesize that recombination-based DNA-repair mechanisms are at least partially responsible for the different evolutionary outcomes observed in such closely related nonphotosynthetic algae.

Organelle_PBA, a pipeline for assembling chloroplast and mitochondrial genomes from PacBio DNA sequencing data

Background

The development of long-read sequencing technologies, such as single-molecule real-time (SMRT) sequencing by PacBio, has produced a revolution in the sequencing of small genomes. Sequencing organelle genomes using PacBio long-read data is a cost effective, straightforward approach. Nevertheless, the availability of simple-to-use software to perform the assembly from raw reads is limited at present.

Results

We present Organelle-PBA, a Perl program designed specifically for the assembly of chloroplast and mitochondrial genomes. For chloroplast genomes, the program selects the chloroplast reads from a whole genome sequencing pool, maps the reads to a reference sequence from a closely related species, and then performs read correction and de novo assembly using Sprai. Organelle-PBA completes the assembly process with the additional step of scaffolding by SSPACE-LongRead. The program then detects the chloroplast inverted repeats and reassembles and re-orients the assembly based on the organelle origin of the reference. We have evaluated the performance of the software using PacBio reads from different species, read coverage, and reference genomes. Finally, we present the assembly of two novel chloroplast genomes from the species Picea glauca (Pinaceae) and Sinningia speciosa (Gesneriaceae).

Conclusion

Organelle-PBA is an easy-to-use Perl-based software pipeline that was written specifically to assemble mitochondrial and chloroplast genomes from whole genome PacBio reads. The program is available at https://github.com/aubombarely/Organelle_PBA.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3412-9) contains supplementary material, which is available to authorized users.

Does Cell Size Impact Chloroplast Genome Size?

There is a strong positive relationship between nuclear genome size and cell size across the eukaryotic domain, but the cause and effect of this relationship is unclear. A positive coupling of cell size and DNA content has also been recorded for various bacteria, suggesting that, with some exceptions, this association might be universal throughout the tree of life. However, the link between cell size and genome size has not yet been thoroughly explored with respect to chloroplasts, or organelles as a whole, largely because of a lack data on cell morphology and organelle DNA content. Here, I speculate about a potential positive scaling of cell size and chloroplast genome size within different plastid-bearing protists, including ulvophyte, prasinophyte, and trebouxiophyte green algae. I provide examples in which large and small chloroplast DNAs occur alongside large and small cell sizes, respectively, as well as examples where this trend does not hold. Ultimately, I argue that a relationship between cellular architecture and organelle genome architecture is worth exploring, and encourage researchers to keep an open mind on this front.

Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae)

Background

The chloroplast genome sustained extensive changes in architecture during the evolution of the Chlorophyceae, a morphologically and ecologically diverse class of green algae belonging to the Chlorophyta; however, the forces driving these changes are poorly understood. The five orders recognized in the Chlorophyceae form two major clades: the CS clade consisting of the Chlamydomonadales and Sphaeropleales, and the OCC clade consisting of the Oedogoniales, Chaetophorales, and Chaetopeltidales. In the OCC clade, considerable variations in chloroplast DNA (cpDNA) structure, size, gene order, and intron content have been observed. The large inverted repeat (IR), an ancestral feature characteristic of most green plants, is present in Oedogonium cardiacum (Oedogoniales) but is lacking in the examined members of the Chaetophorales and Chaetopeltidales. Remarkably, the Oedogonium 35.5-kb IR houses genes that were putatively acquired through horizontal DNA transfer. To better understand the dynamics of chloroplast genome evolution in the Oedogoniales, we analyzed the cpDNA of a second representative of this order, Oedocladium carolinianum.

Methods

The Oedocladium cpDNA was sequenced and annotated. The evolutionary distances separating Oedocladium and Oedogonium cpDNAs and two other pairs of chlorophycean cpDNAs were estimated using a 61-gene data set. Phylogenetic analysis of an alignment of group IIA introns from members of the OCC clade was performed. Secondary structures and insertion sites of oedogonialean group IIA introns were analyzed.

Results

The 204,438-bp Oedocladium genome is 7.9 kb larger than the Oedogonium genome, but its repertoire of conserved genes is remarkably similar and gene order differs by only one reversal. Although the 23.7-kb IR is missing the putative foreign genes found in Oedogonium, it contains sequences coding for a putative phage or bacterial DNA primase and a hypothetical protein. Intergenic sequences are 1.5-fold longer and dispersed repeats are more abundant, but a smaller fraction of the Oedocladium genome is occupied by introns. Six additional group II introns are present, five of which lack ORFs and carry highly similar sequences to that of the ORF-less IIA intron shared with Oedogonium. Secondary structure analysis of the group IIA introns disclosed marked differences in the exon-binding sites; however, each intron showed perfect or nearly perfect base pairing interactions with its target site.

Discussion

Our results suggest that chloroplast genes rearrange more slowly in the Oedogoniales than in the Chaetophorales and raise questions as to what was the nature of the foreign coding sequences in the IR of the common ancestor of the Oedogoniales. They provide the first evidence for intragenomic proliferation of group IIA introns in the Viridiplantae, revealing that intron spread in the Oedocladium lineage likely occurred by retrohoming after sequence divergence of the exon-binding sites.

Chloroplast DNA Structural Variation, Phylogeny, and Age of Divergence among Diploid Cotton Species

The cotton genus (Gossypium spp.) contains 8 monophyletic diploid genome groups (A, B, C, D, E, F, G, K) and a single allotetraploid clade (AD). To gain insight into the phylogeny of Gossypium and molecular evolution of the chloroplast genome in this group, we performed a comparative analysis of 19 Gossypium chloroplast genomes, six reported here for the first time. Nucleotide distance in non-coding regions was about three times that of coding regions. As expected, distances were smaller within than among genome groups. Phylogenetic topologies based on nucleotide and indel data support for the resolution of the 8 genome groups into 6 clades. Phylogenetic analysis of indel distribution among the 19 genomes demonstrates contrasting evolutionary dynamics in different clades, with a parallel genome downsizing in two genome groups and a biased accumulation of insertions in the clade containing the cultivated cottons leading to large (for Gossypium) chloroplast genomes. Divergence time estimates derived from the cpDNA sequence suggest that the major diploid clades had diverged approximately 10 to 11 million years ago. The complete nucleotide sequences of 6 cpDNA genomes are provided, offering a resource for cytonuclear studies in Gossypium.

Nothing in Evolution Makes Sense Except in the Light of Genomics: Read-Write Genome Evolution as an Active Biological Process

The 21st century genomics-based analysis of evolutionary variation reveals a number of novel features impossible to predict when Dobzhansky and other evolutionary biologists formulated the neo-Darwinian Modern Synthesis in the middle of the last century. These include three distinct realms of cell evolution; symbiogenetic fusions forming eukaryotic cells with multiple genome compartments; horizontal organelle, virus and DNA transfers; functional organization of proteins as systems of interacting domains subject to rapid evolution by exon shuffling and exonization; distributed genome networks integrated by mobile repetitive regulatory signals; and regulation of multicellular development by non-coding lncRNAs containing repetitive sequence components. Rather than single gene traits, all phenotypes involve coordinated activity by multiple interacting cell molecules. Genomes contain abundant and functional repetitive components in addition to the unique coding sequences envisaged in the early days of molecular biology. Combinatorial coding, plus the biochemical abilities cells possess to rearrange DNA molecules, constitute a powerful toolbox for adaptive genome rewriting. That is, cells possess "Read-Write Genomes" they alter by numerous biochemical processes capable of rapidly restructuring cellular DNA molecules. Rather than viewing genome evolution as a series of accidental modifications, we can now study it as a complex biological process of active self-modification.

Comparative Chloroplast Genome Analyses of Streptophyte Green Algae Uncover Major Structural Alterations in the Klebsormidiophyceae, Coleochaetophyceae and Zygnematophyceae

The Streptophyta comprises all land plants and six main lineages of freshwater green algae: Mesostigmatophyceae, Chlorokybophyceae, Klebsormidiophyceae, Charophyceae, Coleochaetophyceae and Zygnematophyceae. Previous comparisons of the chloroplast genome from nine streptophyte algae (including four zygnematophyceans) revealed that, although land plant chloroplast DNAs (cpDNAs) inherited most of their highly conserved structural features from green algal ancestors, considerable cpDNA changes took place during the evolution of the Zygnematophyceae, the sister group of land plants. To gain deeper insights into the evolutionary dynamics of the chloroplast genome in streptophyte algae, we sequenced the cpDNAs of nine additional taxa: two klebsormidiophyceans (Entransia fimbriata and Klebsormidium sp. SAG 51.86), one coleocheatophycean (Coleochaete scutata) and six zygnematophyceans (Cylindrocystis brebissonii, Netrium digitus, Roya obtusa, Spirogyra maxima, Cosmarium botrytis and Closterium baillyanum). Our comparative analyses of these genomes with their streptophyte algal counterparts indicate that the large inverted repeat (IR) encoding the rDNA operon experienced loss or expansion/contraction in all three sampled classes and that genes were extensively shuffled in both the Klebsormidiophyceae and Zygnematophyceae. The klebsormidiophycean genomes boast greatly expanded IRs, with the Entransia 60,590-bp IR being the largest known among green algae. The 206,025-bp Entransia cpDNA, which is one of the largest genome among streptophytes, encodes 118 standard genes, i.e., four additional genes compared to its Klebsormidium flaccidum homolog. We inferred that seven of the 21 group II introns usually found in land plants were already present in the common ancestor of the Klebsormidiophyceae and its sister lineages. At 107,236 bp and with 117 standard genes, the Coleochaete IR-less genome is both the smallest and most compact among the streptophyte algal cpDNAs analyzed thus far; it lacks eight genes relative to its Chaetosphaeridium globosum homolog, four of which represent unique events in the evolutionary scenario of gene losses we reconstructed for streptophyte algae. The 10 compared zygnematophycean cpDNAs display tremendous variations at all levels, except gene content. During zygnematophycean evolution, the IR disappeared a minimum of five times, the rDNA operon was broken at four distinct sites, group II introns were lost on at least 43 occasions, and putative foreign genes, mainly of phage/viral origin, were gained.

Hiding in plain sight: Koshicola spirodelophila gen. et sp. nov. (Chaetopeltidales, Chlorophyceae), a novel green alga associated with the aquatic angiosperm Spirodela polyrhiza

PREMISE OF THE STUDY

Discovery and morphological characterization of a novel epiphytic aquatic green alga increases our understanding of Chaetopeltidales, a poorly known order in Chlorophyceae. Chloroplast genomic data from this taxon reveals an unusual architecture previously unknown in green algae.

METHODS

Using light and electron microscopy, we characterized the morphology and ultrastructure of a novel taxon of green algae. Bayesian phylogenetic analyses of nuclear and plastid genes were used to test the hypothesized membership of this taxon in order Chaetopeltidales. With next-generation sequence data, we assembled the plastid genome of this novel taxon and compared its gene content and architecture to that of related species to further investigate plastid genome traits.

KEY RESULTS

The morphology and ultrastructure of this alga are consistent with placement in Chaetopeltidales (Chlorophyceae), but a distinct trait combination supports recognition of this alga as a new genus and species-Koshicola spirodelophila gen. et sp. nov. Its placement in the phylogeny as a descendant of a deep division in the Chaetopeltidales is supported by analysis of molecular data sets. The chloroplast genome is among the largest reported in green algae and the genes are distributed on three large (rather than a single) chromosome, in contrast to other studied green algae.

CONCLUSIONS

The discovery of Koshicola spirodelophila gen. et sp. nov. highlights the importance of investigating even commonplace habitats to explore new microalgal diversity. This work expands our understanding of the morphological and chloroplast genomic features of green algae, and in particular those of the poorly studied Chaetopeltidales.

Distinctive Architecture of the Chloroplast Genome in the Chlorodendrophycean Green Algae Scherffelia dubia and Tetraselmis sp. CCMP 881

The Chlorodendrophyceae is a small class of green algae belonging to the core Chlorophyta, an assemblage that also comprises the Pedinophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Here we describe for the first time the chloroplast genomes of chlorodendrophycean algae (Scherffelia dubia, 137,161 bp; Tetraselmis sp. CCMP 881, 100,264 bp). Characterized by a very small single-copy (SSC) region devoid of any gene and an unusually large inverted repeat (IR), the quadripartite structures of the Scherffelia and Tetraselmis genomes are unique among all core chlorophytes examined thus far. The lack of genes in the SSC region is offset by the rich and atypical gene complement of the IR, which includes genes from the SSC and large single-copy regions of prasinophyte and streptophyte chloroplast genomes having retained an ancestral quadripartite structure. Remarkably, seven of the atypical IR-encoded genes have also been observed in the IRs of pedinophycean and trebouxiophycean chloroplast genomes, suggesting that they were already present in the IR of the common ancestor of all core chlorophytes. Considering that the relationships among the main lineages of the core Chlorophyta are still unresolved, we evaluated the impact of including the Chlorodendrophyceae in chloroplast phylogenomic analyses. The trees we inferred using data sets of 79 and 108 genes from 71 chlorophytes indicate that the Chlorodendrophyceae is a deep-diverging lineage of the core Chlorophyta, although the placement of this class relative to the Pedinophyceae remains ambiguous. Interestingly, some of our phylogenomic trees together with our comparative analysis of gene order data support the monophyly of the Trebouxiophyceae, thus offering further evidence that the previously observed affiliation between the Chlorellales and Pedinophyceae is the result of systematic errors in phylogenetic reconstruction.

Increasing genomic diversity and evidence of constrained lifestyle evolution due to insertion sequences in Aeromonas salmonicida

BACKGROUND

Aeromonads make up a group of Gram-negative bacteria that includes human and fish pathogens. The Aeromonas salmonicida species has the peculiarity of including five known subspecies. However, few studies of the genomes of A. salmonicida subspecies have been reported to date.

RESULTS

We sequenced the genomes of additional A. salmonicida isolates, including three from India, using next-generation sequencing in order to gain a better understanding of the genomic and phylogenetic links between A. salmonicida subspecies. Their relative phylogenetic positions were confirmed by a core genome phylogeny based on 1645 gene sequences. The Indian isolates, which formed a sub-group together with A. salmonicida subsp. pectinolytica, were able to grow at either at 18 °C and 37 °C, unlike the A. salmonicida psychrophilic isolates that did not grow at 37 °C. Amino acid frequencies, GC content, tRNA composition, loss and gain of genes during evolution, pseudogenes as well as genes under positive selection and the mobilome were studied to explain this intraspecies dichotomy.

CONCLUSION

Insertion sequences appeared to be an important driving force that locked the psychrophilic strains into their particular lifestyle in order to conserve their genomic integrity. This observation, based on comparative genomics, is in agreement with previous results showing that insertion sequence mobility induced by heat in A. salmonicida subspecies causes genomic plasticity, resulting in a deleterious effect on the virulence of the bacterium. We provide a proof-of-concept that selfish DNAs play a major role in the evolution of bacterial species by modeling genomes.

Chloroplast phylogenomic analysis of chlorophyte green algae identifies a novel lineage sister to the Sphaeropleales (Chlorophyceae)

Background

The class Chlorophyceae (Chlorophyta) includes morphologically and ecologically diverse green algae. Most of the documented species belong to the clade formed by the Chlamydomonadales (also called Volvocales) and Sphaeropleales. Although studies based on the nuclear 18S rRNA gene or a few combined genes have shed light on the diversity and phylogenetic structure of the Chlamydomonadales, the positions of many of the monophyletic groups identified remain uncertain. Here, we used a chloroplast phylogenomic approach to delineate the relationships among these lineages.

Results

To generate the analyzed amino acid and nucleotide data sets, we sequenced the chloroplast DNAs (cpDNAs) of 24 chlorophycean taxa; these included representatives from 16 of the 21 primary clades previously recognized in the Chlamydomonadales, two taxa from a coccoid lineage (Jenufa) that was suspected to be sister to the Golenkiniaceae, and two sphaeroplealeans. Using Bayesian and/or maximum likelihood inference methods, we analyzed an amino acid data set that was assembled from 69 cpDNA-encoded proteins of 73 core chlorophyte (including 33 chlorophyceans), as well as two nucleotide data sets that were generated from the 69 genes coding for these proteins and 29 RNA-coding genes. The protein and gene phylogenies were congruent and robustly resolved the branching order of most of the investigated lineages. Within the Chlamydomonadales, 22 taxa formed an assemblage of five major clades/lineages. The earliest-diverging clade displayed Hafniomonas laevis and the Crucicarteria, and was followed by the Radicarteria and then by the Chloromonadinia. The latter lineage was sister to two superclades, one consisting of the Oogamochlamydinia and Reinhardtinia and the other of the Caudivolvoxa and Xenovolvoxa. To our surprise, the Jenufa species and the two spine-bearing green algae belonging to the Golenkinia and Treubaria genera were recovered in a highly supported monophyletic group that also included three taxa representing distinct families of the Sphaeropleales (Bracteacoccaceae, Mychonastaceae, and Scenedesmaceae).

Conclusions

Our phylogenomic study advances our knowledge regarding the circumscription and internal structure of the Chlamydomonadales, suggesting that a previously unrecognized lineage is sister to the Sphaeropleales. In addition, it offers new insights into the flagellar structures of the founding members of both the Chlamydomonadales and Sphaeropleales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0544-5) contains supplementary material, which is available to authorized users.

Dynamic Evolution of the Chloroplast Genome in the Green Algal Classes Pedinophyceae and Trebouxiophyceae

Previous studies of trebouxiophycean chloroplast genomes revealed little information regarding the evolutionary dynamics of this genome because taxon sampling was too sparse and the relationships between the sampled taxa were unknown. We recently sequenced the chloroplast genomes of 27 trebouxiophycean and 2 pedinophycean green algae to resolve the relationships among the main lineages recognized for the Trebouxiophyceae. These taxa and the previously sampled members of the Pedinophyceae and Trebouxiophyceae are included in the comparative chloroplast genome analysis we report here. The 38 genomes examined display considerable variability at all levels, except gene content. Our results highlight the high propensity of the rDNA-containing large inverted repeat (IR) to vary in size, gene content and gene order as well as the repeated losses it experienced during trebouxiophycean evolution. Of the seven predicted IR losses, one event demarcates a superclade of 11 taxa representing 5 late-diverging lineages. IR expansions/contractions account not only for changes in gene content in this region but also for changes in gene order and gene duplications. Inversions also led to gene rearrangements within the IR, including the reversal or disruption of the rDNA operon in some lineages. Most of the 20 IR-less genomes are more rearranged compared with their IR-containing homologs and tend to show an accelerated rate of sequence evolution. In the IR-less superclade, several ancestral operons were disrupted, a few genes were fragmented, and a subgroup of taxa features a G+C-biased nucleotide composition. Our analyses also unveiled putative cases of gene acquisitions through horizontal transfer.

The Organellar Genomes of Chromera and Vitrella, the Phototrophic Relatives of Apicomplexan Parasites

Apicomplexa are known to contain greatly reduced organellar genomes. Their mitochondrial genome carries only three protein-coding genes, and their plastid genome is reduced to a 35-kb-long circle. The discovery of coral-endosymbiotic algae Chromera velia and Vitrella brassicaformis, which share a common ancestry with Apicomplexa, provided an opportunity to study possibly ancestral forms of organellar genomes, a unique glimpse into the evolutionary history of apicomplexan parasites. The structurally similar mitochondrial genomes of Chromera and Vitrella differ in gene content, which is reflected in the composition of their respiratory chains. Thus, Chromera lacks respiratory complexes I and III, whereas Vitrella and apicomplexan parasites are missing only complex I. Plastid genomes differ substantially between these algae, particularly in structure: The Chromera plastid genome is a linear, 120-kb molecule with large and divergent genes, whereas the plastid genome of Vitrella is a highly compact circle that is only 85 kb long but nonetheless contains more genes than that of Chromera. It appears that organellar genomes have already been reduced in free-living phototrophic ancestors of apicomplexan parasites, and such reduction is not associated with parasitism.

The complete chloroplast and mitochondrial genomes of the green macroalga Ulva sp. UNA00071828 (Ulvophyceae, Chlorophyta)

Sequencing mitochondrial and chloroplast genomes has become an integral part in understanding the genomic machinery and the phylogenetic histories of green algae. Previously, only three chloroplast genomes (Oltmannsiellopsis viridis, Pseudendoclonium akinetum, and Bryopsis hypnoides) and two mitochondrial genomes (O. viridis and P. akinetum) from the class Ulvophyceae have been published. Here, we present the first chloroplast and mitochondrial genomes from the ecologically and economically important marine, green algal genus Ulva. The chloroplast genome of Ulva sp. was 99,983 bp in a circular-mapping molecule that lacked inverted repeats, and thus far, was the smallest ulvophycean plastid genome. This cpDNA was a highly compact, AT-rich genome that contained a total of 102 identified genes (71 protein-coding genes, 28 tRNA genes, and three ribosomal RNA genes). Additionally, five introns were annotated in four genes: atpA (1), petB (1), psbB (2), and rrl (1). The circular-mapping mitochondrial genome of Ulva sp. was 73,493 bp and follows the expanded pattern also seen in other ulvophyceans and trebouxiophyceans. The Ulva sp. mtDNA contained 29 protein-coding genes, 25 tRNA genes, and two rRNA genes for a total of 56 identifiable genes. Ten introns were annotated in this mtDNA: cox1 (4), atp1 (1), nad3 (1), nad5 (1), and rrs (3). Double-cut-and-join (DCJ) values showed that organellar genomes across Chlorophyta are highly rearranged, in contrast to the highly conserved organellar genomes of the red algae (Rhodophyta). A phylogenomic investigation of 51 plastid protein-coding genes showed that Ulvophyceae is not monophyletic, and also placed Oltmannsiellopsis (Oltmannsiellopsidales) and Tetraselmis (Chlorodendrophyceae) closely to Ulva (Ulvales) and Pseudendoclonium (Ulothrichales).

Mitochondrial and plastid genome architecture: Reoccurring themes, but significant differences at the extremes

Mitochondrial and plastid genomes show a wide array of architectures, varying immensely in size, structure, and content. Some organelle DNAs have even developed elaborate eccentricities, such as scrambled coding regions, nonstandard genetic codes, and convoluted modes of posttranscriptional modification and editing. Here, we compare and contrast the breadth of genomic complexity between mitochondrial and plastid chromosomes. Both organelle genomes have independently evolved many of the same features and taken on similar genomic embellishments, often within the same species or lineage. This trend is most likely because the nuclear-encoded proteins mediating these processes eventually leak from one organelle into the other, leading to a high likelihood of processes appearing in both compartments in parallel. However, the complexity and intensity of genomic embellishments are consistently more pronounced for mitochondria than for plastids, even when they are found in both compartments. We explore the evolutionary forces responsible for these patterns and argue that organelle DNA repair processes, mutation rates, and population genetic landscapes are all important factors leading to the observed convergence and divergence in organelle genome architecture.

The chloroplast genomes of Bryopsis plumosa and Tydemania expeditiones (Bryopsidales, Chlorophyta): compact genomes and genes of bacterial origin

Background

Species of Bryopsidales form ecologically important components of seaweed communities worldwide. These siphonous macroalgae are composed of a single giant tubular cell containing millions of nuclei and chloroplasts, and harbor diverse bacterial communities. Little is known about the diversity of chloroplast genomes (cpDNAs) in this group, and about the possible consequences of intracellular bacteria on genome composition of the host. We present the complete cpDNAs of Bryopsis plumosa and Tydemania expeditiones, as well as a re-annotated cpDNA of B. hypnoides, which was shown to contain a higher number of genes than originally published. Chloroplast genomic data were also used to evaluate phylogenetic hypotheses in the Chlorophyta, such as monophyly of the Ulvophyceae (the class in which the order Bryopsidales is currently classified).

Results

Both DNAs are circular and lack a large inverted repeat. The cpDNA of B. plumosa is 106,859 bp long and contains 115 unique genes. A 13 kb region was identified with several freestanding open reading frames (ORFs) of putative bacterial origin, including a large ORF (>8 kb) closely related to bacterial rhs-family genes. The cpDNA of T. expeditiones is 105,200 bp long and contains 125 unique genes. As in B. plumosa, several regions were identified with ORFs of possible bacterial origin, including genes involved in mobile functions (transposases, integrases, phage/plasmid DNA primases), and ORFs showing close similarity with bacterial DNA methyltransferases. The cpDNA of B. hypnoides differs from that of B. plumosa mainly in the presence of long intergenic spacers, and a large tRNA region. Chloroplast phylogenomic analyses were largely inconclusive with respect to monophyly of the Ulvophyceae, and the relationship of the Bryopsidales within the Chlorophyta.

Conclusions

The cpDNAs of B. plumosa and T. expeditiones are amongst the smallest and most gene dense chloroplast genomes in the core Chlorophyta. The presence of bacterial genes, including genes typically found in mobile elements, suggest that these have been acquired through horizontal gene transfer, which may have been facilitated by the occurrence of obligate intracellular bacteria in these siphonous algae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1418-3) contains supplementary material, which is available to authorized users.