Fine Mapping of a New Major QTL-qGLS8 for Gray Leaf Spot Resistance in Maize

Gray leaf spot (GLS), caused by different species of Cercospora, is a fungal, non-soil-borne disease that causes serious reductions in maize yield worldwide. The identification of major quantitative trait loci (QTLs) for GLS resistance in maize is essential for developing marker-assisted selection strategies in maize breeding. Previous research found a significant difference (P < 0.01) in GLS resistance between T32 (highly resistant) and J51 (highly susceptible) genotypes of maize. Initial QTL analysis was conducted in an F2 : 3 population of 189 individuals utilizing genetic maps that were constructed using 181 simple sequence repeat (SSR) markers. One QTL (qGLS8) was detected, defined by the markers umc1130 and umc2354 in three environments. The qGLS8 QTL detected in the initial analysis was located in a 51.96-Mb genomic region of chromosome 8 and explained 7.89–14.71% of the phenotypic variation in GLS resistance in different environments. We also developed a near isogenic line (NIL) BC3F2 population with 1,468 individuals and a BC3F2-Micro population with 180 individuals for fine mapping. High-resolution genetic and physical maps were constructed using six newly developed SSRs. The QTL-qGLS8 was narrowed down to a 124-kb region flanked by the markers ym20 and ym51 and explained up to 17.46% of the phenotypic variation in GLS resistance. The QTL-qGLS8 contained seven candidate genes, such as an MYB-related transcription factor 24 and a C3H transcription factor 347), and long intergenic non-coding RNAs (lincRNAs). The present study aimed to provide a foundation for the identification of candidate genes for GLS resistance in maize.


INTRODUCTION
Gray leaf spot (GLS), caused by different species of Cercospora, is a serious foliar disease of maize (Zea mays L.) globally, especially in the Americas and Africa (Ward et al., 1999;Crous and Braun, 2003;Crous et al., 2006). GLS was first identified in Asia in 2007 in Nepal (Manandhar et al., 2011). It also occurs in the spring production regions of maize in northeast China (Li and Mei, 2008) and has spread to the north, east, and southwest of China over the past 20 years. It has surpassed the occurrence of maize big spots in the Yunnan, Sichuan, and Hubei provinces of China and has become the most important leaf disease (Zhou et al., 2012). GLS is prevalent in production areas characterized by dew formation in the morning, followed by hot, humid afternoons, and relatively cool nights (Shrestha et al., 2019). GLS has become increasingly important and is currently seen as one of the most serious yieldlimiting diseases of maize (Nutter Jr and Jenco, 1992;Ward and Nowell, 1998;Crous et al., 2006;Dhami et al., 2015). Previous studies have reported yield reductions of 260-320 kg/hm 2 due to outbreaks of GLS (Engelsing et al., 2012).
Maize GLS is caused by various species of Cercospora, such as Cercospora zeae-maydis (C. zeae-maydis), Cercospora zeina (C. zeina), and Cercospora sorghi (C. sorghi) var. maydis (Wang et al., 1998). According to Wang et al. (1998) and Crous et al. (2006), the variation is rather low in populations of C. zeae-maydis and C. zeina, which are genetically and phenotypically distinct. In China, there are two Cercospora species causing GLS that are C. zeae-maydis and C. zeina. The former is found in northern China, such as, Heilongjiang, Jilin, Liaoning, Inner Mongolia, and Shandong provinces, and the latter was distributed in southwest China, such as Yunnan and Hubei provinces (Liu et al., 2013). Pathogen spores spread through the air and infect maize leaves, leaf sheaths, and bracts under high-temperature and highhumidity conditions (Zhang et al., 2012). Breeding GLS resistant cultivars is a prominent strategy being used to control this disease (Shrestha et al., 2019).
GLS resistance in maize is controlled by additive genetic effects, which is characteristic of quantitative trait inheritance (Juliatti et al., 2009;Du et al., 2020). In recent years, a number of quantitative trait loci (QTLs) studies on GLS resistance in maize have been conducted. As a result, a major QTL for GLS resistance has been identified on 10 different chromosomes in maize (Bubeck et al., 1993;Clements et al., 2000;Lehmensiek et al., 2001;Gordon et al., 2004;Pozar et al., 2009;Zwonitzer et al., 2010;Zhang et al., 2012). Results on the identification and stability of the QTLs, however, vary widely in different genetic backgrounds and environments. Shi et al. (2007) identified seven hotspot QTL and candidate genes on chromosomes 2, 3, and 8 of maize based on QTL meta-analysis. Wang et al. (2014) identified 11 consistent QTL intervals and three candidate genes for GLS resistance in maize through meta-analysis methods involving the homologous alignment of the maize genome with rice and Arabidopsis thaliana genomes.
Previous fine-mapping studies of QTL for GLS resistance have narrowed several QTL regions (Zhang et al., 2012(Zhang et al., , 2017Xu et al., 2014;Benson et al., 2015;Lv et al., 2020;Sun et al., 2021). Zhang et al. (2012) identified two major QTLs, qRgls1 and qRgls2, for GLS resistance on chromosomes 8 and 5, respectively, and fine mapped qRgls1 to an interval of 1.4 Mb. In a subsequent study, Zhong (2018) narrowed the region of qRgls1 to 120 kb Abbreviations: GLS, Gray leaf spot; QTLs, Quantitative trait loci; SSR, Simple sequence repeat; NIL, Near isogenic line; MAS, Marker-assisted selection; JD, Jiangdong; MG, Mengga; MW, Mengwen; cM, Centimorgans; QTL-qGLS8, An anti-disease position on chromosome 8 in maize. on chromosome 8 in a 9 MB physical location of the B73 genome. This region was cloned by screening a bacterial artificial chromosome (BAC) library of the resistant parent. QTL-qRgls2 was narrowed from an initial ∼110-Mb region to an interval of ∼1 Mb on chromosome 5, which harbored 15 predicted genes (Xu et al., 2014). Another QTL was finely mapped to an ∼130-Kb region on chromosome 8 of a 152-Mb region of the B73 genome (Zhang et al., 2017). A previous study also narrowed QTL for GLS resistance to 5.2 and 6.5 Mb on chromosome 1 (Benson et al., 2015). Although the above studies have greatly advanced the identification of genetic markers for GLS resistance in maize, the discrepancies among different genetic backgrounds and environments underscore the need for further studies.
Molecular marker-assisted selection (MAS) represents an important strategy for breeding disease-resistant varieties. In the present study, QTL analyses and fine mapping in a Jiao 51 (J51, an elite inbred line) (Yang and Yang, 2007) near isogenic line (NIL) background was conducted and identified a quantitative trait locus for resistance to maize gray leaf spot (QTL-qGLS8) derived from T32 (an inbred line resistant to GLS derived from "Suwan") (Wu et al., 2019). The present study provides new information on genetic markers for GLS resistance in maize that can be utilized for MAS breeding and provides a foundation for studying the functional role of candidate genes associated with GLS resistance in maize.

Plant Material and Field Experiments
The initial mapping population comprised 189 F 2 : 3 families that originated from a cross between inbred T32 (a highly resistant genotype) and J51 (a highly susceptible genotype). F 2 individuals were developed by planting F 1 seeds in the winter of 2014 in Hainan, the southernmost province of China. Next, an F 2 : 3 segregation population developed from selfing each F 2 individual was planted in the spring of 2015 in Jiangdong (JD), Mengga (MG), and Mengwen (MW), in Yunnan, a southwest province of China. The F 2 and the F 2 : 3 families were used to conduct a QTL mapping analysis, which resulted in identifying an anti-disease position (named QTL-qGLS8) on chromosome 8 of maize (Yang, 2016).
Fine QTL mapping was performed by developing 20 BC 3 F 2 lines, in which the recombination occurred within the region containing the target QTL. The original cross T32 (donor parent) × J51 (recurrent parent) was made in Hainan during the winter of 2013. Randomly selected F 1 ears of the cross were pollinated by the recurrent parent J51 in the spring of 2014 to develop a backcross population (BC 1 F 1 ). BC 1 F 1 individuals were continuously backcrossed with the recurrent parent J51 using MAS until the generation of the BC 3 F 1 population in the spring of 2015. Twenty individuals were selected from the BC 3 F 1 population using simple sequence repeat (SSR) markers, denoted as D27, D128, D187, D221, D352, D355, D405, D428, D441, D493, D544, D550, D552, D600, D640, D647, D661, D707, D752, and D1123. The selected individuals formed a NIL   (Figure 1).

GLS Disease Evaluation
A total of 189 F 2 : 3 families were evaluated for GLS resistance under natural conditions in two replicates at each of three locations, JD, MG, and MW. A scoring system with five levels (1, 3, 5, 7, and 9) was used to rank GLS resistance. A completely randomized block design with two replications and 20 plants in each plot was established in each of the three sites. All plants were self-pollinated, and 13 plants from the middle of each plot were sampled at 2-and 4-week post-pollination and ranked for GLS resistance ( Table 1) (Zhang et al., 2012). The phenotypic evaluation of the BC 3 F 2 population was the same as was used for the F 2 : 3 population.

Statistical Analysis of the Phenotypic Data
Statistical analysis of the phenotypic data was performed using IBM SPSS Statistical 20.0 software. The variances for genotype, environment, and interaction between genotype and environment were estimated using Microsoft Office Excel based on a random model. The broad-sense heritability (H 2 ) was calculated as follows: H 2 = δ 2 G /(δ 2 G + δ 2 GE n + δ 2 E /nr), where δ 2 G represents the genetic variance, δ 2 GE represents the genotype and environment interaction, δ 2 E represents the error variance,  and n and r are the numbers environments and replications, respectively (Knapp et al., 1985).

Linkage Map Construction for the F 2 Population
Genomic DNA was extracted from samples of young leaves of each plant using the cationic detergent cetyl-trimethyl ammonium bromide (CTAB) method (Chen and Ronald, 1999). The genotype of an individual was identified using a 10% nondenaturing polyacrylamide gel electrophoresis. A total of 286 polymorphic markers between T32 and J51 were used to develop the genetic map in JoinMap version 4 (Jacobs et al., 1995) based on publicly available SSR markers retrieved from the Maize Genetics and Genomics Database (http://www.maizegdb. org) (Portwood et al., 2019). The linkage map was constructed from 181 markers across 10 chromosomes. The recombination frequency between linked loci was transformed into the genetic distance (centimorgans, cM) through Kosambi's function.

SSR Marker Development
For the QTL mapping of the BC 3 F 2 -micro population, QTL-qGLS8 was detected in a CI between SSR markers CSU329 and umc1034 with a physical distance of 1,102 kb according to the maize B73 whole-genome physical map (http://www. maizesequence.org/). Bacterial artificial chromosome sequences within the qGLS8 CI were downloaded from the website. These sequences were used to design new markers between the T32 and the J51 parent lines. The sequences were first scanned using SSRHunter1.3 software (Li and Wan, 2005) to identify SSRs. Primers were then designed using Primer Premier 5.0 software with the following criteria: primers should be ∼20 nucleotides in length with a 40-60% guanine-cytosine content, there should be no consecutive tracts of a single nucleotide and no secondary structure. The primer sequences, mapping position, and the amplified length of the eight new markers developed in this study are listed in Table 2.

QTL Analysis
Inclusive composite interval mapping was performed using QTL IciMapping software  to identify the QTL and estimate their effect. Parameters for forward regression analysis were set with window size and walking speed of 10 and 1 cM, respectively. The significance threshold for QTL detection included 1,000 random permutations of the phenotypic data at the 5% level. The gene action mode of each significant QTL was estimated in accordance with the rate ( (Edwards et al., 1987;Tuberosa et al., 1998). QTLs were named using a standardized approach, for example, in JDGLS4a, "JD" represents the environment (site) in which the QTL was identified (JD, MG, and MW are abbreviations for the Jiangdong, Mengga, and Mengwen planting sites, respectively); "GLS" represents the name of the trait; the number "4" represents the serial number of the chromosome on which the QTL is located; and "a" represents the serial alpha code of the detected QTL. A total of 1,468 BC 3 F 2 plants were genotyped using 10 SSR markers to construct a small-scale linkage map. Eight of the SSR markers were newly developed. The small-scale genetic mapping of the qGLS8 genes was performed using JoinMap version 4 software to combine the genotypic data of molecular markers with QTL-qGLS8 genes.

Phenotypic Variation of GLS Resistance in the F 2 : 3 Families
The t-tests of the differences in GLS resistance in the two parents within a site were highly significant (P < 0.01; Table 3). GLS resistance in the T32 parent was significantly higher than it was in the J51 parent. The scores for GLS resistance exhibited continuous and approximately normal distributions across the three sites with low skewness and kurtosis, indicating that GLS resistance is a typical quantitative trait (Table 3; Supplementary Figure 1). A transgressive segregation of the resistance score was observed in the populations across the three sites. In addition, variance analysis indicated that the effect of genotype (F 2 : 3 families), environment, and a genotype/environment interaction on GLS resistance was highly significant (P < 0.01), and that heritability was >80% (Table 4).

Linkage Map of the F 2 Population
The length of the QTL map constructed using the 181 polymorphic SSR markers (Figure 2) was 1153.497 cM, and the average distance between markers was 6.37 cM. The average number of markers per chromosome was 18.1, ranging from 16 on chromosome 2 to 22 on chromosome 6. A few differences in marker order on a chromosome were compared to the IBM 2008 Neighbors Frame 6 (Supplementary Figure 2). These variations in distance and order may be due to the different materials used in the population.

Single Environment Analysis of the GLS Resistance Main-Effect QTL in the F 2 : 3 Families
Fourteen QTLs for GLS resistance were detected in F 2 : 3 families, which comprised five QTLs at JD, four QTLs at MG, and five QTLs at MW ( Table 5). Two of the loci were detected in different environments. JDGLS4a, MGGLS4a, and MWGLS4a were located in a common genetic region on chromosome 4 within the Phi076-bnlg1890 marker interval (bin4.11) and explained 5.86, 6.75, and 6.20% of the phenotypic variation in GLS resistance at the JD, MG, and MW sites, respectively. Additive effects of three QTLs on increased GLS resistance were contributed by T32, whose effect values were −0.34, −0.32, and −0.36, respectively. JDGLS8a, MGGLS8a, and MWGLS8a were also located in a common genetic region on chromosome 8 in the umc1130-umc2354 marker interval (bin8.03) and explained 7.89, 14.71, and 11.51% of the phenotypic variation in GLS resistance at the JD, MG, and MW sites, respectively. Additive effects of three QTL on increased GLS resistance were also contributed by T32, whose effect values were −0.22, −0.61, and −0.46 respectively. Compared to qGLS8, the effect value of qGLS4 was lower. Therefore, qGLS8 was used as the main locus for further study. The additive effects of the T32 allele in 86% of the QTL increased the GLS resistance level by 0.22-0.61, compared to the additive effect of the J51 allele. The modes of gene action were A, PD, and OD.

Joint Analysis of QTL Detection in the F 2 : 3 Families Across Three Environments
Twelve QTLs were detected for GLS resistance in the joint analysis combining data from all three sites ( Table 6). The two common genetic regions detected in the single environment (site) analysis were also detected in the joint analysis within the same marker intervals with the relatively higher log of odd (LOD) values of 8.25 and 10.11. The two QTLs detected in only one environment and the one QTL detected only in two environments (sites) in the single environment analysis were also detected in the joint analysis within the same marker intervals. One common genetic regions detected in the MG and MW sites in the single environment analysis was also detected in the joint analysis, located near the genetic region of QTL-JDGLS7a, which was detected only at the JD site with the same marker (umc1944). Some QTLs detected in the joint analysis across all three environments (sites) were not detected in the single environment analysis. These QTLs probably have a minor effect and are detected when the analysis includes a greater number of environments.

Fine Mapping of qGLS8
Randomly selected F 1 ears of the original T32 (donor parent) × J51 (recurrent parent) cross were pollinated by the recurrent parent J51 to develop a backcross population (BC 1 F 1 ). BC 1 F 1 individuals were continuously backcrossed with the recurrent parent J51 using MAS until the BC 3 F 1 generation was obtained. Twenty individuals were selected from the BC 3 F 1 population using SSR markers (designated as umc1130 and umc2354). The BC 3 F 2 population, developed from selfing each BC 3 F 1 individual, was used to identify recombinants for fine mapping.
In 20 individuals, the T32 chromosomal segment harboring qGLS8 was present in the genetic background of J51. A total of 180 BC 3 F 2 -Micro individuals were genotyped using 12 SSR markers located in the target region and evaluated for GLS resistance. Fine mapping confirmed a QTL peak between markers CSU329 and umc1034 with an LOD score of 4.55 and an R 2 of 11.11% (Table 7). This peak differed from the results of the original F 2 : 3 mapping (Tables 5, 6), which may have been due to the increased SSR marker density used in the fine mapping and the elimination of interference from the genetic background. The phenotypic means of BC 3 F 2 -Micro individuals were compared among the three genotypic classes defined by the allele constitution of markers in the CSU329-umc1034 interval to confirm the gene action of the QTL-qGLS8 identified in this study. The difference in the level of GLS resistance between plants having a homozygous T32 introgression in the target region and those having a homozygous introgression of J51 DNA was 0.61 (a difference of 8.24%). These results confirmed the  identification of the primary QTL in the F 2 : 3 population and demonstrated that the transgressive segregant could be identified in subsequent generations. The mean value of GLS resistance in the homozygous T32 class (7.40) was not significantly different from that of the heterozygous class (7.37); however, both were significantly different from that of the homozygous J51 class (8.01, P < 0.01). These results indicate that the T32 allele was completely dominant over the J51 allele. This analysis implied that the homozygous T32 class and the heterozygotes could be combined for use in the phenotypic analysis during fine mapping, which provided a larger sample size and greater statistical power. Four newly developed polymorphic SSR markers were added to the existing SSR markers to genotype BC 3 F 2 -Micro individuals to precisely define the location of the critical breakpoints. A high-resolution genetic linkage map was constructed in the target region. QTL analysis confirmed the highly significant peak between markers ym20 and ym31 with an LOD score of 4.34 and an R 2 value of 10.51% (Table 7).
In addition to the QTL analysis, an extreme sampling strategy was used to clarify on which side of ym20-ym31 the qGLS8 locus was located. A total of 180 BC 3 F 2 -Micro individuals were compared genotypically using six of the markers in the target region. Based on these data, 35 recombinants were identified ( Figure 3A). The mean phenotypic performance of individuals with and without the T32 introgression was compared, and results indicated that groups B and C were significantly different (P < 0.05). Thus, the ym20-ym31 interval represents a critical region of recombination, with the critical breakpoint falling  between ym20 and ym31. This finding further supported the results of the QTL linkage analysis. A total of 1,468 BC 3 F 2 plants were genotyped for the fine mapping of qGLS8 and included the use of four new markers developed to further subdivide the target interval. The QTL analysis confirmed the presence of an extremely significant peak between markers ym20 and ym51 with LOD and R 2 values of 58.61 and 17.46%, respectively ( Table 7). A total of 28 recombinants were identified and grouped into six genotypic classes consisting of 5, 6, 3, 2, 4, and 8 individuals. When the phenotypic means for each class were compared with the T32 and J51 parents (Figure 3B), the most informative comparison was between classes B and C, where the phenotypes (between T32 and J51) were significantly different (P < 0.05). Based on these results, the location of qGLS8 was concluded to be in the recombination area between ym20 and ym51. This conclusion was further supported by the comparison between classes B and E. Consequently, the location of qGLS8 was narrowed down to a 124 kb fragment defined by ym20 and ym51 markers on chromosome 8.
The eight newly developed markers were authentically mapped on bin 8.03 with two of the markers (ym20 and ym51) within and six of the markers located outside of the resistance qGLS8 region. The newly developed markers will greatly facilitate both map-based cloning and MAS breeding for GLS resistance in maize (Figure 4). The 124-kb region harbors seven genes based on the B73 Reference-Gramene_v4 (Tello-Ruiz et al., 2018) and the Maize GDB, including five genes and two lincRNAs ( Table 8).

DISCUSSION
Maize, one of the primary cereal crops produced for both human and animal consumption, suffers from numerous biotic and abiotic stresses, which may result in considerable reductions in yield (Dhami et al., 2015). GLS is a major threat to maize production worldwide (Lv et al., 2020). GLS resistance in maize is a typical quantitative trait (Zhang et al., 2017) and the identification of major QTL and candidate genes is crucial to the use of marker-assisted breeding for GLS resistance.
Previous studies have reported more than 100 QTLs for GLS resistance distributed on all 10 chromosomes of maize (Du et al., 2020). The QTL on chromosome 8, QTL8, is defined by markers npi590A-umc93 (4.7 cM) and umc48-umc30 (6.9 cM) in F 2 and F 2 : 3 populations and explains 7.7-11.0% of the phenotypic variance in GLS resistance (Maroof et al., 1996). Shi et al. (2007) detected a consensus QTL (446.14 cM) in Bin 8.06 using a meta-analysis of five previous studies. In addition, two of the loci of QTL8 are defined by FIGURE 4 | The application of SSR marker ym51 in genotyping the BC 5 F 1 population by electrophoresis on a PAGE gel. Polymorphic bands were observed in the two parental lines, "T32" (lane 1), "J51" (lane 2), and both representative bands in the "F1" (lane 3). Variations on the banding pattern were observed among the BC 5 F 1 individuals (lanes 4-33). Lanes 4-33 are randomly selected BC 5 F 1 individuals that are either homozygous (one band having the same size as that of "J51") or heterozygous (two bands corresponding to "T32" and "J51," respectively). SSR, simple sequence repeat. markers PZA01470.1 (Bin 8.03) and PZA03651.1 (Bin 8.06) using nested association mapping, and six candidate genes have been identified encoding the Rust resistance protein rp3-1, Heat shock factor-binding protein 1, GapC2, T cytoplasm male sterility restorer factor 2, Chloroplast phytoene synthase, and B transcriptional activator (Benson et al., 2015). The fine mapping of qGLS8 has indicated that qGLS8 is defined by the markers ctg358-32 and ctg358-01, located in an ∼130-kb region in Bin 8.06, and contains five genes encoding the Receptorlike protein kinase 2, Phosphoethanolamine N-methyltransferase 3, Receptor-like protein kinase 2, ABC transporter, and an uncharacterized protein (Zhang et al., 2017). In the present study, two populations of maize were developed to identify major QTL for GLS resistance in the F 2 : 3 population derived from the T32 parent and to conduct fine mapping in a BC 3 F 2 population with a J51 NIL background. The T32 and J51 parents may provide novel, major resistance genes. In our study, qGLS8 was detected as the only overlapping QTL, defined by the markers umc1130 and umc2354 (Bin 8.03), in different environments (sites) and explains up to 12.98% of the phenotypic variation for GLS resistance in the F 2 : 3 population. In the BC 3 F 2 population, qGLS8, defined by the ym20 and ym51 (Bin 8.03) markers, explains up to 17.46% of the phenotypic variation in GLS resistance. The genetic effect of qGLS8 is significant and stable in the F 2 : 3 and BC 3 F 2 populations across multiple environments, suggesting that it is a major QTL for GLS resistance in maize.
In the present study, QTL-qGLS8 was narrowed down to a 124-kb fragment defined by markers ym20 and ym51 with a physical location of 21.23 Mb and 21.36 Mbp on chromosome 8, respectively, according to the B73 reference genome. The QTL detected in our study is different from the QTL previously reported by Zhang et al. (2012), with a physical location on chromosome 8 from 8.61 Mb to 10.07 Mbp, according to the B73 reference genome. Moreover, the QTL detected in the physical location from 23.77 Mbp to 101.18 Mbp by Benson et al. (2015) does not harbor the QTL-qGLS8 identified in the present study. Notably, the interval of the QTL identified in the present study does not overlap the map position for GSL resistance reported in other studies, indicating that it represents the discovery of a new major locus for GLS resistance, potentially resulting from the introgression of T32. The major QTL-qGLS8 in bin 8.03 represents a good choice for resistance gene cloning and for use in marker-assisted selection in breeding programs focusing on GSL resistance.
Seven candidate genes were identified in the QTL-qGLS8 region, including two uncharacterized lincRNAs and five genes, according to the B73 Reference-Gramene_v4. Zm00001d008808, encoding the MYB-related transcription factor 24 is significantly upregulated in response to drought (Zenda et al., 2019) and longday conditions (Ku et al., 2016), suggesting that it may be a stressresponsive gene. The identified MYB-related transcription factor has also been shown to function in regulating the circadian clock in Arabidopsis (Nguyen and Lee, 2016) and promoting flowering in rice (Zhang et al., 2016), which corresponds to previous research that found a strong correlation between flowering time and GLS resistance (Zwonitzer et al., 2010). Zm00001d008812, encoding the C3H-transcription factor 347, has been reported to be involved in cell wall biogenesis as a differentially expressed gene in elongating vs. non-elongating maize internodes (Bosch et al., 2011). C3H-transcription factors also contribute to salt stress tolerance in many plants (Jiang et al., 2017;Liu et al., 2020), which may reflect its role in the adaptation to stress response. The potential function of the other candidate genes, however, is unclear. The function of candidate lincRNAs in maize is also unclear; however, lincRNAs have been shown to play an important role in plant development and response to stress in Arabidopsis (Li et al., 2016). The function and the regulatory mechanisms of candidate genes and lincRNAs in GLS resistance require further study. Our future studies will focus on the functional role of the candidate genes Zm00001d008808 and Zm00001d008812 and involve gene cloning, expression analysis, overexpression studies, and other molecular approaches. The fine mapping of the newly identified QTL QTL-qGSL8 provides further insight into GLS resistance in maize and could play an important role in developing maize varieties with increased resistance to GLS.

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 author/s.