Transcription factor FveMYB117a inhibits axillary bud outgrowth by regulating cytokinin homeostasis in woodland strawberry.

Shoot branching affects plant architecture. In strawberry (Fragaria L.), short branches (crowns) develop from dormant axillary buds to form inflorescences and flowers. While this developmental transition contributes greatly to perenniality and yield in strawberry, its regulatory mechanism remains unclear and understudied. In the woodland strawberry (Fragaria vesca), we identified and characterized 2 independent mutants showing more crowns. Both mutant alleles reside in FveMYB117a, a R2R3-MYB transcription factor gene highly expressed in shoot apical meristems, axillary buds, and young leaves. Transcriptome analysis revealed that the expression of several cytokinin pathway genes was altered in the fvemyb117a mutant. Consistently, active cytokinins were significantly increased in the axillary buds of the fvemyb117a mutant. Exogenous application of cytokinin enhanced crown outgrowth in the wild type, whereas the cytokinin inhibitors suppressed crown outgrowth in the fvemyb117a mutant. FveMYB117a binds directly to the promoters of the cytokinin homeostasis genes FveIPT2 encoding an isopentenyltransferase and FveCKX1 encoding a cytokinin oxidase to regulate their expression. Conversely, the type-B Arabidopsis response regulators FveARR1 and FveARR2b can directly inhibit the expression of FveMYB117a, indicative of a negative feedback regulation. In conclusion, we identified FveMYB117a as a key repressor of crown outgrowth by inhibiting cytokinin accumulation and provide a mechanistic basis for bud fate transition in an herbaceous perennial plant.

The red/far-red light photoreceptor FvePhyB regulates tissue elongation and anthocyanin accumulation in woodland strawberry.

Light is an important environmental signal that influences plant growth and development. Among the photoreceptors, phytochromes can sense red/far-red light to coordinate various biological processes. However, their functions in strawberry are not yet known. In this study, we identified an EMS mutant, named P8, in woodland strawberry (Fragaria vesca) that showed greatly increased plant height and reduced anthocyanin content. Mapping-by-sequencing revealed that the causal mutation in FvePhyB leads to premature termination of translation. The light treatment assay revealed that FvePhyB is a bona fide red/far-red light photoreceptor, as it specifically inhibits hypocotyl length under red light. Transcriptome analysis showed that the FvePhyB mutation affects the expression levels of genes involved in hormone synthesis and signaling and anthocyanin biosynthesis in petioles and fruits. The srl mutant with a longer internode is caused by a mutation in the DELLA gene FveRGA1 (Repressor of GA1) in the gibberellin pathway. We found that the P8 srl double mutant has much longer internodes than srl, suggesting a synergistic role of FvePhyB and FveRGA1 in this process. Taken together, these results demonstrate the important role of FvePhyB in regulating plant architecture and anthocyanin content in woodland strawberry.

Brassinosteroid catabolic enzyme CYP734A129 regulates the morphologies of leaves and floral organs in woodland strawberry.

Brassinosteroids (BRs) play critical roles in plant growth and development and regulate many important agronomic traits. However, the functions of BRs in strawberry are unclear. This study identified two mutants, named P6 and R87, in woodland strawberry (Fragaria vesca) from EMS mutagenesis populations that exhibit narrow leaves, petals and sepals. Mapping by sequencing and genetic studies revealed that the F. vesca CYP734A129, encoding a putative BR catabolic enzyme, is the causative gene for both P6 and R87. Overexpression of CYP734A129 in both F. vesca and Arabidopsis causes a severe dwarf phenotype, and the BRI1-EMS-SUPPRESSOR 1 (BES1) protein is less abundant in the CYP734A129-overexpressing Arabidopsis seedlings. This suggests that CYP734A129 is functionally conserved with CYP734A1, as a BR-inactivating enzyme. Transcriptome analysis of young leaves revealed that four BR biosynthetic genes were significantly downregulated in P6 (cyp734a129), and photosynthesis-related genes were highly enriched among the up-regulated genes in P6 compared to the wild type. This further supports that CYP734A129 inactivates BRs in F. vesca. Furthermore, we showed that mutations in CYP734A129 do not affect fruit shape and color during ripening in strawberry. Overall, our results suggest that F. vesca CYP734A129 is a BR catabolic enzyme, and provide insights into the roles of CYP734A129 in strawberry.

Molecular bases of strawberry fruit quality traits: Advances, challenges, and opportunities.

The strawberry is one of the world's most popular fruits, providing humans with vitamins, fibers, and antioxidants. Cultivated strawberry (Fragaria × ananassa) is an allo-octoploid and highly heterozygous, making it a challenge for breeding, quantitative trait locus (QTL) mapping, and gene discovery. Some wild strawberry relatives, such as Fragaria vesca, have diploid genomes and are becoming laboratory models for the cultivated strawberry. Recent advances in genome sequencing and CRISPR-mediated genome editing have greatly improved the understanding of various aspects of strawberry growth and development in both cultivated and wild strawberries. This review focuses on fruit quality traits that are most relevant to the consumers, including fruit aroma, sweetness, color, firmness, and shape. Recently available phased-haplotype genomes, single nucleotide polymorphism (SNP) arrays, extensive fruit transcriptomes, and other big data have made it possible to locate key genomic regions or pinpoint specific genes that underlie volatile synthesis, anthocyanin accumulation for fruit color, and sweetness intensity or perception. These new advances will greatly facilitate marker-assisted breeding, the introgression of missing genes into modern varieties, and precise genome editing of selected genes and pathways. Strawberries are poised to benefit from these recent advances, providing consumers with fruit that is tastier, longer-lasting, healthier, and more beautiful.

Auxin biosynthesis gene FveYUC4 is critical for leaf and flower morphogenesis in woodland strawberry.

Auxin plays an essential role in plant growth and development, particularly in fruit development. The YUCCA (YUC) genes encode flavin monooxygenases that catalyze a rate-limiting step in auxin biosynthesis. Mutations that disrupt YUC gene function provide useful tools for dissecting general and specific functions of auxin during plant development. In woodland strawberry (Fragaria vesca), two ethyl methanesulfonate mutants, Y422 and Y1011, have been identified that exhibit severe defects in leaves and flowers. In particular, the width of the leaf blade is greatly reduced, and each leaflet in the mutants has fewer and deeper serrations. In addition, the number and shape of the floral organs are altered, resulting in smaller fruits. Mapping by sequencing revealed that both mutations reside in the FveYUC4 gene, and were therefore renamed as yuc4-1 and yuc4-2. Consistent with a role for FveYUC4 in auxin synthesis, free auxin and its metabolites are significantly reduced in the yuc4 leaves and flowers. This role of FveYUC4 in leaf and flower development is supported by its high and specific expression in young leaves and flower buds using GUS reporters. Furthermore, germline transformation of pYUC4::YUC4, which resulted in elevated expression of FveYUC4 in yuc4 mutants, not only rescued the leaf and flower defects but also produced parthenocarpic fruits. Taken together, our data demonstrate that FveYUC4 is essential for leaf and flower morphogenesis in woodland strawberry by providing auxin hormone at the proper time and in the right tissues.

The trithorax group factor ULTRAPETALA1 controls flower and leaf development in woodland strawberry.

The trithorax group (TrxG) factors play a critical role in the regulation of gene transcription by modulating histone methylation. However, the biological functions of the TrxG components are poorly characterized in different plant species. In this work, we identified three allelic ethyl methane-sulfonate-induced mutants P7, R67 and M3 in the woodland strawberry Fragaria vesca. These mutants show an increased number of floral organs, a lower pollination rate, raised achenes on the surface of the receptacle and increased leaf complexity. The causative gene is FvH4_6g44900, which contains severe mutations leading to premature stop codons or alternative splicing in each mutant. This gene encodes a protein with high similarity to ULTRAPETALA1, a component of the TrxG complex, and is therefore named as FveULT1. Yeast-two-hybrid and split-luciferase assays revealed that FveULT1 can physically interact with the TrxG factor FveATX1 and the PcG repressive complex 2 (PRC2) accessory protein FveEMF1. Transcriptome analysis revealed that several MADS-box genes, FveLFY and FveUFO were significantly up-regulated in fveult1 flower buds. The leaf development genes FveKNOXs, FveLFYa and SIMPLE LEAF1 were strongly induced in fveult1 leaves, and their promoter regions showed increased H3K4me3 levels and decreased H3K27me3 levels in fveult1 compared to WT. Taken together, our results demonstrate that FveULT1 is important for flower, fruit and leaf development and highlight the potential regulatory functions of histone methylation in strawberry.

Roles and evolution of four LEAFY homologs in floral patterning and leaf development in woodland strawberry.

The plant-specific transcription factor LEAFY (LFY), generally maintained as a single-copy gene in most angiosperm species, plays critical roles in flower development. The woodland strawberry (Fragaria vesca) possesses four LFY homologs in the genome; however, their respective functions and evolution remain unknown. Here, we identified and validated that mutations in one of the four LFY homologs, FveLFYa, cause homeotic conversion of floral organs and reiterative outgrowth of ectopic flowers. In contrast to FveLFYa, FveLFYb/c/d appear dispensable under normal growth conditions, as fvelfyc mutants are indistinguishable from wild type and FveLFYb and FveLFYd are barely expressed. Transgenic analysis and yeast one-hybrid assay showed that FveLFYa and FveLFYb, but not FveLFYc and FveLFYd, are functionally conserved with AtLFY in Arabidopsis (Arabidopsis thaliana). Unexpectedly, LFY-binding site prediction and yeast one-hybrid assay revealed that the transcriptional links between LFY and the APETALA1 (AP1) promoter/the large AGAMOUS (AG) intron are missing in F. vesca, which is due to the loss of LFY-binding sites. The data indicate that mutations in cis-regulatory elements could contribute to LFY evolution. Moreover, we showed that FveLFYa is involved in leaf development, as approximately 30% of mature leaves have smaller or fewer leaflets in fvelfya. Phylogenetic analysis indicated that LFY homologs in Fragaria species may arise from recent duplication events in their common ancestor and are undergoing convergent gene loss. Together, these results provide insight into the role of LFY in flower and leaf development in strawberry and have important implications for the evolution of LFY.

A GT-1 and PKc domain-containing transcription regulator SIMPLE LEAF1 controls compound leaf development in woodland strawberry.

Leaves are strikingly diverse in terms of shapes and complexity. The wild and cultivated strawberry plants mostly develop trifoliate compound leaves, yet the underlying genetic basis remains unclear in this important fruit crop in Rosaceae. Here, we identified two EMS mutants designated simple leaf1 (sl1-1 and sl1-2) and one natural simple-leafed mutant monophylla in Fragaria vesca. Their causative mutations all reside in SL1 (FvH4_7g28640) causing premature stop codon at different positions in sl1-1 and sl1-2 and an eight-nucleotide insertion (GTTCATCA) in monophylla. SL1 encodes a transcription regulator with the conserved DNA-binding domain GT-1 and the catalytic domain of protein kinases PKc. Expression of SL1pro::SL1 in sl1-1 completely restored compound leaf formation. The 35S::SL1 lines developed palmate-like leaves with four or five leaflets at a low penetrance. However, overexpressing the truncated SL1ΔPK caused no phenotypes, probably due to the disruption of homodimerization. SL1 is preferentially expressed at the tips of leaflets and serrations. Moreover, SL1 is closely associated with the auxin pathway and works synergistically with FveLFYa in leaf morphogenesis. Overall, our work uncovered a new type of transcription regulator that promotes compound leaf formation in the woodland strawberry and shed new lights on the diversity of leaf complexity control.

Factor of DNA methylation 1 affects woodland strawberry plant stature and organ size via DNA methylation.

RNA-directed DNA methylation (RdDM) is an epigenetic process that directs silencing to specific genomic regions and loci. The biological functions of RdDM are not well studied in horticultural plants. Here, we isolated the ethyl methane-sulfonate-induced mutant reduced organ size (ros) producing small leaves, flowers, and fruits in woodland strawberry (Fragaria vesca) due to reduced cell numbers compared with that in the wild-type (WT). The candidate mutation causes a premature stop codon in FvH4_6g28780, which shares high similarity to Arabidopsis (Arabidopsis thaliana) Factor of DNA Methylation1 (FDM1) encoding an RdDM pathway component and was named FveFDM1. Consistently, the fvefdm1CR mutants generated by CRISPR/Cas9 also produced smaller organs. Overexpressing FveFDM1 in an Arabidopsis fdm1-1 fdm2-1 double mutant restored DNA methylation at the RdDM target loci. FveFDM1 acts in a protein complex with its homolog Involved in De Novo 2 (FveIDN2). Furthermore, whole-genome bisulfite sequencing revealed that DNA methylation, especially in the CHH context, was remarkably reduced throughout the genome in fvefdm1. Common and specific differentially expressed genes were identified in different tissues of fvefdm1 compared to in WT tissues. DNA methylation and expression levels of several gibberellic acid (GA) biosynthesis and cell cycle genes were validated. Moreover, the contents of GA and auxin were substantially reduced in the young leaves of fvefdm1 compared to in the WT. However, exogenous application of GA and auxin could not recover the organ size of fvefdm1. In addition, expression levels of FveFDM1, FveIDN2, Nuclear RNA Polymerase D1 (FveNRPD1), Domains Rearranged Methylase 2 (FveDRM2), and cell cycle genes were greatly induced by GA treatment. Overall, our work demonstrated the critical roles of FveFDM1 in plant growth and development via RdDM-mediated DNA methylation in horticultural crops.

Anomalous Network Traffic Detection Method Based on an Elevated Harris Hawks Optimization Method and Gated Recurrent Unit Classifier.

In recent years, network traffic contains a lot of feature information. If there are too many redundant features, the computational cost of the algorithm will be greatly increased. This paper proposes an anomalous network traffic detection method based on Elevated Harris Hawks optimization. This method is easier to identify redundant features in anomalous network traffic, reduces computational overhead, and improves the performance of anomalous traffic detection methods. By enhancing the random jump distance function, escape energy function, and designing a unique fitness function, there is a unique anomalous traffic detection method built using the algorithm and the neural network for anomalous traffic detection. This method is tested on three public network traffic datasets, namely the UNSW-NB15, NSL-KDD, and CICIDS2018. The experimental results show that the proposed method does not only significantly reduce the number of features in the dataset and computational overhead, but also gives better indicators for every test.

Variation burst during dedifferentiation and increased CHH-type DNA methylation after 30 years of in vitro culture of sweet orange.

Somaclonal variation arising from tissue culture may provide a valuable resource for the selection of new germplasm, but may not preserve true-to-type characteristics, which is a major concern for germplasm conservation or genome editing. The genomic changes associated with dedifferentiation and somaclonal variation during long-term in vitro culture are largely unknown. Sweet orange was one of the earliest plant species to be cultured in vitro and induced via somatic embryogenesis. We compared four sweet orange callus lines after 30 years of constant tissue culture with newly induced calli by comprehensively determining the single-nucleotide polymorphisms, copy number variations, transposable element insertions, methylomic and transcriptomic changes. We identified a burst of variation during early dedifferentiation, including a retrotransposon outbreak, followed by a variation purge during long-term in vitro culture. Notably, CHH methylation showed a dynamic pattern, initially disappearing during dedifferentiation and then more than recovering after 30 years of in vitro culture. We also analyzed the effects of somaclonal variation on transcriptional reprogramming, and indicated subgenome dominance was evident in the tetraploid callus. We identified a retrotransposon insertion and DNA modification alternations in the potential regeneration-related gene CLAVATA3/EMBRYO SURROUNDING REGION-RELATED 16. This study provides the foundation to harness in vitro variation and offers a deeper understanding of the variation introduced by tissue culture during germplasm conservation, somatic embryogenesis, gene editing, and breeding programs.

The MADS-box gene FveSEP3 plays essential roles in flower organogenesis and fruit development in woodland strawberry.

Flower and fruit development are two key steps for plant reproduction. The ABCE model for flower development has been well established in model plant species; however, the functions of ABCE genes in fruit crops are less understood. In this work, we identified an EMS mutant named R27 in woodland strawberry (Fragaria vesca), showing the conversion of petals, stamens, and carpels to sepaloid organs in a semidominant inheritance fashion. Mapping by sequencing revealed that the class E gene homolog FveSEP3 (FvH4_4g23530) possessed the causative mutation in R27 due to a G to E amino acid change in the conserved MADS domain. Additional fvesep3CR mutants generated by CRISPR/Cas9 displayed similar phenotypes to fvesep3-R27. Overexpressing wild-type or mutated FveSEP3 in Arabidopsis suggested that the mutation in R27 might cause a dominant-negative effect. Further analyses indicated that FveSEP3 physically interacted with each of the ABCE proteins in strawberry. Moreover, both R27 and fvesep3CR mutants exhibited parthenocarpic fruit growth and delayed fruit ripening. Transcriptome analysis revealed that both common and specific differentially expressed genes were identified in young fruit at 6-7 days post anthesis (DPA) of fvesep3 and pollinated wild type when compared to unpollinated wild type, especially those in the auxin pathway, a key hormone regulating fruit set in strawberry. Together, we provided compelling evidence that FveSEP3 plays predominant E functions compared to other E gene homologs in flower development and that FveSEP3 represses fruit growth in the absence of pollination and promotes fruit ripening in strawberry.

Research on Network Security Situation Awareness Based on the LSTM-DT Model.

To better understand the behavior of attackers and describe the network state, we construct an LSTM-DT model for network security situation awareness, which provides risk assessment indicators and quantitative methods. This paper introduces the concept of attack probability, making prediction results more consistent with the actual network situation. The model is focused on the problem of the time sequence of network security situation assessment by using the decision tree algorithm (DT) and long short-term memory(LSTM) network. The biggest innovation of this paper is to change the description of the network situation in the original dataset. The original label only has attack and normal. We put forward a new idea which regards attack as a possibility, obtaining the probability of each attack, and describing the network situation by combining the occurrence probability and attack impact. Firstly, we determine the network risk assessment indicators through the dataset feature distribution, and we give the network risk assessment index a corresponding weight based on the analytic hierarchy process (AHP). Then, the stack sparse auto-encoder (SSAE) is used to learn the characteristics of the original dataset. The attack probability can be predicted by the processed dataset by using the LSTM network. At the same time, the DT algorithm is applied to identify attack types. Finally, we draw the corresponding curve according to the network security situation value at each time. Experiments show that the accuracy of the network situation awareness method proposed in this paper can reach 95%, and the accuracy of attack recognition can reach 87%. Compared with the former research results, the effect is better in describing complex network environment problems.

GRAS transcription factor LOSS OF AXILLARY MERISTEMS is essential for stamen and runner formation in wild strawberry.

Axillary bud development is a major factor that impacts plant architecture. A runner is an elongated shoot that develops from axillary bud and is frequently used for clonal propagation of strawberry. However, the genetic control underlying runner production is largely unknown. Here, we identified and characterized loss of axillary meristems (lam), an ethyl methanesulfonate-induced mutant of the diploid woodland strawberry (Fragaria vesca) that lacked stamens in flowers and had reduced numbers of branch crowns and runners. The reduced branch crown and runner phenotypes were caused by a failure of axillary meristem initiation. The causative mutation of lam was located in FvH4_3g41310, which encodes a GRAS transcription factor, and was validated by a complementation test. lamCR mutants generated by CRISPR/Cas9 produced flowers without stamens and had fewer runners than the wild-type. LAM was broadly expressed in meristematic tissues. Gibberellic acid (GA) application induced runner outgrowth from the remaining buds in lam, but failed to do so at the empty axils of lam. In contrast, treatment with the GA biosynthesis inhibitor paclobutrazol converted the runners into branch crowns. Moreover, genetic studies indicated that lam is epistatic to suppressor of runnerless (srl), a mutant of FveRGA1 in the GA pathway, during runner formation. Our results demonstrate that LAM is required for stamen and runner formation and acts sequentially with GA from bud initiation to runner outgrowth, providing insights into the molecular regulation of these economically important organs in strawberry.

The Roles of Floral Organ Genes in Regulating Rosaceae Fruit Development.

The function of floral organ identity genes, APETALA1/2/3, PISTILLATA, AGAMOUS, and SEPALLATA1/2/3, in flower development is highly conserved across angiosperms. Emerging evidence shows that these genes also play important roles in the development of the fruit that originates from floral organs following pollination and fertilization. However, their roles in fruit development may vary significantly between species depending on the floral organ types contributing to the fruit tissues. Fruits of the Rosaceae family develop from different floral organ types depending on the species, for example, peach fruit flesh develops from carpellary tissues, whereas apple and strawberry fruit flesh develop from extra-carpellary tissues, the hypanthium and receptacle, respectively. In this review, we summarize recent advances in understanding floral organ gene function in Rosaceae fruit development and analyze the similarities and diversities within this family as well as between Rosaceae and the model plant species Arabidopsis and tomato. We conclude by suggesting future research opportunities using genomics resources to rapidly dissect gene function in this family of perennial plants.

Roles of the Brassica napus DELLA Protein BnaA6.RGA, in Modulating Drought Tolerance by Interacting With the ABA Signaling Component BnaA10.ABF2.

Drought is a major threat to plant growth and crop productivity. Reduced level of the gibberellin would result in increased drought tolerance, but the underlying mechanism is still unclear. In Brassica napus, there are four BnaRGA genes that code for DELLA proteins, negative regulators of GA signaling. Among them, expression of BnaA6.RGA was greatly induced by drought and abscisic acid (ABA). Previously, we created the gain-of-function mutant of BnaA6.RGA, bnaa6.rga-D, and the loss-of-function quadruple mutant, bnarga by CRISPR/Cas9, respectively. Here we show that bnaa6.rga-D displayed enhanced drought tolerance, and its stomatal closure was hypersensitive to ABA treatment. By contrast, bnarga displayed reduced drought tolerance and was less sensitive to ABA treatment, but there is no difference in drought tolerance between single BnaRGA mutant and WT, suggesting a functional redundancy between the BnaRGA genes in this process. Furthermore, we found that BnaRGAs were able to interact physically with BnaA10.ABF2, an essential transcription factor in ABA signaling. The BnaA10.ABF2-BnaA6.RGA protein complex greatly increased the expression level of the drought responsive gene BnaC9.RAB18. Taken together, this work highlighted the fundamental roles of DELLA proteins in drought tolerance in B. napus, and provide desirable germplasm for further breeding of drought tolerance in rapeseed.

Genetic modulation of RAP alters fruit coloration in both wild and cultivated strawberry.

Fruit colour affects consumer preference and is an important trait for breeding in strawberry. Previously, we isolated the Reduced Anthocyanins in Petioles (RAP) gene encoding a glutathione S-transferase (GST) that binds anthocyanins to facilitate their transport from cytosol to vacuole in the diploid strawberry Fragaria vesca. The parent of rap was the F. vesca variety 'Yellow Wonder' that develops white fruit due to a natural mutation in the FveMYB10 gene. Here, we investigated the application potential of RAP in modulating fruit colours by overexpression of RAP in F. vesca and knockout of RAP in the cultivated strawberry Fragaria × ananassa. Unexpectedly, the RAP overexpression in Yellow Wonder background caused formation of red fruit. In addition, the red coloration occurs precociously at floral stage 10 and continues in the receptacle during early fruit development. Transcriptome analysis revealed that the anthocyanin biosynthesis genes were not up-regulated in RAP-ox; rap myb10 flowers at anthesis and largely inhibited at the turning stage in fruit, suggesting a coloration mechanism independent of FveMYB10. Moreover, we used CRISPR/Cas9 to knockout RAP in cultivated strawberry which is octoploid. Six copies of RAP were simultaneously knocked out in the T0 generation leading to the green stem and white-fruited phenotype. Several T1 progeny have segregated away the CRISPR/Cas9 transgene but maintain the green stem trait. Our results indicate that enhancing the anthocyanin transport could redirect the metabolic flux from proanthocyanidin to anthocyanin production at early developmental stages of fruit and that RAP is one promising candidate gene in fruit colour breeding of strawberry.

Erratum: Author Correction: Updated annotation of the wild strawberry Fragaria vesca V4 genome.

[This corrects the article DOI: 10.1038/s41438-019-0142-6.].

Retrotransposon promoter of Ruby1 controls both light- and cold-induced accumulation of anthocyanins in blood orange.

Blood orange is generally recognized to accumulate anthocyanins in its fruit pulp in a cold-inducible manner. We observed that the fruit peel of blood orange can also accumulate anthocyanins under ample light conditions. Interestingly, purple pummelo can accumulate anthocyanins only in its fruit peel but not in its pulp. The mechanism underlying the tissue specificity of anthocyanin accumulation in citrus is unknown. Here, we show that the active promoter of Ruby1, a key activator of anthocyanin biosynthesis, is also light inducible in addition to its already known cold inducibility in blood orange. Electrophoretic mobility shift assays and transient expression assays showed that HY5 positively regulated the transcription of Ruby1 by binding to the G-box motif (CACGTC). The tissue specificity of anthocyanin accumulation in the peel of purple pummelo may be due to the lack of a low temperature responsive element and a MYC binding site, which were shown to be involved in cold inducibility of CsRuby1 in blood orange by insertion of a long terminal repeat type retrotransposon in the promoter. These results bring new insights into the regulatory mechanism of anthocyanin biosynthesis in response to environmental stimuli and provide cis-elements for genetic improvement of anthocyanin-stable fruits rich in antioxidant metabolites.