Organelle genome complexity scales positively with organism size in volvocine green algae

One Ring does not rule them all: Linear mtDNA in Metazoa

Recent advances in genome sequencing technologies have facilitated the exploration of the architecture of genomes, including mitochondrial genomes (mtDNA). In particular, whole genome sequencing has provided easier access to mitochondrial genomes with unusual organizations, which were difficult to obtain using traditional PCR-based approaches. As a consequence, there has been a steep increase in complete mtDNA sequences, particularly for Metazoa. The popular view of metazoan mtDNA is that of a small gene-dense circular chromosome. This view clashes with discoveries of a number of linear mtDNAs, particularly in non-bilaterian animals. Here, we review the distribution of linear mtDNA in Metazoa, namely in isopods, cnidarians, and sponges. We discuss the multiple origins of linear mitogenomes in these clades, where linearity has been linked to the likely insertion of a linear plasmid in cnidarians and the demosponge Acanthella acuta, while fixation of a heteroplasmy in the anticodon site of a tRNA might be responsible for the monolinear form of the mtDNA in some isopods. We also summarize our current knowledge of mechanisms that maintain the integrity of linear mitochromosomes, where a recurrent theme is the presence of terminal repeats that likely play the role of telomeres. We caution in defining a linear chromosome as complete, particularly when coding sequences and key features of linear DNA are missing. Finally, we encourage authors interested in mitogenome science to utilize all available data for linear mtDNA, including those tagged as "incomplete" or "unverified" in public databases, as they can still provide useful information such as phylogenetic characters and gene order.

The genome of the Arctic snow alga Limnomonas spitsbergensis (Chlamydomonadales)

Snow algae are a diverse group of extremophilic microeukaryotes found on melting polar and alpine snowfields. They play an important role in the microbial ecology of the cryosphere, and their propagation on snow and ice surfaces may in part accelerate climate-induced melting of these systems. High-quality snow algae genomes are needed for studies on their unique physiology, adaptive mechanisms, and genome evolution under multiple forms of stress, including cold temperatures and intense sunlight. Here, we assembled and annotated the genome of Limnomonas spitsbergensis, a cryophilic biciliate green alga originally isolated from melting snow on Svalbard, in the Arctic. The L. spitsbergensis genome assembly is based primarily on the use of PacBio long reads and secondly Illumina short reads, with an assembly size of 260.248 Mb in 124 contigs. A combination of 3 alternative annotation strategies was used including protein homology, RNA-seq evidence, and PacBio full-length transcript isoforms. The best merged set of annotations identified 18,277 protein-coding genes, which were 95.2% complete based on Benchmarking Universal Single-Copy Orthologs analysis. We also provide the annotated mitogenome, which is a relatively large 77.942 kb circular mapping sequence containing extensive repeats. The L. spitsbergensis genome will provide a new resource for research on snow algae adaptation, behavior, and natural selection in unique, low-temperature terrestrial environments that are under threat from climate change.

Characterization of the mitochondrial genome of Chlorolobion braunii ITBB-AG6, an azolla-associated green alga isolated from sanitary sewage

Sphaeropleales have the characteristics of rapid growth, high oil content, and efficient removal rates of nitrogen and phosphorus in sewage waters, and is potentially valuable in biodiesel production and environmental remediation. In this study, we isolated a strain of Sphaeropleales, Chlorolobion braunii strain ITBB-AG6 from an azolla community in a sewage pond. Its mitochondrial genome contains 110,124 bp and harbors at least 40 genes, including 15 protein-coding genes, 20 tRNA genes, and three rRNA genes. The protein-coding genes include two for ATP synthases, seven for NAD(P)H-quinone oxidoreductases (nad), three for cytochrome c oxidase subunits (coxs), and one for cytochrome b (cob). Transfer RNA genes for 18 amino acids were identified, in which the tRNA genes for leucine and serine are doubled, but the tRNA genes for threonine and valine are not annotated. Phylogenetic analysis using the mitochondrial genomes of seven families of Sphaeropleales indicated that ITBB-AG6 is closely related to Monoraphidium neglectum, and falls in the family Selenastraceae with 100% bootstrap support. Two species in the family Neochloridaceae are separated by a species in Hydrodictyaceae, indicating a polyphyletic nature. These findings revealed the complicated phylogenetic relationships of the Sphaeropleales and the necessity of genome sequences in the taxonomy of microalgae.

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.

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Mitochondrial mRNA Processing in the Chlorophyte Alga Pediastrum duplex and Streptophyte Alga Chara vulgaris Reveals an Evolutionary Branch in Mitochondrial mRNA Processing

Mitochondria carry the remnant of an ancestral bacterial chromosome and express those genes with a system separate and distinct from the nucleus. Mitochondrial genes are transcribed as poly-cistronic primary transcripts which are post-transcriptionally processed to create individual translationally competent mRNAs. Algae post-transcriptional processing has only been explored in Chlamydomonas reinhardtii (Class: Chlorophyceae) and the mature mRNAs are different than higher plants, having no 5' UnTranslated Regions (UTRs), much shorter and more variable 3' UTRs and polycytidylated mature mRNAs. In this study, we analyzed transcript termini using circular RT-PCR and PacBio Iso-Seq to survey the 3' and 5' UTRs and termini for two green algae, Pediastrum duplex (Class: Chlorophyceae) and Chara vulgaris (Class: Charophyceae). This enabled the comparison of processing in the chlorophyte and charophyte clades of green algae to determine if the differences in mitochondrial mRNA processing pre-date the invasion of land by embryophytes. We report that the 5' mRNA termini and non-template 3' termini additions in P. duplex resemble those of C. reinhardtii, suggesting a conservation of mRNA processing among the chlorophyceae. We also report that C. vulgaris mRNA UTRs are much longer than chlorophytic examples, lack polycytidylation, and are polyadenylated similar to embryophytes. This demonstrates that some mitochondrial mRNA processing events diverged with the split between chlorophytic and streptophytic algae.

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.

Genome analyses provide insights into the evolution and adaptation of the eukaryotic Picophytoplankton Mychonastes homosphaera

Background

Picophytoplankton are abundant and can contribute greatly to primary production in eutrophic lakes. Mychonastes species are among the common eukaryotic picophytoplankton in eutrophic lakes. We used third-generation sequencing technology to sequence the whole genome of Mychonastes homosphaera isolated from Lake Chaohu, a eutrophic freshwater lake in China.

Result

The 24.23 Mbp nuclear genome of M.homosphaera, harboring 6649 protein-coding genes, is more compact than the genomes of the closely related Sphaeropleales species. This genome streamlining may be caused by a reduction in gene family number, intergenic size and introns. The genome sequence of M.homosphaera reveals the strategies adopted by this organism for environmental adaptation in the eutrophic lake. Analysis of cultures and the protein complement highlight the metabolic flexibility of M.homosphaera, the genome of which encodes genes involved in light harvesting, carbohydrate metabolism, and nitrogen and microelement metabolism, many of which form functional gene clusters. Reconstruction of the bioenergetic metabolic pathways of M.homosphaera, such as the lipid, starch and isoprenoid pathways, reveals characteristics that make this species suitable for biofuel production.

Conclusion

The analysis of the whole genome of M. homosphaera provides insights into the genome streamlining, the high lipid yield, the environmental adaptation and phytoplankton evolution.

Evolution of simple multicellular life cycles in dynamic environments

The mode of reproduction is a critical characteristic of any species, as it has a strong effect on its evolution. As any other trait, the reproduction mode is subject to natural selection and may adapt to the environment. When the environment varies over time, different reproduction modes could be optimal at different times. The natural response to a dynamic environment seems to be bet hedging, where multiple reproductive strategies are stochastically executed. Here, we develop a framework for the evolution of simple multicellular life cycles in a dynamic environment. We use a matrix population model of undifferentiated multicellular groups undergoing fragmentation and ask which mode maximizes the population growth rate. Counterintuitively, we find that natural selection in dynamic environments generally tends to promote deterministic, not stochastic, reproduction modes.

Phylogenetic Analysis and Substitution Rate Estimation of Colonial Volvocine Algae Based on Mitochondrial Genomes

We sequenced the mitochondrial genome of six colonial volvocine algae, namely: Pandorina morum, Pandorina colemaniae, Volvulina compacta, Colemanosphaera angeleri, Colemanosphaera charkowiensi, and Yamagishiella unicocca. Previous studies have typically reconstructed the phylogenetic relationship between colonial volvocine algae based on chloroplast or nuclear genes. Here, we explore the validity of phylogenetic analysis based on mitochondrial protein-coding genes. We found phylogenetic incongruence of the genera Yamagishiella and Colemanosphaera. In Yamagishiella, the stochastic error and linkage group formed by the mitochondrial protein-coding genes prevent phylogenetic analyses from reflecting the true relationship. In Colemanosphaera, a different reconstruction approach revealed a different phylogenetic relationship. This incongruence may be because of the influence of biological factors, such as incomplete lineage sorting or horizontal gene transfer. We also analyzed the substitution rates in the mitochondrial and chloroplast genomes between colonial volvocine algae. Our results showed that all volvocine species showed significantly higher substitution rates for the mitochondrial genome compared with the chloroplast genome. The nonsynonymous substitution (dN)/synonymous substitution (dS) ratio is similar in the genomes of both organelles in most volvocine species, suggesting that the two counterparts are under a similar selection pressure. We also identified a few chloroplast protein-coding genes that showed high dN/dS ratios in some species, resulting in a significant dN/dS ratio difference between the mitochondrial and chloroplast genomes.

Complete mitogenomes of the chlorophycean green algae Bulbochaete rectangularis var. hiloensis (Oedogoniales) and Stigeoclonium helveticum (Chaetophorales) provide insight into the sequence of events that led to the acquisition of a reduced-derived pattern of evolution in the Chlamydomonadales and Sphaeropleales

Mitogenome evolution in the Chlorophyceae is characterized by the acquisition of a reduced-derived pattern by the Chlamydomonadales + Sphaeropleales clade. Because no mitogenomes are available for the sister clade Oedogoniales + Chaetophorales + Chaetopeltidales, it remains unclear whether the common ancestor of chlorophycean green algae harbored a reduced-derived or ancestral-type mitogenome. The 70,191 and 46,765-bp mitogenomes reported here for Bulbochaete rectangularis var. hiloensis (Oedogoniales) and Stigeoclonium helveticum (Chaetophorales), respectively, shed light on this question. Both contain the same set of 41 conserved genes, a repertoire lacking numerous protein-coding genes but featuring all 27 tRNA genes typically found in ancestral-type mitogenomes.

Assessing Combinability of Phylogenomic Data Using Bayes Factors

With the rapid reduction in sequencing costs of high-throughput genomic data, it has become commonplace to use hundreds of genes to infer phylogeny of any study system. While sampling a large number of genes has given us a tremendous opportunity to uncover previously unknown relationships and improve phylogenetic resolution, it also presents us with new challenges when the phylogenetic signal is confused by differences in the evolutionary histories of sampled genes. Given the incorporation of accurate marginal likelihood estimation methods into popular Bayesian software programs, it is natural to consider using the Bayes Factor (BF) to compare different partition models in which genes within any given partition subset share both tree topology and edge lengths. We explore using marginal likelihood to assess data subset combinability when data subsets have varying levels of phylogenetic discordance due to deep coalescence events among genes (simulated within a species tree), and compare the results with our recently described phylogenetic informational dissonance index (D) estimated for each data set. BF effectively detects phylogenetic incongruence and provides a way to assess the statistical significance of D values. We use BFs to assess data combinability using an empirical data set comprising 56 plastid genes from the green algal order Volvocales. We also discuss the potential need for calibrating BFs and demonstrate that BFs used in this study are correctly calibrated.

The inflated mitochondrial genomes of siphonous green algae reflect processes driving expansion of noncoding DNA and proliferation of introns

Within the siphonous green algal order Bryopsidales, the size and gene arrangement of chloroplast genomes has been examined extensively, while mitochondrial genomes have been mostly overlooked. The recently published mitochondrial genome of Caulerpa lentillifera is large with expanded noncoding DNA, but it remains unclear if this is characteristic of the entire order. Our study aims to evaluate the evolutionary forces shaping organelle genome dynamics in the Bryopsidales based on the C. lentillifera and Ostreobium quekettii mitochondrial genomes. In this study, the mitochondrial genome of O. quekettii was characterised using a combination of long and short read sequencing, and bioinformatic tools for annotation and sequence analyses. We compared the mitochondrial and chloroplast genomes of O. quekettii and C. lentillifera to examine hypotheses related to genome evolution. The O. quekettii mitochondrial genome is the largest green algal mitochondrial genome sequenced (241,739 bp), considerably larger than its chloroplast genome. As with the mtDNA of C. lentillifera, most of this excess size is from the expansion of intergenic DNA and proliferation of introns. Inflated mitochondrial genomes in the Bryopsidales suggest effective population size, recombination and/or mutation rate, influenced by nuclear-encoded proteins, differ between the genomes of mitochondria and chloroplasts, reducing the strength of selection to influence evolution of their mitochondrial genomes.

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

Evolutionary Analysis of Unicellular Species in Chlamydomonadales Through Chloroplast Genome Comparison With the Colonial Volvocine Algae

This study is the first determination of six chloroplast genomes of colonial volvocine algae, Colemanosphaera charkowiensis, Volvulina compacta, Pandorina colemaniae, Pandorina morum, Colemanosphaera angeleri, and Yamagishiella unicocca. Based on 55 chloroplast protein-coding genes, we compared the nonsynonymous (dN) and synonymous (dS) substitution rates between colonial volvocine algae and the other unicellular Chlamydomonadales species. When refer to the dN, we found 27 genes were significantly different, among them, 19 genes were significant higher in unicellular species (FDR-adjusted P < 0.05). When refer to the dS, we found 10 genes were significantly different, among them, 6 genes were significant higher in unicellular species (FDR-adjusted P < 0.05). Then we identified 14 putative fast-evolving genes and 11 putative positively selected genes of unicellular species, we analyzed the function of positively selected sites of the overlap genes of putative fast-evolving and positively selected genes, and found some sites were close to the important functional region of the proteins. Photosynthesis is the process to transform and store solar energy by chloroplast, it plays a vital role in the survival of algae, this study is the first to use the chloroplast genomes to analysis the evolutionary relationship between colonial and unicellular species in Chlamydomonadales. We found more genes have higher substitution rates in unicellular species and proposed that the fast-evolving and positively selected two genes, psbA and psbC, may help to improve the photosynthetic efficiency of unicellular species in Chlamydomonadales.

The 4-Celled Tetrabaena socialis Nuclear Genome Reveals the Essential Components for Genetic Control of Cell Number at the Origin of Multicellularity in the Volvocine Lineage

Multicellularity is the premier example of a major evolutionary transition in individuality and was a foundational event in the evolution of macroscopic biodiversity. The volvocine chlorophyte lineage is well suited for studying this process. Extant members span unicellular, simple colonial, and obligate multicellular taxa with germ-soma differentiation. Here, we report the nuclear genome sequence of one of the most morphologically simple organisms in this lineage-the 4-celled colonial Tetrabaena socialis and compare this to the three other complete volvocine nuclear genomes. Using conservative estimates of gene family expansions a minimal set of expanded gene families was identified that associate with the origin of multicellularity. These families are rich in genes related to developmental processes. A subset of these families is lineage specific, which suggests that at a genomic level the evolution of multicellularity also includes lineage-specific molecular developments. Multiple points of evidence associate modifications to the ubiquitin proteasomal pathway (UPP) with the beginning of coloniality. Genes undergoing positive or accelerating selection in the multicellular volvocines were found to be enriched in components of the UPP and gene families gained at the origin of multicellularity include components of the UPP. A defining feature of colonial/multicellular life cycles is the genetic control of cell number. The genomic data presented here, which includes diversification of cell cycle genes and modifications to the UPP, align the genetic components with the evolution of this trait.

Enlarged and highly repetitive plastome of Lagarostrobos and plastid phylogenomics of Podocarpaceae

Podocarpaceae is the largest family in cupressophytes (conifers II), but its plastid genomes (plastomes) are poorly studied, with plastome data currently existing for only four of the 19 Podocarpaceous genera. In this study, we sequenced and assembled the complete plastomes from representatives of eight additional genera, including Afrocarpus, Dacrydium, Lagarostrobos, Lepidothamnus, Pherosphaera, Phyllocladus, Prumnopitys, and Saxegothaea. We found that Lagarostrobos, a monotypic genus native to Tasmania, has the largest plastome (151,496 bp) among any cupressophytes studied to date. Plastome enlargement in Lagarostrobos coincides with increased intergenic spacers, repeats, and duplicated genes. Among the Podocarpaceae, Lagarostrobos has the most rearranged plastome, but its substitution rates are modest. Plastid phylogenomic analyses based on 81 plastid genes clarify the positions of previously conflicting Podocarpaceous genera. Tree topologies firmly support the division of Podocarpaceae into two sister clades: (1) the Prumnopityoid clade and (2) the clade containing Podocarpoid, Dacrydioid, Pherosphaera, and Saxegothaea. The Phyllocladus is nested within the Podocarpaceae, thus familial status of the monotypic Phyllocladaceae is not supported.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

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.

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.

Repeated evolution and reversibility of self-fertilization in the volvocine green algae

Outcrossing and self-fertilization are fundamental strategies of sexual reproduction, each with different evolutionary costs and benefits. Self-fertilization is thought to be an evolutionary "dead-end" strategy, beneficial in the short term but costly in the long term, resulting in self-fertilizing species that occupy only the tips of phylogenetic trees. Here, we use volvocine green algae to investigate the evolution of self-fertilization. We use ancestral-state reconstructions to show that self-fertilization has repeatedly evolved from outcrossing ancestors and that multiple reversals from selfing to outcrossing have occurred. We use three phylogenetic metrics to show that self-fertilization is not restricted to the tips of the phylogenetic tree, a finding inconsistent with the view of self-fertilization as a dead-end strategy. We also find no evidence for higher extinction rates or lower speciation rates in selfing lineages. We find that self-fertilizing species have significantly larger colonies than outcrossing species, suggesting the benefits of selfing may counteract the costs of increased size. We speculate that our macroevolutionary results on self-fertilization (i.e., non-tippy distribution, no decreased diversification rates) may be explained by the haploid-dominant life cycle that occurs in volvocine algae, which may alter the costs and benefits of selfing.

Multiple Independent Changes in Mitochondrial Genome Conformation in Chlamydomonadalean Algae

Chlamydomonadalean green algae are no stranger to linear mitochondrial genomes, particularly members of the Reinhardtinia clade. At least nine different Reinhardtinia species are known to have linear mitochondrial DNAs (mtDNAs), including the model species Chlamydomonas reinhardtii. Thus, it is no surprise that some have suggested that the most recent common ancestor of the Reinhardtinia clade had a linear mtDNA. But the recent uncovering of circular-mapping mtDNAs in a range of Reinhardtinia algae, such as Volvox carteri and Tetrabaena socialis, has shed doubt on this hypothesis. Here, we explore mtDNA sequence and structure within the colonial Reinhardtinia algae Yamagishiella unicocca and Eudorina sp. NIES-3984, which occupy phylogenetically intermediate positions between species with opposing mtDNA mapping structures. Sequencing and gel electrophoresis data indicate that Y. unicocca has a linear monomeric mitochondrial genome with long (3 kb) palindromic telomeres. Conversely, the mtDNA of Eudorina sp., despite having an identical gene order to that of Y. unicocca, assembled as a circular-mapping molecule. Restriction digests of Eudorina sp. mtDNA supported its circular map, but also revealed a linear monomeric form with a matching architecture and gene order to the Y. unicocca mtDNA. Based on these data, we suggest that there have been at least three separate shifts in mtDNA conformation in the Reinhardtinia, and that the common ancestor of this clade had a linear monomeric mitochondrial genome with palindromic telomeres.

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.

The mutational hazard hypothesis of organelle genome evolution: 10 years on

Why is there such a large variation in size and noncoding DNA content among organelle genomes? One explanation is that this genomic variation results from differences in the rates of organelle mutation and random genetic drift, as opposed to being the direct product of natural selection. Along these lines, the mutational hazard hypothesis (MHH) holds that 'excess' DNA is a mutational liability (because it increases the potential for harmful mutations) and, thus, has a greater tendency to accumulate in an organelle system with a low mutation rate as opposed to one with a high rate of mutation. Various studies have explored this hypothesis and, more generally, the relationship between organelle genome architecture and the mode and efficiency of organelle DNA repair. Although some of these investigations are in agreement with the MHH, others have contradicted it; nevertheless, they support a central role of mutation, DNA maintenance pathways and random genetic drift in fashioning organelle chromosomes. Arguably, one of the most important contributions of the MHH is that it has sparked crucial, widespread discussions about the importance of nonadaptive processes in genome evolution.

The Gonium pectorale genome demonstrates co-option of cell cycle regulation during the evolution of multicellularity

The transition to multicellularity has occurred numerous times in all domains of life, yet its initial steps are poorly understood. The volvocine green algae are a tractable system for understanding the genetic basis of multicellularity including the initial formation of cooperative cell groups. Here we report the genome sequence of the undifferentiated colonial alga, Gonium pectorale, where group formation evolved by co-option of the retinoblastoma cell cycle regulatory pathway. Significantly, expression of the Gonium retinoblastoma cell cycle regulator in unicellular Chlamydomonas causes it to become colonial. The presence of these changes in undifferentiated Gonium indicates extensive group-level adaptation during the initial step in the evolution of multicellularity. These results emphasize an early and formative step in the evolution of multicellularity, the evolution of cell cycle regulation, one that may shed light on the evolutionary history of other multicellular innovations and evolutionary transitions.

Phylogenetic patterns of gene rearrangements in four mitochondrial genomes from the green algal family Hydrodictyaceae (Sphaeropleales, Chlorophyceae)

Background

The variability in gene organization and architecture of green algal mitochondrial genomes is only recently being studied on a finer taxonomic scale. Sequenced mt genomes from the chlorophycean orders Volvocales and Sphaeropleales exhibit considerable variation in size, content, and structure, even among closely related genera. However, sampling of mt genomes on a within-family scale is still poor and the sparsity of information precludes a thorough understanding of genome evolution in the green algae.

Methods

Genomic DNA of representative taxa were sequenced on an Illumina HiSeq2500 to produce 2x100 bp paired reads, and mitochondrial genomes were assembled and annotated using Geneious v.6.1.5. Phylogenetic analysis of 13 protein-coding mitochondrial genes spanning the Sphaeropleales was performed.

Results

This study presents one of the first within-family comparisons of mt genome diversity, and is the first to report complete mt genomes for the family Hydrodictyaceae (order Sphaeropleales). Four complete mt genomes representing three taxa and four phylogenetic groups, Stauridium tetras, Pseudopediastrum boryanum, and Pediastrum duplex, range in size from 37,723 to 53,560 bp. The size variability is primarily due to intergenic region expansion, and intron content is generally low compared with other mt genomes of Sphaeropleales.

Conclusions

Certain gene rearrangements appear to follow a phylogenetic pattern, and with a more thorough taxon sampling genome-level sequence may be useful in resolving systematic conundrums that plague this morphologically diverse family.

Electronic supplementary material

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

Gamete Dialogs in Green Lineages

Gamete fusion is a core process of sexual reproduction and, in both plants and animals, different sex gametes fuse within species. Although most of the molecular factors involved in gamete interaction are still unknown in various sex-possessing eukaryotes, reports of such factors in algae and land plants have been increasing in the past decade. In particular, knowledge of gamete interaction in flowering plants and green algae has increased since the identification of the conserved gamete fusion factor generative cell specific 1/hapless 2 (GCS1/HAP2). GCS1 was first identified as a pollen generative cell-specific transmembrane protein in the lily (Lilium longiflorum), and was then shown to function not only in flowering plant gamete fusion but also in various eukaryotes, including unicellular protists and metazoans. In addition, although initially restricted to Chlamydomonas, knowledge of gamete attachment in flowering plants was also acquired. This review focuses on recent progress in the study of gamete interaction in volvocine green algae and flowering plants and discusses conserved mechanisms of gamete recognition, attachment, and fusion leading to zygote formation.

Highly Conserved Mitochondrial Genomes among Multicellular Red Algae of the Florideophyceae

Two red algal classes, the Florideophyceae (approximately 7,100 spp.) and Bangiophyceae (approximately 193 spp.), comprise 98% of red algal diversity in marine and freshwater habitats. These two classes form well-supported monophyletic groups in most phylogenetic analyses. Nonetheless, the interordinal relationships remain largely unresolved, in particular in the largest subclass Rhodymeniophycidae that includes 70% of all species. To elucidate red algal phylogenetic relationships and study organelle evolution, we determined the sequence of 11 mitochondrial genomes (mtDNA) from 5 florideophycean subclasses. These mtDNAs were combined with existing data, resulting in a database of 25 florideophytes and 12 bangiophytes (including cyanidiophycean species). A concatenated alignment of mt proteins was used to resolve ordinal relationships in the Rhodymeniophycidae. Red algal mtDNA genome comparisons showed 47 instances of gene rearrangement including 12 that distinguish Bangiophyceae from Hildenbrandiophycidae, and 5 that distinguish Hildenbrandiophycidae from Nemaliophycidae. These organelle data support a rapid radiation and surprisingly high conservation of mtDNA gene syntheny among the morphologically divergent multicellular lineages of Rhodymeniophycidae. In contrast, we find extensive mitochondrial gene rearrangements when comparing Bangiophyceae and Florideophyceae and multiple examples of gene loss among the different red algal lineages.

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.

Genome expansion via lineage splitting and genome reduction in the cicada endosymbiont Hodgkinia

Comparative genomics from mitochondria, plastids, and mutualistic endosymbiotic bacteria has shown that the stable establishment of a bacterium in a host cell results in genome reduction. Although many highly reduced genomes from endosymbiotic bacteria are stable in gene content and genome structure, organelle genomes are sometimes characterized by dramatic structural diversity. Previous results from Candidatus Hodgkinia cicadicola, an endosymbiont of cicadas, revealed that some lineages of this bacterium had split into two new cytologically distinct yet genetically interdependent species. It was hypothesized that the long life cycle of cicadas in part enabled this unusual lineage-splitting event. Here we test this hypothesis by investigating the structure of the Ca. Hodgkinia genome in one of the longest-lived cicadas, Magicicada tredecim. We show that the Ca. Hodgkinia genome from M. tredecim has fragmented into multiple new chromosomes or genomes, with at least some remaining partitioned into discrete cells. We also show that this lineage-splitting process has resulted in a complex of Ca. Hodgkinia genomes that are 1.1-Mb pairs in length when considered together, an almost 10-fold increase in size from the hypothetical single-genome ancestor. These results parallel some examples of genome fragmentation and expansion in organelles, although the mechanisms that give rise to these extreme genome instabilities are likely different.

When the lights go out: the evolutionary fate of free-living colorless green algae

The endosymbiotic origin of plastids was a launching point for eukaryotic evolution. The autotrophic abilities bestowed by plastids are responsible for much of the eukaryotic diversity we observe today. But despite its many advantages, photosynthesis has been lost numerous times and in disparate lineages throughout eukaryote evolution. For example, among green algae, several groups have lost photosynthesis independently and in response to different selective pressures; these include the parasitic/pathogenic trebouxiophyte genera Helicosporidium and Prototheca, and the free-living chlamydomonadalean genera Polytomella and Polytoma. Here, we examine the published data on colorless green algae and argue that investigations into the different evolutionary routes leading to their current nonphotosynthetic lifestyles provide exceptional opportunities to understand the ecological and genomic factors involved in the loss of photosynthesis.

Massive and widespread organelle genomic expansion in the green algal genus Dunaliella

The mitochondrial genomes of chlamydomonadalean green algae are renowned for their highly reduced and conserved gene repertoires, which are almost fixed at 12 genes across the entire lineage. The sizes of these genomes, however, are much more variable, with some species having small, compact mitochondrial DNAs (mtDNAs) and others having expanded ones. Earlier work demonstrated that the halophilic genus Dunaliella contains extremely inflated organelle genomes, but to date the mtDNA of only one isolate has been explored. Here, by surveying mtDNA architecture across the Chlamydomonadales, we show that various Dunaliella species have undergone massive levels of mitochondrial genomic expansion, harboring the most inflated, intron-dense mtDNAs available from chlorophyte green algae. The same also appears to be true for their plastid genomes, which are potentially among the largest of all plastid-containing eukaryotes. Genetic divergence data are used to investigate the underlying causes of such extreme organelle genomic architectures, and ultimately reveal order-of-magnitude differences in mitochondrial versus plastid mutation rates within Dunaliella.

Sympatric speciation in a bacterial endosymbiont results in two genomes with the functionality of one

Mutualisms that become evolutionarily stable give rise to organismal interdependencies. Some insects have developed intracellular associations with communities of bacteria, where the interdependencies are manifest in patterns of complementary gene loss and retention among members of the symbiosis. Here, using comparative genomics and microscopy, we show that a three-member symbiotic community has become a four-way assemblage through a novel bacterial lineage-splitting event. In some but not all cicada species of the genus Tettigades, the endosymbiont Candidatus Hodgkinia cicadicola has split into two new cytologically distinct but metabolically interdependent species. Although these new bacterial genomes are partitioned into discrete cell types, the intergenome patterns of gene loss and retention are almost perfectly complementary. These results defy easy classification: they show genomic patterns consistent with those observed after both speciation and whole-genome duplication. We suggest that our results highlight the potential power of nonadaptive forces in shaping organismal complexity.

Gene arrangement convergence, diverse intron content, and genetic code modifications in mitochondrial genomes of sphaeropleales (chlorophyta)

The majority of our knowledge about mitochondrial genomes of Viridiplantae comes from land plants, but much less is known about their green algal relatives. In the green algal order Sphaeropleales (Chlorophyta), only one representative mitochondrial genome is currently available—that of Acutodesmus obliquus. Our study adds nine completely sequenced and three partially sequenced mitochondrial genomes spanning the phylogenetic diversity of Sphaeropleales. We show not only a size range of 25–53 kb and variation in intron content (0–11) and gene order but also conservation of 13 core respiratory genes and fragmented ribosomal RNA genes. We also report an unusual case of gene arrangement convergence in Neochloris aquatica, where the two rns fragments were secondarily placed in close proximity. Finally, we report the unprecedented usage of UCG as stop codon in Pseudomuriella schumacherensis. In addition, phylogenetic analyses of the mitochondrial protein-coding genes yield a fully resolved, well-supported phylogeny, showing promise for addressing systematic challenges in green algae.