Ancestral chloroplast genome in Mesostigma viride reveals an early branch of green plant evolution

Cytoplasmic genome contributions to domestication and improvement of modern maize

BACKGROUND

Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated.

RESULTS

Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes.

CONCLUSIONS

Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.

Chloroplast genome evolution and phylogeny of the early-diverging charophycean green algae with a focus on the Klebsormidiophyceae and Streptofilum

The Klebsormidiophyceae are a class of green microalgae observed globally in both freshwater and terrestrial habitats. Morphology-based classification schemes of this class have been shown to be inadequate due to the simple morphology of these algae, the tendency of morphology to vary in culture versus field conditions, and rampant morphological homoplasy. Molecular studies revealing cryptic diversity have renewed interest in this group. We sequenced the complete chloroplast genomes of a broad series of taxa spanning the known taxonomic breadth of this class. We also sequenced the chloroplast genomes of three strains of Streptofilum, a recently discovered green algal lineage with close affinity to the Klebsormidiophyceae. Our results affirm the previously hypothesized polyphyly of the genus Klebsormidium as well as the polyphyly of the nominal species in this genus, K. flaccidum. Furthermore, plastome sequences strongly support the status of Streptofilum as a distinct, early-diverging lineage of charophytic algae sister to a clade comprising Klebsormidiophyceae plus Phragmoplastophyta. We also uncovered major structural alterations in the chloroplast genomes of species in Klebsormidium that have broad implications regarding the underlying mechanisms of chloroplast genome evolution.

Assembly and analysis of the complete mitochondrial genome of Forsythia suspensa (Thunb.) Vahl

BACKGROUND

Forsythia suspensa (Thunb.) Vahl is a valuable ornamental and medicinal plant. Although the nuclear and chloroplast genomes of F. suspensa have been published, its complete mitochondrial genome sequence has yet to be reported. In this study, the genomic DNA of F. suspensa yellowish leaf material was extracted, sequenced by using a mixture of Illumina Novaseq6000 short reads and Oxford Nanopore PromethION long reads, and the sequencing data were assembled and annotated.

RESULT

The F. suspensa mitochondrial genome was obtained in the length of 535,692 bp with a circular structure, and the GC content was 44.90%. The genome contains 60 genes, including 36 protein-coding genes, 21 tRNA genes, and three rRNA genes. We further analyzed RNA editing of the protein-coding genes, relative synonymous codon usage, and sequence repeats based on the genomic data. There were 25 homologous sequences between F. suspensa mitochondria and chloroplast genome, which involved the transfer of 8 mitochondrial genes, and 9473 homologous sequences between mitochondrial and nuclear genomes. Analysis of the nucleic acid substitution rate, nucleic acid diversity, and collinearity of protein-coding genes of the F. suspensa mitochondrial genome revealed that the majority of genes may have undergone purifying selection, exhibiting a slower rate of evolution and a relatively conserved structure. Analysis of the phylogenetic relationships among different species revealed that F. suspensa was most closely related to Olea europaea subsp. Europaea.

CONCLUSION

In this study, we sequenced, assembled, and annotated a high-quality F. suspensa mitochondrial genome. The results of this study will enrich the mitochondrial genome data of Forsythia, lay a foundation for the phylogenetic development of Forsythia, and promote the evolutionary analysis of Oleaceae species.

Plastid Genome Evolution of Two Colony-Forming Benthic Ochrosphaera neapolitana Strains (Coccolithales, Haptophyta)

Coccolithophores are well-known haptophytes that produce small calcium carbonate coccoliths, which in turn contribute to carbon sequestration in the marine environment. Despite their important ecological role, only two of eleven haptophyte plastid genomes are from coccolithophores, and those two belong to the order Isochrysidales. Here, we report the plastid genomes of two strains of Ochrosphaera neapolitana (Coccolithales) from Spain (CCAC 3688 B) and the USA (A15,280). The newly constructed plastid genomes are the largest in size (116,906 bp and 113,686 bp, respectively) among all the available haptophyte plastid genomes, primarily due to the increased intergenic regions. These two plastid genomes possess a conventional quadripartite structure with a long single copy and short single copy separated by two inverted ribosomal repeats. These two plastid genomes share 110 core genes, six rRNAs, and 29 tRNAs, but CCAC 3688 B has an additional CDS (ycf55) and one tRNA (trnL-UAG). Two large insertions at the intergenic regions (2 kb insertion between ycf35 and ycf45; 0.5 kb insertion in the middle of trnM and trnY) were detected in the strain CCAC 3688 B. We found the genes of light-independent protochlorophyllide oxidoreductase (chlB, chlN, and chlL), which convert protochlorophyllide to chlorophyllide during chlorophyll biosynthesis, in the plastid genomes of O. neapolitana as well as in other benthic Isochrysidales and Coccolithales species, putatively suggesting an evolutionary adaptation to benthic habitats.

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.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnP^GGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Chloroplast genomes of five Oedogonium species: genome structure, phylogenetic analysis and adaptive evolution

Background

The order Oedogoniales within the single family Oedogoniaceae comprised of three genera, Oedogonium, Oedocladium, and Bulbochaete based on traditional morphological criteria. While several molecular phylogenetic studies have suggested that both Oedogonium and Oedocladium may not be monophyletic, broader taxon sampling and large amounts of molecular data acquisition could help to resolve the phylogeny and evolutionary problems of this order. This study determined five chloroplast (cp) genomes of Oedogonium species and aimed to provide further information on cp genome for a better understanding of the phylogenetic and evolutionary relationships of the order Oedogoniales.

Results

The five Oedogonium cp genomes showed typical quadripartite and circular structures, and were relatively conserved in their structure, gene synteny, and inverted repeats boundaries in general, except for small variation in genome sizes, AT contents, introns, and repeats. Phylogenetic analyses based on 54 cp protein-coding genes examined by maximum likelihood and Bayesian analyses using amino acid and nucleotide datasets indicated that both Oedocladium and Oedogonium are polyphyletic groups. A positively selected gene (psbA) was identified in the two Oedocladium species and the terrestrial Oedogonium species, indicating that terrestrial Oedogoniales taxa may have undergone adaptive evolution to adjust to the difference in light intensity between aquatic and terrestrial habitats.

Conclusions

Our results enrich the data on cp genomes of the genus Oedogonium. The availability of these cp genomes can help in understanding the cp genome characteristics and resolve phylogenetic and evolutionary relationships of the order Oedogoniales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08006-1.

Finding Species-Specific Extracellular Surface-Facing Proteomes in Toxic Dinoflagellates

As a sequel to our previous report of the existence of species-specific protein/peptide expression profiles (PEPs) acquired by mass spectrometry in some dinoflagellates, we established, with the help of a plasma-membrane-impermeable labeling agent, a surface amphiesmal protein extraction method (SAPE) to label and capture species-specific surface proteins (SSSPs) as well as saxitoxins-producing-species-specific surface proteins (Stx-SSPs) that face the extracellular space (i.e., SSSPsEf and Stx-SSPsEf). Five selected toxic dinoflagellates, Alexandrium minutum, A. lusitanicum, A. tamarense, Gymnodinium catenatum, and Karenia mikimotoi, were used in this study. Transcriptomic databases of these five species were also constructed. With the aid of liquid chromatography linked-tandem mass spectrometry (LC-MS/MS) and the transcriptomic databases of these species, extracellularly facing membrane proteomes of the five different species were identified. Within these proteomes, 16 extracellular-facing and functionally significant transport proteins were found. Furthermore, 10 SSSPs and 6 Stx-SSPs were identified as amphiesmal proteins but not facing outward to the extracellular environment. We also found SSSPsEf and Stx-SSPsEf in the proteomes. The potential functional correlation of these proteins towards the production of saxitoxins in dinoflagellates and the degree of species specificity were discussed accordingly.

Unique Peripheral Antennas in the Photosystems of the Streptophyte Alga Mesostigma viride

Land plants evolved from a single group of streptophyte algae. One of the key factors needed for adaptation to a land environment is the modification in the peripheral antenna systems of photosystems (PSs). Here, the PSs of Mesostigma viride, one of the earliest-branching streptophyte algae, were analyzed to gain insight into their evolution. Isoform sequencing and phylogenetic analyses of light-harvesting complexes (LHCs) revealed that M. viride possesses three algae-specific LHCs, including algae-type LHCA2, LHCA9 and LHCP, while the streptophyte-specific LHCB6 was not identified. These data suggest that the acquisition of LHCB6 and the loss of algae-type LHCs occurred after the M. viride lineage branched off from other streptophytes. Clear-native (CN)-polyacrylamide gel electrophoresis (PAGE) resolved the photosynthetic complexes, including the PSI-PSII megacomplex, PSII-LHCII, two PSI-LHCI-LHCIIs, PSI-LHCI and the LHCII trimer. Results indicated that the higher-molecular weight PSI-LHCI-LHCII likely had more LHCII than the lower-molecular weight one, a unique feature of M. viride PSs. CN-PAGE coupled with mass spectrometry strongly suggested that the LHCP was bound to PSII-LHCII, while the algae-type LHCA2 and LHCA9 were bound to PSI-LHCI, both of which are different from those in land plants. Results of the present study strongly suggest that M. viride PSs possess unique features that were inherited from a common ancestor of streptophyte and chlorophyte algae.

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ-stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNATyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences "U-U-x2-C" and "U-x-G-x2-T," thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.

The genome of Prasinoderma coloniale unveils the existence of a third phylum within green plants

Genome analysis of the pico-eukaryotic marine green alga Prasinoderma coloniale CCMP 1413 unveils the existence of a novel phylum within green plants (Viridiplantae), the Prasinodermophyta, which diverged before the split of Chlorophyta and Streptophyta. Structural features of the genome and gene family comparisons revealed an intermediate position of the P. coloniale genome (25.3 Mb) between the extremely compact, small genomes of picoplanktonic Mamiellophyceae (Chlorophyta) and the larger, more complex genomes of early-diverging streptophyte algae. Reconstruction of the minimal core genome of Viridiplantae allowed identification of an ancestral toolkit of transcription factors and flagellar proteins. Adaptations of P. coloniale to its deep-water, oligotrophic environment involved expansion of light-harvesting proteins, reduction of early light-induced proteins, evolution of a distinct type of C4 photosynthesis and carbon-concentrating mechanism, synthesis of the metal-complexing metabolite picolinic acid, and vitamin B1, B7 and B12 auxotrophy. The P. coloniale genome provides first insights into the dawn of green plant evolution.

Identification of High Molecular Variation Loci in Complete Chloroplast Genomes of Mammillaria (Cactaceae, Caryophyllales)

In plants, partial DNA sequences of chloroplasts have been widely used in evolutionary studies. However, the Cactaceae family (1500-1800 species) lacks molecular markers that allow a phylogenetic resolution between species and genera. In order to identify sequences with high variation levels, we compared previously reported complete chloroplast genomes of seven species of Mammillaria. We identified repeated sequences (RSs) and two types of DNA variation: short sequence repeats (SSRs) and divergent homologous loci. The species with the highest number of RSs was M. solisioides (256), whereas M. pectinifera contained the highest amount of SSRs (84). In contrast, M. zephyranthoides contained the lowest number (35) of both RSs and SSRs. In addition, five of the SSRs were found in the seven species, but only three of them showed variation. A total of 180 homologous loci were identified among the seven species. Out of these, 20 loci showed a molecular variation of 5% to 31%, and 12 had a length within the range of 150 to 1000 bp. We conclude that the high levels of variation at the reported loci represent valuable knowledge that may help to resolve phylogenetic relationships and that may potentially be convenient as molecular markers for population genetics and phylogeographic studies.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

Comparative analysis of plastid genomes within the Campanulaceae and phylogenetic implications

The conflicts exist between the phylogeny of Campanulaceae based on nuclear ITS sequence and plastid markers, particularly in the subdivision of Cyanantheae (Campanulaceae). Besides, various and complicated plastid genome structures can be found in species of the Campanulaceae. However, limited availability of genomic information largely hinders the studies of molecular evolution and phylogeny of Campanulaceae. We reported the complete plastid genomes of three Cyanantheae species, compared them to eight published Campanulaceae plastomes, and shed light on a deeper understanding of the applicability of plastomes. We found that there were obvious differences among gene order, GC content, gene compositions and IR junctions of LSC/IRa. Almost all protein-coding genes and amino acid sequences showed obvious codon preferences. We identified 14 genes with highly positively selected sites and branch-site model displayed 96 sites under potentially positive selection on the three lineages of phylogenetic tree. Phylogenetic analyses showed that Cyananthus was more closely related to Codonopsis compared with Cyclocodon and also clearly illustrated the relationship among the Cyanantheae species. We also found six coding regions having high nucleotide divergence value. Hotpot regions were considered to be useful molecular markers for resolving phylogenetic relationships and species authentication of Campanulaceae.

Mesostigma viride Genome and Transcriptome Provide Insights into the Origin and Evolution of Streptophyta

The Streptophyta include unicellular and multicellular charophyte green algae and land plants. Colonization of the terrestrial habitat by land plants is a major evolutionary event that has transformed the planet. So far, lack of genome information on unicellular charophyte algae hinders the understanding of the origin and the evolution from unicellular to multicellular life in Streptophyta. This work reports the high-quality reference genome and transcriptome of Mesostigma viride, a single-celled charophyte alga with a position at the base of Streptophyta. There are abundant segmental duplications and transposable elements in M. viride, which contribute to a relatively large genome with high gene content compared to other algae and early diverging land plants. This work identifies the origin of genetic tools that multicellular Streptophyta have inherited and key genetic innovations required for the evolution of land plants from unicellular aquatic ancestors. The findings shed light on the age-old questions of the evolution of multicellularity and the origin of land plants.

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.

Comparative and Phylogenetic Analyses of Ginger (Zingiber officinale) in the Family Zingiberaceae Based on the Complete Chloroplast Genome

Zingiber officinale, commonly known as ginger, is an important plant of the family Zingiberaceae and is widely used as an herbal medicine and condiment. The lack of chloroplast genomic information hinders molecular research and phylogenetic analysis on ginger. We introduced the complete chloroplast genome of Z. officinale and identified its phylogenetic position in Zingiberaceae. The chloroplast genome of Z. officinale is 162,621 bp with a four-part circular structure and 36.1% GC content. All 113 unique genes were annotated. A total of 78 simple sequence repeats (SSRs) and 42 long repeat sequences, which are potential areas for species authentication, were found. Comparative analysis revealed some highly variable regions, including rps16-trnQ-UUG, atpH-atpI, trnT-UGU-trnL-UAA, ycf1, and psaC-ndhE. Moreover, the small single-copy (SSC) region was the most variable region in all four shared regions, indicating that it may be undergoing rapid nucleotide substitution in the family Zingiberaceae. Phylogenetic analysis based on all available chloroplasts of Zingiberales in the National Center for Biotechnology Information indicated that Zingiber is a sister branch to Kaempferia species. The availability of the Z. officinale chloroplast genome provided invaluable data for species-level authentication and phylogenetic analysis and can thus benefit further investigations on species in the family Zingiberaceae.

The Carboxy Terminus of YCF1 Contains a Motif Conserved throughout >500 Myr of Streptophyte Evolution

Plastids evolved from cyanobacteria by endosymbiosis. During the course of evolution, the coding capacity of plastid genomes shrinks due to gene loss or transfer to the nucleus. In the green lineage, however, there were apparent gene gains including that of ycf1. Although its function is still debated, YCF1 has proven to be a useful marker for plastid evolution. YCF1 sequence and predicted structural features unite the plastid genomes of land plants with those of their closest algal relatives, the higher streptophyte algae; YCF1 appears to have undergone pronounced changes during the course of streptophyte algal evolution. Using new data, we show that YCF1 underwent divergent evolution in the common ancestor of higher streptophyte algae and Klebsormidiophycae. This divergence resulted in the origin of an extreme, klebsormidiophycean-specific YCF1 and the higher streptophyte Ste-YCF1. Most importantly, our analysis uncovers a conserved carboxy-terminal sequence stretch within YCF1 that is unique to higher streptophytes and hints at an important, yet unexplored function.

Phylogenomics of tomato chloroplasts using assembly and alignment-free method

Phylogenetics and population genetics are central disciplines in evolutionary biology. Both are based on the comparison of single DNA sequences, or a concatenation of a number of these. However, with the advent of next-generation DNA sequencing technologies, the approaches that consider large genomic data sets are of growing importance for the elucidation of evolutionary relationships among species. Among these approaches, the assembly and alignment-free methods which allow an efficient distance computation and phylogeny reconstruction are of great importance. However, it is not yet clear under what quality conditions and abundance of genomic data such methods are able to infer phylogenies accurately. In the present study we assess the method originally proposed by Fan et al. for whole genome data, in the elucidation of Tomatoes' chloroplast phylogenetics using short read sequences. We find that this assembly and alignment-free method is capable of reproducing previous results under conditions of high coverage, given that low frequency k-mers (i.e. error prone data) are effectively filtered out. Finally, we present a complete chloroplast phylogeny for the best data quality candidates of the recently published 360 tomato genomes.

Evolutionary Dynamics of Cryptophyte Plastid Genomes

Cryptophytes are an ecologically important group of largely photosynthetic unicellular eukaryotes. This lineage is of great interest to evolutionary biologists because their plastids are of red algal secondary endosymbiotic origin and the host cell retains four different genomes (host nuclear, mitochondrial, plastid, and red algal nucleomorph). Here, we report a comparative analysis of plastid genomes from six representative cryptophyte genera. Four newly sequenced cryptophyte plastid genomes of Chroomonas mesostigmatica, Ch. placoidea, Cryptomonas curvata, and Storeatula sp. CCMP1868 share a number of features including synteny and gene content with the previously sequenced genomes of Cryptomonas paramecium, Rhodomonas salina, Teleaulax amphioxeia, and Guillardia theta. Our analysis of these plastid genomes reveals examples of gene loss and intron insertion. In particular, the chlB/chlL/chlN genes, which encode light-independent (dark active) protochlorophyllide oxidoreductase (LIPOR) proteins have undergone recent gene loss and pseudogenization in cryptophytes. Comparison of phylogenetic trees based on plastid and nuclear genome data sets show the introduction, via secondary endosymbiosis, of a red algal derived plastid in a lineage of chlorophyll-c containing algae. This event was followed by additional rounds of eukaryotic endosymbioses that spread the red lineage plastid to diverse groups such as haptophytes and stramenopiles.

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.

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.

Evolutionary Dynamics of Chloroplast Genomes in Low Light: A Case Study of the Endolithic Green Alga Ostreobium quekettii

Some photosynthetic organisms live in extremely low light environments. Light limitation is associated with selective forces as well as reduced exposure to mutagens, and over evolutionary timescales it can leave a footprint on species’ genomes. Here, we present the chloroplast genomes of four green algae (Bryopsidales, Ulvophyceae), including the endolithic (limestone-boring) alga Ostreobium quekettii, which is a low light specialist. We use phylogenetic models and comparative genomic tools to investigate whether the chloroplast genome of Ostreobium corresponds to our expectations of how low light would affect genome evolution. Ostreobium has the smallest and most gene-dense chloroplast genome among Ulvophyceae reported to date, matching our expectation that light limitation would impose resource constraints reflected in the chloroplast genome architecture. Rates of molecular evolution are significantly slower along the phylogenetic branch leading to Ostreobium, in agreement with the expected effects of low light and energy levels on molecular evolution. We expected the ability of Ostreobium to perform photosynthesis in very low light to be associated with positive selection in genes related to the photosynthetic machinery, but instead, we observed that these genes may be under stronger purifying selection. Besides shedding light on the genome dynamics associated with a low light lifestyle, this study helps to resolve the role of environmental factors in shaping the diversity of genome architectures observed in nature.

Mitochondrion-to-Chloroplast DNA Transfers and Intragenomic Proliferation of Chloroplast Group II Introns in Gloeotilopsis Green Algae (Ulotrichales, Ulvophyceae)

To probe organelle genome evolution in the Ulvales/Ulotrichales clade, the newly sequenced chloroplast and mitochondrial genomes of Gloeotilopsis planctonica and Gloeotilopsis sarcinoidea (Ulotrichales) were compared with those of Pseudendoclonium akinetum (Ulotrichales) and of the few other green algae previously sampled in the Ulvophyceae. At 105,236 bp, the G planctonica mitochondrial DNA (mtDNA) is the largest mitochondrial genome reported so far among chlorophytes, whereas the 221,431-bp G planctonica and 262,888-bp G sarcinoidea chloroplast DNAs (cpDNAs) are the largest chloroplast genomes analyzed among the Ulvophyceae. Gains of non-coding sequences largely account for the expansion of these genomes. Both Gloeotilopsis cpDNAs lack the inverted repeat (IR) typically found in green plants, indicating that two independent IR losses occurred in the Ulvales/Ulotrichales. Our comparison of the Pseudendoclonium and Gloeotilopsis cpDNAs offered clues regarding the mechanism of IR loss in the Ulotrichales, suggesting that internal sequences from the rDNA operon were differentially lost from the two original IR copies during this process. Our analyses also unveiled a number of genetic novelties. Short mtDNA fragments were discovered in two distinct regions of the G sarcinoidea cpDNA, providing the first evidence for intracellular inter-organelle gene migration in green algae. We identified for the first time in green algal organelles, group II introns with LAGLIDADG ORFs as well as group II introns inserted into untranslated gene regions. We discovered many group II introns occupying sites not previously documented for the chloroplast genome and demonstrated that a number of them arose by intragenomic proliferation, most likely through retrohoming.

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.

Chloroplast phylogenomic analyses reveal the deepest-branching lineage of the Chlorophyta, Palmophyllophyceae class. nov

The green plants (Viridiplantae) are an ancient group of eukaryotes comprising two main clades: the Chlorophyta, which includes a wide diversity of green algae, and the Streptophyta, which consists of freshwater green algae and the land plants. The early-diverging lineages of the Viridiplantae comprise unicellular algae, and multicellularity has evolved independently in the two clades. Recent molecular data have revealed an unrecognized early-diverging lineage of green plants, the Palmophyllales, with a unique form of multicellularity, and typically found in deep water. The phylogenetic position of this enigmatic group, however, remained uncertain. Here we elucidate the evolutionary affinity of the Palmophyllales using chloroplast genomic, and nuclear rDNA data. Phylogenetic analyses firmly place the palmophyllalean Verdigellas peltata along with species of Prasinococcales (prasinophyte clade VI) in the deepest-branching clade of the Chlorophyta. The small, compact and intronless chloroplast genome (cpDNA) of V. peltata shows striking similarities in gene content and organization with the cpDNAs of Prasinococcales and the streptophyte Mesostigma viride, indicating that cpDNA architecture has been extremely well conserved in these deep-branching lineages of green plants. The phylogenetic distinctness of the Palmophyllales-Prasinococcales clade, characterized by unique ultrastructural features, warrants recognition of a new class of green plants, Palmophyllophyceae class. nov.

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.

Complete Chloroplast Genome Sequence of Aquilaria sinensis (Lour.) Gilg and Evolution Analysis within the Malvales Order

Aquilaria sinensis (Lour.) Gilg is an important medicinal woody plant producing agarwood, which is widely used in traditional Chinese medicine. High-throughput sequencing of chloroplast (cp) genomes enhanced the understanding about evolutionary relationships within plant families. In this study, we determined the complete cp genome sequences for A. sinensis. The size of the A. sinensis cp genome was 159,565 bp. This genome included a large single-copy region of 87,482 bp, a small single-copy region of 19,857 bp, and a pair of inverted repeats (IRa and IRb) of 26,113 bp each. The GC content of the genome was 37.11%. The A. sinensis cp genome encoded 113 functional genes, including 82 protein-coding genes, 27 tRNA genes, and 4 rRNA genes. Seven genes were duplicated in the protein-coding genes, whereas 11 genes were duplicated in the RNA genes. A total of 45 polymorphic simple-sequence repeat loci and 60 pairs of large repeats were identified. Most simple-sequence repeats were located in the noncoding sections of the large single-copy/small single-copy region and exhibited high A/T content. Moreover, 33 pairs of large repeat sequences were located in the protein-coding genes, whereas 27 pairs were located in the intergenic regions. Aquilaria sinensis cp genome bias ended with A/T on the basis of codon usage. The distribution of codon usage in A. sinensis cp genome was most similar to that in the Gonystylus bancanus cp genome. Comparative results of 82 protein-coding genes from 29 species of cp genomes demonstrated that A. sinensis was a sister species to G. bancanus within the Malvales order. Aquilaria sinensis cp genome presented the highest sequence similarity of >90% with the G. bancanus cp genome by using CGView Comparison Tool. This finding strongly supports the placement of A. sinensis as a sister to G. bancanus within the Malvales order. The complete A. sinensis cp genome information will be highly beneficial for further studies on this traditional medicinal plant. Moreover, the results will enhance our understanding about the evolution of cp genomes of the Malvales order, particularly with regard to the role of A. sinensis in plant systematics and evolution.

Complete Plastid Genome Sequence of the Brown Alga Undaria pinnatifida

In this study, we fully sequenced the circular plastid genome of a brown alga, Undaria pinnatifida. The genome is 130,383 base pairs (bp) in size; it contains a large single-copy (LSC, 76,598 bp) and a small single-copy region (SSC, 42,977 bp), separated by two inverted repeats (IRa and IRb: 5,404 bp). The genome contains 139 protein-coding, 28 tRNA, and 6 rRNA genes; none of these genes contains introns. Organization and gene contents of the U. pinnatifida plastid genome were similar to those of Saccharina japonica. There is a co-linear relationship between the plastid genome of U. pinnatifida and that of three previously sequenced large brown algal species. Phylogenetic analyses of 43 taxa based on 23 plastid protein-coding genes grouped all plastids into a red or green lineage. In the large brown algae branch, U. pinnatifida and S. japonica formed a sister clade with much closer relationship to Ectocarpus siliculosus than to Fucus vesiculosus. For the first time, the start codon ATT was identified in the plastid genome of large brown algae, in the atpA gene of U. pinnatifida. In addition, we found a gene-length change induced by a 3-bp repetitive DNA in ycf35 and ilvB genes of the U. pinnatifida plastid genome.

Complete chloroplast genome sequence of poisonous and medicinal plant Datura stramonium: organizations and implications for genetic engineering

Datura stramonium is a widely used poisonous plant with great medicinal and economic value. Its chloroplast (cp) genome is 155,871 bp in length with a typical quadripartite structure of the large (LSC, 86,302 bp) and small (SSC, 18,367 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,601 bp). The genome contains 113 unique genes, including 80 protein-coding genes, 29 tRNAs and four rRNAs. A total of 11 forward, 9 palindromic and 13 tandem repeats were detected in the D. stramonium cp genome. Most simple sequence repeats (SSR) are AT-rich and are less abundant in coding regions than in non-coding regions. Both SSRs and GC content were unevenly distributed in the entire cp genome. All preferred synonymous codons were found to use A/T ending codons. The difference in GC contents of entire genomes and of the three-codon positions suggests that the D. stramonium cp genome might possess different genomic organization, in part due to different mutational pressures. The five most divergent coding regions and four non-coding regions (trnH-psbA, rps4-trnS, ndhD-ccsA, and ndhI-ndhG) were identified using whole plastome alignment, which can be used to develop molecular markers for phylogenetics and barcoding studies within the Solanaceae. Phylogenetic analysis based on 68 protein-coding genes supported Datura as a sister to Solanum. This study provides valuable information for phylogenetic and cp genetic engineering studies of this poisonous and medicinal plant.

Green algae and the origins of multicellularity in the plant kingdom

The green lineage of chlorophyte algae and streptophytes form a large and diverse clade with multiple independent transitions to produce multicellular and/or macroscopically complex organization. In this review, I focus on two of the best-studied multicellular groups of green algae: charophytes and volvocines. Charophyte algae are the closest relatives of land plants and encompass the transition from unicellularity to simple multicellularity. Many of the innovations present in land plants have their roots in the cell and developmental biology of charophyte algae. Volvocine algae evolved an independent route to multicellularity that is captured by a graded series of increasing cell-type specialization and developmental complexity. The study of volvocine algae has provided unprecedented insights into the innovations required to achieve multicellularity.

Six newly sequenced chloroplast genomes from prasinophyte green algae provide insights into the relationships among prasinophyte lineages and the diversity of streamlined genome architecture in picoplanktonic species

Background

Because they represent the earliest divergences of the Chlorophyta, the morphologically diverse unicellular green algae making up the prasinophytes hold the key to understanding the nature of the first viridiplants and the evolutionary patterns that accompanied the radiation of chlorophytes. Nuclear-encoded 18S rDNA phylogenies unveiled nine prasinophyte clades (clades I through IX) but their branching order is still uncertain. We present here the newly sequenced chloroplast genomes of Nephroselmis astigmatica (clade III) and of five picoplanktonic species from clade VI (Prasinococcus sp. CCMP 1194, Prasinophyceae sp. MBIC 106222 and Prasinoderma coloniale) and clade VII (Picocystis salinarum and Prasinophyceae sp. CCMP 1205). These chloroplast DNAs (cpDNAs) were compared with those of the six previously sampled prasinophytes (clades I, II, III and V) in order to gain information both on the relationships among prasinophyte lineages and on chloroplast genome evolution.

Results

Varying from 64.3 to 85.6 kb in size and encoding 100 to 115 conserved genes, the cpDNAs of the newly investigated picoplanktonic species are substantially smaller than those observed for larger-size prasinophytes, are economically packed and contain a reduced gene content. Although the Nephroselmis and Picocystis cpDNAs feature a large inverted repeat encoding the rRNA operon, gene partitioning among the single copy regions is remarkably different. Unexpectedly, we found that all three species from clade VI (Prasinococcales) harbor chloroplast genes not previously documented for chlorophytes (ndhJ, rbcR, rpl21, rps15, rps16 and ycf66) and that Picocystis contains a trans-spliced group II intron. The phylogenies inferred from cpDNA-encoded proteins are essentially congruent with 18S rDNA trees, resolving with robust support all six examined prasinophyte lineages, with the exception of the Pycnococcaceae.

Conclusions

Our results underscore the high variability in genome architecture among prasinophyte lineages, highlighting the strong pressure to maintain a small and compact chloroplast genome in picoplanktonic species. The unique set of six chloroplast genes found in the Prasinococcales supports the ancestral status of this lineage within the prasinophytes. The widely diverging traits uncovered for the clade-VII members (Picocystis and Prasinophyceae sp. CCMP 1205) are consistent with their resolution as separate lineages in the chloroplast phylogeny.

Overly simplistic substitution models obscure green plant phylogeny

Phylogenetic analysis is an increasingly common and valuable component of plant science. Knowledge of the phylogenetic relationships between plant groups is a prerequisite for understanding the origin and evolution of important plant features, and phylogenetic analysis of individual genes and gene families provides fundamental insights into how those genes and their functions evolved. However, despite an active research community exploring and improving phylogenetic methods, the analytical methods commonly used, and the phylogenetic results they produce, are accorded far more confidence than they warrant. In this opinion article, I emphasise that important parts of the green plant phylogeny are inconsistently resolved and I argue that the lack of consistency arises due to inadequate modelling of changes in the substitution process.

Applications of next-generation sequencing to unravelling the evolutionary history of algae

First-generation Sanger DNA sequencing revolutionized science over the past three decades and the current next-generation sequencing (NGS) technology has opened the doors to the next phase in the sequencing revolution. Using NGS, scientists are able to sequence entire genomes and to generate extensive transcriptome data from diverse photosynthetic eukaryotes in a timely and cost-effective manner. Genome data in particular shed light on the complicated evolutionary history of algae that form the basis of the food chain in many environments. In the Eukaryotic Tree of Life, the fact that photosynthetic lineages are positioned in four supergroups has important evolutionary consequences. We now know that the story of eukaryotic photosynthesis unfolds with a primary endosymbiosis between an ancestral heterotrophic protist and a captured cyanobacterium that gave rise to the glaucophytes, red algae and Viridiplantae (green algae and land plants). These primary plastids were then transferred to other eukaryotic groups through secondary endosymbiosis. A red alga was captured by the ancestor(s) of the stramenopiles, alveolates (dinoflagellates, apicomplexa, chromeridae), cryptophytes and haptophytes, whereas green algae were captured independently by the common ancestors of the euglenophytes and chlorarachniophytes. A separate case of primary endosymbiosis is found in the filose amoeba Paulinella chromatophora, which has at least nine heterotrophic sister species. Paulinella genome data provide detailed insights into the early stages of plastid establishment. Therefore, genome data produced by NGS have provided many novel insights into the taxonomy, phylogeny and evolutionary history of photosynthetic eukaryotes.

From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

BACKGROUND

Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data.

RESULTS

We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa.

CONCLUSIONS

Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.

Heterotrimeric G proteins in green algae: an early innovation in the evolution of the plant lineage

Heterotrimeric G-proteins (G-proteins, hereafter) are important signaling components in all eukaryotes. The absence of these proteins in the sequenced genomes of Chlorophyaceaen green algae has raised questions about their evolutionary origin and prevalence in the plant lineage. The existence of G-proteins has often been correlated with the acquisition of embryophytic life-cycle and/or terrestrial habitats of plants which occurred around 450 million years ago. Our discovery of functional G-proteins in Chara braunii, a representative of the Charophycean green algae, establishes the existence of this conserved signaling pathway in the most basal plants and dates it even further back to 1-1.5 billion years ago. We have now identified the sequence homologs of G-proteins in additional algal families and propose that green algae represent a model system for one of the most basal forms of G-protein signaling known to exist to date. Given the possible differences that exist between plant and metazoan G-protein signaling mechanisms, such basal organisms will serve as important resources to trace the evolutionary origin of proposed mechanistic differences between the systems as well as their plant-specific functions.

Agrigenomics for microalgal biofuel production: an overview of various bioinformatics resources and recent studies to link OMICS to bioenergy and bioeconomy

Microalgal biofuels offer great promise in contributing to the growing global demand for alternative sources of renewable energy. However, to make algae-based fuels cost competitive with petroleum, lipid production capabilities of microalgae need to improve substantially. Recent progress in algal genomics, in conjunction with other "omic" approaches, has accelerated the ability to identify metabolic pathways and genes that are potential targets in the development of genetically engineered microalgal strains with optimum lipid content. In this review, we summarize the current bioeconomic status of global biofuel feedstocks with particular reference to the role of "omics" in optimizing sustainable biofuel production. We also provide an overview of the various databases and bioinformatics resources available to gain a more complete understanding of lipid metabolism across algal species, along with the recent contributions of "omic" approaches in the metabolic pathway studies for microalgal biofuel production.

The wheat chloroplastic proteome

With the availability of plant genome sequencing, analysis of plant proteins with mass spectrometry has become promising and admired. Determining the proteome of a cell is still a challenging assignment, which is convoluted by proteome dynamics and convolution. Chloroplast is fastidious curiosity for plant biologists due to their intricate biochemical pathways for indispensable metabolite functions. In this review, an overview on proteomic studies conducted in wheat with a special focus on subcellular proteomics of chloroplast, salt and water stress. In recent years, we and other groups have attempted to understand the photosynthesis in wheat and abiotic stress under salt imposed and water deficit during vegetative stage. Those studies provide interesting results leading to better understanding of the photosynthesis and identifying the stress-responsive proteins. Indeed, recent studies aimed at resolving the photosynthesis pathway in wheat. Proteomic analysis combining two complementary approaches such as 2-DE and shotgun methods couple to high through put mass spectrometry (LTQ-FTICR and MALDI-TOF/TOF) in order to better understand the responsible proteins in photosynthesis and abiotic stress (salt and water) in wheat chloroplast will be focused.

BIOLOGICAL SIGNIFICANCE

In this review we discussed the identification of the most abundant protein in wheat chloroplast and stress-responsive under salt and water stress in chloroplast of wheat seedlings, thus providing the proteomic view of the events during the development of this seedling under stress conditions. Chloroplast is fastidious curiosity for plant biologists due to their intricate biochemical pathways for indispensable metabolite functions. An overview on proteomic studies conducted in wheat with a special focus on subcellular proteomics of chloroplast, salt and water stress. We have attempted to understand the photosynthesis in wheat and abiotic stress under salt imposed and water deficit during seedling stage. Those studies provide interesting results leading to a better understanding of the photosynthesis and identifying the stress-responsive proteins. In reality, our studies aspired at resolving the photosynthesis pathway in wheat. Proteomic analysis united two complementary approaches such as Tricine SDS-PAGE and 2-DE methods couple to high through put mass spectrometry (LTQ-FTICR and MALDI-TOF/TOF) in order to better understand the responsible proteins in photosynthesis and abiotic stress (salt and water) in wheat chloroplast will be highlighted. This article is part of a Special Issue entitled: Translational Plant Proteomics.

RNA processing and decay in plastids

Plastids were derived through endosymbiosis from a cyanobacterial ancestor, whose uptake was followed by massive gene transfer to the nucleus, resulting in the compact size and modest coding capacity of the extant plastid genome. Plastid gene expression is essential for plant development, but depends on nucleus-encoded proteins recruited from cyanobacterial or host-cell origins. The plastid genome is heavily transcribed from numerous promoters, giving posttranscriptional events a critical role in determining the quantity and sizes of accumulating RNA species. The major events reviewed here are RNA editing, which restores protein conservation or creates correct open reading frames by converting C residues to U, RNA splicing, which occurs both in cis and trans, and RNA cleavage, which relies on a variety of exoribonucleases and endoribonucleases. Because the RNases have little sequence specificity, they are collectively able to remove extraneous RNAs whose ends are not protected by RNA secondary structures or sequence-specific RNA-binding proteins (RBPs). Other plastid RBPs, largely members of the helical-repeat superfamily, confer specificity to editing and splicing reactions. The enzymes that catalyze RNA processing are also the main actors in RNA decay, implying that these antagonistic roles are optimally balanced. We place the actions of RBPs and RNases in the context of a recent proteomic analysis that identifies components of the plastid nucleoid, a protein-DNA complex with multiple roles in gene expression. These results suggest that sublocalization and/or concentration gradients of plastid proteins could underpin the regulation of RNA maturation and degradation.

Uncovering the protein translocon at the chloroplast inner envelope membrane

Chloroplasts require protein translocons at the outer and inner envelope membranes, termed TOC and TIC, respectively, to import thousands of cytoplasmically synthesized preproteins. However, the molecular identity of the TIC translocon remains controversial. Tic20 forms a 1-megadalton complex at the inner membrane and directly interacts with translocating preproteins. We purified the 1-megadalton complex from Arabidopsis, comprising Tic20 and three other essential components, one of which is encoded by the enigmatic open reading frame ycf1 in the chloroplast genome. All four components, together with well-known TOC components, were found stoichiometrically associated with different translocating preproteins. When reconstituted into planar lipid bilayers, the purified complex formed a preprotein-sensitive channel. Thus, this complex constitutes a general TIC translocon.

Complex RNA metabolism in the chloroplast: an update on the psbB operon

Expression of most plastid genes involves multiple post-transcriptional processing events, such as splicing, editing, and intercistronic processing. The latter involves the formation of mono-, di-, and multicistronic transcripts, which can further be regulated by differential stability and expression. The plastid pentacistronic psbB transcription unit has been well characterized in vascular plants. It encodes the subunits CP47 (psbB), T (psbT), and H (psbH) of photosystem II as well as cytochrome b (6) (petB) and subunit IV (petD) of the cytochrome b (6) f complex. Each of the petB and petD genes contains a group II intron, which is spliced during post-transcriptional modification. The small subunit of photosystem II, PsbN, is encoded in the intercistronic region between psbH and psbT but is transcribed in the opposite direction. Expression of the psbB gene cluster necessitates different processing events along with numerous newly evolved specificity factors conferring stability to many of the processed RNA transcripts, and thus exemplarily shows the complexity of RNA metabolism in the chloroplast.

Nested in the Chlorellales or independent class? Phylogeny and classification of the Pedinophyceae (Viridiplantae) revealed by molecular phylogenetic analyses of complete nuclear and plastid-encoded rRNA operons

The class Pedinophyceae was established for asymmetric uniflagellate green algae, and was originally considered as an ancestral lineage of viridiplants. However, analyses of 71 concatenated plastid proteins [Turmel et al. (2009): Mol. Biol. Evol. 26: 2317-2331] recovered Pedinomonas within the Chlorellales (Trebouxiophyceae), thereby questioning the Pedinophyceae as an independent class. For the present study, complete nuclear and plastid-encoded rRNA operon sequences have been determined for 37 taxa of green algae including 6 members of the Pedinophyceae, providing 9272 aligned nucleotide positions. Phylogenies using both rRNA operons consistently rejected any relationship between Pedinophyceae and the Chlorellales. Instead, the Pedinophyceae were significantly resolved as sister of all phycoplast-containing 'core' chlorophytes, i.e. Chlorodendrophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Reinvestigation of plastid proteins discovered biased phylogenetic signal among protein partitions, indicating the published Pedinomonas + Chlorellales association as likely artificial. Marine pedinophytes (Resultomonas and Marsupiomonas; Marsupiomonadales ord. nov.), formed a sister clade to the order Pedinomonadales, occurring in freshwater and soil habitats. Synapomorphies in rRNA secondary structures were integrated in taxonomic diagnoses of the Pedinophyceae and were also used for BLAST searches targeting environmental sequence databases. The latter approach revealed conserved habitat preferences for the Marsupiomonadales and Pedinomonadales, and identified several novel pedinophyte lineages.

Mutational patterns in RNA secondary structure evolution examined in three RNA families

The goal of this work was to study mutational patterns in the evolution of RNA secondary structure. We analyzed bacterial tmRNA, RNaseP and eukaryotic telomerase RNA secondary structures, mapping structural variability onto phylogenetic trees constructed primarily from rRNA sequences. We found that secondary structures evolve both by whole stem insertion/deletion, and by mutations that create or disrupt stem base pairing. We analyzed the evolution of stem lengths and constructed substitution matrices describing the changes responsible for the variation in the RNA stem length. In addition, we used principal component analysis of the stem length data to determine the most variable stems in different families of RNA. This data provides new insights into the evolution of RNA secondary structures and patterns of variation in the lengths of double helical regions of RNA molecules. Our findings will facilitate design of improved mutational models for RNA structure evolution.

Evolutionary origins of Brassicaceae specific genes in Arabidopsis thaliana

Background

All sequenced genomes contain a proportion of lineage-specific genes, which exhibit no sequence similarity to any genes outside the lineage. Despite their prevalence, the origins and functions of most lineage-specific genes remain largely unknown. As more genomes are sequenced opportunities for understanding evolutionary origins and functions of lineage-specific genes are increasing.

Results

This study provides a comprehensive analysis of the origins of lineage-specific genes (LSGs) in Arabidopsis thaliana that are restricted to the Brassicaceae family. In this study, lineage-specific genes within the nuclear (1761 genes) and mitochondrial (28 genes) genomes are identified. The evolutionary origins of two thirds of the lineage-specific genes within the Arabidopsis thaliana genome are also identified. Almost a quarter of lineage-specific genes originate from non-lineage-specific paralogs, while the origins of ~10% of lineage-specific genes are partly derived from DNA exapted from transposable elements (twice the proportion observed for non-lineage-specific genes). Lineage-specific genes are also enriched in genes that have overlapping CDS, which is consistent with such novel genes arising from overprinting. Over half of the subset of the 958 lineage-specific genes found only in Arabidopsis thaliana have alignments to intergenic regions in Arabidopsis lyrata, consistent with either de novo origination or differential gene loss and retention, with both evolutionary scenarios explaining the lineage-specific status of these genes. A smaller number of lineage-specific genes with an incomplete open reading frame across different Arabidopsis thaliana accessions are further identified as accession-specific genes, most likely of recent origin in Arabidopsis thaliana. Putative de novo origination for two of the Arabidopsis thaliana-only genes is identified via additional sequencing across accessions of Arabidopsis thaliana and closely related sister species lineages. We demonstrate that lineage-specific genes have high tissue specificity and low expression levels across multiple tissues and developmental stages. Finally, stress responsiveness is identified as a distinct feature of Brassicaceae-specific genes; where these LSGs are enriched for genes responsive to a wide range of abiotic stresses.

Conclusion

Improving our understanding of the origins of lineage-specific genes is key to gaining insights regarding how novel genes can arise and acquire functionality in different lineages. This study comprehensively identifies all of the Brassicaceae-specific genes in Arabidopsis thaliana and identifies how the majority of such lineage-specific genes have arisen. The analysis allows the relative importance (and prevalence) of different evolutionary routes to the genesis of novel ORFs within lineages to be assessed. Insights regarding the functional roles of lineage-specific genes are further advanced through identification of enrichment for stress responsiveness in lineage-specific genes, highlighting their likely importance for environmental adaptation strategies.

Internal transcribed spacer 2 (nu ITS2 rRNA) sequence-structure phylogenetics: towards an automated reconstruction of the green algal tree of life

BACKGROUND

Chloroplast-encoded genes (matK and rbcL) have been formally proposed for use in DNA barcoding efforts targeting embryophytes. Extending such a protocol to chlorophytan green algae, though, is fraught with problems including non homology (matK) and heterogeneity that prevents the creation of a universal PCR toolkit (rbcL). Some have advocated the use of the nuclear-encoded, internal transcribed spacer two (ITS2) as an alternative to the traditional chloroplast markers. However, the ITS2 is broadly perceived to be insufficiently conserved or to be confounded by introgression or biparental inheritance patterns, precluding its broad use in phylogenetic reconstruction or as a DNA barcode. A growing body of evidence has shown that simultaneous analysis of nucleotide data with secondary structure information can overcome at least some of the limitations of ITS2. The goal of this investigation was to assess the feasibility of an automated, sequence-structure approach for analysis of IT2 data from a large sampling of phylum Chlorophyta.

METHODOLOGY/PRINCIPAL FINDINGS

Sequences and secondary structures from 591 chlorophycean, 741 trebouxiophycean and 938 ulvophycean algae, all obtained from the ITS2 Database, were aligned using a sequence structure-specific scoring matrix. Phylogenetic relationships were reconstructed by Profile Neighbor-Joining coupled with a sequence structure-specific, general time reversible substitution model. Results from analyses of the ITS2 data were robust at multiple nodes and showed considerable congruence with results from published phylogenetic analyses.

CONCLUSIONS/SIGNIFICANCE

Our observations on the power of automated, sequence-structure analyses of ITS2 to reconstruct phylum-level phylogenies of the green algae validate this approach to assessing diversity for large sets of chlorophytan taxa. Moreover, our results indicate that objections to the use of ITS2 for DNA barcoding should be weighed against the utility of an automated, data analysis approach with demonstrated power to reconstruct evolutionary patterns for highly divergent lineages.