Chloroplast genomic data provide new and robust insights into the phylogeny and evolution of the Ranunculaceae

Genetic characterization and phylogenetic analysis of the Nigella sativa (black seed) plastome

In this study, the complete plastome sequence of Nigella sativa (black seed), was analyzed for the first time. The plastome spans approximately 154,120 bp, comprising four sections: the Large Single-Copy (LSC) (85,538 bp), the Small Single-Copy (SSC) (17,984 bp), and two Inverted Repeat (IR) regions (25,299 bp). A comparative study of N. sativa's plastome with ten other species from various genera in the Ranunculaceae family reveals substantial structural variations. The contraction of the inverted repeat region in N. sativa influences the boundaries of single-copy regions, resulting in a shorter plastome size than other species. When comparing the plastome of N. sativa with those of its related species, significant divergence is observed, particularly except for N. damascena. Among these, the plastome of A. glaucifolium displays the highest average pairwise sequence divergence (0.2851) with N. sativa, followed by A. raddeana (0.2290) and A. coerulea (0.1222). Furthermore, the study identified 12 distinct hotspot regions characterized by elevated Pi values (> 0.1). These regions include trnH-GUG-psbA, matK-trnQ-UUG, psbK-trnR-UCU, atpF-atpI, rpoB-psbD, ycf3-ndhJ, ndhC-cemA, petA-psaJ, trnN-GUU-ndhF, trnV-GAC-rps12, ycf2-trnI-CAU, and ndhA-ycf1. Approximately, 24 tandem and 48 palindromic and forward repeats were detected in N. sativa plastome. The analysis revealed 32 microsatellites with the majority being mononucleotide repeats. In the N. sativa plastome, phenylalanine had the highest number of codons (1982 codons), while alanine was the least common amino acid with 260 codons. A phylogenetic tree, constructed using protein-coding genes, revealed a distinct monophyletic clade comprising N. sativa and N. damascene, closely aligned with the Cimicifugeae tribe and exhibiting robust support. This plastome provides valuable genetic information for precise species identification, phylogenetic resolution, and evolutionary studies of N. sativa.

Spur development and evolution: An update

Floral spurs, widely recognized as a classic example of key morphological and functional innovation and thought to have promoted the origin and adaptive evolution of many flowering plant lineages, have attracted the attention of researchers for centuries. Despite this, the mechanisms underlying the development and evolution of these structures remain poorly understood. Recent studies have discovered the phytohormones and transcription factor genes that play key roles in regulating patterns of cell division and cell expansion during spur morphogenesis. Spur morphogenesis was also found to be tightly linked with the programs specifying floral zygomorphy, floral organ identity determination, and nectary development. Independent origins and losses of spurs in different flowering plant lineages, therefore, may be attributed to changes in the spur program and/or its upstream ones.

Genomic incongruence accompanies the evolution of flower symmetry in Eudicots: a case study in the poppy family (Papaveraceae, Ranunculales)

Reconstructing evolutionary trajectories and transitions that have shaped floral diversity relies heavily on the phylogenetic framework on which traits are modelled. In this study, we focus on the angiosperm order Ranunculales, sister to all other eudicots, to unravel higher-level relationships, especially those tied to evolutionary transitions in flower symmetry within the family Papaveraceae. This family presents an astonishing array of floral diversity, with actinomorphic, disymmetric (two perpendicular symmetry axes), and zygomorphic flowers. We generated nuclear and plastid datasets using the Angiosperms353 universal probe set for target capture sequencing (of 353 single-copy nuclear ortholog genes), together with publicly available transcriptome and plastome data mined from open-access online repositories. We relied on the fossil record of the order Ranunculales to date our phylogenies and to establish a timeline of events. Our phylogenomic workflow shows that nuclear-plastid incongruence accompanies topological uncertainties in Ranunculales. A cocktail of incomplete lineage sorting, post-hybridization introgression, and extinction following rapid speciation most likely explain the observed knots in the topology. These knots coincide with major floral symmetry transitions and thus obscure the order of evolutionary events.

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification

Introduction

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety.

Methods

The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum.

Results and discussion

The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T. fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium.

Conclusion

Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

A cornucopia of diversity-Ranunculales as a model lineage

The Ranunculales are a hyperdiverse lineage in many aspects of their phenotype, including growth habit, floral and leaf morphology, reproductive mode, and specialized metabolism. Many Ranunculales species, such as opium poppy and goldenseal, have a high medicinal value. In addition, the order includes a large number of commercially important ornamental plants, such as columbines and larkspurs. The phylogenetic position of the order with respect to monocots and core eudicots and the diversity within this lineage make the Ranunculales an excellent group for studying evolutionary processes by comparative studies. Lately, the phylogeny of Ranunculales was revised, and genetic and genomic resources were developed for many species, allowing comparative analyses at the molecular scale. Here, we review the literature on the resources for genetic manipulation and genome sequencing, the recent phylogeny reconstruction of this order, and its fossil record. Further, we explain their habitat range and delve into the diversity in their floral morphology, focusing on perianth organ identity, floral symmetry, occurrences of spurs and nectaries, sexual and pollination systems, and fruit and dehiscence types. The Ranunculales order offers a wealth of opportunities for scientific exploration across various disciplines and scales, to gain novel insights into plant biology for researchers and plant enthusiasts alike.

Comparative plastome analysis of the sister genera Ceratocephala and Myosurus (Ranunculaceae) reveals signals of adaptive evolution to arid and aquatic environments

BACKGROUND

Expansion and contraction of inverted repeats can cause considerable variation of plastid genomes (plastomes) in angiosperms. However, little is known about whether structural variations of plastomes are associated with adaptation to or occupancy of new environments. Moreover, adaptive evolution of angiosperm plastid genes remains poorly understood. Here, we sequenced the complete plastomes for four species of xerophytic Ceratocephala and hydrophytic Myosurus, as well as Ficaria verna. By an integration of phylogenomic, comparative genomic, and selection pressure analyses, we investigated evolutionary patterns of plastomes in Ranunculeae and their relationships with adaptation to dry and aquatic habitats.

RESULTS

Owing to the significant contraction of the boundary of IRA/LSC towards the IRA, plastome sizes and IR lengths of Myosurus and Ceratocephala are smaller within Ranunculeae. Compared to other Ranunculeae, the Myosurus plastome lost clpP and rps16, one copy of rpl2 and rpl23, and one intron of rpoC1 and rpl16, and the Ceratocephala plastome added an infA gene and lost one copy of rpl2 and two introns of clpP. A total of 11 plastid genes (14%) showed positive selection, two genes common to Myosurus and Ceratocephala, seven in Ceratocephala only, and two in Myosurus only. Four genes showed strong signals of episodic positive selection. The rps7 gene of Ceratocephala and the rpl32 and ycf4 genes of Myosurus showed an increase in the rate of variation close to 3.3 Ma.

CONCLUSIONS

The plastomic structure variations as well as the positive selection of two plastid genes might be related to the colonization of new environments by the common ancestor of Ceratocephala and Myosurus. The seven and two genes under positive selection might be related to the adaptation to dry and aquatic habitats in Ceratocephala and Myosurus, respectively. Moreover, intensified aridity and frequent sea-level fluctuations, as well as global cooling, might have favored an increased rate of change in some genes at about 3.3 Ma, associated with adaptation to dry and aquatic environments, respectively. These findings suggest that changing environments might have influenced structural variations of plastomes and fixed new mutations arising on some plastid genes owing to adaptation to specific habitats.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Characterizing complete mitochondrial genome of Aquilegia amurensis and its evolutionary implications

BACKGROUND

Aquilegia is a model system for studying the evolution of adaptive radiation. However, very few studies have been conducted on the Aquilegia mitochondrial genome. Since mitochondria play a key role in plant adaptation to abiotic stress, analyzing the mitochondrial genome may provide a new perspective for understanding adaptive evolution.

RESULTS

The Aquilegia amurensis mitochondrial genome was characterized by a circular chromosome and two linear chromosomes, with a total length of 538,736 bp; the genes included 33 protein-coding genes, 24 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. We subsequently conducted a phylogenetic analysis based on single nucleotide polymorphisms (SNPs) in the mitochondrial genomes of 18 Aquilegia species, which were roughly divided into two clades: the European-Asian clade and the North American clade. Moreover, the genes mttB and rpl5 were shown to be positively selected in European-Asian species, and they may help European and Asian species adapt to environmental changes.

CONCLUSIONS

In this study, we assembled and annotated the first mitochondrial genome of the adaptive evolution model plant Aquilegia. The subsequent analysis provided us with a basis for further molecular studies on Aquilegia mitochondrial genomes and valuable information on adaptive evolution in Aquilegia.

The mechanism underlying asymmetric bending of lateral petals in Delphinium (Ranunculaceae)

During the process of flower opening, most petals move downward in the direction of the pedicel (i.e., epinastic movement). In most Delphinium flowers, however, their two lateral petals display a very peculiar movement, the mirrored helical rotation, which requires the twist of the petal stalk. However, in some lineages, their lateral petals also exhibit asymmetric bending that increases the degree of mirrored helical rotation, facilitating the formation of a 3D final shape. Notably, petal asymmetric bending is a novel trait that has not been noticed yet, so its morphological nature, developmental process, and molecular mechanisms remain largely unknown. Here, by using D. anthriscifolium as a model, we determined that petal asymmetric bending was caused by the localized expansion of cell width, accompanied by the specialized array of cell wall nano-structure, on the adaxial epidermis. Digital gene analyses, gene expression, and functional studies revealed that a class I homeodomain-leucine zipper family transcription factor gene, DeanLATE MERISTEM IDENTITY1 (DeanLMI1), contributes to petal asymmetric bending; knockdown of it led to the formation of explanate 2D petals. Specifically, DeanLMI1 promotes cell expansion in width and influences the arrangement of cell wall nano-structure on the localized adaxial epidermis. These results not only provide a comprehensive portrait of petal asymmetric bending for the first time but also shed some new insights into the mechanisms of flower opening and helical movement in plants.

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.

The complete chloroplast genome of Pulsatilla chinensis f. alba D. K. Zang (Ranunculaceae, Pulsatilla Miller)

Pulsatilla chinensis f. alba D. K. Zang 1993 is a forma of Pulsatilla chinensis (Bge.) Regel, the root of P. chinensis is traditional Chinese medicine called Pulsatillae radix. The biggest difference between P. chinensis f. alba and P. chinensis is that P. chinensis f. alba sepals is white. The complete chloroplast genome of P. chinensis f. alba was sequenced using the Illumina NovaSeq platform for the first time. The lengths of the genome, large single-copy (LSC), small single-copy (SSC), two inverted repeats (IRs), and GC content were 163,654 bp, 82,355 bp, 19,069 bp, 31,115 bp, and 37.2%, respectively. It had 134 genes, including 90 protein-coding genes, 36 tRNA genes, and eight rRNA genes. The maximum-likelihood tree indicated that P. chinensis f. alba had a closer relationship with P. chinensis. This study would provide a theoretical basis for the further study of Pulsatilla plants genetics phylogenetic research.

Comparative Phylogenetic Analysis of Ancient Korean Tea "Hadong Cheon-Nyeon Cha (Camellia sinensis var. sinensis)" Using Complete Chloroplast Genome Sequences

Wild teas are valuable genetic resources for studying evolution and breeding. Here, we report the complete chloroplast genome of the ancient Korean tea 'Hadong Cheon-nyeon Cha' (C. sinensis var. sinensis), which is known as the oldest tea tree in Korea. This study determined seven Camellia sinensis var. sinenesis, including Hadong Cheon-nyeon Cha (HCNC) chloroplast genome sequences, using Illumina sequencing technology via de novo assembly. The chloroplast genome sizes ranged from 157,019 to 157,114 bp and were organized into quadripartite regions with the typical chloroplast genomes. Further, differences in SNPs and InDels were detected across the seven chloroplast genomes through variance analysis. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These genomic resources provide evolutionary insight into Korean tea plant cultivars and lay the foundation for a better understanding of the ancient Korean tea plant HCNC.

Insight into chloroplast genome structural variation of the Mongolian endemic species Adonis mongolica (Ranunculaceae) in the Adonideae tribe

Adonis mongolica is a threatened species that is endemic to Mongolia. It is a medicinal plant from the Adonis genus and has been used to treat heart diseases. However, the genomics and evolution of this species have not been thoroughly studied. We sequenced the first complete plastome of A. mongolica and compared it with ten Adonideae species to describe the plastome structure and infer phylogenetic relationships. The complete plastome of A. mongolica was 157,521 bp long and had a typical quadripartite structure with numerous divergent regions. The plastomes of Adonideae had relatively constant genome structures and sizes, except for those of Adonis. The plastome structure was consistent across Adonis. We identified a 44.8 kb large-scale inversion within the large single-copy region and rpl32 gene loss in the Adonis plastomes compared to other members of the Adonideae tribe. Additionally, Adonis had a smaller plastome size (156,917-157,603 bp) than the other genera within the tribe (159,666-160,940 bp), which was attributed to deletions of intergenic regions and partial and complete gene losses. These results suggested that an intramolecular mutation occurred in the ancestor of the Adonis genus. Based on the phylogenetic results, Adonis separated earlier than the other genera within the Adonideae tribe. The genome structures and divergences of specific regions in the Adonis genus were unique to the Adonideae tribe. This study provides fundamental knowledge for further genomic research in Mongolia and a better understanding of the evolutionary history of endemic plants.

Complete chloroplast genome sequences of the ornamental plant Prunus cistena and comparative and phylogenetic analyses with its closely related species

BACKGROUND

Prunus cistena is an excellent color leaf configuration tree for urban landscaping in the world, which has purplish red leaves, light pink flowers, plant shape and high ornamental value. Genomic resources for P. cistena are scarce, and a clear phylogenetic and evolutionary history for this species has yet to be elucidated. Here, we sequenced and analyzed the complete chloroplast genome of P. cistena and compared it with related species of the genus Prunus based on the chloroplast genome.

RESULTS

The complete chloroplast genome of P. cistena is a 157,935 bp long typical tetrad structure, with an overall GC content of 36.72% and higher GC content in the in the inverted repeats (IR) regions than in the large single-copy (LSC) and small single-copy (SSC) regions. It contains 130 genes, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The ycf3 and clpP genes have two introns, with the longest intron in the trnK-UUU gene in the LSC region. Moreover, the genome has a total of 253SSRs, with the mononucleotide SSRs being the most abundant. The chloroplast sequences and gene arrangements of P. cistena are highly conserved, with the overall structure and gene order similar to other Prunus species. The atpE, ccsA, petA, rps8, and matK genes have undergone significant positive selection in Prunus species. P. cistena has a close evolutionary relationship with P. jamasakura. The coding and IR regions are more conserved than the noncoding regions, and the chloroplast DNA sequences are highly conserved throughout the genus Prunus.

CONCLUSIONS

The current genomic datasets provide valuable information for further species identification, evolution, and phylogenetic research of the genus Prunus.

Complete Plastid Genomes of Nine Species of Ranunculeae (Ranunculaceae) and Their Phylogenetic Inferences

The tribe Ranunculeae, Ranunculaceae, comprising 19 genera widely distributed all over the world. Although a large number of Sanger sequencing-based molecular phylogenetic studies have been published, very few studies have been performed on using genomic data to infer phylogenetic relationships within Ranunculeae. In this study, the complete plastid genomes of nine species (eleven samples) from Ceratocephala, Halerpestes, and Ranunculus were de novo assembled using a next-generation sequencing method. Previously published plastomes of Oxygraphis and other related genera of the family were downloaded from GenBank for comparative analysis. The complete plastome of each Ranunculeae species has 112 genes in total, including 78 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The plastome structure of Ranunculeae samples is conserved in gene order and arrangement. There are no inverted repeat (IR) region expansions and only one IR contraction was found in the tested samples. This study also compared plastome sequences across all the samples in gene collinearity, codon usage, RNA editing sites, nucleotide variability, simple sequence repeats, and positive selection sites. Phylogeny of the available Ranunculeae species was inferred by the plastome data using maximum-likelihood and Bayesian inference methods, and data partitioning strategies were tested. The phylogenetic relationships were better resolved compared to previous studies based on Sanger sequencing methods, showing the potential value of the plastome data in inferring the phylogeny of the tribe.

Comparative plastome genomics and phylogenetic relationships of the genus Trollius

Trollius, a genus in the Ranunculaceae family, has significant medicinal and ornamental value. It is widely distributed in China with 16 different species accepted. However, due to the lack of enough samples and information sites, the molecular phylogenetic relationships of Trollius have been unresolved till now. Here we sequenced, assembled and annotated the plastomes of 16 Trollius species to investigate their genomic characteristics, inverted repeat (IR) boundaries, sequence repeats, and hypervariable loci. In addition, the phylogenetic relationships of this genus was reconstructed based on the whole plastomes and the protein-coding sequences data-sets. The plastomes of Trollius ranged between 159,597 bp and 160,202 bp in length, and contained 113 unique genes, including 79 protein coding, 30 tRNA, and 4 rRNA. The IR boundaries were relatively conserved within the genus Trollius. 959 simple sequence repeats and 657 long sequence repeats were detected in the Trollius plastomes. We identified 12 highly polymorphic loci (Pi > 0.0115) that can be used as plastid markers in molecular identification and phylogenetic investigation of the genus. Besides, Trollius was a monophyletic group with the earliest divergence clade being Trollius lilacinus Bunge, and the remaining species were divided into two strongly-supported clades. The phylogeny in our study supported the traditional classification systems based on the color of sepal, but not the previous classification system based on the types and relative lengths of the nectaries, and distribution. The genomic resources provided in our study can be used in the taxonomy of the genus Trollius, promoting the development and utilization of this genus.

Ancient allopatry and ecological divergence act together to promote plant diversity in mountainous regions: evidence from comparative phylogeography of two genera in the Sino-Himalayan region

BACKGROUND

How geographical isolation and ecological divergence act together to promote plant diversity in mountainous regions remains largely unknown. In this study, we chose two genera comprising a small number of species distributed in the Sino-Himalayan region, Megacodon (Gentianaceae) and Beesia (Ranunculaceae), which both exhibit a fragmented distribution pattern and are found across a wide range of elevations. By summarizing their common patterns of speciation and/or divergence processes, we aim to understand how environmental changes accelerated lineage diversification in the Sino-Himalayan region through ancient allopatry and ecological divergence.

RESULTS

Using ddRAD-seq, chloroplast genome sequences, and specific molecular markers, we studied the phylogenetic relationships, population structure, and historical biogeography of Beesia and Megacodon. Both genera began to diverge from the late Miocene onwards, with ancient allopatry at lower elevations formed narrow-range species or relict populations. Mantel tests between genetic distance and climatic, elevational, or geographic distance revealed an isolation-by-distance pattern in Beesia and Megacodon stylophorus. Megacodon showed two clades occupying entirely different altitudinal ranges, whereas Beesia calthifolia exhibited a genetic divergence pattern along an elevation gradient. Furthermore, we conducted morphological measurements on Beesia calthifolia and found that different elevational groups had distinct leaf shapes.

CONCLUSIONS

The regional disjunctions of plant groups in the Sino-Himalayan region are drastic and closely related to several biogeographic boundaries. As a consequence of major geological and climate change, ecological divergence when different elevations are colonized often happens simultaneously within plant groups. Although habitat fragmentation and parapatric ecological divergence each spur speciation to different extents, a combined effect of these two factors is a common phenomenon in the Sino-Himalayan region.

Complete plastid genome structure of 13 Asian Justicia (Acanthaceae) species: comparative genomics and phylogenetic analyses

BACKGROUND

Justicia L. is the largest genus in Acanthaceae Juss. and widely distributed in tropical and subtropical regions of the world. Previous phylogenetic studies have proposed a general phylogenetic framework for Justicia based on several molecular markers. However, their studies were mainly focused on resolution of phylogenetic issues of Justicia in Africa, Australia and South America due to limited sampling from Asia. Additionally, although Justicia plants are of high medical and ornamental values, little research on its genetics was reported. Therefore, to improve the understanding of its genomic structure and relationships among Asian Justicia plants, we sequenced complete chloroplast (cp.) genomes of 12 Asian plants and combined with the previously published cp. genome of Justicia leptostachya Hemsl. for further comparative genomics and phylogenetic analyses.

RESULTS

All the cp. genomes exhibit a typical quadripartite structure without genomic rearrangement and gene loss. Their sizes range from 148,374 to 151,739 bp, including a large single copy (LSC, 81,434-83,676 bp), a small single copy (SSC, 16,833-17,507 bp) and two inverted repeats (IR, 24,947-25,549 bp). GC contents range from 38.1 to 38.4%. All the plastomes contain 114 genes, including 80 protein-coding genes, 30 tRNAs and 4 rRNAs. IR variation and repetitive sequences analyses both indicated that Justicia grossa C. B. Clarke is different from other Justicia species because its lengths of ndhF and ycf1 in IRs are shorter than others and it is richest in SSRs and dispersed repeats. The ycf1 gene was identified as the candidate DNA barcode for the genus Justicia. Our phylogenetic results showed that Justicia is a polyphyletic group, which is consistent with previous studies. Among them, J. grossa belongs to subtribe Tetramerinae of tribe Justicieae while the other Justicia members belong to subtribe Justiciinae. Therefore, based on morphological and molecular evidence, J. grossa should be undoubtedly recognized as a new genus. Interestingly, the evolutionary history of Justicia was discovered to be congruent with the morphology evolution.

CONCLUSION

Our study not only elucidates basic features of Justicia whole plastomes, but also sheds light on interspecific relationships of Asian Justicia plants for the first time.

The complete chloroplast genome of Ranunculus ternatus (Ranunculaceae)

Ranunculus ternatus Thunb. 1784 is a plant with important medicinal values. Here we report its chloroplast genome. This chloroplast genome was 156,003 bp in length with a GC content of 37.86%. It is composed of a large single copy (LSC) of 85,397 bp and a small single copy (SSC) of 19,856 bp, which are separated by a pair of inverted repeats (IR) of 25,375 bp each. The chloroplast genome contained 110 unique genes, including 77 protein-coding genes, 4 rRNA genes, and 29 tRNA genes. Phylogenetic analysis indicated that R. ternatus was closely related to R. cassubicifolius. This chloroplast genome not only enriches the genome information of Ranunculus but also will be useful for the evolution study of the family Ranunculaceae.

Phylogenomic and population genomic analyses reveal the spatial-temporal dynamics of diversification of the Nigella arvensis complex (Ranunculaceae) in the Aegean archipelago

The continental-shelf islands of the Aegean Sea provide an ideal geographical setting for evolutionary-biogeographical studies but disentangling the relationships between palaeogeographical history and the times, orders of modes of taxon divergence is not straightforward. Here, we used phylogenomic and population genomic approaches, based on orthologous gene sequences and transcriptome-derived SNP data, to reconstruct the spatial-temporal evolution of the Aegean Nigella arvensis complex (Ranunculaceae; 11 out of 12 taxa). The group's early diversification in the Early/Mid-Pliocene (c. 3.77 Mya) resulted in three main lineages (Greek mainland vs. central Aegean + Turkish mainland/eastern Aegean islands), while all extant taxa are of Late Plio-/Early Pleistocene origin (c. 3.30-1.59 Mya). Demographic modelling of the outcrossing taxa uncovered disparate modes of (sub)speciation, including divergence with gene flow on the Greek mainland, para- or peripatric diversification across eastern Aegean islands, and a 'mixing-isolation-mixing (MIM)' mode of subspeciation in the Cyclades. The two selfing species (N. stricta, N. doerfleri) evolved independently from the outcrossers. Present-day island configurations are clearly insufficient to explain the spatial-temporal history of lineage diversification and modes of (sub)speciation in Aegean Nigella. Moreover, our identification of positively selected genes in almost all taxa calls into question that this plant group represents a case of 'non-adaptive' radiation. Our study revealed an episodic diversification history of the N. arvensis complex, giving new insight into the modes and drivers of island speciation and adaption across multiple spatiotemporal scales.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Comparative analysis of the complete chloroplast genome of Papaveraceae to identify rearrangements within the Corydalis chloroplast genome

Chloroplast genomes are valuable for inferring evolutionary relationships. We report the complete chloroplast genomes of 36 Corydalis spp. and one Fumaria species. We compared these genomes with 22 other taxa and investigated the genome structure, gene content, and evolutionary dynamics of the chloroplast genomes of 58 species, explored the structure, size, repeat sequences, and divergent hotspots of these genomes, conducted phylogenetic analysis, and identified nine types of chloroplast genome structures among Corydalis spp. The ndh gene family suffered inversion and rearrangement or was lost or pseudogenized throughout the chloroplast genomes of various Corydalis species. Analysis of five protein-coding genes revealed simple sequence repeats and repetitive sequences that can be potential molecular markers for species identification. Phylogenetic analysis revealed three subgenera in Corydalis. Subgenera Cremnocapnos and Sophorocapnos represented the Type 2 and 3 genome structures, respectively. Subgenus Corydalis included all types except type 3, suggesting that chloroplast genome structural diversity increased during its differentiation. Despite the explosive diversification of this subgenus, most endemic species collected from the Korean Peninsula shared only one type of genome structure, suggesting recent divergence. These findings will greatly improve our understanding of the chloroplast genome of Corydalis and may help develop effective molecular markers.

Complete chloroplast genome of Adonis pseudoamurensis W.T.Wang (Ranunculaceae)

Adonis pseudoamurensis W.T. Wang 1980 is an important traditional medicinal plant used for the treatment of cardiac diseases. The complete chloroplast (cp) genome of Adonis pseudoamurensis is reported for the first time in this study. The circular cp genome is 156,917 bp in length, consisting of a large single-copy region (86,262 bp), a small single-copy region (18,067 bp), and two inverted repeat regions (26,294 bp). The genome encodes 129 genes, comprising 84 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. Phylogenetic analysis showed that A. pseudoamurensis is closely related to A. amurensis.

Complete chloroplast genome sequence of Camellia sinensis: genome structure, adaptive evolution, and phylogenetic relationships

The chloroplast (cp) genome holds immense potential for a variety of applications including species identification, phylogenetic analysis, and evolutionary studies. In this study, we utilized Illumina NovaSeq 6000 to sequence the DNA of Camellia sinensis L. cultivar 'Zhuyeqi', followed by the assembly of its chloroplast genome using SPAdes v3.10.1, with subsequent analysis of its features and phylogenetic placement. The results showed that the cp genome of 'Zhuyeqi' was 157,072 bp, with a large single-copy region (LSC, 86,628 bp), a small single-copy region (SSC,18,282 bp), and two inverted repeat regions (IR, 26,081 bp). The total AT and GC contents of the cp genome of 'Zhuyeqi' were observed to be 62.21% and 37.29%, respectively. The cp genome encoded 135 unique genes, including 90 protein-coding genes (CDS), 37 tRNA genes, and 8 rRNA genes. Moreover, 31 codons and 247 simple sequence repeats (SSRs) were identified. The cp genomes of 'Zhuyeqi' were found to be relatively conserved, with conservation observed in the IR region, which showed no evidence of inversions or rearrangements. The five regions with the largest variations were identified, with four regions (rps12, rps19, rps16, and rpl33) located in the LSC region and one divergent region (trnI-GAU) in the IR region. Phylogenetic analysis revealed that Camellia sinensis (KJ996106.1) was closely related to 'Zhuyeqi', indicating a close phylogenetic relationship between these two species. These findings could provide important genetic information for further research into breeding of tea tree, phylogeny, and evolution of Camellia sinensis.

The First Complete Chloroplast Genome of Campanula carpatica: Genome Characterization and Phylogenetic Diversity

Campanula carpatica is an ornamental flowering plant belonging to the family Campanulaceae. The complete chloroplast genome of C. carpatica was obtained using Illumina HiSeq X and Oxford Nanopore (Nanopore GridION) platforms. The chloroplast genome exhibited a typical circular structure with a total length of 169,341 bp, comprising a large single-copy region of 102,323 bp, a small single-copy region of 7744 bp, and a pair of inverted repeats (IRa/IRb) of 29,637 bp each. Out of a total 120 genes, 76 were protein-coding genes, 36 were transfer RNA genes, and eight were ribosomal RNA genes. The genomic characteristics of C. carpatica are similar to those of other Campanula species in terms of repetitive sequences, sequence divergence, and contraction/expansion events in the inverted repeat regions. A phylogenetic analysis of 63 shared genes in 16 plant species revealed that Campanula zangezura is the closest relative of C. carpatica. Phylogenetic analysis indicated that C. carpatica was within the Campanula clade, and C. pallida occupied the outermost position of that clade.

Biogeographic diversification of Actaea (Ranunculaceae): Insights into the historical assembly of deciduous broad-leaved forests in the Northern Hemisphere

The deciduous broad-leaved forests (DBLFs) cover large temperate and subtropical high-altitude regions in the Northern Hemisphere. They are home to rich biodiversity, especially to numerous endemic and relict species. However, we know little about how this vegetation in the Northern Hemisphere has developed through time. Here, we used Actaea (Ranunculaceae), an herbaceous genus almost exclusively growing in the understory of the Northern Hemisphere DBLFs, to shed light on the historical assembly of this biome in the Northern Hemisphere. We present a complete species-level phylogenetic analysis of Actaea based on five plastid and nuclear loci. Using the phylogenetic framework, we estimated divergence times, ancestral ranges, and diversification rates. Phylogenetic analyses strongly support Actaea as monophyletic. Sections Podocarpae and Oligocarpae compose a clade, sister to all other Actaea. The sister relationship between sections Chloranthae and Souliea is strongly supported. Section Dichanthera is not monophyletic unless section Cimicifuga is included. Actaea originated in East Asia, likely the Qinghai-Tibet Plateau, in the late Paleocene (c. 57 Ma), and subsequently dispersed into North America in the middle Eocene (c. 43 Ma) via the Thulean bridge. Actaea reached Europe twice, Japan twice, and Taiwan once, and all these five colonization events occurred in the late Miocene-early Pliocene, a period when sea level dropped. Actaea began to diversify at c. 43 Ma. The section-level diversification took place at c. 27-37 Ma and the species-level diversification experienced accelerations twice, which occurred at c. 15 Ma and c. 5 Ma, respectively. Our findings suggest that the Northern Hemisphere DBLFs might have risen in the middle Eocene and further diversified in the late Eocene-Oligocene, middle Miocene and early Pliocene, in association with climatic deterioration during these four periods.

Diversity analyses in two ornamental and large-genome Ranunculaceae species based on a low-cost Klenow NGS-based protocol

Persian buttercup (Ranunculus asiaticus L.) and poppy anemone (Anemone coronaria L.) are ornamental, outcrossing, perennial species belonging to the Ranunculaceae family, characterized by large and highly repetitive genomes. We applied K-seq protocol in both species to generate high-throughput sequencing data and produce a large number of genetic polymorphisms. The technique entails the application of Klenow polymerase-based PCR using short primers designed by analyzing k-mer sets in the genome sequence. To date the genome sequence of both species has not been released, thus we designed primer sets based on the reference the genome sequence of the related species Aquilegia oxysepala var. kansuensis (Brühl). A whole of 11,542 SNPs were selected for assessing genetic diversity of eighteen commercial varieties of R. asiaticus, while 1,752 SNPs for assessing genetic diversity in six cultivars of A. coronaria. UPGMA dendrograms were constructed and in R. asiaticus integrated in with PCA analysis. This study reports the first molecular fingerprinting within Persian buttercup, while the results obtained in poppy anemone were compared with a previously published SSR-based fingerprinting, proving K-seq to be an efficient protocol for the genotyping of complex genetic backgrounds.

Pinellia genus: A systematic review of active ingredients, pharmacological effects and action mechanism, toxicological evaluation, and multi-omics application

The dried tuber of Pinellia ternata (Thunb.) Breit, Pinelliae Rhizoma (PR, also named 'Banxia' in Chinese), is widely used in traditional medicine. This review aims to provide detail summary of active ingredients, pharmacological effects, toxic ingredients, detoxification strategies, and omic researches, etc. Pharmacological ingredients from PR are mainly classified into six categories: alkaloids, amino acids, polysaccharides, phenylpropanoids, essential oils, and glucocerebrosides. Diversity of chemical composition determines the broad-spectrum efficacy and gives a foundation for the comprehensive utilization of P. ternata germplasm resources. The pharmacological compounds are involved in inhibition of cancer cells by targeting various pathways, including activation of immune system, inhibition of proliferation and cycle, induction of apoptosis, and inhibition of angiogenesis. The pharmacological components of PR act on nervous system by targeting neurotransmitters, activating immune system, decreasing apoptosis, and increasing redox system. Lectins, one major class of the toxic ingredients extracted from raw PR, possess significant toxic effects on human cells. Inflammatory factors, cytochrome P450 proteins (CYP) family enzymes, mammalian target of rapamycin (mTOR) signaling factors, transforming growth factor-β (TGF-β) signaling factors, and nervous system, are considered to be the target sites of lectins. Recently, omic analysis is widely applied in Pinellia genus studies. Plastome genome-based molecular markers are deeply used for identifying and resolving phylogeny of Pinellia genus plants. Various omic works revealed and functional identified a series of environmental stress responsive factors and active component biosynthesis-related genes. Our review summarizes the recent progress in active and toxic ingredient evaluation, pharmacological effects, detoxification strategies, and functional gene identification and accelerates efficient utilization of this traditional herb.

The diversity of elaborate petals in Isopyreae (Ranunculaceae): a special focus on nectary structure

Elaborate petals are highly diverse in morphology, structure, and epidermal differentiation and play a key role in attracting pollinators. There have been few studies on the elaborate structure of petals in the tribe Isopyreae (Ranunculaceae). Seven genera in Isopyreae (Aquilegia, Semiaquilegia, Urophysa, Isopyrum, Paraquilegia, Dichocarpum, and Leptopyrum) have petals that vary in morphology, and two genera (Enemion and Thalictrum) have no petals. The petals of nine species belonged to 7 genera in the tribe were studied to reveal their nectary structure, epidermal micromorphology and ancestral traits. The petal nectaries of Isopyreae examined in this study were located at the tip of spurs (Aquilegia yabeana and A. rockii), or the bottom of shallow sacs (Semiaquilegia adoxoides, Urophysa henryi, Isopyrum manshuricum, and Paraquilegia microphylla), a cup-shaped structure (Dichocarpum fargesii) and a bilabiate structure (Leptopyrum fumarioides). The petal nectary of eight species in Isopyreae (except A. ecalcarata) was composed of secretory epidermis, nectary parenchyma, and vascular tissues, and some sieve tubes reached the secretory parenchyma cells. Among the eight species with nectaries examined in the present study, A. yabeana had the most developed nectaries, with 10-15 layers of secretory parenchyma cells. The epidermal cells of mature petals of the nine species were divided into 11 types. Among these 11 types, there were two types of secretory cells and two types of trichomes. Aquilegia yabeana and A. rockii had the highest number of cell types (eight types), and I. manshuricum and L. fumarioides had the lowest number of cell types (three types). Aquilegia ecalcarata had no secretory cells, and the papillose conical polygonal secretory cells of D. fargesii were different from those of the other seven species with nectaries. Trichomes were found only in Aquilegia, Semiaquilegia, Urophysa, and Paraquilegia. The ancestral mode of nectar presentation in Isopyreae was petals with hidden nectar (70.58%). The different modes of nectar presentation in petals may reflect adaptations to different pollinators in Isopyreae.

Phylogenomics and phylogeography of Menispermum (Menispermaceae)

INTRODUCTION

Phylogenomics have been widely used to resolve ambiguous and controversial evolutionary relationships among plant species and genera, and the identification of unique indels in plastomes may even help to understand the evolution of some plant families. Menispermum L. (Menispermaceae) consists of three species, M. dauricum DC., M. canadense L., and M. mexicanum Rose, which are disjuncly distributed among East Asia, Eastern North America and Mexico. Taxonomists continue to debate whether M. mexicanum is a distinct species, a variety of M. dauricum, or simply a synonym of M. canadense. To date, no molecular systematics studies have included this doubtful species in phylogenetic analyses.

METHODS

In this study, we examined phylogenomics and phylogeography of Menispermum across its entire range using 29 whole plastomes of Menispermaceae and 18 ITS1&ITS2 sequences of Menispermeae. We reconstructed interspecific relationships of Menispermum and explored plastome evolution in Menispermaceae, revealing several genomic hotspot regions for the family.

RESULTS AND DISCUSSION

Phylogenetic and network analyses based on whole plastome and ITS1&ITS2 sequences show that Menispermum clusters into two clades with high support values, Clade A (M. dauricum) and Clade B (M. canadense + M. mexicanum). However, M. mexicanum is nested within M. canadense and, as a result, we support that M. mexicanum is a synonym of M. canadense. We also identified important molecular variations in the plastomes of Menispermaceae. Several indels and consequently premature terminations of genes occur in Menispermaceae. A total of 54 regions were identified as the most highly variable plastome regions, with nucleotide diversity (Pi) values > 0.05, including two coding genes (matK, ycf1), four introns (trnK intron, rpl16 intron, rps16 intron, ndhA intron), and 48 intergenic spacer (IGS) regions. Of these, four informative hotspot regions (trnH-psbA, ndhF-rpl32, trnK-rps16, and trnP-psaJ) should be especially useful for future studies of phylogeny, phylogeography and conservation genetics of Menispermaceae.

Phylogeny and adaptive evolution of subgenus Rhizirideum (Amaryllidaceae, Allium) based on plastid genomes

The subgenus Rhizirideum in the genus Allium consists of 38 species worldwide and forms five sections (A. sect. Rhizomatosa, A. sect. Tenuissima, A. sect. Rhizirideum, A. sect. Eduardia, and A. sect. Caespitosoprason), A. sect. Caespitosoprason being merged into A. sect. Rhizomatosa recently. Previous studies on this subgenus mainly focused on separate sections. To investigate the inter-section and inter-subgenera phylogenetic relationships and adaptive evolution of A. subg. Rhizirideum, we selected thirteen representative species, which cover five sections of this subgenus and can represent four typical phenotypes of it. We conducted the comparative plastome analysis with our thirteen plastomes. And phylogenetic inferences with CDSs and complete sequences of plastomes of our thirteen species and another fifty-four related species were also performed. As a result, the A. subg. Rhizirideum plastomes were relatively conservative in structure, IR/SC borders, codon usage, and repeat sequence. In phylogenetic results, the inter-subgenera relationships among A. subg. Rhizirideum and other genus Allium subgenera were generally similar to the previous reports. In contrast, the inter-section relationships within our subgenus A. subg. Rhizirideum were newly resolved in this study. A. sect. Rhizomatosa and A. sect. Tenuissima were sister branches, which were then clustered with A. sect. Rhizirideum and A. sect. Eduardia successively. However, Allium Polyrhizum Turcz. ex Regel, type species of A. sect. Caespitosoprason, was resolved as the basal taxon of A. subg. Rhizirideum. Allium siphonanthum J. M. Xu was also found in clade A. subg. Cyathophora instead of clade A. subg. Rhizirideum. The selective pressure analysis was also conducted, and most protein-coding genes were under purifying selection. At the same time, just one gene, ycf2, was found under positive selection, and another three genes (rbcL, ycf1a, ycf1b) presented relaxed selection, which were all involved in the photosynthesis. The low temperature, dry climate, and high altitude of the extreme habitats where A. subg. Rhizirideum species grow might impose intense natural selection forces on their plastome genes for photosynthesis. In summary, our research provides new insights into the phylogeny and adaptive evolution of A. subg. Rhizirideum. Moreover, we suggest that the positions of the A. subg. Rhizirideum species A. polyrhizum and A. siphonanthum should be reconsidered.

Chloroplast genomes of four Carex species: Long repetitive sequences trigger dramatic changes in chloroplast genome structure

The chloroplast genomes of angiosperms usually have a stable circular quadripartite structure that exhibits high consistency in genome size and gene order. As one of the most diverse genera of angiosperms, Carex is of great value for the study of evolutionary relationships and speciation within its genus, but the study of the structure of its chloroplast genome is limited due to its highly expanded and restructured genome with a large number of repeats. In this study, we provided a more detailed account of the chloroplast genomes of Carex using a hybrid assembly of second- and third-generation sequencing and examined structural variation within this genus. The study revealed that chloroplast genomes of four Carex species are significantly longer than that of most angiosperms and are characterized by high sequence rearrangement rates, low GC content and gene density, and increased repetitive sequences. The location of chloroplast genome structural variation in the species of Carex studied is closely related to the positions of long repeat sequences; this genus provides a typical example of chloroplast structural variation and expansion caused by long repeats. Phylogenetic relationships constructed based on the chloroplast protein-coding genes support the latest taxonomic system of Carex, while revealing that structural variation in the chloroplast genome of Carex may have some phylogenetic significance. Moreover, this study demonstrated a hybrid assembly approach based on long and short reads to analyze complex chloroplast genome assembly and also provided an important reference for the analysis of structural rearrangements of chloroplast genomes in other taxa.

Comparative analysis of chloroplast genomes of endangered heterostylous species Primula wilsonii and its closely related species

Primula, well known for its heterostyly, is the largest genus in the family Primulaceae with more than 500 species. The considerable species number has introduced a huge challenge for taxonomy. The phylogenetic relationships among Primula still maintain unresolved due to frequent hybridization and introgression between closely related species. In this study, we sequenced and assembled the complete chloroplast genomes of Primula wilsonii Dunn, which is a PSESP (plant species with extremely small populations), using Illumina sequencing and compared its genomic sequences with those of four related Primula species. The chloroplast genomes of Primula species were similar in the basic structure, gene order, and GC content. The detected 38 SSRs (simple sequence repeats) loci and 17 hypervariable regions had many similarities in P. wilsonii, P. anisodora, P. miyabeana, and P. poissonii, but showed a significant difference compared with those in P. secundiflora. Slight variations were observed among Primula chloroplast genomes, in consideration of the relatively stable patterns of IR (inverted repeats) contraction and expansion. Phylogenetic analysis based on chloroplast genomes and protein-coding genes confirmed three major clades in Chinese Primula, but the infrageneric sections were not in accordance with morphological traits. The P. poissonii complex was confirmed here and P. anisodora was the most closely related species to P. wilsonii. Overall, the chloroplast genome sequences provided useful genetic and evolutionary information for phylogeny and population genetics on Chinese Primula species.

Complete chloroplast genome sequence of Eranthis byunsanensis B.Y. Sun (Ranunculaceae), an endemic species in Korea

The Korean endemic Eranthis byunsanensis B.Y. Sun, 1993 (Ranunculaceae) is a rare plant distributed in the southwestern part of the Korean Peninsula. The complete chloroplast (cp) genome of E. byunsanensis was sequenced by next-generation sequencing (NGS) using an Illumina HiSeq X platform. The cp genome of E. byunsanensis is 160,324 bp in length with 37.9% GC content. It showed a typical quadripartite structure consisting of a pair of inverted repeats (IRs; 28,356 bp), a large single-copy region (LSC; 87,671 bp), and a small single-copy region (SSC; 15,941 bp). The cp genome comprises 130 genes including 85 protein-coding genes (PCGs), 37 tRNA genes, and eight rRNA genes. The molecular phylogenetic analysis indicates that E. byunsanensis is closely related to Eranthis stellata, both of which belong to the genus Eranthis.

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

Evolutionary history of genus Coptis and its dynamic changes in the potential suitable distribution area

The genus Coptis belongs to the Ranunculaceae family, containing 15 recognized species highly diverse in morphology. It is a conspicuous taxon with special evolutionary position, distribution pattern and medicinal value, which makes it to be of great research and conservation significance. In order to better understand the evolutionary dynamics of Coptis and promote more practical conservation measures, we performed plastome sequencing and used the sequencing data in combination with worldwide occurrence data of Coptis to estimate genetic diversity and divergence times, rebuild biogeographic history and predict its potential suitable distribution area. The average nucleotide diversity of Coptis was 0.0067 and the hotspot regions with the highest hypermutation levels were located in the ycf1 gene. Coptis is most likely to have originated in North America and Japanese archipelago and has a typical Eastern Asian and North American disjunct distribution pattern, while the species diversity center is located in Mid-West China and Japan. The crown age of the genus is estimated at around 8.49 Mya. The most suitable climatic conditions for Coptis were as follows: precipitation of driest quarter > 25.5 mm, annual precipitation > 844.9 mm and annual mean temperature -3.1 to 19 °C. The global and China suitable area shows an upward trend in the future when emission of greenhouse gases is well controlled, but the area, especially in China, decreases significantly without greenhouse gas policy interventions. The results of this study provide a comprehensive insight into the Coptis evolutionary dynamics and will facilitate future conservation efforts.

The phylogeny of Seseli (Apiaceae, Apioideae): insights from molecular and morphological data

BACKGROUND

The genus Seseli L., which consists of 125-140 species distributed in the Old World from western Europe and northwestern Africa to China and Japan, is one of the largest and most taxonomically difficult genera of Apiaceae Lindl. Although several previous studies have been conducted on Seseli based on limited morphological characteristics and molecular fragments, a robust and comprehensive phylogeny of Seseli remains elusive. Plastomes provide abundant genetic information and have been widely used in studying plant phylogeny and evolution. Consequently, we newly generated the complete plastomes of eleven Seseli taxa. We combined plastome data and morphological characteristics to investigate the phylogeny of Seseli.

RESULTS

In our study, we observed that the genome length, gene numbers, IR/SC borders, and repeat composition of the eleven Seseli plastomes were variable. Several appropriate mutation hotspot regions may be developed as candidate DNA barcodes for evolution, phylogeny, and species identification of Seseli. The phylogenetic results identified that Seseli was not a monophyletic group. Moreover, the eleven newly sequenced Seseli taxa did not cluster with S. tortuosum (the type species of Seseli, belonging to the tribe Selineae), where S. delavayi clustered with Eriocycla belonging to the tribe Echinophoreae and the other ten belonged to Selineae. The comparative plastome and morphological characteristics analyses confirmed the reliability of the phylogenetic analyses and implied the complex evolution of Seseli.

CONCLUSION

Combining molecular and morphological data is efficient and useful for studying the phylogeny of Seseli. We suggest that "a narrow sense" of Seseli will be meaningful for further study and the current taxonomic system of Seseli needs to be revised. In summary, our study can provide new insights into the phylogenetic relationships and taxonomic framework of Seseli.

Genome-partitioning strategy, plastid and nuclear phylogenomic discordance, and its evolutionary implications of Clematis (Ranunculaceae)

Clematis is one of the largest genera of Ranunculaceae with many phylogenetic problems left to be resolved. Clematis species have considerable genome size of more than 7 Gbp, and there was no whole-genome reference sequence published in this genus. This raises difficulties in acquiring nuclear genome data for its phylogenetic analysis. Previous studies based on Sanger sequencing data, plastid genome data, and nrDNA sequences did not well resolve the phylogeny of Clematis. In this study, we used genome skimming and transcriptome data to assemble the plastid genome sequences, nuclear single nucleotide polymorphisms (SNPs) datasets, and single-copy nuclear orthologous genes (SCOGs) to reconstruct the phylogenetic backbone of Clematis, and test effectiveness of these genome partitioning methods. We also further analyzed the discordance among nuclear gene trees and between plastid and nuclear phylogenies. The results showed that the SCOGs datasets, assembled from transcriptome method, well resolved the phylogenetic backbone of Clematis. The nuclear SNPs datasets from genome skimming method can also produce similar results with the SCOGs data. In contrast to the plastid phylogeny, the phylogeny resolved by nuclear genome data is more robust and better corresponds to morphological characters. Our results suggested that rapid species radiation may have generated high level of incomplete lineage sorting, which was the major cause of nuclear gene discordance. Our simulation also showed that there may have been frequent interspecific hybridization events, which led to some of the cyto-nuclear discordances in Clematis. This study not only provides the first robust phylogenetic backbone of Clematis based on nuclear genome data, but also provides suggestions of genome partitioning strategies for the phylogenomic study of other plant taxa.

Genotyping-by-sequencing reveals range expansion of Adonis vernalis (Ranunculaceae) from Southeastern Europe into the zonal Euro-Siberian steppe

The Euro-Siberian steppe flora consists of warm- and cold-adapted species, which may have responded differently to Pleistocene glacials and interglacials. Genotyping-by-sequencing individuals from across the distribution range of the pheasant's eye (Adonis vernalis), we aimed to gain insight into steppe florogenesis based on the species' evolutionary history. Although the primary area of origin of the species group comprising A. vernalis, A. villosa and A. volgensis is in Asia, our results indicate that recent populations of A. vernalis are not of Asian origin but evolved in the southern part of Europe during the Pleistocene, with Spanish populations clearly genetically distinct from the Southeastern European populations. We inferred that A. vernalis migrated eastwards from the sub-Mediterranean forest-steppes of Southeastern Europe into the continental forest-steppe zone. Eastern European populations had the highest private allelic richness, indicating long-term large population sizes in this region. As a thermophilic species, A. vernalis seems unlikely to have survived in the cold deserts of the Last Glacial Maximum in Western Siberia, so this region was likely (re)colonized postglacially. Overall, our results reinforce the importance of identifying the area of origin and the corresponding ecological requirements of steppe plants in order to understand the composition of today's steppe flora.

Application of Analytical Technologies in the Discrimination and Authentication of Herbs from Fritillaria: A Review

Medicinal plants of Fritillaria are widely distributed in numerous countries around the world and possess excellent antitussive and expectorant effects. In particular, Fritillariae Bulbus (FB) as a precious traditional medicine has thousands of years of medical history in China. Herbs of Fritillaria have a high market value and demand while limited by harsh growing circumstances and scarce wild resources. As a consequence, fraudulent behaviors are regularly engaged by the unscrupulous merchants in an attempt to reap greater profits. It is of an urgent need to evaluate the quality of Fritillaria herbs and their products using various analytical instruments and techniques. This review has scrutinized approximately 160 articles from 1995 to 2022 published on the investigation of Fritillaria herbs and related herbal products. The botanical classification of genus Fritillaria, types of counterfeits, technologies applied for differentiating Fritillaria species were comprehensively summarized and discussed in the current review. Molecular and chromatographic identification were the dominant technologies in the authentication of Fritillaria herbs. Additionally, we brought some potential and promising technologies and analytical strategies into attention, which are worthy attempting in the future researches. This review could conduce to excellent reference value for further investigations of the authenticity assessment of Fritillaria species.

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.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Nonadaptive molecular evolution of plastome during the speciation of Actaea purpurea and its relatives

We have seen an explosive increase of plant plastid genome (plastome) sequences in the last decade, and the view that sequence variation in plastomes is maintained by the mutation-drift balance has been challenged by new evidence. Although comparative genomic and population-level studies provided us with evidence for positive evolution of plastid genes at both the macro- and micro-evolution levels, less studies have systematically investigated how plastomes have evolved during the speciation process. We here sequenced 13 plastomes of Actaea purpurea (P.K. Hsiao) J. Compton, and its closest relatives, and conducted a systematic survey of positive selection in their plastid genes using the McDonald-Kreitman test and codon-based methods using maximum likelihood to estimate the ratio of nonsynonymous to synonymous substitutions (ω) across a phylogeny. We found that during the speciation of A. purpurea and its relatives, all plastid genes evolved neutrally or were under purifying selection. Genome size, gene order, and number were highly conserved. Comparing to A. purpurea, plastomes of Actaea japonica and Actaea biternata had low genetic diversity, consistent with previous studies. Our work not only sheds important light on the evolutionary history of A. purpurea and its kin, but also on the evolution of plastomes during plant speciation.

New Insights Into the Backbone Phylogeny and Character Evolution of Corydalis (Papaveraceae) Based on Plastome Data

A robust backbone phylogeny is fundamental for developing a stable classification and is instructive for further research. However, it was still not available for Corydalis DC., a species-rich (> 500 species), ecologically and medically important, but taxonomically notoriously difficult genus. Here, we constructed backbone phylogeny and estimated the divergence of Corydalis based on the plastome data from 39 Corydalis species (32 newly sequenced), which represent ca. 80% of sections and series across this genus. Our phylogenetic analyses recovered six fully supported main clades (I-VI) and provided full support for the majority of lineages within Corydalis. Section Archaeocapnos was unexpectedly turned out to be sister to the rest of the subg. Corydalis s. l. (clades IV-VI), thus treating as a distinct clade (clade III) to render all the main clades monophyletic. Additionally, some unusual plastome structural rearrangements were constantly detected within Corydalis and were proven to be lineage-specific in this study, which, in turn, provided further support to our phylogeny. A segment containing five genes (trnV-UAC-rbcL) in the plastome's LSC region was either normally located downstream of the ndhC gene in clade I species or translocated downstream of the atpH gene in clade II species or translocated to downstream of the trnK-UUU gene in clade III-VI species. The unique large inversion (ca. 50 kb) in the plastome LSC region of clade III species, representing an intermediate stage of the above translocation in clades IV-VI, firmly supported clade III as a distinct and early diverged clade within this large lineage (clades III-VI). Our phylogeny contradicted substantially with the morphology-based taxonomy, rejected the treatment of tuberous species as an independent evolutionary group, and proved that some commonly used diagnostic characters (e.g., root and rhizome) were results of convergent evolution, suggestive of unreliability in Corydalis. We dated the origin of crown Corydalis to the early Eocene (crown age 49.08 Ma) and revealed possible explosive radiation around 25 Ma, coinciding with the drastic uplift of the Qinghai-Tibetan Plateau in Oligocene and Miocene. This study provided the most reliable and robust backbone phylogeny of Corydalis to date and shed some new insights on the evolution of Corydalis.

Comparative chloroplast genomes and phylogenetic analyses of Pinellia

Background

Pinellia Tenore (Araceae) is a genus of perennial herbaceous plants, all of which have medicinal value. The chloroplast (cp) genome data of Pinellia are scarce, and the phylogenetic relationship and gene evolution remain unclear.

Methods and results

We sequenced and annotated the Pinellia pedatisecta cp genome and combined it with previously published genomes for other Pinellia species. We used bioinformatics methods to analyse the genomic structure, repetitive sequences, interspecific variation, divergence hotspots, phylogenetic relationships, divergence time estimation and selective pressure of four Pinellia plastomes. Results showed that the cp genomes of Pinellia varied in length between 168,178 (P. pedatisecta MN046890) and 164,013 bp (P. ternata KR270823). A total of 68–111 SSR loci were identified as candidate molecular markers for further genetic diversity study. Eight mutational hotspot regions were determined, including psbI-trnG-UCC, psbM-rpoB, ndhJ-trnT-UGU, trnP-UGG-trnW-CCA, ndhF-trnN-GUU, ndhG-ndhE, ycf1-rps15 and trnR-ycf1. Gene selection pressure suggested that four genes were subjected to positive selection. Phylogenetic inferences based on the complete cp genomes revealed a sister relationship between Pinellia and Arisaema plants whose divergence was estimated to occur around 22.48 million years ago. All Pinellia species formed a monophyletic evolutionary clade in which P. peltata, rather than P. pedatisecta, earlier diverged, indicating that P. pedatisecta is not the basal taxon of Pinellia but P. peltata may be.

Conclusions

The cp genomes of Pinellia will provide valuable information for species classification, identification, molecular breeding and evolutionary exploration of the genus Pinellia.

Supplementary Information

The online version of this article (10.1007/s11033-022-07617-5) contains supplementary material, which is available to authorized users.

A phylotranscriptome study using silica gel-dried leaf tissues produces an updated robust phylogeny of Ranunculaceae

The utility of transcriptome data in plant phylogenetics has gained popularity in recent years. However, because RNA degrades much more easily than DNA, the logistics of obtaining fresh tissues has become a major limiting factor for widely applying this method. Here, we used Ranunculaceae to test whether silica-dried plant tissues could be used for RNA extraction and subsequent phylogenomic studies. We sequenced 27 transcriptomes, 21 from silica gel-dried (SD-samples) and six from liquid nitrogen-preserved (LN-samples) leaf tissues, and downloaded 27 additional transcriptomes from GenBank. Our results showed that although the LN-samples produced slightly better reads than the SD-samples, there were no significant differences in RNA quality and quantity, assembled contig lengths and numbers, and BUSCO comparisons between two treatments. Using these data, we conducted phylogenomic analyses, including concatenated- and coalescent-based phylogenetic reconstruction, molecular dating, coalescent simulation, phylogenetic network estimation, and whole genome duplication (WGD) inference. The resulting phylogeny was consistent with previous studies with higher resolution and statistical support. The 11 core Ranunculaceae tribes grouped into two chromosome type clades (T- and R-types), with high support. Discordance among gene trees is likely due to hybridization and introgression, ancient genetic polymorphism and incomplete lineage sorting. Our results strongly support one ancient hybridization event within the R-type clade and three WGD events in Ranunculales. Evolution of the three Ranunculaceae chromosome types is likely not directly related to WGD events. By clearly resolving the Ranunculaceae phylogeny, we demonstrated that SD-samples can be used for RNA-seq and phylotranscriptomic studies of angiosperms.

Petal development and elaboration

Petals can be simple or elaborate, depending on whether they have complex basic structures and/or highly specialized epidermal modifications. It has been proposed that the independent origin and diversification of elaborate petals have promoted plant-animal interactions and, therefore, the evolutionary radiation of corresponding plant groups. Recent advances in floral development and evolution have greatly improved our understanding of the processes, patterns, and mechanisms underlying petal elaboration. In this review, we compare the developmental processes of simple and elaborate petals, concluding that elaborate petals can be achieved through four main paths of modifications (i.e. marginal elaboration, ventral elaboration, dorsal elaboration, and surface elaboration). Although different types of elaborate petals were formed through different types of modifications, they are all results of changes in the expression patterns of genes involved in organ polarity establishment and/or the proliferation, expansion, and differentiation of cells. The deployment of existing genetic materials to perform a new function was also shown to be a key to making elaborate petals during evolution.

Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences

BACKGROUND

Polystachya Hook. is a large pantropical orchid genus (c. 240 species) distributed in Africa, southern Asia and the Americas, with the center of diversity in Africa. Previous studies on species of this genus have not obtained the complete chloroplast genomes, structures and variations. Additionally, the phylogenetic position of the genus in the Orchidaceae is still controversial and uncertain. Therefore, in this study, we sequenced the complete plastomes of six Kenya Polystachya species based on genome skimming, subjected them to comparative genomic analysis, and reconstructed the phylogenetic relationships with other Orchidaceae species.

RESULTS

The results exhibited that the chloroplast genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The plastomes of the six Polystachya species ranged from 145,484 bp to 149,274 bp in length and had an almost similar GC content of 36.9-37.0%. Gene annotation revealed 106-109 single-copy genes. In addition, 19 genes are duplicated in the inverted regions, and 16 genes each possessd one or more introns. Although no large structural variations were observed among the Polystachya plastomes, about 1 kb inversion was found in Polystachya modesta and all 11 ndh genes in the Polystachya plastomes were lost or pseudogenized. Comparative analysis of the overall sequence identity among six complete chloroplast genomes confirmed that for both coding and non-coding regions in Polystachya, SC regions exhibit higher sequence variation than IRs. Furthermore, there were various amplifications in the IR regions among the six Polystachya species. Most of the protein-coding genes of these species had a high degree of codon preference. We screened out SSRs and found seven relatively highly variable loci. Moreover, 13 genes were discovered with significant positive selection. Phylogenetic analysis showed that the six Polystachya species formed a monophyletic clade and were more closely related to the tribe Vandeae. Phylogenetic relationships of the family Orchidaceae inferred from the 85 chloroplast genome sequences were generally consistent with previous studies and robust.

CONCLUSIONS

Our study is the initial report of the complete chloroplast genomes of the six Polystachya species, elucidates the structural characteristics of the chloroplast genome of Polystachya, and filters out highly variable sequences that can contribute to the development of DNA markers for use in the study of genetic variability and evolutionary studies in Polystachya. In addition, the phylogenetic results strongly support that the genus of Polystachya is a part of the tribe Vandeae.