Tight genetic linkage of prezygotic barrier loci creates a multifunctional speciation island in Petunia

The genetic architecture of floral trait divergence between hummingbird- and self-pollinated monkeyflower (Mimulus) species

Pollination syndromes are a key component of flowering plant diversification, prompting questions about the architecture of single traits and genetic coordination among traits. Here, we investigate the genetics of extreme floral divergence between naturally hybridizing monkeyflowers, Mimulus parishii (self-pollinated) and M. cardinalis (hummingbird-pollinated). We mapped quantitative trait loci (QTLs) for 18 pigment, pollinator reward/handling, and dimensional traits in parallel sets of F2 hybrids plus recombinant inbred lines and generated nearly isogenic lines (NILs) for two dimensional traits, pistil length and corolla size. Our multi-population approach revealed a highly polygenic basis (n = 190 QTLs total) for pollination syndrome divergence, capturing minor QTLs even for pigment traits with leading major loci. There was significant QTL overlap within pigment and dimensional categories. Nectar volume QTLs clustered with those for floral dimensions, suggesting a partially shared module. The NILs refined two pistil length QTLs, only one of which has tightly correlated effects on other dimensional traits. An overall polygenic architecture of floral divergence is partially coordinated by genetic modules formed by linkage (pigments) and likely pleiotropy (dimensions plus nectar). This work illuminates pollinator syndrome diversification in a model radiation and generates a robust framework for molecular and ecological genomics.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Rapid genomic evolution in Brassica rapa with bumblebee selection in experimental evolution

BACKGROUND

Insect pollinators shape rapid phenotypic evolution of traits related to floral attractiveness and plant reproductive success. However, the underlying genomic changes remain largely unknown despite their importance in predicting adaptive responses to natural or to artificial selection. Based on a nine-generation experimental evolution study with fast cycling Brassica rapa plants adapting to bumblebees, we investigate the genomic evolution associated with the previously observed parallel phenotypic evolution. In this current evolve and resequencing (E&R) study, we conduct a genomic scan of the allele frequency changes along the genome in bumblebee-pollinated and hand-pollinated plants and perform a genomic principal component analysis (PCA).

RESULTS

We highlight rapid genomic evolution associated with the observed phenotypic evolution mediated by bumblebees. Controlling for genetic drift, we observe significant changes in allelic frequencies at multiple loci. However, this pattern differs according to the replicate of bumblebee-pollinated plants, suggesting putative non-parallel genomic evolution. Finally, our study underlines an increase in genomic variance implying the putative involvement of multiple loci in short-term pollinator adaptation.

CONCLUSIONS

Overall, our study enhances our understanding of the complex interactions between pollinator and plants, providing a stepping stone towards unravelling the genetic basis of plant genomic adaptation to biotic factors in the environment.

How the switch to hummingbird pollination has greatly contributed to our understanding of evolutionary processes

The evolutionary switch to hummingbird pollination exemplifies complex adaptation, requiring evolutionary change in multiple component traits. Despite this complexity, diverse lineages have converged on hummingbird-adapted flowers on a relatively short evolutionary timescale. Here, I review how features of the genetic basis of adaptation contribute to this remarkable evolutionary lability. Large-effect substitutions, large mutational targets for adaptation, adaptive introgression, and concentrated architecture all contribute to the origin and maintenance of hummingbird-adapted flowers. The genetic features of adaptation are likely shaped by the ecological and geographic context of the switch to hummingbird pollination, with implications for future evolutionary trajectories.

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.

Genetic architecture of a pollinator shift and its fate in secondary hybrid zones of two Petunia species

BACKGROUND

Theory suggests that the genetic architecture of traits under divergent natural selection influences how easily reproductive barriers evolve and are maintained between species. Divergently selected traits with a simple genetic architecture (few loci with major phenotypic effects) should facilitate the establishment and maintenance of reproductive isolation between species that are still connected by some gene flow. While empirical support for this idea appears to be mixed, most studies test the influence of trait architectures on reproductive isolation only indirectly. Petunia plant species are, in part, reproductively isolated by their different pollinators. To investigate the genetic causes and consequences of this ecological isolation, we deciphered the genetic architecture of three floral pollination syndrome traits in naturally occurring hybrids between the widespread Petunia axillaris and the highly endemic and endangered P. exserta.

RESULTS

Using population genetics, Bayesian linear mixed modelling and genome-wide association studies, we found that the three pollination syndrome traits vary in genetic architecture. Few genome regions explain a majority of the variation in flavonol content (defining UV floral colour) and strongly predict the trait value in hybrids irrespective of interspecific admixture in the rest of their genomes. In contrast, variation in pistil exsertion and anthocyanin content (defining visible floral colour) is controlled by many genome-wide loci. Opposite to flavonol content, the genome-wide proportion of admixture between the two species predicts trait values in their hybrids. Finally, the genome regions strongly associated with the traits do not show extreme divergence between individuals representing the two species, suggesting that divergent selection on these genome regions is relatively weak within their contact zones.

CONCLUSIONS

Among the traits analysed, those with a more complex genetic architecture are best maintained in association with the species upon their secondary contact. We propose that this maintained genotype-phenotype association is a coincidental consequence of the complex genetic architectures of these traits: some of their many underlying small-effect loci are likely to be coincidentally linked with the actual barrier loci keeping these species partially isolated upon secondary contact. Hence, the genetic architecture of a trait seems to matter for the outcome of hybridization not only then when the trait itself is under selection.

Tight genetic linkage of genes causing hybrid necrosis and pollinator isolation between young species

The mechanisms of reproductive isolation that cause phenotypic diversification and eventually speciation are a major topic of evolutionary research. Hybrid necrosis is a post-zygotic isolation mechanism in which cell death develops in the absence of pathogens. It is often due to the incompatibility between proteins from two parents. Here we describe a unique case of hybrid necrosis due to an incompatibility between loci on chromosomes 2 and 7 between two pollinator-isolated Petunia species. Typical immune responses as well as endoplasmic reticulum stress responses are induced in the necrotic line. The locus on chromosome 2 encodes ChiA1, a bifunctional GH18 chitinase/lysozyme. The enzymatic activity of ChiA1 is dispensable for the development of necrosis. We propose that the extremely high expression of ChiA1 involves a positive feedback loop between the loci on chromosomes 2 and 7. ChiA1 is tightly linked to major genes involved in the adaptation to different pollinators, a form of pre-zygotic isolation. This linkage of pre- and post-zygotic barriers strengthens reproductive isolation and probably contributes to rapid diversification and speciation.

Single gene mutation in a plant MYB transcription factor causes a major shift in pollinator preference

Understanding the molecular basis of reproductive isolation and speciation is a key goal of evolutionary genetics. In the South American genus Petunia, the R2R3-MYB transcription factor MYB-FL regulates the biosynthesis of UV-absorbing flavonol pigments, a major determinant of pollinator preference. MYB-FL is highly expressed in the hawkmoth-pollinated P. axillaris, but independent losses of its activity in sister taxa P. secreta and P. exserta led to UV-reflective flowers and associated pollinator shifts in each lineage (bees and hummingbirds, respectively). We created a myb-fl CRISPR mutant in P. axillaris and studied the effect of this single gene on innate pollinator preference. The mutation strongly reduced the expression of the two key flavonol-related biosynthetic genes but only affected the expression of few other genes. The mutant flowers were UV reflective as expected but additionally contained low levels of visible anthocyanin pigments. Hawkmoths strongly preferred the wild-type P. axillaris over the myb-fl mutant, whereas both social and solitary bee preference depended on the level of visible color of the mutants. MYB-FL, with its specific expression pattern, small number of target genes, and key position at the nexus of flavonol and anthocyanin biosynthetic pathways, provides a striking example of evolution by single mutations of large phenotypic effect.

Convergence without divergence in North American red-flowering Silene

Combinations of correlated floral traits have arisen repeatedly across angiosperms through convergent evolution in response to pollinator selection to optimize reproduction. While some plant groups exhibit very distinct combinations of traits adapted to specific pollinators (so-called pollination syndromes), others do not. Determining how floral traits diverge across clades and whether floral traits show predictable correlations in diverse groups of flowering plants is key to determining the extent to which pollinator-mediated selection drives diversification. The North American Silene section Physolychnis is an ideal group to investigate patterns of floral evolution because it is characterized by the evolution of novel red floral color, extensive floral morphological variation, polyploidy, and exposure to a novel group of pollinators (hummingbirds). We test for correlated patterns of trait evolution that would be consistent with convergent responses to selection in the key floral traits of color and morphology. We also consider both the role of phylogenic distance and geographic overlap in explaining patterns of floral trait variation. Inconsistent with phenotypically divergent pollination syndromes, we find very little clustering of North American Silene into distinct floral morphospace. We also find little evidence that phylogenetic history or geographic overlap explains patterns of floral diversity in this group. White- and pink-flowering species show extensive phenotypic diversity but are entirely overlapping in morphological variation. However, red-flowering species have much less phenotypic disparity and cluster tightly in floral morphospace. We find that red-flowering species have evolved floral traits that align with a traditional hummingbird syndrome, but that these trait values overlap with several white and pink species as well. Our findings support the hypothesis that convergent evolution does not always proceed through comparative phenotypic divergence, but possibly through sorting of standing ancestral variation.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

Testing candidate genes linked to corolla shape variation of a pollinator shift in Rhytidophyllum (Gesneriaceae)

Floral adaptations to specific pollinators like corolla shape variation often result in reproductive isolation and thus speciation. But despite their ecological importance, the genetic bases of corolla shape transitions are still poorly understood, especially outside model species. Hence, our goal was to identify candidate genes potentially involved in corolla shape variation between two closely related species of the Rhytidophyllum genus (Gesneriaceae family) from the Antilles with contrasting pollination strategies. Rhytidophyllum rupincola has a tubular corolla and is strictly pollinated by hummingbirds, whereas R. auriculatum has more open flowers and is pollinated by hummingbirds, bats, and insects. We surveyed the literature and used a comparative transcriptome sequence analysis of synonymous and non-synonymous nucleotide substitutions to obtain a list of genes that could explain floral variation between R. auriculatum and R. rupincola. We then tested their association with corolla shape variation using QTL mapping in a F2 hybrid population. Out of 28 genes tested, three were found to be good candidates because of a strong association with corolla shape: RADIALIS, GLOBOSA, and JAGGED. Although the role of these genes in Rhytidophyllum corolla shape variation remains to be confirmed, these findings are a first step towards identifying the genes that have been under selection by pollinators and thus involved in reproductive isolation and speciation in this genus.

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait's genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor "concentrated" genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.

The genetic basis of floral mechanical isolation between two hummingbird-pollinated Neotropical understorey herbs

Floral divergence can contribute to reproductive isolation among plant lineages, and thus provides an opportunity to study the genetics of speciation, including the number, effect size, mode of action and interactions of quantitative trait loci (QTL). Moreover, flowers represent suites of functionally interrelated traits, but it is unclear to what extent the phenotypic integration of the flower is underlain by a shared genetic architecture, which could facilitate or constrain correlated evolution of floral traits. Here, we examine the genetic architecture of floral morphological traits involved in an evolutionary switch from bill to forehead pollen placement between two species of hummingbird-pollinated Neotropical understorey herbs that are reproductively isolated by these floral differences. For the majority of traits, we find multiple QTL of relatively small effect spread throughout the genome. We also find substantial colocalization and alignment of effects of QTL underlying different floral traits that function together to promote outcrossing and reduce heterospecific pollen transfer. Our results are consistent with adaptive pleiotropy or linkage of many co-adapted genes, either of which could have facilitated a response to correlated selection and helped to stabilize divergent phenotypes in the face of low levels of hybridization. Moreover, our results indicate that floral mechanical isolation can be consistent with an infinitesimal model of adaptation.

Complex evolution of novel red floral color in Petunia

Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.

Do floral and ecogeographic isolation allow the co-occurrence of two ecotypes of Anacamptis papilionacea (Orchidaceae)?

Ecotypes are relatively frequent in flowering plants and considered central in ecological speciation as local adaptation can promote the insurgence of reproductive isolation. Without geographic isolation, gene flow usually homogenizes the allopatrically generated phenotypic and ecological divergences, unless other forms of reproductive isolation keep them separated. Here, we investigated two orchid ecotypes with marked phenotypic floral divergence that coexist in contact zones. We found that the two ecotypes show different ecological habitat preferences with one being more climatically restricted than the other. The ecotypes remain clearly morphologically differentiated both in allopatry and in sympatry and differed in diverse floral traits. Despite only slightly different flowering times, the two ecotypes achieved floral isolation thanks to different pollination strategies. We found that both ecotypes attract a wide range of insects, but the ratio of male/female attracted by the two ecotypes was significantly different, with one ecotype mainly attracts male pollinators, while the other mainly attracts female pollinators. As a potential consequence, the two ecotypes show different pollen transfer efficiency. Experimental plots with pollen staining showed a higher proportion of intra- than interecotype movements confirming floral isolation between ecotypes in sympatry while crossing experiments excluded evident postmating barriers. Even if not completely halting the interecotypes pollen flow in sympatry, such incipient switch in pollination strategy between ecotypes may represent a first step on the path toward evolution of sexual mimicry in Orchidinae.

Combining biotechnology and evolution for understanding the mechanisms of pollinator attraction

Many flowering plants rely on pollinators for their reproductive success. Plant-pollinator interactions usually depend on a complex combination of traits based on a fine-tuned biosynthetic machinery, with many structural and regulatory genes involved. Yet, the physiological mechanisms in plants are the product of evolutionary processes. While evolution has been modifying flowers through millions of years, it is also a rapid process that can change plant traits within few generations. Here we discuss both mechanistic and evolutionary aspects of pollinator attraction. We also propose how latest advances in biotechnology and evolutionary studies, and their combination, will improve the elucidation of molecular mechanisms and evolutionary dynamics of pollinator attraction in changing environments.

Rational Design of a Novel Hawkmoth Pollinator Interaction in Mimulus Section Erythranthe

Diversification of theca. 275,000 extant flowering plant species has been driven in large part by coevolution with animal pollinators. A recurring pattern of pollinator shifts from hummingbird to hawkmoth pollination has characterized plant speciation in many western North American plant taxa, but in the genusMimulus(monkeyflowers) sectionErythranthethe evolution of hawkmoth pollination from hummingbird-pollinated ancestors has not occurred. We manipulated two flower color loci and tested the attractiveness of the resulting four color phenotypes (red, yellow, pink, and white) to naïve hawkmoths (Manduca sexta). Hawkmoths strongly prefer derived colors (yellow, pink, white) over the ancestral red when choosing an initial flower to visit, and generally preferred derived colors when total visits and total visit time were considered, with no hawkmoth preferring ancestral red over derived colors. The simple flower color genetics underlying this innate pollinator preference suggests a potential path for speciation into an unfilled hawkmoth-pollinated niche inMimulussectionErythranthe, and the deliberate design of a hawkmoth-pollinated flower demonstrates a new, predictive method for studying pollination syndrome evolution.

Pollination syndromes in the 21st century: where do we stand and where may we go?

Pollination syndromes, recurring suites of floral traits appearing in connection with specific functional pollinator groups, have served for decades to organise floral diversity under a functional-ecological perspective. Some potential caveats, such as over-simplification of complex plant-animal interactions or lack of empirical observations, have been identified and discussed in recent years. Which of these caveats do indeed cause problems, which have been solved and where do future possibilities lie? I address these questions in a review of the pollination-syndrome literature of 2010 to 2019. I show that the majority of studies was based on detailed empirical pollinator observations and could reliably predict pollinators based on a few floral traits such as colour, shape or reward. Some traits (i.e. colour) were less reliable in predicting pollinators than others (i.e. reward, corolla width), however. I stress that future studies should consider floral traits beyond those traditionally recorded to expand our understanding of mechanisms of floral evolution. I discuss statistical methods suitable for objectively analysing the interplay of system-specific evolutionary constraints, pollinator-mediated selection and adaptive trade-offs at microecological and macroecological scales. I exemplify my arguments on an empirical dataset of floral traits of a neotropical plant radiation in the family Melastomataceae.

Dissecting genetic diversity and genomic background of Petunia cultivars with contrasting growth habits

The cultivated petunia (Petunia ×hybrida) is derived from the progenitor species P. axillaris and P. integrifolia. The hybridization dates back only to the 1830s, though intensive breeding efforts have yielded cultivars exhibiting incredible diversity for many traits, including growth habit, flower color, and flower size. Until now, little is known about the genetic diversity and genomic background of modern cultivars. Here we selected a panel of 13 cultivars with contrasting growth habits and three wild species (the progenitors and P. exserta) to estimate the genomic contribution from the ancestral species and to study whether the variation of the genetic origin could be associated with different breeding programs or morphological variability. Transcriptome sequencing identified 1,164,566 SNPs representing 98.4% (32,451) of the transcripts that cover 99.2% (of 52,697,361 bp) of the P. axillaris transcriptome. Cultivars with an upright growth habit had more homozygous alleles and more P. axillaris-derived alleles than trailing cultivars, while mounded cultivars had intermediate heterozygosity. Unlike previous studies, we found the proportions of alleles derived from each progenitor species varied across cultivars but overall were not biased toward one progenitor species, suggesting diverse selection during cultivar development. For trailing cultivars, alleles potentially introgressed from other wild species ("out" alleles) were enriched. The "out" alleles were clustered in particular regions of chromosomes, suggesting that these regions may be hotspots of introgression. Transcripts in these regions were enriched with gene ontology terms associated with growth habit. This study provides novel insight into the contributions of progenitor species to the genomic background of modern petunia cultivars and identifies genome regions that may harbor genes conferring the trailing growth habit for further exploration.

Identification of transcription factors controlling floral morphology in wild Petunia species with contrasting pollination syndromes

Adaptation to different pollinators is an important driver of speciation in the angiosperms. Genetic approaches such as QTL mapping have been successfully used to identify the underlying speciation genes. However, these methods are often limited by widespread suppression of recombination due to divergence between species. While the mutations that caused the interspecific differences in floral color and scent have been elucidated in a variety of plant genera, the genes that are responsible for morphological differences remain mostly unknown. Differences in floral organ length determine the pollination efficiency of hawkmoths and hummingbirds, and therefore the genes that control these differences are potential speciation genes. Identifying such genes is challenging, especially in non-model species and when studying complex traits for which little prior genetic and biochemical knowledge is available. Here we combine transcriptomics with detailed growth analysis to identify candidate transcription factors underlying interspecific variation in the styles of Petunia flowers. Starting from a set of 2284 genes, stepwise filtering for expression in styles, differential expression between species, correlation with growth-related traits, allele-specific expression in interspecific hybrids, and/or high-impact polymorphisms resulted in a set of 43 candidate speciation genes. Validation by virus-induced gene silencing identified two MYB transcription factors, EOBI and EOBII, that were previously shown to regulate floral scent emission, a trait associated with pollination by hawkmoths.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

What could be the fate of secondary contact zones between closely related plant species?

Interspecific hybridization has been fundamental in plant evolution. Nevertheless, the fate of hybrid zones throughout the generations remains poorly addressed. We analyzed a pair of recently diverged, interfertile, and sympatric Petunia species to ask what fate the interspecific hybrid population has met over time. We analyzed the genetic diversity in two generations from two contact sites and evaluated the effect of introgression. To do this, we collected all adult plants from the contact zones, including canonicals and intermediary colored individuals, and compared them with purebred representatives of both species based on seven highly informative microsatellite loci. We compared the genetic diversity observed in the contact zones with what is seen in isolated populations of each species, considering two generations of these annual species. Our results have confirmed the genetic differentiation between the species and the hybrid origin of the majority of the intermediary colored individuals. We also observed a differentiation related to genetic variability and inbreeding levels among the populations. Over time, there were no significant differences per site related to genetic diversity or phenotype composition. We found two stable populations kept by high inbreeding and backcross rates that influence the genetic diversity of their parental species through introgression.

Pleiotropy facilitates local adaptation to distant optima in common ragweed (Ambrosia artemisiifolia)

Pleiotropy, the control of multiple phenotypes by a single locus, is expected to slow the rate of adaptation by increasing the chance that beneficial alleles also have deleterious effects. However, a prediction arising from classical theory of quantitative trait evolution states that pleiotropic alleles may have a selective advantage when phenotypes are distant from their selective optima. We examine the role of pleiotropy in regulating adaptive differentiation among populations of common ragweed (Ambrosia artemisiifolia); a species that has recently expanded its North American range due to human-mediated habitat change. We employ a phenotype-free approach by using connectivity in gene networks as a proxy for pleiotropy. First, we identify loci bearing footprints of local adaptation, and then use genotype-expression mapping and co-expression networks to infer the connectivity of the genes. Our results indicate that the putatively adaptive loci are highly pleiotropic, as they are more likely than expected to affect the expression of other genes, and they reside in central positions within the gene networks. We propose that the conditionally advantageous alleles at these loci avoid the cost of pleiotropy by having large phenotypic effects that are beneficial when populations are far from their selective optima. We further use evolutionary simulations to show that these patterns are in agreement with a model where populations face novel selective pressures, as expected during a range expansion. Overall, our results suggest that highly connected genes may be targets of positive selection during environmental change, even though they likely experience strong purifying selection in stable selective environments.

Geography is essential for reproductive isolation between florally diversified morning glory species from Amazon canga savannahs

The variety, relative importance and eco-evolutionary stability of reproductive barriers are critical to understanding the processes of speciation and species persistence. Here we evaluated the strength of the biotic prezygotic and postzygotic isolation barriers between closely related morning glory species from Amazon canga savannahs. The flower geometry and flower visitor assemblage analyses supported pollination by the bees in lavender-flowered Ipomoea marabaensis and recruitment of hummingbirds as pollinators in red-flowered Ipomoea cavalcantei. Nevertheless, native bee species and alien honeybees foraged on flowers of both species. Real-time interspecific hybridization underscored functionality of the overlap in flower visitor assemblages, questioning the strength of prezygotic isolation underpinned by diversification in flower colour and geometry. Interspecific hybrids were fertile and produced offspring in nature. No significant asymmetry in interspecific hybridization and hybrid incompatibilities among offspring were found, indicating weak postmating and postzygotic isolation. The results suggested that despite floral diversification, the insular-type geographic isolation remains a major barrier to gene flow. Findings set a framework for the future analysis of contemporary evolution of plant-pollinator networks at the population, community, and ecosystem levels in tropical ecosystems that are known to be distinct from the more familiar temperate climate models.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata

Background

Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata.

Results

We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns.

Conclusions

This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Electronic supplementary material

The online version of this article (10.1186/s13227-019-0129-2) contains supplementary material, which is available to authorized users.

A simple genetic architecture and low constraint allow rapid floral evolution in a diverse and recently radiating plant genus

Genetic correlations among different components of phenotypes, especially those resulting from pleiotropy, can constrain or facilitate trait evolution. These factors could especially influence the evolution of traits that are functionally integrated, such as those comprising the flower. Indeed, pleiotropy is proposed as a main driver of repeated convergent trait transitions, including the evolution of phenotypically similar pollinator syndromes. We assessed the role of pleiotropy in the differentiation of floral and other reproductive traits between two species - Jaltomata sinuosa and J. umbellata (Solanaceae) - that have divergent suites of floral traits consistent with bee and hummingbird pollination, respectively. To do so, we generated a hybrid population and examined the genetic architecture (trait segregation and quantitative trait locus (QTL) distribution) underlying 25 floral and fertility traits. We found that most floral traits had a relatively simple genetic basis (few, predominantly additive, QTLs of moderate to large effect), as well as little evidence of antagonistic pleiotropy (few trait correlations and QTL colocalization, particularly between traits of different classes). However, we did detect a potential case of adaptive pleiotropy among floral size and nectar traits. These mechanisms may have facilitated the rapid floral trait evolution observed within Jaltomata, and may be a common component of rapid phenotypic change more broadly.

Repurposing population genetics data to discern genomic architecture: A case study of linkage cohort detection in mountain pine beetle (Dendroctonus ponderosae)

Genetic surveys of the population structure of species can be used as resources for exploring their genomic architecture. By adjusting filtering assumptions, genome-wide single-nucleotide polymorphism (SNP) datasets can be reused to give new insights into the genetic basis of divergence and speciation without targeted resampling of specimens. Filtering only for missing data and minor allele frequency, we used a combination of principal components analysis and linkage disequilibrium network analysis to distinguish three cohorts of variable SNPs in the mountain pine beetle in western Canada, including one that was sex-linked and one that was geographically associated. These marker cohorts indicate genomically localized differentiation, and their detection demonstrates an accessible and intuitive method for discovering potential islands of genomic divergence without a priori knowledge of a species' genomic architecture. Thus, this method has utility for directly addressing the genomic architecture of species and generating new hypotheses for functional research.

The Genetics of a Behavioral Speciation Phenotype in an Island System

Mating behavior divergence can make significant contributions to reproductive isolation and speciation in various biogeographic contexts. However, whether the genetic architecture underlying mating behavior divergence is related to the biogeographic history and the tempo and mode of speciation remains poorly understood. Here, we use quantitative trait locus (QTL) mapping to infer the number, distribution, and effect size of mating song rhythm variations in the crickets Laupala eukolea and Laupala cerasina, which occur on different islands (Maui and Hawaii). We then compare these results with a similar study of an independently evolving species pair that diverged within the same island. Finally, we annotate the L. cerasina transcriptome and test whether the QTL fall in functionally enriched genomic regions. We document a polygenic architecture behind the song rhythm divergence in the inter-island species pair that is remarkably similar to that previously found for an intra-island species pair in the same genus. Importantly, the QTL regions were significantly enriched for potential homologs of the genes involved in pathways that may be modulating the cricket song rhythm. These clusters of loci could constrain the spatial genomic distribution of the genetic variation underlying the cricket song variation and harbor several candidate genes that merit further study.

Varied Genomic Responses to Maladaptive Gene Flow and Their Evidence

Adaptation to a local environment often occurs in the face of maladaptive gene flow. In this perspective, I discuss several ideas on how a genome may respond to maladaptive gene flow during adaptation. On the one hand, selection can build clusters of locally adaptive alleles at fortuitously co-localized loci within a genome, thereby facilitating local adaptation with gene flow ('allele-only clustering'). On the other hand, the selective pressure to link adaptive alleles may drive co-localization of the actual loci relevant for local adaptation within a genome through structural genome changes or an evolving intra-genomic crossover rate ('locus clustering'). While the expected outcome is, in both cases, a higher frequency of locally adaptive alleles in some genome regions than others, the molecular units evolving in response to gene flow differ (i.e., alleles versus loci). I argue that, although making this distinction is important, we commonly lack the critical empirical evidence to do so. This is mainly because many current approaches are biased towards detecting local adaptation in genome regions with low crossover rates. The importance of low-crossover genome regions for adaptation with gene flow, such as in co-localizing relevant loci within a genome, thus remains unclear. Future empirical investigations should address these questions by making use of comparative genomics, where multiple de novo genome assemblies from species evolved under different degrees of genetic exchange are compared. This research promises to advance our understanding of how a genome adapts to maladaptive gene flow, thereby promoting adaptive divergence and reproductive isolation.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

Selfing Can Facilitate Transitions between Pollination Syndromes

Pollinator-mediated selection on plants can favor transitions to a new pollinator depending on the relative abundances and efficiencies of pollinators present in the community. A frequently observed example is the transition from bee pollination to hummingbird pollination. We present a population genetic model that examines whether the ability to inbreed can influence evolutionary change in traits that underlie pollinator attraction. We find that a transition to a more efficient but less abundant pollinator is favored under a broadened set of ecological conditions if plants are capable of delayed selfing rather than obligately outcrossing. Delayed selfing allows plants carrying an allele that attracts the novel pollinator to reproduce even when this pollinator is rare, providing reproductive assurance. In addition, delayed selfing weakens the effects of Haldane's sieve by increasing the fixation probability for recessive alleles that confer adaptation to the new pollinator. Our model provides novel insight into the paradoxical abundance of recessive mutations in adaptation to hummingbird attraction. It further predicts that transitions to efficient but less abundant pollinators (such as hummingbirds in certain communities) should disproportionately occur in self-compatible lineages. Currently available mating system data sets are consistent with this prediction, and we suggest future areas of research that will enable a rigorous test of this theory.

Genome-wide identification of physically clustered genes suggests chromatin-level co-regulation in male reproductive development in Arabidopsis thaliana

Co-expression of physically linked genes occurs surprisingly frequently in eukaryotes. Such chromosomal clustering may confer a selective advantage as it enables coordinated gene regulation at the chromatin level. We studied the chromosomal organization of genes involved in male reproductive development in Arabidopsis thaliana. We developed an in-silico tool to identify physical clusters of co-regulated genes from gene expression data. We identified 17 clusters (96 genes) involved in stamen development and acting downstream of the transcriptional activator MS1 (MALE STERILITY 1), which contains a PHD domain associated with chromatin re-organization. The clusters exhibited little gene homology or promoter element similarity, and largely overlapped with reported repressive histone marks. Experiments on a subset of the clusters suggested a link between expression activation and chromatin conformation: qRT-PCR and mRNA in situ hybridization showed that the clustered genes were up-regulated within 48 h after MS1 induction; out of 14 chromatin-remodeling mutants studied, expression of clustered genes was consistently down-regulated only in hta9/hta11, previously associated with metabolic cluster activation; DNA fluorescence in situ hybridization confirmed that transcriptional activation of the clustered genes was correlated with open chromatin conformation. Stamen development thus appears to involve transcriptional activation of physically clustered genes through chromatin de-condensation.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Trait Correlations in the Genomics Era

Thinking about the evolutionary causes and consequences of trait correlations has been dominated by quantitative genetics theory that is focused on hypothetical loci. Since this theory was initially developed, technology has enabled the identification of specific genetic variants that contribute to trait correlations. Here, we review studies of the genetic basis of trait correlations to ask: What has this new information taught us? We find that causal variants can be pleiotropic and/or linked in different ways, indicating that pleiotropy and linkage are not alternative genetic mechanisms. Further, many trait correlations have a polygenic basis, suggesting that both pleiotropy and linkage likely contribute. We discuss implications of these findings for the evolutionary causes and consequences of trait correlations.

Genetic Analysis of Natural Variation in Antirrhinum Scent Profiles Identifies BENZOIC ACID CARBOXYMETHYL TRANSFERASE As the Major Locus Controlling Methyl Benzoate Synthesis

The Antirrhinum genus has a considerable complexity in the scent profiles produced by different species. We have analyzed the genetic differences between A. majus and A. linkianum, two species divergent in the emission of methyl benzoate, methyl cinnamate, acetophenone, and ocimene. The genetic analysis showed that all compounds segregated in a Mendelian fashion attributable to one or two loci with simple or epistatic interactions. Several lines lacked methyl benzoate, a major Volatile Organic Compound emitted by A. majus but missing in A. linkianum. Using a candidate gene approach, we found that the BENZOIC ACID CARBOXYMETHYL TRANSFERASE from A. linkianum appeared to be a null allele as we could not detect mRNA expression. The coding region did not show significant differences that could explain the loss of expression. The intron-exon boundaries was also conserved indicating that there is no alternative splicing in A. linkianum as compared to A. majus. However, it showed multiple polymorphisms in the 5' promoter region including two insertions, one harboring an IDLE MITE transposon with additional sequences with high homology to the PLENA locus and a second one with somewhat lower homology to the regulatory region of the VENOSA locus. It also had a 778 bp deletion as compared to the A. majus BAMT promoter region. Our results show that the differences in scent emission between A. majus and A. linkianum may be traced back to single genes involved in discrete biosynthetic reactions such as benzoic acid methylation. Thus, natural variation of this complex trait maybe the result of combinations of wild type, and loss of function alleles in different genes involved in discrete VOCs biosynthesis. Furthermore, the presence of active transposable elements in the genus may account for rapid evolution and instability, raising the possibility of adaptation to local pollinators.

Gain and Loss of Floral Scent Production through Changes in Structural Genes during Pollinator-Mediated Speciation

The interactions of plants with their pollinators are thought to be a driving force in the evolution of angiosperms. Adaptation to a new pollinator involves coordinated changes in multiple floral traits controlled by multiple genes. Surprisingly, such complex genetic shifts have happened numerous times during evolution. Here we report on the genetic basis of the changes in one such trait, floral scent emission, in the genus Petunia (Solanaceae). The increase in the quantity and complexity of the volatiles during the shift from bee to hawkmoth pollination was due to de novo expression of the genes encoding benzoic acid/salicylic acid carboxyl methyltransferase (BSMT) and benzoyl-CoA:benzylalcohol/2-phenylethanol benzoyltransferase (BPBT) together with moderately increased transcript levels for most enzymes of the phenylpropanoid/benzenoid pathway. Loss of cinnamate-CoA ligase (CNL) function as well as a reduction in the expression of the MYB transcription factor ODO1 explain the loss of scent during the transition from moth to hummingbird pollination. The CNL gene in the hummingbird-adapted species is inactive due to a stop codon, but also appears to have undergone further degradation over time. Therefore, we propose that loss of scent happened relatively early in the transition toward hummingbird pollination, and probably preceded the loss of UV-absorbing flavonols. The discovery that CNL is also involved in the loss of scent during the transition from outcrossing to selfing in Capsella (Brassicaceae) (see the accompanying paper) raises interesting questions about the possible causes of deep evolutionary conservation of the targets of evolutionary change.

The molecular bases of floral scent evolution under artificial selection: insights from a transcriptome analysis in Brassica rapa

In an artificial selection experiment using fast-cycling Brassica rapa plants it was recently shown that floral VOCs respond rapidly to selection for increased amounts. Here we carried out transcriptome analysis in these plants to explore the molecular bases of the augmentation in the artificially selected scent compound, phenylacetaldehyde (PAA), as well as other compounds that increased through pleiotropy. In the transcriptome data, we found up-regulation of genes likely underlying PAA synthesis, but also several genes of the shikimate pathway and the related phenylalanine metabolism. As phenylalanine is the precursor of many aromatic volatiles that showed increased emission, this result could explain some of the pleiotropic evolutionary responses. In addition, we found that ribosomal protein genes were up-regulated in “high” (high PAA amount) selection line plants, a mechanism that might further augment the effect of elevated gene expression at the proteomic level. Our study shows that selection on an individual trait can impose changes in the expression of several different genes, which could explain pleiotropic responses in the biosynthetic network of floral volatiles.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

Accessibility, constraint, and repetition in adaptive floral evolution

Adaptive phenotypic evolution is shaped by natural selection on multiple organismal traits as well as by genetic correlations among traits. Genetic correlations can arise through pleiotropy and can bias the production of phenotypic variation to certain combinations of traits. This phenomenon is referred to as developmental bias or constraint. Developmental bias may accelerate or constrain phenotypic evolution, depending on whether selection acts parallel or in opposition to genetic correlations among traits. We discuss examples from floral evolution where genetic correlations among floral traits contribute to rapid, coordinated evolution in multiple floral organ phenotypes and suggest future research directions that will explore the relationship between the genetic basis of adaptation and the pre-existing structure of genetic correlations. On the other hand, natural selection may act perpendicular to a strong genetic correlation, for example when two traits are encoded by a subset of the same genes and natural selection favors change in one trait and stability in the second trait. In such cases, adaptation is constrained by the availability of genetic variation that can influence the focal trait with minimal pleiotropic effects. Examples from plant diversification suggest that the origin of certain adaptations depends on the prior evolution of a gene copy with reduced pleiotropic effects, generated through the process of gene duplication followed by subfunctionalization or neofunctionalization. A history of gene duplication in some developmental pathways appears to have allowed particular flowering plant linages to have repeatedly evolved adaptations that might otherwise have been developmentally constrained.

A Segregating Inversion Generates Fitness Variation in Yellow Monkeyflower (Mimulus guttatus)

Polymorphic chromosomal rearrangements can bind hundreds of genes into single genetic loci with diverse effects. Rearrangements are often associated with local adaptation and speciation and may also be an important component of genetic variation within populations. We genetically and phenotypically characterize a segregating inversion (inv6) in the Iron Mountain (IM) population of Mimulus guttatus (yellow monkeyflower). We initially mapped inv6 as a region of recombination suppression in three F2 populations resulting from crosses among IM plants. In each case, the F1 parent was heterozygous for a derived haplotype, homogenous across markers spanning over 5 Mb of chromsome 6. In the three F2 populations, inv6 reduced male and female fitness components. In addition,i nv6 carriers suffered an ∼30% loss of pollen viability in the field. Despite these costs, inv6 exists at moderate frequency (∼8%) in the natural population, suggesting counterbalancing fitness benefits that maintain the polymorphism. Across 4 years of monitoring in the field, inv6 had an overall significant positive effect on seed production (lifetime female fitness) of carriers. This benefit was particularly strong in harsh years and may be mediated (in part) by strong positive effects on flower production. These data suggest that opposing fitness effects maintain an intermediate frequency, and as a consequence, inv6 generates inbreeding depression and high genetic variance. We discuss these findings in relation to the theory of inbreeding depression and the maintenance of fitness variation.