Skip to main content

Genetic analysis and QTL mapping for pericarp thickness in maize (Zea mays L.)

Abstract

Proper pericarp thickness protects the maize kernel against pests and diseases, moreover, thinner pericarp improves the eating quality in fresh corn. In this study, we aimed to investigate the dynamic changes in maize pericarp during kernel development and identified the major quantitative trait loci (QTLs) for maize pericarp thickness. It was observed that maize pericarp thickness first increased and then decreased. During the growth and formation stages, the pericarp thickness gradually increased and reached the maximum, after which it gradually decreased and reached the minimum during maturity. To identify the QTLs for pericarp thickness, a BC4F4 population was constructed using maize inbred lines B73 (recurrent parent with thick pericarp) and Baimaya (donor parent with thin pericarp). In addition, a high-density genetic map was constructed using maize 10 K SNP microarray. A total of 17 QTLs related to pericarp thickness were identified in combination with the phenotypic data. The results revealed that the heritability of the thickness of upper germinal side of pericarp (UG) was 0.63. The major QTL controlling UG was qPT1-1, which was located on chromosome 1 (212,215,145–212,948,882). The heritability of the thickness of upper abgerminal side of pericarp (UA) was 0.70. The major QTL controlling UA was qPT2-1, which was located on chromosome 2 (2,550,197–14,732,993). In addition, a combination of functional annotation, DNA sequencing analysis and quantitative real-time PCR (qPCR) screened two candidate genes, Zm00001d001964 and Zm00001d002283, that could potentially control maize pericarp thickness. This study provides valuable insights into the improvement of maize pericarp thickness during breeding.

Peer Review reports

Introduction

Pericarp is the outer protective tissue of maize kernel. It protects the embryo and endosperm from mechanical damage or pests. Moreover, it regulates the rate of dehydration and maturation of the kernel and water absorption for germination [1, 2]. Maize pericarp thickness is reported to be associated with seed quality [3], popcorn quality [4], resistance to pathogens [5], seed maturity, and moisture content [6]. Furthermore, pericarp thickness is a quality-related factor for sweet corn breeding. Thin pericarp is a desirable trait for sweet corn because it improves the tenderness [7]. Pericarp thickness negatively correlates with tenderness and influences the resistance to mastication [8]. Therefore, the selection of appropriately thin pericarp is crucial for the quality breeding of maize.

The development process of maize pericarp is divided into four stages: undeveloped, thickening, expansion, and strengthening using the hierarchical clustering and PCA analysis of high temporal resolution transcriptome [9]. Then, it was found that from 15 to 23 day after pollination (DAP), the pericarp thickness of M03 and M08 showed a decreasing trend, and then gradually stabilized from the dynamic measurements of pericarp thickness with these two lines [10]. Certainly, the differences in pericarp thickness of various maize materials were also investigated. The pericarp thickness of 19 maize materials (genotypes Su, su, sh2, and suse) was analyzed and showed that SuSu had the thickest pericarp and suse had the thinnest pericarp [11]. The molecular mechanisms of genes that control pericarp thickness vary. pericarp color 1 gene increases maize pericarp thickness through the accumulation of phlobaphene [12]. And silencing the SlPRE2 gene affected the plant's response to gibberellins resulting in reduced tomato fruit size and pericarp thickness [13]. Gan et al. found that the reduction in tomato pericarp thickness was mainly due to a reduction in cell number, which led to a reduction in cell division. In addition, endogenous cytokinins were demonstrated to regulate pericarp cell division and fruit size using transcriptome analysis [14].

Quantitative trait loci (QTL) mapping is a key step in exploring the genetic mechanisms of the complex quantitative trait of maize pericarp thickness. The genes controlling maize pericarp thickness ranged from 1.4 to 5.9, with an average narrow-sense heritability of 55.2%, and pericarp thickness is a highly heritable trait that is not strongly influenced by the environment [15, 16]. Helm et al. classified the pericarp into six positions: upper, middle and lower part of the germinal and non-germinal surfaces. And they found that the narrow-sense heritability of pericarp thickness was estimated to be up to 80% [17]. Choe et al. measured pericarp thickness at five different positions on maize and reported that the thickness at all five different positions exhibited high heritability and correlation and most QTLs were linked to more than one pericarp thickness position. This suggested that the pericarp thickness at different positions on maize is controlled by common genes with multiple effects [16]. Wang et al. used different molecular markers to investigate the QTLs that control pericarp thickness in maize and reported that three genes were associated with pericarp thickness on chromosomes 1, 2, and 6 [18]. Park et al. found four QTLs associated with pericarp thickness, which were located on chromosomes 4, 5, 8, and 9 [19].

In this study, we first investigated the dynamics of maize pericarp thickness during kernel development in normal corn, sweet corn, and waxy corn inbred lines. In addition, we constructed the BC4F4 population using two maize inbred lines with different pericarp thicknesses [B73 (thick pericarp) and Baimaya (thin pericarp)] and identified the QTLs for pericarp thickness based on inclusive composite interval mapping. The results of our study may provide a theoretical basis for changing the thickness of maize pericarp during maize breeding.

Materials and methods

Plant materials and field experiments

The materials used in this study was provided by the maize group of South China Agricultural University. To investigate the dynamic changes of maize pericarp thickness, sweet corn inbred line ZF1, waxy corn inbred line TN113, and normal corn inbred line N75 were used as materials. They were planted in Zengcheng District, Guangzhou City, Guangdong Province, China (ZC, approximately 113°E and 23°N) in August 2022. In terms of QTL mapping for pericarp thickness, two parentlines B73 and Baimaya were selected to construct the BC4F4 population. B73 (thick pericarp) was used as a recurrent parent, and Baimaya (thin pericarp) and was choosed as donor parent. In February 2022, 318 BC4F4 families and their parents were planted in Huanggang City, Hubei Province, China (HG, approximately 115°E and 30°N), and Gucheng County, Xiangyang City, Hubei Province, China (GC, approximately 111°E and 31°N). The 318 BC4F4 lines and their parents were also planted in August 2022 in Zengcheng District, Guangzhou City, Guangdong Province, China (approximately 113°E and 23°N). In total, 10 plants of each family line or parent were grown with 2 rows. The length of the rows was 3 m. The distance between rows was 70 cm, and the distance between plants was 25 cm. The ears were covered with sulfuric acid paper bags to ensure self-fertilization before the silking stage. Locally recommended herbicide and pesticide programs were followed for crop management.

Pericarp collection and phenotypic analysis

To monitor the changes in maize pericarp, the successfully fertilized kernels of ZF1, TN113, and N75 were collected every 2 days from the 5th day after self-pollination. For QTL mapping, the ears from the BC4F4 population and parents were collected at 21 DAP. Three self-pollinated ears were sampled from each genotype and stored at -20°C.

Ten kernels were collected from the middle of each ear from the above three self-pollinated ears. One-third of the head of the corn kernel was transversely cut off with a double-edged razor blade, followed by slicing the samples with a thickness of < 1 mm on a slide and pipetting 30 μL Safranin O solution (Diluted to 0.05 g/L. Roles-Bio, Guangzhou, China) into the center of the samples. After staining for 10 to 15 s, excess staining solution was removed using a filter paper, and the stained sections were observed under a Zeiss light microscope (Zeiss, Shanghai, China). The thickness of upper germinal side and upper abgerminal side of pericarp and the number of pericarp cell layers were recorded. Dividing pericarp thickness by the number of cell layers as a reference value for individual pericarp cell thickness. The observed values of pericarp thickness and the number of pericarp cell layers for each genotype were determined as the average values obtained from three ears, and these averaged values were utilized for subsequent analysis.

IBM SPSS Statistics 26 software (IBM Corp. Released 2019; IBM SPSS Statistics for Windows, Version 26.0. Armonk, NY: IBM Corp) was used to perform descriptive statistics, normal distribution plotting, normality tests (P-P plots), and Pearson’s correlation analysis [20, 21]. The broad-sense heritability (\({h}_{b}^{2}\)) was estimated as \({h}_{b}^{2}={V}_{ g}^{2}/\left({V}_{ g}^{2}+{V}_{ ge}^{2}/e+{V}_{ e}^{2}/re\right)\). \({V}_{ g}^{2}\) is the genetic variance. \({V}_{ e}^{2}\) is the residual variance. \({V}_{ ge}^{2}\) is the genotype-by-environment interaction. e represents the number of environments and r represents the number of replicates in each environment [22].

Construction of linkage maps and QTL mapping for pericarp thickness in maize

The DNA of the two parents and the BC4F4 populations was extracted using the cetyltrimethylammonium bromide method [23], and genotyping of the DNA was performed using the maize 10 K liquid chip (MolBreeding, Shijiazhuang, China). A total of 10,910 SNP markers were obtained after excluding markers with heterozygosity greater than 0.1, biased segregation ratio greater than 0.3, and minimum allele frequency (MAF) less than 5%. 2,834 high quality SNPs were used for further research. The parameter recombination frequency of QTL IciMapping (version 4.1) [24] was set to be less than 0.35, and the Kosambi mapping function was used for the calculation of genetic map distances, and QTL for pericarp thickness were analyzed using inclusive composite interval mapping. The genetic distance information was imported into R4.2.1 (R Core Team (2023). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/.) and plotted as a 2D genetic map using LinkageMapView’s “lmv.linkage.plot” command [25]. QTL nomenclature was in the form of “q + abbreviation of trait name + chromosome number” [26].

Functional structural domain prediction and DNA sequencing analysis of candidate genes

Prediction of functional structural domains using the MaizeGDB database (B73 RefGen_V4) [27] and InterPro website [28]. The primers of candidate genes for PCR amplification were designed using NCBI primer-BLAST [29] and synthesized by Tsingke Biotech Co., Ltd (Tsingke, Guangzhou, China). Table S1 provides the primer sequences for the candidate genes. The PCR reaction volume was 25.0 μl containing 200 ng of genomic DNA, 1 μl of primer (10 μM) and 10 μl of 2 × plus Taq HiFi PCR mix (Mikx, Guangzhou, China). The quality of the amplification products was assessed by electrophoresis on a 1.2% agarose gel. Subsequently, the PCR products meeting the criteria were forwarded to Tsingke Biotech Co., Ltd. (Tsingke, Guangzhou, China) for sequencing. Sequence comparisons and analysis were performed using SnapGene software (www.snapgene.com). Schematic representation of gene structure and protein functional structural domains done by IBS software [30].

Quantitative real-time PCR

RNA from the seedling leaf of B73 and Baimaya was extracted with RNA extraction reagent (Surbiopure, Guangzhou, China) and reverse transcribed into cDNA with All-in-one First-Strand Synthesis MasterMix with dsDNase (Xinkailai, Guangzhou, China). The cDNA was diluted 5 times and used as the template for quantitative real-time PCR (qPCR). DUF was used as a housekeeping gene [31]. qPCR was performed using CFX Connect™ Real-Time PCR Detection System (Bio-Rad, USA) and ChamQ SYBR qPCR Master Mix without ROX (Vazyme, Guangzhou, China). Quantification of qPCR results using the 2CT method [32]. Primer sequences for related genes are listed in Table S2. Three biological replicates were carried out.

Results

Dynamic changes in pericarp thickness in sweet, waxy, and normal maize during kernel development

There were significant differences in the thickness of the pericarp of sweet, waxy, and normal maize, but the trends in thickness showed consistency, with all showing an increasing and then gradually decreasing trend. Throughout the entire process of maize pericarp development, the number of mesocarp cells is the highest, and their morphological changes are the most dramatic. They constitute the major component that influences pericarp thickness, typically consisting of 5 to 20 cell layers (Fig. 1). In the initial stages following fertilization, pericarp thickness rapidly increases, which can be divided into two parts based on cell morphology: cell division and cell elongation. Before DAP 5, the increase in pericarp thickness is mainly due to the continuous division of mesocarp cells, resulting in a rapid increase in cell number. From DAP 5 to DAP 9, the number of cell layers stabilizes between 19–22 layers. Meanwhile, the cells start elongating both horizontally and vertically. During this stage, the increase in pericarp thickness is primarily caused by cell expansion (Fig. S1). After reaching its peak, there is a sharp decrease in pericarp thickness, which is also a result of rapid changes in the number of mesocarp cells. From DAP 11 to DAP 13, mesocarp cells undergo autophagy from the inside out, causing a significant disappearance of cells and the formation of cavities in the pericarp. The number of cell layers in the pericarp decreases from the initial 19–22 layers to 9–11 layers. As the cell number decreases and the substances in the cavities are absorbed by the endosperm, pericarp thickness rapidly decreases. In the subsequent period, changes in mesocarp cells tend to stabilize, and pericarp thickness no longer exhibits drastic changes but shows a gradual downward trend.

Fig. 1
figure 1

Changes in maize pericarp thickness during the kernel developmental process. A Microscopic observation of the pericarp development in maize. EP: exocarp. EN: endosperm. CA: the cavities formed by apoptotic cells. B Variation in pericarp thickness after self-pollination of three different corn types. C Variation in the number of cell layers after self-pollination. ZF1: sweet corn inbred line; N75: normal corn inbred line; TN113: waxy corn inbred line

Studying the changes in the developmental process of the maize pericarp revealed that the process can be divided into four stages: growth, formation, autophagy, and maturation (Figs. 1 and 2). (1) Growth stage: From the time of self-pollination till DAP 5, the number of dividing pericarp cells, mainly mesocarp cells, rapidly increases. At the same time, these newly generated pericarp cells begin to synthesize starch. (2) Formation stage: from DAP 5 to DAP 9, the pericarp cells stop dividing and start to elongate. At the same time, large amounts of starch are present in the pericarp cells. The primary role of the pericarp cell at this stage is to synthesize and store starch. (3) Autophagy stage: from DAP 10 to 26, the pericarp cells stop synthesizing starch. Autophagy starts to occur in the mesocarp cells, the cells break down, and the contents flow out. This creates a cavity in the pericarp, and the starch stored in the cells is gradually absorbed by the endosperm. (4) Maturation stage: After DAP 27, the nucleus and contents of the pericarp cells mostly disappear, indicating that the pericarp cells have died. The major function of the pericarp at this stage is to protect the embryo and endosperm.

Fig. 2
figure 2

Changes in maize pericarp cells during various developmental stages. EP: exocarp. SG: starch grain. AL: aleurone layer

In summary, this study revealed that the differences in pericarp thickness among the three selected maize inbred lines and the differences in pericarp thickness at different developmental stages of individual kernels are mainly determined by the individual pericarp cell thickness and number of mesocarp cells.

Phenotypic and genetic heritability analyses of pericarp thickness in maize

The phenotypic correlations of pericarp thickness in three different environments were positive and highly significant (Fig. 3). The correlation coefficients for the thickness of upper germinal side of pericarp (UG) and the thickness of upper abgerminal side of pericarp (UA) were 0.235–0.415 and 0.354–0.508, respectively. UG and UA (Fig. S2) was significantly correlated (P < 0.01) in the three different environments, indicating that both UG and UA can reflect the overall level of pericarp thickness in each individual line.

Fig. 3
figure 3

Pearson correlation coefficients for pericarp thickness in three different environments. The diagonal line represents the frequency histogram of peel thickness. The lower left of the diagonal (lower triangle) is a scatter plot of the expression values of the two samples. The greater the slope of the red fitted curve, the stronger the correlation between the two samples. In the upper right of the diagonal (upper triangle), the numbers indicate the correlation values between the two samples and * indicates the level of significance (* P < 0.05, ** P < 0.01, *** P < 0.001). HG: Huanggang, GC: Gucheng, ZC: Zengcheng, UG: The thickness of upper germinal side of pericarp, UA: The thickness of upper abgerminal side of pericarp

There was a significant difference in pericarp thickness between the two parents, with the pericarp of the paternal Baimaya being thinner than that of the maternal B73 (Figs. 4 and 5). The main reason for this difference was the different thickness of individual pericarp cells (Fig. S3). The thickness of the UG was slightly more than that of the UA, which was consistent across the three environments.

Fig. 4
figure 4

Sectioning of B73 and Baimaya to determine pericarp thickness. The thickness of the upper germinal side of the pericarp of (A) B73 and (B) Baimaya. The thickness of the upper abgerminal side of the pericarp of (C) B73 and (D) Baimaya

Fig. 5
figure 5

Analysis of significant differences in pericarp thickness between B73 and Baimaya. A The thickness of upper germinal side of pericarp. B The thickness of upper germinal side of pericarp. Different letters indicate significant differences at the P < 0.05 level

The broad-sense heritability (h) of the UG and UA were 0.63 and 0.70, respectively (Table 1). The segregation of pericarp thickness conformed to a normal distribution (Fig. S4) and was a typical quantitative trait under polygenic control. The pericarp thickness distribution in the BC4F4 population was continuous, and considerable transgressive segregation was observed, indicating that both parents contributed to alleles for pericarp thickness.

Table 1 Pericarp thickness in the two parents and BC4F4 population

Identification of primary QTLs for pericarp thickness in the BC4F4 family lines

A genetic map was constructed using 2,824 polymorphic SNP markers (Fig. S5). There were 10 linkage groups in the map with a total distance of 4506.87 centimorgans (cM), with the longest linkage group being 749.56 cM and the shortest being 292.17 cM with an average distance of 1.57 cM. The maximum number of SNP markers was 414 in linkage group 1, and the minimum number was 199 in linkage group 7. The proportion of marker gaps of < 5 cM was calculated for each cluster to measure the degree of linkage between markers, with higher proportions indicating a more even distribution of markers in that cluster (Table S3). The highest proportion of gap of < 5 cM was 99.64% for linkage group 9, and the lowest was 93.89% for linkage group 2. Additionally, the largest gap of this map was located on linkage group 2 at 39.82 cM.

The QTL mapping for pericarp thickness in maize was based on the phenotypic data of the two replicates and three locations. The results of correlation analysis of pericarp thickness revealed that UG was significantly correlated with that of UA (Fig. 3). This indicated that the pericarp thickness of the two sides might be controlled by different QTLs, therefore, the data of the two sides were analyzed using QTL mapping. A total of 17 QTLs for pericarp thickness were detected in the BC4F4 population, distributed on maize chromosomes 1, 2, 4, 5, 8, 9, and 10 (Table 2).

Table 2 QTLs detected for pericarp thickness in the BC4F4 population

A total of 10 QTLs for UG were identified, and the phenotypic variation explained by a single QTL ranged from 2.15% to 10.69%. qPT1-1, qPT2-1, qPT10-1, and qPT8-2 exhibited negative additive effects. This indicated that the alleles of these QTLs that reduce pericarp thickness originated from the B73 parent, whereas the alleles carried by the Baimaya parent acted to reduce pericarp thickness in the QTLs with positive additive effects. Among them, qPT8-1 exhibited the highest additive effect of 15.56. qPT1-1 was identified in three different environments. This QTL is located at the physical position 212,215,145–212,948,882 on chromosome 1, explained 4.93%–10.69% of the phenotypic variance, and had high LOD values ranging from 7.96 to 10.46.

For UA, eight QTLs were identified. qPT2-1 identified at Gucheng exhibited the highest LOD value of 8.20. The additive effect was − 13.13, representing that the allele for this QTL that reduces pericarp thickness originated from the B73 parent. Furthermore, qPT2-1 was identified in both UG and UA. It is speculated that the pericarp thickness on both sides of the maize kernel may be controlled by this QTL. In summary, qPT1-1 and qPT2-1 were speculated to be stable QTLs controlling maize pericarp thickness (Fig. 6, Table 2).

Fig. 6
figure 6

QTL mapping of maize pericarp thickness. The horizontal dotted line in the figure is the threshold line for a LOD value of 2.5, and the reddish labelled font is the major QTL

DNA sequencing of candidate genes and prediction of functional structural domains

The B73 reference sequence (B73 RefGen_V4) from the Maize Reference Genome Database was used to obtain genetic and functional annotation information within QTL intervals. In total, 14 and 662 genes were found within qPT1-1 and qPT2-1, respectively (Tables S4 and S5). Through functional annotation and DNA sequence analysis, candidate genes (Zm00001d001964 and Zm00001d002283) potentially controlling pericarp thickness were identified.

Through the prediction of functional domains, Zm00001d001964 (ZmSAUR15) has an Auxin_inducible domain (Fig. 7A). It is well known that auxin influences cell division and cell elongation, thereby affecting plant growth and development, and consequently the thickness of the maize pericarp. Many members of the SMALL AUXIN UP RNA (SAUR) gene family are among the most quickly and significantly triggered by auxin [33]. The results of the DNA sequencing analysis indicate that the Zm00001d001964 gene of Baimaya differs from B73 by three base substitutions, located at positions 382 bp, 400 bp, and 417 bp within the coding sequence (CDS). These substitutions result in the conversion of Serine to Alanine, Glycine to Arginine, and Lysine to Asparagine.

Fig. 7
figure 7

Schematic gene structure and protein functional structure domains (A) Zm00001d001964 (B) Zm00001d002283

Furthermore, the analysis of functional domains revealed that Zm00001d002283 has a PKinase domain (Fig. 7B). Protein phosphorylation plays a vital role in regulating cell cycle progression, especially through the activity of cyclin-dependent protein kinases, which are key players in controlling cell division [34]. The DNA sequencing analysis of the Zm00001d002283 gene reveals eleven base substitutions, two deletions, and one insertion in Baimaya compared to B73 (Table S6). Among these changes, four base substitutions were identified within the CDS at positions 2224 bp, 3262 bp, 3268 bp, and 3270 bp. These substitutions cause Valine to change to Alanine, Valine to Isoleucine, and Glutamine to Asparagine. The change from Valine to Alanine occurs within the PKinase domain. The sequencing peaks of some SNPs are shown in Fig. S6.

In addition, the results of qPCR showed that the gene expression of Zm00001d001964 was higher in B73 than that in Baimaya, while the gene expression of Zm00001d002283 was higher in Baimaya (Fig. S7). This suggests that high expression of Zm00001d001964 promotes pericarp thickness, while high expression of Zm00001d002283 reduces pericarp thickness.

Discussion

Cytological changes in the pericarp of different varieties of maize during kernel development

Generally, four methods are used for determining the thickness of maize pericarp: (1) Microscopic method in which sections of frozen seeds are stained with a dye and observed under a calibrated microscope [35]. (2) Weighing the whole pericarp and reporting it as a percentage of the total dry kernel weight [36].(3) Micrometer method, in which the pericarp is first peeled from the kernels and dried at room temperature, and finally, the thickness is measured using a micrometer [37]. (4) Method using penetrometer (an instrument used to assess the hardness of the fruit), which quantifies the thickness of the pericarp by calculating the ratio of spring pressure to surface area [38]. Microscopic method was used in this study. The pericarp thickness of different maize varieties exhibited significant variation. At DAP 27, reaching a minimum of 59 μm in sweet maize and a maximum of 139 μm in waxy maize. Notably, there was an observable trend of initial increase followed by decrease in pericarp thickness across different maize varieties. The trend of pericarp thickness was similar to that of Zhang et al. [9]. Furthermore, by analyzing the changes in pericarp thickness and disparities in pericarp cell structure, maize pericarp development has been classified into four stages: growth, formation, autophagy, and maturation. These findings are consistent with the research on other grass species, such as wheat [39].

During growth and formation, the thickness of the pericarp gradually increases. The main reason for the thickening of the pericarp during the growth stage is the dramatic increase in cell number due to the continuous division of the pericarp cells. The main physiological activity of the pericarp cells during the formation stage changes from cell division to cell elongation, and the main reason for the continued increase in pericarp thickness at this time is the change in the thickness of individual cells. Large amounts of starch are commonly accumulated in the pericarp cells of cereals [40,41,42]. In this study, large number of starch granules were observed to be synthesized and stored in the pericarp cells during the formative stages of maize pericarp development. From the autophagy stage onward, maize pericarp thickness continued to decrease and eventually stabilized. Programmed cell death (PCD) occurs in pericarp cells during cereal development and endosperm development is also affected by this process [43, 44]. During the autophagy stage of maize pericarp development, mesocarp cells enter into the PCD process. A large number of cells begin to break down, and the contents flow out of the pericarp to form a cavity. Further, the starch and other nutrients in the cavity are transported to the endosperm and stored there. This phenomenon is an important reason for the sharp decrease in the thickness of the pericarp. The maturation stage is the final stage of maize pericarp development. After the autophagy stage, most mesocarp cells die, and the pericarp thickness is close to the minimum value. At this time, almost all the starch in the pericarp cells is transported to the endosperm. With the dehydration of the kernel, the water content of the pericarp cells also gradually decreases, and the overall thickness of the pericarp remains unchanged.

Comparisons of QTLs controlling pericarp thickness in maize

Previous studies have demonstrated a significant negative correlation between pericarp thickness and pericarp tenderness, and that pericarp thickness is controlled by multiple genes [3, 15]. Table 1 shows that the broad-sense heritability of the UG and UA were 0.63 and 0.70, respectively. It indicated that genetic effects had a significant influence on this trait, and environmental effects had little effect on pericarp thickness.

Yu et al. observed eight QTL associated with pericarp thickness on chromosomes 2, 3, 5, 6 and 8 using a BC1F2 population containing 190 families. Among them, the QTLs located in bin2.01, bin5.06 and bin8.05 were in similar regions to qPT2-1, qPT5-1 and qPT8-2 in this study, respectively [45]. In addition, a study involving 264 F2:3 family lines generated from crosses between Korean glutinous maize inbred lines BH 20 and BH 30 detected a total of 33 QTLs for pericarp thickness. These QTLs were found to be located on chromosomes 1, 2, 3, 4, 8, 9, and 10. The marker intervals umc1278 to umc1833 (bin:1.07), umc1329 to dupssr34 (bin:4.06), and umc1691 to umc1771 (bin:9.03) are in similar regions to qPT1-1, qPT4-2, and qPT9-2, respectively [16]. Using pericarp weight as a reference for measuring pericarp thickness, Wanlayaporn et al. identified QTLs associated with maize pericarp thickness on chromosomes 1, 2, 4, 9 and 10 [36]. Wu et al. conducted a study where they developed a BC4F3 population comprising 148 families. They constructed a genetic map with 3876 specific length amplified fragment tags and identified 14 QTLs associated with pericarp thickness. Among these QTLs, qPT10-5 (located at position 144,631,242–145,532,401) was identified as the major locus. Within the qPT10-5 region, they found 42 candidate genes. Through transcriptome analysis, they identified five genes that exhibited differential expression between the two parents. Based on gene annotation information, three of these genes were identified as potential candidate genes for pericarp thickness. Specifically, GRMZM2G143352 was annotated as an AUX/IAA transcription factor, GRMZM2G143402 as a ZIM transcription factor, and GRMZM2G143389 as a protein named FATTY ACID EXPORT 2 chloroplastic [46].

In this study, the pericarp thickness of the upper abgerminal and upper germinal was used as a representative to determine the maize pericarp thickness. A total of 17 QTL were finally identified by combining the phenotypic data on maize pericarp thickness with the genotypic data obtained from the maize 10 K SNP microarray. Two main effect QTLs were identified on chromosomes 1 and 2, respectively, with qPT1-1 explaining 4.93%–12.97% of the phenotypic variance and qPT2-1 explaining 4.04%–8.63% of the phenotypic variance. Furthermore, qPT4-2, qPT5-1, qPT8-2, and qPT9-2 were consistent with the QTL results in other articles, although these loci were not repeatedly detected in this study.

Candidate genes for maize pericarp thickness

In maize, the expression of Zm00001d001964 is influenced by the presence of auxin. Auxin, a vital plant hormone that regulates the growth and development of plants, can induce the expression of auxin-responsive genes, including AUX/IAAs, GH3s, and SAURs [33]. The down-regulation of Sl-IAA17, which acts as an active inhibitor of gene transcription regulated by auxin, results in an increase in tomato pericarp cell size [47]. SAUR19 has been found to have a pivotal role in facilitating cell elongation in Arabidopsis. This effect is achieved by stimulating the activity of the plasma membrane H-ATPase while suppressing the phosphatase type 2C-D subfamily [48]. In addition to SAUR19, there are other SAURs in Arabidopsis that play a positive role in cell expansion, such as SAUR36, SAUR41, and SAUR63 [49,50,51]. Li et al. overexpressed the VvSAUR041 gene from tomato and found that the average cell area of the pericarp tissue was greater in the OE lines [52]. Functional structural domain prediction of Zm00001d001964 revealed an Auxin-inducible structural domain. In contrast to the B73 reference gene, Baimaya exhibits three base substitutions, all occurring within the coding sequence. Auxin regulates cell division and expansion, thus controlling plant growth and development, which may impact maize pericarp thickness. These findings suggest that ZmSAUR15 likely plays a crucial role in regulating pericarp thickness development.

The homologue of Zm00001d002283 in Arabidopsis is AT1G16670, and its encoded product belongs to the protein kinase family. Protein kinases play a role in regulating the phosphorylation state of proteins and affect protein activity through catalysis of substrate proteins, which in turn affects a wide range of plant activities [53]. Protein phosphorylation regulates various biological processes such as signal transduction, growth and development, substance metabolism, and environmental responses by altering the conformation of proteins and affecting their activity and protein–protein interactions [54, 55]. Sequencing analysis of Zm00001d002283 showed that Baimaya had eleven base substitutions, two deletions, and one insertion when compared to the reference gene for B73. Four of these base substitutions were in the CDS. Functional structural domain prediction revealed that Zm00001d002283 has a PKinase domain, and there is a change from Valine to Alanine within the PKinase domain. Protein kinases have a variety of important functions in living organisms. It can regulate the function and activity of proteins through phosphorylation reactions. In living organisms, protein kinases are involved in important biological processes such as cell signaling, cell cycle regulation, cell differentiation and proliferation. Czerednik et al. found that cyclin-dependent kinase (CDK) affect the size and number of cell layers in tomato pericarp cells [56]. A previous study found that overexpression of a zinc finger gene, Solanum lycopersicum PERICARP-ASSOCIATED ZINC FINGER PROTEIN 1 (SlPZF1), enhanced the expression of cyclin-dependent protein kinase (SlCDKB1 and SlCDKB2) in the pericarp tissue and reduces pericarp cell thickness [57].

In summary, analysis of the dynamics of pericarp thickness showed that changes in thickness at different developmental stages of the pericarp were mainly attributed to expansion and cell division of mesocarp cells. Our sequencing results indicate that the Zm00001d001964 and Zm00001d002283 genes of Baimaya show multiple base substitutions at the CDS compared to the B73 reference genome. These substitutions result in amino acid changes, which are associated with specific protein functional domains. In tomato, overexpression of the VvSAUR041 gene increases the average cell area of the pericarp tissue and plays a positive role in fruit development [52]. Whereas, enhanced expression of cell cycle protein-dependent protein kinase in pericarp tissues decreased pericarp thickness and mesocarp cell size [57]. These are consistent with the qPCR results of this experiment (Fig. S7).

Conclusions

This study investigated changes in maize pericarp development and identified major QTLs for maize pericarp thickness. During the development of maize kernel, pericarp thickness mainly exhibits a pattern of first increasing and then decreasing. By observing the changes in the structural characteristics of pericarp cells, it was found that the differences in pericarp thickness of individual kernels at different developmental stages were mainly determined by the thickness of individual pericarp cells and the number of mesocarp cells, and that pericarp development was divided into four stages: growth, formation, autophagy, and maturation. The pericarp thickness rapidly increases during the growth stage and reaches the maximum during the formation stage. Further, it begins to decrease during the autophagy stage and stabilizes during the maturation stage. Additionally, a high-density genetic map with a total length of 4,467.48 cM and 10 linkage groups was constructed using the BC4F4 population in combination with a maize 10 K SNP microarray, and the average genetic distance between adjacent markers was 1.57 cM. A total of 17 QTLs related to pericarp thickness were identified by combining phenotypic data and genetic mapping. The results revealed that qPT1-1 was the main QTL controlling UG and explaining 4.93%–12.97% of phenotypic variation, whereas qPT2-1 was the main QTL controlling UA and explaining 4.04%–8.63% of phenotypic variation. In addition, two candidate genes, Zm00001d001964 and Zm00001d002283, were screened for possible control of maize pericarp thickness, in combination with functional annotation, DNA sequencing and qPCR analysis. Among them, Zm00001d001964 has an Auxin_inducible domain, while Zm00001d002283 has a PKinase domain. The application of these candidate genes to the breeding process is expected to help improve maize pericarp thickness and provide valuable guidance and insights to breeders. However, further studies and experiments are needed to confirm the functions and effects of these genes.

Availability of data and materials

All data generated or analyzed during this study are included in this article and its supplementary information files or are available from the corresponding author on reasonable request.

References

  1. García-Lara S, Bergvinson DJ, Burt AJ, Ramputh AI, Díaz-Pontones DM, Arnason JT. The role of pericarp cell wall components in maize weevil resistance. Crop Sci. 2004;44(5):1546–52.

    Article  Google Scholar 

  2. Hu XW, Wang YR, Wu YP. Effects of the pericarp on imbibition, seed germination, and seedling establishment in seeds of Hedysarum scoparium Fisch. et Mey. Ecolog Res. 2009;24(3):559–64.

    Article  Google Scholar 

  3. Bailey D, Bailey R. The relation of the pericarp to tenderness in sweet corn. Proc Am Soc Hort Sci. 1938;36:555–9.

    Google Scholar 

  4. Richardson D. Pericarp thickness in popcorn. Agron J. 1960;52:77–80.

    Article  Google Scholar 

  5. Wolf M, Buzan C, MacMasters MM, Rist C. Structure of the mature corn kernel: III  microscopic structure of the endosperm of dent corn. Cereal Chem. 1952;29:349–61.

    CAS  Google Scholar 

  6. Helm J, Zuber M. Effect of harvest date on pericarp thickness in dent corn. Can J Plant Sci. 1970;50(4):411–3.

    Article  Google Scholar 

  7. Ito GM, Brewbaker JL. Genetic advance through mass selection for tenderness in sweet corn. J Am Soc Hortic Sci. 1981;106(4):496–9.

    Article  Google Scholar 

  8. Huelson WA. Sweet corn: Interscience. 1954.

    Google Scholar 

  9. Zhang J, Yue Y, Hu M, Yi F, Chen J, Lai J, Xin B. Dynamic transcriptome landscape of maize pericarp development. Plant J. 2024; 117(5):1574–91.

  10. Xiong C, Pei H, Zhang Y, Ren W, Ma Z, Tang Y, Huang J. Integrative analysis of transcriptome and miRNAome reveals molecular mechanisms regulating pericarp thickness in sweet corn during kernel development. Front Plant Sci. 2022;13:945379.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Yunian H. Measurement and analysis of pericarp thickness on different genotypes of sweet corn (Zea mays L.). Shanghai Nongye Xuebao. 1995;11(4):51–4.

    Google Scholar 

  12. Landoni M, Puglisi D, Cassani E, Brunoldi G, Comaschi C, Pilu R. Phlobaphenes modify pericarp thickness in maize and accumulation of the fumonisin mycotoxins. Sci Rep. 2020;10:1417.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zhu Z, Liang H, Chen G, Li F, Wang Y, Liao C, Hu Z. The bHLH transcription factor SlPRE2 regulates tomato fruit development and modulates plant response to gibberellin. Plant Cell Rep. 2019;38:1053–64.

    Article  CAS  PubMed  Google Scholar 

  14. Gan L, Song M, Wang X, Yang N, Li H, Liu X, Li Y. Cytokinins are involved in regulation of tomato pericarp thickness and fruit size. Horticulture Res. 2022;9:uhab041.

    Article  Google Scholar 

  15. Ito GM, Brewbaker JL. Genetic analysis of pericarp thickness in progenies of eight corn hybrids. J Am Soc Hortic Sci. 1991;116(6):1072–7.

    Article  Google Scholar 

  16. Choe E, Rocheford TR. Genetic and QTL analysis of pericarp thickness and ear architecture traits of Korean waxy corn germplasm. Euphytica. 2012;183:243–60.

    Article  Google Scholar 

  17. Helm J, Zuber M. Inheritance of pericarp thickness in corn belt maize. Crop Sci. 1972;12(4):428–30.

    Article  Google Scholar 

  18. Wang B, Brewbaker J. Quantitative trait loci affecting pericarp thickness of corn kernels. Maydica. 2001;46:159–65.

    Google Scholar 

  19. Park KJ, Sa KJ, Koh H-J, Lee JK. QTL analysis for eating quality-related traits in an F2: 3 population derived from waxy corn× sweet corn cross. Breed Sci. 2013;63(3):325–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ghasemi A, Zahediasl S. Normality tests for statistical analysis: a guide for non-statisticians. Int J Endocrinol Metab. 2012;10(2):486.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Obilor EI, Amadi EC. Test for significance of Pearson’s correlation coefficient. Int J Innov Mathematics Statistics Energy Policies. 2018;6(1):11–23.

    Google Scholar 

  22. Knapp S, Stroup W, Ross W. Exact confidence intervals for heritability on a progeny mean basis 1. Crop Sci. 1985;25(1):192–4.

    Article  Google Scholar 

  23. Irfan M, Ting ZT, Yang W, Chunyu Z, Qing M, Lijun Z, Feng L. Modification of CTAB protocol for maize genomic DNA extraction. Res J Biotechnol. 2013;8(1):41–5.

    CAS  Google Scholar 

  24. Meng L, Li H, Zhang L, Wang J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3(3):269–83.

    Article  Google Scholar 

  25. Ouellette LA, Reid RW, Blanchard SG, Brouwer CR. LinkageMapView—rendering high-resolution linkage and QTL maps. Bioinformatics. 2018;34(2):306–7.

    Article  CAS  PubMed  Google Scholar 

  26. McCouch SR. Report on QTL nomenclature. Rice Genet Newsl. 1997;14:11–3.

    Google Scholar 

  27. Portwood JL, Woodhouse MR, Cannon EK, Gardiner JM, Harper LC, Schaeffer ML, Walsh JR, Sen TZ, Cho KT, Schott DA. MaizeGDB 2018: the maize multi-genome genetics and genomics database. Nucleic Acids Res. 2019;47(D1):D1146–54.

    Article  PubMed  Google Scholar 

  28. Paysan-Lafosse T, Blum M, Chuguransky S, Grego T, Pinto BL, Salazar GA, Bileschi ML, Bork P, Bridge A, Colwell L. InterPro in 2022. Nucleic Acids Res. 2023;51(D1):D418–27.

    Article  CAS  PubMed  Google Scholar 

  29. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:1–11.

    Article  Google Scholar 

  30. Liu W, Xie Y, Ma J, Luo X, Nie P, Zuo Z, Lahrmann U, Zhao Q, Zheng Y, Zhao Y. IBS: an illustrator for the presentation and visualization of biological sequences. Bioinformatics. 2015;31(20):3359–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Lin F, Jiang L, Liu Y, Lv Y, Dai H, Zhao H. Genome-wide identification of housekeeping genes in maize. Plant Mol Biol. 2014;86:543–54.

    Article  CAS  PubMed  Google Scholar 

  32. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  33. Hagen G, Guilfoyle T. Auxin-responsive gene expression: genes, promoters and regulatory factors. Plant Mol Biol. 2002;49:373–85.

    Article  CAS  PubMed  Google Scholar 

  34. Inzé D, De Veylder L. Cell cycle regulation in plant development. Annu Rev Genet. 2006;40:77–105.

    Article  PubMed  Google Scholar 

  35. Lux A, Morita S, Abe J, Ito K. An improved method for clearing and staining free-hand sections and whole-mount samples. Ann Bot. 2005;96(6):989–96.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Wanlayaporn K, Somyong S, Pootakham W, Shearman J, Vanavichit A, Kumar P, Chee PW, Tragoonrung S. QTL mapping of pericarp thickness in immature and mature stages in thai tropical sweet corn (Zea mays var. saccharata). Chiang Mai J Sci. 2018;45:177–87.

    CAS  Google Scholar 

  37. Wolf M, Cull IM, Helm J, Zuber M. Measuring thickness of excised mature corn pericarp 1. Agron J. 1969;61(5):777–9.

    Article  Google Scholar 

  38. Hale TA, Hassell RL, Phillips T, Halpin E. Penetrometer and taste panel perception of pericarp tenderness in su, se, and sh2 sweet corn at three maturities. HortTechnology. 2004;14(4):521–4.

    Article  Google Scholar 

  39. Xiong F, Yu X, Zhou L, Wang F, Xiong A. Structural and physiological characterization during wheat pericarp development. Plant Cell Rep. 2013;32:1309–20.

    Article  CAS  PubMed  Google Scholar 

  40. Earp C, McDonough C, Rooney L. Microscopy of pericarp development in the caryopsis of Sorghum bicolor (L.) Moench. J Cereal Sci. 2004;39(1):21–7.

    Article  Google Scholar 

  41. Duffus C, Rosie R. Some enzyme activities associated with the chlorophyll containing layers of the immature barley pericarp. Planta. 1973;114:219–26.

    Article  CAS  PubMed  Google Scholar 

  42. Krishnan S, Dayanandan P. Structural and histochemical studies on grain-filling in the caryopsis of rice (Oryza sativa L.). J Biosciences. 2003;28:455–69.

    Article  CAS  Google Scholar 

  43. Thorne JH. Phloem unloading of C and N assimilates in developing seeds. Annu Rev Plant Physiol. 1985;36(1):317–43.

    Article  CAS  Google Scholar 

  44. Zhou Z, Wang L, Li J, Song X, Yang C. Study on programmed cell death and dynamic changes of starch accumulation in pericarp cells of Triticum aestivum L. Protoplasma. 2009;236:49–58.

    Article  CAS  PubMed  Google Scholar 

  45. Yu Y, Li G, Qi X, Li C, Mao J, Hu J. Mapping and epistatic interactions of QTLs for pericarp thickness in sweet corn. Acta Agron Sin. 2015;41(3):359–66.

    Article  CAS  Google Scholar 

  46. Wu X, Wang B, Xie F, Zhang L, Gong J, Zhu W, Li X, Feng F, Huang J. QTL mapping and transcriptome analysis identify candidate genes regulating pericarp thickness in sweet corn. BMC Plant Biol. 2020;20:1–13.

    Article  Google Scholar 

  47. Su L, Bassa C, Audran C, Mila I, Cheniclet C, Chevalier C, Bouzayen M, Roustan J-P, Chervin C. The auxin Sl-IAA17 transcriptional repressor controls fruit size via the regulation of endoreduplication-related cell expansion. Plant Cell Physiol. 2014;55(11):1969–76.

    Article  CAS  PubMed  Google Scholar 

  48. Spartz AK, Lee SH, Wenger JP, Gonzalez N, Itoh H, Inzé D, Peer WA, Murphy AS, Overvoorde PJ, Gray WM. The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J. 2012;70(6):978–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Stamm P, Kumar PP. Auxin and gibberellin responsive Arabidopsis SMALL AUXIN UP RNA36 regulates hypocotyl elongation in the light. Plant Cell Rep. 2013;32:759–69.

    Article  CAS  PubMed  Google Scholar 

  50. Kong Y, Zhu Y, Gao C, She W, Lin W, Chen Y, Han N, Bian H, Zhu M, Wang J. Tissue-specific expression of SMALL AUXIN UP RNA41 differentially regulates cell expansion and root meristem patterning in Arabidopsis. Plant Cell Physiol. 2013;54(4):609–21.

    Article  CAS  PubMed  Google Scholar 

  51. Chae K, Isaacs CG, Reeves PH, Maloney GS, Muday GK, Nagpal P, Reed JW. Arabidopsis SMALL AUXIN UP RNA63 promotes hypocotyl and stamen filament elongation. Plant J. 2012;71(4):684–97.

    Article  CAS  PubMed  Google Scholar 

  52. Li M, Chen R, Gu H, Cheng D, Guo X, Shi C, Li L, Xu G, Gu S, Wu Z. Grape small auxin upregulated RNA (SAUR) 041 is a candidate regulator of berry size in grape. Int J Mol Sci. 2021;22(21):11818.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Silva-Sanchez C, Li H, Chen S. Recent advances and challenges in plant phosphoproteomics. Proteomics. 2015;15(5–6):1127–41.

    Article  CAS  PubMed  Google Scholar 

  54. Mithoe SC, Menke FL. Phosphoproteomics perspective on plant signal transduction and tyrosine phosphorylation. Phytochemistry. 2011;72(10):997–1006.

    Article  CAS  PubMed  Google Scholar 

  55. van Wijk KJ, Friso G, Walther D, Schulze WX. Meta-analysis of Arabidopsis thaliana phospho-proteomics data reveals compartmentalization of phosphorylation motifs. Plant Cell. 2014;26(6):2367–89.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Czerednik A, Busscher M, Bielen BA, Wolters-Arts M, de Maagd RA, Angenent GC. Regulation of tomato fruit pericarp development by an interplay between CDKB and CDKA1 cell cycle genes. J Exp Bot. 2012;63(7):2605–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Zhao F, Zhang J, Weng L, Li M, Wang Q, Xiao H. Fruit size control by a zinc finger protein regulating pericarp cell size in tomato. Mol Hortic. 2021;1:1–16.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all the reviewers who participated in the review, as well as MJEditor (www.mjeditor.com) for providing English editing services during the preparation of this manuscript.

Funding

This research was funded by the Natural Science Foundation of Guangdong Province, China (2023A1515030066), the Natural Science Foundation of Hubei Province, China (2022CFB395), the Science and Technology Innovation Team Project of Hubei University of Arts and Science (No. 2021kptd01), Xiangyang Youth Science and Technology Talent Development Plan (2019).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, J.H.; methodology, G.G., H.P. and Y.T.; software, G.G and H.J.; investigation, G.G., Y.T. and H.P.; resources, J.H.; data curation, G.G. and H.P.; writing—original draft preparation, G.G and H.J.; writing—review and editing, L.Z. and J.H.; supervision J.H.; funding acquisition, L.Z. and J.H. All authors have read and agreed to the published version of the manuscript. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Lihong Zhai or Jun Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gong, G., Jia, H., Tang, Y. et al. Genetic analysis and QTL mapping for pericarp thickness in maize (Zea mays L.). BMC Plant Biol 24, 338 (2024). https://doi.org/10.1186/s12870-024-05052-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05052-1

Keywords