Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

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

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

Extensive genome-wide phylogenetic discordance is due to incomplete lineage sorting in the rapidly radiated East Asian genus Nekemias (Vitaceae)

BACKGROUND AND AIMS

Nekemias is a small genus of the grape family, with nine species discontinuously distributed in temperate to subtropical zones of the Northern Hemisphere but mostly in East Asia. Previous phylogenetic studies on Nekemias have mainly been based on a few chloroplast markers, and the phylogenetic framework and systematic relationships are still highly contested.

METHODS

We carried out a systematic framework reconstruction of Nekemias and intra-generic reticulate evolutionary analyses based on extensive single-copy nuclear and chloroplast genomic data obtained by the Hyb-Seq approach, combining genome skimming and target enrichment.

KEY RESULTS

Both nuclear and chloroplast genomic data strongly support the monophyly of Nekemias with its division into two major lineages from East Asia and North America, respectively. There are strong and extensive topological conflicts among nuclear gene trees and between nuclear and chloroplast topologies within the genus, especially within the East Asian clade.

CONCLUSIONS

Rapid radiation through predominant incomplete lineage sorting (ILS) throughout the evolutionary history of the East Asian taxa is supported to explain the relatively high species diversity of Nekemias in East Asia. This study highlights the important role of short periods of rapid evolutionary radiations accompanied by ILS as a mechanism for the complex and fast species diversifications in the grape family as well as potentially in many other plant lineages in East Asia and beyond.

Disentangling a genome-wide mosaic of conflicting phylogenetic signals in Western Rattlesnakes

Species tree inference is often assumed to be more accurate as datasets increase in size, with whole genomes representing the best-case-scenario for estimating a single, most-likely speciation history with high confidence. However, genomes may harbor a complex mixture of evolutionary histories among loci, which amplifies the opportunity for model misspecification and impacts phylogenetic inference. Accordingly, multiple distinct and well-supported phylogenetic trees are often recovered from genome-scale data, and approaches for biologically interpreting these distinct signatures are a major challenge for evolutionary biology in the age of genomics. Here, we analyze 32 whole genomes of nine taxa and two outgroups from the Western Rattlesnake species complex. Using concordance factors, topology weighting, and concatenated and species tree analyses with a chromosome-level reference genome, we characterize the distribution of phylogenetic signal across the genomic landscape. We find that concatenated and species tree analyses of autosomes, the Z (sex) chromosome, and mitochondrial genome yield distinct, yet strongly supported phylogenies. Analyses of site-specific likelihoods show additional patterns consistent with rampant model misspecification, a likely consequence of several evolutionary processes. Together, our results suggest that a combination of historic and recent introgression, along with natural selection, recombination rate variation, and cytonuclear co-evolution of nuclear-encoded mitochondrial genes, underlie genome-wide variation in phylogenetic signal. Our results highlight both the power and complexity of interpreting whole genomes in a phylogenetic context and illustrate how patterns of phylogenetic discordance can reveal the impacts of different evolutionary processes that contribute to genome-wide variation in phylogenetic signal.

Unraveling the extensive phylogenetic discordance and evolutionary history of spurless taxa within the Aquilegia ecalcarata complex

Parallel evolution of the same, or at least very similar, phenotype(s) in different lineages is often interpreted as evidence for the action of natural selection. However, caution is required when inferring parallel evolution based on uncertain or potentially incorrect phylogenetic frameworks. Here, by conducting extensive phylogenomic and population genetic analyses, we aim to clarify the evolutionary history of spurless taxa within the Aquilegia ecalcarata complex. We observed substantial discordance in the phylogenetic patterns across the entire genome, primarily attributed to ancient introgression and incomplete lineage sorting. Additionally, we identified several spurless lineages whose phylogenetic positions were distorted by admixture events. Using a backbone tree and demographic modeling, we determined that these spurless taxa independently originated twice within this group. Intriguingly, our investigation revealed that the spurless taxa experienced population expansion during global cooling, while their spurred sister groups underwent population contraction. The parallel losses of petal spurs, therefore, may be linked to adaptations for low-temperature conditions. These findings emphasize the importance of comprehensive population-level analyses in phylogenetic inference and provide valuable insights into the dynamics of trait loss and its implications for the adaptive strategies.

Intergeneric and interspecific relationships in tribe Ricineae revealed by phylogenomics of the plastome and transcriptome

Introduction

The taxonomy of Euphorbiaceae is extremely difficult, especially the phylogeny of closely related genera. In Ricinus, which embraces an important non-food oil-seed crop worldwide, Discocleidion and Speranskia are closely related genera based on molecular evidence (tribe Ricineae), however the intergeneric and interspecific relationship of the tribe is not well-resolved.

Methods

Plastome and transcriptome were sequenced and assembled before maximum likelihood and Bayesian inference phylogenetic trees were reconstructed. Plastome features and comparative analyses were conducted. Morphological traits of the tribe were explored as supplement to the molecular data.

Results

The newly sequenced plastomes ranged from 167,327 to 190,093 bp with typical circular quadripartite structures. The longest genome of S. tuberculata may due to higher number of simple sequence repeats. Natural selection pressure on chloroplast genes was relatively small and the tribe likely experienced a population contraction. The transcriptome assembly contig N50 of the tribe ranged from 1506 (D. rufescens) to 2489 bp (S. tuberculata). A total of 50,513 genes (S. cantonensis) to 78,048 genes (D. ulmifolium) were detected, and the GC content varied between 38.17% (S. cantonensis) and 40.01% (R. communis). The three genera formed a well-supported monophyletic lineage, confirmed by different genomic data using different methods. Discocleidion and Ricinus were supported to be closely related. In Speranskia, S. yunnanensis diverged first and the divergence of S. tuberculata and S. cantonensis was followed. Further, morphological similarities supported the monophyletic lineage and intergeneric and interspecific relationship.

Discussion

The relationship in the tribe Ricineae is clearly revealed by genomic and morphological data, providing a genetic basis for future comparative genomic investigations and phylogeny reconstruction of Euphorbiaceae.

Standing genetic variation and introgression shape the cryptic radiation of Aquilegia in the mountains of Southwest China

Cryptic diversity in evolutionary radiation offers an excellent system for investigating the intricacies of evolutionary progress. Understanding the evolution of cryptic diversity is imperative for unraveling the hidden complexities of biodiversity. However, empirical evidence elucidating the mechanisms behind cryptic radiation remains limited, particularly in plants. Here, we focus on a monophyletic group of Aquilegia species mainly distributed in the mountains of Southwest China, one of the world's biodiversity hotspots. Using whole-genome resequencing of 158 individuals from 23 natural populations, we identify three to four paraphyletic lineages within each morphological species. Our findings reveal that 39 out of 43 detected instances of introgression occurred post-lineage formation. Identifying shared genomic regions indicates that the divergence of fixed singletons in lineages from morphological species A. kansuensis and A. rockii predates lineage formation, supporting a scenario where incomplete lineage sorting of standing variation contributes to morphological parallelism. Furthermore, strong positive correlations among genomic differentiation, divergence, and introgression suggest that standing variations and introgression from non-sister lineages contribute to the rapid genetic divergence. Our study illuminates the important roles of standing variations and introgression in plant cryptic radiation, advancing our understanding of the complex mechanisms behind the evolution of biodiversity in recent radiation events.

Phylogenomics of Paragymnopteris (Cheilanthoideae, Pteridaceae): Insights from plastome, mitochondrial, and nuclear datasets

Previous studies have shown that at least six genera of the Cheilanthoideae, a subfamily of the fern family Pteridaceae, may not be monophyletic. In these non-monophyletic genera, the Old-World genus Paragymnopteris including approximately five species have long been controversial. In this study, with an extensive taxon sampling of Paragymnopteris, we assembled 19 complete plastomes of all recognized Paragymnopteris species, plastomes of Pellaea (3 species) and Argyrochosma (1 species), as well as transcriptomes from Paragymnopteris (6 species) and Argyrochosma (1 species). We conducted a comprehensive and systematic phylogenomic analysis focusing on the contentious relationships among the genus of Paragymnopteris through 9 plastid makers, the plastomes, mitochondria, nuclear ribosomal cistron genomes, and single-copy nuclear genes. Moreover, we further combined distribution, ploidy, and morphological features to investigate the evolution of Paragymnopteris. The backbone of Paragymnopteris was resolved consistently in the nuclear and plastid phylogenies. Our major results include: (1) Paragymnopteris is not monophyletic including two fully supported clades; (2) confirming that Paragymnopteris delavayi var. intermedia is a close relative of P. delavayi instead of P. marantae var. marantae; (3) the chromosome base number may not be a stable trait which has previously been used as an important character to divide Paragymnopteris into two groups; and (4) gene flow or introgression might be the main reason for the gene trees conflict of Paragymnopteris, but both gene flow and ILS might simultaneously and/or cumulatively act on the conflict of core pellaeids. The robust phylogeny of Paragymnopteris presented here will help us for the future studies of the arid to semi-arid ferns of Cheilanthoideae at the evolutionary, physiological, developmental, and omics-based levels.

Phylotranscriptomic Analyses Resolve Evolutionary History of Eremopyrum (Triticeae; Poaceae)

Disentangling the phylogenetic relationship of polyploid species is essential for understanding how such polyploid species evolved following their origin. To investigate the speciation and evolutionary history of Eremopyrum, we analyzed 36 transcriptomes from 9 polyploid accessions of Eremopyrum and 27 diploid taxa representing 12 basic genomes in Triticeae. Phylogenetic reconstruction, divergence time, and introgression event demonstrated that (1) Eremopyrum and Agropyron shared a common ancestor; (2) Eremopyrum has undergone ongoing evolutionary diversification since its origin in Late Miocene; (3) the diploid E. triticeum and E. distans were the genome donors of the tetraploid species of Eremopyrum; (4) both Eremopyrum and Agropyron contribute to the nonmonophyletic origin of tetraploid E. orientale via introgression events. Our results shed new light on our understanding of the diversity and ecological adaptation of the species in Eremopyrum.

Subfamily evolution analysis using nuclear and chloroplast data from the same reads

The chloroplast (cp) genome is a widely used tool for exploring plant evolutionary relationships, yet its effectiveness in fully resolving these relationships remains uncertain. Integrating cp genome data with nuclear DNA information offers a more comprehensive view but often requires separate datasets. In response, we employed the same raw read sequencing data to construct cp genome-based trees and nuclear DNA phylogenetic trees using Read2Tree, a cost-efficient method for extracting conserved nuclear gene sequences from raw read data, focusing on the Aurantioideae subfamily, which includes Citrus and its relatives. The resulting nuclear DNA trees were consistent with existing nuclear evolutionary relationships derived from high-throughput sequencing, but diverged from cp genome-based trees. To elucidate the underlying complex evolutionary processes causing these discordances, we implemented an integrative workflow that utilized multiple alignments of each gene generated by Read2Tree, in conjunction with other phylogenomic methods. Our analysis revealed that incomplete lineage sorting predominantly drives these discordances, while introgression and ancient introgression also contribute to topological discrepancies within certain clades. This study underscores the cost-effectiveness of using the same raw sequencing data for both cp and nuclear DNA analyses in understanding plant evolutionary relationships.

Chromosome-Level Genome Assembly and Whole-Genome Resequencing Revealed Contrasting Population Genetic Differentiation of Black Bream (Megalobrama skolkovii) (Teleostei: Cyprinidae) Allopatric and Sympatric to Its Kin Species

The black bream (Megalobrama skolkovii) is an economically important species widely distributed in China, with its geographic populations potentially having undergone differentiations and local adaptations. In this study, we presented a chromosome-level genome assembly of this species and investigated genetic differentiations of its populations that are allopatric (the northern one) and sympatric (the Poyang Lake) to its kin species, the blunt-snout bream (M. amblycephala), using whole genome resequencing analysis. The results showed that the genome size of black bream was 1.13 Gb, very similar to its kin species but larger than its close relatives, the four Chinese major carps. By resequencing individuals from the northern and Poyang Lake populations, we found that the northern population showed lower genetic diversity, larger genetic differentiation, and two sharp historical declines in population size through demographic analysis, indicating the possible bottlenecks after the allopatric isolation. In contrast, the Poyang Lake population, with its higher genetic diversity, higher Tajima's D value, and lower levels of linkage disequilibrium, reflects the ancestral state of black bream. In addition, we also found that the northern population shared more alleles with its kin species, indicating it may retain more ancestral variations. This was further analyzed to be caused by incomplete lineage sorting and ancient introgression. Some key genes related to reproductive processes, body size development, and muscle metabolism were found under selection in the northern population, possibly responsible for its local adaptation. Our findings that the black bream allopatric population had a loss of genetic diversity but retained more ancestral variations can expand our knowledge on population genetic differentiation and give us hints for future genetic conservation.

Contrast and Genomic Characterisation of Ancient and Recent Interspecific Introgression Between Deeply Diverged Moustache Toads (Leptobrachium)

Recent genomic analyses have provided new insights into the process of interspecific introgression and its consequences on species evolution. Most recent studies, however, focused on hybridization between recently radiated species, with few examining the genomic outcomes of ancient hybridization across deeply diverged species. Using whole genome data of moustache toads (Leptobrachium), we identified signals of three hybridization events among nine species that diverged at the Eocene. An ancient introgression from L. leishanense to the ancestral branch (C1) of L. liui introduced adaptive variants. The highly introgressed regions include genes with important functions in odorant detection and immune responses. These genes are preserved in all three descendent populations of L. liui_C1, and these regions likely have been positively selected over a long filtering process. A recent introgression occurred from L. huashen to L. tengchongense, with the introgressed regions being mostly neutral. Furthermore, one F1 hybrid individual was detected between sympatric L. ailaonicum and L. promustache. The signals of introgression largely disappeared after removing the hybrid individual, indicating an occasional hybridization but minimal introgression. Further examination of highly divergent but low introgressed genomic regions revealed both pre-mating isolation and genetic incompatibility as potential mechanisms of resisting introgression and maintaining species boundaries. Additionally, no large X-effect was found in these introgression events. Hybridization between deeply diverged amphibian species may be common, but detectable introgressions are likely less so, with recent introgression being mostly neutral and the rare ancient one potentially adaptive. Our findings complement recent genomic work, and together they provide a better understanding of the genomic characteristics of interspecific introgression and its significance in species adaptation and evolution.

Deciphering the evolution and biogeography of ant-ferns Lecanopteris s.s

Southeast Asia is a biodiversity hotspot characterized by a complex paleogeography, and its Polypodiopsida flora is particularly diverse. While hybridization is recognized as common in ferns, further research is needed to investigate the relationship between hybridization events and fern diversity. Lecanopteris s.s., an ant-associated fern, has been subject to debate regarding species delimitations primarily due to limited DNA markers and species sampling. Our study integrates 22 newly generated plastomes, 22 transcriptomes, and flow cytometry of all native species along with two cultivated hybrids. Our objective is to elucidate the reticulate evolutionary history within Lecanopteris s.s. through the integration of phylobiogeographic reconstruction, gene flow inference, and genome size estimation. Key findings of our study include: (1) An enlarged plastome size (178-187 Kb) in Lecanopteris s.s., attributed to extreme expansion of the Inverted Repeat (IR) regions; (2) The traditional 'pumila' and 'crustacea' groups are paraphyletic; (3) Significant cytonuclear discordance attributed to gene flow; (4) Natural hybridization and introgression in the 'pumila' and 'darnaedii' groups; (5) L. luzonensis is the maternal parent of L. 'Yellow Tip', with L. pumila suggested as a possible paternal parent; (6) L. 'Tatsuta' is a hybrid between L. luzonensis and L. crustacea; (7) Lecanopteris s.s. first diverged during the Neogene and then during the middle Miocene climatic optimum in the Indochina and Sundaic regions. In conclusion, the biogeographic history and speciation of Lecanopteris have been profoundly shaped by past climate changes and geodynamics of Southeast Asia. Dispersals, hybridization and introgression between species act as pivotal factors in the evolutionary trajectory of Lecanopteris s.s.. This research provides a robust framework for further exploration and understanding of the complex dynamics driving the diversification and distribution patterns within Polypodiaceae subfamily Microsoroideae.

Phylogenomics and Pervasive Genome-Wide Phylogenetic Discordance Among Fin Whales (Balaenoptera physalus)

Phylogenomics has the power to uncover complex phylogenetic scenarios across the genome. In most cases, no single topology is reflected across the entire genome as the phylogenetic signal differs among genomic regions due to processes, such as introgression and incomplete lineage sorting. Baleen whales are among the largest vertebrates on Earth with a high dispersal potential in a relatively unrestricted habitat, the oceans. The fin whale (Balaenoptera physalus) is one of the most enigmatic baleen whale species, currently divided into four subspecies. It has been a matter of debate whether phylogeographic patterns explain taxonomic variation in fin whales. Here we present a chromosome-level whole genome analysis of the phylogenetic relationships among fin whales from multiple ocean basins. First, we estimated concatenated and consensus phylogenies for both the mitochondrial and nuclear genomes. The consensus phylogenies based upon the autosomal genome uncovered monophyletic clades associated with each ocean basin, aligning with the current understanding of subspecies division. Nevertheless, discordances were detected in the phylogenies based on the Y chromosome, mitochondrial genome, autosomal genome and X chromosome. Furthermore, we detected signs of introgression and pervasive phylogenetic discordance across the autosomal genome. This complex phylogenetic scenario could be explained by a puzzle of introgressive events, not yet documented in fin whales. Similarly, incomplete lineage sorting and low phylogenetic signal could lead to such phylogenetic discordances. Our study reinforces the pitfalls of relying on concatenated or single locus phylogenies to determine taxonomic relationships below the species level by illustrating the underlying nuances that some phylogenetic approaches may fail to capture. We emphasize the significance of accurate taxonomic delineation in fin whales by exploring crucial information revealed through genome-wide assessments.

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

Phylogeography of Pleurospermum foetens (Apiaceae) From the Sky Islands of Southwest China

Sky islands provide insights on how glacial-interglacial cycles have shaped species distribution and help for predicting species' responses to climate warming. The alpine subnival belt of southwest China, especially in the Hengduan Mountains and adjacent areas, is sky island-like. Among them, the Yunnan-Kweichow Plateau harbors several isolated mountains with well-developed alpine subnival vegetation, sharing a similar species composition with the Hengduan Mountains. However, the relationship between the sky islands of the Hengduan Mountains and the Yunnan-Kweichow Plateau remains insufficiently explored. Pleurospermum foetens (Apiaceae) is a species endemic to the alpine screes of the Yunnan-Kweichow Plateau and the Hengduan Mountains. We used DNA sequence data from 59 individuals across 9 populations, combined with ecological niche modeling, to investigate the evolution history and future distribution of P. foetens within this sky island region. The results indicate the following: (1) P. foetens exhibits a significant phylogeographic structure and can be classified into three nrDNA clades and two cpDNA clades, respectively, (2) a nuclear-plastid discordance observed in P. foetens and its relatives based on phylogenetic analysis. P. foetens is monophyletic in the nrDNA phylogeny, while two major clades (HDM and YGP) are present in the cpDNA phylogeny, each forming a clade with other congeneric species. (3) Ecological niche modeling of P. foetens indicated that the species had the most extensive suitable habitat during the last glacial maximum (LGM). However, anticipated climate warming in the coming decades is expected to reduce the suitable range of P. foetens, posing a significant threat to isolated marginal populations (e.g., Shizi Mountain) with restricted alpine scree habitats. In conclusion, our study highlights the substantial effect of sky island and glacial-interglacial cycles on the population divergence of P. foetens. Conservation efforts for marginal populations of alpine plants in the Yunnan-Kweichow Plateau require increased attention and prioritization.

Whole Genomes Reveal Evolutionary Relationships and Mechanisms Underlying Gene-Tree Discordance in Neodiprion Sawflies

Rapidly evolving taxa are excellent models for understanding the mechanisms that give rise to biodiversity. However, developing an accurate historical framework for comparative analysis of such lineages remains a challenge due to ubiquitous incomplete lineage sorting (ILS) and introgression. Here, we use a whole-genome alignment, multiple locus-sampling strategies, and summary-tree and single nucleotide polymorphism-based species-tree methods to infer a species tree for eastern North American Neodiprion species, a clade of pine-feeding sawflies (Order: Hymenopteran; Family: Diprionidae). We recovered a well-supported species tree that-except for three uncertain relationships-was robust to different strategies for analyzing whole-genome data. Nevertheless, underlying gene-tree discordance was high. To understand this genealogical variation, we used multiple linear regression to model site concordance factors estimated in 50-kb windows as a function of several genomic predictor variables. We found that site concordance factors tended to be higher in regions of the genome with more parsimony-informative sites, fewer singletons, less missing data, lower GC content, more genes, lower recombination rates, and lower D-statistics (less introgression). Together, these results suggest that ILS, introgression, and genotyping error all shape the genomic landscape of gene-tree discordance in Neodiprion. More generally, our findings demonstrate how combining phylogenomic analysis with knowledge of local genomic features can reveal mechanisms that produce topological heterogeneity across genomes.

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Characterization of genome-wide phylogenetic conflict uncovers evolutionary modes of carnivorous fungi

Mass extinction has often paved the way for rapid evolutionary radiation, resulting in the emergence of diverse taxa within specific lineages. The emergence and diversification of carnivorous nematode-trapping fungi (NTF) in Ascomycota have been linked to the Permian-Triassic (PT) extinction, but the processes underlying NTF radiation remain unclear. We conducted phylogenomic analyses using 23 genomes that represent three NTF lineages, each employing distinct nematode traps-mechanical traps (Drechslerella spp.), three-dimensional (3D) adhesive traps (Arthrobotrys spp.), and two-dimensional (2D) adhesive traps (Dactylellina spp.), and the genome of one non-NTF species as the outgroup. These analyses revealed multiple mechanisms that likely contributed to the tempo of the NTF evolution and rapid radiation. The species tree of NTFs based on 2,944 single-copy orthologous genes suggested that Drechslerella emerged earlier than Arthrobotrys and Dactylellina. Extensive genome-wide phylogenetic discordance was observed, mainly due to incomplete lineage sorting (ILS) between lineages. Two modes of non-vertical evolution (introgression and horizontal gene transfer) also contributed to phylogenetic discordance. The ILS genes that are associated with hyphal growth and trap morphogenesis (e.g., those associated with the cell membrane system and polarized cell division) exhibited signs of positive selection.IMPORTANCEBy conducting a comprehensive phylogenomic analysis of 23 genomes across three NTF lineages, the research reveals how diverse evolutionary mechanisms, including ILS and non-vertical evolution (introgression and horizontal gene transfer), contribute to the swift diversification of NTFs. These findings highlight the complex evolutionary dynamics that drive the rapid radiation of NTFs, providing valuable insights into the processes underlying their diversity and adaptation.

Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae)

Obtaining a robust phylogeny proves challenging due to the intricate evolutionary history of species, where processes such as hybridization and incomplete lineage sorting can introduce conflicting signals, thereby complicating phylogenetic inference. In this study, we conducted comprehensive sampling of Elsholtzieae, with a particular focus on its largest genus, Elsholtzia. We utilized 503 nuclear loci and complete plastome sequences obtained from 99 whole-genome sequencing datasets to elucidate the interspecific relationships within the Elsholtzieae. Additionally, we explored various sources of conflicts between gene trees and species trees. Fully supported backbone phylogenies were recovered, and the monophyly of Elsholtzia and Keiskea was not supported. Significant gene tree heterogeneity was observed at numerous nodes, particularly regarding the placement of Vuhuangia and the E. densa clade. Further investigations into potential causes of this discordance revealed that incomplete lineage sorting (ILS), coupled with hybridization events, has given rise to substantial gene tree discordance. Several species, represented by multiple samples, exhibited a closer association with geographical distribution rather than following a strictly monophyletic pattern in plastid trees, suggesting chloroplast capture within Elsholtzieae and providing evidence of hybridization. In conclusion, this study provides phylogenomic insights to untangle taxonomic problems in the tribe Elsholtzieae, especially the genus Elsholtzia.

Linked Selection and Gene Density Shape Genome-Wide Patterns of Diversification in Peatmosses

Genome evolution under speciation is poorly understood in nonmodel and nonvascular plants, such as bryophytes-the largest group of nonvascular land plants. Their genomes are structurally different from angiosperms and likely subjected to stronger linked selection pressure, which may have profound consequences on genome evolution in diversifying lineages, even more so when their genome architecture is conserved. We use the highly diverse, rapidly radiated group of peatmosses (Sphagnum) to characterize the processes affecting genome diversification in bryophytes. Using whole-genome sequencing data from populations of 12 species sampled at different phylogenetic and geographical scales, we describe high correlation of the genomic landscapes of differentiation, divergence, and diversity in Sphagnum. Coupled with evidence from the patterns of covariation among different measures of genetic diversity, phylogenetic discordance, and gene density, this provides strong support that peatmoss genome evolution has been shaped by the long-term effects of linked selection, constrained by distribution of selection targets in the genome. Thus, peatmosses join the growing number of animal and plant groups where functional features of the genome, such as gene density, and linked selection drive genome evolution along predetermined and highly similar routes in different species. Our findings demonstrate the great potential of bryophytes for studying the genomics of speciation and highlight the urgent need to expand the genomic resources in this remarkable group of plants.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Morphological, genetic and ecological divergence in near-cryptic bryophyte species widespread in the Holarctic: the Dicranum acutifolium complex (Dicranales) revisited in the Alps

There is mounting evidence that reproductively isolated, but morphologically weakly differentiated species (so-called cryptic species) represent a substantial part of biological diversity, especially in bryophytes. We assessed the evolutionary history and ecological differentiation of a species pair, Dicranum brevifolium and D. septentrionale, which have overlapping ranges in the Holarctic. Despite their morphological similarity, we found similar genetic differentiation as between morphologically well-differentiated Dicranum species. Moreover, we detected gene tree discordance between plastid and nuclear markers, but neither of the two datasets resolved the two as sister species. The signal in trnL-trnF better reflects the morphological and ecological affinities and indicates a close relationship while ITS sequence data resolved the two taxa as phylogenetically distantly related. The discordance is probably unrelated to the ecological differentiation of D. septentrionale to colonise subneutral to alkaline substrates (vs. acidic in D. brevifolium), because this ability is rare in the genus and shared with D. acutifolium. This taxon is the closest relative of D. septentrionale according to the trnL-trnF data and does not share the discordance in ITS. We furthermore demonstrate that beside D. acutifolium, both D. septentrionale and D. brevifolium occur in the Alps but D. brevifolium is most likely rarer. Based on morphological analyses including factor analysis for mixed data of 45 traits we suggest treating the latter two as near-cryptic species and we recommend verifying morphological determinations molecularly.

A roadmap of phylogenomic methods for studying polyploid plant genera

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

A new nuclear phylogeny of the tea family (Theaceae) unravels rapid radiations in genus Camellia

Molecular analyses of rapidly radiating groups often reveal incongruence between gene trees. This mainly results from incomplete lineage sorting, introgression, and gene tree estimation error, which complicate the estimation of phylogenetic relationships. In this study, we reconstruct the phylogeny of Theaceae using 348 nuclear loci from 68 individuals and two outgroup taxa. Sequence data were obtained by target enrichment using the recently released Angiosperm 353 universal probe set applied to herbarium specimens. The robustness of the topologies to variation in data quality was established under a range of different filtering schemes, using both coalescent and concatenation approaches. Our results confirmed most of the previously hypothesized relationships among tribes and genera, while clarifying additional interspecific relationships within the rapidly radiating genus Camellia. We recovered a remarkably high degree of gene tree heterogeneity indicative of rapid radiation in the group and observed cytonuclear conflicts, especially within Camellia. This was especially pronounced around short branches, which we primarily associate with gene tree estimation error. Our analysis also indicates that incomplete lineage sorting (ILS) contributed to gene-tree conflicts and accounted for approximately 14 % of the explained variation, whereas inferred introgression levels were low. Our study advances the understanding of the evolution of this important plant family and provides guidance on the application of target capture methods and the evaluation of key processes that influence phylogenetic discordances.

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

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

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

The phylogenetic relationship and demographic history of rhesus macaques (Macaca mulatta) in subtropical and temperate regions, China

Pleistocene climatic oscillations exerted significant influences on the genetic structure and demography of rhesus macaque (Macaca mulatta) in eastern China. However, the evolutionary history of rhesus macaques in subtropical and temperate China remained unclear and/or controversial. Herein, we analyzed the autosomes, mitochondrial genomes, and Y-chromosomes from 84 individuals of Chinese rhesus macaque. The results revealed that (1) all individuals were clustered into pan-west and pan-east genetic groups, which exhibited Shaanxi Province as the northernmost region of western dispersal route of rhesus macaques in China; (2) in subtropical and temperate China, rhesus macaques were divided into four lineages (TH, DB, HS, and QL), and their divergence times corresponded to the Penultimate Glaciation (300-130 kya) and Last Glaciation (70-10 kya), respectively; (3) the individuals from Mt. Taihangshan (TH) are closely related to individuals from Mt. Dabashan (DB) in the autosomal tree, rather than individuals from Mt. Huangshan (HS) as indicated by the mitogenome tree, which supports the hypothesis that the ancestral rhesus macaques radiated into Mt. Taihangshan from Mt. Huangshan via Mt. Dabashan; and (4) the demographic scenario of the four lineages showed the ancestral rhesus macaques bottleneck and expansion corresponding to the suitable habitat reduction and expansion, which confirmed they had experienced northward recolonization and southward retreat events from Mt. Huangshan area via Northern China Plain to Northernmost China along with Pleistocene glacial cycles. This study provides a new insight into understanding how Pleistocene glaciation has influenced faunal diversity in subtropical and temperate China, especially for those exhibiting differential patterns of sex dispersal.

Towards a species-level phylogeny for Neotropical Myrtaceae: Notes on topology and resources for future studies

PREMISE

Increasingly complete phylogenies underpin studies in systematics, ecology, and evolution. Myrteae (Myrtaceae), with ~2700 species, is a key component of the exceptionally diverse Neotropical flora, but given its complicated taxonomy, automated assembling of molecular supermatrices from public databases often lead to unreliable topologies due to poor species identification.

METHODS

Here, we build a taxonomically verified molecular supermatrix of Neotropical Myrteae by assembling 3909 published and 1004 unpublished sequences from two nuclear and seven plastid molecular markers. We infer a time-calibrated phylogenetic tree that covers 712 species of Myrteae (~28% of the total diversity in the clade) and evaluate geographic and taxonomic gaps in sampling.

RESULTS

The tree inferred from the fully concatenated matrix mostly reflects the topology of the plastid data set and there is a moderate to strong incongruence between trees inferred from nuclear and plastid partitions. Large, species-rich genera are still the poorest sampled within the group. Eastern South America is the best-represented area in proportion to its species diversity, while Western Amazon, Mesoamerica, and the Caribbean are the least represented.

CONCLUSIONS

We provide a time-calibrated tree that can be more reliably used to address finer-scale eco-evolutionary questions that involve this group in the Neotropics. Gaps to be filled by future studies include improving representation of taxa and areas that remain poorly sampled, investigating causes of conflict between nuclear and plastid partitions, and the role of hybridization and incomplete lineage sorting in relationships that are poorly supported.

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.

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

Abundant incongruence in a clade endemic to a biodiversity hotspot: Phylogenetics of the scrub mint clade (Lamiaceae)

The Scrub Mint clade(Lamiaceae) provides a unique system for investigating the evolutionary processes driving diversification in the North American Coastal Plain from both a systematic and biogeographic context. The clade comprisesDicerandra, Conradina, Piloblephis, Stachydeoma, and four species of the broadly defined genus Clinopodium(Mentheae; Lamiaceae), almost all of which are endemic to the North American Eastern Coastal Plain. Most species of this clade are threatened or endangered and restricted to sandhill or a mosaic of scrub habitats. We analyzed relationships in this clade to understand the evolution of the group and identify evolutionary mechanisms acting on the clade, with important implications for conservation. We used a target-capture method to sequence and analyze 238 nuclear loci across all species of scrub mints, reconstructed the phylogeny, and calculated gene tree concordance, gene tree estimation error, and reticulation indices for every node in the tree using ML methods. Phylogenetic networks were used to determine reticulation events. Our nuclear phylogenetic estimates were consistent with previous results, while greatly increasing the robustness of taxon sampling. The phylogeny resolved the full relationship between Dicerandra and Conradina and the less-studied members of the clade (Piloblephis, Stachydeoma, Clinopodium spp.). We found hotspots of gene tree discordance and reticulation throughout the tree, especially in perennial Dicerandra. Several instances of reticulation events were uncovered between annual and perennial Dicerandra, and within the Conradina + allies clade. Incomplete lineage sorting also likely contributed to phylogenetic discordance. These results clarify phylogenetic relationships in the clade and provide insight on important evolutionary drivers in the clade, such as hybridization. General relationships in the group were confirmed, while the large amount of gene tree discordance is likely due to reticulation across the phylogeny.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

A phylogenomic perspective on interspecific competition

Evolutionary processes may have substantial impacts on community assembly, but evidence for phylogenetic relatedness as a determinant of interspecific interaction strength remains mixed. In this perspective, we consider a possible role for discordance between gene trees and species trees in the interpretation of phylogenetic signal in studies of community ecology. Modern genomic data show that the evolutionary histories of many taxa are better described by a patchwork of histories that vary along the genome rather than a single species tree. If a subset of genomic loci harbour trait-related genetic variation, then the phylogeny at these loci may be more informative of interspecific trait differences than the genome background. We develop a simple method to detect loci harbouring phylogenetic signal and demonstrate its application through a proof-of-principle analysis of Penicillium genomes and pairwise interaction strength. Our results show that phylogenetic signal that may be masked genome-wide could be detectable using phylogenomic techniques and may provide a window into the genetic basis for interspecific interactions.

Reticulate evolution of the tertiary relict Osmanthus

When interspecific gene flow is common, species relationships are more accurately represented by a phylogenetic network than by a bifurcating tree. This study aimed to uncover the role of introgression in the evolution of Osmanthus, the only genus of the subtribe Oleinae (Oleaceae) with its distribution center in East Asia. We built species trees, detected introgression, and constructed networks using multiple kinds of sequencing data (whole genome resequencing, transcriptome sequencing, and Sanger sequencing of nrDNA) combined with concatenation and coalescence approaches. Then, based on well-understood species relationships, historical biogeographic analyses and diversification rate estimates were employed to reveal the history of Osmanthus. Osmanthus originated in mid-Miocene Europe and dispersed to the eastern Tibetan Plateau in the late Miocene. Thereafter, it continued to spread eastwards. Phylogenetic conflict is common within the 'Core Osmanthus' clade and is seen at both early and late stages of diversification, leading to hypotheses of net-like species relationships. Incomplete lineage sorting proved ineffective in explaining phylogenetic conflicts and thus supported introgression as the main cause of conflicts. This study elucidates the diversification history of a relict genus in the subtropical regions of eastern Asia and reveals that introgression had profound effects on its evolutionary history.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Multiple paths toward repeated phenotypic evolution in the spiny-leg adaptive radiation (Tetragnatha; Hawai'i)

The repeated evolution of phenotypes provides clear evidence for the role of natural selection in driving evolutionary change. However, the evolutionary origin of repeated phenotypes can be difficult to disentangle as it can arise from a combination of factors such as gene flow, shared ancestral polymorphisms or mutation. Here, we investigate the presence of these evolutionary processes in the Hawaiian spiny-leg Tetragnatha adaptive radiation, which includes four microhabitat-specialists or ecomorphs, with different body pigmentation and size (Green, Large Brown, Maroon, and Small Brown). We investigated the evolutionary history of this radiation using 76 newly generated low-coverage, whole-genome resequenced samples, along with phylogenetic and population genomic tools. Considering the Green ecomorph as the ancestral state, our results suggest that the Green ecomorph likely re-evolved once, the Large Brown and Maroon ecomorphs evolved twice and the Small Brown evolved three times. We found that the evolution of the Maroon and Small Brown ecomorphs likely involved ancestral hybridization events, while the Green and Large Brown ecomorphs likely evolved through novel mutations, despite a high rate of incomplete lineage sorting in the dataset. Our findings demonstrate that the repeated evolution of ecomorphs in the Hawaiian spiny-leg Tetragnatha is influenced by multiple evolutionary processes.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Diversification and introgression in four chromosomal taxa of the Pearson's horseshoe bat (Rhinolophus pearsoni) group

Chromosomal variation among closely related taxa is common in both plants and animals, and can reduce rates of introgression as well as promote reproductive isolation and speciation. In mammals, studies relating introgression to chromosomal variation have tended to focus on a few model systems and typically characterized levels of introgression using small numbers of loci. Here we took a genome-wide approach to examine how introgression rates vary among four closely related horseshoe bats (Rhinolophus pearsoni group) that possess different diploid chromosome numbers (2n = 42, 44, 46, and 60) resulting from Robertsonian (Rb) changes (fissions/fusions). Using a sequence capture we obtained orthologous loci for thousands of nuclear loci, as well as mitogenomes, and performed phylogenetic and population genetic analyses. We found that the taxon with 2n = 60 was the first to diverge in this group, and that the relationships among the three other taxa (2n = 42, 44 and 46) showed discordance across our different analyses. Our results revealed signatures of multiple ancient introgression events between the four taxa, with evidence of mitonuclar discordance in phylogenetic trees and reticulation events in their evolutionary history. Despite this, we found no evidence of recent and/or ongoing introgression between taxa. Overall, our results indicate that the effects of Rb changes on the reduction of introgression are complicated and that these may contribute to reproductive isolation and speciation in concert with other factors (e.g. phenotypic and genic divergence).

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

De Novo Transcriptome Assembly and Comparative Analysis of Differentially Expressed Genes Involved in Cold Acclimation and Freezing Tolerance of the Arctic Moss Aulacomnium turgidum (Wahlenb.) Schwaegr

Cold acclimation refers to a phenomenon in which plants become more tolerant to freezing after exposure to non-lethal low temperatures. Aulacomnium turgidum (Wahlenb.) Schwaegr is a moss found in the Arctic that can be used to study the freezing tolerance of bryophytes. To improve our understanding of the cold acclimation effect on the freezing tolerance of A. turgidum, we compared the electrolyte leakage of protonema grown at 25 °C (non-acclimation; NA) and at 4 °C (cold acclimation; CA). Freezing damage was significantly lower in CA plants frozen at -12 °C (CA-12) than in NA plants frozen at -12 °C (NA-12). During recovery at 25 °C, CA-12 demonstrated a more rapid and greater level of the maximum photochemical efficiency of photosystem II than NA-12, indicating a greater recovery capacity for CA-12 compared to NA-12. For the comparative analysis of the transcriptome between NA-12 and CA-12, six cDNA libraries were constructed in triplicate, and RNA-seq reads were assembled into 45,796 unigenes. The differential gene expression analysis showed that a significant number of AP2 transcription factor genes and pentatricopeptide repeat protein-coding genes related to abiotic stress and the sugar metabolism pathway were upregulated in CA-12. Furthermore, starch and maltose concentrations increased in CA-12, suggesting that cold acclimation increases freezing tolerance and protects photosynthetic efficiency through the accumulation of starch and maltose in A. turgidum. A de novo assembled transcriptome can be used to explore genetic sources in non-model organisms.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

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.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

Phylogenomic analyses based on the plastid genome and concatenated nrDNA sequence data reveal cytonuclear discordance in genus Atractylodes (Asteraceae: Carduoideae)

Atractylodes species are widely distributed across East Asia and are cultivated as medicinal herbs in China, Japan, and Korea. Their unclear morphological characteristics and low levels of genetic divergence obscure the taxonomic relationships among these species. In this study, 24 plant samples were collected representing five species of Atractylodes located in China; of these, 23 belonged to members of the A. lancea complex. High-throughput sequencing was used to obtain the concatenated nrDNA sequences (18S-ITS1-5.8S-ITS2-28S) and plastid genomes. The concatenated nrDNA sequence lengths for all the Atractylodes species were 5,849 bp, and the GC content was 55%. The lengths of the whole plastid genome sequences ranged from 152,138 bp (A. chinensis) to 153,268 bp (A. lancea), while their insertion/deletion sites were mainly distributed in the intergenic regions. Furthermore, 33, 34, 36, 31, and 32 tandem repeat sequences, as well as 30, 30, 29, 30, and 30 SSR loci, were detected in A. chinensis, A. koreana, A. lancea, A. japonica, and A. macrocephala, respectively. In addition to these findings, a considerable number of heteroplasmic variations were detected in the plastid genomes, implying a complicated phylogenetic history for Atractylodes. The results of the phylogenetic analysis involving concatenated nrDNA sequences showed that A. lancea and A. japonica formed two separate clades, with A. chinensis and A. koreana constituting their sister clade, while A. lancea, A. koreana, A. chinensis, and A. japonica were found based on plastid datasets to represent a mixed clade on the phylogenetic tree. Phylogenetic network analysis suggested that A. lancea may have hybridized with the common ancestor of A. chinensis and A. japonica, while ABBA-BABA tests of SNPs in the plastid genomes showed that A. chinensis was more closely related to A. japonica than to A. lancea. This study reveals the extensive discordance and complexity of the relationships across the members of the A. lancea complex (A. lancea, A. chinensis, A. koreana, and A. japonica) according to cytonuclear genomic data; this may be caused by interspecific hybridization or gene introgression.

Integrating morphological and genetic limits in the taxonomic delimitation of the Cuban taxa of Magnoliasubsect.Talauma (Magnoliaceae)

An accurate taxa delimitation, based on a full understanding of evolutionary processes involved in taxa differentiation, can be gained from a combination of ecological, morphological, and molecular approaches. The taxonomy of Magnoliasubsect.Talauma in Cuba has long been debated and exclusively based on traditional morphological study of a limited number of individuals. A more accurate description of leaf morphology variation using geometric morphometrics combined with genetic data could bring consistency to taxa delimitation in this group. Leaf samples for the morphological (243) and genetic (461) analyses were collected throughout the entire distribution range. The variability of each taxon was analyzed through multivariate and geometric morphometry, and 21 genetic markers (SSR). The observed leaf morphological variability was higher than previously described. Morphological and genetic classifications were highly congruent in two out of four taxa. Our data brought evidence that Magnoliaorbiculata can be considered a true species with very clear genetic and morphological limits. The main taxonomic issues concern the north-eastern Cuban populations of Magnoliasubsect.Talauma. The data supported the existence of two clear groups: corresponding mainly to M.minor-M.oblongifolia and T.ophiticola. However, these two groups cannot be considered fully delimited since genetic markers provided evidence of genetic admixture between them. Due to the likely absence of, at least strong, reproductive barriers between these three taxa, we propose therefore to consider them as a species complex.

Demographic history and gene flow in the peatmosses Sphagnum recurvum and Sphagnum flexuosum (Bryophyta: Sphagnaceae)

Population size changes and gene flow are processes that can have significant impacts on evolution. The aim of this study was to investigate the relationship of geography to patterns of gene flow and population size changes in a pair of closely related Sphagnum (peatmoss) species: S. recurvum and S. flexuosum. Both species occur in eastern North America, and S. flexuosum also occurs in Europe. Genetic data from restriction-site-associated DNA sequencing (RAD-seq) were used in this study. Analyses of gene flow were accomplished using coalescent simulations of site frequency spectra (SFSs). Signatures of gene flow were confirmed by f 4 statistics. For S. flexuosum, genetic diversity of plants in glaciated areas appeared to be lower than that in unglaciated areas, suggesting that glaciation can have an impact on effective population sizes. There is asymmetric gene flow from eastern North America to Europe, suggesting that Europe might have been colonized by plants from eastern North America after the last glacial maximum. The rate of gene flow between S. flexuosum and S. recurvum is lower than that between geographically disjunct S. flexuosum populations. The rate of gene flow between species is higher among sympatric plants of the two species than between currently allopatric S. flexuosum populations. There was also gene flow from S. recurvum to the ancestor S. flexuosum on both continents which occurred through secondary contact. These results illustrate a complex history of interspecific gene flow between S. flexuosum and S. recurvum, which occurred in at least two phases: between ancestral populations after secondary contact and between currently sympatric plants.

OrthoSNAP: A tree splitting and pruning algorithm for retrieving single-copy orthologs from gene family trees

Molecular evolution studies, such as phylogenomic studies and genome-wide surveys of selection, often rely on gene families of single-copy orthologs (SC-OGs). Large gene families with multiple homologs in 1 or more species-a phenomenon observed among several important families of genes such as transporters and transcription factors-are often ignored because identifying and retrieving SC-OGs nested within them is challenging. To address this issue and increase the number of markers used in molecular evolution studies, we developed OrthoSNAP, a software that uses a phylogenetic framework to simultaneously split gene families into SC-OGs and prune species-specific inparalogs. We term SC-OGs identified by OrthoSNAP as SNAP-OGs because they are identified using a splitting and pruning procedure analogous to snapping branches on a tree. From 415,129 orthologous groups of genes inferred across 7 eukaryotic phylogenomic datasets, we identified 9,821 SC-OGs; using OrthoSNAP on the remaining 405,308 orthologous groups of genes, we identified an additional 10,704 SNAP-OGs. Comparison of SNAP-OGs and SC-OGs revealed that their phylogenetic information content was similar, even in complex datasets that contain a whole-genome duplication, complex patterns of duplication and loss, transcriptome data where each gene typically has multiple transcripts, and contentious branches in the tree of life. OrthoSNAP is useful for increasing the number of markers used in molecular evolution data matrices, a critical step for robustly inferring and exploring the tree of life.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.