Skip to main content

Genome-wide association study of trace elements in maize kernels

Abstract

Maize (Zea mays L.), a staple food and significant economic crop, is enriched with riboflavin, micronutrients and other compounds that are beneficial for human health. As emphasis on the nutritional quality of crops increases maize research has expanded to focus on both yield and quality. This study exploreed the genetic factors influencing micronutrient levels in maize kernels through a comprehensive genome-wide association study (GWAS). We utilized a diverse panel of 244 inbred maize lines and approximately 3 million single nucleotide polymorphisms (SNPs) to investigate the accumulation of essential and trace elements including cadmium (Cd), cobalt (Co), copper (Cu), nickel (Ni), selenium (Se) and zinc (Zn). Our analysis identified 842 quantitative trait loci (QTLs), with 12 QTLs shared across multiple elements and pinpointed 524 potential genes within a 100 kb radius of these QTLs. Notably ZmHMA3 has emerged as a key candidate gene previously reported to influence the Cd accumulation. We highlighted ten pivotal genes associated with trace element transport including those encoding heavy metal ATPases, MYB transcription factors, ABC transporters and other crucial proteins involved in metal handling. Additionally, haplotype analysis revealed that eight inbred linesaccumulated relatively high levels of beneficial elements while harmful elements were minimized. These findings elucidate the genetic mechanisms underlying trace element accumulation in maize kernels and provide a foundation for the breeding of nutritionally enhanced maize varieties.

Peer Review reports

Introduction

Since the onset of the 21st century, trace element deficiency has emerged as a significant global health concern and a phenomenon often described as “hidden hunger” [1,2,3]. Deficiencies in trace elements in humans can lead to poor growth, intellectual impairment, and increased risk of morbidity and mortality, thus, nutritional imbalances and hidden hunger are pressing global health issues [1, 4].

Various trace elements, including zinc (Zn), iron (Fe), selenium (Se), and copper (Cu), are indispensable for human health and play crucial roles in numerous biological processes [5]. An imbalance in trace element intake whether insufficient or excessive can adversely affect human health [5, 6]. For instance Cu is essential for hemoglobin production and maintaining normal blood levels helping to prevent iron deficiency anemia and related diseases. Conversely, excessive Cu can denature hemoglobin, leading to hemolytic anemia and bile excretion dysfunction [6, 7]. Zn is crucial for a variety of transcription factors that regulate metabolism-related genes [8, 9]. Se, which is vital for human thyroid gland metabolism and antioxidant systems, and is also beneficial to the crops [10, 11]. Deficiencies in cobalt (Co) can impair vitamin B12 synthesis and red blood cell regeneration leading to anemia and other serious conditions [12, 13]. Similarly, Nickel (Ni) is essential for the function of various enzymes and the deficiency can result in metabolic disorders and growth retardation [14, 15]. Cadmium (Cd), although important for fat metabolism and insulin regulation, can be toxic in excess causing acute intoxication [16, 17]. In plants, Cu, Co, Ni, and Se participate in photosynthesis and increase photosynthetic efficiency, enhancing plant disease resistance and stress resistance [18,19,20,21]. Co can increase nitrogen uptake and utilization [21]. Zinc is an important component of plant enzymes, which can promote the growth of plant [22]. Se, Cd, Co is a non-essential trace element in plants, among which Se is a beneficial trace element, Cd and Co are toxic heavy metals, and excessive amounts of these elements can inhibit photosynthesis [20, 23, 24]. Therefore, the accumulation of trace elements in plants not only affects the growth of plants, but also has a certain impact on humans who use them as food.

Maize (Zea mays L.) is a major global food crop that plays a crucial role in ensuring food security and is a primary dietary source of calories, protein and trace elements worldwide [25, 26]. The absorption rates of trace elements such as Zn and Se vary among maize varieties [27, 28]. This highlights the importance of developing high-quality maize varieties with elevated beneficial element content and reduced harmful element content for food security.

Trace elements from the soil are absorbed by plant roots and transported to accumulate in kernels through a complex genetic network regulated by multiple genes [29]. In plants, elements reach the kernels via transporters and metal chelators involving protein families such as the natural resistance-associated macrophage protein (NRAMP), heavy-metal ATPase (HMA), yellow stripe-like (YSL), zinc-regulated transporter (ZRT)/iron-regulated transporter (IRT)-like protein (ZIP), ATP-binding cassette (ABC) transporter, metal tolerance protein (MTP), nicotianamine synthase (NAS) and multidrug and toxic compound extrusion (MATE) [30,31,32,33]. Trace elements affect various transcription factors that regulate plant genes, including WRKY, bZIP, MYB, ERF, HSF, ARF13, SPLs and NACs [32, 34, 35].

Selenium uptake in plants can occur in forms such as selenate (SeO42) and selenite (SeO32), as well as in organic forms such as selenocysteine (SeCys) and selenomethionine (SeMet), with different absorption pathways for each [36,37,38]. For example selenite is metabolized into organic selenium forms within roots and then transported to kernels primarily via the phloem facilitated by transporters such as the nitrate transporter NRT1.1B, which increases the Se content in kernels [39,40,41].

Genome-wide association studies (GWAS) have been instrumental in identifying genetic variations associated with quantitative traits in plants, leading to the identification of critical loci and candidate genes related to mineral element regulation in species like rice and maize [42,43,44,45,46,47,48,49]. Such studies have revealed significant associations between genetic variants and trace metal ion accumulation in maize, contributing to the understanding of the regulatory networks involved.

In this study, we determined six trace element contents of 244 maize inbred lines and the GWAS was performed with 3 million single nucleotide polymorphisms (SNPs) covering the whole genome. Significant SNPs associated with trace elements were identified by GWAS, and key candidate genes were screened to provide a new reference for further understanding of the genetic mechanisms. Several excellent germplasm resources with higher beneficial trace elements and lower harmful metal elements have been identified by analyzing the genetic variation in trace elements in maize kernels. This study aimed to reveal the genetic mechanism underlying trace element accumulation in maize kernels and provide conditions for breeding new varieties of maize with high nutritional quality.

Results

Phenotypic analysis of trace elements

Descriptive statistical analysis was performed on the trace element content data from the populations in the two environments (Table 1; Fig. S1), and the 6 trace element contents in kernels from Jiangmen were higher than those in kernels from Hainan. The average contents of Cd, Co, Cu, Ni, Se, and Zn in the kernels from Hainan were 0.0145, 0.0318, 2.1368, 0.612, 0.0456 and 15.9197 mg/kg, respectively (Table 1). The average contents of the 6 trace elements in Jiangmen_1 and Jiangmen_2 were 0.0382/0.0311, 0.1996/0.23, 5.0306/4.3635, 5.4285/6.4784, 0.0736/0.1406, and 26.0946/26.6259, respectively (Table 1). The phenotypic variation coefficient of the 6 trace elements in the population was large, indicating that there was wide phenotypic diversity (Table 1). The 6 trace element contents showed skewed or incomplete normal distributions (Table 1; Fig. 1). In addition, the Pearson correlation showed that there was a significant correlation between Ni and Co, and a significant negative correlation existed between Ni and Se (Fig. 1; Fig. S2), indicating that the accumulation of trace elements was affected by the environment.

Table 1 The contents of Cd, Co, Cu, Ni, Se and Zn in the maize kernels of Hainan (HN), Jiangmen_1 (JM1) and Jiangmen_2 (JM2)
Fig. 1
figure 1

Correlation analysis and frequency distributions of Cd, Co, Cu, Ni, Se and Zn. Asterisks represent a significant difference. * p < 0.05, ** p < 0.01, *** p < 0.001

Genome wide association study of trace elements

The Q + K model was used for GWAS of Cd, Co, Cu, Ni, Se, and Zn in Hainan, Jiangmen_1 and Jiangmen_2, and the QQ plots showed that the model was more reliable (Fig. 2; Fig. S3; Fig. S4). The SNPs significantly associated with trace elements were also identified. There were 5621, 198 94, 90, 22, and 248 SNPs of Cd, Co, Cu, Ni, Se, and Zn identified in Hainan, respectively, and 176/2343, 9/51, 57/304, 9/53, 53/363 and 63/16 SNPs identified in Jiangmen_1 and Jiangmen_2, respectively (Table S1). The three sets of environmental phenotypic data were analyzed to obtain BLUP values, and the GWAS was performed. A total of 5653, 16, 107, 37, 423 and 155 SNPs were associated with Cd, Co, Cu, Ni, Se, and Zn, respectively (Table S1). Among the four GWASs, the numbers of colocalized SNPs of Cd, Cu, Zn, Co, Ni and Se (HN, JM_1, JM_2, and BLUP) were 4312, 2, 24, 10, 4 and 17, respectively (Table S1).

Fig. 2
figure 2

Distribution of the significant SNPs on chromosomes. The vertical line represents the position of the SNP on chromosome, and the color represents significance

By integrating the SNPs into QTLs, a total of 392 QTLs were identified in the Hainan environment, and 248, 42, 16, 27, 14 and 47 QTLs were obtained for Cd, Co, Cu, Ni, Se, and Zn, respectively. A total of 108 and 230 QTLs were detected in JM1 and JM2, respectively. A total of 280 QTLs were obtained by GWAS of BLUP values (Table S2). A total of 842 QTLs were detected in the three environments and BLUP, with 386, 69, 81, 68, 152 and 97 QTLs related to Cd, Co, Cu, Ni, Se, and Zn, respectively (Table S2). Among them, 71 QTLs were detected under two different GWAS conditions, 22 were detected three times, and 12 were detected four times (Fig. 3; Table S3). There are twelve QTLs were colocalized with two elements (Table 2).

Fig. 3
figure 3

Number of significant QTLs of the six traits co-located by BLUP and multiple environments. HN: Hainan, JM1: Jiangmen_1, JM2: Jiangmen_2

Table 2 The list of significant SNPs identified by two traits

Candidate genes analysis

A total of 524 genes with functional annotations were retrieved by searching for candidate genes within 100 kb upstream and downstream of the QTLs, and 265, 25, 43, 52, 65, and 91 candidate genes were associated with Cd, Co, Cu, Ni, Se, and Zn, respectively, among which 17 were colocated by two elements (Table S4). Zm00001d005190 (ZmHMA3) is a candidate gene, that is a member of the heavy metal transporting ATPase family and is a key gene known by pioneers to regulate Cd transport and accumulation. Eight significant SNPs were found in ZmHMA3, three of which were located on the exon of ZmHMA3.

Eleven key candidate genes were identified by gene function annotations and previous studies; among them were 4 candidate genes related to Cd. A QTL consisting of 6 SNPs was identified by Cd on chr2. This QTL was located 95 kb from Zm00001d005062 and encodes a vacuolar-sorting receptor 1 (VSR) protein that is closely linked to plant root growth and the absorption of ion elements. Furthermore, s2_163018473 was distributed in the exon of heavy metal ATPase 4 (Zm00001d005189, HMA4). Zm00001d012521 (Chx17) is located 35.6 kb from rs8_176060449 and belongs to the cation/H + exchanger (CAX) family. Zm00001d018500 is a member of the ABC transporter family (Table 3). The key candidate gene associated with Co is Zm00001d030738 (Mybr42), which encodes a MYB-related transcription factor 42. Zm00001d009492 is also an ABC transporter G family member and is a candidate gene associated with Cu. Zm00001d032644 encodes the vacuolar protein sorting-associated protein 52 A, which is a candidate gene related to Ni (Table 3). Two key candidate genes related to Se, Zm00001d004136 and Zm00001d043240 encode the glutathione S-transferase family protein and sulfite exporter TauE/SafE family protein, respectively. Two candidate genes associated with Zn content, Zm00001d002343 (Zip3), which encodes a zinc/iron-regulated transporter - like protein 3, and Zm00001d035918 (Myb1), which is a MYB transcription factor (Table 3).

Table 3 Information on the key candidate genes

Haplotype analysis

The Hap2 haplotype of Zm00001d005062 had a significantly higher Cd content than other haplotypes; this was an inferior haplotype that included 45 resources (Fig. 4A; Table S5). The Cd content in Hap3 of HMA4 was significantly lower than that in the other haplotypes, and Hap3 contained 20 resources (Fig. 4B; Table S5). In addition, we analyzed the haplotypes of key genes associated with Cd, including ZmHMA3. The results showed that the Cd content was significantly lower in Hap3 than in the other haplotypes. Hap3 was considered an excellent haplotype, and there were 21 resources of this haplotype (Fig. 4C; Table S5). According to the colocalized QTL data, the key gene Zm00001d038225 was colocalized with both Co and Ni. The haplotype analysis of Zm00001d038225 showed that the contents of Co and Ni in Hap3 were significantly lower than those in the other haplotypes (Fig. 4D, E). Therefore, Hap3 was an excellent haplotype, and contained 23 resources in the germplasm resource (Table S5). Nas5 is located approximately 130 kb away from rs7_180095813 and is associated with Zn. Hap2 of Nas5 contained 36 resources, and the Zn content was significantly higher than that of the other haplotypes (Fig. 4F; Table S5). The haplotype with higher Zn content was considered an excellent haplotype. By haplotype analysis, three excellent germplasm resources (CAU280, CAU41 and CAU510) with lower Co, Ni and higher Zn content were obtained; and five excellent germplasm resources (CAU238, CAU318, CAU390, CAU440 and CAU91) with lower Cd and higher Zn content were obtained (Table S6).

Fig. 4
figure 4

Haplotype analysis of candidate genes. A: Zm00001d005062, B: HMA4 (Zm00001d005189), C: ZmHMA3 (Zm00001d005190), D: Zm00001d038225, E: Zm00001d0038225, F: Nas5 (Zm00001d022557). The asterisk represents a significant difference. * p < 0.05, ** p < 0.01, *** p < 0.001

Discussion

Phenotypic variation of the 6 trace elements in very different environments

Although many studies have been conducted on trace elements in plants, the genes that regulate trace elements in maize remain largely undetermined. In this study we detected 6 trace elements in 244 maize inbred lines and found that they had rich genetic diversity in this population (Table 1; Fig. S1). These results indicate that this population is an excellent resource for increasing the genetic diversity of maize germplasm and contains a large number of favorable haplotypes, which are necessary for successful detection of the genetic basis of target traits. Pearson correlation analysis revealed significant differences in the contents of trace elements among the different environments (Fig. 1) which may be due to differences in the content of mineral elements in soil and climate differences in different regions. Previous studies have shown that the amount of fertilizer applied to soil is significantly related to the accumulation of trace elements in plants [29, 50]. The content of trace elements in maize kernels is influenced by both the environment and heredity. Correlation analysis of the different elements revealed that Co, Ni, Cu and Zn were significantly correlated (Fig. S2). Co and Ni are toxic heavy metals, Cu and Zn are essential trace elements which means that the uptake and transport of these two related metal elements in maize kernels may be regulated by the same pathway [30, 51]. Ni is negatively correlated with Zn, Cu and Se, and the uptake of toxic heavy metals by plants may have an antagonistic effect on the uptake of beneficial trace elements, which is in accordance with the results for Se and arsenic (As) [52]. These results showed that there was little correlation between the uptake of various elements in plants, and it was necessary to select a certain element or two antagonistic elements, beneficial and harmful elements, as specific traits to improve and cultivate new varieties.

The QTLs identified by 6 trace elements in maize kernels

A total of 386 QTLs were associated with Cd (Table S2), accounting for 45.8% of the total number of QTLs. More than 5,000 strong and dense SNPs associated with Cd were identified on chr2: 149–171 Mb (Fig. 2; Table S1). The results for ZmHMA3 and HMA4 in this region, are consistent with the findings of previous studies, indicating the reliability of the results of this study [53].

Co is a beneficial element for plant growth [12]. However, the mechanism underlying its tolerance and accumulation in maize is still poorly understood. Similarly, little is known about Ni and Cu in maize [30]. In this study, we identified new QTLs associated with these trace elements that were not found in previous studies. The uptake, transport and accumulation of trace elements in plants are complex processes regulated by multiple genes. This study provides a new scientific reference for studies on the genetic basis of trace elements in maize. Among the QTLs associated with Zn, three QTLs located at 180 Mb on chr7 were consistent with previous studies, and new QTLs were also identified (Table S1; Table S2) [54].

Se is an essential trace element for plants and a beneficial element for humans [10]. Se and heavy metals are antagonistic and can combine with heavy metals such as Cd and Ni to form selenomeric protein complexes that are expelled from the body [55]. In this study, we not only identified new QTLs and key genes associated with Se, but also identified QTLs associated with Se and Cd (Table 2). Therefore, this study provides not only a basis for breeding Se-rich maize, but also new ideas for cultivating maize varieties with low concentrations of harmful heavy metal elements [52, 56].

In addition, there were some colocalized QTLs. Rs2_10959451 identified by BLUP, Jiangmen_1, and Jiangmen_2, was associated with Ni. Therefore, this key locus is related to the Ni content in maize and can be used as a key site for breeding new maize varieties with lower levels of the harmful element Ni. Six QTLs were colocalized by Ni and Co (Table 2), which was consistent with the results of the element correlation analysis. In previous studies, the transport of Ni and Co in plants was shown to be regulated by the same transporter, such as IRON REGULATED1/Ferroportin1 (IREG1/FPN1) and IREG2/FPN2 have been identified in Arabidopsis thaliana, which are believed to transport Co and Ni within vacuoles [57]. In rice, ARG1 is an ABC transporter that transports Ni2+ and Co2+ in chloroplasts [58]. The QTLs colocalized with the two elements identified in this study provide new insights into the study of the specific transport of Co and Ni transporters.

Two QTLs colocalized with Cd and Se (Table 2) Se can reduce the accumulation of Cd in plants effectively and alleviate Cd poisoning. After the addition of Se to hydroponic wheat the genes involved in Cd uptake are downregulated, resulting in a decrease in Cd accumulation [56]. In addition, Cd, Cu, Co, and Ni colocalized with QTLs, and they are metallic elements. Many metal transporters are involved in the transport of various metal elements [30, 33]. IRT1 transports Fe2+/Mn2+/Zn2+/Co2+/Ni2+ to the plasma membrane [57]. AtMGT1 is a plasma membrane protein that can transport Mg2+/Co2+/Ni2+ [59]. Therefore, the colocalization QTLs obtained in this study are reasonable.

Key candidate genes for cultivating/developing new maize varieties with higher beneficial elements and lower harmful elements

We identified 11 key candidate genes in this study, among which were two genes encoding the ATP-binding cassette transporter (ABC) protein: Zm00001d018500 and Zm00001d009492 (Table 3). The ABC transporter superfamily has many functions and transports a wide range of substrates, such as heavy metal ions and lipids in plants [60]. In maize, the ABC transporter AtPDR8 is an efflux pump for Cd conjugates at the plasma membrane involved in heavy metal resistance [61]. In barley, transcriptome analysis revealed that ABC transporters actively participate in the detoxification of Cu and revealed the mechanism underlying Cu tolerance [62]. However, the ABC transporter is not involved in Cu transport in maize. Therefore, the regulation of Cu transport by ABC transporter in maize was first reported in this study. Zm00001d005062 and Zm00001d032644 are vacuolar sorting receptors and vacuolar protein sorting associated proteins, respectively (Table 3). Vacuoles also play an important role in ion transport [63]. OsHMA3 is considered the primary determinant of Cd accumulation in rice grains, and plays a role in the vacuolar membrane [64]. ZmHMA3 is involved in mediating the vacuolar sequestration of Cd and affects its transport from roots to stems [53]. Cation/H(+) antiporters (CAXs) are a particular class of vacuolar transporters involved in Ca2+, Mg2+ and Mn2+ transport [65]. StCAX1/4 transport proteins have a unique tolerance to Cd stress in potatoes. In poplar, CAX family genes are involved in Cd stress responses [34]. Zm00001d012521 encodes a CAX transport protein that is associated with Cd, and this is a discovery in maize. The proteins encoded by Zm00001d005189 (HMA4) and Zm00001d002343 (Zip3) are typical metal ion transport proteins in plants (Table 3) [32, 33]. Zm00001d030738 (Myb42) and Zm00001d030738 (Myb1) are specifically expressed genes. Myb42 is expressed in the primary roots and grains, while Myb1 is only expressed in the roots (Table 3). The accumulated trace elements in the grains are taken up by the roots from the soil and then transported to the grains. The genes involved in ion uptake and transport are also regulated by transcription factors [32, 34].

In addition we identified two key candidate genes associated with Se (Table 3). Zm00001d004136 encodes a glutathione S-transferase protein. The process of selenite assimilation is influenced by glutathione. Ren et al. identified 29 glutathione S-transferases after tea plants were treated with selenite [66]. Zm00001d043240 encodes a sulfite exporter, the TauE/SafE family protein, which is a sulfate transporter. Sulfate transporters are involved in selenite absorption in plants [36]. The Arabidopsis thaliana mutant deficient in sulfate transport sel1-9 has a T-DNA insertion in the Sultr1;2 promoter and is a selenate-resistant mutant [67]. Therefore, the key candidate genes identified in this study provide genetic resources for improving trace element accumulation in maize kernels and can promote the breeding and development of new varieties based on endogenous genetic variation and genome editing.

Haplotype analysis was used to select germplasm resources

Good germplasm resources are the basis of maize variety improvement and genetic breeding. Exploring and utilizing maize germplasm resources through various means is important for variety improvement and genetic breeding [68]. There are abundant variations in the contents of trace elements in germplasm resources, which meets the requirements of haplotype analysis [69, 70]. In this study, we used haplotype analysis to select excellent germplasm resources. Haplotype analysis of ZmHMA3 and HMA4 not only revealed excellent germplasm resources with lower Cd content, but also revealed resources with higher Cd content, (Fig. 4). The resources with excellent haplotypes of these two genes were almost identical (Table S5), likely because the two genes have similar functions and act simultaneously in the same regulatory pathway. The haplotype analysis of Zm00001d005062 revealed that the germplasm resources with higher Cd content (Fig. 4; Table S5). Cd is a poisonous heavy metal, and its enrichment is harmful to both plants and humans. Therefore these results offer new insights into the breeding of new maize varieties with lower Cd content. Zm00001d022557 is a nicotianamine synthase 5 (Nas5) gene that encodes a class II NAS purportedly involved in the synthesis of the metal ion chelator nicotianamine, which has also been shown to be associated with Zn in previous studies [54]. The maize germplasm resources with higher Zn content were identified through haplotype analysis of Nas5. Zm00001d 038225 was associated with both Ni and Co (Fig. 4; Table 2) and haplotype analysis revealed that the same haplotype had significantly lower levels of these two elements than did the other haplotypes (Fig. 4), which provides an important reference for improving multiple traits at the same time. These excellent haplotypes could be useful for molecular breeding of maize in the future.

Materials and methods

Plant materials and field experiments

This study utilized a natural population of 244 inbred maize lines provided by the laboratory of Professor Lai Jinsheng at China Agricultural University [71]. This collection includes crucial genetic resources such as Huang C, Chang7-2, B73, and Mo17. The experimental crops were cultivated at two locations: Ledong, Hainan (HN, 18.75°N, 109.18°E) in 2020, and Jiangmen, Guangdong (JM, 22.35°N, 113.05°E) in 2021. The two different areas planted in Jiangmen called Jiangmen_1 and Jiangmen_2. The plants were arranged in single rows with a spacing of 25 cm × 60 cm. Each inbred line was self-pollinated, and upon reaching maturity, the plants were harvested individually and preserved for further analysis.

Trace element contents detected in maize kernels

The ears were harvested from different environments and the maize kernels were thrashed by hand. The kernels were subsequently dried in an oven at 80 °C for 3 days and ground into powder. Five grams of each sample was weighed, added to 5 ml nitric acid, incubated for 1 h or overnight, digested with 65% nitric acid in an Anton Paar Multiwave PRO microwave for 45 min, and left to set for 1 h with a temperature gradient of 120–180 °C. It was then heated at 100 ℃ for 30 min, diluted with deionized water to 10 ml, and mixed for subsequent use. Inductively coupled plasma mass spectrometry (ICP-MS) was used to analyze the Cd, Co, Cu, Se, Ni, and Zn contents in maize kernels. Three technical replicates were performed for each digestion.

Statistical analysis of phenotypic and genotypic data

IBM SPSS Statistics 19 was used to calculate the average, range variation, standard deviation, skewness, kurtosis and variation coefficient of the trace elements in inbred lines from different environments and to perform variance analysis. Origin was used to construct images for each element, and the GGally package of the R language was used to analyze the Pearson correlation coefficient [72]. The best linear unbiased prediction (BLUP) value of trace elements was calculated using the R package lme4, and environmental variables, genotype variables, and repetitions were all taken as random effects. Y = µ + a * env + b * geno + e, and Y is the phenotypic value, environment variable (env) and genotype variable (geno) are random effects, a is the environmental effect, b is the genotype effect, that is the breeding value [73].

Genome-wide association study

The whole genome sequencing was performed for the natural population used in this study and 3 million SNPs were obtained. The Q + K mixed linear model in the mixed line model (MLM) was used to reduce false positives for GWAS via gemma 0.98.3 [74]. The 100 kb distance was taken as the LD distance by Plink [75]. Admixture structure software was used to analyze the population (v1.3.0) [76]. A total of 106,009 SNPs were obtained by plink, parameter r2 > = 0.2, window size = 50 and step size = 50, so the SNPs significantly associated with trace element-related traits were identified by the significance threshold of association (p ≤ 1 × 10− 5). Manhattan plots and quantile - quantile (QQ) plots were created using the CMplot package in R [77]. Significant SNPs were intervened into QTLs [44, 45, 78].

Prediction of candidate genes

All candidate genes within 100 kb upstream and downstream of each significant locus were downloaded from the maize database MaizeGDB (https://www.maizegdb.org) based on the B73 AGPv4 reference genome. Candidate genes were annotated and analyzed by NCBI (https://www.ncbi.nlm.nih.gov) and MaizeGDB, which were screened according to gene function annotation.

Data availability

All relevant data are within the manuscript and its additional files.

References

  1. Gödecke T, Stein AJ, Qaim M. The global burden of chronic and hidden hunger: trends and determinants. Global food Secur. 2018;17:21–9.

    Article  Google Scholar 

  2. Stevens GA, Beal T, Mbuya MN, Luo H, Neufeld LM, Addo OY, Adu-Afarwuah S, Alayón S, Bhutta Z, Brown KH. Micronutrient deficiencies among preschool-aged children and women of reproductive age worldwide: a pooled analysis of individual-level data from population-representative surveys. Lancet Global Health. 2022;10(11):e1590–9.

    Article  CAS  PubMed  Google Scholar 

  3. Beal T, Massiot E, Arsenault JE, Smith MR, Hijmans RJ. Global trends in dietary micronutrient supplies and estimated prevalence of inadequate intakes. PLoS ONE. 2017;12(4):e175554.

    Article  Google Scholar 

  4. Bailey RL, West KP Jr, Black RE. The epidemiology of global micronutrient deficiencies. ANN NUTR METAB. 2015;66(Suppl 2):22–33.

    Article  CAS  PubMed  Google Scholar 

  5. Prashanth L, Kattapagari KK, Chitturi RT, Baddam VR, Prasad LK. A review on role of essential trace elements in health and disease. J NTR Univ Health Sci. 2015;4(2):75–85.

    Article  Google Scholar 

  6. Fraga CG. Relevance, essentiality and toxicity of trace elements in human health. MOL ASPECTS MED. 2005;26(4–5):235–44.

    Article  CAS  PubMed  Google Scholar 

  7. Karim N. Copper and human health-a review. J Bahria Univ Med Dent Coll. 2018;8(2):117–22.

    Article  Google Scholar 

  8. Prasad AS. Discovery of human zinc deficiency: its impact on human health and disease. ADV NUTR. 2013;4(2):176–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Tong J, Sun M, Wang Y, Zhang Y, Rasheed A, Li M, Xia X, He Z, Hao Y. Dissection of molecular processes and genetic architecture underlying iron and zinc homeostasis for biofortification: from model plants to common wheat. Int J Mol Sci. 2020;21(23):9280.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Duntas LH, Benvenga S. Selenium: an element for life. Endocrine. 2015;48:756–75.

    Article  CAS  PubMed  Google Scholar 

  11. Gorini F, Sabatino L, Pingitore A, Vassalle C. Selenium: an element of life essential for thyroid function. Molecules. 2021;26(23):7084.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Yamada K. Cobalt: its role in health and disease. Interrelations between Essent Metal ions Hum Dis 2013:295–320.

  13. González-Montaña J, Escalera-Valente F, Alonso AJ, Lomillos JM, Robles R, Alonso ME. Relationship between vitamin B12 and cobalt metabolism in domestic ruminant: an update. ANIMALS-BASEL. 2020;10(10):1855.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Genchi G, Carocci A, Lauria G, Sinicropi MS, Catalano A. Nickel: human health and environmental toxicology. Int J Environ Res Public Health. 2020;17(3):679.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kamboj N, Malik RS, Dhanker P, Kumar A. Importance of nickel in crops. J Pharmacognosy Phytochemistry. 2018;7(3):3470–5.

    CAS  Google Scholar 

  16. Attia SM, Varadharajan K, Shanmugakonar M, Das SC, Al-Naemi HA. Cadmium: an emerging role in adipose tissue dysfunction. EXPOS HEALTH. 2022;14(1):171–83.

    Article  CAS  Google Scholar 

  17. Larregle EV, Varas SM, Oliveros LB, Martinez LD, Antón R, Marchevsky E, Giménez MS. Lipid metabolism in liver of rat exposed to cadmium. FOOD CHEM TOXICOL. 2008;46(5):1786–92.

    Article  CAS  PubMed  Google Scholar 

  18. Yruela I. Copper in plants: acquisition, transport and interactions. FUNCT PLANT BIOL. 2009;36(5):409–30.

    Article  CAS  PubMed  Google Scholar 

  19. Shahzad B, Tanveer M, Rehman A, Cheema SA, Fahad S, Rehman S, Sharma A. Nickel; whether toxic or essential for plants and environment-A review. PLANT PHYSIOL BIOCH. 2018;132:641–51.

    Article  CAS  Google Scholar 

  20. El-Ramady H, Abdalla N, Taha HS, Alshaal T, El-Henawy A, Faizy SEA, Shams MS, Youssef SM, Shalaby T, Bayoumi Y. Selenium and nano-selenium in plant nutrition. ENVIRON CHEM LETT. 2016;14:123–47.

    Article  CAS  Google Scholar 

  21. Hu X, Wei X, Ling J, Chen J. Cobalt: an essential micronutrient for plant growth? FRONT PLANT SCI. 2021;12:768523.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hamzah Saleem M, Usman K, Rizwan M, Al Jabri H, Alsafran M. Functions and strategies for enhancing zinc availability in plants for sustainable agriculture. FRONT PLANT SCI. 2022;13:1033092.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Dos Reis LL, Alho LDOG, de Abreu CB, Gebara RC, Da Silva Mansano A, Melão MDGG. Effects of cadmium and cobalt mixtures on growth and photosynthesis of Raphidocelis subcapitata (Chlorophyceae). AQUAT TOXICOL. 2022;244:106077.

    Article  CAS  PubMed  Google Scholar 

  24. Dias MC, Monteiro C, Moutinho-Pereira J, Correia C, Gonçalves B, Santos C. Cadmium toxicity affects photosynthesis and plant growth at different levels. ACTA PHYSIOL PLANT. 2013;35:1281–9.

    Article  CAS  Google Scholar 

  25. Shiferaw B, Prasanna BM, Hellin J, Bänziger M. Crops that feed the world 6. Past successes and future challenges to the role played by maize in global food security. FOOD SECUR. 2011;3:307–27.

    Article  Google Scholar 

  26. Nuss ET, Tanumihardjo SA. Maize: a paramount staple crop in the context of global nutrition. COMPR REV FOOD SCI F. 2010;9(4):417–36.

    Article  CAS  Google Scholar 

  27. Molina M, Escudey M, Chang AC, Chen W, Arancibia-Miranda N. Trace element uptake dynamics for maize (Zea mays L.) grown under field conditions. PLANT SOIL. 2013;370:471–83.

    Article  CAS  Google Scholar 

  28. Suganya A, Saravanan A, Manivannan N. Role of zinc nutrition for increasing zinc availability, uptake, yield, and quality of maize (Zea mays L.) grains: an overview. Commun Soil Sci Plant Anal. 2020;51(15):2001–21.

    Article  CAS  Google Scholar 

  29. Singh PK, Pratap SG, Tandon PK. The mechanisms of trace element uptake and transport up to grains of crop plants. Sustainable Solutions Elemental Defic Excess Crop Plants 2020:119–33.

  30. Lange B, van Der Ent A, Baker AJM, Echevarria G, Mahy G, Malaisse F, Meerts P, Pourret O, Verbruggen N, Faucon M. Copper and cobalt accumulation in plants: a critical assessment of the current state of knowledge. NEW PHYTOL. 2017;213(2):537–51.

    Article  CAS  PubMed  Google Scholar 

  31. Van der Pas L, Ingle RA. Towards an understanding of the molecular basis of nickel hyperaccumulation in plants. Plants. 2019;8(1):11.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Mei S, Lin K, Williams DV, Liu Y, Dai H, Cao F. Cadmium accumulation in cereal crops and tobacco: a review. Agronomy. 2022;12(8):1952.

    Article  CAS  Google Scholar 

  33. Zeng H, Wu H, Yan F, Yi K, Zhu Y. Molecular regulation of zinc deficiency responses in plants. J PLANT PHYSIOL. 2021;261:153419.

    Article  CAS  PubMed  Google Scholar 

  34. Liu Y, He G, He Y, Tang Y, Zhao F, He T. Discovery of cadmium-tolerant biomacromolecule (StCAX1/4 transportproteins) in potato and its potential regulatory relationship with WRKY transcription factors. INT J BIOL MACROMOL. 2023;228:385–99.

    Article  CAS  PubMed  Google Scholar 

  35. Gu L, Hou Y, Sun Y, Chen X, Wang G, Wang H, Zhu B, Du X. The maize WRKY transcription factor ZmWRKY64 confers cadmium tolerance in Arabidopsis and maize (Zea mays L). PLANT CELL REP. 2024;43(2):44.

    Article  CAS  PubMed  Google Scholar 

  36. Gupta M, Gupta S. An overview of selenium uptake, metabolism, and toxicity in plants. FRONT PLANT SCI 2017, 7:2074.

  37. Raina M, Sharma A, Nazir M, Kumari P, Rustagi A, Hami A, Bhau BS, Zargar SM, Kumar D. Exploring the new dimensions of selenium research to understand the underlying mechanism of its uptake, translocation, and accumulation. PHYSIOL Plant. 2021;171(4):882–95.

    Article  CAS  PubMed  Google Scholar 

  38. Schiavon M, Pilon-Smits EA. The fascinating facets of plant selenium accumulation-biochemistry, physiology, evolution and ecology. NEW PHYTOL. 2017;213(4):1582–96.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang L, Hu B, Li W, Che R, Deng K, Li H, Yu F, Ling H, Li Y, Chu C. OsPT2, a phosphate transporter, is involved in the active uptake of selenite in rice. NEW PHYTOL. 2014;201(4):1183–91.

    Article  CAS  PubMed  Google Scholar 

  40. Zhang L, Hu B, Deng K, Gao X, Sun G, Zhang Z, Li P, Wang W, Li H, Zhang Z. NRT1.1B improves selenium concentrations in rice grains by facilitating selenomethinone translocation. PLANT BIOTECHNOL J. 2019;17(6):1058–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhao XQ, Mitani N, Yamaji N, Shen RF, Ma JF. Involvement of silicon influx transporter OsNIP2; 1 in selenite uptake in rice. PLANT PHYSIOL. 2010;153(4):1871–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yang M, Lu K, Zhao F, Xie W, Ramakrishna P, Wang G, Du Q, Liang L, Sun C, Zhao H. Genome-wide association studies reveal the genetic basis of ionomic variation in rice. Plant Cell. 2018;30(11):2720–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Genome-wide association studies of grain yield and quality traits under optimum and low-nitrogen stress in tropical maize (Zea mays L.). THEOR APPL GENET 2022, 135(12):4351–4370.

  44. Chen W, Cui F, Zhu H, Zhang X, Lu S, Lu C, Chang H, Fan L, Lin H, Fang J. Genome-wide association study of kernel colour traits and mining of elite alleles from the major loci in maize. BMC PLANT BIOL. 2024;24(1):1–11.

    CAS  Google Scholar 

  45. Zhu H, Lai R, Chen W, Lu C, Chachar Z, Lu S, Lin H, Fan L, Hu Y, An Y. Genetic dissection of maize (Zea maysL.) Trace element traits using genome-wide association studies. BMC PLANT BIOL. 2023;23(1):1–14.

    Article  Google Scholar 

  46. Alomari DZ, Eggert K, Von Wirén N, Polley A, Plieske J, Ganal MW, Liu F, Pillen K, Röder MS. Whole-genome association mapping and genomic prediction for iron concentration in wheat grains. Int J Mol Sci. 2018;20(1):76.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Alqudah AM, Sallam A, Baenziger PS, Börner A. GWAS: fast-forwarding gene identification and characterization in temperate cereals: lessons from barley–a review. J ADV RES. 2020;22:119–35.

    Article  PubMed  Google Scholar 

  48. Alomari DZ, Eggert K, Von Wiren N, Alqudah AM, Polley A, Plieske J, Ganal MW, Pillen K, Röder MS. Identifying candidate genes for enhancing grain zn concentration in wheat. FRONT PLANT SCI. 2018;9:1313.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Alomari DZ, Alqudah AM, Pillen K, Von Wirén N, Röder MS. Toward identification of a putative candidate gene for nutrient mineral accumulation in wheat grains for human nutrition purposes. J EXP BOT. 2021;72(18):6305–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. He ZL, Yang XE, Stoffella PJ. Trace elements in agroecosystems and impacts on the environment. J TRACE ELEM MED BIO. 2005;19(2–3):125–40.

    Article  CAS  Google Scholar 

  51. Persans MW, Salt DE. Possible molecular mechanisms involved in nickel, zinc and selenium hyperaccumulation in plants. Biotechnol Genet Eng Rev. 2000;17(1):389–416.

    Article  CAS  PubMed  Google Scholar 

  52. Chauhan R, Awasthi S, Indoliya Y, Chauhan AS, Mishra S, Agrawal L, Srivastava S, Dwivedi S, Singh PC, Mallick S. Transcriptome and proteome analyses reveal selenium mediated amelioration of arsenic toxicity in rice (Oryza sativa L). J HAZARD MATER. 2020;390:122122.

    Article  CAS  PubMed  Google Scholar 

  53. Tang B, Luo M, Zhang Y, Guo H, Li J, Song W, Zhang R, Feng Z, Kong M, Li H. Natural variations in the P-type ATPase heavy metal transporter gene ZmHMA3 control cadmium accumulation in maize grains. J EXP BOT. 2021;72(18):6230–46.

    Article  CAS  PubMed  Google Scholar 

  54. Wu D, Tanaka R, Li X, Ramstein GP, Cu S, Hamilton JP, Buell CR, Stangoulis J, Rocheford T, Gore MA. High-resolution genome-wide association study pinpoints metal transporter and chelator genes involved in the genetic control of element levels in maize grain. G3. 2021;11(4):jkab59.

    Article  Google Scholar 

  55. Zhao B, Zhao J, Zhou S, Wu X, Xu X, Yang R, Yuan Z. Selenium and toxic metals in human hair of the Dashan Region, China: concentrations, sources, and antagonism effect. ECOTOX ENVIRON SAFE. 2023;250:114479.

    Article  CAS  Google Scholar 

  56. Di X, Jing R, Qin X, Wei Y, Liang X, Wang L, Xu Y, Sun Y, Huang Q. Transcriptome analysis reveals the molecular mechanism of different forms of selenium in reducing cadmium uptake and accumulation in wheat seedlings. Chemosphere. 2023;340:139888.

    Article  CAS  PubMed  Google Scholar 

  57. Morrissey J, Baxter IR, Lee J, Li L, Lahner B, Grotz N, Kaplan J, Salt DE, Guerinot ML. The ferroportin metal efflux proteins function in iron and cobalt homeostasis in Arabidopsis. Plant Cell. 2009;21(10):3326–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Li H, Liu Y, Qin H, Lin X, Tang D, Wu Z, Luo W, Shen Y, Dong F, Wang Y. A rice chloroplast-localized ABC transporter ARG1 modulates cobalt and nickel homeostasis and contributes to photosynthetic capacity. NEW PHYTOL. 2020;228(1):163–78.

    Article  CAS  PubMed  Google Scholar 

  59. Li L, Tutone AF, Drummond RS, Gardner RC, Luan S. A novel family of magnesium transport genes in Arabidopsis. Plant Cell. 2001;13(12):2761–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Jalmi SK. The role of ABC transporters in metal transport in plants. Plant Metal and Metalloid transporters. Springer; 2022. pp. 55–71.

  61. Kim D, Bovet L, Maeshima M, Martinoia E, Lee Y. The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2007;50(2):207–18.

    Article  CAS  PubMed  Google Scholar 

  62. Wa Lwalaba JL, Zvobgo G, Gai Y, Issaka JH, Mwamba TM, Louis LT, Fu L, Nazir MM, Kirika BA, Tshibangu AK. Transcriptome analysis reveals the tolerant mechanisms to cobalt and copper in barley. ECOTOX ENVIRON SAFE. 2021;209:111761.

    Article  Google Scholar 

  63. Rekha K, Usha B, Keeran NS. Role of ABC transporters and other vacuolar transporters during heavy metal stress in plants. Metal and nutrient transporters in Abiotic Stress. Elsevier; 2021. pp. 55–76.

  64. Miyadate H, Adachi S, Hiraizumi A, Tezuka K, Nakazawa N, Kawamoto T, Katou K, Kodama I, Sakurai K, Takahashi H. OsHMA3, a P1B-type of ATPase affects root‐to‐shoot cadmium translocation in rice by mediating efflux into vacuoles. NEW PHYTOL. 2011;189(1):190–9.

    Article  CAS  PubMed  Google Scholar 

  65. He F, Shi Y, Li J, Lin T, Zhao K, Chen L, Mi J, Zhang F, Zhong Y, Lu M. Genome-wide analysis and expression profiling of Cation/H + exchanger (CAX) family genes reveal likely functions in cadmium stress responses in poplar. INT J BIOL MACROMOL. 2022;204:76–88.

    Article  CAS  PubMed  Google Scholar 

  66. Ren H, Li X, Guo L, Wang L, Hao X, Zeng J. Integrative transcriptome and proteome analysis reveals the absorption and metabolism of selenium in tea plants [Camellia sinensis (L.) O. Kuntze]. FRONT PLANT SCI. 2022;13:848349.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Shibagaki N, Rose A, McDermott JP, Fujiwara T, Hayashi H, Yoneyama T, Davies JP. Selenate-resistant mutants of Arabidopsis thaliana identify Sultr1; 2, a sulfate transporter required for efficient transport of sulfate into roots. Plant J. 2002;29(4):475–86.

    Article  CAS  PubMed  Google Scholar 

  68. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.

    Article  CAS  PubMed  Google Scholar 

  69. Yu S, Wu J, Wang M, Shi W, Xia G, Jia J, Kang Z, Han D. Haplotype variations in QTL for salt tolerance in Chinese wheat accessions identified by marker-based and pedigree-based kinship analyses. Crop J. 2020;8(6):1011–24.

    Article  Google Scholar 

  70. He Q, Tang S, Zhi H, Chen J, Zhang J, Liang H, Alam O, Li H, Zhang H, Xing L. A graph-based genome and pan-genome variation of the model plant Setaria. NAT GENET 2023:1–11.

  71. Liu H, Shi J, Sun C, Gong H, Fan X, Qiu F, Huang X, Feng Q, Zheng X, Yuan N. Gene duplication confers enhanced expression of 27-kDa γ-zein for endosperm modification in quality protein maize. Proc Natl Acad Sci. 2016;113(18):4964–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Schloerke B, Cook D, Larmarange J, Briatte F, Marbach M, Thoen E, Elberg A, Toomet O, Crowley J, Hofmann H. GGally: Extension to’ggplot2’. 2021.

  73. Douglas Bates MM, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J STAT SOFTW. 2015;67(1):1–48.

    Google Scholar 

  74. Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. NAT GENET. 2012;44(7):821–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, De Bakker PI, Daly MJ. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. GENOME RES. 2009;19(9):1655–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Li LC. Circle Manhattan plot. R Package Version 2018, 3(2).

  78. A SNP-based GWAS and functional haplotype-based GWAS of flag leaf-related traits and their influence on the yield of bread wheat (Triticum aestivum L.). THEOR APPL GENET 2021, 134(12):3895–3909.

Download references

Acknowledgements

We thank Dr. Jinsheng Lai from China Agricultural University for generously providing the inbred maize lines and their genotypes for our research.

Funding

This work was supported by the GDAS’ Project of Science and Technology Development (2022GDASZH-2022010102), the National Natural Science Foundation of China (32072027), the Guangdong Province Special Projects in Key Fields of Ordinary Colleges and Universities, the Guangdong Province Key Construction Discipline Research Ability Enhancement Project (2022ZDJS023), the Special Project for Rural Revitalization Strategy in Guangdong Province (2022-NPY-00-023-5), the Basic and Applied Basic Research Fund of Guangdong Province (2022A1515110760), the stable support project of Guangdong Academy of Sciences, and germplasm innovation and new variety breeding of heat resistant fresh-eating maize.

Author information

Authors and Affiliations

Authors

Contributions

Y.Q and XL designed the research; W.C, X.Z and X.L analyzed data; W.C wrote the paper; C.L and H.C carried out the experiments. Z.C and Q.W reviewed the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Hailong Chang or Qinnan Wang.

Ethics declarations

Ethics approval and consent to participate

All experimental studies on plants complied with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it.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-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, W., Li, X., Zhang, X. et al. Genome-wide association study of trace elements in maize kernels. BMC Plant Biol 24, 724 (2024). https://doi.org/10.1186/s12870-024-05419-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05419-4

Keywords