Key message: Intensive public sector breeding efforts and public-private partnerships have led to the increase in genetic gains, and deployment of elite climate-resilient maize cultivars for the stress-prone environments in the tropics. Maize (Zea mays L.) plays a critical role in ensuring food and nutritional security, and livelihoods of millions of resource-constrained smallholders. However, maize yields in the tropical rainfed environments are now increasingly vulnerable to various climate-induced stresses, especially drought, heat, waterlogging, salinity, cold, diseases, and insect pests, which often come in combinations to severely impact maize crops. The International Maize and Wheat Improvement Center (CIMMYT), in partnership with several public and private sector institutions, has been intensively engaged over the last four decades in breeding elite tropical maize germplasm with tolerance to key abiotic and biotic stresses, using an extensive managed stress screening network and on-farm testing system. This has led to the successful development and deployment of an array of elite stress-tolerant maize cultivars across sub-Saharan Africa, Asia, and Latin America. Further increasing genetic gains in the tropical maize breeding programs demands judicious integration of doubled haploidy, high-throughput and precise phenotyping, genomics-assisted breeding, breeding data management, and more effective decision support tools. Multi-institutional efforts, especially public–private alliances, are key to ensure that the improved maize varieties effectively reach the climate-vulnerable farming communities in the tropics, including accelerated replacement of old/obsolete varieties.

Genomic selection is being used increasingly in plant breeding to accelerate genetic gain per unit time, one of the most important genomic selection applications in maize breeding is to predict and select the best un-phenotyped lines in bi-parental populations based on genomic estimated breeding values. In the present study, 22 bi-parental tropical maize populations genotyped with low density SNPs were used to evaluate the genomic prediction accuracy (rMG) of the six trait-environment combinations under various levels of training population size (TPS) and marker density (MD), and assess the effect of trait heritability (h2), TPS and MD on rMG estimation. Our results showed that: (1) moderate rMG values obtained for different trait-environment combinations, when 50% of the total genotypes was used as training population and approximate 200 SNPs were used for prediction; (2) rMG increased as the increase of h2, TPS and MD, both correlation and variance analyses showed that h2 is the most important factor and MD is the least important factor on rMG estimation for most of the trait-environment combinations; (3) the trend observed in rMG differed with the trend observed in rMG/h, and h is the square root of heritability of the predicted trait, it indicated that both rMG and rMG/h values should be presented in GS study to show the accuracy of genomic selection and the relative accuracy of genomic selection compared with phenotypic selection, respectively. This study provides useful information to maize breeders to design genomic selection workflow in their breeding programs.

Drought stress (DS) is a major constraint to maize yield production. Heat stress (HS) alone and in combination with DS are likely to become the increasing constraints. Association mapping and genomic prediction (GP) analyses were conducted in a collection of 300 tropical and subtropical maize inbred lines to reveal the genetic architecture of grain yield and flowering time under well-watered (WW), DS, HS, and combined DS and HS conditions. Out of the 381,165 genotyping-by-sequencing SNPs, 1549 SNPs were significantly associated with all the 12 trait-environment combinations, the average PVE (phenotypic variation explained) by these SNPs was 4.33%, and 541 of them had a PVE value greater than 5%. These significant associations were clustered into 446 genomic regions with a window size of 20 Mb per region, and 673 candidate genes containing the significantly associated SNPs were identified. In addition, 33 hotspots were identified for 12 trait-environment combinations and most were located on chromosomes 1 and 8. Compared with single SNP-based association mapping, the haplotype-based associated mapping detected fewer number of significant associations and candidate genes with higher PVE values. All the 688 candidate genes were enriched into 15 gene ontology terms, and 46 candidate genes showed significant differential expression under the WW and DS conditions. Association mapping results identified few overlapped significant markers and candidate genes for the same traits evaluated under different managements, indicating the genetic divergence between the individual stress tolerance and the combined drought and HS tolerance. The GP accuracies obtained from the marker-trait associated SNPs were relatively higher than those obtained from the genome-wide SNPs for most of the target traits. The genetic architecture information of the grain yield and flowering time revealed in this study, and the genomic regions identified for the different trait-environment combinations are useful in accelerating the efforts on rapid development of the stress-tolerant maize germplasm through marker-assisted selection and/or genomic selection.

Key message Molecular characterization information on genetic diversity, population structure and genetic relationships provided by this research will help maize breeders to better understand how to utilize the current CML collection. Abstract CIMMYT maize inbred lines (CMLs) have been widely used all over the world and have contributed greatly to both tropical and temperate maize improvement. Genetic diversity and population structure of the current CML collection and of six temperate inbred lines were assessed and relationships among all lines were determined with genotyping-by-sequencing SNPs. Results indicated that: (1) wider genetic distance and low kinship coefficients among most pairs of lines reflected the uniqueness of most lines in the current CML collection; (2) the population structure and genetic divergence between the Temperate subgroup and Tropical subgroups were clear; three major environmental adaptation groups (Lowland Tropical, Subtropical/Mid-altitude and Highland Tropical subgroups) were clearly present in the current CML collection; (3) the genetic diversity of the three Tropical subgroups was similar and greater than that of the Temperate subgroup; the average genetic distance between the Temperate and Tropical subgroups was greater than among Tropical subgroups; and (4) heterotic patterns in each environmental adaptation group estimated using GBS SNPs were only partially consistent with patterns estimated based on combining ability tests and pedigree information. Combining current heterotic information based on combining ability tests and the genetic relationships inferred from molecular marker analyses may be the best strategy to define heterotic groups for future tropical maize improvement. Information resulting from this research will help breeders to better understand how to utilize all the CMLs to select parental lines, replace testers, assign heterotic groups and create a core set of breeding germplasm.

Genotyping-by-Sequencing (GBS) is a low-cost, high-throughput genotyping method that relies on restriction enzymes to reduce genome complexity. GBS is being widely used for various genetic and breeding applications. In the present study, 2240 individuals from eight maize populations, including two association populations (AM), backcross first generation (BC1), BC1F2, F2, double haploid (DH), intermated B73 × Mo17 (IBM), and a recombinant inbred line (RIL) population, were genotyped using GBS. A total of 955,120 of raw data for SNPs was obtained for each individual, with an average genotyping error of 0.70%. The rate of missing genotypic data for these SNPs was related to the level of multiplex sequencing: ~ 25% missing data for 96-plex and ~ 55% for 384-plex. Imputation can greatly reduce the rate of missing genotypes to 12.65% and 3.72% for AM populations and bi-parental populations, respectively, although it increases total genotyping error. For analysis of genetic diversity and linkage mapping, unimputed data with a low rate of genotyping error is beneficial, whereas, for association mapping, imputed data would result in higher marker density and would improve map resolution. Because imputation does not influence the prediction accuracy, both unimputed and imputed data can be used for genomic prediction. In summary, GBS is a versatile and efficient SNP discovery approach for homozygous materials and can be effectively applied for various purposes in maize genetics and breeding.

To increase genetic gain for tolerance to drought, we aimed to identify environmentally stable QTL in per se and testcross combination under well-watered (WW) and drought stressed (DS) conditions and evaluate the possible deployment of QTL using marker assisted and/or genomic selection (QTL/GS-MAS). A total of 169 doubled haploid lines derived from the cross between CML495 and LPSC7F64 and 190 testcrosses (tester CML494) were evaluated in a total of 11 treatment-by-population combinations under WW and DS conditions. In response to DS, grain yield (GY) and plant height (PHT) were reduced while time to anthesis and the anthesis silking interval (ASI) increased for both lines and hybrids. Forty-eight QTL were detected for a total of nine traits. The allele derived from CML495 generally increased trait values for anthesis, ASI, PHT, the normalized difference vegetative index (NDVI) and the green leaf area duration (GLAD; a composite trait of NDVI, PHT and senescence) while it reduced trait values for leaf rolling and senescence. The LOD scores for all detected QTL ranged from 2.0 to 7.2 explaining 4.4 to 19.4% of the observed phenotypic variance with R² ranging from 0 (GY, DS, lines) to 37.3% (PHT, WW, lines). Prediction accuracy of the model used for genomic selection was generally higher than phenotypic variance explained by the sum of QTL for individual traits indicative of the polygenic control of traits evaluated here. We therefore propose to use QTL-MAS in forward breeding to enrich the allelic frequency for a few desired traits with strong additive QTL in early selection cycles while GS-MAS could be used in more mature breeding programs to additionally capture alleles with smaller additive effects.