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 ten 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.

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 ofinvertase(49-fold) andhexokinase(36-fold) in roots explained the ATP requirement in waterlogging condition. Also, high up-regulation ofexpansins(42-fold),plant aspartic protease A3(19-fold),polygalacturonases(16-fold),respiratory burst oxidase homolog(12-fold), andhydrolases(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.

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.

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.

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.

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.