- Research article
- Open Access
Construction of a high-density genetic map and identification of loci controlling purple sepal trait of flower head in Brassica oleracea L. italica
BMC Plant Biologyvolume 19, Article number: 228 (2019)
Some broccoli (Brassica oleracea L. italic) accessions have purple sepals and cold weather would deepen the purple color, while the sepals of other broccoli lines are always green even in cold winter. The related locus or gene is still unknown. In this study, a high-density genetic map was constructed based on specific locus amplified fragment (SLAF) sequencing in a doubled-haploid segregation population with 127 individuals. And mapping of the purple sepal trait in flower heads based on phenotypic data collected during three seasons was performed.
A genetic map was constructed, which contained 6694 SLAF markers with an average sequencing depth of 81.37-fold in the maternal line, 84-fold in the paternal line, and 15.76-fold in each individual population studied. In all of the annual data recorded, three quantitative trait loci (QTLs) were identified that were all distributed within the linkage group (LG) 1. Among them, a major locus, qPH.C01–2, located at 36.393 cM LG1, was consistently detected in all analysis. Besides this locus, another two minor loci, qPH.C01–4 and qPH.C01–5, were identified near qPH.C01–2, based on the phenotypic data from spring of 2018.
The purple sepal trait could be controlled by a major single locus and two minor loci. The genetic map and location of the purple sepal trait of flower heads provide an important foundation for mapping other compound traits and the identification of the genes related to purple sepal trait in broccoli.
Flavonoid, carotenoid and betalain together with chlorophyll in plants endow them with all kinds of colors [1,2,3,4,5]. Purple color is the pigment display of betalain or anthocyanin (one kind of Flavonoid). Betalain is a tyrosine-derived red-purple and yellow pigments which exist exclusively in Caryophyllales [2, 6], while anthocyanin is an abundant pigment in many different plant species and they can change color from red to purple and blue . In Brassica plants, purple color which was caused by anthocyanin accumulation was associated with the induction of its biosynthesis genes or transcription factors [7,8,9,10,11]. The mutation in the upstream regulatory region of R2R3 MYB transcription factor led to increased expression levels of the transcription factor gene, which then up regulated the expression of the structural genes involved in anthocyanin biosynthesis and endowed cauliflower with vivid purple color . The promoter substitution or deletion of BoMYBL2–1 resulted in a purple coloration of cabbage (Brassica oleracea var. capitata F. rubra) . Purple leaf genes in ornamental kale and purple stem genes in Chinese kale have been fine mapped [13, 14]. However, it is not reported that a gene or locus controlling purple sepals which are affected by low temperature in broccoli (Brassica oleracea var. italica, 2n = 2x = 18).
Some broccoli accessions have purple sepals, and the purple color intensifies during cold weather. Other accessions have green sepals even in cold winter. The purple sepals bring out dull color to flower heads of broccoli. So the broccoli flower heads with green sepals are more welcome and have higher price than those with purple sepals in Chinese market. And breeders tend to breed the broccoli new cultivars which have flower heads with ever-green sepals. It is important to discover the genetic rule, the locus and the gene of the purple sepals for broccoli breeding. We tried to map the locus controlling purple sepals in broccoli by a previous high-dense genetic map which was constructed by referring to TO1000 whole genomic sequences . But it failed to detect a locus. In this study, we mapped the loci controlling the purple sepals by reconstructing a no-reference genetic map based on the previous sequence data submitted to the National Center of Biotechnology Information (NCBI) (the BioProject ID: PRJNA449775) .
SLAF markers were developed based on alignment and clustering of all clean sequence reads using BLAT (BLAST-Like Alignment Tool) software. A total of 182,813 SLAFs were obtained with 61.20-fold average sequencing depth in the maternal line, 65.98-fold in the paternal line, and 15.82-fold in each individual offspring (Table 1). Based on the allele numbers and the sequence differences, SLAF markers could be classified into three types: polymorphic, non-polymorphic and repetitive (Table 2). Of these, 20.77% (37,969 SLAF markers) of the markers were polymorphic. As the segregation group contained DH (doubled haploid) lines, the genotypes that were not aa × bb genotype, and those that lacked parental information were abandoned. This provided us with a high-quality collection of SLAF markers. Together 18,295 SLAF markers out of 37,969 belonged to the aa×bb genotype. Polymorphic SLAF markers with lower-quality or severe partial separation or more than 3 SNPs (single nucleotide polymorphism) or covering less 70% separate individuals were filtered. A final set with 6694 markers was used to construct the genetic map with an average sequencing depth of 81.37-fold in the maternal line, 84-fold in the paternal line, and 15.76-fold in each individual population (Table 1).
Construction of the genetic map
A genetic map with 9 linkage groups (LGs) was obtained from 6694 SLAFs containing 12,980 SNPs, with a total genetic distance 880.78 cM, an average genetic distance of 0.16 cM, a maximum gap of 8.41 cM on the LG8. We identified 1358 markers that indicate segregation distortion, with 0.04% of singletons and 23.41% of miss (Fig. 1; Table 3). The individual integrity of the markers in the map was 97.88% (Additional file 1).
Collinearity analysis between the genetic map and the reference genome
The sequences of SLAF markers on the map and the sequences of reference genome TO1000 were aligned. Approximately, 95.35% (6383) of the markers from a total of 6694 SLAF markers were mapped to corresponding positions in the reference genome. Only 311 SLAF markers (4.65% of 6694 SLAFs) present in the map could not be assigned to positions in the reference genome. Therefore, the markers in the map show high collinearity to the reference genome (Fig. 2).
Purple sepal trait of flower head and its inheritance model
The purple sepal trait of flower heads in the DH population and their parents were surveyed in the autumn of 2015 and 2017 and in the spring of 2018. The sepal color of the maternal line is purple, and low temperature intensifies the purple coloration; while the sepal color of the paternal line never turned purple, even at sub-zero temperatures (Fig. 3). The sepal color of the F1 hybrid was purple. Of the 309 genotypes in the F2 population, the sepals of 236 individuals were purple, and sepals of 73 were not purple. A χ2 test indicated that the fitness of the segregation to the expected ratio of 3:1. From the results, we can infer that the purple sepal trait is likely controlled by a single locus.
Location of the purple sepal allele
Based on the genetic map and using the phenotype data collected during the autumn of 2015 and 2017, a single locus on LG1 was identified. Both the loci were qPH.C01–1 (30.093–39.543 cM) and qPH.C01–2 (36.393–36.393 cM), respectively. The locus qPH.C01–1 contains qPH.C01–2. That’s, the same single locus region was found based on the data of 2015 and 2017. While three loci, qPH.C01–4 (24.58–24.58 cM), qPH.C01–5 (32.455–32.455 cM), and qPH.C01–3 (36.393–40.331 cM) were found in LG1, based on the data from spring of 2018 (Table 4). These three loci were close to each other in LG1. Apparently, the locus at 36.393 cM on LG1, which accounts for 10.3% of the phenotypic variation, was common in all the three analyses performed.
Broccoli is one varietas in B. oleracea . The genomes of other two varietas in this species, B. oleracea var. capitata line 02–12 and kale-like TO1000 DH, have been sequenced and successfully assembled [17, 18]. They are very important for genetic research in B. oleracea crops, and their genomes have been referenced in the construction of a genetic map in cauliflower, and for fine mapping in ornamental kale and Chinese kale [13, 14, 19,20,21]. However, there are assembly errors and mis-anchored sequence scaffolds in the cabbage genome ‘02–12’ and kale-like genome TO1000 [14, 22]. Actually, such gaps and misassemblies are common in whole-genome assembly of many species [13, 22,23,24]. Biological complexities, especially the Brassicaceae-specific triplication events and continuous amplification of transposable elements in B. oleracea, have complicated the assembly of the genomes, resulting in gaps and misassemblies [25,26,27]. These errors lead to incorrect mapping, failure to accurately identify and clone gene(s) based on a related genetic map, inability to identify a related locus, and inefficient marker-assisted selection [13, 14]. Moreover, one or even several genotypes sequenced at a high resolution may not contain all the sequences for that particular species. For example, there were huge structural variations including inversions, translocations, and presence/absence variations between the two genomes of elite ‘indica’ rice varieties’ ZS97RS1 and MH63RS1. Besides, there were some genes present in variety ZS97RS1 that were absent in MH63RS1, and vice-versa . In our previous study, no locus related to the purple sepal trait was found in that genetic mapping exercise, which was constructed based on the reference genome TO1000. Broccoli and the kale-like line TO1000, between which the main phenotypic difference is that the former has a big flower head, are two different varietas of B. oleracea. It is possible that sequences differences and structural variations exist between the two varietas, making the identification of the purple sepal allele in the genetic map difficult, particularly, if based on the TO1000 reference genome. So, it was necessary to reconstruct a non-reference genetic map for mapping the purple sepal loci.
One locus related to the purple sepal trait on the map, which was located at 36.393 cM in LG1 was identified based on the phenotypic data collected during two autumn seasons, while the three loci located at intervals of 15.751 cM, 24.58 cM, and 40.331 cM from LG1, were found based on the phenotypic data from the spring of 2018. The identification of three loci seems inconsistent with the results from the χ2 test, which, indicates that the purple sepal trait may be controlled by a single locus. As the purple sepal color in the maternal line intensifies during low temperature, potential regulatory genes related to the purple color could be induced by cold. In addition, the proximity of all the three loci in LG1, may indicate the presence of a regulatory network involved in controlling the purple coloration, and perhaps is active and sensitive to low temperature, yet the related trait is difficult for the human eye to visualize. In autumn, the temperature reaches 30 °C, coinciding with the beginning of growth, after which the temperature decreases. The temperature during the whole growth period stays above zero. In contrast, the temperature during spring increases gradually and reaches a high point. Assuming that the purple color is regulated by low temperature, it may be easier to observe the purple sepal color of broccoli flower heads in spring. These two loci could be related to low temperature inducing purple coloration of sepals in the flower heads of broccoli. So the purple sepal might be controlled by a single major locus and two minor loci, which would be confirmed by fine mapping, for the other both loci had low LOD values and low phenotypic variation explained by the QTLs.
Anthocyanins are catalytically synthesized by a series of enzymes encoded by corresponding genes, such as chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavanone 3′-hydroxylase (F3’H), flavonol synthase (FLS), dihydroflavonol 4-reductase (DFR), encoding glutathione-S-transferases (GST), leucoanthocyaniidin oxgenase (LDOX), anthocyanidin reductase (ANR), UDP-glucose: flavonoid 3-O-glucosyltransferase (UD3GT) [12, 29, 30]. While, anthocyanin biosynthesis also is regulated by transcription factors such as MYB, bHLH, WD40, and ZAT6 [7, 30,31,32,33]. Anthocyanins biosynthesis and regulation not only involve many structural genes and transcription factors, but also are influenced by environmental factors, such as light, temperature and other substances [12, 34,35,36]. The expression of structural genes of anthocyanins biosynthesis such as CHS, DFR, and GST were higher in low temperature conditions [37,38,39]. Light could induce the expression of the transcription factor genes (such as MYB75 and MYB90), the structural genes (such as CHS, DFR, F3H, and LDOX), the kinase gene (such as MAP KINASE4) and ethylene response factor genes which are involved in anthocyanin biosynthesis [40,41,42]. A PA1-type MYB transcription factor, MdMYBPA1 was identified from red-fleshed apple and it redirected the flavonoid biosynthetic pathway by its promoter’s low-temperature-responsive (LTR) cis-element directly binging MdbHLH33 in low temperature conditions . In purple cabbage, a MYB transcript factors, BoMYBL2–1, negatively regulated anthocyanin synthesis . Substitution or deletion of its promoter resulted in different degree of purple coloration and different sensitivity to low temperature . A R2R3 MYB transcription factor, BoPAP1, might be responsible for purple leaves and up regulated by low temperature in purple kale . These researches show that different plants have different molecular regulation mechanism in anthocyanins accumulation induced by low temperature. Anthocyanins are obviously accumulated in leaves and induced by low temperature in the cabbage and the kale [10, 12]. In this study, no obvious purple color are observed in leaves of the broccoli line DH16–2, but its sepals are purple and become deep purple in cold conditions. What’s more, the genes for the purple traits in kale, cabbage and kohlrabi were all mapped to the chromosome 6 which was the same locus as the BoMYB2 gene in cauliflower . While the loci controlling the purple sepals in broccoli was located on the LG1 in this study. It is still unknown which genes are responsible for that. Further study will be exerted for better understanding of anthyocanins biosynthesis and regulation in purple sepals of broccoli.
In this study, we constructed a high-density genetic map using a DH segregation population of a cross between DH16–2 and DH28–4. Based on this map, a major locus and two minor loci of the purple sepal using phenotypic data collected from three seasons was detected on LG1.
The maternal line DH16–2 (purple sepals), paternal line DH28–4 (green sepals), their hybrid and the segregating population with 127 DH lines was generated using the microspore culturing methods. The parental lines and segregating population individuals were planted in the Yangdu Experimental Greenhouse of the Zhejiang Academy of Agricultural Sciences in August 2014, August 2017 and January 2018. The purple sepal trait was surveyed when the flower heads were mature during the crop season. During the period of cultivation, every plant was checked to make sure no leave covering on the flower head in order to avoid the effect of light on the sepals color.
SLAF markers and genotyping
Based on the raw data from BioProject (PRJNA449775), SLAF markers were identified, and genotyping was performed as described in . Initially, all the low-quality reads were removed and high-quality raw reads were sorted to each progeny according to the duplex barcodes. Clean reads were then clustered based on similarity over 90% after the barcodes and the terminal 5 bp positions were trimmed from each high-quality read. The sequences that were clustered together were defined as one SLAF locus . SNP loci for each SLAF locus was identified between the parents. SLAF markers with more than 3 SNPs were deleted. Alleles of every SLAF locus were recognized on the basis of the reads from the parents. Diploid broccoli plants were used in this study, so, each SLAF locus could include a maximum of 4 genotypes. Only SLAF loci with two to four alleles were reserved as potential markers and those SLAF loci with more than 4 alleles were categorized as repetitive SLAF markers and rejected. All polymorphic SLAF marker loci consistent both, in the parental and offspring SNP loci were genotyped. The marker codes of the polymorphic SLAF markers were analyzed on the basis of the population type, such as DH with only one segregation type (aa × bb).
Genotype scoring was done using a Bayesian approach to guarantee consistent genotyping quality . First, a posterior conditional probability was calculated based on the coverage of each allele and the number of SNPs. Genotyping quality score based on the probability was used to select qualified markers for the succeeding analysis. Low-quality markers or individuals were discarded in the dynamic process. When the average genotype quality values of all SLAF markers reached the critical value, the process stopped. High-quality SLAF markers for constructing the genetic map were screened using the following criteria: Sequences with more than 10-fold depth of the parents, markers covering more than 70% genotypes of all offspring, and the segregation distortion as examined by the χ2 test.
Genetic map construction and QTL mapping
The HighMap software developed by Beijing Biomarker Technologies Corporation, was utilized to construct a high-density and high-quality map . Recombinant frequencies and the maximum likelihood method (MLOD) scores between markers were used to deduce the linkage phases and two-point analysis was used. Molecular markers were divided into different linkage groups based on the MLOD score, and each linkage group was regarded as a chromosome. A genetic map was preliminarily constructed, and the initial sequence of the markers was obtained based on the MLOD score. As the genotyping results contain certain errors related to the molecular markers, genotype correction was performed based on the sequence of the markers in the map. A high-quality genetic map was finally constructed using the Kosambi mapping function after several rounds of correction . MapQTL 5.0 was used to analyze the phenotypic data for interval mapping.
Comparison between the genetic map and the reference genome
In order to compare the linearity between markers on the genetic map and the reference genome, SLAF markers sequences of the map were positioned on the reference genome, kale-like line TO1000. The physical locations of the SLAF markers in the map were identified. These physical locations in the map were compared to their positions in the reference genome and a linear analysis was conducted.
Quantitative trail locus
Specific locus amplified fragment
Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54:733–49.
Polturak G, Aharoni A. "La vie en rose": biosynthesis, sources, and applications of Betalain pigments. Mol Plant. 2018;11(1):7–22.
Llorente B, Martinez-Garcia JF, Stange C, Rodriguez-Concepcion M. Illuminating colors: regulation of carotenoid biosynthesis and accumulation by light. Curr Opin Plant Biol. 2017;37:49–55.
Outchkourov NS, Karlova R, Hölscher M, Schrama X, Blilou I, Jongedijk E, et al. Transcription factor-mediated control of anthocyanin biosynthesis in vegetative tissues. J Plant Physiol. 2018;176(2):1862–78.
Han F, Cui H, Zhang B, Liu X, YangL ZM, et al. Map-based cloning and characterization of BoCCD4, a gene responsible for white/yellow petal color in B. oleracea. BMC Genomics. 2019;20:242.
Lightfoot DJ, Jarvis DE, Ramaraj T, Lee R, Jellen EN, Maughan PJ. Single-molecule sequencing and hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution. BMC Biol. 2017;15(1):74.
Chiu LW, Zhou XJ, Burke S, Wu X, Prior RL, Li L. The purple cauliflower arises from activation of a MYB transcription factor. Plant Physiol. 2010;154(3):1470–80.
Zhang Y, Chen G, Dong T, Pan Y, Zhao Z, Tian S, et al. Anthocyanin accumulation and transcriptional regulation of anthocyanin biosynthesis in purple box choy (Brassica rapa var. chinensis). J Agric Food Chem. 2014;62:12366–76.
Xie Q, Hu Z, Zhang Y, Taian S, Wang Z, Zhao Z, et al. Accumulation and molecular regulation of anthocyanin in purple timorous stem mustard (Brassica juncea var. tumida Tsen et lee). J Agric Food Chem. 2014;62:7813–21.
Zhang B, Hu Z, Zhang Y, Li Y, Zhou S, Chen G. A putative functional MYB transcription factor induced by low temperature regulates anthocyanin biosynthesis in purple kale (Brassica oleracea var. acephala f. tricolor). Plant Cell Rep. 2012;31:281–9.
Zhang Y, Hu Z, Zhu M, Zhu Z, Wang Z, Tian S, et al. Anthocyanin accumulation and molecular analysis of correlated genes in purple kohlrabi (Brassica oleracea var. gongylodes L.). J Agr Food Chem. 2015;63:4160–9.
Song H, Yi H, Lee M, Han CT, Lee J, Kim HR, et al. Purple Brassica oleracea var. capitata F. rubra is due to the loss of BoMYBL2-1 expression. BMC Plant Biol. 2018;18:82–97.
Liu X, Gao BZ, Han FQ, Fang ZY, Yang LM, Zhuang M, et al. Genetics and fine mapping of a purple leaf gene, BoPr, in ornamental kale (Brassica oleracea L. var. acephala). BMC Genomics. 2017;18:230.
Tang Q, Tian M, An G, Zhang W, Chen J, Yan C. Rapid identification of the purple stem (Ps) gene of Chinese kale (Brassica oleracea var. alboglabra) in a segregation distortion population by bulked segregant analysis and RNA sequencing. Mol Breed. 2017;37(153) https://doi.org/10.1007/s11032-017-0752-3.
Yu H, Wang J, Zhao Z, Sheng X, Shen Y, Branca F, et al. Construction of a high-density genetic map and identification of loci related to hollow stem trait in broccoli (Brassic oleracea L. italica). Front Plant Sci. 2019. https://doi.org/10.3389/fpls.2019.00045.
Cheng F, Sun R, Hou X, Zheng H, Zhang F, Zhang Y, et al. Subgenome parallel selection is associated with morphotype diversification and convergent crop domestication in Brassica rapa and Brassica oleracea. Nat Genet. 2016;48(10):1218–24.
Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IA, et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat Commun. 2014;5:3930.
Parkin IA, Koh C, Tang H, Robinson SJ, Kagale S, Clarke WE, et al. Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 2014;15:R77.
Zhao Z, Gu H, Sheng X, Yu H, Wang J, Huang R, et al. Genome-wide single-nucleotide polymorphisms discovery and high-density genetic map construction in cauliflower using specific-locus amplified fragment sequencing. Front Plant Sci. 2016;7(393):334.
Mei JQ, Wang JH, Li YH, Tian S, Wei DY, Shao CG, et al. Mapping of genetic locus for leaf trichome in Brassica oleracea. Theor Appl Genet. 2017;130(9):1953–9.
Liang JL, Ma Y, Wu J, Cheng F, Liu B, Wang XW. Map-based cloning of the dominant genic male sterile Ms-cd1 gene in cabbage (Brassica oleracea). Theor Appl Genet. 2017;130:71–9.
Lee J, Izzah NK, Choi BS, Joh HJ, Lee SC, Perumal S, et al. Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L). DNA Res. 2015;14:S113.
Steinberg KM, Schneider VA, Graves-Lindsay TA, Fulton RS, Agarwala R, Huddleston J, et al. Single haplotype assembly of the human genome from a hydatidiform mole. Genome Res. 2014;24:2066–76.
Ni XM, Xia QJ, Zhang HB, Cheng S, Li H, Fan GY, et al. Updated foxtail millet genome assembly and gene mapping of nine key agronomic traits by resequencing a RIL population. GigaScience. 2017;6:1–8.
Lysak MA, Koch MA, Pecinka A, Schubert I. Chromosome triplication found across the tribe Brassicaceae. Genome Res. 2005;15:516–25.
Wang XW, Wang HZ, Wang J, Sun RF, Wu J, Liu SY, et al. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43:1035–9.
Liu D, Ma C, Hong W, Huang L, Liu M, Liu H, et al. Construction and analysis of high-density linkage map using high-throughput sequencing data. PLoS One. 2014b;9:e98855.
Zhang JW, Chen LL, Xing F, Kudrna DA, Yao W, Copetti D, et al. Extensive sequence divergence between the reference genomes of two elite indica rice varieties Zhenshan 97 and Minghui 63. Proc Natl Acad Sci U S A. 2016;113(35):E5163–71.
Xu W, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015;20(3):176–85.
Lepiniec L, Debeaujon I, Routaboul JM, Baudry A, Pourcel L, Nesi N, et al. Genetics and biochemistry of seed flavonoids. Annu Rev Plant Biol. 2006;57:405–30.
Shi H, Liu G, Wei Y, Chan Z. The zinc-finger transcription factor ZAT6 is essential for hydrogen peroxide induction of anthocyanin synthesis in Arabidopsis. Plant Mol Biol. 2018;97:165–76.
Zhu ZX, Wang HL, Wang YT, Guan S, Wang F, Tang JY, et al. Characterization of the cis elements in the proximal promoter regions of the anthocyanin pathway genes reveals a common regulatory logic that governs pathway regulation. J Exp Bot. 2015;66:3775–89.
Xie Y, Tan H, Ma Z, Huang J. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressors of the MYB/bHLH/WD40 complex in Arabidopsis thaliana. Mol Plant. 2016;9(5):711–21.
Sakuta M. Diversity in plant red pigments: anthocyanins and betacyanins. Plant Biotech Rep. 2014;8:37–48.
Solfanelli C, Poggi A, Loreti E, Alpi A, Perata P. Sucrose-specific induction of the anthocyanin biosynthetic pathway in Arabidopsis. Plant Physiol. 2006;140(2):637–46.
Li J, He YJ, Zhou L, Liu Y, Jiang M, Ren L, et al. Transcriptome profiling of genes related to light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.) before purple color becomes evident. BMC Genomics. 2018;19:201–12.
Leyva A, Jarillo JA, Salinas J, Martinez-Zapater JM. Low temperature induces the accumulation of phenylalanine ammonia-lyase and chalcone synthase mRNAs of Arabidopsis thaliana in a light-dependent manner. Plant Physiol. 1995;108(1):39–46.
Shvarts M, Borochov A, Weiss D. Low temperature enhances petunia flower pigmentation and induces chalcone synthase gene expression. Physiol Plantarum. 1997;99:67–72.
Carmona L, Alquézar B, Marques VV, Peña L. Anthocyanin biosynthesis and accumulation in blood oranges during postharvest storage at different low temperatures. Food Chem. 2017;237:7–14.
Cominelli E, Gusmaroli G, Allegra D, Galbiati M, Wade HK, Jenkins GI, et al. Expression analysis of anthocyanin regulatory genes in response to different light qualities in Arabidopsis thaliana. J Plant Physiol. 2008;165(8):886–94.
Li S, Wang W, Gao J, Yin K, Wang R, Wang CC, et al. Myb75 phosphorylation by mpk4 is required for light-induced anthocyanin accumulation in Arabidopsis. Plant Cell. 2016;28(11):2866–83.
Ni J, Bai S, Zhao Y, Qian M, Tao R, Yin L, et al. Ethylene response factors Pp4ERF24 and Pp12ERF96 regulate blue light-induced anthocyanin biosynthesis in ‘red Zaosu’ pear fruits by interacting with MYB114. Plant Mol Biol. 2019;99:67–78.
Wang N, Qu C, Jiang S, Chen Z, Xu H, Fang H, et al. The proanthocyanidin-specific transcription factor MdMYBPA1 initiates anthocyanin synthesis under low-temperature conditions in red-fleshed apples. Plant J. 2018;96:39–55.
Yan C, An G, Zhu T, Zhang W, Zhang L, Peng L, et al. Independent activation of the BoMYB2 gene leading to purple traits in Brassica oleracea. Theor Appl Genet. 2019;132:895–906.
Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One. 2013;8(3):e58700.
Roman LT, Michael YG, Darren AN, Eugene VK. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucl Acids Res. 2000;28(1):33–6.
Zhang J, Zhang Q, Cheng T, Yang W, Pan H, Zhong J, et al. High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. et Zucc). DNA Res. 2015;22(3):183–91.
Huang XH, Zhao Y, Wei XH, Li CY, Wang AH, Zhao Q, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44:32–9.
The authors thank Editage, a company dedicated to helping international researchers publish their findings in the best English language, for editing the language.
This work was supported by the Ministry of Science and Technology of the People’s Republic of China for key research and development plan (2017YFD0101805) and key projects of international scientific and technological innovation cooperation between governments (2017YFE0114500); Science and Technology Department of Zhejiang Province for Key Agriculture Development Project (2016C02051–5) for Public Welfare Project (2017C32076); the Science and Technology Department of Ningbo for the project (2015C110008); Zhejiang Academy of Agricultural Sciences for subject construction (A1) and the European Union’s Horizon 2020 Research and Innovation Program under grant agreement No. 774244. The funders had no role in study design, data collection, data analysis, data interpretation, the writing of the manuscript or decision to publish.
Availability of data and materials
All the raw data have been submitted to the National Center of Biotechnology Information (NCBI) (the BioProject ID: PRJNA449775).
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.
Individual integrity of specific locus amplified fragment (SLAF) markers in the genetic map. (PNG 15 kb)