Genome-Wide Analysis and Prediction of Resistance to Goss's Wilt in Maize

A genome-wide association study identifies genes associated with cuticular wax metabolism in maize

The plant cuticle is essential in plant defense against biotic and abiotic stresses. To systematically elucidate the genetic architecture of maize (Zea mays L.) cuticular wax metabolism, 2 cuticular wax-related traits, the chlorophyll extraction rate (CER) and water loss rate (WLR) of 389 maize inbred lines, were investigated and a genome-wide association study (GWAS) was performed using 1.25 million single nucleotide polymorphisms (SNPs). In total, 57 nonredundant quantitative trait loci (QTL) explaining 5.57% to 15.07% of the phenotypic variation for each QTL were identified. These QTLs contained 183 genes, among which 21 strong candidates were identified based on functional annotations and previous publications. Remarkably, 3 candidate genes that express differentially during cuticle development encode β-ketoacyl-CoA synthase (KCS). While ZmKCS19 was known to be involved in cuticle wax metabolism, ZmKCS12 and ZmKCS3 functions were not reported. The association between ZmKCS12 and WLR was confirmed by resequencing 106 inbred lines, and the variation of WLR was significant between different haplotypes of ZmKCS12. In this study, the loss-of-function mutant of ZmKCS12 exhibited wrinkled leaf morphology, altered wax crystal morphology, and decreased C32 wax monomer levels, causing an increased WLR and sensitivity to drought. These results confirm that ZmKCS12 plays a vital role in maize C32 wax monomer synthesis and is critical for drought tolerance. In sum, through GWAS of 2 cuticular wax-associated traits, this study reveals comprehensively the genetic architecture in maize cuticular wax metabolism and provides a valuable reference for the genetic improvement of stress tolerance in maize.

Identification of loci conferring resistance to 4 foliar diseases of maize

Foliar diseases of maize are among the most important diseases of maize worldwide. This study focused on 4 major foliar diseases of maize: Goss's wilt, gray leaf spot, northern corn leaf blight, and southern corn leaf blight. QTL mapping for resistance to Goss's wilt was conducted in 4 disease resistance introgression line populations with Oh7B as the common recurrent parent and Ki3, NC262, NC304, and NC344 as recurrent donor parents. Mapping results for Goss's wilt resistance were combined with previous studies for gray leaf spot, northern corn leaf blight, and southern corn leaf blight resistance in the same 4 populations. We conducted (1) individual linkage mapping analysis to identify QTL specific to each disease and population; (2) Mahalanobis distance analysis to identify putative multiple disease resistance regions for each population; and 3) joint linkage mapping to identify QTL across the 4 populations for each disease. We identified 3 lines that were resistant to all 4 diseases. We mapped 13 Goss's wilt QTLs in the individual populations and an additional 6 using joint linkage mapping. All Goss's wilt QTL had small effects, confirming that resistance to Goss's wilt is highly quantitative. We report several potentially important chromosomal bins associated with multiple disease resistance including 1.02, 1.03, 3.04, 4.06, 4.08, and 9.03. Together, these findings indicate that disease QTL distribution is not random and that there are locations in the genome that confer resistance to multiple diseases. Furthermore, resistance to bacterial and fungal diseases is not entirely distinct, and we identified lines resistant to both fungi and bacteria, as well as loci that confer resistance to both bacterial and fungal diseases.

Advancements and Prospects of Genome-Wide Association Studies (GWAS) in Maize

Genome-wide association studies (GWAS) have emerged as a powerful tool for unraveling intricate genotype-phenotype association across various species. Maize (Zea mays L.), renowned for its extensive genetic diversity and rapid linkage disequilibrium (LD), stands as an exemplary candidate for GWAS. In maize, GWAS has made significant advancements by pinpointing numerous genetic loci and potential genes associated with complex traits, including responses to both abiotic and biotic stress. These discoveries hold the promise of enhancing adaptability and yield through effective breeding strategies. Nevertheless, the impact of environmental stress on crop growth and yield is evident in various agronomic traits. Therefore, understanding the complex genetic basis of these traits becomes paramount. This review delves into current and future prospectives aimed at yield, quality, and environmental stress resilience in maize and also addresses the challenges encountered during genomic selection and molecular breeding, all facilitated by the utilization of GWAS. Furthermore, the integration of omics, including genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics has enriched our understanding of intricate traits in maize, thereby enhancing environmental stress tolerance and boosting maize production. Collectively, these insights not only advance our understanding of the genetic mechanism regulating complex traits but also propel the utilization of marker-assisted selection in maize molecular breeding programs, where GWAS plays a pivotal role. Therefore, GWAS provides robust support for delving into the genetic mechanism underlying complex traits in maize and enhancing breeding strategies.

Integrating genome-wide association study into genomic selection for the prediction of agronomic traits in rice (Oryza sativa L.)

Accurately identifying varieties with targeted agronomic traits was thought to contribute to genetic selection and accelerate rice breeding progress. Genomic selection (GS) is a promising technique that uses markers covering the whole genome to predict the genomic-estimated breeding values (GEBV), with the ability to select before phenotypes are measured. To choose the appropriate GS models for breeding work, we analyzed the predictability of nine agronomic traits measured from a population of 459 diverse rice varieties. By the comparison of eight representative GS models, we found that the prediction accuracies ranged from 0.407 to 0.896, with reproducing kernel Hilbert space (RKHS) having the highest predictive ability in most traits. Further results demonstrated the predictivity of GS is altered by several factors. Moreover, we assessed the method of integrating genome-wide association study (GWAS) into various GS models. The predictabilities of GS combined peak-associated markers generated from six different GWAS models were significantly different; a recommendation of Mixed Linear Model (MLM)-RKHS was given for the GWAS-GS-integrated prediction. Finally, based on the above result, we experimented with applying the P-values obtained from optimal GWAS models into ridge regression best linear unbiased prediction (rrBLUP), which benefited the low predictive traits in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01423-y.

Genetic and transcriptomic dissection of host defense to Goss's bacterial wilt and leaf blight of maize

Goss's wilt, caused by the Gram-positive actinobacterium Clavibacter nebraskensis, is an important bacterial disease of maize. The molecular and genetic mechanisms of resistance to the bacterium, or, in general, Gram-positive bacteria causing plant diseases, remain poorly understood. Here, we examined the genetic basis of Goss's wilt through differential gene expression, standard genome-wide association mapping (GWAS), extreme phenotype (XP) GWAS using highly resistant (R) and highly susceptible (S) lines, and quantitative trait locus (QTL) mapping using 3 bi-parental populations, identifying 11 disease association loci. Three loci were validated using near-isogenic lines or recombinant inbred lines. Our analysis indicates that Goss's wilt resistance is highly complex and major resistance genes are not commonly present. RNA sequencing of samples separately pooled from R and S lines with or without bacterial inoculation was performed, enabling identification of common and differential gene responses in R and S lines. Based on expression, in both R and S lines, the photosynthesis pathway was silenced upon infection, while stress-responsive pathways and phytohormone pathways, namely, abscisic acid, auxin, ethylene, jasmonate, and gibberellin, were markedly activated. In addition, 65 genes showed differential responses (up- or down-regulated) to infection in R and S lines. Combining genetic mapping and transcriptional data, individual candidate genes conferring Goss's wilt resistance were identified. Collectively, aspects of the genetic architecture of Goss's wilt resistance were revealed, providing foundational data for mechanistic studies.

Goss's Wilt Resistance in Corn Is Mediated via Salicylic Acid and Programmed Cell Death but Not Jasmonic Acid Pathways

A highly aggressive strain (CMN14-5-1) of Clavibacter nebraskensis bacteria, which causes Goss's wilt in corn, induced severe symptoms in a susceptible corn line (CO447), resulting in water-soaked lesions followed by necrosis within a few days. A tolerant line (CO450) inoculated with the same strain exhibited only mild symptoms such as chlorosis, freckling, and necrosis that did not progress after the first six days following infection. Both lesion length and disease severity were measured using the area under the disease progression curve (AUDPC), and significant differences were found between treatments. We analyzed the expression of key genes related to plant defense in both corn lines challenged with the CMN14-5-1 strain. Allene oxide synthase (ZmAOS), a gene responsible for the production of jasmonic acid (JA), was induced in the CO447 line in response to CMN14-5-1. Following inoculation with CMN14-5-1, the CO450 line demonstrated a higher expression of salicylic acid (SA)-related genes, ZmPAL and ZmPR-1, compared to the CO447 line. In the CO450 line, four genes related to programmed cell death (PCD) were upregulated: respiratory burst oxidase homolog protein D (ZmrbohD), polyphenol oxidase (ZmPPO1), ras-related protein 7 (ZmRab7), and peptidyl-prolyl cis-trans isomerase (ZmPPI). The differential gene expression in response to CMN14-5-1 between the two corn lines provided an indication that SA and PCD are involved in the regulation of corn defense responses against Goss's wilt disease, whereas JA may be contributing to disease susceptibility.

Genome-wide association and genomic prediction for resistance to southern corn rust in DH and testcross populations

Southern corn rust (SCR), caused by Puccinia polysora Underw, is a destructive disease that can severely reduce grain yield in maize (Zea mays L.). Owing to P. polysora being multi-racial, it is very important to explore more resistance genes and develop more efficient selection approaches in maize breeding programs. Here, four Doubled Haploid (DH) populations with 384 accessions originated from selected parents and their 903 testcross hybrids were used to perform genome-wide association (GWAS). Three GWAS processes included the additive model in the DH panel, additive and dominant models in the hybrid panel. As a result, five loci were detected on chromosomes 1, 7, 8, 8, and 10, with P-values ranging from 4.83×10^-7 to 2.46×10^-41. In all association analyses, a highly significant locus on chromosome 10 was detected, which was tight chained with the known SCR resistance gene RPPC and RPPK. Genomic prediction (GP), has been proven to be effective in plant breeding. In our study, several models were performed to explore predictive ability in hybrid populations for SCR resistance, including extended GBLUP with different genetic matrices, maker based prediction models, and mixed models with QTL as fixed factors. For GBLUP models, the prediction accuracies ranged from 0.56-0.60. Compared with traditional prediction only with additive effect, prediction ability was significantly improved by adding additive-by-additive effect (P-value< 0.05). For maker based models, the accuracy of BayesA and BayesB was 0.65, 8% higher than other models (i.e., RRBLUP, BRR, BL, BayesC). Finally, by adding QTL into the mixed linear prediction model, the accuracy can be further improved to 0.67, especially for the G_A model, the prediction performance can be increased by 11.67%. The prediction accuracy of the BayesB model can be further improved significantly by adding QTL information (P-value< 0.05). This study will provide important valuable information for understanding the genetic architecture and the application of GP for SCR in maize breeding.

Characterization and genetic dissection of maize ear leaf midrib acquired by 3D digital technology

The spatial morphological structure of plant leaves is an important index to evaluate crop ideotype. In this study, we characterized the three-dimensional (3D) data of the ear leaf midrib of maize at the grain-filling stage using the 3D digitization technology and obtained the phenotypic values of 15 traits covering four different dimensions of the ear leaf midrib, of which 13 phenotypic traits were firstly proposed for featuring plant leaf spatial structure. Cluster analysis results showed that the 13 traits could be divided into four groups, Group I, -II, -III and -IV. Group I contains HorizontalLength, OutwardGrowthMeasure, LeafAngle and DeviationTip; Group II contains DeviationAngle, MaxCurvature and CurvaturePos; Group III contains LeafLength and ProjectionArea; Group IV contains TipTop, VerticalHeight, UpwardGrowthMeasure, and CurvatureRatio. To investigate the genetic basis of the ear leaf midrib curve, 13 traits with high repeatability were subjected to genome-wide association study (GWAS) analysis. A total of 828 significantly related SNPs were identified and 1365 candidate genes were annotated. Among these, 29 candidate genes with the highest significant and multi-method validation were regarded as the key findings. In addition, pathway enrichment analysis was performed on the candidate genes of traits to explore the potential genetic mechanism of leaf midrib curve phenotype formation. These results not only contribute to further understanding of maize leaf spatial structure traits but also provide new genetic loci for maize leaf spatial structure to improve the plant type of maize varieties.

Dissecting the Genetic Structure of Maize Leaf Sheaths at Seedling Stage by Image-Based High-Throughput Phenotypic Acquisition and Characterization

The rapid development of high-throughput phenotypic detection techniques makes it possible to obtain a large number of crop phenotypic information quickly, efficiently, and accurately. Among them, image-based phenotypic acquisition method has been widely used in crop phenotypic identification and characteristic research due to its characteristics of automation, non-invasive, non-destructive and high throughput. In this study, we proposed a method to define and analyze the traits related to leaf sheaths including morphology-related, color-related and biomass-related traits at V6 stage. Next, we analyzed the phenotypic variation of leaf sheaths of 418 maize inbred lines based on 87 leaf sheath-related phenotypic traits. In order to further analyze the mechanism of leaf sheath phenotype formation, 25 key traits (2 biomass-related, 19 morphology-related and 4 color-related traits) with heritability greater than 0.3 were analyzed by genome-wide association studies (GWAS). And 1816 candidate genes of 17 whole plant leaf sheath traits and 1,297 candidate genes of 8 sixth leaf sheath traits were obtained, respectively. Among them, 46 genes with clear functional descriptions were annotated by single nucleotide polymorphism (SNPs) that both Top1 and multi-method validated. Functional enrichment analysis results showed that candidate genes of leaf sheath traits were enriched into multiple pathways related to cellular component assembly and organization, cell proliferation and epidermal cell differentiation, and response to hunger, nutrition and extracellular stimulation. The results presented here are helpful to further understand phenotypic traits of maize leaf sheath and provide a reference for revealing the genetic mechanism of maize leaf sheath phenotype formation.

Colonization and Movement of Green Fluorescent Protein-Labeled Clavibacter nebraskensis in Maize

Clavibacter nebraskensis causes Goss's bacterial wilt and leaf blight, a major disease of maize. Infected crop residue is the primary inoculum source and infection can occur via wounds or natural openings, such as stomata or hydathodes. The use of resistant hybrids is the primary control method for Goss's wilt. In this study, colonization and movement patterns of C. nebraskensis during infection were examined using green fluorescent protein (GFP)-labeled bacterial strains. We successfully introduced a plasmid to C. nebraskensis via electroporation, which resulted in GFP accumulation. Fluorescence microscopy revealed that in the absence of wounding, bacteria colonize leaf tissue via entry through the hydathodes when guttation droplets are present. Stomatal penetration was not observed under natural conditions. Bacteria initially colonize the xylem and subsequently the mesophyll, which creates the freckles that are characteristic of the disease. Bacteria infiltrated into the mesophyll did not cause disease symptoms, could not enter the vasculature, and did not spread from the initial inoculation point. Bacteria were observed exuding through stomata onto the leaf surface, resulting in the characteristic sheen of diseased leaves. Resistant maize lines exhibited decreased bacterial spread in the vasculature and the mesophyll. These tools to examine C. nebraskensis movement offer opportunities and new insights into the pathogenesis process and can form the basis for improved Goss's wilt management through host resistance.

Genome-wide association mapping in maize: status and prospects

Genome-wide association study (GWAS) provides a robust and potent tool to retrieve complex phenotypic traits back to their underlying genetics. Maize is an excellent crop for performing GWAS due to diverse genetic variability, rapid decay of linkage disequilibrium, availability of distinct sub-populations and abundant SNP information. The application of GWAS in maize has resulted in successful identification of thousands of genomic regions associated with many abiotic and biotic stresses. Many agronomic and quality traits of maize are severely affected by such stresses and, significantly affecting its growth and productivity. To improve productivity of maize crop in countries like India which contribute only 2% to the world's total production in 2019-2020, it is essential to understand genetic complexity of underlying traits. Various DNA markers and trait associations have been revealed using conventional linkage mapping methods. However, it has achieved limited success in improving polygenic complex traits due to lower resolution of trait mapping. The present review explores the prospects of GWAS in improving yield, quality and stress tolerance in maize besides, strengths and challenges of using GWAS for molecular breeding and genomic selection. The information gathered will facilitate elucidation of genetic mechanisms of complex traits and improve efficiency of marker-assisted selection in maize breeding.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02799-4.

Genetic Diversity and Population Structure of Maize Inbred Lines with Varying Levels of Resistance to Striga Hermonthica Using Agronomic Trait-Based and SNP Markers

Striga hermonthica is a serious biotic stress limiting maize production in sub-Saharan Africa. The limited information on the patterns of genetic diversity among maize inbred lines derived from source germplasm with mixed genetic backgrounds limits the development of inbred lines, hybrids, and synthetics with durable resistance to S. hermonthica. This study was conducted to assess the level of genetic diversity in a panel of 150 diverse maize inbred lines using agronomic and molecular data and also to infer the population structure among the inbred lines. Ten Striga-resistance-related traits were used for the phenotypic characterization, and 16,735 high-quality single-nucleotide polymorphisms (SNPs), identified by genotyping-by-sequencing (GBS), were used for molecular diversity. The phenotypic and molecular hierarchical cluster analyses grouped the inbred lines into five clusters, respectively. However, the grouping patterns between the phenotypic and molecular hierarchical cluster analyses were inconsistent due to non-overlapping information between the phenotypic and molecular data. The correlation between the phenotypic and molecular diversity matrices was very low (0.001), which is in agreement with the inconsistencies observed between the clusters formed by the phenotypic and molecular diversity analyses. The joint phenotypic and genotypic diversity matrices grouped the inbred lines into three groups based on their reaction patterns to S. hermonthica, and this was able to exploit a broad estimate of the actual diversity among the inbred lines. The joint analysis shows an invaluable insight for measuring genetic diversity in the evaluated materials. The result indicates that wide genetic variability exists among the inbred lines and that the joint diversity analysis can be utilized to reliably assign the inbred lines into heterotic groups and also to enhance the level of resistance to Striga in new maize varieties.

Identification of Loci That Confer Resistance to Bacterial and Fungal Diseases of Maize

Crops are hosts to numerous plant pathogenic microorganisms. Maize has several major disease issues; thus, breeding multiple disease resistant (MDR) varieties is critical. While the genetic basis of resistance to multiple fungal pathogens has been studied in maize, less is known about the relationship between fungal and bacterial resistance. In this study, we evaluated a disease resistance introgression line (DRIL) population for the foliar disease Goss’s bacterial wilt and blight (GW) and conducted quantitative trait locus (QTL) mapping. We identified a total of ten QTL across multiple environments. We then combined our GW data with data on four additional foliar diseases (northern corn leaf blight, southern corn leaf blight, gray leaf spot, and bacterial leaf streak) and conducted multivariate analysis to identify regions conferring resistance to multiple diseases. We identified 20 chromosomal bins with putative multiple disease effects. We examined the five chromosomal regions (bins 1.05, 3.04, 4.06, 8.03, and 9.02) with the strongest statistical support. By examining how each haplotype effected each disease, we identified several regions associated with increased resistance to multiple diseases and three regions associated with opposite effects for bacterial and fungal diseases. In summary, we identified several promising candidate regions for multiple disease resistance in maize and specific DRILs to expedite interrogation.