Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae).

BACKGROUND: The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed. RESULTS: Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness. CONCLUSIONS: Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae).

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants.

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae).

BACKGROUND AND AIMS: Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse. METHODS: We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains. KEY RESULTS: The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification. CONCLUSIONS: The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

Plastid phylogenomics and biogeographic analysis support a trans-Tethyan origin and rapid early radiation of Cornales in the Mid-Cretaceous.

The Cornales is a relatively small but morphologically diverse order in the basal position of the Asterids clade. Previous study hypothesized that the order might have undergone ancient rapid radiation during the Cretaceous when major angiosperm lineages were established. We conducted the phylogenomic analysis of Cornales using 81 plastid genome sequences with 67 newly generated in this study to test the hypothesis. This sampling represents all the families and 31 out of 48 genera in the order. Phylogenetic analyses were conducted using different datasets to examine the effects of different coding positions and character coding methods. We further conducted divergence time, diversification rate, and biogeographic analyses to understand the early evolutionary history of Cornales in space and time. Our phylogenetic analyses of four datasets (the amino acid characters, the 1st and 2nd codon positions of protein coding genes, nucleotide characters with degenerated coding method, and noncoding regions) resulted in a robust phylogeny congruent with results of previous studies, showing (((Cornaceae-Alangiaceae)-(Curtisiaceae-Grubbiaceae))-(((Nyssaceae-Davidiaceae)-Mastixiaceae)-((Hydrostachyaceae-(Hydrangeaceae-Loasaceae)))). Phylogenetic relationships within families were also well resolved. Conflicts in the placement of Hydrostachyaceae were found from analyses of two datasets, the nucleotide characters of all codon position and the 3rd codon positions, where the family was united with Loasaceae, but not strongly supported. Results from divergence time analyses suggested a mid-Cretaceous origin of Cornales followed by rapid early diversification into major clades/families within 10 million years. The early diversification of Cornales may have been facilitated by divergence in habitat and morphology following geographic dispersals. The ancestral distribution of the order was inferred as a widespread range covering Asia, Europe, North America, and Africa when including fossils in the analyses, suggesting an origin of the order likely along the Tethys Seaway where the areas were connected in the mid-Cretaceous. Inferred geographic origins of each family differed to some extent between analyses including fossils vs excluding fossils. In the analysis with extant and fossil species, the origins of the African Hydrostachyaceae and Grubbiaceae-Curtisiaceae clade were inferred to have involved two independent events, an intercontinental dispersal from the northern hemisphere to Africa and an intercontinental vicariance between the northern hemisphere and Africa, respectively. Other families were inferred to have evolved in the northern hemisphere with subsequent intercontinental dispersal(s) to other areas including to Central and South America, during their subsequent diversification. Net diversification rate analysis based on treePL dated phylogeny using MEDUSA detected a nearly 5-fold decrease in the African endemic Curtisiaceae-Grubbiaceae (CuG) clade and an increase of rate in the Hydrangeaceae-Loasaceae (HL) clade. Within HL, a decrease in the Fendlera-Jamesia clade and an increase in the Philadelphus clade were also detected. The findings are also consistent with the level of present species diversity in these lineages. Our study demonstrated the value of plastid genome in phylogenomic study, but posed an old challenge of biogeographic study with fossil data and raised caution for the synonymous substitution sites of plastid genome in phylogenomics studies.

Niche convergence and biogeographic history shape elevational tree community assembly in a subtropical mountain forest.

Niche convergence or conservatism have been proposed as essential mechanisms underlying elevational plant community assembly in tropical mountain ecosystems. Subtropical mountains, compared to tropical mountains, are likely to be shaped by a mixing of different geographic affinities of species and remain somehow unclear. Here, we used 31 0.1-ha permanent plots distributed in subtropical forests on the eastern and western aspects of the Gaoligong Mountains, southwest China between 1498 m and 3204 m a.sl. to evaluate how niche-based and biogeographic processes shape tree community assembly along elevational gradients. We analyzed the elevational patterns of taxonomic, phylogenetic and functional diversity, as well as of individual traits, and assessed the relative importance of environmental effects on these diversity measures. We then classified tree species as being either tropical affiliated or temperate affiliated and estimated their contribution to the composition of biogeographic affinities. Species richness decreased with elevation, and species composition showed apparent turnover across the aspects and elevations. Most traits exhibited convergent patterns across the entire elevational gradient. Phylogenetic and functional diversity showed opposing patterns, with phylogenetic diversity increasing and functional diversity decreasing with elevation. Soil nutrients, especially phosphorus and nitrogen, appeared to be the main abiotic variables driving the elevational diversity patterns. Communities at lower elevations were occupied by tropical genera, while highlands contained species of tropical and temperate biogeographic affinities. Moreover, the high phylogenetic diversity at high elevations were likely due to differences in evolutionary history between temperate and tropical species. Our results highlight the importance of niche convergence of tropical species and the legacy of biogeographic history on the composition and structure of subtropical mountain forests. Furthermore, limited soil phosphorus caused traits divergence and the partitioning for different forms of phosphorus may explain the high biodiversity found in phosphorus-limited subtropical forests.

Amentotaxus×hybridia (Taxaceae), a new natural Amentotaxus hybrid from southeast Yunnan province, China.

During floristic surveys of Taxaceae in Hekou County, Yunnan Province, China, a putative natural hybrid between A.yunnanensis H.L. Li and A.hekouensis L.M. Gao was collected. Morphological and molecular evidence confirms its status as a natural hybrid. Amentotaxus×hybridia L.M. Gao has linear or linear-lanceolate leaves 6-13 cm × 1.0-1.5 cm, white stomatal bands with 34-40 rows on abaxial side, 2.5-3.5 mm, slightly wider than leaf margins; 3-6 seeds borne at the base of the branchlet, peduncle 1.3-1.6 cm long with 3-4 rows of persistent basal bracts.

Testing the Complete Plastome for Species Discrimination, Cryptic Species Discovery and Phylogenetic Resolution in Cephalotaxus (Cephalotaxaceae).

Species of Cephalotaxus have great economic and ecological values. However, the taxonomy and interspecific phylogenetic relationships within the genus have been controversial and remained not fully resolved until now. To date, no study examined the efficiency of the complete plastome as super-barcode across Cephalotaxus species with multiple samples per taxon. In this study, we have evaluated the complete plastome in species discrimination and phylogenetic resolution in Cephalotaxus by including 32 individuals of all eight recognized species and five varieties following Farjon's classification (2010) with multiple samples per taxon. Our results indicated that not all species recognized in recent taxonomic revisions of Cephalotaxus could be distinguished and not all were monophyletic. Based on the plastome phylogeny, a new taxonomic classification for the genus comprising nine species and two varieties, including a cryptic species, was proposed. The phylogeny also resolved all interspecific relationships. Compared to the plastome based classification, standard DNA barcodes, alone or in combination, only recognized a maximum of seven out of the nine species. Moreover, two highly variable single loci, ycf1 and rps16, each alone achieved full species discrimination. With the moderate length of 1079 bp, rps16 is proposed as a specific barcode to discriminate Cephalotaxus species. The super-barcodes and specific barcode candidates will aid in the identification of endangered Cephalotaxus species, and to help focus conservation measures.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron.

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

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.