The massive mitochondrial genome of the angiosperm Silene noctiflora is evolving by gain or loss of entire chromosomes

Assemble and comparative analysis of the mitochondrial genome of Rhododendron delavayi: Insights into phylogenetic relationships and genomic variations

Rhododendron delavayi, a notable ornamental plant primarily found in regions of China like Yunnan and Guizhou provinces, holds substantial horticultural value. To elucidate the systematic phylogenetic relationships and organelle genomic differences within R. delavayi and related Rhododendron species, we conducted sequencing and assembly of the complete mitochondrial genome of R. delavayi. The full-length mitochondrial genome of it was a singular circular molecule spanning 1,009,263 bp, comprising 53 protein-coding genes, including 18 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 32 protein-coding genes. A total of 1,182 simple sequence repeats (SSRs) loci were identified in the R. delavayi mitochondrial genome, primarily consisting of single nucleotide, dinucleotide, and trinucleotide repeats. Nucleotide diversity analysis highlighted five genes (atp6, atp9, cox2, nad1, and rpl10) with the highest diversity within the mitochondrial genomes of Rhododendron genus. Comparative analysis of the mitochondrial genome of R. delavayi with those of four other Rhododendron species indicated complex rearrangements in 21 genes, including rps4, nad6, rps3, atp6, cob, atp9, nad7, among others. The mitochondrial phylogenetic tree revealed a close relationship between R. delavayi and R. decorum, forming a sister clade to R. × pulchrum and R. simsii. Furthermore, 126 plastid-to-mitochondrial gene transfers in R. delavayi were identified, ranging from 30 bp to 19,385 bp. These fragments collectively constituted 47.54 % and 9.52 % of the chloroplast and mitochondrial genomes (202,169 bp), respectively. Complex mitochondrial-to-mitochondrial transfers were also observed, with 843 identified fragments totaling 312,036 bp (30.92 % of the mitochondrial genome). Segments exceeding 10 kb may mediate homologous recombination within the mitochondrial molecules. Remarkably, our study underscores that the mitochondrial genome of R. delavayi was the largest reported within the Rhododendron genus to date. The intricate rearrangements observed in the mitochondrial genomes of Rhododendron species, alone with the identification of five potential molecular marker sites, provided valuable insights for species classification and parentage identification within the Rhododendron genus.

Horizontally transferred mitochondrial DNA tracts become circular by microhomology-mediated repair pathways

The holoparasitic plant Lophophytum mirabile exhibits remarkable levels of mitochondrial horizontal gene transfer (HGT). Gathering comparative data from other individuals and host plants can provide insights into the HGT process. We sequenced the mitochondrial genome (mtDNA) from individuals of two species of Lophophytum and from mimosoid hosts. We applied a stringent phylogenomic approach to elucidate the origin of the whole mtDNAs, estimate the timing of the transfers, and understand the molecular mechanisms involved. Ancestral and recent HGT events replaced and enlarged the multichromosomal mtDNA of Lophophytum spp., with the foreign DNA ascending to 74%. A total of 14 foreign mitochondrial chromosomes originated from continuous regions in the host mtDNA flanked by short direct repeats. These foreign tracts are circularized by microhomology-mediated repair pathways and replicate independently until they are lost or they eventually recombine with other chromosomes. The foreign noncoding chromosomes are variably present in the population and likely evolve by genetic drift. We present the 'circle-mediated HGT' model in which foreign mitochondrial DNA tracts become circular and are maintained as plasmid-like molecules. This model challenges the conventional belief that foreign DNA must be integrated into the recipient genome for successful HGT.

Assembly and comparative analysis of the complete mitochondrial genome of Fritillaria ussuriensis Maxim. (Liliales: Liliaceae), an endangered medicinal plant

BACKGROUND

Fritillaria ussuriensis is an endangered medicinal plant known for its notable therapeutic properties. Unfortunately, its population has drastically declined due to the destruction of forest habitats. Thus, effectively protecting F. ussuriensis from extinction poses a significant challenge. A profound understanding of its genetic foundation is crucial. To date, research on the complete mitochondrial genome of F. ussuriensis has not yet been reported.

RESULTS

The complete mitochondrial genome of F. ussuriensis was sequenced and assembled by integrating PacBio and Illumina sequencing technologies, revealing 13 circular chromosomes totaling 737,569 bp with an average GC content of 45.41%. A total of 55 genes were annotated in this mitogenome, including 2 rRNA genes, 12 tRNA genes, and 41 PCGs. The mitochondrial genome of F. ussuriensis contained 192 SSRs and 4,027 dispersed repeats. In the PCGs of F. ussuriensis mitogenome, 90.00% of the RSCU values exceeding 1 exhibited a preference for A-ended or U-ended codons. In addition, 505 RNA editing sites were predicted across these PCGs. Selective pressure analysis suggested negative selection on most PCGs to preserve mitochondrial functionality, as the notable exception of the gene nad3 showed positive selection. Comparison between the mitochondrial and chloroplast genomes of F. ussuriensis revealed 20 homologous fragments totaling 8,954 bp. Nucleotide diversity analysis revealed the variation among genes, and gene atp9 was the most notable. Despite the conservation of GC content, mitogenome sizes varied significantly among six closely related species, and colinear analysis confirmed the lack of conservation in their genomic structures. Phylogenetic analysis indicated a close relationship between F. ussuriensis and Lilium tsingtauense.

CONCLUSIONS

In this study, we sequenced and annotated the mitogenome of F. ussuriensis and compared it with the mitogenomes of other closely related species. In addition to genomic features and evolutionary position, this study also provides valuable genomic resources to further understand and utilize this medicinal plant.

A Systematic Review and Developmental Perspective on Origin of CMS Genes in Crops

Cytoplasmic male sterility (CMS) arises from the incompatibility between the nucleus and cytoplasm as typical representatives of the chimeric structures in the mitochondrial genome (mitogenome), which has been extensively applied for hybrid seed production in various crops. The frequent occurrence of chimeric mitochondrial genes leading to CMS is consistent with the mitochondrial DNA (mtDNA) evolution. The sequence conservation resulting from faithfully maternal inheritance and the chimeric structure caused by frequent sequence recombination have been defined as two major features of the mitogenome. However, when and how these chimeric mitochondrial genes appear in the context of the highly conserved reproduction of mitochondria is an enigma. This review, therefore, presents the critical view of the research on CMS in plants to elucidate the mechanisms of this phenomenon. Generally, distant hybridization is the main mechanism to generate an original CMS source in natural populations and in breeding. Mitochondria and mitogenomes show pleomorphic and dynamic changes at key stages of the life cycle. The promitochondria in dry seeds develop into fully functioning mitochondria during seed imbibition, followed by massive mitochondria or mitogenome fusion and fission in the germination stage along with changes in the mtDNA structure and quantity. The mitogenome stability is controlled by nuclear loci, such as the nuclear gene Msh1. Its suppression leads to the rearrangement of mtDNA and the production of heritable CMS genes. An abundant recombination of mtDNA is also often found in distant hybrids and somatic/cybrid hybrids. Since mtDNA recombination is ubiquitous in distant hybridization, we put forward a hypothesis that the original CMS genes originated from mtDNA recombination during the germination of the hybrid seeds produced from distant hybridizations to solve the nucleo-cytoplasmic incompatibility resulting from the allogenic nuclear genome during seed germination.

The Complete Mitogenome of Apostasia fujianica Y.Li & S.Lan and Comparative Analysis of Mitogenomes across Orchidaceae

Apostasia fujianica belongs to the genus Apostasia and is part of the basal lineage in the phylogenetic tree of the Orchidaceae. Currently, there are only ten reported complete mitochondrial genomes in orchids, which greatly hinders the understanding of mitochondrial evolution in Orchidaceae. Therefore, we assembled and annotated the mitochondrial genome of A. fujianica, which has a length of 573,612 bp and a GC content of 44.5%. We annotated a total of 44 genes, including 30 protein-coding genes, 12 tRNA genes, and two rRNA genes. We also performed relative synonymous codon usage (RSCU) analysis, repeat sequence analysis, intergenomic transfer (IGT) analysis, and Ka/Ks analysis for A. fujianica and conducted RNA editing site analysis on the mitochondrial genomes of eight orchid species. We found that most protein-coding genes are under purifying selection, but nad6 is under positive selection, with a Ka/Ks value of 1.35. During the IGT event in A. fujianica's mitogenome, the trnN-GUU, trnD-GUC, trnW-CCA, trnP-UGG, and psaJ genes were identified as having transferred from the plastid to the mitochondrion. Compared to other monocots, the family Orchidaceae appears to have lost the rpl10, rpl14, sdh3, and sdh4 genes. Additionally, to further elucidate the evolutionary relationships among monocots, we constructed a phylogenetic tree based on the complete mitogenomes of monocots. Our study results provide valuable data on the mitogenome of A. fujianica and lay the groundwork for future research on genetic variation, evolutionary relationships, and breeding of Orchidaceae.

Assembly and comparative genome analysis of four mitochondrial genomes from Saccharum complex species

Saccharum complex includes genera Saccharum, Miscanthus, Erianthus, Narenga, and Tripidium. Since the Saccharum complex/Saccharinae constitutes the gene pool used by sugarcane breeders to introduce useful traits into sugarcane, studying the genomic characterization of the Saccharum complex has become particularly important. Here, we assembled graph-based mitochondrial genomes (mitogenomes) of four Saccharinae species (T. arundinaceum, E. rockii, M. sinensis, and N. porphyrocoma) using Illumina and PacBio sequencing data. The total lengths of the mitogenomes of T. arundinaceum, M. sinensis, E. rockii and N. porphyrocoma were 549,593 bp, 514,248 bp, 481,576 bp and 513,095 bp, respectively. Then, we performed a comparative mitogenomes analysis of Saccharinae species, including characterization, organelles transfer sequence, collinear sequence, phylogenetics analysis, and gene duplicated/loss. Our results provided the mitogenomes of four species closely related to sugarcane breeding, enriching the mitochondrial genomic resources of the Saccharinae. Additionally, our study offered new insights into the evolution of mitogenomes at the family and genus levels and enhanced our understanding of organelle evolution in the highly polyploid Saccharum genus.

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata

BACKGROUND

Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species.

RESULTS

In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade.

CONCLUSIONS

We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

Graph-based mitochondrial genomes of three foundation species in the Saccharum genus

KEY MESSAGE

We reported the graph-based mitochondrial genomes of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) for the first time. The results revealed pan-structural variation and evolutionary processes in the mitochondrial genomes within Saccharum. Saccharum belongs to the Andropogoneae, and cultivars species in Saccharum contribute nearly 80% of sugar production in the world. To explore the genomic studies in Saccharum, we assembled 15 complete mitochondrial genomes (mitogenome) of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) using Illumina and Oxford Nanopore Technologies sequencing data. The mitogenomes of the three species were divided into a total of eight types based on contig numbers and linkages. All mitogenomes in the three species encoded 51 unique genes, including 32 protein-coding, 3 ribosomal RNA (rRNA) and 16 transfer RNA (tRNA) genes. The existence of long and short-repeat-mediated recombinations in the mitogenome of S. officinarum and S. robustum was revealed and confirmed through PCR validation. Furthermore, employing comparative genomics and phylogenetic analyses of the organelle genomes, we unveiled the evolutionary relationships and history of the major interspecific lineages in Saccharum genus. Phylogenetic analyses of homologous fragments between S. officinarum and S. robustum showed that S. officinarum and S. robustum are phylogenetically distinct and that they were likely parallel rather than domesticated. The variations between ancient (S. sinense and S. barberi) and modern cultivated species (S. hybrid) possibly resulted from hybridization involving different S. officinarum accessions. Lastly, this project reported the first graph-based mitogenomes of three Saccharum species, and a systematic comparison of the structural organization, evolutionary processes, and pan-structural variation of the Saccharum mitogenomes revealed the differential features of the Saccharum mitogenomes.

De Novo Hybrid Assembly Unveils Multi-Chromosomal Mitochondrial Genomes in Ludwigia Species, Highlighting Genomic Recombination, Gene Transfer, and RNA Editing Events

Biological invasions have been identified as the fifth cause of biodiversity loss, and their subsequent dispersal represents a major ecological challenge. The aquatic invasive species Ludwigia grandiflora subsp. hexapetala (Lgh) and Ludwigia peploides subsp. montevidensis (Lpm) are largely distributed in aquatic environments in North America and in Europe. However, they also present worrying terrestrial forms that are able to colonize wet meadows. To comprehend the mechanisms of the terrestrial adaptation of Lgh and Lpm, it is necessary to develop their genomic resources, which are currently poorly documented. We performed de novo assembly of the mitogenomes of Lgh and Lpm through hybrid assemblies, combining short reads (SR) and/or long reads (LR) before annotating both mitogenomes. We successfully assembled the mitogenomes of Lgh and Lpm into two circular molecules each, resulting in a combined total length of 711,578 bp and 722,518 bp, respectively. Notably, both the Lgh and Lpm molecules contained plastome-origin sequences, comprising 7.8% of the mitochondrial genome length. Additionally, we identified recombinations that were mediated by large repeats, suggesting the presence of multiple alternative conformations. In conclusion, our study presents the first high-quality mitogenomes of Lpm and Lgh, which are the only ones in the Myrtales order found as two circular molecules.

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BACKGROUND

Orchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.

RESULTS

Here, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.

CONCLUSION

In this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

Retention of an Endosymbiont for the Production of a Single Molecule

Sap-feeding insects often maintain two or more nutritional endosymbionts that act in concert to produce compounds essential for insect survival. Many mealybugs have endosymbionts in a nested configuration: one or two bacterial species reside within the cytoplasm of another bacterium, and together, these bacteria have genomes that encode interdependent sets of genes needed to produce key nutritional molecules. Here, we show that the mealybug Pseudococcus viburni has three endosymbionts, one of which contributes only two unique genes that produce the host nutrition-related molecule chorismate. All three bacterial endosymbionts have tiny genomes, suggesting that they have been coevolving inside their insect host for millions of years.

Assembly and comparative analysis of the complete mitochondrial genome of Pinellia ternata

Pinellia ternata is an important natural medicinal herb in China. However, it is susceptible to withering when exposed to high temperatures during growth, which limits its tuber production. Mitochondria usually function in stress response. The P . ternata mitochondrial (mt) genome has yet to be explored. Therefore, we integrated PacBio and Illumina sequencing reads to assemble and annotate the mt genome of P . ternata . The circular mt genome of P . ternata is 876 608bp in length and contains 38 protein-coding genes (PCGs), 20 tRNA genes and three rRNA genes. Codon usage, sequence repeats, RNA editing and gene migration from chloroplast (cp) to mt were also examined. Phylogenetic analysis based on the mt genomes of P . ternata and 36 other taxa revealed the taxonomic and evolutionary status of P . ternata . Furthermore, we investigated the mt genome size and GC content by comparing P . ternata with the other 35 species. An evaluation of non-synonymous substitutions and synonymous substitutions indicated that most PCGs in the mt genome underwent negative selection. Our results provide comprehensive information on the P . ternata mt genome, which may facilitate future research on the high-temperature response of P . ternata and provide new molecular insights on the Araceae family.

Insights into the multi-chromosomal mitochondrial genome structure of the xero-halophytic plant Haloxylon Ammodendron (C.A.Mey.) Bunge ex Fenzl

BACKGROUND

Haloxylon ammodendron holds significance as an ecological plant, showcasing remarkable adaptability to desert conditions, halophytic environments, and sand fixation. With its potential for carbon sequestration, it emerges as a promising candidate for environmental sustainability. Furthermore, it serves as a valuable C4 plant model, offering insights into the genetic foundations of extreme drought tolerance. Despite the availability of plastid and nuclear genomes, the absence of a mitochondrial genome (mitogenome or mtDNA) hinders a comprehensive understanding of its its mtDNA structure, organization, and phylogenetic implications.

RESULTS

In the present study, the mitochondrial genome of H. ammodendron was assembled and annotated, resulting in a multi-chromosomal configuration with two circular chromosomes. The mtDNA measured 210,149 bp in length and contained 31 protein-coding genes, 18 tRNA and three rRNA. Our analysis identified a total of 66 simple sequence repeats along with 27 tandem repeats, 312 forward repeats, and 303 palindromic repeats were found. Notably, 17 sequence fragments displayed homology between the mtDNA and chloroplast genome (cpDNA), spanning 5233 bp, accounting for 2.49% of the total mitogenome size. Additionally, we predicted 337 RNA editing sites, all of the C-to-U conversion type. Phylogenetic inference confidently placed H. ammodendron in the Amaranthacea family and its close relative, Suaeda glacum.

CONCLUSIONS

H. ammodendron mtDNA showed a multi-chromosomal structure with two fully circularized molecules. This newly characterized mtDNA represents a valuable resource for gaining insights into the basis of mtDNA structure variation within Caryophyllales and the evolution of land plants, contributing to their identification, and classification.

microRNAs in Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma) possibly involved in the mitochondrial function

Mitochondria are essential organelles for maintaining vital cellular functions, and microRNAs (miRNAs) regulate gene expression posttranscriptionally. miRNAs exhibit tissue and time-specific patterns in mitochondria and specifically mitochondrial miRNAs (mitomiRs) can regulate the mRNA expression both originating from mitochondrial and nuclear transcription which affect mitochondrial metabolic activity and cell homeostasis. In this study, miRNAs of two insect species, Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma), were investigated for the first time. The known and possible novel miRNAs were predicted and characterized and their potential effects on mitochondrial transcription were investigated in these insect species using deep sequencing. The previously reported mitomiRs were also investigated and housekeeping miRNAs were characterized. miRNAs that are involved in mitochondrial processes such as apoptosis and signaling and that affect genes encoding the subunits of OXPHOS complexes have been identified in each species. Here, 81 and 161 novel mature miRNA candidates were bioinformatically predicted and 9 and 24 of those were aligned with reference mitogenomes of S. parreyssi and L. saccharina, respectively. As a result of RNAHybrid analysis, 51 and 69 potential targets of miRNAs were found in the mitogenome of S. parreyssi and L. saccharina, respectively. cox1 gene was the most targeted gene and cytB, rrnS, and rrnL genes were highly targeted in both of the species by novel miRNAs, hypothetically. We speculate that these novel miRNAs, originating from or targeting mitochondria, influence on rRNA genes or positively selected mitochondrial protein-coding genes. These findings may provide a new perspective in evaluating miRNAs for maintaining mitochondrial function and transcription.

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.

Complete mitochondrial genome of Mentha spicata L. reveals multiple chromosomal configurations and RNA editing events

Mentha spicata L. is a valuable plant that yields spearmint oil, widely utilized in the pharmaceutical, chemical, and cosmetic industries. The mitochondrial genome (mitogenome) is an essential material for molecular breeding and evolution studies. Here, the mitogenome of M. spicata was assembled by combining Nanopore and Illumina reads. It consisted of a linear chromosome (Ch1) and two circular chromosomes (Ch2 and Ch3). Furthermore, we showed two pairs of repeats (R1 and R2) mediated recombinations resulting in multiple chromosomal configurations. The R1-mediated-recombination generated a large molecule formed by joining Ch2 and Ch1. Similarly, the R2-mediated-recombination generated a large molecule formed by joining Ch3 and Ch1. Then, we identified 17 mitochondrial plastid DNAs (MTPTs) by comparing the mitogenome and cpgenome. The MTPT14 was conserved in multiple species, which has undergone the same evolutionary process as the two organellar genomes among M. spicata, Hesperelaea palmeri and Castilleja paramensis. Based on the RNA-seq reads, 246 RNA editing sites were predicted, resulting in the conversion of cytosine to uracil bases. Furthermore, we successfully validated 40 out of 43 predicted sites. This project reported a complex structure of the M. spicata mitogenome resulting from repeat-mediated recombinations, which will provide valuable information for gene function study and the breeding of different varieties.

Assembly and phylogenetic analysis of the mitochondrial genome of endangered medicinal plant Huperzia crispata

Huperzia crispata is a traditional Chinese herb plant and has attracted special attention in recent years for its products Hup A can serve as an acetylcholinesterase inhibitor (AChEI). Although the chloroplast (cp) genome of H. crispata has been studied, there are no reports regarding the Huperzia mitochondrial (mt) genome since the previously reported H. squarrosa has been revised as Phlegmariurus squarrosus. The mt genome of H. crispata was sequenced using a combination of long-read nanopore and Illumina sequencing platforms. The entire H. crispata mt genome was assembled in a circular with a length of 412,594 bp and a total of 91 genes, including 45 tRNAs, 6 rRNAs, 37 protein-coding genes (PCGs), and 3 pseudogenes. Notably, the rps8 gene was present in P. squarrosus and a pseudogene rps8 was presented in H. crispata, which was lacking in most of Pteridophyta and Gymnospermae. Intron-encoded maturase (mat-atp9i85 and mat-cobi787) genes were present in H. crispata and P. squarrosus, but lost in other examined lycophytes, ferns, and Gymnospermae plants. Collinearity analysis showed that the mt genome of H. crispata and P. squarrossus is highly conservative compared to other ferns. Relative synonymous codon usage (RSCU) analysis showed that the amino acids most frequently found were phenylalanine (Phe) (4.77%), isoleucine (Ile) (4.71%), lysine (Lys) (4.26%), while arginine (Arg) (0.32%), and histidine (His) (0.42%) were rarely found. Simple sequence repeats (SSR) analysis revealed that a total of 114 SSRs were identified in the mt genome of H. crispata and account for 0.35% of the whole mt genome. Monomer repeats were the majority types of SSRs and represent 91.89% of the total SSRs. In addition, a total of 1948 interspersed repeats (158 forward, 147 palindromic, and 5 reverse repeats) with a length ranging from 30 bp to 14,945 bp were identified in the H. crispata mt genome and the 30-39-bp repeats were the most abundant type. Gene transfer analysis indicated that a total of 12 homologous fragments were discovered between the cp and mt genomes of H. crispata, accounting for 0.93% and 2.48% of the total cp and mt genomes, respectively. The phylogenetic trees revealed that H. crispata was the sister of P. squarrosus. The Ka/Ks analysis results suggested that most PCGs, except atp6 gene, were subject to purification selection during evolution. Our study provides extensive information on the features of the H. crispata mt genome and will help unravel evolutionary relationships, and molecular identification within lycophytes.

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

The Mitochondrial Genome of the Holoparasitic Plant Thonningia sanguinea Provides Insights into the Evolution of the Multichromosomal Structure

Multichromosomal mitochondrial genome (mitogenome) structures have repeatedly evolved in many lineages of angiosperms. However, the underlying mechanism remains unclear. The mitogenomes of three genera of Balanophoraceae, namely Lophophytum, Ombrophytum, and Rhopalocnemis, have already been sequenced and assembled, all showing a highly multichromosomal structure, albeit with different genome and chromosome sizes. It is expected that characterization of additional lineages of this family may expand the knowledge of mitogenome diversity and provide insights into the evolution of the plant mitogenome structure and size. Here, we assembled and characterized the mitogenome of Thonningia sanguinea, which, together with Balanophora, forms a clade sister to the clade comprising Lophophytum, Ombrophytum, and Rhopalocnemis. The mitogenome of T. sanguinea possesses a multichromosomal structure of 18 circular chromosomes of 8.7-19.2 kb, with a total size of 246,247 bp. There are very limited shared regions and poor chromosomal correspondence between T. sanguinea and other Balanophoraceae species, suggesting frequent rearrangements and rapid sequence turnover. Numerous medium- and small-sized repeats were identified in the T. sanguinea mitogenome; however, no repeat-mediated recombination was detected, which was verified by Illumina reads mapping and PCR experiments. Intraspecific mitogenome variations in T. sanguinea are mostly insertions and deletions, some of which can lead to degradation of perfect repeats in one or two accessions. Based on the mitogenome features of T. sanguinea, we propose a mechanism to explain the evolution of a multichromosomal mitogenome from a master circle, which involves mutation in organellar DNA replication, recombination and repair genes, decrease of recombination, and repeat degradation.

Structural and gene composition variation of the complete mitochondrial genome of Mammillaria huitzilopochtli (Cactaceae, Caryophyllales), revealed by de novo assembly

BACKGROUND

Structural descriptions of complete genomes have elucidated evolutionary processes in angiosperms. In Cactaceae (Caryophyllales), a high structural diversity of the chloroplast genome has been identified within and among genera. In this study, we assembled the first mitochondrial genome (mtDNA) for the short-globose cactus Mammillaria huitzilopochtli. For comparative purposes, we used the published genomes of 19 different angiosperms and the gymnosperm Cycas taitungensis as an external group for phylogenetic issues.

RESULTS

The mtDNA of M. huitzilopochtli was assembled into one linear chromosome of 2,052,004 bp, in which 65 genes were annotated. These genes account for 57,606 bp including 34 protein-coding genes (PCGs), 27 tRNAs, and three rRNAs. In the non-coding sequences, repeats were abundant, with a total of 4,550 (179,215 bp). In addition, five complete genes (psaC and four tRNAs) of chloroplast origin were documented. Negative selection was estimated for most (23) of the PCGs. The phylogenetic tree showed a topology consistent with previous analyses based on the chloroplast genome.

CONCLUSIONS

The number and type of genes contained in the mtDNA of M. huitzilopochtli were similar to those reported in 19 other angiosperm species, regardless of their phylogenetic relationships. Although other Caryophyllids exhibit strong differences in structural arrangement and total size of mtDNA, these differences do not result in an increase in the typical number and types of genes found in M. huitzilopochtli. We concluded that the total size of mtDNA in angiosperms increases by the lengthening of the non-coding sequences rather than a significant gain of coding genes.

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.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Two Species of Calla Lilies (Zantedeschia, Araceae)

In this study, the mitochondrial genomes of two calla species, Zantedeschia aethiopica Spreng. and Zantedeschia odorata Perry., were assembled and compared for the first time. The Z. aethiopica mt genome was assembled into a single circular chromosome, measuring 675,575 bp in length with a 45.85% GC content. In contrast, the Z. odorata mt genome consisted of bicyclic chromosomes (chromosomes 1 and 2), measuring 719,764 bp and exhibiting a 45.79% GC content. Both mitogenomes harbored similar gene compositions, with 56 and 58 genes identified in Z. aethiopica and Z. odorata, respectively. Analyses of codon usage, sequence repeats, gene migration from chloroplast to mitochondrial, and RNA editing were conducted for both Z. aethiopica and Z. odorata mt genomes. Phylogenetic examination based on the mt genomes of these two species and 30 other taxa provided insights into their evolutionary relationships. Additionally, the core genes in the gynoecium, stamens, and mature pollen grains of the Z. aethiopica mt genome were investigated, which revealed maternal mitochondrial inheritance in this species. In summary, this study offers valuable genomic resources for future research on mitogenome evolution and the molecular breeding of calla lily.

When and why are mitochondria paternally inherited?

In contrast with nuclear genes that are passed on through both parents, mitochondrial genes are maternally inherited in most species, most of the time. The genetic conflict stemming from this transmission asymmetry is well-documented, and there is an abundance of population-genetic theory associated with it. While occasional or aberrant paternal inheritance occurs, there are only a few cases where exclusive paternal inheritance of mitochondrial genomes is the evolved state. Why this is remains poorly understood. By examining commonalities between species with exclusive paternal inheritance, we discuss what they may tell us about the evolutionary forces influencing mitochondrial inheritance patterns. We end by discussing recent technological advances that make exploring the causes and consequences of paternal inheritance feasible.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

Mitochondrial Genome Sequence of Salvia officinalis (Lamiales: Lamiaceae) Suggests Diverse Genome Structures in Cogeneric Species and Finds the Stop Gain of Genes through RNA Editing Events

Our previous study was the first to confirm that the predominant conformation of mitochondrial genome (mitogenome) sequence of Salvia species contains two circular chromosomes. To further understand the organization, variation, and evolution of Salvia mitogenomes, we characterized the mitogenome of Salvia officinalis. The mitogenome of S. officinalis was sequenced using Illumina short reads and Nanopore long reads and assembled using a hybrid assembly strategy. We found that the predominant conformation of the S. officinalis mitogenome also had two circular chromosomes that were 268,341 bp (MC1) and 39,827 bp (MC2) in length. The S. officinalis mitogenome encoded an angiosperm-typical set of 24 core genes, 9 variable genes, 3 rRNA genes, and 16 tRNA genes. We found many rearrangements of the Salvia mitogenome through inter- and intra-specific comparisons. A phylogenetic analysis of the coding sequences (CDs) of 26 common protein-coding genes (PCGs) of 11 Lamiales species and 2 outgroup taxa strongly indicated that the S. officinalis was a sister taxon to S. miltiorrhiza, consistent with the results obtained using concatenated CDs of common plastid genes. The mapping of RNA-seq data to the CDs of PCGs led to the identification of 451 C-to-U RNA editing sites from 31 PCGs of the S. officinalis mitogenome. Using PCR amplification and Sanger sequencing methods, we successfully validated 113 of the 126 RNA editing sites from 11 PCGs. The results of this study suggest that the predominant conformation of the S. officinalis mitogenome are two circular chromosomes, and the stop gain of rpl5 was found through RNA editing events of the Salvia mitogenome.

De Novo Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Chaenomeles speciosa (Sweet) Nakai Revealed the Existence of Two Structural Isomers

As a valuable Chinese traditional medicinal species, Chaenomeles speciosa (Sweet) Nakai (C. speciosa) is a natural resource with significant economic and ornamental value. However, its genetic information is not well understood. In this study, the complete mitochondrial genome of C. speciosa was assembled and characterized to explore the repeat sequences, recombination events, rearrangements, and IGT, to predict RNA editing sites, and to clarify the phylogenetic and evolutionary relationship. The C. speciosa mitochondrial genome was found to have two circular chromosomes as its major conformation, with a total length of 436,464 bp and 45.2% GC content. The mitochondrial genome contained 54 genes, including 33 unique protein-coding genes, 18 tRNAs, and 3 rRNA genes. Seven pairs of repeat sequences involving recombination events were analyzed. Both the repeat pairs, R1 and R2, played significant roles in mediating the major and minor conformations. In total, 18 MTPTs were identified, 6 of which were complete tRNA genes. There were 454 RNA editing sites in the 33 protein-coding sequences predicted by the PREPACT3 program. A phylogenetic analysis based on 22 species of mitochondrial genomes was constructed and indicated highly conserved PCG sequences. Synteny analyses showed extensive genomic rearrangements in the mitochondrial genome of C. speciosa and closely related species. This work is the first to report the C. speciosa mitochondrial genome, which is of great significance for conducting additional genetic studies on this organism.

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.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution

BACKGROUND

Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism.

RESULTS

Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa.

CONCLUSIONS

We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.

Mitochondrial genome sequencing and analysis of the invasive Microstegium vimineum: a resource for systematics, invasion history, and management

Premise of the Research

Plants remain underrepresented among species with sequenced mitochondrial genomes (mitogenomes), due to the difficulty in assembly with short-read technology. Invasive species lag behind crops and other economically important species in this respect, representing a lack of tools for management and land conservation efforts.

Methodology

The mitogenome of Microstegium vimineum, one of the most damaging invasive plant species in North America, was sequenced and analyzed using long-read data, providing a resource for biologists and managers. We conducted analyses of genome content, phylogenomic analyses among grasses and relatives based on mitochondrial coding regions, and an analysis of mitochondrial single nucleotide polymorphism in this invasive grass species.

Pivotal Results

The assembly is 478,010 bp in length and characterized by two large, inverted repeats, and a large, direct repeat. However, the genome could not be circularized, arguing against a "master circle" structure. Long-read assemblies with data subsets revealed several alternative genomic conformations, predominantly associated with large repeats. Plastid-like sequences comprise 2.4% of the genome, with further evidence of Class I and Class II transposable element-like sequences. Phylogenetic analysis placed M. vimineum with other Microstegium species, excluding M. nudum, but with weak support. Analysis of polymorphic sites across 112 accessions of M. vimineum from the native and invasive ranges revealed a complex invasion history.

Conclusions

We present an in-depth analysis of mitogenome structure, content, phylogenetic relationships, and range-wide genomic variation in M. vimineum's invasive US range. The mitogenome of M. vimineum is typical of other andropogonoid grasses, yet mitochondrial sequence variation across the invasive and native ranges is extensive. Our findings suggest multiple introductions to the US over the last century, with subsequent spread, secondary contact, long-distance dispersal, and possibly post-invasion selection on awn phenotypes. Efforts to produce genomic resources for invasive species, including sequenced mitochondrial genomes, will continue to provide tools for their effective management, and to help predict and prevent future invasions.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

Assembly of the Complete Mitochondrial Genome of Gelsemium elegans Revealed the Existence of Homologous Conformations Generated by a Repeat Mediated Recombination

Gelsemium elegans (G. elegans) is a Chinese medicinal plant with substantial economic and feeding values. There is a lack of detailed studies on the mitochondrial genome of G. elegans. In this study, the mitochondrial genome of G. elegans was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of G. elegans is represented by two circular chromosomes of 406,009 bp in length with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected 145 pairs of repeats and found that four pairs of repeats could mediate the homologous recombination into one major conformation and five minor conformations, and the presence of conformations was verified by PCR amplification and Sanger sequencing. A total of 124 SSRs were identified in the G. elegans mitochondrial genome. The homologous segments between the chloroplast and mitochondrial genomes accounted for 5.85% of the mitochondrial genome. We also predicted 477 RNA potential editing sites and found that the nad4 gene was edited 38 times, which was the most frequent occurrence. Taken together, the mitochondrial genome of G. elegans was assembled and annotated. We gained a more comprehensive understanding on the genome of this medicinal plant, which is vital for its effective utilization and genetic improvement, especially for cytoplasmic male sterility breeding and evolution analysis in G. elegans.

Comparative analysis of mitochondrial genomes of Broussonetia spp. (Moraceae) reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

The genus Broussonetia (Moraceae) is comprised of three non-hybrid recognized species that all produce high quality fiber essential in the development of papermaking and barkcloth-making technology. In addition, these species also have medicinal value in several countries. Despite their important economical, medicinal, and ecological values, the complete mitogenome of Broussonetia has not been reported and investigated, which would greatly facilitate molecular phylogenetics, species identification and understanding evolutionary processes. Here, we assembled the first-reported three complete Broussonetia (B. papyrifera, B. kaempferi, and B. monoica) mitochondrial genomes (mitogenome) based on a hybrid strategy using Illumina and Oxford Nanopore Technology sequencing data, and performed comprehensive comparisons in terms of their structure, gene content, synteny, intercellular gene transfer, phylogeny, and RNA editing. Our results showed their huge heterogeneities among the three species. Interestingly, the mitogenomes of B. monoica and B. papyrifera consisted of a single circular structure, whereas the B. kaempferi mitogenome was unique and consisted of a double circular structure. Gene content was consistent except for a few transfer RNA (tRNA) genes. The Broussonetia spp. mitogenomes had high sequence conservation but B. monoica with B. kaempferi contained more synteny blocks and were more related, a finding that was well-supported in organellar phylogeny. Fragments that had been transferred between mitogenomes were detected at plastome hotspots that had integrated under potential mediation of tRNA genes. In addition, RNA editing sites showed great differences in abundance, type, location and efficiency among species and tissues. The availability of these complete gap-free mitogenomes of Broussonetia spp. will provide a valuable genetic resource for evolutionary research and understanding the communications between the two organelle genomes.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Highly active repeat-mediated recombination in the mitogenome of the holoparasitic plant Aeginetia indica

Mitogenomes of most flowering plants evolve slowly in sequence, but rapidly in structure. The rearrangements in structure are mainly caused by repeat-mediated recombination. However, patterns of repeat-mediated recombination vary substantially among plants, and to provide a comprehensive picture, characterization of repeat-mediated recombination should extend to more plant species, including parasitic plants with a distinct heterotrophic lifestyle. Here we assembled the mitogenome of the holoparasitic plant Aeginetia indica (Orobanchaceae) using Illumina sequencing reads. The mitogenome was assembled into a circular chromosome of 420,362 bp, 18,734 bp longer than that of another individual of A. indica which was assembled before as a linear molecule. Synteny analysis between the two mitogenomes revealed numerous rearrangements, unique regions of each individual and 0.2% sequence divergence in their syntenic regions. The A. indica mitogenome contains a gene content typical of flowering plants (33 protein-coding, 3 rRNA, and 17 tRNA genes). Repetitive sequences >30 bp in size totals 57,060 bp, representing 13.6% of the mitogenome. We examined recombination mediated by repeats >100 bp in size and found highly active recombination for all the repeats, including a very large repeat of ~16 kb. Recombination between these repeats can form much smaller subgenomic circular chromosomes, which may lead to rapid replication of mitochondrial DNA and thus be advantageous for A. indica with a parasitic lifestyle. In addition, unlike some other parasitic plants, A. indica shows no evidence for horizontal gene transfer of protein-coding genes in its mitogenome.

Characterization, comparative phylogenetic, and gene transfer analyses of organelle genomes of Rhododendron × pulchrum

Rhododendron × pulchrum, an important horticultural species, is widely distributed in Europe, Asia, and North America. To analyze the phylogenetic and organelle genome information of R. × pulchrum and its related species, the organelle genome of R. × pulchrum was sequenced and assembled. The complete mitochondrial genome showed lineage DNA molecules, which were 816,410 bp long and contained 64 genes, namely 24 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 37 protein-coding genes. The chloroplast genome of R. × pulchrum was reassembled and re-annotated; the results were different from those of previous studies. There were 42 and 46 simple sequence repeats (SSR) identified from the mitochondrial and chloroplast genomes of R. × pulchrum, respectively. Five genes (nad1, nad2, nad4, nad7, and rps3) were potentially useful molecular markers. The R. × pulchrum mitochondrial genome collinear alignment among five species of the Ericaceae showed that the mitochondrial genomes of these related species have a high degree of homology with R. × pulchrum in this gene region, and the most conservative genes were trnC-GCA, trnD-GUC, trnM-CAU, trnN-GUU, trnY-GUA, atp4, nad4, nad2, nad5, ccmC, and rrn26. The phylogenetic trees of mitochondrial genome showed that R. simsii was a sister to R. × pulchrum. The results verified that there was gene rearrangement between R. × pulchrum and R. simsii mitochondrial genomes. The codon usage bias of 10 Ericaceae mitochondrial genes and 7 Rhododendron chloroplast genes were influenced by mutation, while other genes codon usages had undergone selection. The study identified 13 homologous fragments containing gene sequences between the chloroplast and mitochondrial genomes of R. × pulchrum. Overall, our results illustrate the organelle genome information could explain the phylogenetics of plants and could be used to develop molecular markers and genetic evolution. Our study will facilitate the study of population genetics and evolution in Rhododendron and other genera in Ericaceae.

Rearrangement and domestication as drivers of Rosaceae mitogenome plasticity

Background

The mitochondrion is an important cellular component in plants and that functions in producing vital energy for the cell. However, the evolution and structure of mitochondrial genomes (mitogenomes) remain unclear in the Rosaceae family. In this study, we assembled 34 Rosaceae mitogenomes and characterized genome variation, rearrangement rate, and selection signal variation within these mitogenomes.

Results

Comparative analysis of six genera from the Amygdaloideae and five from the Rosoideae subfamilies of Rosaceae revealed that three protein-coding genes were absent from the mitogenomes of five Rosoideae genera. Positive correlations between genome size and repeat content were identified in 38 Rosaceae mitogenomes. Twenty repeats with high recombination frequency (> 50%) provided evidence for predominant substoichiometric conformation of the mitogenomes. Variations in rearrangement rates were identified between eleven genera, and within the Pyrus, Malus, Prunus, and Fragaria genera. Based on population data, phylogenetic inferences from Pyrus mitogenomes supported two distinct maternal lineages of Asian cultivated pears. A Pyrus-specific deletion (DEL-D) in selective sweeps was identified based on the assembled genomes and population data. After the DEL-D sequence fragments originally arose, they may have experienced a subsequent doubling event via homologous recombination and sequence transfer in the Amygdaloideae; afterwards, this variant sequence may have significantly expanded to cultivated groups, thereby improving adaptation during the domestication process.

Conclusions

This study characterizes the variations in gene content, genome size, rearrangement rate, and the impact of domestication in Rosaceae mitogenomes and provides insights into their structural variation patterns and phylogenetic relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01383-3.

Deciphering the Multi-Chromosomal Mitochondrial Genome of Populus simonii

Mitochondria, inherited maternally, are energy metabolism organelles that generate most of the chemical energy needed to power cellular various biochemical reactions. Deciphering mitochondrial genome (mitogenome) is important for elucidating vital activities of species. The complete chloroplast (cp) and nuclear genome sequences of Populus simonii (P. simonii) have been reported, but there has been little progress in its mitogenome. Here, we assemble the complete P. simonii mitogenome into three circular-mapping molecules (lengths 312.5, 283, and 186 kb) with the total length of 781.5 kb. All three molecules of the P. simonii mitogenome had protein-coding capability. Whole-genome alignment analyses of four Populus species revealed the fission of poplar mitogenome in P. simonii. Comparative repeat analyses of four Populus mitogenomes showed that there were no repeats longer than 350 bp in Populus mitogenomes, contributing to the stability of genome sizes and gene contents in the genus Populus. As the first reported multi-circular mitogenome in Populus, this study of P. simonii mitogenome are imperative for better elucidating their biological functions, replication and recombination mechanisms, and their unique evolutionary trajectories in Populus.

Characterization and Phylogenetic Analyses of the Complete Mitochondrial Genome of Sugarcane (Saccharum spp. Hybrids) Line A1

Modern sugarcane cultivars are highly polyploid with complex nuclear genomic genetic background, while their mitochondrion (mt) genomes are much simpler, smaller and more manageable and could provide useful phylogenetic information. In this study, the mt genome of a modern commercial cultivar A1 was sequenced via Illumina Hiseq XTen and PacBio Sequel platform. The assembled and annotated mitochondrial genomes of A1 were composed of two circular DNA molecules, one large and one small, which were named Chromosome 1 and Chromosome 2. The two distinct circular chromosomes of mitogenome construct is consisted with other sugarcane cultivars i.e., Saccharum officinarum Khon Kaen 3 and Saccharum spp. hybrids ROC22 and FN15. The Chromosome 1 of A1 mitogenome is 300,822 bp in length with the GC content of 43.94%, and 7.14% of Chromosome 1 sequences (21,468 nucleotides) are protein coding genes (PCGs) while 92.86% (279,354 nucleotides) are intergenic region. The length of Chromosome 2 is 144,744 bp with the GC content of 43.57%, and 8.20% of Chromosome 2 sequences (11,865 nucleotides) are PCGs while 91.80% (132,879 nucleotides) are intergenic region. A total of 43 genes are located on Chromosome 1, which contains 22 PCGs (six nad genes, four rps genes, four atp genes, three ccm genes, three cox genes, one mat gene and one mtt gene) and 21 non-coding genes including 15 tRNAs and 6 rRNAs. Chromosome 2 includes 18 genes in total, which contains 13 PCGs (four nad genes, three rps genes, two atp genes, one ccm gene, one cob gene, one cox gene and one rpl gene) and five non-coding genes (tRNA genes). Analysis of codon usage of 35 PCGs showed that codon ending in A/U was preferred. Investigation of gene composition indicated that the types and copy numbers of CDS genes, tRNAs and rRNAs of A1 and FN15 were identical. The cox1 gene has two copies and the trnP gene has one copy in A1, FN15 and ROC22 three lines, while there is only one copy of cox1 and two copies of trnP in S. officinarum Khon Kaen 3. In addition, S. officinarum Khon Kaen 3 have no nad1 gene and rps7 gene. 100 sequence repeats, 38 SSRs and 444 RNA editing sites in A1 mt genome were detected. Moreover, the maximum likelihood phylogenetic analysis found that A1 were more closely related to S. spp. hybrid (ROC22 and FN15) and S. officinarum (Khon Kaen 3). Herein, the complete mt genome of A1 will provide essential DNA molecular information for further phylogenetic and evolutionary analysis for Saccharum and Poaceae.

Mitochondrial Phylogenomics of Fagales Provides Insights Into Plant Mitogenome Mosaic Evolution

Fagales are an order of woody plants and comprise more than 1,100 species, most of which produce economically important timbers, nuts, and fruits. Their nuclear and plastid genomes are well-sequenced and provided valuable resources to study their phylogeny, breeding, resistance, etc. However, little is known about the mitochondrial genomes (mitogenomes), which hinder a full understanding of their genome evolution. In this study, we assembled complete mitogenomes of 23 species, covering five of the seven families of Fagales. These mitogenomes had similar gene sets but varied 2.4 times in size. The mitochondrial genes were highly conserved, and their capacity in phylogeny was challenging. The mitogenomic structure was extremely dynamic, and synteny among species was poor. Further analyses of the Fagales mitogenomes revealed extremely mosaic characteristics, with horizontal transfer (HGT)-like sequences from almost all seed plant taxa and even mitoviruses. The largest mitogenome, Carpinus cordata, did not have large amounts of specific sequences but instead contained a high proportion of sequences homologous to other Fagales. Independent and unequal transfers of third-party DNA, including nuclear genome and other resources, may partially account for the HGT-like fragments and unbalanced size expansions observed in Fagales mitogenomes. Supporting this, a mitochondrial plasmid-like of nuclear origin was found in Carpinus. Overall, we deciphered the last genetic materials of Fagales, and our large-scale analyses provide new insights into plant mitogenome evolution and size variation.

A complete mitochondrial genome for fragrant Chinese rosewood (Dalbergia odorifera, Fabaceae) with comparative analyses of genome structure and intergenomic sequence transfers

BACKGROUND

Dalbergia odorifera is an economically and culturally important species in the Fabaceae because of the high-quality lumber and traditional Chinese medicines made from this plant, however, overexploitation has increased the scarcity of D. odorifera. Given the rarity and the multiple uses of this species, it is important to expand the genomic resources for utilizing in applications such as tracking illegal logging, determining effective population size of wild stands, delineating pedigrees in marker assisted breeding programs, and resolving gene networks in functional genomics studies. Even the nuclear and chloroplast genomes have been published for D. odorifera, the complete mitochondrial genome has not been assembled or assessed for sequence transfer to other genomic compartments until now. Such work is essential in understanding structural and functional genome evolution in a lineage (Fabaceae) with frequent intergenomic sequence transfers.

RESULTS

We integrated Illumina short-reads and PacBio CLR long-reads to assemble and annotate the complete mitochondrial genome of D. odorifera. The mitochondrial genome was organized as a single circular structure of 435 Kb in length containing 33 protein coding genes, 4 rRNA and 17 tRNA genes. Nearly 4.0% (17,386 bp) of the genome was annotated as repetitive DNA. From the sequence transfer analysis, it was found that 114 Kb of DNA originating from the mitochondrial genome has been transferred to the nuclear genome, with most of the transfer events having taken place relatively recently. The high frequency of sequence transfers from the mitochondria to the nuclear genome was similar to that of sequence transfer from the chloroplast to the nuclear genome.

CONCLUSION

For the first-time, the complete mitochondrial genome of D. odorifera was assembled in this study, which will provide a baseline resource in understanding genomic evolution in the highly specious Fabaceae. In particular, the assessment of intergenomic sequence transfer suggests that transfers have been common and recent indicating a possible role in environmental adaptation as has been found in other lineages. The high turnover rate of genomic colinearly and large differences in mitochondrial genome size found in the comparative analyses herein providing evidence for the rapid evolution of mitochondrial genome structure compared to chloroplasts in Faboideae. While phylogenetic analyses using functional genes indicate that mitochondrial genes are very slowly evolving compared to chloroplast genes.

Rapid shifts in mitochondrial tRNA import in a plant lineage with extensive mitochondrial tRNA gene loss

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence.

Mitochondrial RNase H1 activity regulates R-loop homeostasis to maintain genome integrity and enable early embryogenesis in Arabidopsis

Plant mitochondrial genomes undergo frequent homologous recombination (HR). Ectopic HR activity is inhibited by the HR surveillance pathway, but the underlying regulatory mechanism is unclear. Here, we show that the mitochondrial RNase H1 AtRNH1B impairs the formation of RNA:DNA hybrids (R-loops) and participates in the HR surveillance pathway in Arabidopsis thaliana. AtRNH1B suppresses ectopic HR at intermediate-sized repeats (IRs) and thus maintains mitochondrial DNA (mtDNA) replication. The RNase H1 AtRNH1C is restricted to the chloroplast; however, when cells lack AtRNH1B, transport of chloroplast AtRNH1C into the mitochondria secures HR surveillance, thus ensuring the integrity of the mitochondrial genome and allowing embryogenesis to proceed. HR surveillance is further regulated by the single-stranded DNA-binding protein ORGANELLAR SINGLE-STRANDED DNA BINDING PROTEIN1 (OSB1), which decreases the formation of R-loops. This study uncovers a facultative dual targeting mechanism between organelles and sheds light on the roles of RNase H1 in organellar genome maintenance and embryogenesis.

This study clarifies the function of mitochondrial RNase H1 in genome stability and early embryogenesis in plants, and shows that mitochondrial R-loops are involved in homologous recombination surveillance of mtDNA. Facultative re-targeting of the chloroplast RNase H1 protein to mitochondria, in response to cellular conditions, can help guarantee mitochondrial RNase H1 activity.

Detecting de novo mitochondrial mutations in angiosperms with highly divergent evolutionary rates

Although plant mitochondrial genomes typically show low rates of sequence evolution, levels of divergence in certain angiosperm lineages suggest anomalously high mitochondrial mutation rates. However, de novo mutations have never been directly analyzed in such lineages. Recent advances in high-fidelity DNA sequencing technologies have enabled detection of mitochondrial mutations when still present at low heteroplasmic frequencies. To date, these approaches have only been performed on a single plant species (Arabidopsis thaliana). Here, we apply a high-fidelity technique (Duplex Sequencing) to multiple angiosperms from the genus Silene, which exhibits extreme heterogeneity in rates of mitochondrial sequence evolution among close relatives. Consistent with phylogenetic evidence, we found that Silene latifolia maintains low mitochondrial variant frequencies that are comparable with previous measurements in Arabidopsis. Silene noctiflora also exhibited low variant frequencies despite high levels of historical sequence divergence, which supports other lines of evidence that this species has reverted to lower mitochondrial mutation rates after a past episode of acceleration. In contrast, S. conica showed much higher variant frequencies in mitochondrial (but not in plastid) DNA, consistent with an ongoing bout of elevated mitochondrial mutation rates. Moreover, we found an altered mutational spectrum in S. conica heavily biased towards AT→GC transitions. We also observed an unusually low number of mitochondrial genome copies per cell in S. conica, potentially pointing to reduced opportunities for homologous recombination to accurately repair mismatches in this species. Overall, these results suggest that historical fluctuations in mutation rates are driving extreme variation in rates of plant mitochondrial sequence evolution.

Cyanophora paradoxa mitochondrial tRNAs play a double game

Present-day mitochondria derive from a single endosymbiosis of an α-proteobacterium into a proto-eukaryotic cell. Since this monophyletic event, mitochondria have evolved considerably, and unique traits have been independently acquired in the different eukaryotic kingdoms. Mitochondrial genome expression and RNA metabolism have diverged greatly. Here, Cyanophora paradoxa, a freshwater alga considered as a living fossil among photosynthetic organisms, represents an exciting model for studying the evolution of mitochondrial gene expression. As expected, fully mature tRNAs are released from primary transcripts to function in mitochondrial translation. We also show that these tRNAs take part in an mRNA processing punctuation mechanism in a non-conventional manner, leading to mRNA-tRNA hybrids with a CCA triplet at their 3'-extremities. In this case, tRNAs are probably used as stabilizing structures impeding the degradation of mRNA by exonucleases. From our data we propose that the present-day tRNA-like elements (t-elements) found at the 3'-terminals of mitochondrial mRNAs in land plants originate from true tRNAs like those observed in the mitochondria of this basal photosynthetic glaucophyte.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

Plant Organellar DNA Polymerases Evolved Multifunctionality through the Acquisition of Novel Amino Acid Insertions

The majority of DNA polymerases (DNAPs) are specialized enzymes with specific roles in DNA replication, translesion DNA synthesis (TLS), or DNA repair. The enzymatic characteristics to perform accurate DNA replication are in apparent contradiction with TLS or DNA repair abilities. For instance, replicative DNAPs incorporate nucleotides with high fidelity and processivity, whereas TLS DNAPs are low-fidelity polymerases with distributive nucleotide incorporation. Plant organelles (mitochondria and chloroplast) are replicated by family-A DNA polymerases that are both replicative and TLS DNAPs. Furthermore, plant organellar DNA polymerases from the plant model Arabidopsis thaliana (AtPOLIs) execute repair of double-stranded breaks by microhomology-mediated end-joining and perform Base Excision Repair (BER) using lyase and strand-displacement activities. AtPOLIs harbor three unique insertions in their polymerization domain that are associated with TLS, microhomology-mediated end-joining (MMEJ), strand-displacement, and lyase activities. We postulate that AtPOLIs are able to execute those different functions through the acquisition of these novel amino acid insertions, making them multifunctional enzymes able to participate in DNA replication and DNA repair.

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 Inheritance in Phytopathogenic Fungi-Everything Is Known, or Is It?

Mitochondria are important organelles in eukaryotes that provide energy for cellular processes. Their function is highly conserved and depends on the expression of nuclear encoded genes and genes encoded in the organellar genome. Mitochondrial DNA replication is independent of the replication control of nuclear DNA and as such, mitochondria may behave as selfish elements, so they need to be controlled, maintained and reliably inherited to progeny. Phytopathogenic fungi meet with special environmental challenges within the plant host that might depend on and influence mitochondrial functions and services. We find that this topic is basically unexplored in the literature, so this review largely depends on work published in other systems. In trying to answer elemental questions on mitochondrial functioning, we aim to introduce the aspect of mitochondrial functions and services to the study of plant-microbe-interactions and stimulate phytopathologists to consider research on this important organelle in their future projects.

Genome skimming and microsatellite analysis reveal contrasting patterns of genetic diversity in a rare sandhill endemic (Erysimum teretifolium, Brassicaceae)

Barriers between islands often inhibit gene flow creating patterns of isolation by distance. In island species, the majority of genetic diversity should be distributed among isolated populations. However, a self-incompatible mating system leads to higher genetic variation within populations and very little between-population subdivision. We examine these two contrasting predictions in Erysimum teretifolium, a rare self-incompatible plant endemic to island-like sandhill habitats in Santa Cruz County, California. We used genome skimming and nuclear microsatellites to assess the distribution of genetic diversity within and among eight of the 13 remaining populations. Phylogenetic analyses of the chloroplast genomes revealed a deep separation of three of the eight populations. The nuclear ribosomal DNA cistron showed no genetic subdivision. Nuclear microsatellites suggest 83% of genetic variation resides within populations. Despite this, 18 of 28 between-population comparisons exhibited significant population structure (mean F_ST = 0.153). No isolation by distance existed among all populations, however when one outlier population was removed from the analysis due to uncertain provenance, significant isolation by distance emerged (r² = 0.5611, p = 0.005). Population census size did not correlate with allelic richness as predicted on islands. Bayesian population assignment detected six genetic groupings with substantial admixture. Unique genetic clusters were concentrated at the periphery of the species’ range. Since the overall distribution of nuclear genetic diversity reflects E. tereifolium’s self-incompatible mating system, the vast majority of genetic variation could be sampled within any individual population. Yet, the chloroplast genome results suggest a deep split and some of the nuclear microsatellite analyses indicate some island-like patterns of genetic diversity. Restoration efforts intending to maximize genetic variation should include representatives from both lineages of the chloroplast genome and, for maximum nuclear genetic diversity, should include representatives of the smaller, peripheral populations.