Foreign Plastid Sequences in Plant Mitochondria are Frequently Acquired Via Mitochondrion-to-Mitochondrion Horizontal Transfer

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Comparative genomics and phylogenetic analysis of mitochondrial genomes of Neocinnamomum

BACKGROUND

Neocinnamomum plants are considered a promising feedstock for biodiesel in China, due to the richness in long-chain fatty acids (LCFAs) found in their seeds. However, the mitochondrial genome (mitogenome) of this genus has not yet been systematically described, and the exploration of species relationships within this genus using mitogenome sequences is also an uncharted territory. This has hindered our understanding of mitogenome diversity and the evolutionary relationships within Neocinnamomum.

RESULTS

In this study, a total of 24 individuals representing seven distinct taxa from the genus Neocinnamomum were subjected to Illumina sequencing, and the species N. delavayi was sequenced using Oxford Nanopore sequencing technology. We successfully assembled the mitogenome of N. delavayi, which is 778,066 bp in size and exhibits a single circular structure. The analysis identified 659 dispersed repeats, 211 simple sequence repeats (SSRs), and 30 tandem repeats within the mitogenome. Additionally, 37 homologous fragments, totaling 9929 bp, were found between the mitogenome and the plastid genome (plastome). The codons of 41 protein-coding genes (PCGs) had a preference for ending in A/T, and the codon usage bias of the majority of these genes was influenced by natural selection pressures. Comparative genomic analysis revealed low collinearity and significant gene rearrangements between species. Phylogenetic analysis resulted in the classification of Neocinnamomum into six distinct clades, contradicting previous findings which based on complete plastomes and nuclear ribosomal cistron (nrDNA). In the PCGs of 24 individuals, 86 mutation events were identified, which included three indels and 83 SNPs. Notably, the ccmC gene underwent positive selection in pairwise comparisons of three species pairs. Furthermore, 748 RNA editing sites were predicted within the PCGs of the N. delavayi mitogenome.

CONCLUSIONS

This study enriches our knowledge of the mitogenomes in the family Lauraceae, and provides valuable data and a foundation for genomic evolution research, genetic resource conservation, and molecular breeding in Neocinnamomum.

Features and evolutionary adaptations of the mitochondrial genome of Garuga forrestii W. W. Sm

Introduction

Garuga forrestii W. W. Sm. is a tree species of the Burseraceae family, endemic to China, found in hot/warm-dry valleys. This species plays a crucial role in maintaining biodiversity in these ecosystems.

Methods

We performed de novo assembly of the Garuga forrestii mitochondrial genome using PMAT (v.1.5.4), resulting in a typical circular molecule of 606,853 bp. The genome consists of 31 tRNA genes, 3 rRNA genes, 35 protein-coding genes, and 1 pseudogene. The study also investigates RNA editing sites and evolutionary patterns.

Results

The mitochondrial genome exhibits a low proportion of repetitive sequences (3.30%), suggesting a highly conserved structure. A high copy number of the trnM-CAT gene (4 copies) is noted, which may contribute to genomic rearrangement and adaptive evolution. Among the 476 RNA editing sites, hydrophilic-hydrophobic and hydrophobic-hydrophobic editing events are most common, accounting for 77.10%. Negative selection predominates among most genes (Ka/Ks < 1), while a few genes (e.g., matR, nad3, rps1, rps12, and rps4) show signs of positive selection (Ka/Ks > 1), potentially conferring evolutionary advantages. Additionally, a significant A/T bias is observed at the third codon position. Phylogenomic analysis supports the APG IV classification, with no evidence of horizontal gene transfer.

Discussion

This mitochondrial genome offers valuable insights into the adaptive mechanisms and evolutionary processes of Garuga forrestii. It enhances our understanding of the species' biogeography in tropical Southeast Asia and Southwest China, providing key information on the evolutionary history of this genus.

Array of metabolic pathways in a kleptoplastidic foraminiferan protist supports chemoautotrophy in dark, euxinic seafloor sediments

Investigations of the metabolic capabilities of anaerobic protists advances our understanding of the evolution of eukaryotic life on Earth and for uncovering analogous extraterrestrial complex microbial life. Certain species of foraminiferan protists live in environments analogous to early Earth conditions when eukaryotes evolved, including sulfidic, anoxic and hypoxic sediment porewaters. Foraminifera are known to form symbioses as well as to harbor organelles from other eukaryotes (chloroplasts), possibly bolstering the host's independence from oxygen. The full extent of foraminiferal physiological capabilities is not fully understood. To date, evidence for foraminiferal anaerobiosis was gleaned from specimens first subjected to stresses associated with removal from in situ conditions. Here, we report comprehensive gene expression analysis of benthic foraminiferal populations preserved in situ on the euxinic (anoxic and sulfidic) bathyal seafloor, thus avoiding environmental alterations associated with sample recovery, including pressure reduction, sunlight exposure, warming, and oxygenation. Metatranscriptomics, metagenome-assembled genomes, and measurements of substrate uptake were used to study the kleptoplastidic foraminifer Nonionella stella inhabiting sulfur-oxidizing bacterial mats of the Santa Barbara Basin, off California. We show N. stella energy generation under dark euxinia is unusual because it orchestrates complex metabolic pathways for ATP production and carbon fixation through the Calvin cycle. These pathways include extended glycolysis, anaerobic fermentation, sulfide oxidation, and the presence of a membrane-bound inorganic pyrophosphatase, an enzyme that hydrolyzes inorganic pyrophosphate to actively pump protons across the mitochondrial membrane.

Comprehensive analysis of the first complete mitogenome and plastome of a traditional Chinese medicine Viola diffusa

BACKGROUND

Viola diffusa is used in the formulation of various Traditional Chinese Medicines (TCMs), including antiviral, antimicrobial, antitussive, and anti-inflammatory drugs, due to its richness in flavonoids and triterpenoids. The biosynthesis of these compounds is largely mediated by cytochrome P450 enzymes, which are primarily located in the membranes of mitochondria and the endoplasmic reticulum.

RESULTS

This study presents the complete assembly of the mitogenome and plastome of Viola diffusa. The circular mitogenome spans 474,721 bp with a GC content of 44.17% and encodes 36 unique protein-coding genes, 21 tRNA, and 3 rRNA. Except for the RSCU values of 1 observed for the start codon (AUG) and tryptophan (UGG), the mitochondrial protein-coding genes exhibited a codon usage bias, with most estimates deviating from 1, similar to patterns observed in closely related species. Analysis of repetitive sequences in the mitogenome demonstrated potential homologous recombination mediated by these repeats. Sequence transfer analysis revealed 24 homologous sequences shared between the mitogenome and plastome, including nine full-length genes. Collinearity was observed among Viola diffusa species within the other members of Malpighiales order, indicated by the presence of homologous fragments. The length and arrangement of collinear blocks varied, and the mitogenome exhibited a high frequency of gene rearrangement.

CONCLUSIONS

We present the first complete assembly of the mitogenome and plastome of Viola diffusa, highlighting its implications for pharmacological, evolutionary, and taxonomic studies. Our research underscores the multifaceted importance of comprehensive mitogenome analysis.

Methods for assembling complex mitochondrial genomes in land plants

The large size and complex structural rearrangements inherent in the mitochondrial genomes of land plants pose challenges for their sequencing. Originally, the assembly of these genomes required the cloning of mitochondrial DNA fragments followed by Sanger sequencing. Subsequently, the advent of next-generation sequencing significantly expedited the process. This review highlights examples of plant mitochondrial genome assembly employing various technologies, including 454 sequencing, Illumina short sequencing reads, and Pacific Biosciences or Oxford Nanopore Technology long sequencing reads. The combination of short and long reads in hybrid assembly has proven to be the most efficient approach for achieving reliable assemblies of land plant mitochondrial genomes.

New Insights into the Diversity of Mitochondrial Plastid DNA

The mitochondrial plastid DNAs (MTPTs) in seed plants were reported more than 40 years ago and exhibited a high diversity regarding gene content, quantity, and size. However, the mechanism that resulted in the current diversity of MTPTs in angiosperms has not been fully discovered. In this study, we sequenced and characterized the complete organelle genomes of Limonia acidissima L., a monotypic species of Rutaceae. The newly generated and previously published organelle genomes of 42 species were used to explore the diversity of MTPTs regarding quantity, gene content, size, and coverage of chloroplast genome (cpDNA) regions. The results showed that the number of MTPTs ranged from three to 74, of which the lengths were from 100 to 53,731 bp. The highest coverage of MTPTs was found in the inverted repeat region, whereas the small single repeat region had the lowest coverage. Based on the previous data and current results, we propose a scenario for the diversity of MTPTs in angiosperms. In the first stage, the whole cpDNA might migrate to the mitogenome. Then, different genomic events, such as duplication, deletion, substitution, and inversion, have occurred continuously and independently and resulted in extremely variable profiles of mitogenomes among angiosperms. Our hypothesis provides a new and possibly reliable scenario for explaining the present circumstances of MTPTs in angiosperms. However, more genomic data should be mined, and more studies should be conducted to clarify this natural phenomenon in plants.

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes

BACKGROUND

Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus.

RESULTS

We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT.

CONCLUSIONS

Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

Organelle Genomes of Epipogium roseum Provide Insight into the Evolution of Mycoheterotrophic Orchids

Epipogium roseum, commonly known as one of the ghost orchids due to its rarity and almost transparent color, is a non-photosynthetic and fully mycoheterotrophic plant. Given its special nutritional strategies and evolutionary significance, the mitogenome was first characterized, and three plastomes sampled from Asia were assembled. The plastomes were found to be the smallest among Orchidaceae, with lengths ranging from 18,339 to 19,047 bp, and exhibited high sequence variety. For the mitogenome, a total of 414,552 bp in length, comprising 26 circular chromosomes, were identified. A total of 54 genes, including 38 protein-coding genes, 13 tRNA genes, and 3 rRNA genes, were annotated. Multiple repeat sequences spanning a length of 203,423 bp (45.47%) were discovered. Intriguingly, six plastid regions via intracellular gene transfer and four plastid regions via horizontal gene transfer to the mitogenome were observed. The phylogenomics, incorporating 90 plastomes and 56 mitogenomes, consistently revealed the sister relationship of Epipogium and Gastrodia, with a bootstrap percentage of 100%. These findings shed light on the organelle evolution of Orchidaceae and non-photosynthetic plants.

Phylogenomics and topological conflicts in the tribe Anthospermeae (Rubiaceae)

Genome skimming (shallow whole-genome sequencing) offers time- and cost-efficient production of large amounts of DNA data that can be used to address unsolved evolutionary questions. Here we address phylogenetic relationships and topological incongruence in the tribe Anthospermeae (Rubiaceae), using phylogenomic data from the mitochondrion, the nuclear ribosomal cistron, and the plastome. All three genomic compartments resolve relationships in the Anthospermeae; the tribe is monophyletic and consists of three major subclades. Carpacoce Sond. is sister to the remaining clade, which comprises an African subclade and a Pacific subclade. Most results, from all three genomic compartments, are statistically well supported; however, not fully consistent. Intergenomic topological incongruence is most notable in the Pacific subclade but present also in the African subclade. Hybridization and introgression followed by organelle capture may explain these conflicts but other processes, such as incomplete lineage sorting (ILS), can yield similar patterns and cannot be ruled out based on the results. Whereas the null hypothesis of congruence among all sequenced loci in the individual genomes could not be rejected for nuclear and mitochondrial data, it was rejected for plastid data. Phylogenetic analyses of three subsets of plastid loci identified using the hierarchical likelihood ratio test demonstrated statistically supported intragenomic topological incongruence. Given that plastid genes are thought to be fully linked, this result is surprising and may suggest modeling or sampling error. However, biological processes such as biparental inheritance and inter-plastome recombination have been reported and may be responsible for the observed intragenomic incongruence. Mitochondrial insertions into the plastome are rarely documented in angiosperms. Our results indicate that a mitochondrial insertion event in the plastid trnS GGA - rps4 IGS region occurred in the common ancestor of the Pacific clade of Anthospermeae. Exclusion/inclusion of this locus in phylogenetic analyses had a strong impact on topological results in the Pacific clade.

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Comparative analysis of the complete mitogenome of Geoffroea decorticans: a native tree surviving in the Atacama Desert

Chañar (Geoffroea decorticans (Gill., ex Hook. & Arn.) Burkart) has been highly significant for indigenous people in the Atacama Desert for over 3,000 years. Through evolutionary processes, the G. decorticans mitogenome likely underwent changes facilitating its adaptation to the extreme conditions of the Atacama Desert. Here, we compare the mitochondrial genome of G. decorticans with those of other Papilionoideae family species. The complete mitogenome of G. decorticans was sequenced and assembled, making it the first in the genus Geoffroea. The mitogenome contained 383,963 base pairs, consisting of 33 protein coding genes, 21 transfer RNA genes, and 3 ribosomal RNA genes. The Chañar mitogenome is relatively compact, and has two intact genes (sdh4 and nad1) which were not observed in most other species. Additionally, Chañar possessed the highest amount of mitochondrial DNA of plastid origin among angiosperm species. The phylogenetic analysis of the mitogenomes of Chañar and 12 other taxa displayed a high level of consistency in taxonomic classification, when compared to those of the plastid genome. Atp8 was subjected to positive selection, while the ccmFc and rps1 were subjected to neutral selection. This study provides valuable information regarding its ability to survive the extreme environmental conditions of the Atacama Desert.

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

Factors contributing to mitogenome size variation and a recurrent intracellular DNA transfer in Melastoma

BACKGROUND

Mitogenome sizes of seed plants vary substantially even among closely related species, which are often related to horizontal or intracellular DNA transfer (HDT or IDT) events. However, the mechanisms of this size variation have not been well characterized.

RESULTS

Here we assembled and characterized the mitogenomes of three species of Melastoma, a tropical shrub genus experiencing rapid speciation. The mitogenomes of M. candidum (Mc), M. sanguineum (Ms) and M. dodecandrum (Md) were assembled to a circular mapping chromosome of 391,595 bp, 395,542 bp and 412,026 bp, respectively. While the mitogenomes of Mc and Ms showed good collinearity except for a large inversion of ~ 150 kb, there were many rearrangements in the mitogenomes between Md and either Mc or Ms. Most non-alignable sequences (> 80%) between Mc and Ms are from gain or loss of mitochondrial sequences. Whereas, between Md and either Mc or Ms, non-alignable sequences in Md are mainly chloroplast derived sequences (> 30%) and from putative horizontal DNA transfers (> 30%), and those in both Mc and Ms are from gain or loss of mitochondrial sequences (> 80%). We also identified a recurrent IDT event in another congeneric species, M. penicillatum, which has not been fixed as it is only found in one of the three examined populations.

CONCLUSIONS

By characterizing mitochondrial genome sequences of Melastoma, our study not only helps understand mitogenome size evolution in closely related species, but also cautions different evolutionary histories of mitochondrial regions due to potential recurrent IDT events in some populations or species.

REFMAKER: make your own reference to target nuclear loci in low coverage genome skimming libraries. Phylogenomic application in Sapotaceae

Genome skimming approach is widely used in plant systematics to infer phylogenies mostly from organelle genomes. However, organelles represent only 10% of the produced libraries, and the low coverage associated with these libraries (< 3X) prevents the capture of nuclear sequences, which are not always available in non-model organisms or limited to the ribosomal regions. We developed REFMAKER, a user-friendly pipeline, to create specific sets of nuclear loci that can next be extracted directly from the genome skimming libraries. For this, a catalogue is built from the meta-assembly of each library contigs and cleaned by selecting the nuclear regions and removing duplicates from clustering steps. Libraries are next mapped onto this catalogue and consensus sequences are generated to produce a ready-to-use phylogenetic matrix following different filtering parameters aiming at removing putative errors and paralogous sequences. REFMAKER allowed us to infer a well resolved phylogeny in Capurodendron (Sapotaceae) on 67 nuclear loci from low-coverage libraries (<1X). The resulting phylogeny is concomitant with one previously inferred on 638 nuclear genes from target enrichment libraries. While it remains preliminary because of this low sequencing depth, REFMAKER therefore opens perspectives in phylogenomics by allowing nuclear phylogeny reconstructions with genome skimming datasets.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

Multichromosomal Mitochondrial Genome of Paphiopedilum micranthum: Compact and Fragmented Genome, and Rampant Intracellular Gene Transfer

Orchidaceae is one of the largest families of angiosperms. Considering the large number of species in this family and its symbiotic relationship with fungi, Orchidaceae provide an ideal model to study the evolution of plant mitogenomes. However, to date, there is only one draft mitochondrial genome of this family available. Here, we present a fully assembled and annotated sequence of the mitochondrial genome (mitogenome) of Paphiopedilum micranthum, a species with high economic and ornamental value. The mitogenome of P. micranthum was 447,368 bp in length and comprised 26 circular subgenomes ranging in size from 5973 bp to 32,281 bp. The genome encoded for 39 mitochondrial-origin, protein-coding genes; 16 tRNAs (three of plastome origin); three rRNAs; and 16 ORFs, while rpl10 and sdh3 were lost from the mitogenome. Moreover, interorganellar DNA transfer was identified in 14 of the 26 chromosomes. These plastid-derived DNA fragments represented 28.32% (46,273 bp) of the P. micranthum plastome, including 12 intact plastome origin genes. Remarkably, the mitogenome of P. micranthum and Gastrodia elata shared 18% (about 81 kb) of their mitochondrial DNA sequences. Additionally, we found a positive correlation between repeat length and recombination frequency. The mitogenome of P. micranthum had more compact and fragmented chromosomes compared to other species with multichromosomal structures. We suggest that repeat-mediated homologous recombination enables the dynamic structure of mitochondrial genomes in Orchidaceae.

The mitochondrial genome sequence of Abies alba Mill. reveals a high structural and combinatorial variation

BACKGROUND

Plant mitogenomes vary widely in size and genomic architecture. Although hundreds of plant mitogenomes of angiosperm species have already been sequence-characterized, only a few mitogenomes are available from gymnosperms. Silver fir (Abies alba) is an economically important gymnosperm species that is widely distributed in Europe and occupies a large range of environmental conditions. Reference sequences of the nuclear and chloroplast genome of A. alba are available, however, the mitogenome has not yet been assembled and studied.

RESULTS

Here, we used paired-end Illumina short reads generated from a single haploid megagametophyte in combination with PacBio long reads from high molecular weight DNA of needles to assemble the first mitogenome sequence of A. alba. Assembly and scaffolding resulted in 11 mitogenome scaffolds, with the largest scaffold being 0.25 Mbp long. Two of the scaffolds displayed a potential circular structure supported by PCR. The total size of the A. alba mitogenome was estimated at 1.43 Mbp, similar to the size (1.33 Mbp) of a draft assembly of the Abies firma mitogenome. In total, 53 distinct genes of known function were annotated in the A. alba mitogenome, comprising 41 protein-coding genes, nine tRNA, and three rRNA genes. The proportion of highly repetitive elements (REs) was 0.168. The mitogenome seems to have a complex and dynamic structure featured by high combinatorial variation, which was specifically confirmed by PCR for the contig with the highest mapping coverage. Comparative analysis of all sequenced mitogenomes of gymnosperms revealed a moderate, but significant positive correlation between mitogenome size and proportion of REs.

CONCLUSIONS

The A. alba mitogenome provides a basis for new comparative studies and will allow to answer important structural, phylogenetic and other evolutionary questions. Future long-read sequencing with higher coverage of the A. alba mitogenome will be the key to further resolve its physical structure. The observed positive correlation between mitogenome size and proportion of REs will be further validated once available mitogenomes of gymnosperms would become more numerous. To test whether a higher proportion of REs in a mitogenome leads to an increased recombination and higher structural complexity and variability is a prospective avenue for future research.

Comparative analyses of three complete Primula mitogenomes with insights into mitogenome size variation in Ericales

BACKGROUND

Although knowledge of the sizes, contents, and forms of plant mitochondrial genomes (mitogenomes) is increasing, little is known about the mechanisms underlying their structural diversity. Evolutionary information on the mitogenomes of Primula, an important ornamental taxon, is more limited than the information on their nuclear and plastid counterparts, which has hindered the comprehensive understanding of Primula mitogenomic diversity and evolution. The present study reported and compared three Primula mitogenomes and discussed the size expansion of mitogenomes in Ericales.

RESULTS

Mitogenome master circles were sequenced and successfully assembled for three Primula taxa and were compared with publicly available Ericales mitogenomes. The three mitogenomes contained similar gene contents and varied primarily in their structures. The Primula mitogenomes possessed relatively high nucleotide diversity among all examined plant lineages. In addition, high nucleotide diversity was found among Primula species between the Mediterranean and Himalaya-Hengduan Mountains. Most predicted RNA editing sites appeared in the second amino acid codon, increasing the hydrophobic character of the protein. An early stop in atp6 caused by RNA editing was conserved across all examined Ericales species. The interfamilial relationships within Ericales and interspecific relationships within Primula could be well resolved based on mitochondrial data. Transfer of the two longest mitochondrial plastid sequences (MTPTs) occurred before the divergence of Primula and its close relatives, and multiple independent transfers could also occur in a single MTPT sequence. Foreign sequence [MTPTs and mitochondrial nuclear DNA sequences (NUMTs)] uptake and repeats were to some extent associated with changes in Ericales mitogenome size, although none of these relationships were significant overall.

CONCLUSIONS

The present study revealed relatively conserved gene contents, gene clusters, RNA editing, and MTPTs but considerable structural variation in Primula mitogenomes. Relatively high nucleotide diversity was found in the Primula mitogenomes. In addition, mitogenomic genes, collinear gene clusters, and locally collinear blocks (LCBs) all showed phylogenetic signals. The evolutionary history of MTPTs in Primula was complicated, even in a single MTPT sequence. Various reasons for the size variation observed in Ericales mitogenomes were found.

The Mitogenome of Sedum plumbizincicola (Crassulaceae): Insights into RNA Editing, Lateral Gene Transfer, and Phylogenetic Implications

As the largest family within the order Saxifragales, Crassulaceae contains about 34 genera with 1400 species. Mitochondria play a critical role in cellular energy production. Since the first land plant mitogenome was reported in Arabidopsis, more than 400 mitogenomic sequences have been deposited in a public database. However, no entire mitogenome data have been available for species of Crassulaceae to date. To better understand the evolutionary history of the organelles of Crassulaceae, we sequenced and performed comprehensive analyses on the mitogenome of Sedum plumbizincicola. The master mitogenomic circle is 212,159 bp in length, including 31 protein-coding genes (PCGs), 14 tRNA genes, and 3 rRNA genes. We further identified totally 508 RNA editing sites in PCGs, and demonstrated that the second codon positions of mitochondrial genes are most prone to RNA editing events. Notably, by neutrality plot analyses, we observed that the mitochondrial RNA editing events have large effects on the driving forces of plant evolution. Additionally, 4 MTPTs and 686 NUMTs were detected in the mitochondrial and nuclear genomes of S. plumbizincicola, respectively. Additionally, we conducted further analyses on gene transfer, secondary structures of mitochondrial RNAs, and phylogenetic implications. Therefore, the findings presented here will be helpful for future investigations on plant mitogenomes.

Analysis of the complete plastomes and nuclear ribosomal DNAs from Euonymus hamiltonianus and its relatives sheds light on their diversity and evolution

Euonymus hamiltonianus and its relatives (Celastraceae family) are used for ornamental and medicinal purposes. However, species identification in Euonymus is difficult due to their morphological diversity. Using plastid genome (plastome) data, we attempt to reveal phylogenetic relationship among Euonymus species and develop useful markers for molecular identification. We assembled the plastome and nuclear ribosomal DNA (nrDNA) sequences from five Euonymus lines collected from South Korea: three Euonymus hamiltonianus accessions, E. europaeus, and E. japonicus. We conducted an in-depth comparative analysis using ten plastomes, including other publicly available plastome data for this genus. The genome structures, gene contents, and gene orders were similar in all Euonymus plastomes in this study. Analysis of nucleotide diversity revealed six divergence hotspots in their plastomes. We identified 339 single nucleotide polymorphisms and 293 insertion or deletions among the four E. hamiltonianus plastomes, pointing to abundant diversity even within the same species. Among 77 commonly shared genes, 9 and 33 were identified as conserved genes in the genus Euonymus and E. hamiltonianus, respectively. Phylogenetic analysis based on plastome and nrDNA sequences revealed the overall consensus and relationships between plastomes and nrDNAs. Finally, we developed six barcoding markers and successfully applied them to 31 E. hamiltonianus lines collected from South Korea. Our findings provide the molecular basis for the classification and molecular taxonomic criteria for the genus Euonymus (at least in Korea), which should aid in more objective classification within this genus. Moreover, the newly developed markers will be useful for understanding the species delimitation of E. hamiltonianus and closely related species.

Mitochondrial genome recombination in somatic hybrids of Solanum commersonii and S. tuberosum

Interspecific somatic hybridization has been performed in potato breeding experiments to increase plant resistance against biotic and abiotic stress conditions. We analyzed the mitochondrial and plastid genomes and 45S nuclear ribosomal DNA (45S rDNA) for the cultivated potato (S. tuberosum, St), wild potato (S. commersonii, Sc), and their somatic hybrid (StSc). Complex genome components and structure, such as the hybrid form of 45S rDNA in StSc, unique plastome in Sc, and recombinant mitogenome were identified. However, the mitogenome exhibited dynamic multipartite structures in both species as well as in the somatic hybrid. In St, the mitogenome is 756,058 bp and is composed of five subgenomes ranging from 297,014 to 49,171 bp. In Sc, it is 552,103 bp long and is composed of two sub-genomes of 338,427 and 213,676 bp length. StSc has 447,645 bp long mitogenome with two subgenomes of length 398,439 and 49,206 bp. The mitogenome structure exhibited dynamic recombination mediated by tandem repeats; however, it contained highly conserved genes in the three species. Among the 35 protein-coding genes of the StSc mitogenome, 21 were identical for all the three species, and 12 and 2 were unique in Sc and St, respectively. The recombinant mitogenome might be derived from homologous recombination between both species during somatic hybrid development.

Breaking the limits - multichromosomal structure of an early eudicot Pulsatilla patens mitogenome reveals extensive RNA-editing, longest repeats and chloroplast derived regions among sequenced land plant mitogenomes

BACKGROUND

The mitogenomes of vascular plants are one of the most structurally diverse molecules. In the present study we characterize mitogenomes of a rare and endangered species Pulsatilla patens. We investigated the gene content and its RNA editing potential, repeats distribution and plastid derived sequences.

RESULTS

The mitogenome structure of early divergent eudicot, endangered Pulsatilla patens does not support the master chromosome hypothesis, revealing the presence of three linear chromosomes of total length 986 613 bp. The molecules are shaped by the presence of extremely long, exceeding 87 kbp repeats and multiple chloroplast-derived regions including nearly complete inverted repeat. Since the plastid IR content of Ranunculales is very characteristic, the incorporation into mitogenome could be explained rather by intracellular transfer than mitochondrial HGT. The mitogenome contains almost a complete set of genes known from other vascular plants with exception of rps10 and sdh3, the latter being present but pseudogenized. Analysis of long ORFs enabled the identification of genes which are rarely present in plant mitogenomes, including RNA and DNA polymerases, albeit their presence even at species level is variable. Mitochondrial transcripts of P. patens were edited with a high frequency, which exceeded the level known in other analyzed angiosperms, despite the strict qualification criteria of counting the editing events and taking into analysis generally less frequently edited leaf transcriptome. The total number of edited sites was 902 and nad4 was identified as the most heavily edited gene with 65 C to U changes. Non-canonical, reverse U to C editing was not detected. Comparative analysis of mitochondrial genes of three Pulsatilla species revealed a level of variation comparable to chloroplast CDS dataset and much higher infrageneric differentiation than in other known angiosperm genera. The variation found in CDS of mitochondrial genes is comparable to values found among Pulsatilla plastomes. Despite the complicated mitogenome structure, 14 single copy regions of 329 kbp, not splitted by repeats or plastid-derived sequences (MTPT), revealed the potential for phylogenetic, phylogeographic and population genetics studies by revealing intra- and interspecific collinearity.

CONCLUSIONS

This study provides valuable new information about mitochondrial genome of early divergent eudicots, Pulsatilla patens, revealed multi-chromosomal structure and shed new light on mitogenomics of early eudicots.

Unprecedented organelle genomic variations in morning glories reveal independent evolutionary scenarios of parasitic plants and the diversification of plant mitochondrial complexes

Background

The morning glories (Convolvulaceae) are distributed worldwide and produce economically important crops, medicinal herbs, and ornamentals. Members of this family are diverse in morphological characteristics and trophic modes, including the leafless parasitic Cuscuta (dodders). Organelle genomes were generally used for studying plant phylogeny and genomic variations. Notably, plastomes in parasitic plants always show non-canonical features, such as reduced size and accelerated rates. However, few organelle genomes of this group have been sequenced, hindering our understanding of their evolution, and dodder mitogenome in particular.

Results

We assembled 22 new mitogenomes and 12 new plastomes in Convolvulaceae. Alongside previously known ones, we totally analyzed organelle genomes of 23 species in the family. Our sampling includes 16 leafy autotrophic species and 7 leafless parasitic dodders, covering 8 of the 12 tribes. Both the plastid and mitochondrial genomes of these plants have encountered variations that were rarely observed in other angiosperms. All of the plastomes possessed atypical IR boundaries. Besides the gene and IR losses in dodders, some leafy species also showed gene and intron losses, duplications, structural variations, and insertions of foreign DNAs. The phylogeny reconstructed by plastid protein coding sequences confirmed the previous relationship of the tribes. However, the monophyly of ‘Merremieae’ and the sister group of Cuscuta remained uncertain. The mitogenome was significantly inflated in Cuscuta japonica, which has exceeded over 800 kb and integrated massive DNAs from other species. In other dodders, mitogenomes were maintained in small size, revealing divergent evolutionary strategies. Mutations unique to plants were detected in the mitochondrial gene ccmFc, which has broken into three fragments through gene fission and splicing shift. The unusual changes likely initially happened to the common ancestor of the family and were caused by a foreign insertion from rosids followed by double-strand breaks and imprecise DNA repairs. The coding regions of ccmFc expanded at both sides after the fission, which may have altered the protein structure.

Conclusions

Our family-scale analyses uncovered unusual scenarios for both organelle genomes in Convolvulaceae, especially in parasitic plants. The data provided valuable genetic resources for studying the evolution of Convolvulaceae and plant parasitism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01250-1.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes

Background

The intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

Results

The mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with 10 missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmF_C gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

Conclusions

The lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08105-z.

Plant mitochondria import DNA via alternative membrane complexes involving various VDAC isoforms

Mitochondria possess transport mechanisms for import of RNA and DNA. Based on import into isolated Solanum tuberosum mitochondria in the presence of competitors, inhibitors or effectors, we show that DNA fragments of different size classes are taken up into plant organelles through distinct channels. Alternative channels can also be activated according to the amount of DNA substrate of a given size class. Analyses of Arabidopsis thaliana knockout lines pointed out a differential involvement of individual voltage-dependent anion channel (VDAC) isoforms in the formation of alternative channels. We propose several outer and inner membrane proteins as VDAC partners in these pathways.

Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.

Complete Plastid and Mitochondrial Genomes of Aeginetia indica Reveal Intracellular Gene Transfer (IGT), Horizontal Gene Transfer (HGT), and Cytoplasmic Male Sterility (CMS)

Orobanchaceae have become a model group for studies on the evolution of parasitic flowering plants, and Aeginetia indica, a holoparasitic plant, is a member of this family. In this study, we assembled the complete chloroplast and mitochondrial genomes of A. indica. The chloroplast and mitochondrial genomes were 56,381 bp and 401,628 bp long, respectively. The chloroplast genome of A. indica shows massive plastid genes and the loss of one IR (inverted repeat). A comparison of the A. indica chloroplast genome sequence with that of a previous study demonstrated that the two chloroplast genomes encode a similar number of proteins (except atpH) but differ greatly in length. The A. indica mitochondrial genome has 53 genes, including 35 protein-coding genes (34 native mitochondrial genes and one chloroplast gene), 15 tRNA (11 native mitochondrial genes and four chloroplast genes) genes, and three rRNA genes. Evidence for intracellular gene transfer (IGT) and horizontal gene transfer (HGT) was obtained for plastid and mitochondrial genomes. ψndhB and ψcemA in the A. indica mitogenome were transferred from the plastid genome of A. indica. The atpH gene in the plastid of A. indica was transferred from another plastid angiosperm plastid and the atpI gene in mitogenome A. indica was transferred from a host plant like Miscanthus siensis. Cox2 (orf43) encodes proteins containing a membrane domain, making ORF (Open Reading Frame) the most likely candidate gene for CMS development in A. indica.

Comparative analysis of nuclear, chloroplast, and mitochondrial genomes of watermelon and melon provides evidence of gene transfer

During plant evolution, there is genetic communication between organelle and nuclear genomes. A comparative analysis was performed on the organelle and nuclear genomes of the watermelon and melon. In the watermelon, chloroplast-derived sequences accounted for 7.6% of the total length of the mitochondrial genome. In the melon, chloroplast-derived sequences accounted for approximately 2.73% of the total mitochondrial genome. In watermelon and melon, the chloroplast-derived small-fragment sequences are either a subset of large-fragment sequences or appeared multiple times in the mitochondrial genome, indicating that these fragments may have undergone multiple independent migration integrations or emerged in the mitochondrial genome after migration, replication, and reorganization. There was no evidence of migration from the mitochondria to chloroplast genome. A sequence with a total length of about 73 kb (47%) in the watermelon chloroplast genome was homologous to a sequence of about 313 kb in the nuclear genome. About 33% of sequences in the watermelon mitochondrial genome was homologous with a 260 kb sequence in the nuclear genome. A sequence with a total length of about 38 kb (25%) in the melon chloroplast genome was homologous with 461 sequences in the nuclear genome, with a total length of about 301 kb. A 3.4 Mb sequence in the nuclear genome was homologous with a melon mitochondrial sequence. These results indicate that, during the evolution of watermelon and melon, a large amount of genetic material was exchanged between the nuclear genome and the two organelle genomes in the cytoplasm.

Utility of the Mitochondrial Genome in Plant Taxonomic Studies

Size, structure, and sequence content lability of plant mitochondrial genome (mtDNA) across species has sharply limited its use in taxonomic studies. Historically, mtDNA variation has been first investigated with RFLPs, while the development of universal primers then allowed studying sequence polymorphisms within short genomic regions (<3 kb). The recent advent of NGS technologies now offers new opportunities by greatly facilitating the assembly of longer mtDNA regions, and even full mitogenomes. Phylogenetic works aiming at comparing signals from different genomic compartments (i.e., nucleus, chloroplast, and mitochondria) have been developed on a few plant lineages, and have been shown especially relevant in groups with contrasted inheritance of organelle genomes. This chapter first reviews the main characteristics of mtDNA and the application offered in taxonomic studies. It then presents tips for best sequencing protocol based on NGS data to be routinely used in mtDNA-based phylogenetic studies.

Adaptive innovation of green plants by horizontal gene transfer

Horizontal gene transfer (HGT) refers to the movement of genetic material between distinct species by means other than sexual reproduction. HGT has contributed tremendously to the genome plasticity and adaptive evolution of prokaryotes and certain unicellular eukaryotes. The evolution of green plants from chlorophyte algae to angiosperms and from water to land represents a process of adaptation to diverse environments, which has been facilitated by acquisition of genetic material from other organisms. In this article, we review the occurrence of HGT in major lineages of green plants, including chlorophyte and charophyte green algae, bryophytes, lycophytes, ferns, and seed plants. In addition, we discuss the significance of horizontally acquired genes in the adaptive innovations of green plants and their potential applications to crop breeding and improvement.

Siberian larch (Larix sibirica Ledeb.) mitochondrial genome assembled using both short and long nucleotide sequence reads is currently the largest known mitogenome

BACKGROUND

Plant mitochondrial genomes (mitogenomes) can be structurally complex while their size can vary from ~ 222 Kbp in Brassica napus to 11.3 Mbp in Silene conica. To date, in comparison with the number of plant species, only a few plant mitogenomes have been sequenced and released, particularly for conifers (the Pinaceae family). Conifers cover an ancient group of land plants that includes about 600 species, and which are of great ecological and economical value. Among them, Siberian larch (Larix sibirica Ledeb.) represents one of the keystone species in Siberian boreal forests. Yet, despite its importance for evolutionary and population studies, the mitogenome of Siberian larch has not yet been assembled and studied.

RESULTS

Two sources of DNA sequences were used to search for mitochondrial DNA (mtDNA) sequences: mtDNA enriched samples and nucleotide reads generated in the de novo whole genome sequencing project, respectively. The assembly of the Siberian larch mitogenome contained nine contigs, with the shortest and the largest contigs being 24,767 bp and 4,008,762 bp, respectively. The total size of the genome was estimated at 11.7 Mbp. In total, 40 protein-coding, 34 tRNA, and 3 rRNA genes and numerous repetitive elements (REs) were annotated in this mitogenome. In total, 864 C-to-U RNA editing sites were found for 38 out of 40 protein-coding genes. The immense size of this genome, currently the largest reported, can be partly explained by variable numbers of mobile genetic elements, and introns, but unlikely by plasmid-related sequences. We found few plasmid-like insertions representing only 0.11% of the entire Siberian larch mitogenome.

CONCLUSIONS

Our study showed that the size of the Siberian larch mitogenome is much larger than in other so far studied Gymnosperms, and in the same range as for the annual flowering plant Silene conica (11.3 Mbp). Similar to other species, the Siberian larch mitogenome contains relatively few genes, and despite its huge size, the repeated and low complexity regions cover only 14.46% of the mitogenome sequence.

Multichromosomal structure and foreign tracts in the Ombrophytum subterraneum (Balanophoraceae) mitochondrial genome

Horizontal gene transfer (HGT) is frequent in parasitic plant mitochondria as a result of vascular connections established in host-parasite relationships. Recent studies of the holoparasitic plant Lophophytum mirabile (Balanophoraceae) revealed the unprecedented acquisition of a large amount of mitochondrial sequences from its legume host. We focused on a close relative, the generalist holoparasite Ombrophytum subterraneum, to examine the incidence of HGT events in the mitochondrial genome (mtDNA). The mtDNA of O. subterraneum assembles into 54 circular chromosomes, only 34 of which contain the 51 full-length coding regions. Numerous foreign tracts (totaling almost 100 kb, ~ 14% of the mtDNA), including 12 intact genes, were acquired by HGT from the Asteraceae hosts. Nine chromosomes concentrate most of those regions and eight are almost entirely foreign. Native homologs of each foreign gene coexist in the mtDNA and are potentially functional. A large proportion of shorter regions were related to the Fabaceae (a total of ~ 110 kb, 15.4%), some of which were shared with L. mirabile. We also found evidence of foreign sequences donated by angiosperm lineages not reported as hosts (Apocynaceae, Euphorbiaceae, Lamiaceae, and Malvales). We propose an evolutionary hypothesis that involves ancient transfers from legume hosts in the common ancestor of Ombrophytum and Lophophytum followed by more recent transfer events in L. mirabile. Besides, the O. subterraneum mtDNA was also subjected to additional HGT events from diverse angiosperm lineages, including large and recent transfers from the Asteraceae, and also from Lamiaceae.

Mitochondrial genome of the nonphotosynthetic mycoheterotrophic plant Hypopitys monotropa, its structure, gene expression and RNA editing

Heterotrophic plants—plants that have lost the ability to photosynthesize—are characterized by a number of changes at all levels of organization. Heterotrophic plants are divided into two large categories—parasitic and mycoheterotrophic (MHT). The question of to what extent such changes are similar in these two categories is still open. The plastid genomes of nonphotosynthetic plants are well characterized, and they exhibit similar patterns of reduction in the two groups. In contrast, little is known about the mitochondrial genomes of MHT plants. We report the structure of the mitochondrial genome of Hypopitys monotropa, a MHT member of Ericaceae, and the expression of its genes. In contrast to its highly reduced plastid genome, the mitochondrial genome of H. monotropa is larger than that of its photosynthetic relative Vaccinium macrocarpon, and its complete size is ~810 Kb. We observed an unusually long repeat-rich structure of the genome that suggests the existence of linear fragments. Despite this unique feature, the gene content of the H. monotropa mitogenome is typical of flowering plants. No acceleration of substitution rates is observed in mitochondrial genes, in contrast to previous observations in parasitic non-photosynthetic plants. Transcriptome sequencing revealed the trans-splicing of several genes and RNA editing in 33 of 38 genes. Notably, we did not find any traces of horizontal gene transfer from fungi, in contrast to plant parasites, which extensively integrate genetic material from their hosts.

Why so Complex? The Intricacy of Genome Structure and Gene Expression, Associated with Angiosperm Mitochondria, May Relate to the Regulation of Embryo Quiescence or Dormancy-Intrinsic Blocks to Early Plant Life

Mitochondria play key roles in cellular-energy metabolism and are vital for plant-life, such as for successful germination and early-seedling establishment. Most mitochondria contain their own genetic system (mtDNA, mitogenome), with an intrinsic protein-synthesis machinery. Although the challenges of maintaining prokaryotic-type structures and functions are common to Eukarya, land plants possess some of the most complex organelle composition of all known organisms. Angiosperms mtDNAs are characteristically the largest and least gene-dense among the eukaryotes. They often contain highly-variable intergenic regions of endogenous or foreign origins and undergo frequent recombination events, which result in different mtDNA configurations, even between closely-related species. The expression of the mitogenome in angiosperms involves extensive mtRNA processing steps, including numerous editing and splicing events. Why do land-plant's mitochondria have to be so complex? The answer to this remains a matter of speculation. We propose that this complexity may have arisen throughout the terrestrialization of plants, as a means to control embryonic mitochondrial functions -a critical adaptive trait to optimize seed germination. The unique characteristics of plant mtDNA may play pivotal roles in the nuclear-regulation of organellar biogenesis and metabolism, possibly to control embryos quiescence or dormancy, essential determinants for the establishment of viable plantlets that can survive post-germination.

Mitochondrial plastid DNA can cause DNA barcoding paradox in plants

The transfer of ancestral plastid genomes into mitochondrial genomes to generate mitochondrial plastid DNA (MTPT) is known to occur in plants, but its impacts on mitochondrial genome complexity and the potential for causing a false-positive DNA barcoding paradox have been underestimated. Here, we assembled the organelle genomes of Cynanchum wilfordii and C. auriculatum, which are indigenous medicinal herbs in Korea and China, respectively. In both species, it is estimated that 35% of the ancestral plastid genomes were transferred to mitochondrial genomes over the past 10 million years and remain conserved in these genomes. Some plastid barcoding markers co-amplified the conserved MTPTs and caused a barcoding paradox, resulting in mis-authentication of botanical ingredients and/or taxonomic mis-positioning. We identified dynamic and lineage-specific MTPTs that have contributed to mitochondrial genome complexity and might cause a putative barcoding paradox across 81 plant species. We suggest that a DNA barcoding guidelines should be developed involving the use of multiple markers to help regulate economically motivated adulteration.

Interchangeable parts: The evolutionarily dynamic tRNA population in plant mitochondria

Transfer RNAs (tRNAs) remain one of the very few classes of genes still encoded in the mitochondrial genome. These key components of the protein translation system must interact with a large enzymatic network of nuclear-encoded gene products to maintain mitochondrial function. Plants have an evolutionarily dynamic mitochondrial tRNA population, including ongoing tRNA gene loss and replacement by both horizontal gene transfer from diverse sources and import of nuclear-expressed tRNAs from the cytosol. Thus, plant mitochondria represent an excellent model for understanding how anciently divergent genes can act as "interchangeable parts" during the evolution of complex molecular systems. In particular, understanding the integration of the mitochondrial translation system with elements of the corresponding machinery used in cytosolic protein synthesis is a key area for eukaryotic cellular evolution. Here, we review the increasingly detailed phylogenetic data about the evolutionary history of mitochondrial tRNA gene loss, transfer, and functional replacement that has created extreme variation in mitochondrial tRNA populations across plant species. We describe emerging tRNA-seq methods with promise for refining our understanding of the expression and subcellular localization of tRNAs. Finally, we summarize current evidence and identify open questions related to coevolutionary changes in nuclear-encoded enzymes that have accompanied turnover in mitochondrial tRNA populations.

Chlorella vulgaris genome assembly and annotation reveals the molecular basis for metabolic acclimation to high light conditions

Chlorella vulgaris is a fast-growing fresh-water microalga cultivated on the industrial scale for applications ranging from food to biofuel production. To advance our understanding of its biology and to establish genetics tools for biotechnological manipulation, we sequenced the nuclear and organelle genomes of Chlorella vulgaris 211/11P by combining next generation sequencing and optical mapping of isolated DNA molecules. This hybrid approach allowed us to assemble the nuclear genome in 14 pseudo-molecules with an N50 of 2.8 Mb and 98.9% of scaffolded genome. The integration of RNA-seq data obtained at two different irradiances of growth (high light, HL versus low light, LL) enabled us to identify 10 724 nuclear genes, coding for 11 082 transcripts. Moreover, 121 and 48 genes, respectively, were found in the chloroplast and mitochondrial genome. Functional annotation and expression analysis of nuclear, chloroplast and mitochondrial genome sequences revealed particular features of Chlorella vulgaris. Evidence of horizontal gene transfers from chloroplast to mitochondrial genome was observed. Furthermore, comparative transcriptomic analyses of LL versus HL provided insights into the molecular basis for metabolic rearrangement under HL versus LL conditions leading to enhanced de novo fatty acid biosynthesis and triacylglycerol accumulation. The occurrence of a cytosolic fatty acid biosynthetic pathway could be predicted and its upregulation upon HL exposure was observed, consistent with the increased lipid amount under HL conditions. These data provide a rich genetic resource for future genome editing studies, and potential targets for biotechnological manipulation of Chlorella vulgaris or other microalgae species to improve biomass and lipid productivity.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

Fluctuations in Fabaceae mitochondrial genome size and content are both ancient and recent

BACKGROUND

Organelle genome studies of Fabaceae, an economically and ecologically important plant family, have been biased towards the plastid genome (plastome). Thus far, less than 15 mitochondrial genome (mitogenome) sequences of Fabaceae have been published, all but four of which belong to the subfamily Papilionoideae, limiting the understanding of size variation and content across the family. To address this, four mitogenomes were sequenced and assembled from three different subfamilies (Cercidoideae, Detarioideae and Caesalpinioideae).

RESULTS

Phylogenetic analysis based on shared mitochondrial protein coding regions produced a fully resolved and well-supported phylogeny that was completely congruent with the plastome tree. Comparative analyses suggest that two kinds of mitogenome expansions have occurred in Fabaceae. Size expansion of four genera (Tamarindus, Libidibia, Haematoxylum, and Leucaena) in two subfamilies (Detarioideae and Caesalpinioideae) occurred in relatively deep nodes, and was mainly caused by intercellular gene transfer and/or interspecific horizontal gene transfer (HGT). The second, more recent expansion occurred in the Papilionoideae as a result of duplication of native mitochondrial sequences. Family-wide gene content analysis revealed 11 gene losses, four (rps2, 7, 11 and 13) of which occurred in the ancestor of Fabaceae. Losses of the remaining seven genes (cox2, rpl2, rpl10, rps1, rps19, sdh3, sdh4) were restricted to specific lineages or occurred independently in different clades. Introns of three genes (cox2, ccmFc and rps10) showed extensive lineage-specific length variation due to large sequence insertions and deletions. Shared DNA analysis among Fabaceae mitogenomes demonstrated a substantial decay of intergenic spacers and provided further insight into HGT between the mimosoid clade of Caesalpinioideae and the holoparasitic Lophophytum (Balanophoraceae).

CONCLUSION

This study represents the most exhaustive analysis of Fabaceae mitogenomes so far, and extends the understanding the dynamic variation in size and gene/intron content. The four newly sequenced mitogenomes reported here expands the phylogenetic coverage to four subfamilies. The family has experienced multiple mitogenome size fluctuations in both ancient and recent times. The causes of these size variations are distinct in different lineages. Fabaceae mitogenomes experienced extensive size fluctuation by recruitment of exogenous DNA and duplication of native mitochondrial DNA.

The plastomes of Astrocaryum aculeatum G. Mey. and A. murumuru Mart. show a flip-flop recombination between two short inverted repeats

The plastomes of Astrocaryum murumuru and A. aculeatum revealed a lineage-specific structural feature originated by flip-flop recombination, non-synonymous substitutions in conserved genes and several molecular markers. Astrocaryum murumuru Mart. and A. aculeatum G.Mey. are two palm species of Amazon forest that are economically important as source of food, oil and raw material for several applications. Genetic studies aiming to establish strategies for conservation and domestication of both species are still in the beginning given that the exploitation is mostly by extractive activity. The identification and characterization of molecular markers are essential to assess the genetic diversity of natural populations of both species. Therefore, we sequenced and characterized in detail the plastome of both species. We compared both species and identified 32 polymorphic SSR loci, 150 SNPs, 46 indels and eight hotspots of nucleotide diversity. Additionally, we reported a specific RNA editing site found in the ccsA gene, which is exclusive to A. murumuru. Moreover, the structural analysis in the plastomes of both species revealed a 4.6-kb inversion encompassing a set of genes involved in chlororespiration and plastid translation. This 4.6-kb inversion is a lineage-specific structural feature of the genus Astrocaryum originated by flip-flop recombination between two short inverted repeats. Furthermore, our phylogenetic analysis using whole plastomes of 39 Arecaceae species placed the Astrocaryum species sister to Acrocomia within the tribe Cocoseae. Finally, our data indicated substantial changes in the plastome structure and sequence of both species of the genus Astrocaryum, bringing new molecular markers, several structural and evolving features, which can be applied in several areas such as genetic, evolution, breeding, phylogeny and conservation strategies for both species.

The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes

Short repeats (SR) play an important role in shaping seed plant mitochondrial genomes (mtDNAs). However, their origin, distribution, and relationships across the different plant lineages remain unresolved. We focus on the angiosperm family Solanaceae that shows great variation in repeat content and extend the study to a wide diversity of seed plants. We determined the complete nucleotide sequences of the organellar genomes of the medicinal plant Physochlaina orientalis (Solanaceae), member of the tribe Hyoscyameae. To understand the evolution of the P. orientalis mtDNA we made comparisons with those of five other Solanaceae. P. orientalis mtDNA presents the largest mitogenome (∼685 kb in size) among the Solanaceae and has an unprecedented 8-copy repeat family of ∼8.2 kb in length and a great number of SR arranged in tandem-like structures. We found that the SR in the Solanaceae share a common origin, but these only expanded in members of the tribe Hyoscyameae. We discuss a mechanism that could explain SR formation and expansion in P. orientalis and Hyoscyamus niger. Finally, the great increase in plant mitochondrial data allowed us to systematically extend our repeat analysis to a total of 136 seed plants to characterize and analyze for the first time families of SR among seed plant mtDNAs.

Comparative plastid genomics of Synurophyceae: inverted repeat dynamics and gene content variation

BACKGROUND

The Synurophyceae is one of most important photosynthetic stramenopile algal lineages in freshwater ecosystems. They are characterized by siliceous scales covering the cell or colony surface and possess plastids of red-algal secondary or tertiary endosymbiotic origin. Despite their ecological and evolutionary significance, the relationships amongst extant Synurophyceae are unclear, as is their relationship to most other stramenopiles.

RESULTS

Here we report a comparative analysis of plastid genomes sequenced from five representative synurophycean algae. Most of these plastid genomes are highly conserved with respect to genome structure and coding capacity, with the exception of gene re-arrangements and partial duplications at the boundary of the inverted repeat and single-copy regions. Several lineage-specific gene loss/gain events and intron insertions were detected (e.g., cemA, dnaB, syfB, and trnL).

CONCLUSIONS

Unexpectedly, the cemA gene of Synurophyceae shows a strong relationship with sequences from members of the green-algal lineage, suggesting the occurrence of a lateral gene transfer event. Using a molecular clock approach based on silica fossil record data, we infer the timing of genome re-arrangement and gene gain/loss events in the plastid genomes of Synurophyceae.

Evolution of six novel ORFs in the plastome of Mankyua chejuense and phylogeny of eusporangiate ferns

In this paper, three plastomes of Mankyua chejuense, Helminthostachys zeylanica, and Botrychium ternatum in Ophioglossaceae were completely sequenced in order to investigate the plastome evolution and phylogeny of eusporangiate ferns. They were similar to each other in terms of length and the gene orders; however, six unknown open reading frames (ORFs) were found between rps4 and trnL-UAA genes in M. chejuense. Similar sequence regions of six ORFs of M. chejuense were found at the plastomes of Ophioglossum californicum and H. zeylanica, as well as the mitochondrial genome (mitogenome) of H. zeylanica, but not in B. ternatum. Interestingly, the translated amino acid sequences of three ORFs were more similar to the proteins of distantly related taxa such as algae and bacteria than they were to proteins in land plants. It is likely that the six ORFs region arose from endosymbiotic gene transfer (EGT) or horizontal gene transfer (HGT), but further study is needed to verify this. Phylogenetic analyses suggested that Mankyua was resolved as the earliest diverging lineage and that Ophioglossum was subsequently diverged in Ophioglossaceae. This result supports why the plastome of M. chejuense have contained the most ancestral six ORFs in the family.