Next Article in Journal
Extra Virgin Olive Oil Phenolic Compounds Modulate the Gene Expression of Biomarkers Involved in Fibroblast Proliferation and Differentiation
Next Article in Special Issue
Progress in Rice Breeding Based on Genomic Research
Previous Article in Journal
Progressive Myoclonus Epilepsy: A Scoping Review of Diagnostic, Phenotypic and Therapeutic Advances
Previous Article in Special Issue
Genome-Wide Identification and Characterization of the VQ Motif-Containing Gene Family Based on Their Evolution and Expression Analysis under Abiotic Stress and Hormone Treatments in Foxtail Millet (Setaria italica L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach

1
MARA Key Laboratory of Sustainable Crop Production in the Middle Reaches of the Yangtze River (Co-Construction by Ministry and Province), College of Agriculture, Yangtze University, Jingzhou 434025, China
2
State Key Laboratory of Vegetable Biobreeding, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
3
Zhongyuan Research Center, Chinese Academy of Agricultural Sciences, Xinxiang 453000, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2024, 15(2), 172; https://doi.org/10.3390/genes15020172
Submission received: 21 December 2023 / Revised: 19 January 2024 / Accepted: 23 January 2024 / Published: 28 January 2024
(This article belongs to the Special Issue Genomic Studies of Plant Breeding)

Abstract

:
Spinach is a significant source of vitamins, minerals, and antioxidants. These nutrients make it delicious and beneficial for human health. However, the genetic mechanism underlying the accumulation of nutrients in spinach remains unclear. In this study, we analyzed the content of chlorophyll a, chlorophyll b, oxalate, nitrate, crude fiber, soluble sugars, manganese, copper, and iron in 62 different spinach accessions. Additionally, 3,356,182 high-quality, single-nucleotide polymorphisms were found using resequencing and used in a genome-wide association study. A total of 2077 loci were discovered that significantly correlated with the concentrations of the nutritional elements. Data mining identified key genes in these intervals for four traits: chlorophyll, oxalate, soluble sugar, and Fe. Our study provides insights into the genetic architecture of nutrient variation and facilitates spinach breeding for good nutrition.

1. Introduction

The Spinach (Spinacia oleracea L.), belonging to the Amaranthaceae family, is primarily valued for its leaves as the main edible organ [1]. It is considered one of the healthiest vegetables in the human diet due to its richness in nutrients and health-promoting compounds, such as chlorophyll, vitamins, soluble sugars, and mineral elements, although it also contains some anti-nutritional factors, such as oxalate and nitrate [2]. The Food and Agriculture Organization of the United Nations (FAO) reports that in 2021, Asia produced 30.86 million tons of spinach, or 95.55% of the world’s total, with China producing 96.69% of that amount (https://www.fao.org/faostat (accessed on 25 October 2023)).
Spinach is a good source of nutrients, which is why it is popular among consumers. Mineral nutrients, for example, are required for plant and animal growth and development, and humans can only receive them via plant-based meals [3]. Among these elements, trace elements such as Fe, Mn, and Cu are important cofactors in the structure of certain enzymes and are essential in many biochemical pathways [4]. Crude fiber is a kind of dietary fiber made mostly of cellulose, hemicellulose, and lignin. It improves fecal volume, promotes intestinal peristalsis, prevents constipation and rectal cancer, and reduces obesity [5]. Additionally, cellulose, as a primary component of plant cell walls, directly influences the texture, palatability, and degree of tenderness of vegetables [6]. In addition to being rich in nutrients, spinach also contains more oxalate than most other crops [7]. Oxalate has been shown to exist in two forms: soluble oxalate and insoluble oxalate, with soluble oxalate being the predominant type found in spinach [8]. When absorbed by the digestive tract, oxalate is easily chelated by calcium ions, magnesium ions, and other metal ions to form insoluble oxalates, thus negatively affecting human health by reducing mineral absorption and promoting the formation of kidney stones [9,10].
Previous studies have discovered large variances in nutritional content within spinach; the reasons for these variations may include factors such as variety, growth conditions, and the environment [11,12,13,14]. In recent years, many studies have been conducted to evaluate spinach’s nutritional content and molecular basis. Fourteen mineral elements (Fe, Cu, Mn, and more) were examined in a panel of 292 spinach germplasm resources. Mineral element levels varied significantly, ranging from 50.43 to 138.58 μg/g for Fe, 2.22 to 10.97 μg/g for Cu, and 29.78 to 415.76 μg/g for Mn [14]. In one investigation, oxalate levels varied from 647.2 to 1056.2 mg/100 g (fresh weight) or from 53.4 to 79.9 mg/100 g (dry weight) [15]. Various methodologies have been used to identify the governing genes and relevant loci that influence these variances. Fourteen loci linked with Mn, Fe, and Cu contents in spinach were detected using a genome-wide association study (GWAS) with sequencing genotyping (GBS) data [14]. The GWAS study identified a potential area for oxalate at 915,444–1,022,853 bp on chromosome 5 [16]. Due to a scarcity of markers and the lack of a high-quality reference genome, potential genes related to these mineral elements are difficult to uncover.
GWAS is the study of phenotype–genotype relationships at the population level. It does not require the construction of genetic linkage maps; it can detect the genetic structure of phenotypes and identify numerous candidate genes [17]. When compared to approaches such as quantitative trait loci (QTL) mapping, it is faster, less expensive, and provides more precise localization. GWAS has been widely utilized to identify genes associated with quantitative traits in various crops, such as rice [18], maize [19], snap bean [20], tomato [21], cucumber [22], pepper [23], melon [24], and carrot [25]. For spinach, GWAS has been used to identify markers and genes associated with features such as disease resistance [26,27,28,29,30], sex determination [31], oxalate concentration [15], and mineral element content [14]. However, few GWASs have been reported on the comprehensive nutritional quality of spinach. In this study, the concentrations of chlorophyll a, chlorophyll b, oxalate, nitrate, crude fiber, soluble sugar, manganese, copper, and iron were determined in 62 inbred lines of spinach. These spinach lines were sequenced with around 10-fold genomic coverage. Using the resulting data, the GWAS for nine nutritional traits identified a substantial number of loci, including key genes regulating biosynthesis and the transportation of mineral elements, indicating that the GWAS of spinach could provide an effective approach for gene identification, pathway elucidation, and spinach improvement.

2. Materials and Methods

2.1. Plant Materials

A total of 62 spinach inbred lines were used in this experiment which were obtained from the Institute of Vegetables and Flowers at the Chinese Academy of Agriculture Sciences (CAAS). All accessions were planted in an open field in Beijing (40° N, 116° E) in the spring of 2019, with the same management practices as conventional cultivation. Young leaves were collected 25 days after sowing for resequencing. For nutrient content determination, 10–20 plants with uniform growth were selected from each accession 50 days after sowing.

2.2. Extraction and Determination of Nutrients

The extraction method was performed using the same protocol as before [32]. The chlorophyll content was determined by spectrophotometry and the organic acid content was determined using ion chromatography. Acid-base treatment was used to measure the crude fiber content. The soluble sugar content was determined using the copper reduction iodometric method. The concentration of mineral elements was determined by inductively coupled plasma emission spectrometry.

2.3. Statistical Analysis of Nutrient Content

IBM SPSS Statistics 23 software was used to calculate the mean value, standard deviation, and coefficient of variation of each component. Origin Pro 2022 software was used to determine the frequency distribution and correlation analysis and to generate figures.

2.4. DNA Extraction

Genomic DNA was extracted from the spinach leaves using the cetyltrimethylammonium bromide (CTAB) method [33]. The quality and concentration of DNA were detected using 1.0% agarose gel electrophoresis and an ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, NC, USA), respectively. The DNA solution was diluted to 20–100 ng/µL as a working solution and stored at −20 °C for subsequent testing.

2.5. Resequencing and Genotyping

After determining the extraction quality, the DNA was sent to the Berry Genomics Company (Beijing, China) for resequencing using the Illumina platform. The raw data were filtered using fastp software (v0.20.0) (parameter q 30, the rest were default) [34] to obtain the clean data. The clean reads were then aligned to the spinach reference genome Monoe-Viroflay [35] using Burrows–Wheeler Aligner (BWA) software (v0.7.17-r1188) (http://bio-bwa.sourceforge.net/ (accessed on 14 July 2023)). SAMtools software (v1.6-3-g200708f) [36] was used to sort the data, mark the duplicate sequences, and convert the alignment file into a binary alignment map (BAM) file. Single-nucleotide polymorphism (SNP) calling was performed using the genome analysis toolkit software (GATK4.3.0.0) [37]. Finally, raw SNPs were filtered using VCFtools software (v0.1.15) [38] based on the following criteria: deletion call rate > 0.8, minor allele frequency (MAF) < 0.05, loci with only two alleles to be retained, and information in all INFO fields to be retained without filtering. The high-quality SNPs obtained through these filtering criteria were used for GWAS analysis.
The functional annotation of variants was performed using ANNOVAR software (version: 2017-07-17) [39].

2.6. Linkage Disequilibrium Analysis

PopLDdecay software (v3.41) [40] was used to determine the linkage disequilibrium (LD) coefficient (r2) values for each pair of SNPs, perform genome-wide LD analysis, and visualize the results. LDBlockShow software (v1.40) [41] was utilized to identify LD blocks within the candidate regions.

2.7. Genome-Wide Association Analysis

The GWAS was performed using two statistical models in rMVP (https://github.com/xiaolei-lab/rMVP (accessed on 20 July 2023)): Generalized Linear Model (GLM) and Fixed and Random Model Circulating Probability Unification (FarmCPU) [42]. The kinship matrix (K) was calculated using the default method. Significance thresholds (−log10P) were corrected using Bonferroni correction. In this study, −log10P = 6.52 (P = 1/3356182) was used as the threshold value to determine the significance of SNP markers for loci associated with nutrient concentrations. The Manhattan plot was plotted using rMVp to visualize the GWAS analysis.

2.8. Candidate Gene Identification

Using Monoe-Viroflay as the reference genome, genes within the candidate region were extracted. The Blast Zone of TBtools II (https://github.com/CJ-Chen/TBtools-II/releases (accessed on 24 September 2023)) [43] was used to identify the corresponding homologous genes in Arabidopsis thaliana, and the one with the best score was taken as the homologous gene. The functions of the candidate genes were analyzed by combining the annotation information of the spinach genome and A. thaliana homologous genes.

3. Results

3.1. Statistical Analysis of Phenotyping Results

In this study, the content of nine nutrients was investigated in different spinach accessions to reveal the change rule of the nutritional quality of spinach (Table 1 and Table S1). In terms of chlorophyll, the results showed that chlorophyll a and chlorophyll b in 62 spinach accessions were significantly correlated, and the correlation coefficient was as high as 0.95 (Figure 1). Specifically, the average chlorophyll a content was 0.749 mg/g, with a variance of 0.480–1.187 mg/g and a coefficient of variation of 23.541%. The average content of chlorophyll b was 0.308 mg/g, with values ranging from 0.170 to 0.466 mg/g, and the coefficient of variation was 26.739%. There was a significant variation in oxalate content among the test accessions. The average concentration was 4.614 g/kg, with a range of variation from 1.86 to 8.89 g/kg and a coefficient of variation of 29.509%. For nitrate, the average content was 0.996 g/kg, with a variance of 0.26 to 2.24 g/kg, exhibiting a relatively high coefficient of variation at 57.716%. Additionally, nitrate showed a significant negative correlation with chlorophyll a, chlorophyll b, and soluble sugars (Figure 1). For soluble sugars, the average content was 0.647%, with a variation range of 0.14–1.64%, and the coefficient of variation was 64.849%. For crude fiber, the average content was 0.728%, and it had the lowest coefficient of variation among all qualities at 8.714%. For mineral elements, the Fe content in spinach was much higher than that of Mn and Cu. The average Fe content was 135.524 mg/kg, and the average content of Mn and Cu was 6.560 and 1.235 mg/kg, respectively. Moreover, Fe and Mn showed a significant positive correlation, and the correlation coefficient was 0.69 (Figure 1). The content of all nine nutrients showed a normal distribution or an approximately normal distribution (Figure S1), aligning with the quantitative trait characteristics. This indicates that the content of these qualities is influenced or controlled by multiple genes.

3.2. Resequencing of Spinach Varieties

Resequencing was conducted on 62 spinach samples, resulting in a total of 787.7 GB of raw data. After filtering with the fastp software, 696.4 GB of clean data was retained. In alignment with the reference genome, Monoe-Viroflay identified a total of 7,728,915 SNPs. Following that, further filtering was performed to obtain 3,356,182 high-quality SNPs. Among these variations, 3,347,533 SNPs were anchored to chromosomes, with the majority located on chromosome 4 (734,008), while chromosome 2 had the fewest distributed SNPs (439,501) (Figure 2 and Table S2).

3.3. Genome-Wide Linkage Disequilibrium Analysis

LD refers to the non-random association of alleles at different loci [44], reflecting the population’s recombination strength. It determines the marker density required and the accuracy of association analysis. During the genome-wide LD decay analysis of high-quality SNPs, when the r2 value decreased from its maximum of 0.7 to half of the maximum value (0.35), the corresponding physical distance was 10.47 kb (Figure 3). Therefore, the LD decay value for the entire spinach genome was 10.47 kb, which was faster than that in another study (5 kb) [45].

3.4. Genome-Wide Association Analysis of Nutrients

Intending to identify the genetic loci controlling these traits, we utilized 3,356,182 SNPs based on two models (GLM and FarmCPU) in rMVP to conduct a GWAS for nine nutrient concentrations in spinach. The GWAS significance threshold was determined to be −log10(p) = 6.52 using the Bonferroni correction. Based on GLM and FarmCPU, we found 1996 and 104 linked loci, respectively, with 23 loci recognized by both models (Table 2 and Table S3; Figure 4, Figures S2 and S3). All in all, there were 2077 loci with −log10(p) > 6.52, mainly distributed on chromosomes 1, 2, and 4. As for chlorophyll a and chlorophyll b, the peak signal was located at the same locus, SOVchr4-65150089; in addition, they co-localized on chromosomes 1, 4, and 6. For oxalate, nitrate, and Fe, their lead SNPs were distributed at 50,807,332, 64,960,432, and 43,602,473 bp on chromosome 1, respectively. Strong association signals were located on chromosome 2 for soluble sugar and Cu and on chromosome 5 for Mn. However, no associated loci were identified for crude fiber.

3.5. Candidate Genes

To find candidate genes, we conducted a haplotype analysis of SNPs throughout the whole genome using LDBlockShow. The region associated with the significant locus was identified as a candidate interval. Subsequently, we screened candidate genes related to the target traits based on the annotation information of the spinach genome Monoe-Viroflay. Among the candidate regions, we identified 127 genes, including 19 unknown genes and 108 known annotation genes, of which 85 had homologous genes in A. thaliana (Table S4).
Regarding chlorophyll, within its candidate regions, we identified a total of 62 predicted genes. Among them, 50 genes were connected with chlorophyll a content, 8 with chlorophyll b, and 4 with both. Among the 62 genes, based on gene annotation, 4 genes were discovered to be associated with chlorophyll content. SOV1g020220 was located on chromosome 1: 93,828,109–93,830,147 bp, and was associated with both chlorophyll a and chlorophyll b. Within its coding region, we identified four non-synonymous SNPs, leading to amino acid changes (Table S5). Additionally, it served as a homologous gene to ATFRD3, which encoded a multidrug and toxin efflux (MATE) family member and played a critical role in regulating cellular iron ion balance; another study suggested that iron plays a role in the biosynthesis of chlorophyll [4]. SOV4g012830 was found in the candidate region of chlorophyll a, was homologous to ATFYH3 in A. thaliana, and was a component of the PHYTOCHROME A (PHYA) signaling network, mediating the Far-Red High-Irradiance Response (FR-HIR) in collaboration with FAR-RED IMPAIRED RESPONSE1 (FAR1). SOV4g012890 encoded a FAR1 domain-containing protein, and two non-synonymous SNPs were identified in its coding region (Table S5). The last gene was SOV4g004990, correlated with chlorophyll b concentration, and its homologous gene AT2G01050 regulated zinc ion binding in A. thaliana. Similarly, zinc ions have been demonstrated to participate in the synthesis of chlorophyll [46].
For oxalate, we found a total of 12 genes associated with its content, among which SOV1g009820 was considered a key candidate gene. SOV1g009820, which encodes β-glucuronosyltransferase GlcAT14A, is found on SOVchr1 at 50,879,476–50,883,262 bp. Furthermore, its coding region exhibited extensive variation, including 11 non-synonymous SNPs, 3 synonymous SNPs, 1 stop-gain SNP, two frameshift insertions, and one non-frameshift deletion. This stop-gain caused the gene to terminate prematurely, resulting in a protein with just 215 amino acids (Table S5).
For soluble sugar, within the range of SOVchr2: 108,223,175–108,323,175 bp, we identified eight genes. Among them, SOV2g035100 was located 12 kb upstream of the SOVchr2-108273175 locus and was thought to be an essential gene that regulates soluble sugar concentration. SOV2g035100 was homologous to ATSTP10, which encoded a hexose-H (+) symporter in A. thaliana.
With regard to Fe, we discovered three genes in their LD blocks, two of which were located on chromosome 1 and one was located on chromosome 5. SOV1g008480 was considered to be a potential gene affecting iron content in spinach. It was located within the range of SOVchr1: 43,699,103–43,701,699 bp and was homologous to ATARD1. It encodes an acireductone dioxygenase that is activated by contact with a heterotrimeric G protein β subunit, acting as a metal enzyme that may regulate iron ion binding.

4. Discussion

With the improvement of people’s living standards, growing attention has gradually turned to health and nutrition. Spinach, as a representation of comprehensive nutrition, has progressively gained recognition. As a result, researchers are concerned not only with the yield of spinach but also with increasing the concentrations of beneficial nutrients while decreasing the concentrations of the harmful components of spinach. To better understand the genetic mechanism, several efforts have been made to assess dietary features using various methodologies [32].
In order to identify genes regulating nutrient components in spinach, our study used individual nutrient content and resequencing data, conducting a GWAS analysis of nine nutritional elements of 62 spinach accessions. GWAS has been employed for a variety of crops to uncover linked loci or genes for nutritional elements such as soluble sugars [47] and mineral elements [48]. This shows that GWAS is appropriate for studying nutritional elements. Moreover, past studies have shown that GWAS is not only successful in varied sample investigations but also produces important results, even in studies with less than 100 samples [49]. As a consequence, the number of accessions has no bearing on the credibility of this study’s findings. Spinach is a hybrid species with a high effective recombination rate, which results in a substantial number of heterozygotes and a low LD. Our calculated LD decay distance for the whole spinach genome was 10.47 Kb, with an r-squared value of 0.35. This decay distance is remarkably shorter than that of other crops, such as snap bean [50] and rice [51]. Rapid LD decay in spinach is advantageous for determining the correlation between linked loci, thereby reducing the risk of false positives in candidate gene identification during GWAS.
Chlorophyll is a natural pigment found in plants that plays an important role in the photosynthesis process. It has a variety of activities, including anticancer, antibacterial, and antioxidant properties. Additionally, chlorophyll plays a protective role in photosynthesis by shielding plants from the damage caused by light during complex cellular processes [52]. In plants, chlorophyll a can synthesize chlorophyll b, which may then be transformed back into chlorophyll a [53]. In previous studies, a gene in Chlamydomonas called CAO (chlorophyll a oxygenase) has been discovered that codes for an enzyme that oxidizes chlorophyll a to produce chlorophyll b [54]. This study discovered four chlorophyll-related candidate genes. Among them, SpFRD3 (SOV1g020220) is a homolog of ATFRD3, which is implicated in iron ion transport and balance in A. thaliana [55]. Previous research has indicated that iron, as a redox-active metal, is involved in photosynthesis and the formation of chlorophyll in plants [4]. Another potential gene, SOV4g004990, is homologous to AT2G01050 and is involved in binding zinc ions. Zinc ions, which act as cofactors, are essential components of numerous enzymes and participate in chlorophyll synthesis [46]. Therefore, we suggest that SpFRD3 and SOV4g004990 may influence the chlorophyll content in spinach by regulating the transport or binding of metal ions within the plant. SpFHY3 (SOV4g012830) is homologous to ATFHY3 (AT3G22170), which is a component of the PHYA signaling network, mediating the FR-HIR response to far-red light in conjunction with FAR1. In A. thaliana, ATFHY3 positively regulates chlorophyll biosynthesis by activating the expression of the HEMB1 gene (HEMB1 encodes 5-aminolevulinic acid dehydratase in the chlorophyll biosynthetic pathway) [56]. Additionally, SOV4g012890 encodes a FAR1 domain-containing protein. The Far1 nuclear protein is derived from the plant MuDR transposase, involved in phytochrome signaling [57], and it has been reported that photosensitive pigments enter the nucleus of the cell when exposed to light and interact directly with transcription factors to regulate gene expression and indirectly affect chlorophyll biosynthesis [58,59]. Therefore, we hypothesize that SpFHY3 and SOV4g012890 may influence the chlorophyll content in spinach by regulating the transmission of photosensitive pigment signals within the plant. In future research, we will use virus-induced gene silencing (VIGS) to confirm the functions of these two genes.
In spinach, the occurrence of a substantial amount of oxalate has an influence on both taste and health. The distribution of oxalate is affected by a variety of factors, including the environment, season, growing period, and genetics. In spinach, researchers have observed higher oxalate concentrations in older leaves [60]; the oxalate content in spinach petioles is lower than that in leaves, thus recommending the selection and breeding of cultivars with high leaf/petiole ratios [61]. Fall-planted spinach has the lowest oxalate content, followed by spring- and summer-planted spinach, while winter-planted spinach has the highest. [62]. Previous studies have identified two genes on chromosome 5 encoding metal ion transporter proteins that affect the levels of soluble oxalate in spinach by influencing the levels of calcium and other ions in spinach cells, thereby affecting the levels of soluble oxalate in spinach [16]. Despite the fact that spinach oxalate has been extensively studied by other researchers, the identified loci or genes have not been validated, and the results of previous studies are not consistence, indicating that oxalate is controlled by multilocus genes and is highly influenced by genetic background and environmental differences. The study identified 12 genes connected to oxalate, and SOV1g009820 encodes β-glucuronosyltransferase, which transfers glucuronic acid to Arabinogalactan-proteins (AGPs). Additionally, glucuronic acid can bind to calcium ions [63]. Moreover, there is one stop-gain SNP in the coding region of this gene that contributes to the premature termination of the expression of this gene and may decrease the activity of glucuronosyltransferase, which, in turn, increases the binding of glucuronide and calcium ions in spinach and decreases the binding of oxalate and calcium ions, leading to a higher soluble oxalate content in spinach. Therefore, we consider this gene as an important candidate gene that is related to the oxalate content of spinach.
Soluble sugars play a crucial role in regulating the taste of fruits and vegetables [17]. They provide energy and structural material and are significant regulators of plant growth, development, and gene expression [64]. It has been demonstrated that soluble sugar synthesis and transporter genes are associated with their concentration accumulation in crops. For example, in Chinese cabbage, there are six genes associated with soluble sugar content, BraA05gAOP1, BraA04gAOP4, BraA03gHT7, and BraA01gHT4 (hexose transporter protein), which are structural genes in soluble sugar biosynthesis, and BraA01gCHR11 and BraA07gSCL1, which are two important transcription factors regulating soluble sugar biosynthesis [65]. There are few reports on loci or genes related to soluble sugar content in spinach. Among the four genes identified in this study, the A. thaliana homolog of SpSTP10 (SOV2g035100), ATSTP10, encodes a hexose-H (+) symporter and catalyzes the high-affinity uptake of glucose, galactose, and mannose [66]. This is consistent with earlier studies that suggested that this gene may be involved in the biosynthesis of soluble sugars.
Iron plays an important role in both humans and plants. First, in the human body, iron is a component of hemoglobin, myoglobin, cytochromes, and some respiratory enzymes, and is involved in oxygen and carbon dioxide transport and exchange and tissue respiration in the body [67]. Second, for plants, iron is a catalyst for chlorophyll synthesis, participates in photosynthesis, and affects root development and overall growth [68]. This study found SpARD1 (SOV1g008480) to be a putative iron-related gene. The homologous gene of this gene is ATARD1 (AT4G14716), a metalloenzyme that helps A. thaliana bind iron ions and acireductone dioxygenase (ARD). Within the iron ion-related pathway in A. thaliana, ATARD1 catalyzes the conversion of acid-reduced ketones to formate and 2-keto-4-methylthiobutyrate as part of the iron ion-related pathway [69]. Thus, based on findings from previous research, we hypothesize that SpARD1 may play a similar role in spinach, participating in the binding of iron ions. We will apply the genetic transformation to the candidate gene in order to assess it more thoroughly.
Molecular marker-assisted selection (MAS) breeding makes use of the fact that molecular markers are closely linked to the genes controlling the target traits and selects excellent individuals with the target traits by detecting molecular markers in the progeny. Compared with traditional breeding methods, this method is faster, more precise, and less affected by environmental circumstances. It may be employed in a number of breeding processes, including identifying parentage, transmitting quantitative and recessive traits in backcrossing, selecting hybrid progeny, predicting hybrid advantage, and determining variety purity [70]. Except for crude fiber, this analysis discovered 2077 SNP loci connected with the remaining eight characteristics. Particularly for Fe, we have identified two associated loci, located at 43,602,473 bp on chromosome 1 and 91,913,409 bp on chromosome 5. These loci accounted for 42.03% and 48.63% of the phenotypic variance, respectively. Both SNPs demonstrated unusually high phenotypic variability and might be utilized to develop molecular markers to aid in the selection of Fe-rich spinach types. In conclusion, the SNP loci associated with nutrient elements identified in this study can be developed as molecular markers, and it is possible to polymerize these markers to screen for superior spinach varieties that are rich in chlorophyll, soluble sugars, and mineral elements, as well as containing small amounts of oxalic acid and nitrate in spinach germplasm.

5. Conclusions

This study identified 2077 SNP loci associated with the concentrations of chlorophyll a, chlorophyll b, oxalate, nitrate, soluble sugars, Fe, Mn, and Cu in spinach. We predicted seven candidate genes based on these loci. Among them, SpFRD3, SOV4g004990, SpFHY3, and SOV4g012890 were considered to be associated with chlorophyll. They influence the chlorophyll content in spinach through the regulation of iron ion balance, zinc ion binding, and photosensitive pigment signal transduction. Additionally, SOV1g009820 encodes β-glucuronosyltransferase, SpSTP10 encodes a hexose-H (+) symporter, and SpARD1 is involved in iron ion binding. These genes are associated with oxalate, soluble sugar, and iron content, respectively. Therefore, this study provides SNP locus information and potential candidate genes for enhancing spinach’s positive traits and diminishing its negative ones. These SNP loci should be used in subsequent research to provide molecular markers for molecular-assisted breeding. Additionally, further research is needed to explore these genes and gain a deeper understanding of their molecular mechanisms in altering the nutritional quality of spinach. Such investigations contribute to a better understanding of the regulatory mechanisms of plant nutritional quality, offering a scientific basis for optimizing spinach’s nutritional characteristics.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes15020172/s1. Figure S1: Frequency distribution of different nutrient elements in 62 spinach accessions; Figure S2: Manhattan plot for genome-wide association analysis (GWAS) of nutrient (a: chlorophyll a; b: chlorophyll b; c: nitrate; d: soluble sugar; e: Mn; f: Cu; and g: Fe) content in spinach using GLM and FarmCPU models. The horizontal and vertical axes represent the genomic position of the SNP and the association power for each SNP, respectively, with the trait expressed as −log10(p). The dashed line shows the Bonferroni-corrected genome-wide threshold; Figure S3: Map of the distribution of SNPs associated with nutrient content on each chromosome. Table S1: Comparative analysis of nutritional quality of different spinach accessions; Table S2: SNP distribution on chromosomes; Table S3: Information on SNPs associated with the nutrient elements detected using the GLM model; Table S4: Putative genes associated with the nutritional composition of spinach within candidate regions based on associated SNPs; Table S5: Functional annotation information for exonic regions of candidate genes.

Author Contributions

Conceptualization, W.Q. and Z.F.; methodology, Z.L., H.S. and N.J.; software, Z.L. and N.J.; validation, Z.L., H.S. and N.J.; formal analysis, H.S. and N.J.; investigation, H.S. and N.J.; resources, Z.X. and H.Z.; data curation, N.J.; writing—original draft preparation, N.J. and Z.L.; writing—review and editing, Z.L.; visualization, Z.L. and N.J.; supervision, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Chinese Academy of Agricultural Sciences Innovation Project (CAAS-ASTIP-IVFCAAS, CAAS-ZDRW202103), China Agricultural Research System (CARS-23-A-17), National Key R&D Program of China (2023YFD1200102), and Beijing Joint Research Program for Germplasm Innovation and New Variety Breeding (G20220628003-03).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. van Treuren, R.; Coquin, P.; Lohwasser, U. Genetic resources collections of leafy vegetables (lettuce, spinach, chicory, artichoke, asparagus, lamb’s lettuce, rhubarb and rocket salad): Composition and gaps. Genet. Resour. Crop Evol. 2012, 59, 981–997. [Google Scholar] [CrossRef]
  2. Lester, G.E.; Makus, D.J.; Hodges, D.M.; Jifon, J.L. Summer (subarctic) versus winter (subtropic) production affects spinach (Spinacia oleracea L.) leaf bionutrients: Vitamins (C, E, Folate, K1, provitamin A), lutein, phenolics, and antioxidants. J. Agric. Food Chem. 2013, 61, 7019–7027. [Google Scholar] [CrossRef]
  3. Fletcher, R.; Bell, I.; Lambert, J. Public health aspects of food fortification: A question of balance. Proc. Nutr. Soc. 2004, 63, 605–614. [Google Scholar] [CrossRef]
  4. Soetan, K.O.; Olaiya, C.O.; Oyewole, O.E. The importance of mineral elements for humans, domestic animals and plants: A review. Afr. J. Food Sci. 2010, 4, 200–222. [Google Scholar] [CrossRef]
  5. Ye, S.; Shah, B.R.; Li, J.; Liang, H.; Zhan, F.; Geng, F.; Li, B. A critical review on interplay between dietary fibers and gut microbiota. Trends Food Sci. Technol. 2022, 124, 237–249. [Google Scholar] [CrossRef]
  6. Zhu, N.; Zhao, C.; Wei, Y.; Sun, C.; Wu, D.; Chen, K. Biosynthetic labeling with 3-O-propargylcaffeyl alcohol reveals in vivo cell-specific patterned lignification in loquat fruits during development and postharvest storage. Hortic. Res. 2021, 8, 61. [Google Scholar] [CrossRef]
  7. Mou, B. Evaluation of oxalate concentration in the U.S. spinach germplasm collection. HortScience 2008, 43, 1690–1693. [Google Scholar] [CrossRef]
  8. Bong, W.C.; Vanhanen, L.P.; Savage, G.P. Addition of calcium compounds to reduce soluble oxalate in a high oxalate food system. Food Chem. 2017, 221, 54–57. [Google Scholar] [CrossRef] [PubMed]
  9. Franceschi, V.R.; Nakata, P.A. Calcium oxalate in plants: Formation and function. Annu. Rev. Plant Biol. 2005, 56, 41–71. [Google Scholar] [CrossRef] [PubMed]
  10. Noonan, S.C.; Savage, G.P. Oxalate content of foods and its effect on humans. Asia Pac. J. Clin. Nutr. 1999, 8, 64–74. [Google Scholar] [CrossRef]
  11. Breimer, T. Environmental factors and cultural measures affecting the nitrate content in spinach. Fertil. Res. 1982, 3, 191–292. [Google Scholar] [CrossRef]
  12. Kawazu, Y.; Okimura, M.; Ishii, T.; Yui, S. Varietal and seasonal differences in oxalate content of spinach. Sci. Hortic. 2003, 97, 203–210. [Google Scholar] [CrossRef]
  13. Lee, M.-H.; Han, J.-S.; Kozukue, N. Changes of chlorophyll contents in spinach by growth periods and storage. Korean J. Food Cook. Sci. 2005, 21, 339–345. [Google Scholar]
  14. Qin, J.; Shi, A.; Mou, B.; Grusak, M.A.; Weng, Y.; Ravelombola, W.; Bhattarai, G.; Dong, L.; Yang, W. Genetic diversity and association mapping of mineral element concentrations in spinach leaves. BMC Genom. 2017, 18, 941. [Google Scholar] [CrossRef] [PubMed]
  15. Shi, A.; Mou, B.; Correll, J.C. Association analysis for oxalate concentration in spinach. Euphytica 2016, 212, 17–28. [Google Scholar] [CrossRef]
  16. Cai, X.; Sun, X.; Xu, C.; Sun, H.; Wang, X.; Ge, C.; Zhang, Z.; Wang, Q.; Fei, Z.; Jiao, C.; et al. Genomic analyses provide insights into spinach domestication and the genetic basis of agronomic traits. Nat. Commun. 2021, 12, 7246. [Google Scholar] [CrossRef] [PubMed]
  17. Wu, X.; Chen, F.; Zhao, X.; Pang, C.; Shi, R.; Liu, C.; Sun, C.; Zhang, W.; Wang, X.; Zhang, J. QTL mapping and GWAS reveal the genetic mechanism controlling soluble solids content in Brassica napus shoots. Foods 2021, 10, 2400. [Google Scholar] [CrossRef] [PubMed]
  18. Huang, X.; Wei, X.; Sang, T.; Zhao, Q.; Feng, Q.; Zhao, Y.; Li, C.; Zhu, C.; Lu, T.; Zhang, Z.; et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 2010, 42, 961–967. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, X.; Guan, Z.; Wang, L.; Fu, J.; Zhang, Y.; Li, Z.; Ma, L.; Liu, P.; Zhang, Y.; Liu, M.; et al. Combined GWAS and QTL analysis for dissecting the genetic architecture of kernel test weight in maize. Mol. Genet. Genom. 2020, 295, 409–420. [Google Scholar] [CrossRef] [PubMed]
  20. Li, M.; Wu, X.; Wang, B.; Wu, X.; Wang, Y.; Wang, J.; Dong, J.; Wu, J.; Lu, Z.; Sun, Y.; et al. Genome-wide association analysis reveals the optimal genomic regions for pod size in bean. Front. Plant Sci. 2023, 14, 1138988. [Google Scholar] [CrossRef]
  21. Zhao, J.; Sauvage, C.; Zhao, J.; Bitton, F.; Bauchet, G.; Liu, D.; Huang, S.; Tieman, D.M.; Klee, H.J.; Causse, M. Meta-analysis of genome-wide association studies provides insights into genetic control of tomato flavor. Nat. Commun. 2019, 10, 1534. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, Y.; Dong, S.; Miao, H.; Liu, X.; Dai, Z.; Li, X.; Gu, X.; Zhang, S. Genome-wide association studies reveal candidate genes related to stem diameter in cucumber (Cucumis sativus L.). Genes 2022, 13, 1095. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, L.; Wang, P.; Wang, Y.; Cheng, Q.; Lu, Q.; Liu, J.; Li, T.; Ai, Y.; Yang, W.; Sun, L.; et al. Genome-wide correlation of 36 agronomic traits in the 287 pepper (Capsicum) accessions obtained from the SLAF-seq-based GWAS. Int. J. Mol. Sci. 2019, 20, 5675. [Google Scholar] [CrossRef] [PubMed]
  24. Yan, H.; Wang, K.; Wang, M.; Feng, L.; Zhang, H.; Wei, X. QTL mapping and genome-wide association study reveal genetic loci and candidate genes related to soluble solids content in melon. Curr. Issues Mol. Biol. 2023, 45, 7110–7129. [Google Scholar] [CrossRef] [PubMed]
  25. Brainard, S.H.; Ellison, S.L.; Simon, P.W.; Dawson, J.C.; Goldman, I.L. Genetic characterization of carrot root shape and size using genome-wide association analysis and genomic-estimated breeding values. Theor. Appl. Genet. 2022, 135, 605–622. [Google Scholar] [CrossRef] [PubMed]
  26. Awika, H.; Cochran, K.; Joshi, V.; Bedre, R.; Mandadi, K.; Avila, C. Single-marker and haplotype-based association analysis of anthracnose (Colletotrichum dematium) resistance in spinach (Spinacia oleracea). Plant Breed. 2019, 139, 402–418. [Google Scholar] [CrossRef]
  27. Awika, H.O.; Marconi, T.G.; Bedre, R.; Mandadi, K.K.; Avila, C.A. Minor alleles are associated with white rust (Albugo occidentalis) susceptibility in spinach (Spinacia oleracea). Hortic. Res. 2019, 6, 129. [Google Scholar] [CrossRef]
  28. Bhattarai, G.; Yang, W.; Shi, A.; Feng, C.; Dhillon, B.; Correll, J.C.; Mou, B. High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach. BMC Genom. 2021, 22, 478. [Google Scholar] [CrossRef]
  29. Shi, A.; Bhattarai, G.; Xiong, H.; Avila, C.A.; Feng, C.; Liu, B.; Joshi, V.; Stein, L.; Mou, B.; du Toit, L.J.; et al. Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm. Hortic. Res. 2022, 9, uhac069. [Google Scholar] [CrossRef]
  30. Shi, A.; Mou, B.; Correll, J.C.; Motes, D.; Weng, Y.; Qin, J.; Yang, W. SNP association analysis of resistance to Verticillium wilt (‘Verticillium dahliae’ Kleb.) in spinach. Aust. J. Crop Sci. 2016, 10, 1188–1196. [Google Scholar] [CrossRef]
  31. Wang, H.; Liu, Z.; Wang, X.; Wu, J.; Zhang, H.; Xia, Z.; Xu, Z.; Qian, W. Genome-wide association studies for monoecism in spinach. Acta Hortic. Sin. 2019, 46, 1495–1502. [Google Scholar] [CrossRef]
  32. She, H.; Liu, Z.; Xu, Z.; Zhang, H.; Qian, W. Nutritional quality analysis and evaluation of different spinach accessions. J. Chang. Veg. 2021, 24, 45–51. [Google Scholar] [CrossRef]
  33. Murray, M.G.; Thompson, W.F. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  35. Collins, K.; Zhao, K.; Jiao, C.; Xu, C.; Cai, X.; Wang, X.; Ge, C.; Dai, S.; Wang, Q.; Wang, Q.; et al. SpinachBase: A central portal for spinach genomics. Database J. Biol. Databases Curation 2019, 2019, baz072. [Google Scholar] [CrossRef] [PubMed]
  36. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  37. Van der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ data to high confidence variant calls: The Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  38. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef]
  39. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  40. 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, 1786–1788. [Google Scholar] [CrossRef]
  41. Dong, S.-S.; He, W.-M.; Ji, J.-J.; Zhang, C.; Guo, Y.; Yang, T.-L. LDBlockShow: A fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Brief. Bioinform. 2020, 22, bbaa227. [Google Scholar] [CrossRef] [PubMed]
  42. Yin, L.; Zhang, H.; Tang, Z.; Xu, J.; Yin, D.; Zhang, Z.; Yuan, X.; Zhu, M.; Zhao, S.; Li, X.; et al. rMVP: A memory-efficient, visualization-enhanced, and parallel-accelerated tool for genome-wide association study. Genom. Proteom. Bioinform. 2021, 19, 619–628. [Google Scholar] [CrossRef] [PubMed]
  43. Chen, C.; Wu, Y.; Li, J.; Wang, X.; Zeng, Z.; Xu, J.; Liu, Y.; Feng, J.; Chen, H.; He, Y.; et al. TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Mol. Plant 2023, 16, 1733–1742. [Google Scholar] [CrossRef] [PubMed]
  44. Kim, S.; Plagnol, V.; Hu, T.T.; Toomajian, C.; Clark, R.M.; Ossowski, S.; Ecker, J.R.; Weigel, D.; Nordborg, M. Recombination and linkage disequilibrium in Arabidopsis thaliana. Nat. Genet. 2007, 39, 1151–1155. [Google Scholar] [CrossRef] [PubMed]
  45. Xu, C.; Jiao, C.; Sun, H.; Cai, X.; Wang, X.; Ge, C.; Zheng, Y.; Liu, W.; Sun, X.; Xu, Y.; et al. Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions. Nat. Commun. 2017, 8, 15275. [Google Scholar] [CrossRef]
  46. Sadeghzadeh, B. A review of zinc nutrition and plant breeding. J. Soil Sci. Plant Nutr. 2013, 13, 905–927. [Google Scholar] [CrossRef]
  47. Xu, W.; Liu, H.; Li, S.; Zhang, W.; Wang, Q.; Zhang, H.; Liu, X.; Cui, X.; Chen, X.; Tang, W. GWAS and identification of candidate genes associated with seed soluble sugar content in vegetable soybean. Agronomy 2022, 12, 1470. [Google Scholar] [CrossRef]
  48. Gyawali, S.; Otte, M.L.; Chao, S.; Jilal, A.; Jacob, D.L.; Amezrou, R.; Verma, R.P.S. Genome wide association studies (GWAS) of element contents in grain with a special focus on zinc and iron in a world collection of barley (Hordeum vulgare L.). J. Cereal Sci. 2017, 77, 266–274. [Google Scholar] [CrossRef]
  49. Korte, A.; Farlow, A. The advantages and limitations of trait analysis with GWAS: A review. Plant Methods 2013, 9, 29. [Google Scholar] [CrossRef]
  50. Liu, Z.; Gao, S.; Zhang, H.; Xu, Z.; Qian, W. Genome-wide association study reveals that PvGUX1_1 is associated with pod stringlessness in snap bean (Phaseolus vulgaris L.). Biology 2022, 11, 611. [Google Scholar] [CrossRef]
  51. Mather, K.A.; Caicedo, A.L.; Polato, N.R.; Olsen, K.M.; McCouch, S.; Purugganan, M.D. The extent of linkage disequilibrium in rice (Oryza sativa L.). Genetics 2007, 177, 2223–2232. [Google Scholar] [CrossRef] [PubMed]
  52. Willows, R.D. Chlorophylls. Plant Pigment. Their Manip. 2004, 14, 23–56. [Google Scholar]
  53. Rüdiger, W. Biosynthesis of chlorophyll b and the chlorophyll cycle. Photosynth. Res. 2002, 74, 187–193. [Google Scholar] [CrossRef] [PubMed]
  54. Tanaka, A.; Ito, H.; Tanaka, R.; Tanaka, N.K.; Yoshida, K.; Okada, K. Chlorophyll a oxygenase (CAO) is involved in chlorophyll b formation from chlorophyll a. Proc. Natl. Acad. Sci. USA 1998, 95, 12719–12723. [Google Scholar] [CrossRef]
  55. Rogers, E.E.; Guerinot, M.L. FRD3, a member of the multidrug and toxin efflux family, controls iron deficiency responses in Arabidopsis. Plant Cell 2002, 14, 1787–1799. [Google Scholar] [CrossRef]
  56. Tang, W.; Wang, W.; Chen, D.; Ji, Q.; Jing, Y.; Wang, H.; Lin, R. Transposase-derived proteins FHY3/FAR1 interact with PHYTOCHROME-INTERACTING FACTOR1 to regulate chlorophyll biosynthesis by modulating HEMB1 during deetiolation in Arabidopsis. Plant Cell 2012, 24, 1984–2000. [Google Scholar] [CrossRef]
  57. Makarova, K.S.; Aravind, L.; Koonin, E.V. SWIM, a novel Zn-chelating domain present in bacteria, archaea and eukaryotes. Trends Biochem. Sci. 2002, 27, 384–386. [Google Scholar] [CrossRef] [PubMed]
  58. Davies, T.G.; Thomas, H.; Thomas, B.J.; Rogers, L.J. Leaf senescence in a nonyellowing mutant of festuca pratensis: Metabolism of cytochrome f. Plant Physiol. 1990, 93, 588–595. [Google Scholar] [CrossRef]
  59. López-Juez, E.; Paul Jarvis, R.; Takeuchi, A.; Page, A.M.; Chory, J. New Arabidopsis cue mutants suggest a close connection between plastid-and phytochrome regulation of nuclear gene expression. Plant Physiol. 1998, 118, 803–815. [Google Scholar] [CrossRef]
  60. Okutani, I.; Sugiyama, N. Relationship between oxalate concentration and leaf position in various spinach cultivars. HortScience 1994, 29, 1019–1021. [Google Scholar] [CrossRef]
  61. Oguchi, Y.; Weerakkody, W.A.P.; Tanaka, A.; Nakazawa, S.; Ando, T. Varietal differences of quality-related compounds in leaves and petioles of spinach grown at two locations. Bull. Hiroshima Prefect. Agric. Res. Cent. 1996, 64, 1–9. [Google Scholar]
  62. Kaminishi, A.; Kita, N. Seasonal change of nitrate and oxalate concentration in relation to the growth rate of spinach cultivars. HortScience 2006, 41, 1589–1595. [Google Scholar] [CrossRef]
  63. Ajayi, O.O.; Held, M.A.; Showalter, A.M. Three β-glucuronosyltransferase genes involved in arabinogalactan biosynthesis function in arabidopsis growth and development. Plants 2021, 10, 1172. [Google Scholar] [CrossRef] [PubMed]
  64. Kato-Noguchi, H.; Yasuda, Y.; Sasaki, R. Soluble sugar availability of aerobically germinated barley, oat and rice coleoptiles in anoxia. J. Plant Physiol. 2010, 167, 1571–1576. [Google Scholar] [CrossRef]
  65. Wang, L.; Zhang, S.; Li, J.; Zhang, Y.; Zhou, D.; Li, C.; He, L.; Li, H.; Wang, F.; Gao, J. Identification of key genes controlling soluble sugar and glucosinolate biosynthesis in Chinese cabbage by integrating metabolome and genome-wide transcriptome analysis. Front. Plant Sci. 2022, 13, 1043489. [Google Scholar] [CrossRef]
  66. Reimann, T.; Zierer, W.; Subert, C.; Sauer, N.; Stadler, R. STP10 encodes a high-affinity monosaccharide transporter and is induced under low-glucose conditions in pollen tubes of Arabidopsis. J. Exp. Bot. 2016, 67, erw048. [Google Scholar] [CrossRef]
  67. Abbaspour, N.; Hurrell, R.; Kelishadi, R. Review on iron and its importance for human health. J. Res. Med. Sci. 2014, 19, 164. [Google Scholar] [PubMed]
  68. Rout, G.R.; Sahoo, S. Role of iron in plant growth and metabolism. Rev. Agric. Sci. 2015, 3, 1–24. [Google Scholar] [CrossRef]
  69. Friedman, E.J.; Wang, H.X.; Jiang, K.; Perovic, I.; Deshpande, A.R.; Pochapsky, T.C.; Temple, B.R.S.; Hicks, S.N.; Harden, T.K.; Jones, A.M. Acireductone dioxygenase 1 (ARD1) is an effector of the heterotrimeric G protein β subunit in Arabidopsis*. J. Biol. Chem. 2011, 286, 30107–30118. [Google Scholar] [CrossRef]
  70. Hasan, N.; Choudhary, S.; Naaz, N.; Sharma, N.; Laskar, R.A. Recent advancements in molecular marker-assisted selection and applications in plant breeding programmes. J. Genet. Eng. Biotechnol. 2021, 19, 128. [Google Scholar] [CrossRef]
Figure 1. Correlation analysis of different nutritional elements in spinach.
Figure 1. Correlation analysis of different nutritional elements in spinach.
Genes 15 00172 g001
Figure 2. Distribution of 3,356,182 single-nucleotide polymorphisms (SNPs) among the six chromosomes of spinach within 1 Mb.
Figure 2. Distribution of 3,356,182 single-nucleotide polymorphisms (SNPs) among the six chromosomes of spinach within 1 Mb.
Genes 15 00172 g002
Figure 3. Linkage disequilibrium (LD) decay in the genome-wide association study (GWAS) population. The red dotted lines indicate the attenuation distance when the r-squared value of 0.35.
Figure 3. Linkage disequilibrium (LD) decay in the genome-wide association study (GWAS) population. The red dotted lines indicate the attenuation distance when the r-squared value of 0.35.
Genes 15 00172 g003
Figure 4. Manhattan plot for genome-wide association study (GWAS) of oxalate content in spinach using GLM and FarmCPU models. The horizontal and vertical axes represent the genomic position of the SNP and the association power for each SNP, respectively, with the trait expressed as −log10(p). The dashed line shows the Bonferroni-corrected genome-wide threshold.
Figure 4. Manhattan plot for genome-wide association study (GWAS) of oxalate content in spinach using GLM and FarmCPU models. The horizontal and vertical axes represent the genomic position of the SNP and the association power for each SNP, respectively, with the trait expressed as −log10(p). The dashed line shows the Bonferroni-corrected genome-wide threshold.
Genes 15 00172 g004
Table 1. The average, range, standard deviation, and coefficient of variation of nine nutrients in spinach.
Table 1. The average, range, standard deviation, and coefficient of variation of nine nutrients in spinach.
TraitAverageMinimumMaximumRangeStandard DeviationCoefficient of Variation (%)
Chlorophyll A (mg/g)0.7490.480 1.1870.7070.17623.541
Chlorophyll B (mg/g)0.3080.170 0.4660.2960.08226.739
Oxalate (g/kg)4.6141.860 8.890 7.031.360 29.509
Nitrate (g/kg)0.9960.260 2.240 1.980.57557.716
Crude Fiber (%)0.7280.590 0.870 0.280.0638.714
Soluble Sugar (%)0.6470.140 1.640 1.50.420 64.849
Mn (mg/kg)6.564.00 10.70 6.71.36920.864
Cu (mg/kg)1.2350.913.30 2.390.32926.626
Fe (mg/kg)135.5427726518832.20123.757
Table 2. Information on SNPs associated with nutrient elements detected using the FarmCPU model.
Table 2. Information on SNPs associated with nutrient elements detected using the FarmCPU model.
TraitSNPChr.PositionAlleles−Log10(P)R-Square(%)
Chlorophyll ASOVchr4-65150089465150089CA10.93 0.14
SOVchr4-613436346134363GT10.91 0.11
SOVchr4-63008113463008113GC10.54 0.12
SOVchr1-96376944196376944CA9.55 0.08
SOVchr6-1458806626145880662TG8.13 0.06
SOVchr4-26424208426424208GA6.44 0.07
Chlorophyll BSOVchr4-65150089465150089CA10.93 0.06
SOVchr4-1059254374105925437GT8.82 0.04
SOVchr4-1456203884145620388GT7.75 0.02
SOVchr4-1799512054179951205AG7.30 0.03
SOVchr1-93842468193842468AC7.16 0.03
SOVchr6-87982689687982689TA6.94 0.05
SOVchr4-613436346134363GT6.80 0.03
SOVchr4-1638326534163832653GA6.70 0.03
OxalateSOVchr1-50807332150807332GT7.59 2.59
SOVchr2-79590343279590343TC7.23 1.32
SOVchr1-1190695921119069592TC7.22 1.17
SOVchr3-43424516343424516CT7.22 2.41
SOVchr1-50996180150996180GT6.97 1.23
SOVchr1-50996185150996185TC6.97 1.23
SOVchr1-50934532150934532CT6.88 1.25
SOVchr1-50565605150565605CT6.87 1.44
SOVchr1-14238191114238191AG6.75 1.49
SOVchr1-50453692150453692CT6.67 1.20
SOVchr1-50453768150453768GA6.67 1.20
SOVchr1-51019015151019015TC6.54 1.22
SOVchr1-51024100151024100AC6.54 1.22
SOVchr1-50731214150731214CT6.53 1.23
SOVchr3-1364661833136466183TC6.53 1.76
NitrateSOVchr1-64960432164960432TG7.41 0.65
SOVchr1-64744070164744070CT7.34 0.59
SOVchr4-1489469154148946915GT7.32 0.63
SOVchr4-1489470174148947017AG7.32 0.63
SOVchr1-64639785164639785AC7.28 0.56
SOVchr1-64726231164726231AG7.28 0.56
SOVchr1-64873483164873483TC7.21 0.60
SOVchr1-64971908164971908TG7.21 0.60
SOVchr1-64636577164636577AT7.21 0.60
SOVchr1-64255238164255238AT6.99 0.62
SOVchr1-1227819011122781901TC6.95 0.40
SOVchr2-66337447266337447CA6.85 1.04
SOVchr1-1227799341122779934CT6.84 0.38
SOVchr1-1227799631122779963TC6.84 0.38
SOVchr1-1227799671122779967CT6.84 0.38
SOVchr1-64826287164826287TC6.80 0.52
SOVchr1-64490849164490849AG6.75 0.54
SOVchr1-64740737164740737TC6.75 0.54
SOVchr1-64781014164781014TC6.75 0.54
SOVchr1-64983706164983706CA6.75 0.54
SOVchr1-65008639165008639TG6.75 0.54
SOVchr1-1227818861122781886AT6.71 0.39
SOVchr1-94072908194072908CT6.69 0.43
SOVchr4-1489443984148944398TG6.69 0.67
SOVchr1-64667101164667101AC6.67 0.61
SOVchr1-64675881164675881GC6.67 0.61
SOVchr1-64686189164686189TA6.67 0.61
SOVchr1-64772753164772753AG6.67 0.61
SOVchr1-64772980164772980TC6.67 0.61
SOVchr1-64782044164782044TG6.67 0.61
SOVchr1-64782074164782074GA6.67 0.61
SOVchr1-64858451164858451TG6.67 0.61
SOVchr1-64859462164859462CT6.67 0.61
SOVchr1-64871632164871632TG6.67 0.61
SOVchr1-64879489164879489TC6.67 0.61
SOVchr1-64894852164894852AG6.67 0.61
SOVchr1-64898345164898345AG6.67 0.61
SOVchr1-64899576164899576TC6.67 0.61
SOVchr1-64902305164902305TC6.67 0.61
SOVchr1-64905439164905439GA6.67 0.61
SOVchr1-64916919164916919GC6.67 0.61
SOVchr1-64924310164924310AG6.67 0.61
SOVchr1-64949052164949052TG6.67 0.61
SOVchr1-64972031164972031TG6.67 0.61
SOVchr1-64993886164993886AC6.67 0.61
SOVchr1-65009190165009190TC6.67 0.61
SOVchr1-65011906165011906GC6.67 0.61
SOVchr1-65439679165439679TC6.67 0.61
SOVchr1-1227809681122780968GA6.65 0.38
SOVchr1-36836579136836579CT6.63 0.47
SOVchr1-64730231164730231TG6.56 0.56
SOVchr1-64820621164820621TC6.56 0.56
SOVchr1-64906356164906356AC6.56 0.56
SOVchr1-932356519323565AG6.54 0.47
SOVchr4-1618578444161857844CT6.53 0.57
SOVchr4-1618578534161857853GA6.53 0.57
SOVchr4-1618578644161857864GT6.53 0.57
SOVchr4-1618578714161857871GA6.53 0.57
SOVchr1-65248580165248580AC6.53 0.46
Soluble SugarSOVchr4-62896056462896056GA7.15 0.16
SOVchr1-1227795501122779550TA6.89 0.09
SOVchr2-78543061278543061CT6.87 0.18
MnSOVchr5-42208975542208975AG7.06 1.27
SOVchr5-40744482540744482GA6.80 1.67
SOVchr4-54452300454452300CT6.58 1.92
CuSOVchr2-26477848226477848AG31.54 1.36
SOVchr1-622521116225211CA26.58 0.32
SOVchr3-33311145333311145TC10.51 0.12
SOVchr3-1369970573136997057TC9.49 0.10
SOVchr4-25830607425830607AT8.73 0.05
SOVchr5-91037611591037611AC7.81 0.09
SOVchr6-66852190666852190TC6.93 0.08
SOVchr4-521353845213538CT6.75 0.04
FeSOVchr1-43602473143602473CA7.14 42.40
SOVchr5-91913409591913409TC6.71 48.32
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ji, N.; Liu, Z.; She, H.; Xu, Z.; Zhang, H.; Fang, Z.; Qian, W. A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach. Genes 2024, 15, 172. https://doi.org/10.3390/genes15020172

AMA Style

Ji N, Liu Z, She H, Xu Z, Zhang H, Fang Z, Qian W. A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach. Genes. 2024; 15(2):172. https://doi.org/10.3390/genes15020172

Chicago/Turabian Style

Ji, Ni, Zhiyuan Liu, Hongbing She, Zhaosheng Xu, Helong Zhang, Zhengwu Fang, and Wei Qian. 2024. "A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach" Genes 15, no. 2: 172. https://doi.org/10.3390/genes15020172

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