Comparative analyses of plastid genomes from fourteen Cornales species: inferences for phylogenetic relationships and genome evolution

Complete Organelle Genome of the Desiccation-Tolerant (DT) Moss Tortula atrovirens and Comparative Analysis of the Pottiaceae Family

Tortula atrovirens (Sm.) Lindb. is an important component of biological soil crusts and possesses an extraordinary tolerance against desiccation in dryland habitats. However, knowledge of the organelle genome of this desiccation-tolerant (DT) moss is still lacking. Here, we assembled the first reported Tortula organelle genome and conducted a comprehensive analysis within the Pottiaceae family. T. atrovirens exhibited the second largest chloroplast genome (129,646 bp) within the Pottiaceae, whereas its mitogenome (105,877 bp) and those of other mosses were smaller in size compared to other land plants. The chloroplast and mitochondrial genomes of T. atrovirens were characterized by the expansion of IR boundaries and the absence of homologous recombination-mediated by large repeats. A total of 57 RNA editing sites were detected through mapping RNA-seq data. Moreover, the gene content and order were highly conserved among the Pottiaceae organelle genomes. Phylogenetic analysis showed that bryophytes are paraphyletic, with their three lineages (hornworts, mosses, and liverworts) and vascular plants forming successive sister clades. Timmiella anomala is clearly separated from the monophyletic Pottiaceae, and T. atrovirens is closely related to Syntrichia filaris within the Pottioideae. In addition, we detected four hypervariable regions for candidate-molecular markers. Our findings provide valuable insights into the organelle genomes of T. atrovirens and the evolutionary relationships within the Pottiaceae family, facilitating future discovery of DT genetic resources from bryophytes.

Characteristics and phylogenetic analysis of the complete chloroplast genome of Cornus hongkongensis Hemsl. 1888 (Cornaceae)

Cornus hongkongensis Hemsl. 1888, native to Hong Kong, belongs to the subgenus Syncarpea within the Cornus genus of the Cornaceae family. The complete chloroplast genome of C. hongkongensis spans 156,954 bp, comprising four subregions: a large single-copy region (86,290 bp), a small single-copy region (18,394 bp), and a pair of inverted repeats (26,135 bp). Within the chloroplast genome of C. hongkongensis, we identified 113 unique genes, including 80 protein-encoding genes, four ribosomal RNA (rRNA) genes, and 30 transfer RNA (tRNA) genes. Phylogenetic analysis based on the complete chloroplast genome of 30 related taxa of the Cornus genus indicates that C. hongkongensis has not formed a monophyletic lineage. Analyses of sequence divergence found three intergenic regions including rps19-rpl22, ccsA-ndhD, and atpH-atpI, exhibiting a high degree of variations. The first chloroplast genome of C. hongkongensis was reported in this work contributes to the enrichment of genomic data for the genus Cornus.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae)

BACKGROUND AND AIMS

Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes.

METHODS

To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes).

KEY RESULTS

Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic.

CONCLUSIONS

This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

First record of the genus Camptotheca (Nyssaceae) in Vietnam and the lectotypification of C.acuminata

As a primary source of anticancer camptothecin, Camptotheca (Nyssaceae) is an economically valuable genus and has long been recorded as endemic to China. Here, Camptotheca is reported as a new record to the flora of Vietnam with the discovery of a wild population of C.acuminata from Lai Chau Province of northern Vietnam. Based on the consultation of relevant literature and type specimens of C.acuminata, a lectotype of the species is designated. Photographic illustrations, morphological description and a distribution map of C.acuminata is provided, and a key to all known species of Nyssaceae in Vietnam is presented, too.

Extensive characterization of 28 complete chloroplast genomes of Hydrangea species: A perspective view of their organization and phylogenetic and evolutionary relationships

The tribe Hydrangeeae displays a unique, distinctive disjunct distribution encompassing East Asia, North America and Hawaii. Despite its complex trait variations and polyphyletic nature, comprehensive phylogenomic and biogeographical studies on this tribe have been lacking. To address this gap, we sequenced and characterized 28 plastomes of Hydrangeeae. Our study highlights the highly conserved nature of Hydrangeaceae chloroplast (cp) genomes in terms of gene content and arrangement. Notably, synapomorphic characteristics of tandem repeats in the conserved domain of accD were observed in the Macrophyllae, Chinenses, and Dichroa sections within the Hydrangeeae tribe. Additionally, we found lower expression of accD in these sections using structure prediction and quantitative real-time PCR analysis. Phylogenomic analyses revealed the subdivision of the Hydrangeeae tribe into two clades with robust support values. Consistent with polyphyletic relationships, sect. Broussaisia was identified as the basal group in the tribe Hydrangeeae. Our study also provides insights into the phylogenetic relationships of Hydrangea petiolaris in the Jeju and Ulleung Island populations, suggesting the need for further studies with more samples and molecular data. Divergence time estimation and biogeographical analyses suggested that the common ancestors of the tribe Hydrangeeae likely originated from North America and East Asia during the Paleocene period via the Bering Land Bridge, potentially facilitating migration within the tribe between these regions. In conclusion, this study enhances our understanding of the evolutionary history and biogeography of the tribe Hydrangeeae, shedding light on the dispersal patterns and origins of this intriguing plant group with its unique disjunct distribution.

The complete plastid genome and characteristics analysis of Achillea millefolium

Achillea is a crop with Chinese herbal characteristics and horticultural values. Its leaves and flowers contain aromatic oil, and the ripe herb can also be used as medicine to induce sweat and relieve rheumatic pains. It is seen cultivated in gardens all over China. Currently, the most comprehensive chloroplast genome sample involved in the study refers to New World clades of Achillea, which are used for marker selection and phylogenetic research. We completely sequenced the chloroplast genomes of Achillea millefolium. These sequencing results showed that the plastid genome is 149,078 bp in size and possesses a typical quadripartite structure containing one large single copy (LSC) with 82,352 bp, one small single copy (SSC) with 18,426 bp, and a pair of inverted repeat (IR) regions with 24,150 bp in Achillea millefolium. The chloroplast genome encodes a common number of genes, of which 88 are protein-coding genes, 37 transfer ribonucleic acid genes, and 8 ribosomal ribonucleic acid genes, which are highly similar in overall size, genome structure, gene content, and sequence. The exact similarity was observed when compared to other Asteraceae species. However, there were structural differences due to the restriction or extension of the inverted repeat (IR) regions-the palindromic repeats being the most prevalent form. Based on 12 whole-plastomes, 3 hypervariable regions (rpoB, rbcL, and petL-trnP-UGG) were discovered, which could be used as potential molecular markers.

Phylogenomics and phylogeography of Menispermum (Menispermaceae)

INTRODUCTION

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

METHODS

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

RESULTS AND DISCUSSION

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

Plastid Phylogenomics and Plastome Evolution of Nandinoideae (Berberidaceae)

Subfamily Nandinoideae Heintze (Berberidaceae), comprising four genera and ca. 19 species, is disjunctively distributed in eastern North America vs. Eurasia (eastern Asia, Central Asia, Middle East, and southeastern Europe), and represents an ideal taxon to explore plastid phylogenomics and plastome evolution in Berberidaceae. Many species of this subfamily have been listed as national or international rare and endangered plants. In this study, we sequenced and assembled 20 complete plastomes, representing three genera and 13 species of Nandinoideae. Together with six plastomes from GenBank, a total of 26 plastomes, representing all four genera and 16 species of Nandinoideae, were used for comparative genomic and phylogenomic analyses. These plastomes showed significant differences in overall size (156,626-161,406 bp), which is mainly due to the expansion in inverted repeat (IR) regions and/or insertion/deletion (indel) events in intergenic spacer (IGS) regions. A 75-bp deletion in the ndhF gene occurred in Leontice and Gymnospermium when compared with Nandina and Caulophyllum. We found a severe truncation at the 5' end of ycf1 in three G. altaicum plastomes, and a premature termination of ropC1 in G. microrrhynchum. Our phylogenomic results support the topology of {Nandina, [Caulophyllum, (Leontice, Gymnospermium)]}. Within the core genus Gymnospermium, we identified G. microrrhynchum from northeastern Asia (Clade A) as the earliest diverging species, followed by G. kiangnanense from eastern China (Clade B), while the rest species clustered into the two sister clades (C and D). Clade C included three species from West Tianshan (G. albertii, G. darwasicum, G. vitellinum). Clade D consisted of G. altaicum from northern Central Asia, plus one species from the Caucasus Mountains (G. smirnovii) and three from southeastern Europe (G. odessanum, G. peloponnesiacum, G. scipetarum). Overall, we identified 21 highly variable plastome regions, including two coding genes (rpl22, ycf1) and 19 intergenic spacer (IGS) regions, all with nucleotide diversity (Pi) values > 0.02. These molecular markers should serve as powerful tools (including DNA barcodes) for future phylogenetic, phylogeographic and conservation genetic studies.

The chloroplast genomes of four Bupleurum (Apiaceae) species endemic to Southwestern China, a diversity center of the genus, as well as their evolutionary implications and phylogenetic inferences

Background

As one of the largest genera in Apiaceae, Bupleurum L. is well known for its high medicinal value. The genus has frequently attracted the attention of evolutionary biologist and taxonomist for its distinctive characteristics in the Apiaceae family. Although some chloroplast genomes data have been now available, the changes in the structure of chloroplast genomes and selective pressure in the genus have not been fully understood. In addition, few of the species are endemic to Southwest China, a distribution and diversity center of Chinese Bupleurum. Endemic species are key components of biodiversity and ecosystems, and investigation of the chloroplast genomes features of endemic species in Bupleurum will be helpful to develop a better understanding of evolutionary process and phylogeny of the genus. In this study, we analyzed the sequences of whole chloroplast genomes of 4 Southwest China endemic Bupleurum species in comparison with the published data of 17 Bupleurum species to determine the evolutionary characteristics of the genus and the phylogenetic relationships of Asian Bupleurum.

Results

The complete chloroplast genome sequences of the 4 endemic Bupleurum species are 155,025 bp to 155,323 bp in length including a SSC and a LSC region separated by a pair of IRs. Comparative analysis revealed an identical chloroplast gene content across the 21 Bupleurum species, including a total of 114 unique genes (30 tRNA genes, 4 rRNA genes and 80 protein-coding genes). Chloroplast genomes of the 21 Bupleurum species showed no rearrangements and a high sequence identity (96.4–99.2%). They also shared a similar tendency of SDRs and SSRs, but differed in number (59–83). In spite of their high conservation, they contained some mutational hotspots, which can be potentially exploited as high-resolution DNA barcodes for species discrimination. Selective pressure analysis showed that four genes were under positive selection. Phylogenetic analysis revealed that the 21 Bupleurum formed two major clades, which are likely to correspond to their geographical distribution.

Conclusions

The chloroplast genome data of the four endemic Bupleurum species provide important insights into the characteristics and evolution of chloroplast genomes of this genu, and the phylogeny of Bupleurum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08008-z.

Comparative Chloroplast Genomics and Phylogenetic Analysis of Zygophyllum (Zygophyllaceae) of China

The genus Zygophyllum comprises over 150 species within the plant family Zygophyllaceae. These species predominantly grow in arid and semiarid areas, and about 20 occur in northwestern China. In this study, we sampled 24 individuals of Zygophyllum representing 15 species and sequenced their complete chloroplast (cp) genomes. For comparison, we also sequenced cp genomes of two species of Peganum from China representing the closely allied family, Nitrariaceae. The 24 cp genomes of Zygophyllum were smaller and ranged in size from 104,221 to 106,286 bp, each containing a large single-copy (LSC) region (79,245-80,439 bp), a small single-copy (SSC) region (16,285-17,146 bp), and a pair of inverted repeat (IR) regions (3,792-4,466 bp). These cp genomes contained 111-112 genes each, including 74-75 protein-coding genes (PCGs), four ribosomal RNA genes, and 33 transfer RNA genes, and all cp genomes showed similar gene order, content, and structure. The cp genomes of Zygophyllum appeared to lose some genes such as ndh genes and rRNA genes, of which four rRNA genes were in the SSC region, not in the IR regions. However, the SC and IR regions had greater similarity within Zygophyllum than between the genus and Peganum. We detected nine highly variable intergenic spacers: matK-trnQ, psaC-rps15, psbZ-trnG, rps7-trnL, rps15-trnN, trnE-trnT, trnL-rpl32, trnQ-psbK, and trnS-trnG. Additionally, we identified 156 simple sequence repeat (cpSSR) markers shared among the genomes of the 24 Zygophyllum samples and seven cpSSRs that were unique to the species of Zygophyllum. These markers may be useful in future studies on genetic diversity and relationships of Zygophyllum and closely related taxa. Using the sequenced cp genomes, we reconstructed a phylogeny that strongly supported the division of Chinese Zygophyllum into herbaceous and shrubby clades. We utilized our phylogenetic results along with prior morphological studies to address several remaining taxonomic questions within Zygophyllum. Specifically, we found that Zygophyllum kaschgaricum is included within Zygophyllum xanthoxylon supporting the present treatment of the former genus Sarcozygium as a subgenus within Zygophyllum. Our results provide a foundation for future research on the genetic resources of Zygophyllum.

Comparative genome analysis revealed gene inversions, boundary expansions and contractions, and gene loss in the Stemona sessilifolia (Miq.) Miq. chloroplast genome

Stemona sessilifolia (Miq.) Miq., commonly known as Baibu, is one of the most popular herbal medicines in Asia. In the Chinese Pharmacopoeia, Baibu has multiple authentic sources and there are many similar herbs sold as Baibu in herbal medicine markets. The existence of counterfeits of Baibu brings challenges to its identification. To assist in its accurate identification, we sequenced and analyzed the complete chloroplast genome of S. sessilifolia using next-generation sequencing technology. The genome was found to be 154,037 bp in length, possessing a typical quadripartite structure consisting of a pair of inverted repeats (IRs: 27,090 bp) separated by a large single copy (LSC: 81,949 bp) and a small single copy (SSC: 17,908 bp). A total of 112 unique genes were identified, including 80 protein-coding, 28 transfer RNA and four ribosomal RNA genes. In addition, 45 tandem, 27 forward, 23 palindromic and 104 simple sequence repeats were detected in the genome by repeated analysis. Compared with its counterfeits (Asparagus officinalis and Carludovica palmata) we found that IR expansion and SSC contraction events of S. sessilifolia resulted in two copies of the rpl22 gene in the IR regions and a partial duplication of the ndhF gene in the SSC region. An approximately 3-kb-long inversion was also identified in the LSC region, leading to the petA and cemA genes being presented in the complementary strand of the chloroplast DNA molecule. Comparative analysis revealed some highly variable regions, including trnF-GAA_ndhJ, atpB_rbcL, rps15_ycf1, trnG-UCC_trnR-UCU, ndhF_rpl32, accD_psaI, rps2_rpoC2, trnS-GCU_trnG-UCC, trnT-UGU_trnL-UAA and rps16_trnQ-UUG. Finally, gene loss events were investigated in the context of phylogenetic relationships. In summary, the complete plastome of S. sessilifolia will provide valuable information for the distinction between Baibu and its counterfeits and assist in elucidating the evolution of S. sessilifolia.

Targeted enrichment of novel chloroplast-based probes reveals a large-scale phylogeny of 412 bamboos

BACKGROUND

The subfamily Bambusoideae belongs to the grass family Poaceae and has significant roles in culture, economy, and ecology. However, the phylogenetic relationships based on large-scale chloroplast genomes (CpGenomes) were elusive. Moreover, most of the chloroplast DNA sequencing methods cannot meet the requirements of large-scale CpGenome sequencing, which greatly limits and impedes the in-depth research of plant genetics and evolution.

RESULTS

To develop a set of bamboo probes, we used 99 high-quality CpGenomes with 6 bamboo CpGenomes as representative species for the probe design, and assembled 15 M unique sequences as the final pan-chloroplast genome. A total of 180,519 probes for chloroplast DNA fragments were designed and synthesized by a novel hybridization-based targeted enrichment approach. Another 468 CpGenomes were selected as test data to verify the quality of the newly synthesized probes and the efficiency of the probes for chloroplast capture. We then successfully applied the probes to synthesize, enrich, and assemble 358 non-redundant CpGenomes of woody bamboo in China. Evaluation analysis showed the probes may be applicable to chloroplasts in Magnoliales, Pinales, Poales et al. Moreover, we reconstructed a phylogenetic tree of 412 bamboos (358 in-house and 54 published), supporting a non-monophyletic lineage of the genus Phyllostachys. Additionally, we shared our data by uploading a dataset of bamboo CpGenome into CNGB ( https://db.cngb.org/search/project/CNP0000502/ ) to enrich resources and promote the development of bamboo phylogenetics.

CONCLUSIONS

The development of the CpGenome enrichment pipeline and its performance on bamboos recommended an inexpensive, high-throughput, time-saving and efficient CpGenome sequencing strategy, which can be applied to facilitate the phylogenetics analysis of most green plants.

Chloroplast Genome Sequences and Comparative Analyses of Combretaceae Mangroves with Related Species

In the Combretaceae family, only two species of Lumnitzera and one species of Laguncularia belong to mangroves. Among them, Lumnitzera littorea (Jack) Voigt. is an endangered mangrove plant in China for the limited occurrence and seed abortion. In contrast, Lumnitzera racemosa Willd. is known as the most widespread mangrove plant in China. Laguncularia racemosa C. F. Gaertn., an exotic mangrove in China, has the fast growth and high adaptation ability. To better understand the phylogenetic positions of these mangroves in Combretaceae and in Myrtales and to provide information for studies on evolutionary adaptation for intertidal habitat, the complete chloroplast (cp) genomes of Lu. racemosa and La. racemosa were sequenced. Furthermore, we present here the results from the assembly and annotation of the two cp genomes, which were further subjected to the comparative analysis with Lu. littorea cp genomes we published before and other eleven closely related species within Myrtales. The chloroplast genomes of the three Combretaceae mangrove species: Lu. littorea, Lu. racemosa, and La. racemosa are 159,687 bp, 159,473 bp, and 158,311 bp in size. All three cp genomes host 130 genes including 85 protein-coding genes, 37 tRNAs, and 4 rRNAs. A comparative analysis of those three genomes revealed the high similarity of genes in coding-regions and conserved gene order in the IR and LSC/SSC regions. The differences between Lumnitzera and Laguncularia cp genomes are the locations of rps19 and rpl2 genes in the IR/SC boundary regions. Investigating the effects of selection events on shared protein-coding genes showed a relaxed selection had acted on the ycf2, ycf1, and matK genes of Combretaceae mangroves compared to the nonmangrove species Eucalyptus aromaphloia. The phylogenetic analysis based on the whole chloroplast genome sequence with one outgroup species strongly supported three Combretaceae mangroves together with other two Combretaceae species formed a cluster in Combretaceae. This study is the first report on the comparative analysis of three Combretaceae mangrove chloroplast genomes, which will provide the significant information for understanding photosynthesis and evolution in Combretaceae mangrove plants.

The large single-copy (LSC) region functions as a highly effective and efficient molecular marker for accurate authentication of medicinal Dendrobium species

Having great medicinal values, Dendrobium species of "Fengdou" (DSFs) are a taxonomically complex group in Dendrobium genus including many closely related and recently diverged species. Traditionally used DNA markers have been proved to be insufficient in authenticating many species of this group. Here, we investigated 101 complete plastomes from 23 DSFs, comprising 72 newly sequenced and 29 documented, which all exhibited well-conserved genomic organization and gene order. Plastome-wide comparison showed the co-occurrence of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels), which can be explained by both the repeat-associated and indel-associated mutation hypotheses. Moreover, guanine-cytosine (GC) content was found to be negatively correlated with the three divergence variables (SNPs, indels and repeats), indicating that GC content may reflect the level of the local sequence divergence. Our species authentication analyses revealed that the relaxed filtering strategies of sequence alignment had no negative impact on species identification. By assessing the maximum likelihood (ML) trees inferred from different datasets, we found that the complete plastome and large single-copy (LSC) datasets both successfully identified all 23 DSFs with the maximum bootstrap values. However, owing to the high efficiency of LSC in species identification, we recommend using LSC for accurate authentication of DSFs.

Plastome comparative genomics in maples resolves the infrageneric backbone relationships

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deep-level relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Leveraging plastomes for comparative analysis and phylogenomic inference within Scutellarioideae (Lamiaceae)

Scutellaria, or skullcaps, are medicinally important herbs in China, India, Japan, and elsewhere. Though Scutellaria is the second largest and one of the more taxonomically challenging genera within Lamiaceae, few molecular systematic studies have been undertaken within the genus; in part due to a paucity of available informative markers. The lack of informative molecular markers for Scutellaria hinders our ability to accurately and robustly reconstruct phylogenetic relationships, which hampers our understanding of the diversity, phylogeny, and evolutionary history of this cosmopolitan genus. Comparative analyses of 15 plastomes, representing 14 species of subfamily Scutellarioideae, indicate that plastomes within Scutellarioideae contain about 151,000 nucleotides, and possess a typical quadripartite structure. In total, 590 simple sequence repeats, 489 longer repeats, and 16 hyper-variable regions were identified from the 15 plastomes. Phylogenetic relationships among the 14 species representing four of the five genera of Scutellarioideae were resolved with high support values, but the current infrageneric classification of Scutellaria was not supported in all analyses. Complete plastome sequences provide better resolution at an interspecific level than using few to several plastid markers in phylogenetic reconstruction. The data presented here will serve as a foundation to facilitate DNA barcoding, species identification, and systematic research within Scutellaria, which is an important medicinal plant resource worldwide.

Comparative analyses of chloroplast genomes from 13 Lagerstroemia (Lythraceae) species: identification of highly divergent regions and inference of phylogenetic relationships

Seven divergence hotspots as plastid markers for DNA barcoding was selected, and the phylogeny of 13 Lagerstroemia species based on the cp genome data was reconstructed within Myrtales. The Lagerstroemia species used in this study originated in China and have high economic and ecological value. The shared interspecific morphological characteristics and intraspecific morphological variation resulting from hybridization among Lagerstroemia taxa have made resolving their classification problems and phylogenetic relationships difficult. Systematic comparative genomic analysis has been shown to resolve phylogenetic relationships. We sequenced and annotated 6 Lagerstroemia cp genomes (Lagerstroemia excelsa, Lagerstroemia limii, Lagerstroemia siamica, Lagerstroemia tomentosa, Lagerstroemia venusta, and Lagerstroemia calyculata) for the first time and combined them with previously published genomes for Lagerstroemia species. Bioinformatics was used to analyse the 13 cp genomes in terms of gene structure and organization, codon usage, contraction and expansion of inverted repeat regions, repeat structure, divergence hotspots, species pairwise Ka/Ks ratios and phylogenetic relationships. The length varied between 152,049 bp in Lagerstroemia subcostata and 152,521 bp in L. venusta. We selected seven divergence hotspots in the cp genomes that had the potential to act as plastid markers to distinguish Lagerstroemia species. The phylogenetic relationships within Myrtales inferred from the cp genomes of 13 Lagerstroemia species and 27 other Myrtales species were highly supported, which illustrated several novel relationships within Myrtales. Taken together, our results provide comprehensive chloroplast genomic resources, which can be used further for species identification and molecular breeding of Lagerstroemia species.

Comparative Chloroplast Genomics of Fritillaria (Liliaceae), Inferences for Phylogenetic Relationships between Fritillaria and Lilium and Plastome Evolution

Fritillaria is a genus that has important medicinal and horticultural values. The study involved the most comprehensive chloroplast genome samples referring to Old and New World clades of Fritillaria for marker selection and phylogenetic studies. We reported and compared eleven newly sequenced whole-plastome sequences of Fritillaria which proved highly similar in overall size (151,652–152,434 bp), genome structure, gene content, and order. Comparing them with other species of Liliales (6 out of 10 families) indicated the same similarity but showed some structural variations due to the contraction or expansion of the inverted repeat (IR) regions. A/T mononucleotides, palindromic, and forward repeats were the most common types. Six hypervariable regions (rps16-trnQ, rbcL-accD, accD-psaI, psaJ-rpl33, petD-rpoA, and rpl32-trnL) were discovered based on 26 Fritillaria whole-plastomes to be potential molecular markers. Based on the plastome data that were collected from 26 Fritillaria and 21 Lilium species, a phylogenomic study was carried out with three Cardiocrinum species as outgroups. Fritillaria was sister to Lilium with a high support value, and the interspecies relationships within subgenus Fritillaria were resolved very well. The six hypervariable regions can be used as candidate DNA barcodes of Fritillaria and the phylogenomic framework can guide extensive genomic sampling for further phylogenetic analyses.

Conservation and innovation: Plastome evolution during rapid radiation of Rhodiola on the Qinghai-Tibetan Plateau

The amount of plastome sequence data available has soared in the last decade, but the nature of plastome evolution during rapid radiations is largely unknown. Moreover, although there is increasing evidence showing that plastomes may have undergone adaptive evolution in order to allow adaptation to various environments, few studies have systematically investigated the role of the plastome in alpine adaptation. To address these questions, we sequenced and analyzed 12 representative species of Rhodiola, a genus which includes ca. 70 perennial herbs mainly growing in alpine habitats in the Qinghai-Tibet Plateau and the Hengduan Mountains. Rapid radiation in this genus was triggered by the uplift of the Qinghai-Tibet Plateau. We also included nine species of Crassulaceae as the outgroups. All plastomes were conserved with respect to size, structure, and gene content and order, with few variations: each contained 134 genes, including 85 protein-coding genes, 37 tRNAs, 8 rRNAs, and 4 potential pseudogenes. Four types of repeat sequence were detected. Slight contraction and expansion of the inverted repeats were also revealed. Both the genome-wide alignment and sequence polymorphism analyses showed that the inverted repeats and coding regions were more conserved than the single-copy regions and the non-coding regions. Positive selection analyses identified three genes containing sites of positive selection (rpl16, ndhA, ndhH), and one gene with a faster than average rate of evolution (psaA). The products of these genes may be involved in the adaptation of Rhodiola to alpine environments such as low CO₂ concentration and high-intensity light.

The complete chloroplast genome sequence of Cornus sunhangii (Cornaceae)

Cornus sunhangii is a newly discovered species of Cornaceae and is endemic to Tibet, China. In this study, the complete chloroplast genome of C. sunhangii is described. The length of the chloroplast genome is 157,446 bp. The genome has quadripartite structure, which includes a large single-copy (LSC) region (86,861 bp), a small single-copy (SSC) region (18,331 bp), and a pair of inverted repeat (IR) regions (26,127 bp each). The complete chloroplast genome contains a total of 113 unique genes, of which 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes.

Chloroplast Genomes and Comparative Analyses among Thirteen Taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae)

The Myrsinaceae s.str. clade is a tropical woody representative in Myrsinoideae of Primulaceae and has ca. 1300 species. The generic limits and alignments of this clade are unclear due to the limited number of genetic markers and/or taxon samplings in previous studies. Here, the chloroplast (cp) genomes of 13 taxa within the Myrsinaceae s.str. clade are sequenced and characterized. These cp genomes are typical quadripartite circle molecules and are highly conserved in size and gene content. Three pseudogenes are identified, of which ycf15 is totally absent from five taxa. Noncoding and large single copy region (LSC) exhibit higher levels of nucleotide diversity (Pi) than other regions. A total of ten hotspot fragments and 796 chloroplast simple sequence repeats (SSR) loci are found across all cp genomes. The results of phylogenetic analysis support the notion that the monophyletic Myrsinaceae s.str. clade has two subclades. Non-synonymous substitution rates (d_N) are higher in housekeeping (HK) genes than photosynthetic (PS) genes, but both groups have a nearly identical synonymous substitution rate (d_S). The results indicate that the PS genes are under stronger functional constraints compared with the HK genes. Overall, the study provides hypervariable molecular markers for phylogenetic reconstruction and contributes to a better understanding of plastid gene evolution in Myrsinaceae s.str. clade.

Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastid genomes

BACKGROUND AND AIMS

Phylogenetic relationships within tribe Shoreeae, containing the main elements of tropical forests in Southeast Asia, present a long-standing problem in the systematics of Dipterocarpaceae. Sequencing whole plastomes using next-generation sequencing- (NGS) based genome skimming is increasingly employed for investigating phylogenetic relationships of plants. Here, the usefulness of complete plastid genome sequences in resolving phylogenetic relationships within Shoreeae is evaluated.

METHODS

A pipeline to obtain alignments of whole plastid genome sequences across individuals with different amounts of available data is presented. In total, 48 individuals, representing 37 species and four genera of the ecologically and economically important tribe Shoreeae sensu Ashton, were investigated. Phylogenetic trees were reconstructed using maximum parsimony, maximum likelihood and Bayesian inference.

KEY RESULTS

Here, the first fully sequenced plastid genomes for the tribe Shoreeae are presented. Their size, GC content and gene order are comparable with those of other members of Malvales. Phylogenomic analyses demonstrate that whole plastid genomes are useful for inferring phylogenetic relationships among genera and groups of Shorea (Shoreeae) but fail to provide well-supported phylogenetic relationships among some of the most closely related species. Discordance in placement of Parashorea was observed between phylogenetic trees obtained from plastome analyses and those obtained from nuclear single nucleotide polymorphism (SNP) data sets identified in restriction-site associated sequencing (RADseq).

CONCLUSIONS

Phylogenomic analyses of the entire plastid genomes are useful for inferring phylogenetic relationships at lower taxonomic levels, but are not sufficient for detailed phylogenetic reconstructions of closely related species groups in Shoreeae. Discordance in placement of Parashorea was further investigated for evidence of ancient hybridization.

A systematic comparison of eight new plastome sequences from Ipomoea L

Background

Ipomoea is the largest genus in the family Convolvulaceae. The species in this genus have been widely used in many fields, such as agriculture, nutrition, and medicine. With the development of next-generation sequencing, more than 50 chloroplast genomes of Ipomoea species have been sequenced. However, the repeats and divergence regions in Ipomoea have not been well investigated. In the present study, we sequenced and assembled eight chloroplast genomes from sweet potato's close wild relatives. By combining these with 32 published chloroplast genomes, we conducted a detailed comparative analysis of a broad range of Ipomoea species.

Methods

Eight chloroplast genomes were assembled using short DNA sequences generated by next-generation sequencing technology. By combining these chloroplast genomes with 32 other published Ipomoea chloroplast genomes downloaded from GenBank and the Oxford Research Archive, we conducted a comparative analysis of the repeat sequences and divergence regions across the Ipomoea genus. In addition, separate analyses of the Batatas group and Quamoclit group were also performed.

Results

The eight newly sequenced chloroplast genomes ranged from 161,225 to 161,721 bp in length and displayed the typical circular quadripartite structure, consisting of a pair of inverted repeat (IR) regions (30,798-30,910 bp each) separated by a large single copy (LSC) region (87,575-88,004 bp) and a small single copy (SSC) region (12,018-12,051 bp). The average guanine-cytosine (GC) content was approximately 40.5% in the IR region, 36.1% in the LSC region, 32.2% in the SSC regions, and 37.5% in complete sequence for all the generated plastomes. The eight chloroplast genome sequences from this study included 80 protein-coding genes, four rRNAs (rrn23, rrn16, rrn5, and rrn4.5), and 37 tRNAs. The boundaries of single copy regions and IR regions were highly conserved in the eight chloroplast genomes. In Ipomoea, 57-89 pairs of repetitive sequences and 39-64 simple sequence repeats were found. By conducting a sliding window analysis, we found six relatively high variable regions (ndhA intron, ndhH-ndhF, ndhF-rpl32, rpl32-trnL, rps16-trnQ, and ndhF) in the Ipomoea genus, eight (trnG, rpl32-trnL, ndhA intron, ndhF-rpl32, ndhH-ndhF, ccsA-ndhD, trnG-trnR, and pasA-ycf3) in the Batatas group, and eight (ndhA intron, petN-psbM, rpl32-trnL, trnG-trnR, trnK-rps16, ndhC-trnV, rps16-trnQ, and trnG) in the Quamoclit group. Our maximum-likelihood tree based on whole chloroplast genomes confirmed the phylogenetic topology reported in previous studies.

Conclusions

The chloroplast genome sequence and structure were highly conserved in the eight newly-sequenced Ipomoea species. Our comparative analysis included a broad range of Ipomoea chloroplast genomes, providing valuable information for Ipomoea species identification and enhancing the understanding of Ipomoea genetic resources.

Integrative taxonomy of commercially important deep water penaeoid shrimps from India

The deep water penaeoid shrimp is an important commercial crustacean resource along the Indian coast. The molecular and morphological information of this group from the Indian coast is scarcely known. In this study, we investigated the identification and phylogenetic relationships of the deep water penaeoid shrimps using three mitochondrial (cytochrome oxidase subunit I (COI), cytochrome b, 16S rRNA) genes, which were compared with 54 morphological characters and further used to evaluate character evolution. Our study revealed remarkable molecular divergence (3.3-33.0%) in nine species from three genera of Solenoceridae, four species from three genera of Penaeidae and one species from Aristeidae using COI. Phylogenetic analysis using maximum likelihood and Bayesian approaches revealed that all species from these families are monophyletic. The present analysis revealed the existence of subgroups in the genus Solenocera suggesting the slow reduction of postrostral carina which corresponds to the increase in distributional depth during the evolutionary process which further indicates the origin of the genus in the continental shelf and extending up to the continental slope. In addition, we generated the DNA barcode database involving these species which can help further to investigate the detailed evolution and biogeography of these valuable crustacean resources.

Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide

Taxus (yew) is both the most species-rich and taxonomically difficult genus in Taxaceae. To date, no study has elucidated the complexities of the plastid genome (plastome) or examined the possibility of whole plastomes as super-barcodes across yew species worldwide. In this study, we sequenced plastomes from two to three individuals for each of the 16 recognized yew species (including three potential cryptics) and Pseudotaxus chienii. Our comparative analyses uncovered several gene loss events that independently occurred in yews, resulting in a lower plastid gene number than other Taxaceous genera. In Pseudotaxus and Taxus, we found two isomeric arrangements that differ by the orientation of a 35 kb fragment flanked by "trnQ-IRs". These two arrangements exist in different ratios within each sampled individual, and intraspecific shifts in major isomeric arrangements are first reported here in Taxus. Moreover, we demonstrate that entire plastomes can be used to successfully discriminate all Taxus species with 100% support, suggesting that they are useful as super-barcodes for species identification. We also propose that accD and rrn16-rrn23 are promising special barcodes to discriminate yew species. Our newly developed Taxus plastomic sequences provide a resource for super-barcodes and conservation genetics of several endangered yews and serve as comprehensive data to improve models of plastome complexity in Taxaceae as a whole and authenticate Taxus species.

First reported chloroplast genome sequence of Punica granatum (cultivar Helow) from Jabal Al-Akhdar, Oman: phylogenetic comparative assortment with Lagerstroemia

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It has grown in popularity and is a profitable fruit crop due to its attractive features including a bright red appearance and its biological activities. Scientific exploration of the genetics and evolution of these beneficial traits has been hampered by limited genomic information. In this study, we sequenced the complete chloroplast (cp) genome of the native P. granatum (cultivar Helow) cultivated in the mountains of Jabal Al-Akhdar, Oman. The results revealed a P. granatum cp genome length of 158,630 bp, characterized by a relatively conserved structure containing 2 inverted repeat regions of 25,466 bp, an 18,686 bp small single copy regions, and an 89,015 bp large single copy region. The 86 protein-coding genes included 37 transfer RNA genes and 8 ribosomal RNA genes. Comparison of the P. granatum whole cp genome with seven Lagerstroemia species revealed an overall high degree of sequence similarity with divergence among intergenic spacers. The location, distribution, and divergence of repeat sequences and shared genes of the Punica and Lagerstroemia species were highly similar. Analyses of nucleotide substitution, insertion/deletions, and highly variable regions in these cp genomes identified potential plastid markers for taxonomic and phylogenetic studies in Myrtales. A phylogenetic study of the cp genomes and 76 shared coding regions generated similar cladograms. The complete cp genome of P. granatum will aid in taxonomical studies of the family Lythraceae.

Phylogenomic and Comparative Analyses of Complete Plastomes of Croomia and Stemona (Stemonaceae)

The monocot genus Croomia (Stemonaceae) comprises three herbaceous perennial species that exhibit EA (Eastern Asian)–ENA (Eastern North American) disjunct distribution. However, due to the lack of effective genomic resources, its evolutionary history is still weakly resolved. In the present study, we conducted comparative analysis of the complete chloroplast (cp) genomes of three Croomia species and two Stemona species. These five cp genomes proved highly similar in overall size (154,407–155,261 bp), structure, gene order and content. All five cp genomes contained the same 114 unique genes consisting of 80 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Gene content, gene order, AT content and IR/SC boundary structures were almost the same among the five Stemonaceae cp genomes, except that the Stemona cp genome was found to contain an inversion in cemA and petA. The lengths of five genomes varied due to contraction/expansion of the IR/SC borders. A/T mononucleotides were the richest Simple Sequence Repeats (SSRs). A total of 46, 48, 47, 61 and 60 repeats were identified in C. japonica, C. heterosepala, C. pauciflora, S. japonica and S. mairei, respectively. A comparison of pairwise sequence divergence values across all introns and intergenic spacers revealed that the ndhF–rpl32, psbM–trnD and trnS–trnG regions are the fastest-evolving regions. These regions are therefore likely to be the best choices for molecular evolutionary and systematic studies at low taxonomic levels in Stemonaceae. Phylogenetic analyses of the complete cp genomes and 78 protein-coding genes strongly supported the monophyly of Croomia. Two Asian species were identified as sisters that likely diverged in the Early Pleistocene (1.62 Mya, 95% HPD: 1.125–2.251 Mya), whereas the divergence of C. pauciflora dated back to the Late Miocene (4.77 Mya, 95% HPD: 3.626–6.162 Mya). The availability of these cp genomes will provide valuable genetic resources for further population genetics and phylogeographic studies on Croomia.

Plant phylogenomics based on genome-partitioning strategies: Progress and prospects

The rapid expansion of next-generation sequencing (NGS) has generated a powerful array of approaches to address fundamental questions in biology. Several genome-partitioning strategies to sequence selected subsets of the genome have emerged in the fields of phylogenomics and evolutionary genomics. In this review, we summarize the applications, advantages and limitations of four NGS-based genome-partitioning approaches in plant phylogenomics: genome skimming, transcriptome sequencing (RNA-seq), restriction site associated DNA sequencing (RAD-Seq), and targeted capture (Hyb-seq). Of these four genome-partitioning approaches, targeted capture (especially Hyb-seq) shows the greatest promise for plant phylogenetics over the next few years. This review will aid researchers in their selection of appropriate genome-partitioning approaches to address questions of evolutionary scale, where we anticipate continued development and expansion of whole-genome sequencing strategies in the fields of plant phylogenomics and evolutionary biology research.

Floral morphology and morphogenesis in Camptotheca (Nyssaceae), and its systematic significance

BACKGROUND AND AIMS

Camptotheca is endemic to China and there are limited data about the breeding system and morphogenesis of the flowers. Camptotheca is thought to be related to Nyssa and Davidia in Nyssaceae, which has sometimes been included in Cornaceae. However, molecular phylogenetic studies confirmed the inclusion of Camptotheca in Nyssaceae and its exclusion from Cornaceae. The aim of this study was to reveal developmental features of the inflorescence and flowers in Camptotheca to compare with related taxa in Cornales.

METHODS

Inflorescences and flowers of Camptotheca acuminata at all developmental stages were collected and studied with a scanning electron microscope and stereo microscope.

KEY RESULTS

Camptotheca has botryoids which are composed of several capitate floral units (FUs) that are initiated acropetally. On each FU, flowers are grouped in dyads that are initiated acropetally. All floral organs are initiated centripetally. Calyx lobes are restricted to five teeth. The hypanthium, with five toothed calyx lobes, is adnate to the ovary. The five petals are free and valvate. Ten stamens are inserted in two whorls around the central depression, in which the style is immersed. Three carpels are initiated independently but the ovary is syncarpous and unilocular. The ovule is unitegmic and heterotropous. Inflorescences are functionally andromonoecious varying with the position of the FUs on the inflorescence system. Flowers on the upper FU often have robust styles and fully developed ovules. Flowers on the lower FU have undeveloped styles and aborted ovules, and the flowers on the middle FU are transitional.

CONCLUSIONS

Camptotheca possesses several traits that unify it with Nyssa, Mastixia and Diplopanax. Inflorescence and floral characters support a close relationship with Nyssaceae and Mastixiaceae but a distant relationship with Cornus. Our results corroborate molecular inferences and support a separate family Nyssaceae.