Past climate changes facilitated homoploid speciation in three mountain spiny fescues (Festuca, Poaceae)

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.

Hybridization and introgression events in cooccurring populations of closely related grasses (Poaceae: Stipa) in high mountain steppes of Central Asia

Stipa is a genus comprising ca. 150 species found in warm temperate regions of the Old World and around 30% of its representatives are of hybrid origin. In this study, using integrative taxonomy approach, we tested the hypothesis that hybridization and introgression are the explanations of the morphological intermediacy in species belonging to Stipa sect. Smirnovia, one of the species-rich sections in the mountains of Central Asia. Two novel nothospecies, S. magnifica × S. caucasica subsp. nikolai and S. lingua × S. caucasica subsp. nikolai, were identified based on a combination of morphological characters and SNPs markers. SNPs marker revealed that all S. lingua × S. caucasica samples were F1 hybrids, whereas most of S. magnifica × S. caucasica samples were backcross hybrids. Furthermore, the above mentioned hybrids exhibit transgressive morphological characters to each of their parental species. These findings have implications for understanding the process of hybridization in the genus Stipa, particularly in the sect. Smirnovia. As a taxonomic conclusion, we describe the two new nothospecies S. × muksuensis (from Tajikistan) and S. × ochyrae (from Kyrgyzstan) and present an identification key to species morphologically similar to the taxa mentioned above.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Identifying Climatic Drivers of Hybridization with a New Ancestral Niche Reconstruction Method

Applications of molecular phylogenetic approaches have uncovered evidence of hybridization across numerous clades of life, yet the environmental factors responsible for driving opportunities for hybridization remain obscure. Verbal models implicating geographic range shifts that brought species together during the Pleistocene have often been invoked, but quantitative tests using paleoclimatic data are needed to validate these models. Here, we produce a phylogeny for Heuchereae, a clade of 15 genera and 83 species in Saxifragaceae, with complete sampling of recognized species, using 277 nuclear loci and nearly complete chloroplast genomes. We then employ an improved framework with a coalescent simulation approach to test and confirm previous hybridization hypotheses and identify one new intergeneric hybridization event. Focusing on the North American distribution of Heuchereae, we introduce and implement a newly developed approach to reconstruct potential past distributions for ancestral lineages across all species in the clade and across a paleoclimatic record extending from the late Pliocene. Time calibration based on both nuclear and chloroplast trees recovers a mid- to late-Pleistocene date for most inferred hybridization events, a timeframe concomitant with repeated geographic range restriction into overlapping refugia. Our results indicate an important role for past episodes of climate change, and the contrasting responses of species with differing ecological strategies, in generating novel patterns of range contact among plant communities and therefore new opportunities for hybridization. The new ancestral niche method flexibly models the shape of niche while incorporating diverse sources of uncertainty and will be an important addition to the current comparative methods toolkit. [Ancestral niche reconstruction; hybridization; paleoclimate; pleistocene.].

Phylogenomics reveals widespread hybridization and polyploidization in Henckelia (Gesneriaceae)

BACKGROUND AND AIMS

Hybridization has long been recognized as an important process for plant evolution and is often accompanied by polyploidization, another prominent force in generating biodiversity. Despite its pivotal importance in evolution, the actual prevalence and distribution of hybridization across the tree of life remain unclear.

METHODS

We used whole-genome shotgun (WGS) sequencing and cytological data to investigate the evolutionary history of Henckelia, a large genus in the family Gesneriaceae with a high frequency of suspected hybridization and polyploidization events. We generated WGS sequencing data at about 10× coverage for 26 Chinese Henckelia species plus one Sri Lankan species. To untangle the hybridization history, we separately extracted whole plastomes and thousands of single-copy nuclear genes from the sequencing data, and reconstructed phylogenies based on both nuclear and plastid data. We also explored sources of both genealogical and cytonuclear conflicts and identified signals of hybridization and introgression within our phylogenomic dataset using several statistical methods. Additionally, to test the polyploidization history, we evaluated chromosome counts for 45 populations of the 27 Henckelia species studied.

KEY RESULTS

We obtained well-supported phylogenetic relationships using both concatenation- and coalescent-based methods. However, the nuclear phylogenies were highly inconsistent with the plastid phylogeny, and we observed intensive discordance among nuclear gene trees. Further analyses suggested that both incomplete lineage sorting and gene flow contributed to the observed cytonuclear and genealogical discordance. Our analyses of introgression and phylogenetic networks revealed a complex history of hybridization within the genus Henckelia. In addition, based on chromosome counts for 27 Henckelia species, we found independent polyploidization events occurred within Henckelia after different hybridization events.

CONCLUSIONS

Our findings demonstrated that hybridization and polyploidization are common in Henckelia. Furthermore, our results revealed that H. oblongifolia is not a member of the redefined Henckelia and they suggested several other taxonomic treatments in this genus.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

5S-IGS rDNA in wind-pollinated trees (Fagus L.) encapsulates 55 million years of reticulate evolution and hybrid origins of modern species

Standard models of plant speciation assume strictly dichotomous genealogies in which a species, the ancestor, is replaced by two offspring species. The reality in wind-pollinated trees with long evolutionary histories is more complex: species evolve from other species through isolation when genetic drift exceeds gene flow; lineage mixing can give rise to new species (hybrid taxa such as nothospecies and allopolyploids). The multi-copy, potentially multi-locus 5S rDNA is one of few gene regions conserving signal from dichotomous and reticulate evolutionary processes down to the level of intra-genomic recombination. Therefore, it can provide unique insights into the dynamic speciation processes of lineages that diversified tens of millions of years ago. Here, we provide the first high-throughput sequencing (HTS) of the 5S intergenic spacers (5S-IGS) for a lineage of wind-pollinated subtropical to temperate trees, the Fagus crenata - F. sylvatica s.l. lineage, and its distant relative F. japonica. The observed 4963 unique 5S-IGS variants reflect a complex history of hybrid origins, lineage sorting, mixing via secondary gene flow, and intra-genomic competition between two or more paralogous-homoeologous 5S rDNA lineages. We show that modern species are genetic mosaics and represent a striking case of ongoing reticulate evolution during the past 55 million years.

Genesis, Evolution, and Genetic Diversity of the Hexaploid, Narrow Endemic Centaurea tentudaica

Within the genus Centaurea L., polyploidy is very common, and it is believed that, as to all angiosperms, it was key in the history of its diversification and evolution. Centaurea tentudaica is a hexaploid from subsect. Chamaecyanus of unknown origin. In this study, we examined the possible autopolyploid or allopolyploid origin using allozymes and sequences of three molecular markers: nuclear-ribosomic region ETS, and low-copy genes AGT1 and PgiC. We also included three species geographically and morphologically close to C. tentudaica: C. amblensis, C. galianoi, and C. ornata. Neighbor-Net and Bayesian analyses show a close relationship between C. amblensis and C. tentudaica and no relationship to any of the other species, which suggest that C. tentudaica is an autopolyploid of C. amblensis. Allozyme banding pattern also supports the autopolyploidy hypothesis and shows high levels of genetic diversity in the polyploid, which could suggest multiple origins by recurrent crosses of tetraploid and diploid cytotypes of C. amblensis. Environmental niche modeling was used to analyze the distribution of the possible parental species during the present, Last Glacial Maximum (LGM), Last Interglacial Period (LIG), and Penultimate Glacial Maximum (PGM) environmental conditions. Supporting the molecular suggestions that C. tentudaica originated from C. amblensis, environmental niche modeling confirms that past distribution of C. amblensis overlapped with the distribution of C. tentudaica.

Angiosperms at the edge: Extremity, diversity, and phylogeny

A hallmark of flowering plants is their ability to invade some of the most extreme and dynamic habitats, including cold and dry biomes, to a far greater extent than other land plants. Recent work has provided insight to the phylogenetic distribution and evolutionary mechanisms which have enabled this success, yet needed is a synthesis of evolutionary perspectives with plant physiological traits, morphology, and genomic diversity. Linking these disparate components will not only lead to better understand the evolutionary parallelism and diversification of plants with these two strategies, but also to provide the framework needed for directing future research. We summarize the primary physiological and structural traits involved in response to cold- and drought stress, outline the phylogenetic distribution of these adaptations, and describe the recurring association of these changes with rapid diversification events that occurred in multiple lineages over the past 15 million years. Across these threefold facets of dry-cold correlation (traits, phylogeny, and time) we stress the contrast between (a) the amazing diversity of solutions flowering plants have developed in the face of extreme environments and (b) a broad correlation between cold and dry adaptations that in some cases may hint at deep common origins.

The protected tree Dimorphandra wilsonii (Fabaceae) is a population of inter-specific hybrids: recommendations for conservation in the Brazilian Cerrado/Atlantic Forest ecotone

BACKGROUNDS AND AIMS

Dimorphandra wilsonii Rizzini, a critically endangered and protected tree, has a restricted distribution in the ecotone between the Cerrado and the Atlantic Forest in south-eastern Brazil. In this area, it co-occurs with D. mollis Benth., a common tree from the Cerrado, and D. exaltata Schott., a rare tree from the Atlantic Forest. Previous studies of D. wilsonii indicated heterozygosity excess at the individual level. Field observation of some intermediate phenotypes between D. wilsonii and both congeners suggests hybridization of D. wilsonii with D. mollis and/or D. exaltata. Here, we tested the hypothesis that D. wilsonii may have originated from hybridization between D. exaltata and D. mollis. We also performed cytogenetic analysis to examine if the heterozygosity excess could be explained by polyploidy in D. wilsonii.

METHODS

We evaluated the genetic diversity and population structure of D. wilsonii using 11 nuclear simple sequence repeats (SSRs) genotyped in 152 individuals sampled across the taxon's range. We performed comparative genetic analyses using overlapping SSR markers between D. wilsonii and previously published SSR data in D. mollis and D. exaltata to subsequently perform a series of allelic comparisons, multivariate and Bayesian analysis.

KEY RESULTS

Our results suggest that D. wilsonii individuals are most likely to correspond to F1 hybrids between D. exaltata and D. mollis. Cytogenetic analysis indicated that D. wilsonii is diploid with the same chromosome number as D. mollis (2n = 2x = 28).

CONCLUSIONS

Our study raises questions about the taxonomic status and the evolutionary future of D. wilsonii. We suggest that the conservation and management strategy for D. wilsonii should be revised and that it should take into account both parental Dimorphandra species in the ecotone, with special emphasis on the threatened D. exaltata. Finally, this study highlights the value of genetic information for the design of conservation strategies.

The transferability of microsatellite loci from a homoploid to a polyploid hybrid complex: an example from fine-leaved Festuca species (Poaceae)

Background

Microsatellite loci, or single sequence repeats (SSR), are widely used as powerful markers in population genetics. They represent an attractive tool for studying plants such as grasses, whose evolution is driven by hybridisation and polyploidization. However, the development of microsatellite markers has been challenging and time-consuming, especially for non-model organisms lacking available genome-wide sequence data. One straightforward and low-cost approach is to transfer the SSR loci developed for one species, or complex, to another closely-related one. This work evaluates the transferability of microsatellite loci from homoploid to allopolyploid complexes of fine-leaved Festuca species and to assess their use in two new species. The studied complex (F. amethystina-F. tatrae) is a useful model for research on the local adaptability of grasses with different ploidy levels. Since both species can be considered as rare or threatened (F. tatrae-as a mountain and narrow endemic species and F. amethystina-a mountain species with relict lowland populations), any tool enabling studies on genetic diversity and population genetics, such as SSR markers, could also be very useful in a conservation context.

Methods

The ploidy level within populations was estimated using flow cytometry. One diploid and one tetraploid population of F. amethystina and a diploid population of F. tatrae were chosen to test the transferability of SSR loci. Because our work describes the transfer of SSR nuclear markers designed originally for F. gautieri, a phylogenetic tree was prepared based on the ITS marker to assess the genetic distance between the studied complexes. The PCR products were separated on a high-resolution agarose gel, intended for SSR marker analysis. Appropriate solutions for the allotetraploid population and whole mixed-ploidy complex were implemented.

Results

Flow cytometry confirmed earlier data regarding DNA content in the investigated species and cytotypes. The phylogenetic ITS tree indicated a small genetic distance between F. gautieri complexes and the studied species. Ten microsatellite markers were successfully transferred. All markers were polymorphic. In total, 163 different alleles were scored from the 10 SSR loci. PCoA of accessions revealed well-separated groups corresponding to studied populations. Over 60% of the total variance is explained by differentiation within populations and one third among them.

Conclusions

The transferred markers are valid tools for the study of population genetics and inheritance relationships within cytotypes and species and between them. The presented markers can be used to study inbreeding depression in the Festuca species, and variations in the degrees of genetic diversity between different cytotypes in mountain and lowland areas. Our findings can also be applied to study conservation strategies for ensuring biodiversity at the genetic level in polyploid complexes.

Hybridization Between Yuccas From Baja California: Genomic and Environmental Patterns

Hybridization can occur when two geographically isolated species are reproductively compatible and have come into sympatry due to range shifts. Yucca and yucca moths exhibit obligate pollination mutualism; yucca moths are responsible for the gene flow mediated by pollen among yucca populations. In the Baja California Peninsula, there are two yucca sister species, Y. capensis and Y. valida, that have coevolved with the same pollinator, Tegeticula baja. Both yucca species are endemic to the peninsula, and their current distributions are allopatric. Based on their morphological characteristics, it has been suggested that some plants growing in the southern part of the Magdalena flatland, a spatially disjunct part of Yucca valida's range, have hybrid origins. We conducted genomic and climatic analyses of the two yucca species as well as the putative hybrid populations. We genotyped 3,423 single nucleotide polymorphisms in 120 individuals sampled from 35 localities. We applied Bayesian tests and geographic cline analyses to the genomic data. Using climatic information from the occurrence sites, we projected species distribution models in different periods to assess changes in the distributional range, and we performed a statistical test to define the niche divergence between the paternal species and the putative hybrid populations. Structure analysis revealed mixed ancestry in the genome of hybrid populations, and the Bayesian models supported a scenario of post-divergence gene flow between the yucca species. Our species distribution models reveal that the geographical ranges of the parental species overlapped mainly during the Last Glacial Maximum, which could facilitate genetic admixture between those species. Finally, we found that most of the assessed environmental axes between the parents and hybrid populations are divergent, indicating that the climatic niche of the hybrid populations is shifting from that of the populations' progenitors. Our results show that the populations in the southern part of the Magdalena flatland are the result of combination of the genetic components of two species. Hybrid individuals with this novel genomic combination arose in a different habitat than their parental species, and they exhibit ecological divergence, which contributes to reproductive isolation through spatial and temporal barriers.

Additive Traits Lead to Feeding Advantage and Reproductive Isolation, Promoting Homoploid Hybrid Speciation

Speciation through homoploid hybridization (HHS) is considered extremely rare in animals. This is mainly because the establishment of reproductive isolation as a product of hybridization is uncommon. Additionally, many traits are underpinned by polygeny and/or incomplete dominance, where the hybrid phenotype is an additive blend of parental characteristics. Phenotypically intermediate hybrids are usually at a fitness disadvantage compared with parental species and tend to vanish through backcrossing with parental population(s). It is therefore unknown whether the additive nature of hybrid traits in itself could lead successfully to HHS. Using a multi-marker genetic data set and a meta-analysis of diet and morphology, we investigated a potential case of HHS in the prions (Pachyptila spp.), seabirds distinguished by their bills, prey choice, and timing of breeding. Using approximate Bayesian computation, we show that the medium-billed Salvin's prion (Pachyptila salvini) could be a hybrid between the narrow-billed Antarctic prion (Pachyptila desolata) and broad-billed prion (Pachyptila vittata). Remarkably, P. salvini's intermediate bill width has given it a feeding advantage with respect to the other Pachyptila species, allowing it to consume a broader range of prey, potentially increasing its fitness. Available metadata showed that P. salvini is also intermediate in breeding phenology and, with no overlap in breeding times, it is effectively reproductively isolated from either parental species through allochrony. These results provide evidence for a case of HHS in nature, and show for the first time that additivity of divergent parental traits alone can lead directly to increased hybrid fitness and reproductive isolation.

'At the crossroads towards polyploidy': genomic divergence and extent of homoploid hybridization are drivers for the formation of the ox-eye daisy polyploid complex (Leucanthemum, Compositae-Anthemideae)

Polyploidy plays a paramount role in phytodiversity, but the causes of this evolutionary pathway require further study. Here, we use phylogenetic methods to examine possible polyploidy-promoting factors by comparing diploid representatives of the comprehensive European polyploid complex Leucanthemum with members of its strictly diploid North African counterpart Rhodanthemum. We investigate genetic divergence and gene flow among all diploid lineages of both genera to evaluate the role of genomic differentiation and hybridization for polyploid speciation. To test whether hybridization in Leucanthemum has been triggered by the geological conditions during its diversification, we additionally generate a time-calibrated phylogeny of 46 species of the subtribe Leucantheminae. Leucanthemum shows a significantly higher genetic divergence and hybridization signal among diploid lineages compared with Rhodanthemum, in spite of a similar crown age and diversification pattern during the Quaternary. Our study demonstrates the importance of genetic differentiation among diploid progenitors and their concurrent affinity for natural hybridization for the formation of a polyploid complex. Furthermore, the role of climate-induced range overlaps on hybridization and polyploid speciation during the Quaternary is discussed.

New prospects in the detection and comparative analysis of hybridization in the tree of life

Assessing the relative importance of the various pathways to diversification is a central goal of biodiversity researchers. For plant biologists, and increasingly across the spectrum of biological sciences, among these pathways of interest is hybridization. New methodological developments are moving the field away from questions of whether natural hybridization occurs or hybrids can persist and toward more direct assessments of the long-term impact of hybridization on diversification and genome organization. Advances in theory and new data, especially phylogenomic data, have changed the face of this field, revealing extensive occurrences of hybridization at both shallow and deep levels, but lacking is a synthesis of these advancements. Here we provide an overview of methods that have been proposed for detecting hybridization with molecular data and advocate a time-extended, comparative view of reticulate evolution. In particular, we pose three overarching questions, newly placed within reach, that are critical for advancing our understanding of hybridization pattern and process: (1) How often is introgression biased toward certain genomes and loci, and is this bias selectively neutral? (2) What are the relative rates of formation of hybrid species and introgressants, and how does this compare to their subsequent fates? (3) Has the frequency of hybridization increased under historical periods of greater dynamism in climate and geographic range, such as the Pleistocene?

Geography and Environment Shape Landscape Genetics of Mediterranean Alpine Species Silene ciliata Poiret. (Caryophyllaceae)

The study of the drivers that shape spatial genetic structure across heterogeneous landscapes is one of the main approaches used to understand population dynamics and responses in changing environments. While the Isolation-by-Distance model (IBD) assumes that genetic differentiation increases among populations with geographical distance, the Isolation-by-Resistance model (IBR) also considers geographical barriers and other landscape features that impede gene flow. On the other hand, the Isolation-by-Environment model (IBE) explains genetic differentiation through environmental differences between populations. Although spatial genetic studies have increased significantly in recent years, plants from alpine ecosystems are highly underrepresented, even though they are great suitable systems to disentangle the role of the different factors that structure genetic variation across environmental gradients. Here, we studied the spatial genetic structure of the Mediterranean alpine specialist Silene ciliata across its southernmost distribution limit. We sampled three populations across an altitudinal gradient from 1850 to 2400 m, and we replicated this sample over three mountain ranges aligned across an E-W axis in the central part of the Iberian Peninsula. We genotyped 20 individuals per population based on eight microsatellite markers and used different landscape genetic tools to infer the role of topographic and environmental factors in shaping observed patterns along the altitudinal gradient. We found a significant genetic structure among the studied Silene ciliata populations which was related to the orography and E-W configuration of the mountain ranges. IBD pattern arose as the main factor shaping population genetic differentiation. Geographical barriers between mountain ranges also affected the spatial genetic structure (IBR pattern). Although environmental variables had a significant effect on population genetic diversity parameters, no IBE pattern was found on genetic structure. Our study reveals that IBD was the driver that best explained the genetic structure, whereas environmental factors also played a role in determining genetic diversity values of this dominant plant of Mediterranean alpine environments.