- Research article
- Open Access
A highly mutable GST is essential for bract colouration in Euphorbia pulcherrima Willd. Ex Klotsch
BMC Genomics volume 22, Article number: 208 (2021)
Mutation breeding is an extraordinary tool in plant breeding to increase the genetic variability, where mutations in anthocyanin biosynthesis are targets to generate distinctive phenotypes in ornamental species. In poinsettia, ionizing radiation is routinely applied in breeding programs to obtaining a range of colours, with nearly all pink and white varieties being obtained after γ- or X-ray mutagenesis of red varieties. In the present study we performed a thorough characterization of a potential mutagenesis target gene as the main responsible for the ‘white paradox’ in poinsettia.
We identified a GST gene in poinsettia (Bract1) as an essential factor for the expression of anthocyanin-based red colouration of bracts, which presents a high phylogenetic similarity to known anthocyanin-related GSTs. Red poinsettia varieties and white mutants generated from these varieties by X-ray were analysed for polymorphisms related to the ‘white paradox’ in the species. A 4 bp mutation in a short repeat within the coding region of Bract1 is most likely responsible for the appearance of white phenotypes upon irradiation treatment. The polymorphism between wild-type and mutant alleles co-segregates with the phenotype in progeny from heterozygous red and white parents. Moreover, overexpression of Bract1 wild-type allele in Arabidopsis tt19 mutants restored the anthocyanin phenotype, while the Bract1 mutated allele showed to be non-functional.
The identified repeat seems to be highly unstable, since mutated plants can be easily detected among fewer than 200 shoots derived from 10 mutated plants. Our data indicate that particular short repeat sequences, similar to microsatellite sequences or so-called dynamic mutations, might be hot spots for genetic variability. Moreover, the identification of the Bract1 mutation fills a gap on the understanding on the molecular mechanism of colour formation in poinsettia.
Poinsettia, Euphorbia pulcherrima Willd. ex Klotsch, commonly known as Christmas Star, is an important ornamental crop, especially due to its association with Christmas time in North America, Europe, and Asia, with annual sales reaching nearly 150 million dollars in the USA . Its ornamental value is based on its intensely coloured bracts, which can be red, white, pink, or yellow or even have dual, scattered, or marbled colourations. Nonetheless, poinsettia breeding still focuses on red- and white-coloured varieties due to higher acceptance by consumers. In 2018, in Germany, approximately 80% of the poinsettias grown were red, 11% were white, and 9% were pink or had dual/scattered colouration .
Ionizing radiation is an important tool in mutation breeding for new colour variations in poinsettia, with nearly all pink and white varieties being obtained after gamma or X-ray mutagenesis of shoots of red varieties. Poinsettia mutation breeding is usually performed on cuttings that are irradiated with moderate dosages (~ 20 Gy), and mutants are selected on side shoots of the originally irradiated shoots. Flowering induction in the species occurs under short-day conditions and is accompanied by the development and colouration of bracts. Therefore, green leaves and red bracts occur concomitantly and accumulate different groups of pigments, i.e., chlorophylls and anthocyanins [53, 61]. Several anthocyanin types have been identified in poinsettia bracts and are responsible for its colouration range [3, 55, 66]; however, molecular information is still limited for the species [28, 72]. Nonetheless, genes responsible for the biosynthesis of the anthocyanin pathway have been intensively characterized in a range of species, with its regulation being highly dependent on R2R3-MYB regulatory genes and MYB-bHLH-WD40 (MBW) regulatory complexes [16, 58, 76].
Once synthesized on the cytoplasmic surface of the endoplasmic reticulum (ER), anthocyanin molecules need to be stored in the vacuole to prevent oxidation and loss of colour . Two main models of anthocyanin transport have been proposed: i) a vesicle trafficking-mediated model, where vesicle-like structures filled with anthocyanins are imported into the central vacuole via vesicle fusion [23, 27, 62]; and ii) a transporter-mediated model, where anthocyanins are carried across the vacuolar membrane by transport proteins (e.g., ABC and MATE transporters) with the help of glutathione S-transferase (GST) enzymes [26, 63, 78]. GSTs can bind to anthocyanin molecules to form a complex, thus escorting them from the ER to the vacuole, preventing oxidation [13, 54, 67, 78]. Anthocyanin-related GSTs play major roles in anthocyanin transport, since loss of function of these proteins leads to phenotypes with a lack of pigmentation, such as bz2 (Bronze-2) in maize, an9 (Anthocyanin 9) in petunia, tt19 (Transparent Testa 19) in Arabidopsis, fl3 (Flavonoid3) in carnation, riant (regulator involved in anthocyanin transport) in peach, and rap (reduced anthocyanin in petioles) in strawberry [2, 8, 38, 43, 48]. In our previous study, an anthocyanin-related GST-like gene showed higher expression in a red poinsettia variety than in the white counterpart, thus making it a promising candidate responsible for the so-called ‘white paradox’, e.g. appearance of acyanic (uncolored) phenotype despite the detection of expression of all structural genes and the related enzyme activities involved in the formation of red anthocyanin pigments .
In our current study, we identified an anthocyanin-related GST as the most likely target of the radiation-induced mutation of red poinsettias in white bract sports. Using different approaches, this study demonstrates the functionality of the poinsettia GST as an anthocyanin transporter. Most importantly, we show that a short repeat motif within the coding region of the gene is highly unstable upon mutation treatment, which leads to the high frequency of anthocyanin mutations observed in commercial mutation breeding. In addition to facilitating mutation breeding for bract colours, these results may be a starting point for analysing the genetic instability of short repeat sequences in plants.
Identification and characterization of Bract1
In a previous study , we observed higher expression of an anthocyanin-related GST-like gene (termed Bract1 hereafter) in the red poinsettia variety ‘Christmas Feelings’ than in its white counterpart ‘Christmas Feelings Pearl’. To investigate whether a similar phenomenon is observed in other red and white poinsettia pairs, we performed RT-qPCRs for six pairs of red-bracted poinsettia varieties and their independently generated white mutants. Normalized relative quantity (NRQ) values were calculated relative to one of the biological replicates of the ‘Chr. Glory’ variety according to the Pffafl method and equations . The levels of Bract1 expression varied among all varieties, with the varieties ‘Christmas Feelings’, ‘Titan’ and ‘SK130’ showing the highest relative expression. Although no lack of expression was observed in any of the white varieties, all red varieties showed significantly higher expression of Bract1 than their white counterparts (Fig. 1).
To further characterize the anthocyanin-related GST in poinsettia we sequenced the complete coding and intronic regions of the gene for the ‘Vintage’ variety. The final full-length (from first ATG start codon to stop codon) Bract1 sequence contains 2314 bp, with three exons (147 bp, 48 bp and 450 bp) and two introns (455 bp and 1214 bp) (Fig. 2a). The coding sequences (CDSs) of Bract1 from 14 red- and white-bracted poinsettia varieties (‘Noel’, ‘Valentino’, ‘Christmas Feelings’, ‘Christmas Feelings Pearl’, ‘Christmas Glory’, ‘Christmas Glory White’, ‘Joy’, ‘Joy White’, ‘Titan’, ‘Titan White’, ‘Bravo’, ‘Bravo White’, ‘SK130’ and ‘SK130 White’) were further sequenced. The sequence alignment showed high similarity in the CDS for all varieties, except for six single-nucleotide polymorphisms (SNPs) that were identified in both the ‘Chr. Glory’ and ‘Bravo’ varieties (Table 1). This result shows the presence of at least two allelic forms of the Bract1 gene. Additionally, a 4 bp deletion located 8 bp upstream of the first exon-intron junction was observed in all white varieties (Fig. 2b). The deletion is located in a short repeat locus, resembling a short simple sequence repeat (SSR), with a tetranucleotide motif ((CTTC)3) composition. The exact location of the (CTTC)3 motif is shown in Fig. 1a. The full-length gene sequence and CDS are available in Additional File S1.
Bract1 encodes a putative functional protein of 215 amino acids (aa) and a mass of 24.6 kDa, with distinctive GST components: a conserved GSH-binding site (G-site) located in the N-terminal domain and a C-terminal substrate-binding domain (H-site) . The predicted protein from the CDS containing the 4 bp deletion is a putative truncated protein with an early stop codon at position aa52 due to a frameshift in the mRNA. The full-length amino acid sequence and the truncated version are available in Additional File S1.
Bract colouration associated with a deletion in the Bract1 gene
The colour range in poinsettia varieties is obtained either through classic breeding (crossing) or mutagenic breeding (radiation), thus generating a spectrum of bract colours, such as pink, marble, orange and white/creamy. The white varieties are often obtained through radiation mutagenesis of the red varieties, followed by shoot development and trait selection. Therefore, red and white poinsettias from the same variety are referred to as ‘pairs’ due to their highly similar genetic background. However, not all red varieties can produce white sports through radiation. Therefore, red poinsettia varieties are distinguished into ‘heterozygous’ and ‘homozygous’ for the colouration locus according to their ability to generate white sports and according to the segregation of red and white phenotypes in progeny of crosses with white genotypes.
Since the 4 bp indel in the SSR locus of Bract1 had shown indications of polymorphism among the different poinsettia varieties—and a correlation with bract colouration—we used a genotyping approach based on the fluorescent labelling of PCR fragments. We genotyped 22 different poinsettia varieties bearing red and white bracts (Fig. 3a, Additional File S2). All the red heterozygous varieties showed two distinct copies of the allele (with and without the 4 bp deletion), while their white counterparts showed only the copy with the deletion. On the other hand, homozygous red varieties (i.e., those unable to generate white sports) showed only the copy without the deletion.
We further genotyped a segregating population with 190 progeny from ‘Joy’ (Rr) x ‘Joy White’ (rr) containing 36 white and 154 red plants (Fig. 3b, Additional File S2). Contrary to expectation, we observed a deviation in the segregation ratio, which was approximately 4:1 (red:white), instead of the expected 1:1 ratio for this crossing. This may be explained by the fact that seeds from white varieties are less vital than those from red varieties (von Tubeuf, Selecta One, pers. comm.). In addition, white varieties also exhibit lower pollen fertility, thus increasing the chances of self-pollination when red varieties are used as a female parent (von Tubeuf, Selecta One, pers. comm.). In fact, 17 red progeny showed only the wild-type copy of the allele (data not shown), which can be attributed only to self-pollination. Nonetheless, all the white progeny showed only the allele copy containing the deletion, thus reinforcing our hypothesis that the presence of the allele containing the deletion in a homozygous recessive state is correlated with the white phenotype.
Bract1 is the anthocyanin-related GST orthologue in poinsettia
As GST genes occur in large gene families, we wanted to analyse whether the poinsettia GST gene was related to other GST genes involved in anthocyanin transport to the vacuole. Therefore, we computed a phylogenetic tree from the deducted amino acid sequences of 95 GST family members from our previously assembled poinsettia transcriptome , as well as the Bract1 and anthocyanin-related GSTs from other species (CkmGST3, LcGST4, VvGST4, PhAN9, PpRiant1, PpRiant2, AtGSTF11 and AtTT19). Nine GST classes were identified among the poinsettia GSTs: Tau, Theta, Lambda, Zeta, Phi, tetrachlorohydroquinone dehalogenase (TCHQD), glutathionyl hydroquinone reductase (GHR), dehydroascorbate reductase (DHAR) and eukaryotic translation elongation factor 1B-γ (Ef1Bγ). Except Tau and Ef1Bγ, all other GST classes showed a single cluster (Fig. 4). All anthocyanin-related GSTs belong to the Phi class and clustered together in the phylogenetic tree, with Bract1 showing high similarity with these GSTs.
By aligning the Bract1 nucleotide CDSs with those of anthocyanin-related GSTs from other species, an overall nucleotide similarity of 61.9% was observed (Additional File S3). Protein alignment of BRACT1 with the other anthocyanin-related GSTs resulted in an overall similarity of 58.3%, with the peach RIANT1 protein showing the highest similarity (66.5%) (Additional File S3). Interestingly, we identified seven amino acid residues, previously reported as specific to anthocyanin-related GSTs [32, 37, 40], that are conserved in the protein alignment, except in AtGSTF11: 2Val, 11Ala, 13Cys, 62Phe, 90Leu, 91Glu and 152Ser (Fig. 5). In summary, these results indicate that Bract1 is the anthocyanin-related GST orthologue in poinsettia.
Bract1 functionally complements the Arabidopsis tt19 mutant phenotype
To examine the in vivo function of Bract1 as an anthocyanin transporter, we tested the ability of Bract1 cDNA to functionally complement the Arabidopsis GST mutant tt19, which is defective in the expression of anthocyanins in aboveground organs and seeds. Two constructs containing the Bract1 cDNA (with and without the 4 bp deletion) under the cauliflower mosaic virus (CaMV) 35S promoter were introduced into the tt19 mutant by the floral-dip method [11, 75]. Although the constructs contained a GFP marker for the selection of transgenic events, we genotyped 10 independent biological replicates from the T2 progeny of tt19/35S::Bract1 and tt19/35S::Bract1_mut transgenic plants. All progeny contained the correct allele from the Bract1 gene, thus confirming the correct integration of the transgenic construct (Fig. 6a, Additional File S2).
Upon stimulation of anthocyanin accumulation in seedlings by irradiation with red/blue LEDs, the tt19/35S::Bract1 transgenic lines displayed a purple hypocotyl phenotype at the seedling stage, similar to the Columbia (Col-0) line but not the tt19 mutant (Fig. 6b). On the other hand, tt19/35S::Bract1 transgenic lines did not show complementation of the anthocyanin phenotype. The ProAtUbi::GFP construct, used as a control for infiltration, did not result in any phenotypic changes. Moreover, transgenic plants harbouring Bract1 did not complement the seed colour of tt19, as the seed colour at the ripening stage remained the same as that of the mutant tt19 in transgenic plants (Fig. 6c). This finding suggests that Bract1 may have distinct functions from TT19 during seed coat pigmentation. Taken together, these results not only emphasize the role of Bract1 in anthocyanin transport in poinsettia but also demonstrate that a deletion in its coding region leads to a colourless phenotype.
De novo mutations occur with high frequency and include deletion of the 4 bp repeat
To study the stability of the 4 bp repeat within the first exon of the Bract1 gene, we analysed DNA samples from mutation experiments conducted over the last 4 years at Selecta One. In brief, 10 cuttings from the varieties ‘Aurora’, ‘SK159 Dark Pink’, ‘Aurora Jingle’ and ‘SK183’ were X-ray irradiated with 20 Gy (30 Gy for ‘SK183), and subsequently, side shoots from those cuttings were further propagated. DNA was extracted and analysed as previously described from 377, 191, 188 and 186 of the propagated side shoots. Table 2 shows the results indicating that out of 942 samples, 9 mutations could be detected. Three mutated progeny were identified in both ‘SK159 Dark Pink’ and ‘Aurora Jingle’ individuals and two in the ‘SK183’ individuals, and only one mutated individual was identified in the ‘Aurora’ progeny. Unfortunately, as this was part of a commercial breeding programme, individual shoots were not labelled in a way that would allow tracing them back to one of the original shoots that were irradiated. However, even if all the mutations detected in each of the separate mutation treatments were redundant and originated from one original mutational event, the frequency was extraordinarily high.
Microsatellite repeats are not an anthocyanin-related feature
The microsatellite repeat present in the Bract1 gene shows signs of instability upon irradiation treatment. To identify whether such repeats are a common feature for anthocyanin-related GSTs in Euphorbiaceae or related taxa or are a family-specific feature, we first computed a phylogenetic tree from the CDSs of Bract1, known anthocyanin-related GSTs (CkmGST3, LcGST4, VvGST4, PhAN9, PpRiant1, PpRiant2, and AtTT19) and GST-like orthologues from Euphorbiaceae species (Euphorbia esula, Euphorbia pekinensis, Ricinus communis, Jatropha curcas, Hevea brasiliensis and Manihot esculenta). Figure 7 shows that Bract1 shared high similarity with the GSTs from the two Euphorbia species (E. esula and E. pekinensis) but also closely clustered with the GSTs from the other Euphorbiaceae species. Although none of the GST-like genes from other Euphorbiaceae have been investigated as putative anthocyanin transporters, they may perform a similar function due to their homology with Bract1 and the other anthocyanin-related GSTs.
Furthermore, we wanted to assess the distribution of the microsatellite repeat and possibly investigate its origin. The sequence alignment shows that none of the GSTs analysed contain the same (CTTC)3 SSR motif observed in Bract1 (Fig. 8). When analysing the GSTs from Euphorbiaceae species, two CTTC repeat units were present in E. pekinensis, with a single nucleotide substitution responsible for the loss of the third repeat. J. curcas, H. brasiliensis and M. esculenta showed two CTTC repeat units, with two nucleotide substitutions related to the loss of the third repeat. Last, E. esula and R. communis showed a single CTTC unit and three nucleotide substitutions in the microsatellite region. The anthocyanin-related GSTs showed greater distinction in the microsatellite region. Three of the genes contained a single CTTC repeat unit (CkmGST3, AtTT19 and LcGST4), while the others did not contain any repeat unit (VvGST4, PpRiant1, PpRiant2 and PhAN9), which was due to the higher number of nucleotide substitutions (ranging from 4 to 7 substitutions). In conclusion, even though some of the known anthocyanin-related GSTs contain a CTTC sequence, CTTC repeats seem to be a common feature of Euphorbiaceae GSTs, and the three repeats from the E. pulcherrima GST (Bract1) behave in an unstable manner upon being subjected to ionizing irradiation.
Anthocyanins, a class of flavonoid secondary metabolite compounds , are responsible for providing orange to blue colours in plant tissues, and their biosynthetic and regulatory mechanisms have been widely characterized . However, there is still debate on the mechanism of anthocyanin transport from the ER to the vacuole [62, 63]. Strong evidence for the involvement of transport proteins [21, 24], with a special role of GST enzymes, has been reported in several plant species (Alfenito et al. [38, 43];). In the present study, we demonstrated that the Bract1 gene functions as an anthocyanin transporter in poinsettia and that a highly mutable repeat in its coding region leads to frequent deletions and therefore to a colour-deficient phenotype.
Bract1 is a functional GST gene related to anthocyanin transport in poinsettia
GSTs are a large and diverse group of enzymes with multifunctional roles, especially in the detoxification of xenobiotics as well as in responses to biotic and abiotic stresses [1, 15]. The classification of GSTs is based on sequence conservation, genomic organization, and physiochemical properties, among other features [18, 33, 42]. Based on our previous study , we identified 95 GST genes in poinsettia and phylogenetically classified them into nine different classes based on their similarity with known Arabidopsis GSTs (Fig. 4). To date, 14 GST classes have been identified in plants: tau (U), phi (F), lambda (L), DHAR, theta (T), zeta (Z), EF1Bγ, TCHQD, microsomal prostaglandin E-synthase type 2 (mPGES-2), GHR, metaxin, Ure2p, hemerythrin (H) and iota (I) (reviewed by ).
A large number of GSTs have been identified in plant species, such as 49 in Capsella rubella , 55 in Arabidopsis (), 61 in Citrus  and 139 in L. chinensis . Bract1 clusters with high bootstrap support with anthocyanin-related GSTs from other species (e.g., AtTT19, PhAN9 and VvGST4), with all of these GSTs belonging to the phi class. Known anthocyanin-related GSTs belong almost exclusively to the phi class, except for Bronze-2 from maize, which belongs to the tau class . Further support for Bract1 being a member of the phi-type plant GST genes is provided by the presence of two introns as a characteristic of this group of genes, such as AN9 from petunia and TT19 from Arabidopsis (Alfenito et al. ;).
Complementation studies using Arabidopsis tt19 mutants have been widely applied as proof of concept for the function of GSTs as anthocyanin transporters (Alfenito et al. [32, 34, 37, 40, 46, 54, 57];). Due to the high amino acid conservation of GST enzymes involved in flavonoid accumulation among species , they can complement each other’s anthocyanin-deficient mutants (Alfenito et al. ;). However, similar to our observation for Bract1, not all of these genes complemented both the shoot and seed phenotypes [34, 40, 48]. TT19 is involved in both anthocyanin accumulation in vegetative tissues and proanthocyanidin (PA) accumulation in Arabidopsis seed coats, which provides its brown colouration . Transgenic tt19 mutants overexpressing the petunia AN9 orthologue showed AN9 mRNA expression in developing siliques, but the seed colour still remained the same as the wild-type mutant . Altogether, these results suggest that GST orthologues from different species may have distinct functions from TT19 during seed-coat pigmentation.
A direct complementation of poinsettia white mutants with the functional Bract1 would ultimately prove its function in bract colouration. However, neither Agrobacterium-mediated infiltration nor biolistic particle delivery system (a.k.a. gene gun) were successful for transient expression studies (data not shown). Stable transformation in poinsettia have been done using electrophoresis-based methods [9, 71], but no stable transgenic poinsettia was obtained. Successful stable transformation via Agrobacterium-mediated infiltration has been previously achieved, but the process is time-consuming . Attempts to perform stable transformation of poinsettia with Bract1 alleles will bridge the current knowledge gap but are out of the scope of the present study.
A loss-of-function mutation in Bract1 is the cause of the “white paradox” in poinsettia
Based on our results, we hypothesize that deletion of one unit of the repeat in the Bract1 gene is responsible for most of the white genotypes in poinsettia. This hypothesis is strongly supported by the evidence that the tt19/35S::Bract1_mut mutant was not able to complement the anthocyanin phenotype in the Arabidopsis tt19 mutant, unlike the tt19/35S::Bract1 mutant. Mutations in GSTs leading to colourless phenotypes have been previously reported. A mutation in the fl3 gene in carnation leads to a light pink phenotype, but a brighter phenotype is observed upon complementation by petunia AN9 and maize Bz2 . In peach, four alleles of a GST gene (Riant) were identified, with two of them containing frameshift mutations and unable to complement the Arabidopsis tt19 phenotype. Varieties containing copies of the mutated alleles in a homozygous state showed flowers with white variegated phenotypes . Last, a single-nucleotide polymorphism (SNP) in the strawberry RAP gene, leading to a premature stop codon, results in a mutant with green petioles and leaves. The non-functional rap gene was not able to complement Arabidopsis tt19, while wild-type RAP was successful .
In our analyses, all six independently generated white mutants of red varieties displayed the same deletion of a 4 bp repeat in Bract1, whereas the original varieties all contained a fully functional copy of the gene. In addition, co-segregation of the deletion with the white phenotype was observed in a segregated population of 190 progeny. Furthermore, a novel mutation leading to a homozygous recessive allele of Bract1 among 184 samples obtained from irradiated cuttings of the heterozygous line SK183 led to a white phenotype (Table 3). Altogether, the results of this study present strong evidence that the four-base deletion in Bract1 is the cause of the red-to-white shift in the poinsettia varieties analysed here. However, as anthocyanin biosynthesis involves several steps, other regulatory and structural genes might give rise to white mutants as well, as has been shown in numerous other examples [5, 39, 49, 52]. We did not detect these genes in our current plant material perhaps due to the much higher mutation rate of the Bract1 gene than of less mutable genes.
Bract1 contains a short highly mutable four-base repeat
Upon X-ray treatment, red poinsettia plants produce progeny bearing white phenotypes with high frequencies, often based on only 10 irradiated cuttings (von Tubeuf, Selecta One, pers. comm., Selecta One). This phenomenon is associated with a deletion in a short repeat in the Bract1 gene of white mutants in a homozygous state. The mutations in all six independent mutant pairs that we detected are exactly identical, which indicates that the X-ray treatment did not directly cause the mutation but rather led to changes indirectly by stimulating the DNA repair mechanisms via replication errors, by increasing recombination or by the other mechanisms discussed for mutations in repeat sequences . The possible involvement of replication-based errors is supported by our observation that upon amplification of the repeat via standard PCR from cloned Bract1 wild-type or mutant allele, a low level of variants carrying four-base indels can always be detected (data not shown).
Radiation is frequently used as a tool for mutagenic breeding in poinsettia. In contrast to ethyl-methanesulphonate (EMS)-based chemical mutagenesis, which produces point mutations with high frequency , ionizing radiation (e.g., X-rays and γ-rays) induces DNA oxidative damage, such as double-strand breaks (DSBs), base substitutions, deletions and chromosomal alterations, at a lower frequency, frequently resulting in loss of gene function [35, 36, 51]. SSRs are among the most variable types of repetitive sequences in the genome . Studies have shown that SSR instability increases with plant development  and abiotic stress . This might be another explanation for the frequent observation of repeat changes in the Bract1 gene after X-ray irradiation, although the small number of repeats (i.e., three) of four base pairs each does not fit the most widely applied criteria used to define SSRs, which usually focus on sequences with a larger number of repeats.
However, little information about the genetics and dynamics is available for short repeats. A majority of studies compared historical events for mostly shorter SSRs (2 and 3 bp repeats with larger repeat numbers) in present-day populations or the dynamic repeats responsible for human diseases (mostly trinucleotide repeats), which usually display effects beyond those of large numbers of repeats (> 30 repeats ).
Our observation that a large number of mutation events could be observed in the side shoots of ten irradiated plants indicates an unusually high mutation rate, which is in contrast to the few reports in which exact mutation rates have been reported for vegetatively propagated crops . In one example, the woody ornamental plant Tibouchina urvelliana was irradiated three independent times with a 45 Gy dose, resulting in 0.06% dwarf mutants each time . However, several authors reported that the radiosensitivity of vegetative tissues varies greatly among species and tissues , so exact comparative estimations of mutation frequencies have a very limited accuracy among species and conditions. However, experiments with transgenic Arabidopsis lines harbouring constructs designed to analyse restoration of GUS open reading frames by either recombination or by restoring in-frame translation by mutations in SSRs demonstrate the occurrence of easily detectable numbers of somatic mutation events [22, 74]. Together with the careful selection of side shoots after X-ray irradiation of poinsettia, this finding may explain the high rate of recessive mutations detected here.
In this study, we showed that the poinsettia Bract1 gene is an active GST gene involved in the expression of anthocyanins in poinsettia bracts. Furthermore, a 4 bp deletion in a short repeat within the coding region of Bract1 is the most likely cause of many mutations that lead to a white bract colour. This mutation occurs with an unusually high frequency and is presumably an indirect effect of X-ray mutagenesis. Future analyses using mutagenesis in transgenic Arabidopsis lines harbouring Bract1 might help elucidate the causes of the high instability of this repeat. Moreover, this result might also serve as a reference for the study of other repeat-containing structural genes as potential mutational hot spots in plant genomes.
In this study, we showed that the poinsettia Bract1 gene is an active GST gene involved in the expression of anthocyanins in poinsettia bracts. Furthermore, a four base pair deletion in a short repeat within the coding region of Bract1 is the most likely cause of many white mutations for bract colour. This mutation occurs with an unusually high frequency and is presumably an indirect effect of the x-ray mutagenesis. Future analyses using mutagenesis in transgenic Arabidopsis lines harbouring the Bract1 might help to elucidate the causes of the high instability of this repeat. Moreover, it might also serve as an example for other repeat containing structural genes in plant genomes as potential mutational hot spots.
A range of red- and white-bracted varieties of poinsettia was used in this study for the different analyses (Table 3). Among the varieties, we used so-called pairs of red varieties and white mutants generated from these varieties by γ-ray or X-ray mutagenesis by the company Selecta One (Stuttgart, Germany). In addition, a segregating population containing 190 progeny from the “Joy x Joy white” cross, one such pair, was also used for analysis. The plant material used are exclusively cultivated varieties of the species Euphorbia pulcherrima which were provided by the company Selecta One. Identification of the material was conducted by Dr. Robert Boehm and Mr. Guido von Tubeuf. There were no vouchers taken and deposited.
Leaf and bract samples for DNA and RNA isolation were harvested at Selecta One, immediately frozen in liquid nitrogen, shipped on dry ice and stored at − 80 °C at the Institute for Plant Genetics of the Leibniz Universität Hannover (Hannover, Germany). Moreover, rooted cuttings of three red (‘Aurora’, ‘Aurora Jingle’ and ‘SK183’) and one pink (‘SK159 Dark Pink’) variety were irradiated with γ-rays (20 Gy) and further developed, and new cuttings were generated by the company Selecta One. The progeny were used for fragment analysis (please refer to section “Poinsettia genotyping and fragment analysis”).
Seeds of the Arabidopsis thaliana Columbia (col-0) genotype were available at the Institute for Plant Genetics of the Leibniz Universität Hannover; these seeds were originally obtained from the Arabidopsis information service in Frankfurt, Germany , and subsequently propagated in isolated greenhouse compartments. Seeds of the mutant line tt19–8 (stock number: CS2105587) were obtained from the Arabidopsis Biological Resource Center (ABRC). Seeds were sown in Einheitserde P substrate, and seedlings were placed in long-day conditions (16 h light/8 h dark, 22 °C) for 2 weeks to induce flowering. Seedlings of the wild-type and mutants analysed for the expression of anthocyanins were grown under supplemental red/blue LED light (GP LED production DR/B 120 LB, Philips, Germany) to stimulate anthocyanin development.
DNA and RNA isolation
For the poinsettia samples, DNA was isolated from approximately 100 mg of leaf tissue using the NucleoSpin® Plant II Kit (Macherey–Nagel GmbH & Co., KG, Düren, Germany) according to the manufacturer’s instructions. Total RNA was isolated from approximately 100 mg of bract tissue using the mirPremier™ miRNA Isolation Kit (Sigma-Aldrich, St. Louis, USA). For Arabidopsis samples, total RNA was isolated from approximately 50 mg of leaf tissue using the Quick-RNA Plant Kit (Zymo Research, Irvine, USA). cDNA synthesis was performed using the FastGene Scriptase Basic cDNA Kit (Nippon Genetics Europe GmbH, Düren, Germany) according to the manufacturer’s recommendations. The DNA and total RNA concentrations and quality were assessed using a NanoDrop™ 2000 (Thermo Fisher Scientific, Wilmington, USA) and gel electrophoresis.
GST expression by RT-qPCR
Two endogenous reference genes (Translation elongation factor 1 beta (EF1B) and Translation elongation factor 1 alpha (EF1A)) were used to normalize the Bract1 expression data. Primer sequences are available in Additional File S4. The amplification efficiency for all primers was obtained from relative standard curves. Three independent biological replicates were used for each of the varieties. RT-qPCRs were performed using the qPCRBIO SyGreen Mix Lo-ROX Kit (Nippon Genetics Europe GmbH) according to the manufacturer’s recommendations. Briefly, reactions were carried out in technical triplicates in a volume of 10 μL containing 5 μL of qPCRBIO SyGreen Mix Lo-ROX, 10 μmol of gene-specific forward and reverse primers, and 4 μL of a 1:50 cDNA dilution. RT-qPCRs were performed using a StepOne™ Real-Time PCR System (Applied Biosystems, Singapore, Singapore). The normalized relative quantity (NRQ) was calculated according to the Pfaffl equations . The real-time data for this study are provided according to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines . Pairwise statistical analysis between each red variety and its white counterpart was performed using the Relative Expression Software Tool (REST) v2.0.13 .
Bract1 gene sequencing
The poinsettia variety ‘Vintage’ was used for full-length sequencing of the GST-like gene (hereafter named Bract1). PCRs were performed in a 50 μL reaction containing 50 ng of DNA template, 1X PrimeSTAR® Buffer (Mg2+ plus), 0.2 mM each dNTP, 0.25 μM forward and reverse primers and 1.25 U of PrimeSTAR® HS DNA Polymerase (Takara Bio Inc., Kusatsu, Japan). The cycling conditions were 95 °C for 3 min; 30 cycles of 95 °C for 30 s, 60 °C for 30 s and 72 °C for 2 min; and a final extension of 10 min at 72 °C. The PCR products were resolved in a 1% (w/v) agarose gel by horizontal electrophoresis for 90 min at 100 V. The correct bands were excised from the gel and purified using the NucleoSpin® Gel and PCR Clean-up Kit (Macherey–Nagel) following the manufacturer’s recommendations. Finally, the purified PCR fragments were sent to Eurofins Genomics (Ebersberg, Germany) for Sanger sequencing. The generated sequences were aligned using the ClustalW function in the BioEdit Sequence Alignment Editor v7.2.5 , and a final full-length gene sequence for Bract1 was generated.
The coding sequences (CDSs) of Bract1 from 14 red- and white-bracted poinsettia varieties (‘Noel’, ‘Valentino’, ‘Christmas Feelings’, ‘Christmas Feelings Pearl’, ‘Christmas Glory’, ‘Christmas Glory White’, ‘Joy’, ‘Joy White’, ‘Titan’, ‘Titan White’, ‘Bravo’, ‘Bravo White’, ‘SK130’ and ‘SK130 White’) were further sequenced. PCRs were performed in a 50 μL reaction containing 1 μL of undiluted cDNA, 1X PrimeSTAR® Buffer (Mg2+ plus), 0.2 mM each dNTP, 0.25 μM forward and reverse primers and 1.25 U of PrimeSTAR® HS DNA Polymerase (Takara). The cycling conditions were 95 °C for 3 min; 30 cycles of 95 °C for 30 s, 60 °C for 30 s and 72 °C for 45 s; and a final extension of 10 min at 72 °C. The sequencing strategy was the same as that used for the full-length sequencing analysis. The generated sequences were aligned using the ClustalW function in the BioEdit Sequence Alignment Editor v7.2.5. Primer sequences are available in Additional File S4.
Poinsettia genotyping and fragment analysis
To detect changes in the repeat structure of the Bract1 gene, a genotyping approach based on the fluorescent labelling of PCR fragments  was applied. DNA samples were PCR amplified in a 20 μL reaction containing 50 ng of DNA template, 1X Williams buffer, 0.15 mM each dNTP, 0.0125 μM forward primer, 0.07 μM universal FAM-labelled M13 primer, 0.25 μM reverse primer and 1 U of DCSPol DNA Polymerase (DNA Cloning Service, Hamburg, Germany). The cycling conditions were 94 °C for 3 min; 24 cycles of 94 °C for 45 s, 59 °C for 1 min and 72 °C for 1 min; 6 cycles of 94 °C for 30 s, 52 °C for 45 s and 72 °C for 1 min; and a final extension of 10 min at 72 °C. Fifty microliters of formamide loading dye (98% formamide, 10 mM EDTA, 0.05% pararosanilin) was added to each reaction, which was then incubated at 95 °C for 5 min. The PCR products were resolved in a 6% (w/v) acrylamide gel via vertical electrophoresis using a LI-COR Gene Readir 4200 DNA Analyser (LI-COR Biosciences, Nebraska, USA). The varieties from Table 3 and the progeny of the irradiated variety ‘Aurora’ were genotyped using fluorescent labelling of PCR fragments. The progeny of the irradiated varieties ‘Aurora Jingle’ and ‘SK159 Dark Pink’ were analysed by capillary electrophoresis on an ABI 3730 XL system at Microsynth AG (Balgach, Switzerland). Primer sequences are available in Additional File S4.
Protein sequences for Bract1 and 95 different poinsettia GSTs, retrieved from our previous study , were predicted with TransDecoder  and used for the construction of a phylogenetic tree. Moreover, protein sequences for known anthocyanin-related GSTs from other species were included in the analysis: CkmGST3 (Cyclamen persicum x Cyclamen purpurascens, GenBank - AB682678.1), LcGST4 (Litchi chinensis, GenBank - KT946768.1), VvGST4 (Vitis vinifera, GenBank - AY971515.1), PhAN9 (Petunia hybrida, GenBank - Y07721.1), PpRiant1 (Prunus persica, GenBank - KT312847.1), PpRiant2 (P. persica, GenBank - KT312848.1), AtGSTF11 (Arabidopsis thaliana, GenBank - NM_111189.3) and AtTT19 (A. thaliana, GenBank - NM_121728.4). The putative protein sequences of all the GSTs are available in Additional File S5.
Sequence alignment was performed using MUSCLE , and the phylogenetic tree was constructed with MEGA X v10.0.5  using the maximum likelihood (ML) method with the Whelan and Goldman matrix-based model using a discrete gamma distribution (WAG+G) . The best model was estimated using MEGAX. The tree topology was tested via a bootstrap analysis with 1000 replicates. For better visualization of the phylogenetic tree, Tree Of Life (iTOL) software, version 4.2.3  (https://itol.embl.de/), was used.
Plasmid construction and Agrobacterium-mediated infiltration
The coding sequence of Bract1 was amplified from the poinsettia varieties ‘Vintage’ and ‘Christmas Feelings Pearl’ to capture both wild-type and mutated alleles. The primers used for amplification are available in Additional File S4. The PCR fragments were inserted in the sense orientation into the BamHI-HindIII site of the C757pGFPU10–35 s-ocs-LH (Additional File S6) binary vector (DNA Cloning Service, Hamburg, Germany) using the In-Fusion® HD Cloning Kit (Takara Bio Inc., Kusatsu, Japan). The vector contains a GFP gene under the control of the Arabidopsis ubiquitin promoter (ProAtUbi::GFP) and a 35S promoter upstream of the multiple cloning site into which inserts were cloned. The final expression vectors contained either the wild-type allele (35S::Bract1) or the mutated allele (35S::Bract1_mut). The vector containing only GFP was also used for transformation as a negative control (only ProAtUbi::GFP).
The expression vectors were introduced into Agrobacterium tumefaciens strain GV3101 via electroporation. Transformation of the A. thaliana mutant line tt19–8 was performed using the floral dip method [11, 75]. For transgenic plant selection, T0 seeds were sown in soil, and GFP-expressing seedlings were selected to produce T1 and subsequently T2 progeny to achieve GST homozygosity. T2 seedlings of tt19/35S::Bract1 and tt19/35S::Bract1_mut transgenic plants were used for phenotypic analysis. Non-transformed tt19–8, Col-0 and tt19/ProAtUbi::GFP seedlings were used as controls. Seedlings were placed under red light to stimulate anthocyanin biosynthesis. To confirm correct GST integration, 10 independent biological replicates, each representing an independently selected transgenic line from the T2 progeny of tt19/35S::Bract1 and tt19/35S::Bract1_mut transgenic plants, were used for GST genotyping. The protocol was the same as that used in the section “Poinsettia genotyping and fragment analysis”.
Analysis of the GST repeat
To understand the origin of the microsatellite-like repeat in the Bract1 gene, orthologous GST genes from Euphorbiaceae species were retrieved by BLASTN against the Euphorbiaceae (taxid: 3977) nucleotide database. GST-like genes from the Euphorbiaceae species Ricinus communis (GenBank - XM_002532928.3), Manihot esculenta (GenBank - XM_021748071.1), Jatropha curcas (GenBank - XM_012219312.2), Hevea brasiliensis (GenBank - XM_021787187.1), Euphorbia esula (GenBank - PJAE01736713.1) and Euphorbia pekinensis [7, 72], as well the anthocyanin-related GSTs CkmGST3 (C. persicum x C. purpurascens, GenBank - AB682678.1), LcGST4 (L. chinensis, GenBank - KT946768.1), VvGST4 (V. vinifera, GenBank - AY971515.1), PhAN9 (P. hybrida, GenBank - Y07721.1), PpRiant1 (P. persica, GenBank - KT312847.1), PpRiant2 (P. persica, GenBank - KT312848.1) and AtTT19 (A. thaliana, GenBank - NM_121728.4) were used to construct a phylogenetic tree.
Sequence alignment was performed using MUSCLE , and the phylogenetic tree was constructed with MEGA X v10.0.5  using the ML method with the Tamura 3-parameter matrix-based model  using a discrete gamma distribution with invariant sites (T92 + G + I). The best model was estimated using MEGAX. The tree topology was tested via a bootstrap analysis with 1000 replicates. For better visualization of the phylogenetic tree, Tree Of Life (iTOL) software, version 4.2.3  (https://itol.embl.de/), was used.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files. DNA sequence from the Bract1 gene is available in the GenBank repository under accession number MW718861. The CDS sequences from several poinsettia genotypes generated during the current study are available in the GenBank repository under accession numbers MW718847 - MW718860.
ATP binding cassette transporter
Cauliflower mosaic virus
Green fluorescent protein
Gray (unit of ionizing radiation)
Insertion and deletion variations
Multi-antimicrobial extrusion protein
- MBW complex:
Normalized relative quantity
Polymerase chain reaction
Quantitative reverse transcription PCR
Single nucleotide polymorphism
Single sequence repeat
Transparent testa 19 Arabidopsis mutant
Agrawal GK, Jwa N-S, Rakwal R. A pathogen-induced novel rice (Oryza sativa L.) gene encodes a putative protein homologous to type II glutathione S-transferases. Plant Sci. 2002;163(6):1153–60. https://doi.org/10.1016/S0168-9452(02)00331-X.
Alfenito MR, Souer E, Goodman CD, Buell R, Mol J, Koes R, et al. Functional Complementation of Anthocyanin Sequestration in the Vacuole by Widely Divergent Glutathione S-Transferases. Plant Cell. 1998;10(7):1135-49.
Asen S. Anthocyanins in bracts of Euphorbia Pulcherrima as revealed by paper chromatographic and spectrophotometric methods. Plant Physiol. 1958;33(1):14–7. https://doi.org/10.1104/pp.33.1.14.
Bassham DC. Pigments on the move. Nature. 2015;526(7575):644–5. https://doi.org/10.1038/526644a.
Ben-Simhon Z, Judeinstein S, Trainin T, Harel-Beja R, Bar-Ya’akov I, Borochov-Neori H, et al. A “white” anthocyanin-less pomegranate (Punica granatum L.) caused by an insertion in the coding region of the Leucoanthocyanidin Dioxygenase (LDOX; ANS) gene. PLOS ONE. 2015;10:e0142777. https://doi.org/10.1371/journal.pone.0142777.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22. https://doi.org/10.1373/clinchem.2008.112797.
Cao X, Zhang F, Yuan B, Meng L, Yang X, Jiang J. De novo transcriptome sequencing and analysis of Euphorbia pekinensis Rupr. And identification of genes involved in diterpenoid biosynthesis. Plant Gene. 2017;12:33–42. https://doi.org/10.1016/j.plgene.2017.07.001.
Cheng J, Liao L, Zhou H, Gu C, Wang L, Han Y. A small indel mutation in an anthocyanin transporter causes variegated colouration of peach flowers. J Exp Bot. 2015;66(22):7227–39. https://doi.org/10.1093/jxb/erv419.
Clarke JL, Klemsdal SS, Floistad E, Hvoslef-Eide AK, Haugslien S, Moe R, et al. Genetic engineering of poinsettia with the aim of enhancing its resistance to poinsettia mosaic virus. Acta Hortic. 2006:321–5. https://doi.org/10.17660/ActaHortic.2006.722.40.
Clarke JL, Spetz C, Haugslien S, Xing S, Dees MW, Moe R, Blystad DR. Agrobacterium tumefaciens-mediated transformation of poinsettia, Euphorbia pulcherrima, with virus-derived hairpin RNA constructs confers resistance to poinsettia mosaic virus. Plant Cell Rep. 2008;27(6):1027–38. https://doi.org/10.1007/s00299-008-0526-9.
Clough SJ, Bent AF. Floral dip: a simplified method forAgrobacterium-mediated transformation ofArabidopsis thaliana: floral dip transformation of Arabidopsis. Plant J. 1998;16(6):735–43. https://doi.org/10.1046/j.1365-313x.1998.00343.x.
Debener T, Lehnackers H, Arnold M, Dangl JL. Identification and molecular mapping of a single Arabidopsis thaliana locus determining resistance to a phytopathogenic Pseudomonas syringae isolate. Plant J Cell Mol Biol. 1991;1(3):289–302. doi:10.1046/j.1365-313X.1991.t01-7-00999.x.
Dixon DP, Edwards R. Glutathione Transferases. Arab Book. 2010;8:e0131 doi:10.1199/tab.0131.
Dixon DP, Lapthorn A, Edwards R. Plant glutathione transferases. Genome Biol. 2002;3, reviews3004.1. https://doi.org/10.1186/gb-2002-3-3-reviews3004.
Dixon DP, Skipsey M, Edwards R. Roles for glutathione transferases in plant secondary metabolism. Phytochemistry. 2010;71(4):338–50. https://doi.org/10.1016/j.phytochem.2009.12.012.
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15(10):573–81. https://doi.org/10.1016/j.tplants.2010.06.005.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. https://doi.org/10.1093/nar/gkh340.
Edwards R, Dixon DP, Walbot V. Plant glutathione S-transferases: enzymes with multiple functions in sickness and in health. Trends Plant Sci. 2000;5(5):193–8. https://doi.org/10.1016/s1360-1385(00)01601-0.
Ellegren H. Microsatellites: simple sequences with complex evolution. Nat. Rev. Genet. 2004;5(6):435–45. https://doi.org/10.1038/nrg1348.
Esnault M-A, Legue F, Chenal C. Ionizing radiation: advances in plant response. Environ Exp Bot. 2010;68(3):231–7. https://doi.org/10.1016/j.envexpbot.2010.01.007.
Francisco RM, Regalado A, Ageorges A, Burla BJ, Bassin B, Eisenach C, Zarrouk O, Vialet S, Marlin T, Chaves MM, Martinoia E, Nagy R. ABCC1, an ATP binding cassette protein from grape berry, transports anthocyanidin 3-O-Glucosides. Plant Cell. 2013;25(5):1840–54. https://doi.org/10.1105/tpc.112.102152.
Golubov A, Yao Y, Maheshwari P, Bilichak A, Boyko A, Belzile F, Kovalchuk I. Microsatellite instability in Arabidopsis increases with plant Development1[W][OA]. Plant Physiol. 2010;154(3):1415–27. https://doi.org/10.1104/pp.110.162933.
Gomez C, Conejero G, Torregrosa L, Cheynier V, Terrier N, Ageorges A. In vivo grapevine anthocyanin transport involves vesicle-mediated trafficking and the contribution of anthoMATE transporters and GST: anthocyanin trafficking in grapevine. Plant J. 2011;67(6):960–70. https://doi.org/10.1111/j.1365-313X.2011.04648.x.
Goodman CD, Casati P, Walbot V. A multidrug resistance–associated protein involved in anthocyanin transport in Zea mays. Plant Cell. 2004;16(7):1812–26. https://doi.org/10.1105/tpc.022574.
Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, Reynolds SH, Enns LC, Burtner C, Johnson JE, Odden AR, Comai L, Henikoff S. Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics. 2003;164(2):731–40.
Grotewold E, editor. The science of flavonoids. New York: Springer; 2006. https://doi.org/10.1007/978-0-387-28822-2.
Grotewold E, Davies K. Trafficking and sequestration of Anthocyanins. Nat Prod Commun. 2008;3, 1934578X0800300. https://doi.org/10.1177/1934578X0800300806.
Gu Z, Chen H, Yang R, Ran M. Identification of DFR as a promoter of anthocyanin accumulation in poinsettia ( Euphorbia pulcherrima , willd. ex Klotzsch) bracts under short-day conditions. Sci Hortic. 2018;236:158–65. https://doi.org/10.1016/j.scienta.2018.03.032.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, MacManes MD, Ott M, Orvis J, Pochet N, Strozzi F, Weeks N, Westerman R, William T, Dewey CN, Henschel R, LeDuc RD, Friedman N, Regev A. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512. https://doi.org/10.1038/nprot.2013.084.
Hall, T. A., Hall, T., Hall, T. A., Hall, T., Hall, T. A., Hall, T. A., et al. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Doi: 10.14601/Phytopathol_Mediterr-14998u1.29.
He G, Guan C-N, Chen Q-X, Gou X-J, Liu W, Zeng Q-Y, Lan T. Genome-wide analysis of the glutathione S-Transferase gene family in Capsella rubella: identification, expression, and biochemical functions. Front Plant Sci. 2016;7. https://doi.org/10.3389/fpls.2016.01325.
Hu B, Zhao J, Lai B, Qin Y, Wang H, Hu G. LcGST4 is an anthocyanin-related glutathione S-transferase gene in Litchi chinensis Sonn. Plant Cell Rep. 2016;35(4):831–43. https://doi.org/10.1007/s00299-015-1924-4.
Islam S, Sajib SD, Jui ZS, Arabia S, Islam T, Ghosh A. Genome-wide identification of glutathione S-transferase gene family in pepper, its classification, and expression profiling under different anatomical and environmental conditions. Sci Rep. 2019;9(1):9101. https://doi.org/10.1038/s41598-019-45320-x.
Jiang S, Chen M, He N, Chen X, Wang N, Sun Q, Zhang T, Xu H, Fang H, Wang Y, Zhang Z, Wu S, Chen X. MdGSTF6, activated by MdMYB1, plays an essential role in anthocyanin accumulation in apple. Hortic Res. 2019;6(1):40. https://doi.org/10.1038/s41438-019-0118-6.
Jo YD, Kim J-B. Frequency and Spectrum of radiation-induced mutations revealed by whole-genome sequencing analyses of plants. Quantum Beam Sci. 2019;3(2):7. https://doi.org/10.3390/qubs3020007.
Kazama Y, Hirano T, Saito H, Liu Y, Ohbu S, Hayashi Y, Abe T. Characterization of highly efficient heavy-ion mutagenesis in Arabidopsis thaliana. BMC Plant Biol. 2011;11(1):161. https://doi.org/10.1186/1471-2229-11-161.
Kitamura S, Akita Y, Ishizaka H, Narumi I, Tanaka A. Molecular characterization of an anthocyanin-related glutathione S-transferase gene in cyclamen. J Plant Physiol. 2012;169(6):636–42. https://doi.org/10.1016/j.jplph.2011.12.011.
Kitamura S, Shikazono N, Tanaka A. TRANSPARENT TESTA 19 is involved in the accumulation of both anthocyanins and proanthocyanidins in Arabidopsis. Plant J. 2004;37(1):104–14. https://doi.org/10.1046/j.1365-313X.2003.01943.x.
Koseki M, Goto K, Masuta C, Kanazawa A. The star-type color pattern in Petunia hybrida ‘red star’ flowers is induced by sequence-specific degradation of Chalcone synthase RNA. Plant Cell Physiol. 2005;46(11):1879–83. https://doi.org/10.1093/pcp/pci192.
Kou M, Liu Y, Li Z, Zhang Y, Tang W, Yan H, Wang X, Chen XG, Su ZX, Arisha MH, Li Q, Ma DF. A novel glutathione S-transferase gene from sweetpotato, IbGSTF4, is involved in anthocyanin sequestration. Plant Physiol Biochem. 2019;135:395–403. https://doi.org/10.1016/j.plaphy.2018.12.028.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9. https://doi.org/10.1093/molbev/msy096.
Lallement P-A, Brouwer B, Keech O, Hecker A, Rouhier N. The still mysterious roles of cysteine-containing glutathione transferases in plants. Front Pharmacol. 2014;5. https://doi.org/10.3389/fphar.2014.00192.
Larsen ES, Alfenito MR, Briggs WR, Walbot V. A carnation anthocyanin mutant is complemented by the glutathione S-transferases encoded by maize Bz2 and petunia An9. Plant Cell Rep. 2003;21(9):900–4. https://doi.org/10.1007/s00299-002-0545-x.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5. https://doi.org/10.1093/nar/gkw290.
Licciardello C, D’Agostino N, Traini A, Recupero G, Frusciante L, Chiusano M. Characterization of the glutathione S-transferase gene family through ESTs and expression analyses within common and pigmented cultivars of Citrus sinensis (L) Osbeck. BMC Plant Biol. 2014;14(1):39. https://doi.org/10.1186/1471-2229-14-39, Characterization of the glutathione S-transferase gene family through ESTs and expression analyses within common and pigmented cultivars of Citrus sinensis (L.) Osbeck.
Liu Y, Qi Y, Zhang A, Wu H, Liu Z, Ren X. Molecular cloning and functional characterization of AcGST1, an anthocyanin-related glutathione S-transferase gene in kiwifruit (Actinidia chinensis). Plant Mol Biol. 2019;100(4-5):451–65. https://doi.org/10.1007/s11103-019-00870-6.
Liu Y, Tikunov Y, Schouten RE, Marcelis LFM, Visser RGF, Bovy A. Anthocyanin biosynthesis and degradation mechanisms in Solanaceous vegetables: a review. Front Chem. 2018;6:52. https://doi.org/10.3389/fchem.2018.00052.
Luo H, Dai C, Li Y, Feng J, Liu Z, Kang C. Reduced Anthocyanins in petioles codes for a GST anthocyanin transporter that is essential for the foliage and fruit coloration in strawberry. J Exp Bot. 2018;69(10):2595–608. https://doi.org/10.1093/jxb/ery096.
Luo P, Ning G, Wang Z, Shen Y, Jin H, Li P, et al. Disequilibrium of Flavonol synthase and Dihydroflavonol-4-Reductase expression associated tightly to white vs. red color flower formation in plants. Front. Plant Sci. 6. 2016. https://doi.org/10.3389/fpls.2015.01257.
Marrs KA, Alfenito MR, Lloyd AM, Walbot V. A glutathione S-transferase involved in vacuolar transfer encoded by the maize gene Bronze-2. Nature. 1995;375(6530):397–400. https://doi.org/10.1038/375397a0.
Morita R, Kusaba M, Iida S, Yamaguchi H, Nishio T, Nishimura M. Molecular characterization of mutations induced by gamma irradiation in rice. Genes Genet Syst. 2009;84(5):361–70. https://doi.org/10.1266/ggs.84.361.
Morita Y, Saito R, Ban Y, Tanikawa N, Kuchitsu K, Ando T, Yoshikawa M, Habu Y, Ozeki Y, Nakayama M. Tandemly arranged chalcone synthase a genes contribute to the spatially regulated expression of siRNA and the natural bicolor floral phenotype in Petunia hybrida: endogenous siRNA and bicolored flowers of Petunia. Plant J. 2012;70(5):739–49. https://doi.org/10.1111/j.1365-313X.2012.04908.x.
Moustaka J, Panteris E, Adamakis I-DS, Tanou G, Giannakoula A, Eleftheriou EP, Moustakas M. High anthocyanin accumulation in poinsettia leaves is accompanied by thylakoid membrane unstacking, acting as a photoprotective mechanism, to prevent ROS formation. Environ Exp Bot. 2018;154:44–55. https://doi.org/10.1016/j.envexpbot.2018.01.006.
Mueller LA, Goodman CD, Silady RA, Walbot V. AN9, a Petunia glutathione S -Transferase required for anthocyanin sequestration, Is a Flavonoid-Binding Protein. Plant Physiol. 2000;123(4):1561–70. https://doi.org/10.1104/pp.123.4.1561.
Nitarska D, Stefanini C, Haselmair-Gosch C, Miosic S, Walliser B, Mikulic-Petkovsek M, Regos I, Slatnar A, Debener T, Terefe-Ayana D, Vilperte V, Hadersdorfer J, Stich K, Halbwirth H. The rare orange-red colored Euphorbia pulcherrima cultivar “harvest Orange” shows a nonsense mutation in a flavonoid 3′-hydroxylase allele expressed in the bracts. BMC Plant Biol. 2018;18(1):216. https://doi.org/10.1186/s12870-018-1424-0.
Pearson CE, Edamura KN, Cleary JD. Repeat instability: mechanisms of dynamic mutations. Nat Rev Genet. 2005;6(10):729–42. https://doi.org/10.1038/nrg1689.
Pérez-Díaz R, Madrid-Espinoza J, Salinas-Cornejo J, González-Villanueva E, Ruiz-Lara S. Differential roles for VviGST1, VviGST3, and VviGST4 in Proanthocyanidin and anthocyanin transport in Vitis vinífera. Front Plant Sci. 2016;7. https://doi.org/10.3389/fpls.2016.01166.
Petroni K, Tonelli C. Recent advances on the regulation of anthocyanin synthesis in reproductive organs. Plant Sci. 2011;181(3):219–29. https://doi.org/10.1016/j.plantsci.2011.05.009.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:45e–45. https://doi.org/10.1093/nar/29.9.e45.
Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30:10.
Pomar F, Ros Barceló A. Are red leaves photosynthetically active? Biol Plant. 2007;51(4):799–800. https://doi.org/10.1007/s10535-007-0164-z.
Poustka F, Irani NG, Feller A, Lu Y, Pourcel L, Frame K, Grotewold E. A trafficking pathway for Anthocyanins overlaps with the endoplasmic reticulum-to-vacuole protein-sorting route in Arabidopsis and contributes to the formation of vacuolar inclusions. Plant Physiol. 2007;145(4):1323–35. https://doi.org/10.1104/pp.107.105064.
Saito K, Yonekura-Sakakibara K, Nakabayashi R, Higashi Y, Yamazaki M, Tohge T, Fernie AR. The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol Biochem. 2013;72:21–34. https://doi.org/10.1016/j.plaphy.2013.02.001.
Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000;18(2):233–4. https://doi.org/10.1038/72708.
Schum, A., and Preil, W. (1998). “Induced Mutations in Ornamental Plants,” in Somaclonal Variation and Induced Mutations in Crop Improvement Current Plant Science and Biotechnology in Agriculture. S. M. Jain, D. S. Brar, B. S. Ahloowalia (Dordrecht: Springer Netherlands), 333–366. doi:https://doi.org/10.1007/978-94-015-9125-6_17.
Slatnar A, Mikulic-Petkovsek M, Veberic R, Stampar F, Schmitzer V. Anthocyanin and chlorophyll content during poinsettia bract development. Sci Hortic. 2013;150:142–5. https://doi.org/10.1016/j.scienta.2012.10.014.
Sun Y, Li H, Huang J-R. Arabidopsis TT19 functions as a carrier to transport anthocyanin from the cytosol to Tonoplasts. Mol Plant. 2012;5(2):387–400. https://doi.org/10.1093/mp/ssr110.
Tamura K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C-content biases. Mol Biol Evol. 1992;9:678–87. https://doi.org/10.1093/oxfordjournals.molbev.a040752.
USDA (2018). Floriculture crops 2018 summary. National Agricultural Statistics Service. Available at: https://www.nass.usda.gov/Publications/Todays_Reports/reports/floran19.pdf.
van der Ploeg, R. (2018). FloraCulture international. Int. Assoc. Hortic. Prod. IAHP. Available at: https://www.floraculture.eu/digital-magazine/fci-november-december-2018/ [Accessed March 18, 2020].
Vik NI, Hvoslef-Eide AK, Gjerde H, Bakke K. Stable transformation of pionsettia via electrophoresis. Acta Hortic. 2001:101–3. https://doi.org/10.17660/ActaHortic.2001.560.12.
Vilperte V, Lucaciu CR, Halbwirth H, Boehm R, Rattei T, Debener T. Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts. BMC Genomics. 2019;20:900. https://doi.org/10.1186/s12864-019-6247-3.
Whelan S, Goldman N. A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001;18(5):691–9. https://doi.org/10.1093/oxfordjournals.molbev.a003851.
Yao Y, Kovalchuk I. Abiotic stress leads to somatic and heritable changes in homologous recombination frequency, point mutation frequency and microsatellite stability in Arabidopsis plants. Mutat Res. 2011;707(1-2):61–6. https://doi.org/10.1016/j.mrfmmm.2010.12.013.
Zhang X, Henriques R, Lin S-S, Niu Q-W, Chua N-H. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nat Protoc. 2006;1(2):641–6. https://doi.org/10.1038/nprot.2006.97.
Zhao D, Tao J. Recent advances on the development and regulation of flower color in ornamental plants. Front Plant Sci. 2015;6. https://doi.org/10.3389/fpls.2015.00261.
Zhao J. Flavonoid transport mechanisms: how to go, and with whom. Trends Plant Sci. 2015;20(9):576–85. https://doi.org/10.1016/j.tplants.2015.06.007.
Zhao J, Dixon RA. The ‘ins’ and ‘outs’ of flavonoid transport. Trends Plant Sci. 2010;15(2):72–80. https://doi.org/10.1016/j.tplants.2009.11.006.
We would like to thank Guido von Tubeuf for providing the poinsettia varieties.
This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 675657 Flower Power. This funding body had no role in the design of the study, collection, analysis, or interpretation of data, or in writing the manuscript. The publication of this article was funded by the Open Access fund of Leibniz Universität Hannover. Open Access funding enabled and organized by Projekt DEAL.
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.
Additional file 1
Full-length sequence of the Bract1 gene, CDS sequences of red and white poinsettia varieties, full sequence and truncated version of the BRACT1 protein.
Additional file 2.
Additional file 3
Sequence similarity of Bract1 with other anthocyanin-related GSTs.
Additional file 4.
List of primers used in each of the analyses in the present study.
Additional file 5
Deducted protein sequences from the Bract1 and 95 GSTs from E. pulcherrima GSTs, as well as anthocyanin-related GSTs from other species.
Additional file 6
Schematic representation from the C757pGFPU10–35 s-ocs-LH binary vector used for the constructions of transformation plasmids containing either the wild-type allele (35S::Bract1) or the mutated allele (35S::Bract1_mut) from the poinsettia GST.
Additional file 7
Sequence alignment of Bract1, anthocyanin-related GSTs and orthologue GSTs from Euphorbiaceae species.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Vilperte, V., Boehm, R. & Debener, T. A highly mutable GST is essential for bract colouration in Euphorbia pulcherrima Willd. Ex Klotsch. BMC Genomics 22, 208 (2021). https://doi.org/10.1186/s12864-021-07527-z
- Euphorbia pulcherrima
- Ionizing radiation
- Glutathione S-transferase
- Mutation breeding
- Short repeat sequences