Evidence of mitochondrial DNA in the chloroplast genome of Convallaria keiskei and its subsequent evolution in the Asparagales

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes

BACKGROUND

Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus.

RESULTS

We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT.

CONCLUSIONS

Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.

Characteristics of chloroplast and mitochondrial genomes and intracellular gene transfer in the Korean endemic shrub, Sophora koreensis Nakai (Fabaceae)

Sophora koreensis Nakai, an endemic species distributed only in the Korean Peninsula, is of great geographical, economic, and taxonomic importance. Although its complete chloroplast (cp) genome sequence has been reported, its mitochondrial (mt) genome sequence has not yet been studied. Therefore, in this study, we aimed to investigate its mt genome sequence and compare it with those reported for other Fabaceae species. Total genomic DNA was extracted from fresh S. koreensis leaves collected from natural habitats in Gangwon-do Province, South Korea. This was followed by polymerase chain reaction (PCR) amplification of cpDNA insertions in the mt genome and the detection of microsatellites and dispersed repeats in the cp and mt genomes. Finally, the cp and mt genomes of S. koreensis were compared with those reported for other Fabaceae species. The cp sequence of S. koreensis showed identical gene orders and contents as those previously reported. Only six substitutions and one deletion were detected with 99 % homology. Conversely, the complete mt genome sequence, which was 517,845 bp in length and encoded 61 genes, including 43 protein-coding, 15 transfer RNAs, and 3 ribosomal RNA genes, was considerably different from that of S. japonica in terms of gene order and composition. Further, the mt genome of S. koreensis included ca. 7 and 3 kb insertions, representing an intracellular gene transfer (IGT) event, and the regions with these insertions were determined to be originally present in the cp genome. This IGT event was also confirmed via PCR amplification. IGT events can be induced via biological gene expression control or the use of repetitive sequences, and they provide important insights into the evolutionary lineage of S. koreensis. However, further studies are needed to clarify the gene transfer mechanisms between the two organelles.

Population Structure and Genetic Diversity Analyses Provide New Insight into the Endemic Species Aster spathulifolius Maxim. and Its Evolutionary History

Aster spathulifolius, an ecologically significant plant species native to the coastal regions of Korea and Japan, remains understudied in terms of its genetic structure and evolutionary history. In this study, we employed four chloroplast markers and the nuclear ITS region from 15 populations of A. spathulifolius from both Korea and Japan, including their islands, to unravel the spatial genetic structure, differentiation, gene flow, phylogenetic, and biogeographical relationships. Analysis based on multiple methods identified a low level of genetic diversity, genetic differentiation and gene flow among A. spathulifolius populations. Network analysis and principal coordinates analysis showed that 15 populations could be divided into two groups: mainland and island. Furthermore, UPGMA, neighbor-net, maximum-likelihood and Bayesian inference-based phylogenetic tree confirmed that these populations formed two distinct clades. Therefore, the island populations might be treated as A. spathulifolius populations rather than A. oharai populations. Divergence time analysis estimated the divergence of A. spathulifolius lineages approximately 23.09 million years ago, while ancestral area reconstruction analysis suggested Korea as the potential origin, conflicting with alternative scenarios. These findings contribute to a comprehensive understanding of the evolutionary history, genetic structure, and adaptive strategies of A. spathulifolius in coastal environments. Our study challenges previous assumptions and underscores the necessity for further population studies to elucidate the intricate dynamics of this distinctive plant species.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Comprehensive Comparative Analyses of Aspidistra Chloroplast Genomes: Insights into Interspecific Plastid Diversity and Phylogeny

Limestone karsts are renowned for extremely high species richness and endemism. Aspidistra (Asparagaceae) is among the highly diversified genera distributed in karst areas, making it an ideal group for studying the evolutionary mechanisms of karst plants. The taxonomy and identification of Aspidistra species are mainly based on their specialized and diverse floral structures. Aspidistra plants have inconspicuous flowers, and the similarity in vegetative morphology often leads to difficulties in species discrimination. Chloroplast genomes possess variable genetic information and offer the potential for interspecies identification. However, as yet there is little information about the interspecific diversity and evolution of the plastid genomes of Aspidistra. In this study, we reported chloroplast (cp) genomes of seven Aspidistra species (A. crassifila, A. dolichanthera, A. erecta, A. longgangensis, A. minutiflora, A. nankunshanensis, and A. retusa). These seven highly-conserved plastid genomes all have a typical quartile structure and include a total of 113 unique genes, comprising 79 protein-coding genes, 4 rRNA genes and 30 tRNA genes. Additionally, we conducted a comprehensive comparative analysis of Aspidistra cp genomes. We identified eight divergent hotspot regions (trnC-GCA-petN, trnE-UUC-psbD, accD-psaI, petA-psbJ, rpl20-rps12, rpl36-rps8, ccsA-ndhD and rps15-ycf1) that serve as potential molecular markers. Our newly generated Aspidistra plastomes enrich the resources of plastid genomes of karst plants, and an investigation into the plastome diversity offers novel perspectives on the taxonomy, phylogeny and evolution of Aspidistra species.

Complete mitochondrial genome of Agrostis stolonifera: insights into structure, Codon usage, repeats, and RNA editing

BACKGROUND

Plants possess mitochondrial genomes that are large and complex compared to animals. Despite their size, plant mitochondrial genomes do not contain significantly more genes than their animal counterparts. Studies into the sequence and structure of plant mitochondrial genomes heavily imply that the main mechanism driving replication of plant mtDNA, and offer valuable insights into plant evolution, energy production, and environmental adaptation.

RESULTS

This study presents the first comprehensive analysis of Agrostis stolonifera's mitochondrial genome, characterized by a branched structure comprising three contiguous chromosomes, totaling 560,800 bp with a GC content of 44.07%. Annotations reveal 33 unique protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The predominant codons for alanine and glutamine are GCU and CAA, respectively, while cysteine and phenylalanine exhibit weaker codon usage biases. The mitogenome contains 73, 34, and 23 simple sequence repeats (SSRs) on chromosomes 1, 2, and 3, respectively. Chromosome 1 exhibits the most frequent A-repeat monomeric SSR, whereas chromosome 2 displays the most common U-repeat monomeric SSR. DNA transformation analysis identifies 48 homologous fragments between the mitogenome and chloroplast genome, representing 3.41% of the mitogenome's total length. The PREP suite detects 460 C-U RNA editing events across 33 mitochondrial PCGs, with the highest count in the ccmFn gene and the lowest in the rps7 gene. Phylogenetic analysis confirms A. stolonifera's placement within the Pooideae subfamily, showing a close relationship to Lolium perenne, consistent with the APG IV classification system. Numerous homologous co-linear blocks are observed in A. stolonifera's mitogenomes and those of related species, while certain regions lack homology.

CONCLUSIONS

The unique features and complexities of the A. stolonifera mitochondrial genome, along with its similarities and differences to related species, provide valuable insights into plant evolution, energy production, and environmental adaptation. The findings from this study significantly contribute to the growing body of knowledge on plant mitochondrial genomes and their role in plant biology.

Comparative Analysis and Identification of Terpene Synthase Genes in Convallaria keiskei Leaf, Flower and Root Using RNA-Sequencing Profiling

The 'Lilly of the Valley' species, Convallaria, is renowned for its fragrant white flowers and distinctive fresh and green floral scent, attributed to a rich composition of volatile organic compounds (VOCs). However, the molecular mechanisms underlying the biosynthesis of this floral scent remain poorly understood due to a lack of transcriptomic data. In this study, we conducted the first comparative transcriptome analysis of C. keiskei, encompassing the leaf, flower, and root tissues. Our aim was to investigate the terpene synthase (TPS) genes and differential gene expression (DEG) patterns associated with essential oil biosynthesis. Through de novo assembly, we generated a substantial number of unigenes, with the highest count in the root (146,550), followed by the flower (116,434) and the leaf (72,044). Among the identified unigenes, we focused on fifteen putative ckTPS genes, which are involved in the synthesis of mono- and sesquiterpenes, the key aromatic compounds responsible for the essential oil biosynthesis in C. keiskei. The expression of these genes was validated using quantitative PCR analysis. Both DEG and qPCR analyses revealed the presence of ckTPS genes in the flower transcriptome, responsible for the synthesis of various compounds such as geraniol, germacrene, kaurene, linalool, nerolidol, trans-ocimene and valencene. The leaf transcriptome exhibited genes related to the biosynthesis of kaurene and trans-ocimene. In the root, the identified unigenes were associated with synthesizing kaurene, trans-ocimene and valencene. Both analyses indicated that the genes involved in mono- and sesquiterpene biosynthesis are more highly expressed in the flower compared to the leaf and root. This comprehensive study provides valuable resources for future investigations aiming to unravel the essential oil-biosynthesis-related genes in the Convallaria genus.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BACKGROUND

Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.

RESULTS

The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.

CONCLUSION

The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Phylogeny and evolution of Asparagaceae subfamily Nolinoideae: new insights from plastid phylogenomics

BACKGROUND AND AIMS

Asparagaceae subfamily Nolinoideae is an economically important plant group, but the deep relationships and evolutionary history of the lineage remain poorly understood. Based on a large data set including 37 newly sequenced samples and publicly available plastomes, this study aims to better resolve the inter-tribal relationships of Nolinoideae, and to rigorously examine the tribe-level monophyly of Convallarieae, Ophiopogoneae and Polygonateae.

METHODS

Maximum likelihood (ML) and Bayesian inference (BI) methods were used to infer phylogenetic relationships of Nolinoideae at the genus level and above. The diversification history of Nolinoideae was explored using molecular dating.

KEY RESULTS

Both ML and BI analyses identically recovered five clades within Nolinoideae, respectively corresponding to Dracaeneae + Rusceae, Polygonateae + Theropogon, Ophiopogoneae, Nolineae, and Convallarieae excluding Theropogon, and most deep nodes were well supported. As Theropogon was embedded in Polygonateae, the plastome phylogeny failed to resolve Convallarieae and Polygonateae as reciprocally monophyletic. Divergence time estimation showed that the origins of most Nolinoideae genera were dated to the Miocene and Pliocene. The youthfulness of Nolinoideae genera is well represented in the three herbaceous tribes (Convallarieae, Ophiopogoneae and Polygonateae) chiefly distributed in temperate areas of the Northern Hemisphere, as the median stem ages of all 14 genera currently belonging to them were estimated at <12.37 Ma.

CONCLUSIONS

This study recovered a robust backbone phylogeny, providing new insights for better understanding the evolution and classification of Nolinoideae. Compared with the deep relationships recovered by a previous study based on transcriptomic data, our data suggest that ancient hybridization or incomplete lineage sorting may have occurred in the early diversification of Nolinoideae. Our findings will provide important reference for further study of the evolutionary complexity of Nolinoideae using nuclear genomic data. The recent origin of these herbaceous genera currently belonging to Convallarieae, Ophiopogoneae and Polygonateae provides new evidence to support the hypothesis that the global expansion of temperate habitats caused by the climate cooling over the past 15 million years may have dramatically driven lineage diversification and speciation in the Northern Hemisphere temperate flora.

Plastid phylogenomics and plastome evolution in the morning glory family (Convolvulaceae)

Convolvulaceae, the morning glories or bindweeds, is a large family containing species of economic value, including crops, traditional medicines, ornamentals, and vegetables. However, not only are the phylogenetic relationships within this group still debated at the intertribal and intergeneric levels, but also plastid genome (plastome) complexity within Convolvulaceae is not well surveyed. We gathered 78 plastomes representing 17 genera across nine of the 12 Convolvulaceae tribes. Our plastid phylogenomic trees confirm the monophyly of Convolvulaceae, place the genus Jacquemontia within the subfamily Dicranostyloideae, and suggest that the tribe Merremieae is paraphyletic. In contrast, positions of the two genera Cuscuta and Erycibe are uncertain as the bootstrap support of the branches leading to them is moderate to weak. We show that nucleotide substitution rates are extremely variable among Convolvulaceae taxa and likely responsible for the topological uncertainty. Numerous plastomic rearrangements are detected in Convolvulaceae, including inversions, duplications, contraction and expansion of inverted repeats (IRs), and losses of genes and introns. Moreover, integrated foreign DNA of mitochondrial origin was found in the Jacquemontia plastome, adding a rare example of gene transfer from mitochondria to plastids in angiosperms. In the IR of Dichondra, we discovered an extra copy of rpl16 containing a direct repeat of ca. 200 bp long. This repeat was experimentally demonstrated to trigger effective homologous recombination, resulting in the coexistence of intron-containing and -lacking rpl16 duplicates. Therefore, we propose a hypothetical model to interpret intron loss accompanied by invasion of direct repeats at appropriate positions. Our model complements the intron loss model driven by retroprocessing when genes have lost introns but contain abundant RNA editing sites adjacent to former splicing sites.

Comparative Analyses of Plastomes of Four Anubias (Araceae) Taxa, Tropical Aquatic Plants Endemic to Africa

Anubias Schott (Araceae) have high ornamental properties as aquarium plants. However, the genus has difficulties in species identification, and the mechanism of its adaptation to the aquatic environment is unknown. To better identify species and understand the evolutionary history of Anubias, the plastomes of Anubias barteri Schott, A. barteri var. nana (Engl.) Crusio, and A. hastifolia Engl., were sequenced. The sizes of the plastomes of Anubias ranged from 169,841 bp to 170,037 bp. These plastomes were composed of conserved quadripartite circular structures and comprised 112 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. The comparative analysis of genome structure, repeat sequences, codon usage and RNA editing sites revealed high similarities among the Anubias plastomes, indicating the conservation of plastomes of Anubias. Three spacer regions with relatively high nucleotide diversity, trnL-CAA-ndhB, ycf1-ndhF, and rps15-ycf1, were found within the plastomes of Anubias. Phylogenetic analysis, based on 75 protein-coding genes, showed that Anubias was sister to Montrichardia arborescens (L.) Schott (BS = 99). In addition, four genes (ccsA, matK, ndhF, and ycf4) that contain sites undergoing positive selection were identified within the Anubias plastomes. These genes may play an important role in the adaptation of Anubias to the aquatic environment. The present study provides a valuable resource for further studies on species identification and the evolutionary history of Anubias.

Evolutionary Comparison of the Complete Chloroplast Genomes in Convallaria Species and Phylogenetic Study of Asparagaceae

The genus Convallaria (Asparagaceae) comprises three herbaceous perennial species that are widely distributed in the understory of temperate deciduous forests in the Northern Hemisphere. Although Convallaria species have high medicinal and horticultural values, studies related to the phylogenetic analysis of this genus are few. In the present study, we assembled and reported five complete chloroplast (cp) sequences of three Convallaria species (two of C. keiskei Miq., two of C. majalis L., and one of C. montana Raf.) using Illumina paired-end sequencing data. The cp genomes were highly similar in overall size (161,365–162,972 bp), and all consisted of a pair of inverted repeats (IR) regions (29,140–29,486 bp) separated by a large single-copy (LSC) (85,183–85,521 bp) and a small single-copy (SSC) region (17,877–18,502 bp). Each cp genome contained the same 113 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene content, gene order, AT content and IR/SC boundary structure were nearly identical among all of the Convallaria cp genomes. However, their lengths varied due to contraction/expansion at the IR/LSC borders. Simple sequence repeat (SSR) analyses indicated that the richest SSRs are A/T mononucleotides. Three highly variable regions (petA-psbJ, psbI-trnS and ccsA-ndhD) were identified as valuable molecular markers. Phylogenetic analysis of the family Asparagaceae using 48 cp genome sequences supported the monophyly of Convallaria, which formed a sister clade to the genus Rohdea. Our study provides a robust phylogeny of the Asparagaceae family. The complete cp genome sequences will contribute to further studies in the molecular identification, genetic diversity, and phylogeny of Convallaria.

The entire chloroplast genome sequence of Asparagus cochinchinensis and genetic comparison to Asparagus species

Asparagus cochinchinensis is a traditional Chinese medicinal plant. The chloroplast (cp) genome study on A. cochinchinensis is poorly understood. In this research, we collected the data from the cp genome assembly and gene annotation of A. cochinchinensis, followed by further comparative analysis with six species in the genus Asparagus. The cp genome of A. cochinchinensis showed a circular quadripartite structure in the size of 157,095 bp, comprising a large single-copy (LSC), a small single-copy (SSC), and two inverted repeat (IR) regions. A total of 137 genes were annotated, consisting of 86 protein-coding genes, 8 ribosomal RNAs, 38 transfer RNAs, and 5 pseudo-genes. Forty scattered repetitive sequences and 247 simple sequence repeats loci were marked out. In addition, A/T-ending codons were shown to have a basis in the codon analysis. A cp genome comparative analysis revealed that a similar gene composition was detected in the IR and LSC/SSC regions with Asparagus species. Based on the complete cp genome sequence in Asparagaceae, the result showed that A. cochinchinensis was closely related to A. racemosus by phylogenetic analysis. Therefore, our study providing A. cochinchinensis genomic resources could effectively contribute to the phylogenetic analysis and molecular identification of the genus Asparagus.

Characterization of the Complete Chloroplast Genome Sequence of the Socotra Dragon`s Blood Tree (Dracaena cinnabari Balf.)

The Socotra dragon`s blood tree (Dracaena cinnabari Balf.) is endemic to the island of Socotra in Yemen. This iconic species plays an essential role in the survival of associated organisms, acting as an umbrella tree. Overexploitation, overgrazing by livestock, global climate change, and insufficient regeneration mean that the populations of this valuable species are declining in the wild. Although there are many studies on the morphology, anatomy, and physiology of D. cinnabari, no genomic analysis of this endangered species has been performed so far. Therefore, the main aim of this study was to characterize the complete chloroplast sequence genome of D. cinnabari for conservation purposes. The D. cinnabari chloroplast genome is 155,371 bp with a total GC content of 37.5%. It has a quadripartite plastid genome structure composed of one large single-copy region of 83,870 bp, one small single-copy region of 18,471 bp, and two inverted repeat regions of 26,515 bp each. One hundred and thirty-two genes were annotated, 86 of which are protein-coding genes, 38 are transfer RNAs, and eight are ribosomal RNAs. Forty simple sequence repeats have also been identified in this chloroplast genome. Comparative analysis of complete sequences of D. cinnabari chloroplast genomes with other species of the genus Dracaena showed a very high conservativeness of their structure and organization. Phylogenetic inference showed that D. cinnabari is much closer to D. draco, D. cochinchinensis, and D. cambodiana than to D. terniflora, D. angustifolia, D. hokouensis, and D. elliptica. The results obtained in this study provide new and valuable omics data for further phylogenetic studies of the genus Dracaena as well as enable the protection of genetic resources of highly endangered D. cinnabari.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Complete Plastid and Mitochondrial Genomes of Aeginetia indica Reveal Intracellular Gene Transfer (IGT), Horizontal Gene Transfer (HGT), and Cytoplasmic Male Sterility (CMS)

Orobanchaceae have become a model group for studies on the evolution of parasitic flowering plants, and Aeginetia indica, a holoparasitic plant, is a member of this family. In this study, we assembled the complete chloroplast and mitochondrial genomes of A. indica. The chloroplast and mitochondrial genomes were 56,381 bp and 401,628 bp long, respectively. The chloroplast genome of A. indica shows massive plastid genes and the loss of one IR (inverted repeat). A comparison of the A. indica chloroplast genome sequence with that of a previous study demonstrated that the two chloroplast genomes encode a similar number of proteins (except atpH) but differ greatly in length. The A. indica mitochondrial genome has 53 genes, including 35 protein-coding genes (34 native mitochondrial genes and one chloroplast gene), 15 tRNA (11 native mitochondrial genes and four chloroplast genes) genes, and three rRNA genes. Evidence for intracellular gene transfer (IGT) and horizontal gene transfer (HGT) was obtained for plastid and mitochondrial genomes. ψndhB and ψcemA in the A. indica mitogenome were transferred from the plastid genome of A. indica. The atpH gene in the plastid of A. indica was transferred from another plastid angiosperm plastid and the atpI gene in mitogenome A. indica was transferred from a host plant like Miscanthus siensis. Cox2 (orf43) encodes proteins containing a membrane domain, making ORF (Open Reading Frame) the most likely candidate gene for CMS development in A. indica.

Lineage-specific plastid degradation in subtribe Gentianinae (Gentianaceae)

The structure and sequence of plastid genomes is highly conserved across most land plants, except for a minority of lineages that show gene loss and genome degradation. Understanding the early stages of plastome degradation may provide crucial insights into the repeatability and predictability of genomic evolutionary trends. We investigated these trends in subtribe Gentianinae of the Gentianaceae, which encompasses ca. 450 species distributed around the world, particularly in alpine and subalpine environments. We sequenced, assembled, and annotated the plastomes of 41 species, representing all six genera in subtribe Gentianinae and all main sections of the species-rich genus Gentiana L. We reconstructed the phylogeny, estimated divergence times, investigated the phylogenetic distribution of putative gene losses, and related these to substitution rate shifts and species' habitats. We obtained a strongly supported topology consistent with earlier studies, with all six genera in Gentianinae recovered as monophyletic and all main sections of Gentiana having full support. While closely related species have very similar plastomes in terms of size and structure, independent gene losses, particularly of the ndh complex, have occurred in multiple clades across the phylogeny. Gene loss was usually associated with a shift in the boundaries of the small single-copy and inverted repeat regions. Substitution rates were variable between clades, with evidence for both elevated and decelerated rate shifts. Independent lineage-specific loss of ndh genes occurred at a wide range of times, from Eocene to Pliocene. Our study illustrates that diverse degradation patterns shape the evolution of the plastid in this species-rich plant group.

Chloroplast Genome of Rambutan and Comparative Analyses in Sapindaceae

Rambutan (Nephelium lappaceum L.) is an important fruit tree that belongs to the family Sapindaceae and is widely cultivated in Southeast Asia. We sequenced its chloroplast genome for the first time and assembled 161,321 bp circular DNA. It is characterized by a typical quadripartite structure composed of a large (86,068 bp) and small (18,153 bp) single-copy region interspersed by two identical inverted repeats (IRs) (28,550 bp). We identified 132 genes including 78 protein-coding genes, 29 tRNA and 4 rRNA genes, with 21 genes duplicated in the IRs. Sixty-three simple sequence repeats (SSRs) and 98 repetitive sequences were detected. Twenty-nine codons showed biased usage and 49 potential RNA editing sites were predicted across 18 protein-coding genes in the rambutan chloroplast genome. In addition, coding gene sequence divergence analysis suggested that ccsA, clpP, rpoA, rps12, psbJ and rps19 were under positive selection, which might reflect specific adaptations of N. lappaceum to its particular living environment. Comparative chloroplast genome analyses from nine species in Sapindaceae revealed that a higher similarity was conserved in the IR regions than in the large single-copy (LSC) and small single-copy (SSC) regions. The phylogenetic analysis showed that N. lappaceum chloroplast genome has the closest relationship with that of Pometia tomentosa. The understanding of the chloroplast genomics of rambutan and comparative analysis of Sapindaceae species would provide insight into future research on the breeding of rambutan and Sapindaceae evolutionary studies.

Impact of parasitic lifestyle and different types of centromere organization on chromosome and genome evolution in the plant genus Cuscuta

The parasitic genus Cuscuta (Convolvulaceae) is exceptional among plants with respect to centromere organization, including both monocentric and holocentric chromosomes, and substantial variation in genome size and chromosome number. We investigated 12 species representing the diversity of the genus in a phylogenetic context to reveal the molecular and evolutionary processes leading to diversification of their genomes. We measured genome sizes and investigated karyotypes and centromere organization using molecular cytogenetic techniques. We also performed low-pass whole genome sequencing and comparative analysis of repetitive DNA composition. A remarkable 102-fold variation in genome sizes (342-34 734 Mbp/1C) was detected for monocentric Cuscuta species, while genomes of holocentric species were of moderate sizes (533-1545 Mbp/1C). The genome size variation was primarily driven by the differential accumulation of LTR-retrotransposons and satellite DNA. The transition to holocentric chromosomes in the subgenus Cuscuta was associated with loss of histone H2A phosphorylation and elimination of centromeric retrotransposons. In addition, basic chromosome number of holocentric species (x = 7) was smaller than in monocentrics (x = 15 or 16). We demonstrated that the transition to holocentricity in Cuscuta was accompanied by significant changes in epigenetic marks, chromosome number and the repetitive DNA sequence composition.

Dissection for Floral Micromorphology and Plastid Genome of Valuable Medicinal Borages Arnebia and Lithospermum (Boraginaceae)

The genera Arnebia and Lithospermum (Lithospermeae-Boraginaceae) comprise 25-30 and 50-60 species, respectively. Some of them are economically valuable, as their roots frequently contain a purple-red dye used in the cosmetic industry. Furthermore, dried roots of Arnebia euchroma, A. guttata, and Lithospermum erythrorhizon, which have been designated Lithospermi Radix, are used as traditional Korean herbal medicine. This study is the first report on the floral micromorphology and complete chloroplast (cp) genome sequences of A. guttata (including A. tibetana), A. euchroma, and L. erythrorhizon. We reveal great diversity in floral epidermal cell patterns, gynoecium, and structure of trichomes. The cp genomes were 149,361-150,465 bp in length, with conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 tRNA genes, and four rRNA genes. Gene order, content, and orientation were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the four cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (atpH-atpI, petN-psbM, rbcL-psaI, ycf4-cemA, ndhF-rpl32, and ndhC-trnV-UAC). To facilitate species identification, we developed molecular markers psaA- ycf3 (PSY), trnI-CAU- ycf2 (TCY), and ndhC-trnV-UAC (NCTV) based on divergence hotspots. High-resolution phylogenetic analysis revealed clear clustering and a close relationship of Arnebia to its Lithospermum sister group, which was supported by strong bootstrap values and posterior probabilities. Overall, gynoecium characteristics and genetic distance of cp genomes suggest that A. tibetana, might be recognized as an independent species rather than a synonym of A. guttata. The present morphological and cp genomic results provide useful information for future studies, such as taxonomic, phylogenetic, and evolutionary analysis of Boraginaceae.

Comparative Plastid Genomics of Neotropical Bulbophyllum (Orchidaceae; Epidendroideae)

Pantropical Bulbophyllum, with ∼2,200 species, is one of the largest genera in Orchidaceae. Although phylogenetics and taxonomy of the ∼60 American species in the genus are generally well understood, some species complexes need more study to clearly delimit their component species and provide information about their evolutionary history. Previous research has suggested that the plastid genome includes phylogenetic markers capable of providing resolution at low taxonomic levels, and thus it could be an effective tool if these divergent regions can be identified. In this study, we sequenced the complete plastid genome of eight Bulbophyllum species, representing five of six Neotropical taxonomic sections. All plastomes conserve the typical quadripartite structure, and, although the general structure of plastid genomes is conserved, differences in ndh-gene composition and total length were detected. Total length was determined by contraction and expansion of the small single-copy region, a result of an independent loss of the seven ndh genes. Selection analyses indicated that protein-coding genes were generally well conserved, but in four genes, we identified 95 putative sites under positive selection. Furthermore, a total of 54 polymorphic simple sequence repeats were identified, for which we developed amplification primers. In addition, we propose 10 regions with potential to improve phylogenetic analyses of Neotropical Bulbophyllum species.

Comparative Analysis of the Complete Chloroplast Genome of Mainland Aster spathulifolius and Other Aster Species

Aster spathulifolius, a common ornamental and medicinal plant, is widely distributed in Korea and Japan, and is genetically classified into mainland and island types. Here, we sequenced the whole chloroplast genome of mainland A. spathulifolius and compared it with those of the island type and other Aster species. The chloroplast genome of mainland A. spathulifolius is 152,732 bp with a conserved quadripartite structure, has 37.28% guanine-cytosine (GC) content, and contains 114 non-redundant genes. Comparison of the chloroplast genomes between the two A. spathulifolius lines and the other Aster species revealed that their sequences, GC contents, gene contents and orders, and exon-intron structure were well conserved; however, differences were observed in their lengths, repeat sequences, and the contraction and expansion of the inverted repeats. The variations were mostly in the single-copy regions and non-coding regions, which, together with the detected simple sequence repeats, could be used for the development of molecular markers to distinguish between these plants. All Aster species clustered into a monophyletic group, but the chloroplast genome of mainland A. spathulifolius was more similar to the other Aster species than to that of the island A. spathulifolius. The accD and ndhF genes were detected to be under positive selection within the Aster lineage compared to other related taxa. The complete chloroplast genome of mainland A. spathulifolius presented in this study will be helpful for species identification and the analysis of the genetic diversity, evolution, and phylogenetic relationships in the Aster genus and the Asteraceae.