Phylogenomic analyses reveal a deep history of hybridization and polyploidy in the Neotropical genus Lachemilla (Rosaceae)

Phylogenomics provides new insight into the phylogeny and diversification of Asian Lappula (Boraginaceae)

The application of omics data serves as a powerful tool for investigating the roles of incomplete lineage sorting (ILS) and hybridization in shaping genomic diversity, offering deeper insights into complex evolutionary processes. In this study, we utilized deep genome sequencing data from 76 individuals of Lappula and its closely allied genera, collected from China and Central Asia. By employing the HybPiper and Easy353 pipelines, we recovered 262-279 single-copy nuclear genes (SCNs) and 352-353 Angiosperms353 genes, respectively. We analyzed multiple datasets, including complete chloroplast genomes and a filtered set of 475 SCNs, to conduct phylogenetic analyses using both concatenated and coalescent-based methods. Furthermore, we employed Quartet Sampling (QS), coalescent simulations, MSCquartets, HyDe, and reticulate network analyses to investigate the sources of phylogenetic discordance. Our results confirm that Lappula is polyphyletic, with L. mogoltavica clustering with Pseudolappula sinaica and forming a sister relationship with other taxa included in this study. Additionally, three Lepechiniella taxa nested within distinct clades of Lappula. Significant gene tree discordance was observed at several nodes within Lappula. Coalescent simulations and hybrid detection analyses suggest that both ILS and hybridization contribute to these discrepancies. Flow cytometry (FCM) analyses confirmed the presence of both diploid and tetraploid taxa within Lappula. Phylogenetic network analyses further revealed that Clades IV and VII likely originated through hybridization, with the tetraploids in Clade IV arising from two independent hybridization events. Additionally, the "ghost lineage" identified as sister to Lappula redowskii serves as one of the donors in allopolyploidization. In conclusion, our study provides new insights into the deep phylogenetic relationships of Asian Lappula and its closely allied genera, contributing to a more comprehensive understanding of the evolution and diversification of Lappula.

An integrative framework reveals widespread gene flow during the early radiation of oaks and relatives in Quercoideae (Fagaceae)

Although the frequency of ancient hybridization across the Tree of Life is greater than previously thought, little work has been devoted to uncovering the extent, timeline, and geographic and ecological context of ancient hybridization. Using an expansive new dataset of nuclear and chloroplast DNA sequences, we conducted a multifaceted phylogenomic investigation to identify ancient reticulation in the early evolution of oaks (Quercus). We document extensive nuclear gene tree and cytonuclear discordance among major lineages of Quercus and relatives in Quercoideae. Our analyses recovered clear signatures of gene flow against a backdrop of rampant incomplete lineage sorting, with gene flow most prevalent among major lineages of Quercus and relatives in Quercoideae during their initial radiation, dated to the Early-Middle Eocene. Ancestral reconstructions including fossils suggest ancestors of Castanea + Castanopsis, Lithocarpus, and the Old World oak clade probably co-occurred in North America and Eurasia, while the ancestors of Chrysolepis, Notholithocarpus, and the New World oak clade co-occurred in North America, offering ample opportunity for hybridization in each region. Our study shows that hybridization-perhaps in the form of ancient syngameons like those seen today-has been a common and important process throughout the evolutionary history of oaks and their relatives. Concomitantly, this study provides a methodological framework for detecting ancient hybridization in other groups.

Extensive cytonuclear discordance revealed by phylogenomic analyses suggests complex evolutionary history in the holly genus Ilex (Aquifoliaceae)

Ilex L., the exclusive genus of Aquifoliaceae, encompasses over 600 dioecious wood species with a highly irregular distribution, predominantly found in South America and Asia. The phylogeny and classification of this genus remain enigmatic due to significant early extinctions, constrained morphological diversity, recent hybridization/introgression, and conflicting signals from previously utilized markers. This study presents phylogenetic reconstructions based on complete chloroplast genome sequences and single nucleotide polymorphisms (SNPs) derived from genome resequencing data. A total of 116 accessions of Ilex, representing approximately 108 taxa, were included as the ingroup, with five accessions of two species serving as outgroups. Analysis of the chloroplast genome and nuclear SNP data individually resulted in two robust phylogenetic trees, revealing substantial discrepancies between the chloroplast genome and nuclear SNP phylogenies at both the species and clade levels. The chloroplast genome sequences successfully resolved relationships within this genus into eight strongly supported major clades, while the nuclear SNPs resolved relationships into seven highly supported major clades. Our nuclear SNP phylogenetic tree, in comparison to the chloroplast genome tree, aligns more closely with the recently updated classification of Ilex in multiple instances. The extensive cytonuclear discordance identified may be attributed to recent hybridization events and incomplete lineage sorting (ILS).

Polyploidy linked with species richness but not diversification rates or niche breadth in Australian Pomaderreae (Rhamnaceae)

BACKGROUND AND AIMS

Polyploidy is an important evolutionary driver for plants and has been linked with higher species richness and increases in diversification rate. These correlations between ploidy and plant radiations could be the result of polyploid lineages exploiting broader niche space and novel niches due to their enhanced adaptability. The evolution of ploidy and its link to plant diversification across the Australian continent is not well understood. Here, we focus on the ploidy evolution of the Australasian Rhamnaceae tribe Pomaderreae.

METHODS

We generated a densely sampled phylogeny (90 %, 215/240 species) of the tribe and used it to test for the evolution of ploidy. We obtained 30 orthologous nuclear loci per sample and dated the phylogeny using treePL. Ploidy estimates for each sequenced species were obtained using nQuire, based on phased sequence data. We used MiSSE to obtain tip diversification rates and tested for significant relationships between diversification rates and ploidy. We also assessed for relationships between ploidy level and niche breadth, using distributional records, species distributional modelling and WorldClim data.

KEY RESULTS

Polyploidy is extensive across the tribe, with almost half (45 %) of species and the majority of genera exhibiting this trait. We found a significant positive relationship between polyploidy and genus size (i.e. species richness), but a non-significant positive relationship between polyploidy and diversification rates. Polyploidy did not result in significantly wider niche space occupancy for Pomaderreae; however, polyploidy did allow transitions into novel wetter niches. Spatially, eastern Australia is the diversification hotspot for Pomaderreae in contrast to the species hotspot of south-west Western Australia.

CONCLUSIONS

The relationship between polyploidy and diversification is complex. Ancient polyploidization events likely played an important role in the diversification of species-rich genera. A lag time effect may explain the uncoupling of tip diversification rates and polyploidy of extant lineages. Further studies on other groups are required to validate these hypotheses.

Phylotranscriptomic analyses reveal the evolutionary complexity of Paris L. (Melanthiaceae), a morphologically distinctive genus with significant pharmaceutical importance

BACKGROUND AND AIMS

Previous phylogenetic studies on the pharmaceutically significant genus Paris (Melanthiaceae) have consistently revealed substantial cytonuclear discordance, yet the underlying mechanism responsible for this phenomenon remains elusive. This study aims to reconstruct a robust nuclear backbone phylogeny and elucidate the potential evolutionarily complex events contributing to previously observed cytonuclear discordance within Paris.

METHODS

Based on a comprehensive set of nuclear low-copy orthologous genes obtained from transcriptomic data, the intrageneric phylogeny of Paris, along with its phylogenetic relationships to allied genera, were inferred using coalescent and concatenated approaches. The analysis of gene tree discordance and reticulate evolution, in conjunction with an incomplete lineage sorting (ILS) simulation, was conducted to explore potential hybridization and ILS events in the evolutionary history of Paris and assess their contribution to the discordance of gene trees.

KEY RESULTS

The nuclear phylogeny unequivocally confirmed the monophyly of Paris and its sister relationship with Trillium, while widespread incongruences in gene trees were observed at the majority of internal nodes within Paris. The reticulate evolution analysis identified five instances of hybridization events in Paris, indicating that hybridization events might have occurred recurrently throughout the evolutionary history of Paris. In contrast, the ILS simulations revealed that only two internal nodes within section Euthyra experienced ILS events.

CONCLUSIONS

Our data suggest that the previously observed cytonuclear discordance in the phylogeny of Paris can primarily be attributed to recurrent hybridization events, with secondary contributions from infrequent ILS events. The recurrent hybridization events in the evolutionary history of Paris not only drove lineage diversification and speciation but also facilitated morphological innovation, and enhanced ecological adaptability. Therefore, artificial hybridization has great potential for breeding medicinal Paris species. These findings significantly contribute to our comprehensive understanding of the evolutionary complexity of this pharmaceutically significant plant lineage, thereby facilitating effective exploitation and conservation efforts.

Comparative Analysis of the Chloroplast Genomes of the Melliodendron (Styracaceae) Species: Providing Insights into Molecular Evolution and Phylogenetic Relationships

Melliodendron xylocarpum is a member of the Styracaceae family, which is well-known for its remarkable ornamental and medicinal properties. In this research, we conducted comparative analysis of the chloroplast genomes from four samples of M. xylocarpum, representing Melliodendron. The results demonstrated that the chloroplast genome of four M. xylocarpum samples ranging from 157,103 bp to 158,357 bp exhibited a typical quadripartite structure, including one large single-copy (LSC) region (90,131 bp to 90,342 bp), one small single-copy (SSC) region (18,467 bp to 18,785 bp), and two inverted repeat regions (IRs) (24,115 bp to 24,261 bp). Different levels of expansion and contraction were observed in the IR region of four M. xylocarpum samples. Besides, accD and ycf1 have been identified under positive selection, potentially linked to the adaptive response of Melliodendron to various environmental changes. Conflicting phylogenetic relationships were identified among various genera within the Styracaceae family in the phylogenetic tree constructed using CDS sequences and complete chloroplast genomes. Furthermore, the significance of a large sample size was also highlighted in this study for enhancing the accuracy of findings from phylogenetic analyses. The findings of this research will provide significant insights for future investigations into the evolutionary trends and conservation of the Melliodendron species.

Prezygotic barriers effectively limit hybridization in a rapid evolutionary radiation

Hybridization is increasingly recognized as an important evolutionary process across the tree of life. In many clades, phylogenomic approaches have permitted unparalleled insight into the extent and frequency of hybridization. However, we continue to lack a deep understanding of the factors that limit and shape patterns of hybridization, especially in evolutionary radiations. In this study, we characterized patterns of introgression across Costus (Costaceae), a young evolutionary radiation of tropical understory plants that maintain widespread interfertility despite exhibiting strong prezygotic reproductive isolation. We analyzed a phylogenomic dataset of 756 genes from 54 Costus species using multiple complementary approaches - D-statistics, gene-tree-based tests, and phylogenetic network analyses - to detect and characterize introgression events throughout the evolutionary history of the radiation. Our results identified a moderate number of introgression events, including a particularly ancient, well-supported event spanning one of the deepest divergences in the clade. Most introgression events occurred between taxa or ancestral lineages that shared the same pollination syndrome (bee-pollinated or hummingbird-pollinated). These findings suggest that prezygotic barriers, including pollinator specialization, have been key to the balance between introgression and reproductive isolation in Costus.

Phylogenomics and phylogeographic model testing using convolutional neural networks reveal a history of recent admixture in the Canarian Kleinia neriifolia

Multiple-island endemics (MIE) are considered ideal natural subjects to study patterns of island colonization that involve recent population-level genetic processes. Kleinia neriifolia is a Canarian MIE widespread across the archipelago, which exhibits a close phylogenetic relationship with species in northwest Africa and at the other side of the Sahara Desert. Here, we used target sequencing with plastid skimming (Hyb-Seq), a dense population-level sampling of K. neriifolia, and representatives of its African-southern Arabian relatives to infer phylogenetic relationships and divergence times at the species and population levels. Using population genetic techniques and machine learning (convolutional neural networks [CNNs]), we reconstructed phylogeographic relationships and patterns of genetic admixture based on a multilocus SNP nuclear dataset. Phylogenomic analysis based on the nuclear dataset identifies the northwestern African Kleinia anteuphorbium as the sister species of K. neriifolia, with divergence starting in the early Pliocene. Divergence from its sister clade, comprising species from the Horn of Africa and southern Arabia, is dated to the arid Messinian period, lending support to the climatic vicariance origin of the Rand Flora. Phylogeographic model testing with CNNs supports an initial colonization of the central island of Tenerife followed by eastward and westward migration across the archipelago, which resulted in the observed east/west phylogeographic split. Subsequent population extinctions linked to aridification events, and recolonization from Tenerife, are proposed to explain the patterns of genetic admixture in the eastern Canary Islands. We demonstrate that CNNs based on SNPs can be used to discriminate among complex scenarios of island migration and colonization.

Apocynaceae wood evolution matches key morphological innovations

PREMISE

This paper provides an overview of the wood anatomy of the dogbane family (Apocynaceae), reconstructs wood anatomical trait evolution, and links this evolution with woody growth-form transitions and floral and seed trait innovations across the family.

METHODS

Over 200 published wood anatomical descriptions were revised, and original light microscopic sections were made and described for another 50 species. Changes in wood anatomical characters through time were visualized with ancestral state reconstructions. Tests for correlated evolution were performed using a combined data set of anatomical and key floral and seed traits to identify potential synnovations and traits associated with growth-form adaptations.

RESULTS

There was a shift toward a suite of wood anatomical traits that separate the rauvolfioids and early-branching apocynoids from the core apocynoids, including an increased presence of vessel multiples, vessel dimorphism, laticifers, vascular (cambial) variants, and paratracheal axial parenchyma. The presence of this trait suite, which continues in Periplocoideae, Secamonoideae, and Asclepiadoideae, coincides with a progression of floral morphological innovations that evolved on consecutive nodes in the family, and also relates to more frequent transitions toward the climbing and herbaceous habits. In addition, a considerable shortening of vessel elements and fibers along the phylogenetic backbone of the family is correlated with a general reduction in plant size.

CONCLUSIONS

There are clear evolutionary transitions in the wood anatomy of Apocynaceae representing structural adaptations across the family that are associated with a quick succession of evolutionary changes of the floral bauplan.

Phylogenomic insights into Neotropical Magnolia relationships

Despite extensive research into the phylogenetic relationships of the genus Magnolia, Neotropical taxa have been neglected. This is partly because their numbers have recently doubled and now account for almost half of the global richness. Therefore, by sampling one-third of all Neotropical taxa their relationships were studied using morphological, nuclear, and plastome data. Two major clades were identified: Clade I, comprising Magnolia sect. Talauma, Magnolia sect. Splendentes from the Neotropics, and the Asian Magnolia sect. Gwillimia; while Clade II included the Neotropical clades Magnolia sect. Macrophylla and Magnolia sect. Magnolia, along with the remaining non-Neotropical sections. Within Clade I, Magnolia sect. Talauma was geographically divided into a northern subclade grouping Mexican and Central American taxa, and a southern subclade comprising South American and Caribbean taxa. Magnolia sect. Splendentes was also dichotomously divided, corresponding to the former Magnolia sect. Talauma subsect. Cubenses and Magnolia sect. Talauma subsect. Dugandiodendron. In Clade II, the relationships within Magnolia sect. Macrophylla and Magnolia sect. Magnolia were unclear, suggesting a species complex in all Magnolia sect. Macrophylla taxa. In total, 25 morphological traits were assessed, and ancestral state reconstructions were carried out. Only the joined clustering of mature follicles was a synapomorphy for the southern subclade of Magnolia sect. Talauma. In conclusion, this highlights the need to re-assess the taxonomic delimitation of certain groups, to update the infrageneric classification of Neotropical clades and to explore morphological traits to support them.

Benefits and Limits of Phasing Alleles for Network Inference of Allopolyploid Complexes

Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared with haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where the depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared with using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical nonidentifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes.

The Vitextrifolia complex (Lamiaceae) in the Philippines

The Vitextrifolia complex in the Philippines comprises economically important medicinal species, but its taxonomic status has become ambiguous due to numerous historical nomenclatural revisions. We assembled the complete chloroplast genomes of five species belonging to this species complex to provide additional evidence for their species delimitation. Based on a detailed analysis of specimens which combined molecular and morphological data, we propose reinstating V.elmeri Moldenke and delineating V.arvensis Gentallan, Sengun & M.B.Bartolome as a new endemic species belonging to this complex. The new species is a putative hybrid between V.bicolor Willd. and V.elmeri Moldenke. The specific epithet arvensis reflects its predominantly cultivated nature, both on a commercial scale and in home gardens, as a valued medicinal plant. We also provided a key to identify the five species belonging to the V.trifolia complex in the Philippines.

Recurrent hybridization and gene flow shaped Norway and Siberian spruce evolutionary history over multiple glacial cycles

Most tree species underwent cycles of contraction and expansion during the Quaternary. These cycles led to an ancient and complex genetic structure that has since been affected by extensive gene flow and by strong local adaptation. The extent to which hybridization played a role in this multi-layered genetic structure is important to be investigated. To study the effect of hybridization on the joint population genetic structure of two dominant species of the Eurasian boreal forest, Picea abies and P. obovata, we used targeted resequencing and obtained around 480 K nuclear SNPs and 87 chloroplast SNPs in 542 individuals sampled across most of their distribution ranges. Despite extensive gene flow and a clear pattern of Isolation-by-Distance, distinct genetic clusters emerged, indicating the presence of barriers and corridors to migration. Two cryptic refugia located in the large hybrid zone between the two species played a critical role in shaping their current distributions. The two species repeatedly hybridized during the Pleistocene and the direction of introgression depended on latitude. Our study suggests that hybridization helped both species to overcome main shifts in their distribution ranges during glacial cycles and highlights the importance of considering whole species complex instead of separate entities to retrieve complex demographic histories.

Morphology and ITS sequences provide insights into the phylogeny of Tongoloa (Apiaceae) from China

BACKGROUND

Tongoloa is a genus comprising approximately 20 species, primarily distributed in the mountainous regions of southwest China. The insufficiency of specimen materials and morphological similarities among species render it a taxonomically challenging genus within the Apiaceae family. To elucidate the phylogenetic relationships and taxonomy of Chinese Tongoloa, this study utilized a total of 115 nrITS sequences, including 47 recently obtained sequences, for phylogenetic reconstruction.

RESULTS

Phylogenetic relationships reconstructed from ITS sequences indicate that the East Asia Clade and the Komarovia Clade are sister groups, and Tongoloa belongs to the East Asia Clade. Species of Tongoloa are subdivided into 3 distinct groups, all sharing similar fruit morphologies and are clearly differentiated from related taxa. Several Tongoloa-like members classified under other genera are interpreted to be closely related to Tongoloa. Morphological and molecular data indicate that Tongoloa, Sinolimprichtia subclade and Chinese Trachydium subclade are separate yet genetically contiguous taxa. It is confirmed that Tongoloa zhongdianensis belongs to the Hymenidium Clade, while Sinocarum is classified within the Acronema Clade. Two new taxa are found in the Hengduan Mountains.

CONCLUSION

Tongoloa is a genus within the East Asia Clade of Apiaceae, and the phylogeny reconstructed based on ITS sequences divides it into 3 main groups. By integrating fruit morphology and molecular phylogenetic analyses, we preliminary clarified the intricate taxonomic relationships among Tongoloa and related taxa. These results provide valuable opportunities for a deeper understanding of the phylogeny of Tongoloa.

Artifactual Orthologs and the Need for Diligent Data Exploration in Complex Phylogenomic Datasets: A Museomic Case Study from the Andean Flora

The Andes mountains of western South America are a globally important biodiversity hotspot, yet there is a paucity of resolved phylogenies for plant clades from this region. Filling an important gap in our understanding of the World's richest flora, we present the first phylogeny of Freziera (Pentaphylacaceae), an Andean-centered, cloud forest radiation. Our dataset was obtained via hybrid-enriched target sequence capture of Angiosperms353 universal loci for 50 of the ca. 75 spp., obtained almost entirely from herbarium specimens. We identify high phylogenomic complexity in Freziera, including the presence of data artifacts. Via by-eye observation of gene trees, detailed examination of warnings from recently improved assembly pipelines, and gene tree filtering, we identified that artifactual orthologs (i.e., the presence of only one copy of a multicopy gene due to differential assembly) were an important source of gene tree heterogeneity that had a negative impact on phylogenetic inference and support. These artifactual orthologs may be common in plant phylogenomic datasets, where multiple instances of genome duplication are common. After accounting for artifactual orthologs as source of gene tree error, we identified a significant, but nonspecific signal of introgression using Patterson's D and f4 statistics. Despite phylogenomic complexity, we were able to resolve Freziera into 9 well-supported subclades whose evolution has been shaped by multiple evolutionary processes, including incomplete lineage sorting, historical gene flow, and gene duplication. Our results highlight the complexities of plant phylogenomics, which are heightened in Andean radiations, and show the impact of filtering data processing artifacts and standard filtering approaches on phylogenetic inference.

The Complex Genetic Legacy of Hybridization and Introgression between the Rare Ocotea loxensis van der Werff and the Widespread O. infrafoveolata van der Werff (Lauraceae)

Hybridization and introgression are complex evolutionary mechanisms that can increase species diversity and lead to speciation, but may also lead to species extinction. In this study, we tested the presence and genetic consequences of hybridization between the rare and Ecuadorian endemic O. loxensis van der Werff and the widespread O. infrafoveolata van der Werff (Lauraceae). Phenotypically, some trees are difficult to identify, and we expect that some might in fact be cryptic hybrids. Thus, we developed nuclear microsatellites to assess the existence of hybrids, as well as the patterns of genetic diversity and population structure in allopatric and sympatric populations. The results revealed high levels of genetic diversity, even in the rare O. loxensis, being usually significantly higher in sympatric than in allopatric populations. The Bayesian assignment of individuals into different genetic classes revealed a complex scenario with different hybrid generations occurring in all sympatric populations, but also in allopatric ones. The absence of some backcrossed hybrids suggests the existence of asymmetric gene flow, and that some hybrids might be more fitted than others might. The existence of current and past interspecific gene flow also explains the blurring of species boundaries in these species and could be linked to the high rates of species found in Ocotea.

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

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

Testing plastomes and nuclear ribosomal DNA sequences as the next-generation DNA barcodes for species identification and phylogenetic analysis in Acer

BACKGROUND

Acer is a taxonomically intractable and speciose genus that contains over 150 species. It is challenging to distinguish Acer species only by morphological method due to their abundant variations. Plastome and nuclear ribosomal DNA (nrDNA) sequences are recommended as powerful next-generation DNA barcodes for species discrimination. However, their efficacies were still poorly studied. The current study will evaluate the application of plastome and nrDNA in species identification and perform phylogenetic analyses for Acer.

RESULT

Based on a collection of 83 individuals representing 55 species (c. 55% of Chinese species) from 13 sections, our barcoding analyses demonstrated that plastomes exhibited the highest (90.47%) species discriminatory power among all plastid DNA markers, such as the standard plastid barcodes matK + rbcL + trnH-psbA (61.90%) and ycf1 (76.19%). And the nrDNA (80.95%) revealed higher species resolution than ITS (71.43%). Acer plastomes show abundant interspecific variations, however, species identification failure may be due to the incomplete lineage sorting (ILS) and chloroplast capture resulting from hybridization. We found that the usage of nrDNA contributed to identifying those species that were unidentified by plastomes, implying its capability to some extent to mitigate the impact of hybridization and ILS on species discrimination. However, combining plastome and nrDNA is not recommended given the cytonuclear conflict caused by potential hybridization. Our phylogenetic analysis covering 19 sections (95% sections of Acer) and 128 species (over 80% species of this genus) revealed pervasive inter- and intra-section cytonuclear discordances, hinting that hybridization has played an important role in the evolution of Acer.

CONCLUSION

Plastomes and nrDNA can significantly improve the species resolution in Acer. Our phylogenetic analysis uncovered the scope and depth of cytonuclear conflict in Acer, providing important insights into its evolution.

Historical biogeography and evolutionary diversification of Lilium (Liliaceae): New insights from plastome phylogenomics

Here, we infer the historical biogeography and evolutionary diversification of the genus Lilium. For this purpose, we used the complete plastomes of 64 currently accepted species in the genus Lilium (14 plastomes were newly sequenced) to recover the phylogenetic backbone of the genus and a time-calibrated phylogenetic framework to estimate biogeographical history scenarios and evolutionary diversification rates of Lilium. Our results suggest that ancient climatic changes and geological tectonic activities jointly shaped the distribution range and drove evolutionary radiation of Lilium, including the Middle Miocene Climate Optimum (MMCO), the late Miocene global cooling, as well as the successive uplift of the Qinghai-Tibet Plateau (QTP) and the strengthening of the monsoon climate in East Asia during the late Miocene and the Pliocene. This case study suggests that the unique geological and climatic events in the Neogene of East Asia, in particular the uplift of QTP and the enhancement of monsoonal climate, may have played an essential role in formation of uneven distribution of plant diversity in the Northern Hemisphere.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

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

Plastid phylogenomics of the tribe potentilleae (Rosaceae)

The tribe Potentilleae comprises approximately 1700 species in 13 genera, making it one of the largest of the 16 tribes in Rosaceae. Our understanding of the composition and relationships among members of Potentilleae has advanced dramatically with the application of molecular markers in the last two decades. Yet there is still much work remaining toward a robust phylogenetic framework for the entire Potentilleae and a comprehensive genus-level dating framework for the tribe. The goals of the present study were to establish a phylogenetic framework for Potentilleae, infer the origin and diversification of the tribe using a temporal framework, and explore the taxonomic implications in light of the updated phylogenetic framework. We used the plastome sequences from 158 accessions representing 139 taxa covering all 13 recognized genera of the tribe to reconstruct the Potentilleae phylogeny. High phylogenetic resolution was recovered along the Potentilleae backbone. Two major clades were recovered within Potentilleae, corresponding to the two subtribes Fragariinae and Potentillinae. Within Fragariinae, two subclades were recovered. In one subclade, Sibbaldia sensu stricto is sister to a clade containing Sibbaldianthe, Comarum, Farinopsis, and Alchemilla sensu lato. In the other subclade, Fragaria is sister to a clade comprising Chamaerhodos, Chamaecallis, Drymocallis, Dasiphora, and Potaninia. Within Potentillinae, Argentina is sister to Potentilla sensu stricto. Within Potentilla sensu stricto, clade Himalaya is sister to Alba, and the Himalaya-Alba clade together is sister to a clade comprising Reptans, Potentilla ancistrifolia Bunge, Fragarioides, Ivesioid, and Argentea. Divergence time estimates indicated that tribe Potentilleae originated during the middle Eocene, and subtribes Fragariinae and Potentillinae diverged around the Eocene-Oligocene transition, and divergence times dated for Potentilleae genera ranged from the early Miocene to the late Pleistocene.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems

Premise

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization.

Methods

We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support.

Results

We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061.

Discussion

Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

The low level of plastome differentiation observed in some lineages of Poales hinders molecular species identification

Chloroplast genomes are a source of information successfully used in various fields of plant genetics, including molecular species identification. However, recent studies indicate an extremely low level of interspecific variability in the plastomes of some taxonomic groups of plants, including the genus Stipa L., which is a representative of the grass family. In this study we aimed to analyze the level of chloroplast genome diversity within particular genera as well as the effectiveness of identifying plant species in the Poaceae family and the other representatives of Poales order. Analysis of complete plastid genome alignments created for 96 genera comprising 793 species and 1707 specimens obtained from the GenBank database allowed defining and categorizing molecular diagnostic characters distinguishing the analyzed species from the other representatives of the genus. The results also demonstrate which species do not have any species-specific mutations, thereby they cannot be identified on the basis of differences between the complete chloroplast genomes. Our research showed a huge diversity of the analyzed species in terms of the number of molecular diagnostic characters and indicated which genera pose a particular challenge in terms of molecular species identification. The results show that a very low level of genetic diversity between plastomes is not uncommon in Poales. This is the first extensive research on super-barcoding that tests this method on a large data set and illustrates its effectiveness against the background of phylogenetic relationships.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

Generic concepts and species diversity within the Gynoxyoid clade (Senecioneae, Compositae)

The Gynoxyoid clade of the Senecioneae (Asteraceae) until now included the five genera Aequatorium, Gynoxys, Nordenstamia, Paracalia and Paragynoxys as diagnosed using selected morphological characters. In their pre-phylogenetic circumscription, the genera Aequatorium and Paragynoxys were considered to inhabit the northern Andes in contrast to Nordenstamia and Paracalia that occur in the central Andes. The most species-rich genus, Gynoxys, was believed to be distributed throughout the Andes. We use a recently established plastid phylogenomic framework that rendered Gynoxys paraphyletic to further evaluate the delimitation of genera in the Gynoxyoid clade. We examine the morphological variation of all members of the Gynoxyoid to identify characters potentially informative at genus level. This results in a matrix of eleven, mostly multistate characters, including those originally used to diagnose these genera. The ancestral character state inference displays a high level of homoplasy, but nevertheless supports the recognition of four genera. Aequatorium is characterised by white radiate capitula. Paracalia and Paragynoxys share white flowers and floral characteristics, such as flower opening and length of disc flowers lobes, as plesiomorphic states, but differ in habit (scandent shrubs vs. trees). Paracalia also retained white flowers, but its two species are characterised by the absence of outer phyllaries. The genera Gynoxys and Nordenstamia comprise species with yellow capitula which appear to be a derived feature in the Gynoxyoids. The genus Nordenstamia, with eight species, is synonymised under Gynoxys since molecular evidence shows its species nested within various parts of the Gynoxys subclade and the morphological variation of Nordenstamia falls well within that of Gynoxys. With the goal to assign all species to four genera (Aequatorium, Gynoxys, Paracalia and Paragynoxys), we assess the states for the eleven characters for all members of the Gynoxyoids and generate new ETS and ITS sequences for 171 specimens belonging to 49 species to further support their generic placement. We provide a taxonomic treatment for the four genera recognised here including amended diagnoses and morphological descriptions. Furthermore, a species-level taxonomic backbone is elaborated for all genera using electronic tools that list 158 currently accepted names and synonyms (209 names in total) with the respective protologue and type information, as well as notes on the current understanding of species limits. Eleven names are newly synonymised, two are lectotypified and eight are newly transferred to other genera.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

Phylogenomic inference of the African tribe Monodoreae (Annonaceae) and taxonomic revision of Dennettia, Uvariodendron and Uvariopsis

Monodoreae (Annonaceae) is a tribe composed of 11 genera and 90 species restricted to the tropical African rain forests. All the genera are taxonomically well circumscribed except the species rich genera Uvariodendron and Uvariopsis which lack a recent taxonomic revision. Here, we used a robust phylogenomic approach, including all the 90 currently accepted species, with several specimens per species, and based on more than 300 Annonaceae-specific nuclear genes, to infer the phylogenetic tree of the Monodoreae and test the limits between the genera and species. We recover all the genera as monophyletic, except the genus Uvariopsis for which the species Uvariopsistripetala falls outside this clade. We thus reinstate the monotypic genus Dennettia for its single species Dennettiatripetala. We also erect a new tribe, Ophrypetaleae trib. nov., to accommodate the genera Ophrypetalum and Sanrafaelia, as we recover them excluded from the Monodoreae tribe with good support. Below the genus level, the genera Isolona, Monodora, Uvariastrum, Uvariodendron and Uvariopsis show weakly supported nodes and phylogenetic conflicts, suggesting that population level processes of evolution might occur in these clades. Our results also support, at the molecular level, the description of several new species of Uvariodendron and Uvariopsis, as well as several new synonymies. Finally, we present a taxonomic revision of the genera Dennettia, Uvariodendron and Uvariopsis, which contain one, 18 and 17 species respectively. We provide a key to the 11 genera of the Monodoraeae and describe four new species to science: Uvariodendronkimbozaense Dagallier & Couvreur, sp. nov., Uvariodendronmossambicense Robson ex Dagallier & Couvreur, sp. nov., Uvariodendronpilosicarpum Dagallier & Couvreur, sp. nov. and Uvariopsisoligocarpa Dagallier & Couvreur, sp. nov., and provide provisional descriptions of three putatively new species. We also present lectotypifications and nomenclatural changes implying synonymies and new combinations (Uvariodendroncitriodorum (Le Thomas) Dagallier & Couvreur, comb. et stat. nov., Uvariodendronfuscumvar.magnificum (Verdc.) Dagallier & Couvreur, comb. et stat. nov., Uvariopsiscongensisvar.angustifolia Dagallier & Couvreur, var. nov., Uvariopsisguineensisvar.globiflora (Keay) Dagallier & Couvreur, comb. et stat. nov., and Uvariopsissolheidiivar.letestui (Pellegr.) Dagallier & Couvreur, comb. et stat. nov.).

Phylotranscriptomics unveil a Paleoproterozoic-Mesoproterozoic origin and deep relationships of the Viridiplantae

The Viridiplantae comprise two main clades, the Chlorophyta (including a diverse array of marine and freshwater green algae) and the Streptophyta (consisting of the freshwater charophytes and the land plants). Lineages sister to core Chlorophyta, informally refer to as prasinophytes, form a grade of mainly planktonic green algae. Recently, one of these lineages, Prasinodermophyta, which is previously grouped with prasinophytes, has been identified as the sister lineage to both Chlorophyta and Streptophyta. Resolving the deep relationships among green plants is crucial for understanding the historical impact of green algal diversity on marine ecology and geochemistry, but has been proven difficult given the ancient timing of the diversification events. Through extensive taxon and gene sampling, we conduct large-scale phylogenomic analyses to resolve deep relationships and reveal the Prasinodermophyta as the lineage sister to Chlorophyta, raising questions about the necessity of classifying the Prasinodermophyta as a distinct phylum. We unveil that incomplete lineage sorting is the main cause of discordance regarding the placement of Prasinodermophyta. Molecular dating analyses suggest that crown-group green plants and crown-group Prasinodermophyta date back to the Paleoproterozoic-Mesoproterozoic. Our study establishes a plausible link between oxygen levels in the Paleoproterozoic-Mesoproterozoic and the origin of Viridiplantae.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Pervasive Introgression During Rapid Diversification of the European Mountain Genus Soldanella (L.) (Primulaceae)

Hybridization is a key mechanism involved in lineage diversification and speciation, especially in ecosystems that experienced repeated environmental oscillations. Recently radiated plant groups, which have evolved in mountain ecosystems impacted by historical climate change provide an excellent model system for studying the impact of gene flow on speciation. We combined organellar (whole-plastome) and nuclear genomic data (RAD-seq) with a cytogenetic approach (rDNA FISH) to investigate the effects of hybridization and introgression on evolution and speciation in the genus Soldanella (snowbells, Primulaceae). Pervasive introgression has already occurred among ancestral lineages of snowbells and has persisted throughout the entire evolutionary history of the genus, regardless of the ecology, cytotype, or distribution range size of the affected species. The highest extent of introgression has been detected in the Carpathian species, which is also reflected in their extensive karyotype variation. Introgression occurred even between species with dysploid and euploid cytotypes, which were considered to be reproductively isolated. The magnitude of introgression detected in snowbells is unprecedented in other mountain genera of the European Alpine System investigated hitherto. Our study stresses the prominent evolutionary role of hybridization in facilitating speciation and diversification on the one hand, but also enriching previously isolated genetic pools. [chloroplast capture; diversification; dysploidy; European Alpine system; introgression; nuclear-cytoplasmic discordance; ribosomal DNA.].

Allopatric speciation is more prevalent than parapatric ecological divergence in a recent high-Andean diversification (Linochilus: Asteraceae)

Elucidating how species accumulate in diversity hotspots is an ongoing debate in evolutionary biology. The páramo, in the Northern Andes, has remarkably high indices of plant diversity, endemicity, and diversification rates. A hypothesis for explaining such indices is that allopatric speciation is high in the páramo given its island-like distribution. An alternative hypothesis is that the altitudinal gradient of the Andean topography provides a variety of niches that drive vertical parapatric ecological speciation. A formal test for evaluating the relative roles of allopatric and parapatric ecological speciation is lacking. The main aim of our study is to test which kind of speciation is more common in an endemic páramo genus. We developed a framework incorporating phylogenetics, species' distributions, and a morpho-ecological trait (leaf area) to compare sister species and infer whether allopatric or parapatric ecological divergence caused their speciation. We applied our framework to the species-rich genus Linochilus (63 spp.) and found that the majority of recent speciation events in it (12 events, 80%) have been driven by allopatric speciation, while a smaller fraction (one event, 6.7%) is attributed to parapatric ecological speciation; two pairs of sister species produced inconclusive results (13.3%). We conclude that páramo autochthonous (in-situ) diversification has been primarily driven by allopatric speciation.

The organellar genomes of Pellidae (Marchantiophyta): the evidence of cryptic speciation, conflicting phylogenies and extraordinary reduction of mitogenomes in simple thalloid liverwort lineage

Organellar genomes of liverworts are considered as one of the most stable among plants, with rare events of gene loss and structural rearrangements. However, not all lineages of liverworts are equally explored in the field of organellar genomics, and subclass Pellidae is one of the less known. Hybrid assembly, using both short- and long-read technologies enabled the assembly of repeat-rich mitogenomes of Pellia and Apopellia revealing extraordinary reduction of length in the latter which impacts only intergenic spacers. The mitogenomes of Apopellia were revealed to be the smallest among all known liverworts-109 k bp, despite retaining all introns. The study also showed the loss of one tRNA gene in Apopellia mitogenome, although it had no impact on the codon usage pattern of mitochondrial protein coding genes. Moreover, it was revealed that Apopellia and Pellia differ in codon usage by plastome CDSs, despite identical tRNA gene content. Molecular identification of species is especially important where traditional taxonomic methods fail, especially within Pellidae where cryptic speciation is well recognized. The simple morphology of these species and a tendency towards environmental plasticity make them complicated in identification. Application of super-barcodes, based on complete mitochondrial or plastid genomes sequences enable identification of all cryptic lineages within Apopellia and Pellia genera, however in some particular cases, mitogenomes were more efficient in species delimitation than plastomes.

Target sequence capture data shed light on the deeper evolutionary relationships of subgenus Chamaecerasus in Lonicera (Caprifoliaceae)

The genus Lonicera L. is widely distributed in the north temperate zone and is well-known for its high species richness and morphological diversity. Previous studies have suggested that many sections of Lonicera are not monophyletic and phylogenetic relationships within the genus are still poorly resolved. In this study, we sampled 37 accessions of Lonicera, covering four sections of subgenus Chamaecerasus plus six outgroup taxa, to recover the main clades of Lonicera based on sequences of nuclear loci generated by target enrichment and cpDNA from genome skimming. We found extensive cytonuclear discordance across the subgenus. Both nuclear and plastid phylogenetic analyses supported subgenus Chamaecerasus sister to subgenus Lonicera. Within subgenus Chamaecerasus, sections Isika and Nintooa were each polyphyletic. Based on the nuclear and chloroplast phylogenies, we propose to merge Lonicera korolkowii into section Coeloxylosteum and Lonicera caerulea into section Nintooa. In addition, Lonicera is estimated to have originated in the mid Oligocene (26.45 Ma). The stem age of section Nintooa was estimated to be 17.09 Ma (95% HPD: 13.30-24.45). The stem age of subgenus Lonicera was estimated to be 16.35 Ma (95% HPD: 14.12-23.66). Ancestral area reconstruction analyses indicate that subgenus Chamaecerasus originated in East Asia and Central Asia. In addition, sections Coeloxylosteum and Nintooa originated in East Asia, with subsequent dispersals into other areas. The aridification of the Asian interior likely promoted the rapid radiation of sections Coeloxylosteum and Nintooa within this region. Moreover, our biogeographic analysis fully supports the Bering and the North Atlantic Land Bridge hypotheses for the intercontinental migrations in the Northern Hemisphere. Overall, this study provides new insights into the taxonomically complex lineages of subgenus Chamaecerasus and the process of speciation.

Climate-influenced boreotropical survival and rampant introgressions explain the thriving of New World grapes in the north temperate zone

The north temperate region was characterized by a warm climate and a rich thermophilic flora before the Eocene, but early diversifications of the temperate biome under global climate change and biome shift remain uncertain. Moreover, it is becoming clear that hybridization/introgression is an important driving force of speciation in plant diversity. Here, we applied analyses from biogeography and phylogenetic networks to account for both introgression and incomplete lineage sorting based on genomic data from the New World Vitis, a charismatic component of the temperate North American flora with known and suspected gene flow among species. Biogeographic inference and fossil evidence suggest that the grapes were widely distributed from North America to Europe during the Paleocene to the Eocene, followed by widespread extinction and survival of relicts in the tropical New World. During the climate warming in the early Miocene, a Vitis ancestor migrated northward from the refugia with subsequent diversification in the North American region. We found strong evidence for widespread incongruence and reticulate evolution among nuclear genes within both recent and ancient lineages of the New World Vitis. Furthermore, the organellar genomes showed strong conflicts with the inferred species tree from the nuclear genomes. Our phylogenomic analyses provided an important assessment of the wide occurrence of reticulate introgression in the New World Vitis, which potentially represents one of the most important mechanisms for the diversification of Vitis species in temperate North America and even the entire temperate Northern Hemisphere. The scenario we report here may be a common model of temperate diversification of flowering plants adapted to the global climate cooling and fluctuation in the Neogene.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Evolution of hygrophytic plant species in the Anatolia-Caucasus region: insights from phylogenomic analyses of Cardamine perennials

BACKGROUND AND AIMS

Southwestern Asia is a significant centre of biodiversity and a cradle of diversification for many plant groups, especially xerophytic elements. In contrast, little is known about the evolution and diversification of its hygrophytic flora. To fill this gap, we focus on Cardamine (Brassicaceae) species that grow in wetlands over a wide altitudinal range. We aimed to elucidate their evolution, assess the extent of presumed historical gene flow between species, and draw inferences about intraspecific structure.

METHODS

We applied the phylogenomic Hyb-Seq approach, ecological niche analyses and multivariate morphometrics to a total of 85 Cardamine populations from the target region of Anatolia-Caucasus, usually treated as four to six species, and supplemented them with close relatives from Europe.

KEY RESULTS

Five diploids are recognized in the focus area, three of which occur in regions adjacent to the Black and/or Caspian Sea (C. penzesii, C. tenera, C. lazica), one species widely distributed from the Caucasus to Lebanon and Iran (C. uliginosa), and one western Anatolian entity (provisionally C. cf. uliginosa). Phylogenomic data suggest recent speciation during the Pleistocene, likely driven by both geographic separation (allopatry) and ecological divergence. With the exception of a single hybrid (allotetraploid) speciation event proven for C. wiedemanniana, an endemic of southern Turkey, no significant traces of past or present interspecific gene flow were observed. Genetic variation within the studied species is spatially structured, suggesting reduced gene flow due to geographic and ecological barriers, but also glacial survival in different refugia.

CONCLUSIONS

This study highlights the importance of the refugial regions of the Black and Caspian Seas for both harbouring and generating hygrophytic species diversity in Southwestern Asia. It also supports the significance of evolutionary links between Anatolia and the Balkan Peninsula. Reticulation and polyploidization played a minor evolutionary role here in contrast to the European relatives.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae

Hybridization is common and considered as an important evolutionary force to increase intraspecific genetic diversity. Detecting hybridization events is crucial for understanding the evolutionary history of species and further improving molecular breeding. The studies on identifying hybridization events through the phylogenomic approach are still limited. We proposed the conception and method of identifying allopolyploidy events by phylogenomics. The reconciliation and summary of nuclear multi-labeled gene family trees were adopted to untangle hybridization events from next-generation data in our novel phylogenomic approach. Given horticulturalists’ relatively clear cultivated crossbreeding history, the water lily family is a suitable case for examining recent allopolyploidy events. Here, we reconstructed and confirmed the well-resolved nuclear phylogeny for the Nymphaeales family in the context of geological time as a framework for identifying hybridization signals. We successfully identified two possible allopolyploidy events with the parental lineages for the hybrids in the family Nymphaeaceae based on summarization from multi-labeled gene family trees of Nymphaeales. The lineages where species Nymphaea colorata and Nymphaea caerulea are located may be the progenitors of horticultural cultivated species Nymphaea ‘midnight’ and Nymphaea ‘Woods blue goddess’. The proposed hybridization hypothesis is also supported by horticultural breeding records. Our methodology can be widely applied to identify hybridization events and theoretically facilitate the genome breeding design of hybrid plants.

Target capture data resolve recalcitrant relationships in the coffee family (Rubioideae, Rubiaceae)

Subfamily Rubioideae is the largest of the main lineages in the coffee family (Rubiaceae), with over 8,000 species and 29 tribes. Phylogenetic relationships among tribes and other major clades within this group of plants are still only partly resolved despite considerable efforts. While previous studies have mainly utilized data from the organellar genomes and nuclear ribosomal DNA, we here use a large number of low-copy nuclear genes obtained via a target capture approach to infer phylogenetic relationships within Rubioideae. We included 101 Rubioideae species representing all but two (the monogeneric tribes Foonchewieae and Aitchinsonieae) of the currently recognized tribes, and all but one non-monogeneric tribe were represented by more than one genus. Using data from the 353 genes targeted with the universal Angiosperms353 probe set we investigated the impact of data type, analytical approach, and potential paralogs on phylogenetic reconstruction. We inferred a robust phylogenetic hypothesis of Rubioideae with the vast majority (or all) nodes being highly supported across all analyses and datasets and few incongruences between the inferred topologies. The results were similar to those of previous studies but novel relationships were also identified. We found that supercontigs [coding sequence (CDS) + non-coding sequence] clearly outperformed CDS data in levels of support and gene tree congruence. The full datasets (353 genes) outperformed the datasets with potentially paralogous genes removed (186 genes) in levels of support but increased gene tree incongruence slightly. The pattern of gene tree conflict at short internal branches were often consistent with high levels of incomplete lineage sorting (ILS) due to rapid speciation in the group. While concatenation- and coalescence-based trees mainly agreed, the observed phylogenetic discordance between the two approaches may be best explained by their differences in accounting for ILS. The use of target capture data greatly improved our confidence and understanding of the Rubioideae phylogeny, highlighted by the increased support for previously uncertain relationships and the increased possibility to explore sources of underlying phylogenetic discordance.

Genomic resolution of cryptic species diversity in chipmunks

Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species diversity within an ecologically well-known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence-capture data for 21,551 autosomal and 493 X-linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow-pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population-level sequence data revealed that the so-called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow-pine chipmunk but does show evidence for recurrent hybridization with this and other species.