Genome-Wide Association Study Dissects the Genetic Architecture of Maize Husk Tightness

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.

Utilizing Two Populations Derived from Tropical Maize for Genome-Wide Association Analysis of Banded Leaf and Sheath Blight Resistance

Banded leaf and sheath blight (BLSB) in maize is a soil-borne fungal disease caused by Rhizoctonia solani Kühn, resulting in significant yield losses. Investigating the genes responsible for regulating resistance to BLSB is crucial for yield enhancement. In this study, a multiparent maize population was developed, comprising two recombinant inbred line (RIL) populations totaling 442 F8RILs. The populations were generated by crossing two tropical inbred lines, CML444 and NK40-1, known for their BLSB resistance, as female parents, with the high-yielding but BLSB-susceptible inbred line Ye107 serving as the common male parent. Subsequently, we utilized 562,212 high-quality single nucleotide polymorphisms (SNPs) generated through genotyping-by-sequencing (GBS) for a comprehensive genome-wide association study (GWAS) aimed at identifying genes responsible for BLSB resistance. The objectives of this study were to (1) identify SNPs associated with BLSB resistance through genome-wide association analyses, (2) explore candidate genes regulating BLSB resistance in maize, and (3) investigate pathways involved in BLSB resistance and discover key candidate genes through Gene Ontology (GO) analysis. The GWAS analysis revealed nineteen SNPs significantly associated with BLSB that were consistently identified across four environments in the GWAS, with phenotypic variation explained (PVE) ranging from 2.48% to 11.71%. Screening a 40 kb region upstream and downstream of the significant SNPs revealed several potential candidate genes. By integrating information from maize GDB and the NCBI, we identified five novel candidate genes, namely, Zm00001d009723, Zm00001d009975, Zm00001d009566, Zm00001d009567, located on chromosome 8, and Zm00001d026376, on chromosome 10, related to BLSB resistance. These candidate genes exhibit association with various aspects, including maize cell membrane proteins and cell immune proteins, as well as connections to cell metabolism, transport, transcriptional regulation, and structural proteins. These proteins and biochemical processes play crucial roles in maize defense against BLSB. When Rhizoctonia solani invades maize plants, it induces the expression of genes encoding specific proteins and regulates corresponding metabolic pathways to thwart the invasion of this fungus. The present study significantly contributes to our understanding of the genetic basis of BLSB resistance in maize, offering valuable insights into novel candidate genes that could be instrumental in future breeding efforts to develop maize varieties with enhanced BLSB resistance.

Reveal the kernel dehydration mechanisms in maize based on proteomic and metabolomic analysis

BACKGROUND

Kernel dehydration is an important factor for the mechanized harvest in maize. Kernel moisture content (KMC) and kernel dehydration rate (KDR) are important indicators for kernel dehydration. Although quantitative trait loci and genes related to KMC have been identified, where most of them only focus on the KMC at harvest, these are still far from sufficient to explain all genetic variations, and the relevant regulatory mechanisms are still unclear. In this study, we tried to reveal the key proteins and metabolites related to kernel dehydration in proteome and metabolome levels. Moreover, we preliminarily explored the relevant metabolic pathways that affect kernel dehydration combined proteome and metabolome. These results could accelerate the development of further mechanized maize technologies.

RESULTS

In this study, three maize inbred lines (KB182, KB207, and KB020) with different KMC and KDR were subjected to proteomic analysis 35, 42, and 49 days after pollination (DAP). In total, 8,358 proteins were quantified, and 2,779 of them were differentially expressed proteins in different inbred lines or at different stages. By comparative analysis, K-means cluster, and weighted gene co-expression network analysis based on the proteome data, some important proteins were identified, which are involved in carbohydrate metabolism, stress and defense response, lipid metabolism, and seed development. Through metabolomics analysis of KB182 and KB020 kernels at 42 DAP, 18 significantly different metabolites, including glucose, fructose, proline, and glycerol, were identified.

CONCLUSIONS

In sum, we inferred that kernel dehydration could be regulated through carbohydrate metabolism, antioxidant systems, and late embryogenesis abundant protein and heat shock protein expression, all of which were considered as important regulatory factors during kernel dehydration process. These results shed light on kernel dehydration and provide new insights into developing cultivars with low moisture content.

Genomic selection to improve husk tightness based on genomic molecular markers in maize

Introduction

The husk tightness (HTI) in maize plays a crucial role in regulating the water content of ears during the maturity stage, thereby influencing the quality of mechanical grain harvesting in China. Genomic selection (GS), which employs molecular markers, offers a promising approach for identifying and selecting inbred lines with the desired HTI trait in maize breeding. However, the effectiveness of GS is contingent upon various factors, including the genetic architecture of breeding populations, sequencing platforms, and statistical models.

Methods

An association panel of maize inbred lines was grown across three sites over two years, divided into four subgroups. GS analysis for HTI prediction was performed using marker data from three sequencing platforms and six marker densities with six statistical methods.

Results

The findings indicate that a loosely attached husk can aid in the dissipation of water from kernels in temperate maize germplasms across most environments but not nessarily for tropical-origin maize. Considering the balance between GS prediction accuracy and breeding cost, the optimal prediction strategy is the rrBLUP model, the 50K sequencing platform, a 30% proportion of the test population, and a marker density of r2=0.1. Additionally, selecting a specific SS subgroup for sampling the testing set significantly enhances the predictive capacity for husk tightness.

Discussion

The determination of the optimal GS prediction strategy for HTI provides an economically feasible reference for the practice of molecular breeding. It also serves as a reference method for GS breeding of other agronomic traits.

The RHW1-ZCN4 regulatory pathway confers natural variation of husk leaf width in maize

Maize husk leaf - the outer leafy layers covering the ear - modulates kernel yield and quality. Despite its importance, however, the genetic controls underlying husk leaf development remain elusive. Our previous genome-wide association study identified a single nucleotide polymorphism located in the gene RHW1 (Regulator of Husk leaf Width) that is significantly associated with husk leaf-width diversity in maize. Here, we further demonstrate that a polymorphic 18-bp InDel (insertion/deletion) variant in the 3' untranslated region of RHW1 alters its protein abundance and accounts for husk leaf width variation. RHW1 encodes a putative MYB-like transcriptional repressor. Disruption of RHW1 altered cell proliferation and resulted in a narrower husk leaf, whereas RHW1 overexpression yielded a wider husk leaf. RHW1 positively regulated the expression of ZCN4, a well-known TFL1-like protein involved in maize ear development. Dysfunction of ZCN4 reduced husk leaf width even in the context of RHW1 overexpression. The InDel variant in RHW1 is subject to selection and is associated with maize husk leaf adaption from tropical to temperate regions. Overall, our results identify that RHW1-ZCN4 regulates a pathway conferring husk leaf width variation at a very early stage of husk leaf development in maize.

Investigating the genetic basis of maize ear characteristics: a comprehensive genome-wide study utilizing high-throughput phenotypic measurement method and system

The morphology of maize ears plays a critical role in the breeding of new varieties and increasing yield. However, the study of traditional ear-related traits alone can no longer meet the requirements of breeding. In this study, 20 ear-related traits, including size, shape, number, and color, were obtained in 407 maize inbred lines at two sites using a high-throughput phenotypic measurement method and system. Significant correlations were found among these traits, particularly the novel trait ear shape (ES), which was correlated with traditional traits: kernel number per row and kernel number per ear. Pairwise comparison tests revealed that the inbred lines of tropical-subtropical were significantly different from other subpopulations in row numbers per ear, kernel numbers per ear, and ear color. A genome-wide association study identified 275, 434, and 362 Single nucleotide polymorphisms (SNPs) for Beijing, Sanya, and best linear unbiased prediction scenarios, respectively, explaining 3.78% to 24.17% of the phenotypic variance. Furthermore, 58 candidate genes with detailed functional descriptions common to more than two scenarios were discovered, with 40 genes being associated with color traits on chromosome 1. After analysis of haplotypes, gene expression, and annotated information, several candidate genes with high reliability were identified, including Zm00001d051328 for ear perimeter and width, zma-MIR159f for ear shape, Zm00001d053080 for kernel width and row number per ear, and Zm00001d048373 for the blue color channel of maize kernels in the red-green-blue color model. This study emphasizes the importance of researching novel phenotypic traits in maize by utilizing high-throughput phenotypic measurements. The identified genetic loci enrich the existing genetic studies related to maize ears.

Genome-Wide Association Study Identified Novel SNPs Associated with Chlorophyll Content in Maize

Chlorophyll is an essential component that captures light energy to drive photosynthesis. Chlorophyll content can affect photosynthetic activity and thus yield. Therefore, mining candidate genes of chlorophyll content will help increase maize production. Here, we performed a genome-wide association study (GWAS) on chlorophyll content and its dynamic changes in 378 maize inbred lines with extensive natural variation. Our phenotypic assessment showed that chlorophyll content and its dynamic changes were natural variations with a moderate genetic level of 0.66/0.67. A total of 19 single-nucleotide polymorphisms (SNPs) were found associated with 76 candidate genes, of which one SNP, 2376873-7-G, co-localized in chlorophyll content and area under the chlorophyll content curve (AUCCC). Zm00001d026568 and Zm00001d026569 were highly associated with SNP 2376873-7-G and encoded pentatricopeptide repeat-containing protein and chloroplastic palmitoyl-acyl carrier protein thioesterase, respectively. As expected, higher expression levels of these two genes are associated with higher chlorophyll contents. These results provide a certain experimental basis for discovering the candidate genes of chlorophyll content and finally provide new insights for cultivating high-yield and excellent maize suitable for planting environment.

Genome-wide association analysis of sucrose and alanine contents in edamame beans

The sucrose and Alanine (Ala) content in edamame beans significantly impacts the sweetness flavor of edamame-derived products as an important attribute to consumers' acceptance. Unlike grain-type soybeans, edamame beans are harvested as fresh beans at the R6 to R7 growth stages when beans are filled 80-90% of the pod capacity. The genetic basis of sucrose and Ala contents in fresh edamame beans may differ from those in dry seeds. To date, there is no report on the genetic basis of sucrose and Ala contents in the edamame beans. In this study, a genome-wide association study was conducted to identify single nucleotide polymorphisms (SNPs) related to sucrose and Ala levels in edamame beans using an association mapping panel of 189 edamame accessions genotyped with a SoySNP50K BeadChip. A total of 43 and 25 SNPs was associated with sucrose content and Ala content in the edamame beans, respectively. Four genes (Glyma.10g270800, Glyma.08g137500, Glyma.10g268500, and Glyma.18g193600) with known effects on the process of sucrose biosynthesis and 37 novel sucrose-related genes were characterized. Three genes (Gm17g070500, Glyma.14g201100 and Glyma.18g269600) with likely relevant effects in regulating Ala content and 22 novel Ala-related genes were identified. In addition, by summarizing the phenotypic data of edamame beans from three locations in two years, three PI accessions (PI 532469, PI 243551, and PI 407748) were selected as the high sucrose and high Ala parental lines for the perspective breeding of sweet edamame varieties. Thus, the beneficial alleles, candidate genes, and selected PI accessions identified in this study will be fundamental to develop edamame varieties with improved consumers' acceptance, and eventually promote edamame production as a specialty crop in the United States.

Mapping Genetic Variation in Arabidopsis in Response to Plant Growth-Promoting Bacterium Azoarcus olearius DQS-4T

Plant growth-promoting bacteria (PGPB) can enhance plant health by facilitating nutrient uptake, nitrogen fixation, protection from pathogens, stress tolerance and/or boosting plant productivity. The genetic determinants that drive the plant-bacteria association remain understudied. To identify genetic loci highly correlated with traits responsive to PGPB, we performed a genome-wide association study (GWAS) using an Arabidopsis thaliana population treated with Azoarcus olearius DQS-4^T. Phenotypically, the 305 Arabidopsis accessions tested responded differently to bacterial treatment by improving, inhibiting, or not affecting root system or shoot traits. GWA mapping analysis identified several predicted loci associated with primary root length or root fresh weight. Two statistical analyses were performed to narrow down potential gene candidates followed by haplotype block analysis, resulting in the identification of 11 loci associated with the responsiveness of Arabidopsis root fresh weight to bacterial inoculation. Our results showed considerable variation in the ability of plants to respond to inoculation by A. olearius DQS-4^T while revealing considerable complexity regarding statistically associated loci with the growth traits measured. This investigation is a promising starting point for sustainable breeding strategies for future cropping practices that may employ beneficial microbes and/or modifications of the root microbiome.

Meta-analysis and co-expression analysis revealed stable QTL and candidate genes conferring resistances to Fusarium and Gibberella ear rots while reducing mycotoxin contamination in maize

Fusarium (FER) and Gibberella ear rots (GER) are the two most devastating diseases of maize (Zea mays L.) which reduce yield and affect grain quality worldwide, especially by contamination with mycotoxins. Genetic improvement of host resistance to effectively tackle FER and GER diseases requires the identification of stable quantitative trait loci (QTL) to facilitate the application of genomics-assisted breeding for improving selection efficiency in breeding programs. We applied improved meta-analysis algorithms to re-analyze 224 QTL identified in 15 studies based on dense genome-wide single nucleotide polymorphisms (SNP) in order to identify meta-QTL (MQTL) and colocalized genomic loci for fumonisin (FUM) and deoxynivalenol (DON) accumulation, silk (SR) and kernel (KR) resistances of both FER and GER, kernel dry-down rate (KDD) and husk coverage (HC). A high-resolution genetic consensus map with 36,243 loci was constructed and enabled the projection of 164 of the 224 collected QTL. Candidate genes (CG) mining was performed within the most refined MQTL, and identified CG were cross-validated using publicly available transcriptomic data of maize under Fusarium graminearum infection. The meta-analysis revealed 40 MQTL, of which 29 were associated each with 2-5 FER- and/or GER-related traits. Twenty-eight of the 40 MQTL were common to both FER and GER resistances and 19 MQTL were common to silk and kernel resistances. Fourteen most refined MQTL on chromosomes 1, 2, 3, 4, 7 and 9 harbored a total of 2,272 CG. Cross-validation identified 59 of these CG as responsive to FER and/or GER diseases. MQTL ZmMQTL2.2, ZmMQTL9.2 and ZmMQTL9.4 harbored promising resistance genes, of which GRMZM2G011151 and GRMZM2G093092 were specific to the resistant line for both diseases and encoded "terpene synthase21 (tps21)" and "flavonoid O-methyltransferase2 (fomt2)", respectively. Our findings revealed stable refined MQTL harboring promising candidate genes for use in breeding programs for improving FER and GER resistances with reduced mycotoxin accumulation. These candidate genes can be transferred into elite cultivars by integrating refined MQTL into genomics-assisted backcross breeding strategies.

Genome-wide association study reveals the genetic basis of brace root angle and diameter in maize

Brace roots are the main organ to support the above-ground part of maize plant. It involves in plant growth and development by water absorption and lodging resistance. The bracing root angle (BRA) and diameter (BRD) are important components of brace root traits. Illuminating the genetic basis of BRA and BRD will contribute the improvement for mechanized harvest and increasing production. A GWAS of BRA and BRD was conducted using an associated panel composed of 508 inbred lines of maize. The broad-sense heritability of BRA and BRD was estimated to be respectively 71% ± 0.19 and 52% ± 0.14. The phenotypic variation of BRA and BRD in the non-stiff stalk subgroup (NSS) and the stiff stalk subgroup (SS) subgroups are significantly higher than that in the tropical/subtropical subgroup (TST) subgroups. In addition, BRA and BRD are significantly positive with plant height (PH), ear length (EL), and kernel number per row (KNPR). GWAS revealed 27 candidate genes within the threshold of p < 1.84 × 10⁻⁶ by both MLM and BLINK models. Among them, three genes, GRMZM2G174736, GRMZM2G445169 and GRMZM2G479243 were involved in cell wall function, and GRMZM2G038073 encoded the NAC transcription factor family proteins. These results provide theoretical support for clarifying the genetic basis of brace roots traits.

A Genome-Wide Association Study Dissects the Genetic Architecture of the Metaxylem Vessel Number in Maize Brace Roots

Metaxylem vessels in maize brace roots are key tissue, and their number (MVN) affects plant water and inorganic salt transportation and lodging resistance. Dissecting the genetic basis of MVN in maize brace roots can help guide the genetic improvement of maize drought resistance and lodging resistance during late developmental stages. In this study, we used 508 inbred lines with tropical, subtropical, and temperate backgrounds to analyze the genetic architecture of MVN in maize brace roots. The phenotypic variation in MVN in brace roots was evaluated in three environments, which revealed broad natural variation and relative low levels of heritability (h ² = 0.42). Stiff-stalk lines with a temperate background tended to have higher MVNs than plants in other genetic backgrounds. MVN was significantly positively correlated with plant height, tassel maximum axis length, ear length, and kernel number per row, which indicates that MVN may affect plant morphological development and yield. In addition, MVN was extremely significantly negatively correlated with brace root radius, but significantly positively correlated with brace root angle (BRA), diameter, and number, thus suggesting that the morphological function of some brace root traits may be essentially determined by MVN. Association analysis of MVN in brace roots combined 1,253,814 single nucleotide polymorphisms (SNPs) using FarmCPU revealed a total of nine SNPs significantly associated with MVN at P < 7.96 × 10^-7. Five candidate genes for MVN that may participate in secondary wall formation (GRMZM2G168365, GRMZM2G470499, and GRMZM2G028982) and regulate flowering time (GRMZM2G381691 and GRMZM2G449165). These results provide useful information for understanding the genetic basis of MVN in brace root development. Further functional studies of identified candidate genes should help elucidate the molecular pathways that regulate MVN in maize brace roots.

Genetic analysis of three maize husk traits by QTL mapping in a maize-teosinte population

BACKGROUND

The maize husk consists of numerous leafy layers and plays vital roles in protecting the ear from pathogen infection and dehydration. Teosinte, the wild ancestor of maize, has about three layers of small husk outer covering the ear. Although several quantitative trait loci (QTL) underlying husk morphology variation have been reported, the genetic basis of husk traits between teosinte and maize remains unclear.

RESULTS

A linkage population including 191 BC₂F₈ inbred lines generated from the maize line Mo17 and the teosinte line X26-4 was used to identify QTL associated with three husk traits: i.e., husk length (HL), husk width (HW) and the number of husk layers (HN). The best linear unbiased predictor (BLUP) depicted wide phenotypic variation and high heritability of all three traits. The HL exhibited greater correlation with HW than HN. A total of 4 QTLs were identified including 1, 1, 2, which are associated with HL, HW and HN, respectively. The proportion of phenotypic variation explained by these QTLs was 9.6, 8.9 and 8.1% for HL, HN and HW, respectively.

CONCLUSIONS

The QTLs identified in this study will pave a path to explore candidate genes regulating husk growth and development, and benefit the molecular breeding program based on molecular marker-assisted selection to cultivate maize varieties with an ideal husk morphology.

Genome-Wide Association Analysis Identifies Resistance Loci for Bacterial Leaf Streak Resistance in Rice (Oryza sativa L.)

Bacterial leaf streak (BLS) caused by Xanthomonas oryzae pv. oryzicola (Xoc) is one of the most devastating diseases in rice production areas, especially in humid tropical and subtropical zones throughout Asia and worldwide. A genome-wide association study (GWAS) analysis conducted on a collection of 236 diverse rice accessions, mainly indica varieties, identified 12 quantitative trait loci (QTLs) on chromosomes 1, 2, 3, 4, 5, 8, 9 and 11, conferring resistance to five representative isolates of Thai Xoc. Of these, five QTLs conferred resistance to more than one Xoc isolates. Two QTLs, qBLS5.1 and qBLS2.3, were considered promising QTLs for broad-spectrum resistance to BLS. The xa5 gene was proposed as a potential candidate gene for qBLS5.1 and three genes, encoding pectinesterase inhibitor (OsPEI), eukaryotic zinc-binding protein (OsRAR1), and NDP epimerase function, were proposed as candidate genes for qBLS2.3. Results from this study provide an insight into the potential QTLs and candidate genes for BLS resistance in rice. The recessive xa5 gene is suggested as a potential candidate for strong influence on broad-spectrum resistance and as a focal target in rice breeding programs for BLS resistance.