Abstract
Pigmented rice is desirable in healthy human diets because of potential health benefits compared to white rice. Based on the high diversity of Vietnamese black rice landrace, mapping the loci and predicting putative candidate genes for anthocyanin and flavonoid contents in these landraces would be useful to support the improvement of pigmented rice. A panel of 94 black rice accessions from the National Plant Genebank of Vietnam was analyzed to determine the anthocyanin and flavonoid levels in the whole grain and genotyped using an improved 7 K SNP array. Phenotypic and genotypic data were integrated in a genome wide association study. This identified 32 SNPs significantly associated with total anthocyanin content and 16 SNPs associated with total flavonoid content located on Chromosomes (Chrs.) 1, 2, 3, 6, 7, 8, 11 and 12 and on Chrs. 1, 3, 4, 6, 7, 8, 10, 11 and 12, respectively. A total of 72 candidate genes were co-located with the significant SNPs, of which 44 were predicted to be related to the pathways of biosynthesis of anthocyanin and flavonoid. Twenty candidate genes first reported in this study including the structural genes, LOC_Os03g62300.1 coding for chalcone isomerase and LOC_Os06g41800.1 coding for dihydroflavonol reductase and the regulatory gene, LOC_Os06g43090.1 coding for a MYB transcription factor, were strongly associated with anthocyanin synthesis. These candidate genes and significant SNPs could be used to analyze the function of genes and factors involved in anthocyanin synthesis and to improve pigmented rice by marker-assisted selection.
Similar content being viewed by others
Introduction
Pigmented rice is considered to have a high nutritional quality compared to white rice because it contains more bioactive and nutritional compounds. These compounds found in pigmented rice include, phytosterols, carotenoids, vitamins, micronutrients, and flavonoids (Mbanjo et al. 2020; Shao et al. 2018). Based on different accumulation of flavonoids, pigmented rice is classified into black and red rice varieties with the presence in the rice bran of anthocyanins and proanthocyanidins respectively (Min et al. 2011). Anthocyanins and other flavonoids in pigmented rice grain are known to have strong antioxidant activities that may be beneficial to human health (Semmarath et al. 2022; Shao et al. 2018; Guo et al. 2007; Dwiningsih and Al-Kahtani 2022; Kocic et al. 2011). Anthocyanins are a major subclass of flavonoids that are important for plants in responding to biotic and abiotic stress (Landi et al. 2015; Zaidi et al. 2019).
Anthocyanins are produced through a special branch of the flavonoid pathway and there are two gene groups, structural and regulatory genes, related to the biosynthesis that have been isolated had their functions confirmed (Winkel-Shirley 2001; Xia et al. 2021). The structural genes of rice related to enzymes and factors of anthocyanin pathway include the genes encoding chalcone synthase (OsCHS1 and OsCHS2) (Reddy et al. 1996; Shih et al. 2008), chalcone isomerase (OsCHI/gh) (Druka et al. 2003), flavanone 3-hydroxylase (OsF3H, OsF3H-1, OsF3H-2, and OsF3H-3) (Kim et al. 2008; Xia et al. 2021), flavonoid 3’-hydroxylase (OsF3’H) (Shih et al. 2008), Dihydroflavonol 4-reductase (OsDfr/Rd) (Furukawa et al. 2007; Xia et al. 2021), leucoanthocyanidin dioxygenase (OsANS1 and OsANS2) (Shih et al. 2008; Xia et al. 2021), and UDP-dependent glucosyltransferase (OsUGT) (Ko et al. 2008; Xia et al. 2021). The complex of regulatory genes included R2R3-MYB transcription factor (OsC1, OsP1 and Kala3), basic helix-loop-helix transcription factor (OsB1/Ra1, OsB2/Kala4, Rc and Rb), WD40 repeat (tryptophan-aspartate repeats) (Akhter et al. 2019; Liu et al. 2015; Koes et al. 2005; Hichri et al. 2011; Kim et al. 2011). In the anthocyanin pathway, the structural genes encode functional enzymes and the regulatory genes regulate the expression of structural genes.
The genetics of anthocyanin content and related traits in rice has been studied and documented by a large number of studies. The pigmentation of the pericarp has bene found to be controlled by two dominant complementary genes (Pp and Pb) that are located on rice chromosomes 1 and 4, respectively. A segregation ratio of 9:3:4 was reported for the three phenotypes, black pericarp, brown pericarp, and white pericarp in F2 populations (Tae-Ho et al. 2015; Rahman et al. 2013; Hsien and Chang 1964). Fine mapping of the Pb gene revealed that two other genes (Ra and bhlh16) closely linked with the Pb gene. Ra has been reported to have some common functions in anthocyanin synthesis and the deletion of 2-bp (GT) the Ra gene was found in all purple pericarp compared to white pericarp (Wang and Shu 2007). In addition, Rd and Rc genes were identified and related to red and brown rice on chromosome 1 and chromosome 7, respectively (Nagao et al. 1957; Yu et al. 1995; Sweeney et al. 2006; Nagao and Takahashi 1963). Rd and Rc genes were also identified as being involved in proanthocyanidin synthesis in rice pericarp at the LOC_Os07g11020 and LOC_Os01g39560 loci respectively (Furukawa et al. 2007; Sweeney et al. 2006; Xia et al. 2021). In addition to the two loci on chromosome 1 and 7, two other genomic regions having minor effects on the degree of red coloration in rice pericarp were found on chromosomes 9 and 11 (Dong et al. 2008). In another study, the black pericarp pigmentation was governed by the key loci for anthocyanin (KALA), named Kala1, Kala3, and Kala4 that were identified on chromosome 1, chromosome 3 and chromosome 4 by QTL mapping respectively and these three loci possibly relate to previously reported loci such as Pp, Rd, A, P, Pb, Pl and C (Maeda et al. 2014). Thus, many loci/genes for regulating anthocyanin/flavonoid synthesis were found and characterized by basic genetic analysis and QTL mapping.
Recently, genome-wide association studies (GWAS) using large number of SNPs, and high density of physical maps have been shown to be an efficient method to map QTLs and pinpoint the candidate genes for traits of interest in plants. For anthocyanin and flavonoid contents in rice, both structural genes and regulatory genes involved anthocyanin and flavonoid synthesis have been confirmed and characterized by GWAS and post GWAS. According to previous studies, a number of significant QTLs/SNPs associated with pericarp color in rice have been reported and validated by using the GWAS approach (Shao et al. 2011; Wang et al. 2016, 2020; Yang et al. 2018; Huang et al. 2010, 2012). Two genes, Os02g0650900 and Os08g0301500, encoding glutamate dehydrogenase and sucrose phosphate synthase, were suggested to be candidate genes underlying QTLs associated with pericarp color on chromosomes 2 and 8, respectively (Huang et al. 2010, 2012). Genes encoding chalcone isomerase, chalcone synthase, leucoanthocyanidin reductase, bHLH, MYB, WDR, and B-box underlying 30 significant chromosome positions from association mapping analysis were identified as candidate genes in proanthocyanidin synthesis in the rice grain pericarp (Haghi et al. 2021). Based on analysis of a whole population consisting of pigmented rice accessions and non-pigmented rice accessions and subpopulations of pigmented rice accessions, 26 and 18 significant SNP peaks associated with pigmented traits were identified for the whole population and the subpopulation, respectively. Of these, six peaks of significant SNPs were co-localized with the genes related to secondary metabolites, such as Os01g44260 (Rd), Os03g60509 (OsCHI), Os07g40570 (OsWRKY96), Os10g01480 (OsIPTK), and Os11g32650 (OsCHS1) for flavonoid pathway, and Os12g38400 (OsMYB91) for the phenolic pathway (Wang et al. 2023). Other, post GWAS studies used techniques such as gene expression, re-sequencing and RNA-sequencing to validate putative candidate genes from GWAS and some genes involved in anthocyanin pathway were confirmed in rice (Haghi et al. 2021; Yang et al. 2022; Wang et al. 2023). Consequently, most of the genes/factors involved anthocyanin and flavonoid synthesis have been reported by use of GWAS and post GWAS approaches.
In this study, an improved 7 K SNP dataset combined the top-performing SNPs from several previous rice arrays, including 4,007 SNPs from the C6AIR, 2,056 SNPs from the High Density Rice Array (HDRA), 910 SNPs from the 384-SNP GoldenGate sets, 189 SNPs from the 44 K array and 8 trait specific SNPs (Morales et al. 2020) was used for GWAS analysis of anthocyanin and flavonoid contents among black rice accessions from the National Plant Genebank of Vietnam. The study was carried out to map additional loci associated with anthocyanin and flavonoid contents and to predict putative candidate genes for genetic study and identify promising black rice accessions for breeding purposes.
Materials and methods
Plant materials
A set of 94 black rice accessions, collected from Northern and Central Vietnam, was selected from the National Plant Genebank of Vietnam. This diverse panel was divided into 75 Indica and 19 Japonica accessions based on phenol reaction. The seed coat color was purple and variable purple (S1 Table). The rice plants were grown in experimental fields in 2021 and seeds of each accession were obtained for analysis of total anthocyanin and total flavonoid contents.
Genotypic dataset
DNA of the 94 black rice accessions was extracted by using a CTAB protocol (Doyle and Doyle 1987) and genotyped with a Infinium 7k rice SNP chip (Morales et al. 2020). Over 7k SNP markers were obtained by running on the Illumina Infinium technology to create an SNP dataset for GWAS. This SNP dataset was reduced by exclusion of SNPs with greater than 5% missing data and a minor allele frequency (MAF) of less than 5% to obtain 3744 qualified SNP using TASSEL software (Bradbury et al. 2007).
Analyses of total anthocyanin and flavonoid contents
Total anthocyanin content of brown rice was measured by a UV-visible spectroscopy method (Giusti and Wrolstad 2001) with minor modifications. Briefly, 10 g of powder of brown rice was kept at − 15 °C and mixed with 400 ml solution of 1 V ethanol, 1 V distilled water and 1% HCl for 60 min. The mixture was centrifuged 10 min at 5000 rpm to obtain a supernatant. The absorbance of each dilution was measured at the λvis−max and at 700 nm, against a blank cell filled with distilled water. The absorbance of the diluted sample (A) was calculated as follows:
Then, the anthocyanin concentration in the original sample was calculated by using the following formula:
where MW is the molecular weight, DF is the dilution, ε is the molar absorptivity and l is the thickness of cuvette.
The powder of brown rice was also used to measure total flavonoid content by a spectrophotometric method based on the formation of aluminium-flavonoid complexes (Pękal and Pyrzynska 2014) with minor modifications. Briefly, 1 g powder of brown rice was mixed in 10 ml methanol and filtered to obtain solution. The test solution (1 ml of standard or sample), was mixed with 0.3 ml of NaNO2 (5%) and after 5 min, 0.3 ml of AlCl3 (10%) was added. The sample was mixed and 6 min later was neutralized with 2 ml of 1 N NaOH solution and made-up to 10 ml. The mixture was left for 45 min at room temperature and then subjected to spectral analysis at 510 nm against the blank, where AlCl3 solution was substituted by water. Quercetin (in the 20–100 mg/ml concentration range) was the standard of choice for the expression of results at 510 nm. The flavonoid concentration in the original sample was calculated by using the following formula:
where c is the constant from the standard evaluation, V is the volume of extracted solution and m is weight of powder.
Genome-wide association study (GWAS)
Linkage disequilibrium (LD) and correlation coefficients (r2) were measured on pairs within a sliding window and a fitted curve was computed based on nonlinear regression of LD (r2) versus distance (kb) as previously defined (Remington et al. 2001). Principle components (PCs) and kinship matrix (relatedness) were applied to correct for population structure and relatedness in mixed linear models. Identity by state (IBS) matrices indicating relatedness among rice accessions were calculated using TASSEL 5 (Li and Zhu 2013) to construct phylogenic trees in MEGA 7 (Kumar et al. 2016) and to apply as kinship matrix in GWAS. Genome-wide associations between SNPs and total anthocyanin and flavonoid contents were identified using the mixed linear model with correction for population structure and relatedness in TASSEL 5 (Li and Zhu 2013). False positives were controlled by multiple test correction with false discovery rate (FDR) ≤ 0.05 and the threshold of -log10(p–value) for identifying significant associations was calculated at FDR = 0.05 (Qu et al. 2010).
Predicting candidate genes related total anthocyanin and flavonoid contents
The search for candidate genes was conducted by zooming in at positions of the significant SNPs associated with total anthocyanin and flavonoid contents. Information about the candidate genes was displayed by searching in the Genome Browser of phytozome website (https://phytozome.jgi.doe.gov/jbrowse/) for Oryza sativa v7.0 and Gene loci search in Rice SNP-Seek Database (https://snp-seek.irri.org/_locus.zul) (Mansueto et al. 2016). The putative candidate genes related to interesting traits including total anthocyanin and flavonoid contents were identified by functional annotation of genes that were located near significant SNPs in the physical map (Haghi et al. 2022; Do et al. 2019).
Results
Phenotypic variation of total anthocyanin and total flavonoid contents in 94 black rice accessions
Based on evaluation of total anthocyanin and total flavonoid contents, the phenotypic variation among 94 black rice accessions was significant with ranging from 3.6 to 759.5 mg/100 g and 0.6–11.4 mg/100 g for total anthocyanin content and total flavonoid content, respectively. Total anthocyanin and total flavonoid contents in the panel of black rice accessions were likely a continuous distribution with the graphs being approximately bell‑shaped and symmetric about the mean. However, most accessions were distributed left of the bell-shaped curve with values lower than mean for both total anthocyanin content and total flavonoid content (Fig. 1).
Population structure and relatedness among black rice accessions
Linkage disequilibrium (LD) decay, population structure and relatedness were analyzed for the panel of 94 black rice accessions using 3744 SNPs. LD decay on adjacent SNP pairs were in nonlinear curves (Fig. S1) and the LD blocks at 1/2r2 of 0.1 were an average of 7.5 Mb covering the genome-wide haplotype blocks for the SNP dataset. The principal component analysis (PCA) showed that the cumulative eigenvalue of the first three PCs were 60% of variances for the black rice panel. Based on the first three PCs, 94 black rice accessions were separated into groups roughly corresponding to Indica and Japonica groups (Fig. 2A). The cryptic relatedness among rice accessions was evaluated by kinship matrix from identity by state (IBS) each from paired accessions. The matrix was calculated using the 3744 SNPs to construct a phylogenetic tree showing the relationship among 94 accessions (Fig. 2B). In this dendrogram, four main groups, I, II, II and IV groups were clustered at the lowest IBS value and there were subgroups at higher IBS values. Thus, adding PCs and the kinship matrix into the GWAS model was efficient in control of false positives because of population structure or cryptic relatedness.
GWAS for total anthocyanin content of the black rice collection
A total of 3744 polymorphic SNPs were selected after marker quality control and assurance for further analysis of total anthocyanin in the GWAS. The quantile-quantile (Q-Q) distribution showed most SNPs matched what was expected (Fig. 3A). The number of SNPs significantly associated with total anthocyanin content was 32 based on correction Benjamini-Hochberg of False Discovery Rate (FDR) ≤ 0.05. The significant SNPs associated with total anthocyanin content were located on Chrs. 1, 2, 3, 6, 7, 8, 11 and 12 and the number of significant SNPs were 4, 1, 5, 8, 1, 10, 2 and 1, respectively (Table 1; Fig. 3B). The most significant SNP was id3003846 with – log10(P) of 6.4 that was found on Chr. 3, however, the most significant segment is probably that located on Chr. 6 from position of ss6873320 to ss6893305’s.
GWAS for total anthocyanin content of the black rice collection
The association between SNPs and total flavonoid content was also carried out by using the mixed linear model in TASSEL 5. The result showed that most of the expected SNPs were close to the observed SNPs in the quantile-quantile (Q-Q) plot (Fig. 4A). Based on correction for Benjamini–Hochberg of False Discovery Rate (FDR) ≤ 0.05, – log(P) the threshold was 3.6 for GWAS of total flavonoid content. A total of 16 significant SNPs associated with total flavonoid content were located on Chrs. 1, 3, 4, 6, 7, 8, 10, 11 and 12 (Table 2; Fig. 4B). All significant SNPs were identified with –log(P) values ranging from 3.6 to 7.4 of which the most significant was – log(P) of 7.4 on Chr. 10. GWASs of total anthocyanin and total flavonoid contents revealed significant SNPs found on most of the chromosomes in the rice genome.
Putative candidate genes related total anthocyanin and flavonoid contents
Phytozome genome browsers (https://phytozome-next.jgi.doe.gov/jbrowse/) and Gene loci search in Rice SNP-Seek Database (https://snp-seek.irri.org/_locus.zul) (Mansueto et al. 2016) were used to detect putative candidate genes in the significant genomic regions associated with total anthocyanin and total flavonoid contents. The genomic intervals surrounding the significant SNPs were considered significant genomic regions to search for putative candidate genes. A total of 72 candidate genes were found in the significant genomic intervals on each chromosome (Table 2). Among these, 20 genes included the significant SNPs and other genes located near the significant SNPs associated with the total anthocyanin content or the total flavonoid content. In addition, 44 genes had predicted functions related to biosynthesis pathways of anthocyanin in the Phytozome database/Rice SNP-Seek Database (Table S2). Based on functional annotation, the candidate genes were grouped into structural genes (29), regulatory genes (9) and genes related to malonyl-CoA/coumaroyl-CoA (6) (Table S2). Notably, LOC_Os06g46920.1 and LOC_Os06g41800.1 coding for dihydroflavonol reductase, the key enzyme in the metabolic pathway of anthocyanin were found by GWAS of total anthocyanin and total flavonoid contents, respectively (Fig. 5).
Discussion
The Infinium 7k rice SNP chip was shown to be useful for quantitative trait loci (QTL) and association mapping in diverse materials and 20 QTLs contributing to plant height were identified by GWAS analysis among 189 rice accessions (Morales et al. 2020). Other successful applications of the 7k rice SNP chip for association mapping have been reported, such as 51 QTLs for the different chilling indices identified by GWAS analysis of a panel of 257 rice accessions (Thapa et al. 2020), and 17 SNPs, significantly associated with root plasticity traits under soil moisture stress conditions located on chromosomes 2, 5, 7, 9, 10, and 12 (Niones et al. 2021). In this study, the Infinium 7k rice SNP chip was also validated to be efficient in GWAS analysis of anthocyanin and flavonoid contents with LD blocks averaging 7.5 Mb at 1/2r2 of 0.1, the first three PCs reflected 60% of the variances for the black rice panel, kinship matrix showed the relationship among 94 accessions and number of significant SNPs associated with both traits.
Rice pericarp color and related traits including anthocyanin and flavonoid contents are controlled by multiple genes (Furukawa et al. 2007; Dong et al. 2008; Maeda et al. 2014). By association mapping, more loci for anthocyanin and flavonoid contents in rice have been identified that located on most of the chromosomes such as 30 significant chromosome positions associated with anthocyanin pigmentation in rice (Haghi et al. 2021), and 11 loci for pericarp color (Yang et al. 2022). In a more recent study, 26 significant SNPs associated with rice pericarp color were mapped on chromosomes (Chrs.) 1, 2, 3, 4, 6, 7, 9, 10, 11, and 12 (Wang et al. 2023). By using a mixed model for GWAS and FDR tests, 32 SNPs significantly associated with total anthocyanin content were found to be contributed from Chrs. 1, 2, 3, 6, 7, 8, 11 and 12 (Table 3; Fig. 3), while 16 significant SNPs associated with total flavonoid content were located on Chrs. 1, 3, 4, 6, 7, 8, 10, 11 and 12 (Table 1; Fig. 4). No significant SNPs were found on Chrs. 5 and 9, however significant SNPs identified were identified for the first time on Chr. 8 by GWAS for both total anthocyanin and total flavonoid contents compared to previous studies.
Candidate genes involved in anthocyanin and flavonoid synthesis were predicted by some approaches including GWAS for different tissues of rice. Though genome-wide identification and expression profiles, 85 key structural gene candidates belong to 13 families were suggested for flavonoid biosynthesis in rice (Wang et al. 2022). Based on 8 significant SNPs associated with anthocyanin pigmentation in rice stem, 19 putative candidate genes were highlighted for anthocyanin synthesis (Haghi et al. 2022).The genes involved in the anthocyanin pathway have been identified and classified as structural and regulatory genes in which structural genes encode functional enzymes and regulatory genes encode mainly transcription factors (TFs). The genes involved in anthocyanin synthesis are well described in plants (Yonekura-Sakakibara et al. 2019; Xia et al. 2021; Haghi et al. 2022; Mackon et al. 2021). By Predicting candidate genes related total anthocyanin and flavonoid contents in this study, most of the genes related to anthocyanin pathway were identified, however some of them were located in new genomic regions (Table S2) that could be considered as new candidate genes for anthocyanin pigmentation in rice.
The first step of anthocyanin biosynthesis that branches from the general flavonoid pathway is the convert of coumaroyl-CoA and malonyl-CoA to nargenine chalcone (Mackon et al. 2021; Haghi et al. 2021, 2022; Xia et al. 2021). According to functional annotation, some putative candidate genes for anthocyanin synthesis (LOC_Os01g55590.1, LOC_Os06g44620.1, LOC_Os07g06660.1, LOC_Os08g07720.1, LOC_Os08g07730.1 and LOC_Os12g41250.1) (Table S2) encoding enzymes that are related to coumaroyl-CoA and malonyl-CoA could be involved in anthocyanin synthesis. Structural genes coding chalcone isomerase, dihydroflavonol reductase and UDP-glycosyl transferase, and regulatory genes coding for MYB transcription factor are suggested as important genes for the anthocyanin pathway (Kim et al. 2021; Xia et al. 2021; Mackon et al. 2021). Of the first reported genes, the structural genes (LOC_Os03g62300.1 and LOC_Os06g41800.1) coding for chalcone isomerase and dihydroflavonol reductase and the regulatory gene (LOC_Os06g43090.1) coding for a MYB transcription factor were strongly associated with anthocyanin synthesis.
In conclusion, GWAS results showed that the using the Infinium 7k rice SNP in association mapping was an efficient way to detect significant SNPs associated with anthocyanin and flavonoid contents in rice. Based on applying mix model and correction tests, the significant SNPs from GWAS for anthocyanin and flavonoid contents were significant and accurate for the development of SNP markers for marker-assisted selection. LOC_Os03g62300.1, LOC_Os06g41800.1 (chalcone isomerase) and LOC_Os06g43090.1 (dihydroflavonol reductase) were strongly assocuiated with anthocyanin synthesis. These genes and associated SNPs should be considered for us in studies of pigmented rice in the future.
References
Akhter D, Qin R, Nath UK, Eshag J, Jin X, Shi C (2019) A rice gene, OsPL, encoding a MYB family transcription factor confers anthocyanin synthesis, heat stress response and hormonal signaling. Gene 699:62–72. https://doi.org/10.1016/j.gene.2019.03.013
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES (2007) TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23(19):2633–2635. https://doi.org/10.1093/bioinformatics/btm308
Do TD, Vuong TD, Dunn D, Clubb M, Valliyodan B, Patil G, Chen P, Xu D, Nguyen HT, Shannon JG (2019) Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping. BMC Genomics 20(1):318. https://doi.org/10.1186/s12864-019-5662-9
Dong Y, Xu J, Xiao K, Zhang Y, Zhang J, Luo L, Matsuo M (2008) Genomic regions associated with the degree of red coloration in pericarp of rice (Oryza sativa L). J Cereal Sci 48(2):556–560. https://doi.org/10.1016/j.jcs.2007.11.011
Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bull 19:11–15
Druka A, Kudrna D, Rostoks N, Brueggeman R, von Wettstein D, Kleinhofs A (2003) Chalcone isomerase gene from rice (Oryza sativa) and barley (Hordeum vulgare): physical, genetic and mutation mapping. Gene 302(1–2):171–178. https://doi.org/10.1016/s0378-1119(02)01105-8
Dwiningsih Y, Al-Kahtani J (2022) Genetics, biochemistry and biophysical analysis of anthocyanin in rice (Oryza sativa L.). Adv Sustain Sci Eng Technol 4(1):1–19. https://doi.org/10.26877/asset.v4i1.11659
Furukawa T, Maekawa M, Oki T, Suda I, Iida S, Shimada H, Takamure I, Kadowaki K (2007) The rc and rd genes are involved in proanthocyanidin synthesis in rice pericarp. Plant J 49(1):91–102. https://doi.org/10.1111/j.1365-313X.2006.02958.x
Giusti MM, Wrolstad RE (2001) Characterization and measurement of anthocyanins by UV–visible spectroscopy. Curr Protoc Food Anal Chem. https://doi.org/10.1002/0471142913.faf0102s00
Guo H, Ling W, Wang Q, Liu C, Hu Y, Xia M, Feng X, Xia X (2007) Effect of anthocyanin-rich extract from black rice (Oryza sativa L. indica) on hyperlipidemia and insulin resistance in fructose-fed rats. Plant Foods Hum Nutr 62(1):1–6. https://doi.org/10.1007/s11130-006-0031-7
Haghi R, fazeli A, Ahmadikhah A, Shariati V (2021) Identification of novel genes involved in anthocyanin and proanthocyanidin pigments accumulation in rice tissues using genome-wide association study (GWAS). J Plant Res (Iran J Biol) 34(2):369–384
Haghi R, Ahmadikhah A, Fazeli A, Shariati V (2022) Candidate genes for anthocyanin pigmentation in rice stem revealed by GWAS and whole-genome resequencing. Plant Genome-Us 15(3):e20224. https://doi.org/10.1002/tpg2.20224
Ham TH, Kwon SW, Ryu SN, Koh HJ (2015) correlation analysis between grain color and cyanidin-3-glucoside content of rice grain in segregate population. Plant Breed Biotechnol 3(2):160–166. https://doi.org/10.9787/PBB.2015.3.2.160
Hichri I, Barrieu F, Bogs J, Kappel C, Delrot S, Lauvergeat V (2011) Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway. J Exp Bot 62(8):2465–2483. https://doi.org/10.1093/jxb/erq442
Hsien S-C, Chang T-M (1964) Genic analysis in Rice. IV. Genes for purple pericarp and other characters. Jpn J Breed 14(3):141–149
Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, Li C, Zhu C, Lu T, Zhang Z, Li M, Fan D, Guo Y, Wang A, Wang L, Deng L, Li W, Lu Y, Weng Q, Liu K, Huang T, Zhou T, Jing Y, Li W, Lin Z, Buckler ES, Qian Q, Zhang Q-F, Li J, Han B (2010) Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet 42:961. https://doi.org/10.1038/ng.695
Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, Li W, Guo Y, Deng L, Zhu C, Fan D, Lu Y, Weng Q, Liu K, Zhou T, Jing Y, Si L, Dong G, Huang T, Lu T, Feng Q, Qian Q, Li J, Han B (2012) Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet 44(1):32–39. https://doi.org/10.1038/ng.1018
Kim JH, Lee YJ, Kim BG, Lim Y, Ahn JH (2008) Flavanone 3beta-hydroxylases from rice: key enzymes for favonol and anthocyanin biosynthesis. Mol Cells 25(2):312–316
Kim CK, Cho MA, Choi YH, Kim JA, Kim YH, Kim YK, Park SH (2011) Identification and characterization of seed-specific transcription factors regulating anthocyanin biosynthesis in black rice. J Appl Genet 52(2):161–169. https://doi.org/10.1007/s13353-011-0027-3
Kim D-H, Yang J, Ha S-H, Kim JK, Lee J-Y, Lim S-H (2021) An OsKala3, R2R3 MYB TF, is a common key player for black rice pericarp as main partner of an OsKala4, bHLH TF. Front Plant Sci 12:765049. https://doi.org/10.3389/fpls.2021.765049
Ko JH, Kim BG, Kim JH, Kim H, Lim CE, Lim J, Lee C, Lim Y, Ahn JH (2008) Four glucosyltransferases from rice: cDNA cloning, expression, and characterization. J Plant Physiol 165(4):435–444. https://doi.org/10.1016/j.jplph.2007.01.006
Kocic B, Filipovic S, Nikolic M, Petrovic B (2011) Effects of anthocyanins and anthocyanin-rich extracts on the risk for cancers of the gastrointestinal tract. J Buon 16(4):602–608
Koes R, Verweij W, Quattrocchio FM (2005) Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci 10(5):236–242
Kumar S, Stecher G, Tamura K (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 33(7):1870–1874. https://doi.org/10.1093/molbev/msw054
Landi M, Tattini M, Gould KS (2015) Multiple functional roles of anthocyanins in plant-environment interactions. Environ Exp Bot 119:4–17. https://doi.org/10.1016/j.envexpbot.2015.05.012
Li G, Zhu H (2013) Genetic studies: the linear mixed models in genome-wide association studies. Open Bioinform J 7:27–33
Liu J, Osbourn A, Ma P (2015) MYB transcription factors as regulators of phenylpropanoid metabolism in plants. Mol Plant 8(5):689–708. https://doi.org/10.1016/j.molp.2015.03.012
Mackon E, Jeazet Dongho Epse Mackon GC, Ma Y, Haneef Kashif M, Ali N, Usman B, Liu P (2021) Recent insights into anthocyanin pigmentation, synthesis, trafficking, and regulatory mechanisms in rice (Oryza sativa L.) caryopsis. Biomolecules. https://doi.org/10.3390/biom11030394
Maeda H, Yamaguchi T, Omoteno M, Takarada T, Fujita K, Murata K, Iyama Y, Kojima Y, Morikawa M, Ozaki H, Mukaino N, Kidani Y, Ebitani T (2014) Genetic dissection of black grain rice by the development of a near isogenic line. Breed Sci 64(2):134–141. https://doi.org/10.1270/jsbbs.64.134
Mansueto L, Fuentes RR, Borja FN, Detras J, Abriol-Santos JM, Chebotarov D, Sanciangco M, Palis K, Copetti D, Poliakov A, Dubchak I, Solovyev V, Wing RA, Hamilton RS, Mauleon R, McNally KL, Alexandrov N (2016) Rice SNP-seek database update: new SNPs, indels, and queries. Nucleic Acids Res 45(D1):D1075–D1081. https://doi.org/10.1093/nar/gkw1135
Mbanjo EGN, Kretzschmar T, Jones H, Ereful N, Blanchard C, Boyd LA, Sreenivasulu N (2020) The genetic basis and nutritional benefits of pigmented rice grain. Front Genet 11:229. https://doi.org/10.3389/fgene.2020.00229
Min B, McClung AM, Chen MH (2011) Phytochemicals and antioxidant capacities in rice brans of different color. J Food Sci 76(1):C117–C126. https://doi.org/10.1111/j.1750-3841.2010.01929.x
Morales KY, Singh N, Perez FA, Ignacio JC, Thapa R, Arbelaez JD, Tabien RE, Famoso A, Wang DR, Septiningsih EM, Shi Y, Kretzschmar T, McCouch SR, Thomson MJ (2020) An improved 7K SNP array, the C7AIR, provides a wealth of validated SNP markers for rice breeding and genetics studies. PLoS ONE 15(5):e0232479. https://doi.org/10.1371/journal.pone.0232479
Nagao S, Takahashi M-e (1963) Trial construction of twelve linkage groups in Japanese Rice: (genetical studies on rice plant, XXVII). J Fac Agric Hokkaido Univ 53(1):72–130
Nagao S, Takahashi M-e, Miyamoto T (1957) Genetical studies on rice plant, XXI biochemical studies on red rice pigmentation. Jpn J Genet 32(4):124–128. https://doi.org/10.1266/jjg.32.124
Niones JM, Lipio PLG, Cruz AS, Cabral MCJ, Hautea DM, Lucob-Agustin N, Suralta RR (2021) Genome-wide association mapping for the identification of SNPs controlling lateral root plasticity in selected rice germplasms of the Philippines. Philippine J Sci 150(3):663–674
Pękal A, Pyrzynska K (2014) Evaluation of aluminium complexation reaction for flavonoid content assay. Food Anal Methods 7(9):1776–1782. https://doi.org/10.1007/s12161-014-9814-x
Qu HQ, Tien M, Polychronakos C (2010) Statistical significance in genetic association studies. Clin Invest Med 33(5):E266–E270
Rahman MM, Lee KE, Lee ES, Matin MN, Lee DS, Yun JS, Kim JB, Kang SG (2013) The genetic constitutions of complementary genes pp and pb determine the purple color variation in pericarps with cyanidin-3-O-glucoside depositions in black rice. J Plant Biol 56(1):24–31. https://doi.org/10.1007/s12374-012-0043-9
Reddy AR, Scheffler B, Madhuri G, Srivastava MN, Kumar A, Sathyanarayanan PV, Nair S, Mohan M (1996) Chalcone synthase in rice (Oryza sativa L.): detection of the CHS protein in seedlings and molecular mapping of the chs locus. Plant Mol Biol 32(4):735–743. https://doi.org/10.1007/bf00020214
Remington DL, Thornsberry JM, Matsuoka Y, Wilson LM, Whitt SR, Doebley J, Kresovich S, Goodman MM, Buckler ES (2001) Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc Natl Acad Sci USA 98(20):11479–11484. https://doi.org/10.1073/pnas.201394398
Semmarath W, Mapoung S, Umsumarng S, Arjsri P, Srisawad K, Thippraphan P, Yodkeeree S, Dejkriengkraikul P (2022) Cyanidin-3-O-glucoside and peonidin-3-o-glucoside-rich fraction of black rice germ and bran suppresses inflammatory responses from SARS-CoV-2 spike glycoprotein S1-induction in vitro in A549 lung cells and THP-1 macrophages via Inhibition of the NLRP3 inflammasome pathway. Nutrients 14(13):2738. https://doi.org/10.3390/nu14132738
Shao Y, Hu Z, Yu Y, Mou R, Zhu Z, Beta T (2018) Phenolic acids, anthocyanins, proanthocyanidins, antioxidant activity, minerals and their correlations in non-pigmented, red, and black rice. Food Chem 239:733–741. https://doi.org/10.1016/j.foodchem.2017.07.009
Shao Y, Jin L, Zhang G, Lu Y, Shen Y, Bao J (2011) Association mapping of grain color, phenolic content, flavonoid content and antioxidant capacity in dehulled rice. Theor Appl Genet 122(5):1005–1016. https://doi.org/10.1007/s00122-010-1505-4
Shih CH, Chu H, Tang LK, Sakamoto W, Maekawa M, Chu IK, Wang M, Lo C (2008) Functional characterization of key structural genes in rice flavonoid biosynthesis. Planta 228(6):1043–1054. https://doi.org/10.1007/s00425-008-0806-1
Sweeney MT, Thomson MJ, Pfeil BE, McCouch S (2006) Caught red-handed: Rc encodes a basic helix-loop-helix protein conditioning red pericarp in rice. Plant Cell 18(2):283–294. https://doi.org/10.1105/tpc.105.038430
Thapa R, Tabien RE, Thomson MJ, Septiningsih EM (2020) Genome-wide association mapping to identify genetic loci for cold tolerance and cold recovery during germination in rice. Front Genet 11:22. https://doi.org/10.3389/fgene.2020.00022
Wang C, Shu Q (2007) Fine mapping and candidate gene analysis of purple pericarp gene pb in rice (Oryza sativa L). Chin Sci Bull 52(22):3097–3104. https://doi.org/10.1007/s11434-007-0472-x
Wang H, Xu X, Vieira FG, Xiao Y, Li Z, Wang J, Nielsen R, Chu C (2016) The power of inbreeding: NGS-based GWAS of Rice reveals convergent evolution during rice domestication. Mol Plant 9(7):975–985. https://doi.org/10.1016/j.molp.2016.04.018
Wang W, Zhao M, Zhang G, Liu Z, Hua Y, Jia X, Song J, Ma D, Sun J (2020) Weedy rice as a novel gene resource: a genome-wide association study of anthocyanin biosynthesis and an evaluation of nutritional quality. Front Plant Sci 11:878. https://doi.org/10.3389/fpls.2020.00878
Wang J, Zhang C, Li Y (2022) Genome-wide identification and expression profiles of 13 key structural gene families involved in the biosynthesis of rice flavonoid scaffolds. Genes (Basel) 13(3):410. https://doi.org/10.3390/genes13030410
Wang W, Qiu X, Wang Z, Xie T, Sun W, Xu J, Zhang F, Yu S (2023) Deciphering the genetic architecture of color variation in whole grain rice by genome-wide association. Plants. https://doi.org/10.3390/plants12040927
Winkel-Shirley B (2001) Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol 126(2):485–493. https://doi.org/10.1104/pp.126.2.485
Xia D, Zhou H, Wang Y, Li P, Fu P, Wu B, He Y (2021) How rice organs are colored: the genetic basis of anthocyanin biosynthesis in rice. Crop J 9(3):598–608. https://doi.org/10.1016/j.cj.2021.03.013
Yang X, Xia X, Zeng Y, Nong B, Zhang Z, Wu Y, Xiong F, Zhang Y, Liang H, Deng G, Li D (2018) Identification of candidate genes for gelatinization temperature, gel consistency and pericarp color by GWAS in rice based on SLAF-sequencing. PLoS ONE 13(5):e0196690. https://doi.org/10.1371/journal.pone.0196690
Yang W, Chen L, Zhao J, Wang J, Li W, Yang T, Dong J, Ma Y, Zhou L, Chen J, Wu W, Zhang S, Liu B (2022) Genome-wide association study of pericarp color in rice using different germplasm and phenotyping methods reveals different genetic architectures. Front Plant Sci 13:841191. https://doi.org/10.3389/fpls.2022.841191
Yonekura-Sakakibara K, Higashi Y, Nakabayashi R (2019) The origin and evolution of plant flavonoid metabolism. Front Plant Sci 10:943. https://doi.org/10.3389/fpls.2019.00943
Yu ZH, McCouch SR, Tanksley SD, Kinoshita T, Sato S (1995) Association of morphological and RFLP markers in rice (Oryza sativa L). Genome 38(3):566–574. https://doi.org/10.1139/g95-073
Zaidi SHR, Zakari SA, Zhao Q, Khan AR, Shah JM, Cheng F (2019) Anthocyanin accumulation in black kernel mutant Rice and its contribution to ROS detoxification in response to high temperature at the filling stage. Antioxidants 8(11):510
Acknowledgements
The article authors would like to thank the Ministry of Science and Technology (MOST) for funding to carry out the research project code NĐT.77.AU/20.
Author information
Authors and Affiliations
Contributions
CHP conceived and designed the study, reviewed and edited the manuscript; TDD analysed data and wrote the draft version of the manuscript; HTLN undertook the acquisition and analysis of data for the article; NTH acquired and interpreted the morphological data; TDT made a substantial contribution to revising critically for important intellectual content for the article; MTTV analysed the data for the article; HHD made the acquisition and analysis of morphological data; TGTB analyzed SNP data; RJH reviewed and edited the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
10681_2023_3268_MOESM2_ESM.jpg
Figure S1: Linkage disequilibrium (LD) decay plot of coefficient of correlation (r2) between adjacent marker pairs and genomic distance (kb) (JPG 746 kb)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pham, C.H., Do, T.D., Nguyen, H.T.L. et al. Genome-wide association mapping of genes for anthocyanin and flavonoid contents in Vietnamese landraces of black rice. Euphytica 220, 11 (2024). https://doi.org/10.1007/s10681-023-03268-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10681-023-03268-0