An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus)

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

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

Biased Gene Introgression and Adaptation in the Face of Chloroplast Capture in Aquilegia amurensis

-Chloroplast capture, a phenomenon that can occur through interspecific hybridization and introgression, is frequently invoked to explain cytonuclear discordance in plants. However, relatively few studies have documented the mechanisms of cytonuclear coevolution and its potential for driving species differentiation and possible functional differences in the context of chloroplast capture. To address this crucial question, we chose the Aquilegia genus, which is known for having minimal sterility among species, and inferred that A. amurensis captured the plastome of A. parviflora based on cytonuclear discordance and gene flow between the 2 species. We focused on the introgression region and its differentiation from corresponding regions in closely related species, especially its composition in a chloroplast capture scenario. We found that nuclear genes encoding cytonuclear enzyme complexes (CECs; i.e., organelle-targeted genes) of chloroplast donor species were selectively retained and displaced the original CEC genes in chloroplast-receiving species due to cytonuclear interactions during introgression. Notably, the intrinsic correlation of CEC introgression was a greater degree of evolutionary distance for these CECs between A. amurensis and A. parviflora. Terpene synthase activity genes (GO: 0010333) were overrepresented among the introgressed genes, and more than 30% of these genes were CEC genes. These findings support our observations that floral terpene release pattern is similar between A. amurensis and A. parviflora compared with A. japonica. Our study clarifies the mechanisms of cytonuclear coevolution, species differentiation, and functional differences in the context of chloroplast capture and highlights the potential role of chloroplast capture in adaptation.

Demography-driven and adaptive introgression in a hybrid zone of the Armeria syngameon

Syngameons represent networks of otherwise distinct species connected by limited gene exchange. Although most studies have focused on how species maintain their cohesiveness despite gene flow, there are additional relevant questions regarding the evolutionary dynamics of syngameons and their drivers, as well as the success of their members and the network as a whole. Using a ddRADseq approach, we analysed the genetic structure, genomic clines and demographic history of a coastal hybrid zone involving two species of the Armeria (Plumbaginaceae) syngameon in southern Spain. We inferred that a peripheral population of the sand dune-adapted A. pungens diverged from the rest of the conspecific populations and subsequently hybridized with a locally more abundant pinewood congener, A. macrophylla. Both species display extensive plastid DNA haplotype sharing. Genomic cline analysis identified bidirectional introgression, but more outlier loci with excess A. pungens than A. macrophylla ancestry, suggesting the possibility of selection for A. pungens alleles. This is consistent with the finding that the A. pungens phenotype is selected for in open habitats, and with the strong correlation found between ancestry and phenotype. Taken together, our analyses suggest an intriguing scenario in which bidirectional introgression may, on the one hand, help to avoid reduced levels of genetic diversity due to the small size and isolated location of the A. pungens range-edge population, thereby minimizing demographic risks of stochastic extinction. On the other hand, the data also suggest that introgression into A. macrophylla may allow individuals to grow in open, highly irradiated, deep sandy, salt-exposed habitats.

Phylogenetic incongruence in an Asiatic species complex of the genus Caryodaphnopsis (Lauraceae)

BACKGROUND

Caryodaphnopsis, a group of tropical trees (ca. 20 spp.) in the family Lauraceae, has an amphi-Pacific disjunct distribution: ten species are distributed in Southeast Asia, while eight species are restricted to tropical rainforests in South America. Previously, phylogenetic analyses using two nuclear markers resolved the relationships among the five species from Latin America. However, the phylogenetic relationships between the species in Asia remain poorly known.

RESULTS

Here, we first determined the complete mitochondrial genome (mitogenome), plastome, and the nuclear ribosomal cistron (nrDNA) sequences of C. henryi with lengths of 1,168,029 bp, 154,938 bp, and 6495 bp, respectively. We found 2233 repeats and 368 potential SSRs in the mitogenome of C. henryi and 50 homologous DNA fragments between its mitogenome and plastome. Gene synteny analysis revealed a mass of rearrangements in the mitogenomes of Magnolia biondii, Hernandia nymphaeifolia, and C. henryi and only six conserved clustered genes among them. In order to reconstruct relationships for the ten Caryodaphnopsis species in Asia, we created three datasets: one for the mitogenome (coding genes and ten intergenic regions), another for the plastome (whole genome), and the other for the nuclear ribosomal cistron. All of the 22 Caryodaphnopsis individuals were divided into four, five, and six different clades in the phylogenies based on mitogenome, plastome, and nrDNA datasets, respectively.

CONCLUSIONS

The study showed phylogenetic conflicts within and between nuclear and organellar genome data of Caryodaphnopsis species. The sympatric Caryodaphnopsis species in Hekou and Malipo SW China may be related to the incomplete lineage sorting, chloroplast capture, and/or hybridization, which mixed the species as a complex in their evolutionary history.

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

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.

Historical climate change and vicariance events contributed to the intercontinental disjunct distribution pattern of ash species (Fraxinus, Oleaceae)

The Northern Hemisphere temperate forests exhibit a disjunct distributional pattern in Europe, North America, and East Asia. Here, to reveal the promoter of intercontinental disjunct distribution, Fraxinus was used as a model organism to integrate abundant fossil evidence with high-resolution phylogenies in a phytogeographic analysis. We constructed a robust phylogenetic tree using genomic data, reconstructed the geographic ancestral areas, and evaluated the effect of incorporating fossil information on the reconstructed biogeographic history. The phylogenetic relationships of Fraxinus were highly resolved and divided into seven clades. Fraxinus originated in western North America during Eocene, and six intercontinental dispersal events and five intercontinental vicariance events were occured. Results suggest that climate change and vicariance contributed to the intercontinental disjunct distribution pattern of Fraxinus. Moreover, results highlight the necessity of integrating phylogenetic relationship and fossil to improve the reliability of inferred biogeographic events and our understanding of the processes underlying disjunct distributions.

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.

Evolutionary history of an Irano-Turanian cushion-forming legume (Onobrychis cornuta)

The Irano-Turanian region is one of the largest floristic regions in the world and harbors a high percentage of endemics, including cushion-like and dwarf-shrubby taxa. Onobrychis cornuta is an important cushion-forming element of the subalpine/alpine flora of the Irano-Turanian floristic region. To specify the genetic diversity among the populations of this species (including individuals of O. elymaitica), we employed nrDNA ITS and two noncoding regions of plastid DNA (rpl32-trnL(UAG) and trnT(UGU)-trnL(UAA)). The most striking feature of O. cornuta assemblages was the unexpectedly high nucleotide diversity in both the nDNA and cpDNA dataset. In the analyses of nuclear and plastid regions, 25 ribotypes and 42 haplotypes were found among 77 and 59 accessions, respectively, from Iran, Turkey, and Afghanistan. Network analysis of the datasets demonstrated geographic differentiation within the species. Phylogenetic analyses of all dataset retrieved O. cornuta as a non-monophyletic species due to the inclusion of O. elymaitica, comprising four distinct lineages. In addition, our analyses showed cytonuclear discordance between both nuclear and plastid topologies regarding the position of some O. cornuta individuals. The underlying causes of this inconsistency remain unclear. However, we speculate that chloroplast capture, incomplete lineage sorting, and introgression were the main reasons for this event. Furthermore, molecular dating analysis indicated that O. cornuta originated in the early Pliocene (around 4.8 Mya) and started to diversify throughout the Pliocene and in particular the Pleistocene. Moreover, O. elymaitica was reduced to a subspecific rank within the species.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

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

Development of molecular identification methods for Dryophytes suweonensis and D. japonicus, and their hybrids

Background

As hybridization can reduce biodiversity or cause extinction, it is important to identify both purebred parental species and their hybrids prior to conserving them. The Suwon tree frog, Dryophytes suweonensis, is an endangered wildlife species in Korea that shares its habitat and often hybridizes with the Japanese tree frog, D. japonicus. In particular, D. suweonensis, D. japonicus, and their hybrids often have abnormal ovaries and gonads, which are known causes that could threaten their existence.

Methods

We collected 57 individuals from six localities where D. suweonensis is known to be present. High-resolution melting curve (HRM) analysis of the mitochondrial 12S ribosomal RNA gene was performed to determine the maternal species. Thereafter, the DNA sequences of five nuclear genes (SIAH, TYR, POMC, RAG1, and C-MYC) were analyzed to determine their parental species and hybrid status.

Results

The HRM analysis showed that the melting temperature of D. suweonensis was in the range of 79.0-79.3 °C, and that of D. japonicus was 77.7-78.0 °C, which clearly distinguished the two tree frog species. DNA sequencing of the five nuclear genes revealed 37 single-nucleotide polymorphism (SNP) sites, and STRUCTURE analysis showed a two-group structure as the most likely grouping solution. No heterozygous position in the purebred parental sequences with Q values ≥ 0.995 were found, which clearly distinguished the two treefrog species from their hybrids; 11 individuals were found to be D. suweonensis, eight were found to be D. japonicus, and the remaining 38 individuals were found to be hybrids.

Conclusion

Thus, it was possible to unambiguously identify the parental species and their hybrids using HRM analysis and DNA sequencing methods. This study provided fundamental information for D. suweonensis conservation and restoration research.

Assessing the phylogenetic relationship among varieties of Toona ciliata (Meliaceae) in sympatry with chloroplast genomes

Toona ciliata is an endangered species due to over-cutting and low natural regeneration in China. Its genetic conservation is of an increasing concern. However, several varieties are recognized according to the leaf and flower traits, which complicates genetic conservation of T. ciliata. Here, we sequenced the whole chloroplast genome sequences of three samples for each of four varieties (T. ciliata var. ciliata, T. ciliata var. yunnanensis, T. ciliata var. pubescens, and T. ciliata var. henryi) in sympatry and assessed their phylogenetic relationship at a fine spatial scale. The four varieties had genome sizes ranged from 159,546 to 159,617 bp and had small variations in genome structure. Phylogenomic analysis indicated that the four varieties were genetically well-mixed in branch groups. Genetic diversity from the whole chloroplast genome sequences of 12 samples was low among varieties (average π = 0.0003). Besides, we investigated genetic variation of 58 samples of the four varieties in sympatry using two markers (psaA and trnL-trnF) and showed that genetic differentiation was generally insignificant among varieties (Ф st = 0%-5%). Purifying selection occurred in all protein-coding genes except for the ycf2 gene that was under weak positive selection. Most amino acid sites in all protein-coding genes were under purifying selection except for a few sites that were under positive selection. The chloroplast genome-based phylogeny did not support the morphology-based classification. The overall results implicated that a conservation strategy based on the T. ciliata complex rather than on intraspecific taxon was more appropriate.

Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains

Mountainous regions provide a multitude of habitats and opportunities for complex speciation scenarios. Hybridization leading to chloroplast capture, which can be revealed by incongruent phylogenetic trees, is one possible outcome. Four allopatric Taxus lineages (three species and an undescribed lineage) from the Hengduan Mountains, southwest China, exhibit conflicting phylogenetic relationships between nuclear and chloroplast phylogenies. Here, we use multi-omic data at the population level to investigate their historical speciation processes. Population genomic analysis based on ddRAD-seq data revealed limited contemporary inter-specific gene flow involving only populations located close to another species. In a historical context, chloroplast and nuclear data (transcriptome) consistently showed conflicting phylogenetic relationships for T. florinii and the Emei type lineage. ILS and chloroplast recombination were excluded as possible causes, and transcriptome and ddRAD-seq data revealed an absence of the mosaic nuclear genomes that characterize hybrid origin scenarios. Therefore, T. florinii appears to have originated when a lineage of T. florinii captured the T. chinensis plastid type, whereas plastid introgression in the opposite direction generated the Emei Type. All four species have distinct ecological niche based on community investigations and ecological niche analyses. We propose that the origins of both species represent very rare examples of chloroplast capture events despite the paternal cpDNA inheritance of gymnosperms. Specifically, allopatrically and/or ecologically diverged parental species experienced a rare secondary contact, subsequent hybridization and reciprocal chloroplast capture, generating two new lineages, each of which acquired a unique ecological niche. These events might have been triggered by orogenic activities of the Hengduan Mountains and an intensification of the Asian monsoon in the late Miocene, and may represent a scenario more common in these mountains than presently known.

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

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

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

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

Ixeridiumsagittarioides (Asteraceae-Cichorieae) revisited: range extension and molecular evidence for its systematic position in the Lactuca alliance

Our first record of the rare and scatteredly distributed Ixeridiumsagittarioides for Guizhou, China, triggered a study to assess its systematic position. The species was placed in four different genera in the course of its taxonomic history and was recently treated with doubts as a member of Ixeridium in the Flora of China. Comparative morphological investigation and phylogenetic analyses based on the nuclear ribosomal DNA internal transcribed spacer (nrITS) and five non-coding plastid DNA regions (petD region, psbA-trnH, trnL-trnF, rpl32-trnL (UAG) and 5´rps16-trnQ (UUG) spacers) provided evidence that the species is not a member of Ixeridium and the Crepidinae but has evolved by ancient hybridisation of members of the Lactuca alliance (Lactucinae). It is reinstated as Lactucasagittarioides and a comprehensive morphological description is provided, based on material from its entire range of distribution.

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.

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

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

A phylogenomic analysis of Lonicera and its bearing on the evolution of organ fusion

PREMISE

The ~140 species of Lonicera are characterized by variously fused leaves, bracteoles, and ovaries, making it a model system for studying the evolution and development of organ fusion. However, previous phylogenetic analyses, based mainly on chloroplast DNA markers, have yielded uncertain and conflicting results. A well-supported phylogeny of Lonicera will allow us to trace the evolutionary history of organ fusion.

METHODS

We inferred the phylogeny of Lonicera using restriction site-associated DNA sequencing (RADSeq), sampling all major clades and 18 of the 23 subsections. This provided the basis for inferring the evolution of five fusion-related traits.

RESULTS

RADSeq data yielded a well-resolved and well-supported phylogeny. The two traditionally recognized subgenera (Periclymenum and Chamaecerasus), three of the four sections (Isoxylosteum, Coeloxylosteum, and Nintooa), and half of the subsections sampled were recovered as monophyletic. However, the large and heterogeneous section Isika was strongly supported as paraphyletic. Nintooa, a clade of ~22 mostly vine-forming species, including L. japonica, was recovered in a novel position, raising the possibility of cytonuclear discordance. We document the parallel evolution of fused leaves, bracteoles, and ovaries, with rare reversals. Most strikingly, complete cupules, in which four fused bracteoles completely enclose two unfused ovaries, arose at least three times. Surprisingly, these appear to have evolved directly from ancestors with free bracteoles instead of partial cupules.

CONCLUSIONS

We provide the most comprehensive and well-supported phylogeny of Lonicera to date. Our inference of multiple evolutionary shifts in organ fusion provides a solid foundation for in depth developmental and functional analyses.

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.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species

BACKGROUND AND AIMS

Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes.

METHODS

We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression.

KEY RESULTS

Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases.

CONCLUSIONS

The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

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.

Re-evaluating Homoploid Reticulate Evolution in Helianthus Sunflowers

Sunflowers of the genus Helianthus are models for hybridization research and contain three of the best-studied examples of homoploid hybrid speciation. To understand a broader picture of hybridization within the annual sunflowers, we used whole-genome resequencing to conduct a phylogenomic analysis and test for gene flow between lineages. We find that all annual sunflower species tested have evidence of admixture, suggesting hybridization was common during the radiation of the genus. Support for the major species tree decreases with increasing recombination rate, consistent with hybridization and introgression contributing to discordant topologies. Admixture graphs found hybridization to be associated with the origins of the three putative hybrid species (Helianthus anomalus, Helianthus deserticola, and Helianthus paradoxus). However, the hybridization events are more ancient than suggested by previous work. Furthermore, H. anomalus and H. deserticola appear to have arisen from a single hybridization event involving an unexpected donor, rather than through multiple independent events as previously proposed. This means our results are consistent with, but not definitive proof of, two ancient independent homoploid hybrid speciation events in the genus. Using a broader data set that covers the whole Helianthus genus, including perennial species, we find that signals of introgression span the genus and beyond, suggesting highly divergent introgression and/or the sorting of ancient haplotypes. Thus, Helianthus can be viewed as a syngameon in which largely reproductively isolated species are linked together by occasional or frequent gene flow.

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots.

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.

Genome-partitioning strategy, plastid and nuclear phylogenomic discordance, and its evolutionary implications of Clematis (Ranunculaceae)

Clematis is one of the largest genera of Ranunculaceae with many phylogenetic problems left to be resolved. Clematis species have considerable genome size of more than 7 Gbp, and there was no whole-genome reference sequence published in this genus. This raises difficulties in acquiring nuclear genome data for its phylogenetic analysis. Previous studies based on Sanger sequencing data, plastid genome data, and nrDNA sequences did not well resolve the phylogeny of Clematis. In this study, we used genome skimming and transcriptome data to assemble the plastid genome sequences, nuclear single nucleotide polymorphisms (SNPs) datasets, and single-copy nuclear orthologous genes (SCOGs) to reconstruct the phylogenetic backbone of Clematis, and test effectiveness of these genome partitioning methods. We also further analyzed the discordance among nuclear gene trees and between plastid and nuclear phylogenies. The results showed that the SCOGs datasets, assembled from transcriptome method, well resolved the phylogenetic backbone of Clematis. The nuclear SNPs datasets from genome skimming method can also produce similar results with the SCOGs data. In contrast to the plastid phylogeny, the phylogeny resolved by nuclear genome data is more robust and better corresponds to morphological characters. Our results suggested that rapid species radiation may have generated high level of incomplete lineage sorting, which was the major cause of nuclear gene discordance. Our simulation also showed that there may have been frequent interspecific hybridization events, which led to some of the cyto-nuclear discordances in Clematis. This study not only provides the first robust phylogenetic backbone of Clematis based on nuclear genome data, but also provides suggestions of genome partitioning strategies for the phylogenomic study of other plant taxa.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Plastome sequences fail to resolve shallow level relationships within the rapidly radiated genus Isodon (Lamiaceae)

As one of the largest genera of Lamiaceae and of great medicinal importance, Isodon is also phylogenetically and taxonomically recalcitrant largely ascribed to its recent rapid radiation in the Hengduan Mountains. Previous molecular phylogenetic studies using limited loci have only successfully resolved the backbone topology of the genus, but the interspecific relationships suffered from low resolution, especially within the largest clade (Clade IV) which comprises over 80% species. In this study, we attempted to further elucidate the phylogenetic relationships within Isodon especially Clade IV using plastome sequences with a broad taxon sampling of ca. 80% species of the genus. To reduce systematic errors, twelve different plastome data sets (coding and non-coding regions with ambiguously aligned regions and saturated loci removed or not) were employed to reconstruct phylogeny using maximum likelihood and Bayesian inference. Our results revealed largely congruent topologies of the 12 data sets and recovered major lineages of Isodon consistent with previous studies, but several incongruences are also found among these data sets and among single plastid loci. Most of the shallow nodes within Clade IV were resolved with high support but extremely short branch lengths in plastid trees, and showed tremendous conflicts with the nrDNA tree, morphology and geographic distribution. These incongruences may largely result from stochasticity (due to insufficient phylogenetic signal) and hybridization and plastid capture. Therefore, the uniparental-inherited plastome sequences are insufficient to disentangle relationships within a genus which has undergone recent rapid diversification. Our findings highlight a need for additional data from nuclear genome to resolve the relationships within Clade IV and more focused studies to assess the influences of multiple processes in the evolutionary history of Isodon. Nevertheless, the morphology of the shape and surface sculpture/indumentum of nutlets is of systematic importance that they can distinguish the four major clades of Isodon.

Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L

Background

Ferula L. is one of the largest and most taxonomically complicated genera as well as being an important medicinal plant resource in the family Apiaceae. To investigate the plastome features and phylogenetic relationships of Ferula and its neighboring genera Soranthus Ledeb., Schumannia Kuntze., and Talassia Korovin, we sequenced 14 complete plastomes of 12 species. 

Results

The size of the 14 complete chloroplast genomes ranged from 165,607 to 167,013 base pairs (bp) encoding 132 distinct genes (87 protein-coding, 37 tRNA, and 8 rRNA genes), and showed a typical quadripartite structure with a pair of inverted repeats (IR) regions. Based on comparative analysis, we found that the 14 plastomes were similar in codon usage, repeat sequence, simple sequence repeats (SSRs), and IR borders, and had significant collinearity. Based on our phylogenetic analyses, Soranthus, Schumannia, and Talassia should be considered synonymous with Ferula. Six highly divergent regions (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, and ycf1) were also detected, which may represent potential molecular markers, and combined with selective pressure analysis, the weak positive selection gene ccsA may be a discriminating DNA barcode for Ferula species.

Conclusion

Plastids contain abundant informative sites for resolving phylogenetic relationships. Combined with previous studies, we suggest that there is still much room for improvement in the classification of Ferula. Overall, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of this genus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08868-z.

Nonadaptive molecular evolution of plastome during the speciation of Actaea purpurea and its relatives

We have seen an explosive increase of plant plastid genome (plastome) sequences in the last decade, and the view that sequence variation in plastomes is maintained by the mutation-drift balance has been challenged by new evidence. Although comparative genomic and population-level studies provided us with evidence for positive evolution of plastid genes at both the macro- and micro-evolution levels, less studies have systematically investigated how plastomes have evolved during the speciation process. We here sequenced 13 plastomes of Actaea purpurea (P.K. Hsiao) J. Compton, and its closest relatives, and conducted a systematic survey of positive selection in their plastid genes using the McDonald-Kreitman test and codon-based methods using maximum likelihood to estimate the ratio of nonsynonymous to synonymous substitutions (ω) across a phylogeny. We found that during the speciation of A. purpurea and its relatives, all plastid genes evolved neutrally or were under purifying selection. Genome size, gene order, and number were highly conserved. Comparing to A. purpurea, plastomes of Actaea japonica and Actaea biternata had low genetic diversity, consistent with previous studies. Our work not only sheds important light on the evolutionary history of A. purpurea and its kin, but also on the evolution of plastomes during plant speciation.

Phylogenomic analysis of 997 nuclear genes reveals the need for extensive generic re-delimitation in Caesalpinioideae (Leguminosae)

Subfamily Caesalpinioideae with ca. 4,600 species in 152 genera is the second-largest subfamily of legumes (Leguminosae) and forms an ecologically and economically important group of trees, shrubs and lianas with a pantropical distribution. Despite major advances in the last few decades towards aligning genera with clades across Caesalpinioideae, generic delimitation remains in a state of considerable flux, especially across the mimosoid clade. We test the monophyly of genera across Caesalpinioideae via phylogenomic analysis of 997 nuclear genes sequenced via targeted enrichment (Hybseq) for 420 species and 147 of the 152 genera currently recognised in the subfamily. We show that 22 genera are non-monophyletic or nested in other genera and that non-monophyly is concentrated in the mimosoid clade where ca. 25% of the 90 genera are found to be non-monophyletic. We suggest two main reasons for this pervasive generic non-monophyly: (i) extensive morphological homoplasy that we document here for a handful of important traits and, particularly, the repeated evolution of distinctive fruit types that were historically emphasised in delimiting genera and (ii) this is an artefact of the lack of pantropical taxonomic syntheses and sampling in previous phylogenies and the consequent failure to identify clades that span the Old World and New World or conversely amphi-Atlantic genera that are non-monophyletic, both of which are critical for delimiting genera across this large pantropical clade. Finally, we discuss taxon delimitation in the phylogenomic era and especially how assessing patterns of gene tree conflict can provide additional insights into generic delimitation. This new phylogenomic framework provides the foundations for a series of papers reclassifying genera that are presented here in Advances in Legume Systematics (ALS) 14 Part 1, for establishing a new higher-level phylogenetic tribal and clade-based classification of Caesalpinioideae that is the focus of ALS14 Part 2 and for downstream analyses of evolutionary diversification and biogeography of this important group of legumes which are presented elsewhere.

Rampant nuclear-mitochondrial-plastid phylogenomic discordance in globally distributed calcifying microalgae

Incongruent phylogenies have been widely observed between nuclear and plastid or mitochondrial genomes in terrestrial plants and animals. However, few studies have examined these patterns in microalgae or the discordance between the two organelles. Here we investigated the nuclear-mitochondrial-plastid phylogenomic incongruence in Emiliania-Gephyrocapsa, a group of cosmopolitan calcifying phytoplankton with enormous populations and recent speciations. We assembled mitochondrial and plastid genomes of 27 strains from across global oceans and temperature regimes, and analyzed the phylogenomic histories of the three compartments using concatenation and coalescence methods. Six major clades with varying morphology and distribution are well recognized in the nuclear phylogeny, but such relationships are absent in the mitochondrial and plastid phylogenies, which also differ substantially from each other. The rampant phylogenomic discordance is due to a combination of organellar capture (introgression), organellar genome recombination, and incomplete lineage sorting of ancient polymorphic organellar genomes. Hybridization can lead to replacements of whole organellar genomes without introgression of nuclear genes and the two organelles are not inherited as a single cytoplasmic unit. This study illustrates the convoluted evolution and inheritance of organellar genomes in isogamous haplodiplontic microalgae and provides a window into the phylogenomic complexity of marine unicellular eukaryotes.

New Insights Into the Relationships Within Subtribe Scorzonerinae (Cichorieae, Asteraceae) Using Hybrid Capture Phylogenomics (Hyb-Seq)

Subtribe Scorzonerinae (Cichorieae, Asteraceae) contains 12 main lineages and approximately 300 species. Relationships within the subtribe, either at inter- or intrageneric levels, were largely unresolved in phylogenetic studies to date, due to the lack of phylogenetic signal provided by traditional Sanger sequencing markers. In this study, we employed a phylogenomics approach (Hyb-Seq) that targets 1,061 nuclear-conserved ortholog loci designed for Asteraceae and obtained chloroplast coding regions as a by-product of off-target reads. Our objectives were to evaluate the potential of the Hyb-Seq approach in resolving the phylogenetic relationships across the subtribe at deep and shallow nodes, investigate the relationships of major lineages at inter- and intrageneric levels, and examine the impact of the different datasets and approaches on the robustness of phylogenetic inferences. We analyzed three nuclear datasets: exon only, excluding all potentially paralogous loci; exon only, including loci that were only potentially paralogous in 1-3 samples; exon plus intron regions (supercontigs); and the plastome CDS region. Phylogenetic relationships were reconstructed using both multispecies coalescent and concatenation (Maximum Likelihood and Bayesian analyses) approaches. Overall, our phylogenetic reconstructions recovered the same monophyletic major lineages found in previous studies and were successful in fully resolving the backbone phylogeny of the subtribe, while the internal resolution of the lineages was comparatively poor. The backbone topologies were largely congruent among all inferences, but some incongruent relationships were recovered between nuclear and plastome datasets, which are discussed and assumed to represent cases of cytonuclear discordance. Considering the newly resolved phylogenies, a new infrageneric classification of Scorzonera in its revised circumscription is proposed.

Pervasive Phylogenomic Incongruence Underlies Evolutionary Relationships in Eyebrights (Euphrasia, Orobanchaceae)

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (Euphrasia), a group renowned for taxonomic complexity, with the aim of documenting the extent of phylogenetic discordance at both deep and at shallow phylogenetic scales. We generate whole-genome sequencing data and integrate this with prior genomic data to perform a comprehensive analysis of nuclear genomic, nuclear ribosomal (nrDNA), and complete plastid genomes from 57 individuals representing 36 Euphrasia species. The species tree analysis of 3,454 conserved nuclear scaffolds (46 Mb) reveals that at shallow phylogenetic scales postglacial colonization of North Western Europe occurred in multiple waves from discrete source populations, with most species not being monophyletic, and instead combining genomic variants from across clades. At a deeper phylogenetic scale, the Euphrasia phylogeny is structured by geography and ploidy, and partially by taxonomy. Comparative analyses show Southern Hemisphere tetraploids include a distinct subgenome indicative of independent polyploidy events from Northern Hemisphere taxa. In contrast to the nuclear genome analyses, the plastid genome phylogeny reveals limited geographic structure, while the nrDNA phylogeny is informative of some geographic and taxonomic affinities but more thorough phylogenetic inference is impeded by the retention of ancestral polymorphisms in the polyploids. Overall our results reveal extensive phylogenetic discordance at both deeper and shallower nodes, with broad-scale geographic structure of genomic variation but a lack of definitive taxonomic signal. This suggests that Euphrasia species either have polytopic origins or are maintained by narrow genomic regions in the face of extensive homogenizing gene flow. Moreover, these results suggest genome skimming will not be an effective extended barcode to identify species in groups such as Euphrasia, or many other postglacial species groups.

Increased resolution in the face of conflict: phylogenomics of the Neotropical bellflowers (Campanulaceae: Lobelioideae), a rapid plant radiation

BACKGROUND AND AIMS

The centropogonid clade (Lobelioideae: Campanulaceae) is an Andean-centred rapid radiation characterized by repeated convergent evolution of morphological traits, including fruit type and pollination syndromes. While previous studies have resolved relationships of lineages with fleshy fruits into subclades, relationships among capsular species remain unresolved. This lack of resolution has impeded reclassification of non-monophyletic genera, whose current taxonomy relies heavily on traits that have undergone convergent evolution.

METHODS

Targeted sequence capture using a probe-set recently developed for the centropogonid clade was used to obtain phylogenomic data from DNA extracted from both silica-dried and herbarium leaf tissue. These data were used to infer relationships among species using concatenated and partitioned species tree methods, and to quantify gene tree discordance.

KEY RESULTS

While silica-dried leaf tissue resulted in longer assembled sequence data, the inclusion of herbarium samples improved taxonomic representation. Relationships among baccate lineages are similar to those inferred in previous studies, although they differ for lineages within and among capsular clades. We improve the phylogenetic resolution of Siphocampylus, which forms ten groups of closely related species which we informally name. Two subclades of Siphocampylus and two individual species are rogue taxa whose placement differs widely across analyses. Gene tree discordance (including cytonuclear discordance) is rampant.

CONCLUSIONS

This first phylogenomic study of the centropogonid clade considerably improves our understanding of relationships in this rapid radiation. Differences across analyses and the possibility of additional lineage discoveries still hamper a solid and stable reclassification. Rapid morphological innovation corresponds with a high degree of phylogenomic complexity, including cytonuclear discordance, nuclear gene tree conflict and well-supported differences between analyses based on different nuclear loci. Together, these results point to a potential role of hemiplasy underlying repeated convergent evolution. This hallmark of rapid radiations is probably present in many other species-rich Andean plant radiations.

Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family

BACKGROUND

Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes.

RESULTS

We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae.

CONCLUSIONS

This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae)

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Phylogeny and biogeography of the northern temperate genus Dracocephalum s.l. (Lamiaceae)

The northern temperate genus Dracocephalum consists of approximately 70 species mainly distributed in the steppe-desert biomes of Central and West Asia and the alpine region of the Qinghai-Tibetan Plateau (QTP). Previous work has shown that Dracocephalum is not monophyletic and might include Hyssopus and Lallemantia. This study attempts to clarify the phylogenetic relationships, diversification patterns, and the biogeographical history of the three genera (defined as Dracocephalum s.l.). Based on a sampling of 66 taxa comprising more than 80% from extant species of Dracocephalum s.l., morphological, phylogenetic (maximum parsimony, likelihood, and Bayesian inference based on nuclear ITS and ETS, plastid rpl32-trnL, trnL-trnF, ycf1, and ycf1-rps15, and two low-copy nuclear markers AT3G09060 and AT1G09680), molecular dating, diversification, and ancestral range estimation analyses were carried out. Our results demonstrate that both Hyssopus and Lallemantia are embedded within Dracocephalum and nine well-supported clades can be recognized within Dracocephalum s.l. Analyses of divergence times suggest that the genus experienced an early rapid radiation during the middle to late Miocene with major lineages diversifying within a relatively narrow timescale. Ancestral area reconstruction analyses indicate that Dracocephalum s.l. originated in Central and West Asia and southern Siberia, and dispersed from Central and West Asia into the QTP and adjacent areas twice independently during the Pliocene. The aridification of the Asian interior possibly promoted the rapid radiation of Dracocephalum within this region, and the uplift of the QTP appears to have triggered the dispersal and recent rapid diversification of the genus in the QTP and adjacent regions. Combining molecular phylogenetic and morphological evidence, a revised infrageneric classification of Dracocephalum s.l. is proposed, which recognizes nine sections within the genus.

Comparative Plastome Analysis of Three Amaryllidaceae Subfamilies: Insights into Variation of Genome Characteristics, Phylogeny, and Adaptive Evolution

In the latest APG IV classification system, Amaryllidaceae is placed under the order of Asparagus and includes three subfamilies: Agapanthoideae, Allioideae, and Amaryllidoideae, which include many economically important crops. With the development of molecular phylogeny, research on the phylogenetic relationship of Amaryllidaceae has become more convenient. However, the current comparative analysis of Amaryllidaceae at the whole chloroplast genome level is still lacking. In this study, we sequenced 18 Allioideae plastomes and combined them with publicly available data (a total of 41 plastomes), including 21 Allioideae species, 1 Agapanthoideae species, 14 Amaryllidoideae species, and 5 Asparagaceae species. Comparative analyses were performed including basic characteristics of genome structure, codon usage, repeat elements, IR boundary, and genome divergence. Phylogenetic relationships were detected using single-copy genes (SCGs) and ribosomal internal transcribed spacer sequences (ITS), and the branch-site model was also employed to conduct the positive selection analysis. The results indicated that all Amaryllidaceae species showed a highly conserved typical tetrad structure. The GC content and five codon usage indexes in Allioideae species were lower than those in the other two subfamilies. Comparison analysis of Bayesian and ML phylogeny based on SCGs strongly supports the monophyly of three subfamilies and the sisterhood among them. Besides, positively selected genes (PSGs) were detected in each of the three subfamilies. Almost all genes with significant posterior probabilities for codon sites were associated with self-replication and photosynthesis. Our study investigated the three subfamilies of Amaryllidaceae at the whole chloroplast genome level and suggested the key role of selective pressure in the adaptation and evolution of Amaryllidaceae.

Reassessment of the Phylogeny and Systematics of Chinese Parnassia (Celastraceae): A Thorough Investigation Using Whole Plastomes and Nuclear Ribosomal DNA

Parnassia L., a perennial herbaceous genus in the family Celastraceae, consists of about 60 species and is mainly distributed in the Pan-Himalayan and surrounding mountainous regions. The taxonomic position and phylogenetic relationships of the genus are still controversial. Herein, we reassessed the taxonomic status of Parnassia and its intra- and inter-generic phylogeny within Celastraceae. To that end, we sequenced and assembled the whole plastid genomes and nuclear ribosomal DNA (nrDNA) of 48 species (74 individuals), including 25 species of Parnassia and 23 species from other genera of Celastraceae. We integrated high throughput sequence data with advanced statistical toolkits and performed the analyses. Our results supported the Angiosperm Phylogeny Group IV (APG IV) taxonomy which kept the genus to the family Celastraceae. Although there were topological conflicts between plastid and nrDNA phylogenetic trees, Parnassia was fully supported as a monophyletic group in all cases. We presented a first attempt to estimate the divergence of Parnassia, and molecular clock analysis indicated that the diversification occurred during the Eocene. The molecular phylogenetic results confirmed numerous taxonomic revisions, revealing that the morphological characters used in Parnassia taxonomy and systematics might have evolved multiple times. In addition, we speculated that hybridization/introgression might exist during genus evolution, which needs to be further studied. Similarly, more in-depth studies will clarify the diversification of characters and species evolution models of this genus.

Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridization

BACKGROUND AND AIMS

Abelia (Caprifoliaceae) is a small genus with five species, including one artificial hybrid and several natural hybrids. The genus has a discontinuous distribution in Mainland China, Taiwan Island and the Ryukyu Islands, providing a model system to explore the mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain.

METHODS

We reconstructed the phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction and ecological niche modelling (ENM) were used to examine the diversification history of Abelia.

KEY RESULTS

We found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times and morphology, we propose to merge Abelia macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been multiple widespread hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in the high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing Taiwan Island during the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region and Taiwan Island is probably the consequence of topographical isolation and postglacial contraction.

CONCLUSIONS

We used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results suggest that shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction between Mainland China and Taiwan Island. This study provides important new insights into the speciation process and taxonomy of Abelia.

Synthesis of Nuclear and Chloroplast Data Combined With Network Analyses Supports the Polyploid Origin of the Apple Tribe and the Hybrid Origin of the Maleae-Gillenieae Clade

Plant biologists have debated the evolutionary origin of the apple tribe (Maleae; Rosaceae) for over a century. The "wide-hybridization hypothesis" posits that the pome-bearing members of Maleae (base chromosome number x = 17) resulted from a hybridization and/or allopolyploid event between progenitors of other tribes in the subfamily Amygdaloideae with x = 8 and x = 9, respectively. An alternative "spiraeoid hypothesis" proposed that the x = 17 of Maleae arose via the genome doubling of x = 9 ancestors to x = 18, and subsequent aneuploidy resulting in x = 17. We use publicly available genomic data-448 nuclear genes and complete plastomes-from 27 species representing all major tribes within the Amygdaloideae to investigate evolutionary relationships within the subfamily containing the apple tribe. Specifically, we use network analyses and multi-labeled trees to test the competing wide-hybridization and spiraeoid hypotheses. Hybridization occurred between an ancestor of the tribe Spiraeeae (x = 9) and an ancestor of the clade Sorbarieae (x = 9) + Exochordeae (x = 8) + Kerrieae (x = 9), giving rise to the clade Gillenieae (x = 9) + Maleae (x = 17). The ancestor of the Maleae + Gillenieae arose via hybridization between distantly related tribes in the Amygdaloideae (i.e., supporting the wide hybridization hypothesis). However, some evidence supports an aspect of the spiraeoid hypothesis-the ancestors involved in the hybridization event were likely both x = 9, so genome doubling was followed by aneuploidy to result in x = 17 observed in Maleae. By synthesizing existing genomic data with novel analyses, we resolve the nearly century-old mystery regarding the origin of the apple tribe. Our results also indicate that nuclear gene tree-species tree conflict and/or cytonuclear conflict are pervasive at several other nodes in subfamily Amygdaloideae of Rosaceae.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.

Mojiangia oreophila (Crepidinae, Cichorieae, Asteraceae), a new species and genus from Mojiang County, SW Yunnan, China, and putative successor of the maternal Faberia ancestor

A single small population of chasmophytic plants is described as M ojiangia oreophila, a monotypic genus in the subtribe Crepidinae, characterised by a unique combination of morphological features, in particular densely long-papillose homomorphic achenes with five main ribs each accompanied by two secondary ribs, coarse brownish pappus bristles, moderately many-flowered capitula, a small involucre with numerous outer phyllaries, perennial rosette herb growth and brown-woolly caudex and leaf axils. Molecular phylogenetic analysis detected that in the nrITS phylogeny M. oreophila forms a clade of its own in the Crepidinae; in the plastid DNA phylogeny it is nested in the clade formed by the hybridogenous genus Faberia, the maternal ancestor of which comes from the Crepidinae and the paternal ancestor from the Lactucinae, where Faberia is placed in nrITS phylogenies. M. oreophila shares several morphological features with Faberia and also shares the expected chromosome number of 2n = 16 with its hitherto unknown maternal ancestor. M. oreophila may therefore be a successor of the maternal ancestor of Faberia. Alternatively, cytonuclear discordance is to be assumed in Mojiangia, caused by chloroplast capture as a result of hybridisation and introgression with Faberia.

Contrasting gene-level signatures of selection with reproductive fitness

Selection leaves signatures in the DNA sequence of genes, with many test statistics devised to detect its action. While these statistics are frequently used to support hypotheses about the adaptive significance of particular genes, the effect these genes have on reproductive fitness is rarely quantified experimentally. Consequently, it is unclear how gene-level signatures of selection are associated with empirical estimates of gene effect on fitness. Eukaryotic datasets that permit this comparison are very limited. Using the model plant Arabidopsis thaliana, for which these resources are available, we calculated seven gene-level substitution and polymorphism-based statistics commonly used to infer selection (dN/dS, NI, DOS, Tajima's D, Fu and Li's D*, Fay and Wu's H, and Zeng's E) and, using knockout lines, compared these to gene-level estimates of effect on fitness. We found that consistent with expectations, essential genes were more likely to be classified as negatively selected. By contrast, using 379 Arabidopsis genes for which data was available, we found no evidence that genes predicted to be positively selected had a significantly different effect on fitness than genes evolving more neutrally. We discuss these results in the context of the analytic challenges posed by Arabidopsis, one of the only systems in which this study could be conducted, and advocate for examination in additional systems. These results are relevant to the evaluation of genome-wide studies across species where experimental fitness data is unavailable, as well as highlighting an increasing need for the latter.

Hybridization Dynamics and Extensive Introgression in the Daphnia longispina Species Complex: New Insights from a High-Quality Daphnia galeata Reference Genome

Hybridization and introgression are recognized as an important source of variation that influence adaptive processes; both phenomena are frequent in the genus Daphnia, a keystone zooplankton taxon in freshwater ecosystems that comprises several species complexes. To investigate genome-wide consequences of introgression between species, we provide here the first high-quality genome assembly for a member of the Daphnia longispina species complex, Daphnia galeata. We further resequenced 49 whole genomes of three species of the complex and their interspecific hybrids both from genotypes sampled in the water column and from single resting eggs extracted from sediment cores. Populations from habitats with diverse ecological conditions offered an opportunity to study the dynamics of hybridization linked to ecological changes and revealed a high prevalence of hybrids. Using phylogenetic and population genomic approaches, we provide first insights into the intra- and interspecific genome-wide variability in this species complex and identify regions of high divergence. Finally, we assess the length of ancestry tracts in hybrids to characterize introgression patterns across the genome. Our analyses uncover a complex history of hybridization and introgression reflecting multiple generations of hybridization and backcrossing in the Daphnia longispina species complex. Overall, this study and the new resources presented here pave the way for a better understanding of ancient and contemporary gene flow in the species complex and facilitate future studies on resting egg banks accumulating in lake sediment.

Sage Insights Into the Phylogeny of Salvia: Dealing With Sources of Discordance Within and Across Genomes

Next-generation sequencing technologies have facilitated new phylogenomic approaches to help clarify previously intractable relationships while simultaneously highlighting the pervasive nature of incongruence within and among genomes that can complicate definitive taxonomic conclusions. Salvia L., with ∼1,000 species, makes up nearly 15% of the species diversity in the mint family and has attracted great interest from biologists across subdisciplines. Despite the great progress that has been achieved in discerning the placement of Salvia within Lamiaceae and in clarifying its infrageneric relationships through plastid, nuclear ribosomal, and nuclear single-copy genes, the incomplete resolution has left open major questions regarding the phylogenetic relationships among and within the subgenera, as well as to what extent the infrageneric relationships differ across genomes. We expanded a previously published anchored hybrid enrichment dataset of 35 exemplars of Salvia to 179 terminals. We also reconstructed nearly complete plastomes for these samples from off-target reads. We used these data to examine the concordance and discordance among the nuclear loci and between the nuclear and plastid genomes in detail, elucidating both broad-scale and species-level relationships within Salvia. We found that despite the widespread gene tree discordance, nuclear phylogenies reconstructed using concatenated, coalescent, and network-based approaches recover a common backbone topology. Moreover, all subgenera, except for Audibertia, are strongly supported as monophyletic in all analyses. The plastome genealogy is largely resolved and is congruent with the nuclear backbone. However, multiple analyses suggest that incomplete lineage sorting does not fully explain the gene tree discordance. Instead, horizontal gene flow has been important in both the deep and more recent history of Salvia. Our results provide a robust species tree of Salvia across phylogenetic scales and genomes. Future comparative analyses in the genus will need to account for the impacts of hybridization/introgression and incomplete lineage sorting in topology and divergence time estimation.