Altitudinal gradients, plant hybrid zones and evolutionary novelty

Hybrid zones in the European Alps impact the phylogeography of alpine vicariant willow species (Salix L.)

Introduction

In the European Alps, Pleistocene climate oscillations resulted in geographical range expansions and restrictions of species. Postglacial recolonizations often result in secondary contact hybridization of vicariant species, thereby creating hybrid zones with patterns of introgression. Here, we compare the genetic structure of two secondary contact hybrid zones between two vicariant willow species pairs occurring in the European Alpine System. Supplemented by morphological and ecological data, we try to understand the factors shaping the hybrid zones and their influence on geographical range filling patterns.

Methods

RAD sequencing and morphometric data were used to characterize biogeographical history, genetic diversity and the hybrid zone of each species pair. Vegetation relevés and species distribution models provided ecological context and support.

Key results

Results suggest that recolonization of the Alps happened from peripheral glacial refugia, resulting in broad secondary contact zones in the Eastern Alps in both species pairs. Both hybrid zones show introgression, but differ in symmetry and intensity of gene flow, in the type of introgressed loci, and in the geographical range. Habitat preferences and species distribution models do not indicate ecological barriers to recolonization.

Conclusions

Hybrid zones do not only affect the genetic structure of species by gene flow and introgression, but also appear to impact the biogeographical patterns of species.

Genetic Diversity and Ecogeographical Niche Overlap Among Hybridising Ox-Eye Daisies (Leucanthemum, Asteraceae) in the Carpathian Mountains: The Impact of Anthropogenic Disturbances

Climate change and human influence are transforming mountain ecosystems, significantly impacting species distributions and biodiversity. Among these changes, the upward migration of lowland species into mountain regions stands out. This study examines the ecogeographical niche overlap and genetic diversity among three Leucanthemum species distributed along an altitudinal gradient in the Carpathian Mountains: the lowland L. ircutianum (4x), the montane L. rotundifolium (2x) and the alpine L. gaudinii (2x). By genotyping over 600 individuals using SNP analysis, followed by principal coordinate analysis (PCoA), Neighbour-Net Network and Structure clustering, we reveal not just distinct genetic groups but also hybridisation across all species, suggesting the potential for triple hybrids. Genetic admixture is further supported by environmental background and niche overlap analyses that reveal substantial overlap among species, particularly in line with their vertical distribution. Climate envelope plots indicate a likely reduction in available habitat for mountainous species due to climate change, leading to an increase in competition and an intensification of hybridisation. Anthropogenic influences are further intensifying these hybridisation trends. Among the studied species, L. gaudinii is most at risk of overwhelming hybridisation, whereas L. ircutianum may experience habitat expansion. By providing a comprehensive genetic and ecological overview, our research highlights the significance of hybridisation in biodiversity conservation and the challenges posed by environmental changes and anthropogenic activities in mountain environments. This study not only contributes to the understanding of genetic diversity in the Carpathians but also underscores the broader implications for molecular ecology and conservation strategies in mountain ecosystems.

Genomic Introgression in the Hybrid zones at the Margins of the Species' Range Between Ecologically Distinct Rubus Species

Populations in extreme environments at the margins of a species' range are often the most vulnerable to climate change, but they may also experience novel evolutionary processes, such as secondary contact and hybridization with their relatives. The range overlap resulting from secondary contact with related species that have adapted to different climatic zones may act as corridors for adaptive introgression. To test this hypothesis, we examined the hybrid zones along the altitude of two closely related Rubus species, one temperate and the other subtropical species, at their southern and northern limits on Yakushima Island, Japan. Genomic cline analysis revealed non-neutral introgression throughout the genome in both directions in the two species. Some of these genomic regions involve gene ontology terms related to the regulation of several biological processes. Our niche modeling suggests that, assuming niche conservatism, the temperate species are likely to lose their suitable habitat, and the backcrossed hybrids with the subtropical species are already expanding upslope on the island. Adaptive introgression through the hybrid zone may contribute to the persistence and expansion of the species in the southernmost and northernmost populations.

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.

Weak Genetic Isolation and Putative Phenotypic Selection in the Wild Carnation Dianthus virgineus (Caryophyllaceae)

By relating genetic divergence at neutral loci, phenotypic variation, and geographic and environmental distances, it is possible to dissect micro-evolutionary scenarios involving natural selection and neutral evolution. In this work, we tested the patterns of intraspecific genetic and phenotypic variation along an elevational gradient, using Dianthus virgineus as study system. We genotyped genome-wide SNPs through ddRAD sequencing and quantified phenotypic variation through multivariate morphological variation. We assessed patterns of variation by testing the statistical association between genetic, phenotypic, geographic, and elevational distances and explored the role of genetic drift and selection by comparing the Fst and Pst of morphometric traits. We revealed a weak genetic structure related to geographic distance among populations, but we excluded the predominant role of genetic drift acting on phenotypic traits. A high degree of phenotypic differentiation with respect to genetic divergence at neutral loci allowed us to hypothesize the effect of selection, putatively fuelled by changing conditions at different sites, on morphological traits. Thus, natural selection acting despite low genetic divergence at neutral loci can be hypothesized as a putative driver explaining the observed patterns of variation.

Demographic responses of hybridizing cinquefoils to changing climate in the Colorado Rocky Mountains

Hybridization between taxa generates new pools of genetic variation that can lead to different environmental responses and demographic trajectories over time than seen in parental lineages. The potential for hybrids to have novel environmental tolerances may be increasingly important in mountainous regions, which are rapidly warming and drying due to climate change. Demographic analysis makes it possible to quantify within- and among-species responses to variation in climate and to predict population growth rates as those conditions change. We estimated vital rates and population growth in 13 natural populations of two cinquefoil taxa (Potentilla hippiana and P. pulcherrima) and their hybrid across elevation gradients in the Southern Rockies. Using three consecutive years of environmental and demographic data, we compared the demographic responses of hybrid and parental taxa to environmental variation across space and time. All three taxa had lower predicted population growth rates under warm, dry conditions. However, the magnitude of these responses varied among taxa and populations. Hybrids had consistently lower predicted population growth rates than P. hippiana. In contrast, hybrid performance relative to P. pulcherrima varied with population and climate, with the hybrid maintaining relatively stable growth rates while populations of P. pulcherrima shrank under warm, dry conditions. Our findings demonstrate that hybrids in this system are neither intrinsically unfit nor universally more vigorous than parents, suggesting that the demographic consequences of hybridization are context-dependent. Our results also imply that shifts to warmer and drier conditions could have particularly negative repercussions for P. pulcherrima, which is currently the most abundant taxon in the study area, possibly as a legacy of more favorable historical climates. More broadly, the distributions of these long-lived taxa are lagging behind their demographic trajectories, such that the currently less common P. hippiana could become the most abundant of the Potentilla taxa as this region continues to warm and dry.

Population structure and hybridization under contemporary and future climates in a heteroploid foundational shrub species (Artemisia tridentata)

Current and past climatic changes can shift plant climatic niches, which may cause spatial overlap or separation between related taxa. The former often leads to hybridization and introgression, which may generate novel variation and influence the adaptive capacity of plants. An additional mechanism facilitating adaptations to novel environments and an important evolutionary driver in plants is polyploidy as the result of whole genome duplication. Artemisia tridentata (big sagebrush) is a landscape-dominating foundational shrub in the western United States which occupies distinct ecological niches, exhibiting diploid and tetraploid cytotypes. Tetraploids have a large impact on the species' landscape dominance as they occupy a preponderance of the arid spectrum of A. tridentata range. Three distinct subspecies are recognized, which co-occur in ecotones - the transition zone between two or more distinct ecological niches - allowing for hybridization and introgression. Here we assess the genomic distinctiveness and extent of hybridization among subspecies at different ploidies under both contemporary and predicted future climates. We sampled five transects throughout the western United States where a subspecies overlap was predicted using subspecies-specific climate niche models. Along each transect, we sampled multiple plots representing the parental and the potential hybrid habitats. We performed reduced representation sequencing and processed the data using a ploidy-informed genotyping approach. Population genomic analyses revealed distinct diploid subspecies and at least two distinct tetraploid gene pools, indicating independent origins of the tetraploid populations. We detected low levels of hybridization (2.5%) between the diploid subspecies, while we found evidence for increased admixture between ploidy levels (18%), indicating hybridization has an important role in the formation of tetraploids. Our analyses highlight the importance of subspecies co-occurrence within these ecotones to maintain gene exchange and potential formation of tetraploid populations. Genomic confirmations of subspecies in the ecotones support the subspecies overlap predicted by the contemporary climate niche models. However, future mid-century projections of subspecies niches predict a substantial loss in range and subspecies overlap. Thus, reductions in hybridization potential could affect new recruitment of genetically variable tetraploids that are vital to this species' ecological role. Our results underscore the importance of ecotone conservation and restoration.

Variation in gene expression along an elevation gradient of Rhododendron sanguineum var. haemaleum assessed in a comparative transcriptomic analysis

Selection along environmental gradients may play a vital role in driving adaptive evolution. Nevertheless, genomic variation and genetic adaptation along environmental clines remains largely unknown in plants in alpine ecosystems. To close this knowledge gap, we assayed transcriptomic profiles of late flower bud and early leaf bud of Rhododendron sanguineum var. haemaleum from four different elevational belts between 3,000 m and 3,800 m in the Gaoligong Mountains. By comparing differences in gene expression of these samples, a gene co-expression network (WGCNA) was constructed to identify candidate genes related to elevation. We found that the overall gene expression patterns are organ-specific for the flower and leaf. Differentially expressed unigenes were identified in these organs. In flowers, these were mainly related to terpenoid metabolism (RsHMGR, RsTPS), while in leaves mainly related to anthocyanin biosynthesis (RsCHS, RsF3'5'H). Terpenoids are the main components of flower scent (fragrance) likely attracting insects for pollination. In response to fewer pollinators at higher elevation zone, it seems relatively less scent is produced in flower organs to reduce energy consumption. Secondary metabolites in leaves such as anthocyanins determine the plants' alternative adaptive strategy to extreme environments, such as selective pressures of insect herbivory from environmental changes and substrate competition in biosynthesis pathways at high elevations. Our findings indicated that the gene expression profiles generated from flower and leaf organs showed parallel expression shifts but with different functionality, suggesting the existence of flexibility in response strategies of plants exposed to heterogeneous environments across elevational gradients. The genes identified here are likely to be involved in the adaptation of the plants to these varying mountainous environments. This study thus contributes to our understanding of the molecular mechanisms of adaptation in response to environmental change.

Rapid evolution of hybrid breakdown following recent divergence with gene flow in Senecio species on Mount Etna, Sicily

How do nascent species evolve reproductive isolation during speciation with on-going gene flow? How do hybrid lineages become stabilised hybrid species? While commonly used genomic approaches provide an indirect way to identify species incompatibility factors, synthetic hybrids generated from interspecific crosses allow direct pinpointing of phenotypic traits involved in incompatibilities and the traits that are potentially adaptive in hybrid species. Here we report the analysis of phenotypic variation and hybrid breakdown in crosses between closely-related Senecio aethnensis and S. chrysanthemifolius, and their homoploid hybrid species, S. squalidus. The two former species represent a likely case of recent (<200 ky) speciation with gene flow driven by adaptation to contrasting conditions of high- and low-elevations on Mount Etna, Sicily. As these species form viable and fertile hybrids, it remains unclear whether they have started to evolve reproductive incompatibility. Our analysis represents the first study of phenotypic variation and hybrid breakdown involving multiple Senecio hybrid families. It revealed wide range of variation in multiple traits, including the traits previously unrecorded in synthetic hybrids. Leaf shape, highly distinct between S. aethnensis and S. chrysanthemifolius, was extremely variable in F₂ hybrids, but more consistent in S. squalidus. Our study demonstrates that interspecific incompatibilities can evolve rapidly despite on-going gene flow between the species. Further work is necessary to understand the genetic bases of these incompatibilities and their role in speciation with gene flow.

Evidence of hybrid breakdown among invasive hybrid cattails (Typha × glauca)

Interspecific hybridization has varied consequences for offspring fitness, with implications for the maintenance of species integrity. Hybrid vigour, when it occurs, can peak in first-generation (F1) hybrids and then decline in advanced-generation (F2+) hybrids. This hybrid breakdown, together with the processes affecting patterns of hybridization and hybrid fitness, determine the evolutionary stability of hybrid zones. An extensive hybrid zone in North America involving the cattails Typha latifolia, T. angustifolia, and their invasive hybrid T. × glauca is characterized by hybrid vigour among F1s, but the fitness of advanced-generation hybrids has not been studied. We compared seed germination and plant growth of T. latifolia (parental L), F1 T. × glauca (F1), hybrid backcrosses to T. angustifolia (bcA) and T. latifolia (bcL), and advanced-generation (F2) hybrids. Consistent with expectations under hybrid breakdown, we found reduced plant growth for F2 hybrids in comparison with F1s (plant height and above-ground biomass) and parental Ls (above-ground biomass). Backcrossed hybrids had intermediate measures of plant growth and bcLs were characterized by reduced seed germination in comparison with parental Ls. Hybrid breakdown could make the formation of F1s in North America finite because (1) hybridization among cattails is asymmetric, with T. angustifolia but not T. latifolia subject to genetic swamping, and (2) T. angustifolia is less common and subject to competitive displacement by F1s. Hybrid breakdown is therefore expected to reduce hybrid frequencies over time, contributing to the long-term maintenance of T. latifolia - the only native cattail in the study region.

Transcriptome Analysis Provides Insights into Potentilla bifurca Adaptation to High Altitude

Potentilla bifurca is widely distributed in Eurasia, including the Tibetan Plateau. It is a valuable medicinal plant in the Tibetan traditional medicine system, especially for the treatment of diabetes. This study investigated the functional gene profile of Potentilla bifurca at different altitudes by RNA-sequencing technology, including de novo assembly of 222,619 unigenes from 405 million clean reads, 57.64% of which were annotated in Nr, GO, KEGG, Pfam, and Swiss-Prot databases. The most significantly differentially expressed top 50 genes in the high-altitude samples were derived from plants that responded to abiotic stress, such as peroxidase, superoxide dismutase protein, and the ubiquitin-conjugating enzyme. Pathway analysis revealed that a large number of DEGs encode key enzymes involved in secondary metabolites, including phenylpropane and flavonoids. In addition, a total of 298 potential genomic SSRs were identified in this study, which provides information on the development of functional molecular markers for genetic diversity assessment. In conclusion, this study provides the first comprehensive assessment of the Potentilla bifurca transcriptome. This provides new insights into coping mechanisms for non-model organisms surviving in harsh environments at high altitudes, as well as molecular evidence for the selection of superior medicinal plants.

Extensive Sympatry and Frequent Hybridization of Ecologically Divergent Aquatic Plants on the Qinghai-Tibetan Plateau

Hybridization has fascinated biologists in recent centuries for its evolutionary importance, especially in plants. Hybrid zones are commonly located in regions across environmental gradients due to more opportunities to contact and ecological heterogeneity. For aquatic taxa, intrazonal character makes broad overlapping regions in intermediate environments between related species. However, we have limited information on the hybridization pattern of aquatic taxa in alpines, especially submerged macrophytes. In this study, we aimed to test the hypotheses that niche overlap and hybridization might be extensive in related aquatic plants across an altitudinal gradient. We evaluated the niche overlap in three related species pairs on the Qinghai-Tibetan Plateau and assessed the spatial pattern of hybrid populations. Obvious niche overlap and common hybridization were revealed in all three pairs of related aquatic plants. The plateau edge and river basins were broad areas for the sympatry of divergent taxa, where a large proportion of hybrid populations occurred. Hybrids are also discretely distributed in diverse habitats on the plateau. Differences in the extent of niche overlap, genetic incompatibility and phylogeographic history might lead to variation differences in hybridization patterns among the three species pairs. Our results suggested that plateau areas are a hotspot for ecologically divergent aquatic species to contact and mate and implied that hybridization may be important for the freshwater biodiversity of highlands.

Do habitat and elevation promote hybridization during secondary contact between three genetically distinct groups of warbling vireo (Vireo gilvus)?

Following postglacial expansion, secondary contact can occur between genetically distinct lineages. These genetic lineages may be associated with specific habitat or environmental variables and therefore, their distributions in secondary contact could reflect such conditions within these areas. Here we used mtDNA, microsatellite, and morphological data to study three genetically distinct groups of warbling vireo (Vireo gilvus) and investigate the role that elevation and habitat play in their distributions. We studied two main contact zones and within each contact zone, we examined two separate transects. Across the Great Plains contact zone, we found that hybridization between eastern and western groups occurs along a habitat and elevational gradient, whereas hybridization across the Rocky Mountain contact zone was not as closely associated with habitat or elevation. Hybrids in the Great Plains contact zone were more common in transitional areas between deciduous and mixed-wood forests, and at lower elevations (<1000 m). Hybridization patterns were similar along both Great Plains transects indicating that habitat and elevation play a role in hybridization between distinct eastern and western genetic groups. The observed patterns suggest adaptation to different habitats, perhaps originating during isolation in multiple Pleistocene refugia, is facilitating hybridization in areas where habitat types overlap.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

The Role of Interspecific Hybridisation in Adaptation and Speciation: Insights From Studies in Senecio

Hybridisation is well documented in many species, especially plants. Although hybrid populations might be short-lived and do not evolve into new lineages, hybridisaiton could lead to evolutionary novelty, promoting adaptation and speciation. The genus Senecio (Asteraceae) has been actively used to unravel the role of hybridisation in adaptation and speciation. In this article, we first briefly describe the process of hybridisation and the state of hybridisation research over the years. We then discuss various roles of hybridisation in plant adaptation and speciation illustrated with examples from different Senecio species, but also mention other groups of organisms whenever necessary. In particular, we focus on the genomic and transcriptomic consequences of hybridisation, as well as the ecological and physiological aspects from the hybrids' point of view. Overall, this article aims to showcase the roles of hybridisation in speciation and adaptation, and the research potential of Senecio, which is part of the ecologically and economically important family, Asteraceae.

Demographic decline and lineage-specific adaptations characterize New Zealand kiwi

Small and fragmented populations may become rapidly differentiated due to genetic drift, making it difficult to distinguish whether neutral genetic structure is a signature of recent demographic events, or of long-term evolutionary processes that could have allowed populations to adaptively diverge. We sequenced 52 whole genomes to examine Holocene demographic history and patterns of adaptation in kiwi (Apteryx), and recovered 11 strongly differentiated genetic clusters corresponding to previously recognized lineages. Demographic models suggest that all 11 lineages experienced dramatic population crashes relative to early- or mid-Holocene levels. Small population size is associated with low genetic diversity and elevated genetic differentiation (FST), suggesting that population declines have strengthened genetic structure and led to the loss of genetic diversity. However, population size is not correlated with inbreeding rates. Eight lineages show signatures of lineage-specific selective sweeps (284 sweeps total) that are unlikely to have been caused by demographic stochasticity. Overall, these results suggest that despite strong genetic drift associated with recent bottlenecks, most kiwi lineages possess unique adaptations and should be recognized as separate adaptive units in conservation contexts. Our work highlights how whole-genome datasets can address longstanding uncertainty about the evolutionary and conservation significance of small and fragmented populations of threatened species.

Natural hybridization among three Rhododendron species (Ericaceae) revealed by morphological and genomic evidence

BACKGROUND

Natural hybridization can influence the adaptive response to selection and accelerate species diversification. Understanding the composition and structure of hybrid zones may elucidate patterns of hybridization processes that are important to the formation and maintenance of species, especially for taxa that have experienced rapidly adaptive radiation. Here, we used morphological traits, ddRAD-seq and plastid DNA sequence data to investigate the structure of a Rhododendron hybrid zone and uncover the hybridization patterns among three sympatric and closely related species.

RESULTS

Our results show that the hybrid zone is complex, where bi-directional hybridization takes place among the three sympatric parental species: R. spinuliferum, R. scabrifolium, and R. spiciferum. Hybrids between R. spinuliferum and R. spiciferum (R. ×duclouxii) comprise multiple hybrid classes and a high proportion of F₁ generation hybrids, while a novel hybrid taxon between R. spinuliferum and R. scabrifolium dominated the F₂ generation, but no backcross individuals were detected. The hybrid zone showed basically coincident patterns of population structure between genomic and morphological data.

CONCLUSIONS

Natural hybridization exists among the three Rhododendron species in the hybrid zone, although patterns of hybrid formation vary between hybrid taxa, which may result in different evolutionary outcomes. This study represents a unique opportunity to dissect the ecological and evolutionary mechanisms associated with adaptive radiation of Rhododendron species in a biodiversity hotspot.

Exploring potentialities of avian genomic research in Nepalese Himalayas

Nepal, a small landlocked country in South Asia, holds about 800 km of Himalayan Mountain range including the Earth's highest mountain. Within such a mountain range in the north and plain lowlands in the south, Nepal provides a habitat for about 9% of global avian fauna. However, this diversity is underrated because of the lack of enough studies, especially using molecular tools to quantify and understand the distribution patterns of diversity. In this study, we reviewed the studies in the last two decades (2000‒2019) that used molecular methods to study the biodiversity in Nepal to examine the ongoing research trend and focus. Although Nepalese Himalaya has many opportunities for cutting-edge molecular research, our results indicated that the rate of genetic/genomic studies is much slower compared to the regional trends. We found that genetic research in Nepal heavily relies on resources from international institutes and that too is mostly limited to research on species monitoring, distribution, and taxonomic validations. Local infrastructures to carry out cutting-edge genomic research in Nepal are still in their infancy and there is a strong need for support from national/international scientists, universities, and governmental agencies to expand such genomic infrastructures in Nepal. We particularly highlight avian fauna as a potential future study system in this region that can be an excellent resource to explore key biological questions such as understanding eco-physiology and molecular basis of organismal persistence to changing environment, evolutionary processes underlying divergence and speciation, or mechanisms of endemism and restrictive distribution of species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s40657-021-00290-5.

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones

PREMISE

The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data.

METHODS

We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape.

RESULTS

Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations.

CONCLUSIONS

The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

Spatiotemporal maintenance of flora in the Himalaya biodiversity hotspot: Current knowledge and future perspectives

Mountain ecosystems support a significant one-third of all terrestrial biodiversity, but our understanding of the spatiotemporal maintenance of this high biodiversity remains poor, or at best controversial. The Himalaya hosts a complex mountain ecosystem with high topographic and climatic heterogeneity and harbors one of the world's richest floras. The high species endemism, together with increasing anthropogenic threats, has qualified the Himalaya as one of the most significant global biodiversity hotspots. The topographic and climatic complexity of the Himalaya makes it an ideal natural laboratory for studying the mechanisms of floral exchange, diversification, and spatiotemporal distributions. Here, we review literature pertaining to the Himalaya in order to generate a concise synthesis of the origin, distribution, and climate change responses of the Himalayan flora. We found that the Himalaya supports a rich biodiversity and that the Hengduan Mountains supplied the majority of the Himalayan floral elements, which subsequently diversified from the late Miocene onward, to create today's relatively high endemicity in the Himalaya. Further, we uncover links between this Miocene diversification and the joint effect of geological and climatic upheavals in the Himalaya. There is marked variance regarding species dispersal, elevational gradients, and impact of climate change among plant species in the Himalaya, and our review highlights some of the general trends and recent advances on these aspects. Finally, we provide some recommendations for conservation planning and future research. Our work could be useful in guiding future research in this important ecosystem and will also provide new insights into the maintenance mechanisms underpinning other mountain systems.

The Evolutionary History, Diversity, and Ecology of Willows (Salix L.) in the European Alps

The genus Salix (willows), with 33 species, represents the most diverse genus of woody plants in the European Alps. Many species dominate subalpine and alpine types of vegetation. Despite a long history of research on willows, the evolutionary and ecological factors for this species richness are poorly known. Here we will review recent progress in research on phylogenetic relationships, evolution, ecology, and speciation in alpine willows. Phylogenomic reconstructions suggest multiple colonization of the Alps, probably from the late Miocene onward, and reject hypotheses of a single radiation. Relatives occur in the Arctic and in temperate Eurasia. Most species are widespread in the European mountain systems or in the European lowlands. Within the Alps, species differ ecologically according to different elevational zones and habitat preferences. Homoploid hybridization is a frequent process in willows and happens mostly after climatic fluctuations and secondary contact. Breakdown of the ecological crossing barriers of species is followed by introgressive hybridization. Polyploidy is an important speciation mechanism, as 40% of species are polyploid, including the four endemic species of the Alps. Phylogenomic data suggest an allopolyploid origin for all taxa analyzed so far. Further studies are needed to specifically analyze biogeographical history, character evolution, and genome evolution of polyploids.

Hybridization, ecogeographical displacement and the emergence of new lineages - A genotyping-by-sequencing and ecological niche and species distribution modelling study of Sempervivum tectorum L. (Houseleek)

Ecogeographical displacement of homoploid hybrid lineages from their parents is well documented and considered an important mechanism to achieve reproductive isolation. In this study, we investigated the origin of the flowering plant species Sempervivum tectorum in the Massif Central (France) through homoploid hybridization between lineages of the species from the Rhine Gorge area (Germany) and the Pyrenees (France). We used genotyping-by-sequencing genetic data as evidence for the hybrid origin of the Massif Central lineage, and WorldClim climatic data and soil pH and soil temperature data collected by us for ecological niche and species distribution modelling. We could show that the Massif Central lineage shows hybrid admixture and that the niche of this lineage is significantly different from those of the parental lineages. In comparison with the parental niches, different variables of the niche of the hybrid lineage are intermediate, parental-combined or extreme. The different niche of the Massif Central populations thus can plausibly be interpreted as hybridization-derived. Our species distribution modelling for the Last Glacial Maximum and Mid-Holocene showed that the potential distribution of the hybrid lineage at the likely time of its origin in the Quaternary possibly was parapatric in relation to the largely sympatric distributions of the parental lineages. We hypothesize that reproductive isolation of the hybrid lineage from the parental lineages resulted from the segregation of distribution ranges by a differential response of the three lineages to a warming climate.

Biased Hybridization and Its Impact on Adaptive Introgression

Gene exchange between species can influence ecological and evolutionary processes ranging from population rescue to adaptive radiation. Genomic tools have provided new insights into the prevalence and nature of gene exchange between species. However, much remains unknown of how ecological, behavioral, and evolutionary factors determine what genetic variation moves between species in the first place. In particular, more research is needed that evaluates whether such factors bias gene flow from one species to another, and whether any such biases affect how genetic variation from another species is ultimately retained in the genome of a given species. Addressing this issue is crucial in a changing world where hybridization and introgression might determine which species succeed and which become extinct.

Recurrent hybridization underlies the evolution of novelty in Gentiana (Gentianaceae) in the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) and adjacent areas are centres of diversity for several alpine groups. Although it is known that the QTP acted as a source area for diversification of the alpine genus Gentiana, the evolutionary processes underlying diversity in this genus, especially the formation of narrow endemics, are still poorly understood. Hybridization has been proposed as a driver of plant endemism in the QTP but few cases have been documented with genetic data. Here, we describe a new endemic species in Gentiana section Cruciata as G. hoae sp. nov., and explore its evolutionary history with complete plastid genomes and nuclear ribosomal internal transcribed spacer sequence data. Genetic divergence within G. hoae ~3 million years ago was followed by postglacial expansion on the QTP, suggesting Pleistocene glaciations as a key factor shaping the population history of G. hoae. Furthermore, a mismatch between plastid and nuclear data suggest that G. hoae participated in historical hybridization, while population sequencing show this species continues to hybridize with the co-occurring congener G. straminea in three locations. Our results indicate that hybridization may be a common process in the evolution of Gentiana and may be widespread among recently diverged taxa of the QTP.

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST -FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST ) and quantitative genetic (Q_ST ) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass-related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.

Analyses of Hybrid Viability across a Hybrid Zone between Two Alnus Species Using Microsatellites and cpDNA Markers

Diploid Alnus glutinosa s. str. and autotetraploid A. rohlenae form a narrow hybrid zone in a study area in southern Serbia, which results in triploid hybrid formation. The vast majority of previous studies have been focused on studies of maternal plants, but the offspring resulting from their crossing have not been much studied. Here, we use the variability of microsatellites and chloroplast DNA between these species and their putative hybrids to create an overall picture of the development of the hybrid zone and its predicted type. To elucidate the gene transfer within both species, the origins of individual ploidies and especially the role of triploid hybrids, a germination experiment was carried out linked with a flow cytometry study of the resulting seedlings. The tension zone model seems to offer the most adequate explanation of our observations, with selection against triploid hybrids and the spatial positioning of the hybrid zone. Despite selection against them, the triploid hybrids play an important role in the exchange of genes between the two species and therefore serve as a bridge for introgression. The presence of fertile triploids is essential for enriching the haplotype diversity between these species and for the development of new genetic lineages.

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental-like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.

Fine-scale species distribution modelling and genotyping by sequencing to examine hybridisation between two narrow endemic plant species

Hybridization has an important and often positive role in plant evolution. However, it can also have negative consequences for species. Two closely related species of Ornduffia are endemic to the Porongurup Range in the South West Australian Global Biodiversity Hotspot. The rare Ornduffia calthifolia is found exclusively on the summits, while O. marchantii is more widely dispersed across a greater range of elevation and is not considered threatened. Hybridisation in suitable overlapping habitat has been suspected between them for decades. Here we combine genotyping by sequencing to verify hybridisation genetically, and fine scale (2 m resolution) species distribution modelling (SDM) to test if hybrids occur in suitable intersecting habitat. From a study area of c. 4700 ha, SDM identified c. 275 ha and c. 322 ha of suitable habitat for O. calthifolia and O. marchantii, respectively. We identified range overlap between species of c. 59 ha), which enveloped 32 individuals confirmed to be hybrids. While the hybrids were at the margin of suitable habitat for O. marchantii, their preference for elevated habitat was closer to the more narrowly distributed O. calthifolia. The combination of genetic data and fine scale spatial modelling approaches enabled a better understanding of hybridisation among taxa of conservation significance. However, the level to which hybrid proliferation and competition for habitat presents as a threat to O. calthifolia is currently unknown and requires priority in conservation management given the threats from global warming and disturbance by tourism.

Biogeographical, ecological, and phylogenetic analyses clarifying the evolutionary history of Calibrachoa in South American grasslands

Calibrachoa is a charismatic South American genus of Solanaceae, closely related to Petunia, which encompasses approximately 30 species. Studies that were based solely on plastid molecular markers indicated the monophyly of the genus and distributed its species in two subgenera; to date no phylogeny has included a broad morphological variants and nuclear markers. Here, we present a phylogenetic analysis based on eight plastid and eight nuclear markers that cover the most extensive geographic distribution for the genus. We use this phylogeny to infer the biogeographic history of the genus and to understand the primary drivers for species diversification. Our results yield a fully supported tree where monophyly is confirmed to genus and subgenera. The species of Stimomphis subgenus that were previously considered uncertain, here emerge in four highly supported clades. The hypothesis of niche conservatism is confirmed, and adaptive radiation explains the species diversification. The lowlands are the most likely ancestral area of the genus, subgenera, and two clades of Stimomphis subgenus. Our results constitute an excellent starting point for further evolutionary and taxonomic studies and explain several uncertain evolutionary relationships in the group and the evolution of their distribution.

Influence of a climatic gradient on genetic exchange between two oak species

PREMISE

In plant groups with limited intrinsic barriers to gene flow, it is thought that environmental conditions can modulate interspecific genetic exchange. Oaks are known for limited barriers to gene flow among closely related species. Here, we use Quercus as a living laboratory in which to pursue a fundamental question in plant evolution: Do environmental gradients restrict or promote genetic exchange between species?

METHODS

We focused on two North American oaks, the rare Quercus dumosa and the widespread Q. berberidifolia. We sampled intensively along a contact zone in California, USA. We sequenced restriction site-associated DNA markers and measured vegetative phenotype. We tested for genetic exchange, the association with climate, and the effect on phenotype.

RESULTS

There is evidence for genetic exchange between the species. Admixed plants are found in areas of intermediate climate, while less admixed plants are found at the extremes of the climatic gradient. Genetic and phenotypic patterns are out of phase in the contact zone; some plants display the phenotype of one species but are genetically associated with another.

CONCLUSIONS

Our results support the hypothesis that a strong climatic gradient can promote genetic exchange between species. The overall weak correlation between genotype and phenotype in the contact zone between the species suggests that genetic exchange can lead to the breakdown of trait combinations used to define species. This incongruency predicts ongoing problems for conservation of Q. dumosa, with implications for conservation of other oaks.

Incomplete reproductive isolation between Rhododendron taxa enables hybrid formation and persistence

The evolutionary consequences of hybridization ultimately depend on the magnitude of reproductive isolation between hybrids and their parents. We evaluated the relative contributions of pre- and post-zygotic barriers to reproduction for hybrid formation, hybrid persistence and potential for reproductive isolation of hybrids formed between two Rhododendron species, R. spiciferum and R. spinuliferum. Our study established that incomplete reproductive isolation promotes hybrid formation and persistence and delays hybrid speciation. All pre-zygotic barriers to reproduction leading to hybrid formation are incomplete: parental species have overlapping flowering; they share the same pollinators; reciprocal assessments of pollen tube germination and growth do not differ among parents. The absence of post-zygotic barriers between parental taxa indicates that the persistence of hybrids is likely. Reproductive isolation was incomplete between hybrids and parents in all cases studied, although asymmetric differences in reproductive fitness were prevalent and possibly explain the genetic structure of natural hybrid swarms where hybridization is known to be bidirectional but asymmetric. Introgression, rather than speciation, is a probable evolutionary outcome of hybridization between the two Rhododendron taxa. Our study provides insights into understanding the evolutionary implications of natural hybridization in woody plants.

Species Boundaries and Parapatric Speciation in the Complex of Alpine Shrubs, Rosa sericea (Rosaceae), Based on Population Genetics and Ecological Tolerances

Discerning species boundaries among closely related taxa is fundamental to studying evolution and biodiversity. However, species boundaries can be difficult to access in plants because ongoing divergence and speciation may leave an evolutionary footprint similar to introgression, which occurs frequently among species and genera. In this study, we sought to determine species boundaries between two closely related alpine shrubs, Rosa sericea and Rosa omeiensis, using population genetics, environmental data and ecological niche modeling, and morphological traits. We analyzed populations of R. sericea and R. omeiensis using genetic markers comprising a fragment of the single-copy nuclear gene, LEAFY, micro-satellites (EST-SSR), and plastid DNA sequences. The DNA sequence data suggested clusters of populations consistent with geography but not with previously proposed species boundaries based on morphology. Nevertheless, we found that the ecological niches of the previously proposed species only partially overlap. Thus, we suspect that these species are in the process of parapatric speciation; that is, differentiating along an ecological gradient, so that they exhibit differing morphology. Morphology has previously been the basis of recognizing the species R. sericea and R. omeiensis, which are the most widely distributed species within a broader R. sericea complex that includes several other narrow endemics. Here, we recognize R. sericea and R. omeiensis as independent species based on morphological and ecological data under the unified species concept, which emphasizes that these data types are of equal value to DNA for determining species boundaries and refining taxonomic treatments. While the DNA data did not delimit species within the R. sericea complex, we expect to develop and utilize new, robust DNA tools for understanding speciation within this group in future studies.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

Coadapted genomes and selection on hybrids: Fisher's geometric model explains a variety of empirical patterns

Natural selection plays a variety of roles in hybridization, speciation, and admixture. Most research has focused on two extreme cases: crosses between closely related inbred lines, where hybrids are fitter than their parents, or crosses between effectively isolated species, where hybrids suffer severe breakdown. But many natural populations must fall into intermediate regimes, with multiple types of gene interaction, and these are more difficult to study. Here, we develop a simple fitness landscape model, and show that it naturally interpolates between previous modeling approaches, which were designed for the extreme cases, and invoke either mildly deleterious recessives, or discrete hybrid incompatibilities. Our model yields several new predictions, which we test with genomic data from Mytilus mussels, and published data from plants (Zea, Populus, and Senecio) and animals (Mus, Teleogryllus, and Drosophila). The predictions are generally supported, and the model explains a number of surprising empirical patterns. Our approach enables novel and complementary uses of genome-wide datasets, which do not depend on identifying outlier loci, or "speciation genes" with anomalous effects. Given its simplicity and flexibility, and its predictive successes with a wide range of data, the approach should be readily extendable to other outstanding questions in the study of hybridization.

Phylogeography of two closely related species of Allium endemic to East Asia: Population evolution in response to climate oscillations

This study investigated the effects of climate oscillations on the evolution of two closely related Allium species, A. neriniflorum and A. tubiflorum. We sequenced three cp DNA (cpDNA) fragments (rps16, rpl32-trnL, and trnD-trnT, together approximately 2,500 bp in length) of two closely related Allium species, with samples from 367 individuals in 47 populations distributed across the total range of these species. The interspecific and intraspecific divergence times of the two species were in the Quaternary glaciation. The population divergence was high for the cpDNA variation, suggesting a significant phylogeographic structure (NST = 0.844, GST = 0.798, p < 0.05). Remarkable ecological differentiation was also revealed by Niche models and statistical analyses. Our results suggest the speciation event of the two species was triggered by violent climatic changes during the Quaternary glaciation.

Fitness of natural willow hybrids in a pioneer mosaic hybrid zone

Hybrid fitness is an important parameter to predict the evolutionary consequences of a hybridization event and to characterize hybrid zones. We studied fitness parameters of F1 and later-generation hybrids between the lowland species Salix purpurea and the alpine S. helvetica that have recently emerged during colonization of an alpine glacier forefield. Fruit production (number of capsules per catkin and fruit set) did not differ between hybrids and parents, but the number of seeds per capsule of F1 hybrids was slightly lower than that of later-generation hybrids and of the parents. Germination rates and seedling growth were tested on three substrates (pH 4.5, 7.0, and 8.0). Germination rates of seeds collected from F1 hybrids were lower on acid and neutral substrates, but equal at pH 8.0 compared to all other groups, while the seeds from later-generation hybrids performed as well as the parents on all three substrates. In seedling growth, the colonizer S. purpurea performed better than all other taxa on all three substrates, while hybrids resembled the subalpine species S. helvetica. Results suggest that endogenous selection acts against F1 hybrids, but favors fitter genotypes in later-generation hybrids. Exogenous selection via soil pH appears to be weak during seedling establishment. The pioneer vegetation on the glacier forefield may offer sufficient niche space for hybrid seedlings. Owing to the relatively high fitness of the hybrids and the scattered distribution of hybrids and parental individuals on the glacier forefield, this hybrid zone can be assigned to a mosaic model, probably facilitating gene flow and introgression between the parental species. As establishment of the hybrid zone appears to be linked to a colonization process, we propose to call it a pioneer mosaic hybrid zone.

Genomics of Natural Populations: How Differentially Expressed Genes Shape the Evolution of Chromosomal Inversions in Drosophila pseudoobscura

Chromosomal rearrangements can shape the structure of genetic variation in the genome directly through alteration of genes at breakpoints or indirectly by holding combinations of genetic variants together due to reduced recombination. The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for >30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described >80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. In this study, we tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, we asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. We find evidence for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. We find the majority of loci show additivity in heterokaryotypes. Our results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds.

Genomic architecture of phenotypic divergence between two hybridizing plant species along an elevational gradient

Knowledge of the genetic basis of phenotypic divergence between species and how such divergence is caused and maintained is crucial to an understanding of speciation and the generation of biodiversity. The hybrid zone between Senecio aethnensis and S. chrysanthemifolius on Mount Etna, Sicily, provides a well-studied example of species divergence in response to conditions at different elevations, despite hybridization and gene flow. Here, we investigate the genetic architecture of divergence between these two species using a combination of quantitative trait locus (QTL) mapping and genetic differentiation measures based on genetic marker analysis. A QTL architecture characterized by physical QTL clustering, epistatic interactions between QTLs, and pleiotropy was identified, and is consistent with the presence of divergent QTL complexes resistant to gene flow. A role for divergent selection between species was indicated by significant negative associations between levels of interspecific genetic differentiation at mapped marker gene loci and map distance from QTLs and hybrid incompatibility loci. Within-species selection contributing to interspecific differentiation was evidenced by negative associations between interspecific genetic differentiation and genetic diversity within species. These results show that the two Senecio species, while subject to gene flow, maintain divergent genomic regions consistent with local selection within species and selection against hybrids between species which, in turn, contribute to the maintenance of their distinct phenotypic differences.

Demographic history of speciation in a Senecio altitudinal hybrid zone on Mt. Etna

Hybrid zones typically form as a result of species coming into secondary contact, but can also be established in situ as an ecotonal hybrid zone, a situation which has been reported far less frequently. An altitudinal hybrid zone on Mount Etna between two ragwort species (the low elevation Senecio chrysanthemifolius and high elevation S. aethnensis) could potentially represent either of these possibilities. However, a scenario of secondary contact vs. speciation with gene flow has not been explicitly tested. Here, we test these alternatives and demonstrate that the data do not support secondary contact. Furthermore, we report that the previous analyses of speciation history of these species were based on admixed populations, which has led to inflated estimates of ongoing, interspecific gene flow. Our new analyses, based on 'pure' S. aethnensis and S. chrysanthemifolius populations, reveal gene exchange of less than one effective migrant per generation, a level low enough to allow the species to accumulate neutral, genomewide differences. Overall, our results are consistent with a scenario of speciation with gene flow and a divergence time which coincides with the rise of Mt. Etna to altitudes above 2000 m (~150 KY). Further work to quantify the role of adaptation to contrasting environments of high and low altitudes will be needed to support the scenario of recent ecological speciation in this system.

Genetic origin and composition of a natural hybrid poplar Populus × jrtyschensis from two distantly related species

BACKGROUND

The factors that contribute to and maintain hybrid zones between distinct species are highly variable, depending on hybrid origins, frequencies and fitness. In this study, we aimed to examine genetic origins, compositions and possible maintenance of Populus × jrtyschensis, an assumed natural hybrid between two distantly related species. This hybrid poplar occurs mainly on the floodplains along the river valleys between the overlapping distributions of the two putative parents.

RESULTS

We collected 566 individuals from 45 typical populations of P. × jrtyschensis, P. nigra and P. laurifolia. We genotyped them based on the sequence variations of one maternally inherited chloroplast DNA (cpDNA) fragment and genetic polymorphisms at 20 SSR loci. We further sequenced eight nuclear genes for 168 individuals from 31 populations. Two groups of cpDNA haplotypes characteristic of P. nigra and P. laurifolia respectively were both recovered for P. × jrtyschensis. Genetic structures and coalescent tests of two sets of nuclear population genetic data suggested that P. × jrtyschensis originated from hybridizations between the two assumed parental species. All examined populations of P. × jrtyschensis comprise mainly F1 hybrids from interspecific hybridizations between P. nigra and P. laurifolia. In the habitats of P. × jrtyschensis, there are lower concentrations of soil nitrogen than in the habitats occupied by the other two species.

CONCLUSIONS

Our extensive examination of the genetic composition of P. × jrtyschensis suggested that it is typical of F1-dominated hybrid zones. This finding plus the low concentration of soil nitrogen in the floodplain soils support the F1-dominated bounded hybrid superiority hypothesis of hybrid zone maintenance for this particular hybrid poplar.