Genetic dissection of quantitative trait loci for flag leaf size in bread wheat (Triticum aestivum L.)

Flag leaf size is a crucial trait influencing plant architecture and yield potential in wheat. A recombinant inbred line (RIL) population derived from the cross of W7268 and Chuanyu 12 was employed to identify quantitative trait loci (QTL) controlling flag leaf length (FLL), flag leaf width (FLW), and flag leaf area (FLA) in six environments and the best linear unbiased estimator (BLUE) datasets. Using a 55 K SNP-based genetic map, six major and stable QTL were detected with 6.33–53.12% of explained phenotypic variation. Except for QFlw.cib-4B.3, the other five major QTL were co-located within two intervals on chromosomes 2B and 2D, namely QFll/Fla.cib-2B and QFll/Flw/Fla.cib-2D, respectively. Their interactions and effects on the corresponding traits and yield-related traits were also assessed based on flanking markers. QFll/Fla.cib-2B showed pleiotropic effects on spikelet number per spike (SNS). QFlw.cib-4B.3 and QFll/Flw/Fla.cib-2D had effects on grain number per spike (GNS) and thousand-grain weight (TGW). Comparison analysis suggested that QFll/Fla.cib-2B was likely a new locus. Two candidate genes, TraesCS2B03G0222800 and TraesCS2B03G0230000, associated with leaf development within the interval of QFll/Fla.cib-2B were identified based on expression-pattern analysis, gene annotation, ortholog analysis, and sequence variation. The major QTL and markers reported here provide valuable information for understanding the genetic mechanism underlying flag leaf size as well as breeding utilization in wheat.

Flag leaf size is a crucial trait influencing plant architecture and yield potential in wheat. A recombinant inbred line (RIL) population derived from the cross of W7268 and Chuanyu 12 was employed to identify quantitative trait loci (QTL) controlling flag leaf length (FLL), flag leaf width (FLW), and flag leaf area (FLA) in six environments and the best linear unbiased estimator (BLUE) datasets. Using a 55 K SNP-based genetic map, six major and stable QTL were detected with 6.33-53.12% of explained phenotypic variation. Except for QFlw.cib-4B.3, the other five major QTL were co-located within two intervals on chromosomes 2B and 2D, namely QFll/Fla.cib-2B and QFll/Flw/Fla.cib-2D, respectively. Their interactions and effects on the corresponding traits and yield-related traits were also assessed based on flanking markers. QFll/Fla.cib-2B showed pleiotropic effects on spikelet number per spike (SNS). QFlw.cib-4B.3 and QFll/ Flw/Fla.cib-2D had effects on grain number per spike (GNS) and thousand-grain weight (TGW). Comparison analysis suggested that QFll/Fla.cib-2B was likely a new locus. Two candidate genes, TraesCS2B03G0222800 and TraesCS2B03G0230000, associated with leaf development within the interval of QFll/Fla.cib-2B were identified based on expression-pattern analysis, gene annotation, ortholog analysis, and sequence variation. The major QTL and markers reported here provide valuable information for understanding the genetic mechanism underlying flag leaf size as well as breeding utilization in wheat.

Introduction
Bread wheat (Triticum aestivum L.), one of the most widely adapted food crops, provides about a quarter of the calories consumed by humans (Curtis and Halford, 2014). The shape, size, and posture of the leaves together with ear and awn, play a decisive role in the photosynthetic capacity of plants and also regulate many important agronomic traits, such as yield and biotic and abiotic stress responses (Sourdille et al., 2002;Liu et al., 2015).
Flag leaf size, estimated by flag leaf length (FLL), width (FLW), and area (FLA), is an important control of plant structure and is correlated with yield-related traits Niu et al., 2022;Zhao et al., 2022). In cereals, flag leaves are the main organ for photosynthesis and plays a crucial role in grain development, such as enhanced proteostasis, lipid remodeling, and nitrogen remobilization (Cohen et al., 2022). Therefore, breeding wheat with optimal leaf morphology has been regarded as an effective method to improve grain yield.
Previously studies attempting to uncover the genetic mechanism of flag leaf morphology in crops showed that flag leaf size was determined by quantitative trait loci (QTL) and significantly influenced by the environment (Coleman et al., 2001;Kobayashi et al., 2003). In rice, the genes that control flag leaf size have been extensively investigated and several major types of signaling pathways have been identified, including the transcription factor signaling, cell expansion pathway, microRNA regulator, and others Tian et al., 2015;Alamin et al., 2017;He et al., 2018;Zhou et al., 2020;Uzair et al., 2021;Huang et al., 2022). In rice, for example, OsCKX3 mediated the accumulation of cytokinins, the mutant osckx3 exhibited a larger flag leaf size. (Huang et al., 2022). Mapbased cloning in rice identified a GATA family transcription factor was a candidate gene for SNFL1 and the mutant snfl1 exhibited a reduction in flag leaf epidermal cell length (He et al., 2018). The MIR319 gene family contains two members, Osa-MIR319a and Osa-MIR319b. Overexpression of MIR319 in rice and MiR319/TaGAMYB3 module in wheat regulate the number of longitudinal small veins in the leaf, which led to an increase in leaf blade width, and improves grain yield. (Yang et al., 2013;Wang et al., 2014;Jian et al., 2022). NRL1 encodes the cellulose synthase-like protein D4 (OsCslD4) and plays a critical role in leaf morphogenesis by regulating longitudinal veins and adaxial bulliform cells development (Hu et al., 2010).
The molecular cloning of genes related to flag leaf size in wheat falls behind rice owing to its huge genome. So far, QTL for flag leaf size have been detected on almost all wheat chromosomes based on various genetic populations and environments (Fan et al., 2015;Wu et al., 2016;Liu et al., 2018;Zhao et al., 2018;Jin et al., 2020;Ma et al., 2020;Tu et al., 2021). For example, TaFLW1, a major QTL for FLW, was fine-mapped into a 0.2-cM interval on chromosome 5A, which is tightly linked to Fhb5 (Xue et al., 2013). QFlw-6A was fine-mapped to a small interval on chromosome 6A and 10 genes were predicted in this region (Yan et al., 2020). More recently, QFlw-5B was narrowed to a 2.5 Mb region and contained 27 predicted genes . To date, there is no report of map-based cloning of the gene controlling flag leaf size in wheat.
The present study was undertaken to (i) evaluate the performance of flag leaf size in a recombinant inbred line (RIL) population in multiple environments; (ii) identify QTL for flag leaf size using a wheat 55 K SNP-based genetic map and analyze their effects; (iii) assess relationships between flag leaf size and yield-related traits; (iv) predict candidate genes for major QTL.

Plant materials
A RIL population (WC12, 180 F 9 lines) derived from the cross of W7268 and Chuanyu 12 (CY12) by the single-seed descend method was used for gene mapping in this study. The wheat line W7268 was selected by our lab. It was characterized by desirable agronomic traits including high SNS, GNS, and flag leaf size. Thus, it has been widely used in wheat breeding programs and several elite varieties have been selected during the past decade. However, the genetic control of flag leaf size in W7268 is uncovered. CY12 is a commercial cultivar with a smaller flag leaf size. In addition, 135 wheat accessions (including 60 modern cultivars and 75 landraces) were used to genotyping. The population was constructed and the accessions were conserved by our laboratory.

Phenotyping and statistical analysis
RILs of WC12 were evaluated at two experimental sites in three growing seasons: 2018-2019 in Shuangliu (103°52′ E, 30°3 4′ N) (E1); 2018-2019 in Shifang (104°11′ E, 31°6′ N) (E2); 2019-2020 in Shuangliu (E3); 2019-2020 in Shifang (E4); 2020-2021 in Shuangliu (E5); 2020-2021 in Shifang (E6). A completely randomized block design was used for all of the trials in each environment. Each line was planted in a one-row plot with a row length of 1.2 m, a row spacing of 0.2 m, and 11 seeds per row. Two replicates were employed in each environment. Fertilizer (N: 25%, P 2 O 5 : 10%, K 2 O: 10%) was applied at sowing time at a rate of 450 kg/ha. Field management and disease control were performed in accordance with conventional practices in wheat production. After anthesis, the main tillers of ten representative plants from each line were selected for measuring flag leaf length (FLL), flag leaf width (FLW), and flag leaf area (FLA). FLL was measured as the distance from the base to the tip of the leaf; FLW was the width of the widest part of the leaf; FLA was derived from the FLL and FLW and estimated as FLL × FLW × 0.75 (Zeuli and Qualset, 1990). Also, the phenotypic values of some yield-related traits were measured. At maturity, ten representative plants were randomly selected to measure agronomic traits, including spikelet number per spike (SNS), spike compactness (SC), grain number per spike (GNS), and fertile tiller number (FTN). Spikelet number per spike (SNS) was determined by counting the number of spikelets in main spikes; SC was calculated by dividing the spike length (SL) by SNS. The main spikes of target plants were then harvested and threshed manually. GNS was then counted manually, and the thousandgrain weight (TGW) and grain yield per plant (GYP) was assessed with SC-G software (Wanshen Detection Technology Co., Ltd., Hangzhou, China).
Basic statistical analysis, frequency distribution analysis, and correlation coefficients analysis among traits were conducted on the phenotypic data using software SPSS25 (Chicago, IL, USA), R 4.1.2, and QTL IciMapping v4.2 (Meng et al., 2015). The best linear unbiased estimator (BLUE) was calculated using the R package "lme4" and used for combined QTL detection, correlation analysis, and effect analysis. Estimation of the broad-sense heritability (H 2 ) of each trait according to the method described by Smith (Smith et al., 1998). The significance of the difference was measured by the Student's t test (P < 0.05), Welch's two-sample t test (P < 0.05), and Wilcoxon's symbol rank-sum (P < 0.05) using SPSS25 and R 4.1.2, respectively.

Linkage map construction and QTL detection
A high-density genetic map containing 2,186 bin markers was constructed using the wheat 55 K SNP array according to a previous study (Liao et al., 2022). This genetic map spanned 2,398.67 cM across 21 chromosomes with an average interval of 1.10 cM/marker. There are 1,598 SNP markers on the map, of which genome A contains 428 markers, genome D contains 467 markers, and genome B contains the most markers (703). In addition, markers derived from three reported genes Ppd-B1, Rht-B1, and Ppd-D1 were integrated into the genetic map for QTL detection (Ellis et al., 2002;Guo et al., 2010;Dıáz et al., 2012). They were added to the corresponding chr2B, chr4B, and chr2D genetic map by JoinMap 4.1 (Voorrips et al., 2006). QTL analysis was conducted using the inclusive composite interval mapping (ICIM) function of QTL IciMapping v4.2 (Meng et al., 2015). Individual environmental QTL values were measured using the bi-parental populations (BIP) module with walking step = 1.00 cM, PIN = 0.001, and LOD score values ≥ 3. QTL with PVE value greater than 20% in at least one environment and could be stably detected in more than four environments (including the BLUE dataset) were considered as the major ones. QTL with common flanking markers or less than 1 cM apart were considered identical. QTL were named according to the International Rules of Genetic Nomenclature (https://wheat.pw.usda.gov/ggpages/wgc/98/Intro.htm). The "cib" represents "Chengdu Institute of Biology".

Marker development and genotyping
To develop new markers within the mapped interval, variants between parental lines were detected using exome sequencing. Exon capture, sequencing, and analysis were performed by Bioacme Biotechnology Co., Ltd (Wuhan, China, http://www.whbioacme.com).
On the basis of the preliminary QTL mapping result, new Kompetitive Allele Specific PCR (KASP) markers were developed within the interval of major QTL. According to the exome sequencing result, the variant sites (including SNPs and InDels) between W7268 and CY12 were screened. Only SNPs with differences between W7268 and CY12 at the major QTL, by referring to the Chinese Spring (CS) reference genome IWGSC RefSeq v2.1 (http://www.wheatgenome.org), were converted into KASP markers. KASP assay reaction procedure and data analysis implemented according to the method depicted by Ji (Ji et al., 2020).

Conditional QTL analysis for FLA
Conditional QTL analysis is an excellent tool for interpreting the relationship between complex traits and components at the QTL level (Cui et al., 2011). We carried out conditional QTL analysis as described by Liu (Liu et al., 2018). QTL IciMapping v4.2 (Meng et al., 2015) was used to identify the conditional loci with conditional phenotypic values. The conditional phenotypic values (T1|T2), which means the value of Trait 1 conditional on Trait 2, were obtained using QGAStation2.0 (Chen et al., 2012). Here, 'FLA|FLL' and 'FLA|FLW' refer to the value of FLA excluding the influences of FLL and FLW, respectively. In the analysis of QTL mapping results, if a QTL was detected only in unconditional QTL analysis, the locus was considered to have a large contribution to the corresponding trait; whereas a QTL was considered to be unassociated with the corresponding trait if it was detected in both conditional and unconditional QTL analysis.

Prediction of candidate gene
Physical intervals of the major QTL detected in this study were obtained by blasting against (E-value of 1e-5) their flanking markers sequences to genome sequences of CS reference genome IWGSC RefSeq v2.1 (http://www.wheatgenome.org). The annotations and functions of genes between flanking markers were further analyzed using UniProt (https://www.uniprot.org/). The expression pattern analysis of candidate genes was performed using Gene Expression of Triticeae Multi-omics Center (http: //202.194.139.32/expression/wheat.html) (Ma et al., 2021) and Wheat Expression Browser (http://www. wheat-expression.com) (Ramıŕez-Gonzaĺez et al., 2018). The circle graph of expression values was drawn using the ggplot2 package in R 4.1.2. Analysis of orthologous genes in wheat and rice was carried out using the Triticeae-Gene Tribe (http:// wheat.cau.edu.cn/TGT/) . Furthermore, to analyze the potential candidate genes, SNPs in the target regions were collected using the exome sequencing result.

Phenotypic variation of flag leaf-related traits
The phenotypic means of flag leaf-related traits (including FLL, FLW and FLA) for the parents and the population were listed in Table 1, as well as basic statistics from plants grown in six environments and the BLUE datasets. Significant differences in FLL, FLW, and FLA between W7268 and CY12 were observed ( Figure 1 and Table 1). W7268 had significantly higher values for FLL (except E3 and E4), FLW, and FLA than those of CY12. The flag leaf-related traits of the WC12 population ranged from 19.50 to 33.80 cm for FLL, 1.62 to 2.95 cm for FLW, and 27.80 to 73.20 cm 2 for FLA, in BLUE datasets, respectively. The estimated H 2 for FLL, FLW, and FLA were 0.88, 0.92, and 0.89, respectively. These results indicated that FLL, FLW, and FLA were environmentally stable and were mainly determined by genetic factors. A pattern of continuous distributions for FLL, FLW, and FLA was observed in each environment and the BLUE dataset of the WC12 population, suggesting that they were common quantitative traits and controlled by multi-genes ( Figure 2). Furthermore, significant and positive correlations among FLL, FLW, and FLA were detected among six environments, with Pearson's correlations (r) of 0.21-0.93 in the WC12 population (Supplementary Figure 1).
Two major QTL associated with FLA were detected. QFla.cib-2B.2 was detected in four environments and the BLUE dataset with LOD values ranging from 3.51 to 11.27. It explained 6.33-25.33% of the phenotypic variance. QFla.cib-2D, detected in three environments and the BLUE dataset, had LOD values of 12.88-18.05 and accounted for 28.6-37.19% of the phenotypic variance. The positive alleles of QFla.cib-2B.2 and QFla.cib-2D were contributed by CY12 and W7268, respectively ( Figure 4 and Table 2).
According to the interval of QFll/Fla.cib-2B, KASP markers (KA01 to KA29) were developed and integrated into the genetic map based on the exome sequencing result of the two parents (Supplementary Table 3). Among them, four markers, KA12, KA13, KA15, and KA17, were found to be closely linked to QFll/ Fla.cib-2B ( Figure 4A).

Conditional QTL analysis for FLA
FLA is a complex trait composed of FLL and FLW. We performed conditional QTL analysis of FLA in the WC12 population to further evaluate the effect of FLL and FLW on FLA. The results showed that when FLA was conditional on FLL, the LOD values of QFla.cib-2B.2 and QFla.cib-2D were significantly reduced; when FLW was conditioned, the LOD values were lower than those in the unconditional analysis but remained at a high level (Supplementary Figure 2). These results indicated that FLL was primarily responsible for FLA in the WC12 population.

Effects of major QTL on corresponding traits
In the WC12 population, we identified seven major QTL (Table 2). Their effects on corresponding traits were assessed on the basis of flanking markers. As expected, lines in the WC12 population with the positive alleles at the three loci showed Coefficients of pairwise Pearson's correlations between flag leaf length (A), flag leaf width (B), and flag leaf area (C) and yield-related traits in the WC12 population. *, *** and **** represent significant at P < 0.05, P < 0.005 and P < 0.001, respectively.
significantly higher values of the corresponding traits in all environments and the combined data than those with the negative alleles (except for E3 of QFll.cib-2D, E4 of QFlw.cib-2D.1, and E3 of QFla.cib-2D) (Figure 4).

Effects of major QTL on yield-related traits in the mapping population
We further evaluated the effects of major QTL on yieldrelated traits using the BLUE datasets. Compared with lines harboring the alleles from CY12 at QFll/Fla.cib-2B, lines with the alleles from W7268 had lower SNS (P < 0.05), but higher SC (P < 0.005). QFll/Fla.cib-2B had no significant effects on GNS, TGW, and FTN ( Figure 5A). For QFlw.cib-4B.3, lines containing the alleles from W7268 had higher GNS (P < 0.001) but lower TGW (P < 0.001) than lines containing the alleles from CY12 ( Figure 5B). No significant difference in SNS, SC, and FTN was observed between the two groups. For QFll/Flw/Fla.cib-2D, lines possessing the alleles from W7268 had significantly higher SNS (P < 0.001), GNS (P < 0.001), and FTN (P < 0.005) than lines containing the alleles from CY12, but had significantly lower TGW (P < 0.001) ( Figure 5C).

Potential candidate genes for QFll/ Fla.cib-2B
According to the CS reference genome (IWGSC RefSeq v2.1), there were 56 annotated high-confidence genes in the interval of QFll/Fla.cib-2B (Supplementary Table 4). Expression pattern analysis showed that 36 (Triticeae Multi-omics Center) and 23 (Wheat Expression Browser) genes were expressed in the leaf, respectively (Supplementary Figure 4). Gene annotation and orthologous gene analyses (Supplementary Table 4), combined with previous expression pattern analysis, suggested that TraesCS2B03G0222800 and TraesCS2B03G0230000 were likely to be associated with flag leaf development and growth. Six SNPs, two located in introns and four in the coding region, were detected in TraesCS2B03G0230000 between W7268 and CY12 by sequence analysis. Among the four SNPs in the exons, there were two synonymous SNPs and two non-synonymous SNPs (Arg to Trp, and Ala to Val) (Supplementary Table 5).

Discussion
Comparison of the major QTL to those reported previously In this study, seven major and stable QTL were identified in multiple environments, explaining 6.33-53.12% of phenotypic variations. These QTL will benefit cloning and marker-assisted selection (MAS) in wheat breeding.
Because Ppd-B1 has a similar physical interval (63.36-63.37 Mb on chromosome 2B) to QFll/Fla.cib-2B, we used the functional marker of Ppd-B1 to perform the comparison. By genotyping the WC12 lines, we integrated it into the genetic map and the result showed that Ppd-B1 was genetically separate from QFll/Fla.cib-2B ( Figure 4A). As a result, QFll/Fla.cib-2B is likely a novel locus.
Relationships between flag leaf size and yield-related traits and pleiotropic effects of major QTL Optimizing flag leaf size, including length, width, and area, plays an important role in increasing grain yield (Zhao et al., 2018). In wheat, the flag leaf size was significantly correlated with yield-related traits (Cui et al., 2003;Wang et al., 2011;Fan et al., 2015;Liu et al., 2018;Zhao et al., 2022). In the present study, we also found that FLL, FLW, and FLA showed a positive correlation to SNS and GNS, and a negative correlation to SC, FTN, and TGW ( Figure 3). This result suggested that larger flag leaf benefits from forming more SNS and GNS but hinders the tillering and grain weight increase.
W7268 is an elite line characterized by high SNS and GNS in the Yangzi River region of China. It has been widely used in wheat breeding programs and several varieties have been selected using W7268 as one of the parents. In the present study, we found that QFlw.cib-4B.3 and QFll/Flw/Fla.cib-2D have various pleiotropic effects on SNS and GNS, which may contribute to the higher yield of W7268. Moreover, QFll/Fla.cib-2B from CY12 also had additive effects on SNS and GNS and may be an unreported locus with potential to increase yield ( Figure 5). Thus, cloning and utilization of these QTL will be valuable for grain yield improvement by optimizing flag leaf size.
To analyze the utilization of the positive alleles of the major QTL during artificial selection, we used their flanking markers to genotype 135 wheat accessions (75 of landraces and 60 of cultivars), and the results were shown in Supplementary Based on our results, QFlw.cib-4B.3 and QFll/Flw/Fla.cib-2D may be allelic to Rht-B1 and Ppd-D1, respectively. And the positive allele for QFlw.cib-4B.3 may be Rht-B1b and for QFll/Flw/Fla.cib-2D may be Ppd-D1b. According to previous reports (Yang et al., 2009;Guo et al., 2010;Tang et al., 2012;Wang et al., 2013;Xu et al., 2014;Bai et al., 2015;Zhang et al., 2016), the distribution of Rht-B1a and Ppd-D1a was enriched during the breeding process compared to Rht-B1b and Ppd-D1b in China, which was consistent with our study. The dwarf gene Rht-B1b could reduce the thousand-grain weight, which may be the reason why breeders did not prefer this genotype (Xu et al., 2014;Zhang et al., 2016). Ppd-D1a, insensitive to photoperiod, facilitated flowering and allowed wheat to finish filling before the onset of summer heat. (Yang et al., 2009;Guo et al., 2010;Bai et al., 2015).
As shown in Supplementary Figure 8, significant differences in GYP were identified between lines with positive and negative alleles, with one, three, and four environments for QFlw.cib-4B.3, QFll/Fla.cib-2B, and QFll/Flw/Fla.cib-2D, respectively. As expected, lines with the negative alleles of QFlw.cib-4B.3 and QFll/Flw/Fla.cib-2D showed higher GYP values, consistent with their distribution changes during the artificial selection. For QFll/Fla.cib-2B, lines with the negative allele had a higher GYP than the positive allele.
The significant effects of these flag leaf-related QTL on different traits, their distribution in different accessions, and various impacts on yield, suggested that utilization of them will help to improve grain yield by optimizing flag leaf size.

Potential candidate genes for QFll/ Fla.cib-2B
In the physical interval of QFll/Fla.cib-2B in the CS genome, 36 (Triticeae Multi-omics Center) and 23 (Wheat Expression Browser) genes expressed in leaves were screened (Supplementary Figure 4 and Supplementary Table 4). According to expression-pattern analysis, gene annotation, and ortholog analysis, two candidate genes TraesCS2B03G0222800 and TraesCS2B03G0230000 were preliminarily selected to be associated with flag leaf development. Among them, TraesCS2B03G0222800 encodes a pseudo-response regulator and is an ortholog of rice Os07g0695100, which is also known as DTH7 (Days to heading 7)/Ghd7.1(Grain number, plant height, and heading date 7), a pleiotropic gene controlling tassel stage, plant height and the number of glumes per spike in rice (Yan et al., 2013). A recent study has shown that Ghd7.1 mutations resulted in reduced leaf size in rice, and allelic variance analysis has verified that Ghd7.1 is a functional candidate gene for leaf size (Tang et al., 2018). TraesCS2B03G0230000 encodes an omega-3 fatty acid desaturase, which is associated with the production of alphalinolenic acid, a major class of fatty acids found in the membrane lipid cells of higher plants. (John et al., 1995). Plastid omega-3 desaturase activity is necessary to increase the level of alphalinolenic acid in extra plastid lipids during leaf development (Horiguchi et al., 1998). Sequence analysis also revealed abundant sequence polymorphisms in TraesCS2B03G0230000 between W7268 and CY12 (Supplementary Table 5). Therefore, TraesCS2B03G0222800 and TraesCS2B03G0230000 may be potential candidates for map-based cloning in the future.

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
LC carried out the entire study and wrote the draft manuscript. ZX constructed the mapping population. XF and QZ assisted in field trials. QY, SL, XL, CJ, DL, and FM participated in phenotyping and data analysis. TW guided the entire study and discussed the results. BF designed the experiments, participated in data analysis, discussed results, and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding
This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Precision Seed Design and Breeding, XDA24030402) and the CAS "Light of West China" (2023XBZG_XBQNXZ_A), the Sichuan Science and Technology Program (2022ZDZX0014). on FLL, FLW, and FLA, respectively. + and − represent lines with and without the positive alleles of the target QTL based on the flanking markers of the corresponding QTL, respectively. **** represents significant at P < 0.001 SUPPLEMENTARY FIGURE 4 Expression pattern of genes within the QFll/Fla.cib-2B interval. (A-D) represent leaf, grain, spike, and root by Triticeae Multi-omics Center, respectively; (E-H) represent leaves/shoots, grain, spike, and roots by Wheat Expression Browser, respectively. Effect analysis of QFll/Fla.cib-2B (A), QFlw.cib-4B.3 (B), and QFll/Flw/ Fla.cib-2D (C) on grain yield per plant (GYP) in WC12 population. W7268 and CY12 represent lines with alleles from W7268 and CY12, respectively; *, **, ***, and **** represent significant at P < 0.05, P < 0.01, P < 0.005, and P < 0.001, respectively.