Inflated organelle genomes and a circular-mapping mtDNA probably existed at the origin of coloniality in volvocine green algae

The volvocine lineage is a monophyletic grouping of unicellular, colonial and multicellular algae, and a model for studying the evolution of multicellularity. In addition to being morphologically diverse, volvocine algae boast a surprising amount of organelle genomic variation. Moreover, volvocine organelle genome complexity appears to scale positively with organismal complexity. However, the organelle DNA architecture at the origin of colonial living is not known. To examine this issue, we sequenced the plastid and mitochondrial DNAs (ptDNA and mtDNA) of the 4-celled alga Tetrabaena socialis, which is basal to the colonial and multicellular volvocines.

Tetrabaena socialis has a circular-mapping mitochondrial genome, contrasting with the linear mtDNA architecture of its relative Chlamydomonas reinhardtii. This suggests that a circular-mapping mtDNA conformation emerged at or near the transition to group living in the volvocines, or represents the ancestral state of the lineage as a whole. The T. socialis ptDNA is very large (>405 kb) and dense with repeats, supporting the idea that a shift from a unicellular to a colonial existence coincided with organelle genomic expansion, potentially as a result of increased random genetic drift. These data reinforce the idea that volvocine algae harbour some of the most expanded plastid chromosomes from the eukaryotic tree of life. Circular-mapping mtDNAs are turning out to be more common within volvocines than originally thought, particularly for colonial and multicellular species. Altogether, volvocine organelle genomes became markedly more inflated during the evolution of multicellularity, but complex organelle genomes appear to have existed at the very beginning of colonial living.     

Unparalleled GC content in the plastid DNA of Selaginella

One of the more conspicuous features of plastid DNA (ptDNA) is its low guanine and cytosine (GC) content. As of February 2009, all completely-sequenced plastid genomes have a GC content below 43% except for the ptDNA of the lycophyte Selaginella uncinata, which is 55% GC. The forces driving the S. uncinata ptDNA towards G and C are undetermined, and it is unknown if other Selaginella species have GC-biased plastid genomes. This study presents the complete ptDNA sequence of Selaginella moellendorffii and compares it with the previously reported S. uncinata plastid genome. Partial ptDNA sequences from 103 different Selaginella species are also described as well as a significant proportion of the S. moellendorffii mitochondrial genome. Moreover, S. moellendorffii express sequence tags are data-mined to estimate levels of plastid and mitochondrial RNA editing. Overall, these data are used to show that: (1) there is a genus-wide GC bias in Selaginella ptDNA, which is most pronounced in South American articulate species; (2) within the Lycopsida class (and among plants in general), GC-biased ptDNA is restricted to the Selaginella genus; (3) the cause of this GC bias is arguably a combination of reduced AT-mutation pressure relative to other plastid genomes and a large number of C-to-U RNA editing sites; and (4) the mitochondrial DNA (mtDNA) of S. moellendorffii is also GC biased (even more so than the ptDNA) and is arguably the most GC-rich organelle genome observed to date—the high GC content of the mtDNA also appears to be influenced by RNA editing. Ultimately, these findings provide convincing support for the earlier proposed theory that the GC content of land-plant organelle DNA is positively correlated and directly connected to levels of organelle RNA editing.     

Two Distinct Plastid Genome Configurations and Unprecedented Intraspecies Length Variation in the accD Coding Region in Medicago truncatula

We fully sequenced four and partially sequenced six additional plastid genomes of the model legume Medicago truncatula. Three accessions, Jemalong 2HA, Borung and Paraggio, belong to ssp. truncatula, and R108 to ssp. tricycla. We report here that the R108 ptDNA has a ∼45-kb inversion compared with the ptDNA in ssp. truncatula, mediated by a short, imperfect repeat. DNA gel blot analyses of seven additional ssp. tricycla accessions detected only one of the two alternative genome arrangements, represented by three and four accessions each. Furthermore, we found a variable number of repeats in the essential accD and ycf1 coding regions. The repeats within accD are recombinationally active, yielding variable-length insertions and deletions in the central part of the coding region. The length of ACCD was distinct in each of the 10 sequenced ecotypes, ranging between 650 and 796 amino acids. The repeats in the ycf1 coding region are also recombinationally active, yielding short indels in 10 regions of the reading frames. Thus, the plastid genome variability we report here could be linked to repeat-mediated genome rearrangements. However, the rate of recombination was sufficiently low, so that no heterogeneity of ptDNA could be observed in populations maintained by single-seed descent.     

High Sequence Divergence but Limited Architectural Rearrangements in Organelle Genomes of Cyanophora (Glaucophyta) Species

Cyanophora is the glaucophyte model taxon. Following the sequencing of the nuclear genome of C. paradoxa, studies based on single organelle and nuclear molecular markers revealed previously unrecognized species diversity within this glaucophyte genus. Here, we present the complete plastid (ptDNA) and mitochondrial (mtDNA) genomes of C. kugrensii, C. sudae, and C. biloba. The respective sizes and coding capacities of both ptDNAs and mtDNAs are conserved among Cyanophora species with only minor differences due to specific gene duplications. Organelle phylogenomic analyses consistently recover the species C. kugrensii and C. paradoxa as a clade and C. sudae and C. biloba as a separate group. The phylogenetic affiliations of the four Cyanophora species are consistent with architectural similarities shared at the organelle genomic level. Genetic distance estimations from both organelle sequences are also consistent with phylogenetic and architecture evidence. Comparative analyses confirm that the Cyanophora mitochondrial genes accumulate substitutions at 3-fold higher rates than plastid counterparts, suggesting that mtDNA markers are more appropriate to investigate glaucophyte diversity and evolutionary events that occur at a population level. The study of complete organelle genomes is becoming the standard for species delimitation and is particularly relevant to study cryptic diversity in microbial groups.     

SecA is plastid-encoded in a red alga: implications for the evolution of plastid genomes and the thylakoid protein import apparatus

Summary Partial sequence analysis of the plastid DNA (ptDNA) from a red alga, Antithamnion sp., revealed the presence of a homologue to the Escherichia coli SecA gene as well as two open reading frames (ORF 510, ORF 179). In addition a sec Y homologue has been detected on the plastid genome by heterologous hybridization. None of these genes has been found in completely sequenced chlorophytic plastid genomes. SecA and secY gene copies were also detected in the ptDNA of a chromophytic alga, indicating that secAY may be ubiquitous in rhodophytes and chromophytes. The significance of these findings for the evolution of plastid genomes and the thylakoid protein import mechanism is discussed.     

Nucleotide Diversity of the Colorless Green Alga Polytomella parva (Chlorophyceae,Chlorophyta): High for the Mitochondrial Telomeres,Surprisingly Low Everywhere Else*

ABSTRACT. Silent‐site nucleotide diversity data (π_silent) can provide insights into the forces driving genome evolution. Here we present π_silent statistics for the mitochondrial and nuclear DNAs of Polytomella parva, a nonphotosynthetic green alga with a highly reduced, linear fragmented mitochondrial genome. We show that this species harbors very little genetic diversity, with the exception of the mitochondrial telomeres, which have an excess of polymorphic sites. These data are compared with previously published π_silent values from the mitochondrial and nuclear genomes of the model species Chlamydomonas reinhardtii and Volvox carteri, which are close relatives of P. parva, and are used to understand the modes and tempos of genome evolution within green algae.     

Effective Extraction and Assembly Methods for Simultaneously Obtaining Plastid and Mitochondrial Genomes

Background

In conventional approaches to plastid and mitochondrial genome sequencing, the sequencing steps are performed separately; thus, plastid DNA (ptDNA) and mitochondrial DNA (mtDNA) should be prepared independently. However, it is difficult to extract pure ptDNA and mtDNA from plant tissue. Following the development of high-throughput sequencing technology, many researchers have attempted to obtain plastid genomes or mitochondrial genomes using high-throughput sequencing data from total DNA. Unfortunately, the huge datasets generated consume massive computing and storage resources and cost a great deal, and even more importantly, excessive pollution reads affect the accuracy of the assembly. Therefore, it is necessary to develop an effective method that can generate base sequences from plant tissue and that is suitable for all plant species. Here, we describe a highly effective, low-cost method for obtaining plastid and mitochondrial genomes simultaneously.

Results

First, we obtained high-quality DNA employing Partial Concentration Extraction. Second, we evaluated the purity of the DNA sample and determined the sequencing dataset size employing Vector Control Quantitative Analysis. Third, paired-end reads were obtained using a high-throughput sequencing platform. Fourth, we obtained scaffolds employing Two-step Assembly. Finally, we filled in gaps using specific methods and obtained complete plastid and mitochondrial genomes. To ensure the accuracy of plastid and mitochondrial genomes, we validated the assembly using PCR and Sanger sequencing. Using this method,we obtained complete plastid and mitochondrial genomes with lengths of 153,533 nt and 223,412 nt separately.

Conclusion

A simple method for extracting, evaluating, sequencing and assembling plastid and mitochondrial genomes was developed. This method has many advantages: it is timesaving, inexpensive and reproducible and produces high-quality sequence. Furthermore, this method can produce plastid and mitochondrial genomes simultaneously and be used for other plant species. Due to its simplicity and extensive applicability, this method will support research on plant cytoplasmic genomes.     

Patterns of Genomic Integration of Nuclear Chloroplast DNA Fragments in Plant Species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.     

Plastid Genome of <Emphasis Type="Italic">Dictyopteris divaricata</Emphasis> (Dictyotales,Phaeophyceae): Understanding the Evolution of Plastid Genomes in Brown Algae

Dictyotophycidae is a subclass of brown algae containing 395 species that are distributed worldwide. A complete plastid (chloroplast) genome (ptDNA or cpDNA) had not previously been sequenced from this group. In this study, the complete plastid genome of Dictyopteris divaricata (Okamura) Okamura (Dictyotales, Phaeophyceae) was characterized and compared to other brown algal ptDNAs. This plastid genome was 126,099 bp in size with two inverted repeats (IRs) of 6026 bp. The D. divaricata IRs contained rpl21, making its IRs larger than representatives from the orders Fucales and Laminariales, but was smaller than that from Ectocarpales. The G + C content of D. divaricata (31.19%) was the highest of the known ptDNAs of brown algae (28.94–31.05%). Two protein-coding genes, rbcR and rpl32, were present in ptDNAs of Laminariales, Ectocarpales (Ectocarpus siliculosus), and Fucales (LEF) but were absent in D. divaricata. Reduced intergenic space (13.11%) and eight pairs of overlapping genes in D. divaricata ptDNA made it the most compact plastid genome in brown algae so far. The architecture of D. divaricata ptDNA showed higher similarity to that of Laminariales compared with Fucales and Ectocarpales. The difference in general features, gene content, and architecture among the ptDNAs of D. divaricata and LEF clade revealed the diversity and evolutionary trends of plastid genomes in brown algae.     

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

The mitochondrial genome of grape (Vitis vinifera), the largestorganelle genome sequenced so far, is presented. The genomeis 773,279 nt long and has the highest coding capacity amongknown angiosperm mitochondrial DNAs (mtDNAs). The proportionof promiscuous DNA of plastid origin in the genome is also thelargest ever reported for an angiosperm mtDNA, both in absoluteand relative terms. In all, 42.4% of chloroplast genome of Vitishas been incorporated into its mitochondrial genome. In orderto test if horizontal gene transfer (HGT) has also contributedto the gene content of the grape mtDNA, we built phylogenetictrees with the coding sequences of mitochondrial genes of grapeand their homologs from plant mitochondrial genomes. Many incongruentgene tree topologies were obtained. However, the extent of incongruencebetween these gene trees is not significantly greater than thatobserved among optimal trees for chloroplast genes, the commonancestry of which has never been in doubt. In both cases, weattribute this incongruence to artifacts of tree reconstruction,insufficient numbers of characters, and gene paralogy. Thisfinding leads us to question the recent phylogenetic interpretationof Bergthorsson et al. (2003, 2004) and Richardson and Palmer(2007) that rampant HGT into the mtDNA of Amborella best explainsphylogenetic incongruence between mitochondrial gene trees forangiosperms. The only evidence for HGT into the Vitis mtDNAfound involves fragments of two coding sequences stemming fromtwo closteroviruses that cause the leaf roll disease of thisplant. We also report that analysis of sequences shared by bothchloroplast and mitochondrial genomes provides evidence fora previously unknown gene transfer route from the mitochondrionto the chloroplast.     

Mitochondrial Genome of Phlebia radiata Is the Second Largest (156 kbp) among Fungi and Features Signs of Genome Flexibility and Recent Recombination Events

Mitochondria are eukaryotic organelles supporting individual life-style via generation of proton motive force and cellular energy, and indispensable metabolic pathways. As part of genome sequencing of the white rot Basidiomycota species Phlebia radiata, we first assembled its mitochondrial genome (mtDNA). So far, the 156 348 bp mtDNA is the second largest described for fungi, and of considerable size among eukaryotes. The P. radiata mtDNA assembled as single circular dsDNA molecule containing genes for the large and small ribosomal RNAs, 28 transfer RNAs, and over 100 open reading frames encoding the 14 fungal conserved protein subunits of the mitochondrial complexes I, III, IV, and V. Two genes (atp6 and tRNA-Ile^GAU) were duplicated within 6.1 kbp inverted region, which is a unique feature of the genome. The large mtDNA size, however, is explained by the dominance of intronic and intergenic regions (sum 80% of mtDNA sequence). The intergenic DNA stretches harness short (≤200 nt) repetitive, dispersed and overlapping sequence elements in abundance. Long self-splicing introns of types I and II interrupt eleven of the conserved genes (cox1,2,3; cob; nad1,2,4,4L,5; rnl; rns). The introns embrace a total of 57 homing endonucleases with LAGLIDADGD and GYI-YIG core motifs, which makes P. radiata mtDNA to one of the largest known reservoirs of intron-homing endonucleases. The inverted duplication, intergenic stretches, and intronic features are indications of dynamics and genetic flexibility of the mtDNA, not fully recognized to this extent in fungal mitochondrial genomes previously, thus giving new insights for the evolution of organelle genomes in eukaryotes.     

MITOCHONDRIAL GENOME CONFORMATION AMONG CW‐GROUP CHLOROPHYCEAN ALGAE1

Most green algal taxa have circular‐mapping mitochondrial genomes, whereas some have linear genome‐ or subgenomic‐sized mitochondrial DNAs (mtDNA). It is not clear, however, if the circular‐mapping genomes represent genome‐sized circular molecules, if such circular molecules and the linear forms are the predominant in vivo mtDNA structures, or if the linear forms arose only once or multiple times among extant green algal lineages. We therefore examined the DNA components detected with homologous mtDNA probes after pulsed‐field gel electrophoresis of total cellular DNA from the chlorophycean basal bodies displaced clockwise(CW)‐group taxa Chlamydomonas reinhardtii and Chlamydomonas moewusii. For C. reinhardtii, the 15.8‐kb linear mtDNA was the only DNA component detected, and there was no evidence of circular or large linear precursors of this DNA. In the case of C. moewusii, which is known to have a circular‐mapping 22.9‐kb mitochondrial genome, three DNA components were detected; these appeared to be circular (relaxed and supercoiled) and genome‐sized linear DNA molecules, the latter of which likely resulted from random double‐strand breaks in the circular forms during DNA isolation. In further studies, DNA from additional CW‐group taxa was examined using conventional gel electrophoresis and DNA‐filter blot analysis with C. reinhardtii and C. moewusii mtDNA probes. We conclude that all taxa from the “Volvox clade” (sensu Nakayama et al. 1996 of the CW‐group have genome‐ or subgenomic‐sized linear mtDNAs as their predominant mtDNA form and that these arose from a genome‐sized circular form in an ancestor that existed near the base of this clade.     

The <Emphasis Type="Italic">Dunaliella salina</Emphasis> organelle genomes: large sequences,inflated with intronic and intergenic DNA

Background  

Dunaliella salina Teodoresco, a unicellular, halophilic green alga belonging to the Chlorophyceae, is among the most industrially important microalgae. This is because D. salina can produce massive amounts of β-carotene, which can be collected for commercial purposes, and because of its potential as a feedstock for biofuels production. Although the biochemistry and physiology of D. salina have been studied in great detail, virtually nothing is known about the genomes it carries, especially those within its mitochondrion and plastid. This study presents the complete mitochondrial and plastid genome sequences of D. salina and compares them with those of the model green algae Chlamydomonas reinhardtii and Volvox carteri.     

The Complete Sequence of the Mitochondrial Genome of Butomus umbellatus – A Member of an Early Branching Lineage of Monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.     

Genome skimming by shotgun sequencing helps resolve the phylogeny of a pantropical tree family

Whole genome sequencing is helping generate robust phylogenetic hypotheses for a range of taxonomic groups that were previously recalcitrant to classical molecular phylogenetic approaches. As a case study, we performed a shallow shotgun sequencing of eight species in the tropical tree family Chrysobalanaceae to retrieve large fragments of high‐copy number DNA regions and test the potential of these regions for phylogeny reconstruction. We were able to assemble the nuclear ribosomal cluster (nrDNA), the complete plastid genome (ptDNA) and a large fraction of the mitochondrial genome (mtDNA) with approximately 1000×, 450× and 120× sequencing depth respectively. The phylogenetic tree obtained with ptDNA resolved five of the seven internal nodes. In contrast, the tree obtained with mtDNA and nrDNA data were largely unresolved. This study demonstrates that genome skimming is a cost‐effective approach and shows potential in plant molecular systematics within Chrysobalanaceae and other under‐studied groups.     

Comparative genomics of large mitochondria in placozoans

The first sequenced mitochondrial genome of a placozoan, Trichoplax adhaerens, challenged the conventional wisdom that a compact mitochondrial genome is a common feature among all animals. Three additional placozoan mitochondrial genomes representing highly divergent clades have been sequenced to determine whether the large Trichoplax mtDNA is a shared feature among members of the phylum Placozoa or a uniquely derived condition. All three mitochondrial genomes were found to be very large, 32- to 37-kb, circular molecules, having the typical 12 respiratory chain genes, 24 tRNAs, rnS, and rnL. They share with the Trichoplax mitochondrial genome the absence of atp8, atp9, and all ribosomal protein genes, the presence of several cox1 introns, and a large open reading frame containing an intron group I LAGLIDADG endonuclease domain. The differences in mtDNA size within Placozoa are due to variation in intergenic spacer regions and the presence or absence of long open reading frames of unknown function. Phylogenetic analyses of the 12 respiratory chain genes support the monophyly of Placozoa. The similarities in composition and structure between the three mitochondrial genomes reported here and that of Trichoplax's mtDNA suggest that their uncompacted state is a shared ancestral feature to other nonmetazoans while their gene content is a derived feature shared only among the Metazoa.     

Nucleotide diversity of the Chlamydomonas reinhardtii plastid genome: addressing the mutational-hazard hypothesis

Background  

The mutational-hazard hypothesis argues that the noncoding-DNA content of a genome is a consequence of the mutation rate (μ) and the effective number of genes per locus in the population (N _g ). The hypothesis predicts that genomes with a high N _g μ will be more compact than those with a small N _g μ. Approximations of N _g μ can be gained by measuring the nucleotide diversity at silent sites (π_silent). We addressed the mutation-hazard hypothesis apropos plastid-genome evolution by measuring π_silent of the Chlamydomonas reinhardtii plastid DNA (ptDNA), the most noncoding-DNA-dense plastid genome observed to date. The data presented here in conjunction with previously published values of π_silent for the C. reinhardtii mitochondrial and nuclear genomes, which are respectively compact and bloated, allow for a complete analysis of nucleotide diversity and genome compactness in all three genetic compartments of this model organism.     

Identification of a highly conserved hydroxyproline-rich glycoprotein in the cell walls of Chlamydomonas reinhardtii and two other Volvocales

The unicellular alga Chlamydomonas reinhardtii Dang, has a cell wall made entirely from hydroxyproline-rich glycoproteins (HRGPs). We recently employed a quantiative in vitro reconstitution system (Adair et al. 1987, J. Cell Biol. 105, 2373–2382) to assign outer-wall HRGPs of C. reinhardtii to specific sublayers, and describe the major interactions responsible for their assembly. Some of these interactions appear to involve relatively conserved HRGP domains, as evidenced by interspecific cell-wall reconstitution between C. reinhardtii and two multicellular Volvocales (Volvoxcarteri lyengar and Gonium pectorale Müller). In the present report we provide biochemical and immunological evidence that the outer cell-walls of V. carteri and G. pectorale both contain prominent HRGPs closely related to C. reinhardtii GP2. Identification of conserved GP2 homologues indicates a molecular basis for interspecific reconstitution and provides a useful avenue for characterization of HRGP domains mediating cell-wall formation in these algae.Abbreviations GP1, 2, 3 outer-cell wall glycoproteins 1, 2, and 3 - GP2_dg deglycosylated GP2 - HRGP hydroxyprolinerich glycoprotein - SDS-PAGE sodium docecyl sulfate polyacrylamide gel electrophoresis     

The GC-rich mitochondrial and plastid genomes of the green alga Coccomyxa give insight into the evolution of organelle DNA nucleotide landscape

Most of the available mitochondrial and plastid genome sequences are biased towards adenine and thymine (AT) over guanine and cytosine (GC). Examples of GC-rich organelle DNAs are limited to a small but eclectic list of species, including certain green algae. Here, to gain insight in the evolution of organelle nucleotide landscape, we present the GC-rich mitochondrial and plastid DNAs from the trebouxiophyte green alga Coccomyxa sp. C-169. We compare these sequences with other GC-rich organelle DNAs and argue that the forces biasing them towards G and C are nonadaptive and linked to the metabolic and/or life history features of this species. The Coccomyxa organelle genomes are also used for phylogenetic analyses, which highlight the complexities in trying to resolve the interrelationships among the core chlorophyte green algae, but ultimately favour a sister relationship between the Ulvophyceae and Chlorophyceae, with the Trebouxiophyceae branching at the base of the chlorophyte crown.