Deep phylogenetic incongruence in the angiosperm clade Rosidae

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

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

Plastid phylogenomics and fossil evidence provide new insights into the evolutionary complexity of the 'woody clade' in Saxifragales

BACKGROUND

The "woody clade" in Saxifragales (WCS), encompassing four woody families (Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae), is a phylogenetically recalcitrant node in the angiosperm tree of life, as the interfamilial relationships of the WCS remain contentious. Based on a comprehensive sampling of WCS genera, this study aims to recover a robust maternal backbone phylogeny of the WCS by analyzing plastid genome (plastome) sequence data using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and to explore the possible causes of the phylogenetic recalcitrance with respect to deep relationships within the WCS, in combination with molecular and fossil evidence.

RESULTS

Although the four WCS families were identically resolved as monophyletic, the MP analysis recovered different tree topologies for the relationships among Altingiaceae, Cercidiphyllaceae, and Daphniphyllaceae from the ML and BI phylogenies. The fossil-calibrated plastome phylogeny showed that the WCS underwent a rapid divergence of crown groups in the early Cretaceous (between 104.79 and 100.23 Ma), leading to the origin of the stem lineage ancestors of Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae within a very short time span (∼4.56 Ma). Compared with the tree topology recovered in a previous study based on nuclear genome data, cytonuclear discordance regarding the interfamilial relationships of the WCS was detected.

CONCLUSIONS

Molecular and fossil evidence imply that the early divergence of the WCS might have experienced radiative diversification of crown groups, extensive extinctions at the genus and species levels around the Cretaceous/Paleocene boundary, and ancient hybridization. Such evolutionarily complex events may introduce biases in topological estimations within the WCS due to incomplete lineage sorting, cytonuclear discordance, and long-branch attraction, potentially impacting the accurate reconstruction of deep relationships.

Anatomy of a mega-radiation: Biogeography and niche evolution in Astragalus

PREMISE

Astragalus (Fabaceae), with more than 3000 species, represents a globally successful radiation of morphologically highly similar species predominant across the northern hemisphere. It has attracted attention from systematists and biogeographers, who have asked what factors might be behind the extraordinary diversity of this important arid-adapted clade and what sets it apart from close relatives with far less species richness.

METHODS

Here, for the first time using extensive phylogenetic sampling, we asked whether (1) Astragalus is uniquely characterized by bursts of radiation or whether diversification instead is uniform and no different from closely related taxa. Then we tested whether the species diversity of Astragalus is attributable specifically to its predilection for (2) cold and arid habitats, (3) particular soils, or to (4) chromosome evolution. Finally, we tested (5) whether Astragalus originated in central Asia as proposed and (6) whether niche evolutionary shifts were subsequently associated with the colonization of other continents.

RESULTS

Our results point to the importance of heterogeneity in the diversification of Astragalus, with upshifts associated with the earliest divergences but not strongly tied to any abiotic factor or biogeographic regionalization tested here. The only potential correlate with diversification we identified was chromosome number. Biogeographic shifts have a strong association with the abiotic environment and highlight the importance of central Asia as a biogeographic gateway.

CONCLUSIONS

Our investigation shows the importance of phylogenetic and evolutionary studies of logistically challenging "mega-radiations." Our findings reject any simple key innovation behind high diversity and underline the often nuanced, multifactorial processes leading to species-rich clades.

Characterization of the Complete Mitochondrial Genome of Wintersweet (Chimonanthus praecox) and Comparative Analysis within Magnoliids

Mitochondrial genome sequencing is a valuable tool for investigating mitogenome evolution, species phylogeny, and population genetics. Chimonanthus praecox (L.) Link, also known as "La Mei" in Chinese, is a famous ornamental and medical shrub belonging to the order Laurales of the Calycanthaceae family. Although the nuclear genomes and chloroplast genomes of certain Laurales representatives, such as Lindera glauca, Laurus nobilis, and Piper nigrum, have been sequenced, the mitochondrial genome of Laurales members remains unknown. Here, we reported the first complete mitogenome of C. praecox. The mitogenome was 972,347 bp in length and comprised 60 unique coding genes, including 40 protein-coding genes (PCGs), 17 tRNA genes, and three rRNA genes. The skewness of the PCGs showed that the AT skew (-0.0096233) was negative, while the GC skew (0.031656) was positive, indicating higher contents of T's and G's in the mitochondrial genome of C. praecox. The Ka/Ks ratio analysis showed that the Ka/Ks values of most genes were less than one, suggesting that these genes were under purifying selection. Furthermore, there is a substantial abundance of dispersed repeats in C. praecox, constituting 16.98% of the total mitochondrial genome. A total of 731 SSR repeats were identified in the mitogenome, the highest number among the eleven available magnoliids mitogenomes. The mitochondrial phylogenetic analysis based on 29 conserved PCGs placed the C. praecox in Lauraceae, and supported the sister relationship of Laurales with Magnoliales, which was congruent with the nuclear genome evidence. The present study enriches the mitogenome data of C. praecox and promotes further studies on phylogeny and plastid evolution.

Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: A case study in the cinquefoil genus Potentilla sensu lato (Potentilleae, Rosaceae)

Phylogenomic conflicts are widespread among genomic data, with most previous studies primarily focusing on nuclear datasets instead of organellar genomes. In this study, we investigate phylogenetic conflict analyses within and between plastid and mitochondrial genomes using Potentilla as a case study. We generated three plastid datasets (coding, noncoding, and all-region) and one mitochondrial dataset (coding regions) to infer phylogenies based on concatenated and multispecies coalescent (MSC) methods. Conflict analyses were then performed using PhyParts and Quartet Sampling (QS). Both plastid and mitochondrial genomes divided the Potentilla into eight highly supported clades, two of which were newly identified in this study. While most organellar loci were uninformative for the majority of nodes (bootstrap value < 70%), PhyParts and QS detected conflicting signals within the two organellar genomes. Regression analyses revealed that conflict signals mainly occurred among shorter loci, whereas longer loci tended to be more concordant with the species tree. In addition, two significant disagreements between the two organellar genomes were detected, likely attributed to hybridization and/or incomplete lineage sorting. Our results demonstrate that mitochondrial genes can fully resolve the phylogenetic relationships among eight major clades of Potentilla and are not always linked with plastome in evolutionary history. Stochastic inferences appear to be the primary source of observed conflicts among the gene trees. We recommend that the loci with short sequence length or containing limited informative sites should be used cautiously in MSC analysis, and suggest the joint application of concatenated and MSC methods for phylogenetic inference using organellar genomes.

Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications?

Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.

Assembly and comparative analysis of the complete mitochondrial genome of Isopyrum anemonoides (Ranunculaceae)

Ranunculaceae is a large family of angiosperms comprising 2500 known species-a few with medicinal and ornamental values. Despite this, only two mitochondrial genomes (mitogenomes) of the family have been released in GenBank. Isopyrum anemonoides is a medicinal plant belonging to the family Ranunculaceae, and its chloroplast genome has recently been reported; however, its mitogenome remains unexplored. In this study, we assembled and analyzed the complete mitochondrial genome of I. anemonoides and performed a comparative analysis against different Ranunculaceae species, reconstructing the phylogenetic framework of Isopyrum. The circular mitogenome of I. anemonoides has a length of 206,722 bp, with a nucleotide composition of A (26.4%), T (26.4%), C (23.6%), and G (23.6%), and contains 62 genes, comprising 37 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and three ribosomal RNA (rRNA) genes. Abundantly interspersed repetitive and simple sequence repeat (SSR) loci were detected in the I. anemonoides mitogenome, with tetranucleotide repeats accounting for the highest proportion of SSRs. By detecting gene migration, we observed gene exchange between the chloroplast and mitogenome in I. anemonoides, including six intact tRNA genes, six PCG fragments, and fragments from two rRNA genes. Comparative mitogenome analysis of three Ranunculaceae species indicated that the PCG contents were conserved and the GC contents were similar. Selective pressure analysis revealed that only two genes (nad1 and rpl5) were under positive selection during their evolution in Ranunculales, and two specific RNA editing sites (atp6 and mttB) were detected in the I. anemonoides mitogenome. Moreover, a phylogenetic analysis based on the mitogenomes of I. anemonoides and the other 15 taxa accurately reflected the evolutionary and taxonomic status of I. anemonoides. Overall, this study provides new insights into the genetics, systematics, and evolution of mitochondrial evolution in Ranunculaceae, particularly I. anemonoides.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.

Multiple Outgroups Can Cause Random Rooting in Phylogenomics

Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.

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.

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.

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.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

Incomplete lineage sorting and reticulate evolution mask species relationships in Brunelliaceae, an Andean family with rapid, recent diversification

PREMISE

To date, phylogenetic relationships within the monogeneric Brunelliaceae have been based on morphological evidence, which does not provide sufficient phylogenetic resolution. Here we use target-enriched nuclear data to improve our understanding of phylogenetic relationships in the family.

METHODS

We used the Angiosperms353 toolkit for targeted recovery of exonic regions and supercontigs (exons + introns) from low copy nuclear genes from 53 of 70 species in Brunellia, and several outgroup taxa. We removed loci that indicated biased inference of relationships and applied concatenated and coalescent methods to infer Brunellia phylogeny. We identified conflicts among gene trees that may reflect hybridization or incomplete lineage sorting events and assessed their impact on phylogenetic inference. Finally, we performed ancestral-state reconstructions of morphological traits and assessed the homology of character states used to define sections and subsections in Brunellia.

RESULTS

Brunellia comprises two major clades and several subclades. Most of these clades/subclades do not correspond to previous infrageneric taxa. There is high topological incongruence among the subclades across analyses.

CONCLUSIONS

Phylogenetic reconstructions point to rapid species diversification in Brunelliaceae, reflected in very short branches between successive species splits. The removal of putatively biased loci slightly improves phylogenetic support for individual clades. Reticulate evolution due to hybridization and/or incomplete lineage sorting likely both contribute to gene-tree discordance. Morphological characters used to define taxa in current classification schemes are homoplastic in the ancestral character-state reconstructions. While target enrichment data allows us to broaden our understanding of diversification in Brunellia, the relationships among subclades remain incompletely understood.

High-quality chromosome-level genome assembly of Litsea coreana L. provides insights into Magnoliids evolution and flavonoid biosynthesis

The magnoliid Litsea coreana has been the subject of a substantial amount of research owing to its production of many flavonoid metabolites, high food processing value, and a controversial phylogenetic position. For this study, we assembled a high-grade genome at the chromosome scale and annotation of L. coreana that was anchored to 12 chromosomes. The total genome was 1139.45 Mb, while the N50 scaffold was 97.18 Mb long. The analysis of phylogenetic trees constructed by different methods show that the phylogeny of Magnoliids is inconsistent, indicating that the differentiation process of monocots, eudicots, and Magnoliids still remains in dispute. An ancient whole-genome duplication (WGD) event was shown to have occurred before the Magnoliales and Laurels had differentiated. Subsequently, an independent WGD appeared in the Lauralean lineage. A total of 27 types of flavonoids were detected in all five tissues of L. coreana. Chalcone synthases (CHSs) that are responsible for production of flavonoids have been validated at the bioinformatics level. The retention of comparative genomic analyses of the CHS gene family showed that this family had contracted significantly in L. coreana. Our research further elaborated the evolution of Lauraceae and perfected the genetic basis of flavonoid biosynthesis in L. coreana. SIGNIFICANCE STATEMENT: Provides evidence that determines the evolutionary status of Magnoliids. The chalcone synthase gene family was significantly contracted in Litsea coreana.

Incompatibility Phylogenetic Signals between Double-Digest Restriction Site-Associated DNA Sequencing and Plastid Genomes in Chinese Curcuma (Zingiberaceae)—A Recent Qinghai–Tibetan Plateau Diversification Genera

Curcuma is of high economic value, credited to its medicinal, edible, and ornamental properties, which possess all signatures of adaptability, and rapid radiation, especially species of Curcuma (Chinese Curcuma, a recent Qinghai–Tibetan Plateau diversification genera) scattered in China. However, little is known about the incongruent phylogenetic signals within this genera from different inheritance patterns that will militate against the further development of this genera. In this research, we applied complete chloroplast genome data together with double-digest restriction site-associated DNA sequencing data (ddRAD-seq) strategy to investigate phylogenetic signals of Chinese Curcuma species, clustering using two RAD analysis pipelines (STACKS and pyRAD). Phylogenetic trees were obtained from each locus based on the maximum likelihood (ML) and multispecies coalescent (BEAST) methods. For visual comparison, multi-method and different datasets were used to infer the phylogeny. We discovered inconsistent relationships for the Chinese Curcuma with varying degrees of support using different methods and datasets.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

Plastome Diversity and Phylogenomic Relationships in Asteraceae

Plastid genomes are in general highly conserved given their slow evolutionary rate, and thus large changes in their structure are unusual. However, when specific rearrangements are present, they are often phylogenetically informative. Asteraceae is a highly diverse family whose evolution is long driven by polyploidy (up to 48x) and hybridization, both processes usually complicating systematic inferences. In this study, we generated one of the most comprehensive plastome-based phylogenies of family Asteraceae, providing information about the structure, genetic diversity and repeat composition of these sequences. By comparing the whole-plastome sequences obtained, we confirmed the double inversion located in the long single-copy region, for most of the species analyzed (with the exception of basal tribes), a well-known feature for Asteraceae plastomes. We also showed that genome size, gene order and gene content are highly conserved along the family. However, species representative of the basal subfamily Barnadesioideae-as well as in the sister family Calyceraceae-lack the pseudogene rps19 located in one inverted repeat. The phylogenomic analysis conducted here, based on 63 protein-coding genes, 30 transfer RNA genes and 21 ribosomal RNA genes from 36 species of Asteraceae, were overall consistent with the general consensus for the family's phylogeny while resolving the position of tribe Senecioneae and revealing some incongruences at tribe level between reconstructions based on nuclear and plastid DNA data.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

Comparative analysis of plastomes in Oxalidaceae: Phylogenetic relationships and potential molecular markers

The wood sorrel family, Oxalidaceae, is mainly composed of annual or perennial herbs, a few shrubs, and trees distributed from temperate to tropical zones. Members of Oxalidaceae are of high medicinal, ornamental, and economic value. Despite the rich diversity and value of Oxalidaceae, few molecular markers or plastomes are available for phylogenetic analysis of the family. Here, we reported four new whole plastomes of Oxalidaceae and compared them with plastomes of three species in the family, as well as the plastome of Rourea microphylla in the closely related family Connaraceae. The eight plastomes ranged in length from 150,673 bp (Biophytum sensitivum) to 156,609 bp (R. microphylla). Genome annotations revealed a total of 129-131 genes, including 83-84 protein-coding genes, eight rRNA genes, 37 tRNA genes, and two to three pseudogenes. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. The smaller plastomes of herbs B. sensitivum and three Oxalis species are associated with variations in IR region sizes, intergenic region variation, and gene or intron loss. We identified sequences with high variation that may serve as molecular markers in taxonomic studies of Oxalidaceae. The phylogenetic trees of selected superrosid representatives based on 76 protein-coding genes corroborated the Oxalidaceae position in Oxalidales and supported it as a sister to Connaraceae. Our research also supported the monophyly of the COM (Celastrales, Oxalidales, and Malpighiales) clade.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Phylogenomics and biogeography of Cunoniaceae (Oxalidales) with complete generic sampling and taxonomic realignments

PREMISE

Cunoniaceae are a family of shrubs and trees with 27 genera and ca. 335 species, mostly confined to tropical and wet temperate zones of the southern hemisphere. There are several known issues regarding generic limits, and the family also displays a number of intriguing long-range disjunctions.

METHODS

We performed a phylogenomic study using the universal Angiosperms353 probe set for targeted sequence capture. We sampled 37 species covering all genera in the Cunoniaceae, and those in the three closely related families of the crown Oxalidales (Brunelliaceae, Cephalotaceae, and Elaeocarpaceae). We also performed analyses for molecular dating and ancestral area reconstruction.

RESULTS

We recovered the topology (Cunoniaceae, (Cephalotaceae, (Brunelliaceae, Elaeocarpaceae))) and a well-resolved genus-level phylogeny of Cunoniaceae with strongly supported clades corresponding to all previously recognized tribes. As previously suspected, the genera Ackama and Weinmannia were recovered as paraphyletic. Australasia was inferred as the likely ancestral area for the family.

CONCLUSIONS

The current distribution of Cunoniaceae is best explained by long-distance dispersal with a few possible cases of Australasian-American vicariance events. Extinctions may have been important in determining the mostly Oceanian distribution of this family while some genera in the tribe Cunonieae and in New Caledonia have undergone recent bursts of diversification. New generic diagnoses, 80 new combinations, and one new name are provided for a recircumscribed Ackama (including Spiraeopsis), a much smaller Weinmannia (mostly New World), and a resurrected Pterophylla to accommodate Old World taxa previously in Weinmannia.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

Biogeographical patterns and speciation of the genus Pinguicula (Lentibulariaceae) inferred by phylogenetic analyses

Earlier phylogenetic studies in the genus Pinguicua (Lentibulariaceae) suggested that the species within a geographical region was rather monophyletic, although the sampling was limited or was restricted to specific regions. Those results conflicted with the floral morphology-based classification, which has been widely accepted to date. In the current study, one nuclear ribosomal DNA (internal transcribed spacer; ITS) and two regions of chloroplast DNA (matK and rpl32-trnL), from up to ca. 80% of the taxa in the genus Pinguicula, covering all three subgenera, were sequenced to demonstrate the inconsistency and explore a possible evolutionary history of the genus. Some incongruence was observed between nuclear and chloroplast topologies and the results from each of the three DNA analyses conflicted with the morphology-based subgeneric divisions. Both the ITS tree and network, however, corresponded with the biogeographical patterns of the genus supported by life-forms (winter rosette or hibernaculum formation) and basic chromosome numbers (haploidy). The dormant strategy evolved in a specific geographical region is a phylogenetic constraint and a synapomorphic characteristic within a lineage. Therefore, the results denied the idea that the Mexican group, morphologically divided into the three subgenera, independently acquired winter rosette formations. Topological incongruence among the trees or reticulations, indicated by parallel edges in phylogenetic networks, implied that some taxa originated by introgressive hybridisation. Although there are exceptions, species within the same geographical region arose from a common ancestor. Therefore, the classification by the floral characteristics is rather unreliable. The results obtained from this study suggest that evolution within the genus Pinguicula has involved; 1) ancient expansions to geographical regions with gene flow and subsequent vicariance with genetic drift, 2) acquirement of a common dormant strategy within a specific lineage to adapt a local climate (i.e., synapomorphic characteristic), 3) recent speciation in a short time span linked to introgressive hybridisation or multiplying the ploidy level (i.e., divergence), and 4) parallel evolution in floral traits among lineages found in different geographical regions (i.e., convergence). As such, the floral morphology masks and obscures the phylogenetic relationships among species in the genus.

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Genomic evidence of prevalent hybridization throughout the evolutionary history of the fig-wasp pollination mutualism

Ficus (figs) and their agaonid wasp pollinators present an ecologically important mutualism that also provides a rich comparative system for studying functional co-diversification throughout its coevolutionary history (~75 million years). We obtained entire nuclear, mitochondrial, and chloroplast genomes for 15 species representing all major clades of Ficus. Multiple analyses of these genomic data suggest that hybridization events have occurred throughout Ficus evolutionary history. Furthermore, cophylogenetic reconciliation analyses detect significant incongruence among all nuclear, chloroplast, and mitochondrial-based phylogenies, none of which correspond with any published phylogenies of the associated pollinator wasps. These findings are most consistent with frequent host-switching by the pollinators, leading to fig hybridization, even between distantly related clades. Here, we suggest that these pollinator host-switches and fig hybridization events are a dominant feature of fig/wasp coevolutionary history, and by generating novel genomic combinations in the figs have likely contributed to the remarkable diversity exhibited by this mutualism.

The Deep Evolutionary Relationships of the Morphologically Heterogeneous Nolinoideae (Asparagaceae) Revealed by Transcriptome Data

The subfamily Nolinoideae of Asparagaceae is an extremely morphologically heterogeneous group, which is comprised of seven lineages, formerly known as Eriospermaceae, Polygonateae, Ophiopogoneae, Convallarieae, Ruscaceae s.s., Dracaenaceae, and Nolinaceae from different families or even orders. Their drastically divergent morphologies and low level of molecular resolution have hindered our understanding on their evolutionary history. To resolve reliable and clear phylogenetic relationships of the Nolinoideae, a phylogenetic study was conducted based on transcriptomic sequencing of 15 species representing all the seven lineages. A dataset containing up to 2,850,331 sites across 2,126 genes was analyzed using both concatenated and coalescent methods. Except for Eriospermum as outgroup, the transcriptomic data strongly resolved the remaining six lineages into two groups, one is a paraphyletic grade including the woody lineages of dracaenoids, ruscoids, and nolinoids and a monophyletic herbaceous clade. Within the herbaceous group, the Ophiopogoneae + Theropogon is sister to a clade that is composed of Convallarieae and the monophyletic Polygonateae. Our work provides a first robust deep relationship of the highly heterogeneous Nolinoideae and paves the way for further investigations of its complex evolution.

Comparative Performance of Popular Methods for Hybrid Detection using Genomic Data

Interspecific hybridization is an important evolutionary phenomenon that generates genetic variability in a population and fosters species diversity in nature. The availability of large genome scale datasets has revolutionized hybridization studies to shift from the observation of the presence or absence of hybrids to the investigation of the genomic constitution of hybrids and their genome-specific evolutionary dynamics. Although a handful of methods have been proposed in an attempt to identify hybrids, accurate detection of hybridization from genomic data remains a challenging task. In addition to methods that infer phylogenetic networks or that utilize pairwise divergence, site pattern frequency based and population genetic clustering approaches are popularly used in practice, though the performance of these methods under different hybridization scenarios has not been extensively examined. Here, we use simulated data to comparatively evaluate the performance of four tools that are commonly used to infer hybridization events: the site pattern frequency based methods HyDe and the D-statistic (i.e., the ABBA-BABA test) and the population clustering approaches structure and ADMIXTURE. We consider single hybridization scenarios that vary in the time of hybridization and the amount of incomplete lineage sorting (ILS) for different proportions of parental contributions (γ); introgressive hybridization; multiple hybridization scenarios; and a mixture of ancestral and recent hybridization scenarios. We focus on the statistical power to detect hybridization and the false discovery rate (FDR) for comparisons of the D-statistic and HyDe, and the accuracy of the estimates of γ as measured by the mean squared error for HyDe, structure, and ADMIXTURE. Both HyDe and the D-statistic are powerful for detecting hybridization in all scenarios except those with high ILS, although the D-statistic often has an unacceptably high FDR. The estimates of γ in HyDe are impressively robust and accurate whereas structure and ADMIXTURE sometimes fail to identify hybrids, particularly when the proportional parental contributions are asymmetric (i.e., when γ is close to 0). Moreover, the posterior distribution estimated using structure exhibits multimodality in many scenarios, making interpretation difficult. Our results provide guidance in selecting appropriate methods for identifying hybrid populations from genomic data.

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

How to Tackle Phylogenetic Discordance in Recent and Rapidly Radiating Groups? Developing a Workflow Using Loricaria (Asteraceae) as an Example

A major challenge in phylogenetics and -genomics is to resolve young rapidly radiating groups. The fast succession of species increases the probability of incomplete lineage sorting (ILS), and different topologies of the gene trees are expected, leading to gene tree discordance, i.e., not all gene trees represent the species tree. Phylogenetic discordance is common in phylogenomic datasets, and apart from ILS, additional sources include hybridization, whole-genome duplication, and methodological artifacts. Despite a high degree of gene tree discordance, species trees are often well supported and the sources of discordance are not further addressed in phylogenomic studies, which can eventually lead to incorrect phylogenetic hypotheses, especially in rapidly radiating groups. We chose the high-Andean Asteraceae genus Loricaria to shed light on the potential sources of phylogenetic discordance and generated a phylogenetic hypothesis. By accounting for paralogy during gene tree inference, we generated a species tree based on hundreds of nuclear loci, using Hyb-Seq, and a plastome phylogeny obtained from off-target reads during target enrichment. We observed a high degree of gene tree discordance, which we found implausible at first sight, because the genus did not show evidence of hybridization in previous studies. We used various phylogenomic analyses (trees and networks) as well as the D-statistics to test for ILS and hybridization, which we developed into a workflow on how to tackle phylogenetic discordance in recent radiations. We found strong evidence for ILS and hybridization within the genus Loricaria. Low genetic differentiation was evident between species located in different Andean cordilleras, which could be indicative of substantial introgression between populations, promoted during Pleistocene glaciations, when alpine habitats shifted creating opportunities for secondary contact and hybridization.

Resolving the Phylogeny of the Olive Family (Oleaceae): Confronting Information from Organellar and Nuclear Genomes

The olive family, Oleaceae, is a group of woody plants comprising 28 genera and ca. 700 species, distributed on all continents (except Antarctica) in both temperate and tropical environments. It includes several genera of major economic and ecological importance such as olives, ash trees, jasmines, forsythias, osmanthuses, privets and lilacs. The natural history of the group is not completely understood yet, but its diversification seems to be associated with polyploidisation events and the evolution of various reproductive and dispersal strategies. In addition, some taxonomical issues still need to be resolved, particularly in the paleopolyploid tribe Oleeae. Reconstructing a robust phylogenetic hypothesis is thus an important step toward a better comprehension of Oleaceae's diversity. Here, we reconstructed phylogenies of the olive family using 80 plastid coding sequences, 37 mitochondrial genes, the complete nuclear ribosomal cluster and a small multigene family encoding phytochromes (phyB and phyE) of 61 representative species. Tribes and subtribes were strongly supported by all phylogenetic reconstructions, while a few Oleeae genera are still polyphyletic (Chionanthus, Olea, Osmanthus, Nestegis) or paraphyletic (Schrebera, Syringa). Some phylogenetic relationships among tribes remain poorly resolved with conflicts between topologies reconstructed from different genomic regions. The use of nuclear data remains an important challenge especially in a group with ploidy changes (both paleo- and neo-polyploids). This work provides new genomic datasets that will assist the study of the biogeography and taxonomy of the whole Oleaceae.

The Perfect Storm: Gene Tree Estimation Error, Incomplete Lineage Sorting, and Ancient Gene Flow Explain the Most Recalcitrant Ancient Angiosperm Clade, Malpighiales

The genomic revolution offers renewed hope of resolving rapid radiations in the Tree of Life. The development of the multispecies coalescent model and improved gene tree estimation methods can better accommodate gene tree heterogeneity caused by incomplete lineage sorting (ILS) and gene tree estimation error stemming from the short internal branches. However, the relative influence of these factors in species tree inference is not well understood. Using anchored hybrid enrichment, we generated a data set including 423 single-copy loci from 64 taxa representing 39 families to infer the species tree of the flowering plant order Malpighiales. This order includes 9 of the top 10 most unstable nodes in angiosperms, which have been hypothesized to arise from the rapid radiation during the Cretaceous. Here, we show that coalescent-based methods do not resolve the backbone of Malpighiales and concatenation methods yield inconsistent estimations, providing evidence that gene tree heterogeneity is high in this clade. Despite high levels of ILS and gene tree estimation error, our simulations demonstrate that these two factors alone are insufficient to explain the lack of resolution in this order. To explore this further, we examined triplet frequencies among empirical gene trees and discovered some of them deviated significantly from those attributed to ILS and estimation error, suggesting gene flow as an additional and previously unappreciated phenomenon promoting gene tree variation in Malpighiales. Finally, we applied a novel method to quantify the relative contribution of these three primary sources of gene tree heterogeneity and demonstrated that ILS, gene tree estimation error, and gene flow contributed to 10.0$\%$, 34.8$\%$, and 21.4$\%$ of the variation, respectively. Together, our results suggest that a perfect storm of factors likely influence this lack of resolution, and further indicate that recalcitrant phylogenetic relationships like the backbone of Malpighiales may be better represented as phylogenetic networks. Thus, reducing such groups solely to existing models that adhere strictly to bifurcating trees greatly oversimplifies reality, and obscures our ability to more clearly discern the process of evolution. [Coalescent; concatenation; flanking region; hybrid enrichment, introgression; phylogenomics; rapid radiation, triplet frequency.].

Conflicting phylogenetic signals in plastomes of the tribe Laureae (Lauraceae)

Background

Gene tree discordance is common in phylogenetic analyses. Many phylogenetic studies have excluded non-coding regions of the plastome without evaluating their impact on tree topology. In general, plastid loci have often been treated as a single unit, and tree discordance among these loci has seldom been examined. Using samples of Laureae (Lauraceae) plastomes, we explored plastome variation among the tribe, examined the influence of non-coding regions on tree topology, and quantified intra-plastome conflict.

Results

We found that the plastomes of Laureae have low inter-specific variation and are highly similar in structure, size, and gene content. Laureae was divided into three groups, subclades I, II and III. The inclusion of non-coding regions changed the phylogenetic relationship among the three subclades. Topologies based on coding and non-coding regions were largely congruent except for the relationship among subclades I, II and III. By measuring the distribution of phylogenetic signal across loci that supported different topologies, we found that nine loci (two coding regions, two introns and five intergenic spacers) played a critical role at the contentious node.

Conclusions

Our results suggest that subclade III and subclade II are successively sister to subclade I. Conflicting phylogenetic signals exist between coding and non-coding regions of Laureae plastomes. Our study highlights the importance of evaluating the influence of non-coding regions on tree topology and emphasizes the necessity of examining discordance among different plastid loci in phylogenetic studies.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

Complex evolutionary history of coffees revealed by full plastid genomes and 28,800 nuclear SNP analyses, with particular emphasis on Coffea canephora (Robusta coffee)

For decades coffees were associated with the genus Coffea. In 2011, the closely related genus Psilanthus was subsumed into Coffea. However, results obtained in 2017-based on 28,800 nuclear SNPs-indicated that there is not substantial phylogenetic support for this incorporation. In addition, a recent study of 16 plastid full-genome sequences highlighted an incongruous placement of Coffea canephora (Robusta coffee) between maternal and nuclear trees. In this study, similar global features of the plastid genomes of Psilanthus and Coffea are observed. In agreement with morphological and physiological traits, the nuclear phylogenetic tree clearly separates Psilanthus from Coffea (with exception to C. rhamnifolia, closer to Psilanthus than to Coffea). In contrast, the maternal molecular tree was incongruent with both morphological and nuclear differentiation, with four main clades observed, two of which include both Psilanthus and Coffea species, and two with either Psilanthus or Coffea species. Interestingly, Coffea and Psilanthus taxa sampled in West and Central Africa are members of the same group. Several mechanisms such as the retention of ancestral polymorphisms due to incomplete lineage sorting, hybridization leading to homoploidy (without chromosome doubling) and alloploidy (for C. arabica) are involved in the evolutionary history of the coffee species. While sharing similar morphological characteristics, the genetic relationships within C. canephora have shown that some populations are well differentiated and genetically isolated. Given the position of its closely-related species, we may also consider C. canephora to be undergoing a long process of speciation with an intermediate step of (sub-)speciation.

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

The Litsea genome and the evolution of the laurel family

The laurel family within the Magnoliids has attracted attentions owing to its scents, variable inflorescences, and controversial phylogenetic position. Here, we present a chromosome-level assembly of the Litsea cubeba genome, together with low-coverage genomic and transcriptomic data for many other Lauraceae. Phylogenomic analyses show phylogenetic discordance at the position of Magnoliids, suggesting incomplete lineage sorting during the divergence of monocots, eudicots, and Magnoliids. An ancient whole-genome duplication (WGD) event occurred just before the divergence of Laurales and Magnoliales; subsequently, independent WGDs occurred almost simultaneously in the three Lauralean lineages. The phylogenetic relationships within Lauraceae correspond to the divergence of inflorescences, as evidenced by the phylogeny of FUWA, a conserved gene involved in determining panicle architecture in Lauraceae. Monoterpene synthases responsible for production of specific volatile compounds in Lauraceae are functionally verified. Our work sheds light on the evolution of the Lauraceae, the genetic basis for floral evolution and specific scents.

Litsea cubeba belongs to the Lauraceae family within the Magnoliids clade. Here, the authors assemble its genome and reveal divergence of inflorescence and sexual differentiation, the phylogenetic relationships across the Lauraceae and related species, and biosynthetic genes related to essential oil synthesis.

Prickly waterlily and rigid hornwort genomes shed light on early angiosperm evolution

Angiosperms represent one of the most spectacular terrestrial radiations on the planet1, but their early diversification and phylogenetic relationships remain uncertain2-5. A key reason for this impasse is the paucity of complete genomes representing early-diverging angiosperms. Here, we present high-quality, chromosomal-level genome assemblies of two aquatic species-prickly waterlily (Euryale ferox; Nymphaeales) and the rigid hornwort (Ceratophyllum demersum; Ceratophyllales)-and expand the genomic representation for key sectors of the angiosperm tree of life. We identify multiple independent polyploidization events in each of the five major clades (that is, Nymphaeales, magnoliids, monocots, Ceratophyllales and eudicots). Furthermore, our phylogenomic analyses, which spanned multiple datasets and diverse methods, confirm that Amborella and Nymphaeales are successively sister to all other angiosperms. Furthermore, these genomes help to elucidate relationships among the major subclades within Mesangiospermae, which contain about 350,000 species. In particular, the species-poor lineage Ceratophyllales is supported as sister to eudicots, and monocots and magnoliids are placed as successively sister to Ceratophyllales and eudicots. Finally, our analyses indicate that incomplete lineage sorting may account for the incongruent phylogenetic placement of magnoliids between nuclear and plastid genomes.

-Complete plastid genome sequences of two species of the Neotropical genus Brunellia (Brunelliaceae)

Here we present the first two complete plastid genomes for Brunelliaceae, a Neotropical family with a single genus, Brunellia. We surveyed the entire plastid genome in order to find variable cpDNA regions for further phylogenetic analyses across the family. We sampled morphologically different species, B. antioquensis and B. trianae, and found that the plastid genomes are 157,685 and 157,775 bp in length and display the typical quadripartite structure found in angiosperms. Despite the clear morphological distinction between both species, the molecular data show a very low level of divergence. The amount of nucleotide substitutions per site is one of the lowest reported to date among published congeneric studies (π = 0.00025). The plastid genomes have gene order and content coincident with other COM (Celastrales, Oxalidales, Malpighiales) relatives. Phylogenetic analyses of selected superrosid representatives show high bootstrap support for the ((C,M)O) topology. The N-fixing clade appears as the sister group of the COM clade and Zygophyllales as the sister to the rest of the fabids group.

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BACKGROUND

Paris (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.

RESULTS

All 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.

CONCLUSIONS

Our results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

A Universal Probe Set for Targeted Sequencing of 353 Nuclear Genes from Any Flowering Plant Designed Using k-Medoids Clustering

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5-15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself.

Multi-gene incongruence consistent with hybridisation in Cladocopium (Symbiodiniaceae), an ecologically important genus of coral reef symbionts

Coral reefs rely on their intracellular dinoflagellate symbionts (family Symbiodiniaceae) for nutritional provision in nutrient-poor waters, yet this association is threatened by thermally stressful conditions. Despite this, the evolutionary potential of these symbionts remains poorly characterised. In this study, we tested the potential for divergent Symbiodiniaceae types to sexually reproduce (i.e. hybridise) within Cladocopium, the most ecologically prevalent genus in this family. With sequence data from three organelles (cob gene, mitochondrion; psbA^ncr region, chloroplast; and ITS2 region, nucleus), we utilised the Incongruence Length Difference test, Approximately Unbiased test, tree hybridisation analyses and visual inspection of raw data in stepwise fashion to highlight incongruences between organelles, and thus provide evidence of reticulate evolution. Using this approach, we identified three putative hybrid Cladocopium samples among the 158 analysed, at two of the seven sites sampled. These samples were identified as the common Cladocopium types C40 or C1 with respect to the mitochondria and chloroplasts, but the rarer types C3z, C3u and C1# with respect to their nuclear identity. These five Cladocopium types have previously been confirmed as evolutionarily distinct and were also recovered in non-incongruent samples multiple times, which is strongly suggestive that they sexually reproduced to produce the incongruent samples. A concomitant inspection of next generation sequencing data for these samples suggests that other plausible explanations, such as incomplete lineage sorting or the presence of co-dominance, are much less likely. The approach taken in this study allows incongruences between gene regions to be identified with confidence, and brings new light to the evolutionary potential within Symbiodiniaceae.