Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone).

MAIN CONCLUSION: The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Genome-wide identification and characterization of NBLRR genes in finger millet (Eleusine coracana L.) and their expression in response to Magnaporthe grisea infection.

BACKGROUND: The nucleotide binding site leucine rich repeat (NBLRR) genes significantly regulate defences against phytopathogens in plants. The genome-wide identification and analysis of NBLRR genes have been performed in several species. However, the detailed evolution, structure, expression of NBLRRs and functional response to Magnaporthe grisea are unknown in finger millet (Eleusine coracana (L.) Gaertn.). RESULTS: The genome-wide scanning of the finger millet genome resulted in 116 NBLRR (EcNBLRRs1-116) encompassing 64 CC-NB-LRR, 47 NB-LRR and 5 CCR-NB-LRR types. The evolutionary studies among the NBLRRs of five Gramineae species, viz., purple false brome (Brachypodium distachyon (L.) P.Beauv.), finger millet (E. coracana), rice (Oryza sativa L.), sorghum (Sorghum bicolor L. (Moench)) and foxtail millet (Setaria italica (L.) P.Beauv.) showed the evolution of NBLRRs in the ancestral lineage of the target species and subsequent divergence through gene-loss events. The purifying selection (Ka/Ks < 1) shaped the expansions of NBLRRs paralogs in finger millet and orthologs among the target Gramineae species. The promoter sequence analysis showed various stress- and phytohormone-responsive cis-acting elements besides growth and development, indicating their potential role in disease defence and regulatory mechanisms. The expression analysis of 22 EcNBLRRs in the genotypes showing contrasting responses to Magnaporthe grisea infection revealed four and five EcNBLRRs in early and late infection stages, respectively. The six of these nine candidate EcNBLRRs proteins, viz., EcNBLRR21, EcNBLRR26, EcNBLRR30, EcNBLRR45, EcNBLRR55 and EcNBLRR76 showed CC, NB and LRR domains, whereas the EcNBLRR23, EcNBLRR32 and EcNBLRR83 showed NB and LRR somains. CONCLUSION: The identification and expression analysis of EcNBLRRs showed the role of EcNBLRR genes in assigning blast resistance in finger millet. These results pave the foundation for in-depth and targeted functional analysis of EcNBLRRs through genome editing and transgenic approaches.

Genome-wide identification of potassium channels in maize showed evolutionary patterns and variable functional responses to abiotic stresses.

Potassium (K) channels are essential components of plant biology, mediating not only K ion (K+) homeostasis but also regulating several physiological processes and stress tolerance. In the current investigation, we identified 27 K+ channels in maize and deciphered the evolution and divergence pattern with four monocots and five dicot species. Chromosomal localization and expansion of K+ channel genes showed uneven distribution and were independent of genome size. The dispersed duplication is the major force in expanding K+ channels in the target genomes. The mean Ka/Ks ratio of <0.5 in paralogs and orthologs indicates horizontal and vertical expansions of K+ channel genes under strong purifying selection. The one-to-one K+ channel orthologs were prominent among the closely related species, with higher synteny between maize and the rest of the monocots. Comprehensive K+ channels promoter analysis revealed various cis-regulatory elements mediating stress tolerance with the predominance of MYB and STRE binding sites. The regulatory network showed AP2-EREBP TFs, miR164 and miR399 are prominent regulatory elements of K+ channels. The qRT-PCR analysis of K+ channels and regulatory miRNAs showed significant expressions in response to drought and waterlogging stresses. The present study expanded the knowledge on K+ channels in maize and will serve as a basis for an in-depth functional analysis.

Development of bread wheat (Triticum aestivum L) variety HD3411 following marker-assisted backcross breeding for drought tolerance.

Marker-assisted backcross breeding enables selective insertion of targeted traits into the genome to improve yield, quality, and stress resistance in wheat. In the current investigation, we transferred four drought tolerance quantitative trait loci (QTLs) controlling traits, viz canopy temperature, normalized difference vegetative index, chlorophyll content, and grain yield from the drought-tolerant donor line, C306, into a popular high-yielding, drought-sensitive variety, HD2733. Marker-assisted selection coupled with stringent phenotypic screening was used to advance each generation. This study resulted in 23 improved lines carrying combinations of four drought tolerance QTLs with a range of 85.35%-95.79% background recovery. The backcross-derived lines gave a higher yield under moisture-deficit stress conditions compared with the recipient parent. They also showed higher phenotypic mean values for physiological traits and stability characteristics of HD2733. A promising genotype, HD3411, derived from this cross was identified for release after national multi-location coordinating trials under the All India Coordinated Wheat Improvement Project. Our study is a prime example of the advantages of precision breeding using integrating markers and phenotypic selection to develop new cultivars with desirable traits like drought tolerance.

Unravelling structural, functional, evolutionary and genetic basis of SWEET transporters regulating abiotic stress tolerance in maize.

Sugars Will Eventually be Exported Transporters (SWEETs) are the novel sugar transporters widely distributed among living systems. SWEETs play a crucial role in various bio-physiological processes, viz., plant developmental, nectar secretion, pollen development, and regulation of biotic and abiotic stresses, in addition to their prime sugar-transporting activity. Thus, in-depth structural, evolutionary, and functional characterization of maize SWEET transporters was performed for their utility in maize improvement. The mining of SWEET genes in the latest maize genome release (v.5) showed an uneven distribution of 20 ZmSWEETs. The comprehensive structural analyses and docking of ZmSWEETs with four sugars, viz., fructose, galactose, glucose, and sucrose, revealed frequent amino acid residues forming hydrogen (asparagine, valine, serine) and hydrophobic (tryptophan, glycine, and phenylalanine) interactions. Evolutionary analyses of SWEETs showed a mixed lineage with 50-100 % commonality of ortho-groups and -sequences evolved under strong purifying selection (Ka/Ks < 0.5). The duplication analysis showed non-functionalization (ZmSWEET18 in B73) and neo- and sub-functionalization (ZmSWEET3, ZmSWEET6, ZmSWEET9, ZmSWEET19, and ZmSWEET20) events in maize. Functional analyses of ZmSWEET genes through co-expression, in silico expression and qRT-PCR assays showed the relevance of ZmSWEETs expression in regulating drought, heat, and waterlogging stress tolerances in maize. The first ever ZmSWEET-regulatory network revealed 286 direct (ZmSWEET-TF: 140 ZmSWEET-miRNA: 146) and 1226 indirect (TF-TF: 597; TF-miRNA: 629) edges. The present investigation has given new insights into the complex transcriptional and post-transcriptional regulation and the regulatory and functional relevance of ZmSWEETs in assigning stress tolerance in maize.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi.

Zinc binuclear cluster proteins (ZBC; Zn(II)2Cys6) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)2Cys6 coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)2Cys6 genes. The perseverance of GC-rich codons of Zn(II)2Cys6 in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

Evolutionary and functional characterisation of glutathione peroxidases showed splicing mediated stress responses in Maize.

Maize (Zea mays L) is an important cereal with extensive adaptability and multifaceted usages. However, various abiotic and biotic stresses limit the productivity of maize across the globe. Exposure of plant to stresses disturb the balance between reactive oxygen species (ROS) production and scavenging, which subsequently increases cellular damage and death of plants. Tolerant genotypes have evolved higher output of scavenging antioxidative defence compounds (ADCs) during stresses as one of the protective mechanisms. The glutathione peroxidases (GPXs) are the broad class of ADCs family. The plant GPXs catalyse the reduction of hydrogen peroxide (H2O2), lipid hydroperoxides and organic hydroperoxides to the corresponding alcohol, and facilitate the regulation of stress tolerance mechanisms. The present investigation was framed to study the maize GPXs using evolutionary and functional analyses. Seven GPX genes with thirteen splice-variants and sixty-three types of cis-acting elements were identified through whole-genome scanning in maize. Evolutionary analysis of GPXs in monocots and dicots revealed mixed and lineage-specific grouping patterns in phylogeny. The expression of ZmGPX splice variants was studied in drought and waterlogging tolerant (L1621701) and sensitive (PML10) genotypes in root and shoot tissues. Further, the differential expression of splice variants of ZmGPX1, ZmGPX3, ZmGPX6 and ZmGPX7 and regulatory network analysis suggested the splicing and regulatory elements mediated stress responses. The present investigation suggests targeting the splicing machinery of GPXs as an approach to enhance the stress tolerance in maize.

Comparative Transcriptome Analysis of Iron and Zinc Deficiency in Maize (Zea mays L.).

Globally, one-third of the population is affected by iron (Fe) and zinc (Zn) deficiency, which is severe in developing and underdeveloped countries where cereal-based diets predominate. The genetic biofortification approach is the most sustainable and one of the cost-effective ways to address Fe and Zn malnutrition. Maize is a major source of nutrition in sub-Saharan Africa, South Asia and Latin America. Understanding systems' biology and the identification of genes involved in Fe and Zn homeostasis facilitate the development of Fe- and Zn-enriched maize. We conducted a genome-wide transcriptome assay in maize inbred SKV616, under -Zn, -Fe and -Fe-Zn stresses. The results revealed the differential expression of several genes related to the mugineic acid pathway, metal transporters, photosynthesis, phytohormone and carbohydrate metabolism. We report here Fe and Zn deficiency-mediated changes in the transcriptome, root length, stomatal conductance, transpiration rate and reduced rate of photosynthesis. Furthermore, the presence of multiple regulatory elements and/or the co-factor nature of Fe and Zn in enzymes indicate their association with the differential expression and opposite regulation of several key gene(s). The differentially expressed candidate genes in the present investigation would help in breeding for Fe and Zn efficient and kernel Fe- and Zn-rich maize cultivars through gene editing, transgenics and molecular breeding.

RNAseq revealed the important gene pathways controlling adaptive mechanisms under waterlogged stress in maize.

Waterlogging causes yield penalty in maize-growing countries of subtropical regions. Transcriptome analysis of the roots of a tolerant inbred HKI1105 using RNA sequencing revealed 21,364 differentially expressed genes (DEGs) under waterlogged stress condition. These 21,364 DEGs are known to regulate important pathways including energy-production, programmed cell death (PCD), aerenchyma formation, and ethylene responsiveness. High up-regulation of invertase (49-fold) and hexokinase (36-fold) in roots explained the ATP requirement in waterlogging condition. Also, high up-regulation of expansins (42-fold), plant aspartic protease A3 (19-fold), polygalacturonases (16-fold), respiratory burst oxidase homolog (12-fold), and hydrolases (11-fold) explained the PCD of root cortical cells followed by the formation of aerenchyma tissue during waterlogging stress. We hypothesized that the oxygen transfer in waterlogged roots is promoted by a cross-talk of fermentative, metabolic, and glycolytic pathways that generate ATPs for PCD and aerenchyma formation in root cortical cells. SNPs were mapped to the DEGs regulating aerenchyma formation (12), ethylene-responsive factors (11), and glycolysis (4) under stress. RNAseq derived SNPs can be used in selection approaches to breed tolerant hybrids. Overall, this investigation provided significant evidence of genes operating in the adaptive traits such as ethylene production and aerenchyma formation to cope-up the waterlogging stress.

Variable Level of Dominance of Candidate Genes Controlling Drought Functional Traits in Maize Hybrids.

Breeding maize for drought tolerance necessitates the knowledge on tolerant genotypes, molecular basis of drought tolerance mechanism, action, and expression pattern of genes. Studying the expression pattern and gene action of candidate genes during drought stress in the hybrids will help in choosing target genes for drought tolerance breeding. In the present investigation, a set of five hybrids and their seven parents with a variable level of tolerance to drought stress was selected to study the magnitude and the direction of 52 drought-responsive candidate genes distributed across various biological functions, viz., stomatal regulation, root development, detoxification, hormone signaling, photosynthesis, and sugar metabolism. The tolerant parents, HKI1105 and CML425, and their hybrid, ADWLH2, were physiologically active under drought stress, since vital parameters viz., chlorophyll, root length and relative water content, were on par with the respective well-watered control. All the genes were up-regulated in ADWLH2, many were down-regulated in HM8 and HM9, and most were down-regulated in PMH1 and PMH3 in the shoots and roots. The nature of the gene action was controlled by the parental combination rather than the parent per se. The differentially expressed genes in all five hybrids explained a mostly non-additive gene action over additivity, which was skewed toward any of the parental lines. Tissue-specific gene action was also noticed in many of the genes. The non-additive gene action is driven by genetic diversity, allele polymorphism, events during gene regulation, and small RNAs under the stress condition. Differential regulation and cross-talk of genes controlling various biological functions explained the basis of drought tolerance in subtropical maize hybrids. The nature of the gene action and the direction of the expression play crucial roles in designing introgression and hybrid breeding programmes to breed drought tolerant maize hybrids.

Identification, Characterization, and Functional Validation of Drought-responsive MicroRNAs in Subtropical Maize Inbreds.

MicroRNA-mediated gene regulation plays a crucial role in controlling drought tolerance. In the present investigation, 13 drought-associated miRNA families consisting of 65 members and regulating 42 unique target mRNAs were identified from drought-associated microarray expression data in maize and were subjected to structural and functional characterization. The largest number of members (14) was found in the zma-miR166 and zma-miR395 families, with several targets. However, zma-miR160, zma-miR390, zma-miR393, and zma-miR2275 each showed a single target. Twenty-three major drought-responsive cis-regulatory elements were found in the upstream regions of miRNAs. Many drought-related transcription factors, such as GAMYB, HD-Zip III, and NAC, were associated with the target mRNAs. Furthermore, two contrasting subtropical maize genotypes (tolerant: HKI-1532 and sensitive: V-372) were used to understand the miRNA-assisted regulation of target mRNA under drought stress. Approximately 35 and 31% of miRNAs were up-regulated in HKI-1532 and V-372, respectively. The up-regulation of target mRNAs was as high as 14.2% in HKI-1532 but was only 2.38% in V-372. The expression patterns of miRNA-target mRNA pairs were classified into four different types: Type I- up-regulation, Type II- down-regulation, Type III- neutral regulation, and Type IV- opposite regulation. HKI-1532 displayed 46 Type I, 13 Type II, and 23 Type III patterns, whereas V-372 had mostly Type IV interactions (151). A low level of negative regulations of miRNA associated with a higher level of mRNA activity in the tolerant genotype helped to maintain crucial biological functions such as ABA signaling, the auxin response pathway, the light-responsive pathway and endosperm expression under stress conditions, thereby leading to drought tolerance. Our study identified candidate miRNAs and mRNAs operating in important pathways under drought stress conditions, and these candidates will be useful in the development of drought-tolerant maize hybrids.

Genomic Selection for Drought Tolerance Using Genome-Wide SNPs in Maize.

Traditional breeding strategies for selecting superior genotypes depending on phenotypic traits have proven to be of limited success, as this direct selection is hindered by low heritability, genetic interactions such as epistasis, environmental-genotype interactions, and polygenic effects. With the advent of new genomic tools, breeders have paved a way for selecting superior breeds. Genomic selection (GS) has emerged as one of the most important approaches for predicting genotype performance. Here, we tested the breeding values of 240 maize subtropical lines phenotyped for drought at different environments using 29,619 cured SNPs. Prediction accuracies of seven genomic selection models (ridge regression, LASSO, elastic net, random forest, reproducing kernel Hilbert space, Bayes A and Bayes B) were tested for their agronomic traits. Though prediction accuracies of Bayes B, Bayes A and RKHS were comparable, Bayes B outperformed the other models by predicting highest Pearson correlation coefficient in all three environments. From Bayes B, a set of the top 1053 significant SNPs with higher marker effects was selected across all datasets to validate the genes and QTLs. Out of these 1053 SNPs, 77 SNPs associated with 10 drought-responsive transcription factors. These transcription factors were associated with different physiological and molecular functions (stomatal closure, root development, hormonal signaling and photosynthesis). Of several models, Bayes B has been shown to have the highest level of prediction accuracy for our data sets. Our experiments also highlighted several SNPs based on their performance and relative importance to drought tolerance. The result of our experiments is important for the selection of superior genotypes and candidate genes for breeding drought-tolerant maize hybrids.

Genomewide Expression and Functional Interactions of Genes under Drought Stress in Maize.

A genomewide transcriptome assay of two subtropical genotypes of maize was used to observe the expression of genes at seedling stage of drought stress. The number of genes expressed differentially was greater in HKI1532 (a drought tolerant genotype) than in PC3 (a drought sensitive genotype), indicating primary differences at the transcriptional level in stress tolerance. The global coexpression networks of the two genotypes differed significantly with respect to the number of modules and the coexpression pattern within the modules. A total of 174 drought-responsive genes were selected from HKI1532, and their coexpression network revealed key correlations between different adaptive pathways, each cluster of the network representing a specific biological function. Transcription factors related to ABA-dependent stomatal closure, signalling, and phosphoprotein cascades work in concert to compensate for reduced photosynthesis. Under stress, water balance was maintained by coexpression of the genes involved in osmotic adjustments and transporter proteins. Metabolism was maintained by the coexpression of genes involved in cell wall modification and protein and lipid metabolism. The interaction of genes involved in crucial biological functions during stress was identified and the results will be useful in targeting important gene interactions to understand drought tolerance in greater detail.

Comparative Analysis of CDPK Family in Maize, Arabidopsis, Rice, and Sorghum Revealed Potential Targets for Drought Tolerance Improvement.

Calcium dependent protein kinases (CDPKs) play significant role in regulation of plant growth and development in response to various stresses including drought. A set of 32 CDPK genes identified in maize were further used for searching of orthologs in the model plant Arabidopsis (72) and major food crops such as rice (78) and sorghum (91). We comprehensively studied the phylogenetic relationship, annotations, gene duplications, gene structure, divergence time, 3-D protein structures and tissue-specific drought induced expression of CDPK genes in all four species. Variation in intron frequency in the studied species was one of the reasons for the functional diversity of CDPK genes to various stress responses. Protein kinase and protein kinase C phosphorylation site domains were the most conserved motifs identified in all species. Four groups were identified from the sequence-based phylogenetic analysis, in which maize CDPKs were clustered in group III. Expression data showed that the CDPK genes were highly expressed in leaf of maize, rice, and sorghum whereas in Arabidopsis the maximum expression was observed in root. The expression assay showed 5, 6, 11, and 9 were the commonly and differentially expressed drought-related orthologous genes in maize, Arabidopsis, rice, and sorghum, respectively. 3-D protein structure were predicted for the nine genes (Arabidopsis: 2, maize: 2, rice: 3, and sorghum: 2) showing differential expression in at least three species. The predicted 3-D structures were further evaluated and validated by Ramachandran plot, ANOLEA, ProSA, and Verify-3D. The superimposed 3-D structure of drought-related orthologous proteins retained similar folding pattern owing to their conserved nature. Functional annotation revealed the involvement of CDPK genes in various pathways such as osmotic homeostasis, cell protection, and root growth. The interactions of CDPK genes in various pathways play crucial role in imparting drought tolerance through different ABA and MAPK signaling cascades. These selected candidate genes could be targeted in development of drought tolerant genotypes in maize, rice, and sorghum through appropriate breeding approaches. Our comparative experiments of CDPK genes could also be extended in the drought stress breeding programmes of the related species.

In Silico Characterization and Functional Validation of Cell Wall Modification Genes Imparting Waterlogging Tolerance in Maize.

Cell wall modification (CWM) promotes the formation of aerenchyma in roots under waterlogging conditions as an adaptive mechanism. Lysigenous aerenchyma formation in roots improves oxygen transfer in plants, which highlights the importance of CWM as a focal point in waterlogging stress tolerance. We investigated the structural and functional compositions of CWM genes and their expression patterns under waterlogging conditions in maize. Cell wall modification genes were identified for 3 known waterlogging-responsive cis-acting regulatory elements, namely, GC motif, anaerobic response elements, and G-box, and 2 unnamed elements. Structural motifs mapped in CWM genes were represented in genes regulating waterlogging stress-tolerant pathways, including fermentation, glycolysis, programmed cell death, and reactive oxygen species signaling. The highly aligned regions of characterized and uncharacterized CWM proteins revealed common structural domains amongst them. Membrane spanning regions present in the protein structures revealed transmembrane activity of CWM proteins in the plant cell wall. Cell wall modification proteins had interacted with ethylene-responsive pathway regulating genes (E3 ubiquitin ligases RNG finger and F-box) in a maize protein-protein interaction network. Cell wall modification genes had also coexpressed with energy metabolism, programmed cell death, and reactive oxygen species signaling, regulating genes in a single coexpression cluster. These configurations of CWM genes can be used to modify the protein expression in maize under waterlogging stress condition. Our study established the importance of CWM genes in waterlogging tolerance, and these genes can be used as candidates in introgression breeding and genome editing experiments to impart tolerance in maize hybrids.

Stability Performance of Inductively Coupled Plasma Mass Spectrometry-Phenotyped Kernel Minerals Concentration and Grain Yield in Maize in Different Agro-Climatic Zones.

Deficiency of iron and zinc causes micronutrient malnutrition or hidden hunger, which severely affects ~25% of global population. Genetic biofortification of maize has emerged as cost effective and sustainable approach in addressing malnourishment of iron and zinc deficiency. Therefore, understanding the genetic variation and stability of kernel micronutrients and grain yield of the maize inbreds is a prerequisite in breeding micronutrient-rich high yielding hybrids to alleviate micronutrient malnutrition. We report here, the genetic variability and stability of the kernel micronutrients concentration and grain yield in a set of 50 maize inbred panel selected from the national and the international centres that were raised at six different maize growing regions of India. Phenotyping of kernels using inductively coupled plasma mass spectrometry (ICP-MS) revealed considerable variability for kernel minerals concentration (iron: 18.88 to 47.65 mg kg(-1); zinc: 5.41 to 30.85 mg kg(-1); manganese: 3.30 to 17.73 mg kg(-1); copper: 0.53 to 5.48 mg kg(-1)) and grain yield (826.6 to 5413 kg ha(-1)). Significant positive correlation was observed between kernel iron and zinc within (r = 0.37 to r = 0.52, p < 0.05) and across locations (r = 0.44, p < 0.01). Variance components of the additive main effects and multiplicative interactions (AMMI) model showed significant genotype and genotype × environment interaction for kernel minerals concentration and grain yield. Most of the variation was contributed by genotype main effect for kernel iron (39.6%), manganese (41.34%) and copper (41.12%), and environment main effects for both kernel zinc (40.5%) and grain yield (37.0%). Genotype main effect plus genotype-by-environment interaction (GGE) biplot identified several mega environments for kernel minerals and grain yield. Comparison of stability parameters revealed AMMI stability value (ASV) as the better representative of the AMMI stability parameters. Dynamic stability parameter GGE distance (GGED) showed strong and positive correlation with both mean kernel concentrations and grain yield. Inbreds (CM-501, SKV-775, HUZM-185) identified from the present investigation will be useful in developing micronutrient-rich as well as stable maize hybrids without compromising grain yield.