- 1National Engineering Laboratory for Tree Breeding, Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, The Tree and Ornamental Plant Breeding and Biotechnology Laboratory of National Forestry and Grassland Administration, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 2Department of Ecology and Environmental Science, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
- 3Department of Bioinformatics, Ori (Shandong) Gene Science and Technology Co., Ltd., Weifang, China
Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion (r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time (r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8–12 protein-coding genes, of which 4–7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes.
1 Introduction
In contrast to the ∼8,600 chloroplast (https://www.ncbi.nlm.nih.gov/) and ∼800 nuclear reference genomes assembled (Marks et al., 2021) for land plant taxa thus far, only ∼350 mitogenomes have been released (https://www.ncbi.nlm.nih.gov/; by November 2021). These published plant mitogenomes revealed extensive diversity in genome size, genome rearrangement, nucleotide substitution rate, and gene contents (Palmer and Herbon, 1988; Chaw et al., 2008; Smith and Keeling, 2015; Mower, 2020; Wu et al., 2020). In angiosperms, mitogenomes span a 170-fold range in size, from ∼66 Kb in Viscum scurruloideum to ∼11.3 Mb in Silene conica (Sloan et al., 2012; Skippington et al., 2015). The mechanisms contributing to the size variation are hypothesized to include horizontal transfer, intergenomic transfer, and repetitive sequence proliferation (Goremykin et al., 2012; Rice et al., 2013; Guo et al., 2017). For example, the mitogenome of the basal angiosperm Amborella trichopoda acquired ∼20% of its sequences through horizontal transfer from green algae, mosses, and other plants (Rice et al., 2013), while apple and maize acquired 20–78% of mitochondrial DNA (mtDNA) from their nuclei (Goremykin et al., 2012). The repeat content in 82 angiosperm mitogenomes ranged from 1.3% to 48.9% (Dong et al., 2018). A phylogenetic pattern in these traits is, however, not yet discernable.
Plant mitogenomes are often characterized by a strong dichotomy between rates of sequence and structural evolution (Palmer and Herbon, 1988). Analyses of the available mitogenomes show that while most plant lineages have low substitution rates, it is misleading to categorize all plants in this way; synonymous substitution rates are exceptionally variable, with four orders of magnitude across seed plants (Mower et al., 2007). The reason for such fluid substitution rates remains unclear.
Gymnosperms contain ∼1,000 species representing five of the six main lineages of seed plants, including Cupressophytes, Pinaceae, Gnetophytes, Ginkgo, and Cycads (Wang and Ran, 2014). Up to now, only ∼20 gymnosperm mitogenomes have been released, which are too few to allow a proper understanding of the properties and diversities of this genome. Despite the limited data, interesting and sometimes unexpected patterns of mtDNA variation have emerged. Similar to angiosperms, gymnosperm mitogenomes also show a wide range in size, with Ginkgo biloba (∼347 Kb) and Larix sibirica (∼11.7 Mb) being the extremes (Guo et al., 2016; Putintseva et al., 2020). The driving factors for this size heterogeneity are, however, debated because estimations of repetitive contents and nucleus-derived sequences vary among studies (Guo et al., 2016; Kan et al., 2020; Sullivan et al., 2020).
Fast structural rearrangement of gymnosperm mitogenomes is reported to result in a substantial loss of synteny and shared sequences between closely related species with divergence times as small as 15 million years ago (MYA) (Guo et al., 2016; Cole et al., 2018; Sullivan et al., 2020). On the other hand, conserved structures and high degrees of synteny are found between Cycas and Ginkgo that diversified >300 MYA (Guo et al., 2016). Parallel to this unpredictable structural diversity are the highly variable substitution rates among the taxa (Guo et al., 2016; Sullivan et al., 2020; Kan et al., 2021). In addition, the number of protein-coding and other genes vary between gymnosperm lineages (Guo et al., 2020; Kan et al., 2021). These observations show a largely lineage-specific mode of mitogenome evolution in gymnosperms, but the mechanisms and phylogenetic patterns for this diversity are unclear due to sparse sampling and limited analyses.
In this study, we report a complete mitogenome assembly of Platycladus orientalis, Cupressaceae (Cupressophytes). Cupressaceae is the largest conifer family with 25 genera; Platycladus is a monotypic genus in this family. This new genome thus fills a gap in the still sparse coverage for this clade. We then sampled 17 gymnosperm and two angiosperm mitogenomes, representing all six major lineages of seed plants, to characterize gene content, transposable elements (TEs), dispersed repeats, tandem repeats, sequence turnover, structural rearrangement and substitution rate variation. We provide quantitative estimates of the role of repetitive elements in genome expansion, the mode of repeats accumulation, and the rate of synteny decay among taxa. Our results reinforce some early observations but also shed new light on the dynamics of mitogenome evolution.
2 Materials and Methods
2.1 Mitochondrial Genome Assembly and Annotation
We collected fresh young leaves of P. orientalis from an elite tree in a seed orchard in Jiaxian, Henan Province, China. DNA used for Illumina sequencing was isolated using a cetyl trimethyl ammonium bromide (CTAB)-based method (Doyle and Doyle, 1987). For the Oxford Nanopore Technologies (ONT) sequencing, we prepared DNA using a DNeasy Plant Mini Kit (QIAGEN). A short-insert 2 × 150 bp pair-end library was constructed following the manufacturer’s PCR-free protocol (Illumina) and sequenced on an Illumina NovaSeq 6000 platform. Library for ONT sequencing was prepared following the Nanopore 1D Genomic DNA by ligation (SQK-LSK109)–PromethION protocol, and sequenced on a PromethION platform. The short reads from Illumina were filtered using fastp v0.21.0 (Chen et al., 2018) with parameters of “-q 20 -l 36 --cut_right --n_base_limit 0 -c.” The long reads from Nanopore were corrected using NextDenovo v1.1.0 (https://github.com/Nextomics/NextDenovo). We generated ∼110 Gb Illumina clean short reads and ∼808 Gb (40.3 M reads with a N50 of 25 Kb) corrected ONT long reads.
We applied a long and short reads hybrid strategy to assemble the P. orientalis mitogenome. First, the clean short reads were de novo assembled using GetOrganelle v1.6.4 (Jin et al., 2020) with parameters of “-R 50 -k 67,87,107,127 -F embplant_mt --reduce-reads-for-coverage inf -w 127 -t 12 --max-reads 7.5E8”. The resulting unitig graph was manually edited using Bandage v0.8.1 (Wick et al., 2015) to remove chloroplast- and nuclear-derived unitig nodes. The corrected long reads were then mapped to the graph using the minimap2 v2.24-r1122 (Li, 2018). Finally, repeats on the graph were resolved by aligning with the mapped long reads.
The protein-coding genes, non-coding rRNA, and tRNA genes were annotated using the Organelle Genome Annotation Pipeline (OGAP, https://github.com/zhangrengang/OGAP) with parameters of “-taxon Acrogymnospermae Bryophyta Liliopsida commelinids fabids malvids -mt -trn_struct”. More specifically, protein-coding genes were identified with Exonerate v2.2.0 (Slater and Birney, 2005) and AUGUSTUS v3.3.3 (Stanke et al., 2008), tRNA genes with tRNAscan-SE v2.0.5 (Lowe and Eddy, 1997), and rRNA genes with BLAT v36x2 (Kent, 2002) within the OGAP.
RNA editing sites in protein-coding genes of each species were predicted using the PREP-Mt (Mower, 2005) with a cutoff of 0.2. To validate the accuracy of the predictions, we searched for RNAseq data of P. orientalis (Supplementary Table S1) and mapped them to the protein-coding genes with bwa mem v0.7.17-r1188 (Li, 2013) to determine the edited sites using REDItoolDenovo v1.3 (Lo Giudice et al., 2020) with minimum coverage of 10 and mapping quality score of 30.
2.2 Repetitive Sequences
To estimate the repetitive sequences in mitogenomes, we annotated tandem repeats, TEs, and dispersed repeats in P. orientalis and 19 other taxa representing the six major lineages of seed plants included in this study (Table 1). Tandem repeats were identified using Tandem Repeats Finder v4.10.0 with the parameters of “2 7 7 80 10 50 800 -d” (Benson, 1999). For TEs, we first used the de novo TE family identification and modeling package RepeatModeler v1.0.10 (Smit and Hubley, 2008) to train a TE database from each mitogenome. The RepeatModeler derived library was then combined with the database “RepeatMasker.lib” of RepeatMasker v4.07 (Smit et al., 2013) and the repeat library of Picea abies v1.0 assembly (Nystedt et al., 2013) to generate a joint repeats library. Finally, TEs in each genome were identified by RepeatMasker using this reference database.
Dispersed repeats were identified with a BLASTN search using a word size of seven and an e-value of 1E-6 as in a previous study (Guo et al., 2016). Overlapping regions of dispersed repeats were counted only once when calculating their total size. The dispersed repeats were grouped into small-, medium-, and large-size classes with sequence lengths of <100 bp, 100–1,000 bp, and ≥1,000 bp, respectively. Divergence among repeats was estimated using Kimura two-parameter distance (K2P) (Kimura, 1980) for pairwise repeat sequences.
2.3 Identify Mitochondrial Genes From Transcriptomes
Based on gene annotation results, we found 8–12 genes were absent in mitogenomes of Cupressophytes and Gnetophytes (see Section 3.3). To investigate the possibility of their transfer to nuclei, we first tried to recover these genes from transcriptomes. We searched for RNAseq data of Cupressophytes and Gnetophytes in NCBI (https://www.ncbi.nlm.nih.gov/), and de novo assembled transcriptomes for Cupressus sempervirens, Hesperocyparis glabra, Taxus cuspidata, Podocarpus macrophyllus, Gnetum gnemon, and Welwitschia mirabilis (Supplementary Table S2). The Illumina reads were filtered using fastp v0.21.0 (Chen et al., 2018) with the parameters of “-q 20 -l 36 --cut_right --n_base_limit 0 -c,” generating 6.3–14.7 Gb clean reads per species. We then de novo assembled the transcriptomes using Trinity v2.8.5 with default parameters (Grabherr et al., 2011). The transcriptome of P. orientalis is from our early study (Hu et al., 2016). For Araucaria heterophylla, we only found 212 ESTs/cDNAs in TreeGenes database (Wegrzyn et al., 2008) (Supplementary Table S2). We searched for intact mitochondrial genes in these transcriptomes and ESTs/cDNAs using OGAP with the parameters of “-taxon Acrogymnospermae Bryophyta Liliopsida commelinids fabids malvids -mt -trn_struct -trans”. Potential chloroplast genes were filtered out after BLAST against NCBI’s NT database (E-value < 1E-40, sequence coverage >99%, and sequence identity >99%).
Second, to confirm that the transcriptome-recovered mitochondrial genes (mitogenes) were in nuclear genomes, we used TBLASTN to search for matches of these genes in reference genomes of Cupressophytes and Gnetophytes, including Sequoiadendron giganteum (Scott et al., 2020), Sequoia sempervirens (Neale et al., 2021), Taxus chinensis (Xiong et al., 2021), Taxus wallichiana (Cheng et al., 2021), Gnetum montanum (Wan et al., 2018), and Welwitschia mirabilis (Wan et al., 2021) (Supplementary Table S3).
Third, we used mapping depth of reads of resequencing data to verify the genome origins of the transcriptome-recovered mitogenes. In a plant cell, the abundance of nuclear and organelle genomes can differ by thousands fold with the copy numbers of plastid the highest, mitochondrial intermediate and nuclear the lowest (Petit and Vendramin, 2007). Based on this expectation, we investigated the depth of read coverage of each recovered genes. The mitochondrial genes retrieved from mitogenome and transcriptomes were used as reference. We also include the low-copy nuclear gene LEAFY and chloroplast gene matK as spike-in reference. We downloaded whole-genome resequencing data of Cupressus sempervirens, Hesperocyparis glabra, Taxus cuspidata, Podocarpus macrophyllus, Gnetum gnemon, and Welwitschia mirabilis from Sequence Read Archive (SRA) of NCBI (https://www.ncbi.nlm.nih.gov/sra) (Supplementary Table S4), the resequencing data of P. orientalis is from this study. The Illumina reads were filtered using fastp v0.21.0 (Chen et al., 2018) with the same parameters as in Section 2.1, generating 7.2–31.4 Gb clean reads per species. The clean reads were mapped to the reference genes using bwa mem v0.7.17-r1188 with default parameters. The resulting sequence alignment maps (SAMs) were sorted with SAMtools v1.14 (https://github.com/samtools/samtools), reads with alignment score (AS) <25 were removed by using ngsutilsj (https://github.com/compgen-io/ngsutilsj). We used SAMtools to compute the depth of the aligned reads at each position of the gene.
2.4 Genome Alignments and Shared DNA Blocks
To identify synteny blocks among the mitogenomes of the major gymnosperm lineages, we performed pairwise genome alignment using Mauve v2.4.0 with default parameters (Darling et al., 2010). To calculate the degree of shared DNA between species, each pair of mitogenomes were aligned using BLASTN with a word size of seven and an E-value of 1E-6.
2.5 Phylogenetic Analysis and Estimation of Nucleotide Substitution Rates
We used 28 protein-coding genes shared among all 20 mitogenomes included in this study to reconstruct a phylogenetic tree. Edited protein sequences were individually aligned with MAFFT v7.453 with default parameters and converted into the corresponding edited codon sequence alignments. Poorly aligned regions were trimmed using BMGE v1.1 (Criscuolo and Gribaldo, 2010) with parameters of “-t CODON -g 0.5.” The trimmed alignments of each gene were concatenated into a single alignment of 26,844 bp. This alignment was used to build a maximum-likelihood (ML) phylogenetic tree using IQ-TREE v1.6.12 (Nguyen et al., 2015), using the best-fit model TVM + F + R3 selected by ModelFinder (Kalyaanamoorthy et al., 2017) and with 1,000 replications of ultrafast bootstrap and Shimodaira-Hasegawa-like approximate likelihood-ratio (SH-aLRT) test. The branch lengths in units of synonymous (dS) and nonsynonymous (dN) substitution rates were estimated under the free-ratio branch model using codeml in PAML v4.9j package (Yang, 2007).
To estimate the divergence time between lineages, we run the MCMCTree in PAML v4.9j (Yang, 2007) for 2,000,000 iterations and 100,000 iterations for a burn-in. To calibrate the tree, we used the following time points of diversification from the TimeTree web (http://www.timetree.org/): 289–337 MYA between gymnosperms and angiosperms, 282–324 MYA between Cycads and other gymnosperm clades, 159–246 MYA between Cycas taitungensis and the other Cycads species, 134–197 MYA between Abies sibirica and other species within Pinaceae, 7.6 to 30.7 MYA between Picea abies and Picea glauca, 59–87 MYA between Pinus taeda and Pinus lambertiana, 99–123 MYA between Gnetum gnemon and Welwitschia mirabilis, and 213–272 MYA between the clades consisting of Podocarpus macrophyllus and Araucaria heterophylla and the other species of Cupressophytes.
3 Results
3.1 The Mitogenome of P. orientalis
We assembled the mitogenome of P. orientalis into two circular-mapping chromosomes representing a master ring and a sub-genomic ring with a total size of ∼2.6 Mb (Figure 1A, Table 2). This genome lacks large structural variants but is rich in dispersed repeats (Figures 1B,C). The genome size is similar to C. sempervirens (2.7 Mb) and H. glabra (2.0 Mb) in the family of Cupressaceae, but more than 5-fold larger than T. cuspidata and ∼2–4 fold smaller than Picea and Larix in Pinaceae (Table 1). Across the Cupressophytes and Pinaceae, the mitogenomes vary 6–10 fold in size, in contrast to the relatively stable and small genomes in Cycads and Gnetophytes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-g001.jpg)
FIGURE 1. The mitogenome assembly of Platycladus orientalis. (A) The order, orientation, and size of genes on the two chromosomes of P. orientalis mitogenome. The nad1, nad2, nad4, nad5, and cox2 are trans-splicing genes. Each box is proportional to the size of the gene including introns; genes with <200 nucleotides are shown as 200 nucleotides. Intergenic regions are not to scale. (B) The dot plot shows self-alignment. Dots represent the homology sequences with length >100 bp. (C) Distribution of dispersed repeats. Large-size repeats ≥1,000 bp in length are indicated in red, medium-size repeats in the range of 100–1,000 bp are in orange, and small-size repeats <100 bp are in blue. The numbers on the ring anchor genome coordinates in kilobases.
The mitogenome of P. orientalis contains 32 distinct protein-coding genes (Figure 1A, Table 2), seven tRNA genes (Asp-trnD-GUC, Glu-trnE-UUC, Ile-trnI-CAU, Met-trnM-CAU, Met-trnfM-CAU, Trp-trnW-CCA, and Tyr-trnY-GUA), and three rRNA genes. Among them, ccmFC and cox1 have four copies, rrn5, rrn18, and rrn26 have three copies, nad1, trnD-GUC, trnE-UUC, trnM-CAU, and trnfM-CAU have two copies. We found five cis-spliced introns and 10 trans-spliced introns in the protein-coding genes (Supplementary Table S5). Both the number of protein-coding genes (Figure 2B) and the number of cis- and trans-spliced introns identified in P. orientalis are similar to that of Cupressophytes taxa reported previously (Guo et al., 2020). We also detected seven complete chloroplast genes (although they may not be functional), atpA, ndhA, ndhB, ndhC, ndhD, ndhE, and ndhH, and two chloroplast pseudogenes, petB and petD, in the mitogenome of P. orientalis (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-g002.jpg)
FIGURE 2. Phylogeny and gene content of the 20 mitogenomes included in this study. (A) The phylogenetic tree among major gymnosperm lineages based on 28 shared genes. Numbers next to nodes indicate bootstrap values from 1,000 replicates, and blue bars correspond to the 95% highest posterior density (HPD) of divergence time. Period is shown at the bottom of the tree. Ca., Carboniferous; Pe. Permian; Tr. Triassic; Ju. Jurassic; Cr. Cretaceous; Pa. Paleogene; N. Neogene. (B) Gene contents in the 20 mitogenomes. The 28 shared genes are atp[1, 4, 6, 8, 9], ccm[B, C, Fc, Fn], cob, cox[1–3], matR, mttB, nad[1, 2, 3, 4, 4L, 5, 6, 7, 9], rps[3, 4, 12], and sdh4. Intact genes are indicated by dark blue blocks and missing genes in white. The lost genes that are confirmed as transferred to nucleus are highlighted in orange. Pseudogenes are indicated as ψ. (C) RNA editing sites in the 20 mitogenomes. The two plots in the far right show the number of editing sites for all genes and the 28 shared genes, respectively. (D) Nucleotide substitution rates in the 20 mitogenomes. Branch lengths are proportional to rates of synonymous (dS) and nonsynonymous (dN) substitutions inferred by the maximum-likelihood method implemented in PAML.
3.2 RNA Editing, Mitochondrial Genes-Based Phylogeny, and Variable Substitution Rates in Cupressophytes and Gnetophytes
RNA editing remodels variability in higher plant mitogenomes, and could bias phylogenetic inferences if not accounted for (Szmidt et al., 2001; Mower, 2005). We predicted the RNA editing sites in protein coding genes of all 20 species and found that Ginkgo (1,300 sites) and Gnetophytes (225–336 sites) had the highest and lowest C-to-U editing sites across all protein-coding regions, respectively, and the other species spanned a range of 520–1,220 sites (Figure 2C). This variation in editing frequency reflects a lineage-specific pattern of editing, for example 25% (254) and 44–45% (102–147) of the editing sites in Ginkgo and Gnetophytes were unique (Figure 2C). In contrast, C. sempervirens and H. glabra of Cupressophytes each had only 6% (43–47) unique sites (Figure 2C).
To assess the accuracy of the predicted editing sites, we identified the edited sites from RNAseq data in P. orientalis. A total of 1,068 C-to-U edit sites were detected in 32 genes, of which 745 sites were among the 1,007 predicted sites (Table 2, Supplementary Table S6). This comparison suggests that the prediction method provides a reasonably reliable estimate of RNA editing (74% in P. orientalis), consistent with the results from a previous study (Guo et al., 2016).
We inferred the phylogeny of the major gymnosperm lineages using 28 conserved protein-coding genes common to all 20 species included in this study. The concatenated sequence alignment consisted 26,844 sites. On this phylogenetic tree, Cupressophytes and Gnetophytes are sister clades with a strong bootstrap support of 92% (Figure 2A). Within Cupressophytes, P. orientalis is sister to the group consisting of C. sempervirens and H. glabra (Figure 2A), and diverged from this group 87 MYA in the Paleogene. The divergence of the Cupressaceae family is dated to 145 MYA (Figure 2A), which is similar to the 159 MYA estimate from nuclear data (De La Torre et al., 2017).
To evaluate the substitution rates in mitochondrial genes, we estimated synonymous (dS) and nonsynonymous (dN) substitution rates along the branches of the inferred phylogeny. We found that dS and dN varied up to 46- and 234-fold among species, respectively. This wide range is mainly caused by Cupressophytes and Gnetophytes with their substantially longer branches and larger variation in dS and dN among species (Figure 2D; Table 3) compared to the other clades, all of which have low substitution rates. The overall dS of Cupressophytes and Gnetophytes is more than 5-fold higher than Pinaceae. The same overall pattern is seen with dN, although it is about 20–70% of the dS values, reflecting selective constraints on nonsynonymous mutations (Figure 2D; Table 3). Within Cupressophytes, dS varied by 46-fold and dN varied by 25-fold among lineages, with P. orientalis at the low end for both values (dS = 0.0030, dN = 0.0015; Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-t003.jpg)
TABLE 3. Relative (dS, dN) and absolute (RS, RN) synonymous and nonsynonymous substitution rates. MYA, million years ago.
Moreover, we calculated the absolute substitution rates RS and RN per branch by adjusting the dS and dN by dividing the divergence time of the branch. Again, RS and RN varied up to 118- and 171-fold among species, respectively, with the highest RS and RN in Cupressophytes and Gnetophytes, and the lowest in Pinaceae. Within Cupressophytes, RS and RN varied by 40-fold and 70-fold, with P. orientalis having the lowest rates (Table 3). This wide range of substitution rates indicates that sequence evolution is heterogeneous within and among plant lineages, especially in Cupressophytes.
3.3 Mitochondrial Gene Loss and Transfer to Nucleus in Cupressophytes and Gnetophytes
Early studies suggest that the mitogenome of the common ancestor of seed plants probably contained 41 protein-coding genes (Richardson et al., 2013; Guo et al., 2016). We found all 41 genes in Cycads, Ginkgo and Pinaceae, but 8–12 genes were absent in Cupressophytes and Gnetophytes (white and orange boxes in Figure 2B). All species in Cupressophytes have lost the same set of eight ribosomal protein genes, while Cupressaceae species and T. cuspidata lost an additional sdh3 gene (Figure 2B). In Gnetophytes, both G. gnemon and W. mirabilis have lost 11 genes, including 10 ribosomal protein genes and the sdh3 gene; their rpl16 appeared as degraded pseudogene (Figure 2B).
The reduction in mitochondrial protein-coding genes in Cupressophytes and Gnetophytes has been suggested to be a result of gene transfer to the nucleus prior to their loss from the mitogenome (Guo et al., 2016; Guo et al., 2020). To verify whether the lost genes were transferred to the nucleus, we first searched for them in the transcriptomes of Cupressophytes and Gnetophytes (Supplementary Table S2). We recovered 4–7 intact lost mitochondrial genes in transcriptomes of each respective species, but none in A. heterophylla (Figure 2B). The recovery of few intact mitochondrial genes in A. heterophylla is likely to be a consequence of the limited numbers of ESTs or cDNA sequences in the database, which would hamper detection of these genes in its nuclear genome. We then searched for these recovered genes in their own or related species’ reference genomes and found all of them with high identity (89–100%) and coverage (85–100%) scores (Supplementary Table S3). There were a few exceptions, e.g., rps11 and sdh3 showed lower identity (52–75%) but high coverage (92–100%) scores, likely due to different substitution rates between genes and species. The BLAST hits of Podocarpus macrophyllus genes gave moderate identity scores (48–81%) likely because of the lack of closely related reference. We further analyzed the mapping depth of genome resequencing reads to the mitogenes and transcriptome-recovered genes in each Cupressophytes and Gnetophytes species. In P. orientalis, the median depth of read coverage of the 32 genes in mitogenome ranged 275–1,169×, while the nuclear gene LEAFY had 15× and the chloroplast matK 14,395× (Supplementary Figure S1). The six recovered mitogenes had median depth of read coverage 11–17×, similar to the LEAFY, supporting their location in nucleus. Using the same approach, we confirmed that 4–7 of the lost genes in the seven other species of Cupressophytes and Gnetophytes were transferred to the nucleus (orange boxes in Figure 2B). Our results illustrate the differential abundance of nuclear, mitochondrial and plastid genomes in plant cells, and demonstrate the power of depth of read coverage as a quantitative tool for verification of genome origin of a gene.
3.4 Sequence Turnover and Structural Rearrangement in Gymnosperm Mitogenomes
Rapid structural evolution is thought to be characteristic of plant mitogenomes (Palmer and Herbon, 1988; Cole et al., 2018). To understand the rate of sequence turnover and structural rearrangement in gymnosperms, we examined the degree of synteny and shared DNA among the mitogenomes of the main gymnosperms lineages. We found that syntenic blocks and shared DNA between genomes reduced rapidly with phylogenetic distance (Figures 3A–D; Table 4). For example, P. orientalis and C. sempervirens diversified for 87 million years (MY) shared ∼60% DNA sequence, but this value dropped sharply to only ∼10% between P. macrophyllus and A. heterophylla with a divergence time of 126 MY (Figures 3A,C; Table 4). We examined the patterns of shared DNA between 16 pairs of gymnosperm species, and found a strong negative correlation (r = 0.53, p = 0.03) with time of divergence (Figure 3F). However, as an outlier, exceptionally low DNA turnover was found between G. biloba and C. taitungensis (Table 4), similar to the report by Guo et al. (2016). In addition, high proportions of shared DNA (58–64%) were retained between C. taitungensis and Ceratozamia hildae/Zamia integrifolia with the divergence time of 194 MY (Figure 3E; Table 4). These results illustrate a variable rate of mitogenome turnover in gymnosperms.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-g003.jpg)
FIGURE 3. Mitogenome alignments between gymnosperm species showing the locally collinear blocks at varying divergence time. (A) Platycladus orientalis vs. Cupressus sempervirens. (B) Cupressus sempervirens vs. Hesperocyparis glabra. (C) Podocarpus macrophyllus vs. Araucaria heterophylla. (D) Gnetum gnemon vs. Welwitschia mirabilis. (E) Cycas taitungensis vs. Ceratozamia hildae. The color bars represent corresponding local collinear blocks, heights are proportional to pairwise sequence identity, and blocks below the center line indicate regions that align in the reverse complement (inverse) orientation. MYA, million years ago. (F) Correlation between divergence time and the degree of shared sequences. The blue line represents the regression fitting with a linear model. The grey shadow is the 95% confidence interval of the linear regression.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-t004.jpg)
TABLE 4. Amount of shared mitochondrial DNA in gymnosperms spanning a wide range of divergence times. MYA, million years ago.
3.5 Repetitive Sequences in Gymnosperm Mitogenomes
Repetitive sequences, particularly TEs, are important components of plant nuclear genomes and contribute to genome size variation and functional regulation (Mehrotra and Goyal, 2014). In contrast, TEs and other repetitive elements in plant mitogenomes are poorly characterized. In this study, we annotated TEs, tandem repeats, and dispersed repeats in the 20 mitogenomes, and found that they have a strong correlation with the genome size and are highly variable among species (Figure 4; Table 5). The highest proportion of TE is in A. sibirica (19.26%), followed by Z. integrifolia (15.57%), and the lowest in W. mirabilis (1.06%). TE contents vary widely even in the same clade, e.g. from 1.58% (P. macrophyllus) to 6.80% (H. glabra) in Cupressophytes, from 5.29% (P. lambertiana) to 19.26% (A. sibirica) in Pinaceae, and from 2.52% (C. taitungensis) to 15.57% (Z. integrifolia) in Cycads. Among TEs, LTRs (long terminal repeats) are the most common group followed by LINE (long interspersed nuclear elements) in most species (Figure 4A). For the three Cupressaceae species and the two Pinus species, the most abundant element is the “Unknown” category (Figure 4A). The “Unknown” category of TEs could be due to accumulation of various mutations over long evolutionary time that have degraded the TEs into sequences that cannot be properly classified.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-g004.jpg)
FIGURE 4. Composition of repetitive sequences in gymnosperm mitogenomes. (A) The proportion of TE classes identified in each genome. LTR, long terminal repeat; LINE, long interspersed nuclear element; SINE, short interspersed element; DNA, DNA transposable element; Unknown, unclassified transposable element. (B) Correlation between mitogenome size and tandem repeat size, (C) TE size, (D) dispersed repeat size, and (E) repetitive sequence size. The blue line represents the regression fitting with a linear model. The grey shadow is the 95% confidence interval of the linear regression.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-t005.jpg)
TABLE 5. Repetitive sequences in each mitogenome. Dispersed repeats are divided into small-size (<100 bp), medium-size (100–1,000 bp), and large-size (≥1,000 bp) classes. Note, overlapping repetitive sequences are counted only once when calculating their total size.
In addition to TEs, dispersed repeats, including inverted and direct orientation duplicated sequences, are prevalent and made up to 20–30% of the mitogenomes, with a few exceptions (Table 5). The highest proportion is found in Cycads (22–32%) and A. sibirica in Pinaceae (32%). The two species in Gnetophytes are highly variable in this regard with 6% in W. mirabilis but 27% in G. gnemon; this value was stable (21–22%) among the three taxa in Cupressaceae. Tandem repeats are a smaller component compared to dispersed repeats and TEs, and were observed at the level of 1–5% in most of the mitogenomes (Table 5). The overall repetitive sequences composed 20–38% of the analyzed mitogenomes with a few exceptions of 9–17%, and show a strong correlation with genome size (r = 0.99, p < 0.0001; Figure 4E; Table 5).
3.6 Accumulation and Elimination of Dispersed Repeats in Mitogenomes
To understand the dynamics of dispersed repeats in mitogenomes, we divided the repeats into small-, medium-, and large-size classes with lengths of <100 bp, 100–1,000 bp, and ≥1,000 bp, respectively. We found that small-size repeats are the most abundant class in the 20 mitogenomes, followed by the medium-size repeats, while large-size repeats are rare (Table 5). The median lengths for the three classes are 52 bp, 136 bp, and 1,517 bp, respectively. The number of small-size repeats are 2–20 fold more common than medium-size repeats; G. biloba was extreme having 112-fold more small-size than medium-size repeats (Table 5). In A. sibirica, the two classes are balanced. Large-size dispersed repeats are limited, with 50% of the species containing fewer than 15 large repeats (Table 5). On the other hand, three species A. sibirica, H. glabra and L. sibirica, had more than 100 large repeats (Table 5). Thus, the density plots for small-size and medium-size repeats are similar among species, while it varied widely for the large-size repeats (Figure 5A). The GC contents are relatively conserved for the three classes of dispersed repeats (Figure 5B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/867736/fgene-13-867736-HTML/image_m/fgene-13-867736-g005.jpg)
FIGURE 5. The distribution of repeat size, GC content, and Kimuta 2-parameter (K2P) distance of dispersed repeats in each genome. (A) The density plot showing the length distribution of three classes of dispersed repeats. (B) The distribution of GC content in each size class. (C) The distribution of K2P distances among dispersed repeats. Repeats are grouped into small-size (<100 bp), medium-size (100–1,000 bp) and large-size (≥1,000 bp) classes. The K2P distance of 0 and 0.2 are shown with vertical red dash lines. Note, Taxus cuspidata contains two large-size dispersed repeats with K2P values of 0 and 0.0046, respectively, which failed to show in the density plot (C)
To understand the tempo of evolution in dispersed repeats, we analyzed the distribution of K2P distances among the repeats in the 20 mitogenomes. Small values of K2P represent recent duplications while large values represent more ancient duplication events. We divided the K2P values into three distance classes: ∼0 as recent duplications, >0 and <0.2 as intermediate events, and >0.2 as ancient insertions. A K2P distribution with a dominant narrow peak indicates rapid accumulation of repeats in a short period, otherwise a flattened peak indicates a slow generation of repeats over a prolonged period. As shown in Figure 5C, the K2P distributions of the small-size class is similar among all species with a dominant recent peak and lack of ancient insertions, suggesting they have been accumulated recently over a short time span. In contrast, the K2P distribution of the medium-size and large-size repeats are visibly different among species, often with two or three peaks spanning recent to ancient classes, suggesting multiple and distant expansion events.
Ancient dispersed repeats are easy to identify based on K2P distances, but they are not prevalent in the 20 mitogenomes. This may suggest that the persistance of ancient large-size repeats is governed by the rates of accumulation and elimination. Large repeats can be eliminated or broken down by repeat-mediated recombination, resulting in an increase in small- or medium-size repeats that appear to be more recent expansions (Maréchal and Brisson, 2010; Sullivan et al., 2020).
4 Discussion
Assembly of gymnosperm mitogenomes from genomic sequence data is challenging because of complex mechanisms contributing to structural and sequence variation, and lack of well-developed databases and tools for accurate genome mapping. Up to now, only ∼20 gymnosperm mitogenomes have been released. This study contributes a complete mitogenome sequence for Cupressaceae and expands the coverage in gymnosperms. With this new genome, we performed in-depth comparative analyses among major gymnosperm lineages to understand the tempo and mode of mitogenome evolution.
4.1 Heterogeneity in Repetitive Elements in Mitogenomes
Among the available gymnosperm mitogenomes, genome sizes vary markedly by 30-fold (Table 1). This large variation is hypothesized to be due to lineage- or species-specific accumulation of repetitive sequences, horizontal transfers and intergenomic transfers (Goremykin et al., 2012; Rice et al., 2013; Guo et al., 2017). It has been reported that 29–100% of sequences in gymnosperm mitogenomes are similar to the nuclear genome sequences (Sullivan et al., 2020). In P. orientalis we found 76% undefined sequences (Table 2), which could in part have been derived from the nucleus by intracellular transfer.
Another source of mitogenome size variation is repetitive sequences. Estimations of repetitive contents can be affected by filter and search criteria, e.g. setting the length to 50 bp or 100 bp could underestimate the full size range of repeats. In this study, we searched for dispersed repeats in all size classes. We found that dispersed repeats occupied 6–32% and TEs 1–19% in the sampled mitogenomes. The wide range of dispersed repeats content quantified in this study mirror the results of 1–49% found in 82 angiosperms (Dong et al., 2018). All categories of repetitive sequences are strongly correlated with genome size, giving an overall relationship of r = 0.99. Our results illustrate that repetitive sequences, especially TEs and dispersed repeats, are major components of mitogenome expansion.
Although a major component of mitogenomes, the size distributions of repetitive sequences and the dynamics of their accumulation have rarely been investigated in gymnosperms. We analyzed the frequency and sequence diversity of different size classes of repeats, and found small-size repeats are most common while large-size repeats are rare and vary widely among species. K2P distances among repeats indicate that most of the repeats have accumulated by recent expansions. Large-size repeats usually tend to recombine more frequently than small repeats (Skippington et al., 2015; Guo et al., 2016), which may eliminate large repeats over time. However, in P. abies, most recombinogenic repeats are small- to medium-size, which deviates from the general expectation that, in the absence of other factors, recombination should scale positively with repeat length (Sullivan et al., 2020). We are of the opinion that genome- or region-specific recombination rates could lead to unique sequence turnover rates resulting in the heterogeneous distribution of large repeats in mitogenomes.
4.2 Extensive Rearrangement and Sequence Turnover in Mitogenomes
The plant mitogenome is characterized by rapid structural rearrangement and sequence turnover, driven by a high recombination activity, making mtDNA highly variable even between closely related species (Guo et al., 2016; Cole et al., 2018; Sullivan et al., 2020). In line with this view, extensive rearrangement is found between two spruce species, P. abies and P. glauca (Sullivan et al., 2020). On the other hand, the mitogenomes of G. biloba and C. taitungensis have retained substantial amounts of shared DNA despite their divergence time >300 MY (Table 4; Guo et al., 2016). Even with these extreme cases, we found that the synteny blocks and shared DNA between mitogenomes generally decayed with divergence time, with a negative correlation r = 0.53 (p = 0.03) among the gymnosperms (Figure 3F). Although not exhaustive, our analysis provides a quantitative estimation of the rate of sequence turnover in gymnosperm mitogenomes. Previous studies show that rearrangement in mitogenomes occurs largely through repeat-mediated recombination (Cole et al., 2018; Xia et al., 2020). However, the relationship between inter- and intra-molecule recombination and rearrangement rates among species is unclear (Cole et al., 2018). Recombination dynamics in plant mitogenomes are poorly understood to date and deserve further scrutiny as it drives sequence and structural evolution of plant mitogenomes, and affects genome integrity and biological function.
4.3 mtDNA Evidence Support Cupressophytes and Gnetophytes as Sister Groups
The phylogenetic placement of Gnetophytes differs among studies. Five conflicting hypotheses were supported by different marker systems and include: the Gnetophytes-other seed plant hypothesis, the Gnetophytes-other gymnosperm hypothesis, the Gnetifer (Gnetophytes + Conifers) hypothesis, the Gnepine (Gnetophytes + Pinaceae) hypothesis, and the GneCup (Gnetophytes + Cupressophytes) hypothesis (Ran et al., 2018). Yet another suggestion is that Gnetophytes is sister to or within the conifers (Wan et al., 2018). The Gnepine hypothesis is favored by most studies (Ran et al., 2018), including a recent phylogeny based on mitochondrial protein-coding genes (Kan et al., 2021). On the other hand, plastid DNA strongly supports the GneCup hypothesis (Ruhfel et al., 2014; Gitzendanner et al., 2018; Leebens-Mack et al., 2019). This incongruence between datasets might be due to differences in the substitution rates of mitochondrial, plastid and nuclear genes, as well as the amount of informative sites included in each study.
In this study, we inferred the relationships among major gymnosperm lineages using 28 conserved protein-coding genes. The recovered phylogeny supports the GneCup hypothesis with a strong bootstrap support. We are more in favor of GneCup hypothesis because of the shared gene transfer pattern in these two clades. Over evolutionary time, some segments of mtDNA have diverged so much that groups of genes have been lost. The common ancestor of seed plants is inferred to contain 41 protein-coding genes (Guo et al., 2016; Mower, 2020), and Cycads, Ginkgo, and Pinaceae maintained this same set. The mitogenomes of Cupressophytes and Gnetophytes stand out from other gymnosperms in that both have lost 8–11 ribosomal protein-coding genes and the sdh3 gene, with the exception of Ephedra przewalskii (Gnetophytes) which has lost 19 genes (Guo et al., 2020; Kan et al., 2021). To confirm whether these lost genes have been transferred to nucleus, we first annotated them in transcriptomes, then identified their homologs in nuclear genomes, and further examined their depth of mapped reads using genome resequencing data. All three lines of evidence support the transfer of 4–7 mitogenes to nucleus. Our findings are in good agreement with the recent studies by Kan et al. (2020) and Kan et al. (2021), in which they report finding 6–8 intact or partial homologs of lost genes in the nuclear genomes. In sharp contrast, the three other lineages of gymnosperms, including Pinaceae, Ginkgo, and Cycads, share the same set of genes as the common ancestor more than 300 MYA (Guo et al., 2020). The unrecovered missing mitogenes could either be truly lost or not present in the transcriptomes we analyzed due to tissue- and/or time-specific expression, or incomplete annotation from transcriptomes, which identified the genes as partial homologs instead of intact genes.
Other striking characteristics of Cupressophyte and Gnetophyte mtDNA are their accelerated and hugely variable substitution rates across the 28 conserved protein-coding genes relative to other gymnosperms. The synonymous substitution rate dS in Cupressophytes and Gnetophytes averaged 0.0641 and 0.0536, respectively. After scaling up by divergence time, the absolute synonymous substitution rates RS were similar in Cupressophytes and Gnetophytes with 2.5 × 10−10 and 2.1 × 10−10 site/year, respectively, and are ∼2–7 fold higher than Pinaceae, Ginkgo, and Cycads. Previous work similarly inferred exceptionally low substitution rates in Pinaceae (Pinus and Picea) and higher substitution rates in Gnetophytes species (Guo et al., 2016; Sullivan et al., 2020; Kan et al., 2021). However, the extraordinary variation in RS in Cupressophytes and Gnetophytes is revealed only by denser sampling of taxa from each clade. Within Cupressophytes, RS varied 46-fold among taxa. Our results add to the recognition that while most plant mitogenomes show very low synonymous substitution rates, this rate of divergence is also exceptionally variable for reasons still unclear (Mower et al., 2007). Another unexpected finding is that within Gnetophytes, G. gnemon and W. mirabilis showed distinct species-specific RNA editing patters as 87–91% of their editing sites differ from each other. Such a striking diversification in RNA editing is not observed within other clades. Future studies with broader taxa sampling are required to understand the mechanisms driving these patterns, including gene loss and transfer, TEs, recombination and selection.
5 Conclusion
This study contributes a complete mitogenome assembly of P. orientalis to the still very limited accessions in gymnosperms. Our comparative analyses with 19 other mitogenomes of seed plants characterized the composition and distribution of repetitive elements, the tempo of sequence turnover and structural rearrangement, and the frequency of RNA editing in protein-coding genes. Our study revealed shared patterns of cyto-unclear gene transfer and accelerated substitution rates in Cupressophytes and Gnetophytes, and lend support to their sister group placement within the gymnosperms. Our findings highlight and reinforce the dynamic and enigmatic evolution of mitogenomes in gymnosperms.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, OL703044-OL703045.
Author Contributions
WZ, HL and X-RW designed the study; HL, WZ, R-GZ and J-FM generated and analyzed data; HL, WZ, and X-RW wrote the manuscript with contributions from all authors.
Funding
This study was supported by grants from the National Natural Science Foundation of China (No. 32171816), and the Swedish Research Council (VR 2017-04686).
Conflict of Interest
Author R-GZ was employed by the company Ori (Shandong) Gene Science and Technology Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Dr. Jeffrey Mower for sharing data with us, and Prof. Barbara Giles for valuable comments and linguistic editing. We also thank one of the reviewers for constructive suggestions that improved the clarity of the manuscript. Genomic data processing and analyses were performed using resources provided by the Swedish National Infrastructure for Computing (SNIC), through the High Performance Computing Centre North (HPC2N), Umeå University.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.867736/full#supplementary-material
References
Benson, G. (1999). Tandem Repeats Finder: A Program to Analyze DNA Sequences. Nucleic Acids Res. 27 (2), 573–580. doi:10.1093/nar/27.2.573
Chaw, S.-M., Chun-Chieh Shih, A., Wang, D., Wu, Y.-W., Liu, S.-M., and Chou, T.-Y. (2008). The Mitochondrial Genome of the Gymnosperm Cycas Taitungensis Contains a Novel Family of Short Interspersed Elements, Bpu Sequences, and Abundant RNA Editing Sites. Mol. Biol. Evol. 25 (3), 603–615. doi:10.1093/molbev/msn009
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: An Ultra-fast All-In-One FASTQ Preprocessor. Bioinformatics 34 (17), i884–i890. doi:10.1093/bioinformatics/bty560
Cheng, J., Wang, X., Liu, X., Zhu, X., Li, Z., Chu, H., et al. (2021). Chromosome-level Genome of Himalayan Yew Provides Insights into the Origin and Evolution of the Paclitaxel Biosynthetic Pathway. Mol. Plant 14 (7), 1199–1209. doi:10.1016/j.molp.2021.04.015
Cole, L. W., Guo, W., Mower, J. P., and Palmer, J. D. (2018). High and Variable Rates of Repeat-Mediated Mitochondrial Genome Rearrangement in a Genus of Plants. Mol. Biol. Evol. 35 (11), 2773–2785. doi:10.1093/molbev/msy176
Criscuolo, A., and Gribaldo, S. (2010). BMGE (Block Mapping and Gathering with Entropy): A New Software for Selection of Phylogenetic Informative Regions from Multiple Sequence Alignments. BMC Evol. Biol. 10 (1), 210. doi:10.1186/1471-2148-10-210
Darling, A. E., Mau, B., and Perna, N. T. (2010). ProgressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement. PLoS One 5 (6), e11147. doi:10.1371/journal.pone.0011147
De La Torre, A. R., Li, Z., Van de Peer, Y., and Ingvarsson, P. K. (2017). Contrasting Rates of Molecular Evolution and Patterns of Selection Among Gymnosperms and Flowering Plants. Mol. Biol. Evol. 34 (6), 1363–1377. doi:10.1093/molbev/msx069
Dong, S., Zhao, C., Chen, F., Liu, Y., Zhang, S., Wu, H., et al. (2018). The Complete Mitochondrial Genome of the Early Flowering Plant Nymphaea Colorata Is Highly Repetitive with Low Recombination. BMC Genomics 19 (1), 614. doi:10.1186/s12864-018-4991-4
Doyle, J. J., and Doyle, J. L. (1987). A Rapid DNA Isolation Procedure for Small Quantities of Fresh Leaf Tissue. Phytochem. Bull. 19, 11–15.
Gitzendanner, M. A., Soltis, P. S., Wong, G. K.-S., Ruhfel, B. R., and Soltis, D. E. (2018). Plastid Phylogenomic Analysis of Green Plants: A Billion Years of Evolutionary History. Am. J. Bot. 105 (3), 291–301. doi:10.1002/ajb2.1048
Goremykin, V. V., Lockhart, P. J., Viola, R., and Velasco, R. (2012). The Mitochondrial Genome of Malus Domestica and the Import-Driven Hypothesis of Mitochondrial Genome Expansion in Seed Plants. Plant J. 71 (4), 615–626. doi:10.1111/j.1365-313X.2012.05014.x
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length Transcriptome Assembly from RNA-Seq Data without a Reference Genome. Nat. Biotechnol. 29 (7), 644–652. doi:10.1038/nbt.1883
Guo, W., Grewe, F., Fan, W., Young, G. J., Knoop, V., Palmer, J. D., et al. (2016). Ginkgo and Welwitschia Mitogenome Reveal Extreme Contrasts in Gymnosperm Mitochondrial Evolution. Mol. Biol. Evol. 33 (6), 1448–1460. doi:10.1093/molbev/msw024
Guo, W., Zhu, A., Fan, W., Adams, R. P., and Mower, J. P. (2020). Extensive Shifts from Cis- to Trans-splicing of Gymnosperm Mitochondrial Introns. Mol. Biol. Evol. 37 (6), 1615–1620. doi:10.1093/molbev/msaa029
Guo, W., Zhu, A., Fan, W., and Mower, J. P. (2017). Complete Mitochondrial Genomes from the Ferns Ophioglossum Californicum and Psilotum Nudum Are Highly Repetitive with the Largest Organellar Introns. New Phytol. 213 (1), 391–403. doi:10.1111/nph.14135
Hu, X.-G., Liu, H., Jin, Y., Sun, Y.-Q., Li, Y., Zhao, W., et al. (2016). De Novo transcriptome Assembly and Characterization for the Widespread and Stress-Tolerant Conifer Platycladus Orientalis. PLoS One 11 (2), e0148985. doi:10.1371/journal.pone.0148985
Jackman, S. D., Warren, R. L., Gibb, E. A., Vandervalk, B. P., Mohamadi, H., Chu, J., et al. (2015). Organellar Genomes of White Spruce (Picea glauca): Assembly and Annotation. Genome Biol. Evol. 8 (1), 29–41. doi:10.1093/gbe/evv244
Jin, J.-J., Yu, W.-B., Yang, J.-B., Song, Y., dePamphilis, C. W., Yi, T.-S., et al. (2020). GetOrganelle: A Fast and Versatile Toolkit for Accurate De Novo Assembly of Organelle Genomes. Genome Biol. 21 (1), 241. doi:10.1186/s13059-020-02154-5
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: Fast Model Selection for Accurate Phylogenetic Estimates. Nat. Methods 14 (6), 587–589. doi:10.1038/nmeth.4285
Kan, S.-L., Shen, T.-T., Gong, P., Ran, J.-H., and Wang, X.-Q. (2020). The Complete Mitochondrial Genome of Taxus Cuspidata (Taxaceae): Eight Protein-Coding Genes Have Transferred to the Nuclear Genome. BMC Evol. Biol. 20 (1), 10. doi:10.1186/s12862-020-1582-1
Kan, S.-L., Shen, T.-T., Ran, J.-H., and Wang, X.-Q. (2021). Both Conifer II and Gnetales Are Characterized by a High Frequency of Ancient Mitochondrial Gene Transfer to the Nuclear Genome. BMC Biol. 19 (1), 146. doi:10.1186/s12915-021-01096-z
Kent, W. J. (2002). BLAT-The BLAST-like Alignment Tool. Genome Res. 12 (4), 656–664. doi:10.1101/gr.229202
Kimura, M. (1980). A Simple Method for Estimating Evolutionary Rates of Base Substitutions through Comparative Studies of Nucleotide Sequences. J. Mol. Evol. 16 (2), 111–120. doi:10.1007/BF01731581
Leebens-Mack, J. H., Barker, M. S., Carpenter, E. J., Deyholos, M. K., Gitzendanner, M. A., Graham, S. W., et al. (2019). One Thousand Plant Transcriptomes and the Phylogenomics of Green Plants. Nature 574 (7780), 679–685. doi:10.1038/s41586-019-1693-2
Li, H. (2013). Aligning Sequence Reads, Clone Sequences and Assembly Contigs with BWA-MEM. Available at: https://arxiv.org/abs/1303.3997v2 [q-bio.GN] (Accessed June 21, 2018).
Li, H. (2018). Minimap2: Pairwise Alignment for Nucleotide Sequences. Bioinformatics 34 (18), 3094–3100. doi:10.1093/bioinformatics/bty191
Lo Giudice, C., Tangaro, M. A., Pesole, G., and Picardi, E. (2020). Investigating RNA Editing in Deep Transcriptome Datasets with REDItools and REDIportal. Nat. Protoc. 15 (3), 1098–1131. doi:10.1038/s41596-019-0279-7
Lowe, T. M., and Eddy, S. R. (1997). tRNAscan-SE: A Program for Improved Detection of Transfer RNA Genes in Genomic Sequence. Nucleic Acids Res. 25 (5), 955–964. doi:10.1093/nar/25.5.955
Maréchal, A., and Brisson, N. (2010). Recombination and the Maintenance of Plant Organelle Genome Stability. New Phytol. 186 (2), 299–317. doi:10.1111/j.1469-8137.2010.03195.x
Marks, R. A., Hotaling, S., Frandsen, P. B., and VanBuren, R. (2021). Representation and Participation across 20 Years of Plant Genome Sequencing. Nat. Plants 7, 1571–1578. doi:10.1038/s41477-021-01031-8
Mehrotra, S., and Goyal, V. (2014). Repetitive Sequences in Plant Nuclear DNA: Types, Distribution, Evolution and Function. Genomics, Proteomics Bioinforma. 12 (4), 164–171. doi:10.1016/j.gpb.2014.07.003
Mower, J. P. (2005). PREP-mt: Predictive RNA Editor for Plant Mitochondrial Genes. BMC Bioinforma. 6 (1), 96. doi:10.1186/1471-2105-6-96
Mower, J. P., Touzet, P., Gummow, J. S., Delph, L. F., and Palmer, J. D. (2007). Extensive Variation in Synonymous Substitution Rates in Mitochondrial Genes of Seed Plants. BMC Evol. Biol. 7 (1), 135. doi:10.1186/1471-2148-7-135
Mower, J. P. (2020). Variation in Protein Gene and Intron Content Among Land Plant Mitogenomes. Mitochondrion 53, 203–213. doi:10.1016/j.mito.2020.06.002
Neale, D. B., Zimin, A. V., Zaman, S., Scott, A. D., Shrestha, B., Workman, R. E., et al. (2021). Assembled and Annotated 26.5 Gbp Coast Redwood Genome: A Resource for Estimating Evolutionary Adaptive Potential and Investigating Hexaploid Origin. G3 Genes|Genomes|Genetics 12 (1), jkab380. doi:10.1093/g3journal/jkab380
Nguyen, L.-T., Schmidt, H. A., Von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 32 (1), 268–274. doi:10.1093/molbev/msu300
Nystedt, B., Street, N. R., Wetterbom, A., Zuccolo, A., Lin, Y.-C., Scofield, D. G., et al. (2013). The Norway Spruce Genome Sequence and Conifer Genome Evolution. Nature 497 (7451), 579–584. doi:10.1038/nature12211
Palmer, J. D., and Herbon, L. A. (1988). Plant Mitochondrial DNA Evolved Rapidly in Structure, but Slowly in Sequence. J. Mol. Evol. 28 (1), 87–97. doi:10.1007/BF02143500
Petit, R. J., and Vendramin, G. G. (2007). “Plant Phylogeography Based on Organelle Genes: an Introduction,” in Phylogeography of Southern European Refugia: Evolutionary Perspectives on the Origins and Conservation of European Biodiversity. Editors S. Weiss, and N. Ferrandeds (Dordrecht: Springer Netherlands), 23–97. doi:10.1007/1-4020-4904-8_2
Putintseva, Y. A., Bondar, E. I., Simonov, E. P., Sharov, V. V., Oreshkova, N. V., Kuzmin, D. A., et al. (2020). Siberian Larch (Larix Sibirica Ledeb.) Mitochondrial Genome Assembled Using Both Short and Long Nucleotide Sequence Reads Is Currently the Largest Known Mitogenome. BMC Genomics 21 (1), 654. doi:10.1186/s12864-020-07061-4
Ran, J.-H., Shen, T.-T., Wang, M.-M., and Wang, X.-Q. (2018). Phylogenomics Resolves the Deep Phylogeny of Seed Plants and Indicates Partial Convergent or Homoplastic Evolution between Gnetales and Angiosperms. Proc. R. Soc. B 285 (1881), 20181012. doi:10.1098/rspb.2018.1012
Rice, D. W., Alverson, A. J., Richardson, A. O., Young, G. J., Sanchez-Puerta, M. V., Munzinger, J., et al. (2013). Horizontal Transfer of Entire Genomes via Mitochondrial Fusion in the Angiosperm Amborella. Science 342 (6165), 1468–1473. doi:10.1126/science.1246275
Richardson, A. O., Rice, D. W., Young, G. J., Alverson, A. J., and Palmer, J. D. (2013). The “Fossilized” Mitochondrial Genome of Liriodendron Tulipifera: Ancestral Gene Content and Order, Ancestral Editing Sites, and Extraordinarily Low Mutation Rate. BMC Biol. 11 (1), 1–17. doi:10.1186/1741-7007-11-29
Ruhfel, B. R., Gitzendanner, M. A., Soltis, P. S., Soltis, D. E., and Burleigh, J. (2014). From Algae to Angiosperms-Inferring the Phylogeny of Green Plants (Viridiplantae) from 360 Plastid Genomes. BMC Evol. Biol. 14 (1), 23. doi:10.1186/1471-2148-14-23
Scott, A. D., Zimin, A. V., Puiu, D., Workman, R., Britton, M., Zaman, S., et al. (2020). A Reference Genome Sequence for Giant Sequoia. G3 Genes|Genomes|Genetics 10 (11), 3907–3919. doi:10.1534/g3.120.401612
Skippington, E., Barkman, T. J., Rice, D. W., and Palmer, J. D. (2015). Miniaturized Mitogenome of the Parasitic Plant Viscum Scurruloideum Is Extremely Divergent and Dynamic and Has Lost All Nad Genes. Proc. Natl. Acad. Sci. U.S.A. 112 (27), E3515–E3524. doi:10.1073/pnas.1504491112
Slater, G., and Birney, E. (2005). Automated Generation of Heuristics for Biological Sequence Comparison. BMC Bioinforma. 6, 31. doi:10.1186/1471-2105-6-31
Sloan, D. B., Alverson, A. J., Chuckalovcak, J. P., Wu, M., McCauley, D. E., Palmer, J. D., et al. (2012). Rapid Evolution of Enormous, Multichromosomal Genomes in Flowering Plant Mitochondria with Exceptionally High Mutation Rates. PLoS Biol. 10 (1), e1001241. doi:10.1371/journal.pbio.1001241
Smit, A., Hubley, R., and Green, P. (2013). RepeatMasker Open-4.0. 2013-2015. Available at: http://www.repeatmasker.org (accessed June 01, 2021).
Smit, A., and Hubley, R. (2008). RepeatModeler Open-1.0. 2008-2015. Available at: http://www.repeatmasker.org (accessed June 01, 2021).
Smith, D. R., and Keeling, P. J. (2015). Mitochondrial and Plastid Genome Architecture: Reoccurring Themes, but Significant Differences at the Extremes. Proc. Natl. Acad. Sci. U.S.A. 112 (33), 10177–10184. doi:10.1073/pnas.1422049112
Stanke, M., Diekhans, M., Baertsch, R., and Haussler, D. (2008). Using Native and Syntenically Mapped cDNA Alignments to Improve De Novo Gene Finding. Bioinformatics 24 (5), 637–644. doi:10.1093/bioinformatics/btn013
Sullivan, A. R., Eldfjell, Y., Schiffthaler, B., Delhomme, N., Asp, T., Hebelstrup, K. H., et al. (2020). The Mitogenome of Norway Spruce and a Reappraisal of Mitochondrial Recombination in Plants. Genome Biol. Evol. 12 (1), 3586–3598. doi:10.1093/gbe/evz263
Szmidt, A. E., Lu, M.-Z., and Wang, X.-R. (2001). Effects of RNA Editing on the coxI Evolution and Phylogeny Reconstruction. Euphytica 118 (1), 9–18. doi:10.1023/A:1004046220115
Wan, T., Liu, Z.-M., Li, L.-F., Leitch, A. R., Leitch, I. J., Lohaus, R., et al. (2018). A Genome for Gnetophytes and Early Evolution of Seed Plants. Nat. Plants 4 (2), 82–89. doi:10.1038/s41477-017-0097-2
Wan, T., Liu, Z., Leitch, I. J., Xin, H., Maggs-Kölling, G., Gong, Y., et al. (2021). The Welwitschia Genome Reveals a Unique Biology Underpinning Extreme Longevity in Deserts. Nat. Commun. 12, 4247. doi:10.1038/s41467-021-24528-4
Wang, X.-Q., and Ran, J.-H. (2014). Evolution and Biogeography of Gymnosperms. Mol. Phylogenetics Evol. 75, 24–40. doi:10.1016/j.ympev.2014.02.005
Wegrzyn, J. L., Lee, J. M., Tearse, B. R., and Neale, D. B. (2008). TreeGenes: A Forest Tree Genome Database. Int. J. Plant Genomics 2008, 1–7. doi:10.1155/2008/412875
Wick, R. R., Schultz, M. B., Zobel, J., and Holt, K. E. (2015). Bandage: Interactive Visualization of De Novo Genome Assemblies. Bioinformatics 31 (20), 3350–3352. doi:10.1093/bioinformatics/btv383
Wu, Z. Q., Liao, X. Z., Zhang, X. N., Tembrock, L. R., and Broz, A. (2020). Genomic Architectural Variation of Plant Mitochondria-A Review of Multichromosomal Structuring. J Sytematics Evol. 60, 160–168. doi:10.1111/jse.12655
Xia, H., Zhao, W., Shi, Y., Wang, X.-R., and Wang, B. (2020). Microhomologies Are Associated with Tandem Duplications and Structural Variation in Plant Mitochondrial Genomes. Genome Biol. Evol. 12 (11), 1965–1974. doi:10.1093/gbe/evaa172
Xiong, X., Gou, J., Liao, Q., Li, Y., Zhou, Q., Bi, G., et al. (2021). The Taxus Genome Provides Insights into Paclitaxel Biosynthesis. Nat. Plants 7 (8), 1026–1036. doi:10.1038/s41477-021-00963-5
Keywords: Platycladus orientalis, gene transfer, gymnosperm mitogenome, rearrangement, repetitive sequences, substitution rate, sequence turnover, RNA editing
Citation: Liu H, Zhao W, Zhang R-G, Mao J-F and Wang X-R (2022) Repetitive Elements, Sequence Turnover and Cyto-Nuclear Gene Transfer in Gymnosperm Mitogenomes. Front. Genet. 13:867736. doi: 10.3389/fgene.2022.867736
Received: 01 February 2022; Accepted: 28 April 2022;
Published: 25 May 2022.
Edited by:
Jinhui Chen, Nanjing Forestry University, ChinaReviewed by:
Peng-Fei Ma, Kunming Institute of Botany (CAS), ChinaJia-Yu Xue, Nanjing Agricultural University, China
Copyright © 2022 Liu, Zhao, Zhang, Mao and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Zhao, zhao.wei@umu.se
†These authors have contributed equally to this work