Whole-genome sequence and annotation of Penstemon davidsonii.

Penstemon is the most speciose flowering plant genus endemic to North America. Penstemon species' diverse morphology and adaptation to various environments have made them a valuable model system for studying evolution. Here, we report the first full reference genome assembly and annotation for Penstemon davidsonii. Using PacBio long-read sequencing and Hi-C scaffolding technology, we constructed a de novo reference genome of 437,568,744 bases, with a contig N50 of 40 Mb and L50 of 5. The annotation includes 18,199 gene models, and both the genome and transcriptome assembly contain over 95% complete eudicot BUSCOs. This genome assembly will serve as a valuable reference for studying the evolutionary history and genetic diversity of the Penstemon genus.

A few genetic variants go a long way in differentiating Penstemon species.

The integrity of hybridizing species is usually maintained by genome-wide selection or by selection on a few genomic regions. A study published in PLOS Biology finds a different pattern-60 SNPs spread across the genome differentiate a Penstemon species pair.

A few essential genetic loci distinguish Penstemon species with flowers adapted to pollination by bees or hummingbirds.

In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genus Penstemon displays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from a Penstemon species complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-wide dXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 "species-diagnostic loci," which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.

Improving phylogenetic resolution of the Lamiales using the complete plastome sequences of six Penstemon species.

The North American endemic genus Penstemon (Mitchell) has a recent geologic origin of ca. 3.6 million years ago (MYA) during the Pliocene/Pleistocene transition and has undergone a rapid adaptive evolutionary radiation with ca. 285 species of perennial forbs and sub-shrubs. Penstemon is divided into six subgenera occupying all North American habitats including the Arctic tundra, Central American tropical forests, alpine meadows, arid deserts, and temperate grasslands. Due to the rapid rate of diversification and speciation, previous phylogenetic studies using individual and concatenated chloroplast sequences have failed to resolve many polytomic clades. We investigated the efficacy of utilizing the plastid genomes (plastomes) of 29 species in the Lamiales order, including five newly sequenced Penstemon plastomes, for analyzing phylogenetic relationships and resolving problematic clades. We compared whole-plastome based phylogenies to phylogenies based on individual gene sequences (matK, ndhF, psaA, psbA, rbcL, rpoC2, and rps2) and concatenated sequences. We also We found that our whole-plastome based phylogeny had higher nodal support than all other phylogenies, which suggests that it provides greater accuracy in describing the hierarchal relationships among taxa as compared to other methods. We found that the genus Penstemon forms a monophyletic clade sister to, but separate from, the Old World taxa of the Plantaginaceae family included in our study. Our whole-plastome based phylogeny also supports the rearrangement of the Scrophulariaceae family and improves resolution of major clades and genera of the Lamiales.

SNP-skimming: A fast approach to map loci generating quantitative variation in natural populations.

Genome-wide association mapping (GWAS) is a method to estimate the contribution of segregating genetic loci to trait variation. A major challenge for applying GWAS to nonmodel species has been generating dense genome-wide markers that satisfy the key requirement that marker data are error-free. Here, we present an approach to map loci within natural populations using inexpensive shallow genome sequencing. This "SNP-skimming" approach involves two steps: an initial genome-wide scan to identify putative targets followed by deep sequencing for confirmation of targeted loci. We apply our method to a test data set of floral dimension variation in the plant Penstemon virgatus, a member of a genus that has experienced dynamic floral adaptation that reflects repeated transitions in primary pollinator. The ability to detect SNPs that generate phenotypic variation depends on population genetic factors such as population allele frequency, effect size and epistasis, as well as sampling effects contingent on missing data and genotype uncertainty. However, both simulations and the Penstemon data suggest that the most significant tests from the initial SNP skim are likely to be true positives-loci with subtle but significant quantitative effects on phenotype. We discuss the promise and limitations of this method and consider optimal experimental design for a given sequencing effort. Simulations demonstrate that sampling a larger number of individual at the expense of average read depth per individual maximizes the power to detect loci.

Multiplexed shotgun genotyping resolves species relationships within the North American genus Penstemon.

PREMISE OF THE STUDY: Evolutionary radiations provide opportunities to examine large-scale patterns in diversification and character evolution, yet are often recalcitrant to phylogenetic resolution due to rapid speciation events. The plant genus Penstemon has been difficult to resolve using Sanger sequence-based markers, leading to the hypothesis that it represents a recent North American radiation. The current study demonstrates the utility of multiplexed shotgun genotyping (MSG), a style of restriction site-associated DNA sequencing (RADseq), to infer phylogenetic relationships within a subset of species in this genus and provide insight into evolutionary patterns. METHODS: We sampled genomic DNA, primarily from herbarium material, and subjected it to MSG library preparation and Illumina sequencing. The resultant sequencing reads were clustered into homologous loci, aligned, and concatenated into data matrices that differed according to clustering similarity and amount of missing data. We performed phylogenetic analyses on these matrices using maximum likelihood (RAxML) and a species tree approach (SVDquartets). KEY RESULTS: MSG data provide a highly resolved estimate of species relationships within Penstemon. While most species relationships were highly supported, the position of certain taxa remains ambiguous, suggesting that increased taxonomic sampling or additional methodologies may be required. The data confirm that evolutionary shifts from hymenopteran- to hummingbird-adapted flowers have occurred independently many times. CONCLUSIONS: This study demonstrates that phylogenomic approaches yielding thousands of variable sites can greatly improve species-level resolution of recent and rapid radiations. Similar to other studies, we found that less conservative similarity and missing data thresholds resulted in more highly supported topologies.

Predicting the genetic consequences of future climate change: The power of coupling spatial demography, the coalescent, and historical landscape changes.

PREMISE OF THE STUDY: Climate change is a widely accepted threat to biodiversity. Species distribution models (SDMs) are used to forecast whether and how species distributions may track these changes. Yet, SDMs generally fail to account for genetic and demographic processes, limiting population-level inferences. We still do not understand how predicted environmental shifts will impact the spatial distribution of genetic diversity within taxa. METHODS: We propose a novel method that predicts spatially explicit genetic and demographic landscapes of populations under future climatic conditions. We use carefully parameterized SDMs as estimates of the spatial distribution of suitable habitats and landscape dispersal permeability under present-day, past, and future conditions. We use empirical genetic data and approximate Bayesian computation to estimate unknown demographic parameters. Finally, we employ these parameters to simulate realistic and complex models of responses to future environmental shifts. We contrast parameterized models under current and future landscapes to quantify the expected magnitude of change. KEY RESULTS: We implement this framework on neutral genetic data available from Penstemon deustus. Our results predict that future climate change will result in geographically widespread declines in genetic diversity in this species. The extent of reduction will heavily depend on the continuity of population networks and deme sizes. CONCLUSIONS: To our knowledge, this is the first study to provide spatially explicit predictions of within-species genetic diversity using climatic, demographic, and genetic data. Our approach accounts for climatic, geographic, and biological complexity. This framework is promising for understanding evolutionary consequences of climate change, and guiding conservation planning.

Ecological transition predictably associated with gene degeneration.

Gene degeneration or loss can significantly contribute to phenotypic diversification, but may generate genetic constraints on future evolutionary trajectories, potentially restricting phenotypic reversal. Such constraints may manifest as directional evolutionary trends when parallel phenotypic shifts consistently involve gene degeneration or loss. Here, we demonstrate that widespread parallel evolution in Penstemon from blue to red flowers predictably involves the functional inactivation and degeneration of the enzyme flavonoid 3',5'-hydroxylase (F3'5'H), an anthocyanin pathway enzyme required for the production of blue floral pigments. Other types of genetic mutations do not consistently accompany this phenotypic shift. This pattern may be driven by the relatively large mutational target size of degenerative mutations to this locus and the apparent lack of associated pleiotropic effects. The consistent degeneration of F3'5'H may provide a mechanistic explanation for the observed asymmetry in the direction of flower color evolution in Penstemon: Blue to red transitions are common, but reverse transitions have not been observed. Although phenotypic shifts in this system are likely driven by natural selection, internal constraints may generate predictable genetic outcomes and may restrict future evolutionary trajectories.

Identification of major quantitative trait loci underlying floral pollination syndrome divergence in Penstemon.

Distinct floral pollination syndromes have emerged multiple times during the diversification of flowering plants. For example, in western North America, a hummingbird pollination syndrome has evolved more than 100 times, generally from within insect-pollinated lineages. The hummingbird syndrome is characterized by a suite of floral traits that attracts and facilitates pollen movement by hummingbirds, while at the same time discourages bee visitation. These floral traits generally include large nectar volume, red flower colour, elongated and narrow corolla tubes and reproductive organs that are exerted from the corolla. A handful of studies have examined the genetic architecture of hummingbird pollination syndrome evolution. These studies find that mutations of relatively large effect often explain increased nectar volume and transition to red flower colour. In addition, they suggest that adaptive suites of floral traits may often exhibit a high degree of genetic linkage, which could facilitate their fixation during pollination syndrome evolution. Here, we explore these emerging generalities by investigating the genetic basis of floral pollination syndrome divergence between two related Penstemon species with different pollination syndromes--bee-pollinated P. neomexicanus and closely related hummingbird-pollinated P. barbatus. In an F2 mapping population derived from a cross between these two species, we characterized the effect size of genetic loci underlying floral trait divergence associated with the transition to bird pollination, as well as correlation structure of floral trait variation. We find the effect sizes of quantitative trait loci for adaptive floral traits are in line with patterns observed in previous studies, and find strong evidence that suites of floral traits are genetically linked. This linkage may be due to genetic proximity or pleiotropic effects of single causative loci. Interestingly, our data suggest that the evolution of floral traits critical for hummingbird pollination was not constrained by negative pleiotropy at loci that show co-localization for multiple traits.

Predictability and irreversibility of genetic changes associated with flower color evolution in Penstemon barbatus.

Two outstanding questions in evolutionary biology are whether, and how often, the genetic basis of phenotypic evolution is predictable; and whether genetic change constrains evolutionary reversibility. We address these questions by studying the genetic basis of red flower color in Penstemon barbatus. The production of red flowers often involves the inactivation of one or both of two anthocyanin pathway genes, Flavonoid 3',5'-hydroxylase (F3'5'h) and Flavonoid 3'-hydroxylase (F3'h). We used gene expression and enzyme function assays to determine that redundant inactivating mutations to F3'5'h underlie the evolution of red flowers in P. barbatus. Comparison of our results to previously characterized shifts from blue to red flowers suggests that the genetic change associated with the evolution of red flowers is predictable: when it involves elimination of F3'5'H activity, functional inactivation or deletion of this gene tends to occur; however, when it involves elimination of F3'H activity, tissue-specific regulatory substitutions occur and the gene is not functionally inactivated. This pattern is consistent with emerging data from physiological experiments indicating that F3'h may have pleiotropic effects and is thus subject to purifying selection. The multiple, redundant inactivating mutations to F3'5'h suggest that reversal to blue-purple flowers in this group would be unlikely.

Developing molecular tools and insights into the Penstemon genome using genomic reduction and next-generation sequencing.

BACKGROUND: Penstemon's unique phenotypic diversity, hardiness, and drought-tolerance give it great potential for the xeric landscaping industry. Molecular markers will accelerate the breeding and domestication of drought tolerant Penstemon cultivars by, creating genetic maps, and clarifying of phylogenetic relationships. Our objectives were to identify and validate interspecific molecular markers from four diverse Penstemon species in order to gain specific insights into the Penstemon genome. RESULTS: We used a 454 pyrosequencing and GR-RSC (genome reduction using restriction site conservation) to identify homologous loci across four Penstemon species (P. cyananthus, P. davidsonii, P. dissectus, and P. fruticosus) representing three diverse subgenera with considerable genome size variation. From these genomic data, we identified 133 unique interspecific markers containing SSRs and INDELs of which 51 produced viable PCR-based markers. These markers produced simple banding patterns in 90% of the species × marker interactions (~84% were polymorphic). Twelve of the markers were tested across 93, mostly xeric, Penstemon taxa (72 species), of which ~98% produced reproducible marker data. Additionally, we identified an average of one SNP per 2,890 bp per species and one per 97 bp between any two apparent homologous sequences from the four source species. We selected 192 homologous sequences, meeting stringent parameters, to create SNP markers. Of these, 75 demonstrated repeatable polymorphic marker functionality across the four sequence source species. Finally, sequence analysis indicated that repetitive elements were approximately 70% more prevalent in the P. cyananthus genome, the largest genome in the study, than in the smallest genome surveyed (P. dissectus). CONCLUSIONS: We demonstrated the utility of GR-RSC to identify homologous loci across related Penstemon taxa. Though PCR primer regions were conserved across a broadly sampled survey of Penstemon species (93 taxa), DNA sequence within these amplicons (12 SSR/INDEL markers) was highly diverse. With the continued decline in next-generation sequencing costs, it will soon be feasible to use genomic reduction techniques to simultaneously sequence thousands of homologous loci across dozens of Penstemon species. Such efforts will greatly facilitate our understanding of the phylogenetic structure within this important drought tolerant genus. In the interim, this study identified thousands of SNPs and over 50 SSRs/INDELs which should provide a foundation for future Penstemon phylogenetic studies and breeding efforts.

Natural selection on floral volatile production in Penstemon digitalis: highlighting the role of linalool.

Natural selection is thought to have shaped the evolution of floral scent; however, unlike other floral characters, we have a rudimentary knowledge of how phenotypic selection acts on scent. We found that floral scent was under stronger selection than corolla traits such as flower size and flower color in weakly scented Penstemon digitalis. Our results suggest that to understand evolution in floral phenotypes, including scent in floral selection, studies are crucial. For P. digitalis, linalool was the direct target of selection in the scent bouquet. Therefore, we determined the enantiomeric configuration of linalool because interacting insects may perceive the enantiomers differentially. We found that P. digitalis produces only (S)-(+)-linalool and, more interestingly, it is also taken up into the nectar. Because the nectar is scented and flavored with (S)-(+)-linalool, it may be an important cue for pollinators visiting P. digitalis flowers.

Influences of landscape and pollinators on population genetic structure: examples from three Penstemon (Plantaginaceae) species in the Great Basin.

PREMISE OF THE STUDY: Despite rapid growth in the field of landscape genetics, our understanding of how landscape features interact with life history traits to influence population genetic structure in plant species remains limited. Here, we identify population genetic divergence in three species of Penstemon (Plantaginaceae) similarly distributed throughout the Great Basin region of the western United States but with different pollination syndromes (bee and hummingbird). The Great Basin's mountainous landscape provides an ideal setting to compare the interaction of landscape and dispersal ability in isolating populations of different species. METHODS: We used eight highly polymorphic microsatellite loci to identify neutral population genetic structure between populations within and among mountain ranges for eight populations of P. deustus, 10 populations of P. pachyphyllus, and 10 populations of P. rostriflorus. We applied traditional population genetics approaches as well as spatial and landscape genetics approaches to infer genetic structure and discontinuities among populations. KEY RESULTS: All three species had significant genetic structure and exhibited isolation by distance, ranging from high structure and low inferred gene flow in the bee-pollinated species P. deustus (F(ST) = 0.1330, R(ST) = 0.4076, seven genetic clusters identified) and P. pachyphyllus (F(ST) = 0.1896, R(ST) = 0.2531, four genetic clusters identified) to much lower structure and higher inferred gene flow in the hummingbird-pollinated P. rostriflorus (F(ST) = 0.0638, R(ST) = 0.1116, three genetic clusters identified). CONCLUSIONS: These three Penstemon species have significant yet strikingly different patterns of population genetic structure, findings consistent with different interactions between landscape features and the dispersal capabilities of their pollinators.

A survey of Penstemon's genome size.

Penstemon is the largest genus in North America with more than 270 reported species. However, little is known about its genome size. This information may be useful in developing hybrids for landscape use and for gaining insight into its current taxonomy. Using flow cytometry, we estimated the genome size of approximately 40% of the genus (115 accessions from 105 different species). Genome sizes for both reported and probable diploids range from P. dissectus 2C = 0.94 pg (1C = 462 Mbp) to P. pachyphyllus var. mucronatus 2C = 1.88 pg (1C = 919 Mbp), and the polyploids range from P. attenuatus var. attenuatus 2C = 2.35 pg (1C = 1148 Mbp) to P. digitalis 2C = 6.45 pg (1C = 3152 Mbp). Chromosome counts were done for ten previously published and four previously unreported Penstemon species (P. dissectus, P. navajoa, P. caespitosus var. desertipicti, and P. ramaleyi). These counts were compiled with all previously published chromosome data and compared with the flow cytometry results. Ploidy within this study ranged from diploid to dodecaploid. These data were compared and contrasted with the current taxonomy of Penstemon and previously published internal transcribed spacer and chloroplast DNA phylogenetic work. Based on genome size and previous studies, reassigning P. montanus to the subgenus Penstemon and P. personatus to the subgenus Dasanthera, would better reflect the phylogeny of the genus. Furthermore, our data concur with previous studies suggesting that the subgenus Habroanthus be included in the subgenus Penstemon. The DNA content of subgenus Penstemon exhibits high plasticity and spans a sixfold increase from the smallest to the largest genome (P. linarioides subsp. sileri and P. digitalis, respectively). Our study found flow cytometry to be useful in species identification and verification.

Pollinator-mediated natural selection in Penstemon digitalis.

Measuring the agents of natural selection is important because it allows us to understand not only which traits are expected to evolve but also why they will evolve. Natural selection by pollinators on floral traits is often assumed because in outcrossing animal-pollinated plants flowers are generally thought to function as advertisements of rewards directed at pollinators. We tested the role of bee pollinators in selection on Penstemon digitalis and found that pollinators were driving selection for larger and more flowers. However, what makes our publication unique is the additional information we gained from reviewing the few other studies that also directly tested whether pollinators were agents of selection on floral traits. As we would expect if pollinators are important agents of selection, selection on floral traits was significantly stronger when pollinators were present than when they were experimentally removed. Taken together, these results suggest that pollinators can be important drivers of selection in contemporary populations.

Physiological differences among two Penstemon species and their hybrids in field and common garden environments.

Hybrids can exhibit unique combinations of the physiological traits of their parents. These particular combinations may influence hybrid fitness and the evolutionary trajectory of a hybrid zone. Here, a hybrid zone between Penstemon newberryi and Penstemon davidsonii along an elevational gradient was examined, and physiological traits of parents and hybrids were measured in their native environment and a common garden. Gas exchange rates of nine different crosses were also measured. Alpine P. davidsonii had less negative pre-dawn water potential and lower water use efficiency (WUE) than its montane relative P. newberryi in a common garden and in field measurements. The species difference in WUE was attributable to lower conductance in P. newberryi in the field, but to a higher photosynthetic rate in this species in the common garden. The alpine species took less time to produce mature fruits and reached maximum photosynthetic rate at a lower temperature. Natural hybrids were intermediate for most characters. F(1) hybrids had lower conductance than progeny of natural hybrids. The intermediate WUE of natural hybrids may be one factor that allows them to persist in intermediate environments. Comparisons of different crosses suggest that the genotypic composition of hybrids influences their physiological performance.

Constrained lability in floral evolution: counting convergent origins of hummingbird pollination in Penstemon and Keckiella.

In the clade of Penstemon and segregate genera, pollination syndromes are well defined among the 284 species. Most display combinations of floral characters associated with pollination by Hymenoptera, the ancestral mode of pollination for this clade. Forty-one species present characters associated with hummingbird pollination, although some of these ornithophiles are also visited by insects. The ornithophiles are scattered throughout the traditional taxonomy and across phylogenies estimated from nuclear (internal transcribed spacer (ITS)) and chloroplast DNA (trnCD/TL) sequence data. Here, the number of separate origins of ornithophily is estimated, using bootstrap phylogenies and constrained parsimony searches. Analyses suggest 21 separate origins, with overwhelming support for 10 of these. Because species sampling was incomplete, this is probably an underestimate. Penstemons therefore show great evolutionary lability with respect to acquiring hummingbird pollination; this syndrome acts as an attractor to which species with large sympetalous nectar-rich flowers have frequently been drawn. By contrast, penstemons have not undergone evolutionary shifts backwards or to other pollination syndromes. Thus, they are an example of both striking evolutionary lability and constrained evolution.