Frequent, phylogenetically local horizontal transfer of the cox1 group I Intron in flowering plant mitochondria

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica

MAIN CONCLUSION

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BACKGROUND

Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.

RESULTS

PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.

CONCLUSIONS

We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

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.

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.

Evolution of cox2 introns in angiosperm mitochondria and efficient splicing of an elongated cox2i691 intron

In angiosperms, the mitochondrial cox2 gene harbors up to two introns, commonly referred to as cox2i373 and cox2i691. We studied the cox2 from 222 fully-sequenced mitogenomes from 30 angiosperm orders and analyzed the evolution of their introns. Unlike cox2i373, cox2i691 shows a distribution among plants that is shaped by frequent intron loss events driven by localized retroprocessing. In addition, cox2i691 exhibits sporadic elongations, frequently in domain IV of introns. Such elongations are poorly related to repeat content and two of them showed the presence of LINE transposons, suggesting that increasing intron size is very likely due to nuclear intracelular DNA transfer followed by incorporation into the mitochondrial DNA. Surprisingly, we found that cox2i691 is erroneously annotated as absent in 30 mitogenomes deposited in public databases. Although each of the cox2 introns is ∼1.5 kb in length, a cox2i691 of 4.2 kb has been reported in Acacia ligulata (Fabaceae). It is still unclear whether its unusual length is due to a trans-splicing arrangement or the loss of functionality of the interrupted cox2. Through analyzing short-read RNA sequencing of Acacia with a multi-step computational strategy, we found that the Acacia cox2 is functional and its long intron is spliced in cis in a very efficient manner despite its length.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

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.

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.

Plant mitochondria - past, present and future

The study of plant mitochondria started in earnest around 1950 with the first isolations of mitochondria from animal and plant tissues. The first 35 years were spent establishing the basic properties of plant mitochondria and plant respiration using biochemical and physiological approaches. A number of unique properties (compared to mammalian mitochondria) were observed: (i) the ability to oxidize malate, glycine and cytosolic NAD(P)H at high rates; (ii) the partial insensitivity to rotenone, which turned out to be due to the presence of a second NADH dehydrogenase on the inner surface of the inner mitochondrial membrane in addition to the classical Complex I NADH dehydrogenase; and (iii) the partial insensitivity to cyanide, which turned out to be due to an alternative oxidase, which is also located on the inner surface of the inner mitochondrial membrane, in addition to the classical Complex IV, cytochrome oxidase. With the appearance of molecular biology methods around 1985, followed by genomics, further unique properties were discovered: (iv) plant mitochondrial DNA (mtDNA) is 10-600 times larger than the mammalian mtDNA, yet it only contains approximately 50% more genes; (v) plant mtDNA has kept the standard genetic code, and it has a low divergence rate with respect to point mutations, but a high recombinatorial activity; (vi) mitochondrial mRNA maturation includes a uniquely complex set of activities for processing, splicing and editing (at hundreds of sites); (vii) recombination in mtDNA creates novel reading frames that can produce male sterility; and (viii) plant mitochondria have a large proteome with 2000-3000 different proteins containing many unique proteins such as 200-300 pentatricopeptide repeat proteins. We describe the present and fairly detailed picture of the structure and function of plant mitochondria and how the unique properties make their metabolism more flexible allowing them to be involved in many diverse processes in the plant cell, such as photosynthesis, photorespiration, CAM and C4 metabolism, heat production, temperature control, stress resistance mechanisms, programmed cell death and genomic evolution. However, it is still a challenge to understand how the regulation of metabolism and mtDNA expression works at the cellular level and how retrograde signaling from the mitochondria coordinates all those processes.

Shifts from cis-to trans-splicing of five mitochondrial introns in Tolypanthus maclurei

Shifts from cis-to trans-splicing of mitochondrial introns tend to correlate with relative genome rearrangement rates during vascular plant evolution, as is particularly apparent in some lineages of gymnosperms. However, although many angiosperms have also relatively high mitogenomic rearrangement rates, very few cis-to trans-splicing shifts except for five trans-spliced introns shared in seed plants have been reported. In this study, we sequenced and characterized the mitogenome of Tolypanthus maclurei, a hemiparasitic plant from the family Loranthaceae (Santalales). The mitogenome was assembled into a circular chromosome of 256,961 bp long, relatively small compared with its relatives from Santalales. It possessed a gene content of typical angiosperm mitogenomes, including 33 protein-coding genes, three rRNA genes and ten tRNA genes. Plastid-derived DNA fragments took up 9.1% of the mitogenome. The mitogenome contained one group I intron (cox1i729) and 23 group II introns. We found shifts from cis-to trans-splicing of five additional introns in its mitogenome, of which two are specific in T. maclurei. Moreover, atp1 is a chimeric gene and phylogenetic analysis indicated that a 356 bp region near the 3′ end of atp1 of T. maclurei was acquired from Lamiales via horizontal gene transfer. Our results suggest that shifts to trans-splicing of mitochondrial introns may not be uncommon among angiosperms.

Horizontal Gene Transfers in Plants

In plants, as in all eukaryotes, the vertical transmission of genetic information through reproduction ensures the maintenance of the integrity of species. However, many reports over the past few years have clearly shown that horizontal gene transfers, referred to as HGTs (the interspecific transmission of genetic information across reproductive barriers) are very common in nature and concern all living organisms including plants. The advent of next-generation sequencing technologies (NGS) has opened new perspectives for the study of HGTs through comparative genomic approaches. In this review, we provide an up-to-date view of our current knowledge of HGTs in plants.

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.

Mitochondrial Fostering: The Mitochondrial Genome May Play a Role in Plant Orphan Gene Evolution

Plant mitochondrial genomes exhibit unique evolutionary patterns. They have a high rearrangement but low mutation rate, and a large size. Based on massive mitochondrial DNA transfers to the nucleus as well as the mitochondrial unique evolutionary traits, we propose a "Mitochondrial Fostering" theory where the organelle genome plays an integral role in the arrival and development of orphan genes (genes with no homologs in other lineages). Two approaches were used to test this theory: (1) bioinformatic analysis of nuclear mitochondrial DNA (Numts: mitochondrial originating DNA that migrated to the nucleus) at the genome level, and (2) bioinformatic analysis of particular orphan sequences present in both the mitochondrial genome and the nuclear genome of Arabidopsis thaliana. One study example is given about one orphan sequence that codes for two unique orphan genes: one in the mitochondrial genome and another one in the nuclear genome. DNA alignments show regions of this A. thaliana orphan sequence exist scattered throughout other land plant mitochondrial genomes. This is consistent with the high recombination rates of mitochondrial genomes in land plants. This may also enable the creation of novel coding sequences within the orphan loci, which can then be transferred to the nuclear genome and become exposed to new evolutionary pressures. Our study also reveals a high correlation between the amount of mitochondrial DNA transferred to the nuclear genome and the number of orphan genes in land plants. All the data suggests the mitochondrial genome may play a role in nuclear orphan gene evolution in land plants.

Episodic and guanine-cytosine-biased bursts of intragenomic and interspecific synonymous divergence in Ajugoideae (Lamiaceae) mitogenomes

Synonymous substitution rates in plant mitochondrial genomes vary by orders of magnitude among species, whereas synonymous rates among genes within a genome are generally consistent. Exceptionally, genes within the Ajuga reptans (Lamiaceae) mitochondrial genome exhibit unprecedented intragenomic heterogeneity in synonymous sequence divergence, but the biological mechanisms underlying this rate variation remain unclear. We tracked the origin and evolutionary trajectory of mitochondrial rate variations by dense sampling in Ajugoideae and found differences in the timing and magnitude of rate acceleration for particular genes. The most divergent genes accelerated earlier, retained a high rate across Ajugoideae, and are generally devoid of RNA editing, whereas moderately diverged genes accelerated later and retained relatively higher RNA editing frequency. The acceleration of mutation rates correlates with increased guanine-cytosine (GC) content, suggesting a key role for GC-biased gene conversion and/or repair after the breakage of ancestral gene clusters.

Extensive variation in nucleotide substitution rate and gene/intron loss in mitochondrial genomes of Pelargonium

Geraniaceae organelle genomes have been shown to exhibit several highly unusual features compared to most other photosynthetic angiosperms. This includes massively rearranged plastomes with considerable size variation, extensive gene and intron loss, accelerated rates of nucleotide substitutions in both mitogenomes and plastomes, and biparental inheritance and cytonuclear incompatibility of the plastome. Most previous studies have focused on plastome evolution with mitogenome comparisons limited to only a few taxa or genes. In this study, mitogenomes and transcriptomes were examined for 27 species of Geraniales, including 13 species of Pelargonium. Extensive gene and intron losses were detected across the Geraniales with Pelargonium representing the most gene depauperate lineage in the family. Plotting these events on the Geraniaceae phylogenetic tree showed that gene losses occurred multiple times, whereas intron losses more closely reflected the relationships among taxa. In addition, P. australe acquired an intron by horizontal transfer. Comparisons of nucleotide substitution rates in Pelargonium showed that synonymous changes in nuclear genes were much lower than in mitochondrial genes. This is in contrast to the previously published studies that indicated that nuclear genes have 16 fold higher rates than mitochondrial genes across angiosperms. Elevated synonymous substitutions occurred for each mitochondrial gene in Pelargonium with the highest values 783 and 324 times higher than outgroups and other Geraniaceae, respectively. Pelargonium is one of four unrelated genera of angiosperms (Ajuga, Plantago and Silene) that have experienced highly accelerated nucleotide substitutions in mitogenomes. It is distinct from most angiosperms in also having elevated substitution rates in plastid genes but the cause of rate accelerations in Pelargonium plastomes and mitogenomes may be different.

The insertion of a mitochondrial selfish element into the nuclear genome and its consequences

Homing endonucleases (HE) are enzymes capable of cutting DNA at highly specific target sequences, the repair of the generated double-strand break resulting in the insertion of the HE-encoding gene ("homing" mechanism). HEs are present in all three domains of life and viruses; in eukaryotes, they are mostly found in the genomes of mitochondria and chloroplasts, as well as nuclear ribosomal RNAs. We here report the case of a HE that accidentally integrated into a telomeric region of the nuclear genome of the fungal maize pathogen Ustilago maydis. We show that the gene has a mitochondrial origin, but its original copy is absent from the U. maydis mitochondrial genome, suggesting a subsequent loss or a horizontal transfer from a different species. The telomeric HE underwent mutations in its active site and lost its original start codon. A potential other start codon was retained downstream, but we did not detect any significant transcription of the newly created open reading frame, suggesting that the inserted gene is not functional. Besides, the insertion site is located in a putative RecQ helicase gene, truncating the C-terminal domain of the protein. The truncated helicase is expressed during infection of the host, together with other homologous telomeric helicases. This unusual mutational event altered two genes: The integrated HE gene subsequently lost its homing activity, while its insertion created a truncated version of an existing gene, possibly altering its function. As the insertion is absent in other field isolates, suggesting that it is recent, the U. maydis 521 reference strain offers a snapshot of this singular mutational event.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

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.

Variation in protein gene and intron content among land plant mitogenomes

The functional content of the mitochondrial genome (mitogenome) is highly diverse across eukaryotes. Among land plants, our understanding of the variation in mitochondrial gene and intron content is improving from concerted efforts to densely sample mitogenomes from diverse land plants. Here I review the current state of knowledge regarding the diversity in content of protein genes and introns in the mitogenomes of all major land plant lineages. Mitochondrial protein gene content is largely conserved among mosses and liverworts, but it varies substantially among and within other land plant lineages due to convergent losses of genes encoding ribosomal proteins and, to a lesser extent, genes for proteins involved in cytochrome c maturation and oxidative phosphorylation. Mitochondrial intron content is fairly stable within each major land plant lineage, but highly variable among lineages, resulting from occasional gains and many convergent losses over time. Trans-splicing has evolved dozens of times in various vascular plant lineages, particularly those with relatively higher rates of mitogenomic rearrangement. Across eukaryotes, mitochondrial protein gene and intron content has been shaped massive convergent evolution.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Mitochondrial Genome of Fagopyrum esculentum and the Genetic Diversity of Extranuclear Genomes in Buckwheat

Fagopyrum esculentum (common buckwheat) is an important agricultural non-cereal grain plant. Despite extensive genetic studies, the information on its mitochondrial genome is still lacking. Using long reads generated by single-molecule real-time technology coupled with circular consensus sequencing (CCS) protocol, we assembled the buckwheat mitochondrial genome and detected that its prevalent form consists of 10 circular chromosomes with a total length of 404 Kb. In order to confirm the presence of a multipartite structure, we developed a new targeted assembly tool capable of processing long reads. The mitogenome contains all genes typical for plant mitochondrial genomes and long inserts of plastid origin (~6.4% of the total mitogenome length). Using this new information, we characterized the genetic diversity of mitochondrial and plastid genomes in 11 buckwheat cultivars compared with the ancestral subspecies, F. esculentum ssp. ancestrale. We found it to be surprisingly low within cultivars: Only three to six variations in the mitogenome and one to two in the plastid genome. In contrast, the divergence with F. esculentum ssp. ancestrale is much higher: 220 positions differ in the mitochondrial genome and 159 in the plastid genome. The SNPs in the plastid genome are enriched in non-synonymous substitutions, in particular in the genes involved in photosynthesis: psbA, psbC, and psbH. This presumably reflects the selection for the increased photosynthesis efficiency as a part of the buckwheat breeding program.

Evolution of Fusarium tricinctum and Fusarium avenaceum mitochondrial genomes is driven by mobility of introns and of a new type of palindromic microsatellite repeats

Background

Increased contamination of European and Asian wheat and barley crops with “emerging” mycotoxins such as enniatins or beauvericin, produced by Fusarium avenaceum and Fusarium tricinctum, suggest that these phylogenetically close species could be involved in future food-safety crises.

Results

The mitochondrial genomes of F. tricinctum strain INRA104 and F. avenaceum strain FaLH27 have been annotated. A comparative analysis was carried out then extended to a set of 25 wild strains. Results show that they constitute two distinct species, easily distinguished by their mitochondrial sequences. The mitochondrial genetic variability is mainly located within the intergenic regions. Marks of variations show they have evolved (i) by Single Nucleotide Polymorphisms (SNPs), (ii) by length variations mediated by insertion/deletion sequences (Indels), and (iii) by length mutations generated by DNA sliding events occurring in mononucleotide (A)_n or (T)_n microsatellite type sequences arranged in a peculiar palindromic organization. The optionality of these palindromes between both species argues for their mobility. The presence of Indels and SNPs in palindrome neighbouring regions suggests their involvement in these observed variations. Moreover, the intraspecific and interspecific variations in the presence/absence of group I introns suggest a high mobility, resulting from several events of gain and loss during short evolution periods. Phylogenetic analyses of intron orthologous sequences suggest that most introns could have originated from lateral transfers from phylogenetically close or distant species belonging to various Ascomycota genera and even to the Basidiomycota fungal division.

Conclusions

Mitochondrial genome evolution between F. tricinctum and F. avenaceum is mostly driven by two types of mobile genetic elements, implicated in genome polymorphism. The first one is represented by group I introns. Indeed, both genomes harbour optional (inter- or intra-specifically) group I introns, all carrying putatively functional hegs, arguing for a high mobility of these introns during short evolution periods. The gain events were shown to involve, for most of them, lateral transfers between phylogenetically distant species. This study has also revealed a new type of mobile genetic element constituted by a palindromic arrangement of (A) n and (T) n microsatellite sequences whose presence was related to occurrence of SNPs and Indels in the neighbouring regions.

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.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

Ancient Mitochondrial Gene Transfer between Fungi and the Orchids

The mitochondrial genomes (mitogenomes) of plants are known to incorporate and accumulate DNA from intra- and extracellular donors. Despite the intimate relationships formed between flowing plants (angiosperms) and fungi, lengthy fungal-like sequence has not been identified in angiosperm mitogenomes to date. Here, we present multiple lines of evidence documenting horizontal gene transfer (HGT) between the mitogenomes of fungi and the ancestors of the orchids, plants that are obligate parasites of fungi during their early development. We show that the ancestor of the orchids acquired an ∼270-bp fungal mitogenomic region containing three transfer RNA genes. We propose that the short HGT was later replaced by a second HGT event transferring >8 kb and 14 genes from a fungal mitogenome to that of the ancestor of the largest orchid subfamily, Epidendroideae. Our results represent the first evidence of genomic-scale HGT between fungal and angiosperm mitogenomes and demonstrate that the length intergenic spacer regions of angiosperm mitogenomes can effectively fossilize the genomic remains of ancient, nonplant organisms.

Comprehensive analyses of ZFP gene family and characterization of expression profiles during plant hormone response in cotton

BACKGROUND

Zinc finger proteins (ZFPs) containing only a single zinc finger domain play important roles in the regulation of plant growth and development, as well as in biotic and abiotic stress responses. To date, the evolutionary history and functions of the ZFP gene family have not been identified in cotton.

RESULTS

In this paper, we identified 29 ZFP genes in Gossypium hirsutum. This gene family was divided into seven subfamilies, 22 of which were distributed over 17 chromosomes. Bioinformatic analysis revealed that 20 GhZFP genes originated from whole genome duplications and two originated from dispersed duplication events, indicating that whole genome duplication is the main force in the expansion of the GhZFP gene family. Most GhZFP8 subfamily genes, except for GhZFP8-3, were highly expressed during fiber cell growth, and were induced by brassinosteroids in vitro. Furthermore, we found that a large number of GhZFP genes contained gibberellic acid responsive elements, auxin responsive elements, and E-box elements in their promoter regions. Exogenous application of these hormones significantly stimulated the expression of these genes.

CONCLUSIONS

Our findings reveal that GhZFP8 genes are involved in cotton fiber development and widely induced by auxin, gibberellin and BR, which provides a foundation for the identification of more downstream genes with potential roles in phytohormone stimuli, and a basis for breeding better cotton varieties in the future.

Elucidating genomic patterns and recombination events in plant cybrid mitochondria

KEY MESSAGE

Cybrid plant mitochondria undergo homologous recombination, mainly BIR, keep a single allele for each gene, and maintain exclusive sequences of each parent and a single copy of the homologous regions. The maintenance of a dynamic equilibrium between the mitochondrial and nuclear genomes requires continuous communication and a high level of compatibility between them, so that alterations in one genetic compartment need adjustments in the other. The co-evolution of nuclear and mitochondrial genomes has been poorly studied, even though the consequences and effects of this interaction are highly relevant for human health, as well as for crop improvement programs and for genetic engineering. The mitochondria of plants represent an excellent system to understand the mechanisms of genomic rearrangements, chimeric gene formation, incompatibility between nucleus and cytoplasm, and horizontal gene transfer. We carried out detailed analyses of the mtDNA of a repeated cybrid between the solanaceae Nicotiana tabacum and Hyoscyamus niger. The mtDNA of the cybrid was intermediate between the size of the parental mtDNAs and the sum of them. Noticeably, most of the homologous sequences inherited from both parents were lost. In contrast, the majority of the sequences exclusive of a single parent were maintained. The mitochondrial gene content included a majority of N. tabacum derived genes, but also chimeric, two-parent derived, and H. niger-derived genes in a tobacco nuclear background. Any of these alterations in the gene content could be the cause of CMS in the cybrid. The parental mtDNAs interacted through 28 homologous recombination events and a single case of illegitimate recombination. Three main homologous recombination mechanisms were recognized in the cybrid mitochondria. Break induced replication (BIR) pathway was the most frequent. We propose that BIR could be one of the mechanisms responsible for the loss of the majority of the repeated regions derived from H. niger.

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.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5732-z) contains supplementary material, which is available to authorized users.

Evolution and Diversification of Kiwifruit Mitogenomes through Extensive Whole-Genome Rearrangement and Mosaic Loss of Intergenic Sequences in a Highly Variable Region

Angiosperm mitochondrial genomes (mitogenomes) are notable for their extreme diversity in both size and structure. However, our current understanding of this diversity is limited, and the underlying mechanism contributing to this diversity remains unclear. Here, we completely assembled and compared the mitogenomes of three kiwifruit (Actinidia) species, which represent an early divergent lineage in asterids. We found conserved gene content and fewer genomic repeats, particularly large repeats (>1 kb), in the three mitogenomes. However, sequence transfers such as intracellular events are variable and dynamic, in which both ancestral shared and recently species-specific events as well as complicated transfers of two plastid-derived sequences into the nucleus through the mitogenomic bridge were detected. We identified extensive whole-genome rearrangements among kiwifruit mitogenomes and found a highly variable V region in which fragmentation and frequent mosaic loss of intergenic sequences occurred, resulting in greatly interspecific variations. One example is the fragmentation of the V region into two regions, V1 and V2, giving rise to the two mitochondrial chromosomes of Actinidia chinensis. Finally, we compared the kiwifruit mitogenomes with those of other asterids to characterize their overall mitogenomic diversity, which identified frequent gain/loss of genes/introns across lineages. In addition to repeat-mediated recombination and import-driven hypothesis of genome size expansion reported in previous studies, our results highlight a pattern of dynamic structural variation in plant mitogenomes through global genomic rearrangements and species-specific fragmentation and mosaic loss of intergenic sequences in highly variable regions on the basis of a relatively large ancestral mitogenome.

Horizontally-acquired genetic elements in the mitochondrial genome of a centrohelid Marophrys sp. SRT127

Mitochondrial genomes exhibit diverse features among eukaryotes in the aspect of gene content, genome structure, and the mobile genetic elements such as introns and plasmids. Although the number of published mitochondrial genomes is increasing at tremendous speed, those of several lineages remain unexplored. Here, we sequenced the complete mitochondrial genome of a unicellular heterotrophic eukaryote, Marophrys sp. SRT127 belonging to the Centroheliozoa, as the first report on this lineage. The circular-mapped mitochondrial genome, which is 113,062 bp in length, encodes 69 genes typically found in mitochondrial genomes. In addition, the Marophrys mitochondrial genome contains 19 group I introns. Of these, 11 introns have genes for homing endonuclease (HE) and phylogenetic analyses of HEs have shown that at least five Marophrys HEs are related to those in green algal plastid genomes, suggesting intron transfer between the Marophrys mitochondrion and green algal plastids. We also discovered a putative mitochondrial plasmid in linear form. Two genes encoded in the circular-mapped mitochondrial genome were found to share significant similarities to those in the linear plasmid, suggesting that the plasmid was integrated into the mitochondrial genome. These findings expand our knowledge on the diversity and evolution of the mobile genetic elements in mitochondrial genomes.

Annotation and analysis of the mitochondrial genome of Coniothyrium glycines, causal agent of red leaf blotch of soybean, reveals an abundance of homing endonucleases

Coniothyrium glycines, the causal agent of soybean red leaf blotch, is a USDA APHIS-listed Plant Pathogen Select Agent and potential threat to US agriculture. Sequencing of the C. glycines mt genome revealed a circular 98,533-bp molecule with a mean GC content of 29.01%. It contains twelve of the mitochondrial genes typically involved in oxidative phosphorylation (atp6, cob, cox1-3, nad1-6, and nad4L), one for a ribosomal protein (rps3), four for hypothetical proteins, one for each of the small and large subunit ribosomal RNAs (rns and rnl) and a set of 30 tRNAs. Genes were encoded on both DNA strands with cox1 and cox2 occurring as adjacent genes having no intergenic spacers. Likewise, nad2 and nad3 are adjacent with no intergenic spacers and nad5 is immediately followed by nad4L with an overlap of one base. Thirty-two introns, comprising 54.1% of the total mt genome, were identified within eight protein-coding genes and the rnl. Eighteen of the introns contained putative intronic ORFs with either LAGLIDADG or GIY-YIG homing endonuclease motifs, and an additional eleven introns showed evidence of truncated or degenerate endonuclease motifs. One intron possessed a degenerate N-acetyl-transferase domain. C. glycines shares some conservation of gene order with other members of the Pleosporales, most notably nad6-rnl-atp6 and associated conserved tRNA clusters. Phylogenetic analysis of the twelve shared protein coding genes agrees with commonly accepted fungal taxonomy. C. glycines represents the second largest mt genome from a member of the Pleosporales sequenced to date. This research provides the first genomic information on C. glycines, which may provide targets for rapid diagnostic assays and population studies.

Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta)

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Multiple origins of endosymbionts in Chlorellaceae with no reductive effects on the plastid or mitochondrial genomes

Ancient endosymbiotic relationships have led to extreme genomic reduction in many bacterial and eukaryotic algal endosymbionts. Endosymbionts in more recent and/or facultative relationships can also experience genomic reduction to a lesser extent, but little is known about the effects of the endosymbiotic transition on the organellar genomes of eukaryotes. To understand how the endosymbiotic lifestyle has affected the organellar genomes of photosynthetic green algae, we generated the complete plastid genome (plastome) and mitochondrial genome (mitogenome) sequences from three green algal endosymbionts (Chlorella heliozoae, Chlorella variabilis and Micractinium conductrix). The mitogenomes and plastomes of the three newly sequenced endosymbionts have a standard set of genes compared with free-living trebouxiophytes, providing no evidence for functional genomic reduction. Instead, their organellar genomes are generally larger and more intron rich. Intron content is highly variable among the members of Chlorella, suggesting very high rates of gain and/or loss of introns during evolution. Phylogenetic analysis of plastid and mitochondrial genes demonstrated that the three endosymbionts do not form a monophyletic group, indicating that the endosymbiotic lifestyle has evolved multiple times in Chlorellaceae. In addition, M. conductrix is deeply nested within the Chlorella clade, suggesting that taxonomic revision is needed for one or both genera.

Horizontal transfer of a ß-1,6-glucanase gene from an ancestral species of fungal endophyte to a cool-season grass host

Molecular characterisation has convincingly demonstrated some types of horizontal gene transfer in eukaryotes, but nuclear gene transfer between distantly related eukaryotic groups appears to have been rare. For angiosperms (flowering plants), nuclear gene transfer events identified to date have been confined to genes originating from prokaryotes or other plant species. In this report, evidence for ancient horizontal transfer of a fungal nuclear gene, encoding a ß-1,6-glucanase enzyme for fungal cell wall degradation, into an angiosperm lineage is presented for the first time. The gene was identified from de novo sequencing and assembly of the genome and transcriptome of perennial ryegrass, a cool-season grass species. Molecular analysis confirmed the presence of the complete gene in the genome of perennial ryegrass. No corresponding sequence was found in other plant species, apart from members of the Poeae sub-tribes Loliinae and Dactylidinae. Evidence suggests that a common ancestor of the two sub-tribes acquired the gene from a species ancestral to contemporary grass-associated fungal endophytes around 9-13 million years ago. This first report of horizontal transfer of a nuclear gene from a taxonomically distant eukaryote to modern flowering plants provides evidence for a novel adaptation mechanism in angiosperms.

The mitochondrial genome of the terrestrial carnivorous plant Utricularia reniformis (Lentibulariaceae): Structure, comparative analysis and evolutionary landmarks

The carnivorous plants of the family Lentibulariaceae have attained recent attention not only because of their interesting lifestyle, but also because of their dynamic nuclear genome size. Lentibulariaceae genomes span an order of magnitude and include species with the smallest genomes in angiosperms, making them a powerful system to study the mechanisms of genome expansion and contraction. However, little is known about mitochondrial DNA (mtDNA) sequences of this family, and the evolutionary forces that shape this organellar genome. Here we report the sequencing and assembly of the complete mtDNA from the endemic terrestrial Brazilian species Utricularia reniformis. The 857,234bp master circle mitochondrial genome encodes 70 transcriptionaly active genes (42 protein-coding, 25 tRNAs and 3 rRNAs), covering up to 7% of the mtDNA. A ltrA-like protein related to splicing and mobility and a LAGLIDADG homing endonuclease have been identified in intronic regions, suggesting particular mechanisms of genome maintenance. RNA-seq analysis identified properties with putative diverse and important roles in genome regulation and evolution: 1) 672kbp (78%) of the mtDNA is covered by full-length reads; 2) most of the 243kbp intergenic regions exhibit transcripts; and 3) at least 69 novel RNA editing sites in the protein-coding genes. Additional genomic features are hypothetical ORFs (48%), chloroplast insertions, including truncated plastid genes that have been lost from the chloroplast DNA (5%), repeats (5%), relics of transposable elements mostly related to LTR retrotransposons (5%), and truncated mitovirus sequences (0.4%). Phylogenetic analysis based on 32 different Lamiales mitochondrial genomes corroborate that Lentibulariaceae is a monophyletic group. In summary, the U. reniformis mtDNA represents the eighth largest plant mtDNA described to date, shedding light on the genomic trends and evolutionary characteristics and phylogenetic history of the family Lentibulariaceae.

Group II introns in wheat mitochondria have degenerate structural features and varied splicing pathways

Mitochondrial introns in flowering plant genes are virtually all classified as members of the group II ribozyme family although certain structural features have degenerated to varying degrees over evolutionary time. We are interested in the impact that unconventional intron architecture might have on splicing biochemistry in vivo and we have focused in particular on intronic domains V and VI, which for self-splicing introns provide a key component of the catalytic core and the bulged branchpoint adenosine, respectively. Notably, the two transesterification steps in classical group II splicing are the same as for nuclear spliceosomal introns and release the intron as a lariat. Using RT-PCR and circularized RT-PCR, we had previously demonstrated that several wheat mitochondrial introns which lack a branchpoint adenosine have atypical splicing pathways, and we have now extended this analysis to the full set of wheat introns, namely six trans-splicing and sixteen cis-splicing ones. A number of introns are excised using non-lariat pathways and interestingly, we find that several introns which do have a conventional domain VI also use pathways that appear to exploit other internal or external nucleophiles, with the lariat form being relatively minor. Somewhat surprisingly, several introns with weakly-structured domain V/VI helices still exhibit classical lariat splicing, suggesting that accessory factors aid in restoring a splicing-competent conformation. Our observations illustrate that the loss of conventional group II features during evolution is correlated with altered splicing biochemistry in an intron-distinctive manner.

Unparalleled replacement of native mitochondrial genes by foreign homologs in a holoparasitic plant

Horizontal gene transfer (HGT) among flowering plant mitochondria occurs frequently and, in most cases, leads to nonfunctional transgenes in the recipient genome. Parasitic plants are particularly prone to this phenomenon, but their mitochondrial genomes (mtDNA) have been largely unexplored. We undertook a large-scale mitochondrial genomic study of the holoparasitic plant Lophophytum mirabile (Balanophoraceae). Comprehensive phylogenetic analyses were performed to address the frequency, origin, and impact of HGT. The sequencing of the complete mtDNA of L. mirabile revealed the unprecedented acquisition of host-derived mitochondrial genes, representing 80% of the protein-coding gene content. All but two of these foreign genes replaced the native homologs and are probably functional in energy metabolism. The genome consists of 54 circular-mapping chromosomes, 25 of which carry no intact genes. The likely functional replacement of up to 26 genes in L. mirabile represents a stunning example of the potential effect of rampant HGT on plant mitochondria. The use of host-derived genes may have a positive effect on the host-parasite relationship, but could also be the result of nonadaptive forces.

Evolution of group I introns in Porifera: new evidence for intron mobility and implications for DNA barcoding

BACKGROUND

Mitochondrial introns intermit coding regions of genes and feature characteristic secondary structures and splicing mechanisms. In metazoans, mitochondrial introns have only been detected in sponges, cnidarians, placozoans and one annelid species. Within demosponges, group I and group II introns are present in six families. Based on different insertion sites within the cox1 gene and secondary structures, four types of group I and two types of group II introns are known, which can harbor up to three encoding homing endonuclease genes (HEG) of the LAGLIDADG family (group I) and/or reverse transcriptase (group II). However, only little is known about sponge intron mobility, transmission, and origin due to the lack of a comprehensive dataset. We analyzed the largest dataset on sponge mitochondrial group I introns to date: 95 specimens, from 11 different sponge genera which provided novel insights into the evolution of group I introns.

RESULTS

For the first time group I introns were detected in four genera of the sponge family Scleritodermidae (Scleritoderma, Microscleroderma, Aciculites, Setidium). We demonstrated that group I introns in sponges aggregate in the most conserved regions of cox1. We showed that co-occurrence of two introns in cox1 is unique among metazoans, but not uncommon in sponges. However, this combination always associates an active intron with a degenerating one. Earlier hypotheses of HGT were confirmed and for the first time VGT and secondary losses of introns conclusively demonstrated.

CONCLUSION

This study validates the subclass Spirophorina (Tetractinellida) as an intron hotspot in sponges. Our analyses confirm that most sponge group I introns probably originated from fungi. DNA barcoding is discussed and the application of alternative primers suggested.

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

Angiosperm mitochondrial genomes (mtDNA) exhibit variable quantities of alien sequences. Many of these sequences are acquired by intracellular gene transfer (IGT) from the plastid. In addition, frequent events of horizontal gene transfer (HGT) between mitochondria of different species also contribute to their expanded genomes. In contrast, alien sequences are rarely found in plastid genomes. Most of the plant-to-plant HGT events involve mitochondrion-to-mitochondrion transfers. Occasionally, foreign sequences in mtDNAs are plastid-derived (MTPT), raising questions about their origin, frequency, and mechanism of transfer. The rising number of complete mtDNAs allowed us to address these questions. We identified 15 new foreign MTPTs, increasing significantly the number of those previously reported. One out of five of the angiosperm species analyzed contained at least one foreign MTPT, suggesting a remarkable frequency of HGT among plants. By analyzing the flanking regions of the foreign MTPTs, we found strong evidence for mt-to-mt transfers in 65% of the cases. We hypothesize that plastid sequences were initially acquired by the native mtDNA via IGT and then transferred to a distantly-related plant via mitochondrial HGT, rather than directly from a foreign plastid to the mitochondrial genome. Finally, we describe three novel putative cases of mitochondrial-derived sequences among angiosperm plastomes.

Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns

Currently, complete mitochondrial genomes (mitogenomes) are available from all major land plant lineages except ferns. Sequencing of fern mitogenomes could shed light on the major evolutionary transitions that established mitogenomic diversity among extant lineages. In this study, we generated complete mitogenomes from the adder's tongue fern (Ophioglossum californicum) and the whisk fern (Psilotum nudum). The Psilotum mitogenome (628 kb) contains a rich complement of genes and introns, some of which are the largest of any green plant organellar genome. In the Ophioglossum mitogenome (372 kb), gene and intron content is slightly reduced, including the loss of all four mitochondrial ccm genes. Transcripts of nuclear Ccm genes also were not detected, suggesting loss of the entire mitochondrial cytochrome c maturation pathway from Ophioglossum. Both fern mitogenomes are highly repetitive, yet they show extremely low levels of active recombination. Transcriptomic sequencing uncovered ˜1000 sites of C-to-U RNA editing in both species, plus a small number (< 60) of U-to-C edit sites. Overall, the first mitochondrial genomes of ferns show a mix of features shared with lycophytes and/or seed plants and several novel genomic features, enabling a robust reconstruction of the mitogenome in the common ancestor of vascular plants.

Limited mitogenomic degradation in response to a parasitic lifestyle in Orobanchaceae

In parasitic plants, the reduction in plastid genome (plastome) size and content is driven predominantly by the loss of photosynthetic genes. The first completed mitochondrial genomes (mitogenomes) from parasitic mistletoes also exhibit significant degradation, but the generality of this observation for other parasitic plants is unclear. We sequenced the complete mitogenome and plastome of the hemiparasite Castilleja paramensis (Orobanchaceae) and compared them with additional holoparasitic, hemiparasitic and nonparasitic species from Orobanchaceae. Comparative mitogenomic analysis revealed minimal gene loss among the seven Orobanchaceae species, indicating the retention of typical mitochondrial function among Orobanchaceae species. Phylogenetic analysis demonstrated that the mobile cox1 intron was acquired vertically from a nonparasitic ancestor, arguing against a role for Orobanchaceae parasites in the horizontal acquisition or distribution of this intron. The C. paramensis plastome has retained nearly all genes except for the recent pseudogenization of four subunits of the NAD(P)H dehydrogenase complex, indicating a very early stage of plastome degradation. These results lend support to the notion that loss of ndh gene function is the first step of plastome degradation in the transition to a parasitic lifestyle.

The Reverse Transcriptase/RNA Maturase Protein MatR Is Required for the Splicing of Various Group II Introns in Brassicaceae Mitochondria

Group II introns are large catalytic RNAs that are ancestrally related to nuclear spliceosomal introns. Sequences corresponding to group II RNAs are found in many prokaryotes and are particularly prevalent within plants organellar genomes. Proteins encoded within the introns themselves (maturases) facilitate the splicing of their own host pre-RNAs. Mitochondrial introns in plants have diverged considerably in sequence and have lost their maturases. In angiosperms, only a single maturase has been retained in the mitochondrial DNA: the matR gene found within NADH dehydrogenase 1 (nad1) intron 4. Its conservation across land plants and RNA editing events, which restore conserved amino acids, indicates that matR encodes a functional protein. However, the biological role of MatR remains unclear. Here, we performed an in vivo investigation of the roles of MatR in Brassicaceae. Directed knockdown of matR expression via synthetically designed ribozymes altered the processing of various introns, including nad1 i4. Pull-down experiments further indicated that MatR is associated with nad1 i4 and several other intron-containing pre-mRNAs. MatR may thus represent an intermediate link in the gradual evolutionary transition from the intron-specific maturases in bacteria into their versatile spliceosomal descendants in the nucleus. The similarity between maturases and the core spliceosomal Prp8 protein further supports this intriguing theory.

Loss of a Trans-Splicing nad1 Intron from Geraniaceae and Transfer of the Maturase Gene matR to the Nucleus in Pelargonium

The mitochondrial nad1 gene of seed plants has a complex structure, including four introns in cis or trans configurations and a maturase gene (matR) hosted within the final intron. In the geranium family (Geraniaceae), however, sequencing of representative species revealed that three of the four introns, including one in a trans configuration and another that hosts matR, were lost from the nad1 gene in their common ancestor. Despite the loss of the host intron, matR has been retained as a freestanding gene in most genera of the family, indicating that this maturase has additional functions beyond the splicing of its host intron. In the common ancestor of Pelargonium, matR was transferred to the nuclear genome, where it was split into two unlinked genes that encode either its reverse transcriptase or maturase domain. Both nuclear genes are transcribed and contain predicted mitochondrial targeting signals, suggesting that they express functional proteins that are imported into mitochondria. The nuclear localization and split domain structure of matR in the Pelargonium nuclear genome offers a unique opportunity to assess the function of these two domains using transgenic approaches.