- Research article
- Open Access
Stability of petal color polymorphism: the significance of anthocyanin accumulation in photosynthetic tissues
BMC Plant Biology volume 19, Article number: 496 (2019)
Anthocyanins are the primary source of colour in flowers and also accumulate in vegetative tissues, where they have multiple protective roles traditionally attributed to early compounds of the metabolic pathway (flavonols, flavones, etc.). Petal-specific loss of anthocyanins in petals allows plants to escape from the negative pleiotropic effects of flavonoid and anthocyanins loss in vegetative organs, where they perform a plethora of essential functions. Herein, we investigate the degree of pleiotropy at the biochemical scale in a pink-white flower colour polymorphism in the shore campion, Silene littorea. We report the frequencies of pink and white individuals across 21 populations and underlying biochemical profiles of three flower colour variants: anthocyanins present in all tissues (pink petals), petal-specific loss of anthocyanins (white petals), and loss of anthocyanins in all tissues (white petals).
Individuals lacking anthocyanins only in petals represent a stable polymorphism in two populations at the northern edge of the species range (mean frequency 8–21%). Whereas, individuals lacking anthocyanins in the whole plant were found across the species range, yet always at very low frequencies (< 1%). Biochemically, the flavonoids detected were anthocyanins and flavones; in pigmented individuals, concentrations of flavones were 14–56× higher than anthocyanins across tissues with differences of > 100× detected in leaves. Loss of anthocyanin pigmentation, either in petals or in the whole plant, does not influence the ability of these phenotypes to synthesize flavones, and this pattern was congruent among all sampled populations.
We found that all colour variants showed similar flavone profiles, either in petals or in the whole plant, and only the flower colour variant with anthocyanins in photosynthetic tissues persists as a stable flower colour polymorphism. These findings suggest that anthocyanins in photosynthetic tissues, not flavonoid intermediates, are the targets of non-pollinator mediated selection.
Mutations are the primary source of genetic variation in all organisms and have a key contribution to phenotypic diversity [1, 2], but not all mutations are evolutionarily relevant. Some phenotypic changes are produced through spontaneous mutations with deleterious effects that are consistently eliminated by purifying selection . In contrast, persistent phenotypic changes arise from mutations maintained by balancing selection through frequency-dependent or heterogeneous selection or through the promotion of multiple adaptive peaks [4,5,6,7,8], resulting in a population polymorphism for that trait [9,10,11]. Factors that may determine why some new phenotypes are fleeting and some persists as polymorphisms are still open [12, 13], but the study of flower color is helping to shed light to this issue [14,15,16,17].
Flower color variation has drawn the attention of many naturalists through the history [18,19,20], and nowadays continues to be an important focus of research to evolutionary biologist [14, 21]. Flower color has been considered as an adaptive trait for pollinator attraction , but underlying pigments also have other functions, especially in vegetative tissues. Anthocyanins are the most common plant pigment that color flowers, conferring orange, red, pink and blue colors  that attract diverse functional groups of pollinators [22, 24]. For example, variation in anthocyanin content in monkeyflowers (Mimulus) results in red and pink-flowered species that are visited by hummingbirds and bees, respectively [25, 26]. In vegetative tissues, anthocyanins may show protective roles such as sunscreens, antioxidants or antipathogens, among others [27, 28]. Thus, loss of anthocyanins may affect pollinator activity, but may also have physiological effects depending on whether they are accumulated or not in vegetative tissues. If the loss of anthocyanins is confined to the flowers (usually in petals), the rest of the plant can produce anthocyanins and reduce any negative pleiotropic effects in other tissues [29, 30], whereas anthocyanins-lacking individuals in the whole plant potentially grow and reproduce, but frequently exhibit fitness disadvantages that seem to explain their scarcity in the wild [31,32,33] (see Additional file 1: Table S1).
Petal-specific loss of anthocyanins is frequently induced by regulatory mutations (that is changes in the regulation of gene expression) and shows a mutation bias to Myb transcription factors, the key regulatory factors controlling anthocyanin biosynthesis in plants [30, 34]. Different copies of Myb proteins regulate floral and vegetative anthocyanins, thus the specificity of this regulation is predicted to have low pleiotropic consequences . For example, in Ipomoea purpurea mutations the regulatory IpMyb1 gene are responsible for anthocyanin loss in pigmented flowers . However, these mutations do not affect the fitness of white-flowered plants and show equal or even higher reproductive success than that of the pigmented individuals [15, 29].
Anthocyanin-lacking individuals, on the other hand, are conferred by loss-of-function in any of the structural loci or whole plant regulatory genes of the anthocyanin biosynthetic pathway (hereafter ABP) [30, 34]. Loss-of-function mutations may target a high spectrum of genes since there are more possible loci that could confer the non-pigmented phenotype. Thus, the inactivation of any structural gene of the pathway often limit the flux down the ABP and block the anthocyanin production, but also may affect the synthesis of uncolored/pale-yellow non-anthocyanin flavonoids in the side branches of the pathway . These flavonoids, such as flavones or flavonols, also perform important ecological functions because they show similar or even more protective functions against environmental stressors than anthocyanins themselves . Therefore, the persistence of loss-of-function mutations should be limited by the negative pleiotropic effects associated to the absence of anthocyanins and/or intermediate non-anthocyanin flavonoids [14, 17, 38]. The selection against these variants may depend on where the mutation occurs and the associated negative consequences for the flavonoid loss.
Loss of pigmentation, particularly due to absence of anthocyanins, represents the most frequent cases of flower color polymorphism in plants [39, 40]. White-flowered morphs represent valuable natural genotypes to know the possible selective disadvantages of lack of anthocyanins in the whole plant [33, 39, 41,42,43], but the non-anthocyanin flavonoid composition of such plants is unknown. In addition, the quantities of anthocyanins are usually correlated with those of non-anthocyanin flavonoids, at least in some tissues, and the concentrations of the latter are even higher than that of anthocyanins [44, 45]. Consequently, it is difficult to distinguish which group of flavonoids is responsible for the putative selective disadvantage of anthocyanin lacking plants [31,32,33], and studies that clearly differentiate between flavonoid groups are limited.
The shore campion Silene littorea Brot. (Caryophyllaceae) is an annual pink-flowered species that accumulates anthocyanins and non-anthocyanin flavonoids in petals and in calyces, leaves and stems [44, 46]. The accumulation of both kinds of flavonoids in vegetative tissues is highly variable and seems to respond to light stress . Silene littorea grows along the Iberian coast, and exhibits an anthocyanin-based pink-white flower color polymorphism in two populations of the northwest distribution range, but anthocyanin-lacking individuals are occasionally observed in some populations . In S. littorea, petal-specific polymorphism is likely due to downregulation of the flavanone-3-hydroxylase (F3 h) gene through a downregulation of the SlMyb1a transcription factor , but genetic causes of anthocyanin-lacking plants are still unknown. Flowers of this species are mainly visited by generalist pollinators from the orders Diptera, Hymenoptera, and Lepidoptera; however, they do not seem to show strong pollinator preference for either pink or white flowers (M.L.B. 2019, unpublished data). Thus, the occurrence of petal anthocyanin loss and whole-plant anthocyanin loss individuals in S. littorea (hereafter PAL and WAL phenotypes, respectively) offers an excellent opportunity for understanding the importance of non-pollinator selection due to lack of anthocyanins and/or non-anthocyanin flavonoids.
In this study, we seek to understand the factors that determine the fate of different forms of anthocyanin variation in S. littorea. Thus, we investigated the population frequency of three anthocyanin phenotypes (PAL, WAL and fully pigmented) in 21 populations across the species distribution range over five years. Then, we used high-performance liquid chromatography coupled with diode-array detection and electrospray ionization tandem mass spectrometry (HPLC-DAD-MSn) to study the flavonoid profiles at the whole plant level in the fully pigmented phenotype and compare them to that of white-flowered variants (i.e. PAL and WAL phenotypes); after that, this study was expanded to more individuals and populations using spectrophotometric quantification of flavonoids. Because of the negative consequences of the absence of flavonoids [30, 38], we expect the PAL phenotype to be more common within a population compared to WAL. Thus, loss of anthocyanins and non-anthocyanin flavonoids are expected to be limited to petals in PAL plants, but extended to the whole plant in WAL individuals.
Frequency of PAL, WAL and fully pigmented phenotypes
Our population surveys confirmed that the PAL phenotype is limited to two populations in the northern portion of the species range (Fig. 1). A polymorphism from 8 to 21% of PAL plants has been maintained over the years at these two populations (Additional file 2: Table S2). In contrast, WAL individuals were found in nine of the 21 populations surveyed, including the two polymorphic populations, but always at very low frequencies (< 1% of total plants; Fig. 1 and Additional file 2: Table S2) and without any clear geographic pattern.
Flavonoid identification and composition in each plant tissue
Five anthocyanins and 21 flavones were identified in petals, as well as four anthocyanins and 19 flavones in photosynthetic tissues (Additional file 4: Table S3). The anthocyanins detected were cyanidin derivatives in all cases, but with different substituents in petals and photosynthetic tissues (Fig. 2 and Additional file 4: Table S3). The main anthocyanin present in pigmented petals was a glycosylated cyanidin with two sugars (one rhamnosylglucose and one glucose) and acylated with acetic acid (representing 71.0–74.1% of the total anthocyanin concentrations; peak 3 in Fig. 2c). In photosynthetic tissues, the structures of the predominant anthocyanins were simpler, with only one sugar attached to the aglycone (78.0–99.4%; peaks 6–9 in Fig. 2g).
The flavone composition was also different in petals compared to photosynthetic tissues: isovitexin derivatives were mostly accumulated in petals whereas isoorientin derivatives were the main flavones present in photosynthetic tissues (Table 1 and Additional file 4: Table S3). The primary petal flavone was an isovitexin glycosylated with two pentose sugars (69.5–88.3% of the total flavone concentrations in the three phenotypes; peak 17 in Fig. 2a and b). An isoorientin derivative containing an additional hexose and a caffeoyl residue was the main flavone present in calyces (58.8–63.2%), leaves (42.9–58.1%) and stems (50.6–57.8%; peak 35 in Fig. 2e and f). In all tissues, isoscoparin derivatives were also detected, but at relatively low levels (< 11%; Table 1 and Additional file 4: Table S3).
Flavones were the most abundant flavonoids detected across tissues. In fully pigmented individuals, concentrations of flavones were 14–56× higher than anthocyanins in petals, calyces and stems; leaf tissues showed an ever great bias, producing flavones at rates >100x that of anthocyanins (Table 1). These differences were even more apparent in white petals of PAL plants and anthocyanin-lacking WAL individuals.
Variation in flavonoid content among phenotypes using HPLC-DAD-MSn
The three phenotypes showed significant differences in their anthocyanin concentrations. In petals, PAL and WAL phenotypes accumulated only 1% of the same anthocyanins found in the fully pigmented phenotype (Fig. 2d; Tables 1 and 2). In calyces and stems, WAL phenotypes produced undetectable concentration of anthocyanins (Fig. 2h), whereas fully pigmented and PAL phenotypes showed similar anthocyanin levels. In leaves, anthocyanin concentration was very low and statistically similar for the three phenotypes.
In contrast to the differences found in anthocyanin concentrations, the three phenotypes showed minimal differences in their flavone content (Fig. 2a, b, e and f), and only three petal flavones (~ 1.5% of total flavones) were not present in all phenotypes (see Additional file 4: Table S3). We found differences in flavone composition between the polymorphic and non-polymorphic populations, with five compounds specific to Breña (compounds 16, 18, 21, 33a and 37a in Additional file 4: Table S3). The first three compounds were rare non-acylated O-glycosyl-C-monoglycoside flavones of petals (< 1% of total flavones), whereas the other two were moderately abundant in photosynthetic tissues (4.17–21.2% of total flavones). Thus, PCA based on the flavone composition and concentration showed higher separation between populations than among phenotypes of each population (Additional file 7: Figure S4). When we compared the flavone concentration of each specific group of flavones (i.e. derivatives of isovitexin, isoorientin and isoscoparin, and di-C-glycosides), we found no significant difference among phenotypes neither in the polymorphic nor the non-polymorphic population (Table 2). Similar results were obtained when using the relative proportion of flavones in MANOVAs (Additional file 8: Table S4).
Variation in flavonoid content among phenotypes measured spectrophotometrically
When expanding the sampling to more individuals and populations, spectrophotometric quantification showed similar pattern of anthocyanin and flavone production to that found in the HPLC analysis. In the two polymorphic populations, the three phenotypes showed significant differences for the anthocyanin accumulation in all tissues except for leaves, which showed very low values in all phenotypes (Fig. 3; Table 3). In petals, PAL and WAL phenotypes produced near zero anthocyanin concentration, whereas in photosynthetic tissues only the WAL phenotype lacked anthocyanins. The three phenotypes showed statistically similar flavone concentrations in photosynthetic tissues. In petals, significant differences were found due to the higher flavone content in WAL plants from Barra (Fig. 3). Between populations, significant differences were found for the anthocyanin production in petals and stems, and for the flavone production in all tissues. Anthocyanins and flavones concentrations were, in general, higher in Barra than Louro.
In non-polymorphic populations, anthocyanin production in petals and calyces were significantly different between fully pigmented and WAL phenotypes, and near significant in stems (Fig. 3; Table 3). Flavone concentrations in both phenotypes were similar in all tissues except for calyces, in which WAL plants from Breña showed nearly half concentrations compared to fully pigmented plants (Fig. 3). Between populations, significant differences were found for anthocyanins in all tissues except for the stems, and for flavones in all photosynthetic tissues, showing higher concentration levels in Breña population.
In this study, we found that shore campion accumulates both anthocyanins and flavones, but specific classes of these compounds were differentially produced in petals versus photosynthetic calyces, leaves and stems. Fully pigmented and PAL plants showed similar anthocyanin content in the analyzed tissues, except for the obvious absence in petals, whereas WAL phenotype lacks anthocyanins in the whole plant. In contrast, plants with white petals (both PAL and WAL phenotypes) have similar flavone composition and concentration compared to pink-flowered plants. Thus, the synthesis of flavones in each tissue of both PAL and WAL phenotypes seems to be not influenced by the loss of anthocyanins. This pattern of anthocyanin and flavone production in all phenotypes was congruent in the distant polymorphic and non-polymorphic populations. Together, these results suggest that anthocyanin accumulation in photosynthetic tissues are directly or indirectly involved in petal color polymorphism persistence. Below, we discuss these findings in view of the frequency in which PAL and WAL phenotypes are found in natural populations.
One of the most significant findings reported here is that PAL and WAL plants exhibited similar flavone content as fully pigmented plants, even though they lacked anthocyanins in either their petals or petals and vegetative tissues. In a previous study analyzing the sequences and gene expression of ABP genes in petals of S. littorea, Casimiro-Soriguer et al.  suggested that anthocyanin petal-loss in PAL individuals is caused by a decreased expression of flavanone-3-hydroxylase (F3h) controlled by a petal specific regulatory gene, SlMyb1a. Detection of flavones in PAL petals is consistent with the blockage of the ABP at F3 h since flavones are synthetized from naringenin or eriodictyol, which are produced in the steps immediately preceding F3H (Fig. 4). We suggest that downregulation of F3 h prevents a blockage of flavone production  and redirects flux from anthocyanins to flavones in white petals of S. littorea, as is described in other flower color polymorphic species [49, 50]. Similarly, mutations leading to the complete lack of anthocyanins in WAL individuals may occur later in the pathway to preserve flavone production. In Mimulus lewisii and Iochroma calycinum, mutations in coding regions of a late gene of the ABP, the dihydroflavonol 4-reductase (Dfr), cause the complete loss of anthocyanins of rare white-flowered individuals [51, 52]. A similar downstream blockage of the ABP, but not necessarily caused by inactivation of DFR, could explain the absence of anthocyanins in WAL individuals.
Although several studies have examined the genetic and molecular basis for PAL and WAL phenotypes [11, 29, 48, 51,52,53], this is the first study that has compared the complete flavonoid profile in these two forms of white-flowered individuals. Nevertheless, some studies have approached this goal in wild species. For example, flowers of the rare white-petal phenotype of Iochroma calycinum accumulate the same flavonol (quercetin) that the pigmented phenotype as determined by thin layer chromatography . Although comparisons between natural and horticultural plants should be addressed carefully since they are under different selection when in cultivation, there is some biochemical knowledge available for ornamental plants. Several studies have reported similar flavonol profiles for white-flowered and pigmented lines, as for lisianthus (Eustoma grandiflorum) and gentians (Gentiana triflora) [54, 55]; however, some white-flowered lines showed different flavonoid profiles, probably because of distinct genetic blockage of the ABP. We have found that all WAL plants of the shore campion showed the same flavonoid profile. This lack of flavonoid diversity in WAL phenotype reinforces the assumption that loss-of-function mutations may target specific late genes rather than the early genes of the ABP, which would compromise flavone production. The shore campion does not produce others classes of flavonoids others than flavones and anthocyanins . Thus, any mutation affecting an ABP gene prior to F3 h should preclude any flavone production (Fig. 4), making this mutation selectively disadvantageous due to the decisive role of these compounds to plant development and survival [28, 37, 56].
Loss of anthocyanin pigments is relatively common in nature [39, 40, 57], but its effects on plant fitness will determine the fate of white-flowered individuals. In several species, PAL phenotypes generally show similar, or even higher, fitness than fully pigmented plants ([40, 58, 59]; but see, for instance, ), resulting in stable flower color polymorphism in populations. In the shore campion, the pink-white polymorphism is maintained over the years, and white flowers represents 8–21% of total plants in the two polymorphic populations. Myb-mediated loss of anthocyanins, as for S. littorea, are frequently cell or tissue specific [60, 61] and allow downregulation of petal anthocyanins without hindering anthocyanin accumulation in other tissues. Our biochemical results confirmed that PAL and fully pigmented plants have similar anthocyanin content in photosynthetic tissues, in addition to have similar flavone content, which is expected to have few or no pleiotropic effects that could alter fitness of PAL phenotypes [17, 34]. In fact, snails and caterpillars produced similar herbivory levels in leaves of fully pigmented and PAL plants of Barra, but white petals of PAL plants showed higher hemipteran florivory that petals of fully pigmented plants (M.L.B. 2019, unpublished data).
A key question arising from our findings is, “Why are WAL phenotypes so rare even if they have similar amounts of protective flavones than fully pigmented plants?”. Flavones share many of the numerous protective biological functions attributed to anthocyanins [27, 56, 62], and are at least 14 times more abundant than anthocyanins across tissues of the shore campion, as commonly found in other species . In S. littorea, petal flavones (isovitexins) could be involved in regulation of vacuole homeostasis in epidermal cells and/or act as co-pigments of anthocyanins [64, 65], whereas flavones accumulated in photosynthetic tissues (isoorientins) are effective antioxidants that may play important functions in stress tolerance [28, 62]. Since these protective flavones are accumulated in reproductive and vegetative tissues of WAL plants, it seems plausible that the loss of anthocyanins could be involved in the ecological disadvantage by which WAL phenotype remains scarce in the populations. In fact, it is recently proposed that anthocyanins may play a decisive role in the regulation of signaling cascades responsible for cell growth and differentiation; thus, controlling important developmental processes [27, 28, 66, 67]. In addition, genetic linkage between ABP genes and other loci affecting fitness , as well as the metabolic cross-talk between flavonoid and other metabolic pathways , are other possible explanations for why loss of anthocyanins (not the flavonoid intermediates) seems to restrict the spread of WAL phenotypes. Taken together, our results suggest that the ability to produce anthocyanin pigments in photosynthetic tissues of the shore campion is associated with the ability to generate stable petal color polymorphisms.
In summary, our results show striking differences in the ability to synthesize anthocyanins between fully pigmented and white-petal variants of S. littorea, whereas flavone production is not affected by loss of anthocyanin pigments. Differences in flavonoid profile between PAL and WAL individuals are based in the absence of anthocyanins in petals or the whole plant, respectively. The low frequency of WAL plants in natural populations leads us to consider the negative ecological consequences of anthocyanin loss, or other putative pleiotropically-linked traits, in photosynthetic tissues and suggests the critical role of these compounds to generate stable flower color polymorphism. It will be interesting to analyze the flavonoid profile in other species with anthocyanin loss individuals either in the petals or in whole plant, and linking this information with their population frequency. These data will provide new insights into the flower color evolution, which also may produce new knowledge on microevolutionary processes.
Study species and frequency of anthocyanin-deficient phenotypes
Silene littorea is an annual wild plant that grows in coastal dune ecosystems from the northwestern to the southeastern Iberian Peninsula . Depending on the degree of human disturbance on their habitat, population’s size ranges from approx. 100 (e.g. Algezur, Odiel) to more than 10,000 individuals (e.g. Sines, Barra). This self-compatible species is entomophilous, but may produce up to 20% of fruit and seed set by spontaneous autogamy . Calyces, leaves and stems produce chlorophylls showing photosynthetic activity .
To account for the population frequency of PAL and WAL phenotypes, we sampled 21 populations covering the full distribution range of S. littorea (Fig. 1). In populations with more than 1000 estimated individuals, frequency estimates were carried out following a random sampling; in populations with lower number of individuals, the whole population was carefully sampled. Plants were visually categorized to each phenotype based on the absence or presence of anthocyanins in petals and photosynthetic tissues. PAL and WAL individuals showed a similar morphology (plant height, number of open flowers, flower size) than the fully pigmented phenotype.
Sampling for flavonoid analyses
Flavonoid analyses were conducted in two northwestern petal-color polymorphic populations (Barra and Louro) and in two southern non-polymorphic populations (Breña and Trafalgar; Fig. 1). For complete flavonoid identification, we randomly selected four plants of each fully pigmented, PAL and WAL phenotypes in Barra and four plants from each fully pigmented and WAL phenotypes in Breña. For spectrophotometric flavonoids quantification in these four populations, we extended the sampling to 11–20 plants of each phenotype except for the rare WAL individuals, in which six, five, six and three individuals were analyzed (Barra, Louro, Breña and Trafalgar, respectively). Sampling was carried out from March to April 2016, during the early to mid-flowering period of the species (permissions were not necessary to collect these samples). For each plant, we collected four samples: the petals and calyx of one flower, a leaf and stem sections (1 cm length) from the middle region of the stem. Flavonoids were extracted in 1.5 ml of MeOH containing 1% HCl and stored at − 20 °C in the dark, following the procedure described in Del Valle et al. . Voucher specimens from these populations are deposited at the Pablo de Olavide University Herbarium (UPOS-3954, UPOS-8982, UPOS-8983 and UPOS-8984).
Flavonoid identification and quantification by HPLC-DAD-MSn
A volume of 500 μL of methanolic extracts of each plant tissue (i.e. petals, calyces, leaves and stems) was concentrated in a SpeedVac concentrator (Savant ISS110, Thermo Fisher Scientific, NC, USA) after the addition of 100 μL of ultrapure water (Autwomatic, Wasserlab, Barbatáin, Spain). The volume of the aqueous extracts was then adjusted to 250 μL with acidified water (pH = 1.4, HCl). Aqueous extracts were filtered (Clarinert™ Syring Filters, 0.45 μm, Agela Technologies, DE, USA) prior to the HPLC-DAD-MSn analysis. We followed the procedure described in Alcalde-Eon et al. [71, 72], which has provided satisfactory results in the analyses of anthocyanins in other plant materials and also allows the simultaneous detection of anthocyanins and flavones (see results) in a single run. HPLC analyses were carried out in a Hewlett-Packard 1100 series liquid chromatograph (Agilent Technologies, Waldbronn, Germany). Detection was carried out at 360 nm and 520 nm for flavone and anthocyanin analysis, respectively. Spectra were recorded from 220 to 600 nm.
Mass spectrometric analyses were performed in an API 3200 Qtrap equipped with an ESI source and a triple quadrupole-ion trap mass analyzer that was controlled by Analyst v.5.1 software (Applied Biosystems, Darmstadt, Germany). The HPLC system was connected to the mass spectrometer via the UV cell outlet. Positive mode (ESI+) and specific conditions were selected to allow the simultaneous detection of anthocyanins and flavones (see Additional file 9: Supplementary methods for details).
Identification of the compounds was done considering their retention times, UV-Vis spectra, m/z of the molecular ion (M+) for anthocyanins or m/z of the protoned ion [M + H]+ for flavones, fragment ions and fragmentation patterns of the compounds. These data supplied valuable information concerning the nature of the aglycones and substituents, which was further compared to the features of standards reported in the literature  and to those standards (cyanidin 3-O-glucoside, isovitexin (apigenin 6-C-glucoside) and isoorientin (luteolin 6-C-glucoside)) and samples with known composition on barley (Hordeum vulgare L. ) analyzed in the same conditions. Furthermore, alkaline and acid hydrolyses were performed in all the types of sample extracts (petals and vegetative parts) to verify the presence of acids in the molecules as well as to determine their identity (see Additional file 9: Supplementary methods for details). Some of the major compounds of petals and photosynthetic tissues were isolated and alkaline hydrolysis was also carried out on them.
Anthocyanin quantification was done from the area of the peaks detected in the chromatogram recorded at 520 nm and using a Cyanidin 3-O-glucoside (Polyphenols Labs, Sandnes, Norway) calibration curve. Likewise, flavone quantification was done from the area of the peaks detected in the chromatogram recorded at 360 nm. A calibration curve of isovitexin and isoorientin (Extrasynthese, Genay, France) was employed to quantify the flavones present in petals and photosynthetic tissues, respectively.
HPLC results were analyzed separately for flavones and anthocyanins. In each tissue, exploratory analyses of flavone content of the different phenotypes were performed by principal-component analysis (PCA), using concentrations of all compounds identified. We retained those compounds with the highest principal component loadings. In addition, compounds highly correlated were eliminated to overcome co-linearity. PCA was based on the covariance matrix and without rotation of the extracted component . Confirmatory MANOVAs (multivariate analysis of variance) were carried out when differences between populations were detected . MANOVAs were performed grouping the total flavone composition of plants into four functional groups according to the C-glycoside flavone core (i.e. derivatives of isovitexin, isoorientin and isoscoparin, and di-C-glycosides; see Table 1). Since plants from Barra and Breña locations showed differences in their flavone content (see results) and the PAL phenotype is only present in Barra, both populations were analyzed independently. Because only a few anthocyanin compounds were found in the samples and some phenotypes did not produce anthocyanins (see results), exploratory PCAs could not be performed; instead, differences in the total anthocyanin concentrations among phenotypes were analyzed using ANOVAs with post hoc Tukey test. ANOVAs, PCAs and MANOVAs were carried out in SPSS v.22.0 (Armonk, NY, IBM Corp.).
Spectrophotometric flavonoids quantification
We used a Multiskan GO microplate spectrophotometer (Thermo Fisher Scientific Inc., MA, USA) to quantify global anthocyanin and flavone concentrations in order to expand the study to more individuals and populations. Three replicas of 200 μL were measured for of each tissues of all individuals. Absorbances were read at 350 and 520 nm to determine the concentrations of flavones and anthocyanins respectively , and their concentrations were calculated using calibration curves of standards of the main compounds found (i.e. cyanidin 3-O-glucoside, isovitexin and isoorientin) and expressed as cyanidin-3-glucoside, isovitexin and isoorientin equivalents in fresh weight, respectively.
Generalized linear models (GLMs) with Gaussian or gamma error distribution were performed in R v3.4.0  to test for differences in the accumulation of anthocyanins and flavones in each plant tissue; phenotype and population were considered as fixed factors. Previously, we tested the error distributions that generated the smaller deviance in the model, using the Akaike’s Information Criterion. F-test for analysis of deviance was used to correct for overdispersion . Multiple post hoc comparisons were performed using the “multcomp” R-package . Separated GLMs were performed for polymorphic and non-polymorphic populations, using post hoc comparisons with Bonferroni adjustment to test for differences among phenotypes within each population.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Anthocyanin biosynthetic pathway
Generalized linear models
- HPLC-DAD-MSn :
High-performance liquid chromatography coupled with diode-array detection and electrospray ionization tandem mass spectrometry
Petal anthocyanin loss
Wright S. Evolution in Mendelian populations. Genetics. 1931;16:97–159.
Herron JC, Freeman S. Evolutionary analysis. New York: Pearson; 2013.
Huxley J. Morphism and evolution. Heredity (Edinb). 1955;9:1–51.
Gigord LDB, Macnair MR, Smithson A. Negative frequency-dependent selection maintains a dramatic flower color polymorphism in the rewardless orchid Dactylorhiza sambucina (L.) Soò. Proc Natl Acad Sci. 2001;98:6253–5. https://doi.org/10.1073/pnas.111162598.
Hoballah ME, Gubitz T, Stuurman J, Broger L, Barone M, Mandel T, et al. Single gene-mediated shift in pollinator attraction in Petunia. Plant Cell. 2007;19:779–90. https://doi.org/10.1105/tpc.106.048694.
Mitchell-Olds T, Willis JH, Goldstein DB. Which evolutionary processes influence natural genetic variation for phenotypic traits? Nat Rev Genet. 2007;8:845–56.
Laine AL, Tellier A. Heterogeneous selection promotes maintenance of polymorphism in host-parasite interactions. Oikos. 2008;117:1281–8.
Chouteau M, Arias M, Joron M. Warning signals are under positive frequency-dependent selection in nature. Proc Natl Acad Sci. 2016;113:2164–9. https://doi.org/10.1073/pnas.1519216113.
Ford EB. Ecological Genetics. London: Chapman and Hall; 1975.
Wright S. Evolution and the genetics of populations. Volumen 4. Variability between and among natural populations. Chicago: University of Chicago Press; 1978.
Kellenberger RT, Byers KJRP, De Brito Francisco RM, Staedler YM, LaFountain AM, Schönenberger J, et al. Emergence of a floral colour polymorphism by pollinator-mediated overdominance. Nat Commun. 2019;10:63. https://doi.org/10.1038/s41467-018-07936-x.
Do R, Balick D, Li H, Adzhubei I, Sunyaev S, Reich D. No evidence that selection has been less effective at removing deleterious mutations in Europeans than in Africans. Nat Genet. 2015;47:126–31.
Sharp NP, Agrawal AF. An experimental test of the mutation-selection balance model for the maintenance of genetic variance in fitness components. Proc R Soc B Biol Sci. 2018;285:20181864.
Rausher MD. Evolutionary transitions in floral color. Int J Plant Sci. 2008;169:7–21.
Stern DL, Orgogozo V. The loci of evolution: how predictable is genetic evolution? Evolution. 2008;62:2155–77.
McKinnon JS, Pierotti MER. Colour polymorphism and correlated characters: genetic mechanisms and evolution. Mol Ecol. 2010;19:5101–25.
Streisfeld MA, Rausher MD. Population genetics, pleiotropy, and the preferential fixation of mutations during adaptive evolution. Evolution. 2010;65:629–42.
Hopkirk T. Flora Anomoia. John Smith & Son: London; 1817.
Mendel G. Versuche über Plflanzenhybriden. Verh natf ver Brunn abh. 1866;3:S3. https://doi.org/10.2307/2528198.
Coonen LP. Aristotle’s biology. Bioscience. 1977;27:733–8. https://doi.org/10.2307/1297678.
Clegg MT, Durbin ML. Flower color variation: a model for the experimental study of evolution. Proc Natl Acad Sci. 2000;97:7016–23. https://doi.org/10.1073/pnas.97.13.7016.
Fenster CB, Armbruster WS, Wilson P, Dudash MR, Thomson JD. Pollination syndromes and floral specialization. Annu Rev Ecol Evol Syst. 2004;35:375–403. https://doi.org/10.1146/annurev.ecolsys.34.011802.132347.
Grotewold E. The genetics and biochemistry of floral pigments. Annu Rev Plant Biol. 2006;57:761–80. https://doi.org/10.1146/annurev.arplant.57.032905.105248.
Harbone JB, Grayer RJ. Flavonoids and insects. In: Harborne JB, editor. The flavonoids. Advances in research since 1986. London: Chapman & Hall; 1994. p. 589–618.
Wilbert SM, Schemske DW, Bradshaw HD. Floral anthocyanins from two monkeyflower species with different pollinators. Biochem Syst Ecol. 1997;25:437–43.
Schemske DW, Bradshaw HD. Pollinator preference and the evolution of floral traits in monkeyflowers (Mimulus). Proc Natl Acad Sci. 1999;96:11910–5. https://doi.org/10.1073/pnas.96.21.11910.
Landi M, Tattini M, Gould KS. Multiple functional roles of anthocyanins in plant-environment interactions. Environ Exp Bot. 2015;119:4–17. https://doi.org/10.1016/j.envexpbot.2015.05.012.
Davies KM, Albert NW, Zhou Y, Schwinn KE. Functions of flavonoid and betalain pigments in abiotic stress tolerance in plants. Annu Plant Rev. 2018;1:1–41.
Coberly LC, Rausher MD. Pleiotropic effects of an allele producing white flowers in Ipomoea purpurea. Evolution. 2008;62:1076–85.
Sobel JM, Streisfeld MA. Flower color as a model system for studies of plant evo-devo. Front Plant Sci. 2013;4:321.
Wheldale M. The anthocyanin pigments of plants. Cambridge: Cambridge University Press; 1916.
Burdon JJ, Marshall DR, Brown AHD. Demographic and genetic changes in populations of Echium Plantagineum. J Ecol. 1983;71:667–79. https://doi.org/10.2307/2259584.
Levin DA, Brack ET. Natural selection against white petals in Phlox. Evolution. 1995;49:1017–22. https://doi.org/10.1111/j.1558-5646.1995.tb02336.x.
Wessinger CA, Rausher MD. Lessons from flower colour evolution on targets of selection. J Exp Bot. 2012;63:5741–9. https://doi.org/10.1093/jxb/err267.
Streisfeld MA, Liu D, Rausher MD. Predictable patterns of constraint among anthocyanin-regulating transcription factors in Ipomoea. New Phytol. 2011;191:264–74.
Chang SM, Lu Y, Rausher MD. Neutral evolution of the nonbinding region of the anthocyanin regulatory gene Ipmyb1 in Ipomoea. Genetics. 2005;170:1967–78.
Pollastri S, Tattini M. Flavonols: old compounds for old roles. Ann Bot. 2011;108:1225–33.
Strauss SY, Whittall JB. Non-pollinator agents of selection on floral traits. In: Harder LD, Barret HSC, editors. Ecology and evolution of flowers. Oxford University Press: Oxford; 2006. p. 120–38. https://doi.org/10.1111/j.1439-0418.1973.tb01782.x.
Warren J, Mackenzie S. Why are all colour combinations not equally represented as flower-colour polymorphisms? New Phytol. 2001;151:237–41.
Narbona E, Wang H, Ortiz PL, Arista M, Imbert E. Flower colour polymorphism in the Mediterranean Basin: occurrence, maintenance and implications for speciation. Plant Biol. 2018;20:8–20.
Ernst WHO. Scarcity of flower color polymorphism in field populations of Digitalis Purpurea L. Flora. 1987;179:231–9.
Twyford AD, Caola AM, Choudhary P, Raina R, Friedman J. Loss of color pigmentation is maintained at high frequency in a monkey flower population. Am Nat. 2018;191:135–45. https://doi.org/10.1086/694853.
Vaidya P, McDurmon A, Mattoon E, Keefe M, Carley L, Lee CR, et al. Ecological causes and consequences of flower color polymorphism in a self-pollinating plant (Boechera stricta). New Phytol. 2018;218:380–92.
Del Valle JC, Buide ML, Casimiro-Soriguer I, Whittall JB, Narbona E. On flavonoid accumulation in different plant parts: variation patterns among individuals and populations in the shore campion (Silene littorea). Front Plant Sci. 2015;6:939. https://doi.org/10.3389/fpls.2015.00939.
Berardi AE, Hildreth SB, Helm RF, Winkel BSJ, Smith SD. Evolutionary correlations in flavonoid production across flowers and leaves in the Iochrominae (Solanaceae). Phytochemistry. 2016;130:119–27.
Del Valle JC, Buide ML, Whittall J, Narbona E. Phenotypic plasticity in light-induced flavonoids varies among tissues in Silene littorea (Caryophyllaceae). Environ Exp Bot. 2018;153:100–7. https://doi.org/10.1016/j.envexpbot.2018.05.014.
Casimiro-Soriguer I, Narbona E, Buide ML, Del Valle JC, Whittall JB. Transcriptome and biochemical analysis of a flower color polymorphism in Silene littorea (Caryophyllaceae). Front Plant Sci. 2016;7:204. https://doi.org/10.3389/fpls.2016.00204.
Dick CA, Buenrostro J, Butler T, Carlson ML, Kliebenstein DJ, Whittall JB. Arctic mustard flower color polymorphism controlled by petal-specific downregulation at the threshold of the anthocyanin biosynthetic pathway. PLoS One. 2011;6:e18230.
Van Houwelingen A, Souer E, Spelt K, Kloos D, Mol J, Koes R. Analysis of flower pigmentation mutants generated by random transposon mutagenesis in Petunia hybrida. Plant J. 1998;13:39–50.
Nishihara M, Yamada E, Saito M, Fujita K, Takahashi H, Nakatsuka T. Molecular characterization of mutations in white-flowered torenia plants. BMC Plant Biol. 2014;14:86. https://doi.org/10.1186/1471-2229-14-86.
Wu CA, Streisfeld MA, Nutter LI, Cross KA. The genetic basis of a rare flower color polymorphism in Mimulus lewisii provides insight into the repeatability of evolution. PLoS One. 2013;8:e81173.
Coburn RA, Griffin RH, Smith SD. Genetic basis for a rare floral mutant in an andean species of solanaceae. Am J Bot. 2015;102:264–72.
Coberly LC, Rausher MD. Analysis of a chalcone synthase mutant in Ipomoea purpurea reveals a novel function for flavonoids: amelioration of heat stress. Mol Ecol. 2003;12:1113–24.
Davies KM, Bradley JM, Schwinn KE, Markham KR, Podivinsky E. Flavonoid biosynthesis in flower petals of five lines of lisianthus (Eustoma grandiflorum Grise.). Plant Sci. 1993;95:67–77.
Nakatsuka T, Nishihara M, Mishiba K, Yamamura S. Two different mutations are involved in the formation of white-flowered gentian plants. Plant Sci. 2005;169:949–58.
Jiang N, Doseff A, Grotewold E. Flavones: from biosynthesis to health benefits. Plants. 2016;5:27. https://doi.org/10.3390/plants5020027.
Carlson JE, Holsinger KE. Extrapolating from local ecological processes to genus-wide patterns in colour polymorphism in south African Protea. Proc R Soc B Biol Sci. 2015;282:20150583.
Gómez JM. Phenotypic selection and response to selection in Lobularia maritima: importance of direct and correlational components of natural selection. J Evol Biol. 2000;13:689–99.
Carlson JE, Holsinger KE. Direct and indirect selection on floral pigmentation by pollinators and seed predators in a color polymorphic south African shrub. Oecologia. 2013;171:905–19.
Albert NW, Davies KM, Lewis DH, Zhang HB, Montefiori M, Brendolise C, et al. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell. 2014;26:962–80.
Martins TR, Jiang P, Rausher MD. How petals change their spots: cis-regulatory re-wiring in Clarkia (Onagraceae). New Phytol. 2017;216:510–8.
Agati G, Azzarello E, Pollastri S, Tattini M. Flavonoids as antioxidants in plants: location and functional significance. Plant Sci. 2012;196:67–76. https://doi.org/10.1016/j.plantsci.2012.07.014.
Manetas Y. Why some leaves are anthocyanic and why most anthocyanic leaves are red? Flora. 2006;201:163–77.
van Brederode J, van Genderen HH, Berendsen W. Morphological effects of the flavone isovitexin in a non-glycosylating genotype of Silene pratensis (Caryophyllaceae). Experientia. 1982;38:929–31.
Yabuya T, Saito M, Iwashina T, Yamaguchi M. Stability of flower colors due to anthocyanin-flavone copigmentation in Japanese garden iris, Iris ensata Thunb. Euphytica. 2000;115:1–5.
Taylor LP, Grotewold E. Flavonoids as developmental regulators. Curr Opin Plant Biol. 2005;8:317–23.
Yao P, Deng R, Huang Y, Stael S, Shi J, Shi G, et al. Diverse biological effects of glycosyltransferase genes from Tartary buckwheat. BMC Plant Biol. 2019;19:339.
Mouradov A, Spangenberg G. Flavonoids: a metabolic network mediating plants adaptation to their real estate. Front Plant Sci. 2014;620. https://doi.org/10.3389/fpls.2014.00620.
Casimiro-Soriguer, I., Buide, M. L., and Narbona, E. (2013). The roles of female and hermaphroditic flowers in the gynodioecious-gynomonoecious Silene littorea, insights into the phenology of sex expression. Plant Biol. 15, 941–947. https://doi.org/10.1111/j.1438-8677.2012.00697.x
Narbona, E, Jaca, J., Del Valle, J. C., Valladares, F., and Buide,M. L. (2018). Whole-plant reddening in Silene germana is due to anthocyanin accumulation in response to visible light. Plant Biol. 20, 968-977. https://doi.org/10.1111/plb.12875.
Alcalde-Eon C, Rivas-Gonzalo JC, Muñoz O, Escribano-Bailón MT. Schizanthus grahamii and Schizanthus hookeri. Is there any relationship between their anthocyanin compositions and their different pollination syndromes? Phytochemistry. 2013;85:62–71. https://doi.org/10.1016/j.phytochem.2012.09.013.
Alcalde-Eon C, García-Estévez I, Rivas-Gonzalo JC. Rodríguez De La Cruz D, Escribano-Bailón MT. Anthocyanins of the anthers as chemotaxonomic markers in the genus Populus L.. differentiation between Populus nigra, Populus alba and Populus tremula. Phytochemistry. 2016;128:35–49.
Mabry TJ, Markham KR, Thomas MB. The systematic identification of flavonoids. New York: Springer-Verlag; 1970.
Ferreres F, Andrade PB, Valentão P, Gil-Izquierdo A. Further knowledge on barley (Hordeum vulgare L.) leaves O-glycosyl-C-glycosyl flavones by liquid chromatography-UV diode-array detection-electrospray ionisation mass spectrometry. J Chromatogr A. 2008;1182:56–64.
Valcárcel V, Vargas P. Quantitative morphology and species delimitation under the general lineage concept: optimization for Hedera (Araliaceae). Am J Bot. 2010;97:1555–73.
R Core Team. R: a language and environment for statistical computing. 2017. http://www.r-project.org.
Crawley MJ. The R Book. Chichester: Wiley; 2007.
Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008;50:346–63.
The authors thank Ágata Cardoso, Rebeca Ferreras and Antonio Gallardo for field and laboratory assistance. We also thank reviewers for their useful comments and suggestions.
Financial support was provided by the Spanish Government MINECO projects (CGL2012–37646 and CGL2015–63827-P) and also through a Predoctoral Training Program grant to JCV (BES-2013–062610).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Frequency of petal anthocyanin loss (PAL) and whole-plant anthocyanin loss (WAL) individuals in natural populations of species with polymorphism caused by loss of anthocyanins in petal and spontaneous white mutants, respectively.
Sites, geographical locations, estimated population size and percentage of petal anthocyanin loss (PAL) and whole-plant anthocyanin loss (WAL) phenotypes of S. littorea.
Pictures of white-flowered plants of S. littorea in which differences in anthocyanin accumulation in photosynthetic tissues between the petal anthocyanin loss (PAL) and whole-plant anthocyanin loss (WAL) phenotypes are shown.
Anthocyanins and flavones identified through HPLC-DAD-MSn from methanolic extracts of petals and photosynthetic tissues (calyces, leaves and stems) of S. littorea plants.
Chromatogram of the petal extract of a WAL specimen from Breña population recorded at 360 nm.
Chromatogram of the leaf extract of a fully pigmented specimen from Barra population recorded at 360 nm.
Scatter plot of principal components extracted from PCAs using flavone composition detected in petals, calyces, leaves and stems of S. littorea phenotypes through HPLC-DAD-MSn.
Results from MANOVAs comparing the relative proportion of flavones among phenotypes in Barra (fully pigmented, PAL and WAL) and Breña (fully pigmented and WAL).
Supplementary methods. Details of the mass spectrometry conditions, isolation of flavones, and alkaline and acid hydrolysis.
About this article
Cite this article
Del Valle, J.C., Alcalde-Eon, C., Escribano-Bailón, M.T. et al. Stability of petal color polymorphism: the significance of anthocyanin accumulation in photosynthetic tissues. BMC Plant Biol 19, 496 (2019). https://doi.org/10.1186/s12870-019-2082-6
- Flower color polymorphism
- Loss of pigmentation
- Non-pollinator mediated selection
- Plant secondary metabolites