Natural selection is an important driver of genetic and phenotypic differentiation between species. A powerful way to test the role of natural selection in the formation and maintenance of species is to study species complexes in which potential gene flow is high but realized gene flow is low. For a recent radiation of New World desert shrubs ( Encelia : Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of species hybridize at ecologically divergent parapatric boundaries in the adaptive radiation of Encelia , such mechanisms may maintain species boundaries throughout this group. SIGNIFICANCE STATEMENTIn Baja California, the deserts meet the coastal dunes in a narrow transition visible even from satellite images. We study two species pairs of desert shrubs (Encelia) that occur across this transition. Although these species can interbreed, they remain distinct. Using a combination of genetics, field experiments, 3D-imaging, and physiological measurements, we show that natural selection counteracts the homogenizing effects of gene exchange. The different habitats of these species create multiple mechanisms of selection -drought, salinity, herbivory, and burial, which together maintain these species in their native habitats and their hybrids in intermediate habitats. This study illustrates how Page 4 DiVittorio et al. Extremely strong natural selection environmental factors influence traits and fitness and how this in turn maintain species, highlighting the importance of natural selection in speciation.

There are multiple hypotheses for the spectacular plant diversity found in deserts. We explore how different factors, including the roles of ecological opportunity and selection, promote diversification and disparification in Encelia, a lineage of woody plants in the deserts of the Americas. Using a nearly complete species-level phylogeny based on double-digest restriction-aided sequencing along with a broad set of phenotypic traits, we estimate divergence times and diversification rates, identify instances of hybridization, quantify trait disparity and assess phenotypic divergence across environmental gradients. We show that Encelia originated and diversified recently (mid-Pleistocene) and rapidly, with rates comparable to notable adaptive radiations in plants. Encelia probably originated in the hot deserts of North America, with subsequent diversification across steep environmental gradients. We uncover multiple instances of gene flow between species. The radiation of Encelia is characterized by fast rates of phenotypic evolution, trait lability and extreme disparity across environments and between species pairs with overlapping geographic ranges. Encelia exemplifies how interspecific gene flow in combination with high trait lability can enable exceptionally fast diversification and disparification across steep environmental gradients.

SummaryThere are multiple hypotheses for the spectacular plant diversity found in deserts. We explore how different factors, including the roles of ecological opportunity and selection, promote diversification and disparification in Encelia, a lineage of woody plants in the deserts of the Americas.Using a nearly complete species-level phylogeny along with a broad set of phenotypic traits, we estimate divergence times and diversification rates, identify instances of hybridization, quantify trait disparity, and assess phenotypic divergence across environmental gradients.We show that Encelia originated and diversified recently (mid-Pleistocene) and rapidly, with rates comparable to notable adaptive radiations in plants. Encelia probably originated in the hot deserts of North America, with subsequent diversification across steep environmental gradients. We uncover multiple instances of gene flow between species. The radiation of Encelia is characterized by fast rates of phenotypic evolution, trait lability, and extreme disparity across environments and between species-pairs with overlapping geographic ranges.Encelia exemplifies how interspecific gene flow in combination with high trait lability can enable exceptionally fast diversification and disparification across steep environmental gradients.