Next Article in Journal
Identification of Genes Involved in Fe–S Cluster Biosynthesis of Nitrogenase in Paenibacillus polymyxa WLY78
Next Article in Special Issue
Contribution of Massive Mitochondrial Fusion and Subsequent Fission in the Plant Life Cycle to the Integrity of the Mitochondrion and Its Genome
Previous Article in Journal
Human Amniotic Epithelial Cells as a Tool to Investigate the Effects of Cyanidin 3-O-Glucoside on Cell Differentiation
Previous Article in Special Issue
A PPR Protein ACM1 Is Involved in Chloroplast Gene Expression and Early Plastid Development in Arabidopsis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

1
Institute of Agriculture Science and Technology, Chungbuk National University, Cheongju, Chungbuk 28644, Korea
2
Department of Forest Science, Chungbuk National University, Cheongju, Chungbuk 28644, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(7), 3770; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22073770
Submission received: 3 February 2021 / Revised: 4 March 2021 / Accepted: 1 April 2021 / Published: 5 April 2021
(This article belongs to the Special Issue Organelle Genetics in Plants 2.0)

Abstract

:
Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.

1. Introduction

Caprifoliaceae sensu lato is one of two large clades in the order Dipsacales, which comprises Caprifoliaceae sensu stricto, Diervillaceae, Dipsacaceae, Linnaeaceae, Morinaceae, Valerianaceae, and Zabelia [1]. Among these taxa, the family Valerianaceae contains approximately 350 species of annual and perennial plants, which, with the exception of Australia and New Zealand, have an extensive worldwide distribution [2]. The members of this family are characterized by sympetalous, bilaterally symmetric to strongly asymmetric flowers; inferior tricarpellate ovaries with a single carpel at maturity; a single anatropous ovule; achene fruits; the absence of endosperm; and the presence of valeriana-epoxy-triesters (valepotriates) [2,3,4]. Given the established therapeutic properties of valepotriates, a distinct class of iridoid compounds, members of Valerianaceae have been used in traditional medicine for several hundred years [4].
The genus Valeriana L. is the largest among those comprising the family Valerianaceae, with species distributed throughout Asia, Europe, and America [5]. Among these, Valeriana sambucifolia Mikan is morphologically similar to Valeriana officinalis L., although it can be distinguished from the latter based on the small numbers of leaflets, aboveground runners [6], and the trichomes borne on fruits [2]. V. sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara (Figure 1), a broad-leaved valerian, is endemic to Ulleung Island, considered one of the hot spots of endemism in Korea [7]. Although this infra-specific taxon has for long been studied in a pharmacological context [8,9,10], there are numerous features of this plant that have yet to be sufficiently characterized. In contrast to V. sambucifolia, which is generally distributed in Scandinavia and Central and Atlantic Europe [11], V. sambucifolia f. dageletiana occurs only in Northeast Asia. Moreover, these plants differ with respect to the chromosome number, with that of the former and latter being 2n = 56 [12,13] and 2n = 54 or 55 [14], respectively. Despite it being taxonomically considered a forma of V. sambucifolia, based on the characteristics of a small number of broader leaflets (Figure 1A-B) and belowground runners (Figure 1C), we still have insufficient information about this pharmacologically attractive species [8].
Plants are characterized by two organellar genomes, namely the chloroplast genome (plastome) and the mitochondrial genome (mitogenome) [15], which in terrestrial plants have undergone independent evolution [16]. The plastomes of plants are generally conserved in terms of gene content, gene order, and genome length [17,18,19], although certain lineages, such as those of heterotrophs, have been found to have contracted genomes [20,21]. On the basis of these genomic characteristics, plastomes have been used as resources for resolving the phylogenetic relationships of uncertain lineages [22,23], providing an understanding of the transition from photosynthetic to non-photosynthetic modes of nutrition [24,25,26] and elucidating synapomorphies in the common ancestor of early diverged terrestrial plants [27,28]. In contrast to the relatively conserved plastomes, the mitogenomes of plants have undergone comparatively rapid structural evolution [29] by incorporating foreign DNA via horizontal [30] or intracellular [31,32] gene transfer. Consequently, there tends to be low similarity among the mitogenomes of terrestrial plants, even among closely related taxa [33,34,35].
Although the inter-organelle transfer of DNA from plastomes to mitogenomes has been extensively reported, transfer of DNA in the opposite direction, from mitogenomes to plastomes, appears to be a comparatively rare event that has only recently been discovered in some lineages. Moreover, whereas the former mitochondrial plastid DNAs (MTPTs) have been found to vary in length [16,36,37,38], the latter plastid mitochondrial DNAs (PTMTs) are generally smaller in length [39,40,41], although PTMTs of the Anacardium occidentale plastome appear to be a notable exception in this regard [42].
To date, the reference sequences of 4732 plastomes and 223 mitogenomes from terrestrial plants have been deposited in the NCBI database (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genome/ accessed on 5 April 2021). As these numbers indicate, the mitogenome sequence data are comparatively limited (only 5% compared to the plastome data), and currently there are no reference mitogenomes for plant orders such as Bruniales, Escalloniales, Paracryphiales, and Dipsacales among asterids, thereby hampering studies on mitogenome evolution and inter-organelle DNA transfer in these clades.
In this study, we sequenced the plastome and mitogenome of the Ulleung Island endemic V. sambucifolia f. dageletiana with the following aims: (1) to compare the genetic variation between the whole plastomes of two closely related valerians, V. officinalis and V. sambucifolia f. dageletiana, thereby characterizing this endemic species at the genomic level; (2) to produce a reference mitogenome sequences for the order Dipsacales (to the best of our knowledge, the first reported); and (3) to gain a comprehensive insight into inter-organelle DNA transfer in the order Dipsacales.

2. Results

2.1. Variation in the V. sambucifolia f. dageletiana Plastome

The average coverage of the V. sambucifolia f. dageletiana plastome was 2037.8 ± 233.6 (mean ± SD). The sequenced plastome has a length of 155,179 bp, including 85,334 bp large single-copy (LSC) and 15,243 bp small single-copy (SSC) regions and two inverted repeats (IRs) covering 27,301 bp. The plastome contains 134 genes, encoding 82 protein-coding genes, 8 ribosomal RNAs, 39 transfer RNAs, 2 pseudogenes (accD, ycf1), and 3 partial genes (ndhD and 2 trnR-ACG in the IR region). Compared with the plastome of V. officinalis, the most closely related taxa in the genera, the sequenced plastome differs in that the clpP, rps3, rps18, and ycf2 genes remain intact.
Notably, we found that the trnH gene, which, with the exception of certain Patrinia species, is generally located in the LSC region of Caprifoliaceae plants [43,44], has undergone translocation to the IR region in the plastome of V. sambucifolia f. dageletiana via IR expansion. Similar IR expansion and contraction at the boundaries between the IR and single-copy (SC) regions were found to account for differences in the lengths of the two Valeriana plastomes. With reference to the V. officinalis plastome structure, we assumed that IRb initially expanded to rpl32 over ndhF and that, subsequently, both the ndhF and rpl32 genes were eliminated from IRb and a copy of trnR-ACG was pseudogenized by a deletion in both IR regions. Finally, IRb expanded to within ndhD over ccsA and trnL. As a consequence, the IR region of V. sambucifolia f. dageletiana is now 3502 bp longer than that of V. officinalis (Figure 2).
Of the 155,179 bp of the full-length V. sambucifolia f. dageletiana plastome, 10 regions, comprising 152,753 bp, were found to correspond to the plastid DNAs of other viridiplantae. However, among the unmatched sequences, we identified eight non-plastome regions (NPRs), two of which are duplicated in the IR regions, that did not correspond to any other viridiplantae plastid DNAs: one of them is located in rbcL-accD (LSC region), two are positioned in trnL-trnL (IR regions), and five are found in ycf1 (SSC region). Moreover, we were unable to match six of these NPRs to any other nucleotide sequence. The 183 bp low-complexity region of NPR_4 was found to correspond to certain nuclear genome sequences or mRNAs of eukaryotes (Supplementary Tables S1 and S2). Interestingly, our database search revealed that the 348 bp NPR_7 corresponded primarily with sequences in animal chromosomes, notably those of bony fishes, along with a few matches with fungal and bacterial sequences. To verify the presence of these NPRs in the V. sambucifolia f. dageletiana plastome, we re-examined the occurrence of these regions using other genomic DNA extracted from a different individual of V. sambucifolia f. dageletiana deposited in our laboratory and three individuals of V. sambucifolia f. dageletiana with different voucher specimens acquired from the Plant DNA Bank in Korea (PDBK). The eight NPRs in each of the extracted DNAs were re-sequenced using the primer sets newly designed in this study. With the exception of the number of poly-G bases in NPR_1 and a 12 bp deletion in NPR_7 in one sample, all sequenced regions were found to be identical to the initially determined plastome sequences of V. sambucifolia f. dageletiana. Additionally, the average depths of the 10 NPRs were almost identical to those of their flanking sequences. Consequently, we assumed that these NPRs are true V. sambucifolia f. dageletiana plastome sequences, as opposed to copies of genomic regions originating from the nucleus or the mitogenome.
Insertions or deletions (indels) were found to be common features of the plastomes of both V. sambucifolia f. dageletiana and V. officinalis. Given the difficulty in determining whether the indels in these plastomes represent insertions or deletions, based on the available data, for the purposes of the present study, we simply defined all detected indels as deletions, compared with the sequence of the other plastome. In addition to 79 deletions caused by IR expansion/contraction at IR–LSC junctions, we detected a further 205 deletions 1–762 bp in length (Figure 3A), which we classified into four types (I to IV) (Figure 3B). Type I deletions were caused by simple sequence repeats (SSRs), whereas those in categories II and III included a single pair of repetitive sequences. If there was an intervening sequence of 1 bp or no intervening sequence between two repeats, we considered these as tandem repeats and defined the deletion of one repeating unit as type II. Type III deletions were defined as those including an intervening sequence more than 2 bp in length and a repeating unit. The remaining deletions that were unassociated with repetitive sequences were classified as type IV. We established that a total of 35 type I deletions were caused by length variations in homopolymer repeats of less than 7 bp. Type II deletions were found to occur more frequently within the plastome of V. sambucifolia f. dageletiana than in that of V. officinalis, with those longer than 30 bp being detected exclusively in the V. sambucifolia f. dageletiana plastome. Similarly, type III deletions were detected solely in the plastome of V. sambucifolia f. dageletiana. In contrast to type II and III deletions associated with repetitive sequences and slipped strand mispairing [45], those in the type IV category were detected mainly within the V. officinalis plastome.
We also identified 246 SSRs in the plastome of V. sambucifolia f. dageletiana, among which mono-nucleotide repeats accounted for the overwhelming majority of 95.2% (85% for A/T and 10.2% for C/G) (Supplementary Table S3). Although di-, tetra-, and hexa-nucleotide repeats were also identified, we were unable to detect any tri- or penta-nucleotide repeats. Consistent with these observations, at the genus level, mono-nucleotide repeats appear to be particularly prevalent throughout Caprifoliaceae (Supplementary Figure S1). Interestingly, however, whereas penta-nucleotide repeats were not found in any genera, hexa-nucleotide repeats were identified in all genera within the family.

2.2. Features of the V. sambucifolia f. dageletiana Mitochondrial Genome

The mitogenome of V. sambucifolia f. dageletiana was found to be 1,187,459 bp in length, with an average coverage of 199.4 ± 32.7, a GC content of 42.5%, and 74 genes encoding 32 protein-coding genes, 5 ribosomal RNAs, 32 transfer RNAs, and 5 pseudogenes. Compared with the mitogenomes of Asterales species, the following elements, which play roles in mitochondrial electron transport and oxidative phosphorylation, remained intact: NADH dehydrogenase subunits in complex I, apocytochrome in complex III, cytochrome c oxidase subunits in complex IV, ATP synthase subunits in complex V, and ccm genes in cytochrome c biogenesis system I (Table 1). However, certain genes encoding ribosomal proteins (rpl16, rps1, rps13, and rps19) were identified as pseudogenes, owing to the presence of internal stop codons. Among the 32 tRNAs, 14 are similar to the plastid tRNAs of V. sambucifolia f. dageletiana, namely trnD-GUC, trnF-GAA, trnH-GUG (2), trnM-CAU (2), trnN-GUU, trnP-UGG (2), trnS-GGA, trnT-UGU, trnV-GAC, and trnW-CCA (2).
Among the 1,187,459 bp of the mitogenome, 1313 regions encompassing 749,000 bp were found to correspond to sequences in the mitogenomes of 364 terrestrial plants. However, we were unable to detect any sequences in these plants corresponding to a further 1015 regions with lengths in excess of 30 bp. Excluding the 39 regions duplicated within these 1015 non-mitogenomic regions, the remaining 976 regions cover 418,671 bp and range in length from 30 to 5965 bp, and we found that 83.4% of the sequences showed no correspondence with sequences in the searched nucleotide database, based on the search options 11-word size and 1e-3 expected value (e-value). To determine whether a substantial proportion of large-sized mitogenomes is unique, and if so, whether these unique genomic sequences are a consequence of nucleus-to-mitochondrion gene transfer, we examined five mitogenomes of similar lengths to the Valeriana sambucifolia f. dageletiana mitogenome, characterized by corresponding sequences, three of which were found to show correspondence to whole nuclear genome sequences (Table 2). BLAST results revealed that the proportion of regions showing no matches with sequences in the nucleotide database ranged from 12.0% (Acer yangbiense) to 63.5% (Platycodon grandiflorus), whereas portions of mitogenomes ranging from 1.1% (Acer yangbiense) to 25.0% (Cucurbita pepo subsp. pepo) were identified in the corresponding nuclear genomes.
We established that dispersed repeat regions account for 17.1% of the total Valeriana sambucifolia f. dageletiana mitogenome and that these regions are evenly distributed throughout the genome (Figure 4A). Moreover, we found that for repeats less than 100 bp in length, the sequence length was inversely proportional to the percentage pairwise identity between dispersed repeats. In contrast, however, we detected no significant correlation between length and pairwise identity among the repeats exceeding 200 bp in length (Figure 4B). Furthermore, compared with repeats of Asterales, the lengths of identical repeats were found to be of medium size (Supplementary Figure S2).
We also identified 1205 SSRs in the mitogenome of V. sambucifolia f. dageletiana, 901 of which are mono-nucleotide repeats, comprising 74.8% of the total (Supplementary Table S4). With the exception of tetra-nucleotide repeats, the occurrence of which was found to be more frequent than that of di-, tri-, penta-, or hexa-nucleotide repeats, the length of the motif was inversely correlated with the frequency of repeat.

2.3. Intercompartmental Gene Transfer in the Mitogenome of V. sambucifolia f. dageletiana

Although phylogenetic analysis would ostensibly appear an effective approach for investigating the direction of intercompartmental gene transfer, certain mutations, such as those causing rearrangements, deletions, and contig degradation, tend to prevent alignment between transferred DNAs and the sequences of origin. Compared with mitogenomes, we assumed that the acquisitions of foreign organellar DNA would be relatively uncommon in the plastomes of terrestrial plants [38,41,46,47], even though early terrestrial plant lineages, such as mosses, appear to have relatively conserved mitogenomes [48]. On the basis of this assumption, we reasoned that if DNAs prevalent in the plastomes of these plants are also found in the mitogenomes, then they should be considered to have originated from the plastomes and not vice versa.
In the present study, we detected 162 regions in the mitogenome of V. sambucifolia f. dageletiana that either partially or entirely correspond to sequences in terrestrial plant plastomes, including 13 duplicated regions. To determine the direction of gene transfer between the two organellar compartments, based on how data were categorized by certain minimum thresholds, we used three thresholds (40%, 10%, and 1%), depending on how many two-organellar genomes were hit against queries. Among 149 regions, we found that 68 corresponded to sequences in more than 40% of 9,025 plastomes and in more than 1% of a further 364 mitogenomes (Table 3 and Figure 5). These regions were accordingly considered mitochondrial plastid DNAs (MTPTs) that may have been translocated from the plastome to the mitogenome. Conversely, we identified 46 regions corresponding to sequences in at least 150 mitogenomes but to sequences in less than 1% of plastomes (Table 3 and Figure 5). These regions were thus considered to be plastid mitochondrial DNAs (PTMTs), which may have undergone translocation from mitogenomes to plastomes. However, the complexity of the sequences of other regions precluded an assumption of their origin. Only four regions of the seven that showed comparatively low correspondence to sequences in organellar genomes were similar to sequences in certain plastomes of Orobanchaceae or Caprifoliaceae.

3. Discussion

3.1. Frequent Gene Transfers from the V. sambucifolia f. dageletiana mitogenome

Although the mitogenomes of early diverged terrestrial plants, such as mosses, tend to be stable with respect to size and structure [48], these genomes have undergone comparatively rapid increases in size and structural diversity in the more evolutionarily advanced vascular plants [49,50,51]. The notable differences in the lengths of vascular plant mitogenomes are believed to be attributable to the transfer of genes from the co-occurring nucleus and chloroplasts [16,39,52] or that of foreign DNA derived from other plants [30,53] and duplications within the mitogenome [54,55,56].
With the exception of the mitogenome of Platycodon grandifloras, the length of the V. sambucifolia f. dageletiana mitogenome is three to four times that of the mitogenomes in other plants within Asterales, although without any extensive variation in gene content. In this regard, although dispersed repeat sequences account for 17.1% of the V. sambucifolia f. dageletiana mitogenome, the primary factor contributing to an expansion in genome size, accounting for 31.6% of the entire mitogenome, is the presence of ostensibly unique regions comprising sequences for which we obtained no matches with sequences in the mitogenomes of other vascular plants reported to date. However, far from being unique to V. sambucifolia f. dageletiana, the presence of substantial unmatched portions of the mitogenome appears to be a common phenomenon associated with large-sized plant mitogenomes. Interestingly, however, some of these unmatched sequences were found to correspond to sequences in the nuclear genome of the same organism (Table 2).
On the basis of our understanding of the monophyly of terrestrial plants, it might be assumed that there would not have been high mitogenome sequence diversity among the populations of ancestral plants. As these plants evolved, however, there would have been increasing opportunities for bidirectional intracellular gene transfer between mitogenomes and nuclear genomes [32,46,57], along with the degradation and deletion of foreign DNAs in both nuclear and mitochondrial genomes, during the diversification of plants [16,58] (Figure 6). Accordingly, these events may have led to alterations in mitochondrial-incorporated nuclear DNA such that these sequences diverged from those of the original nuclear DNA, thereby obscuring identification of the origins of mitogenome sequences in terrestrial plants. In addition, the horizontal transfer of genes from other plant genomes to mitogenomes may have contributed to greater sequence and structural complexity [53]. Accordingly, these series of evolutionary events, involving gene transfer and decay, can be assumed to have led to a marked diversification in mitogenome sequences among the different orders, families, and genera of terrestrial plants, with the presence of certain modified sequences becoming intrinsic features of the mitogenomes of these plants.
Unfortunately, apart from the mitogenome of V. sambucifolia f. dageletiana sequenced in the present study, there is currently no information available regarding the sequences of mitogenomes in the order Dipsacales, and to date, there have been no reports of nuclear genome sequences in species of the genus Valeriana. Consequently, at present, we are unable to establish the identity of those sequences transferred directly to the mitogenome of V. sambucifolia f. dageletiana. Nonetheless, it is clear that there has at least been a frequent transfer of genes into mitogenomes within the Dipsacales lineage, and further study of the nuclear genome should accordingly be undertaken to gain a more in-depth insight into the intrinsic portions of incorporated mitogenomic DNA at the species level.

3.2. Origins of Non-Plastomic Regions Found in the V. sambucifolia f. dageletiana Plastome

Although in theory, six types of reciprocal inter-compartmental gene transfer among the three plant genomes are possible, the occurrence of gene flow from either the nucleus or mitochondria to plastids has been considered to be extremely infrequent [46,47], given the structural compactness of plastomes, limitation of non-homologous recombination, lack of an efficient uptake system for exogenous DNA, and the improbability of inter-chloroplast fusion [47]. However, accumulating genomic evidence is providing new insights into the intracellular transfer of genes from the mitogenome to the plastome in a number of different plant lineages: Apiales [40,59], Vitales [39], Spindales [42], Poales [60], and Gentianales [41,61]. Moreover, the horizontal transfer of genes from the genomes of distinctly related species, mainly bacteria, to the plastomes of non-plant organisms [62,63,64,65] and ferns [66,67,68] has also been reported.
In the present study, we identified eight NPRs, including two NPRs duplicated in the IR regions. Among the six excluding the duplicated ones, there was found no correspondence to any nucleotide sequences reported to date (Supplementary Table S1). As mentioned in the previous section, an additional 1.15% to 25% of mitogenomic regions in certain plants is conjectured to be derived from nuclear genome sequences, although, at present, this cannot be verified owing to the current deficiency in sequence data. Thus, subject to further confirmation, we speculated that unidentified sequences in the V. sambucifolia f. dageletiana plastome are nuclear derived and have been incorporated into the plastome in the same manner whereby other nuclear DNA sequences have integrated into plastid genomes. However, it remains conceivable that these unidentified sequences could have originated directly from bacteria or fungi. Indeed, in this regard, it has been established that the plastomes of Mankyua chejuense, an early diverging species of eusporangiate fern, contain an approximate 8 kb insertion between rps4 and trnL-UAA and that certain open reading frames in this insertion are likely to have been acquired from distantly related taxa such as algae or bacteria [66]. Moreover, although these open reading frames have also been identified in other fern organelles [67,68], most of the intergenic spacers between open reading frames detected in M. chejuense do not match any sequence regions in current nucleotide databases. These findings thus tend to imply that the identity of horizontally transferred non-coding sequences would not be ascertained based on nucleotide database searches. Accordingly, it can be anticipated that the accumulation of larger amounts of sequence data, including those of bacteria, fungi, and the nuclear DNA of terrestrial plants, will, in the future, enable us to gain a more comprehensive insight into the origins of these unidentified NPRs.
Nevertheless, the source of the foreign DNA sequence NPR_7 detected in the V. sambucifolia f. dageletiana plastome remains something of a conundrum. Although sequences similar to NPR_7 have infrequently been identified in some lineages of bacteria and fungi, this sequence was found to correspond most closely to certain regions in the DNA of bony fishes. Generally, gene transfer is mediated via physical contact, notably through inter-compartmental gene transfer within the same cell or horizontal gene transfer between endosymbionts and hosts or two physically linked species. In contrast, given their spatial isolation, it is difficult to envisage any similar mode of contact between bony fishes and V. sambucifolia f. dageletiana. On the basis of the terrestrial plant plastome database and BLAST results, we assume that the NPR_7 sequence has recently undergone to be transferred into the V. sambucifolia plastome during the evolution and that the ancestor of bony fishes contained an NPR_7-like region in its nuclear genome.
As mechanisms to account for the similarity of sequences identified in the V. sambucifolia f. dageletiana plastome and bony fish, we propose three conceivable scenarios. The first of these envisages vertical inheritance from common ancestors of eukaryotes and bacteria (Figure 7A). Although most of the BLAST-identified sequences corresponding to NPR_7 are derived from bony fishes, we found that the NPR_7 sequence also rarely showed a certain degree of match to the messenger RNAs and DNAs of mosquitos, primates, apicomplexans, crustaceans, bivalves, and bacteria. It is thus conceivable that regions of NPR_7-like ancient DNA did not need to be maintained for survival. Accordingly, they have been less influenced by purifying selection during evolution and thus more readily degraded or deleted than coding sequences and sequence motifs. Moreover, whereas some of these non-essential DNA sequences harbored in the nuclear genomes of eukaryotes may have disappeared within certain lineages, others may have been retained, as in V. sambucifolia f. dageletiana. In the final event in this scenario, it is probable that the NPR_7 sequence underwent translocation from the nucleus to the plastid in a manner similar to that of other plastid components derived from nuclear DNA.
In the second scenario, we posit that the NPR_7 region of V. sambucifolia f. dageletiana and NPR_7-like regions of bony fish originated from other organisms, such as bacteria (Figure 7B). The horizontal transfer of genes from bacteria to plants [69,70] and animals [71,72,73] is well documented, and the acquisition of genes from bacteria has similarly been detected in protists [74,75] and fungi [76]. Accordingly, the transmission of bacterial DNA into eukaryotes is probably a universal phenomenon. Consequently, the appearance of NPR_7-like sequences in diverse eukaryote lineages might be the result of independent horizontal gene transfer from bacteria to different lineages.
In the third scenario, we speculate that NPR_7 was transferred from bony fishes to an ancestor of V. sambucifolia f. dageletiana via bacteria, fungi, or other vectors (Figure 7C). Although the transfer of genes between animals and terrestrial plants is considered to be an extremely rare evolutionary phenomenon, its likelihood should not be completely dismissed. Repetitive DNA sequences, referred to as transposable elements or transposons, have a tendency to relocate within genomes [77] and even surmount species boundaries [78,79]. For example, Lin et al. [80] revealed that the prominent Penelope-like elements, a group of arthropod retrotransposons, were transferred to a common ancestor of conifers, either directly or via vectors such as bacteria and fungi.
Nevertheless, despite the plausibility of these three scenarios as explanations for the origin of the anomalous NPR_7 sequence, we still require more persuasive experimental evidence to fully comprehend this finding. For example, given the current limitations regarding database resources, we have been unable to establish whether NPR_7-like sequences are found in the genomes of a wider range of terrestrial plants or whether, with the exception of certain taxa such as V. sambucifolia f. dageletiana, most of these sequences have been eliminated during the course of evolution, as envisaged in the vertical inheritance scenario (Figure 7A). Furthermore, if NPR_7-like sequences are of bacterial origin, it is worth exploring as to why these sequences were identified in the genomes of only a few bacteria (Figure 7B) or why, if the sequence was derived from bony fishes via vector mediation, horizontal gene transfer has seemingly occurred only in an endemic plant isolated on a small island in the nearly 250 million years ago since the common ancestor of bony fishes diverged [81] (Figure 7C). Although, given the lack of relevant data, we are currently unable to resolve these questions with any degree of clarity, a maximum-parsimony scenario would tend to suggest that DNA derived from bony fish was transferred to V. sambucifolia f. dageletiana or its most recent common ancestor via intermediary vectors such as bacteria or fungi, whereas independent gene transfer from bony fish may have occurred in other, more distantly related taxa.

3.3. Rapid Structural Mutation of Plastomes in the Genus Valeriana

With the exception of certain lineages, such as those of Campanulaceae [82], Geraniaceae [83,84], Oleaceae [85], Orchidaceae [86], and Orobanchaceae [26], the plastomes of flowering plants tend to be structurally conserved [17]. However, although the structures of plastomes in plants of the family Caprifoliaceae have appeared to be highly conserved and characterized by the elimination of rpl2, rps19, and ycf1 from the typical IR [44], the recently sequenced plastomes of two Dipsacus [87] and a Linnaeae [88] species in this family have revealed retention of the typical IRs of angiosperms. In addition to contracted IR regions with respect to the rpl2, rps19, and ycf1 loci, the boundaries of IR regions in the V. sambucifolia f. dageletiana plastome appear to have undergone a more complex series of shifts. In Caprifoliaceae, transcription of the trnH gene appears to be modified to a certain extent, owing to gene duplication in the LSC region throughout the family, including V. officinalis [44]. Contrastingly, in the V. sambucifolia f. dageletiana plastome, duplication of the trnH gene shows a different pattern, involving complete translocation into the IR region, similar to that seen in certain Patrinia species [43,44]. Furthermore, it is assumed that two expansions and one contraction of the IR have occurred in the plastomes of V. sambucifolia f. dageletiana and V. officinalis, at least in the IR/SSC border regions (Figure 2). Interestingly, we found that a majority of the NPRs detected in the V. sambucifolia f. dageletiana plastome are located in the vicinity of the trnN-GUU locus or within ycf1, which are typically close to each other at adjacent IR/SSC borders in angiosperms. In this regard, the mutation of transfer RNAs has been well documented. For example, the inversion of plastid DNA between distinct tRNA genes has been identified in cereals [89] and extensive rearrangements associated with tRNA genes in the plastome of Trachelium caeruleum have been reported [82]. Moreover, it has been found that tRNA genes in the plastomes of ferns appear to contribute to genome instability, leading to the insertion of mobile elements [68].
Although the association between IR/SSC border instability and the translocation of foreign DNA is ambiguous, given that IR/SSC borders in the plastomes of Caprifoliaceae plants tend to be highly conserved [44] and our finding that NPRs are generally located in the vicinity of tRNA genes, we assume that in the V. sambucifolia f. dageletiana plastome, foreign DNA has become inserted in sequences adjacent to the IR/SSC borders. Such insertion events are presumed to weaken the stability of IR/SSC borders, thereby providing conditions conducive to IR expansions/contractions. Consequently, we speculate that gene duplication, horizontal gene transfer, and IR expansions/contraction have occurred independently and recently in Valeriana, subsequent to the divergence of V. sambucifolia and V. officinalis.

3.4. Genetic Consequence of Anagenetic Speciation on Ulleung Island

Among the species endemic to Ulleung Island, the plastome sequences of those reported to date have been found to be highly conserved with respect to both structure and gene content [90,91,92,93,94,95]. The plastome sequence of V. sambucifolia f. dageletiana reported herein thus appears to be a notable exception. Although Yang et al. [91] proposed the use of four highly variable regions to assess the genetic consequences of anagenetic speciation on Ulleung Island, these four loci are generally highly variable throughout angiosperms [23,96,97,98,99]. Therefore, there is currently no conclusive evidence supporting anagenetic speciation on Ulleung Island based on organelle genome evolution.
In contrast to the plastome data obtained for other plant taxa distributed on Ulleung Island, we identified pronounced structural changes in the organellar genomes of V. sambucifolia f. dageletiana, including rearrangements and gene transfer. This, thus, raises the question as to why such structural variants have been detected only in the plastid genome of V. sambucifolia f. dageletiana and why the plastomes of 10 other species endemic to the island have a typical structure without any notable alterations. To resolve these issues, we recommend extensive comparative genomic analyses, focusing not only on plastomes but also on other genomes (mitogenomes and/or nuclear genomes) among island endemics and closely related mainland taxa.

4. Materials and Methods

4.1. Sampling, DNA Extraction, and Next-Generation Sequencing

Specimens of V. sambucifolia f. dageletiana (Figure 1) were collected from the plant’s natural habitat on Ulleung Island, Korea, and a voucher specimen has been deposited at Kyungpook National University. For the purposes of DNA extraction, we desiccated two leaves using silica gel and then ground these to a powder using a homogenizer. Total genomic DNA was extracted using a DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) and was quantified using a NanoDrop spectrophotometer. Using this quantified DNA, three different libraries for next-generation sequencing (two paired-end libraries and one mate-pair library (5 kb paired distance)) were generated using the services of Macrogen, Seoul, Korea, and sequenced using the Illumina platform (Supplementary Table S5). Raw reads were trimmed using Trimmomatic 0.39 [100], with the options LEADING:10, TRAILING:10, SLIDINGWINDOWS:4:20, and MINLEN:50.

4.2. Assembly and Annotation of the Plastome and Mitogenome

The plastome of V. sambucifolia f. dageletiana was assembled and annotated using the protocol described by Kim and Chase [101], with the exception of read trimming. For the mitogenome, we used slightly modified versions of previously developed baiting and iteration strategies [16,102]. Having trimmed the raw reads, these were initially assembled using the de novo assembly programs Abyss [103,104] and Megahit [105,106] (Supplementary Figure S3). Plastome contigs and contigs less than 10 kb in size were then filtered out, and contigs including mitochondrial fragments were extracted using BLASTn [107]. Thereafter, the reads were mapped to the contigs using the Geneious assembler [108] with zero mismatches and zero gaps for 25 iterations. The consensus contigs thus generated were aligned using de novo assembly, and the resultant contigs were repeatedly re-used as references until all contigs were incorporated into a single circular contig.

4.3. Analysis of Repetitive Sequences in Organellar Genomes

Simple sequence repeats (SSRs) detected in both organellar genomes were identified using MISA-web, with the options 7, 5, 4, 3, 3, and 3 for mono- to hexa-nucleotides, respectively. In addition, dispersed repeats in the mitogenome were identified using BLASTn [107], with a word size of 7 and an expected value (e-value) of 1 × 10−6. To minimize overestimation of repeat values, BLAST hits with identical query start and end positions were removed.

4.4. Identification of Mitochondrial Plastid DNAs in the Mitogenome of V. sambucifolia f. dageletiana

To investigate the occurrence MTPTs in the mitogenome of V. sambucifolia f. dageletiana, we blasted the mitogenome against 9025 plastome sequences of terrestrial plant origin (downloaded from the NCBI database) using BLASTn with an e-value of 1e-5 and a word size of 11. When queries corresponding to sequences in the database overlapped on the mitogenome, they were annotated using a single MTPT on the mitogenome of V. sambucifolia f. dageletiana. In addition, using BLASTn with an e-value of 1e-5 and a word size of 11, the mitogenome was blasted against 364 terrestrial plant mitogenome sequences downloaded from the NCBI database to detect non-mitogenomic regions. The different e-value thresholds applied to the aforementioned BLASTn searches were based on different total database lengths, which directly influenced the e-value.

4.5. Identification of non-plastome-originated DNAs in the plastome of V. sambucifolia f. dageletiana

The plastome of V. sambucifolia f. dageletiana was blasted against Viridiplantae plastid DNAs downloaded from GenBank (1,419,429 sequences) using BLASTn with an e-value of 1e−3 and a word size of 11. The intervening sequences that showed no match with those in the database were extracted and to identify their origin were blasted against the NCBI nucleotide database (https://blast.ncbi.nlm.nih.gov/Blast.cgi accessed on 5 April 2021) using BLASTn with an e-value of 1e−3 and a word size of 11.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms22073770/s1: Figure S1: Simple sequence repeats at the genus level in Caprifoliaceae, Figure S2: The identical repeats in mitogenome of V. sambucifolia f. dageletiana compared with those of Asterales, Figure S3: Assembly strategy for the mitogenome of V. sambucifolia f. dageletian, Table S1: BLAST results for non-plastomic regions of the V. sambucifolia f. dageletiana plastome, Table S2: Eight NPR sequences., Table S3: Simple sequence repeats in the plastome of V. sambucifolia f. dageletiana, Table S4: Simple sequence repeats in the mitogenome of V. sambucifolia f. dageletiana, Table S5: Information of next generation sequencing data.

Author Contributions

Conceptualization, H.T.K. and J.S.K.; methodology, H.T.K.; software and data curation, H.T.K.; resources, J.S.K.; writing—original draft preparation, H.T.K.; writing—review and editing, H.T.K. and J.S.K.; visualization, H.T.K.; supervision, J.S.K.; project administration, H.T.K.; and funding acquisition, J.S.K. All authors have read and agreed to the publication of the final version of this manuscript.

Funding

This research was funded by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1I1A2053517).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the Dryad Digital Repository, https://datadryad.org/stash/share/WB4rygpAm0IuoBDTuIgcMWaoH7I03UNPGy8jvIMKTWc (accessed on 5 April 2021).

Acknowledgments

The authors would like to thank W. Lee, who assisted in collecting the plant material for initial investigation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xiang, C.L.; Dong, H.J.; Landrein, S.; Zhao, F.; Yu, W.B.; Soltis, D.E.; Soltis, P.S.; Backlund, A.; Wang, H.F.; Li, D.Z.; et al. Revisiting the phylogeny of Dipsacales: New insights from phylogenomic analyses of complete plastomic sequences. J. Syst. Evol. 2020, 58, 103–117. [Google Scholar] [CrossRef]
  2. Jacobs, B.; Bell, C.; Smets, E. Fruits and seeds of the valeriana clade (Dipsacales): Diversity and evolution. Int. J. Plant Sci. 2010, 171, 421–434. [Google Scholar] [CrossRef]
  3. Bell, C.D.; Donoghue, M.J. Phylogeny and biogeography of Valerianaceae (Dipsacales) with special reference to the South American valerians. Org. Divers. Evol. 2005, 5, 147–159. [Google Scholar] [CrossRef]
  4. Backlund, A.; Moritz, T. Phylogenetic implications of an expanded valepotriate distribution in the Valerianaceae. Biochem. Syst. Ecol. 1998, 26, 309–335. [Google Scholar] [CrossRef]
  5. eFloras. 2008. Available online: http://www.efloras.org (accessed on 10 April 2020).
  6. Bjørn, G.K.; Olsson, K.; Bladh, K.W. Valeriana officinalis L.(Common valerian). Spice-and Medicinal Plants in the Nordic and Baltic Countries. 2006. Available online: https://www.nordgen.org/ngdoc/plants/publications/SPIMED_report_maj_2006.pdf (accessed on 5 April 2020).
  7. Flora of Korea Editorial Committe. The Genera of Vascular Plants of Korea; Academy Publishing, Co.: Seoul, Korea, 2007. [Google Scholar]
  8. Wang, Z.; Hwang, S.H.; Kim, J.H.; Lim, S.S. Anti-obesity effect of the above-ground part of valeriana dageletiana nakai ex f. maek extract in high-fat diet-induced obese C57BL/6N Mice. Nutrients 2017, 9, 689. [Google Scholar] [CrossRef] [PubMed]
  9. Kwon, C.-Y.; Suh, H.-W.; Choi, E.-J.; Chung, S.-Y.; Kim, J.-W. Complementary and alternative medicine in clinical practice guideline for insomnia. Korean Soc. Orient. Neuropsychiatry 2016, 27, 235–248. [Google Scholar] [CrossRef]
  10. Ryu, K.-S. Monoterpenoid of Korean valerian roots. Korean J. Pharmacogn. 1974, 5, 1–6. [Google Scholar]
  11. Federov, A.A.; Bobrov, A. Flora of Russia, the European Part and Bordering Regions; AA Balkema Publishers: Rotterdam, The Netherlands, 2001. [Google Scholar]
  12. Skalińska, M. Meiosis in a polyhaploid twin plant and a hexaploid hybrid of Valeriana sambucifolia Mikan. Acta Soc. Bot. Poloniae 1954, 23, 359–374. [Google Scholar] [CrossRef]
  13. Lövkvist, B.R. Chromosome numbers in south Swedish vascular plants. Opera Bot. 1999, 137, 1–42. [Google Scholar]
  14. Lee, Y. Chromosome numbers of flowering plants in Korea (2). J. Korean Res. Inst. Better Living 1969, 2, 141–145. [Google Scholar]
  15. Bock, R.; Knoop, V. Genomics of Chloroplasts and Mitochondria; Springer Science & Business Media: Berlin/Heidelberg, Germany; Amsterdam, The Netherlands, 2012; Volume 35. [Google Scholar]
  16. Kim, H.T.; Lee, J.M. Organellar genome analysis reveals endosymbiotic gene transfers in tomato. PLoS ONE 2018, 13, e0202279. [Google Scholar] [CrossRef]
  17. Ruhlman, T.A.; Jansen, R.K. The plastid genomes of flowering plants. Methods Mol. Biol. 2014, 1132, 3–38. [Google Scholar] [CrossRef]
  18. Wolf, P.G.; Roper, J.M.; Duffy, A.M. The evolution of chloroplast genome structure in ferns. Genome 2010, 53, 731–738. [Google Scholar] [CrossRef]
  19. Wolf, P.G.; Karol, K.G.; Mandoli, D.F.; Kuehl, J.; Arumuganathan, K.; Ellis, M.W.; Mishler, B.D.; Kelch, D.G.; Olmstead, R.G.; Boore, J.L. The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae). Gene 2005, 350, 117–128. [Google Scholar] [CrossRef] [PubMed]
  20. De Pamphilis, C.W.; Palmer, J.D. Loss of photosynthetic and chlororespiratory genes from the plastid genome of a parasitic flowering plant. Nature 1990, 348, 337–339. [Google Scholar] [CrossRef] [PubMed]
  21. Barrett, C.F.; Freudenstein, J.V.; Li, J.; Mayfield-Jones, D.R.; Perez, L.; Pires, J.C.; Santos, C. Investigating the path of plastid genome degradation in an early-transitional clade of heterotrophic orchids, and implications for heterotrophic angiosperms. Mol. Biol. Evol. 2014, 31, 3095–3112. [Google Scholar] [CrossRef] [PubMed]
  22. Givnish, T.J.; Spalink, D.; Ames, M.; Lyon, S.P.; Hunter, S.J.; Zuluaga, A.; Iles, W.J.; Clements, M.A.; Arroyo, M.T.; Leebens-Mack, J.; et al. Orchid phylogenomics and multiple drivers of their extraordinary diversification. Proc. Biol. Sci. 2015, 282, 171–180. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, J.B.; Tang, M.; Li, H.T.; Zhang, Z.R.; Li, D.Z. Complete chloroplast genome of the genus Cymbidium: Lights into the species identification, phylogenetic implications and population genetic analyses. BMC Evol. Biol. 2013, 13, 84. [Google Scholar] [CrossRef]
  24. Kim, H.T.; Shin, C.-H.; Sun, H.; Kim, J.-H. Sequencing of the plastome in the leafless green mycoheterotroph Cymbidium macrorhizon helps us to understand an early stage of fully mycoheterotrophic plastome structure. Plant Syst. Evol. 2017, 304, 1–14. [Google Scholar] [CrossRef]
  25. Barrett, C.F.; Davis, J.I. The plastid genome of the mycoheterotrophic Corallorhiza striata (Orchidaceae) is in the relatively early stages of degradation. Am. J. Bot. 2012, 99, 1513–1523. [Google Scholar] [CrossRef]
  26. Wicke, S.; Muller, K.F.; de Pamphilis, C.W.; Quandt, D.; Wickett, N.J.; Zhang, Y.; Renner, S.S.; Schneeweiss, G.M. Mechanisms of functional and physical genome reduction in photosynthetic and nonphotosynthetic parasitic plants of the broomrape family. Plant Cell 2013, 25, 3711–3725. [Google Scholar] [CrossRef] [PubMed]
  27. Karol, K.G.; Arumuganathan, K.; Boore, J.L.; Duffy, A.M.; Everett, K.D.; Hall, J.D.; Hansen, S.K.; Kuehl, J.V.; Mandoli, D.F.; Mishler, B.D.; et al. Complete plastome sequences of Equisetum arvense and Isoetes flaccida: Implications for phylogeny and plastid genome evolution of early land plant lineages. BMC Evol. Biol. 2010, 10, 321. [Google Scholar] [CrossRef]
  28. Raubeson, L.A.; Jansen, R.K. Chloroplast DNA evidence on the ancient evolutionary split in vascular land plants. Science 1992, 255, 1697–1699. [Google Scholar] [CrossRef] [PubMed]
  29. Palmer, J.D.; Herbon, L.A. Plant mitochondrial DNA evolved rapidly in structure, but slowly in sequence. J. Mol. Evol. 1988, 28, 87–97. [Google Scholar] [CrossRef]
  30. Bergthorsson, U.; Adams, K.L.; Thomason, B.; Palmer, J.D. Widespread horizontal transfer of mitochondrial genes in flowering plants. Nature 2003, 424, 197–201. [Google Scholar] [CrossRef]
  31. Wang, D.; Wu, Y.W.; Shih, A.C.; Wu, C.S.; Wang, Y.N.; Chaw, S.M. Transfer of chloroplast genomic DNA to mitochondrial genome occurred at least 300 MYA. Mol. Biol. Evol. 2007, 24, 2040–2048. [Google Scholar] [CrossRef] [PubMed]
  32. Rodriguez-Moreno, L.; Gonzalez, V.M.; Benjak, A.; Marti, M.C.; Puigdomenech, P.; Aranda, M.A.; Garcia-Mas, J. Determination of the melon chloroplast and mitochondrial genome sequences reveals that the largest reported mitochondrial genome in plants contains a significant amount of DNA having a nuclear origin. BMC Genom. 2011, 12, 424. [Google Scholar] [CrossRef] [PubMed]
  33. Naito, K.; Kaga, A.; Tomooka, N.; Kawase, M. De novo assembly of the complete organelle genome sequences of azuki bean (Vigna angularis) using next-generation sequencers. Breed. Sci. 2013, 63, 176–182. [Google Scholar] [CrossRef]
  34. Jo, Y.D.; Choi, Y.; Kim, D.H.; Kim, B.D.; Kang, B.C. Extensive structural variations between mitochondrial genomes of CMS and normal peppers (Capsicum annuum L.) revealed by complete nucleotide sequencing. BMC Genom. 2014, 15, 561. [Google Scholar] [CrossRef]
  35. Chang, S.; Wang, Y.; Lu, J.; Gai, J.; Li, J.; Chu, P.; Guan, R.; Zhao, T. The mitochondrial genome of soybean reveals complex genome structures and gene evolution at intercellular and phylogenetic levels. PLoS ONE 2013, 8, e56502. [Google Scholar] [CrossRef]
  36. Ellis, J. Promiscuous DNA-chloroplast genes inside plant mitochondria. Nature 1982, 299, 678–679. [Google Scholar] [CrossRef]
  37. Stern, D.B.; Lonsdale, D.M. Mitochondrial and chloroplast genomes of maize have a 12-kilobase DNA sequence in common. Nature 1982, 299, 698–702. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, X.C.; Chen, H.; Yang, D.; Liu, C. Diversity of mitochondrial plastid DNAs (MTPTs) in seed plants. Mitochondrial DNA A DNA Mapp. Seq. Anal. 2018, 29, 635–642. [Google Scholar] [CrossRef] [PubMed]
  39. Goremykin, V.V.; Salamini, F.; Velasco, R.; Viola, R. Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer. Mol. Biol. Evol. 2009, 26, 99–110. [Google Scholar] [CrossRef] [PubMed]
  40. Spooner, D.M.; Ruess, H.; Iorizzo, M.; Senalik, D.; Simon, P. Entire plastid phylogeny of the carrot genus (Daucus, Apiaceae): Concordance with nuclear data and mitochondrial and nuclear DNA insertions to the plastid. Am. J. Bot. 2017, 104, 296–312. [Google Scholar] [CrossRef] [PubMed]
  41. Smith, D.R. Mitochondrion-to-plastid DNA transfer: It happens. New Phytol. 2014, 202, 736–738. [Google Scholar] [CrossRef] [PubMed]
  42. Rabah, S.O.; Lee, C.; Hajrah, N.H.; Makki, R.M.; Alharby, H.F.; Alhebshi, A.M.; Sabir, J.S.M.; Jansen, R.K.; Ruhlman, T.A. Plastome sequencing of ten nonmodel crop species uncovers a large insertion of mitochondrial DNA in cashew. Plant Genome 2017, 10. [Google Scholar] [CrossRef]
  43. Jung, E.-H.; Lim, C.E.; Lee, B.Y.; Hong, S.-P. Complete chloroplast genome sequence of Patrinia saniculifolia hemsl.(Disacales: Caprifoliaceae), an endemic plant in Korea. Mitochondrial DNA Part B 2018, 3, 60–61. [Google Scholar] [CrossRef]
  44. Wang, H.-X.; Liu, H.; Moore, M.J.; Landrein, S.; Liu, B.; Zhu, Z.-X.; Wang, H.-F. Plastid phylogenomic insights into the evolution of the Caprifoliaceae s.l. (Dipsacales). Mol. Phylogenet. Evol. 2020, 142, 106641. [Google Scholar] [CrossRef]
  45. Levinson, G.; Gutman, G.A. Slipped-strand mispairing: A major mechanism for DNA sequence evolution. Mol. Biol. Evol. 1987, 4, 203–221. [Google Scholar] [CrossRef]
  46. Timmis, J.N.; Ayliffe, M.A.; Huang, C.Y.; Martin, W. Endosymbiotic gene transfer: Organelle genomes forge eukaryotic chromosomes. Nat. Rev. Genet. 2004, 5, 123–135. [Google Scholar] [CrossRef] [PubMed]
  47. Kleine, T.; Maier, U.G.; Leister, D. DNA transfer from organelles to the nucleus: The idiosyncratic genetics of endosymbiosis. Annu. Rev. Plant Biol. 2009, 60, 115–138. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, Y.; Medina, R.; Goffinet, B. 350 my of mitochondrial genome stasis in mosses, an early land plant lineage. Mol. Biol. Evol. 2014, 31, 2586–2591. [Google Scholar] [CrossRef]
  49. Guo, W.; Zhu, A.; Fan, W.; Mower, J.P. Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns. New Phytol. 2017, 213, 391–403. [Google Scholar] [CrossRef] [PubMed]
  50. Chaw, S.M.; Shih, A.C.; Wang, D.; Wu, Y.W.; Liu, S.M.; Chou, T.Y. 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. 2008, 25, 603–615. [Google Scholar] [CrossRef]
  51. Alverson, A.J.; Wei, X.X.; Rice, D.W.; Stern, D.B.; Barry, K.; Palmer, J.D. Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae). Mol. Biol. Evol. 2010, 27, 1436–1448. [Google Scholar] [CrossRef] [PubMed]
  52. Leister, D. Origin, evolution and genetic effects of nuclear insertions of organelle DNA. Trends Genet. 2005, 21, 655–663. [Google Scholar] [CrossRef]
  53. Rice, D.W.; Alverson, A.J.; Richardson, A.O.; Young, G.J.; Sanchez-Puerta, M.V.; Munzinger, J.; Barry, K.; Boore, J.L.; Zhang, Y.; de Pamphilis, C.W.; et al. Horizontal transfer of entire genomes via mitochondrial fusion in the angiosperm Amborella. Science 2013, 342, 1468–1473. [Google Scholar] [CrossRef] [PubMed]
  54. Dong, S.; Zhao, C.; Chen, F.; Liu, Y.; Zhang, S.; Wu, H.; Zhang, L.; Liu, Y. The complete mitochondrial genome of the early flowering plant Nymphaea colorata is highly repetitive with low recombination. BMC Genom. 2018, 19, 614. [Google Scholar] [CrossRef]
  55. Dong, S.; Chen, L.; Liu, Y.; Wang, Y.; Zhang, S.; Yang, L.; Lang, X.; Zhang, S. The draft mitochondrial genome of Magnolia biondii and mitochondrial phylogenomics of angiosperms. PLoS ONE 2020, 15, e0231020. [Google Scholar] [CrossRef]
  56. Satoh, M.; Kubo, T.; Nishizawa, S.; Estiati, A.; Itchoda, N.; Mikami, T. The cytoplasmic male-sterile type and normal type mitochondrial genomes of sugar beet share the same complement of genes of known function but differ in the content of expressed ORFs. Mol. Genet. Genom. 2004, 272, 247–256. [Google Scholar] [CrossRef]
  57. Michalovova, M.; Vyskot, B.; Kejnovsky, E. Analysis of plastid and mitochondrial DNA insertions in the nucleus (NUPTs and NUMTs) of six plant species: Size, relative age and chromosomal localization. Heredity 2013, 111, 314–320. [Google Scholar] [CrossRef] [PubMed]
  58. Richly, E.; Leister, D. NUPTs in sequenced eukaryotes and their genomic organization in relation to NUMTs. Mol. Biol. Evol. 2004, 21, 1972–1980. [Google Scholar] [CrossRef]
  59. Downie, S.R.; Jansen, R.K. A comparative analysis of whole plastid genomes from the apiales: Expansion and contraction of the inverted repeat, mitochondrial to plastid transfer of DNA, and identification of highly divergent noncoding regions. Syst. Bot. 2015, 40, 336–351. [Google Scholar] [CrossRef]
  60. Ma, P.F.; Zhang, Y.X.; Guo, Z.H.; Li, D.Z. Evidence for horizontal transfer of mitochondrial DNA to the plastid genome in a bamboo genus. Sci. Rep. 2015, 5, 11608. [Google Scholar] [CrossRef]
  61. Straub, S.C.; Cronn, R.C.; Edwards, C.; Fishbein, M.; Liston, A. Horizontal transfer of DNA from the mitochondrial to the plastid genome and its subsequent evolution in milkweeds (Apocynaceae). Genome Biol. Evol. 2013, 5, 1872–1885. [Google Scholar] [CrossRef] [PubMed]
  62. Rice, D.W.; Palmer, J.D. An exceptional horizontal gene transfer in plastids: Gene replacement by a distant bacterial paralog and evidence that haptophyte and cryptophyte plastids are sisters. BMC Biol. 2006, 4, 31. [Google Scholar] [CrossRef]
  63. Sheveleva, E.V.; Hallick, R.B. Recent horizontal intron transfer to a chloroplast genome. Nucleic Acids Res. 2004, 32, 803–810. [Google Scholar] [CrossRef]
  64. Brouard, J.S.; Otis, C.; Lemieux, C.; Turmel, M. Chloroplast DNA sequence of the green alga Oedogonium cardiacum (Chlorophyceae): Unique genome architecture, derived characters shared with the Chaetophorales and novel genes acquired through horizontal transfer. BMC Genom. 2008, 9, 290. [Google Scholar] [CrossRef]
  65. Delwiche, C.F.; Palmer, J.D. Rampant horizontal transfer and duplication of rubisco genes in eubacteria and plastids. Mol. Biol. Evol. 1996, 13, 873–882. [Google Scholar] [CrossRef]
  66. Kim, H.T.; Kim, K.J. Evolution of six novel ORFs in the plastome of Mankyua chejuense and phylogeny of eusporangiate ferns. Sci. Rep. 2018, 8, 16466. [Google Scholar] [CrossRef]
  67. Robison, T.A.; Grusz, A.L.; Wolf, P.G.; Mower, J.P.; Fauskee, B.D.; Sosa, K.; Schuettpelz, E. Mobile elements shape plastome evolution in ferns. Genome Biol. Evol. 2018, 10, 2558–2571. [Google Scholar] [CrossRef] [PubMed]
  68. Kim, H.T.; Kim, J.S. The dynamic evolution of mobile open reading frames in plastomes of Hymenophyllum Sm. and new insight on Hymenophyllum coreanum Nakai. Sci. Rep. 2020, 10, 11059. [Google Scholar] [CrossRef]
  69. Yue, J.; Hu, X.; Sun, H.; Yang, Y.; Huang, J. Widespread impact of horizontal gene transfer on plant colonization of land. Nat. Commun. 2012, 3, 1152. [Google Scholar] [CrossRef] [PubMed]
  70. Broothaerts, W.; Mitchell, H.J.; Weir, B.; Kaines, S.; Smith, L.M.A.; Yang, W.; Mayer, J.E.; Roa-Rodríguez, C.; Jefferson, R.A. Gene transfer to plants by diverse species of bacteria. Nature 2005, 433, 629–633. [Google Scholar] [CrossRef] [PubMed]
  71. Dunning Hotopp, J.C.; Clark, M.E.; Oliveira, D.C.; Foster, J.M.; Fischer, P.; Munoz Torres, M.C.; Giebel, J.D.; Kumar, N.; Ishmael, N.; Wang, S.; et al. Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science 2007, 317, 1753–1756. [Google Scholar] [CrossRef] [PubMed]
  72. Dunning Hotopp, J.C. Horizontal gene transfer between bacteria and animals. Trends Genet. 2011, 27, 157–163. [Google Scholar] [CrossRef] [PubMed]
  73. Wybouw, N.; Dermauw, W.; Tirry, L.; Stevens, C.; Grbic, M.; Feyereisen, R.; Van Leeuwen, T. A gene horizontally transferred from bacteria protects arthropods from host plant cyanide poisoning. eLife 2014, 3, e02365. [Google Scholar] [CrossRef]
  74. Ricard, G.; McEwan, N.R.; Dutilh, B.E.; Jouany, J.P.; Macheboeuf, D.; Mitsumori, M.; McIntosh, F.M.; Michalowski, T.; Nagamine, T.; Nelson, N.; et al. Horizontal gene transfer from Bacteria to rumen Ciliates indicates adaptation to their anaerobic, carbohydrates-rich environment. BMC Genom. 2006, 7, 22. [Google Scholar] [CrossRef]
  75. Bowler, C.; Allen, A.E.; Badger, J.H.; Grimwood, J.; Jabbari, K.; Kuo, A.; Maheswari, U.; Martens, C.; Maumus, F.; Otillar, R.P.; et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature 2008, 456, 239–244. [Google Scholar] [CrossRef]
  76. Richards, T.A.; Leonard, G.; Soanes, D.M.; Talbot, N.J. Gene transfer into the fungi. Fungal Biol. Rev. 2011, 25, 98–110. [Google Scholar] [CrossRef]
  77. Finnegan, D.J. Transposable elements in eukaryotes. In International Review of Cytology; Bourne, G.H., Danielli, J.F., Jeon, K.W., Eds.; Academic Press: Cambridge, MA, USA, 1985; Volume 93, pp. 281–326. [Google Scholar]
  78. Schaack, S.; Gilbert, C.; Feschotte, C. Promiscuous DNA: Horizontal transfer of transposable elements and why it matters for eukaryotic evolution. Trends Ecol. Evol. 2010, 25, 537–546. [Google Scholar] [CrossRef]
  79. El Baidouri, M.; Carpentier, M.C.; Cooke, R.; Gao, D.; Lasserre, E.; Llauro, C.; Mirouze, M.; Picault, N.; Jackson, S.A.; Panaud, O. Widespread and frequent horizontal transfers of transposable elements in plants. Genome Res. 2014, 24, 831–838. [Google Scholar] [CrossRef]
  80. Lin, X.; Faridi, N.; Casola, C. An ancient transkingdom horizontal transfer of penelope-like retroelements from arthropods to conifers. Genome Biol. Evol. 2016, 8, 1252–1266. [Google Scholar] [CrossRef]
  81. Betancur, R.R.; Wiley, E.O.; Arratia, G.; Acero, A.; Bailly, N.; Miya, M.; Lecointre, G.; Orti, G. Phylogenetic classification of bony fishes. BMC Evol. Biol. 2017, 17, 162. [Google Scholar] [CrossRef] [PubMed]
  82. Haberle, R.C.; Fourcade, H.M.; Boore, J.L.; Jansen, R.K. Extensive rearrangements in the chloroplast genome of Trachelium caeruleum are associated with repeats and tRNA genes. J. Mol. Evol. 2008, 66, 350–361. [Google Scholar] [CrossRef] [PubMed]
  83. Weng, M.L.; Blazier, J.C.; Govindu, M.; Jansen, R.K. Reconstruction of the ancestral plastid genome in Geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Mol. Biol. Evol. 2014, 31, 645–659. [Google Scholar] [CrossRef] [PubMed]
  84. Guisinger, M.M.; Kuehl, J.V.; Boore, J.L.; Jansen, R.K. Extreme reconfiguration of plastid genomes in the angiosperm family Geraniaceae: Rearrangements, repeats, and codon usage. Mol. Biol. Evol. 2011, 28, 583–600. [Google Scholar] [CrossRef]
  85. Lee, H.L.; Jansen, R.K.; Chumley, T.W.; Kim, K.J. Gene relocations within chloroplast genomes of Jasminum and Menodora (Oleaceae) are due to multiple, overlapping inversions. Mol. Biol. Evol. 2007, 24, 1161–1180. [Google Scholar] [CrossRef]
  86. Feng, Y.L.; Wicke, S.; Li, J.W.; Han, Y.; Lin, C.S.; Li, D.Z.; Zhou, T.T.; Huang, W.C.; Huang, L.Q.; Jin, X.H. Lineage-specific reductions of plastid genomes in an orchid tribe with partially and fully mycoheterotrophic species. Genome Biol. Evol. 2016, 8, 2164–2175. [Google Scholar] [CrossRef]
  87. Park, I.; Yang, S.; Kim, W.J.; Noh, P.; Lee, H.O.; Moon, B.C. Authentication of herbal medicines Dipsacus asper and Phlomoides umbrosa using DNA barcodes, chloroplast genome, and sequence characterized amplified region (SCAR) marker. Molecules 2018, 23, 1748. [Google Scholar] [CrossRef]
  88. Liu, H.; Xia, M.; Xiao, Q.; Fang, J.; Wang, A.; Chen, S.; Zhang, D. Characterization of the complete chloroplast genome of Linnaea borealis, a rare, clonal self-incompatible plant. Mitochondrial DNA Part B 2020, 5, 200–201. [Google Scholar] [CrossRef] [PubMed]
  89. Hiratsuka, J.; Shimada, H.; Whittier, R.; Ishibashi, T.; Sakamoto, M.; Mori, M.; Kondo, C.; Honji, Y.; Sun, C.-R.; Meng, B.-Y. The complete sequence of the rice (Oryza sativa) chloroplast genome: Intermolecular recombination between distinct tRNA genes accounts for a major plastid DNA inversion during the evolution of the cereals. Mol. General Genet. MGG 1989, 217, 185–194. [Google Scholar] [CrossRef]
  90. Yang, J.Y.; Pak, J.-H.; Kim, S.-C. The complete plastome sequence of Rubus takesimensis endemic to Ulleung Island, Korea: Insights into molecular evolution of anagenetically derived species in Rubus (Rosaceae). Gene 2018, 668, 221–228. [Google Scholar] [CrossRef] [PubMed]
  91. Yang, J.; Kang, G.-H.; Pak, J.-H.; Kim, S.-C. Characterization and comparison of two complete plastomes of Rosaceae species (Potentilla dickinsii var. glabrata and Spiraea insularis) endemic to Ulleung Island, Korea. Int. J. Mol. Sci. 2020, 21, 4933. [Google Scholar] [CrossRef] [PubMed]
  92. Cho, M.-S.; Yang, J.Y.; Kim, S.-C. Complete chloroplast genome of Ulleung Island endemic flowering cherry, Prunus takesimensis (Rosaceae), in Korea. Mitochondrial DNA Part B 2018, 3, 274–275. [Google Scholar] [CrossRef]
  93. Park, J.-S.; Jin, D.-P.; Park, J.-W.; Choi, B.-H. Complete chloroplast genome of Fagus multinervis, a beech species endemic to Ulleung Island in South Korea. Mitochondrial DNA Part B 2019, 4, 1698–1699. [Google Scholar] [CrossRef]
  94. Yang, J.Y.; Pak, J.-H.; Kim, S.-C. Chloroplast genome of critically endangered Cotoneaster wilsonii (Rosaceae) endemic to Ulleung Island, Korea. Mitochondrial DNA Part B 2019, 4, 3892–3893. [Google Scholar] [CrossRef] [PubMed]
  95. Gil, H.-Y.; Kim, S.-C. The plastome sequence of Ulleung Rowan, Sorbus ulleungensis (Rosaceae), a new endemic species on Ulleung Island, Korea. Mitochondrial DNA Part B 2018, 3, 284–285. [Google Scholar] [CrossRef]
  96. Shaw, J.; Lickey, E.B.; Schilling, E.E.; Small, R.L. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. Am. J. Bot. 2007, 94, 275–288. [Google Scholar] [CrossRef]
  97. Shaw, J.; Lickey, E.B.; Beck, J.T.; Farmer, S.B.; Liu, W.; Miller, J.; Siripun, K.C.; Winder, C.T.; Schilling, E.E.; Small, R.L. The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am. J. Bot. 2005, 92, 142–166. [Google Scholar] [CrossRef]
  98. Shaw, J.; Shafer, H.L.; Leonard, O.R.; Kovach, M.J.; Schorr, M.; Morris, A.B. Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: The tortoise and the hare IV. Am. J. Bot. 2014, 101, 1987–2004. [Google Scholar] [CrossRef]
  99. Kim, H.T.; Kim, J.S.; Lee, Y.M.; Mun, J.-H.; Kim, J.-H. Molecular markers for phylogenetic applications derived from comparative plastome analysis of Prunus species. J. Syst. Evol. 2019, 57, 15–22. [Google Scholar] [CrossRef]
  100. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  101. Kim, H.T.; Chase, M.W. Independent degradation in genes of the plastid ndh gene family in species of the orchid genus Cymbidium (Orchidaceae; Epidendroideae). PLoS ONE 2017, 12, e0187318. [Google Scholar] [CrossRef] [PubMed]
  102. Hahn, C.; Bachmann, L.; Chevreux, B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—A baiting and iterative mapping approach. Nucleic Acids Res. 2013, 41, e129. [Google Scholar] [CrossRef] [PubMed]
  103. Jackman, S.D.; Vandervalk, B.P.; Mohamadi, H.; Chu, J.; Yeo, S.; Hammond, S.A.; Jahesh, G.; Khan, H.; Coombe, L.; Warren, R.L. ABySS 2.0: Resource-efficient assembly of large genomes using a Bloom filter. Genome Res. 2017, 27, 768–777. [Google Scholar] [CrossRef] [PubMed]
  104. Simpson, J.T.; Wong, K.; Jackman, S.D.; Schein, J.E.; Jones, S.J.; Birol, I. ABySS: A parallel assembler for short read sequence data. Genome Res. 2009, 19, 1117–1123. [Google Scholar] [CrossRef] [PubMed]
  105. Li, D.; Luo, R.; Liu, C.-M.; Leung, C.-M.; Ting, H.-F.; Sadakane, K.; Yamashita, H.; Lam, T.-W. MEGAHIT v1. 0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods 2016, 102, 3–11. [Google Scholar] [CrossRef]
  106. Li, D.; Liu, C.M.; Luo, R.; Sadakane, K.; Lam, T.W. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015, 31, 1674–1676. [Google Scholar] [CrossRef]
  107. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  108. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Valeriana sambucifolia f. dageletiana in its natural habitat on Ulleung Island, Korea. (A) General feature of entire plant; (B) Leaves and flower buds., (C) The base of stems and runners.
Figure 1. Valeriana sambucifolia f. dageletiana in its natural habitat on Ulleung Island, Korea. (A) General feature of entire plant; (B) Leaves and flower buds., (C) The base of stems and runners.
Ijms 22 03770 g001
Figure 2. Comparison of inverted repeat (IR) expansion and contraction at the boundaries between IR and single-copy (SC) regions in the plastomes of Valeriana sambucifolia f. dageletiana and Valeriana officinalis.
Figure 2. Comparison of inverted repeat (IR) expansion and contraction at the boundaries between IR and single-copy (SC) regions in the plastomes of Valeriana sambucifolia f. dageletiana and Valeriana officinalis.
Ijms 22 03770 g002
Figure 3. Deletions detected in the plastomes of Valeriana sambucifolia f. dageletiana (red bars) and Valeriana officinalis (green bars). (A) Total deletions. (B) The different types of deletions.
Figure 3. Deletions detected in the plastomes of Valeriana sambucifolia f. dageletiana (red bars) and Valeriana officinalis (green bars). (A) Total deletions. (B) The different types of deletions.
Ijms 22 03770 g003
Figure 4. Dispersed repeats in the mitogenome of Valeriana sambucifolia f. dageletiana. Different-colored dots denote different percentage pairwise identities between dispersed repeats. (A) The length of repeats based on position. (B) Correlation between pairwise identity and the length of repeats.
Figure 4. Dispersed repeats in the mitogenome of Valeriana sambucifolia f. dageletiana. Different-colored dots denote different percentage pairwise identities between dispersed repeats. (A) The length of repeats based on position. (B) Correlation between pairwise identity and the length of repeats.
Ijms 22 03770 g004
Figure 5. BLAST results for searches against 9025 plastomes (X axis) and 364 mitogenomes (Y axis) of terrestrial plants. Different-colored dots denote the difference in the percentage matches for subjects.
Figure 5. BLAST results for searches against 9025 plastomes (X axis) and 364 mitogenomes (Y axis) of terrestrial plants. Different-colored dots denote the difference in the percentage matches for subjects.
Ijms 22 03770 g005
Figure 6. Model of the proposed origins of unidentified mitogenome sequences found in Valeriana sambucifolia f. dageletiana.
Figure 6. Model of the proposed origins of unidentified mitogenome sequences found in Valeriana sambucifolia f. dageletiana.
Ijms 22 03770 g006
Figure 7. Three scenarios proposed to explain the occurrence of sequence region similarity between the plastome of Valeriana sambucifolia f. dageletiana and the nuclear genome of bony fishes. (A) Vertical inheritance from a common ancestor of eukaryotes and bacteria. (B) Independent gene transfer from bacteria. (C) Horizontal gene transfer from bony fishes to V. sambucifolia f. dageletiana via bacteria.
Figure 7. Three scenarios proposed to explain the occurrence of sequence region similarity between the plastome of Valeriana sambucifolia f. dageletiana and the nuclear genome of bony fishes. (A) Vertical inheritance from a common ancestor of eukaryotes and bacteria. (B) Independent gene transfer from bacteria. (C) Horizontal gene transfer from bony fishes to V. sambucifolia f. dageletiana via bacteria.
Ijms 22 03770 g007
Table 1. Mitogenome-coding genes in V. sambucifolia f. dageletiana compared with those in species in the order Asterales. The shaded rows indicate conserved genes in the orders Asterles and Dipsacales.
Table 1. Mitogenome-coding genes in V. sambucifolia f. dageletiana compared with those in species in the order Asterales. The shaded rows indicate conserved genes in the orders Asterles and Dipsacales.
GeneAsteralesDipsacales
Chrysanthemum borealeChrysanthemum indicumCodonopsis lanceolataDiplostephium hartwegiiHelianthus annuus 1Helianthus annuus 2Helianthus annuus 3Helianthus annuus 4Helianthus annuus 5Helianthus annuus 6Lactuca salignaLactuca sativaLactuca serriolaParaprenanthes iversifoliaPlatycodon grandiflorusValeriana sambucifolia
f. dageletiana
atp1+ a+++++++++++++++
atp4++++++++++++++++
atp6++++++D b+++++++++
atp8++++++++++++++++
atp9++++++++++++++++
ccmB++++++++++DDD+++
ccmC++++++++++++++++
ccmFc++++++++++++++++
ccmFn++++++++++++++++
cob++++++++++++++++
cox1++++++++++++++++
cox2++++++++++++++++
cox3++++++++++++++++
matR++++++++++++++++
mttB++++++++++++++++
nad1++++++++++++++++
nad2++++++++++++++++
nad3++++++++++++++++
nad4++++++++++++++++
nad4L+D++++++++++++++
nad5++++++++++DDDDP c+
nad6++P+++++++++++++
nad7++++++++++++++++
nad9++++++++++++++++
rpl5++- d+++++++++++-+
rpl10++++++++++DDD+++
rpl16PP+++++++++++++P
rps1P+++---------P+P
rps3++++++++++++++++
rps4++++++++++++++++
rps7--+-----------++
rps10---------------+
rps11-P-P++++++PPPP--
rps12++++++++++++++++
rps13+++++++++++++++P
rps14PP-PPPPPPPPPPP-+
rps19---------------P
shd3--P-----------+-
shd4++-+PPP+PPPPPP+-
a +: present; b D: duplicated; c P: pseudogene or partial gene; d -: absent.
Table 2. The results of BLASTn searches for seven mitogenomes with or without whole-genome sequences.
Table 2. The results of BLASTn searches for seven mitogenomes with or without whole-genome sequences.
TaxaAccession of Mitogenome/Length (bp)Accession of WGS a% of Non-Hit Regions
nt bnt + WGS
Valeriana sambucifolia f. dageletianaPresent study/1,187,459-31.7-
Platycodon grandiflorusNC_035958/1249,593-63.5-
Pinus taedaNC_039746/1,191,054-43.3-
Schisandra sphenantheraNC_042758/1,101,768-19.9-
Cucurbita pepo subsp. pepoNC_014050/982,833NC_036638.1~NC_036657.146.221.2
Mangifera indicaCM021857/871,458CM021837.1~CM021856.119.815.6
Acer yangbienseCM017774/803,281CM017761.1~CM017773.112.010.9
a WGS: whole-genome sequences; b nt: nucleotide database.
Table 3. Number of organelle genomes against 149 regions in the mitogenome of V. sambucifolium f. dageletiana that corresponded to sequences in the plastomes of terrestrial plants.
Table 3. Number of organelle genomes against 149 regions in the mitogenome of V. sambucifolium f. dageletiana that corresponded to sequences in the plastomes of terrestrial plants.
Number of QueriesNumber of Subjects Against QueriesGene Transfer
Percentage of 9025 PlastomesPercentage of 364 Mitogenomes
31≥40≥40Plastome → Mitogenome (mitochondrial plastid DNA (MTPT))
21≥40<40 and ≥10
16≥40<10 and ≥1
1<40 and ≥10<40 and ≥10
1<40 and ≥10<10 and ≥1
4<10 and ≥1≥40Mitogenome → Plastome (plastid mitochondrial DNA (PTMT))
1<10 and ≥1<40 and ≥10
1<10 and ≥1<10 and ≥1
42<1≥40Mitogenome → Plastome (PTMT)
7<1<40 and ≥10
17<1<10 and ≥1
7<1<1Nuclear genome → Organelle genomes
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, H.T.; Kim, J.S. Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer. Int. J. Mol. Sci. 2021, 22, 3770. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22073770

AMA Style

Kim HT, Kim JS. Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer. International Journal of Molecular Sciences. 2021; 22(7):3770. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22073770

Chicago/Turabian Style

Kim, Hyoung Tae, and Jung Sung Kim. 2021. "Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer" International Journal of Molecular Sciences 22, no. 7: 3770. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22073770

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop