Phylogenomics and divergence pattern of Polygonatum (Asparagaceae: Polygonateae) in the north temperate region.

Polygonatum is the largest genus of tribe Polygonateae (Asparagaceae) and is widely distributed in the temperate Northern Hemisphere, especially well diversified in southwestern China to northeastern Asia. Phylogenetic relationships of many species are still controversial. Hence it is necessary to clarify their phylogenetic relationships and infer possible reticulate relationships for the genus. In this study, genome-wide data of 43 species from Polygonatum and its closely related taxa were obtained by Hyb-Seq sequencing. The phylogenetic trees constructed from genome-wide nuclear and chloroplast sequences strongly supported the monophyly of Polygonatum with division into three major clades. A high level of incongruence was detected between nuclear and chloroplast trees as well as among gene trees within the genus, but all occurred within each major clade. However, introgression tests and reticulate evolution analyses revealed low level of gene flow and weak introgression events in the genus, suggesting hybridization and introgression were not dominant during the evolutionary diversification of Polygonatum in the Northern Hemisphere. This study provides important insights into reconstructing evolutionary relationships and speciation pattern of taxa from the north temperate flora.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications.

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

The sweet tabaiba or there and back again: phylogeographical history of the Macaronesian Euphorbia balsamifera.

BACKGROUND AND AIMS: Biogeographical relationships between the Canary Islands and north-west Africa are often explained by oceanic dispersal and geographical proximity. Sister-group relationships between Canarian and eastern African/Arabian taxa, the 'Rand Flora' pattern, are rare among plants and have been attributed to the extinction of north-western African populations. Euphorbia balsamifera is the only representative species of this pattern that is distributed in the Canary Islands and north-west Africa; it is also one of few species present in all seven islands. Previous studies placed African populations of E. balsamifera as sister to the Canarian populations, but this relationship was based on herbarium samples with highly degraded DNA. Here, we test the extinction hypothesis by sampling new continental populations; we also expand the Canarian sampling to examine the dynamics of island colonization and diversification. METHODS: Using target enrichment with genome skimming, we reconstructed phylogenetic relationships within E. balsamifera and between this species and its disjunct relatives. A single nucleotide polymorphism dataset obtained from the target sequences was used to infer population genetic diversity patterns. We used convolutional neural networks to discriminate among alternative Canary Islands colonization scenarios. KEY RESULTS: The results confirmed the Rand Flora sister-group relationship between western E. balsamifera and Euphorbia adenensis in the Eritreo-Arabian region and recovered an eastern-western geographical structure among E. balsamifera Canarian populations. Convolutional neural networks supported a scenario of east-to-west island colonization, followed by population extinctions in Lanzarote and Fuerteventura and recolonization from Tenerife and Gran Canaria; a signal of admixture between the eastern island and north-west African populations was recovered. CONCLUSIONS: Our findings support the Surfing Syngameon Hypothesis for the colonization of the Canary Islands by E. balsamifera, but also a recent back-colonization to the continent. Populations of E. balsamifera from northwest Africa are not the remnants of an ancestral continental stock, but originated from migration events from Lanzarote and Fuerteventura. This is further evidence that oceanic archipelagos are not a sink for biodiversity, but may be a source of new genetic variability.

Next-generation sequencing data show rapid radiation and several long-distance dispersal events in early Costaceae.

The monocot family Costaceae Nakai consists of seven genera but their mutual relationships have not been satisfactorily resolved in previous studies employing classical molecular markers. Phylogenomic analyses of 365 nuclear genes and nearly-complete plastome data provide almost fully resolved insights into their diversification. Paracostus is identified as sister to all other taxa, followed by several very short branches leading to discrete lineages, suggesting an ancient rapid radiation of these early lineages and leaving the exact relationships among them unresolved. Relationships among Chamaecostus, Dimerocostus and Monocostus confirmed earlier findings that these genera form a monophyletic group. The Afro-American Costus is also monophyletic. By contrast, Tapeinochilos appeared as a well-supported crown lineage of Cheilocostus rendering it paraphyletic. As these two genera differ morphologically from one another owing to a shift from insect- to bird-pollination, we propose to keep both names. The divergence time within Costaceae was estimated using penalized likelihood utilizing two fossils within Zingiberales, †Spirematospermum chandlerae and †Ensete oregonense, indicated a relatively recent diversification of Costaceae, between 18 and 9 Mya. Based on these data, the current pantropical distribution of the family is hypothesized to be the result of several long-distance intercontinental dispersal events, which do not correlate with global geoclimatic changes.

Phylogenomics and genome size evolution in Amomum s. s. (Zingiberaceae): Comparison of traditional and modern sequencing methods.

BACKGROUND AND AIMS: A targeted enrichment NGS approach was used to construct the phylogeny of Amomum Roxb. (Zingiberaceae). Phylogenies based on hundreds of nuclear genes, the whole plastome and the rDNA cistron were compared with an ITS-based phylogeny. Trends in genome size (GS) evolution were examined, chromosomes were counted and the geographical distribution of phylogenetic lineages was evaluated. METHODS: In total, 92 accessions of 54 species were analysed. ITS was obtained for 79 accessions, 37 accessions were processed with Hyb-Seq and sequences from 449 nuclear genes, the whole cpDNA, and the rDNA cistron were analysed using concatenation, coalescence and supertree approaches. The evolution of absolute GS was analysed in a phylogenetic and geographical context. The chromosome numbers of 12 accessions were counted. KEY RESULTS: Four groups were recognised in all datasets though their mutual relationships differ among datasets. While group A (A. subulatum and A. petaloideum) is basal to the remaining groups in the nuclear gene phylogeny, in the cpDNA topology it is sister to group B (A. repoeense and related species) and, in the ITS topology, it is sister to group D (the Elettariopsis lineage). The former Elettariopsis makes a monophyletic group. There is an increasing trend in GS during evolution. The largest GS values were found in group D in two tetraploid taxa, A. cinnamomeum and A. aff. biphyllum (both 2n = 96 chromosomes). The rest varied in GS (2C = 3.54-8.78 pg) with a constant chromosome number 2n = 48. There is a weak connection between phylogeny, GS and geography in Amomum. CONCLUSIONS: Amomum consists of four groups, and the former Elettariopsis is monophyletic. Species in this group have the largest GS. Two polyploids were found and GS greatly varied in the rest of Amomum.

Phylogenomic analyses using a new 1013-gene Vitaceae bait-set support major groups of North American Vitis.

A set of newly designed Vitaceae baits targeting 1013 genes was employed to explore phylogenetic relationships among North American Vitis. Eurasian Vitis taxa including Vitis vinifera were found to be nested within North American Vitis subgenus Vitis. North American Vitis subgenus Vitis can be placed into nine main groups: the Monticola group, the Occidentales group, the Californica group, the Vinifera group (introduced from Eurasia), the Mustangensis group, the Palmata group, the Aestivalis group, the Labrusca group, and the Cinerea group. Strong cytonuclear discordances were detected in North American Vitis, with many species non-monophyletic in the plastid phylogeny, while monophyletic in the nuclear phylogeny. The phylogenomic analyses support recognizing four distinct species in the Vitis cinerea complex in North America: V. cinerea, V. baileyana, V. berlandieri, and V. simpsonii. Such treatment will better serve the conservation of wild Vitis diversity in North America. Yet the evolutionary history of Vitis is highly complex, with the concordance analyses indicating conflicting signals across the phylogeny. Cytonuclear discordances and Analyses using the Species Networks applying Quartets (SNaQ) method support extensive hybridizations in North American Vitis. The results further indicate that plastid genomes alone are insufficient for resolving the evolutionary history of plant groups that have undergone rampant hybridization, like the case in North American Vitis. Nuclear gene data are essential for species delimitation, identification and reconstructing evolutionary relationships; therefore, they are imperative for plant phylogenomic studies.

Impact of the climatic changes in the Pliocene-Pleistocene transition on Irano-Turanian species. The radiation of genus Jurinea (Compositae).

The Irano-Turanian region is one of the world's richest floristic regions and the centre of diversity for numerous xerophytic plant lineages. However, we still have limited knowledge on the timing of evolution and biogeographic history of its flora, and potential drivers of diversification remain underexplored. To fill this knowledge gap, we focus on the Eurasian genus Jurinea (ca. 200 species), one of the largest plant radiations that diversified in the region. We applied a macroevolutionary integrative approach to explicitly test diversification hypotheses and investigate the relative roles of geography vs. ecology and niche conservatism vs. niche lability in speciation processes. To do so, we gathered a sample comprising 77% of total genus richness and obtained data about (1) its phylogenetic history, recovering 502 nuclear loci sequences; (2) growth forms; (3) ecological niche, compiling data of 21 variables for more than 2500 occurrences; and (4) paleoclimatic conditions, to estimate climatic stability. Our results revealed that climate was a key factor in the evolutionary dynamics of Jurinea. The main diversification and biogeographic events that occurred during past climate changes, which led to colder and drier conditions, are the following: (1) the origin of the genus (10.7 Ma); (2) long-distance dispersals from the Iranian Plateau to adjacent regions (∼7-4 Ma); and (3) the diversification shift during Pliocene-Pleistocene Transition (ca. 3 Ma), when net diversification rate almost doubled. Our results supported the pre-adaptation hypothesis, i.e., the evolutionary success of Jurinea was linked to the retention of the ancestral niche adapted to aridity. Interestingly, the paleoclimatic analyses revealed that in the Iranian Plateau long-term climatic stability favoured old-lineage persistence, resulting in current high species richness of semi-arid and cold adapted clades; whereas moderate climate oscillations stimulated allopatric diversification in the lineages distributed in the Circumboreal region. In contrast, growth form lability and high niche disparity among closely related species in the Central Asian clade suggest adaptive radiation to mountain habitats. In sum, the radiation of Jurinea is the result of both adaptive and non-adaptive processes influenced by climatic, orogenic and ecological factors.

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.

Genomic, spatial and morphometric data for discrimination of four species in the Mediterranean Tamus clade of yams (Dioscorea, Dioscoreaceae).

BACKGROUND AND AIMS: Among the numerous pantropical species of the yam genus, Dioscorea, only a small group occurs in the Mediterranean basin, including two narrow Pyrenean endemics (Borderea clade) and two Mediterranean-wide species (D. communis and D. orientalis, Tamus clade). However, several currently unrecognized species and infraspecific taxa have been described in the Tamus clade due to significant morphological variation associated with D. communis. Our overarching aim was to investigate taxon delimitation in the Tamus clade using an integrative approach combining phylogenomic, spatial and morphological data. METHODS: We analysed 76 herbarium samples using Hyb-Seq genomic capture to sequence 260 low-copy nuclear genes and plastomes, together with morphometric and environmental modelling approaches. KEY RESULTS: Phylogenomic reconstructions confirmed that the two previously accepted species of the Tamus clade, D. communis and D. orientalis, are monophyletic and form sister clades. Three subclades showing distinctive geographic patterns were identified within D. communis. These subclades were also identifiable from morphometric and climatic data, and introgression patterns were inferred between subclades in the eastern part of the distribution of D. communis. CONCLUSIONS: We propose a taxonomy that maintains D. orientalis, endemic to the eastern Mediterranean region, and splits D. communis sensu lato into three species: D. edulis, endemic to Macaronesia (Canary Islands and Madeira); D. cretica, endemic to the eastern Mediterranean region; and D. communis sensu stricto, widespread across western and central Europe. Introgression inferred between D. communis s.s. and D. cretica is likely to be explained by their relatively recent speciation at the end of the Miocene, disjunct isolation in eastern and western Mediterranean glacial refugia and a subsequent westward recolonization of D. communis s.s. Our study shows that the use of integrated genomic, spatial and morphological approaches allows a more robust definition of species boundaries and the identification of species that previous systematic studies failed to uncover.

A Hyb-Seq phylogeny of Boechera and related genera using a combination of Angiosperms353 and Brassicaceae-specific bait sets.

PREMISE: Although Boechera (Boechereae, Brassicaceae) has become a plant model system for both ecological genomics and evolutionary biology, all previous phylogenetic studies have had limited success in resolving species relationships within the genus. The recent effective application of sequence data from target enrichment approaches to resolve the evolutionary relationships of several other challenging plant groups prompted us to investigate their usefulness in Boechera and Boechereae. METHODS: To resolve the phylogeny of Boechera and closely related genera, we utilized the Hybpiper pipeline to analyze two combined bait sets: Angiosperms353, with broad applicability across flowering plants; and a Brassicaceae-specific bait set designed for use in the mustard family. Relationships for 101 samples representing 81 currently recognized species were inferred from a total of 1114 low-copy nuclear genes using both supermatrix and species coalescence methods. RESULTS: Our analyses resulted in a well-resolved and highly supported phylogeny of the tribe Boechereae. Boechereae is divided into two major clades, one comprising all western North American species of Boechera, the other encompassing the eight other genera of the tribe. Our understanding of relationships within Boechera is enhanced by the recognition of three core clades that are further subdivided into robust regional species complexes. CONCLUSIONS: This study presents the first broadly sampled, well-resolved phylogeny for most known sexual diploid Boechera. This effort provides the foundation for a new phylogenetically informed taxonomy of Boechera that is crucial for its continued use as a model system.

An updated phylogeny, biogeography, and PhyloCode-based classification of Cornaceae based on three sets of genomic data.

PREMISE: A major goal of systematic biology is to uncover the evolutionary history of organisms and translate that knowledge into stable classification systems. Here, we integrate three sets of genome-wide data to resolve phylogenetic relationships in Cornaceae (containing only Cornus s.l.), reconstruct the biogeographic history of the clade, and provide a revised classification using the PhyloCode to stabilize names for this taxonomically controversial group. METHODS: We conducted phylogenetic analyses using 312 single-copy nuclear genes and 70 plastid genes from Angiosperms353 Hyb-Seq, plus numerous loci from RAD-Seq. We integrated fossils using morphological data and produced a dated phylogeny for biogeographical analysis. RESULTS: A well-resolved, strongly supported, comprehensive phylogeny was obtained. Biogeographic analyses support an origin and rapid diversification of Cornus into four morphologically distinct major clades in the Northern Hemisphere (with an eastern Asian ancestor) during the late Cretaceous. Dispersal into Africa from eastern Asia likely occurred along the Tethys Seaway during the Paleogene, whereas dispersal into South America likely occurred during the Neogene. Diversification within the northern hemisphere likely involved repeated independent colonization of new areas during the Paleogene and Neogene along the Bering Land Bridge, the North Atlantic Land Bridge, and the Tethys Seaway. Thirteen strongly supported clades were named following rules of the PhyloCode. CONCLUSIONS: Our study provides an example of integrating genomic and morphological data to produce a robust, explicit species phylogeny that includes fossil taxa, which we translate into an updated classification scheme using the PhyloCode to stabilize names.

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.

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex.

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

Phylogenomics AND biogeography of Castanea (chestnut) and Hamamelis (witch-hazel) - Choosing between RAD-seq and Hyb-Seq approaches.

Hyb-Seq and RAD-seq are well-established high throughput sequencing technologies that have been increasingly used for plant phylogenomic studies. Each method has its own pros and cons. The choice between them is a practical issue for plant systematists studying the evolutionary histories of biodiversity of relatively recent origins. However, few studies have compared the congruence and conflict between results from the two methods within the same group of organisms in plants. In this study, we employed RAD-seq and Hyb-Seq of Angiosperms353 genes in phylogenomic and biogeographic studies of Hamamelis (the witch-hazels) and Castanea (chestnuts), two classic examples exhibiting the well-known eastern Asian (EA) -eastern North American (ENA) disjunct distribution, and compared them side by side. Our results showed congruences in phylogenetic inference and divergence time dating between the two data sets obtained through our customized procedures of library preparation and sequence trimming, although they differed in the number of loci and informative sites, the amount of missing data, and sampling within species. We provide recommendations regarding the selection of the two methods for phylogenomic study at generic level based on fund availability and sampling scale. If funds and time are not constrained, we recommend Hyb-Seq. If funds and time are somewhat limited and sampling is large, we recommend RAD-seq. However, we found greater conflict among gene trees from the RAD-seq data due to the short sequences per locus. Therefore, species tree building and network detecting with the RAD-seq data with short RAD-seq loci (e.g., <150 bp) should avoid using analytical methods relying on gene trees of individual locus, but using site-based methods such as SVDQuartets and D-statistic method. Our phylogenetic analyses of RAD-seq and Hyb-Seq data resulted in well-resolved species relationships. Analyses of the data using the D-statistic test and PhyloNet revealed ancient introgressions in both genera. Biogeographic analyses including fossil data using total evidence-based dated tree and DEC model applying specific inter-area dispersal probabilities revealed a complicated history for each genus, indicating multiple interareal dispersals and local extinctions within and outside areas of the taxa's modern ranges in both the Paleogene and Neogene. The study demonstrates the importance of including fossil taxa for a more accurate reconstruction of biogeographic histories of taxa to understand the EA and ENA floristic disjunction. Our results support a widespread ancestral range in EA-western North America (WNA) followed by early diversification in EA and expansion to North America (NA) and Europe for Castanea and a more widespread ancestral range in EA-ENA-WNA for Hamamelis. The origins of the modern EA-ENA disjunction in both genera were suggested to be the result of vicariance from widespread ancestors in Eurasia-ENA of the mid-Miocene and in EA-NA of the late Oligocene, respectively.

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.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers.

PREMISE: The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods. METHODS: We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs). RESULTS: The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae). CONCLUSIONS: Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Coalescent-based species delimitation in North American pinyon pines using low-copy nuclear genes and plastomes.

PREMISE: Accurate species delimitation is essential for evolutionary biology, conservation, and biodiversity management. We studied species delimitation in North American pinyon pines, Pinus subsection Cembroides, a natural group with high levels of incomplete lineage sorting. METHODS: We used coalescent-based methods and multivariate analyses of low-copy number nuclear genes and nearly complete high-copy number plastomes generated with the Hyb-Seq method. The three coalescent-based species delimitation methods evaluated were the Generalized Mixed Yule Coalescent (GMYC), Poisson Tree Process (PTP), and Trinomial Distribution of Triplets (Tr2). We also measured admixture in populations with possible introgression. RESULTS: Our results show inconsistencies among GMYC, PTP, and Tr2. The single-locus based GMYC analysis of plastid DNA recovered a higher number of species (up to 24 entities, including singleton lineages and clusters) than PTP and the multi-locus coalescent approach. The PTP analysis identified 10 species whereas Tr2 recovered 13, which agreed closely with taxonomic treatments. CONCLUSIONS: We found that PTP and GMYC identified species with low levels of ILS and high morphological divergence (P. maximartinezii, P. pinceana, and P. rzedowskii). However, GMYC method oversplit species by identification of more divergent samples as singletons. Moreover, both PTP and GMYC were incapable of identifying some species that are readily identified morphologically. We suggest that the divergence times between lineages within North American pinyon pines are so disparate that GMYC results are unreliable. Results of the Tr2 method coincided well with previous delimitations based on morphology, DNA, geography, and secondary chemistry.

Phylogeny, origin, and dispersal of Dubyaea (Asteraceae) based on Hyb-Seq data.

Dubyaea DC. is a small genus of Asteraceae that is almost exclusively endemic to the Pan-Himalayan region. Within Dubyaea, phylogenetic relationships remain poorly understood. Here, our well-supported phylogeny based on Hyb-Seq data shows that all samples of Dubyaea in this study belong to a monophyletic group, which is sister to the clade of Soroseris, Syncalathium, and Nabalus. Dubyaea (s. str.) can be divided into three major clades, which are supported by flower color as well as morphological features of the stems and basal leaves. Based on our phylogenetic results, we performed biogeographic analyses and inferred that Dubyaea arose in the late Miocene in Hengduan Mountains and its eastern areas. Following its evolutionary origin, Dubyaea underwent diversification in situ as well as spread to the Himalayas.

Phylogenomic relationships and character evolution of the grape family (Vitaceae).

The grape family consists of 16 genera and ca. 950 species. It is best known for the economically important fruit crop - the grape Vitis vinifera. The deep phylogenetic relationships and character evolution of the grape family have attracted the attention of researchers in recent years. We herein reconstruct the phylogenomic relationships within Vitaceae using nuclear and plastid genes based on the Hyb-Seq approach and test the newly proposed classification system of the family. The five tribes of the grape family, including Ampelopsideae, Cayratieae, Cisseae, Parthenocisseae, and Viteae, are each robustly supported by both nuclear and chloroplast genomic data and the backbone relationships are congruent with previous reports. The cupular floral disc (raised above and free from ovary at the upper part) is an ancestral state of Vitaceae, with the inconspicuous floral disc as derived in the tribe Parthenocisseae, and the state of adnate to the ovary as derived in the tribe Viteae. The 5-merous floral pattern was inferred to be the ancestral in Vitaceae, with the 4-merous flowers evolved at least two times in the family. The compound dichasial cyme (cymose with two secondary axes) is ancestral in Vitaceae and the thyrse inflorescence (a combination of racemose and cymose branching) in tribe Viteae is derived. The ribbon-like trichome only evolved once in Vitaceae, as a synapomorphy for the tribe Viteae.

Settling a family feud: a high-level phylogenomic framework for the Gentianales based on 353 nuclear genes and partial plastomes.

PREMISE: Comprising five families that vastly differ in species richness-ranging from Gelsemiaceae with 13 species to the Rubiaceae with 13,775 species-members of the Gentianales are often among the most species-rich and abundant plants in tropical forests. Despite considerable phylogenetic work within particular families and genera, several alternative topologies for family-level relationships within Gentianales have been presented in previous studies. METHODS: Here we present a phylogenomic analysis based on nuclear genes targeted by the Angiosperms353 probe set for approximately 150 species, representing all families and approximately 85% of the formally recognized tribes. We were able to retrieve partial plastomes from off-target reads for most taxa and infer phylogenetic trees for comparison with the nuclear-derived trees. RESULTS: We recovered high support for over 80% of all nodes. The plastid and nuclear data are largely in agreement, except for some weakly to moderately supported relationships. We discuss the implications of our results for the order's classification, highlighting points of increased support for previously uncertain relationships. Rubiaceae is sister to a clade comprising (Gentianaceae + Gelsemiaceae) + (Apocynaceae + Loganiaceae). CONCLUSIONS: The higher-level phylogenetic relationships within Gentianales are confidently resolved. In contrast to recent studies, our results support the division of Rubiaceae into two subfamilies: Cinchonoideae and Rubioideae. We do not formally recognize Coptosapelteae and Luculieae within any particular subfamily but treat them as incertae sedis. Our framework paves the way for further work on the phylogenetics, biogeography, morphological evolution, and macroecology of this important group of flowering plants.