Complex Hybridization in a Clade of Polytypic Salamanders (Plethodontidae: Desmognathus) Uncovered by Estimating Higher-Level Phylogenetic Networks.

Reticulation between incipient lineages is a common feature of diversification. We examine these phenomena in the Pisgah clade of Desmognathus salamanders from the southern Appalachian Mountains of the eastern United States. The group contains four to seven species exhibiting two discrete phenotypes, aquatic "shovel-nosed" and semi-aquatic "black-bellied" forms. These ecomorphologies are ancient and have apparently been transmitted repeatedly between lineages through introgression. Geographically proximate populations of both phenotypes exhibit admixture, and at least two black-bellied lineages have been produced via reticulations between shovel-nosed parentals, suggesting potential hybrid speciation dynamics. However, computational constraints currently limit our ability to reconstruct network radiations from gene-tree data. Available methods are limited to level-1 networks wherein reticulations do not share edges, and higher-level networks may be non-identifiable in many cases. We present a heuristic approach to recover information from higher-level networks across a range of potentially identifiable empirical scenarios, supported by theory and simulation. When extrinsic information indicates the location and direction of reticulations, our method can successfully estimate a reduced possible set of non-level-1 networks. Phylogenomic data support a single backbone topology with up to five overlapping hybrid edges in the Pisgah clade. These results suggest an unusual mechanism of ecomorphological hybrid speciation, wherein a binary threshold trait causes some hybrid populations to shift between microhabitat niches, promoting ecological divergence between sympatric hybrids and parentals. This contrasts with other well-known systems in which hybrids exhibit intermediate, novel, or transgressive phenotypes. The genetic basis of these phenotypes is unclear and further data are needed to clarify the evolutionary basis of morphological changes with ecological consequences.

Genomic changes and stabilization following homoploid hybrid speciation of the Oxford ragwort Senecio squalidus.

Oxford ragwort (Senecio squalidus) is one of only two homoploid hybrid species known to have originated very recently, so it is a unique model for determining genomic changes and stabilization following homoploid hybrid speciation. Here, we provide a chromosome-level genome assembly of S. squalidus with 95% of the assembly contained in the 10 longest scaffolds, corresponding to its haploid chromosome number. We annotated 30,249 protein-coding genes and estimated that ∼62% of the genome consists of repetitive elements. We then characterized genome-wide patterns of linkage disequilibrium, polymorphism, and divergence in S. squalidus and its two parental species, finding that (1) linkage disequilibrium is highly heterogeneous, with a region on chromosome 4 showing increased values across all three species but especially in S. squalidus; (2) regions harboring genetic incompatibilities between the two parental species tend to be large, show reduced recombination, and have lower polymorphism in S. squalidus; (3) the two parental species have an unequal contribution (70:30) to the genome of S. squalidus, with long blocks of parent-specific ancestry supporting a very rapid stabilization of the hybrid lineage after hybrid formation; and (4) genomic regions with major parent ancestry exhibit an overrepresentation of loci with evidence for divergent selection occurring between the two parental species on Mount Etna. Our results show that both genetic incompatibilities and natural selection play a role in determining genome-wide reorganization following hybrid speciation and that patterns associated with homoploid hybrid speciation-typically seen in much older systems-can evolve very quickly following hybridization.

Phylogeographic analysis reveals multiple origins of the desert shrub Reaumuria songarica in northern Xinjiang, involving homoploid and tetraploid hybrids.

Hybrid speciation plays an important role in species diversification. The establishment of reproductive isolation is crucial for hybrid speciation, and the identification of diverse types of hybrids, particularly homoploid hybrid species, contributes to a comprehensive understanding of this process. Reaumuria songarica is a constructive shrub widespread in arid Central Asia. Previous studies have inferred that the R. songarica populations in the Gurbantunggut Desert (GuD) originated from homoploid hybridizations between its eastern and western lineages and may have evolved into an incipient species. To further elucidate the genetic composition of different hybrid populations and to determine the species boundary of this hybrid lineage, we investigated the overall phylogeographic structure of R. songarica based on variation patterns of five cpDNA and one nrITS sequences across 32 populations. Phylogenetic analyses demonstrated that within the GuD lineage, the Wuerhe population evolved directly from ancestral lineages, whereas the others originated from hybridizations between the eastern and western lineages. PCoA and genetic barrier analysis supported the subdivision of the GuD lineage into the southern (GuD-S) and northern (GuD-N) groups. Populations in the GuD-S group had a consistent genetic composition and the same ancestral female parent, indicating that they belonged to a homoploid hybrid lineage. However, the GuD-N group experienced genetic admixture of the eastern and western lineages on nrITS and cpDNA, with some populations inferred to be allopolyploid based on ploidy data. Based on cpDNA haplotypes, BEAST analyses showed that the GuD-S and GuD-N groups originated after 0.5 Ma. Our results suggest that multiple expansions and contractions of GuD, driven by Quaternary climatic oscillations and the Kunlun-Yellow River tectonic movement, are important causes of the complex origins of R. songarica populations in northern Xinjiang. This study highlights the complex origins of the Junggar Basin flora and the underappreciated role of hybridization in increasing its species diversity.

The role of deep hybridization in fern speciation: Examples from the Thelypteridaceae.

PREMISE: Hybridization is recognized as an important mechanism in fern speciation, with many allopolyploids known among congeners, as well as evidence of ancient genome duplications. Several contemporary instances of deep (intergeneric) hybridization have been noted, invariably resulting in sterile progeny. We chose the christelloid lineage of the family Thelypteridaceae, recognized for its high frequency of both intra- and intergeneric hybrids, to investigate recent hybrid speciation between deeply diverged lineages. We also seek to understand the ecological and evolutionary outcomes of resulting lineages across the landscape. METHODS: By phasing captured reads within a phylogenomic data set of GoFlag 408 nuclear loci using HybPhaser, we investigated candidate hybrids to identify parental lineages. We estimated divergence ages by inferring a dated phylogeny using fossil calibrations with treePL. We investigated ecological niche conservatism between one confirmed intergeneric allotetraploid and its diploid progenitors using the centroid, overlap, unfilling, and expansion (COUE) framework. RESULTS: We provide evidence for at least six instances of intergeneric hybrid speciation within the christelloid clade and estimate up to 45 million years of divergence between progenitors. The niche quantification analysis showed moderate niche overlap between an allopolyploid species and its progenitors, with significant divergence from the niche of one progenitor and conservatism to the other. CONCLUSIONS: The examples provided here highlight the overlooked role that allopolyploidization following intergeneric hybridization may play in fern diversification and range and niche expansions. Applying this approach to other fern taxa may reveal a similar pattern of deep hybridization resulting in highly successful novel lineages.

Topology Testing and Demographic Modeling Illuminate a Novel Speciation Pathway in the Greater Caribbean Sea Following the Formation of the Isthmus of Panama.

Recent genomic analyses have highlighted the prevalence of speciation with gene flow in many taxa and have underscored the importance of accounting for these reticulate evolutionary processes when constructing species trees and generating parameter estimates. This is especially important for deepening our understanding of speciation in the sea where fast-moving ocean currents, expanses of deep water, and periodic episodes of sea level rise and fall act as soft and temporary allopatric barriers that facilitate both divergence and secondary contact. Under these conditions, gene flow is not expected to cease completely while contemporary distributions are expected to differ from historical ones. Here, we conduct range-wide sampling for Pederson's cleaner shrimp (Ancylomenes pedersoni), a species complex from the Greater Caribbean that contains three clearly delimited mitochondrial lineages with both allopatric and sympatric distributions. Using mtDNA barcodes and a genomic ddRADseq approach, we combine classic phylogenetic analyses with extensive topology testing and demographic modeling (10 site frequency replicates × 45 evolutionary models × 50 model simulations/replicate = 22,500 simulations) to test species boundaries and reconstruct the evolutionary history of what was expected to be a simple case study. Instead, our results indicate a history of allopatric divergence, secondary contact, introgression, and endemic hybrid speciation that we hypothesize was driven by the final closure of the Isthmus of Panama and the strengthening of the Gulf Stream Current ~3.5 Ma. The history of this species complex recovered by model-based methods that allow reticulation differs from that recovered by standard phylogenetic analyses and is unexpected given contemporary distributions. The geologically and biologically meaningful insights gained by our model selection analyses illuminate what is likely a novel pathway of species formation not previously documented that resulted from one of the most biogeographically significant events in Earth's history.

Documenting homoploid hybrid speciation.

Homoploid hybrid speciation is challenging to document because hybridization can lead to outcomes other than speciation. Thus, some authors have argued that establishment of homoploid hybrid speciation should include evidence that reproductive barriers isolating the hybrid neo-species from its parental species were derived from hybridization. While this criterion is difficult to satisfy, several recent papers have successfully employed a common pipeline to identify candidate genes underlying such barriers and (in one case) to validate their function. We describe this pipeline, its application to several plant and animal species and what we have learned about homoploid hybrid speciation as a consequence. We argue that - given the ubiquity of admixture and the polygenic basis of reproductive isolation - homoploid hybrid speciation could be much more common and more protracted than suggested by earlier conceptual arguments and theoretical studies.

A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau.

While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.

Niche and phenotypic differentiation in fern hybrid speciation, a case study of Pteris fauriei (Pteridaceae).

BACKGROUND AND AIMS: Niche differentiation is a crucial issue in speciation. Although it has a well-known role in adaptive processes of hybrid angiosperms, it is less understood in hybrid ferns. Here, we investigate whether an intermediate ecological niche of a fern hybrid is a novel adaptation that provides insights into fern hybrid speciation. METHODS: Pteris fauriei (Pteridaceae) is a natural hybrid fern, occurring in environments between its parent species. The maternal Pteris minor is found in sunny areas, but the habitat of the paternal Pteris latipinna is shady. We combined data from morphology, leaf anatomy and photosynthetic traits to explore adaptation and differentiation, along with measuring the environmental features of their niches. We also performed experiments in a common garden to understand ecological plasticity. KEY RESULTS: The hybrid P. fauriei was intermediate between the parent species in stomatal density, leaf anatomical features and photosynthetic characteristics in both natural habitats and a common garden. Interestingly, the maternal P. minor showed significant environmental plasticity and was more similar to the hybrid P. fauriei in the common garden, suggesting that the maternal species experiences stress in its natural habitats but thrives in environments similar to those of the hybrid. CONCLUSIONS: Based on the similar niche preferences of the hybrid and parents, we propose hybrid superiority. Our results indicate that the hybrid P. fauriei exhibits greater fitness and can compete with and occupy the initial niches of the maternal P. minor. Consequently, we suggest that the maternal P. minor has experienced a niche shift, elucidating the pattern of niche differentiation in this hybrid group. These findings offer a potential explanation for the frequent occurrence of hybridization in ferns and provide new insights into fern hybrid speciation, enhancing our understanding of fern diversity.

A model of hybrid speciation process drawn from three new poplar species originating from distant hybridization between sections.

Although numerous studies have been conducted on hybrid speciation, our understanding of this process remains limited. Through an 18-year systematic investigation of all taxa of Populus on the Qinghai-Tibet Plateau, we discovered three new taxa with clear characteristics of sect. Leucoides. Further evidence was gathered from morphology, whole-genome bioinformatics, biogeography, and breeding to demonstrate synthetically that they all originated from distant hybridization between sect. Leucoides and sect. Tacamahaca. P. gonggaensis originated from the hybridization of P. lasiocarpa with P. cathayana, P. butuoensis from the hybridization of P. wilsonii with P. szechuanica, and P. dafengensis from the hybridization of P. lasiocarpa with P. szechuanica. Due to heterosis, the three hybrid taxa possess greater ecological adaptability than their ancestral species. We propose a hybrid speciation process model that incorporates orthogonal, reverse, and backcrossing events. This model can adequately explain some crucial evolutionary concerns, such as the nuclear-cytoplasmic conflict on phylogeny and the extinction of ancestral species within the distribution range of hybrid species.

In search of the Goldilocks zone for hybrid speciation II: hard times for hybrid speciation?

Hybridization opens a unique window for observing speciation mechanisms and is a potential engine of speciation. One controversially discussed outcome of hybridization is homoploid hybrid speciation by reciprocal sorting, where a hybrid population maintains a mixed combination of the parental genetic incompatibilities, preventing further gene exchange between the newly formed population and the two parental sources. Previous work showed that, for specific linkage architectures (i.e., the genomic location and order of hybrid incompatibilities), reciprocal sorting could reliably result in hybrid speciation. Yet, the sorting of incompatibilities creates a risk of population extinction. To understand how the demographic consequences of the purging of incompatibilities interact with the formation of a hybrid species, we model an isolated hybrid population resulting from a single admixture event. We study how population size, linkage architecture, and the strength of the incompatibility affect survival of the hybrid population, resolution/purging of the genetic incompatibilities and the probability of observing hybrid speciation. We demonstrate that the extinction risk is highest for intermediately strong hybrid incompatibilities. In addition, the linkage architecture displaying the highest hybrid speciation probabilities changes drastically with population size. Overall, this indicates that population dynamics can strongly affect the outcome of hybridization and the hybrid speciation probability.

Comparative transcriptomics and gene expression divergence associated with homoploid hybrid speciation in Argyranthemum.

Ecological isolation is increasingly thought to play an important role in speciation, especially for the origin and reproductive isolation of homoploid hybrid species. However, the extent to which divergent and/or transgressive gene expression changes are involved in speciation is not well studied. In this study, we employ comparative transcriptomics to investigate gene expression changes associated with the origin and evolution of two homoploid hybrid plant species, Argyranthemum sundingii and A. lemsii (Asteraceae). As there is no standard methodology for comparative transcriptomics, we examined five different pipelines for data assembly and analysing gene expression across the four species (two hybrid and two parental). We note biases and problems with all pipelines, and the approach used affected the biological interpretation of the data. Using the approach that we found to be optimal, we identify transcripts showing DE between the parental taxa and between the homoploid hybrid species and their parents; in several cases, putative functions of these DE transcripts have a plausible role in ecological adaptation and could be the cause or consequence of ecological speciation. Although independently derived, the homoploid hybrid species have converged on similar expression phenotypes, likely due to adaptation to similar habitats.

Post-hybridization introgression and natural selection promoted genomic divergence of Aegilops speltoides and the four S*-genome diploid species.

Understanding how different driving forces have promoted biological divergence and speciation is one of the central issues in evolutionary biology. The Triticum/Aegilops species complex contains 13 diploid species belonging to the A-, B- and D-lineages and offers an ideal system to address the evolutionary dynamics of lineage fusion and splitting. Here, we sequenced the whole genomes of one S-genome species (Aegilops speltoides) of the B-lineage and four S*-genome diploid species (Aegilops bicornis, Aegilops longissima, Aegilops sharonensis and Aegilops searsii) of the D-lineage at the population level. We performed detailed comparisons of the five species and with the other four representative A-, B- and D-lineage species. Our estimates identified frequent genetic introgressions from A- and B-lineages to the D-lineage species. A remarkable observation is the contrasting distributions of putative introgressed loci by the A- and B-lineages along all the seven chromosomes to the extant D-lineage species. These genetic introgressions resulted in high levels of genetic divergence at centromeric regions between Ae. speltoides (B-lineage) and the other four S*-genome diploid species (D-lineage), while natural selection is a potential contributor to divergence among the four S*-genome species at telomeric regions. Our study provides a genome-wide view on how genetic introgression and natural selection acted together yet chromosome-regionally divided to promote genomic divergence among the five S- and S*-genome diploid species, which provides new and nuanced insights into the evolutionary history of the Triticum/Aegilops species complex.

Independent origins of Spiranthes×kapnosperia (Orchidaceae) and their nomenclatural implications.

Spiranthes Rich. (Orchidaceae) is a commonly encountered but systematically and nomenclaturally challenging component of the North American orchid flora. Here, the evolutionary history and hybrid origin of the recently described S.sheviakii Hough and Young are critically examined. The available molecular data unambiguously support a hybrid origin of S.cernua (L.) Rich. × S.ochroleuca (Rydb.) Rydb. for S.sheviakii, the same parentage as the priority name S.×kapnosperia M.C. Pace. As hybrid formulas can have only one correct name, S.sheviakii is a synonym of S.×kapnosperia. It is likely that S.×kapnosperia evolved independently at least twice in at least two widely disjunct locations.

Transgressive segregation in mating traits drives hybrid speciation.

Hybridization can rapidly generate novel genetic variation, which can promote ecological speciation by creating novel adaptive phenotypes. However, it remains unclear how hybridization, creating novel mating phenotypes (e.g., mating season, genitalia shapes, sexual displays, mate preferences), affects speciation especially when the phenotypes do not confer adaptive advantages. Here, based on individual-based evolutionary simulations, we propose that transgressive segregation of mating traits can drive incipient hybrid speciation. Simulations demonstrated that incipient hybrid speciation occurred most frequently when the hybrid population received moderate continued immigration from parental lineages causing recurrent episodes of hybridization. Recurrent hybridization constantly generated genetic variation, which promoted the rapid stochastic evolution of mating phenotypes in a hybrid population. The stochastic evolution continued until a novel mating phenotype came to dominate the hybrid population, which reproductively isolates the hybrid population from parental lineages. However, too frequent hybridization rather hindered the evolution of reproductive isolation by inflating the variation of mating phenotypes to produce phenotypes allowing mating with parental lineages. Simulations also revealed conditions for the long-term persistence of hybrid species after their incipient emergence. Our results suggest that recurrent transgressive segregation of mating phenotypes can offer a plausible explanation for hybrid speciation and radiations that involved little adaptive ecological divergence.

Genomics of plant speciation.

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Multispecies polyploidization, chromosome shuffling, and genome extraction in Zea/Tripsacum hybrids.

By hybridization and special sexual reproduction, we sequentially aggregated Zea mays, Zea perennis, and Tripsacum dactyloides in an allohexaploid, backcrossed it with maize, derived self-fertile allotetraploids of maize and Z. perennis by natural genome extraction, extended their first six selfed generations, and finally constructed amphitetraploid maize using nascent allotetraploids as a genetic bridge. Transgenerational chromosome inheritance, subgenome stability, chromosome pairings and rearrangements, and their impacts on an organism's fitness were investigated by fertility phenotyping and molecular cytogenetic techniques genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH). Results showed that diversified sexual reproductive methods produced highly differentiated progenies (2n = 35-84) with varying proportions of subgenomic chromosomes, of which one individual (2n = 54, MMMPT) overcame self-incompatibility barriers and produced a self-fertile nascent near-allotetraploid by preferentially eliminating Tripsacum chromosomes. Nascent near-allotetraploid progenies showed persistent chromosome changes, intergenomic translocations, and rDNA variations for at least up to the first six selfed generations; however, the mean chromosome number preferably maintained at the near-tetraploid level (2n = 40) with full integrity of 45S rDNA pairs, and a trend of decreasing variations by advancing generations with an average of 25.53, 14.14, and 0.37 for maize, Z. perennis, and T. dactyloides chromosomes, respectively. The mechanisms for three genome stabilities and karyotype evolution for formatting new polyploid species were discussed.

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.

Hybrid zone or hybrid lineage: a genomic reevaluation of Sibley's classic species conundrum in Pipilo towhees.

Hybrid zones can be studied by modeling clines of trait variation (e.g., morphology, genetics) over a linear transect. Yet, hybrid zones can also be spatially complex, can shift over time, and can even lead to the formation of hybrid lineages with the right combination of dispersal and vicariance. We reassessed Sibley's (1950) gradient between Collared Towhee (Pipilo ocai) and Spotted Towhee (Pipilo maculatus) in Central Mexico to test whether it conformed to a typical tension-zone cline model. By comparing historical and modern data, we found that cline centers for genetic and phenotypic traits have not shifted over the course of 70 years. This equilibrium suggests that secondary contact between these species, which originally diverged over 2 million years ago, likely dates to the Pleistocene. Given the amount of mtDNA divergence, parental ends of the cline have very low autosomal nuclear differentiation (FST = 0.12). Dramatic and coincident cline shifts in mtDNA and throat color suggest the possibility of sexual selection as a factor in differential introgression, while a contrasting cline shift in green back color hints at a role for natural selection. Supporting the idea of a continuum between clinal variation and hybrid lineage formation, the towhee gradient can be analyzed as one population under isolation-by-distance, as a two-population cline, and as three lineages experiencing divergence with gene flow. In the middle of the gradient, a hybrid lineage has become partly isolated, likely due to forested habitat shrinking and fragmenting as it moved upslope after the last glacial maximum and a stark environmental transition. This towhee system offers a window into the potential outcomes of hybridization across a dynamic landscape including the creation of novel genomic and phenotypic combinations and incipient hybrid lineages.

Multiple concordant cytonuclear divergences and potential hybrid speciation within a species complex in Asia.

Various environmental factors impact the distribution, population structure, demography and evolutionary trajectory of a bird species, leading to genetic and morphological divergences between populations across its distribution. The Paradoxornis webbianus species complex is found throughout much of East Asia, where its geographically distinct populations exhibit dramatic morphological variation. This has resulted in a hotly debated taxonomy. This study intended to identify genetic divergence patterns and their underlying contributing factors for this species complex. We collected 243 birds, whose data was combined with those available in GenBank to perform phylogeographic analyses using one mitochondrial and six nuclear loci. Six mitochondrial clades were observed in the species complex, while individual-based Bayesian clustering using nuclear markers showed multiple congruent breaks. Overall, the six molecular lineages could be recognized as independent species under the lineage species concept in view of genetic divergence, clade-specific morphological changes and distribution: P. webbianus, P. w. bulomachus, P. alphonsianus, P. a. ganluoensis, P. brunneus brunneus and P. b. ricketti. The estimated divergence times range from 0.46 to 3.36 million years ago, suggesting it was likely impacted by paleoclimatic changes. Interestingly, P. alphonsianus carries two divergent mitochondrial lineages shared with P. webbianus and P. a. ganluoensis, respectively, and analyses based on nuclear loci found a similar pattern. We discussed the various hypotheses for this pattern and argued that P. alphonsianus was likely the result of hybridization between P. webbianus and P. a. ganluoensis. Further data on genome, transcriptome and breeding ecology are needed to address the hypothesis of hybrid speciation and its underlying mechanisms.

Sequential hybridization may have facilitated ecological transitions in the Southwestern pinyon pine syngameon.

Multispecies interbreeding networks, or syngameons, have been increasingly reported in natural systems. However, the formation, structure, and maintenance of syngameons have received little attention. Through gene flow, syngameons can increase genetic diversity, facilitate the colonization of new environments, and contribute to hybrid speciation. In this study, we evaluated the history, patterns, and consequences of hybridization in a pinyon pine syngameon using morphological and genomic data to assess genetic structure, demographic history, and geographic and climatic data to determine niche differentiation. We demonstrated that Pinus edulis, a dominant species in the Southwestern US and a barometer of climate change, is a core participant in the syngameon, involved in the formation of two drought-adapted hybrid lineages including the parapatric and taxonomically controversial fallax-type. We found that species remain morphologically and genetically distinct at range cores, maintaining species boundaries while undergoing extensive gene flow in areas of sympatry at range peripheries. Our study shows that sequential hybridization may have caused relatively rapid speciation and facilitated the colonization of different niches, resulting in the rapid formation of two new lineages. Participation in the syngameon may allow adaptive traits to be introgressed across species barriers and provide the changes needed to survive future climate scenarios.