Next Article in Journal
Tryptophan Metabolism and Gut-Brain Homeostasis
Next Article in Special Issue
MYB-Mediated Regulation of Anthocyanin Biosynthesis
Previous Article in Journal
Exendin-4 Promotes Schwann Cell Survival/Migration and Myelination In Vitro
Previous Article in Special Issue
Overexpression of Salicylic Acid Carboxyl Methyltransferase (CsSAMT1) Enhances Tolerance to Huanglongbing Disease in Wanjincheng Orange (Citrus sinensis (L.) Osbeck)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 6
10.3390/ijms22062972
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Mapping and Identification of the Candidate Gene for White Seed Coat in Cucurbita maxima

by
Yuzi Shi
,
Meng Zhang
,
Qin Shu
,
Wei Ma
,
Tingzhen Sun
,
Chenggang Xiang
,
Changlin Wang
* and
Ying Duan
*
Key Laboratory of Biology and Genetic Improvement of Horticultural Crops of Ministry of Agriculture and Rural Affairs, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(6), 2972; https://doi.org/10.3390/ijms22062972
Submission received: 7 February 2021 / Revised: 5 March 2021 / Accepted: 8 March 2021 / Published: 15 March 2021
(This article belongs to the Collection Genetics and Molecular Breeding in Plants)
Browse Figures
Versions Notes

Abstract

:
Seed coat color is an important agronomic trait of edible seed pumpkin in Cucurbita maxima. In this study, the development pattern of seed coat was detected in yellow and white seed coat accessions Wuminglv and Agol. Genetic analysis suggested that a single recessive gene white seed coat (wsc) is involved in seed coat color regulation in Cucurbita maxima. An F2 segregating population including 2798 plants was used for fine mapping and a candidate region containing nine genes was identified. Analysis of 54 inbred accessions revealed four main Insertion/Deletion sites in the promoter of CmaCh15G005270 encoding an MYB transcription factor were co-segregated with the phenotype of seed coat color. RNA-seq analysis and qRT-PCR revealed that some genes involved in phenylpropanoid/flavonoid metabolism pathway displayed remarkable distinction in Wuminglv and Agol during the seed coat development. The flanking InDel marker S1548 was developed to predict the seed coat color in the MAS breeding with an accuracy of 100%. The results may provide valuable information for further studies in seed coat color formation and structure development in Cucurbitaceae crops and help the molecular breeding of Cucurbita maxima.

1. Introduction

Cucurbita maxima is a worldwide economic crop in Cucurbita species used mainly for the consumption of the flesh and edible seed [1,2]. The kernel of edible seed pumpkin in Cucurbita maxima is rich in unsaturated fatty acids, amino acid, natural antioxidants, and other nutrients [3,4]. The predominant fatty acids in pumpkin seed are linoleic, oleic, palmitic, and stearic acid with high oxidative stability that would be suitable for food and industrial applications [5,6].
Seed coat color is an important agronomic trait in seed oil crops such as Brassica napus, Glycine max, and Medicago truncatula. The seed coat color and thickness are related to the efficiency of seed processing [7]. Compare with dark-colored seed coat, light-colored seed is a desirable trait with great potential for improving seed quality [8,9]. In the edible seed genetic breeding, the varieties displaying white and thin seed coat are preferred for removing the seed coat conveniently and display a high kernel yield [10].
The seed coat displayed a wide range of colors such as brown, yellow, ivory, and white in Cucurbita species [10,11]. Seed coat color has been extensively studied in Arabidopsis and several oil crops such as Brassica napus, Glycine max, and some species in Cucurbitaceae [12]. In Arabidopsis, the seed coat arises from the integument of the maternal ovule and mainly develops into five layers. The integument layers include endothelium and two adjacent layers. The outer integument layers include subepidermal layers and epidermal cell layers [13,14]. A secondary thickened cell wall is formed in subepidermal cell layers and is deposited around the epidermis layers. Seed coat structure observed in Cucurbita species are referred as chlorenchyma, parenchyma, sclerenchyma, hypodermis, and epidermis [15,16]. The normal seed coat hypodermis, sclerenchyma, and parenchyma tissues are lignified gradually in the development of seed coat [15]. After ripening, the outer side of seed coat will dehydrate and form a high lignified shell to protect the seeds from stress of external environmental factors.
Flavonoids are secondary metabolites in phenylpropanoid pathway. Flavonols, anthocyanins, and proanthocyanidins are the three main flavonoids [17]. Proanthocyanidins are also known as concentrated tannic acid, and the seed coat is brown when the proanthocyanidins are oxidized [13]. Flavonoid metabolisms (predominantly tannins) are produced in endothelium and deposited to form brown seed coat after oxidization in mature seeds [12,14,18,19]. By analyzing the components of different color seed coats in Arabidopsis, Brassica rapus, and other plants, the key enzyme regulation steps from phenylpropanes to proanthocyanidins have been determined clearly [20]. The initial substrates for the synthesis of flavonoids are derived from coumaroyl-CoA of phenylpropane metabolic pathways. Serial components in seed coat pigments synthesis have been identified in flavonoid metabolism pathway including Chalcone synthase (CHS), Chalcone isomerase (CHI), Flavonoid 3′-hydroxylase (FLS), Dihydroflavanol-4-reductase (DFR), Anthocyanin synthase (ANS), and Leucoanthocyanins reductase (LAR) [17,21]. Arabidopsis transparent testa (tt) genes encoding proteins involved in flavonoid metabolism display important roles in seed coat color formation and seed structures [22,23,24]. tt mutants display light seed coat color, defect in seed coat permeability and lower dormancy [25]. TT genes are regulated by MYB-bHLH-WD40 (MBW) protein complex. MBW is composed of R2R3-MYB, bHLH transcription factor, and WD40-repeat family members [12,26,27]. Among them, R2R3-MYB family transcription factors play a critical role in the production of anthocyanins or proanthocyanidins [28].
In Cucurbitaceae species, compared to melon and cucumber, watermelon and pumpkin/squash display rich diversity on seed traits including seed size, seed weight, and seed colors, which seem to depend on the intended uses and under a significant artificial selection in modern breeding [10]. At present, genetic analysis and QTL localization of seed coat color has been studies in Cucurbitaceae crops such as melon and watermelon. In melon, the white seed coat is a quality trait, and it is dominant to yellow and brown [29]. In watermelon, the seed coat color was determined by four loci (r, w, t for red, white, and tan, respectively; d for producing a black and dotted seed coat) which have been identified on chromosome 3, 5, 6, 8 [7,30,31]. The light cream color was recessive to the brown or red seed coat color, and black color was dominant over brown light cream and red color of seed coat separately or independently [32]. Several QTLs associated with black and light yellow seed coat were located on chromosome 1, 3, and 8 [30]. The candidate genes for qsc-c3-1 associated with black seed coat color was narrowed to a 70.2 kb region and a gene encoding Polyphenol oxidase (PPO) was identified to associate with black seed coat [30]. However, the inheritance of seed coat color was not very clear in Cucurbita maxima.
Recently, great efforts have been made in genetic linkage map construction and genome information in Cucurbita species. Molecular markers including random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), simple sequence repeat (SSR), and single-nucleotide polymorphism (SNP) have been developed to construct genetic linkage map in Cucurbita [33,34,35]. Improvements of whole genome sequences and high-density genetic maps have promoted the precision of genomic data in Cucurbita and support a better foundation for the research and application of candidate genes in Cucurbita species [36,37,38,39,40,41]. With effective utilization of molecular markers and genetic map resources, QTLs for many agronomic traits and candidate genes have been identified and applicated accordingly in breeding selection [42,43,44,45,46]. Some QTLs of seed trait such as seed size, seed width, and seed yields have been identified in Cucurbita maxima and candidate genes are precited [10,47]. As one of the most important traits in Cucurbita maxima, studies about genetic and molecular mechanism of seed coat color have not been elucidated clearly.
In this study, two inbred accessions displaying white and yellow seed coat color were used to identify the pigment components of seed coat color formation in Cucurbita maxima. The inheritance of seed coat color was analyzed and a candidate gene encoding MYB transcription factor was predicted by map-based cloning. The molecular marker associated with the phenotype of seed coat color has been developed for MAS. The present results may highlight the understanding of the genetic basis for seed coat color formation and help further functional research in Cucurbita maxima.

2. Results

2.1. The Pigment Deposition of Seed Coat in Cucurbita maxima

Seed coat color is a very important agronomic trait of edible seed pumpkin in Cucurbita maxima. In this study, the seed coat color of 54 inbred accessions were observed and categorized into white or yellow. Twenty-three inbreds displaying white and 31 inbreds displaying yellow seed coat were obtained. Wuminglv (yellow seed coat, P1), Agol (white seed coat, P2), Lizimanhui (yellow seed coat, P3), and Diaogua (white seed coat, P4) were selected as parental lines and crossed to produced genetic populations (Figure 1A,B, Supplement Figure S1). At 0 DPA (Day Post Anthesis) of Wuminglv and Agol, the seed coat colors were both transparent. At 21~28 DPA, the seed coat color of Wuminglv was transferred into yellow, while Agol still remained white, suggesting 21~28 DPA was the important period of seed coat color formation in Cucurbita maxima (Figure 1C). L*, a*, and b* values and the comprehensive chromaticity E were obtained by chromometer to eliminate the deviation of vision subjective understanding of color. The L* color value in white seed coat ranged from 79.6 to 92.79, while that in yellow seed coat ranged from 60.39 to 72.29. The color value in BC1, BC2, and two F2 segregating populations displayed an expected arrangement according to the distribution of separation ratio (Table 1).
The color of petal and of flowers displayed yellow and the difference was not very obvious in Wuminglv and Agol (Supplement Figure S2A), followed by a slight difference in flesh color after fruit maturation (42 DPA) (Supplement Figure S2B,C).
DMACA (4-(Dimethylamino) cinnamaldehyde) is used as an indicator of procyanidins and flavan-3-ol [14,24,48]. To determine the composition of color pigment of seed coat qualitatively, DMACA staining was conducted in Wuminglv, Agol, Diaogua, Lizimanhui, and 76 inbred accessions including 29 white and 47 yellow seed coat accessions. The accessions displaying yellow seed coat were stained dark red and the accessions displaying white seed coat were stained rarely, which was closely related to the seed coat color (Figure 1B, Supplement Figure S3). Dynamic analysis of soluble and insoluble PA in different stages of seed development indicated that the content of procyanidins in the seed coat was significantly different at 21 DPA in Wuminglv and Agol, which is in accordance with the observation of seed coat color (Figure 1D,E).

2.2. The Histochemical Analysis of Seed Coat

To determine the pigment accumulation and differentiation layers of seed coat, anatomic structure was observed at 0, 7, 14, 21, 28, 35 DPA seed coats in Wuminglv and Agol (Figure 2A). At 14 DPA, there were five layers of seed coat: Epidermis (E), Hypodermis (H), Sclerenchyma (S), Parenchyma (P), and Chlorenchyma (C) (Figure 2(Ae,Af)). At 14~21 DPA, the cell differentiation in S layer displayed distinguish in Wuminglv and Agol (Figure 2(Ae–Ah)). The statistical analysis of cell width, cell length, and layers thickness in H, S, and P layers were calculated (Figure 2B). The thickness of seed coat was not obvious at the early stage of development (0~14 DPA), but the structure of H and S layers displayed significant differentiation where the pigment deposition at 21~28 DPA. There were two S-cell layers in Wuminglv but only one S-cell layer in Agol. Toluidine blue staining revealed that the polymeric phenolic compounds in Wuminglv were stained more intensively than that in Agol at 14~21 DPA, followed a dramatic difference at 28~35 DPA (Figure 2(Ai–Al)). As lignin was one of the important polyphenolic compounds in the secondary cell wall formation, the content of seed coat lignin was determined by phloroglucinol staining and thiogylcolic acid method (Supplement Figure S4A,B). The results indicated that more lignin was accumulated at the mature seed stage in Wuminglv, in accordance with a secondary cell wall thickened in S layers, compared to that in Agol.

2.3. Genetic Analysis of Seed Coat Color

Two genetic populations were used to determine the inheritance of seed coat color. For Wuminglv (P1) and Agol (P2), a six-generation population including P1, P2, F1, BC1, and BC2 was constructed and the seed coat color was determined by visual method and chromometer. As visual observation, all F1 plants exhibited the same seed coat color as Wuminglv. In 377 individuals of an F2 segregating population, 282 plants displayed yellow and 95 plants displayed white seed coat. The segregation pattern fit a Mendelian ratio of 3:1 (χ2 = 0.002, p < 0.05). In the 121 individuals of BC2 population, 65 plants displayed yellow seed coat and 56 plants displayed white seed coat. Statistical analysis showed that the segregating pattern of BC1P2 fit 1:1 separation ratio (χ2 = 0.17, p < 0.05) (Table 2). The similar results were displayed in the L*, a*, and b* determination (Supplement Table S1). For Lizimanhui (P3) and Diaogua (P4), similar inheritance was displayed in F2 segregating population (Supplement Table S2). The results indicated that the phenotype of seed coat color was controlled by a recessive single gene and was named as white seed coat (wsc).

2.4. Map-Based Cloning of White Seed Coat

The whole genome resequencing (20× coverage) of Wuminglv and Agol were performed on Illuminar Hiseq platform and 324 InDel markers were designed evenly distributed every 1 Mb on each chromosome. One hundred and twenty-five InDel markers displayed clear and stable bands and polymorphisms consistent with two parents and F1 were selected to analyze an F2 segregating population of 2798 individuals. WSC was located at Chromosome 15 in primary mapping. Using two flanking markers S1517 and S1528 at a distance of 315.4 kb, 34 recombinants were identified and the phenotype of each recombinant was confirmed by F2:3 segregating population (Figure 3A). Forty-two novel InDel markers and 11 dCAPS markers were designed in this region based on the whole genome sequencing (Supplement Table S3). Seventeen InDel markers and nine dCAPS markers displayed polymorphisms among Wuminglv and Agol. The candidate region was ultimately narrowed in 60.4 kb interval between dCAPS marker SD-4 and InDel marker S1546 (Figure 3B).

2.5. Candidate Gene and Expression Analysis

There were nine genes in the candidate region (Table 3) including vesicle-associated membrane protein, MYB transcription factor, pentatricopeptide repeat-containing protein, unknown gene, and some other genes. The function of some genes has been studied in Arabidopsis and other crops, such as CmaCh15G005220/vesicle-associated membrane protein 727 and CmaCh15G005270/AtDIV2 [49,50]. From the gene annotations, few information of genes was related to seed coat color formation or pigment synthesis.
As CmaCh15G005270, CmaCh15G005280, CmaCh15G005290, and CmaCh15G005300 were considered as the candidate gene, the full length of each gene including promoter, 5′-UTR and 3′-UTR was sequenced in Wuminglv and Agol. In CmaCh15G005270, there were two insertions, two deletions, and nine mutations presented in Wuminglv vs. Agol in a region of –1300~–800 bp from the initiation site. Meanwhile, one T insertion and one T–C transition were found in the first intron of the gene (Figure 3A).
To confirm whether these variant sequences were linked to seed coat color, the –1300~–800 bp region of the promoter were sequenced in 54 accessions with different genetic background. Three variant types of promoter sequence were found in 23 white accessions and named Type-Agol (21 accessions), Type 1 (one accession, No. 308), and Type 2 (one accession, No. 264) (Supplement Figure S5).
Gene annotation indicated that CmaCh15G005270 was predicted to encode an R2-MYB transcription factor homolog with DIVARICATA2 (ATDIV2, At5g04760) and involved in the salt and ABA signaling in Arabidopsis. DIV1 displayed a critical role for determining ventral identities in Antirrhinum corolla [51,52]. qRT-PCR showed that the candidate gene CmaCh15G005270 was highly expressed in flowers and flesh, with the highest expression level in the seed coat at 28 DPA, and then gradually decreased (Figure 4). The results showed CmaCh15G005270 may be the candidate gene in regulating the seed coat color in Cucurbita maxima.

2.6. RNA-Seq Analysis

To further determine the metabolic network and downstream regulatory genes involved in seed coat color formation and seed coat development, RNA-seq was conducted in 0, 7, 14, 21, 28, 35, 42 DPA seed coats in Wuminglv and Agol. Three biological duplicates were conducted at each time course and 312.27 G data was obtained with a Q30 value ranging from 90.57% to 94.27% (Supplement Table S4).
At all the development periods, the number of up-regulated DEGs in Wuminglv was ranged from 322 to 3862, and the number of down-regulated DEGs was from 419 to 3405, compared to Agol. The amount of DEGs was most remarkable at 28 DPA, in accordance with the period of seed coat color formation (Figure 5A,B). GO enrichment analysis showed that biological processes such as metabolic process, physiological process, antioxidant activity, and nuclear acid binding transcription factor activity were significantly enriched during the development of seed coat (Figure 5C). KEGG enrichment analysis indicated phenylalanine metabolism and phenylpropanoid biosynthesis were clustered significantly (Figure 5D). The difference of seed coat color in Wuminglv and Agol was closely related to the de novo biosynthesis pathway of phenylpropanoid.
The expression of some genes involved in flavonoid metabolism was confirmed by qRT-PCR in Wuminglv and Agol. The result of expression was consistent with that of RNA-seq data, indicating the reliability of RNA-seq (Figure 6). PAL, C4H, 4CL, CHI, and DFR involved in phenylpropanoid/flavonoid were displayed remarkable relative expression differences in Wuminglv and Agol. RNA-seq data revealed that the transcript level of CmaCh15G005270 increased gradually and arrived at a peak at 28 DPA. Meanwhile, the expression of this gene in Wuminglv was higher than that in Agol at all periods in seed coat development (Figure 5E). The results of RNA-seq illustrated that the phenylpropanoid/flavonoid metabolism pathway was affected in Wuminglv and Agol.

2.7. Validation of Molecular Marker

Fifty-four accessions (31 accessions of yellow seed coat and 23 accessions of white seed coat) and 222 F2 segregating population and RILs plants (107 accessions of yellow seed coat and 115 accessions of white seed coat) derived from P3 and P4 were used to validate the accuracy of the closest flanking InDel marker S1548 (Supplement Table S5). In these accessions, all of the accessions displaying yellow and white seed coat matched Wuminglv and Agol, suggesting that the accuracy of S1548 marker was 100% and exhibited validation in the molecular assistance breeding program in edible seed pumpkin in Cucurbita maxima.

3. Discussion

Seed coat color is one of the most important traits in oil seed and edible seed crops such as rape, soybean, and watermelon. In some crops, there was correlation between the seed quality, flesh appearance, and seed coat color [9]. In Glycine max, the distribution and variation of metabolites in seeds displayed a linkage with seed coat color [12,53,54]. The isoflavones (including malonylglucosides, glucosides, aglycone, and acetylglucosides) were gradually decreased in green, yellow, black, and brown soybeans. Only anthocyanins (cyanidin-3-O-glucoside, delphinidine-3-O-glucoside, and petunidin-3-O-glucoside) were examined in black seed coat cultivars [53]. The isoflavones were associated with seed maturity while fatty acids, glucose, fructose, and sucrose were less linked with seed maturity [55]. In this study, the insoluble PA and soluble PA in yellow seed coat accession Wuminglv was higher than that in white seed coat accession Agol, indicating that flavonoid metabolisms contributed to the seed coat color formation. The correlation between seed coat color, flavonoids, and food nutrition needs to be further demonstrated in Cucurbita maxima. Whether the seed nutrition of yellow seed coat is significantly different from that of white seed coat also provides some new research prospects for variety breeding.
Flavonoids are one of the largest families of polyphenols, with complicated structural and functional diversity [56]. It is derived from phenylalanine and malonyl-coenzyme A (CoA) [57] and the C6-C3-C6 carbon bridge can be further divided by the nature of C3 element, accordingly producing nine major subgroups including chalcones, aurones, isoflavonoids, flavones, flavonols, anthocyanins, condensed tannins, and phlobaphene pigments [57]. Alanine and phenylalanine metabolism in the seed embryo of Brassica napus was found to be related to the embryonal control of seed coat color [58]. In Brassica species, yellow seed coat accessions exhibit a significantly thinner seed coat than that in black seed coat accessions and displayed a lower hull proportion [58]. Besides the major branch pathways of flavonoid biosynthesis, lignin is another branch of the phenylpropanoid pathway and participates in the formation of secondarily thickened cell walls [54]. In this study, a serial of metabolisms involved in flavonoids pathway including lignin, soluble PA, and insoluble PA displayed different distribution patterns in white and yellow seed coat accessions in Cucurbita maxima. Meanwhile, the secondary cell wall thickening of S-layers and the width of seed coat in Agol were obvious thinner than that in yellow seed coat accession Wuminglv (Figure 2). It was suggested that the mechanism of seed coat color determination and secondary cell wall thickening in Cucurbita maxima is associated with flavonoid metabolism and are similar to those in Brassica and Glycine species.
As oxidization compounds of flavonoid, PAs are thought to be the main pigment of seed coat in Arabidopsis, and it exists in the endothelial layer of the inner integument in seed coat [59,60]. In Arabidopsis, PAs synthesis is accompanied by seed coat development, cell differentiation, and seed maturation. In this study, flavonoid compounds in seed coat layer were determined by histochemical analysis during the seed development. At the beginning of pollination, the seed coat color of Wuminglv and Agol were both white. At 21~28 DPA, PAs were accumulated at the H layer and S layer accompanied by a secondary cell wall thickening at S layer, and the seed coat color turned yellow in Wuminglv. The pigment accumulation and the secondary cell structures of seed coat layer are inseparable [16,24], indicating that WSC may displayed a multiple function in seed coat development. Meanwhile, the development of seed coat in Cucurbita maxima possesses similar characteistics typically associated with that in Cucurbita pepo, indicating the regulatory of seed coat development displays a common pattern in Cucurbita species.
In Arabidopsis and many crops, most of the phenotype of light seed coat color displayed recessive and maternal inheritance against dark seed coat color [61]. However, there are quantitative characters that affect seed coat color and some QTLs that have been identified [30,62,63]. In watermelon, the light cream color was recessive to the brown or red seed coat color, and any of these seed coat colors was recessive to the black color independently [31]. Some of the QTLs and candidate genes have been identified [7,30,32]. In this study, genetic analysis of six-generation population and F2 segregating population derived from different genetic background revealed that the phenotype of seed coat color was a recessive trait controlled by a single gene (Table 2 and Supplement Table S2). In the DMACA staining of 76 accessions, although all the accessions of white seed coat remain white or displayed a little pink, the degree of coloration of yellow seed coat ranged from dark red to light red, indicating that there may be some modification loci associated with pigment content, pigment deposition, or layer structure of the seed coat in Cucurbita maxima.
From the analysis of map-based cloning by F2 segregating population, CmaCh15G005270 encoding R2 MYB transcription factor was predicted to be the candidate gene of white seed coat in Cucurbita maxima. In Arabidopsis, the homolog of CmaCh15G005270 was AtDIV1 (DIVARICATA 1) and AtDIV2. It is most similar to LeMYBI from tomato [52]. Interestingly, loss-of-function mutant div2 was tolerant to salt stress and sensitivity to ABA during seed germination. AtDIV2 played a negative role in salt stress and was required in ABA signal pathway [49]. DIV was studied in the ventral petal of Antirrhinum majus flowers [51,64]. An R2R3 MYB transcription factor AmMYBML1 (Antirrhinum majus MYB MIXTA LIKE 1) very similar to MIXTA was regulated by DIV in association with the B-function proteins MADS-domain transcription factors and control the specification of cell types in the ventral petal [64]. AmMYBML was the orthologous to MIXTA, while AmMYBML2 and AmMYBML3 were the orthologous to PhMYB1, AtMYB16, and AtMYB106. These R2R3-MYBs are not only involved in seed coat color formation and pigment accumulation, but also in other complicated physiological and biochemical regulation processes.
Sequence analysis revealed four main Insertion-Deletions of DNA sequences in the promoter regions. As displayed in 54 accessions, all the alleles exhibiting Insertion-Deletions in the promoter of CmaCh15G005270 displayed a strong linkage with the phenotype of white or yellow seed coat. In many plants, it was an efficient way to regulate gene transcription through sequence variation of promoter. In Sesamum indicum, an AACACACAC-motif in single and duplicated copy in different germplasms displayed strong association with root biomass. It can be recognized by SiMYB181 in 5′-UTR of Big root biomass (BRB) and the loci have been selected by modern breeding [65]. In fruit quality trait studies, R2R3-MYB10 transcription factor were involved in anthocyanin biosynthesis in the color determination in peel and flesh and displayed a general regulatory role in the flavonoid/phenylpropanoid pathway during the ripening of fruits in Fragaria × ananassa (strawberry) [66], Malus domestica (apple) [67,68,69], Pyrus pyrifolia (pear) [70,71], Prunus persica (peach) [72], and other fruits. In the fruit of Prunus persica, PpMYB10.1 is the major regulator of anthocyanin accumulation in red-skinned fruit, and the allelic types of MYB10.1 displayed a close relationship with the intensity of red skin coloration [73]. In the leaves of Prunus persica, PpMYB10.4 were identified in anthocyanin accumulation and 18 single nucleotide polymorphisms (SNP) and a 3-bp InDel were identified within a 2.0-kb region upstream the starting codon [72]. In Malus domestica, five direct tandem repeats of a 23-bp sequence forming a minisatellite-like structure contributed to an allelic rearrangement in the region of promoter and lead to ectopic accumulation of anthocyanins in the flesh and leaves of the plant [69]. In Fragaria, a CACTA-like transposon insertion in the MYB10-2 promoter of red-fleshed accessions was associated with enhanced expression of MYB10 and responsible for enhanced anthocyanin biosynthesis [74]. Recent studies on MYB10-D in synthetic hexaploidy wheat revealed that a 2.4 Mb presence-absence variation (PAV) region was involved in the flavonoid and ABA biosynthesis pathways, accordingly, the seed color of MYB10-D overexpression lines was red and that in the wild type was white [75]. The sequence variation of promoter region in MYB transcription has been reported to contribute to pigment deposition. In many vegetable crops, R2R3-MYBs such as MYB2, MYB75/PAP1, and PAP2 were involved in the anthocyanin regulatory in eggplant, radish, lablab pods, kohlrabi, and cabbages, not only in seed coat [26,76,77]. Therefore, the MYBs would be suitable for metabolic genetic engineering for improvement of fruit quality and vegetable nutrients [78].
MBW complex has been studied including R2R3-MYB, bHLH, and WD40 complex [59]. In Phaseolus vulgaris (Common bean), an R2R3 MYB was associated with seed coat darkening and involved in PAs oxidation and accumulation in the seed coat [79]. LACCASE family in Arabidopsis was involved in oxidative lignin polymerization and it was downstream of MYB46 and MYB58 [54]. A basic helix-loop-helix (bHLH) Noemi transcription factor controls the pigments. Lack of anthocyanins in young leaves and flowers is also associated with a lack of proanthocyanidins in seeds and, most notably, with an extreme reduction in fruit acidity [80,81,82].
Some components of regulation are upstream of MBW complex and are involved in the cell layer differentiation. The ovule identity factor SEEDSTICK (STK, a MADS transcription factor) displayed a connection between cell fate determination, seed coat development, and secondary metabolism as it controls PAs accumulation in the inner layer of seed coat by regulating the expression of BANYULS/ANTHOCYANIDIN REDUCTASE (BAN) though H3K9ac chromatin modification [83,84]. TESTA8 and GLABRA3 composited of WD40-bHLH-MYB complex are also downstream of STK [84,85]. In the promoter of Ruby gene, T-to-C transition on TATA and CAAT boxes was identified to determine the efficiency of the expression [80,81]. Cyanidin 3-(6′-malonyl)-β-glucoside and cyanidin 3-glucoside, accounting for most of the anthocyanin biosynthesis, were repressed in ruby accessions. Whether WSC directly participates in the formation of MBW complex like other R2R3-MYBs or indirectly regulates flavonoid synthesis by regulating downstream genes involved in MBW complex indirectly is still being studied.
The RNA-seq data displayed many genes involved in phenylpropanoid metabolism pathway and its downstream anthocyanin and flavonoid pathway. Meanwhile, RNA-seq data and qRT-PCR showed that WSC played a positive role in yellow seed coat formation. The expression pattern of PAL, C4H, 4CL, CHI, and DFR were extremely different in Wuminglv and Agol at the stage of seed coat color formation. It can be confirmed that the color and structure of seed coat are closely related and affected by the same regulatory pathway. In this metabolic network, we have not only discovered genes related to pigment accumulation, but also some genes related to lignin metabolism. The target genes downstream of CmaCh15G005270 are to be further identified.

4. Materials and Methods

4.1. Plant Materials

Inbred accessions Wuminglv (yellow seed coat color, P1), Agol (white seed coat color, P2), Lizimanhui (yellow seed coat color, P3), and Diaogua (white seed coat color, P4) derived from the Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences (Beijing, China) were used as the parents. F1 was obtained through the hybridization of Wuminglv as female parents and Agol as male parents. A six-generation population including P1, P2, F1, F2, BC1, and BC1 was constructed by cross and self-pollination in Wuminglv and Agol. The 2798 plants in F2 segregating population derived from Wuminglv and Agol were used for the fine mapping. The 222 plants in F2 segregating population and RILs derived from Lizimanhui and Diaogua and 54 inbreds displayed different seed coat color were used to examine the accuracy of molecular markers.

4.2. Phenotypic Characterization

The seed coat color was categorized as two types by visual observation directly: the white seed type, or the yellow seed type. Due to the complexity of color concept, the L*, a*, and b* values were obtained by means of chromometer (HP-200, Hanpu, Changzhou), and the comprehensive chromaticity E was obtained by using the formula to carry out the quantitative comparison of seed coat color as E = [(L*)2 + (a*)2 + (b*)2]1/2.

4.3. Histochemical Analysis

Female flowers with conspicuous ovary were artificially pollinated and the date was recorded as 0 DPA (Day post anthesis). Three young fruits were harvested at 0, 7, 14, 21, 28, 35, and 42 DPA. The seeds were immediately collected and embedded in FAA solution at 4 °C (40% formaldehyde/70% glacial acetic acid/50% alcohol = 5:5:90). Paraffin sections were processed in Epsilon Company (Beijing, China) and toluidine blue was used for staining. The paraffin section of seed coat was observed by microscope Olympus BX51 (Tokyo, Japan). The width and length of each cell layer of the seed coat were measured by ImageJ.

4.4. DNA and RNA Extraction

Young leaf samples of P1, P2, F1, and all of the F2 plants with unique number were collected and frozen in liquid nitrogen. DNA extraction was performed as described by Xiang et al. [42]. For RNA extraction, the sample of seed coat were collected after removing the pulp in different development stages (0, 7, 14, 21, 28, 35, 42 DPA). Different tissues of older root, stem, young leaf, petals of anthesis, and ripening flesh were collected. RNA extraction was processed by MiniBEST Universal RNA Extraction Kit (TaKaRa, Japan, 9769).

4.5. Genetic Analysis

A six-generation population (derived from Wuminglv and Agol including P1, P2, F1, F2, BC1, BC2) and another F2 segregating population (derived from Lizimanhui and Diaogua and including P3, P4, F1, F2) were used to study the genetic analysis of seed coat color. The Chi-square test of fit was performed against the expected 3:1 segregation ratio in the F2 segregating population.

4.6. Map-Based Cloning of wsc Gene

The Whole genome sequencing was conductded in Wuminglv and Agol to develop the InDel markers for primary mapping (Biomarker Techonology Company, Beijing, China). Three hundred and twenty-four pairs of primers evenly distributed on each chromosome were designed by Primer 3.0 (http://bioinfo.ut.ee/primer3-0.4.0/, accessed on 6 February 2021) followed the reference genome (http://cucurbitgenomics.org/, accessed on 6 February 2021). The PCR program was conducted as follows: step 1, 95 °C 3 min; step 2, 95 °C 15 s, 54~56 °C 15 s, 72 °C 30 s, repeat 35 cycles; step 3, 72 °C 5 min. The DNA templates were amplified by high-fidelity Polymerase (P115, Vazyme, Nanjing, China) and the polymorphism of PCR products was determined by 8% Polyacrylamide gel electrophoresis (PAGE) in P1, P2, and F1. For fine mapping, DNA templates of 2798 F2 plants were amplified with primary flanking markers S1517 and S1528, followed by 17 InDel markers and nine dCAPS markers designed in the 315.4-kb candidate region for recombinants screening. The phenotype of seed coat color was confirmed in recombinants and the homozygous were selected to confirm the final candidate region.

4.7. Candidate Gene Prediction and Validation of Molecular Marker

The full lengths of four genes in the candidate region were amplified and sequenced from Wuminglv and Agol. The promoter region of candidate gene was amplified in 54 inbred accessions and sequenced in Sangon Biotech (Shanghai, China). The flanking InDel marker S1548 was used for validation of molecular markers in 54 inbred accessions and 222 lines of F2 segregating population and RILs derived from P3 and P4.

4.8. RNA-Seq and Spatiotemporal Expression Analysis

The samples of seed coat were collected after removing the pulp in different development stages (0, 7, 14, 21, 28, 35, 42 DPA). RNA-seq was conducted on Illumina platform with three replicates for each sample in Biomarker Techonology Company (Beijing, China). The sequencing Reads were compared to the reference genome using TopHat2 [86]. The transcriptome assemblies were performed by Cufflflinks [87]. Then, according to the expression level of genes in seed coat samples in different periods, DESeq software was used for differential expression analysis, and differentially expressed genes were screened through indicators FDR = 0.01 and FC = 2 [88]. Finally, the annotation information of each differentially expressed gene set was extracted for bioinformatics analysis.
For gene expression analysis, ChamQ Universal SYBR qPCR Master Mix (Q711, Vazyme, Nanjing, China) was used and the expression of candidate genes was determined by quantitative Real-time PCR (Bio-Rad, C1000 Touch, Foster City, CA, USA). The data was analyzed by Bio-Red CFX Manager (Version 3.1, Bio-Rad).

4.9. DMACA Staining and PA Determination

DMACA staining and PA determination were followed as in Pang et al. [23]. Total soluble PA content was calculated spectrophotometrically after reaction with DMACA reagent (0.2% wt/vol DMACA in methanol-3N HCl) at 640 nm by Spectrophotometer (Unico, UV/VIS 2802PC, Princeton, NJ, USA with (+)-catechin (Solarbio, 154-23-4, Beijing, China) as standard. For total insoluble PA analysis, Absorbance values were converted into PA equivalents using a standard curve of procyanidin B1 (4852-22-6, MYM, Beijing, China).

4.10. Quantitative Determination for Lignin

Quantitative determination for lignin was followed as in [89]. Aliquots of 20 mg of seed coat (six replicates per individual sample) of Wuminglv and Agol were mixed thiogylcolic acid (MYM, Beijing, China) and the standard lignin (Biotapped, Beijing, China) was used to generate the calibration curves.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/22/6/2972/s1.

Author Contributions

Conceptualization, C.W. and Y.D.; methodology, Y.D. and Y.S.; validation, Y.S. and Y.D.; Data curation, Y.S., M.Z., and W.M.; investigation, Y.S., Q.S., T.S., and C.X.; writing—original draft preparation, Y.D. and Y.S.; writing—review and editing, Y.D. and C.W.; project administration, C.W.; funding acquisition, Y.D. and C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 32002037; the Science and Technology Innovation Program of the Chinese Academy of Agricultural Science, grant number CAAS-ASTIP-IVFCAAS; the Central Public-interest Scientific Institution Basal Research Funding; Special Fund for Agro-scientific Research in the Public Interest (201303112).

Acknowledgments

We thank Kailiang Bo (Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences) for providing important advice on sequencing and experiment designing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chomicki, G.; Schaefer, H.; Renner, S.S. Origin and domestication of Cucurbitaceae crops: Insights from phylogenies, genomics and archaeology. New Phytol. 2019, 226, 1240–1255. [Google Scholar] [CrossRef] [Green Version]
  2. Kostecka-Gugała, A.; Kruczek, M.; Ledwożyw-Smoleń, I.; Kaszycki, P. Antioxidants and health-beneficial nutrients in fruits of eighteen Cucurbita cultivars: Analysis of diversity and dietary implications. Molecules 2020, 25, 1792. [Google Scholar] [CrossRef] [Green Version]
  3. Nawirska-Olszańska, A.; Kita, A.; Biesiada, A.; Sokół-Łętowska, A.; Kucharska, A.Z. Characteristics of antioxidant activity and composition of pumpkin seed oils in 12 cultivars. Food Chem. 2013, 139, 155–161. [Google Scholar] [CrossRef]
  4. Amin, M.Z.; Islam, T.; Mostofa, F.; Uddin, M.J.; Rahman, M.; Satter, M.A. Comparative assessment of the physicochemical and biochemical properties of native and hybrid varieties of pumpkin seed and seed oil (Cucurbita maxima Linn.). Heliyon 2019, 5, e02994. [Google Scholar] [CrossRef] [Green Version]
  5. Stevenson, D.G.; Eller, F.J.; Wang, L.; Jane, J.L.; Wang, A.T.; Inglett, G.E. Oil and tocopherol content and composition of pumpkin seed oil in 12 cultivars. J. Agric. Food Chem. 2007, 55, 4005–4013. [Google Scholar] [CrossRef] [Green Version]
  6. Rezig, L.; Chouaibi, M.; Msaada, K.; Hamdi, S. Chemical composition and profile characterisation of pumpkin (Cucurbita maxima) seed oil. Ind. Crops Prod. 2012, 37, 82–87. [Google Scholar] [CrossRef]
  7. Li, B.; Lu, X.; Gebremeskel, H.; Zhao, S.; He, N.; Yuan, P.; Gong, C.; Mohammed, U.; Liu, W. Genetic mapping and discovery of the candidate gene for black seed coat color in Watermelon (Citrullus lanatus). Front. Plant Sci. 2020, 10, 1689. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, C.; Jun, J.H.; Dixon, R.A. MYB5 and MYB14 play pivotal roles in seed coat polymer biosynthesis in Medicago truncatula. Plant Physiol. 2014, 165, 1424–1439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Zhai, Y.; Yu, K.; Cai, S.; Hu, L.; Amoo, O.; Xu, L.; Yang, Y.; Ma, B.; Jiao, Y.; Zhang, C.; et al. Targeted mutagenesis of BnTT8 homologs controls yellow seed coat development for effective oil production in Brassica napus L. Plant Biotechnol. J. 2020, 18, 1153–1168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Guo, Y.; Gao, M.; Liang, X.; Xu, M.; Liu, X.; Zhang, Y.; Liu, X.; Liu, J.; Gao, Y.; Qu, S.; et al. Quantitative trait loci for seed size variation in Cucurbits—A Review. Front. Plant Sci. 2020, 11, 304. [Google Scholar] [CrossRef]
  11. Heneidak, S.; Khalik, K. Seed coat diversity in some tribes of Cucurbitaceae: Implications for taxonomy and species iden-tification. Acta Bot. Bras. 2015, 29, 129–142. [Google Scholar] [CrossRef] [Green Version]
  12. Lee, J.; Hwang, Y.S.; Kim, S.T.; Yoon, W.B.; Han, W.Y.; Kang, I.K.; Choung, M.G. Seed coat color and seed weight contribute differential responses of targeted metabolites in soybean seeds. Food Chem. 2017, 214, 248–258. [Google Scholar] [CrossRef]
  13. Lepiniec, L.; Debeaujon, I.; Routaboul, J.M.; Baudry, A.; Pourcel, L.; Nesi, N.; Caboche, M. Genetics and biochemistry of seed flavonoids. Annu. Rev. Plant Biol. 2006, 57, 405–430. [Google Scholar] [CrossRef]
  14. Demonsais, L.; Utz-Pugin, A.; Loubéry, S.; Lopez-Molina, L. Identification of tannic cell walls at the outer surface of the endosperm upon Arabidopsis seed coat rupture. Plant J. 2020, 104, 567–580. [Google Scholar] [CrossRef]
  15. Bezold, T.N.; Loy, J.; Minocha, S.C. Changes in the cellular content of polyamines in different tissues of seed and fruit of a normal and a hull-less seed variety of pumpkin during development. Plant Sci. 2003, 164, 743–752. [Google Scholar] [CrossRef]
  16. Li, X.; Ge, Y.; Gui, M.; Cui, C.S.; Qu, S.P. Histological analysis of hulled and hull-less squash seed. Coat. Int. J. Agric. Biol. 2013, 15, 1047–1050. [Google Scholar]
  17. Zhao, J. Flavonoid transport mechanisms: How to go, and with whom. Trends Plant Sci. 2015, 20, 576–585. [Google Scholar] [CrossRef]
  18. Malenčić, D.; Popović, M.; Miladinović, J. Phenolic content and antioxidant properties of Soybean (Glycine max (L.) Merr.) Seeds. Molecules 2007, 12, 576–581. [Google Scholar] [CrossRef]
  19. Qu, C.; Fu, F.; Lu, K.; Zhang, K.; Wang, R.; Xu, X.; Wang, M.; Lu, J.; Wan, H.; Zhanglin, T.; et al. Differential accumulation of phenolic compounds and expression of related genes in black- and yellow-seeded Brassica napus. J. Exp. Bot. 2013, 64, 2885–2898. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Bhargava, A.; Mansfield, S.D.; Hall, H.C.; Douglas, C.J.; Ellis, B.E. MYB75 functions in regulation of secondary cell wall formation in the Arabidopsis inflorescence stem. Plant Physiol. 2010, 154, 1428–1438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Dong, N.; Lin, H. Contribution of phenylpropanoid metabolism to plant development and plant–environment interactions. J. Integr. Plant Biol. 2021, 63, 180–209. [Google Scholar] [CrossRef] [PubMed]
  22. Pang, Y.; Peel, G.J.; Sharma, S.B.; Tang, Y.; Dixon, R.A. A transcript profiling approach reveals an epicatechin-specific glucosyltransferase expressed in the seed coat of Medicago truncatula. Proc. Natl. Acad. Sci. USA 2008, 105, 14210–14215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pang, Y.; Wenger, J.P.; Saathoff, K.; Peel, G.J.; Wen, J.; Huhman, D.; Allen, S.N.; Tang, Y.; Cheng, X.; Tadege, M.; et al. A WD40 repeat protein from Medicago truncatula is necessary for tissue-specific anthocyanin and proanthocyanidin biosynthesis but not for trichome development. Plant Physiol. 2009, 151, 1114–1129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. DeBolt, S.; Scheible, W.R.; Schrick, K.; Auer, M.; Beisson, F.; Bischoff, V.; Bouvier-Navé, P.; Carroll, A.; Hematy, K.; Li, Y.; et al. Mutations in UDP-Glucose: Sterol glucosyltransferase in Arabidopsis cause transparent testa phenotype and suberization defect in seeds. Plant Physiol. 2009, 151, 78–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shirley, B.W.; Kubasek, W.L.; Storz, G.; Bruggemann, E.; Koornneef, M.; Ausubel, F.M.; Goodman, H.M. Analysis of Arabidopsis mutants deficient in flavonoid biosynthesis. Plant J. 1995, 8, 659–671. [Google Scholar] [CrossRef] [PubMed]
  26. Teng, S.; Keurentjes, J.; Bentsink, L.; Koornneef, M.; Smeekens, S. Sucrose-specific induction of anthocyanin biosynthesis in Arabidopsis requires the MYB75/PAP1 gene. Plant Physiol. 2005, 139, 1840–1852. [Google Scholar] [CrossRef] [Green Version]
  27. Zheng, T.; Tan, W.; Yang, H.; Zhang, L.; Li, T.; Liu, B.; Zhang, D.; Lin, H. Regulation of anthocyanin accumulation via MYB75/HAT1/TPL-mediated transcriptional repression. PLoS Genet. 2019, 15, e1007993. [Google Scholar] [CrossRef]
  28. Jaakola, L. New insights into the regulation of anthocyanin biosynthesis in fruits. Trends Plant Sci. 2013, 18, 477–483. [Google Scholar] [CrossRef] [Green Version]
  29. McKay, J.W. Factor Interaction in Citrullus. Jpn. J. Crop Sci. 1936, 27, 110–112. [Google Scholar] [CrossRef]
  30. Li, B.; Lu, X.; Dou, J.; Aslam, A.; Gao, L.; Zhao, S.; He, N.; Liu, W. Construction of a high-density genetic map and mapping of fruit traits in Watermelon (Citrullus Lanatus L.) based on Whole-genome resequencing. Int. J. Mol. Sci. 2018, 19, 3268. [Google Scholar] [CrossRef] [Green Version]
  31. Paudel, L.; Clevenger, J.; McGregor, C. Chromosomal locations and interactions of four loci associated with seed coat color in Watermelon. Front. Plant Sci. 2019, 10, 788. [Google Scholar] [CrossRef]
  32. Nath, P.; Khandelwal, R.C. Genetics of seed-coat colour in Citrullus lanatus (Thunb.) Mansf. Theor. Appl. Genet. 1978, 53, 29–31. [Google Scholar] [CrossRef]
  33. Brown, R.N.; Myers, J.R. A genetic map of Squash (Cucurbita sp.) with randomly amplified polymorphic DNA markers and morphological markers. J. Am. Soc. Hortic. Sci. 2002, 127, 568–575. [Google Scholar] [CrossRef] [Green Version]
  34. Zraidi, A.; Lelley, T.; Lebeda, A.; Paris, H.S. Genetic map for pumpkin Cucurbita pepo using random amplified polymorphic DNA markers. In Proceedings of the Meeting on Progress in Cucurbit Genetics and Breeding Research Cucurbitaceae, Olomouc, Czech Republic, 12–17 July 2004. [Google Scholar]
  35. Gong, L.; Stift, G.; Kofler, R.; Pachner, M.; Lelley, T. Microsatellites for the genus Cucurbita and an SSR-based genetic linkage map of Cucurbita pepo L. Theor. Appl. Genet. 2008, 117, 37–48. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, G.; Ren, Y.; Sun, H.; Guo, S.; Zhang, F.; Zhang, J.; Zhang, H.; Jia, Z.; Fei, Z.; Xu, Y.; et al. A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.). BMC Genom. 2015, 16, 1–13. [Google Scholar] [CrossRef] [Green Version]
  37. Montero-Pau, J.; Blanca, J.; Esteras, C.; Martínez-Pérez, E.M.; Gómez, P.; Monforte, A.J.; Cañizares, J.; Picó, B. An SNP-based saturated genetic map and QTL analysis of fruit-related traits in Zucchini using Genotyping-by-sequencing. BMC Genom. 2017, 18, 1–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Sun, H.; Wu, S.; Zhang, G.; Jiao, C.; Guo, S.; Ren, Y.; Zhang, J.; Zhang, H.; Gong, G.; Jia, Z.; et al. Karyotype stability and unbiased fractionation in the Paleo-Allotetraploid Cucurbita genomes. Mol. Plant 2017, 10, 1293–1306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Zhong, Y.J.; Zhou, Y.Y.; Li, J.X.; Yu, T.; Wu, T.Q.; Luo, J.N.; Luo, S.B.; Huang, H.X. A high-density linkage map and QTL mapping of fruit-related traits in pumpkin (Cucurbita moschata Duch.). Sci. Rep. 2017, 7, 12785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Montero-Pau, J.; Blanca, J.; Bombarely, A.; Ziarsolo, P.; Esteras, C.; Martí-Gómez, C.; Ferriol, M.; Gómez, P.; Jamilena, M.; Muelleret, L.; et al. De novo assembly of the zucchini genome reveals a whole-genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol. J. 2018, 16, 1161–1171. [Google Scholar] [CrossRef] [Green Version]
  41. Nguyen, N.N.; Kim, M.; Jung, J.K.; Shim, E.J.; Chung, S.M.; Park, Y.; Lee, G.P.; Sim, S.C. Genome-wide SNP discovery and core marker sets for assessment of genetic variations in cultivated pumpkin (Cucurbita spp.). Hortic. Res. 2020, 7, 121. [Google Scholar] [CrossRef] [PubMed]
  42. Xiang, C.; Duan, Y.; Li, H.; Ma, W.; Huang, S.; Sui, X.; Zhang, Z.; Wang, C. A high-density EST-SSR-based genetic map and QTL analysis of dwarf trait in Cucurbita pepo L. Int. J. Mol. Sci. 2018, 19, 3140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kaźmińska, K.; Hallmann, E.; Korzeniewska, A.; Niemirowicz-Szczytt, K.; Bartoszewski, G. Identification of fruit-associated QTLs in Winter Squash (Cucurbita maxima Duchesne) using recombinant inbred lines. Genes 2020, 11, 419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ramos, A.; Fu, Y.; Michael, V.; Meru, G. QTL-seq for identification of loci associated with resistance to Phytophthora crown rot in squash. Sci. Rep. 2020, 10, 5326. [Google Scholar] [CrossRef]
  45. Sáez, C.; Martínez, C.; Montero-Pau, J.; Esteras, C.; Sifres, A.; Blanca, J.; Ferriol, M.; López, C.; Picó, B. A major QTL located in chromosome 8 of Cucurbita moschata is responsible for resistance to Tomato leaf curl New Delhi virus. Front. Plant Sci. 2020, 11, 207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Vogel, G.; LaPlant, K.E.; Mazourek, M.; Gore, M.A.; Smart, C.D. A combined BSA-Seq and linkage mapping approach identifies genomic regions associated with Phytophthora root and crown rot resistance in squash. Theor. Appl. Genet. 2021, 1–17. [Google Scholar] [CrossRef]
  47. Wang, Y.; Wang, C.; Han, H.; Luo, Y.; Wang, Z.; Yan, C.; Xu, W.; Qu, S. Construction of a high-density genetic map and analysis of seed-related traits using specific length amplified fragment sequencing for Cucurbita maxima. Front. Plant Sci. 2020, 10, 1782. [Google Scholar] [CrossRef]
  48. Abeynayake, S.W.; Panter, S.; Mouradov, A.; Spangenberg, G. A high-resolution method for the localization of proanthocyanidins in plant tissues. Plant Methods 2011, 7, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Fang, Q.; Wang, Q.; Mao, H.; Xu, J.; Wang, Y.; Hu, H.; He, S.; Tu, J.; Cheng, C.; Tian, G.; et al. AtDIV2, an R-R-type MYB transcription factor of Arabidopsis, negatively regulates salt stress by modulating ABA signaling. Plant Cell Rep. 2018, 37, 1499–1511. [Google Scholar] [CrossRef]
  50. Choi, Y.; Kim, S.; Kwon, C.; Yun, H.S. Regulation of cellular VAMP721/722 abundance in Arabidopsis. Plant Signal. Behav. 2019, 14, e1632690. [Google Scholar] [CrossRef]
  51. Galego, L. Role of DIVARICATA in the control of dorsoventral asymmetry in Antirrhinum flowers. Genes Dev. 2002, 16, 880–891. [Google Scholar] [CrossRef] [Green Version]
  52. Machemer, K.; Shaiman, O.; Salts, Y.; Shabtai, S.; Sobolev, I.; Belausov, E.; Grotewold, E.; Barg, R. Interplay of MYB factors in differential cell expansion, and consequences for tomato fruit development. Plant J. 2011, 68, 337–350. [Google Scholar] [CrossRef] [PubMed]
  53. Cho, K.M.; Ha, T.J.; Lee, Y.B.; Seo, W.D.; Kim, J.Y.; Ryu, H.W.; Jeong, S.H.; Kang, Y.M.; Lee, J.H. Soluble phenolics and antioxidant properties of soybean (Glycine max L.) cultivars with varying seed coat colours. J. Funct. Foods 2013, 5, 1065–1076. [Google Scholar] [CrossRef]
  54. Zhao, Q.; Nakashima, J.; Chen, F.; Yin, Y.; Fu, C.; Yun, J.; Shao, H.; Wang, X.; Wang, Z.Y.; Dixon, R.A. LACCASE is necessary and nonredundant with PEROXIDASE for lignin polymerization during vascular development in Arabidopsis. Plant Cell 2013, 25, 3976–3987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lee, J.; Hwang, Y.S.; Chang, W.S.; Moon, J.K.; Choung, M.G. Seed maturity differentially mediates metabolic responses in black soybean. Food Chem. 2013, 141, 2052–2059. [Google Scholar] [CrossRef]
  56. Wen, W.; Alseekh, S.; Fernie, A.R. Conservation and diversification of flavonoid metabolism in the plant kingdom. Curr. Opin. Plant Biol. 2020, 55, 100–108. [Google Scholar] [CrossRef]
  57. Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001, 126, 485–493. [Google Scholar] [CrossRef] [Green Version]
  58. Wang, F.; He, J.; Shi, J.; Zheng, T.; Xu, F.; Wu, G.; Liu, R.; Liu, S. Embryonal control of yellow seed coat locus ECY1 is related to Alanine and Phenylalanine metabolism in the seed embryo of Brassica napus. G3 Genes Genomes Genet. 2016, 6, 1073–1081. [Google Scholar] [CrossRef] [Green Version]
  59. Debeaujon, I.; Léon-Kloosterziel, K.M.; Koornneef, M. Influence of the testa on seed dormancy, germination, and longevity in Arabidopsis. Plant Physiol. 2000, 122, 403–414. [Google Scholar] [CrossRef] [Green Version]
  60. Haughn, G.; Chaudhury, A. Genetic analysis of seed coat development in Arabidopsis. Trends Plant Sci. 2005, 10, 472–477. [Google Scholar] [CrossRef] [PubMed]
  61. Loubéry, S.; De Giorgi, J.; Utz-Pugin, A.; Demonsais, L.; Lopez-Molina, L. A maternally deposited endosperm cuticle contributes to the physiological defects of transparent testa seeds. Plant Physiol. 2018, 177, 1218–1233. [Google Scholar] [CrossRef] [Green Version]
  62. Erfatpour, M.; Navabi, A.; Pauls, K.P. Mapping the non-darkening trait from ‘Wit-rood boontje’ in bean (Phaseolus vulgaris). Theor. Appl. Genet. 2018, 131, 1331–1343. [Google Scholar] [CrossRef] [PubMed]
  63. Junk-Knievel, D.C.; Vandenberg, A.; Bett, K.E. Slow darkening in Pinto Bean (Phaseolus vulgaris L.) seed coats is controlled by a single major gene. Crops Sci. 2008, 48, 189–193. [Google Scholar] [CrossRef]
  64. Perez-Rodriguez, M.; Jaffe, F.W.; Butelli, E.; Glover, B.J.; Martin, C. Development of three different cell types is associated with the activity of a specific MYB transcription factor in the ventral petal of Antirrhinum majus flowers. Development 2005, 132, 359–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Dossa, K.; Zhou, R.; Li, D.; Liu, A.; Qin, L.; Mmadi, M.A.; Su, R.; Zhang, Y.; Wang, J.; Gao, Y.; et al. A novel motif in the 5′-UTR of an orphan gene ‘Big Root Biomass’ modulates root biomass in sesame. Plant Biotechnol. J. 2020. [Google Scholar] [CrossRef]
  66. Medina-Puche, L.; Cumplido-Laso, G.; Amil-Ruiz, F.; Hoffmann, T.; Ring, L.; Rodríguez-Franco, A.; Caballero, J.L.; Schwab, W.; Muñoz-Blanco, J.; Blanco-Portales, R. MYB10 plays a major role in the regulation of flavonoid/phenylpropanoid metabolism during ripening of Fragaria × ananassa fruits. J. Exp. Bot. 2014, 65, 401–417. [Google Scholar] [CrossRef] [Green Version]
  67. Takos, A.M.; Jaffé, F.W.; Jacob, S.R.; Bogs, J.; Robinson, S.P.; Walker, A.R. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 2006, 142, 1216–1232. [Google Scholar] [CrossRef] [Green Version]
  68. Chagné, D.; Carlisle, C.M.; Blond, C.; Volz, R.K.; Whitworth, C.J.; Oraguzie, N.C.; Crowhurst, R.N.; Allan, A.C.; Espley, R.V.; Hellens, R.P.; et al. Mapping a candidate gene (MdMYB10) for red flesh and foliage colour in apple. BMC Genom. 2007, 8, 212. [Google Scholar] [CrossRef] [Green Version]
  69. Espley, R.V.; Brendolise, C.; Chagné, D.; Kutty-Amma, S.; Green, S.; Volz, R.; Putterill, J.; Schouten, H.J.; Gardiner, S.E.; Hellens, R.P.; et al. Multiple repeats of a promoter segment causes transcription factor autoregulation in red apples. Plant Cell 2009, 21, 168–183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Feng, S.; Wang, Y.; Yang, S.; Xu, Y.; Chen, X. Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10. Planta 2010, 232, 245–255. [Google Scholar] [CrossRef]
  71. Bai, S.; Tao, R.; Tang, Y.; Yin, L.; Ma, Y.; Ni, J.; Yan, X.; Yang, Q.; Wu, Z.; Zeng, Y.; et al. BBX16, a B-box protein, positively regulates light-induced anthocyanin accumulation by activating MYB10 in red pear. Plant Biotechnol. J. 2019, 17, 1985–1997. [Google Scholar] [CrossRef] [Green Version]
  72. Zhou, Y.; Zhou, H.; Lin-Wang, K.; Vimolmangkang, S.; Espley, R.V.; Wang, L.; Allan, A.C.; Han, Y. Transcriptome analysis and transient transformation suggest an ancient duplicated MYB transcription factor as a candidate gene for leaf red coloration in peach. BMC Plant Biol. 2014, 14, 388. [Google Scholar] [CrossRef] [Green Version]
  73. Tuan, P.A.; Bai, S.; Yaegaki, H.; Tamura, T.; Hihara, S.; Moriguchi, T.; Oda, K. The crucial role of PpMYB10.1 in anthocyanin accumulation in peach and relationships between its allelic type and skin color phenotype. BMC Plant Biol. 2015, 15, 280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Castillejo, C.; Waurich, V.; Wagner, H.; Ramos, R.; Oiza, N.; Muñoz, P.; Triviño, J.C.; Caruana, J.; Liu, Z.; Cobo, N.; et al. Allelic variation of MYB10 is the major force controlling natural variation of skin and flesh color in Strawberry (Fragaria spp.) fruit. Plant Cell 2020, 32, 3723–3749. [Google Scholar] [CrossRef] [PubMed]
  75. Lang, J.; Fu, Y.; Zhou, Y.; Cheng, M.; Deng, M.; Li, M.; Zhu, T.; Yang, J.; Guo, X.; Gui, L.; et al. Myb10-D confers PHS-3D resistance to pre-harvest sprouting by regulating NCED in ABA biosynthesis pathway of wheat. New Phytol. 2021. [Google Scholar] [CrossRef] [PubMed]
  76. Maier, A.; Schrader, A.; Kokkelink, L.; Falke, C.; Welter, B.; Iniesto, E.; Rubio, V.; Uhrig, J.F.; Hülskamp, M.; Hoecker, U. Light and the E3 ubiquitin ligase COP1/SPA control the protein stability of the MYB transcription factors PAP1 and PAP2 involved in anthocyanin accumulation in Arabidopsis. Plant J. 2013, 74, 638–651. [Google Scholar] [CrossRef] [PubMed]
  77. Jun, J.H.; Liu, C.; Xiao, X.; Dixon, R.A. The transcriptional repressor MYB2 regulates both spatial and temporal patterns of proanthocyandin and anthocyanin pigmentation in Medicago truncatula. Plant Cell 2015, 27, 2860–2879. [Google Scholar] [CrossRef] [Green Version]
  78. Naing, A.H.; Kim, C.K. Roles of R2R3-MYB transcription factors in transcriptional regulation of anthocyanin biosynthesis in horticultural plants. Plant Mol. Biol. 2018, 98, 1–18. [Google Scholar] [CrossRef]
  79. Erfatpour, M.; Pauls, K.P. A R2R3-MYB gene-based marker for the non-darkening seed coat trait in pinto and cranberry beans (Phaseolus vulgaris L.) derived from ‘Wit-rood boontje’. Theor. Appl. Genet. 2020, 133, 1977–1994. [Google Scholar] [CrossRef] [Green Version]
  80. Butelli, E.; Licciardello, C.; Zhang, Y.; Liu, J.; Mackay, S.; Bailey, P.; Reforgiato-Recupero, G.; Martin, C. Retrotransposons control fruit-specific, cold-dependent accumulation of anthocyanins in Blood Oranges. Plant Cell 2012, 24, 1242–1255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Butelli, E.E.; Garcia-Lor, A.; Licciardello, C.; Licciardello, G.; Hill, L.; Recupero, G.R.; Keremane, M.L.; Ramadugu, C.; Krueger, R.; Xu, Q.; et al. Changes in anthocyanin production during domestication of Citrus. Plant Physiol. 2017, 173, 2225–2242. [Google Scholar] [CrossRef] [Green Version]
  82. Butelli, E.; Licciardello, C.; Ramadugu, C.; Durand-Hulak, M.; Celant, A.; Reforgiato, R.G.; Froelicher, Y.; Martin, C. Noemi controls production of flavonoid pigments and fruit acidity and illustrates the domestication routes of modern citrus varieties. Curr. Biol. 2019, 29, 158–164.e2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Xie, D.Y.; Sharma, S.B.; Paiva, N.L.; Ferreira, D.; Dixon, R.A. Role of anthocyanidin reductase, encoded by BANYULS in plant flavonoid biosynthesis. Science 2003, 299, 396–399. [Google Scholar] [CrossRef] [Green Version]
  84. Ezquer, I.; Mizzotti, C.; Nguema-Ona, E.; Gotté, M.; Beauzamy, L.; Viana, V.E.; Dubrulle, N.; Costa de Oliveira, A.; Caporali, E.; Koroney, A.S. The developmental regulator SEEDSTICK controls structural and mechanical properties of the Arabidopsis seed coat. Plant Cell 2016, 28, 2478–2492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Mizzotti, C.; Ezquer, I.; Paolo, D.; Rueda-Romero, P.; Guerra, R.F.; Battaglia, R.; Rogachev, I.; Aharoni, A.; Kater, M.M.; Caporali, E. SEEDSTICK is a master regulator of development and metabolism in the Arabidopsis seed coat. PLoS Genet. 2014, 10, e1004856. [Google Scholar] [CrossRef]
  86. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R13–R36. [Google Scholar] [CrossRef] [Green Version]
  87. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Dyckmans, J.; Flessa, H.; Brinkmann, K.; Mai, C.; Polle, A. Carbon and nitrogen dynamics in acid detergent fibre lignins of beech (Fagus sylvatica L.) during the growth phase. Plant Cell Environ. 2002, 25, 469–478. [Google Scholar] [CrossRef]
Ijms 22 02972 g001 550
Figure 1. The phenotype of seed coat color and pigment accumulation in Cucurbita maxima. (A) Fruit and seed coat color in Wuminglv and Agol. The upper bar = 5 cm, the lower bar = 2 cm. (B) DMACA staining in Wuminglv (P1), Agol (P2), F1, F1′, Lizimanhui (P3), and Diaogua (P4). Bar = 2 cm. (C) Seed coat color formation during the seed development at different stage after pollination. Dynamic changes of seed coat color at 0, 7, 14, 21, 28, 35, 42 DPA (days post anthesis). Bar = 2 cm. (D) Soluble and insoluble PA determination of mature dry seeds in Wuminglv and Agol. (E) Soluble and insoluble PA determination of flesh seeds at 7, 14, 21, 28, 35, 42 DPA in Wuminglv and Agol.
Figure 1. The phenotype of seed coat color and pigment accumulation in Cucurbita maxima. (A) Fruit and seed coat color in Wuminglv and Agol. The upper bar = 5 cm, the lower bar = 2 cm. (B) DMACA staining in Wuminglv (P1), Agol (P2), F1, F1′, Lizimanhui (P3), and Diaogua (P4). Bar = 2 cm. (C) Seed coat color formation during the seed development at different stage after pollination. Dynamic changes of seed coat color at 0, 7, 14, 21, 28, 35, 42 DPA (days post anthesis). Bar = 2 cm. (D) Soluble and insoluble PA determination of mature dry seeds in Wuminglv and Agol. (E) Soluble and insoluble PA determination of flesh seeds at 7, 14, 21, 28, 35, 42 DPA in Wuminglv and Agol.
Ijms 22 02972 g001
Ijms 22 02972 g002 550
Figure 2. Histochemical analysis of the seed coat at different stage of seed development. (A) Morphological structures of seed coat. Dynamic changes of seed coat in Wuminglv and Agol at 0 (a, b), 7 (c, d), 14 (e, f), 21(g, h), 28 (i, j), 35 (k, l) DPA. Different characters represented cell layers as Epidermis (E), Hypodermis (H), Sclerenchyma (S), Parenchyma (P), and Chlorenchyma (C). Bar = 100 μm. (B) The statistics analysis of cell layers of hypodermis (a, b, c), sclerenchyma (d, e, f), and parenchyma (g, h, l).
Figure 2. Histochemical analysis of the seed coat at different stage of seed development. (A) Morphological structures of seed coat. Dynamic changes of seed coat in Wuminglv and Agol at 0 (a, b), 7 (c, d), 14 (e, f), 21(g, h), 28 (i, j), 35 (k, l) DPA. Different characters represented cell layers as Epidermis (E), Hypodermis (H), Sclerenchyma (S), Parenchyma (P), and Chlorenchyma (C). Bar = 100 μm. (B) The statistics analysis of cell layers of hypodermis (a, b, c), sclerenchyma (d, e, f), and parenchyma (g, h, l).
Ijms 22 02972 g002
Ijms 22 02972 g003 550
Figure 3. Map-based cloning of WSC. (A) Map-based cloning of WSC to 60.4 kb on Chromosome 15 in Cucurbita maxima. Sequence analysis in the promoter region of CmaCh15G005270. Red box represented two insertions and two deletions in −1300~−800 bp of the promoter region in Wuminglv. (B) Recombinants analysis for map-based cloning of WSC. White and black columns represented homozygous fragments identity with Wuminglv background (white columns) and Agol background (black columns), respectively. Gray columns represented heterozygous fragments. The numbers in brackets indicated the accession number of the recombined plants derived from F2 segregating population.
Figure 3. Map-based cloning of WSC. (A) Map-based cloning of WSC to 60.4 kb on Chromosome 15 in Cucurbita maxima. Sequence analysis in the promoter region of CmaCh15G005270. Red box represented two insertions and two deletions in −1300~−800 bp of the promoter region in Wuminglv. (B) Recombinants analysis for map-based cloning of WSC. White and black columns represented homozygous fragments identity with Wuminglv background (white columns) and Agol background (black columns), respectively. Gray columns represented heterozygous fragments. The numbers in brackets indicated the accession number of the recombined plants derived from F2 segregating population.
Ijms 22 02972 g003
Ijms 22 02972 g004 550
Figure 4. Transcript expression of CmaCh15G005270 in root, stem, leaf, flower, pulp, and different development stage of seed coat in Wuminglv and Agol. CmaCh15G005270 was highly expressed in leaf, flowers, pulp, and seed coat in Wuminglv.
Figure 4. Transcript expression of CmaCh15G005270 in root, stem, leaf, flower, pulp, and different development stage of seed coat in Wuminglv and Agol. CmaCh15G005270 was highly expressed in leaf, flowers, pulp, and seed coat in Wuminglv.
Ijms 22 02972 g004
Ijms 22 02972 g005 550
Figure 5. RNA-seq data analysis of phenylpropanoid/flavonoid metabolism pathway. (A) The amount of DEGs in different stages of seed coat development. (B) Heatmap of DEGs in phenylpropanoid/flavonoid metabolism pathway at 0~42 DPA seed coat. (C) GO enrichment analysis. (D) KEGG enrichment analysis. (E) Transcript expression levels in nine candidate genes in different stages of seed coat development in Wuminglv and Agol.
Figure 5. RNA-seq data analysis of phenylpropanoid/flavonoid metabolism pathway. (A) The amount of DEGs in different stages of seed coat development. (B) Heatmap of DEGs in phenylpropanoid/flavonoid metabolism pathway at 0~42 DPA seed coat. (C) GO enrichment analysis. (D) KEGG enrichment analysis. (E) Transcript expression levels in nine candidate genes in different stages of seed coat development in Wuminglv and Agol.
Ijms 22 02972 g005
Ijms 22 02972 g006 550
Figure 6. Expression patterns of twelve genes involved in phenylpropanoid/flavonoid metabolism pathway in Wuminglv and Agol seed coat.
Figure 6. Expression patterns of twelve genes involved in phenylpropanoid/flavonoid metabolism pathway in Wuminglv and Agol seed coat.
Ijms 22 02972 g006
Table
Table 1. Seed coat color values of parents and segregating populations obtained by chromometer.
Table 1. Seed coat color values of parents and segregating populations obtained by chromometer.
L*a*b*E
P1 (yellow)68.89 ± 3.4012.45 ± 6.0225.14 ± 6.6875.32 ± 3.95
P2 (white)86.71 ± 6.080.98 ± 0.964.44 ± 3.0586.89 ± 5.90
F1 (P1 × P2) (yellow)66.13 ± 4.6911.46 ± 2.1219.74 ± 5.4170.82 ± 5.03
F1′ (P1 × P2) (yellow)67.03 ± 4.0811.13 ± 2.5819.54 ± 5.0470.83 ± 5.41
F2 (separation)59.63~94.710.04~17.281.67~30.4664.53~94.80
BC1 (yellow)68.71 ± 6.6211.94 ± 4.4021.43 ± 5.2472.41 ± 6.23
BC2 (separation)60.89~96.270.01~15.050.87~21.8364.17~96.27
P3 (yellow)65.35 ± 4.967.59 ± 3.6619.39 ± 3.4668.89 ± 5.28
P4 (white)82.67 ± 3.073.57 ± 1.5815.13 ± 4.5383.79 ± 2.62
F1 (P3 × P4) (yellow)64.66 ± 4.2913.61 ± 2.4521.95 ± 2.5969.58 ± 4.02
F1’ (P3 × P4) (yellow)64.28 ± 5.1113.09 ± 2.4423.82 ± 2.4269.48 ± 4.81
E = [(L*)2 + (a*)2 + (b*)2]1/2.
Table
Table 2. The segregation proportion of yellow seed coat and white seed coat in P1, P2, F1, F1′, F2, BC1, and BC2.
Table 2. The segregation proportion of yellow seed coat and white seed coat in P1, P2, F1, F1′, F2, BC1, and BC2.
PopulationNumber of PlantsExpected RatioChi Square Test
Yellow Seed CoatWhite Seed CoatTotal
P1505
P2066
F113013
F1’13013
F2282953773:1χc2 = 0.002 < 0.052 = 3.84
BC11090109
BC265561211:1χc2 = 0.17 < 0.052 = 3.84
Table
Table 3. Annotation and homologs of candidate genes in Arabidopsis thaliana.
Table 3. Annotation and homologs of candidate genes in Arabidopsis thaliana.
Gene IDLocation in Chr.15Gene Function AnnotationHomologs in AtGene Description in TAIR10
CmaCh15G0052202407955~2413296Vesicle-associated membrane protein 727AT3G54300Synaptobrevin-like protein family
CmaCh15G0052302413306~2414285Unknown protein––
CmaCh15G0052402415047~2420995AMMECR1 familyAT2G38710AMMECR1 family
CmaCh15G0052502421828~2424505Multiple myeloma tumor-associated protein, putativeAT3G52220Leukocyte immunoglobulin-like receptor family protein
CmaCh15G0052602424963~2427834F1F0-ATPase inhibitor proteinAT5G04750F1F0-ATPase inhibitor protein
CmaCh15G0052702428976~2431637MYB transcription factorAT5G04760R-R-type MYB protein which plays negative roles in salt stress and is required for ABA signaling in Arabidopsis
CmaCh15G0052802441111~2445124Cationic amino acid transporter 7AT5G04770Cationic amino acid transporter (CAT) subfamily of amino acid polyamine choline transporters
CmaCh15G0052902449424~2462027Pentatricopeptide repeat-containing proteinAT5G04780Pentatricopeptide repeat (PPR) superfamily protein
CmaCh15G0053002463445~2466291Double Clp-N motif-containing P-loop nucleoside triphosphate hydrolases superfamily protein, putativeAT3G52490SMAX1 and SMAX1-like protein that has weak similarity to AtHSP101
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, Y.; Zhang, M.; Shu, Q.; Ma, W.; Sun, T.; Xiang, C.; Wang, C.; Duan, Y. Genetic Mapping and Identification of the Candidate Gene for White Seed Coat in Cucurbita maxima. Int. J. Mol. Sci. 2021, 22, 2972. https://doi.org/10.3390/ijms22062972

AMA Style

Shi Y, Zhang M, Shu Q, Ma W, Sun T, Xiang C, Wang C, Duan Y. Genetic Mapping and Identification of the Candidate Gene for White Seed Coat in Cucurbita maxima. International Journal of Molecular Sciences. 2021; 22(6):2972. https://doi.org/10.3390/ijms22062972

Chicago/Turabian Style

Shi, Yuzi, Meng Zhang, Qin Shu, Wei Ma, Tingzhen Sun, Chenggang Xiang, Changlin Wang, and Ying Duan. 2021. "Genetic Mapping and Identification of the Candidate Gene for White Seed Coat in Cucurbita maxima" International Journal of Molecular Sciences 22, no. 6: 2972. https://doi.org/10.3390/ijms22062972

APA Style

Shi, Y., Zhang, M., Shu, Q., Ma, W., Sun, T., Xiang, C., Wang, C., & Duan, Y. (2021). Genetic Mapping and Identification of the Candidate Gene for White Seed Coat in Cucurbita maxima. International Journal of Molecular Sciences, 22(6), 2972. https://doi.org/10.3390/ijms22062972

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop