The molecular basis of phenotypic evolution: beyond the usual suspects.

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

Genetic basis of nectar guide trichome variation between bumblebee- and self-pollinated monkeyflowers (Mimulus): role of the MIXTA-like gene GUIDELESS.

Nectar guide trichomes play crucial ecological roles in bee-pollinated flowers, as they serve as footholds and guides for foraging bees to access the floral rewards. However, the genetic basis of natural variation in nectar guide trichomes among species remains poorly understood. In this study, we performed genetic analysis of nectar guide trichome variation between two closely related monkeyflower (Mimulus) species, the bumblebee-pollinated Mimulus lewisii and self-pollinated M. parishii. We demonstrate that a MIXTA-like R2R3-MYB gene, GUIDELESS, is a major contributor to the nectar guide trichome length variation between the two species. The short-haired M. parishii carries a recessive allele due to non-synonymous substitutions in a highly conserved motif among MIXTA-like MYB proteins. Furthermore, our results suggest that besides GUIDELESS, additional loci encoding repressors of trichome elongation also contribute to the transition from bumblebee-pollination to selfing. Taken together, these results suggest that during a pollination syndrome switch, changes in seemingly complex traits such as nectar guide trichomes could have a relatively simple genetic basis, involving just a few genes of large effects.

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe).

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson's D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.

To stripe or not to stripe: the origin of a novel foliar pigmentation pattern in monkeyflowers (Mimulus).

The origin of phenotypic novelty is one of the most challenging problems in evolutionary biology. Although genetic regulatory network rewiring or co-option has been widely recognised as a major contributor, in most cases how such genetic rewiring/co-option happens is completely unknown. We have studied a novel foliar pigmentation pattern that evolved recently in the monkeyflower species Mimulus verbenaceus. Through genome-wide association tests using wild-collected samples, experimental crosses of laboratory inbred lines, gene expression analyses, and functional assays, we identified an anthocyanin-activating R2R3-MYB gene, STRIPY, as the causal gene triggering the emergence of the discrete, mediolateral anthocyanin stripe in the M. verbenaceus leaf. Chemical mutagenesis revealed the existence of upstream activators and repressors that form a 'hidden' prepattern along the leaf proximodistal axis, potentiating the unique expression pattern of STRIPY. Population genomics analyses did not reveal signatures of positive selection, indicating that nonadaptive processes may be responsible for the establishment of this novel trait in the wild. This study demonstrates that the origin of phenotypic novelty requires at least two separate phases, potentiation and actualisation. The foliar stripe pattern of M. verbenaceus provides an excellent platform to probe the molecular details of both processes in future studies.

A Tetratricopeptide Repeat Protein Regulates Carotenoid Biosynthesis and Chromoplast Development in Monkeyflowers (Mimulus).

Little is known about the factors regulating carotenoid biosynthesis in flowers. Here, we characterized the REDUCED CAROTENOID PIGMENTATION2 (RCP2) locus from two monkeyflower (Mimulus) species, the bumblebee-pollinated species Mimulus lewisii and the hummingbird-pollinated species Mimulus verbenaceus We show that loss-of-function mutations of RCP2 cause drastic down-regulation of the entire carotenoid biosynthetic pathway. The causal gene underlying RCP2 encodes a tetratricopeptide repeat protein that is closely related to the Arabidopsis (Arabidopsis thaliana) REDUCED CHLOROPLAST COVERAGE proteins. RCP2 appears to regulate carotenoid biosynthesis independently of RCP1, a previously identified R2R3-MYB master regulator of carotenoid biosynthesis. We show that RCP2 is necessary and sufficient for chromoplast development and carotenoid accumulation in floral tissues. Simultaneous down-regulation of RCP2 and two closely related paralogs, RCP2-L1 and RCP2-L2, yielded plants with pale leaves deficient in chlorophyll and carotenoids and with reduced chloroplast compartment size. Finally, we demonstrate that M. verbenaceus is just as amenable to chemical mutagenesis and in planta transformation as the more extensively studied M. lewisii, making these two species an excellent platform for comparative developmental genetics studies of closely related species with dramatic phenotypic divergence.

Developmental Genetics of Corolla Tube Formation: Role of the tasiRNA-ARF Pathway and a Conceptual Model.

Over 80,000 angiosperm species produce flowers with petals fused into a corolla tube. The corolla tube contributes to the tremendous diversity of flower morphology and plays a critical role in plant reproduction, yet it remains one of the least understood plant structures from a developmental genetics perspective. Through mutant analyses and transgenic experiments, we show that the tasiRNA-ARF pathway is required for corolla tube formation in the monkeyflower species Mimulus lewisii Loss-of-function mutations in the M. lewisii orthologs of ARGONAUTE7 and SUPPRESSOR OF GENE SILENCING3 cause a dramatic decrease in abundance of TAS3-derived small RNAs and a moderate upregulation of AUXIN RESPONSE FACTOR3 (ARF3) and ARF4, which lead to inhibition of lateral expansion of the bases of petal primordia and complete arrest of the upward growth of the interprimordial regions, resulting in unfused corollas. Using the DR5 auxin-responsive promoter, we discovered that auxin signaling is continuous along the petal primordium base and the interprimordial region during the critical stage of corolla tube formation in the wild type, similar to the spatial pattern of MlARF4 expression. Auxin response is much weaker and more restricted in the mutant. Furthermore, exogenous application of a polar auxin transport inhibitor to wild-type floral apices disrupted petal fusion. Together, these results suggest a new conceptual model highlighting the central role of auxin-directed synchronized growth of the petal primordium base and the interprimordial region in corolla tube formation.

Comparative analysis of corolla tube development across three closely related Mimulus species with different pollination syndromes.

Fusion of petals to form a corolla tube is considered a key innovation contributing to the diversification of many flowering plant lineages. Corolla tube length often varies dramatically among species and is a major determinant of pollinator preference. However, our understanding of the developmental dynamics underlying corolla tube length variation is very limited. Here we examined corolla tube growth in the Mimulus lewisii species complex, an emerging model system for studying the developmental genetics and evo-devo of pollinator-associated floral traits. We compared developmental and cellular processes associated with corolla tube length variation among the bee-pollinated M. lewisii, the hummingbird-pollinated Mimulus verbenaceus, and the self-pollinated Mimulus parishii. We found that in all three species, cell size is non-uniformly distributed along the mature tube, with the longest cells just distal to the stamen insertion site. Differences in corolla tube length among the three species are not associated with processes of organogenesis or early development but are associated with variation in multiple processes occurring later in development, including the location and duration of cell division and cell elongation. The tube growth curves of the small-flowered M. parishii and large-flowered M. lewisii are essentially indistinguishable, except that M. parishii tubes stop growing earlier at a smaller size, suggesting a critical role of heterochrony in the shift from outcrossing to selfing. These results not only highlight the developmental process associated with corolla tube variation among species but also provide a baseline reference for future developmental genetic analyses of mutants or transgenic plants with altered corolla tube morphology in this emerging model system.

Repressors of anthocyanin biosynthesis.

Anthocyanins play a variety of adaptive roles in both vegetative tissues and reproductive organs of plants. The broad functionality of these compounds requires sophisticated regulation of the anthocyanin biosynthesis pathway to allow proper localization, timing, and optimal intensity of pigment deposition. While it is well-established that the committed steps of anthocyanin biosynthesis are activated by a highly conserved MYB-bHLH-WDR (MBW) protein complex in virtually all flowering plants, anthocyanin repression seems to be achieved by a wide variety of protein and small RNA families that function in different tissue types and in response to different developmental, environmental, and hormonal cues. In this review, we survey recent progress in the identification of anthocyanin repressors and the characterization of their molecular mechanisms. We find that these seemingly very different repression modules act through a remarkably similar logic, the so-called 'double-negative logic'. Much of the double-negative regulation of anthocyanin production involves signal-induced degradation or sequestration of the repressors from the MBW protein complex. We discuss the functional and evolutionary advantages of this logic design compared with simple or sequential positive regulation. These advantages provide a plausible explanation as to why plants have evolved so many anthocyanin repressors.

The leaf polarity factors SGS3 and YABBYs regulate style elongation through auxin signaling in Mimulus lewisii.

Style length is a major determinant of breeding strategies in flowering plants and can vary dramatically between and within species. However, little is known about the genetic and developmental control of style elongation. We characterized the role of two classes of leaf adaxial-abaxial polarity factors, SUPPRESSOR OF GENE SILENCING3 (SGS3) and the YABBY family transcription factors, in the regulation of style elongation in Mimulus lewisii. We also examined the spatiotemporal patterns of auxin response during style development. Loss of SGS3 function led to reduced style length via limiting cell division, and downregulation of YABBY genes by RNA interference resulted in shorter styles by decreasing both cell division and cell elongation. We discovered an auxin response minimum between the stigma and ovary during the early stages of pistil development that marks style differentiation. Subsequent redistribution of auxin response to this region was correlated with style elongation. Auxin response was substantially altered when both SGS3 and YABBY functions were disrupted. We suggest that auxin signaling plays a central role in style elongation and that the way in which auxin signaling controls the different cell division and elongation patterns underpinning natural style length variation is a major question for future research.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower.

Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.

A Novel R2R3-MYB Transcription Factor Contributes to Petal Blotch Formation by Regulating Organ-Specific Expression of PsCHS in Tree Peony (Paeonia suffruticosa).

Flower color patterns play critical roles in plant-pollinator interactions and represent one of the most common adaptations during angiosperm evolution. However, the molecular mechanisms underlying flower color pattern formation are less understood in non-model organisms. The aim of this study was to identify genes involved in the formation of petal blotches in tree peony (Paeonia suffruticosa) through transcriptome profiling and functional experiments. We identified an R2R3-MYB gene, PsMYB12, representing a distinct R2R3-MYB subgroup, with a spatiotemporal expression pattern tightly associated with petal blotch development. We further demonstrated that PsMYB12 interacts with a basic helix-loop-helix (bHLH) and a WD40 protein in a regulatory complex that directly activates PsCHS expression, which is also specific to the petal blotches. Together, these findings advance our understanding of the molecular mechanisms of pigment pattern formation beyond model plants. They also benefit molecular breeding of tree peony cultivars with novel color patterns and promote germplasm innovation.

Monkeyflowers (Mimulus): new model for plant developmental genetics and evo-devo.

Contents Summary 694 I. Introduction 694 II. The system 695 III. Regulation of carotenoid pigmentation 695 IV. Formation of periodic pigmentation patterns 696 V. Developmental genetics of corolla tube formation and elaboration 697 VI. Molecular basis of floral trait variation underlying pollinator shift 698 VII. Outlook 699 Acknowledgements 699 References 699 SUMMARY: Monkeyflowers (Mimulus) have long been recognized as a classic ecological and evolutionary model system. However, only recently has it been realized that this system also holds great promise for studying the developmental genetics and evo-devo of important plant traits that are not found in well-established model systems such as Arabidopsis. Here, I review recent progress in four different areas of plant research enabled by this new model, including transcriptional regulation of carotenoid biosynthesis, formation of periodic pigmentation patterns, developmental genetics of corolla tube formation and elaboration, and the molecular basis of floral trait divergence underlying pollinator shift. These examples suggest that Mimulus offers ample opportunities to make exciting discoveries in plant development and evolution.

Competition between anthocyanin and flavonol biosynthesis produces spatial pattern variation of floral pigments between Mimulus species.

Flower color patterns have long served as a model for developmental genetics because pigment phenotypes are visually striking, yet generally not required for plant viability, facilitating the genetic analysis of color and pattern mutants. The evolution of novel flower colors and patterns has played a key role in the adaptive radiation of flowering plants via their specialized interactions with different pollinator guilds (e.g., bees, butterflies, birds), motivating the search for allelic differences affecting flower color pattern in closely related plant species with different pollinators. We have identified LIGHT AREAS1 (LAR1), encoding an R2R3-MYB transcription factor, as the causal gene underlying the spatial pattern variation of floral anthocyanin pigmentation between two sister species of monkeyflower: the bumblebee-pollinated Mimulus lewisii and the hummingbird-pollinated Mimulus cardinalis. We demonstrated that LAR1 positively regulates FLAVONOL SYNTHASE (FLS), essentially eliminating anthocyanin biosynthesis in the white region (i.e., light areas) around the corolla throat of M. lewisii flowers by diverting dihydroflavonol into flavonol biosynthesis from the anthocyanin pigment pathway. FLS is preferentially expressed in the light areas of the M. lewisii flower, thus prepatterning the corolla. LAR1 expression in M. cardinalis flowers is much lower than in M. lewisii, explaining the unpatterned phenotype and recessive inheritance of the M. cardinalis allele. Furthermore, our gene-expression analysis and genetic mapping results suggest that cis-regulatory change at the LAR1 gene played a critical role in the evolution of different pigmentation patterns between the two species.

Saltatory rolling circle amplification for sensitive visual detection of Staphylococcus aureus in milk.

Monitoring Staphylococcus aureus with high sensitivity is very important for ensuring milk quality and food safety. In this study, we used a rapid nucleic acid isothermal amplification method, saltatory rolling circle amplification (SRCA), for the detection of Staph. aureus in milk. The results of the SRCA method can be assessed visually by the presence of white precipitate or by fluorescence measurement. Thirteen Staph. aureus strains and 31 non-Staph. aureus strains were used to evaluate the specificity of SRCA. The method exhibited excellent detection of Staph. aureus genomic DNA at a concentration of 7.8 × 101 fg/µL when assessed by visible precipitate, and at 7.8 × 100 fg/µL when detected by fluorescence after addition of the fluorochrome SYBR Green I. In artificially inoculated milk, the detection limits of SRCA were 5.6 × 102 cfu/mL by precipitate and 5.6 × 101 cfu/mL by fluorescence, respectively. Compared with conventional PCR approaches, the SRCA assay achieved at least 100-fold higher sensitivity. Moreover, the sensitivity, specificity, and accuracy of the SRCA-based system were calculated to be 100.00, 97.73, and 97.78%, respectively. These results indicate that SRCA has potential application as a sensitive and visual technique for the detection of Staph. aureus in milk.

A dominant-negative actin mutation alters corolla tube width and pollinator visitation in Mimulus lewisii.

A third of all angiosperm species produce flowers with petals fused into a corolla tube. The various elaborations of corolla tube attributes, such as length, width and curvature, have enabled plants to exploit many specialized pollinator groups. These elaborations often differ dramatically among closely related species, contributing to pollinator shift and pollinator-mediated reproductive isolation and speciation. However, very little is known about the genetic and developmental control of these corolla tube attributes. Here we report the characterization of a semi-dominant mutant in the monkeyflower species Mimulus lewisii, with a substantial decrease in corolla tube width but no change in tube length. This morphological alteration leads to a ˜ 70% decrease in bumblebee visitation rate for the homozygous mutant compared to the wild-type. Through bulk segregant analysis and transgenic experiment, we show that the mutant phenotype is caused by a dominant-negative mutation in an actin gene. This mutation decreases epidermal cell width but not length, and probably also reduces the number of lateral cell divisions. These results suggest a surprising potential role for a 'housekeeping' gene in fine-tuning the development of an ecologically important floral trait.

An R2R3-MYB transcription factor regulates carotenoid pigmentation in Mimulus lewisii flowers.

Carotenoids are yellow, orange, and red pigments that contribute to the beautiful colors and nutritive value of many flowers and fruits. The structural genes in the highly conserved carotenoid biosynthetic pathway have been well characterized in multiple plant systems, but little is known about the transcription factors that control the expression of these structural genes. By analyzing a chemically induced mutant of Mimulus lewisii through bulk segregant analysis and transgenic experiments, we have identified an R2R3-MYB, Reduced Carotenoid Pigmentation 1 (RCP1), as the first transcription factor that positively regulates carotenoid biosynthesis during flower development. Loss-of-function mutations in RCP1 lead to down-regulation of all carotenoid biosynthetic genes and reduced carotenoid content in M. lewisii flowers, a phenotype recapitulated by RNA interference in the wild-type background. Overexpression of this gene in the rcp1 mutant background restores carotenoid production and, unexpectedly, results in simultaneous decrease of anthocyanin production in some transgenic lines by down-regulating the expression of an activator of anthocyanin biosynthesis. Identification of transcriptional regulators of carotenoid biosynthesis provides the 'toolbox' genes for understanding the molecular basis of flower color diversification in nature and for potential enhancement of carotenoid production in crop plants via genetic engineering.

Testing the utility of fluorescent proteins in Mimulus lewisii by an Agrobacterium-mediated transient assay.

KEY MESSAGE: The Agrobacterium -mediated transient expression assay by leaf infiltration in Mimulus lewisii is robust. Fluorescent proteins EGFP, EYFP and DsRed give bright fluorescence signals in the infiltrated tissue. Mimulus lewisii is an emerging developmental genetic model system. Recently developed genomic and genetic resources and a stable transformation protocol have greatly facilitated the identification and functional characterization of genes controlling the development of ecologically important floral traits using this species. To further expedite gene and protein function analyses in M. lewisii, we adopted and simplified the Agrobacterium-mediated transient gene expression method routinely used in tobacco plants. With the validated transient assay, we examined the performance of fluorescent proteins EGFP, EYFP and DsRed in M. lewisii. All three proteins gave bright fluorescence signals when transiently expressed in agroinfiltrated leaves. Furthermore, we demonstrated the utility of fluorescent proteins in M. lewisii by showing the nuclear localization of Reduced Carotenoid Pigmentation 1 (RCP1), a recently discovered R2R3-MYB transcription factor that regulates carotenoid pigmentation during flower development. Both the transient assay and the fluorescent proteins are valuable additions to the M. lewisii toolbox, making this emerging genetic and developmental model system even more powerful.

Fine-scale variation in meiotic recombination in Mimulus inferred from population shotgun sequencing.

Meiotic recombination rates can vary widely across genomes, with hotspots of intense activity interspersed among cold regions. In yeast, hotspots tend to occur in promoter regions of genes, whereas in humans and mice, hotspots are largely defined by binding sites of the positive-regulatory domain zinc finger protein 9. To investigate the detailed recombination pattern in a flowering plant, we use shotgun resequencing of a wild population of the monkeyflower Mimulus guttatus to precisely locate over 400,000 boundaries of historic crossovers or gene conversion tracts. Their distribution defines some 13,000 hotspots of varying strengths, interspersed with cold regions of undetectably low recombination. Average recombination rates peak near starts of genes and fall off sharply, exhibiting polarity. Within genes, recombination tracts are more likely to terminate in exons than in introns. The general pattern is similar to that observed in yeast, as well as in positive-regulatory domain zinc finger protein 9-knockout mice, suggesting that recombination initiation described here in Mimulus may reflect ancient and conserved eukaryotic mechanisms.

Carotenoid composition of the flowers of Mimulus lewisii and related species: Implications regarding the prevalence and origin of two unique, allenic pigments.

The genus Mimulus has been used as a model system in a wide range of ecological and evolutionary studies and contains many species with carotenoid pigmented flowers. However, the detailed carotenoid composition of these flowers has never been reported. In this paper the floral carotenoid composition of 11 Mimulus species are characterized using high-performance liquid chromatography, mass spectrometry and chemical methods with a particular focus on the genetic model species, Mimulus lewisii. M. lewisii flowers have five major carotenoids: antheraxanthin, violaxanthin, neoxanthin, and the unique allenic carotenoids, deepoxyneoxanthin and mimulaxanthin. This carotenoid profile is consistent with the expression levels of putative carotenoid biosynthetic genes in the M. lewisii flower. The other 10 species possess the same five carotenoids or a subset of these. Comparison of the carotenoid profiles among species in a phylogenetic context provides new insights into the biosynthesis and evolution of deepoxyneoxanthin and mimulaxanthin. This work also lays the foundation for future studies regarding transcriptional control of the carotenoid biosynthesis pathway in Mimulus flowers.

Transcriptional control of floral anthocyanin pigmentation in monkeyflowers (Mimulus).

A molecular description of the control of floral pigmentation in a multi-species group displaying various flower color patterns is of great interest for understanding the molecular bases of phenotypic diversification and pollinator-mediated speciation. Through transcriptome profiling, mutant analyses and transgenic experiments, we aim to establish a 'baseline' floral anthocyanin regulation model in Mimulus lewisii and to examine the different ways of tinkering with this model in generating the diversity of floral anthocyanin patterns in other Mimulus species. We find one WD40 and one bHLH gene controlling anthocyanin pigmentation in the entire corolla of M. lewisii and two R2R3-MYB genes, PELAN and NEGAN, controlling anthocyanin production in the petal lobe and nectar guide, respectively. The autoregulation of NEGAN might be a critical property to generate anthocyanin spots. Independent losses of PELAN expression (via different mechanisms) explain two natural yellow-flowered populations of M. cardinalis (typically red-flowered). The NEGAN ortholog is the only anthocyanin-activating MYB expressed in the M. guttatus flowers. The mutant lines and transgenic tools available for M. lewisii will enable gene-by-gene replacement experiments to dissect the genetic and developmental bases of more complex floral color patterns, and to test hypotheses on phenotypic evolution in general.