Plastid phylogenomic insights into the evolution of Caryophyllales

Characteristics and phylogenetic analysis of the complete chloroplast genome of Mesembryanthemum cordifolium L. F. (Aizoaceae)

Mesembryanthemum cordifolium, a perennial plant with crassulacean acid metabolism (CAM) in the Aizoaceae family, has significant ornamental and medicinal values. In this study, we reported the first complete chloroplast genome sequence of this species. The total genome size was 153,734 bp in length, including a large single-copy (LSC) region of 85,692 bp, a small single-copy (SSC) region of 18,212 bp, and a pair of inverted repeat (IR) regions of 24,915 bp by each. The overall GC content of the M. cordifolium chloroplast genome was 37.08%. The genome encodes 131 genes, comprising 87 protein-coding genes (PCGs), 36 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). Phylogenetic analysis shows this species was relatively close to M. crystallinum. This chloroplast genome sequence will be valuable for species discrimination and for understanding phylogenetic relationships within the genus Mesembryanthemum.

Comparative analysis of the complete chloroplast genomes of thirteen Bougainvillea cultivars from South China with implications for their genome structures and phylogenetic relationships

Bougainvillea spp., belonging to the Nyctaginaceae family, have high economic and horticultural value in South China. Despite the high similarity in terms of leaf appearance and hybridization among Bougainvillea species, especially Bougainvillea × buttiana, their phylogenetic relationships are very complicated and controversial. In this study, we sequenced, assembled and analyzed thirteen complete chloroplast genomes of Bougainvillea cultivars from South China, including ten B. × buttiana cultivars and three other Bougainvillea cultivars, and identified their phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family for the first time. These 13 chloroplast genomes had typical quadripartite structures, comprising a large single-copy (LSC) region (85,169-85,695 bp), a small single-copy (SSC) region (18,050-21,789 bp), and a pair of inverted-repeat (IR) regions (25,377-25,426 bp). These genomes each contained 112 different genes, including 79 protein-coding genes, 29 tRNAs and 4 rRNAs. The gene content, codon usage, simple sequence repeats (SSRs), and long repeats were essentially conserved among these 13 genomes. Single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were detected among these 13 genomes. Four divergent regions, namely, trnH-GUG_psbA, trnS-GCU_trnG-UCC-exon1, trnS-GGA_rps4, and ccsA_ndhD, were identified from the comparative analysis of 16 Bougainvillea cultivar genomes. Among the 46 chloroplast genomes of the Nyctaginaceae family, nine genes, namely, rps12, rbcL, ndhF, rpoB, rpoC2, ndhI, psbT, ycf2, and ycf3, were found to be under positive selection at the amino acid site level. Phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family based on complete chloroplast genomes and protein-coding genes revealed that the Bougainvillea genus was a sister to the Belemia genus with strong support and that 35 Bougainvillea individuals were divided into 4 strongly supported clades, namely, Clades Ⅰ, Ⅱ, Ⅲ and Ⅳ. Clade Ⅰ included 6 individuals, which contained 2 cultivars, namely, B. × buttiana 'Gautama's Red' and B. spectabilis 'Flame'. Clades Ⅱ only contained Bougainvillea spinosa. Clade Ⅲ comprised 7 individuals of wild species. Clade Ⅳ included 21 individuals and contained 11 cultivars, namely, B. × buttiana 'Mahara', B. × buttiana 'California Gold', B. × buttiana 'Double Salmon', B. × buttiana 'Double Yellow', B. × buttiana 'Los Banos Beauty', B. × buttiana 'Big Chitra', B. × buttiana 'San Diego Red', B. × buttiana 'Barbara Karst', B. glabra 'White Stripe', B. spectabilis 'Splendens' and B. × buttiana 'Miss Manila' sp. 1. In conclusion, this study not only provided valuable genome resources but also helped to identify Bougainvillea cultivars and understand the chloroplast genome evolution of the Nyctaginaceae family.

Phylogeographic analysis reveals multiple origins of the desert shrub Reaumuria songarica in northern Xinjiang, involving homoploid and tetraploid hybrids

Hybrid speciation plays an important role in species diversification. The establishment of reproductive isolation is crucial for hybrid speciation, and the identification of diverse types of hybrids, particularly homoploid hybrid species, contributes to a comprehensive understanding of this process. Reaumuria songarica is a constructive shrub widespread in arid Central Asia. Previous studies have inferred that the R. songarica populations in the Gurbantunggut Desert (GuD) originated from homoploid hybridizations between its eastern and western lineages and may have evolved into an incipient species. To further elucidate the genetic composition of different hybrid populations and to determine the species boundary of this hybrid lineage, we investigated the overall phylogeographic structure of R. songarica based on variation patterns of five cpDNA and one nrITS sequences across 32 populations. Phylogenetic analyses demonstrated that within the GuD lineage, the Wuerhe population evolved directly from ancestral lineages, whereas the others originated from hybridizations between the eastern and western lineages. PCoA and genetic barrier analysis supported the subdivision of the GuD lineage into the southern (GuD-S) and northern (GuD-N) groups. Populations in the GuD-S group had a consistent genetic composition and the same ancestral female parent, indicating that they belonged to a homoploid hybrid lineage. However, the GuD-N group experienced genetic admixture of the eastern and western lineages on nrITS and cpDNA, with some populations inferred to be allopolyploid based on ploidy data. Based on cpDNA haplotypes, BEAST analyses showed that the GuD-S and GuD-N groups originated after 0.5 Ma. Our results suggest that multiple expansions and contractions of GuD, driven by Quaternary climatic oscillations and the Kunlun-Yellow River tectonic movement, are important causes of the complex origins of R. songarica populations in northern Xinjiang. This study highlights the complex origins of the Junggar Basin flora and the underappreciated role of hybridization in increasing its species diversity.

Characteristics of the chloroplast genome and genetic divergence of Tamarix hispida Willd. 1816 (Tamaricaceae)

Tamarix hispida Willd. 1816, a crucial native plant species in the arid desert region of northwestern China, plays a significant role in maintaining ecological stability. It is instrumental in addressing soil salinity-alkalinity and heavy metal pollution. This research aims to analyze the phylogenetic divergence pattern and evolutionary history of T. hispida by comparing chloroplast genome structures across different populations. Despite the minimal differences in chloroplast genome structure due to conserved genes and junction regions, sequencing was conducted using the Illumina NovaSeq platform to verify the historical evolutionary processes between different populations, followed by assembly and annotation. The results revealed that the T. hispida chloroplast genome is approximately 156,164-156,186 bp in length, with a quadripartite structure and 131 annotated genes. Phylogenetic analysis indicated two lineages within T. hispida, with a divergence time of 3.15 Ma. These findings emphasize the low genetic diversity in T. hispida and offer valuable insights into its evolutionary past. To effectively protect and manage this species, increased scientific research and monitoring of its genetic diversity are necessary. This study underscores the importance of comprehending the genetic mechanisms behind species divergence to develop informed conservation strategies.

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

A taxonomic backbone for the Plumbaginaceae (Caryophyllales)

A taxonomic backbone of the Plumbaginaceae is presented and the current state of knowledge on phylogenetic relationships and taxon limits is reviewed as a basis for the accepted taxon concepts. In total, 4,476 scientific names and designations are treated of which 30 are not in the family Plumbaginaceae. The Plumbaginaceae are subdivided in three tribes with 26 genera and 1,179 accepted species. Two subgenera, 17 sections, two subsections and 187 infraspecific taxa are accepted. At the species and infraspecific level 2,782 synonyms were assigned to accepted taxa, whereas 194 names were excluded from the core checklist (i.e., unplaced taxa, infrageneric subdivisions with still uncertain application, names of verified uncertain application, invalid horticultural names, excluded names from other families, other excluded designations, and unresolved names). The EDIT Platform for Cybertaxonomy was utilized as the tool to compile and manage the names and further taxonomic data under explicit taxon concepts. Secundum references are given in case taxon concepts were taken from the literature, whereas this study serves as reference for newly circumscribed taxa. The family's division into the tribes Aegialitideae, Limonieae, and Plumbagineae departs from earlier two-subfamily classifications, prompted by recent phylogenetic findings that challenge the subfamilial affinity of Aegialitis. The genus Acantholimon was extended to include Gladiolimon, as currently available phylogenetic and morphological data support this merger. In Limonium, all accepted species could be assigned to sections and subsections or the "Mediterranean lineage", respectively, making use of the phylogenetic distribution of their morphological characters and states. A new combination and/or status is proposed for Dyerophytumsocotranum, Limoniumthymoides, Limonium×fraternum, Limonium×rossmaessleri, and Limoniumsect.Jovibarba. Special attention is given to nomenclatural issues, particularly for Staticenomenambiguum to resolve the names under accepted names. The use of artificial groupings like "aggregates", "complexes" and "species groups" in alpha-taxonomic treatments is discussed. The taxonomic backbone will receive continued updates and through the Caryophyllales Taxonomic Expert Network, it contributes the treatment of the Plumbaginaceae for the World Flora Online.

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera genus.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

The complete chloroplast genome of the medicinally important plant Plumbago zeylanica L. (plumbaginaceae) and phylogenetic analysis

Plumbago zeylanica L. 1753 is a medicinally-important herb in family Plumbaginaceae. In this study, we assembled and reported the complete chloroplast genome of P. zeylanica. The plastome of P. zeylanica was 169,178 bp, including a large single-copy region of 92,135 bp, a small single-copy region (SSC) of 13,455 bp and a pair of inverted repeat regions (IRs) of 31,794 bp. It contained 124 genes, including 79 protein-coding genes, 37 tRNA genes and eight rRNA genes. Phylogenetic analysis showed that P. zeylanica formed a close relationship with P. auriculata in Plumbago. The first complete chloroplast genome report of P. zeylanica providing an opportunity to explore the genetic diversity, and would be also helpful in the species identification and conservation.

The complete chloroplast genome sequence of Myricaria wardii Marquand 1929 (Tamaricaceae): a shrub species endemic to the Tibet Plateau

Myricaria wardii Marquand 1929, endemic to the Tibet Plateau, is a perennial shrub with important medicinal and ecological values. In this study, the complete chloroplast (cp) genome of M. wardii was assembled, and the phylogenetic tree was reconstructed to evaluate the phylogenetic location of the species. The results showed that the cp genome size of the M. wardii was 155,299 bp, which contained a pair of inverted repeat (IR) regions (26,150 bp), a large single copy (LSC) region (84,715 bp), and a small single copy (SSC) region (18,284 bp). The total GC content of the cp genome was 36.30%. A total of 128 genes were annotated, consisting of 83 protein-coding genes, 37 tRNA genes and 8 rRNA genes. The phylogenetic analysis showed that M. wardii was closely related to M. prostrata. This study provides useful information for the conservation of this species and the phylogenetic analysis of Tamaricaceae.

Medullary bundles in Caryophyllales: form, function, and evolution

The occurrence of conducting vascular tissue in the pith (CVTP) of tracheophytes is noteworthy. Medullary bundles, one of the remarkable examples of CVTP, evolved multiple times across angiosperms, notably in the Caryophyllales. Yet, information on the occurrence of medullary bundles is fragmented, hampering our understanding of their structure-function relationships, and evolutionary implications. Using three plastid molecular markers (matK, rbcL, and rps16 intron), a phylogeny is constructed for 561 species of Caryophyllales, and anatomical data are assembled for 856 species across 40 families to investigate the diversity of medullary bundles, their function, evolution, and diversification dynamics. Additionally, correlated evolution between medullary bundles and successive cambia was tested. Medullary bundles are ancestrally absent in Caryophyllales and evolved in core and noncore families. They are structurally diverse (e.g. number, arrangement, and types of bundles) and functionally active throughout the plant's lifespan, providing increased hydraulic conductivity, especially in herbaceous plants. Acquisition of medullary bundles does not explain diversification rate heterogeneity but is correlated to a higher diversification rate. Disparate developmental pathways were found leading to rampant convergent evolution of CVTP in Caryophyllales. These findings indicate the diversification of medullary bundles and vascular tissues as another central theme for functional and comparative molecular studies in Caryophyllales.

Phylogeny and diversification of genus Sanicula L. (Apiaceae): novel insights from plastid phylogenomic analyses

BACKGROUND

The genus Sanicula L. is a unique perennial herb that holds important medicinal values. Although the previous studies on Sanicula provided us with a good research basis, its taxonomic system and interspecific relationships have not been satisfactorily resolved, especially for those endemic to China. Moreover, the evolutionary history of this genus also remains inadequately understood. The plastid genomes possessing highly conserved structure and limited evolutionary rate have proved to be an effective tool for studying plant phylogeny and evolution.

RESULTS

In the current study, we newly sequenced and assembled fifteen Sanicula complete plastomes. Combined with two previously reported plastomes, we performed comprehensively plastid phylogenomics analyses to gain novel insights into the evolutionary history of this genus. The comparative results indicated that the seventeen plastomes exhibited a high degree of conservation and similarity in terms of their structure, size, GC content, gene order, IR borders, codon bias patterns and SSRs profiles. Such as all of them displayed a typical quadripartite structure, including a large single copy region (LSC: 85,074-86,197 bp), a small single copy region (SSC: 17,047-17,132 bp) separated by a pair of inverted repeat regions (IRs: 26,176-26,334 bp). And the seventeen plastomes had similar IR boundaries and the adjacent genes were identical. The rps19 gene was located at the junction of the LSC/IRa, the IRa/SSC junction region was located between the trnN gene and ndhF gene, the ycf1 gene appeared in the SSC/IRb junction and the IRb/LSC boundary was located between rpl12 gene and trnH gene. Twelve specific mutation hotspots (atpF, cemA, accD, rpl22, rbcL, matK, ycf1, trnH-psbA, ycf4-cemA, rbcL-accD, trnE-trnT and trnG-trnR) were identified that can serve as potential DNA barcodes for species identification within the genus Sanicula. Furthermore, the plastomes data and Internal Transcribed Spacer (ITS) sequences were performed to reconstruct the phylogeny of Sanicula. Although the tree topologies of them were incongruent, both provided strong evidence supporting the monophyly of Saniculoideae and Apioideae. In addition, the sister groups between Saniculoideae and Apioideae were strongly suggested. The Sanicula species involved in this study were clustered into a clade, and the Eryngium species were also clustered together. However, it was clearly observed that the sections of Sanicula involved in the current study were not respectively recovered as monophyletic group. Molecular dating analysis explored that the origin of this genus was occurred during the late Eocene period, approximately 37.84 Ma (95% HPD: 20.33-52.21 Ma) years ago and the diversification of the genus was occurred in early Miocene 18.38 Ma (95% HPD: 10.68-25.28 Ma).

CONCLUSION

The plastome-based tree and ITS-based tree generated incongruences, which may be attributed to the event of hybridization/introgression, incomplete lineage sorting (ILS) and chloroplast capture. Our study highlighted the power of plastome data to significantly improve the phylogenetic supports and resolutions, and to efficiently explore the evolutionary history of this genus. Molecular dating analysis explored that the diversification of the genus occurred in the early Miocene, which was largely influenced by the prevalence of the East Asian monsoon and the uplift of the Hengduan Mountains (HDM). In summary, our study provides novel insights into the plastome evolution, phylogenetic relationships, taxonomic framework and evolution of genus Sanicula.

Comprehensive Analysis of the Complete Mitochondrial Genome of Rehmannia chingii: An Autotrophic Species in the Orobanchaceae Family

Rehmannia chingii is an important medicinal plant with immense value in scientific research. However, its mitochondrial genome (mitogenome) has not yet been characterized. Herein, based on whole-genome Illumina short reads and PacBio HiFi reads, we obtained the complete mitogenome of R. chingii through a de novo assembly strategy. We carried out comparative genomic analyses and found that, in comparison with the plastid genome (plastome) showing a high degree of structural conservation, the R. chingii mitogenome structure is relatively complex, showing an intricate ring structure with 16 connections, owing to five repetitive sequences. The R. chingii mitogenome was 783,161 bp with a GC content of 44.8% and contained 77 genes, comprising 47 protein-coding genes (CDS), 27 tRNA genes, and 3 rRNA genes. We counted 579 RNA editing events in 47 CDS and 12,828 codons in all CDSs of the R. chingii mitogenome. Furthermore, 24 unique sequence transfer fragments were found between the mitogenome and plastome, comprising 8 mitogenome CDS genes and 16 plastome CDS genes, corresponding to 2.39% of the R. chingii mitogenome. Mitogenomes had shorter but more collinear regions, evidenced by a comparison of the organelles of non-parasitic R. chingii, hemiparasitic Pedicularis chinensis, and holoparasitic Aeginetia indica in the Orobanchaceae family. Moreover, from non-parasitic to holoparasitic species, the genome size in the mitogenomes of Orobanchaceae species did not decrease gradually. Instead, the smallest mitogenome was found in the hemiparasitic species P. chinensis, with a size of 225,612 bp. The findings fill the gap in the mitogenome research of the medicinal plant R. chingii, promote the progress of the organelle genome research of the Orobanchaceae family, and provide clues for molecular breeding.

Distinctive plastome evolution in carnivorous angiosperms

BACKGROUND

Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive.

RESULTS

We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection.

CONCLUSION

Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

The complete chloroplast genome of Spergularia marina (Caryophyllaceae) and its phylogenetic analysis

Spergularia marina (L.) Griseb, 1843 (Caryophyllaceae) is a halophytic plant widely distributed along the southwestern coast of the Korean Peninsula. In this study, the complete chloroplast genome sequence of S. marina was determined using next-generation sequencing (NGS). The chloroplast genome of S. marina is 152,460 bp in length with 36.7% GC content. It comprises a large single-copy (LSC; 83,321), a small single copy (SSC; 17,205 bp), and a pair of inverted repeats (IRs; 25,967 bp) with a typical quadripartite structure. It consists of 131 genes, including 86 protein-coding genes, 8 ribosomal RNAs, and 37 transfer RNAs. Phylogenetic analysis using complete chloroplast genomes showed that among the 17 Caryophyllaceae species, S. marina is most closely related to Spergula arvensis. Since no complete chloroplast genome of the genus Spergularia has been reported to date, our study provides useful genetic information for determining phylogenetic relationships within the Caryophyllaceae.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

Chromosome-scale genome sequence of Suaeda glauca sheds light on salt stress tolerance in halophytes

Soil salinity is a growing concern for global crop production and the sustainable development of humanity. Therefore, it is crucial to comprehend salt tolerance mechanisms and identify salt-tolerance genes to enhance crop tolerance to salt stress. Suaeda glauca, a halophyte species well adapted to the seawater environment, possesses a unique ability to absorb and retain high salt concentrations within its cells, particularly in its leaves, suggesting the presence of a distinct mechanism for salt tolerance. In this study, we performed de novo sequencing of the S. glauca genome. The genome has a size of 1.02 Gb (consisting of two sets of haplotypes) and contains 54 761 annotated genes, including alleles and repeats. Comparative genomic analysis revealed a strong synteny between the genomes of S. glauca and Beta vulgaris. Of the S. glauca genome, 70.56% comprises repeat sequences, with retroelements being the most abundant. Leveraging the allele-aware assembly of the S. glauca genome, we investigated genome-wide allele-specific expression in the analyzed samples. The results indicated that the diversity in promoter sequences might contribute to consistent allele-specific expression. Moreover, a systematic analysis of the ABCE gene families shed light on the formation of S. glauca's flower morphology, suggesting that dysfunction of A-class genes is responsible for the absence of petals in S. glauca. Gene family expansion analysis demonstrated significant enrichment of Gene Ontology (GO) terms associated with DNA repair, chromosome stability, DNA demethylation, cation binding, and red/far-red light signaling pathways in the co-expanded gene families of S. glauca and S. aralocaspica, in comparison with glycophytic species within the chenopodium family. Time-course transcriptome analysis under salt treatments revealed detailed responses of S. glauca to salt tolerance, and the enrichment of the transition-upregulated genes in the leaves associated with DNA repair and chromosome stability, lipid biosynthetic process, and isoprenoid metabolic process. Additionally, genome-wide analysis of transcription factors indicated a significant expansion of FAR1 gene family. However, further investigation is needed to determine the exact role of the FAR1 gene family in salt tolerance in S. glauca.

Structural and gene composition variation of the complete mitochondrial genome of Mammillaria huitzilopochtli (Cactaceae, Caryophyllales), revealed by de novo assembly

BACKGROUND

Structural descriptions of complete genomes have elucidated evolutionary processes in angiosperms. In Cactaceae (Caryophyllales), a high structural diversity of the chloroplast genome has been identified within and among genera. In this study, we assembled the first mitochondrial genome (mtDNA) for the short-globose cactus Mammillaria huitzilopochtli. For comparative purposes, we used the published genomes of 19 different angiosperms and the gymnosperm Cycas taitungensis as an external group for phylogenetic issues.

RESULTS

The mtDNA of M. huitzilopochtli was assembled into one linear chromosome of 2,052,004 bp, in which 65 genes were annotated. These genes account for 57,606 bp including 34 protein-coding genes (PCGs), 27 tRNAs, and three rRNAs. In the non-coding sequences, repeats were abundant, with a total of 4,550 (179,215 bp). In addition, five complete genes (psaC and four tRNAs) of chloroplast origin were documented. Negative selection was estimated for most (23) of the PCGs. The phylogenetic tree showed a topology consistent with previous analyses based on the chloroplast genome.

CONCLUSIONS

The number and type of genes contained in the mtDNA of M. huitzilopochtli were similar to those reported in 19 other angiosperm species, regardless of their phylogenetic relationships. Although other Caryophyllids exhibit strong differences in structural arrangement and total size of mtDNA, these differences do not result in an increase in the typical number and types of genes found in M. huitzilopochtli. We concluded that the total size of mtDNA in angiosperms increases by the lengthening of the non-coding sequences rather than a significant gain of coding genes.

Comparative chloroplast genomics and insights into the molecular evolution of Tanaecium (Bignonieae, Bignoniaceae)

Species of Tanaecium (Bignonieae, Bignoniaceae) are lianas distributed in the Neotropics and centered in the Amazon. Members of the genus exhibit exceptionally diverse flower morphology and pollination systems. Here, we sequenced, assembled, and annotated 12 complete and four partial chloroplast genomes representing 15 Tanaecium species and more than 70% of the known diversity in the genus. Gene content and order were similar in all species of Tanaecium studied, with genome sizes ranging between 158,470 and 160,935 bp. Tanaecium chloroplast genomes have 137 genes, including 80-81 protein-coding genes, 37 tRNA genes, and four rRNA genes. No rearrangements were found in Tanaecium plastomes, but two different patterns of boundaries between regions were recovered. Tanaecium plastomes show nucleotide variability, although only rpoA was hypervariable. Multiple SSRs and repeat regions were detected, and eight genes were found to have signatures of positive selection. Phylogeny reconstruction using 15 Tanaecium plastomes resulted in a strongly supported topology, elucidating several relationships not recovered previously and bringing new insights into the evolution of the genus.

Biosynthesis of Betalains Elicited by Methyl Jasmonate in Two Species of Alternanthera Genus: Antagonistic Regulations Result in Increase of Pigments

Natural pigments are components very important in the dye industry. The betalains are pigments found in plants from Caryophyllales order and are relevant in the food manufacturing. The main source of betalains is beetroot, which has unfavorable aftertaste. Therefore, the demand for alternative species producing betalains has increased. Elicitor molecules such as methyl jasmonate (MeJA) induce metabolic reprogramming acting in the biosynthesis of specialized metabolites and can enhance pigment concentrations. Here, we used this strategy to identify if treatment with MeJA at 100 µM can promote the accumulation of betalains and other bioactive compounds in Alternanthera philoxeroides and Alternanthera sessilis. We performed the gene expression, concentration of betalains, phenols, flavonoids, amino acids (phenylalanine and tyrosine), and antioxidant activity. The results showed that MeJA treatment increased betalains and other bioactive compounds in the two Alternanthera species but A. sessilis had a better performance. One key factor in this pathway is related to the phenylalanine and tyrosine concentration. However, the species have distinct metabolic regulation: in A. philoxeroides, high concentrations of betalain pigments increase the tyrosine concentration and gene expression (include ADH) under MeJA and in A. sessilis, high concentrations of betalain pigments reduce the gene expression and tyrosine concentration after 2 days under MeJA. This study brings new questions about betalain biosynthesis and sheds light on the evolution of this pathway in Caryophyllales.

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

A high-quality Bougainvillea genome provides new insights into evolutionary history and pigment biosynthetic pathways in the Caryophyllales

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana 'Mrs Butt', a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

The genome of the glasshouse plant noble rhubarb (Rheum nobile) provides a window into alpine adaptation

Glasshouse plants are species that trap warmth via specialized morphology and physiology, mimicking a human glasshouse. In the Himalayan alpine region, the highly specialized glasshouse morphology has independently evolved in distinct lineages to adapt to intensive UV radiation and low temperature. Here we demonstrate that the glasshouse structure - specialized cauline leaves - is highly effective in absorbing UV light but transmitting visible and infrared light, creating an optimal microclimate for the development of reproductive organs. We reveal that this glasshouse syndrome has evolved at least three times independently in the rhubarb genus Rheum. We report the genome sequence of the flagship glasshouse plant Rheum nobile and identify key genetic network modules in association with the morphological transition to specialized glasshouse leaves, including active secondary cell wall biogenesis, upregulated cuticular cutin biosynthesis, and suppression of photosynthesis and terpenoid biosynthesis. The distinct cell wall organization and cuticle development might be important for the specialized optical property of glasshouse leaves. We also find that the expansion of LTRs has likely played an important role in noble rhubarb adaptation to high elevation environments. Our study will enable additional comparative analyses to identify the genetic basis underlying the convergent occurrence of glasshouse syndrome.

Phylogenetic analysis and temporal diversification of the tribe Alsineae (Caryophyllaceae) with the description of three new genera, Hesperostellaria, Reniostellaria and Torreyostellaria

Alsineae are one of the most taxonomically difficult tribes in Caryophyllaceae and consist of over 500 species distributed in the northern temperate zone. Recent phylogenetic results have improved our understanding on the evolutionary relationships among Alsineae members. Nevertheless, there are still some unresolved taxonomic and phylogenetic problems at the generic level, and the evolutionary history of major clades within the tribe was unexplored to date. In this study, we carried out phylogenetic analyses and divergence time estimation of Alsineae using the nuclear ribosomal internal transcribed spacer (nrITS) and four plastid regions (matK, rbcL, rps16, trnL-F). The present analyses yielded a robustly supported phylogenetic hypothesis of the tribe. Our results showed that the monophyletic Alsineae are strongly supported to be the sister of Arenarieae, and the inter-generic relationships within Alsineae were mostly resolved with strong support. Both molecular phylogenetic and morphological evidence supported the Asian species Stellaria bistylata and the two North American species Pseudostellaria jamesiana and Stellaria americana all should be recognized as new monotypic genera respectively, and three new genera Reniostellaria, Torreyostellaria, and Hesperostellaria were thereby proposed here. Additionally, molecular and morphological evidence also supported the proposal of the new combination Schizotechium delavayi. Nineteen genera were accepted within Alsineae and a key to these genera was provided. Molecular dating analysis suggested that Alsineae splitted from its sister tribe at ca. 50.2 million-years ago (Ma) during the early Eocene and began to diverge at ca. 37.9 Ma during the late Eocene, and divergent events within Alsineae occurred mainly since the late Oligocene. Results from the present study provide insights into the historical assembly of herbaceous flora in northern temperate regions.

Characterization of the complete chloroplast genome of Gypsophila huashanensis Y. W. Tsui & D. Q. Lu, an endemic herb species in China

Gypsophila huashanensis Y. W. Tsui & D. Q. Lu (Caryophyllaceae) is an endemic herb species to the Qinling Mountains in China. In this study, we characterized its whole plastid genome using the Illumina sequencing platform. The complete plastid genome of G. huashanensis is 152,457 bp in length, including a large single-copy DNA region of 83,476 bp, a small single-copy DNA region of 17,345 bp, and a pair of inverted repeat DNA sequences of 25,818 bp. The genome contains 130 genes comprising 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Evolutionary analysis showed that the non-coding regions of Caryophyllaceae exhibit a higher level of divergence than the exon regions. Gene site selection analysis suggested that 11 coding protein genes (accD, atpF, ndhA, ndhB, petB, petD, rpoCl, rpoC2, rps16, ycfl, and ycf2) have some sites under protein sequence evolution. Phylogenetic analysis showed that G. huashanensis is most closely related to the congeneric species G. oldhamiana. These results are very useful for studying phylogenetic evolution and species divergence in the family Caryophyllaceae.

Characteristics of plastid genomes in the genus Ceratostigma inhabiting arid habitats in China and their phylogenomic implications

BACKGROUND

Ceratostigma, a genus in the Plumbaginaceae, is an ecologically dominant group of shrubs, subshrub and herb mainly distributed in Qinghai-Tibet Plateau and North China. Ceratostigma has been the focal group in several studies, owing to their importance in economic and ecological value and unique breeding styles. Despite this, the genome information is limited and interspecific relationships within the genus Cerotastigma remains unexplored. Here we sequenced, assembled and characterized the 14 plastomes of five species, and conducted phylogenetic analyses of Cerotastigma using plastomes and nuclear ribosomal DNA (nrDNA) data.

RESULTS

Fourteen Cerotastigma plastomes possess typical quadripartite structures with lengths from 164,076 to 168,355 bp that consist of a large single copy, a small single copy and a pair of inverted repeats, and contain 127-128 genes, including 82-83 protein coding genes, 37 transfer RNAs and eight ribosomal RNAs. All plastomes are highly conservative and similar in gene order, simple sequence repeats (SSRs), long repeat repeats and codon usage patterns, but some structural variations in the border of single copy and inverted repeats. Mutation hotspots in coding (Pi values > 0.01: matK, ycf3, rps11, rps3, rpl22 and ndhF) and non-coding regions (Pi values > 0.02: trnH-psbA, rps16-trnQ, ndhF-rpl32 and rpl32-trnL) were identified among plastid genomes that could be served as potential molecular markers for species delimitation and genetic variation studies in Cerotastigma. Gene selective pressure analysis showed that most protein-coding genes have been under purifying selection except two genes. Phylogenetic analyses based on whole plastomes and nrDNA strongly support that the five species formed a monophyletic clade. Moreover, interspecific delimitation was well resolved except C. minus, individuals of which clustered into two main clades corresponding to their geographic distributions. The topology inferred from the nrDNA dataset was not congruent with the tree derived from the analyses of the plastid dataset.

CONCLUSION

These findings represent the first important step in elucidating plastome evolution in this widespread distribution genus Cerotastigma in the Qinghai-Tibet Plateau. The detailed information could provide a valuable resource for understanding the molecular dynamics and phylogenetic relationship in the family Plumbaginaceae. Lineage genetic divergence within C. minus was perhaps promoted by geographic barriers in the Himalaya and Hengduan Mountains region, but introgression or hybridization could not be completely excluded.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

A plastid phylogenomic framework for the palm family (Arecaceae)

BACKGROUND

Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research.

RESULTS

Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support.

CONCLUSIONS

The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

Plastome variations reveal the distinct evolutionary scenarios of plastomes in the subfamily Cereoideae (Cactaceae)

BACKGROUND

The cactus family (Cactaceae) has been reported to have evolved a minimal photosynthetic plastome size, with the loss of inverted-repeat (IR) regions and NDH gene suites. However, there are very limited genomic data on the family, especially Cereoideae, the largest subfamily of cacti.

RESULTS

In the present study, we assembled and annotated 35 plastomes, 33 of which were representatives of Cereoideae, alongside 2 previously published plastomes. We analyzed the organelle genomes of 35 genera in the subfamily. These plastomes have variations rarely observed in those of other angiosperms, including size differences (with ~ 30 kb between the shortest and longest), dramatic dynamic changes in IR boundaries, frequent plastome inversions, and rearrangements. These results suggested that cacti have the most complex plastome evolution among angiosperms.

CONCLUSION

These results provide unique insight into the dynamic evolutionary history of Cereoideae plastomes and refine current knowledge of the relationships within the subfamily.

The molecular phylogenetic position of Harpagocarpus (Polygonaceae) sheds new light on the infrageneric classification of Fagopyrum

In the context of the molecular phylogeny of Polygonaceae, the phylogenetic positions of most genera and their relationships have been resolved. However, the monotypic genus Harpagocarpus has never been included in any published molecular phylogenetic studies. In the present study, we adopt a two-step approach to confirm the phylogenetic placement of Harpagocarpus using two datasets: (1) a concatenated dataset of three chloroplast DNA (cpDNA) regions (matK, rbcL and trnL-F) for Polygonaceae and (2) a combined cpDNA dataset of five sequences (accD, matK, psbA-trnH, rbcL and trnL-F) for Fagopyrum. Our analyses confirm the previous hypothesis based on morphological, anatomical and palynological investigations that Harpagocarpus is congeneric with Fagopyrum and further reveal that H.snowdenii (≡ F.snowdenii) is sister to the woody buckwheat F.tibeticum. Within Fagopyrum, three highly supported clades were discovered and the first sectional classification was proposed to accommodate them: sect. Fagopyrum comprises the two domesticated common buckwheat (F.esculentum and F.tataricum) and their wild relatives (F.esculentumsubsp.ancestrale, F.homotropicum and F.dibotrys) which are characterised by having large corymbose inflorescences and achenes greatly exceeding the perianth; sect. Tibeticum, including F.snowdenii and F.tibeticum, is characterised by the achene having appurtenances along the ribs, greatly exceeding the perianth and the perianth accrescent in fruit; sect. Urophyllum contains all other species of which the achenes were completely enclosed in the perianth. This study is very helpful to understand the phylogeny of the Fagopyrum and sheds light on the future study of taxonomy, biogeography, diversification and character evolution of the genus.

Genomic profiling of dioecious Amaranthus species provides novel insights into species relatedness and sex genes

BACKGROUND

Amaranthus L. is a diverse genus consisting of domesticated, weedy, and non-invasive species distributed around the world. Nine species are dioecious, of which Amaranthus palmeri S. Watson and Amaranthus tuberculatus (Moq.) J.D. Sauer are troublesome weeds of agronomic crops in the USA and elsewhere. Shallow relationships among the dioecious Amaranthus species and the conservation of candidate genes within previously identified A. palmeri and A. tuberculatus male-specific regions of the Y (MSYs) in other dioecious species are poorly understood. In this study, seven genomes of dioecious amaranths were obtained by paired-end short-read sequencing and combined with short reads of seventeen species in the family Amaranthaceae from NCBI database. The species were phylogenomically analyzed to understand their relatedness. Genome characteristics for the dioecious species were evaluated and coverage analysis was used to investigate the conservation of sequences within the MSY regions.

RESULTS

We provide genome size, heterozygosity, and ploidy level inference for seven newly sequenced dioecious Amaranthus species and two additional dioecious species from the NCBI database. We report a pattern of transposable element proliferation in the species, in which seven species had more Ty3 elements than copia elements while A. palmeri and A. watsonii had more copia elements than Ty3 elements, similar to the TE pattern in some monoecious amaranths. Using a Mash-based phylogenomic analysis, we accurately recovered taxonomic relationships among the dioecious Amaranthus species that were previously identified based on comparative morphology. Coverage analysis revealed eleven candidate gene models within the A. palmeri MSY region with male-enriched coverages, as well as regions on scaffold 19 with female-enriched coverage, based on A. watsonii read alignments. A previously reported FLOWERING LOCUS T (FT) within A. tuberculatus MSY contig was also found to exhibit male-enriched coverages for three species closely related to A. tuberculatus but not for A. watsonii reads. Additional characterization of the A. palmeri MSY region revealed that 78% of the region is made of repetitive elements, typical of a sex determination region with reduced recombination.

CONCLUSIONS

The results of this study further increase our understanding of the relationships among the dioecious species of the Amaranthus genus as well as revealed genes with potential roles in sex function in the species.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Comparison of Boraginales Plastomes: Insights into Codon Usage Bias, Adaptive Evolution, and Phylogenetic Relationships

The Boraginales (Boraginaceae a.l.) comprise more than 2450 species worldwide. However, little knowledge exists of the characteristics of the complete plastid genome. In this study, three new sequences representing the first pt genome of Heliotropiaceae and Cordiaceae were assembled and compared with other Boraginales species. The pt genome sizes of Cordia dichotoma, Heliotropium arborescens, and Tournefortia montana were 151,990 bp, 156,243 bp, and 155,891 bp, respectively. Multiple optimal codons were identified, which may provide meaningful information for enhancing the gene expression of Boraginales species. Furthermore, codon usage bias analyses revealed that natural selection and other factors may dominate codon usage patterns in the Boraginales species. The boundaries of the IR/LSC and IR/SSC regions were significantly different, and we also found a signal of obvious IR region expansion in the pt genome of Nonea vesicaria and Arnebia euchroma. Genes with high nucleic acid diversity (pi) values were also calculated, which may be used as potential DNA barcodes to investigate the phylogenetic relationships in Boraginales. psaI, rpl33, rpl36, and rps19 were found to be under positive selection, and these genes play an important role in our understanding of the adaptive evolution of the Boraginales species. Phylogenetic analyses implied that Boraginales can be divided into two groups. The existence of two subfamilies (Lithospermeae and Boragineae) in Boraginaceae is also strongly supported. Our study provides valuable information on pt genome evolution and phylogenetic relationships in the Boraginales species.

Plastid phylogenomics sheds light on divergence time and ecological adaptations of the tribe Persicarieae (Polygonaceae)

Southwestern China, adjacent to the Qinghai-Tibetan Plateau (QTP), is known as a hotspot for plant diversity and endemism, and it is the origin and diversification center of Persicarieae. As one of the major lineages in Polygonaceae, Persicarieae represents a diverse adaptation to various habitats. As a result of morphological plasticity and poorly resolving molecular markers, phylogenetic relationships and infrageneric classification within Persicarieae have long been controversial. In addition, neither plastome phylogenomic studies nor divergence time estimates on a larger sample of Persicarieae species have been made thus far. We sequenced and assembled 74 complete plastomes, including all of the recognized genera within Persicarieae and their relatives. We conducted a comprehensive phylogenetic study of the major clades within Persicarieae and, based on the thus obtained robust phylogeny, also estimated divergence time and the evolution of diagnostic morphological traits. Major relationships found in previous phylogenetic studies were confirmed, including those of the backbone of the tree, which had been a major problem in previous phylogenies of the tribe. Phylogenetic analysis revealed strong support for Koenigia as sister to Bistorta, and together they were sister to the robustly supported Persicaria. Based on the phylogenetic and morphological evidence, we recognize five sections in Persicaria: Persicaria, Amphibia, Tovara, Echinocaulon, and Cephalophilon. It is estimated that the divergence of the Persicarieae began around the late Paleocene, with diversification concentrated in the Eocene and Miocene. In addition, it is suggested that the increasing westerly and monsoon winds in conjunction with the uplift of the QTP may be the driving force for origin and diversification of Persicarieae species. These results provide a valuable evolutionary framework for the study of adaptation in Polygonaceae and insights into plant diversification on the QTP and adjacent areas.

Chromosome-scale Amaranthus tricolor genome provides insights into the evolution of the genus Amaranthus and the mechanism of betalain biosynthesis

Amaranthus tricolor is a vegetable and ornamental amaranth, with high lysine, dietary fibre and squalene content. The red cultivar of A. tricolor possesses a high concentration of betalains, which has been used as natural food colorants. Here, we constructed the genome of A. tricolor, the first reference genome for the subgenus Albersia, combining PacBio HiFi, Nanopore ultra-long and Hi-C data. The contig N50 size was 906 kb, and 99.58% of contig sequence was anchored to the 17 chromosomes, totalling 520 Mb. We annotated 27,813 protein-coding genes with an average 1.3 kb coding sequence and 5.3 exons. We inferred that A. tricolor underwent a whole-genome duplication (WGD) and that the WGD shared by amaranths occurred in the last common ancestor of subfamily Amaranthoideae. Moreover, we comprehensively identified candidate genes in betalain biosynthesis pathway. Among them, DODAα1 and CYP76ADα1, located in one topologically associated domain (TAD) of an active (A) compartment on chromosome 16, were more highly expressed in red leaves than in green leaves, and DODAα1 might be the rate-limiting enzyme gene in betalains biosynthesis. This study presents new genome resources and enriches our understanding of amaranth evolution, betalains production, facilitating molecular breeding improvements and the understanding of C4 plants evolution.

Evidence from Phylogenomics and Morphology Provide Insights into the Phylogeny, Plastome Evolution, and Taxonomy of Kitagawia

Kitagawia Pimenov is one of the segregate genera of Peucedanum sensu lato within the Apiaceae. The phylogenetic position and morphological delimitation of Kitagawia have been controversial. In this study, we used plastid genome (plastome) and nuclear ribosomal DNA (nrDNA) sequences to reconstruct the phylogeny of Kitagawia, along with comparative plastome and morphological analyses between Kitagawia and related taxa. The phylogenetic results identified that all examined Kitagawia species were divided into Subclade I and Subclade II within the tribe Selineae, and they were all distant from the representative members of Peucedanum sensu stricto. The plastomes of Kitagawia and related taxa showed visible differences in the LSC/IRa junction (JLA) and several hypervariable regions, which separated Subclade I and Subclade II from other taxa. Fruit anatomical and micromorphological characteristics, as well as general morphological characteristics, distinguished the four Kitagawia species within Subclade I from Subclade II and other related genera. This study supported the separation of Kitagawia from Peucedanum sensu lato, confirmed that Kitagawia belongs to Selineae, and two species (K. praeruptora and K. formosana) within Subclade II should be placed in a new genus. We believe that the "core" Kitagawia should be limited to Subclade I, and this genus can be distinguished by the association of a series of morphological characteristics. Overall, our study provides new insights into the phylogeny, plastome evolution, and taxonomy of Kitagawia.

A Comparative Phylogenetic Analysis on the Chloroplast Genome in Different Reynoutria japonica Populations

Reynoutria japonica Houtt., a traditional medicine herb of the Polygonaceae family, has been used since ancient times in China due to its various pharmacological effects. Chloroplast genomes are conservative and play an essential role in population diversity analysis. However, there are few studies on the chloroplast genome of R. japonica. In this study, the complete chloroplast genomes of three R. japonica from different regions were performed by next-generation sequencing technology. The results revealed that the lengths of the three chloroplast genomes are between 163,371~163,372 bp, and they have a highly conserved structure with a pair of inverted repeat (IR) regions (31,121 bp), a large single-copy (LSC) region (87,571~87,572 bp), and a small single-copy (SSC) region (13,558 bp). In total, 132 genes were annotated, including 8 rRNA genes, 37 tRNA genes, and 87 protein-coding genes. The phylogenetic analysis strongly revealed that 13 populations of R. japonica form a monophyly, and Fallopia multiflora (Polygonaceae) is its closest species. The two species diverged at ~20.47 million years ago, and R. japonica in China could be further divided into two major groups based on genetic structure analysis. In addition, several potential loci with suitable polymorphism were identified as molecular markers. Our study provides important genetic resources for further development and utilization of R. japonica germplasm, as well as some new insights into the evolutionary characteristics of this medicinal plant.

Plastome phylogenomics and biogeography of the subfam. Polygonoideae (Polygonaceae)

Polygonaceae has a complex taxonomic history, although a few studies using plastid or nuclear DNA fragments have explored relationships within this family, intrafamilial relationships remain controversial. Here, we newly sequenced and annotated 17 plastomes representing 12 genera within Polygonaceae. Combined with previously published data, a total of 49 plastomes representing 22/46 Polygonaceae genera and 16/20 Polygonoideae genera were collected to infer the phylogeny of Polygonaceae, with an emphasis on Polygonoideae. Plastome comparisons revealed high conservation within Polygonoideae in structure and gene order. Phylogenetic analyses using both Maximum Likelihood and Bayesian methods revealed two major clades and seven tribes within Polygonoideae. BEAST and S-DIVA analyses suggested a Paleocene origin of Polygonoideae in Asia. While most genera of Polygonoideae originated and further diversified in Asia, a few genera experienced multiple long-distance dispersal events from Eurasia to North America after the Miocene, with a few dispersal events to the Southern Hemisphere also being detected. Both ancient vicariance and long-distance events have played important roles in shaping the current distribution pattern of Polygonoideae.

Remarkably rapid, recent diversification of Cochemiea and Mammillaria in the Baja California, Mexico region

PREMISE

The Cactaceae of northwestern Mexico and the southwestern United States constitute a major component of the angiosperm biodiversity of the region. The Mammilloid clade, (Cactaceae, tribe Cacteae), composed of the genera Cochemiea, Coryphantha, Cumarinia, Mammillaria, and Pelecyphora is especially species rich. We sought to understand the timing, geographical and climate influences correlated with expansion of the Mammilloid clade, through the Sonoran Desert into Baja California.

METHODS

We reconstructed the historical biogeography of the Mammilloid clade, using Bayesian and maximum likelihood methods, based on a strongly supported molecular phylogeny. We also estimated divergence times, the timing of emergence of key characters, and diversification rates and rate shifts of the Mammilloid clade.

RESULTS

We found that the most recent common ancestor of Cochemiea arrived in the Cape region of Baja California from the Sonoran Desert region approximately 5 million years ago, coinciding with the timing of peninsular rifting from the mainland, suggesting dispersal and vicariance as causes of species richness and endemism. The diversification rate for Cochemiea is estimated to be approximately 12 times that of the mean background diversification rate for angiosperms. Divergence time estimation shows that many of the extant taxa in Cochemiea and Baja California Mammillaria emerged from common ancestors 1 million to 200,000 years ago, having a mid-Pleistocene origin.

CONCLUSIONS

Cochemiea and Mammillaria of the Baja California region are examples of recent, rapid diversification. Geological and climatic forces at multiple spatial and temporal scales are correlated with the western distributions of the Mammilloid clade.

Strong plastid degradation is consistent within section Chondrophyllae, the most speciose lineage of Gentiana

Recovering phylogenetic relationships in lineages experiencing intense diversification has always been a persistent challenge in evolutionary studies, including in Gentiana section Chondrophyllae sensu lato (s.l.). Indeed, this subcosmopolitan taxon encompasses more than 180 mostly annual species distributed around the world. We sequenced and assembled 22 new plastomes representing 21 species in section Chondrophyllae s.l. In addition to previously released plastome data, our study includes all main lineages within the section. We reconstructed their phylogenetic relationships based on protein-coding genes and recombinant DNA (rDNA) cistron sequences, and then investigated plastome structural evolution as well as divergence time. Despite an admittedly humble species cover overall, we recovered a well-supported phylogenetic tree based on plastome data, and found significant discordance between phylogenetic relationships and taxonomic treatments. Our results show that G. capitata and G. leucomelaena diverged early within the section, which is then further divided into two clades. The divergence time estimation showed that section Chondrophyllae s.l. evolved in the second half of the Oligocene. We found that section Chondrophyllae s.l. had the smallest average plastome size (128 KB) in tribe Gentianeae (Gentianaceae), with frequent gene and sequence losses such as the ndh complex and its flanking regions. In addition, we detected both expansion and contraction of the inverted repeat (IR) regions. Our study suggests that plastome degradation parallels the diversification of this group, and illustrates the strong discordance between phylogenetic relationships and taxonomic treatments, which now need to be carefully revised.

Plastome comparison and phylogenomics of Fagopyrum (Polygonaceae): insights into sequence differences between Fagopyrum and its related taxa

BACKGROUND

Fagopyrum (Polygonaceae) is a small plant lineage comprised of more than fifteen economically and medicinally important species. However, the phylogenetic relationships of the genus are not well explored, and the characteristics of Fagopyrum chloroplast genomes (plastomes) remain poorly understood so far. It restricts the comprehension of species diversity in Fagopyrum. Therefore, a comparative plastome analysis and comprehensive phylogenomic analyses are required to reveal the taxonomic relationship among species of Fagopyrum.

RESULTS

In the current study, 12 plastomes were sequenced and assembled from eight species and two varieties of Fagopyrum. In the comparative analysis and phylogenetic analysis, eight previously published plastomes of Fagopyrum were also included. A total of 49 plastomes of other genera in Polygonaceae were retrieved from GenBank and used for comparative analysis with Fagopyrum. The variation of the Fagopyrum plastomes is mainly reflected in the size and boundaries of inverted repeat/single copy (IR/SC) regions. Fagopyrum is a relatively basal taxon in the phylogenomic framework of Polygonaceae comprising a relatively smaller plastome size (158,768-159,985 bp) than another genus of Polygonaceae (158,851-170,232 bp). A few genera of Polygonaceae have nested distribution of the IR/SC boundary variations. Although most species of Fagopyrum show the same IRb/SC boundary with species of Polygonaceae, only a few species show different IRa/SC boundaries. The phylogenomic analyses of Fagopyrum supported the cymosum and urophyllum groups and resolved the systematic position of subclades within the urophyllum group. Moreover, the repeat sequence types and numbers were found different between groups of Fagopyrum. The plastome sequence identity showed significant differences between intra-group and inter-group.

CONCLUSIONS

The deletions of intergenic regions cause a short length of Fagopyrum plastomes, which may be the main reason for plastome size diversity in Polygonaceae species. The phylogenomic reconstruction combined with the characteristics comparison of plastomes supports grouping within Fagopyrum. The outcome of these genome resources may facilitate the taxonomy, germplasm resources identification as well as plant breeding of Fagopyrum.

Phylogeny and Taxonomic Synopsis of the Genus Bougainvillea (Nyctaginaceae)

Bougainvillea Comm. ex Juss. is one of the renowned genera in the Nyctaginaceae, but despite its recognized horticultural value, the taxonomy and phylogeny of the genus is not well-studied. Phylogenetic reconstructions based on plastid genomes showed that B. pachyphylla and B. peruviana are basal taxa, while B. spinosa is sister to two distinct clades: the predominantly cultivated Bougainvillea clade (B. spectabilis, B. glabra, B. arborea, B. cultivar, B. praecox) and the clade containing wild species of Bougainvillea (B. berberidifolia, B. campanulata, B. infesta, B. modesta, B. luteoalba, B. stipitata, and B. stipitata var. grisebachiana). Early divergence of B. peruviana, B. pachyphylla and B. spinosa is highly supported, thus the previously proposed division of Bougainvillea into two subgenera (Bougainvillea and Tricycla) was not reflected in this study. Morphological analysis also revealed that leaf arrangement, size, and indumentum together with the perianth tube and anthocarp shape and indumentum are important characteristics in differentiating the species of Bougainvillea. In the present study, 11 species and one variety are recognized in Bougainvillea. Six names are newly reduced to synonymy, and lectotypes are designated for 27 names. In addition, a revised identification key and illustrations of the distinguishing parts are also provided in the paper.

Phylogenomics Reveals the Evolutionary History of Phytolacca (Phytolaccaceae)

Phytolacca is the largest genus of Phytolaccaceae. Owing to interspecific hybridization, infraspecific variation, and apparent weak genetic control of many qualitative characters, which have obscured boundaries between species, the classification and phylogenetic relationships of this genus are unclear. Native Phytolacca is disjunctly distributed in America, eastern Asia, and Africa, and the biogeographic history of the genus remained unresolved. In this study, we used the whole chloroplast genome and three markers (nrDNA, rbcL, and matK) to reconstruct phylogenetic relationships within Phytolacca, analyze divergence times, and infer biogeographic histories. The phylogenetic results indicate that Phytolacca is monophyletic, which is inconsistent with the infrageneric classification based on morphology. According to the divergence time estimation, Phytolacca began to diversify at approximately 20.30 Ma during the early Miocene. Central America, including Mexico, Costa Rica, and Colombia, is the center of species diversity. Biogeographical analysis indicated five main dispersal events and Phytolacca originated from Central and South America. Birds may be the primary agents of dispersal because of the fleshy fruiting of Phytolacca. This study extended sampling and added more genetic characteristics to infer the evolutionary history of Phytolacca, providing new insights for resolving the classification and elucidating the dispersal events of Phytolacca.

Deep Insights Into the Plastome Evolution and Phylogenetic Relationships of the Tribe Urticeae (Family Urticaceae)

Urticeae s.l., a tribe of Urticaceae well-known for their stinging trichomes, consists of more than 10 genera and approximately 220 species. Relationships within this tribe remain poorly known due to the limited molecular and taxonomic sampling in previous studies, and chloroplast genome (CP genome/plastome) evolution is still largely unaddressed. To address these concerns, we used genome skimming data-CP genome and nuclear ribosomal DNA (18S-ITS1-5.8S-ITS2-26S); 106 accessions-for the very first time to attempt resolving the recalcitrant relationships and to explore chloroplast structural evolution across the group. Furthermore, we assembled a taxon rich two-locus dataset of trnL-F spacer and ITS sequences across 291 accessions to complement our genome skimming dataset. We found that Urticeae plastomes exhibit the tetrad structure typical of angiosperms, with sizes ranging from 145 to 161 kb and encoding a set of 110-112 unique genes. The studied plastomes have also undergone several structural variations, including inverted repeat (IR) expansions and contractions, inversion of the trnN-GUU gene, losses of the rps19 gene, and the rpl2 intron, and the proliferation of multiple repeat types; 11 hypervariable regions were also identified. Our phylogenomic analyses largely resolved major relationships across tribe Urticeae, supporting the monophyly of the tribe and most of its genera except for Laportea, Urera, and Urtica, which were recovered as polyphyletic with strong support. Our analyses also resolved with strong support several previously contentious branches: (1) Girardinia as a sister to the Dendrocnide-Discocnide-Laportea-Nanocnide-Zhengyia-Urtica-Hesperocnide clade and (2) Poikilospermum as sister to the recently transcribed Urera sensu stricto. Analyses of the taxon-rich, two-locus dataset showed lower support but was largely congruent with results from the CP genome and nuclear ribosomal DNA dataset. Collectively, our study highlights the power of genome skimming data to ameliorate phylogenetic resolution and provides new insights into phylogenetic relationships and chloroplast structural evolution in Urticeae.

Complete Chloroplast Genome of Cnidium monnieri (Apiaceae) and Comparisons with Other Tribe Selineae Species

Cnidium monnieri is an economically important traditional Chinese medicinal plant. In this study, the complete chloroplast (cp) genome of C. monnieri was determined using the Illumina paired-end sequencing, the GetOrganelle de novo assembly strategy, as well as the GeSeq annotation method. Our results showed that the cp genome was 147,371 bp in length with 37.4% GC content and included a large single-copy region (94,361 bp) and a small single-copy region (17,552 bp) separated by a pair of inverted repeat regions (17,729 bp). A total of 129 genes were contained in the cp genome, including 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. We also investigated codon usage, RNA editing, repeat sequences, simple sequence repeats (SSRs), IR boundaries, and pairwise Ka/Ks ratios. Four hypervariable regions (trnD-trnY-trnE-trnT, ycf2, ndhF-rpl32-trnL, and ycf1) were identified as candidate molecular markers for species authentication. The phylogenetic analyses supported non-monophyly of Cnidium and C. monnieri located in tribe Selineae based on the cp genome sequences and internal transcribed spacer (ITS) sequences. The incongruence of the phylogenetic position of C. monnieri between ITS and cpDNA phylogenies suggested that C. monnieri might have experienced complex evolutions with hybrid and incomplete lineage sorting. All in all, the results presented herein will provide plentiful chloroplast genomic resources for studies of the taxonomy, phylogeny, and species authentication of C. monnieri. Our study is also conducive to elucidating the phylogenetic relationships and taxonomic position of Cnidium.

Characterization of the complete chloroplast genome of Rumex dentatus L. (Polygonaceae)

Rumex dentatus L. is a flowering plant with promising therapeutic effects. Here, we reported the complete chloroplast genome sequence of R. dentatus. The length of the complete genome is 159,103 bp, with a pair of inverted repeat regions (IRb/IRa, 30,628 bp) separated by a large single-copy (LSC, 84,848 bp) region and a small single-copy (SSC, 12,999 bp) region. The overall GC content of the genome was 37.6%, and the corresponding values in the LSC, SSC, and IR regions were 35.8%, 32.6%, and 41.1%, respectively. There are 129 genes, including 84 protein coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic tree showed that R. dentatus was most relative to the species of R. nepalensis and R. crispus.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives

Background

As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes – traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution.

Results

We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species.

Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated.

Conclusion

In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08336-8.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.