Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 30 October 2025

Sec. Plant Genetics, Epigenetics and Chromosome Biology

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1671460

Genome-wide association studies on leaf midrib architecture in maize

Dongdong DangDongdong Dang1Shuwen JiShuwen Ji1Yubo LiuYubo Liu2Chunming BaiChunming Bai3Yanshu ZhuYanshu Zhu1Xuecai Zhang,Xuecai Zhang1,2Ao Zhang*Ao Zhang1*Yanye Ruan*Yanye Ruan1*
  • 1Shenyang City Key Laboratory of Maize Genomic Selection Breeding, College of Bioscience and Biotechnology, Shenyang Agricultural University, Shenyang, Liaoning, China
  • 2International Maize and Wheat Improvement Center (CIMMYT), El Batan, Texcoco, Mexico
  • 3Seed Industry Innovation Research Institute, Liaoning Academy of Agricultural Sciences, Shenyang, Liaoning, China

Introduction: Leaf midrib architecture (LMA) in maize plays essential roles in supporting leaf structure and facilitating photosynthetic performance. Despite its importance for plant architecture and yield potential, the genetic basis underlying key midrib traits remains poorly understood.

Methods: To dissect the genetic architecture of LMA, we performed a genome-wide association study (GWAS) using a diverse panel of 508 maize inbred lines representing tropical, subtropical, and temperate germplasm. Three traits—midrib width (MW), total midrib thickness (TMT), and midrib basal thickness (MBT)—were measured across two environments in Liaoning Province, China. Phenotypic data were analyzed using mixed linear models with population structure and kinship correction. Significant loci were identified under a Bonferroni-adjusted threshold, and candidate genes were functionally annotated and examined for expression patterns using publicly available transcriptomic data.

Results: All LMA traits exhibited continuous variation and moderate heritability (H^2 = 0.33–0.61), with significant positive correlations among them. GWAS identified six SNPs significantly associated with LMA on chromosomes 5, 7, and 8, corresponding to 97 genes in adjacent genomic intervals. Among these, 27 annotated genes were enriched in functions related to transcriptional regulation, hormone signaling, cytoskeleton organization, and cell development. Key candidate genes included GRMZM2G074124 (YABBY-domain factor), GRMZM2G130953 and GRMZM2G332390 (auxin-related), GRMZM2G079185 (LOB-domain protein), and GRMZM2G407517 (Actin7). Correlation analyses further revealed that LMA traits are significantly associated with yield-related parameters, including cob diameter and grain weight.

Discussion: Our findings demonstrate that LMA in maize is a complex quantitative trait governed by multiple genes with diverse biological functions. Several candidate genes involved in auxin response, cytoskeletal dynamics, and lateral organ development play crucial roles in midrib formation. The observed associations between midrib traits and yield components suggest that LMA can serve as a valuable target for improving canopy structure and productivity in maize breeding. These results provide novel insights and candidate loci for molecular dissection and genetic improvement of leaf midrib traits.

Introduction

Leaf is a photosynthetic organ of plants which provides organic nutrients for whole plant (Sobrado and Turner, 1986; Abadie et al., 2018). The leaf photosynthetic cells and the other parts of plant are connected by vascular bundle system, on which the export of leaf photosynthetic products and input of nutrients and water needed by photosynthesis depend (Bagniewska-Zadworna et al., 2012; Dumlao et al., 2012). Leaf vascular system or venation also control the angle of leaf facing sunlight, which affects photosynthetic efficiency. A midrib is a midvein or primary vein locating at the longitudinal central axis of leaf blade, playing a pivotal role in the transport of substances and in the control on leaf upright angle and canopy architecture (Rolland-Lagan and Prusinkiewicz, 2005). However, the molecular mechanism of venation formation in maize remains unclear.

The development of venation starts at the leaf primordium, which is also the first step of leaf development. DL (drooping leaf) gene is essential for the process of venation development (Eshed et al., 2004). Loss-of-function mutations in the dl locus in rice cause lack of midrib, resulting in the drooping leaf phenotype (Ishikawa et al., 2009). DL gene is expressed in the central region of leaf primordia where the midrib develops subsequently and regulates cell proliferation (Yamaguchi et al., 2004). For Arabidopsis and rice, the DL gene encodes a zinc finger protein belonging to YABBY family of transcription factor, a plant-specific transcription factor family, which play an essential role in plant developmental processes by regulating the expression of many downstream genes, either directly or indirectly (Ohmori et al., 2008). The function of leaf midrib depends on its architecture such as thickness, width and concavity, but, the molecular mechanism that regulates leaf midrib architecture (LMA) is yet to be revealed, which may involve multiple quantitative genes.

Genome-wide association analysis (GWAS) is a powerful approach to identify the quantitative genes, via detecting the statistical association of a specific trait with genetic variants across multiple genomes (Visscher et al., 2017; Luo et al., 2024). Maize has enormous genetic diversity and rapid linkage disequilibrium decay, which is quite suitable for performing GWAS. Several quantitative traits in maize have been dissected by GWAS (Sahito et al., 2024), including low-phosphorus stress tolerance (Xu et al., 2018), field grain drying rate (Dai et al., 2017), husk number and weight (Zhou et al., 2016), aflatoxin accumulation resistance (Tang et al., 2015), water-stressed resistance (Xue et al., 2013), ear rot resistance caused by Fusarium verticillioides (de Jong et al., 2018).

In this study, 508 maize inbred lines with 558,630 high-quality single nucleotide polymorphism (SNP) marker sites were used to perform GWAS on three leaf midrib traits: width, total thickness and base thickness from two locations. Six significant SNPs have been located and 27 candidate genes, encoding many different families of proteins, are selected. These results show that the architecture of leaf midrib in maize is a quantitative trait regulated by multiple genes with different functions.

Materials and methods

Plant materials and field experiments

An association population of 508 maize inbred lines consist of germplasms collected from tropical, subtropical and temperate areas (Yang et al., 2011, 2013). There were 252 tropical/subtropical lines and 256 temperate lines (Supplementary Table S1). The association panel can also be classified into four subgroups (Yang et al., 2011): stiff stalk (SS), non-stiff stalk (NSS), tropical/subtropical (TST) and admixed group (MIX), respectively, with 27, 70, 196 and 215 inbred lines.

The population was grown under two environments in China: Fushun (N41°51′, E123°54′) and Tieling (N42°18′, E123°51′) of Liaoning Province in 2017. The experiments used a randomized complete block design with three blocks, row lengths of 2 m, row spacing of 0.6 m, and plant spacing of 0.3 m. The experiment was conducted in three replicates. Agronomic practices was identical in both environments.

Phenotyping of midrib architecture

Three traits of leaf midrib architecture (LMA): midrib width (MW), total midrib thickness (TMT) and midrib base thickness (MBT) were measured at the one third from the blade base of ear-leaf (Figure 1), at the same time during the maturity stage, using Electronic Vernier Caliper. More than five well-pollinated plants were selected from each line, respectively, in two environments.

Figure 1
Curved geometric shape with a metallic gold finish, featuring three-dimensional measurements labeled as 1, 2, and 3. The dimensions indicate various angles and lengths along the structure.

Figure 1. Architectural traits of leaf midrib in maize. 1 midrib width (MW); 2 midrib base thickness (MBT); 3 total midrib thickness (TMT).

Phenotypic data statistical analyses

Correlation analysis, mapping, and analysis of variance for midrib architecture traits were carried out with TASSEL employing a mixed linear model (MLM). To avoid the occurrence of false-positive results, the MLM incorporated the influence of population structure (Q matrix) and affinities (K matrix), turning into K + Q model of MLM.

The MLM model algorithm is yijk=μ+Geni+Envj+Geni×Envj+Repk+ϵijkWhere yijk is the trait of interest; μ  is the overall mean; Geni, Envj, and are the effects of the i-th genotype, j-th year, and i-th genotype by j-th year interaction, respectively; is the effect of k-th replication; ϵijk is the residual effect of the i-th genotype, j-th year, k-th replication.

Genotype is considered as the fixed effect, whereas all other terms are declared as the random effects. Locations with heritability below 0.05 were excluded from the across environment analysis. Best linear unbiased prediction (BLUP) was done by “PROC MIXED” in SAS software for estimation of leaf midrib trait values of all testers (Kump et al., 2011).

Phenotypic data were subjected to descriptive statistical analysis, analysis of variance, and correlation analysis using SPSS 19.0 (IBM Corp., Armonk, NY, USA). Broad-sense heritability was calculated using the formula proposed by Knapp et al. (1985): h2=σg2/(σg2+σge2nEnvs+σ2nEnvs×nreps), Where σg2, σ2 and are the genotypic variance, error variance, and genotype-by-environment interaction variance, respectively; and nreps and nEnvs are the numbers of replications and environments, respectively.

Genome-wide association analysis

The BLUPs for each trait per inbred line across environments were estimated using PROC MIXED. All 508 lines were genotyped using the MaizeSNP50 BeadChip. The genome-wide association analyses based on these models were conducted with the software in the FarmCPU procedure in R-4.2.1 (Liu et al., 2016). We conducted genome-wide association mapping with 558,630 SNPs (MAF > 0.05) on phenotypic data of each environment. A Bonferroni-corrected threshold probability based on individual tests was calculated to correct for multiple comparisons, using P< 1/N, where N is the number of individual trait-SNP combinations tested, the final threshold is set at P<1.04 × 10-5.

The independent contribution rates of significant SNP loci to midrib traits were determined based on the Bonferroni-corrected threshold and were defined as the phenotypic contribution rates. These calculations were performed using the anova function in R. The algorithm is as follows:

Y=αi=1nXi+βP+ϵ

Where Y represents the three variables of leaf midrib traits; X represents the SNP genotype variable corresponding to each trait; P is the subgroup variable; i is the number of individuals in the population; α is the SNP effect; β is the subgroup effect; and ϵ is the random error term.

Gene function prediction

We used the maize line B73 reference genome (B73 RefGen_v2, https://www.maizegdb.org/) to identify candidate genes that were either included or close to the significantly associated SNPs. Annotated functions and relevant information for the candidate genes were obtained from the Maize Genetics and Genomics Database and the U.S. National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/). Finally, possible candidate genes were selected based on the functional annotations of all genes in the SNP locus and their expression in each tissue of B73.

Heat map of candidate genes

The raw RNA-Seq data, originating from Genevieve M. Hoopes et al. and downloaded from the NCBI database, were aligned to the maize B73 reference genome (RefGen_v2), and the corresponding raw transcriptomic reads are available from NCBI under the following BioProject IDs: PRJNA171684, PRJEB10574, PRJNA226757, PRJNA244661, PRJNA323555, and PRJNA369690 (Hoopes et al., 2019)., and the RNA-Seq reading aligned with the maize B73 reference genome (RefGen_v2). The Cufflinks software keeps only unique map reads for counting FPKMS. The value used in Figure 2 is the normalized FPKM count in the pericarp with respect to other tissues.

Figure 2
Heatmap illustrating gene expression across different plant tissues, labeled on the vertical axis, with specific genes listed on the horizontal axis. The color gradient represents expression levels from low (blue) to high (red). A dendrogram on the left clusters the tissues based on similarity.

Figure 2. Heat map of tissue-specific expression patterns of candidate genes identified by GWAS. The value used in the figure is the PRKM counts in leaf tissue and other tissues, as shown at the bottom of each column. Columns and rows are sorted according to similarity (top and left are hierarchical cluster analysis). Red, white and blue indicate the degree of expression in different tissues.

Results

Trait variation and heritability

Under the two environments, the MW range of the association mapping population is 2.99-11.34mm, the TMT range is 0.99-7.28mm, and the MBT range is 0.23-4.42mm, and the ranges of trait variations are wide (Table 1). Significant variance components for genotype (G) and genotype×environment (G×E) interactions were observed for all three traits, however, G×E interactions represented relatively a small proportion of the total variance. According to the formula proposed by Knapp et al. (1985). The broad-sense heritability’s (H2) estimated for MW, TMT, and MBT were moderate, which were 0.59, 0.61, and 0.33, respectively. The H2 of midrib traits estimated show that leaf midrib architecture traits are controlled by genetic effects, thus these traits are suitable for the association analysis.

Table 1
www.frontiersin.org

Table 1. Variance composition and broad-sense heritability of leaf midrib architecture traits in the association population of maize.

The best unbiased linear predictive value (BLUP) of three midrib traits shows a normal distribution, and the correlation between them are significant positive correlation (p ≤ 0.01) (Figure 3). TMT is highly positive correlated with MBT (r = 0.85, p ≤ 0.01), Pearson correlation coefficient reaching 0.87, the thicker the midrib thickness is, the thicker the midrib base thickness is. The significant correlation is present between leaf midrib architecture traits, indicating that the growth and development of LMA traits are highly coordinated.

Figure 3
Scatterplot matrix showing relationships between variables MW, TMT, and MBT. Each diagonal contains histograms with red density lines. Off-diagonal scatterplots include red trend lines. Correlation coefficients (0.38, 0.28, 0.85) and significance stars are shown.

Figure 3. Frequency distributions and correlations of three leaf midrib architecture traits. Plots on diagonal line show phenotypic distribution of each trait as indicated; values above diagonal line are Pearson’s correlation coefficients between traits; plots below diagonal line are scatter plots of compared traits. MW, midrib width; TMT, total midrib thickness; MBT, midrib base thickness. * Represents the significance level calculated at P ≤0.05; ** represents the significance level calculated at P ≤ 0.01.

The four subgroups, SS, NSS, TST and MIXED, of the maize association panel have different origins (Yang et al., 2011). The SS and NSS subpopulations are of temperate origins, TST subpopulations is of tropical and subtropical origins, and MIXED subpopulations contain other inbred lines that are not classified in the afore mentioned three subpopulations (Yang et al., 2011). We compared the variation of leaf midrib architecture traits among different subgroups through boxplots (Figure 4). The mean value of MW in the TST subgroup was larger compared with SS, NSS and MIXED subgroups, indicating that for the tropical/subtropical origins inbred lines of maize, MW tends to be wider (Figure 4A). For TMT and MBT, no significant differences in the mean value were observed between the subgroups, indicating that the population origins has no significant effect on leaf midrib traits of thickness (Figures 4B, C). However, the difference of leaf midrib architecture traits between inbred lines in each subgroup were observed (Figure 4), which endows the possibility of genetic improvements on LMA.

Figure 4
Violin plots showing data distributions for three variables across subgroups: (a) VW, (b) VT, and (c) VD. Subgroups are labeled as MIXED, NSS, SS, and TST, each represented by unique colors. Box plots within violins indicate median and quartile ranges. Significant differences are marked with asterisks in plot (a).

Figure 4. Boxplots of the distribution of leaf midrib architecture traits in four subpopulations. Analysis of variance (ANOVA) was applied to examine the difference of traits among subpopulations. Different numbers indicate statistically significant difference at P ≤ 0.05. No. of inbred lines included in each subpopulation are 215, 70, 27 and 196 for MIXED (admixed), NSS (non-stiff stalk), SS (stiff stalk) and TST (tropical/subtropical), respectively. (a) MW (midrib width); (b) TMT (total midrib thickness); (c) MBT (midrib base thickness).

Leaf midrib architecture traits related to other agronomic traits

In this study, correlation analysis was conducted on the width, total thickness and base thickness of leaf midrib architecture with 14 agronomic traits, including 7 plant traits, plant height (PH), ear (position) height (EH), ear leaf width (ELW), ear leaf length (ELL), tassel maximum axis length (TMAL), tassel branch number (TBN), leaf number above ear (LNAE); and 4 ear and kernel traits, ear length (EL), ear diameter (ED), cob diameter (CD), cob grain weight (GW); and 3 fertility stage traits, days to anthesis (DTA), days to silking (DTS) and days to heading (DTH). As shown in Figure 5, the three leaf midrib architecture traits are positively and significantly correlated with two their own ear leaf blade traits of ELW and ELL, indicating that the ear leaf blade traits and midrib architecture traits develop synchronously and interplay each other. The midrib traits are also significantly correlated with GW and other grain traits. The MW is significantly correlated with DTS, DTA and DTH.

Figure 5
Heatmap displaying correlations between three variables (MBT, TMT, MW) and multiple factors (PH, EH, ELW, etc.). Values range from -0.09 to 0.32, color-coded from blue (positive) to orange (negative). Asterisks indicate significance levels, with darker shades representing stronger correlations.

Figure 5. Correlation coefficient between leaf midrib architecture traits and other 14 agronomic traits based on BLUP value. MW, midrib width; TMT, total midrib thickness; MBT, midrib base thickness. PH, plant height; EH, ear (position) height; ELW, ear leaf width; ELL, ear leaf length; TMAL, tassel maximum axis length; TBN, tassel branch number; LNAE, leaf number above ear; EL, ear length; ED, ear diameter; CD, cob diameter; GW, cob grain weight; DTA, days to anthesis; DTS, days to silking; DTH, days to heading. *Significant at P ≤ 0.05; **significant at P ≤ 0.01.

GWAS for leaf midrib architecture traits

Genome-wide association analyses were performed using the FarmCPU procedure implemented in R. To minimize the influence of environmental variation, phenotypic BLUP values from two environments were used for GWAS (Figure 6). The results showed that four significant SNPs associated with MBT were identified on chromosomes 5 and 8, one significant SNP associated with TMT was detected on chromosome 7, and one significant SNP associated with MW was identified on chromosome 5.

Figure 6
Manhattan and Q-Q plots showing genetic association data. Plots a, c, and e are Manhattan plots with points colored by chromosome. Significance threshold indicated by red dashed lines. Plots b, d, and f are Q-Q plots displaying observed versus expected -log₁₀(p) values, with data points and a reference line.

Figure 6. Genome-wide association studies of leaf midrib architecture traits. GWAS-derived Manhattan plots and Quantile-quantile plot showing significant P-values associated with leaf midrib architecture traits using MLM. Each dot represents an SNP. The horizontal dashed red line represents the Bonferroni-corrected significant threshold 1.04 × 10-5. (a, b) MW; (c, d) TMT; (e, f) MBT.

Genes co-localized with significant GWAS SNPs

Searching the range of 50 kb for each significant SNP upstream and downstream, a total of 97 genes (Supplementary Table S2) were found in the afore-mentioned 6 sites, based on the B73 reference genome and the LD attenuation distance of each significant site. Of which, 31 genes have functional annotations, and the functions of the remaining 66 genes are unknown. In detail, for the MBT trait, 30 genes were located, of which 12 are annotated; for the TMT, 17 genes were determined, among them 3 are annotated; for the MW, up to 50 genes were detected, 12 of which are annotated (Supplementary Table S3).

To determine whether these candidate genes are specifically expressed in leaf venation traits, the published RNA-seq data set were used to compile expression patterns in different organs/tissues, and a heat map of the tissue-specific expression patterns of candidate genes was produced (Figure 2). Twenty-seven annotated genes were identified as candidate genes associated with the leaf midrib architecture traits. According to NCBI and previous research results, these genes were divided into 6 categories, including metabolism, cell transport, transcriptional regulation, plant hormone regulation, structural proteins, and cell division, which directly or indirectly regulate leaf and midrib architecture trait development in maize.

Discussion

The ultimate goal of maize breeding is not only to improve the quality of grains, but more importantly to enhance the grain yield. The realization of these goals depends on the excellent overall agronomic traits of crop plants. Leaf architecture directly influences canopy structure, consequentially affecting yield. Therefore, it is necessary to clarify the relationship between quality traits, agronomic traits and yield traits. Leaf midrib is responsible for the mutual transport of nutrients and water between the leaf and stalk of plant. It also provides mechanical support to leaf blades and controls the angle of leaf facing sunlight, which influences the efficiency of leaf photosynthesis. The correlation analysis suggests that the ear leaf midrib architecture traits develop synchronously with leaf blade traits and interplay each other; the midrib traits are also significantly correlated with GW and other grain formation traits (Figure 5). The transport capacity and rigidity of leaf midrib are affected by its architecture. Therefore, defining the suitable midrib architecture and dissecting its genetic basis is very important for ideal crop architecture breeding by molecular marker assisted selection.

Leaf midrib architecture is genetically controlled

Current researches have shown that leaf venation is a quantitative trait under polygenic control, with high heritability (Al-Shugeairy, 2016; Rishmawi et al., 2017). In this study, the LMA traits of 508 maize inbred lines manifest continuous distribution, and high broad-sense heritability (Table 1), indicating that the three midrib architecture traits in this panel are conditioned mainly by genetic factors. The value of MBT is stable in two environments, which indicates that the trait has good recurrence in different environments and further demonstrates that it is feasible to map genes related to leaf midrib architecture traits. The larger mean width of leaf midrib appears in the tropical/subtropical origins inbred lines (Figure 4), suggesting that this subgroup possesses abundant or special genes to promote the leaf midrib width development. Although the mean values of TMT and MBT do not show significant difference among subgroups, their significant variations are observed in one single subgroup. Therefore, it is possible to modify leaf midrib architecture using genetic resources inside the subgroups of maize.

Significant correlation of leaf midrib architecture with reproduction traits

The leaf midrib architecture traits are not only significantly correlated with the length (ELL) and width (ELW) of their own ear leaf blade, but also with the reproduction traits, including ear length (EL), ear diameter (ED), cob diameter (CD), cob grain weight (GW), days to anthesis (DTA), days to silking (DTS) and days to heading (DTH) (Figure 5). These correlations show that well-developed leaf midrib architecture traits facilitate the reproductive process in maize. The leaf midrib architecture is a potential aspect of improving whole plant architecture and grain production in maize.

Leaf midrib architecture is regulated by quantitative genes

The leaf midrib architecture is identified as quantitative traits by heritability and phenotypic variation in the population (Table 1; Figure 4). The MBT trait is associated with 4 significant SNPs, the TMT and MW traits are only associated with 1 significant SNP, respectively. However, this does not mean that the TMT and MW traits are not controlled by multiple genes, because many different relative genes may locate in the adjacent region of significant SNP. It has been reported that some similar function or function-related genes locate closely each other on chromosome (Calusinska et al., 2015; Wu et al., 2021). In this study, we screened out several candidate genes regulating LMA traits from the adjacent region of each significant SNP locus (Figure 6), their expression supports their functions in LMA trait development (Figure 2). Twenty-seven from 97 genes identified in 6 significant SNP loci associating LMA are functionally annotated (Figure 2). These candidate genes may participate in modifying gene transcription, regulating auxin functions, small RNA transcription and function, cytoskeleton, cell division and movement, lateral organ growth.

The candidate gene GRMZM2G074124 is associated with MW, presumed to be a YABBY domain transcription factor family protein. YABBY gene family play important roles in the growth and polar establishment of lateral organs such as leaves and floral organs in plants (İnal et al., 2017). The loss-of-function mutations of dl gene belonging to this family cause lack of midrib and the drooping leaf phenotype (Ishikawa et al., 2009) (Yamada et al., 2011). In maize, pleiotropic mutations in drl1 was reported to affect leaf length and width, leaf angle (Strable et al., 2017). This study adds a new function of YABBY domain transcription factor to regulate leaf midrib architecture.

GRMZM2G130953, an associated gene with MBT, encodes IAA27 protein, an auxin-response transcription factor, and GRMZM2G332390 is located in the significant SNP locus of MW, encoding ZmSAUR (SMALL AUXIN UP RNA) 48 protein (Chen et al., 2014), both belonging to early auxin response genes (Hagen and Guilfoyle, 2002). IAA27 protein can regulate the activity of ARF transcription factors (Galli et al., 2015) and SAUR gene is mainly expressed in the hypocotyl and other elongating tissues (Stamm and Kumar, 2013; Favero et al., 2016; Spartz et al., 2017). Therefore, GRMZM2G130953 and GRMZM2G332390 may participate in the development of leaf midrib width and base thickness…. which binds to calmodulin in a calcium-dependent manner, suggesting that they may mediate responses to changes in intracellular calcium level (Yang and Poovaiah, 2000).

GRMZM2G109787 and GRMZM2G481730, both associated with MBT, encode RNA-dependent RNA polymerases (RDR) to synthesize sRNAs (Qian et al., 2011; Kelliher and Walbot, 2014). The two main types of small RNAs are microRNAs (miRNAs) and small interfering RNAs (siRNAs) (Carthew and Sontheimer, 2009). GRMZM2G045977, also associated with MBT, is a target gene of miR396, a growth-regulating factor for maize (Liu et al., 2014), suggesting that miR396 regulates midrib formation (Debernardi et al., 2012).

The candidate gene GRMZM2G078612, associated with MW, encodes γ tubulin -2 chain; GRMZM5G832989, associated with MBT, encodes microtubule-associated protein 70-1. Microtubule (MT) is an important component of cytoskeleton. GRMZM2G407517, associated with MBT, encodes Actin7 is involved in gibberellin and brassinosteroid biosynthesis, auxin transport and cellulose synthesis (Wei et al., 2018). GRMZM2G079185, associated with MW, encodes the LOB (lateral organ boundaries) domain protein, which regulate lateral organ growth in plants (Zhang et al., 2014; Gombos et al., 2017).

Conclusion

The leaf midrib architecture traits width (MW), total thickness (TMT) and base thickness (MBT) have higher heredity and are controlled by quantitative genes in maize. These leaf traits are also significant correlated with reproduction traits such as days to anthesis (DTA), days to silking (DTS), days to heading (DTH), ear length (EL), ear diameter (ED), cob diameter (CD), and cob grain weight (GW). Six SNP loci are significantly associated with LMA, and 97 genes are identified in the SNP location interval.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

DD: Writing – original draft, Writing – review & editing. SJ: Data curation, Writing – original draft. YL: Visualization, Data curation, Writing – original draft. CB: Writing – original draft, Investigation. YZ: Writing – original draft, Project administration. XZ: Writing – review & editing, Writing – original draft. AZ: Writing – review & editing, Writing – original draft. YR: Writing – review & editing, Project administration, Writing – original draft.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Science Foundation of China (31801442), Major Special Sub project of Liaoning Provincial Department of Science and Technology Project (2024JHI/11700007-3), Liaoning Province Germplasm Innovation Grain Storage Special Program (2023JH1/10200009, 2023JH1/1020000903) and Shenyang Key Laboratory for Genome Selection Breeding of Maize.

Acknowledgments

We thank the three reviewers for critical reading of the manuscript and constructive comments. We gratefully acknowledge the invaluable support and contributions of the following SRA samples/projects, which provided essential data for our analyses and conclusions: SRS726671, SRS726672, SRS726811, SRS726802, SRS726799, SRS729984, SRS729963, SRS729962, SRS726704, SRS726705, SRS291020, SRS291021, SRS291022, SRS291073, SRS291074, SRS291075, SRS352445, SRS352446, SRS352447, SRS352451, SRS352452, SRS352453, SRS352454, SRS352455, SRS352456, SRS352457, SRS352458, SRS352459, SRS352460, SRS352461, SRS352462, SRS352507, SRS352508, SRS352511, SRS352514, SRS352517, SRS352518, SRS729955, SRS729960, SRS729961.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1671460/full#supplementary-material

Supplementary Table 1 | Genome-wide diversity and haplotype summary of maize inbred lines.

Supplementary Table 2 | The candidate gene information.

Supplementary Table 3 | Information of candidate genes for leaf midrib architecture in maize. acandidate genes were found in maizegdb. MW, midrib width; TMT, total midrib thickness; MBT, midrib base thickness.

References

Abadie, C., Bathellier, C., and Tcherkez, G. (2018). Carbon allocation to major metabolites in illuminated leaves is not just proportional to photosynthesis when gaseous conditions (CO2 and O2) vary. New Phytol. 218, 94–106. doi: 10.1111/nph.14984

PubMed Abstract | Crossref Full Text | Google Scholar

Al-Shugeairy, Z. (2016). Rice diversity panel mapping for identifying genetics behind leaf vein density trait in rice (Oryza sativa L.). Int. J. Appl. Agric. Sci. 2, 49. doi: 10.11648/j.ijaas.20160204.12

Crossref Full Text | Google Scholar

Bagniewska-Zadworna, A., Byczyk, J., Eissenstat, D. M., Oleksyn, J., and Zadworny, M. (2012). Avoiding transport bottlenecks in an expanding root system: Xylem vessel development in fibrous and pioneer roots under field conditions. Am. J. Bot. 99, 1417–1426. doi: 10.3732/ajb.1100552

PubMed Abstract | Crossref Full Text | Google Scholar

Calusinska, M., Hamilton, C., Monsieurs, P., Mathy, G., Leys, N., Franck, F., et al. (2015). Genome-wide transcriptional analysis suggests hydrogenase- and nitrogenase-mediated hydrogen production in Clostridium butyricum CWBI 1009. Biotechnol. Biofuels 8, 27. doi: 10.1186/s13068-015-0203-5

PubMed Abstract | Crossref Full Text | Google Scholar

Carthew, R. W. and Sontheimer, E. J. (2009). Origins and mechanisms of miRNAs and siRNAs. Cell 136, 642–655. doi: 10.1016/j.cell.2009.01.035

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Y., Hao, X., and Cao, J. (2014). Small auxin upregulated RNA (SAUR) gene family in maize: Identification, evolution, and its phylogenetic comparison with Arabidopsis, rice, and sorghum. J. Integr. Plant Biol. 2, 133–150. doi: 10.1111/jipb.12127

PubMed Abstract | Crossref Full Text | Google Scholar

Dai, L., Wu, L., Dong, Q., Zhang, Z., Wu, N., Song, Y., et al. (2017). Genome-wide association study of field grain drying rate after physiological maturity based on a resequencing approach in elite maize germplasm. Euphytica 213, 1–12. doi: 10.1007/s10681-017-1970-9

Crossref Full Text | Google Scholar

Debernardi, J. M., Rodriguez, R. E., Mecchia, M. A., and Palatnik, J. F. (2012). Functional specialization of the plant miR396 regulatory network through distinct microRNA-target interactions. PLoS Genet. 8, e1002419. doi: 10.1371/journal.pgen.1002419

PubMed Abstract | Crossref Full Text | Google Scholar

de Jong, G., Pamplona, A. K. A., Von Pinho, R. G., and Balestre, M. (2018). Genome-wide association analysis of ear rot resistance caused by Fusarium verticillioides in maize. Genomics 110, 291–303. doi: 10.1016/j.ygeno.2017.12.001

PubMed Abstract | Crossref Full Text | Google Scholar

Dumlao, M. R., Darehshouri, A., Cohu, C. M., Muller, O., Mathias, J., Adams, W. R., et al. (2012). Low temperature acclimation of photosynthetic capacity and leaf morphology in the context of phloem loading type. Photosynth Res. 113, 181–189. doi: 10.1007/s11120-012-9762-5

PubMed Abstract | Crossref Full Text | Google Scholar

Eshed, Y., Izhaki, A., Baum, S. F., Floyd, S. K., and Bowman, J. L. (2004). Asymmetric leaf development and blade expansion in Arabidopsis are mediated by KANADI and YABBY activities. Development 131, 2997–3006. doi: 10.1242/dev.01186

PubMed Abstract | Crossref Full Text | Google Scholar

Favero, D. S., Jacques, C. N., Iwase, A., Le, K. N., Zhao, J., Sugimoto, K., et al. (2016). SUPPRESSOR of PHYTOCHROME b4-3 represses genes associated with auxin signaling to modulate hypocotyl growth. Plant Physiol. 171, 2701–2716. doi: 10.1104/pp.16.00405

PubMed Abstract | Crossref Full Text | Google Scholar

Galli, M., Liu, Q., Moss, B. L., Malcomber, S., Li, W., Gaines, C., et al. (2015). Auxin signaling modules regulate maize inflorescence architecture. Proc. Natl. Acad. Sci. 112, 13372–13377. doi: 10.1073/pnas.1516473112

PubMed Abstract | Crossref Full Text | Google Scholar

Gombos, M., Zombori, Z., Szécsényi, M., Sándor, G., Kovács, H., Györgyey, J., et al. (2017). Characterization of the LBD gene family in Brachypodium: A phylogenetic and transcriptional study. Plant Cell Rep. 36, 61–79. doi: 10.1007/s00299-016-2057-0

PubMed Abstract | Crossref Full Text | Google Scholar

Hagen, G. and Guilfoyle, T. (2002). Auxin-responsive gene expression: Genes, promoters and regulatory factors. Plant Mol. Biol. 49, 373–385. doi: 10.1023/A:1015207114117

PubMed Abstract | Crossref Full Text | Google Scholar

Hoopes, G. M., Hamilton, J. P., Wood, J. C., Esteban, E., Pasha, A., Vaillancourt, B., et al. (2019). An updated gene atlas for maize reveals organ-specific and stress-induced genes. Plant J. 97, 1154–1167. doi: 10.1111/tpj.14184

PubMed Abstract | Crossref Full Text | Google Scholar

İnal, B., Büyük, İ., İlhan, E., and Aras, S. (2017). Genome-wide analysis of Phaseolus vulgaris C2C2-YABBY transcription factors under salt stress conditions. 3 Biotech. 7, 302. doi: 10.1007/s13205-017-0933-0

PubMed Abstract | Crossref Full Text | Google Scholar

Ishikawa, M., Ohmori, Y., Tanaka, W., Hirabayashi, C., Murai, K., Ogihara, Y., et al. (2009). The spatial expression patterns of DROOPING LEAF orthologs suggest a conserved function in grasses. Genes Genet. Syst. 84, 137–146. doi: 10.1266/ggs.84.137

PubMed Abstract | Crossref Full Text | Google Scholar

Kelliher, T. and Walbot, V. (2014). Maize germinal cell initials accommodate hypoxia and precociously express meiotic genes. Plant J. 77, 639–652. doi: 10.1111/tpj.12414

PubMed Abstract | Crossref Full Text | Google Scholar

Kump, K. L., Bradbury, P. J., Wisser, R. J., Buckler, E. S., Belcher, A. R., Oropeza-Rosas, M. A., et al. (2011). Genome-wide association study of quantitative resistance to southern leaf blight in the maize nested association mapping population. Nat. Genet. 43, 163–168. doi: 10.1038/ng.747

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Qin, C., Chen, Z., Zuo, T., Yang, X., Zhou, H., et al. (2014). Identification of miRNAs and their target genes in developing maize ears by combined small RNA and degradome sequencing. BMC Genomics 15, 25. doi: 10.1186/1471-2164-15-25

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Huang, M., Fan, B., Buckler, E. S., and Zhang, Z. (2016). Iterative usage of fixed and random effect models for powerful and efficient Genome-Wide association studies. PLoS Genet. 12, e1005767. doi: 10.1371/journal.pgen.1005767

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, B., Sahito, J. H., Zhang, H., Zhao, J., Yang, G., Wang, W., et al. (2024). SPX family response to low phosphorus stress and the involvement of ZmSPX1 in phosphorus homeostasis in maize. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1385977

PubMed Abstract | Crossref Full Text | Google Scholar

Ohmori, Y., Abiko, M., Horibata, A., and Hirano, H. (2008). A transposon, ping, is integrated into intron 4 of the DROOPING LEAF gene of rice, weakly reducing its expression and causing a mild drooping leaf phenotype. Plant Cell Physiol. 49, 1176–1184. doi: 10.1093/pcp/pcn093

PubMed Abstract | Crossref Full Text | Google Scholar

Qian, Y., Cheng, Y., Cheng, X., Jiang, H., Zhu, S., Cheng, B., et al. (2011). Identification and characterization of Dicer-like, Argonaute and RNA-dependent RNA polymerase gene families in maize. Plant Cell Rep. 30, 1347–1363. doi: 10.1007/s00299-011-1046-6

PubMed Abstract | Crossref Full Text | Google Scholar

Rishmawi, L., Bühler, J., Jaegle, B., Hülskamp, M., and Koornneef, M. (2017). Quantitative trait loci controlling leaf venation inArabidopsis. Plant Cell Environ. 40, 1429–1441. doi: 10.1111/pce.12938

PubMed Abstract | Crossref Full Text | Google Scholar

Rolland-Lagan, A. and Prusinkiewicz, P. (2005). Reviewing models of auxin canalization in the context of leaf vein pattern formation in Arabidopsis. Plant J. 44, 854–865. doi: 10.1111/j.1365-313X.2005.02581.x

PubMed Abstract | Crossref Full Text | Google Scholar

Sahito, J. H., Zhang, H., Gishkori, Z., Ma, C., Wang, Z., Ding, D., et al. (2024). Advancements and prospects of Genome-Wide association studies (GWAS) in maize. Int. J. Mol. Sci. 25, 1918. doi: 10.3390/ijms25031918

PubMed Abstract | Crossref Full Text | Google Scholar

Sobrado, M. A. and Turner, N. C. (1986). Photosynthesis, dry matter accumulation and distribution in the wild sunflower Helianthus petiolaris and the cultivated sunflower Helianthus annuus as influenced by water deficits. Oecologia 69, 181–187. doi: 10.1007/BF00377619

PubMed Abstract | Crossref Full Text | Google Scholar

Spartz, A. K., Lor, V. S., Ren, H., Olszewski, N. E., Miller, N. D., Wu, G., et al. (2017). Constitutive expression of Arabidopsis SMALL AUXIN UP RNA19 (SAUR19) in tomato confers Auxin-Independent hypocotyl elongation. Plant Physiol. 173, 1453–1462. doi: 10.1104/pp.16.01514

PubMed Abstract | Crossref Full Text | Google Scholar

Stamm, P. and Kumar, P. P. (2013). Auxin and gibberellin responsive Arabidopsis SMALL AUXIN UP RNA36 regulates hypocotyl elongation in the light. Plant Cell Rep. 32, 759–769. doi: 10.1007/s00299-013-1406-5

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, J. D., Perkins, A., Williams, W. P., and Warburton, M. L. (2015). Using genome-wide associations to identify metabolic pathways involved in maize aflatoxin accumulation resistance. BMC Genomics 16, 5–22. doi: 10.1186/s12864-015-1874-9

PubMed Abstract | Crossref Full Text | Google Scholar

Visscher, P. M., Wray, N. R., Zhang, Q., Sklar, P., McCarthy, M. I., Brown, M. A., et al. (2017). 10 years of GWAS discovery: biology, function, and translation. Am. J. Hum. Genet. 101, 5–22. doi: 10.1016/j.ajhg.2017.06.005

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, L., Zhang, X., Zhang, Z., Liu, H., and Lin, Z. (2018). A new allele of the Brachytic2 gene in maize can efficiently modify plant architecture. Heredity 121, 75–86. doi: 10.1038/s41437-018-0056-3

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, T., Wei, B., Lin, H., Zhou, B., Lin, T., Liu, Q., et al. (2021). Integral analyses of competing endogenous RNA mechanisms and DNA methylation reveal regulatory mechanisms in osteosarcoma. Front. Cell Dev. Biol. 9. doi: 10.3389/fcell.2021.763347

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, C., Zhang, H., Sun, J., Guo, Z., Zou, C., Li, W., et al. (2018). Genome-wide association study dissects yield components associated with low-phosphorus stress tolerance in maize. Theor. Appl. Genet. 131, 1699–1714. doi: 10.1007/s00122-018-3108-4

PubMed Abstract | Crossref Full Text | Google Scholar

Xue, Y., Warburton, M. L., Sawkins, M., Zhang, X., Setter, T., Xu, Y., et al. (2013). Genome-wide association analysis for nine agronomic traits in maize under well-watered and water-stressed conditions. Theor. Appl. Genet. 126, 2587–2596. doi: 10.1007/s00122-013-2158-x

PubMed Abstract | Crossref Full Text | Google Scholar

Yamada, T., Yokota, S. Y., Hirayama, Y., Imaichi, R., Kato, M., Gasser, C. S., et al. (2011). Ancestral expression patterns and evolutionary diversification of YABBY genes in angiosperms. Plant J. 67, 26–36. doi: 10.1111/j.1365-313X.2011.04570.x

PubMed Abstract | Crossref Full Text | Google Scholar

Yamaguchi, T., Nagasawa, N., Kawasaki, S., Matsuoka, M., Nagato, Y., Hirano, H., et al. (2004). TheYABBY GeneDROOPING LEAF regulates carpel specification and midrib development inOryza sativa. Plant Cell 16, 500–509. doi: 10.1105/tpc.018044

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, X., Gao, S., Xu, S., Zhang, Z., Prasanna, B. M., Li, L., et al. (2011). Characterization of a global germplasm collection and its potential utilization for analysis of complex quantitative traits in maize. Mol. Breed. 28, 511–526. doi: 10.1007/s11032-010-9500-7

Crossref Full Text | Google Scholar

Yang, Q., Li, Z., Li, W., Ku, L., Wang, C., Ye, J., et al. (2013). CACTA-like transposable element in ZmCCT attenuated photoperiod sensitivity and accelerated the postdomestication spread of maize. Proc. Natl. Acad. Sci. 110, 16969–16974. doi: 10.1073/pnas.1310949110

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, T. and Poovaiah, B. W. (2000). Molecular and biochemical evidence for the involvement of calcium/calmodulin in auxin action. J. Biol. Chem. 275, 3137–3143. doi: 10.1074/jbc.275.5.3137

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y. M., Zhang, S. Z., and Zheng, C. C. (2014). Genomewide analysis of LATERAL ORGAN BOUNDARIES Domain gene family in Zea mays. J. Genet. 93, 79–91. doi: 10.1007/s12041-014-0342-7

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, G., Hao, D., Chen, G., Lu, H., Shi, M., Mao, Y., et al. (2016). Genome-wide association study of the husk number and weight in maize (Zea mays L.). Euphytica 210, 195–205. doi: 10.1007/s10681-016-1698-y

Crossref Full Text | Google Scholar

Keywords: leaf venation, leaf vein, vascular bundle, single nucleotide polymorphisms (SNPs), YABBY gene

Citation: Dang D, Ji S, Liu Y, Bai C, Zhu Y, Zhang X, Zhang A and Ruan Y (2025) Genome-wide association studies on leaf midrib architecture in maize. Front. Plant Sci. 16:1671460. doi: 10.3389/fpls.2025.1671460

Received: 23 July 2025; Accepted: 09 October 2025;
Published: 30 October 2025.

Edited by:

Muhammad Saeed, Government College University, Faisalabad, Pakistan

Reviewed by:

Javed Hussain Sahito, Henan Agricultural University, China
Krishna Sai Karnatam, West Virginia State University, United States

Copyright © 2025 Dang, Ji, Liu, Bai, Zhu, Zhang, Zhang and Ruan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ao Zhang, emhhbmdhbzdAc3lhdS5lZHUuY24=; Yanye Ruan, eWFueWVydWFuQHN5YXUuZWR1LmNu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.