1. Introduction
Sweet potato (Ipomoea batatas (L.) Lam) is a hexaploid (2n = 6× = 90) species, rich in high calories, proteins, vitamins, and minerals. It is the seventh most important crop in the world and the fourth most significant crop cultivated in China [1]. Sweet potato breeding programs worldwide are integrating sensory characteristics to provide global markets with improved cultivars [2,3,4]. SRFC is identified as a key driver of consumers’ sensory and hedonic expectations before tasting [5].
SRFC fluctuates in thousands of ways, white, light yellow, yellow, orange, light red, pink, red, fuchsia, purple, deep purple, and dark purple viewed by human eyes. According to SRFC, sweet potato can be divided into white, yellow or orange, and purple cultivars. According to chemical composition, sweet potato can be classified into starch, anthocyanin, and carotenoid cultivars. The major pigments in sweet potato are anthocyanins and carotenoids, which are responsible for sweet potato’s storage root colors [6]. There is a correlation between color and substance content of storage root flesh: the purple-fleshed has the highest anthocyanins content, and orange-fleshed has the highest carotene, while white-fleshed have the lowest [7,8]. Anthocyanin and carotene contents varied widely among the sweet potato varieties with different flesh color.
Anthocyanins are responsible for the shiny orange, pink, red, purple, and blue colors in the flowers and fruits of some plants [9]. Anthocyanins are initially synthesized in the cytoplasm, and transported to vacuoles, where they are sequestered to form colors in different plant organs [10]. Vesicles, glutathione S-transferase, and membrane transporters are involved in the transport of anthocyanins [11]. Anthocyanin biosynthesis is a branch of flavonoid synthesis, and its biosynthetic genes are regulated by the MBW complex comprising myeloblastosis (MYB), basic helix-loop-helix (bHLH), and the WD40 repeat-containing (WD40) protein [12]. IbMYB1, IbMYB340 and ethylene-responsive transcription factor 1 (IbERF1) have been reported to regulate the anthocyanins synthesis of storage root in sweet potato [13,14,15].
Carotenoids are responsible for yellow, orange, and red colors found in nature, and carotenoid synthesis is initiated by geranylgeranyl diphosphate (GGPP) synthesis in the cytosolic mevalonic acid (MVA) pathway and methylerythritol 4-phosphate (MEP) pathway [16]. Plant carotenoid biosynthesis is mostly derived from the MEP route, which is only found in plastids [16]. The diverse carotenoid accumulations, such as α-carotene, β-carotene, lutein and lycopene, are the primary contributors to a variety of carrot root colors, including orange, yellow and red [17,18].The β-carotene, mutatochrome, phytoene, and lutein are the predominant carotenoid components found from high to low in sweet potato storage roots [19]. Carotenoid biosynthetic genes such as geranylgeranyl pyrophosphate synthase (GGPS), phytoene synthase (PSY), phytoene desaturase (PDS), ζ-carotene desaturase (ZDS), carotenoid isomerase (CRTISO), lycopene ε-cyclase (LCY-ε), lycopene β-cyclase (LCY-β), β-carotene hydroxylase (CHY-β), zeaxanthin epoxidase (ZEP), and carotenoid cleavage dioxygenase (CCD) have so far been isolated and characterized in sweet potato [20]. In addition to the carotenoid pathway genes, Orange (Or) gene can also enhance the accumulation of carotenoids in sweet potato [21].
Sweet potato is a hexaploid heterozygote with a complex genetic background, self-pollination infertility and cross incompatibility, which make genetic linkage analysis quite difficult. GWAS has become an important and effective method for crop breeding programs. Risch first proposed GWAS for the genetic studies of complex human disease [22]. GWAS uses high-density SNPs in the mapping group to filter the molecular markers associated with complex trait performance variation based on linkage disequilibrium (LD), then analyzes its genetic effect on phenotype [23]. SNP is widely used for genetic linkage construction of molecular breeding, greatly improving the resolution and precision of the genetic map with the rapid development of next generation sequencing (NGS) techniques [24,25]. At present, GWAS has become a routine strategy for decoding genotype–phenotype associations in many species, more than 1000 studies over the last decade have used GWAS for genetic association [26].
Previously, we used specific-locus amplified fragment sequencing (SLAF-seq) technology to sequence 300 sweet potato accessions [27,28]. In this study, we analyzed the raw data of the 300 sweet potato accessions with the genome of the modern cultivar of sweet potato, ‘Taizhong6’, as the reference genome. Finally, we developed a large number of high-quality SNPs and integrated GWAS and transcriptomic data to identify the main genes driving SRFC differences in sweet potato.
2. Materials and Methods
2.1. Selection of Plant Material and Storage Root Phenotyping
A set of 300 sweet potato accessions, including 76 landraces and 224 modern cultivars, were evaluated in the present study. We recorded the color of the storage root flesh by eye according to the identification standards which were shown in Table 1; SRFC was recorded (Table S1). These accessions were generated from different agro-climatic zones, and were cultivated in 2019 in the Experimental Farm of Hubei Academy of Agricultural Sciences in Ezhou District, China.
2.2. Read Alignment and Variation Calling
Initially, the SLAF-sequencing reads of 300 sweet potato accessions were mapped to the sweet potato reference genome (
2.3. GWAS Analysis
To increase the number of markers that can be tested for association with a trait, we used Beagle to perform genotype imputation [33]. GWAS was performed using the Essicient Mixed-Model Association expedited (EMMAX) software package (
2.4. RNA-Seq Analysis
The transcriptomic data of root ‘Xiangshu99′ (white-fleshed) and ‘Zhezi No1 (purple-fleshed) was found by downloading the raw sequence data from the NCBI Short Read Archive (No. PRJNA721067) [36]; the transcriptomic data of root in ‘Beniharuka (BH)’ (yellow-fleshed) cultivar and its ‘white-fleshed mutants (WH)’ was found by downloading the raw sequence data from the NCBI Short Read Archive (No. PRJDB10052) [37]. Raw data were first processed with Fastpsoftware (
2.5. qRT-PCR
qRT-PCR analyses were carried out to determine the reliability of the RNA-seq results for expression profile analysis. All primers were designed according to the mRNA sequences and were synthesized commercially in the company Tianyi Huiyuan, Wuhan (Table S2). qRT-PCR for mRNAs were carried out in a 10 μL system: 1.0 μL cDNA product, 5 μL 2× qPCR mix, 3.2 μL nucleotide-free water and 0.4 μL 2.5 μM for each of the forward and reverse primers. The reactions were incubated in a BIO-RAD CFX96 for 30 s at 95 °C, followed by 40 cycles of 5 s at 95 °C, then 30 s at 60 °C. Three replications were used for all reactions and β-actin served as the endogenous reference gene. The 2-ΔΔCT method was employed to analyze the relative fold changes of genes [41].
2.6. Statistical Analysis
A total of 567,828 SNPs across 15 linkage groups (LG) remained after quality filtering. The 567,828 SNPs were subsequently used for further analysis. ADMIXTURE was employed to investigate population structure based on the maximum-likelihood method and generated a Q matrix [42]. A phylogenetic tree was constructed using MEGA7 based on the neighbor-joining method, visually edited using FigTree software (
3. Results
3.1. SLAF-Sequencing of 300 Sweet Potato Accessions
In our previous study, we constructed SLAF-seq library of 300 sweet potato accessions using SLAF-seq technology and obtained a total of 498.14 Mb reads [44]. In this study, we mapped the raw data to the sweet potato reference genome with BWA, the rate ranged from 6.46% to 95.11%, with an average rate of 92.61% (Table S3). After filtering, we identified 567,828 high confidence SNPs (Table S4) (missing data < 50%, minor allele frequency (MAF) > 5%) which were used for subsequent analyses. The SNPs ranged from 29,393 (LG10) to 50,397 (LG2), and were distributed across 15 LGs of sweet potato (Figure 1), the SNP number of LG1 to LG5 was 34,856, 50,379, 30,708, 36,340, 36,129, 36,915, 43,424, 30,858, 38,001, 29,393, 47,309, 38,169, 38,107, 391,69, 38,053 in turn. We examined the distribution of polymorphisms in sweet potato genomic regions and found that 109,738, 153,623, 112,695, 99,800, and 100,673 SNPs were located within exons, introns, upstream regions (within 5 kb upstream of the transcription start sites), downstream regions (within 5 kb downstream of transcription stop sites), and intergenic regions, respectively.
3.2. Population Structure and Genomic Variation among the 300 Sweet Potato Accessions
The 300 sweet potato accessions included 76 landraces and 224 modern cultivars, which corresponded to roughly their geographical distributions in Asia and North America, respectively (Figure 2a). Neighbor-joining cluster analysis showed that the 300 accessions were divided into seven groups (Figure 2b); the seven groups contained 63, 23, 24, 21, 66, 56, 47 accessions, respectively. The population structure analysis and principal component analysis (PCA) found that the 300 sweet potato accessions were clustered into 4 groups (K = 4), namely groups 1, 2, 3, and 4 (Figure 2c). Group 1, 2, 3 and 4 contain 69, 71, 83, and 77 accessions, respectively. PCA clustering showed that the accessions of the core set were distributed into two clusters: landraces and modern cultivars. Principal component 1 (PC1) represents 15.51% of the total variation, while PC2 represents 12.0% of the total variation (Figure 2d). LD analysis showed that the delay rate of landraces was faster than that of modern cultivars (Figure 2e), suggesting a high frequency of genetic recombination in landraces, further supporting the richer genetic diversity of landraces.
3.3. GWAS for Identifying Significant SNPs Associated with SRFC
SRFC is one of the most important traits associated with the quality of sweet potato. There is a correlation between the color and chemical composition of sweet potato flesh. We conducted GWAS to identify significant SNP associated with SRFC by utilizing the 567,828 SNP markers (Table S5). GWAS results were visually examined by using Manhattan plots (Figure 3a) and quantile-quantile plots (Figure 3b). The line shows the threshold for genome-wide significance (−log10P > 7). A total of 5 significant SNPs associated with SRFC were identified (Table 2).
Interestingly, g55964, F3′H, was identified within the 4460,608-4466,687 bp in LG14, where 1 SNP (LG14: 467,066) was strongly associated with SRFC. GWAS analysis revealed that the SNP was located at ~379 bp upstream of the start codon of the F3′H which is involved in anthocyanins biosynthesis.
g17506, BAG2-like, was identified within the 5866,570-5867,826 bp in LG5, where one SNP (LG5: 5,867,009) was strongly associated with SRFC. GWAS analysis discovered that the SNP was located in the intron region of the BAG2-like.
g25206, UGT-73D1-like, was identified within the 31,266,628-31,269,316 bp in LG6, where two SNPs (LG6: 31,267,101 and LG6: 31,267,107) were strongly associated with SRFC. GWAS analysis found that the two SNPs were located in the UTR3 of the UGT-73A-like.
g58377, VQ25 isoform X2, was identified within the 21,793,980-21,794,619 bp in LG14, where one SNP (LG14: 21,791,164) was strongly associated with SRFC. GWAS analysis showed that the SNP was located at ~2816 bp upstream of VQ motif-containing protein 25 isoform X2.
3.4. Differentially Expressed Gene between White-Fleshed and Purple-Fleshed Sweet Potato
To study the key genes that regulate purple formation of root in sweet potato, we analyzed the transcriptomic data of root in ‘Xiangshu99′ (white-fleshed) and ‘Zhezi No1 (purple-fleshed). In the transcriptomic analysis, the expression levels of the genes were shown in Table S6a. Compared with ‘Xiangshu99′, 442 upregulated and 117 downregulated DEGs in ‘Zhezi No1’were detected (Table S6b). GO enrichment analysis was conducted to further understand the function of the identified DEGs. A total of 139 (Table S7a) and 64 (Table S7b) terms were significantly enriched in the upregulated and downregulated genes, with p < 0.05, respectively. The 20 most enriched GO terms in the upregulated and downregulated genes are shown in Figure 4. The “flavonoid biosynthetic process” (GO: 0009813) and “flavonoid metabolic process” (GO: 0009812) were the most abundant upregulated DEGs (Figure 4a). The predominantly enriched pathways in the downregulated DEGs were the “Phenylpropanoid biosynthetic process” (GO: 0009699) and “Suberin biosynthetic process” (GO: 0010345) (Figure 4b). The KEGG enrichment analysis found that a total of five (Table S8a) and five (Table S8b) terms were significantly enriched in the upregulated and downregulated genes, with p < 0.05, respectively. The “Flavonoid biosynthesis” (ko00941), and “Anthocyanin biosynthesis” (ko00942) were the most abundant in upregulated DEGs (Figure 5a). The predominantly enriched pathways in the downregulated DEGs were involved in the “Protein processing in endoplasmic reticulum” (ko04141), whereas none of the metabolic pathways were significantly enriched (Figure 5b). The results of GO and KEGG enrichment implied that the flavonoid anthocyanin biosynthetic-related genes regulate the purple formation of roots in sweet potato.
3.5. Differentially Expressed Gene between White Fleshed and Yellow-Fleshed Sweet Potato
To identify the key genes involved in yellow-fleshed, we analyzed the transcriptomic data of the root in ‘Beniharuka (BH)’ (yellow-fleshed) cultivar and its ‘white-fleshed mutants (WH)’. In the transcriptomic analysis, the expression levels of the genes were shown in Table S9a. Compared with ‘WH’, a total of 108 upregulated and 97 downregulated DEGs were detected in ‘Beniharuka’ (Table S9b). GO enrichment analysis was conducted to further understand the function of these DEGs, and a total of 75 (Table S10a) and 56 (Table S10b) terms were significantly enriched in the upregulated and downregulated genes, with a p < 0.05, respectively. The 20 most enriched GO terms in the upregulated and downregulated genes were shown in Figure 6. Among the GO terms in upregulated genes, “abscisic acid transport” (GO: 0080168), “isoprenoid transport” (GO: 0046864), and “terpenoid transporter” (GO: 0046865) were significantly enriched in upregulated genes (Figure 6a). The predominantly enriched pathways in the downregulated DEGs were reported for “cadmium ion transport” (GO: 001569) (Figure 6b). The KEGG enrichment analysis found that a total of 16 (Table S11a) and one (Table S11b) terms were significantly enriched in the upregulated and downregulated genes, with p < 0.05, respectively. The most abundant upregulated DEGs were observed for “Glycolysis/Gluconeogenesis” (ko00010) (Figure 7a). The enriched pathways in the downregulated DEGs were involved in “Isoquinoline alkaloid biosynthesis” (ko00950) (Figure 7b).The results of GO enrichment implied that the terpenes transporter-related genes regulate the yellow formation of sweet potato root.
3.6. Exploring Candidate Genes for SRFC by RNA-Seq
According to the result of RNA-seq, a total of 35 GO terms were found to be associated with the two candidate genes, and four KEGG pathways were found to be associated with two candidate genes (Table 3 and Table 4). GO analysis showed that g55964 belonged to “secondary metabolic process”, “secondary metabolite biosynthetic process”, “metabolic process”, and “biosynthetic process”, whereas g25206 belonged to “glucosyltransferase activity” and “UDP-glycosyltransferase activity”. KEGG analysis suggested that g55964 and g25206 were involved in the “biosynthesis of secondary metabolites”, signifying the involvement of g55964 and g25206 in the biosynthesis of metabolites in sweet potato.
3.7. Candidate Genes Expression Analysis
We quantified the expression of four candidate genes in the flesh of sweet potato varieties that varied in color by analyzing the transcriptomic data (Figure 6, Table 5). The result displayed that only g55964 displayed significant differential expression between ‘Zhezi No1′ and ‘Xiangshu 99′. We found that the expression of g55964 in ‘Zhezi No1′ was 46.8-fold higher than that of ‘Xiangshu99′ (Table S6b), and all candidate genes had no expression in ‘WH’ and ‘BH’ (Figure 8). We performed qRT-PCR to quantify the relative expression of g55964 among eight varieties. The result of qRT-PCR showed that the expression of g55964 in purple-fleshed sweet potato was significantly (p < 0.01) higher than that of non-purple fleshed sweet potato (Figure 9, Table S12).
4. Discussion
4.1. Genetic Variation among the 300 Accessions
In our previous study, we sequenced and characterized the genetic diversities of 300 natural accessions of sweet potato by SLAF-seq, most of which were from major sweet potato production regions (i.e., America, Africa, China, Japan, and Korea) [28]. In this study, population-genetics analyses of the SNPs were generated by mapping reads to hexaploid sweet potato reference genomes, the ‘Taizhong6′. Based on population structure and PCA, 300 sweet potato accessions were clustered into four groups. The faster LD decay rate of landraces was faster than that of modern cultivars, suggesting that the genetic diversity of landraces was higher than modern cultivars. The result was consistent with He et al. [45], who reported that sweet potato landraces possess high genetic diversity and that the genetic diversity of modern cultivars is low.
4.2. GWAS Analysis of SRFC in Sweet Potato
In order to mine the candidate genes for SRFC, we recorded the SRFC of 300 sweet potato accessions, and we used the SNPs site to conduct GWAS of SRFC. We identified five potential SNPs which were associated with SRFC (−log10P > 7), and linked to four candidate genes. g55964 was annotated as F3′H, catalyzing naringenin and dihydokaempferol, and generated eriodictyol and dihydroquercetin which are important intermediates in the biosynthesis of anthocyanins and procyanidins [46]. We found that g25206 was annotated as UGT-73D1-like. It has been reported that uridine diphosphate glycosyltransferase (UGT) regulates the bioactivity, solubility, and transportation of receptors within the cells through modulating glycosyl transfer reaction [47]. UGT can catalyze flavonols and anthocyanins into corresponding glycosides [48,49]. g17506 was annotated as BAG2-like. Previously, it was reported that BCL2-associated athanogene (BAG) family members can regulate plant stress and development [50]. To the best of our knowledge, the BAG family has not been reported to be involved in metabolite synthesis. The g58377 was annotated as VQ25-isoform X2. The VQ motif-containing proteins (VQ) can interact with the WRKY transcription factor through the conserved VQ domain and plays an important role in the signal pathway mediated by the WRKY transcription factor [51,52]. The WRKY transcription factor can control the secondary metabolites synthesis, such as phenols, terpenes, and alkaloids [53]. Among them, anthocyanins and carotenoids are the two main pigments in sweet potato; it is not clear whether VQ can interact with WRKY transcription factors to regulate the synthesis of anthocyanins and carotenoids. Whether BAG and VQ are involved in metabolite synthesis remained to be studied.
4.3. Functional Analysis of DEGs in Different SRFC Varieties
The GO and KEGG enrichment result showed that flavonoid and anthocyanin biosynthetic-related genes regulated the purple formation of the root in sweet potato. Li et al. also reported that flavonoid and anthocyanin biosynthesis were related to the flesh color of storage root [54]. The g55694 was annotated into the flavonoid biosynthesis pathway, and its expression in purple-fleshed sweet potato was much higher than that of non-purple fleshed sweet potato. Thus we speculated that g55694 was the one of the key genes involved in the purple formation of the root in sweet potato.
4.4. Quantitative Trait Loci Analysis of Flesh Color in Sweet Potato
Previously, some studies have reported that QTLs significantly associated with flesh color in sweet potatoes are located on LG3 and LG12. Zhang et al. [55] performed GWAS for flesh color in purple sweet potato and identified one QTL for IbMYB1-2 on LG12. Gemenet et al. [56] identified two major QTLs that affect β-carotene and flesh color; they are phytoene synthase and Orange gene located within QTLs on LG3 and LG12, respectively. Yan et al. [57] identified ten QTLs for flesh color, skin color, and anthocyanin content on LG12 based on the linkages map. In this study, we identified one QTL for flesh color on LG14, the QTL was linked to F3′H which was involved in flavonoid biosynthesis. Anthocyanins and carotenes are responsible for sweet potato storage root colors, with complex components, which are affected by various environmental factors. Thus, further research is required to dissect the genomic mechanisms controlling different types of anthocyanins and carotenes.
5. Conclusions
We recorded SRFC of 300 sweet potato accessions and conducted GWAS in sweet potato for SRFC. The study identified five unique SNPs which were significantly associated with SRFC trait, and linked to four genes. The transcriptome analyses showed that anthocyanin is responsible for the formation of purple roots in sweet potato, while carotenoid is responsible for the formation of yellow roots. We also employed GWAS and transcriptomic analyses to identify one QTL for flesh color on LG14, and the QTL was linked to g55964 which was involved in flavonoid biosynthesis.
C.J. and X.Y. designed the study. Y.L. and R.P. analyzed the data and prepared the manuscript. W.Z., L.W., J.L., S.C. and X.J. provided manuscript revision services. All authors have read and agreed to the published version of the manuscript.
This work was supported by the National Key R&D Program of China (2018YFD1000700, 2018YFD1000705); China Agriculture Research System (CARS-11-C-15); the Characteristic Discipline of Hubei Academy of Agricultural Sciences (2015TSXK06).
Not applicable.
Not applicable.
Not applicable.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 2. Distribution, population structure, PCA, and LD decay of the 300 sweet potato germplasm accessions. (a) The geographic distribution of the 300 sweet potato accessions based on country of origin. (b)The phylogenetic tree of 300 sweet potato accessions based on the analysis of 567,828 SNPs. (c) Population structure analysis using 567,828 SNPs (missing data < 50%, MAF > 5%, r2 < 0.2), K = 4. (d) PCA plot of the first two components (PC1 and PC2) of the 300 sweet potato accessions. (e) Genome-wide average LD decay of the 300 sweet potato accessions.
Enlarge this image.Figure 3. Genome-wide association study analysis of SRFC in 300 sweet potato accessions. (a) Manhattan plots of SNPs association for SRFC. (b) Quantile–quantile plots of expected and observed − log10 (p) values at an FDR cut-off ≤ 0.05. The gray areas are the 95% confidence intervals under the null hypothesis of no association between the SNP and the trait. The red line is the expected line under the null distribution. The blue points are the observed distribution.
Enlarge this image.Figure 4. GO enrichment of DEGs in the root of Xiangshu99 vs. Zhezi No1. (a) upregulated genes; (b) downregulated genes. GO terms belonging to biological processes (BP), cellular components (CC), and molecular functions (MF).
Enlarge this image.Figure 4. GO enrichment of DEGs in the root of Xiangshu99 vs. Zhezi No1. (a) upregulated genes; (b) downregulated genes. GO terms belonging to biological processes (BP), cellular components (CC), and molecular functions (MF).
Enlarge this image.Figure 5. KEGG enrichment of DEGs in the root of Xiangshu99 vs. Zhezi No1. (a) upregulated genes; (b) downregulated genes.
Enlarge this image.Figure 5. KEGG enrichment of DEGs in the root of Xiangshu99 vs. Zhezi No1. (a) upregulated genes; (b) downregulated genes.
Enlarge this image.Figure 6. GO enrichment of DEGs in the root of WH3 vs. BH. (a) upregulated genes; (b) downregulated genes. GO terms belonging to biological processes (BP), cellular components (CC), and molec-ular functions (MF).
Enlarge this image.Figure 7. KEGG enrichment of DEGs in the root of WH3 vs. BH. (a) upregulated genes; (b) downregulated genes.
Enlarge this image.Figure 8. Heatmap of candidate genes in storage root by RNA-seq. (a) The expression of candidate genes in ‘Xiangshu99′ and ‘Zhezi No1′ (No. PRJNA721067); (b) The expression of candidate genes in ’BH’ and ‘WH’ (No. PRJDB10052).
Enlarge this image.Figure 9. g55964 expression comparisons between different colors of flesh by RT-qPCR. The SRFC in ‘Xushu18′and ‘Wansu31′are white, the SRFC in ‘Qinshu No.5′and ‘Eshu11′are light-yellow, the SRFC in ‘Xiangshu99′, ‘Guangshu87′, and ‘Chuanshu294′ are jacinth, the SRFC in ‘Zhezishu No.1′, ‘Ganshu No.1′ and ‘Xuzishu No.8′ are purple.
Morphological traits and their criteria for sweet potato.
| Morphological Traits | Criteria for Recording |
|---|---|
| Storage root flesh color, SRFC | 1 = White, 2 = Light yellow, 3 = Yellow, 4 = Orange, 5 = Jacinth ikaika 6 = Light red, 7 = Red, 8 = Fuchsia, 9 = Purple, 10 = Dark purple, 11 = Black purple |
Summary of SNPs above the significance thresholds.
| Marker Name | Candidate Genes | Annotated | p-Value | R2 | LG | Allele | Position |
|---|---|---|---|---|---|---|---|
| LG14: 4,467,066 | g55964 | F3′H | 3.63 × 10−8 | 0.096366 | LG14 | A/G | 4,467,066 |
| LG5: 5,867,009 | g17506 | BAG2-like | 6.31 × 10−8 | 0.093125 | LG5 | A/G | 5,867,009 |
| LG6: 31,267,101 | g25206 | UGT-73D1-like | 6.41 × 10−8 | 0.093027 | LG6 | C/G | 31,267,101 |
| LG6: 31,267,107 | g25206 | UGT-73D1-like | 7.99 × 10−8 | 0.091735 | LG6 | T/C | 31,267,107 |
| LG14: 21,791,164 | g58377 | VQ25-isoform X2 | 9.43 × 10−8 | 0.090762 | LG14 | C/T | 21,791,164 |
GO analysis of candidate genes related to the RNA-seq results.
| Gene ID | GO |
|---|---|
| g55964 | GO:0016705; GO:0016020; GO:0055114; GO:0016709; GO:0004497; GO:0009733; GO:0042221; GO:0019748; GO:0016491; GO:0050896; GO:0044550; GO:0008150; GO:0008152; GO:0009725; GO:0010033; GO:0009719; GO:0003674; GO:0005575; GO:0009058 |
| g17506 | none |
| g25206 | GO:0046527; GO:0005622; GO:0080044; GO:0043226; GO:0043231; GO:0043229; GO:0044424; GO:0035251; GO:0005623; GO:0016740; GO:0016758; GO:0016757; GO:0044464; GO:0080043; GO:0003674; GO:0043227; GO:0005575; GO:0008194 |
| g58377 | none |
KEGG analysis of candidate genes related to the RNA-seq results.
| Pathway_ID | Pathway_Name | Gene_ID | KO_Entry | EC |
|---|---|---|---|---|
| ko00941 | Flavonoid biosynthesis | g55964 | K05280 | EC:1.14.14.82 |
| ko00944 | Flavone and flavonol biosynthesis | g55964 | K05280 | EC:1.14.14.82 |
| ko01100 | Metabolic pathways | g55964 | K05280 | EC:1.14.14.82 |
| ko01110 | Biosynthesis of secondary metabolites | g55967 | K05280 | EC:1.14.14.82 |
| ko01110 | Biosynthesis of secondary metabolites | g25206 | K13496 | EC:2.4.1. |
The FPKM value of genes in varieties.
| Xiangshu99 | Zhezi No1 | BH | WH | |
|---|---|---|---|---|
| g55964 | 0.46 | 16.86 | 0.00 | 0.00 |
| g17506 | 11.17 | 15.35 | 0.00 | 0.00 |
| g25206 | 0.11 | 0.00 | 0.00 | 0.00 |
| g58377 | 0.12 | 0.00 | 0.00 | 0.00 |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Yang, J.; Moeinzadeh, M.H.; Kuhl, H.; Helmuth, J.; Xiao, P.; Haas, S.; Liu, G.; Zheng, J.; Sun, Z.; Fan, W. et al. Haplotype-resolved sweet potato genome traces back its hexaploidization history. Nat. Plants.; 2017; 3, pp. 696-703. [DOI: https://dx.doi.org/10.1038/s41477-017-0002-z]
2. Leighton, C.S.; Schönfeldt, H.; Kruger, R. Quantitative descriptive sensory analysis of five different cultivars of sweet potato to determine sensory and textural profiles. J. Sens. Stud.; 2010; 25, pp. 2-18. [DOI: https://dx.doi.org/10.1111/j.1745-459X.2008.00188.x]
3. Tomlins, K.; Rwiza, E.; Nyango, A.; Amour, R.; Ngendello, T.; Kapinga, R.; Rees, D.; Jolliffe, F. The use of sensory evaluation and consumer preference for the selection of sweetpotato cultivars in East Africa. J. Sci. Food Agric.; 2004; 84, pp. 791-799. [DOI: https://dx.doi.org/10.1002/jsfa.1712]
4. Leksrisompong, P.P.; Whitson, M.E.; Truong, V.D.; Drake, M.A. Sensory attributes and consumer acceptance of sweet potato cultivars with varying flesh colors. J. Sens. Stud.; 2012; 27, pp. 59-69. [DOI: https://dx.doi.org/10.1111/j.1745-459X.2011.00367.x]
5. Lado, J.; Vicente, E.; Moltini, A.I.; Alcaire, F.; Ares, G. Integrating consumer perception in the selection of purple skin sweet potato cultivars. J. Sci. Food Agric.; 2021; 101, pp. 3620-3629. [DOI: https://dx.doi.org/10.1002/jsfa.10991] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33275810]
6. Yang, Y.; Shi, D.; Wang, Y.; Zhang, L.; Chen, X.; Yang, X.; Xiong, H.; Bhattarai, G.; Ravelombola, W.; Olaoye, D. et al. Transcript profiling for regulation of sweet potato skin color in Sushu8 and its mutant Zhengshu20. Plant Physiol. Biochem.; 2020; 148, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.plaphy.2019.12.035] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31923733]
7. Teow, C.C.; Truong, V.D.; McFeeters, R.F.; Thompson, R.L.; Pecota, K.V.; Yencho, G.C. Antioxidant activities, phenolic and β-carotene contents of sweet potato genotypes with varying flesh colours. Food Chem.; 2007; 103, pp. 829-838. [DOI: https://dx.doi.org/10.1016/j.foodchem.2006.09.033]
8. Azeem, M.; Mu, T.H.; Zhang, M. Effects of high hydrostatic pressure and soaking solution on proximate composition, polyphenols, anthocyanins, β-carotene, and antioxidant activity of white, orange, and purple fleshed sweet potato flour. Food Sci. Technol. Int.; 2020; 26, pp. 388-402. [DOI: https://dx.doi.org/10.1177/1082013219892716]
9. Castañeda-Ovando, A.; Pacheco-Hernandez, M.D.L.; Páez-Hernández, M.E.; Rodríguez, J.A.; Galán-Vidal, C.A. Chemical studies of anthocyanins: A review. Food Chem.; 2009; 113, pp. 859-871. [DOI: https://dx.doi.org/10.1016/j.foodchem.2008.09.001]
10. Zhao, J.; Dixon, R.A. The “ins” and “outs” of flavonoid transport. Trends Plant Sci.; 2010; 15, pp. 72-80. [DOI: https://dx.doi.org/10.1016/j.tplants.2009.11.006]
11. Zhao, J. Flavonoid transport mechanisms: How to go, and with whom. Trends Plant Sci.; 2015; 20, pp. 576-585. [DOI: https://dx.doi.org/10.1016/j.tplants.2015.06.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26205169]
12. Xu, W.; Dubos, C.; Lepiniec, L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci.; 2015; 20, pp. 176-185. [DOI: https://dx.doi.org/10.1016/j.tplants.2014.12.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25577424]
13. Wei, Z.Z.; Hu, K.D.; Zhao, D.L.; Tang, J.; Huang, Z.Q.; Jin, P.; Li, Y.H.; Han, Z.; Hu, L.Y.; Yao, G.F. et al. MYB44 competitively inhibits the formation of the MYB340-bHLH2-NAC56 complex to regulate anthocyanin biosynthesis in purple-fleshed sweet potato. BMC Plant Biol.; 2020; 20, 258. [DOI: https://dx.doi.org/10.1186/s12870-020-02451-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32503504]
14. Park, S.C.; Kim, Y.H.; Kim, S.H.; Jeong, Y.J.; Kim, C.Y.; Lee, J.S.; Bae, J.Y.; Ahn, M.J.; Jeong, J.C.; Lee, H.S. et al. Overexpression of the IbMYB1 gene in an orange-fleshed sweet potato cultivar produces a dual-pigmented transgenic sweet potato with improved antioxidant activity. Physiol. Plant.; 2015; 153, pp. 525-537. [DOI: https://dx.doi.org/10.1111/ppl.12281]
15. Ning, Z.; Hu, K.; Zhou, Z.; Zhao, D.; Tang, J.; Wang, H.; Li, L.; Ding, C.; Chen, X.; Yao, G. et al. IbERF71, with IbMYB340 and IbbHLH2, coregulates anthocyanin accumulation by binding to the IbANS1 promoter in purple-fleshed sweet potato (Ipomoea batatas L.). Plant Cell Rep.; 2021; 40, pp. 157-169. [DOI: https://dx.doi.org/10.1007/s00299-020-02621-0]
16. Rosas-Saavedra, C.; Stange, C. Biosynthesis of Carotenoids in Plants: Enzymes and Color. Subcell Biochem.; 2016; 79, pp. 35-69.
17. Wang, Y.H.; Li, T.; Zhang, R.R.; Khadr, A.; Tian, Y.S.; Xu, Z.S.; Xiong, A.S. Transcript profiling of genes involved in carotenoid biosynthesis among three carrot cultivars with various taproot colors. Protoplasma.; 2020; 257, pp. 949-963. [DOI: https://dx.doi.org/10.1007/s00709-020-01482-4]
18. Jourdan, M.; Gagné, S.; Dubois-Laurent, C.; Maghraoui, M.; Huet, S.; Suel, A.; Hamama, L.; Briard, M.; Peltier, D.; Geoffriau, E. Carotenoid content and root color of cultivated carrot: A candidate-gene association study using an original broad unstructured population. PLoS ONE.; 2015; 10, e0116674. [DOI: https://dx.doi.org/10.1371/journal.pone.0116674]
19. Drapal, M.; Rossel, G.; Heider, B.; Fraser, P.D. Metabolic diversity in sweet potato (Ipomoea batatas, Lam.) leaves and storage roots. Hortic. Res.; 2019; 6, 2. [DOI: https://dx.doi.org/10.1038/s41438-018-0075-5]
20. Kang, L.; Park, S.C.; Ji, C.Y.; Kim, H.S.; Lee, H.S.; Kwak, S.S. Metabolic engineering of carotenoids in transgenic sweetpotato. Breed. Sci.; 2017; 67, pp. 27-34. [DOI: https://dx.doi.org/10.1270/jsbbs.16118]
21. Park, S.C.; Kim, S.H.; Park, S.; Lee, H.U.; Lee, J.S.; Park, W.S.; Ahn, M.J.; Kim, Y.H.; Jeong, J.C.; Lee, H.S. et al. Enhanced accumulation of carotenoids in sweetpotato plants overexpressing IbOr-Ins gene in purple-fleshed sweetpotato cultivar. Plant Physiol. Biochem.; 2015; 86, pp. 82-90. [DOI: https://dx.doi.org/10.1016/j.plaphy.2014.11.017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25438140]
22. Khoury, M.J.; Yang, Q. The future of genetic studies of complex human diseases. Epidemiology.; 1998; 9, pp. 350-354. [DOI: https://dx.doi.org/10.1097/00001648-199805000-00023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9583430]
23. Zhongqu, D.; Jun, Z. Research progress of genome-wide association study. J. Zhejiang Univ.; 2015; 41, pp. 385-393.
24. Pang, Y.; Liu, C.; Wang, D.; St Amand, P.; Bernardo, A.; Li, W.; He, F.; Li, L.; Wang, L.; Yuan, X. et al. High-Resolution Genome-Wide Association Study Identifies Genomic Regions and Candidate Genes for Important Agronomic Traits in Wheat. Mol. Plant; 2020; 13, pp. 1311-1327. [DOI: https://dx.doi.org/10.1016/j.molp.2020.07.008]
25. Kump, K.L.; Bradbury, P.J.; Wisser, R.J.; Buckler, E.S.; Belcher, A.R.; Oropeza-Rosas, M.A.; Zwonitzer, J.C.; Kresovich, S.; McMullen, M.D.; Ware, D. et al. Genome-wide association study of quantitative resistance to southern leaf blight in the maize nested association mapping population. Nat. Genet.; 2011; 43, pp. 163-168. [DOI: https://dx.doi.org/10.1038/ng.747] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21217757]
26. Haijun, L.; Jianbing, Y. Crop genome-wide association study a harvest of biological relevance. Plant J.; 2019; 97, pp. 8-18.
27. Su, W.; Wang, L.; Lei, J.; Chai, S.; Liu, Y.; Yang, Y.; Yang, X.; Jiao, C. Genome-wide assessment of population structure and genetic diversity and development of a core germplasm set for sweet potato based on specific length amplified fragment (SLAF) sequencing. PLoS ONE.; 2017; 12, e0172066. [DOI: https://dx.doi.org/10.1371/journal.pone.0172066]
28. Su, W.J.; Zhao, N.; Lei, J.; Wang, L.J.; Chai, S.S.; Yang, X.S. SNP Sites Developed by Specific Length Amplification Fragment Sequencing (SLAF-seq) in Sweetpotato. Sci. Agric. Sin.; 2016; 49, pp. 27-34.
29. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics.; 2009; 25, pp. 1754-1760. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp324]
30. Schmidt, S. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res.; 2010; 20, pp. 1297-1303.
31. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res.; 2010; 38, e164. [DOI: https://dx.doi.org/10.1093/nar/gkq603] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20601685]
32. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; De Bakker, P.I.W.; Daly, M.J. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet.; 2007; 81, pp. 559-575. [DOI: https://dx.doi.org/10.1086/519795] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17701901]
33. Browning, B.L.; Zhou, Y.; Browning, S.R. A one penny imputed genome from next generation reference panels. Am. J. Hum. Genet.; 2018; 103, pp. 338-348. [DOI: https://dx.doi.org/10.1016/j.ajhg.2018.07.015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30100085]
34. Kang, H.M.; Sul, J.H.; Service, S.K.; Zaitlen, N.A.; Kong, S.Y.; Freimer, N.B.; Sabatti, C.; Eskin, E. Variance component model to account for sample structure in genome-wide genome-wide association studies. Nat. Genet.; 2010; 42, pp. 348-354. [DOI: https://dx.doi.org/10.1038/ng.548]
35. Yu, J.; Pressoir, G.; Briggs, W.H.; Bi, I.V.; Yamasaki, M.; Doebley, J.F.; McMullen, M.D.; Gaut, B.S.; Nielsen, D.M.; Holland, J.B. et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet.; 2006; 38, pp. 203-208. [DOI: https://dx.doi.org/10.1038/ng1702] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16380716]
36. Zhang, D.; Tan, Y.; Dong, F.; Zhang, Y.; Huang, Y.; Zhou, Y.; Zhao, Z.J.; Yin, Q.; Xie, X.; Gao, X. et al. The Expression of IbMYB1 Is Essential to Maintain the Purple Color of Leaf and Storage Root in Sweet Potato [Ipomoea batatas (L.) Lam]. Front. Plant Sci.; 2021; 12, 688707. [DOI: https://dx.doi.org/10.3389/fpls.2021.688707] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34630449]
37. Suematsu, K.; Tanaka, M.; Kurata, R.; Kai, Y. Comparative transcriptome analysis implied a ZEP paralog was a key gene involved in carotenoid accumulation in yellow-fleshed sweetpotato. Sci. Rep.; 2020; 10, 20607. [DOI: https://dx.doi.org/10.1038/s41598-020-77293-7]
38. Kim, D.; Langmead, B.; Salzberg1, S.L. HISAT: A fast spliced aligner with low memory requirements Daehwan HHS Public Access. Nat. Methods.; 2015; 12, pp. 357-360. [DOI: https://dx.doi.org/10.1038/nmeth.3317]
39. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc.; 2016; 11, pp. 1650-1667. [DOI: https://dx.doi.org/10.1038/nprot.2016.095]
40. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol.; 2010; 11, R106. [DOI: https://dx.doi.org/10.1186/gb-2010-11-10-r106]
41. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods.; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11846609]
42. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res.; 2009; 19, pp. 1655-1664. [DOI: https://dx.doi.org/10.1101/gr.094052.109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19648217]
43. Zhang, C.; Dong, S.S.; Xu, J.Y.; He, W.M.; Yang, T.L. PopLDdecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics.; 2019; 35, pp. 1786-1788. [DOI: https://dx.doi.org/10.1093/bioinformatics/bty875] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30321304]
44. Zhang, J.; Song, Q.; Cregan, P.B.; Nelson, R.L.; Wang, X.; Wu, J.; Jiang, G.L. Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm. BMC Genom.; 2015; 16, 217. [DOI: https://dx.doi.org/10.1186/s12864-015-1441-4]
45. He, X.; Liu, Q.; Wang, Y.; Zhai, H.; He, X.Q.; Liu, Q.C.; Wang, Y.P.; Zhai, H. Analysis of genetic diversity among sweetpotato landraces in China. Sci. Agric. Sin.; 2005; 38, pp. 250-257.
46. Holton, T.A.; Cornish, E.C. Genetics and biochemistry of anthocyanin biosynthesis. Plant Cell.; 1995; 7, pp. 1071-1083. [DOI: https://dx.doi.org/10.2307/3870058]
47. Ross, J.; Li, Y.; Lim, E.; Bowles, D.J. Higher plant glycosyltransferases. Genome Biol.; 2001; 2, pp. 3004.1-3004.6. [DOI: https://dx.doi.org/10.1186/gb-2001-2-2-reviews3004]
48. Ono, E.; Homma, Y.; Horikawa, M.; Kunikane-Doi, S.; Imai, H.; Takahashi, S.; Kawai, Y.; Ishiguro, M.; Fukui, Y.; Nakayama, T. Functional differentiation of the glycosyltransferases that contribute to the chemical diversity of bioactive flavonol glycosides in grapevines (Vitis vinifera). Plant Cell.; 2010; 22, pp. 2856-2871. [DOI: https://dx.doi.org/10.1105/tpc.110.074625]
49. Hall, D.; Yuan, X.X.; Murata, J.; De Luca, V. Molecular cloning and biochemical characterization of the UDP-glucose: Flavonoid 3-O-glucosyltransferase from Concord grape (Vitis labrusca). Phytochemistry; 2012; 74, pp. 90-99. [DOI: https://dx.doi.org/10.1016/j.phytochem.2011.10.007]
50. Doukhanina, E.V.; Chen, S.; Van Der Zalm, E.; Godzik, A.; Reed, J.; Dickman, M.B. Identification and functional characterization of the BAG protein family in Arabidopsis thaliana. J. Biol. Chem.; 2006; 281, pp. 18793-18801. [DOI: https://dx.doi.org/10.1074/jbc.M511794200]
51. Chi, Y.; Yang, Y.; Zhou, Y.; Zhou, J.; Fan, B.; Yu, J.Q.; Chen, Z. Protein-protein interactions in the regulation of WRKY transcription factors. Mol. Plant; 2013; 6, pp. 287-300. [DOI: https://dx.doi.org/10.1093/mp/sst026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23455420]
52. Cheng, Y.; Zhou, Y.; Yang, Y.; Chi, Y.J.; Zhou, J.; Chen, J.Y.; Wang, F.; Fan, B.; Shi, K.; Zhou, Y.H. et al. Structural and functional analysis of VQ motif-containing proteins in Arabidopsis as interacting proteins of WRKY transcription factors. Plant Physiol.; 2012; 159, pp. 810-825. [DOI: https://dx.doi.org/10.1104/pp.112.196816] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22535423]
53. Schluttenhofer, C.; Yuan, L. Regulation of specialized metabolism by WRKY transcription factors. Plant Physiol.; 2015; 167, pp. 295-306. [DOI: https://dx.doi.org/10.1104/pp.114.251769] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25501946]
54. Li, Q.; Kou, M.; Li, C.; Zhang, Y.G. Comparative transcriptome analysis reveals candidate genes involved in anthocyanin biosynthesis in sweetpotato (Ipomoea batatas L.). Plant Physiol. Biochem.; 2021; 158, pp. 508-517. [DOI: https://dx.doi.org/10.1016/j.plaphy.2020.11.035]
55. Zhang, L.; Yu, Y.; Shi, T.; Kou, M.; Sun, J.; Xu, T.; Li, Q.; Wu, S.; Cao, Q. Genome-wide analysis of expression quantitative trait loci ( eQTLs ) reveals the regulatory architecture of gene expression variation in the storage roots of sweet potato. Hortic. Res.; 2020; 7, 90. [DOI: https://dx.doi.org/10.1038/s41438-020-0314-4]
56. Gemenet, D.C.; da Silva Pereira, G.; De Boeck, B.; Wood, J.C.; Mollinari, M.; Olukolu, B.A.; Diaz, F.; Mosquera, V.; Ssali, R.T.; David, M. et al. Quantitative trait loci and differential gene expression analyses reveal the genetic basis for negatively associated β-carotene and starch content in hexaploid sweetpotato [Ipomoea batatas (L.) Lam.]. Theor. Appl. Genet.; 2020; 133, pp. 23-36. [DOI: https://dx.doi.org/10.1007/s00122-019-03437-7]
57. Yan, H.; Ma, M.; Ahmad, M.Q.; Arisha, M.H.; Tang, W.; Li, C.; Zhang, Y.; Kou, M.; Wang, X.; Gao, R. et al. High-Density Single Nucleotide Polymorphisms Genetic Map Construction and Quantitative Trait Locus Mapping of Color-Related Traits of Purple Sweet Potato [Ipomoea batatas (L.) Lam.]. Front. Plant Sci.; 2022; 12, 797041. [DOI: https://dx.doi.org/10.3389/fpls.2021.797041]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.