Overexpression of ZmEXPA5 reduces anthesis-silking interval and increases grain yield under drought and well-watered conditions in maize.

Drought is one of the major abiotic stresses affecting the maize production worldwide. As a cross-pollination crop, maize is sensitive to water stress at flowering stage. Drought at this stage leads to asynchronous development of male and female flower organ and increased interval between anthesis and silking, which finally causes failure of pollination and grain yield loss. In the present study, the expansin gene ZmEXPA5 was cloned and its function in drought tolerance was characterized. An indel variant in promoter of ZmEXPA5 is significantly associated with natural variation in drought-induced anthesis-silking interval. The drought susceptible haplotypes showed lower expression level of ZmEXPA5 than tolerant haplotypes and lost the cis-regulatory activity of ZmDOF29. Increasing ZmEXPA5 expression in transgenic maize decreases anthesis-silking interval and improves grain yield under both drought and well-watered environments. In addition, the expression pattern of ZmEXPA5 was analyzed. These findings provide insights into the genetic basis of drought tolerance and a promising gene for drought improvement in maize breeding. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01432-x.

cis-Regulatory variation affecting gene expression contributes to the improvement of maize kernel size.

cis-Regulatory variations contribute to trait evolution and adaptation during crop domestication and improvement. As the most important harvested organ in maize (Zea mays L.), kernel size has undergone intensive selection for size. However, the associations between maize kernel size and cis-regulatory variations remain unclear. We chose two independent association populations to dissect the genetic architecture of maize kernel size together with transcriptomic and genotypic data. The resulting phenotypes reflected a strong influence of population structure on kernel size. Compared with genome-wide association studies (GWASs), which accounted for population structure and relatedness, GWAS based on a naïve or simple linear model revealed additional associated single-nucleotide polymorphisms significantly involved in the conserved pathways controlling seed size in plants. Regulation analyses through expression quantitative trait locus mapping revealed that cis-regulatory variations likely control kernel size by fine-tuning the expression of proximal genes, among which ZmKL1 (GRMZM2G098305) was transgenically validated. We also proved that the pyramiding of the favorable cis-regulatory variations has contributed to the improvement of maize kernel size. Collectively, our results demonstrate that cis-regulatory variations, together with their regulatory genes, provide excellent targets for future maize improvement.

Genomic insights into historical improvement of heterotic groups during modern hybrid maize breeding.

Single-cross maize hybrids display superior heterosis and are produced from crossing two parental inbred lines belonging to genetically different heterotic groups. Here we assembled 1,604 historically utilized maize inbred lines belonging to various female heterotic groups (FHGs) and male heterotic groups (MHGs), and conducted phenotyping and genomic sequencing analyses. We found that the FHGs and MHGs have undergone both convergent and divergent changes for different sets of agronomic traits. Using genome-wide selection scans and association analyses, we identified a large number of candidate genes that contributed to the improvement of agronomic traits of the FHGs and MHGs. Moreover, we observed increased genetic differentiation between the FHGs and MHGs across the breeding eras, and we found a positive correlation between increasing heterozygosity levels in the differentiated genes and heterosis in hybrids. Furthermore, we validated the function of two selected genes and a differentiated gene. This study provides insights into the genomic basis of modern hybrid maize breeding.

Fine mapping qKRN5.04 provides a functional gene negatively regulating maize kernel row number.

KEY MESSAGE: Zm00001d016075 was identified by fine mapping qKRN5.04. The function of Zm00001d016075, negatively modulating maize (Zea Mays L.) kernel row number (KRN), was verified by CRISPR-Cas9. InDel308 located in the promoter of Zm00001d016075 has potential for use as a molecular marker to identify KRN in maize breeding. Kernel row number (KRN), controlled by multiple quantitative trait loci (QTLs), is one of the most important traits that relate to maize production and domestication. Here, fine mapping was conducted to study a major QTL, qKRN5.04, to a 65-kb genomic region using a progeny test strategy in an advanced backcross population, in which Nong531 (N531) and the inbred line of Silunuo (SLN) were employed as the recurrent and the donor parent, respectively. Within this region, there was only one gene (Zm00001d016075) based on the B73 reference genome. Furthermore, we performed regional association mapping using a panel of 236 diverse inbred lines and observed that all significant SNPs were located within Zm00001d016075. The expression of Zm00001d016075 was significantly higher in N531 and qKRN5.04N531 than qKRN5.04SLN, resulting from the different promoter activity of Zm00001d016075. Sequence analysis revealed that InDel308, located in the promoter of Zm00001d016075, was related to the KRN variation in different maize inbred lines. Using the CRISPR-Cas9 strategy, we determined Zm00001d016075 played a role in negatively regulating KRN and had a moderate effect on 10-kernel width, 100-kernel weight, kernels per ear, and grain yield per ear. These results provide critical insights on the genetic basis and quantitative variation for KRN.

Genetic Variation in ZmPAT7 Contributes to Tassel Branch Number in Maize.

Tassel branch number (TBN) is one of the important agronomic traits that contribute to the efficiency of seed production and has been selected strongly during the modern maize breeding process. However, the genetic mechanisms of TBN in maize are not entirely clear. In this study, we used a B73 × CML247 recombination inbred lines (RILs) population to detect quantitative trait loci (QTLs) for TBN. A total of four QTLs (qTBN2a, qTBN2b, qTBN4, and qTBN6) and six candidate genes were identified through expression analysis. Further, one of the candidates (GRMZM2G010011, ZmPAT7) encoding an S-acyltransferase was selected to validate its function by CRISPR-Cas9 technology, and its loss-of-function lines showed a significant increase in TBN. A key SNP(-101) variation in the promoter of ZmPAT7 was significantly associated with TBN. A total of 17 distant eQTLs associated with the expression of ZmPAT7 were identified in expression quantitative trait loci (eQTL) analysis, and ZmNAC3 may be a major factor involved in regulating ZmPAT7. These findings of the present study promote our understanding of the genetic basis of tassel architecture and provide new gene resources for maize breeding improvement.

Genome-wide association studies of leaf angle in maize.

Compact plant-type with small leaf angle has increased canopy light interception, which is conducive to the photosynthesis of the population and higher population yield at high density planting in maize. In this study, a panel of 285 diverse maize inbred lines genotyped with 56,000 SNPs was used to investigate the genetic basis of leaf angle across 3 consecutive years using a genome-wide association study (GWAS). The leaf angle showed broad phenotypic variation and high heritability across different years. Population structure analysis subdivided the panel into four subgroups that correspond to the four major empirical germplasm origins in China, i.e., Tangsipingtou, Reid, Lancaster and P. When tested with the optimal GWAS model, we found that the Q + K model was the best in reducing false positive. In total, 96 SNPs accounting for 5.54-10.44% of phenotypic variation were significantly (P < 0.0001) associated with leaf angle across three years. According to the linkage disequilibrium decay distance, 96 SNPs were binned into 43 QTLs for leaf angle. Seven major QTLs with R2 > 8% stably detected in at least 2 years, and BLUP values were clustered in four genomic regions (bins 2.01, 2.07, 5.06, and 10.04). Seven important candidate genes, Zm00001d001961, Zm00001d006348, Zm00001d006463, Zm00001d017618, Zm00001d024919, Zm00001d025018, and Zm00001d025033 were predicted for the seven stable major QTLs, respectively. The markers identified in this study can be used for molecular breeding for leaf angle, and the candidate genes would contribute to further understanding of the genetic basis of leaf angle. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01241-0.

Genome-wide association studies and whole-genome prediction reveal the genetic architecture of KRN in maize.

BACKGROUND: Kernel row number (KRN) is an important trait for the domestication and improvement of maize. Exploring the genetic basis of KRN has great research significance and can provide valuable information for molecular assisted selection. RESULTS: In this study, one single-locus method (MLM) and six multilocus methods (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB and ISIS EM-BLASSO) of genome-wide association studies (GWASs) were used to identify significant quantitative trait nucleotides (QTNs) for KRN in an association panel including 639 maize inbred lines that were genotyped by the MaizeSNP50 BeadChip. In three phenotyping environments and with best linear unbiased prediction (BLUP) values, the seven GWAS methods revealed different numbers of KRN-associated QTNs, ranging from 11 to 177. Based on these results, seven important regions for KRN located on chromosomes 1, 2, 3, 5, 9, and 10 were identified by at least three methods and in at least two environments. Moreover, 49 genes from the seven regions were expressed in different maize tissues. Among the 49 genes, ARF29 (Zm00001d026540, encoding auxin response factor 29) and CKO4 (Zm00001d043293, encoding cytokinin oxidase protein) were significantly related to KRN, based on expression analysis and candidate gene association mapping. Whole-genome prediction (WGP) of KRN was also performed, and we found that the KRN-associated tagSNPs achieved a high prediction accuracy. The best strategy was to integrate all of the KRN-associated tagSNPs identified by all GWAS models. CONCLUSIONS: These results aid in our understanding of the genetic architecture of KRN and provide useful information for genomic selection for KRN in maize breeding.

Identification of genomic insertion and flanking sequences of the transgenic drought-tolerant maize line "SbSNAC1-382" using the single-molecule real-time (SMRT) sequencing method.

Safety assessment of genetically modified (GM) crops is crucial at the product-development phase before GM crops are placed on the market. Determining characteristics of sequences flanking exogenous insertion sequences is essential for the safety assessment and marketing of transgenic crops. In this study, we used genome walking and whole-genome sequencing (WGS) to identify the flanking sequence characteristics of the SbSNAC1 transgenic drought-tolerant maize line "SbSNAC1-382", but both of the two methods failed. Then, we constructed a genomic fosmid library of the transgenic maize line, which contained 4.18×105 clones with an average insertion fragment of 35 kb, covering 5.85 times the maize genome. Subsequently, three positive clones were screened by pairs of specific primers, and one of the three positive clones was sequenced by using single-molecule real-time (SMRT) sequencing technology. More than 1.95 Gb sequence data (~105× coverage) for the sequenced clone were generated. The junction reads mapped to the boundaries of T-DNA, and the flanking sequences in the transgenic line were identified by comparing all sequencing reads with the maize reference genome and the sequence of the transgenic vector. Furthermore, the putative insertion loci and flanking sequences were confirmed by PCR amplification and Sanger sequencing. The results indicated that two copies of the exogenous T-DNA fragments were inserted at the same genomic site, and the exogenous T-DNA fragments were integrated at the position of Chromosome 5 from 177155650 to 177155696 in the transgenic line 382. In this study, we demonstrated the successful application of the SMRT technology for the characterization of genomic insertion and flanking sequences.

Transcriptome and GWAS analyses reveal candidate gene for seminal root length of maize seedlings under drought stress.

Water deficits are a major constraint on maize growth and yield, and deep roots are one of the major mechanisms of drought tolerance. In this study, four root and shoot traits were evaluated within an association panel consisting of 209 diverse maize accessions under well-watered (WW) and water-stressed (WS) conditions. A significant positive correlation was observed between seminal root length (SRL) under WS treatment and the drought tolerance index (DI) of maize seedlings. The transcriptome profiles of maize seminal roots were compared between four drought-tolerant lines and four drought-sensitive lines under both water conditions to identify genes associated with the drought stress response. After drought stress, 343 and 177 common differentially expressed genes (DEGs) were identified in the drought-tolerant group and drought-sensitive group, respectively. In parallel, a coexpression network underlying SRL was constructed on the basis of transcriptome data, and 10 hub genes involved in two significant associated modules were identified. Additionally, a genome-wide association study (GWAS) of the SRL revealed 62 loci for the two water treatments. By integrating the results of the GWAS, the common DEGs and the coexpression network analysis, 7 promising candidate genes were prioritized for further research. Together, our results provide a foundation for the enhanced understanding of seminal root changes in response to drought stress in maize.

The retromer protein ZmVPS29 regulates maize kernel morphology likely through an auxin-dependent process(es).

Kernel size and morphology are two important yield-determining traits in maize, but their molecular and genetic mechanisms are poorly characterized. Here, we identified a major QTL, qKM4.08, which explains approximately 24.20% of the kernel morphology variance in a recombinant population derived from two elite maize inbred lines, Huangzaosi (HZS, round kernel) and LV28 (slender kernel). Positional cloning and transgenic analysis revealed that qKM4.08 encodes ZmVPS29, a retromer complex component. Compared with the ZmVPS29 HZS allele, the ZmVPS29 LV28 allele showed higher expression in developing kernels. Overexpression of ZmVPS29 conferred a slender kernel morphology and increased the yield per plant in different maize genetic backgrounds. Sequence analysis revealed that ZmVPS29 has been under purifying selection during maize domestication. Association analyses identified two significant kernel morphology-associated polymorphic sites in the ZmVPS29 promoter region that were significantly enriched in modern maize breeding lines. Further study showed that ZmVPS29 increased auxin accumulation during early kernel development by enhancing auxin biosynthesis and transport and reducing auxin degradation and thereby improved kernel development. Our results suggest that ZmVPS29 regulates kernel morphology, most likely through an auxin-dependent process(es).

Characterization and fine mapping of qkrnw4, a major QTL controlling kernel row number in maize.

KEY MESSAGE: A major QTL controlling kernel row number, qkrnw4, was identified by combining linkage analysis and association mapping. Within qkrnw4, on the basis of its expression and bioinformatics analysis, Zm00001d052910 was supposed to be the candidate gene for kernel row number. Kernel row number (KRN) is an important yield-related trait that affects kernel number in maize. Understanding the genetic basis of KRN is important for increasing maize yields. In the present study, by the use of a near-isogenic line (NIL) that has a B73 background and that consistently displays a low KRN across environments, qkrnw4, a major quantitative trait locus (QTL) associated with KRN within a yield trait-related QTL hotspot in bin 4.08, was finely mapped to an ~ 33-kb interval. Regional association analysis of a nested association mapping population comprising 5000 recombinant inbred lines revealed Zm00001d052910, which encodes a protein with an unknown function, as the important candidate gene responsible for qkrnw4. Different expression levels of this candidate gene in immature ears were detected between the NIL and its recurrent parent. Moreover, the expression of several auxin-related genes was consistent with that of the candidate gene. Furthermore, the potential associations of this candidate gene with well-known inflorescence-related genes were discussed. The results of this study provide important information for the genetic elucidation of KRN variation in maize.

Genetic architecture of phenotypic means and plasticities of kernel size and weight in maize.

KEY MESSAGE: Genetic relationships between the phenotypic means and plasticities of kernel size and weight revealed the common genetic control of these traits in maize. Kernel size and weight are crucial components of grain yield in maize, and phenotypic plasticity in these traits facilitates adaptations to changing environments. Elucidating the genetic architecture of the mean phenotypic values and plasticities of kernel size and weight may be essential for breeding climate-robust maize varieties. Here, a maize nested association mapping (CN-NAM) population and association panel were grown in different environments. A joint linkage analysis and genome-wide association mapping were performed for five kernel size and weight phenotypic traits and two phenotypic plasticity measures. The mean phenotypes and plasticities were significantly correlated. The overall results of quantitative trait locus (QTL) and candidate gene analyses indicated moderate and high levels of common genetic control for the two traits. Furthermore, the mean phenotypes or plasticities of the hundred-kernel weight and volume were commonly regulated to a high degree. One pleiotropic locus on chromosome 10 simultaneously controlled the mean phenotypic values and plasticities of kernel size and weight. Therefore, the plasticity of kernel size and weight might be indirectly selected during maize breeding; however, selecting for high or low plasticity in combination with high or low mean phenotypic values of kernel size and weight traits may be difficult.

Genome-wide identification and comparative analysis of drought-related microRNAs in two maize inbred lines with contrasting drought tolerance by deep sequencing.

Drought has become one of the most serious abiotic stresses influencing crop production worldwide. Understanding the molecular regulatory networks underlying drought adaption and tolerance in crops is of great importance for future breeding. microRNAs (miRNAs), as important components of post-transcriptional regulation, play crucial roles in drought response and adaptation in plants. Here, we report a miRNome analysis of two maize inbred lines with contrasting levels of drought tolerance under soil drought in the field. Differential expression analysis showed 11 and 34 miRNAs were uniquely responded to drought in H082183 (drought tolerant) and Lv28 (drought sensitive), respectively, in leaves. In roots, 19 and 23 miRNAs uniquely responded to drought in H082183 and Lv28, respectively. Expression analysis of these drought-responsive miRNA-mRNA modules revealed miR164-MYB, miR164-NAC, miR159-MYB, miR156-SPL and miR160-ARF showed a negative regulatory relationship. Further analysis showed that the miR164-MYB and miR164-NAC modules in the tolerant line modulated the stress response in an ABA (abscisic acid)-dependent manner, while the miR156-SPL and miR160-ARF modules in the sensitive line participated in the inhibition of metabolism in drought-exposed leaves. Together, our results provide new insight into not only drought-tolerance-related miRNA regulation networks in maize but also key miRNAs for further characterization and improvement of maize drought tolerance.

Candidate loci for the kernel row number in maize revealed by a combination of transcriptome analysis and regional association mapping.

BACKGROUND: The kernel row number (KRN) of an ear is an important trait related to yield and domestication in maize. Exploring the underlying genetic mechanisms of KRN has great research significance and application potential. RESULTS: In the present study, N531 with a KRN of 18-22 and SLN with a KRN of 4-6 were used as the recurrent parent and the donor parent, respectively, to develop two introgression lines (ILs), IL_A and IL_B, both of which have common negative-effect alleles from SLN on chromosomes 1, 5 and 10 and significantly reduced inflorescence meristem (IM) diameter and KRN compared with those of N531. We used RNA-Seq to investigate the transcriptome profiles of 5-mm immature ears of N531, IL_A and IL_B. We identified a total of 2872 differentially expressed genes (DEGs) between N531 and IL_A, 2428 DEGs between N531 and IL_B and 1811 DEGs between IL_A and IL_B. A total of 1252 DEGs were detected as overlapping DEGs, and 89 DEGs were located on the common introgression fragments. Furthermore, three DEGs (Zm00001d013277, Zm00001d015310 and Zm00001d015377) containing three SNPs associated with KRN were identified using regional association mapping. CONCLUSIONS: These results will facilitate our understanding of ear development and provide important candidate genes for further study on KRN.

Genome-wide analysis of the pentatricopeptide repeat gene family in different maize genomes and its important role in kernel development.

BACKGROUND: The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants (450 PPR genes in Arabidopsis, 477 PPR genes in rice and 486 PPR genes in foxtail millet) and is important for plant development and growth. Most PPR genes are encoded by plastid and mitochondrial genomes, and the gene products regulate the expression of the related genes in higher plants. However, the functions remain largely unknown, and systematic analysis and comparison of the PPR gene family in different maize genomes have not been performed. RESULTS: In this study, systematic identification and comparison of PPR genes from two elite maize inbred lines, B73 and PH207, were performed. A total of 491 and 456 PPR genes were identified in the B73 and PH207 genomes, respectively. Basic bioinformatics analyses, including of the classification, gene structure, chromosomal location and conserved motifs, were conducted. Examination of PPR gene duplication showed that 12 and 15 segmental duplication gene pairs exist in the B73 and PH207 genomes, respectively, with eight duplication events being shared between the two genomes. Expression analysis suggested that 53 PPR genes exhibit qualitative variations in the different genetic backgrounds. Based on analysis of the correlation between PPR gene expression in kernels and kernel-related traits, four PPR genes are significantly negatively correlated with hundred kernel weight, 12 are significantly negatively correlated with kernel width, and eight are significantly correlated with kernel number. Eight of the 24 PPR genes are also located in metaQTL regions associated with yield and kernel-related traits in maize. Two important PPR genes (GRMZM2G353195 and GRMZM2G141202) might be regarded as important candidate genes associated with maize kernel-related traits. CONCLUSIONS: Our results provide a more comprehensive understanding of PPR genes in different maize inbred lines and identify important candidate genes related to kernel development for subsequent functional validation in maize.

Increased experimental conditions and marker densities identified more genetic loci associated with southern and northern leaf blight resistance in maize.

Southern leaf blight (SLB) and northern leaf blight (NLB) are the two major foliar diseases limiting maize production worldwide. Upon previous study with the nested association mapping (NAM) population, which consist of 5,000 recombinant inbred lines from 25 parents crossed with B73, we expanded the phenotyping environments from the United States (US) to China, and increased the marker densities from 1106 to 7386 SNPs for linkage mapping, and from 1.6 to 28.5 million markers for association mapping. We identified 49 SLB and 48 NLB resistance-related unique QTLs in linkage mapping, and multiple loci in association mapping with candidate genes involved in known plant disease-resistance pathways. Furthermore, an independent natural population with 282 diversified inbred lines were sequenced for four candidate genes selected based on their biological functions. Three of them demonstrated significant associations with disease resistance. These findings provided valuable resources for further implementations to develop varieties with superior resistance for NLB and SLB.

Characterization and fine mapping of qkc7.03: a major locus for kernel cracking in maize.

KEY MESSAGE: A major locus conferring kernel cracking in maize was characterized and fine mapped to an interval of 416.27 kb. Meanwhile, combining the results of transcriptomic analysis, the candidate gene was inferred. Seed development requires a proper structural and physiological balance between the maternal tissues and the internal structures of the seeds. In maize, kernel cracking is a disorder in this balance that seriously limits quality and yield and is characterized by a cracked pericarp at the kernel top and endosperm everting. This study elucidated the genetic basis and characterization of kernel cracking. Primarily, a near isogenic line (NIL) with a B73 background exhibited steady kernel cracking across environments. Therefore, deprived mapping populations were developed from this NIL and its recurrent parent B73. A major locus on chromosome 7, qkc7.03, was identified to be associated with the cracking performance. According to a progeny test of recombination events, qkc7.03 was fine mapped to a physical interval of 416.27 kb. In addition, obvious differences were observed in embryo development and starch granule arrangement within the endosperm between the NIL and its recurrent parent upon the occurrence of kernel cracking. Moreover, compared to its recurrent parent, the transcriptome of the NIL showed a significantly down-regulated expression of genes related to zeins, carbohydrate synthesis and MADS-domain transcription factors. The transcriptomic analysis revealed ten annotated genes within the target region of qkc7.03, and only GRMZM5G899476 was differently expressed between the NIL and its recurrent parent, indicating that this gene might be a candidate gene for kernel cracking. The results of this study facilitate the understanding of the potential mechanism underlying kernel cracking in maize.

Genome-wide identification of gene expression in contrasting maize inbred lines under field drought conditions reveals the significance of transcription factors in drought tolerance.

Drought is a major threat to maize growth and production. Understanding the molecular regulation network of drought tolerance in maize is of great importance. In this study, two maize inbred lines with contrasting drought tolerance were tested in the field under natural soil drought and well-watered conditions. In addition, the transcriptomes of their leaves was analyzed by RNA-Seq. In total, 555 and 2,558 genes were detected to specifically respond to drought in the tolerant and the sensitive line, respectively, with a more positive regulation tendency in the tolerant genotype. Furthermore, 4,700, 4,748, 4,403 and 4,288 genes showed differential expression between the two lines under moderate drought, severe drought and their well-watered controls, respectively. Transcription factors were enriched in both genotypic differentially expressed genes and specifically responsive genes of the tolerant line. It was speculated that the genotype-specific response of 20 transcription factors in the tolerance line and the sustained genotypically differential expression of 22 transcription factors might enhance tolerance to drought in maize. Our results provide new insight into maize drought tolerance-related regulation systems and provide gene resources for subsequent studies and drought tolerance improvement.

Candidate Loci for Yield-Related Traits in Maize Revealed by a Combination of MetaQTL Analysis and Regional Association Mapping.

Maize grain yield and related traits are complex and are controlled by a large number of genes of small effect or quantitative trait loci (QTL). Over the years, a large number of yield-related QTLs have been identified in maize and deposited in public databases. However, integrating and re-analyzing these data and mining candidate loci for yield-related traits has become a major issue in maize. In this study, we collected information on QTLs conferring maize yield-related traits from 33 published studies. Then, 999 of these QTLs were iteratively projected and subjected to meta-analysis to obtain metaQTLs (MQTLs). A total of 76 MQTLs were found across the maize genome. Based on a comparative genomics strategy, several maize orthologs of rice yield-related genes were identified in these MQTL regions. Furthermore, three potential candidate genes (Gene ID: GRMZM2G359974, GRMZM2G301884, and GRMZM2G083894) associated with kernel size and weight within three MQTL regions were identified using regional association mapping, based on the results of the meta-analysis. This strategy, combining MQTL analysis and regional association mapping, is helpful for functional marker development and rapid identification of candidate genes or loci.

Exploring Identity-By-Descent Segments and Putative Functions Using Different Foundation Parents in Maize.

Maize foundation parents (FPs) play no-alternative roles in hybrid breeding because they were widely used in the development of new lines and hybrids. The combination of different identity-by-descent (IBD) segments and genes could account for the formation patterns of different FPs, and knowledge of these IBD regions would provide an extensive foundation for the development of new candidate FP lines in future maize breeding. In this paper, a panel of 304 elite lines derived from FPs, i.e., B73, 207, Mo17, and Huangzaosi (HZS), was collected and analyzed using 43,252 single nucleotide polymorphism (SNP) markers. Most IBD segments specific to particular FP groups were identified, including 116 IBD segments in B73, 105 in Mo17, 111 in 207, and 190 in HZS. In these regions, 423 quantitative trait nucleotides (QTNs) associated with 15 agronomic traits and 804 candidate genes were identified. Some known adaptation-related genes, e.g., dwarf8 and vgt1 in HZS, zcn8 and epc in Mo17, and ZmCCT in 207, were validated as being tightly linked to particular IBD segments. In addition, numerous new candidate genes were also identified. For example, GRMZM2G154278 in HZS, which belongs to the cell cycle control family, was closely linked to a QTN of the ear height/plant height (EH/PH) trait; GRMZM2G051943 in 207, which encodes an endochitinase precursor (EP) chitinase, was closely linked to a QTN for kernel density; and GRMZM2G170586 in Mo17 was closely linked to a QTN for ear diameter. Complex correlations among these genes were also found. Many IBD segments and genes were included in the formation of FP lines, and complex regulatory networks exist among them. These results provide new insights on the genetic basis of complex traits and provide new candidate IBD regions or genes for the improvement of special traits in maize production.