Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey

Detailing the genetic basis of adaptive variation in natural populations is a first step towards understanding the process of adaptive evolution, yet few ecologically relevant traits have been characterized at the genetic level in wild populations. Traits that mediate coevolutionary interactions between species are ideal for studying adaptation because of the intensity of selection and the well-characterized ecological context. We have previously described the ecological context, evolutionary history and partial genetic basis of tetrodotoxin (TTX) resistance in garter snakes (Thamnophis). Derived mutations in a voltage-gated sodium channel gene (Na_v1.4) in three garter snake species are associated with resistance to TTX, the lethal neurotoxin found in their newt prey (Taricha). Here we evaluate the contribution of Na_v1.4 alleles to TTX resistance in two of those species from central coastal California. We measured the phenotypes (TTX resistance) and genotypes (Na_v1.4 and microsatellites) in a local sample of Thamnophis atratus and Thamnophis sirtalis. Allelic variation in Na_v1.4 explains 23 per cent of the variation in TTX resistance in T. atratus while variation in a haphazard sample of the genome (neutral microsatellite markers) shows no association with the phenotype. Similarly, allelic variation in Na_v1.4 correlates almost perfectly with TTX resistance in T. sirtalis, but neutral variation does not. These strong correlations suggest that Na_v1.4 is a major effect locus. The simple genetic architecture of TTX resistance in garter snakes may significantly impact the dynamics of phenotypic coevolution. Fixation of a few alleles of major effect in some garter snake populations may have led to the evolution of extreme phenotypes and an ‘escape’ from the arms race with newts.     

Independent FLC Mutations as Causes of Flowering-Time Variation in Arabidopsis thaliana and Capsella rubella

Capsella rubella is an inbreeding annual forb closely related to Arabidopsis thaliana, a model species widely used for studying natural variation in adaptive traits such as flowering time. Although mutations in dozens of genes can affect flowering of A. thaliana in the laboratory, only a handful of such genes vary in natural populations. Chief among these are FRIGIDA (FRI) and FLOWERING LOCUS C (FLC). Common and rare FRI mutations along with rare FLC mutations explain a large fraction of flowering-time variation in A. thaliana. Here we document flowering time under different conditions in 20 C. rubella accessions from across the species’ range. Similar to A. thaliana, vernalization, long photoperiods and elevated ambient temperature generally promote flowering. In this collection of C. rubella accessions, we did not find any obvious loss-of-function FRI alleles. Using mapping-by-sequencing with two strains that have contrasting flowering behaviors, we identified a splice-site mutation in FLC as the likely cause of early flowering in accession 1408. However, other similarly early C. rubella accessions did not share this mutation. We conclude that the genetic basis of flowering-time variation in C. rubella is complex, despite this very young species having undergone an extreme genetic bottleneck when it split from C. grandiflora a few tens of thousands of years ago.     

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait‐specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high‐throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy‐mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.     

Thyroid hormone responsive QTL and the evolution of paedomorphic salamanders

The transformation of ancestral phenotypes into novel traits is poorly understood for many examples of evolutionary novelty. Ancestrally, salamanders have a biphasic life cycle with an aquatic larval stage, a brief and pronounced metamorphosis, followed by a terrestrial adult stage. Repeatedly during evolution, metamorphic timing has been delayed to exploit growth-permissive environments, resulting in paedomorphic salamanders that retain larval traits as adults. We used thyroid hormone (TH) to rescue metamorphic phenotypes in paedomorphic salamanders and then identified quantitative trait loci (QTL) for life history traits that are associated with amphibian life cycle evolution: metamorphic timing and adult body size. We demonstrate that paedomorphic tiger salamanders (Ambystoma tigrinum complex) carry alleles at three moderate effect QTL (met1–3) that vary in responsiveness to TH and additively affect metamorphic timing. Salamanders that delay metamorphosis attain significantly larger body sizes as adults and met2 explains a significant portion of this variation. Thus, substitution of alleles at TH-responsive loci suggests an adaptive pleiotropic basis for two key life-history traits in amphibians: body size and metamorphic timing. Our study demonstrates a likely pathway for the evolution of novel paedomorphic species from metamorphic ancestors via selection of TH-response alleles that delay metamorphic timing and increase adult body size.     

Exogenous selection shapes germination behaviour and seedling traits of populations at different altitudes in a Senecio hybrid zone

Background and Aims

The Senecio hybrid zone on Mt Etna, Sicily, is characterized by steep altitudinal clines in quantitative traits and genetic variation. Such clines are thought to be maintained by a combination of ‘endogenous’ selection arising from genetic incompatibilities and environment-dependent ‘exogenous’ selection leading to local adaptation. Here, the hypothesis was tested that local adaptation to the altitudinal temperature gradient contributes to maintaining divergence between the parental species, S. chrysanthemifolius and S. aethnensis.

Methods

Intra- and inter-population crosses were performed between five populations from across the hybrid zone and the germination and early seedling growth of the progeny were assessed.

Key Results

Seedlings from higher-altitude populations germinated better under low temperatures (9–13 °C) than those from lower altitude populations. Seedlings from higher-altitude populations had lower survival rates under warm conditions (25/15 °C) than those from lower altitude populations, but also attained greater biomass. There was no altitudinal variation in growth or survival under cold conditions (15/5 °C). Population-level plasticity increased with altitude. Germination, growth and survival of natural hybrids and experimentally generated F₁s generally exceeded the worse-performing parent.

Conclusions

Limited evidence was found for endogenous selection against hybrids but relatively clear evidence was found for divergence in seed and seedling traits, which is probably adaptive. The combination of low-temperature germination and faster growth in warm conditions might enable high-altitude S. aethnensis to maximize its growth during a shorter growing season, while the slower growth of S. chrysanthemifolius may be an adaptation to drought stress at low altitudes. This study indicates that temperature gradients are likely to be an important environmental factor generating and maintaining adaptive divergence across the Senecio hybrid zone on Mt Etna.     

QTL analysis of flowering time inArabidopsis thaliana

Quantitative trait loci (QTL) analyses based on restriction fragment length polymorphism maps have been used to resolve the genetic control of flowering time in a cross between twoArabidopsis thaliana ecotypes H51 and Landsbergerecta, differing widely in flowering time. Five quantitative trait loci affecting flowering time were identified in this cross (RLN1-5), four of which are located in regions containing mutations or loci previously identified as conferring a late-flowering phenotype. One of these loci is coincident with theFRI locus identified as the major determinant for late flowering and vernalization responsiveness in theArabidopsis ecotype Stockholm.RLN5, which maps to the lower half of chromosome five (between markers mi69 and m233), only affected flowering time significantly under short day conditions following a vernalization period. The late-flowering phenotype of H51 compared to Landsbergerecta was due to alleles conferring late flowering at only two of the five loci. At the three other loci, H51 possessed alleles conferring early flowering in comparison to those of Landsbergerecta. Combinations of alleles conferring early and late flowering from both parents accounted for the transgressive segregation of flowering time observed within the F₂ population. Three QTL,RLN1,RLN2 andRLN3 displayed significant genotype-by-environment interactions for flowering time. A significant interaction between alleles atRLN3 andRLN4 was detected.     

Genetic Variation for Life History Sensitivity to Seasonal Warming in Arabidopsis thaliana

Climate change has altered life history events in many plant species; however, little is known about genetic variation underlying seasonal thermal response. In this study, we simulated current and three future warming climates and measured flowering time across a globally diverse set of Arabidopsis thaliana accessions. We found that increased diurnal and seasonal temperature (1°–3°) decreased flowering time in two fall cohorts. The early fall cohort was unique in that both rapid cycling and overwintering life history strategies were revealed; the proportion of rapid cycling plants increased by 3–7% for each 1° temperature increase. We performed genome-wide association studies (GWAS) to identify the underlying genetic basis of thermal sensitivity. GWAS identified five main-effect quantitative trait loci (QTL) controlling flowering time and another five QTL with thermal sensitivity. Candidate genes include known flowering loci; a cochaperone that interacts with heat-shock protein 90; and a flowering hormone, gibberellic acid, a biosynthetic enzyme. The identified genetic architecture allowed accurate prediction of flowering phenotypes (R² > 0.95) that has application for genomic selection of adaptive genotypes for future environments. This work may serve as a reference for breeding and conservation genetic studies under changing environments.     

Natural variation in the genes responsible for maturity loci E1, E2, E3 and E4 in soybean

Background and Aims

The timing of flowering has a direct impact on successful seed production in plants. Flowering of soybean (Glycine max) is controlled by several E loci, and previous studies identified the genes responsible for the flowering loci E1, E2, E3 and E4. However, natural variation in these genes has not been fully elucidated. The aims of this study were the identification of new alleles, establishment of allele diagnoses, examination of allelic combinations for adaptability, and analysis of the integrated effect of these loci on flowering.

Methods

The sequences of these genes and their flanking regions were determined for 39 accessions by primer walking. Systematic discrimination among alleles was performed using DNA markers. Genotypes at the E1–E4 loci were determined for 63 accessions covering several ecological types using DNA markers and sequencing, and flowering times of these accessions at three sowing times were recorded.

Key Results

A new allele with an insertion of a long interspersed nuclear element (LINE) at the promoter of the E1 locus (e1-re) was identified. Insertion and deletion of 36 bases in the eighth intron (E2-in and E2-dl) were observed at the E2 locus. Systematic discrimination among the alleles at the E1–E3 loci was achieved using PCR-based markers. Allelic combinations at the E1–E4 loci were found to be associated with ecological types, and about 62–66 % of variation of flowering time could be attributed to these loci.

Conclusions

The study advances understanding of the combined roles of the E1–E4 loci in flowering and geographic adaptation, and suggests the existence of unidentified genes for flowering in soybean,     

Adaptive basis of geographic variation: genetic,phenotypic and environmental differences among beach mouse populations

A major goal in evolutionary biology is to understand how and why populations differentiate, both genetically and phenotypically, as they invade a novel habitat. A classical example of adaptation is the pale colour of beach mice, relative to their dark mainland ancestors, which colonized the isolated sandy dunes and barrier islands on Florida''s Gulf Coast. However, much less is known about differentiation among the Gulf Coast beach mice, which comprise five subspecies linearly arrayed on Florida''s shoreline. Here, we test the role of selection in maintaining variation among these beach mouse subspecies at multiple levels—phenotype, genotype and the environments they inhabit. While all beach subspecies have light pelage, they differ significantly in colour pattern. These subspecies are also genetically distinct: pair-wise F_st-values range from 0.23 to 0.63 and levels of gene flow are low. However, we did not find a correlation between phenotypic and genetic distance. Instead, we find a significant association between the average ‘lightness’ of each subspecies and the brightness of the substrate it inhabits: the two most genetically divergent subspecies occupy the most similar habitats and have converged on phenotype, whereas the most genetically similar subspecies occupy the most different environments and have divergent phenotypes. Moreover, allelic variation at the pigmentation gene, Mc1r, is statistically correlated with these colour differences but not with variation at other genetic loci. Together, these results suggest that natural selection for camouflage—via changes in Mc1r allele frequency—contributes to pigment differentiation among beach mouse subspecies.     

Female and male genetic contributions to post-mating immune defence in female Drosophila melanogaster

Post-mating reduction in immune defence is common in female insects, and a trade-off between mating and immunity could affect the evolution of immunity. In this work, we tested the capacity of virgin and mated female Drosophila melanogaster to defend against infection by four bacterial pathogens. We found that female D. melanogaster suffer post-mating immunosuppression in a pathogen-dependent manner. The effect of mating was seen after infection with two bacterial pathogens (Providencia rettgeri and Providencia alcalifaciens), though not after infection with two other bacteria (Enterococcus faecalis and Pseudomonas aeruginosa). We then asked whether the evolution of post-mating immunosuppression is primarily a ‘female’ or ‘male’ trait by assaying for genetic variation among females for the degree of post-mating immune suppression they experience and among males for the level of post-mating immunosuppression they elicit in their mates. We also assayed for an interaction between male and female genotypes to test the specific hypothesis that the evolution of a trade-off between mating and immune defence in females might be being driven by sexual conflict. We found that females, but not males, harbour significant genetic variation for post-mating immunosuppression, and we did not detect an interaction between female and male genotypes. We thus conclude that post-mating immune depression is predominantly a ‘female’ trait, and find no evidence that it is evolving under sexual conflict.     

Effects of Infection by Belonolaimus longicaudatus on Rooting Dynamics among St. Augustinegrass and Bermudagrass Genotypes

Understanding rooting dynamics using the minirhizotron technique is useful for cultivar selection and to quantify nematode damage to roots. A 2-yr microplot study including five bermudagrass (‘Tifway’, Belonolaimus longicaudatus susceptible; two commercial cultivars [TifSport and Celebration] and two genotypes [‘BA132’ and ‘PI 291590’], which have been reported to be tolerant to B. longicaudatus) and two St. Augustinegrass (‘FX 313’, susceptible, and ‘Floratam’ that was reported as tolerant to B. longicaudatus) genotypes in a 5 x 2 and 2 x 2 factorial design with four replications, respectively, was initiated in 2012. Two treatments included were uninoculated and B. longicaudatus inoculated. In situ root images were captured each month using a minirhizotron camera system from April to September of 2013 and 2014. Mixed models analysis and comparison of least squares means indicated significant differences in root parameters studied across the genotypes and soil depths of both grass species. ‘Celebration’, ‘TifSport’ and ‘PI 291590’ bermudagrass, and ‘Floratam’ St. Augustinegrass had significantly different root parameters compared to the corresponding susceptible genotypes (P ≤ 0.05). Only ‘TifSport’ had no significant root loss when infested with B. longicaudatus compared to non-infested. ‘Celebration’ and ‘PI 291590’ had significant root loss but retained significantly greater root densities than ‘Tifway’ in B. longicaudatus-infested conditions (P ≤ 0.05). Root lengths were greater at the 0 to 5 cm depth followed by 5 to 10 and 10 to 15 cm of vertical soil depth for both grass species (P ≤ 0.05). ‘Celebration’, ‘TifSport’, and ‘PI 291590’ had better root vigor against B. longicaudatus compared to Tifway.     

Genome-Wide Association Analysis of Adaptation Using Environmentally Predicted Traits

Current methods for studying the genetic basis of adaptation evaluate genetic associations with ecologically relevant traits or single environmental variables, under the implicit assumption that natural selection imposes correlations between phenotypes, environments and genotypes. In practice, observed trait and environmental data are manifestations of unknown selective forces and are only indirectly associated with adaptive genetic variation. In theory, improved estimation of these forces could enable more powerful detection of loci under selection. Here we present an approach in which we approximate adaptive variation by modeling phenotypes as a function of the environment and using the predicted trait in multivariate and univariate genome-wide association analysis (GWAS). Based on computer simulations and published flowering time data from the model plant Arabidopsis thaliana, we find that environmentally predicted traits lead to higher recovery of functional loci in multivariate GWAS and are more strongly correlated to allele frequencies at adaptive loci than individual environmental variables. Our results provide an example of the use of environmental data to obtain independent and meaningful information on adaptive genetic variation.     

Genomic changes under rapid evolution: selection for parasitoid resistance

In this study, we characterize changes in the genome during a swift evolutionary adaptation, by combining experimental selection with high-throughput sequencing. We imposed strong experimental selection on an ecologically relevant trait, parasitoid resistance in Drosophila melanogaster against Asobara tabida. Replicated selection lines rapidly evolved towards enhanced immunity. Larval survival after parasitization increased twofold after just five generations of selection. Whole-genome sequencing revealed that the fast and strong selection response in innate immunity produced multiple, highly localized genomic changes. We identified narrow genomic regions carrying a significant signature of selection, which were present across all chromosomes and covered in total less than 5% of the whole D. melanogaster genome. We identified segregating sites with highly significant changes in frequency between control and selection lines that fell within these narrow ‘selected regions’. These segregating sites were associated with 42 genes that constitute possible targets of selection. A region on chromosome 2R was highly enriched in significant segregating sites and may be of major effect on parasitoid defence. The high genetic variability and small linkage blocks in our base population are likely responsible for allowing this complex trait to evolve without causing widespread erosive effects in the genome, even under such a fast and strong selective regime.     

Global depression in gene expression as a response to rapid thermal changes in vent mussels

Hydrothermal vent mussels belonging to the genus Bathymodiolus are distributed worldwide and dominate communities at shallow Atlantic hydrothermal sites. While organisms inhabiting coastal ecosystems are subjected to predictable oscillations of physical and chemical variables owing to tidal cycles, the vent mussels sustain pronounced temperature changes over short periods of time, correlated to the alternation of oxic/anoxic phases. In this context, we focused on the short-term adaptive response of mussels to temperature change at a molecular level. The mRNA expression of 23 genes involved in various cell functions of the vent mussel Bathymodiolus azoricus was followed after heat shocks for either 30 or 120 min, at 25 and 30°C over a 48 h recovery period at 5°C. Mussels were genotyped at 10 enzyme loci to explore a relationship between natural genetic variation, gene expression and temperature adaptation. Results indicate that the mussel response to increasing temperature is a depression in gene expression, such a response being genotypically correlated at least for the Pgm-1 locus. This suggests that an increase in temperature could be a signal triggering anaerobiosis for B. azoricus or this latter alternatively behaves more like a ‘cold’ stenotherm species, an attribute more related to its phylogenetic history, a cold seeps/wood fall origin.     

Seasonal variation in the mating system of a selfing annual with large floral displays

Background and Aims Flowering plants display considerable variation in mating system, specifically the relative frequency of cross- and self-fertilization. The majority of estimates of outcrossing rate do not account for temporal variation, particularly during the flowering season. Here, we investigated seasonal variation in mating and fertility in Incarvillea sinensis (Bignoniaceae), an annual with showy, insect-pollinated, ‘one-day’ flowers capable of delayed selfing. We examined the influence of several biotic and abiotic environmental factors on day-to-day variation in fruit set, seed set and patterns of mating.Methods We recorded daily flower number and pollinator abundance in nine 3 × 3-m patches in a population at Mu Us Sand land, Inner Mongolia, China. From marked flowers we collected data on daily fruit and seed set and estimated outcrossing rate and biparental inbreeding using six microsatellite loci and 172 open-pollinated families throughout the flowering period.Key Results Flower density increased significantly over most of the 50-d flowering season, but was associated with a decline in levels of pollinator service by bees, particularly on windy days. Fruit and seed set declined over time, especially during the latter third of the flowering period. Multilocus estimates of outcrossing rate were obtained using two methods (the programs MLTR and BORICE) and both indicated high selfing rates of ∼80 %. There was evidence for a significant increase in levels of selfing as the flowering season progressed and pollinator visitation declined. Biparental inbreeding also declined significantly as the flowering season progressed.Conclusions Temporal variation in outcrossing rates may be a common feature of the mating biology of annual, insect-pollinated plants of harsh environments but our study is the first to examine seasonal mating-system dynamics in this context. Despite having large flowers and showy floral displays, I. sinensis attracted relatively few pollinators. Delayed selfing by corolla dragging largely explains the occurrence of mixed mating in I. sinensis, and this mode of self-fertilization probably functions to promote reproductive assurance when pollinator service is limited by windy environmental conditions and competition from co-occurring flowering plants.     

QTL analysis of flowering time inArabidopsis thaliana

Quantitative trait loci (QTL) analyses based on restriction fragment length polymorphism maps have been used to resolve the genetic control of flowering time in a cross between twoArabidopsis thaliana ecotypes H51 and Landsbergerecta, differing widely in flowering time. Five quantitative trait loci affecting flowering time were identified in this cross (RLN1-5), four of which are located in regions containing mutations or loci previously identified as conferring a late-flowering phenotype. One of these loci is coincident with theFRI locus identified as the major determinant for late flowering and vernalization responsiveness in theArabidopsis ecotype Stockholm.RLN5, which maps to the lower half of chromosome five (between markers mi69 and m233), only affected flowering time significantly under short day conditions following a vernalization period. The late-flowering phenotype of H51 compared to Landsbergerecta was due to alleles conferring late flowering at only two of the five loci. At the three other loci, H51 possessed alleles conferring early flowering in comparison to those of Landsbergerecta. Combinations of alleles conferring early and late flowering from both parents accounted for the transgressive segregation of flowering time observed within the F₂ population. Three QTL,RLN1,RLN2 andRLN3 displayed significant genotype-by-environment interactions for flowering time. A significant interaction between alleles atRLN3 andRLN4 was detected.     

Tests for the joint evolution of mating system and drought escape in Mimulus

Background and Aims

Self-fertilizing taxa are often found at the range margins of their progenitors, where sub-optimal habitats may select for alternative physiological strategies. The extent to which self-fertilization is favoured directly vs. arising indirectly through correlations with other adaptive life history traits is unclear. Trait responses to selection depend on genetic variation and covariation, as well as phenotypic and genetic responses to altered environmental conditions. We tested predictions of the hypothesis that self-fertilization in Mimulus arises through direct selection on physiological and developmental traits that allow seasonal drought escape.

Methods

Phenotypic selection on mating system and drought escape traits was estimated in field populations of M. guttatus. In addition, trait phenotype and phenotypic selection were compared between experimental wet and dry soil in two greenhouse populations each of M. guttatus and M. nasutus. Finally, genetic variation and covariation for traits were compared between wet and dry soil treatments in a greenhouse population of M. guttatus.

Key Results

Consistent with predictions, selection for early flowering was generally stronger than for mating system traits, and selection for early flowering was stronger in dry soil. Inconsistent with predictions, selection for water-use efficiency was largely absent; selection for large flowers was stronger than for drought escape in the field; and most drought escape and mating system traits were not genetically correlated. A positive genetic correlation between flowering time and flower size, which opposed the adaptive contour, emerged only in wet soil, suggesting that variation in water availability may maintain variation in these traits. Plastic responses to soil moisture treatments supported the idea that taxonomic divergence could have been facilitated by plasticity in flowering time and selfing.

Conclusions

The hypothesis that plant mating systems may evolve indirectly via selection on correlated life history characteristics is plausible and warrants increased attention.     

Combination of Eight Alleles at Four Quantitative Trait Loci Determines Grain Length in Rice

Grain length is an important quantitative trait in rice (Oryza sativa L.) that influences both grain yield and exterior quality. Although many quantitative trait loci (QTLs) for grain length have been identified, it is still unclear how different alleles from different QTLs regulate grain length coordinately. To explore the mechanisms of QTL combination in the determination of grain length, five mapping populations, including two F₂ populations, an F₃ population, an F₇ recombinant inbred line (RIL) population, and an F₈ RIL population, were developed from the cross between the U.S. tropical japonica variety ‘Lemont’ and the Chinese indica variety ‘Yangdao 4’ and grown under different environmental conditions. Four QTLs (qGL-3-1, qGL-3-2, qGL-4, and qGL-7) for grain length were detected using both composite interval mapping and multiple interval mapping methods in the mapping populations. In each locus, there was an allele from one parent that increased grain length and another allele from another parent that decreased it. The eight alleles in the four QTLs were analyzed to determine whether these alleles act additively across loci, and lead to a linear relationship between the predicted breeding value of QTLs and phenotype. Linear regression analysis suggested that the combination of eight alleles determined grain length. Plants carrying more grain length-increasing alleles had longer grain length than those carrying more grain length-decreasing alleles. This trend was consistent in all five mapping populations and demonstrated the regulation of grain length by the four QTLs. Thus, these QTLs are ideal resources for modifying grain length in rice.     

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 F₂ 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.