Hybrid Speciation and Introgression Both Underlie the Genetic Structures and Evolutionary Relationships of Three Morphologically Distinct Species of (Liliaceae) Forming a Hybrid Zone Along an Elevational Gradient.

We studied hybrid interactions of , , and using an integrative approach combining population genetics, fieldwork, and phenological research. These three species occur along an elevational gradient, with occurring at lower elevations, at higher elevations, and between them. The species show strong morphological differentiation despite there being no clear environmental barriers to gene flow among them. is likely to have a hybrid origin based on our prior work, but its progenitors remain uncertain. We sought to determine whether gene flow occurs among these three parapatric species, and, if so, whether is a hybrid of and/or . We analyzed data from multiple chloroplast genes and spacers, nuclear internal transcribed spacer (ITS), and 18 nuclear Expressed Sequence Tag-Simple Sequence Repeat (EST-SSR) microsatellites for accessions of the three species representing dense population-level sampling. We also inferred phenology by examining species in the field and using herbarium specimens. We found that there are only two types of chloroplast genomes shared among the three species and that forms two distinct groups with closest links to other species of based on ITS. Taken together, is unlikely to be a hybrid species resulting from a cross between and , but gene flow is occurring among the three species. The gene flow is likely to be rare according to evidence from all molecular datasets, and this is corroborated by detection of only one putative hybrid individual in the field and asynchronous phenology. We suspect that the rarity of hybridization events among the species facilitates their continued genetic separation.