Absence of long-term balancing selection on variation in EuMYB3, an R2R3-MYB gene responsible for the anther-color polymorphism in Erythronium umbilicatum.

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

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

Complex cross-incompatibility in morning glories is consistent with a role for mating system in plant speciation.

Reproductive isolation between selfing and outcrossing species can arise through diverse mechanisms, some of which are directly associated with differences in mating system. We dissected cross-incompatibility between the highly selfing morning glory Ipomoea lacunosa and its mixed-mating sister species Ipomoea cordatotriloba. We found that cross-incompatibility is complex, with contributions acting both before and after fertilization. We then investigated whether the transition in mating system may have facilitated the evolution of these reproductive barrier components through mismatched floral morphology, differences in reproductive context, or both. We found evidence that morphological mismatch likely contributes to reproductive isolation in at least one cross-direction and that other pollen-pistil interactions are present. We also identified hybrid seed inviability consistent with the predictions of the weak-inbreeder, strong-outbreeder hypotheses, suggesting endosperm misregulation plays an important role in cross-incompatibility. In contrast, we did not find evidence consistent with the prezygotic weak-inbreeder, strong-outbreeder hypothesis. Our study highlights the complexity of reproductive isolation between outcrossing and selfing species and the extent to which evolutionary consequences of mating system transitions can facilitate speciation.

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation.

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Floral scent divergence across an elevational hybrid zone with varying pollinators.

Divergence in floral traits attractive to different pollinators can promote reproductive isolation in related species. When isolation is incomplete, hybridization may occur, which offers the opportunity to explore mechanisms underlying reproductive isolation. Recent work suggests that divergence in floral scent may frequently contribute to reproductive barriers, although such divergence has seldom been examined in species with generalized pollination. Here, we used two closely related Penstemon species, P. newberryi and P. davidsonii, and their natural hybrids from an elevational gradient with pollinator communities that are predicted to vary in their reliance on floral scent (i.e., primarily hummingbirds at low elevation vs. bees at high elevation). The species vary in a suite of floral traits, but scent is uncharacterized. To address whether scent varies along elevation and potentially contributes to reproductive isolation, we genetically characterized individuals collected at field and identified whether they were parental species or hybrids. We then characterized scent amount and composition. Although the parental species had similar total emissions, some scent characteristics (i.e., scent composition, aromatic emission) diverged between them and may contribute to their isolation. However, the species emitted similar compound sets which could explain hybridization in the contact area. Hybrids were similar to the parents for most scent traits, suggesting that their floral scent would not provide a strong barrier to backcrossing. Our study suggests floral scent may be a trait contributing to species boundaries even in plants with generalized pollination, and reinforces the idea that evolutionary pollinator transitions may involve changes in multiple floral traits.

Morphological and molecular characterization of variation in common bean (Phaseolus vulgaris L.) germplasm from Azad Jammu and Kashmir, Pakistan.

Phaseolus vulgaris, an essential food and source of protein, is cultivated across the world. This study was carried out to investigate the diversity and population structure of 34 P. vulgaris landrace accessions collected from the Azad Jammu and Kashmir (AJ&K) regions of Pakistan. The samples were analyzed both morphologically and using genetic variation identified through RNA sequencing. Our results indicated that most genetic variation occurs among local accessions, with little genetic variation occurring between geographical regions. In addition, the accessions fell into two major genetic groups. Morphological analysis revealed that these two genetic groups differ in a number of quantitative traits, including seed length, seed width, and seed weight. One accession, DUD-11, appears to be a mixture of the two major groups genetically as well as morphologically. Among the other accessions, DUD-8, RWK-2, and NGD-1 depicted particularly high seed weight along with higher seed length, seed width, and seed yield per plant. We suggest focusing on these accessions in future breeding programs. More generally, our results provide baseline data that will be useful for crop improvement and effective cultivation practices in Pakistan.

Modularity and selection of nectar traits in the evolution of the selfing syndrome in Ipomoea lacunosa (Convolvulaceae).

Although the evolution of the selfing syndrome often involves reductions in floral size, pollen and nectar, few studies of selfing syndrome divergence have examined nectar. We investigate whether nectar traits have evolved independently of other floral size traits in the selfing syndrome, whether nectar traits diverged due to drift or selection, and the extent to which quantitative trait locus (QTL) analyses predict genetic correlations. We use F5 recombinant inbred lines (RILs) generated from a cross between Ipomoea cordatotriloba and Ipomoea lacunosa. We calculate genetic correlations to identify evolutionary modules, test whether trait divergence was due to selection, identify QTLs and perform correlation analyses to evaluate how well QTL properties reflect genetic correlations. Nectar and floral size traits form separate evolutionary modules. Selection has acted to reduce nectar traits in the selfing I. lacunosa. Genetic correlations predicted from QTL properties are consistent with observed genetic correlations. Changes in floral traits associated with the selfing syndrome reflect independent evolution of at least two evolutionary modules: nectar and floral size traits. We also demonstrate directional selection on nectar traits, which is likely to be independent of selection on floral size traits. Our study also supports the expected mechanistic link between QTL properties and genetic correlations.

Genetic architecture of divergence: the selfing syndrome in Ipomoea lacunosa.

PREMISE: Highly selfing plant species frequently display a distinctive suite of traits termed the selfing syndrome. Here we tested the hypothesis that these traits are grouped into correlated evolutionary modules and determined the degree of independence between such modules. METHODS: We evaluated phenotypic correlations and QTL overlaps in F2 offspring of a cross between the morning glories Ipomoea lacunosa and I. cordatotriloba and investigated how traits clustered into modules at both the phenotypic and genetic level. We then compared our findings to other QTL studies of the selfing syndrome. RESULTS: In the I. lacunosa selfing syndrome, traits grouped into modules that displayed correlated evolution within but not between modules. QTL overlap predicted phenotypic correlations, and QTLs affecting the same trait module were significantly physically clustered in the genome. The genetic architecture of the selfing syndrome varied across systems, but the pattern of stronger within- than between-module correlation was widespread. CONCLUSIONS: The genetic architecture we observe in the selfing syndrome is consistent with a growing understanding of floral morphological integration achieved via pleiotropy in clustered traits. This view of floral evolution is consistent with resource limitation or predation driving the evolution of the selfing syndrome, but invites further research into both the selective causes of the selfing syndrome and how genetic architecture itself evolves in response to changes in mating system.

Ancient Gene Duplications, Rather Than Polyploidization, Facilitate Diversification of Petal Pigmentation Patterns in Clarkia gracilis (Onagraceae).

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. Clarkia gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, that is, before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.

Morning glory species co-occurrence is associated with asymmetrically decreased and cascading reproductive isolation.

Hybridization between species can affect the strength of the reproductive barriers that separate those species. Two extensions of this effect are (1) the expectation that asymmetric hybridization or gene flow will have asymmetric effects on reproductive barrier strength and (2) the expectation that local hybridization will affect only local reproductive barrier strength and could therefore alter within-species compatibility. We tested these hypotheses in a pair of morning glory species that exhibit asymmetric gene flow from highly selfing Ipomoea lacunosa into mixed-mating Ipomoea cordatotriloba in regions where they co-occur. Because of the direction of this gene flow, we predicted that reproductive barrier strength would be more strongly affected in I. cordatotriloba than I. lacunosa. We also predicted that changes to reproductive barriers in sympatric I. cordatotriloba populations would affect compatibility with allopatric populations of that species. We tested these predictions by measuring the strength of a reproductive barrier to seed set across the species' ranges. Consistent with our first prediction, we found that sympatric and allopatric I. lacunosa produce the same number of seeds in crosses with I. cordatotriloba, whereas crosses between sympatric I. cordatotriloba and I. lacunosa are more successful than crosses between allopatric I. cordatotriloba and I. lacunosa. This difference in compatibility appears to reflect an asymmetric decrease in the strength of the barrier to seed set in sympatric I. cordatotriloba, which could be caused by I. lacunosa alleles that have introgressed into I. cordatotriloba. We further demonstrated that changes to sympatric I. cordatotriloba have decreased its ability to produce seeds with allopatric populations of the same species, in line with our second prediction. Thus, in a manner analogous to cascade reinforcement, we suggest that introgression associated with hybridization not only influences between-species isolation but can also contribute to isolation within a species.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae).

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F2 plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

Host Plant Choice Determined by Reproductive Interference between Closely Related Butterflies.

AbstractA number of empirical studies have concluded that reproductive interference (RI) contributes to parapatric species distributions or sexual exclusion. However, the possibility that divergent host plant use in phytophagous insects is due to sexual exclusion has seldom been considered. Here, we present evidence that RI is responsible for different host plant use by two Pierid butterfly species, Pieris napi and Pieris melete. When a novel host species was introduced about 50 years ago, two Pierid butterfly species at first used both the ancestral host species and the novel one. Subsequently, P. napi shifted to use only the novel host, while P. melete shifted to specialize on the ancestral host. To explain these patterns, we investigated whether the two host species differ in suitability for larval growth and survival. Additionally, we tested whether RI occurred between the two butterfly species using large outdoor field cages. Courtship of females by conspecific and heterospecific males reduces the number of eggs laid by approximately half. However, RI is asymmetric and would generate selection on P. melete females to evolve to avoid the more suitable host species preferred by P. napi. Thus, our study suggests that sexual exclusion can explain the shift in host plant use by these two butterfly species.

Selection favors loss of floral pigmentation in a highly selfing morning glory.

A common evolutionary trend in highly selfing plants is the evolution of the "selfing syndrome", in which traits associated with pollinator attraction are lost or greatly reduced. Limited information is available on whether these trait reductions are favored by natural selection or result from reduced purifying selection coupled with genetic drift. This study attempted to distinguish between these two possibilities for the evolutionary loss of floral pigmentation in the highly selfing species Ipomoea lacunosa. This study also tested the hypothesis that loss of floral pigmentation is caused by downregulation or loss of function in a tissue-specific anthocyanin transcription factor, as has been found in other plants. F2 individuals of a cross between white and pigmented individuals revealed segregation at two epistatically acting loci: one affecting pigmentation in both corolla throat and limbs (Anl1) and one affecting limb pigmentation (Anl2). Individuals that are homozygous for the "white" allele at Anl1 have white throats and limbs regardless of genotype at Anl2. In individuals with pigmented throats, homozygosity of the "white" allele at Anl2 produces white limbs. Flower color variation at Anl1 cosegregates with an R2R3-Myb anthocyanin transcription factor, which is down-regulated in white-flowers but not in pigmented flowers. Differential expression of the two alleles of this gene indicates that down regulation is caused by a cis-regulatory change. Finally, allele-frequency differences at Anl1 were substantially and significantly greater than differences in allele frequencies at four microsatellite loci. These results are consistent with the hypotheses that the identified R2R3-Myb gene corresponds to Anl1 and that evolutionary loss of pigmentation in I. lacunosa was caused by selection. They are also consistent with previous studies demonstrating that loss of floral pigmentation is usually caused by down-regulation or functional inactivation of an R2R3-Myb gene.

A taxonomic monograph of Ipomoea integrated across phylogenetic scales.

Taxonomic monographs have the potential to make a unique contribution to the understanding of global biodiversity. However, such studies, now rare, are often considered too daunting to undertake within a realistic time frame, especially as the world's collections have doubled in size in recent times. Here, we report a global-scale monographic study of morning glories (Ipomoea) that integrated DNA barcodes and high-throughput sequencing with the morphological study of herbarium specimens. Our approach overhauled the taxonomy of this megadiverse group, described 63 new species and uncovered significant increases in net diversification rates comparable to the most iconic evolutionary radiations in the plant kingdom. Finally, we show that more than 60 species of Ipomoea, including sweet potato, independently evolved storage roots in pre-human times, indicating that the storage root is not solely a product of human domestication but a trait that predisposed the species for cultivation. This study demonstrates how the world's natural history collections can contribute to global challenges in the Anthropocene.

Adaptation to hummingbird pollination is associated with reduced diversification in Penstemon.

A striking characteristic of the Western North American flora is the repeated evolution of hummingbird pollination from insect-pollinated ancestors. This pattern has received extensive attention as an opportunity to study repeated trait evolution as well as potential constraints on evolutionary reversibility, with little attention focused on the impact of these transitions on species diversification rates. Yet traits conferring adaptation to divergent pollinators potentially impact speciation and extinction rates, because pollinators facilitate plant reproduction and specify mating patterns between flowering plants. Here, we examine macroevolutionary processes affecting floral pollination syndrome diversity in the largest North American genus of flowering plants, Penstemon. Within Penstemon, transitions from ancestral bee-adapted flowers to hummingbird-adapted flowers have frequently occurred, although hummingbird-adapted species are rare overall within the genus. We inferred macroevolutionary transition and state-dependent diversification rates and found that transitions from ancestral bee-adapted flowers to hummingbird-adapted flowers are associated with reduced net diversification rate, a finding based on an estimated 17 origins of hummingbird pollination in our sample. Although this finding is congruent with hypotheses that hummingbird adaptation in North American Flora is associated with reduced species diversification rates, it contrasts with studies of neotropical plant families where hummingbird pollination has been associated with increased species diversification. We further used the estimated macroevolutionary rates to predict the expected pattern of floral diversity within Penstemon over time, assuming stable diversification and transition rates. Under these assumptions, we find that hummingbird-adapted species are expected to remain rare due to their reduced diversification rates. In fact, current floral diversity in the sampled Penstemon lineage, where less than one-fifth of species are hummingbird adapted, is consistent with predicted levels of diversity under stable macroevolutionary rates.

Multiple aspects of the selfing syndrome of the morning glory Ipomoea lacunosa evolved in response to selection: A Qst-Fst comparison.

The frequent transition from outcrossing to selfing in flowering plants is often accompanied by changes in multiple aspects of floral morphology, termed the "selfing syndrome." While the repeated evolution of these changes suggests a role for natural selection, genetic drift may also be responsible. To determine whether selection or drift shaped different aspects of the pollination syndrome and mating system in the highly selfing morning glory Ipomoea lacunosa, we performed multivariate and univariate Qst-Fst comparisons using a wide sample of populations of I. lacunosa and its mixed-mating sister species Ipomoea cordatotriloba. The two species differ in early growth, floral display, inflorescence traits, corolla size, nectar, and pollen number. Our analyses support a role for natural selection driving trait divergence, specifically in corolla size and nectar traits, but not in early growth, display size, inflorescence length, or pollen traits. We also find evidence of selection for reduced herkogamy in I. lacunosa, consistent with selection driving both the transition in mating system and the correlated floral changes. Our research demonstrates that while some aspects of the selfing syndrome evolved in response to selection, others likely evolved due to drift or correlated selection, and the balance between these forces may vary across selfing species.

Gene flow, divergent selection and resistance to introgression in two species of morning glories (Ipomoea).

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.

Genetic architecture of quantitative flower and leaf traits in a pair of sympatric sister species of Primulina.

Flowers and leaves each represent suites of functionally interrelated traits that are often involved in species divergence and local adaptation. However, a major unresolved issue is how the individual component traits that make up a complex trait such as a flower evolve in a coordinated fashion to retain a high degree of functionality. We use a quantitative trait loci (QTL) approach to elucidate the genetic architecture of divergence in flower and leaf traits between the sister species Primulina depressa and Primulina danxiaensis, which grow sympatrically but in contrasting microhabitats. We found that flower traits were controlled by multiple QTL of small effect, while leaf physiological and morphological traits tended to be controlled by QTL of larger effect. The observed floral integration, manifested by a high degree overlap in both individual trait QTL and QTL for principal component scores (PCA QTL), may have been critical for evolutionary divergence of floral morphology in relation to their pollinators. This overlap suggests that direct selection on only one or a few of the component traits could have caused substantial divergence in other floral traits due to genetic correlations, while the low QTL overlap between floral and vegetative traits suggests that these trait suites are genetically unlinked and can evolve independently in response to different selective pressures corresponding to their distinct functions.

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.