Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes

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 first complete mitochondrial genome of the invasive water hyacinth, Eichhornia crassipes

Eichhornia crassipes is an aquatic plant in tropical and subtropical regions, renowned for its notorious invasive tendencies. In this study, we assembled the complete mitogenome of E. crassipes into a single circle molecule of 397,361 bp. The mitogenome has 58 unique genes, including 37 protein-coding genes (PCGs), 18 tRNA genes, three rRNA genes, and 47 % GC content. Sixteen (6.93 %) homologous fragments, ranging from 31 bp to 8548 bp, were identified, indicating the transfer of genetic material from chloroplasts to mitochondria. In addition, we detected positive selection in six PCGs (ccmB, ccmC, ccmFC, nad3, nad4 and sdh4), along with the identification of 782 RNA editing sites across 37 mt-PCGs. These findings suggest a potential contribution to the robust adaptation of this invasive plant to the stressful environment. Lastly, we inferred that phylogenetic conflicts of E. crassipes between the plastome and mitogenome may be attributed to the difference in nucleotide substitution rates between the two organelle genomes. In conclusion, our study provided vital genomic resources for further understanding the invasive mechanism of this species and exploring the dynamic evolution of mitogenomes within the monocot clade.

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.

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

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

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

Comparative analysis of mitochondrial genomes reveals marine adaptation in seagrasses

BACKGROUND

Seagrasses are higher marine flowering plants that evolved from terrestrial plants, but returned to the sea during the early evolution of monocotyledons through several separate lineages. Thus, they become a good model for studying the adaptation of plants to the marine environment. Sequencing of the mitochondrial (mt) genome of seagrasses is essential for understanding their evolutionary characteristics.

RESULTS

In this study, we sequenced the mt genome of two endangered seagrasses (Zostera japonica and Phyllospadix iwatensis). These data and data on previously sequenced mt genomes from monocotyledons provide new evolutionary evidence of genome size reduction, gene loss, and adaptive evolution in seagrasses. The mt genomes of Z. japonica and P. iwatensis are circular. The sizes of the three seagrasses (including Zostera marine) that have been sequenced to date are smaller than that of other monocotyledons. Additionally, we found a large number of repeat sequences in seagrasses. The most abundant long repeat sequences were 31-40 bp repeats. Our study also found that seagrass species lost extensive ribosomal protein genes during evolution. The rps7 gene and the rpl16 gene of P. iwatensis are exceptions to this trend. The phylogenetic analysis based on the mt genome strongly supports the previous results. Furthermore, we identified five positive selection genes (atp8, nad3, nad6, ccmFn, and matR) in seagrasses that may be associated with their adaptation to the marine environment.

CONCLUSIONS

In this study, we sequenced and annotated the mt genomes of Z. japonica and P. iwatensis and compared them with the genome of other monocotyledons. The results of this study will enhance our understanding of seagrass adaptation to the marine environment and can inform further investigations of the seagrass mt genome.

Mitochondrial genomic data are effective at placing mycoheterotrophic lineages in plant phylogeny

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing

Salix L. (willows) is one of the most taxonomically complex genera of flowering plants, including shrubs, tall trees, bushes, and prostrate plants. Despite the high species diversity, only five mitochondrial genomes (mitogenomes) have been released in this genus. Salix wilsonii is an important ornamental and economic willow tree in section Wilsonia of the genus Salix. In this study, the S. wilsonii mitogenome was assembled into a typical circular structure with a size of 711,456 bp using PacBio HiFi sequencing. A total of 58 genes were annotated in the S. wilsonii mitogenome, including 33 protein-coding genes (PCGs), 22 tRNAs, and 3 rRNAs. In the S. wilsonii mitogenome, four genes (mttB, nad3, nad4, and sdh4) were found to play important roles in its evolution through selection pressure analysis. Collinearity analysis of six Salix mitogenomes revealed high structural variability. To determine the evolutionary position of S. wilsonii, we conducted a phylogenetic analysis of the mitogenomes of S. wilsonii and 12 other species in the order Malpighiales. Results strongly supported the segregation of S. wilsonii and other five Salix species with 100% bootstrap support. The comparative analysis of the S. wilsonii mitogenome not only sheds light on the functional and structural features of S. wilsonii but also provides essential information for genetic studies of the genus Salix.

A High-Quality Genome Assembly of the Mitochondrial Genome of the Oil-Tea Tree Camellia gigantocarpa (Theaceae)

Camellia gigantocarpa is one of the oil-tea trees whose seeds can be used to extract high-quality vegetable oil. To date, there are no data on the mitochondrial genome of the oil-tea tree, in contrast to the tea-tree C. sinensis, which belongs to the same genus. In this paper, we present the first complete mitochondrial genomes of C. gigantocarpa obtained using PacBio Hi-Fi (high-fidelity) and Hi-C sequencing technologies to anchor the 970,410 bp genome assembly into a single sequence. A set of 44 protein-coding genes, 22 non-coding genes, 746 simple sequence repeats (SSRs), and more than 201 kb of repetitive sequences were annotated in the genome assembly. The high percentage of repetitive sequences in the mitochondrial genome of C. gigantocarpa (20.81%) and C.sinensis (22.15%, tea tree) compared to Arabidopsis thaliana (4.96%) significantly increased the mitogenome size in the genus Camellia. The comparison of the mitochondrial genomes between C. gigantocarpa and C. sinensis revealed genes exhibit high variance in gene order and low substitution rate within the genus Camellia. Information on the mitochondrial genome provides a better understanding of the structure and evolution of the genome in Camellia and may contribute to further study of the after-ripening process of oil-tea trees.

Mitochondrial genome of Garcinia mangostana L. variety Mesta

Fruits of Garcinia mangostana L. (mangosteen) are rich in nutrients with xanthones found in the pericarp having great pharmaceutical potential. Mangosteen variety Mesta is only found in Malaysia, which tastes sweeter than the common Manggis variety in Southeast Asia. In this study, we report the complete mitogenome of G. mangostana L. variety Mesta with a total sequence length of 371,235 bp of which 1.7% could be of plastid origin. The overall GC content of the mitogenome is 43.8%, comprising 29 protein-coding genes, 3 rRNA genes, and 21 tRNA genes. Repeat and tandem repeat sequences accounted for 5.8% and 0.15% of the Mesta mitogenome, respectively. There are 333 predicted RNA-editing sites in Mesta mitogenome. These include the RNA-editing events that generated the start codon of nad1 gene and the stop codon of ccmFC gene. Phylogenomic analysis using both maximum likelihood and Bayesian analysis methods showed that the mitogenome of mangosteen variety Mesta was grouped under Malpighiales order. This is the first complete mitogenome from the Garcinia genus for future evolutionary studies.

Pollen in water of unstable salinity: Evolution and function of dynamic apertures in monocot aquatics

PREMISE

The sporoderm of seed-plant pollen grains typically has apertures in which the outer sporopollenin-bearing layer is relatively sparse. The apertures allow regulation of the internal volume of the pollen grain during desiccation and rehydration (harmomegathy) and also serve as sites of pollen germination. A small fraction of angiosperms undergo pollination in water or at the water surface, where desiccation is unlikely. Their pollen grains commonly lack apertures, though with some notable exceptions. We tested a hypothesis that in some angiosperm aquatics that inhabit water of unstable salinity, the pollen apertures accommodate osmotic effects that occur during pollination in such conditions.

METHODS

Pollen grains of the tepaloid clade of the monocot order Alismatales, which contains ecologically diverse aquatic and marshy plants, were examined using light microscopy and scanning electron microscopy. We used Ruppia as a model to test pollen grain response in water of various salinities. Pollen aperture evolution was also analyzed using molecular tree topologies.

RESULTS

Phylogenetic optimizations demonstrated an evolutionary loss and two subsequent regains of the aperturate condition in the tepaloid clade of Alismatales. Both of the taxa that have reverted to aperturate pollen (Ruppia, Ruppiaceae; Althenia, Potamogetonaceae) are adapted to changeable water salinity. Direct experiments with Ruppia showed that the pollen apertures have a role in a harmomegathic response to differences in water salinity.

CONCLUSIONS

Our results showed that the inferred regain of pollen apertures represents an adaptation to changeable water salinity. We invoke a loss-and-regain scenario, prompting questions that are testable using developmental genetics and plant physiology.

Lemnaceae and Orontiaceae Are Phylogenetically and Morphologically Distinct from Araceae

Duckweeds comprise a distinctive clade of pleustophytic monocots that traditionally has been classified as the family Lemnaceae. However, molecular evidence has called into question their phylogenetic independence, with some authors asserting instead that duckweeds should be reclassified as subfamily Lemnoideae of an expanded family Araceae. Although a close phylogenetic relationship of duckweeds with traditional Araceae has been supported by multiple studies, the taxonomic disposition of duckweeds must be evaluated more critically to promote nomenclatural stability and utility. Subsuming duckweeds as a morphologically incongruent lineage of Araceae effectively eliminates the family category of Lemnaceae that has been widely used for many years. Instead, we suggest that Araceae subfamily Orontioideae should be restored to family status as Orontiaceae, which thereby would enable the recognition of three morphologically and phylogenetically distinct lineages: Araceae, Lemnaceae, and Orontiaceae.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Both Conifer II and Gnetales are characterized by a high frequency of ancient mitochondrial gene transfer to the nuclear genome

BACKGROUND

Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes.

RESULTS

Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants.

CONCLUSIONS

This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants.

Comparative Analysis of Mitochondrial Genome Features among Four Clonostachys Species and Insight into Their Systematic Positions in the Order Hypocreales

The mycoparasite fungi of Clonostachys have contributed to the biological control of plant fungal disease and nematodes. The Clonostachys fungi strains were isolated from Ophiocordyceps highlandensis, Ophiocordycepsnigrolla and soil, which identified as Clonostachyscompactiuscula, Clonostachysrogersoniana, Clonostachyssolani and Clonostachys sp. To explore the evolutionary relationship between the mentioned species, the mitochondrial genomes of four Clonostachys species were sequenced and assembled. The four mitogenomes consisted of complete circular DNA molecules, with the total sizes ranging from 27,410 bp to 42,075 bp. The GC contents, GC skews and AT skews of the mitogenomes varied considerably. Mitogenomic synteny analysis indicated that these mitogenomes underwent gene rearrangements. Among the 15 protein-coding genes within the mitogenomes, the nad4L gene exhibited the least genetic distance, demonstrating a high degree of conservation. The selection pressure analysis of these 15 PCGs were all below 1, indicating that PCGs were subject to purifying selection. Based on protein-coding gene calculation of the significantly supported topologies, the four Clonostachys species were divided into a group in the phylogenetic tree. The results supplemented the database of mitogenomes in Hypocreales order, which might be a useful research tool to conduct a phylogenetic analysis of Clonostachys. Additionally, the suitable molecular marker was significant to study phylogenetic relationships in the Bionectriaceae family.

Complete mitochondrial genome of Iniistius trivittatus and unique variation in two observed inserts between rRNA and tRNA genes in wrasses

Background

The family Labridae made up of 519 species in the world. The functional evolution of the feeding-related jaws leaded to differentiation of species, and the pharyngeal jaw apparatus evolved independently, but evolutionary mechanism still remain unaddressed in wrasses. Mitogenomes data can be used to infer genetic diversification and investigate evolutionary history of wrasses, whereas only eight complete mitogenomes in this family have been sequenced to date. Here, we sequenced the complete mitogenomes of Iniistius trivittatus to investigate genetic differentiation among wrasse species.

Results

We sequenced the complete mitogenomes of I. trivittatus using a novel PCR strategy. The I. trivittatus mitogenomes is 16,820 bp in length and includes 13 protein -coding genes, 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region. Compared to eight known mitochondrial genome, 2 additional noncoding regions (lengths of 121 and 107 bp), or so-called inserts, are found in the intergenic regions 12S rRNA - tRNA^Val - 16S rRNA. The presumed origin of the two rare inserts is from tRNA- related retrotransposons. Compared with cytochrome b gene, the two insert sequences are highly conserved at the intraspecies level, but they showed significant variation and low similarity (< 70%) at the interspecies level. The insert events were only observed in I. trivittatus by checking the phylogenetic trees based on the complete mitogenomes of Labrida species. This finding provides evidence that in the mitogenomes, retrotransposon inserts result in intraspecific homoplasmy and interspecific heteroplasmy by natural selection and adaptation to various environments.

Conclusions

This study found additional mitogenome inserts limited in wrasse species. The rRNA genes with inserts might have experienced a selective pressure for adaptation to feeding modes. Such knowledge can enable a better understanding of molecular mechanism underlying morphological evolution in wrasses.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

Mycoheterotrophic Epirixanthes (Polygalaceae) has a typical angiosperm mitogenome but unorthodox plastid genomes

BACKGROUND AND AIMS

Fully mycoheterotrophic plants derive carbon and other nutrients from root-associated fungi and have lost the ability to photosynthesize. While mycoheterotroph plastomes are often degraded compared with green plants, the effect of this unusual symbiosis on mitochondrial genome evolution is unknown. By providing the first complete organelle genome data from Polygalaceae, one of only three eudicot families that developed mycoheterotrophy, we explore how both organellar genomes evolved after loss of photosynthesis.

METHODS

We sequenced and assembled four complete plastid genomes and a mitochondrial genome from species of Polygalaceae, focusing on non-photosynthetic Epirixanthes. We compared these genomes with those of other mycoheterotroph and parasitic plant lineages, and assessed whether organelle genes in Epirixanthes experienced relaxed or intensified selection compared with autotrophic relatives.

KEY RESULTS

Plastomes of two species of Epirixanthes have become substantially degraded compared with that of autotrophic Polygala. Although the lack of photosynthesis is presumably homologous in the genus, the surveyed Epirixanthes species have marked differences in terms of plastome size, structural rearrangements, gene content and substitution rates. Remarkably, both apparently replaced a canonical plastid inverted repeat with large directly repeated sequences. The mitogenome of E. elongata incorporated a considerable number of fossilized plastid genes, by intracellular transfer from an ancestor with a less degraded plastome. Both plastid and mitochondrial genes in E. elongata have increased substitution rates, but the plastid genes of E. pallida do not. Despite this, both species have similar selection patterns operating on plastid housekeeping genes.

CONCLUSIONS

Plastome evolution largely fits with patterns of gene degradation seen in other heterotrophic plants, but includes highly unusual directly duplicated regions. The causes of rate elevation in the sequenced Epirixanthes mitogenome and of rate differences in plastomes of related mycoheterotrophic species are not currently understood.

The alternative reality of plant mitochondrial DNA: One ring does not rule them all

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms.

The complete mitochondrial genome of Paravannella minima (Amoebozoa, Discosea, Vannellida)

We present a complete sequence and describe the organization of the mitochondrial genome of the amoeba Paravannella minima (Amoebooza, Discosea, Vannellida). This tiny species represents a branch at the base of Vannellida tree, to the moment being its earliest-branching lineage. The circular mitochondrial DNA of this species has 53,464 bp in length and contains 30 protein-coding genes, 2 ribosomal RNAs, 23 transfer RNAs, and 15 open reading frames. This genome is significantly longer and contains more protein-coding genes than any yet sequenced mitochondrial genome of vannellid amoebae. Unlike the previously sequenced mitochondrial genomes of Vannellida, which should be translated using the "Table 4" (the mold, protozoan, and coelenterate mitochondrial code), that of P. minima can be properly translated using the universal genetic code.

Assembly of the Mitochondrial Genome in the Campanulaceae Family Using Illumina Low-Coverage Sequencing

Platycodongrandiflorus (balloon flower) and Codonopsislanceolata (bonnet bellflower) are important herbs used in Asian traditional medicine, and both belong to the botanical family Campanulaceae. In this study, we designed and implemented a de novo DNA sequencing and assembly strategy to map the complete mitochondrial genomes of the first two members of the Campanulaceae using low-coverage Illumina DNA sequencing data. We produced a total of 28.9 Gb of paired-end sequencing data from the genomic DNA of P.grandiflorus (20.9 Gb) and C.lanceolata (8.0 Gb). The assembled mitochondrial genome of P.grandiflorus was found to consist of two circular chromosomes; the master circle contains 56 genes, and the minor circle contains 42 genes. The C.lanceolata mitochondrial genome consists of a single circle harboring 54 genes. Using a comparative genome structure and a pattern of repeated sequences, we show that the P.grandiflorus minor circle resulted from a recombination event involving the direct repeats of the master circle. Our dataset will be useful for comparative genomics and for evolutionary studies, and will facilitate further biological and phylogenetic characterization of species in the Campanulaceae.

Practical considerations for plant phylogenomics

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome-scale data from diverse lineages. The development of high-throughput sequencing and long-read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project-dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up-and-coming technologies that may help propel the field even further.