Next Article in Journal
The Colonization and Effect of Isaria cateinannulata on Buckwheat Sprouts
Previous Article in Journal
The Impact of N.I. Vavilov on the Conservation and Use of Plant Genetic Resources in Scandinavia: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Genomic Regions Associated with High Grain Zn Content in Polished Rice Using Genotyping-by-Sequencing (GBS)

by
Goparaju Anurag Uttam
,
Karre Suman
,
Veerendra Jaldhani
,
Pulagam Madhu Babu
,
Durbha Sanjeeva Rao
,
Raman Meenakshi Sundaram
and
Chirravuri Naga Neeraja
*
Department of Biotechnology, Crop Improvement Section, ICAR-Indian Institute of Rice Research (IIRR), Hyderabad 500030, India
*
Author to whom correspondence should be addressed.
Plants 2023, 12(1), 144; https://doi.org/10.3390/plants12010144
Submission received: 8 November 2022 / Revised: 24 November 2022 / Accepted: 29 November 2022 / Published: 28 December 2022
(This article belongs to the Topic Plant Nutrition Biofortification)

Abstract

:
Globally, micronutrient (iron and zinc) enriched rice has been a sustainable and cost-effective solution to overcome malnutrition or hidden hunger. Understanding the genetic basis and identifying the genomic regions for grain zinc (Zn) across diverse genetic backgrounds is an important step to develop biofortified rice varieties. In this case, an RIL population (306 RILs) obtained from a cross between the high-yielding rice variety MTU1010 and the high-zinc rice variety Ranbir Basmati was utilized to pinpoint the genomic region(s) and QTL(s) responsible for grain zinc (Zn) content. A total of 2746 SNP markers spanning a genetic distance of 2445 cM were employed for quantitative trait loci (QTL) analysis, which resulted in the identification of 47 QTLs for mineral (Zn and Fe) and agronomic traits with 3.5–36.0% phenotypic variance explained (PVE) over the seasons. On Chr02, consistent QTLs for grain Zn polished (qZnPR.2.1) and Zn brown (qZnBR.2.2) were identified. On Chr09, two additional reliable QTLs for grain Zn brown (qZnBR.9.1 and qZnBR.9.2) were identified. The major-effect QTLs identified in this study were associated with few key genes related to Zn and Fe transporter activity. The genomic regions, candidate genes, and molecular markers associated with these major QTLs will be useful for genomic-assisted breeding for developing Zn-biofortified varieties.

1. Introduction

Micronutrient malnutrition or hidden hunger has become a major issue for people in developing countries, especially in Asia and Africa [1,2]. Marginal intake of micronutrients has been shown to contribute to increasing mortality rates, affect livelihoods, and result in adverse effects on millions of school-going children and pregnant women [3]. About three billion people worldwide are estimated to be affected and more than 24,000 people die daily owing to micronutrient malnutrition or hidden hunger [4]. Recently, the United Nations (UN) declared that tackling micronutrient deficiencies is one of the sustainable development goals (SDGs) and set a goal of SDG 2 to be achieved by 2035 (https://sustainabledevelopment.un.org/sdgs (20 September 2022)). Among different micronutrients, zinc (Zn) and iron (Fe) are the most essential and are associated with the development of physical and mental health in humans and reducing the risk of micronutrient malnutrition [5].
Microelements, such as iron (Fe), zinc (Zn), copper (Cu), and manganese (Mn), are found in many enzymes and have great importance for maintaining normal metabolic pathways in humans. Zn serves as a major co-factor for more than 300 enzymes and 2000 transcription factors and is involved in the metabolism of carbohydrates, lipids, proteins, and nucleic acids in humans [6,7]. It is the only metal present in all six enzyme classes (oxidoreductase, transferase, hydrolases, lyases, isomerases, and ligases) [8]. One third of the human population, particularly children and women, suffer from Zn deficiency-associated health complications such as stunting, diarrhea, reduced immunity, poor cognitive development, and skin problems [9]. According to World Health Organization (WHO), Zn deficiency is sixth among the top 10 major causes of illness in developing nations, affecting 27–30% of the global population. A large portion of the global population suffers from mineral malnutrition, as they depend on plant-based diets that have a low mineral/calorific density [10].
Rice is a predominant and major staple food crop that provides the caloric needs of half of the global population [11,12]. A total of 483.3 million tons of milled rice is produced globally, according to a recent survey (ICRISAT., 2018). To meet the need, however, production must rise to 800 to 900 million tons by 2025 [13]. Mega varieties under cultivation are a poor source of micronutrients, particularly Zn in polished form (12 ± 14 ppm), and provide only one fifth of daily recommended Zn requirements (15 mg/day) [14,15]. In other words, 11 to 15 mg of Zn must be present in 220 gm of rice, which is the per capita amount of rice identified for India. The present international target as well as that from the AICRIP biofortification trial is 28 ppm, which can only meet half of the RDA. Through the AICRIP rice biofortification trial, medium slender, long slender, and short bold grain type rice varieties were released to some of the states in the country. Therefore, there is a need to develop biofortified varieties suitable for rice consumers under aromatic grain and other grain types. To ameliorate this situation, researchers are working to enhance the grain Zn content in rice by different approaches. Biofortification has emerged as a cost-effective and sustainable strategy to tackle mineral deficiencies by enhancing micronutrient content without compromising the yield achieved through agronomic, conventional, and biotechnological approaches [16,17]. The availability of genetic variability in grain mineral nutrients in the rice germplasm can be employed to develop biofortified high-yielding rice varieties [18,19].
Understanding the genetic potential of the genotypes as well as interaction between genotype and environment is essential for developing successful biofortified rice varieties. Developing Zn-biofortified rice is challenging, as it involves complex genetics of the trait, genetic interactions such as epistasis, and environmental factors such as soil and water [14,20,21,22]. Several genetic studies have also been carried out to identify markers linked with quantitative trait loci (QTLs) for high Zn in grains, which would help in the development of biofortified rice varieties through marker-assisted selection (MAS) [23]. QTL mapping provides information for identification of genomic regions associated with targeted traits of interest (Zn) by combining genetic information with phenotypic data [24]. Several QTLs with moderate to high phenotypic variance were reported for grain Zn on all 12 chromosomes of rice [25]. Recent studies on QTL mapping for mineral elements in rice using bi-parental and multi-parental mapping populations have identified multiple loci and demonstrated the genetic complexity of grain micronutrient traits [26].
Compared with other molecular markers (RFLP, RAPD, SSRs, etc.) single nucleotide polymorphisms (SNPs) are preferred because of their uniform distribution and wide occurrence across the genome (one marker every 100–500 bp, specific to species), thus making them the ideal choice for constructing high-density linkage maps and identification of markers closely associated with the trait of interest [27,28]. Next generation sequencing (NGS) technology with reduced cost provides abundant sequence information along with significant improvements in genome coverage and time [29]. Genotyping by sequencing (GBS), an advancement in NGS technology based on genome target reduction and restriction enzyme use has gained popularity as a cost-effective method in the development of genome-wide markers for genetic studies [30].
Numerous methods for reduced representation genome sequencing have been developed and include genotyping by sequencing (GBS) [31], double digest restriction-site-associated DNA sequencing (ddRAD) [32], and specific-locus amplified fragment sequencing (SLAF) [33]. GBS is being deployed for genetic studies in various crops [34,35], particularly GWAS analysis, diversity studies, genomic selection (GS), marker and gene discovery, and high-resolution QTL mapping in rice [36,37,38].
The present study aims to construct a high-density linkage map and identify QTLs for high grain Zn in polished rice and related agronomic traits using a recombinant inbred line (RIL) population derived from a MTU1010 X Ranbir Basmati cross. Putative candidate genes present in the QTLs identified for high grain Zn in polished rice were retrieved. The QTLs and putative candidate genes identified in the study can be deployed to advance genetic research and breeding applications for enhancing grain Zn content in rice.

2. Results

2.1. Phenotypic Variation in MTU1010 X Ranbir Basmati RIL Population

Wide variations were observed among the RILs across the four seasons for the nine agronomic, yield-, and mineral-related traits, zinc (Zn) in brown rice (ZnBR; 8.5–37.8 ppm), Zn in polished rice (ZnPR; 5.3–31.7 ppm), iron (Fe) in brown rice (FeBR; 1.2–13.3 ppm), Fe in polished rice (FePR; 0.4–6.4 ppm), days to 50% flowering (DFF; 37–131 days), plant height (PH; 8–157 cm), number of tillers per plant (NT; 7–27), panicle length (PL; 16–33 cm), and single plant yield (SPY; 3.6–51 gm) (Table 1). Normal distribution was observed for most of the traits as per the descriptive statistics and histograms across seasons (Supplementary Figures S1 and S2).
In WS 17, G3 RIL showed Zn content above the threshold value of 24 ppm as well as SPY above 26 gm of yield check (IR64) (Figure 1a). In WS 20, G5 RIL showed Zn content above the threshold value (24 ppm) and G3 along with another seven RILs noted SPY. In the DS 18, G1 noted Zn content above the threshold value (24 ppm) and nine lines noted SPY above the yield check (26 gm). In DS 20, G6 noted Zn content above the threshold value and G3 line noted on-par yield with yield check. Based on mean Zn content and SPY of the four seasons, none of the genotypes showed Zn content above the threshold value, whereas seven RILs showed SPY ≥ the yield check (26 gm) (Figure 1b).

2.2. Statistical Analysis

2.2.1. Correlation Analysis

The relation between the mean values of individual traits varied from wet season to dry season in the correlation analysis. In the wet season of 2017–2020, positive correlations were observed between SPY and DFF at 0.05 or 0.01 probability. SPY showed significant negative correlation at 0.05 or 0.01 probability. Panicle length showed significant negative correlation at 0.05 probability with FeBR. NT showed significant negative correlation at 0.05 probability with ZnPR. In the dry season of 2018–2020, DFF showed significant negative correlation with FeBR and FePR at 0.001 and 0.05 probability, respectively. Significant positive correlations were observed among the four mineral traits, namely ZnBR, ZnPR, FeBR, and FePR, except for FeBR and FePR in the dry season of 2018–2020 (see Table 2). Significant positive correlations were observed among the four mineral traits ZnBR, ZnPR, FeBR, and FePR, except for FeBR and FePR in the dry season of 2018–2020.
The combined ANOVA across four seasons or environments indicating significant variance at probability < 0.001 was observed among the genotypes (Supplementary Table S5).

2.2.2. PCA

To explain the relationships among agronomic, yield-related, and mineral traits, PCA (Principal Component Analysis) with 92 RILs was performed using the mean values of all the four seasons.
The nine traits were resolved into nine principal components (PCs). The first four PCs noted an eigenvalue ≥ 1 (Figure 2) with a cumulative proportion of 75.56% (see Table 2). ZnBR, ZnPR, and FeBR were significant in PC1. PH and PL were negatively significant in PC2; DFF, NT, and SPY were significant in PC3. DFF and NT were significant in PC4 (Figure 2 and Table 3).

2.3. Ranking of Genotypes

Ranking of genotypes was conducted using SPY and ZnPR of multi-seasonal data in the “Metan” package of ‘R’. The genotypes G3, G52, G68, G70, G73, G75, and G82 were grouped in WS17 (Figure 3).

2.4. High-Density Genetic Map

The genome-wide polymorphic SNPs were evaluated for expected segregation ratio using chi-square analysis in MTU1010 X Ranbir Basmati RIL population. Out of 32,245 SNPs between the parents, a total of 23,183 SNPs with contrasting alleles between the parents were found to be homozygous polymorphic and displayed segregation within the mapping population. A genetic map was constructed using 2746 SNPs spanning a genetic distance of 2445 cM (Table 4). The total number of polymorphic SNPs mapped across 12 chromosomes ranged from 28 (Chr04) to 418 (Chr06), and the average genetic distance ranged from 122 cM (Chr04) to 371 cM (Chr02) (Table 4). A large variation in marker density was seen across the chromosomes, with the highest marker density on Chr12 (2.05 SNPs/cM) and lowest on Chr04 (0.23 SNPs/cM) (Table 4).

2.5. QTL Mapping for Agronomic and Mineral Traits

A total of 47 QTLs were identified for eight out of nine traits (except PH). These included 24 QTLs for agronomic and yield traits (11 QTLs for DFF, 4 QTLs for NT, 5 QTLs for PL, and 4 QTLs for SPY) and 23 QTLs for mineral traits (5 QTLs each for ZnBR and ZnPR, 4 QTLs for FeBR, and 9 QTLs for FePR) across four seasons (Table 5 and Table 6). The phenotypic variation explained (PVE) ranged from 3.52% to 36.0%. Altogether, 39 major effect QTLs (PVE ≥ 10%) and 8 minor effect QTLs (PVE < 10%) were identified (Table 5 and Table 6). The highest number of QTLs (16) was observed in the wet season (WS) of 2017.

2.5.1. QTLs for Mineral Traits

Grain Zn Brown (ZnBR)

Five major effect QTLs were identified on Chr02 (qZnBR.2.1; PVE 20.80%), Chr05 (qZnBR.5.1; PVE 23.49%), andChr04 (qZnBR.4.1; PVE 16.01%), and one consistent QTL was identified on Chr09 (qZnBR.9.1 PVE 11.40% and qZnBR.9.2 13.88%) in WS20 and DS20 seasons (Table 5). The majority of QTLs showed negative additive value, indicating that the allele contribution was from the donor parent Ranbir Basmati.

Grain Zn Polished (ZnPR)

Five QTLs were identified with PVE ranging from 8.10 to 19.42%. Two major-effect consistent QTLs qZnPR.2.1 (PVE 19.42%) and qZnPR.2.2 (11.78%) were identified on Chr02 in WS17 and WS20. Further, two QTLs, one major (qZnPR.6.1; PVE 12.34%) and one minor (qZnPR.6.2; PVE 8.33%), were identified on Chr06. One minor-effect QTL (qZnPR.7.1; PVE 8.10%) was also identified on Chr07 (Table 5).

Grain Fe Brown (FeBR)

Four major QTLs for FeBR trait were identified on Chr02 (qFeBR.2.1; PVE 15.73%), Chr03 (qFeBR.3.1; PVE 16.78%), and Chr06 (qFeBR.6.1; PVE 11.27%) for WS17 and Chr05 (qFeBR.5.1; PVE 34.15%) for WS20. All the QTLs showed negative additive values, indicating that the alleles were contributed from the donor parent (Ranbir Basmati) (Table 5).

Grain Fe Polished (FePR)

Nine QTLs were identified for the FePR trait, with the PVE ranging from 3.52% to 36.0%. Three major QTLs on Chr01 (qFePR.1.1; PVE 28.73%), Chr05 (qFePR 5.1; PVE 11.97%), and Chr06 (qFePR.6.1; PVE 12.43%) were identified during WS17. During DS18, four QTLs, two major and two minor, were identified each on Chr02 (qFePR 2.1; PVE 3.52%), Chr08 (qFePR.8.1; PVE 18.27%), Chr10 (qFePR.10.1; PVE 5.53%), and Chr12 (qFePR.12.1; PVE 21.55%). One minor QTL was identified on Chr05 (qFePR.5.2; PVE 8.68%) for WS20 and one major QTL was identified on Chr03 (qFePR.3.1; PVE 36.0%) for DS20 (Table 5).

2.5.2. QTLs for Agronomic and Yield-Related Traits

Days to 50% Flowering (DFF)

A total of nine QTLs were identified for the DFF trait across four seasons. The phenotypic variance (PVE) ranged from 6.63% to 21.35%. One QTL each was observed on Chr01 (qDFF 1.1, PVE 8.86%), Chr03 (qDFF.3.1, PVE 9.43%), Chr08 (qDFF.8.1, PVE 20.96%), and Chr11 (qDFF.11.1, PVE 19.25%).
Three major effect QTLs, qDFF.2.1 with PVE 16.05%, were observed in WS17, and another two consistent QTLs (qDFF.2.2, PVE 10.32% and 12.25%; and qDFF.2.3, PVE 14.40% and 14.31%) were observed in DS18 and DS20, respectively. On Chr06, one major (qDFF.6.2; PVE 21.35%) and one minor effect (qDFF.6.1; PVE 6.63%) QTLs were identified (Table 6).

Number of Tillers (NT)

Four major QTLs, two QTLs on Chr09 qNT.9.1 (PVE 17.94%) and qNT.9.2 (PVE 10.37%) were identified. The additive effect of these two QTLs indicates inheritance from a donor parent (Ranbir Basmati). Furthermore, two QTLs each on Chr06, qNT.6.1 (PVE 12.03%) and Chr10, qNT.10.1 (PVE 11.02%) were identified (Table 6).

Panicle Length (PL)

Four major QTLs for panicle length (PL) were identified on Chr02 (qPL.2.1) with PVE 13.29% in WS17. Two major QTLs on Chr11, qPL.11.1 (PVE 12.39%) and qPL.11.2 (PVE 16.51%) were identified for DS18 and DS20. Further, one major effect QTL on Chr12, qPL.12.1 with PVE 15.88% was identified for WS20 and one minor effect QTL qPL.9.1 with PVE 8.03% on Chr09 was observed (Table 6).

Single Plant Yield (SPY)

Four major QTLs were identified for three seasons, with the PVE ranging from 11.26% to 15.32%. Two QTLs, namely qSPY.9.1 (PVE 11.28%) and qSPY.9.2 (PVE 12.99%), were identified on Chr09. Another two QTLs each on Chr03 (qSPY.3.1; PVE 11.26%) and Chr06 (qSPY.6.1; PVE 15.32%) were identified with an additive effect value of −4.98, indicating the allele inheritance from donor parent Ranbir Basmati (Table 6).

2.6. Co-Localization of QTLs

In our present study, the co-localization of QTLs for traits was observed in the genic regions of Chr02, 05, and 06. DFF (days to 50% flowering) on Chr02 (SNP_ 24248265–SNP_ 10686645) was co-localized with four mineral trait QTLs: FeBR, FePR, ZnBR, and ZnPR. Similarly, the region on Chr05 (SNP_14288078–SNP_24423261) FePR QTL was co-localized with ZnBR QTL (Table 5, Table 6 and Figure 4).

2.6.1. Consistent QTLs across Seasons

Two major consistent QTLs, qDFF.2.2 and qDFF.2.3, were identified across two seasons (DS18 and DS20), and four major mineral QTLs, qZnBR.2.1, qZnPR.2.2 on Chr02, were consistent in WS17 and WS20. Another two QTLs, qZnBR.9.1 and qZnBR.9.2 on Chr09, were consistent in DS20 and WS20 (Table 5, Table 6 and Figure 4).

2.6.2. Common QTLs

Two common QTLs for grain Zn were identified on Chr02. The QTLs for grain Zn polished, qZnPR.2.1 (PVE 19.42%), and grain Zn brown, qZnBR.2.1 (PVE 20.80%), on Chr02 were located in the 25.68–24.25 Mb region (Table 5, Table 6 and Figure 4).

2.7. Analysis of Epistatic Interactions

A total of 170 epistatic interactions were detected for agronomic and yield (83) along with grain mineral (87) traits across four seasons (Supplementary Table S1). Eight epistatic interactions for five traits (PH, NT, FeBR, FePR, and ZnPR) were considered as the main effect QTLs, with PVE ranging from 11.20 to 26.83%. For two QTLs of ZnPR, one on Chr05 interacted with Chr10 (PVE 12.15%) and another QTL on Chr07 interacted with Chr11 (PVE 17.24%). Epistatic interaction for FeBR trait was observed on Chr04, which interacted with the locus on Chr05 (PVE 16.16%). Similarly, another epistatic interaction was observed for FePR on Chr02, which interacted with locus on Chr10 (PVE 26.83%). Two interactions were observed for NT on Chr01 with two loci of Chr02 and Chr08, and another QTL for NT on Chr02 interacted with loci on Chr05. Two digenic epistatic QTL interactions within the two loci of Chr02 were observed for PH, with PVE 8.74% and 24.49%, respectively (Figure 5).

2.8. Candidate Genes Underlying QTLs for Zn QTLs

A total of 1637 genes present within the identified QTL region for Zn brown (ZnBR) and polished (ZnPR) on Chr02, Chr06, and Chr09 were retrieved (Supplementary Table S2). They were annotated and categorized broadly on cellular, molecular, and biological processes using WEGO (Web Gene Ontology Plot) software (Supplementary Figure S3). Among several functional categories, 36 putative candidate genes (seven on Chr02, 18 on Chr06, and 11 on Chr09) belonging to transporter activity roles were shortlisted for the targeted QTLs (Supplementary Table S3). Out of 36 genes, two candidate genes on Chr06 showed a putative function of Zn and Fe transport; thus, the role of these candidate genes in Zn transport and metabolism was evaluated using knetminer software.
Os06g0705700 present in the QTL region of qZnPR.6.2 encoded by peptide transporter was found to be linked with 17 genes and three QTLs. Knetminer network analysis showed that this gene has a role in Zn and Fe transport (Figure 6a), and another gene Os06g0706100 present in the same QTL region on Chr06 is encoded by NRAMP4 and has a role in metal ion transport, particularly Zn transport. The two putative candidate genes (Os06g0705700 and Os06g0706100) within the identified QTL region are being further investigated (Figure 6b).

2.9. Expression Analysis of Identified Candidate Genes

To predict the expression pattern of selected Zn- and Fe-related candidate genes (Os06g0705700 and Os06g0706100), expression profiling was carried out using the two parental lines MTU1010 and Ranbir Basmati, their derived inbred lines (I-42, I-44, and I-83), two mega rice varieties (IR64, BPT5204), a set of Zn-biofortified varieties (ZINCORICE, PAUSTIC1, PAUSTIC9) and a high-Zn landrace karuppunel. The recipient parent MTU1010, which has low Zn content, showed low expression of both the genes, whereas Ranbir Basmati showed higher expression levels (Supplementary Figure S4).

3. Discussion

One third of the global population is reported to be suffering from a lack of sufficient Zn nutrition [39]. Biofortification is considered as long-term sustainable strategy to address micronutrient malnutrition by enhancing the grain micronutrient density. Over the last decade, using conventional breeding approaches, a few biofortified rice varieties with high grain Zn have been developed and released globally, as Zn plays a critical role in several cellular and metabolic activities [40,41]. Though wide variability for grain Fe is available in brown rice, because of polishing, around 70 to 80% of grain Fe is lost [42].
Marker-assisted breeding has a vast potential to achieve desirable phenotypic variations in less time through the deployment of molecular markers linked to QTLs for desirable traits [43]. Large QTL regions identified with low-density SSR markers may have undesirable linkages, resulting in unsuccessful introgression [44]. In some of the recent studies, undesirable linkage groups were successfully eliminated by employing NGS-generated markers such as SNPs by identifying the recombinants. Next generation sequencing (NGS) technologies have become potential tools for the discovery of millions of SNPs in a cost-effective manner. The genotyping-by-sequencing (GBS) technique has facilitated the identification of key genomic regions for both complex and simple traits and has accelerated the breeding process [45,46,47,48]. ddRAD sequencing facilitates the identification of SNPs by reducing genome complexity in rice and other crops [49].
In the present study, we developed a high-density genetic linkage map from an RIL population MTU1010 X Ranbir Basmati using a GBS approach to identify QTLs for grain Zn and Fe using multi-seasonal phenotype data of the mapping population. Detailed analysis of phenotype data revealed a wide genetic variability within the RIL population for grain Zn and Fe along with agronomic and yield-related traits and showed normal distribution across four seasons. Previous studies indicated that complex inherited traits and micronutrient values exhibit wide variation, and these can be attributed to various genetic and environmental factors [50,51]. We examined data from the two wet and dry seasons in the current study and confirmed our results based on the previous variability studies in rice grain mineral content across the seasons (Kharif and Rabi) [52,53].

3.1. Correlation

Highly significant positive correlations were obtained among ZnBR, ZnPR, FeBR, and FePR, irrespective of season and location; however, there was a negative correlation with single plant yield (SPY). Several studies have reported positive correlation between Zn and Fe in rice [54,55,56,57]. Swamy et al. (2018b) [58] attributed the positive correlation between Zn and Fe contents to commonalities in their pathways and genetic networks of Zn and Fe uptake, translocation, and loading in rice. Negative correlations between yield and Zn in rice were earlier reported by [20,59] (G3: Mean Zn and Mean SPY). A QTL (qGZn9a) identified in Australian wild rice strain O.meridionalis was reported to be associated with an increase in grain Zn levels also found to be concomitant with fertility reduction [60]. The negative linkages between the yield and Zn must be eliminated for developing a successful biofortified variety by generating a greater number of recombinants per cross for identifying positive combinations for grain Zn and yield.

3.2. Identification of Promising RILs

Six stable RILs were identified for four mineral traits in the present study with >24 ppm (Zn) based on stability and G × E interaction analysis across four environments. From the study, we could also identify promising RILs with high grain Zn in polished rice (35.6 ppm) and promising RILs with single plant yield (28.4 g to 38.4 g), which could be further deployed in the breeding of biofortified rice varieties. For rice grain Zn and Fe, stability and G × E analyses are used for identification of stable donors [47,61,62,63].
Considering the wide variability observed for the breeding lines with high grain Zn and Fe, stability and G × E analyses are being applied for selecting promising RILs in cereals. The contribution of environmental variation to grain Zn and Fe along with other agronomic traits in a RIL population of sorghum was demonstrated through genotype · environment interaction, correlation, and GGE biplot analyses [64]. Stable RILs with higher grain Zn and Fe content were also identified in the RILs of pearl millet using AMMI and GGE biplot analyses [65]. Different stable breeding lines were identified for different environments among eight Zn biofortified lines through stability and G × E analyses [51].

3.3. Identification of Quantitative Trait Loci (QTLs)

Out of 13 QTLs identified for grain Fe and 10 QTLs for grain Zn in brown and polished rice using ICIM, six QTLs were found to be stable/consistent across the seasons, namely qDFF.2.2, qDFF.2.3, qZnPR.2.1, qZnPR.2.2, qZnBR.9.1, and qZnBR.9.2 on Chr02 and Chr09. The lack of stability of QTL effects across the environments/genetic backgrounds has been one of the most limiting factors for successful deployment of QTLs through MAS breeding for various complex traits [66,67,68,69]. The consistent QTLs (qZnPR.2.1, qZnPR.2.2) identified in this study can reliably be deployed in marker-assisted breeding for high grain Zn.
Several researchers have mapped QTLs for grain Zn and Fe traits in rice using various mapping populations, such as RILs, ILs, F2, DH, and MAGIC populations, and have identified multiple loci and clearly demonstrated the genetic complexity of the grain micronutrient traits. The genomic region on Chr02 (SNP_24248265-SNP_25682947) has shown co-localization in the qZPR.2.1 region [42].
Co-localization of QTLs for different element contents in grain is common in rice [70]. In the study, co-localization was observed: the genomic regions of Chr02, 05, and 06 of the DFF trait was co-localized with mineral trait QTLs FeBR, FePR, ZnBR, and ZnPR on Chr02. The Fe QTLs on Chr05 were co-located with ZnBR QTL, suggesting that there may be a possibility of selecting high grain Zn and Fe lines using molecular markers; Ref. [71] reported co-location of Zn and Fe QTLs on Chr02, 03, 08, and 12 in RIL populations.
In the present study, eight significant epistatic interactions were observed for PH, NT, FeBR, FePR, and ZnPR on Chr01, 02, 04, 05, and 07 and accounted for 11.20–26.83% PVE, suggesting the complex genetic variation of the traits. Two digenic interactions were detected for PH on Chr02, with PVE 8.74% and 24.49%, respectively. None of the identified digenic interactions in the present study were found to be involved with the main effect QTLs. Similar observations were earlier reported for epistatic interactions for grain Zn in rice [52,59]. Recent QTL mapping studies suggest that epistasis is considered to be a crucial genetic component underlying complex quantitative traits. Epistasis should be underscored in studying complex traits because it can account for hidden quantitative genetic variations in natural populations [72].
Integration of genomics with conventional breeding efforts can decipher the molecular mechanisms underlying traits of interest. Using the rice genome sequence, it is now possible to identify the putative candidate genes for grain Zn within QTL regions. In the present study, major QTLs involved in grain mineral element content that show consistency across seasons merit further examination by fine mapping and candidate gene analysis for ultimate use in MAS. The identified five candidate genes Os02t0629200 on Chr02 with Zn binding and transport activity, two genes Os09t0268300 and Os09t0297300 on Chr09 with transporter activity, and two genes Os06t0705700 and Os06t0706100 on Chr06 with Zn and Fe transporter activity and metal ion transport function could be promising genes for further validation. The enhanced expression of Os06t0705700 and Os06t0706100 in the flag leaf samples of the donor parent (Ranbir Basmati) and some of the donors/RILs with grain Zn further supported the identified QTLs/candidate genes in the present study.
Out of 92 RILs evaluated, one promising RIL G3 showed the highest SPY and Zn content across 92 RILs. The G3 RIL can be used as donor in biofortification breeding programs. Likewise, three other RILs, G1, G4 and G6, with high Zn content donors, and RILs G73, G49, G70, G82, and G52, with promising yield, can be used for future breeding programs.

4. Materials and Methods

4.1. Development of RIL Population

A recombinant inbred line (RIL) population (n = 306) was developed using the single seed descent (SSD) method by crossing MTU1010 (also known as Cotton Dora Sannalu, pedigree: Krishnaveni/IR-64; Zn: 16.6 ppm in brown rice; 12.9 ppm in polished rice; mean yield of 6.5 tons ha−1; widely cultivated mega variety) as a recipient parent for grain zinc (Zn) and Ranbir Basmati (pure line selection from Basmati 370-90-95; mean yield of 2–2.5 tons ha−1; Zn: 27.5 ppm in brown rice; 23.4 ppm in polished rice) as a donor parent. In this study, two parental lines and a set of selected 92 RILs (F8) were used for genotyping and phenotypic analysis. All 306 lines’ Fe and Zn values are included in Supplementary Table S4.

4.2. Field Experiment Details

A total of 92 RILs along with parents were evaluated during four seasons (Wet season (WS) 2017; Dry season (DS) 2018; Wet season (WS) 2020 and Dry Season (DS) 2020) at ICAR-Indian Institute of Rice Research (IIRR) farm, Hyderabad, India (17.53° N latitude and 78.27° E longitude with 545 mm mean rainfall). Each RIL was planted in four rows. with each row consisting of 15 plants with a spacing of 20 × 15 cm following randomized block design (RCBD). The crop was raised by following recommended field management practices of production and protection.

4.3. Phenotypic Evaluation of Agronomic and Mineral Traits

Across four environments (seasons), five uniform plants were tagged from the center of each row, and phenotype data was collected for the traits. Days to 50% flowering (DFF), plant height (PH) (cm), number of tillers per plant (NT), panicle length (PL) (cm) at heading/early filling stage, single plant yield (SPY) (g), grain zinc (ZnBR, ZnPR), and iron (FeBR, FePR) in brown and polished rice were estimated at the post-harvest stage.

4.4. Estimation of Grain Zn and Fe Content

For estimation of grain Zn and Fe in brown rice (ZnBR, FeBR) and polished rice (ZnPR, FePR), seed samples were dehusked using a JLGJ4.5 rice husker (Jingjian Huayuan International Trade Co., Ltd., Hangzhou, Zhejiang, China). The brown rice was polished using a polisher with non-ferrous and non-zinc components (Krishi International India Ltd., New Delhi, India). After a thorough cleaning, samples of brown and polished rice (5 g) were subjected to energy-dispersive XRF (ED-XRF) (OXFORD Instruments X-Supreme 8000, Abingdon, UK) as per the standardized protocol [73].

4.5. Statistical Analysis

For each season, descriptive statistics such as mean, standard error of the mean (SEM), skewness, kurtosis, and coefficient of variations (CV%) were calculated for five plants using MS-Excel (2010). Trait-wise frequency distribution histograms were generated using R software [74]. Pearson correlation analysis and principal component analysis were computed using R script among nine traits. Different R packages (R core Team, 2018, Vienna, Austria), namely ggplot2, gge, FactoMinerR, and factoextra, were used to generate frequency distribution plots and individual comparison of RILs; the ‘Metan’ package was used for ranking of RILs [75].

4.6. Genotyping by Sequencing (ddRAD-Seq)

Genomic DNA was extracted from pooled young leaves (20–25 days) of both parents and RILs using the DNeasy plant kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. DNA quality of each sample was checked on 0.8% agarose gel. Furthermore, DNA quantification was done using a Qubit 2.0 fluorometer (Themo Fisher Scientific Inc., Waltham, MA, USA). ddRAD–Seq protocol (modified GBS) was followed in the present study.
To perform GBS, genomic DNA was double digested with SphI and MlucI restriction enzymes (NEB, UK) and fractionated in 2% agarose gel to check the product size. The digested fragments were cleaned (Agencourt AMpure XP beads, Invitrogen, Waltham, MA, USA) using standard protocols. The ligation enzyme, T4 DNA ligase (NEB, England) was used to ligate the unique barcode adapters (4–8 nt sequence) at 16 °C for 30 min and heat-inactivated at 65 °C for 10 min.
This was followed by indexing with the addition of index-1 and index-2 (6–8 nt long) for multiplexing sequencing libraries in NGS Illumina. These libraries were PCR amplified (8–12 cycles) using a PhusionTM polymerase kit (Fisher Scientific, Loughborough, UK) followed by AMPure bead cleanup for purification to remove excess adapters and were sequenced on the Illumina HiSeqX platform (Illumina Inc, San Diego, CA, USA).

4.6.1. Genotyping and Filtration

The sequence reads for the parents and RILs were obtained as FASTQ files, which are used for SNP discovery. Raw reads were de-multiplexed according to their barcodes, and the adapters/barcode sequences were removed using standard software [76]. High-quality reads were aligned onto the rice reference genome of Nipponbare (MSU7) using bowtie2-2.2.6. The mapped reads were exported in the form of a Sequence Alignment Map (SAM) file by SAM tools, version 0.1.19 [77]. The alignment file was then processed for SNP calling using a Bayesian approach at specific site. Potential SNPs were filtered using the following criteria: loci with >70% missing data and those that showed distorted segregation of the two parental genotypes were excluded. SNPs with minimum allele frequency (MAF) > 0.05 and 90% call rate were considered for further analysis. The variant annotation was performed based on rice gene models, using in-house pipeline software developed by Agri Genome Labs Pvt., Ltd., Kochi, India.

4.6.2. Genetic Linkage Map Construction

High-quality SNPs obtained after filtering were used for map construction using the linkage mapping function using IciMapping v4.2. [78]. The grouping and ordering of 2746 SNP events between adjacent markers was performed at a minimum logarithm of odds (LOD) value of 2.5. ‘Rippling’ was done for fine-tuning of the ordered markers on their respective chromosomes. The Kosambi map function was used for the construction of the genetic map and to convert the recombination frequencies into map distances in centiMorgans (cM) [79].

4.6.3. QTL Analysis

The main effect QTL analysis was carried out using IciMapping. Deploying the BIP function in IciMapping, additive QTLs were identified by inclusive composite interval mapping (ICIM) based on a 1000—permutation test at a 95% confidence level. The QTLs with >3.0 LOD and phenotypic variance explained (PVE) > 10% were considered as major effect QTLs, those with PVE <10% were considered as minor effect QTLs for a particular trait. Epistatic interactions with the logarithm of odds (LOD) threshold value at 5.0 were considered as significant epistatic QTLs. The QTLs were visualized using MapChart V2.3 [80].

4.7. Mining of Candidate Genes

To identify the candidate genes within the identified QTL region, the genomic position of the flanking markers was used. Using the rice genome annotation project ((http://rice.plantbiology.msu.edu/ (22 August 2022)), RAP-DB genome browser (https://rapdb.dna.affrc.go.jp (22 August 2022)), and Q-TARO (QTL Annotation Rice Online) database (Yonemaru et al., 2010) [81] (http://qtaro.abr.affrc.go.jp/ (22 August 2022)), the candidate genes within the QTL regions were retrieved. Genes were functionally characterized into various categories using WEGO [82]. Network analysis of genes present in the major QTL genomic regions was derived using the Knetminer program (http://knetminer.rothamsted.ac.uk/Oryza_sativa/ (24 August 2022)). The temporal and spatial expression of the identified candidate genes were studied using RiceXpro (www.ricexpro.dna.affrc.go.jp (24 August 2022)).

4.8. Confirmation of Identified QTLs and Markers

Out of mined candidates from QTLs, primers were designed for two candidate genes with known function of Zn transport, namely Os06g0705700 (CS528–F primer: GTGACGGCTTCGGATGAG and R primer: CCCGGTGTAGAAGAAGGTAATG), Os06g0706100 (CS533–F primer: CCAATGCCTGGCCTACTT and R primer: CCACGTGGACACGTTCTT) using the Primer Quest tool (https://eu.idtdna.com/pages/tools/primerquest (2 September 2022)) with qPCR parameters (amplicon size < 150 bp; GC content > 50%; melting temp ~60 °C; no secondary structures) and were synthesized at Integrated DNA Technologies (USA). Based on the RiceXpro data, flag leaf samples during the anthesis stage (after 10 days of flowering) from 11 genotypes (with differential grain Zn) were collected in RNALater™ (Sigma Aldrich, St. Louis, MO, USA). The published protocols for RNA isolation, cDNA synthesis, and checking of integrity of RNA and cDNA and their concentration were followed [83]. To investigate the expression of Os06g0705700 and Os06g0706100 genes, qPCR was performed with three biological replicates. Each qPCR reaction was performed in three technical replicates containing 5 µL SYBR® qPCR Master Mix (Promega Corporation, Madison, WI, USA), 0.6 µM forward and reverse gene-specific primers, and 100 ng template cDNA. The qPCR was performed using QuantStudio5 Applied Biosystems (Thermo Fisher Scientific, Waltham, MA, USA) Real-Time PCR Detection System with the following conditions: 40 PCR cycles of denaturation at 95 °C for 15 s, amplification at 58 °C for 15 s, followed by extension at 72 °C for 10 s. The melting curve (melting—95 °C for 10 s, 65 °C for 1 min, 97 °C for 1 s) was analyzed for checking the amplicon specificity. Two internal control genes, Membrane Protein (Memp) (LOC_Os12g32950.1-F Primer: GAGCGC AAAGTTCCAGAAGAA and R Primer: CGCCACTAGTTGCCGTCCTGAT) and Tumor Protein Homolog (TPH) (LOC_Os11g43900.1-F Primer: CATTGGTGCCAACCCATC and R Primer: AAGGAGGTTGCTCCTGAAGA) were used for the normalization [84]. Relative fold change was calculated using the 2(−ΔΔC(T)) method as described by [85].

5. Conclusions

The developed genetic map in the present study using MTU1010 X Ranbir Basmati cloud be a foundation for fine mapping of grain zinc (Zn) QTLs, and the identified genomic regions could be utilized in future breeding programs. The genes present in the identified QTL regions in this study reportedly play a key role in transporter activity. Further fine mapping, sequencing, and functional validation is required to identify the candidate genes in these QTL regions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12010144/s1, Table S1: List of epistatic QTLs. Table S2: List of genes in selected QTL regions. Table S3: List of candidate genes with transporter activity. Table S4: Fe and Zn data of 306 RILs. Table S5: Analysis of Variance. Figure S1: Frequency distribution of 92 RILs for mineral-related traits for all the four seasons. Figure S2: Frequency distribution of 92 RILs for agronomic and yield-related traits for all the four seasons. Figure S3: WEGO analysis with the identified candidate genes in the QTL regions. Figure S4: Relative expression of (a) Os06g0705700 and (b) Os06g0706100 in the flag leaf of rice genotypes.

Author Contributions

C.N.N. and G.A.U. conceived the idea. C.N.N., D.S.R. and R.M.S. provided critical inputs to the concepts. G.A.U., K.S. and P.M.B. generated genotyping and phenotyping data. G.A.U. and. K.S. conducted phenotyping, genotyping and QTL analysis, including statistical analysis. G.A.U., V.J. and K.S. contributed to analyses, interpretation of data, and writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by DST (Department of Science and Technology); SERB (Science and Engineering Research Board) through DST-SERB National Post-Doctoral fellowship (PDF/2016/003924).

Data Availability Statement

Not applicable.

Acknowledgments

A.U. thanks the DST SERB (Department of Science and Technology-Science and Engineering Research Board), Government of India, for providing funding support through the DST SERB N-PDF scheme. The authors gratefully acknowledge the support extended by the ICAR CRP-Biofortification and ICAR-Indian Institute of Rice Research (ICAR-IIRR) for providing facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cakmak, I.; Kutman, U.B. Agronomic Biofortification of Cereals with Zinc: A Review. Eur. J. Soil Sci. 2018, 69, 172–180. [Google Scholar] [CrossRef] [Green Version]
  2. Garcia-Oliveira, A.L.; Chander, S.; Ortiz, R.; Menkir, A.; Gedil, M. Genetic Basis and Breeding Perspectives of Grain Iron and Zinc Enrichment in Cereals. Front. Plant Sci. 2018, 9, 937. [Google Scholar] [CrossRef] [PubMed]
  3. Sharma, V.; Saini, D.K.; Kumar, A.; Kaushik, P. A Review of Important QTLs for Biofortification Traits in Rice. Preprints 2019, 2019120158. [Google Scholar] [CrossRef] [Green Version]
  4. Roy Choudhury, A.; Review, M. Agronomic and Genetic Biofortification of Rice Grains with Microelements to Assure Human Nutritional Security. SF J. Agric. Crop Manag. 2020, 1, 1005. [Google Scholar]
  5. Welch, R.M.; Graham, R.D. Breeding for Micronutrients in Staple Food Crops from a Human Nutrition Perspective. J. Exp. Bot. 2004, 55, 353–364. [Google Scholar] [CrossRef] [Green Version]
  6. McCall, K.A.; Huang, C.C.; Fierke, C.A. Function and Mechanism of Zinc Metalloenzymes. J. Nutr. 2000, 130, 1437S–1446S. [Google Scholar] [CrossRef] [Green Version]
  7. Roohani, N.; Hurrell, R.; Kelishadi, R.; Schulin, R. Zinc and Its Importance for Human Health: An Integrative Review. J. Res. Med. Sci. 2013, 18, 144–157. [Google Scholar]
  8. Prasad, A.S. Discovery of Human Zinc Deficiency: Its Impact on Human Health and Disease. Adv. Nutr. 2013, 4, 176. [Google Scholar] [CrossRef] [Green Version]
  9. Gibson, R.S. Zinc: The Missing Link in Combating Micronutrient Malnutrition in Developing Countries. Proc. Nutr. Soc. 2006, 65, 51–60. [Google Scholar] [CrossRef] [Green Version]
  10. Waters, B.M.; Grusak, M.A. Quantitative Trait Locus Mapping for Seed Mineral Concentrations in Two Arabidopsis Thaliana Recombinant Inbred Populations. New Phytol. 2008, 179, 1033–1047. [Google Scholar] [CrossRef] [Green Version]
  11. Delseny, M.; Salses, J.; Cooke, R.; Sallaud, C.; Regad, F.; Lagoda, P.; Guiderdoni, E.; Ventelon, M.; Brugidou, C.; Ghesquière, A. Rice Genomics: Present and Future. Plant Physiol. Biochem. 2001, 39, 323–334. [Google Scholar] [CrossRef]
  12. Feng, Y.; Zhai, R.; Lin, Z.; Cao, L.; Wei, X.; Cheng, S. QTL Analysis for Yield Traits in Rice under Two Nitrogen Levels. Chin. J. Rice Sci. 2013, 27, 577–584. [Google Scholar]
  13. Khush, G.S. Strategies for Increasing the Yield Potential of Cereals: Case of Rice as an Example. Plant Breed. 2013, 132, 433–436. [Google Scholar] [CrossRef]
  14. Pradhan, S.K.; Pandit, E.; Pawar, S.; Naveenkumar, R.; Barik, S.R.; Mohanty, S.P.; Nayak, D.K.; Ghritlahre, S.K.; Sanjiba Rao, D.; Reddy, J.N.; et al. Linkage Disequilibrium Mapping for Grain Fe and Zn Enhancing QTLs Useful for Nutrient Dense Rice Breeding. BMC Plant Biol. 2020, 20, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sharma, A.; Patni, B.; Shankhdhar, D.; Shankhdhar, S.C. Zinc—An Indispensable Micronutrient. Physiol. Mol. Biol. Plants 2013, 19, 11–20. [Google Scholar] [CrossRef]
  16. Bouis, H.E.; Saltzman, A. Improving Nutrition through Biofortification: A Review of Evidence from HarvestPlus, 2003 through 2016. Glob. Food Secur. 2017, 12, 49–58. [Google Scholar] [CrossRef]
  17. Bouis, H.E.; Welch, R.M. Biofortification—A Sustainable Agricultural Strategy for Reducing Micronutrient Malnutrition in the Global South. Crop Sci. 2010, 50, S-20–S-32. [Google Scholar] [CrossRef] [Green Version]
  18. Graham, R.; Senadhira, D.; Beebe, S.; Iglesias, C.; Monasterio, I. Breeding for micronutrient density in edible portions of staple food crops: Conventional approaches. Field Crop. Res. 1999, 60, 57–80. [Google Scholar] [CrossRef]
  19. Zhao, M.; Lin, Y.; Chen, H. Improving Nutritional Quality of Rice for Human Health. Theor. Appl. Genet. 2020, 133, 1397–1413. [Google Scholar] [CrossRef]
  20. Wissuwa, M.; Ismail, A.M.; Graham, R.D. Rice Grain Zinc Concentrations as Affected by Genotype, Native Soil-Zinc Availability, and Zinc Fertilization. Plant Soil 2008, 306, 37–48. [Google Scholar] [CrossRef]
  21. Zhang, M.; Pinson, S.R.M.; Tarpley, L.; Huang, X.-Y.; Lahner, B.; Yakubova, E.; Baxter, I.; Guerinot, M.L.; Salt, D. E Mapping and Validation of Quantitative Trait Loci Associated with Concentrations of 16 Elements in Unmilled Rice Grain. Theor. Appl. Genet. 2014, 127, 137–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Swamy, B.P.M.; Descalsota, G.I.L.; Nha, C.T.; Amparado, A.; Inabangan-Asilo, M.A.; Manito, C.; Tesoro, F.; Reinke, R. Identification of Genomic Regions Associated with Agronomic and Biofortification Traits in DH Populations of Rice. PLoS ONE 2018, 13, e0201756. [Google Scholar] [CrossRef] [PubMed]
  23. Neeraja, C.N.; Babu, V.R.; Ram, S.; Hossain, F.; Hariprasanna, K.; Rajpurohit, B.S. Biofortification in Cereals: Progress and Prospects. Curr. Sci. 2017, 113, 1050–1057. [Google Scholar] [CrossRef]
  24. Collard, B.C.Y.; Jahufer, M.Z.Z.; Brouwer, J.B.; Pang, E.C.K. An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: The basic concepts. Euphytica 2005, 142, 169–196. [Google Scholar] [CrossRef]
  25. Calayugan, M.I.C.; Swamy, B.P.; Nha, C.T.; Palanog, A.D.; Biswas, P.S.; Descalsota-Empleo, G.I.; Inabangan-Asilo, M. A Zinc-Biofortified Rice: A Sustainable Food-Based Product for Fighting Zinc Malnutrition. In Rice Improvement; Springer: Cham, Switzerland, 2021. [Google Scholar] [CrossRef]
  26. Kumar, S.; Hash, C.T.; Nepolean, T.; Mahendrakar, M.D.; Satyavathi, C.T.; Singh, G.; Rathore, A.; Yadav, R.S.; Gupta, R.; Srivastava, R.K. Mapping Grain Iron and Zinc Content Quantitative Trait Loci in an Iniadi-Derived Immortal Population of Pearl Millet. Genes 2018, 9, 248. [Google Scholar] [CrossRef] [Green Version]
  27. Ott., J.; Wang, J.; Leal, S.M. Genetic Linkage Analysis in the Age of Whole-Genome Sequencing. Nat. Rev. Genet. 2015, 16, 275–284. [Google Scholar] [CrossRef]
  28. de Leon, T.B.; Linscombe, S.; Subudhi, P.K. Molecular Dissection of Seedling Salinity Tolerance in Rice (Oryza sativa L.) Using a High-Density GBS-Based SNP Linkage Map. Rice 2016, 9, 1–22. [Google Scholar] [CrossRef] [Green Version]
  29. Varshney, R.K.; Singh, V.K.; Kumar, A.; Powell, W.; Sorrells, M.E. Can Genomics Deliver Climate-Change Ready Crops? Curr. Opin. Plant Biol. 2018, 45, 205–211. [Google Scholar] [CrossRef]
  30. Davey, J.W.; Hohenlohe, P.A.; Etter, P.D.; Boone, J.Q.; Catchen, J.M.; Blaxter, M. Genome-Wide Genetic Marker Discovery and Genotyping Using next-Generation Sequencing. Nat. Rev. Genet. 2011, 12, 499–510. [Google Scholar] [CrossRef]
  31. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  32. Peterson, B.K.; Weber, J.N.; Kay, E.H.; Fisher, H.S.; Hoekstra, H.E. Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species. PLoS ONE 2012, 7, 37135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Sun, X.; Liu, D.; Zhang, X.; Li, W.; Liu, H.; Hong, W.; Jiang, C.; Guan, N.; Ma, C.; Zeng, H.; et al. SLAF-Seq: An Efficient Method of Large-Scale De Novo SNP Discovery and Genotyping Using High-Throughput Sequencing. PLoS ONE 2013, 8, e58700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gali, K.K.; Sackville, A.; Tafesse, E.G.; Lachagari, V.B.R.; McPhee, K.; Hybl, M.; Mikić, A.; Smýkal, P.; McGee, R.; Burstin, J.; et al. Genome-Wide Association Mapping for Agronomic and Seed Quality Traits of Field Pea (Pisum sativum L.). Front. Plant Sci. 2019, 10, 1538. [Google Scholar] [CrossRef] [PubMed]
  35. Dissanayaka, D.N.; Gali, K.K.; Jha, A.B.; Lachagari, V.B.R.; Warkentin, T.D. Genome-Wide Association Study to Identify Single Nucleotide Polymorphisms Associated with Fe, Zn, and Se Concentration in Field Pea. Crop Sci. 2020, 60, 2070–2084. [Google Scholar] [CrossRef] [Green Version]
  36. Furuta, T.; Ashikari, M.; Jena, K.K.; Doi, K.; Reuscher, S. Adapting Genotyping-by-Sequencing for Rice F2 Populations. G3 Genes Genomes Genet. 2017, 7, 881–893. [Google Scholar] [CrossRef] [Green Version]
  37. Bhatia, D.; Wing, R.A.; Yu, Y.; Chougule, K.; Kudrna, D.; Lee, S.; Rang, A.; Singh, K. Genotyping by Sequencing of Rice Interspecific Backcross Inbred Lines Identifies QTLs for Grain Weight and Grain Length. Euphytica 2018, 214, 41. [Google Scholar] [CrossRef]
  38. Yadav, S.; Sandhu, N.; Singh, V.K.; Catolos, M.; Kumar, A. Genotyping-by-Sequencing Based QTL Mapping for Rice Grain Yield under Reproductive Stage Drought Stress Tolerance. Sci. Rep. 2019, 9, 14326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. White, P.J.; Broadley, M.R. Physiological Limits to Zinc Biofortification of Edible Crops. Front. Plant Sci. 2011, 2, 80. [Google Scholar] [CrossRef] [Green Version]
  40. Nakandalage, N.; Nicolas, M.; Norton, R.M.; Hirotsu, N.; Milham, P.J.; Seneweera, S. Improving Rice Zinc Biofortification Success Rates through Genetic and Crop Management Approaches in a Changing Environment. Front. Plant Sci. 2016, 7, 764. [Google Scholar] [CrossRef] [Green Version]
  41. Khan, J.A.; Narayana, K.K.; Holla, S.; Shrinivas, S.M.; Dar, Z.A.; Shashidhar, H.E. Micronutrient Productivity: A Comprehensive Parameter for Biofortification in Rice (Oryza sativa L.) Grain. J. Sci. Food Agric. 2019, 99, 1311–1321. [Google Scholar] [CrossRef]
  42. Suman, K.; Neeraja, C.N.; Madhubabu, P.; Rathod, S.; Bej, S.; Jadhav, K.P.; Kumar, J.A.; Chaitanya, U.; Pawar, S.C.; Rani, S.H.; et al. Identification of Promising RILs for High Grain Zinc Through Genotype × Environment Analysis and Stable Grain Zinc QTL Using SSRs and SNPs in Rice (Oryza sativa L.). Front. Plant Sci. 2021, 12, 587482. [Google Scholar] [CrossRef] [PubMed]
  43. Janila, P.; Pandey, M.K.; Shasidhar, Y.; Variath, M.T.; Sriswathi, M.; Khera, P.; Manohar, S.S.; Nagesh, P.; Vishwakarma, M.K.; Mishra, G.P.; et al. Molecular Breeding for Introgression of Fatty Acid Desaturase Mutant Alleles (AhFAD2A and AhFAD2B) Enhances Oil Quality in High and Low Oil Containing Peanut Genotypes. Plant Sci. 2016, 242, 203–213. [Google Scholar] [CrossRef]
  44. Varshney, R.K.; Terauchi, R.; McCouch, S.R. Harvesting the Promising Fruits of Genomics: Applying Genome Sequencing Technologies to Crop Breeding. PLoS Biol. 2014, 12, e1001883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Collard, B.C.Y.; Mackill, D.J. Marker-Assisted Selection: An Approach for Precision Plant Breeding in the Twenty-First Century. Philos. Trans. R. Soc. B Biol. Sci. 2007, 363, 557–572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Soren, K.R.; Madugula, P.; Kumar, N.; Barmukh, R.; Sengar, M.S.; Bharadwaj, C.; Sharma, P.C.; Singh, S.; Bhandari, A.; Singh, J.; et al. Genetic Dissection and Identification of Candidate Genes for Salinity Tolerance Using Axiom® cicerSNP Array in Chickpea. Int. J. Mol. Sci. 2020, 21, 5058. [Google Scholar] [CrossRef] [PubMed]
  47. Babu, P.M.; Neeraja, C.N.; Rathod, S.; Suman, K.; Uttam, G.A.; Chakravartty, N.; Lachagari, V.B.R.; Chaitanya, U.; Rao, L.V.S.; Voleti, S.R. Stable SNP Allele Associations with High Grain Zinc Content in Polished Rice (Oryza sativa L.) Identified Based on DdRAD Sequencing. Front. Genet. 2020, 11, 763. [Google Scholar] [CrossRef]
  48. Heffner, E.L.; Sorrells, M.E.; Jannink, J.-L. Genomic Selection for Crop Improvement. Crop. Sci. 2009, 49, 1–12. [Google Scholar] [CrossRef]
  49. Sudan, J.; Singh, R.; Sharma, S.; Salgotra, R.K.; Sharma, V.; Singh, G.; Sharma, I.; Sharma, S.; Gupta, S.K.; Zargar, S.M. ddRAD sequencing-based identification of inter-genepool SNPs and association analysis in Brassica juncea. BMC Plant Biol. 2019, 19, 594. [Google Scholar] [CrossRef] [Green Version]
  50. Anuradha, K.; Agarwal, S.; Batchu, A.K.; Prasad, B.A.; Mallikarjuna, S.B.P.; Longvah, T. Evaluating rice germplasm for iron and zinc concentration in brown rice and seed dimensions. J. Phytol. 2012, 4, 19–25. [Google Scholar]
  51. Inabangan-Asilo, M.A.; Mallikarjuna Swamy, B.P.; Amparado, A.F.; Descalsota-Empleo, G.I.L.; Arocena, E.C.; Reinke, R. Stability and G × E Analysis of Zinc-Biofortified Rice Genotypes Evaluated in Diverse Environments. Euphytica 2019, 215, 61. [Google Scholar] [CrossRef] [Green Version]
  52. Descalsota-Empleo, G.I.; Amparado, A.; Inabangan-Asilo, M.A.; Tesoro, F.; Stangoulis, J.; Reinke, R.; Swamy, B.P.M. Genetic Mapping of QTL for Agronomic Traits and Grain Mineral Elements in Rice. Crop. J. 2019, 7, 560–572. [Google Scholar] [CrossRef]
  53. Dixit, S.; Singh, U.M.; Abbai, R.; Ram, T.; Singh, V.K.; Paul, A.; Virk, P.S.; Kumar, A. Identification of genomic region(s) responsible for high iron and zinc content in rice. Sci. Rep. 2019, 9, 8136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Anandan, A.; Rajiv, G.; Eswaran, R.; Prakash, M. Genotypic Variation and Relationships between Quality Traits and Trace Elements in Traditional and Improved Rice (Oryza Sativa L.) Genotypes. J. Food Sci. 2011, 76, H122–H130. [Google Scholar] [CrossRef]
  55. Descalsota, G.I.L.; Swamy, B.P.M.; Zaw, H.; Inabangan-Asilo, M.A.; Amparado, A.; Mauleon, R.; Chadha-Mohanty, P.; Arocena, E.C.; Raghavan, C.; Leung, H.; et al. Genome-Wide Association Mapping in a Rice Magic plus Population Detects Qtls and Genes Useful for Biofortification. Front. Plant Sci. 2018, 9, 1347. [Google Scholar] [CrossRef] [PubMed]
  56. Xu, Q.; Zheng, T.Q.; Hu, X.; Cheng, L.R.; Xu, J.L.; Shi, Y.M.; Li, Z.K. Examining Two Sets of Introgression Lines in Rice (Oryza sativa L.) Reveals Favorable Alleles That Improve Grain Zn and Fe Concentrations. PLoS ONE 2015, 10, e0131846. [Google Scholar] [CrossRef] [PubMed]
  57. Stangoulis, J.C.R.; Huynh, B.L.; Welch, R.M.; Choi, E.Y.; Graham, R.D. Quantitative Trait Loci for Phytate in Rice Grain and Their Relationship with Grain Micronutrient Content. Euphytica 2007, 154, 289–294. [Google Scholar] [CrossRef]
  58. Swamy, B.P.M.; Kaladhar, K.; Anuradha, K.; Batchu, A.K.; Longvah, T.; Sarla, N. QTL Analysis for Grain Iron and Zinc Concentrations in Two O. Nivara Derived Backcross Populations. Rice Sci. 2018, 25, 197–207. [Google Scholar] [CrossRef]
  59. Norton, G.J.; Deacon, C.M.; Xiong, L.; Huang, S.; Meharg, A.A.; Price, A.H. Genetic Mapping of the Rice Ionome in Leaves and Grain: Identification of QTLs for 17 Elements Including Arsenic, Cadmium, Iron and Selenium. Plant Soil 2010, 329, 139–153. [Google Scholar] [CrossRef]
  60. Ogasawara, M.; Miyazaki, N.; Monden, G.; Taniko, K.; Lim, S.; Iwata, M.; Ishii, T.; Ma, J.F.; Ishikawa, R. Role of QGZn9a in Controlling Grain Zinc Concentration in Rice, Oryza sativa L. Theor. Appl. Genet. 2021, 134, 3013–3022. [Google Scholar] [CrossRef]
  61. Naik, S.M.; Raman, A.K.; Nagamallika, M.; Venkateshwarlu, C.; Singh, S.P.; Kumar, S.; Singh, S.K.; Tomizuddin, A.; Das, S.P.; Prasad, K.; et al. Genotype  ×  Environment Interactions for Grain Iron and Zinc Content in Rice. J. Sci. Food Agric. 2020, 100, 4150–4164. [Google Scholar] [CrossRef]
  62. Ajmera, S. Studies on Stability Analysis for Grain Yield and Its Attributes in Rice (Oryza sativa L.) Genotypes. Int. J. Pure Appl. Biosci. 2017, 5, 892–908. [Google Scholar] [CrossRef]
  63. Suwarto; Nasrullah. Genotype × Environment Interaction for Iron Concentration of Rice in Central Java of Indonesia. Rice Sci. 2011, 18, 75–78. [Google Scholar] [CrossRef]
  64. Phuke, R.M.; Anuradha, K.; Radhika, K.; Jabeen, F.; Anuradha, G.; Ramesh, T.; Hariprasanna, K.; Mehtre, S.P.; Deshpande, S.P.; Anil, G.; et al. Genetic Variability, Genotype × Environment Interaction, Correlation, and GGE Biplot Analysis for Grain Iron and Zinc Concentration and Other Agronomic Traits in RIL Population of Sorghum (Sorghum bicolor L. Moench). Front. Plant Sci. 2017, 8, 712. [Google Scholar] [CrossRef] [PubMed]
  65. Singhal, T.; Satyavathi, C.T.; Kumar, A.; Sankar, S.M.; Singh, S.P.; Bharadwaj, C.; Aravind, J.; Anuradha, N.; Meena, M.C.; Singh, N.; et al. Genotype × Environment Interaction and Genetic Association of Grain Iron and Zinc Content with Other Agronomic Traits in RIL Population of Pearl Millet. Crop Pasture Sci. 2018, 69, 1092–1102. [Google Scholar] [CrossRef]
  66. Courtois, B.; Shen, L.; Petalcorin, W.; Carandang, S.; Mauleon, R.; Li, Z. Locating QTLs Controlling Constitutive Root Traits in the Rice Population IAC 165 × Co39. Euphytica 2003, 134, 335–345. [Google Scholar] [CrossRef]
  67. Lafitte, H.R.; Price, A.H.; Courtois, B. Yield Response to Water Deficit in an Upland Rice Mapping Population: Associations among Traits and Genetic Markers. Theor. Appl. Genet. 2004, 109, 1237–1246. [Google Scholar] [CrossRef]
  68. Bernier, J.; Atlin, G.N.; Serraj, R.; Kumar, A.; Spaner, D. Breeding upland rice for drought resistance. J. Sci. Food Agric. 2008, 88, 927–939. [Google Scholar] [CrossRef]
  69. Kumar, A.; Dixit, S.; Ram, T.; Yadaw, R.B.; Mishra, K.K.; Mandal, N.P. Breeding High-Yielding Drought-Tolerant Rice: Genetic Variations and Conventional and Molecular Approaches. J. Exp. Bot. 2014, 65, 6265–6278. [Google Scholar] [CrossRef] [Green Version]
  70. Huang, Y.; Sun, C.; Min, J.; Chen, Y.; Tong, C.; Bao, J. Association Mapping of Quantitative Trait Loci for Mineral Element Contents in Whole Grain Rice (Oryza sativa L.). J. Agric. Food Chem. 2015, 63, 10885–10892. [Google Scholar] [CrossRef]
  71. Swamy, B.P.M.; Kaladhar, K.; Anuradha, K.; Batchu, A.K.; Longvah, T.; Viraktamath, B.C.; Sarla, N. Enhancing iron and zinc concentration in rice grains using wild species. In ADNAT Convention and International Symposium on Genomics and Biodiversity; CCMB: Hyderabad, India, 2011; p. 71. [Google Scholar]
  72. Mackay, T.F.C. Epistasis and quantitative traits: Using model organisms to study gene–gene interactions. Nat. Rev. Genet. 2013, 15, 22–33. [Google Scholar] [CrossRef] [Green Version]
  73. Rao, D.S.; Babu, P.M.; Swarnalatha, P.; Kota, S.; Bhadana, V.P.; Varaprasad, G.S.; Babu, V.R. Assessment of Grain Zinc and Iron Variability in Rice Germplasm Using Energy Dispersive X-ray Fluorescence Spectrophotometer. J. Rice Res. 2014, 7, 45–52. [Google Scholar]
  74. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  75. Yaseen, M.; Eskridge, K.M.; Murtaza, G. Stability: Stability Analysis of Genotype by Environment Interaction (GEI), R package version 0.5.0; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  76. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for Association Mapping of Complex Traits in Diverse Samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef] [PubMed]
  77. 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] [Green Version]
  78. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated Software for Genetic Linkage Map Construction and Quantitative Trait Locus Mapping in Biparental Populations. Crop J. 2015, 3, 269–283. [Google Scholar] [CrossRef]
  79. Kosambi, D.D. The estimation of map distances from recombination values. Ann. Eugen. 1943, 12, 172–175. [Google Scholar] [CrossRef]
  80. Voorrips, R.E. Computer Note MapChart: Software for the Graphical Presentation of Linkage Maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Yonemaru, J.-I.; Yamamoto, T.; Fukuoka, S.; Uga, Y.; Hori, K.; Yano, M. Q-TARO: QTL Annotation Rice Online Database. Rice 2010, 3, 194–203. [Google Scholar] [CrossRef] [Green Version]
  82. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  83. Bej, S.; Neeraja, C.; Krishna Kanth, T.; Suman, K.; Barbadikar, K.M.; Voleti, S. Correlation of Expressional Pattern of Ubiquitin Activating Gene with Grain Fe Content in Rice. Oryza-Int. J. Rice 2020, 57, 251–259. [Google Scholar] [CrossRef]
  84. Phule, A.S.; Barbadikar, K.M.; Madhav, M.S.; Senguttuvel, P.; Babu, M.B.B.P.; Ananda Kumar, P. Genes Encoding Membrane Proteins Showed Stable Expression in Rice under Aerobic Condition: Novel Set of Reference Genes for Expression Studies. 3 Biotech 2018, 8, 383. [Google Scholar] [CrossRef]
  85. Schmittgen, T.D.; Livak, K.J. Analyzing Real-Time PCR Data by the Comparative CT Method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Individual comparison of 92 RILs with respect to zinc (Zn) content and single plant yield (SPY) in four seasons; (b) mean zinc (Zn) content vs. mean single plant yield (SPY) across four seasons.
Figure 1. (a) Individual comparison of 92 RILs with respect to zinc (Zn) content and single plant yield (SPY) in four seasons; (b) mean zinc (Zn) content vs. mean single plant yield (SPY) across four seasons.
Plants 12 00144 g001
Figure 2. PCA of 92 RILs for nine traits based on mineral, yield, and agronomic traits.
Figure 2. PCA of 92 RILs for nine traits based on mineral, yield, and agronomic traits.
Plants 12 00144 g002
Figure 3. Ranking of genotypes across all seasons.
Figure 3. Ranking of genotypes across all seasons.
Plants 12 00144 g003
Figure 4. Major quantitative trait loci (QTLs) for various mineral, agronomic, and yield-related traits identified in the MTU1010 X Ranbir Basmati RIL population.
Figure 4. Major quantitative trait loci (QTLs) for various mineral, agronomic, and yield-related traits identified in the MTU1010 X Ranbir Basmati RIL population.
Plants 12 00144 g004aPlants 12 00144 g004b
Figure 5. Epistatic interactions in MTU1010 X Ranbir Basmati RIL population using IciMapping v4.2.
Figure 5. Epistatic interactions in MTU1010 X Ranbir Basmati RIL population using IciMapping v4.2.
Plants 12 00144 g005
Figure 6. (a) Network analysis of candidate gene (Os06g0705700) in a major QTL qZnPR.6.2 using Knetminer analysis; (b) Network analysis of candidate gene (Os06g0706100) in a major QTL qZnPR.6.2 using Knetminer analysis.
Figure 6. (a) Network analysis of candidate gene (Os06g0705700) in a major QTL qZnPR.6.2 using Knetminer analysis; (b) Network analysis of candidate gene (Os06g0706100) in a major QTL qZnPR.6.2 using Knetminer analysis.
Plants 12 00144 g006
Table 1. Descriptive statistics of selected RILs (four seasons).
Table 1. Descriptive statistics of selected RILs (four seasons).
StatisticsYearDFFPHNTPLSPYZnBRZnPRFeBRFePR
MeanWS17105.9126.611.226.520.622.2189.12.8
MTU101097921221.130.2112.46.810.11.9
R. Basmati921151819.721.117.213.311.62.9
DS1899.4110.814.622.417.416.7139.42
MTU101096901322.230.711.76.710.11.9
R. Basmati851141919.414.616.612.511.62.9
WS20104.6119.111.925.617.217.413.89.82.6
MTU1010100103122514.520.114.611.82.9
R. Basmati92101111813.619.116.313.52.3
DS2098.1107.41322.412.816.713.410.73.1
MTU1010103109122313.817.210.98.11.1
R. Basmati100101142215.515.212.610.42.2
VarianceWS1760.9247.14.55108.627.725.12.41.9
DS18236.7327.816.9548.416.715.93.40.9
WS2032.3185.43.24.43916.717.32.10.7
DS2062.467.67.44.811.613.3131.22.4
Std Error (Dev)WS177.815.72.12.210.45.351.61.4
DS1815.418.14.12.274.141.80.9
WS205.713.61.82.16.24.14.21.50.8
DS207.98.22.72.23.43.73.61.11.6
SkewnessWS170.4−0.40.6−0.30.60.60.70.32.4
DS18−0.4−1.60.2−0.20.81.11−0.90.9
WS200.5−0.10.2−0.81.21.21.90.61
DS200.40.90.2010.20.30.45.5
KurtosisWS17−1−0.80.80.9−0.3−0.40−0.59.5
DS181.39.5−0.10.21.23.42.83.31.1
WS200.8−1−0.31.72.92.97.3−0.21
DS20−0.30.50.2−0.22.41.51.4140.3
MinimumWS179289719.76.512.46.86.41
DS183787163.68.55.31.20.4
WS2092927187.57.25.67.41.4
DS2080927176.27.34.58.11.1
MaximumWS1712315718335137.831.713.210.6
DS18131153272742.935.230.513.35.1
WS2012514716304532.935.613.55.7
DS20118135222725.928.426.314.315.5
RangeWS1731681113.344.525.424.96.89.6
DS1894145201139.326.725.212.14.7
WS20335591237.525.7306.14.3
DS203843151019.721.121.86.214.4
Standard ErrorWS170.81.60.220.231.080.540.520.160.14
DS181.61.90.420.230.720.420.410.190.1
WS200.61.40.190.210.640.420.430.150.09
DS200.80.80.280.220.350.380.370.110.16
CV %WS171.63.20.430.452.131.081.030.320.28
DS183.23.70.840.451.420.840.820.380.19
WS201.22.80.370.421.280.840.850.30.17
DS201.61.70.550.440.70.750.740.230.32
WS: Wet Season; DS: Dry Season; DFF: Days to 50% Flowering (days); PH: Plant Height (cm); NT: No. of Tillers; PL: Panicle Length (cm); SPY: Single Plant Yield (g); ZnBR: Zinc Content in Brown Rice; ZnPR (ppm): Zinc Content in Polished Rice (ppm); FeBR: Iron Content in Brown Rice (ppm); FePR: Iron Content in Polished Rice (ppm).
Table 2. Pearson correlation analysis for agronomic, yield-, and mineral-related traits evaluated during two wet seasons (WS) and two dry seasons (DS) (average).
Table 2. Pearson correlation analysis for agronomic, yield-, and mineral-related traits evaluated during two wet seasons (WS) and two dry seasons (DS) (average).
DFFPHNTPLSPYZnBRZnPRFeBRFePR
DFF1.00−0.010.060.020.28 **0.040.02−0.140.06
1.00−0.09−0.120.03−0.01−0.06−0.06−0.35 ***−0.24 *
PH−0.011.000.100.25 *0.04−0.040.03−0.04−0.08
−0.091.000.25 *0.18−0.20 *0.000.050.11−0.05
NT0.060.101.00−0.30 **0.02−0.06−0.190.120.03
−0.120.25 *1.000.00−0.01−0.040.040.180.13
PL0.020.25 *−0.30 **1.000.060.020.15−0.26 *−0.04
0.030.180.001.00−0.16−0.04−0.07−0.020.03
SPY0.28 **0.040.020.061.00−0.26 *−0.21 **−0.27 **−0.19
−0.01−0.20 *−0.01−0.161.00−0.17−0.23 *−0.09−0.27 *
ZnBR0.04−0.04−0.060.02−0.26 *1.000.79 ***0.64 ***0.40 ***
−0.060.00−0.04−0.04−0.171.000.76 ***0.59 ***0.23 *
ZnPR0.020.03−0.190.15−0.21 **0.79 ***1.000.45 ***0.35 ***
−0.060.050.04−0.07−0.23 *0.76 ***1.000.47 ***0.43 ***
FeBR−0.14−0.040.12−0.26 *−0.270.64 ***0.45 ***1.000.34 ***
−0.35 ***0.110.18−0.02−0.090.59 ***0.47 ***1.000.10
FePR0.06−0.080.03−0.04−0.190.40 ***0.35 ***0.34 ***1.00
−0.24 *−0.050.130.03−0.27 *0.23 *0.43 ***0.101.00
*** p < 0.001; ** p < 0.01; * p < 0.05; DFF: Days to 50% Flowering; PH: Plant Height; NT: No. of Tillers; PL: Panicle Length; SPY: Single Plant Yield; ZnBR: Zinc Content in Brown Rice; ZnPR; Zinc Content in Polished Rice; FeBR: Iron Content in Brown Rice; FePR: Iron Content in Polished Rice. Upper value indicates wet season (WS), lower value indicates dry season (DS).
Table 3. Importance of components.
Table 3. Importance of components.
PC1PC2PC3PC4PC5PC6PC7PC8PC9
Standard Deviation1.791.201.100.980.910.790.680.480.23
Proportion of Variance0.350.160.130.110.090.070.050.030.01
Cumulative Proportion0.350.520.650.760.850.920.970.991.00
Table 4. Distribution of markers on the 12 linkage groups (LGs) of rice genetic map for the MTU1010 X Ranbir Basmati RIL population.
Table 4. Distribution of markers on the 12 linkage groups (LGs) of rice genetic map for the MTU1010 X Ranbir Basmati RIL population.
Linkage GrpNo. of SNPs MappedMap Distance (cM)Density (SNPs/cM)
LG012152001.08
LG021713710.46
LG032521941.30
LG04281220.23
LG052612381.10
LG064182721.54
LG072051861.10
LG082661541.73
LG091322650.50
LG103151552.03
LG111941471.32
LG122891412.05
Total274624451.2
Table 5. Summary of major and minor QTLs for mineral traits (zinc (Zn) and iron (Fe)).
Table 5. Summary of major and minor QTLs for mineral traits (zinc (Zn) and iron (Fe)).
S. No.Trait NameQTL NameChrMarker IntervalLODPVE (%)AddAllele
1ZnBR-WS17qZnBR.2.12SNP_25682947-242482654.5320.8−2.56P2
2ZnBR-WS20qZnBR.4.14SNP_34361396-348627893.8416.01−1.81P2
3ZnBR-DS18qZnBR.5.15SNP_23987433-239874402.5923.49−1.3P2
4ZnBR-WS20qZnBR.9.19SNP_3990952-40522672.7211.41.65P1
5ZnBR-DS20qZnBR.9.29SNP_3036200-93876762.5513.88−1.43P2
6ZnPR-WS17qZnPR.2.12SNP_25682947-242482654.7519.42−2.45P2
7ZnPR-WS20qZnPR.2.22SNP_35521280-242482653.4611.78−8.29P2
8ZnPR-WS17qZnPR.6.16SNP_28083826-28240312.8112.34−3.04P2
9ZnPR-DS18qZnPR.6.26SNP_28684145-308950323.198.33−2.19P2
10ZnPR-DS18qZnPR.7.17SNP_14229996-134953153.318.13.52P1
11FeBR-WS17qFeBR.2.12SNP_24248265-106866454.1115.73−0.61P2
12FeBR-WS17qFeBR.3.13SNP_31758461-319085094.0616.78−0.85P2
13FeBR-WS17qFeBR.6.16SNP_2824031-280335552.7311.27−0.76P2
14FeBR-WS20qFeBR.5.15SNP_25431969-261113608.3234.15−0.95P2
15FePR-WS17qFePR.1.11SNP_12900648-3597333210.3828.73−4.04P2
16FePR-DS18qFePR.2.12SNP_25682947-242482652.733.52−0.27P2
17FePR-DS20qFePR.3.13SNP_13716046-2827889126.5936−5.81P2
18FePR-WS17qFePR.5.15SNP_14288078-244232614.5611.97−0.68P2
19FePR-WS20qFePR.5.25SNP_28525048-278592062.918.680.69P1
20FePR-WS17qFePR.6.16SNP_2824031-280335552.8112.43−0.86P2
21FePR-DS18qFePR.8.18SNP_8033030-1061502711.5218.27−0.81P2
22FePR-DS18qFePR.10.110SNP_18122916-179945424.125.530.37P1
23FePR-DS18qFePR.12.112SNP_214104199-21757505713.1721.55−0.64P2
Chr: Chromosome; QTL: Quantitative Trait Loci; cM: Centimorgan; PVE: Phenotypic Variance Explained; LOD: Logarithm of Odds; Add: Additive Effect; P1: Parent 1 (MTU1010); P2: Parent 2 (Ranbir Basmati); ZnBR: Zinc Content in Brown Rice; ZnPR: Zinc Content in Polished Rice; FeBR: Iron Content in Brown Rice; FePR: Iron Content in Polished Rice; WS: Wet Season; DS: Dry Season.
Table 6. Summary of major and minor QTLs for agronomic and yield related traits.
Table 6. Summary of major and minor QTLs for agronomic and yield related traits.
S. No.TraitQTL NameChrMarker IntervalLODPVE (%)AddAllele
1DFF-DS18qDFF.1.11SNP_24168773-45436018.318.86−8.33P2
2DFF-WS17qDFF.2.12SNP_12479844-313587684.1716.054P1
3DFF-DS18qDFF.2.22SNP_18086957- 173527079.910.325.52P1
4DFF-DS20qDFF.2.32SNP_18086957-173527072.8212.252.46P1
5DFF-DS20qDFF.2.42SNP_24248265-106866452.6114.312.97P1
6DFF-DS18qDFF.2.52SNP_24248265-1068664511.6114.47.07P1
7DFF-DS18qDFF.3.13SNP_22060292-228802598.869.43−6.4P2
8DFF-WS17qDFF.6.16SNP_1521959-22114573.046.63−2.98P2
9DFF-WS17qDFF.6.26SNP_6569100-123105777.2621.35−4.68P2
10DFF-DS18qDFF.8.18SNP_14611920-1532988716.5420.9613.22P1
11DFF-WS20qDFF.11.111SNP_11928941-14213113.6819.25−3.6P2
12PL-WS17qPL.2.12SNP_13823910-152783682.5613.29−0.73P2
13PL-DS20qPL.9.19SNP_18030272-179302223.38.030.83P1
14PL-DS18qPL.11.111SNP_116504977-122940272.6512.391.47P1
15PL-DS20qPL.11.211SNP_125186904-1232205545.5116.51−1.26P2
16PL-WS20qPL.12.112SNP_21721652-217976993.215.880.83P1
17SPY-DS18qSPY.3.13SNP_32176877-9622762.6311.264.71P1
18SPY-WS17qSPY.6.16SNP_22324866-16917513.5715.32−4.98P2
19SPY-WS17qSPY.9.19SNP_18053950-180539492.8211.283.74P1
20SPY-DS20qSPY.9.29SNP_18915073-187452002.6912.991.23P1
21NT-DS20qNT.6.16SNP_14518832-156215712.5112.031.59P1
22NT-WS17qNT.9.19SNP_17649556-176546623.2417.94−0.81P2
23NT-WS20qNT.9.29SNP_17887929-187446233.3310.37−0.73P2
24NT-WS20qNT.10.110SNP_22330346-182485112.5711.021.04P1
Chr: Chromosome; QTL: Quantitative Trait Loci; cM: Centimorgan; PVE: Phenotypic Variance Explained; LOD: Logarithm Of Odds; Add: Additive Effect; P1: Parent 1 (MTU1010); P2: Parent 2 (Ranbir Basmati); DFF: Days to 50% Flowering; PL: Panicle Length; SPY: Single Plant Yield; NT: No. of Tillers; WS: Wet Season; DS: Dry Season.
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

Uttam, G.A.; Suman, K.; Jaldhani, V.; Babu, P.M.; Rao, D.S.; Sundaram, R.M.; Neeraja, C.N. Identification of Genomic Regions Associated with High Grain Zn Content in Polished Rice Using Genotyping-by-Sequencing (GBS). Plants 2023, 12, 144. https://doi.org/10.3390/plants12010144

AMA Style

Uttam GA, Suman K, Jaldhani V, Babu PM, Rao DS, Sundaram RM, Neeraja CN. Identification of Genomic Regions Associated with High Grain Zn Content in Polished Rice Using Genotyping-by-Sequencing (GBS). Plants. 2023; 12(1):144. https://doi.org/10.3390/plants12010144

Chicago/Turabian Style

Uttam, Goparaju Anurag, Karre Suman, Veerendra Jaldhani, Pulagam Madhu Babu, Durbha Sanjeeva Rao, Raman Meenakshi Sundaram, and Chirravuri Naga Neeraja. 2023. "Identification of Genomic Regions Associated with High Grain Zn Content in Polished Rice Using Genotyping-by-Sequencing (GBS)" Plants 12, no. 1: 144. https://doi.org/10.3390/plants12010144

APA Style

Uttam, G. A., Suman, K., Jaldhani, V., Babu, P. M., Rao, D. S., Sundaram, R. M., & Neeraja, C. N. (2023). Identification of Genomic Regions Associated with High Grain Zn Content in Polished Rice Using Genotyping-by-Sequencing (GBS). Plants, 12(1), 144. https://doi.org/10.3390/plants12010144

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