Identification of three QTLs that additively affect heading time in bread wheat (Triticum aestivum L.) by QTL-seq approach

Identification of three QTLs that additively affect heading time in bread wheat (Triticum aestivum L.) by QTL-seq approach | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Identification of three QTLs that additively affect heading time in bread wheat (Triticum aestivum L.) by QTL-seq approach Shoya Komura, Fuminori Kobayashi, Youko Oono, Hirokazu Handa, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4923172/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 14 Nov, 2024 Read the published version in Euphytica → Version 1 posted 9 You are reading this latest preprint version Abstract Optimizing the timing of heading is crucial for achieving stable production in bread wheat ( Triticum aestivum L.). We analyzed a breeding line, ‘B1-0393’, that headed 5 days earlier than the Japanese elite cultivar ‘Kitahonami’. To identify quantitative trait loci (QTLs) responsible for this difference, we conducted QTL-seq using F 3 progeny of a cross between ‘Kitahonami’ and ‘B1-0393’. We detected QTLs on chromosomes 2D ( QDth.kupg-2D ), 3D ( QDth.kupg-3D ), and 4A ( QDth.kupg-4A ), explaining 10.25%, 2.00%, and 8.24%, respectively, of phenotypic variance. The QDth.kupg-2D locus corresponded to the major photoperiod-regulating gene Photoperiod-D1 ( Ppd-D1 ). ‘Kitahonami’ had the photoperiod-sensitive Ppd-D1b allele and ‘B1-0393’ had the insensitive Ppd-D1a allele. QDth.kupg-4A overlapped with previously reported QTLs, while QDth.kupg-3D was a novel QTL. Segregation analysis using F 2 and F 3 plants confirmed that alleles from ‘B1-0393’ at all three loci accelerated heading and that the three QTLs had additive effects on days to heading. Our findings could be valuable for regulating heading time to optimize wheat yield. Wheat Heading Quantitative trait locus Bulked segregant analysis Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction The appropriate timing of the transition from vegetative to reproductive growth and subsequent development is crucial for ensuring grain yields and quality. In bread wheat ( Triticum aestivum L.), the optimal heading date is one of the aims in breeding, as it allows growers to avoid winter cold stress while mitigating subsequent heat and drought stresses. Especially, the harvest season is close to the rainy season in Japan. Such conditions can lead to pre-harvest sprouting and Fusarium head blight. Thus, early heading is one of the most essential aims in the breeding program in Japan to avoid losing yield and quality (Hoshino and Seko 1996 ). Vernalization and photoperiod are major wheat heading components (Kato and Shunji 1991 ; Hyles et al. 2020 ). In winter wheat, heading requires prolonged exposure to cold, which ensures the transition to reproductive growth in warm conditions. In contrast, plants with spring growth habit can head without exposure to cold. Vernalization is genetically controlled by VERNALIZATION-1 ( Vrn-1 ), VERNALIZATION-2 ( Vrn-2 ), and Wheat FLOWERING LOCUS T ( WFT ) (Yan et al. 2003 ; Shimada et al. 2009 ; Chen and Dubcovsky 2012 ). The spring growth habit is conferred by recessive loss-of-function alleles of Vrn-2 or by at least one dominant allele of Vrn-1 (Yan et al. 2004 ; Hyles et al. 2020 ). The sensitivity to the duration of daylight, also known as photoperiod sensitivity, is also a crucial element in deciding the timing of heading. Photoperiod sensitivity is regulated mainly by the Photoperiod-1 ( Ppd-1 ) gene, which is homologous to Arabidopsis thaliana PSEUDO RESPONSE REGULATOR 7 ( PRR7 ) (Beales et al. 2007 ; Shaw et al. 2012 ). Ppd-1 alleles are classified into photoperiod sensitive ( Ppd-A1b , Ppd-B1b , and Ppd-D1b ) and photoperiod insensitive ( Ppd-A1a , Ppd-B1a , and Ppd-D1a ) (Beales et al. 2007 ; Shaw et al. 2012 ; Hyles et al. 2020 ). The photoperiod-sensitive alleles show diurnal expression patterns and promote heading under long-day conditions. In contrast, the photoperiod-insensitive alleles are hypermorphic because of deletions in their promoter regions, resulting in dysregulation of Ppd-1 expression (Beales et al. 2007 ; Díaz et al. 2012 ; Nishida et al. 2013 ). The misexpression of Ppd-1 induces an increased WFT expression and promotes heading regardless of daylength (Beales et al. 2007 ; Shaw et al. 2012 ). In the case of Ppd-B1 , photoperiod insensitivity is caused by an increased copy number, which results in high basal expression levels of Ppd-B1 (Díaz et al. 2012 ). Since Vrn-1 and Ppd-1 significantly affect heading timing, allelic variations of these genes have been used to adapt plant heading to various environments (Kiss et al. 2014 ; Chen et al. 2018 ; Mizuno et al. 2022 ). The Ppd-D1a allele is more common in wheat cultivars grown in southern Japan than in northern regions (Seki et al. 2011 ), whereas the Ppd-A1a allele is present only in cultivars in northern Japan (Seki et al. 2013 ). Iwaki et al. ( 2001 ) found that Vrn-D1 is more common in South Asia than other spring habitat alleles. Kiss et al. ( 2014 ) and Mizuno et al. ( 2022 ) revealed that the environment significantly modifies the effects of Vrn-1 and Ppd-1 , and some alleles alter heading dates only in specific environments. Therefore, an appropriate adoption of Vrn-1 and Ppd-1 alleles and their combinations is essential to optimize the timing of heading. Even if the vernalization requirement and day length are satisfied, heading time can vary to some extent. These minor variations are caused by genes referred to as earliness per se ( Eps ) genes. Since Eps genes alter heading dates only slightly, they contribute to local adaptation and can be used for fine-tuning (Maeda and Nakamichi 2022 ). Although numerous Eps quantitative trait loci (QTLs) have been identified, only a limited number of Eps genes have been cloned (Snape et al. 2001 ; Singh K. et al. 2022 ). In wheat, circadian clock genes WHEAT LUX ARRHYTHMO/PHYTOCLOCK 1 ( WPCL1 ), EARLY FLOWERING 3 ( ELF3 ), and WHEAT WD REPEAT 1 ( WWDR1 ) have been identified as Eps genes (Mizuno et al. 2012 ; Alvarez et al. 2016 ; Zikhali et al. 2016 ). Although loss-of-function in one or two homoeologs has a small effect, the null mutants of WPCL1 or WWDR1 show extremely early heading, which results in severe yield reduction (Mizuno et al. 2012 ; Hashimoto et al. 2021 ). The degree of earliness influenced by ELF3 is altered by the prevailing temperature (Ochagavía et al. 2019 ). Therefore, further exploration of loci that could be used to manipulate heading is important to enable flexible breeding to optimize cultivars to the target geographic region or environment. The detection of a QTL is an essential strategy for improving breeding cultivars and identifying causal genes. QTL-seq is a rapid and efficient approach to detect QTLs based on bulked segregant analysis (Takagi et al. 2013 ). Using two parents that show opposite extreme phenotypes and pooled DNA bulks from their progeny, QTL-seq detects genomic regions linked to the phenotype of interest on the basis of ∆SNP-index distribution. This strategy has detected QTLs in rice, pigeon pea, tomato, maize, barley, and oilseed rape (Takagi et al. 2013 ; Illa-Berenguer et al. 2015 ; Hisano et al. 2017 ; Singh V. et al. 2022 ; Wang et al. 2022 ; Ni et al. 2023 ). In bread wheat, because of its large genome (2 n = 6 x = 42, ~ 17 GB), QTL-seq has been used mainly in combination with exome sequencing or RNA sequencing (Li et al. 2020 ; Jiang et al. 2022 ; Ji et al. 2023 ). These approaches have detected QTLs for agronomically important traits, such as grain size, stripe rust resistance, and heading date. However, with the release of pangenome reference sequences for 10 wheat cultivars (Walkowiak et al. 2020 ) and the reduction in cost of next-generation sequencing, whole-genome sequencing of the bread wheat genome has become more accessible. In this study, we aimed to identify candidate genomic regions responsible for the difference in days to heading between the Japanese elite bread wheat ‘Kitahonami’ and the breeding line ‘B1-0393’. Since ‘B1-0393’ headed 5 days earlier than ‘Kitahonami’, the early-heading gene(s) in ‘B1-0393’ may be used to fine-tune heading date. To identify the causal genomic regions, we performed whole-genome sequencing, conducted QTL-seq, and detected three QTLs associated with heading dates, including the Ppd-D1 locus. Segregation tests in the F 2 and F 3 populations confirmed that the three QTLs had additive effects, and the F 3 plants harboring triple early-heading alleles had the earliest heading dates. This study provides valuable insights into QTL mapping in wheat to control heading date. Materials and Methods Plant materials ‘B1-0393’ was found in a gamma-irradiated wheat population and heads earlier than ‘Kitahonami’, whichwhi was released by the Kitami Agricultural Experiment Station, Hokkaido (Yanagisawa et al. 2007 ). The heading dates of ‘Kitahonami’ and ‘B1-0393’ were evaluated in the experimental field of the Institute of Crop Science, NARO (Tsukuba, Japan; lat. 36.1° N, long. 140.6° E), under natural field conditions in the 2018–19 (sown 24 October 2018) and 2019–20 (sown 31 October 2019) cropping seasons. ‘Kitahonami’ and ‘B1-0393’ were crossed, generating F 2 ( n = 130) and F 3 ( n = 598) populations. To measure heading time, ‘Kitahonami’, ‘B1-0393’, and their F 2 progeny were sown on 27 October 2020 (2020–21 season), and their F 3 progeny were sown on 29 October 2021 (2021–22 season) in the same field under natural conditions. The first appearance of the spike tip was recorded as the heading date in both seasons. Genome sequencing and bioinformatics analysis Total DNA was extracted from a leaf of ‘B1-0393’ using a DNeasy Plant Mini Kit (Qiagen, Hilden, Germany). For QTL-seq, 20 individuals with the earliest heading dates and 20 with the latest heading dates were selected from 598 F 3 progeny. The total DNA of each individual was extracted as above. and bulked to produce an early heading bulk and a late heading bulk. DNA libraries of ‘B1-0393’ and each bulk were constructed using a TruSeq DNA library Preparation Kit (Illumina, San Diego, CA, USA). Genome sequencing was performed on an Illumina NovaSeq in paired-end 150-bp mode. The sequencing reads of the wild-type ‘Kitahonami’ were obtained from our previous study (Komura et al. 2022 ). The quality of sequence reads was checked with FASTQC v. 0.11.7 and the reads were filtered by Trimmomatic v. 0.33 (Bolger et al. 2014 ) with an average minimum Phred quality score per 4 bp of < 20 and a length of < 40 bp. We used BWA v. 0.7.17 with default options to align the filtered paired-end reads to the reference genomes of bread wheat cultivars ‘Norin 61’ and ‘Chinese Spring’ (‘CS’) (IWGSC 2018; Walkowiak et al. 2020 ; Shimizu et al. 2021 ). PCR duplications were removed using SAMtools v. 1.9 (Li et al. 2009 ). The average number and genome coverage of aligned reads were calculated using BBMap v. 37.77 (Bushnell 2014 ). SNPs and indels were called using BCFtools v. 1.9 (Li et al. 2009 ) with options “-A” and “-Q 0”. Variants differing between the parents were extracted using BCFtools v. 1.9. We removed variants that did not meet the criteria of (1) mapping quality of reads containing variants of ≥ 40, (2) a homozygous variant detection site, (3) ≥ 5 aligned reads with the altered genotype at the variant detection site but no reads of the reference allele in one parent, and (4) ≥ 5 reads with the reference allele but no variant-type read counts in the other parent. The effects of SNPs and indels were predicted in SnpEff v. 4.3t (Cingolani et al. 2012 ). To calculate SNP density, the number of SNPs per 1 Mbp was counted in Python v. 3.9.0. The R package ‘ggplot2’ (Wickham, 2016 ) was used to visualize the distribution of the moving average of SNP density over the chromosomes. QTL-seq analysis (Takagi et al. 2013 ; Sugihara et al. 2022 ) was performed according to the pipeline of Sugihara et al. ( 2022 ). The reference genome of ‘Norin 61’ was used for QTL-seq analysis. The ΔSNP-index of the two bulks were estimated at the positions where SNPs were detected between the ‘Kitahonami’ and ‘B1-0393’. To narrow down causal variants in the candidate regions, indel-index was also estimated. The SNP-index and indel-index were calculated by dividing the number of reads supporting the alternative sequence by the total number of reads aligned to the position. The ΔSNP-index was calculated as the difference between the SNP-index of the early-heading bulk and that of the late-heading bulk. Marker development and genotyping For segregation analysis of QTL regions, we designed cleaved amplified polymorphic sequence (CAPS) markers on the basis of SNPs between ‘Kitahonami’ and ‘B1-0393’. For CAPS markers, we designed subgenome-specific primers by BLAST searches against the ‘Norin 61’ reference sequence through WheatOmics (Ma et al. 2021 ). The markers used in this study are listed in Supplementary Table S1 . ‘Kitahonami’, ‘B1-0393’, and their 130 F 2 and 598 F 3 progeny were used for genotyping. Total DNA was extracted as above and PCR was performed using Quick-Taq HS DyeMix (Toyobo, Osaka, Japan) as follows: pre-denaturation at 94°C for 2 min; 35 cycles of denaturation at 94°C for 30 s, annealing at 62°C for 30 s, and extension at 68°C for 30 s; and post-extension at 68°C for 1 min. The PCR products were digested with the restriction enzymes listed in Supplementary Table S1 and electrophoresed in 2% agarose gel. To confirm the sequences of Ppd-D1 in ‘Kitahonami’ and ‘B1-0393’, we checked the read alignment using Integrative Genomics Viewer (Robinson et al. 2011 ). Ppd-D1 genotyping was performed using allele-specific primers for Ppd-D1a or Ppd-D1b , as described by Nishida et al. ( 2013 ). Statistical analysis The normality of the distributions in the F 2 and F 3 populations was confirmed using the Shapiro–Wilk test. Variance among genotypes was assessed by an analysis of variance (ANOVA) or the Kruskal–Wallis test. For multiple comparisons, the Tukey–Kramer test or the Steel–Dwass test was used. To detect epistatic interactions between three QTLs, we used ANOVA analysis with the following general linear model (GLM): $$\:{GLM}_{iJl}=\mu\:+{QTL}_{1i}+{QTL}_{2j}+{QTL}_{1}{QTL}_{2ij}+{e}_{iil}$$ where µ is the population mean value; \(\:{QTL}_{1i}\:\) and \(\:{QTL}_{2j}\) indicate the effects of the i th and j th QTLs; and \(\:{e}_{iil}\) is experimental error. The additive and dominant effects and phenotypic variation explained by each locus were determined using the R package ‘R/qtl’ (Broman et al. 2003 ) with the following model: $$\:{HD}_{ijl}=\mu\:+{Ppd\text{-}D1}_{i}+{QDth.kupg\text{-}3D}_{j}+{QDth.kupg\text{-}4A}_{l}+{e}_{ijl}$$ where \(\:\:{Ppd\text{-}D1}_{i}\) , \(\:{QDth.kupg\text{-}3D}_{j}\) , and \(\:{QDth.kupg\text{-}4A}_{l}\) indicate the effects of the Ppd-D1 , QDth.kupg-3D , and QDth.kupg-4A loci; and \(\:{e}_{ijl}\) is experimental error. Additive effect ( \(\:a\) ) and dominance effect ( \(\:d\) ) were calculated as follows: $$\:a=\frac{({A}_{1}{A}_{1}-{A}_{2}{A}_{2})}{2}$$ $$\:d={A}_{1}{A}_{2}-\frac{({A}_{1}{A}_{1}+{A}_{2}{A}_{2})}{2}$$ where \(\:{A}_{1}{A}_{1}\) and \(\:{A}_{2}{A}_{2}\) are the mean values of plants homozygous for the ‘Kitahonami’ and ‘B1-0393’, respectively, alleles; and \(\:{A}_{1}{A}_{2}\) is the mean of heterozygous plants. Estimation of heritability Phenotypic variance in the progeny ( \(\:{V}_{F}\) ) includes genetic and environmental variances, whereas that in the parents ( \(\:{V}_{P1}\) and \(\:{V}_{P2}\) ) includes only environmental variances. Therefore, the broad-sense heritability ( \(\:{H}^{2}\) ) was estimated as (Hartl 2020 ): $$\:{H}^{2}=\frac{{V}_{F}-\raisebox{1ex}{$1$}\!\left/\:\!\raisebox{-1ex}{$2$}\right.({V}_{P1}+{V}_{P2})}{{V}_{F}}$$ The narrow-sense heritability ( \(\:{h}^{2}\) ) was estimated as (Hartl 2020 ): $$\:{M}_{F3\_late}^{”}={M}_{F2}+{h}^{2}({M}_{F2\_late}^{{\prime\:}}-{M}_{F2})$$ $$\:{M}_{F3\_early}^{”}={M}_{F2}-{h}^{2}({M}_{F2}-{M}_{F2\_early}^{{\prime\:}})$$ By rearranging the above equations, the following equation can be obtained: $$\:{h}^{2}=\frac{({M}_{F3\_late}^{”}-{M}_{F3\_early}^{”})}{({M}_{F2\_late}^{{\prime\:}}-{M}_{F2\_early}^{{\prime\:}})}$$ where \(\:{M}_{F2}\) is the mean of the F 2 population; \(\:{M}_{F2\_late}^{{\prime\:}}\) is the mean of F 2 individuals that headed later than a threshold; \(\:{M}_{F2\_early}^{{\prime\:}}\) is the mean of F 2 individuals that headed earlier than a threshold; and \(\:{M}_{F3\_late}^{”}\) and \(\:{M}_{F3\_early}^{”}\) are the means of F 3 individuals generated from the selected F 2 individuals. We set the thresholds for selecting F 2 individuals as > 179 days to heading for late heading and < 174 days to heading for early heading. Results Characterization of early-heading line ‘B1-0393’ ‘B1-0393’ headed 5 days earlier than ‘Kitahonami’ in the 2018–19 season and 3 days earlier in the 2019–20 season (Fig. 1 b). To characterize nucleotide variations between ‘Kitahonami’ and ‘B1-0393’, we resequenced the whole genome of ‘B1-0393’ with short reads and used the sequencing reads of ‘Kitahonami’ from our previous study (Komura et al. 2022 ). After quality control, 2.7 billion reads of ‘Kitahonami’ and 2.1 billion reads of ‘B1-0393’ were filtered, and 2.3 and 1.8 billion reads, respectively, were aligned to the ‘Norin 61’ reference genome (Table 1 ). These reads covered 96.3% and 96.1%, respectively, of the reference genome, with an average depth-of-coverage of 15.08 and 18.56. We detected 255,369 indels and 4,618,884 SNPs between ‘Kitahonami’ and ‘B1-0393’ (Supplementary Table S2 ). To evaluate genetic divergence between ‘Kitahonami’ and ‘B1-0393’, we calculated SNP density per 10 Mbp (Fig. 1 c). The SNPs were unevenly distributed on the chromosomes, with interspersed regions of high and low SNP density. For example, SNP density was high in the 520–640 Mbp interval of chromosome 4A, but it was low in the 0–520 Mbp and 640–720 Mbp intervals, suggesting that the genetic background was similar in these regions. Fewer SNPs were detected in the D genome than in the other genomes, suggesting a closer genetic background of the D genome (Supplementary Fig. S1 ). These results indicate that the genetic background of ‘B1-0393’ is similar to that of ‘Kitahonami’ in some regions but is divergent in others (Fig. 1 c). Table 1 Summary of alignments and coverages of genome sequencing of ‘Kitahonami’ and the early-heading line ‘B1-0393’. Name Total filtered reads a Aligned reads after removal of low-quality reads b Reference bases covered Average depth-of-coverage Kitahonami c 2,724,430,844 (93.63%) 2,334,502,744 (85.69%) 96.32% 15.08 B1-0393 2,151,447,818 (93.26%) 1,846,489,772 (85.83%) 96.10% 18.56 Early heading bulk 2,477,029,056 (85.58%) 2,148,834,659 (86.75%) 96.71% 21.47 Late heading bulk 2,530,545,270 (92.51%) 2,246,339,991 (88.77%) 96.83% 21.66 a Rate of filtered reads = Total filtered reads / Total raw reads × 100. b Rate of filtered reads = Aligned reads after removal of low-quality reads / Total filtered reads × 100. c ‘Kitahonami’ sequencing data are reproduced from Komura et al. ( 2022 ). The F 2 population showed a Gaussian distribution of days to heading ranging from 169 to 185, while ‘Kitahonami’ headed at 179 days and ‘B1-0393’ at 174 days (Fig. 2 , Table 2 ). In the F 3 population, progeny headed from 175 to 191 days, while ‘Kitahonami’ headed at 180 days and ‘B1-0393’ at 177 days. The distribution of days to heading was skewed towards early heading and showed transgressive segregation (Fig. 2 , Table 2 ), suggesting that multiple loci determine the timing of heading. The broad-sense heritability was 0.65 in the F 2 population and 0.57 in the F 3 population, while the narrow-sense heritability was 0.29 (Table 2 ). The former includes additive and dominant variances, while the latter includes only additive variances. The narrow-sense heritability was almost half of the broad-sense heritability, the difference indicating the degree of dominant variances. Therefore, the difference in heading time between ‘Kitahonami’ and ‘B1-0393’ was inferred to be a combination of dominance and additive effects. Table 2 Days to heading in the F 2 and F 3 populations. Population No. of plants Mean SD Skewness Kurtosis Broad-sense heritability Narrow-sense heritability F 2 130 176.13 2.73 0.24 0.23 0.65 0.29 F 3 596 179.47 1.50 1.22 3.91 0.57 Detection of QTLs for heading date by QTL-seq From the F 3 population, we selected 20 plants that headed on 176 or 177 days as the early-heading bulk and 20 plants that headed on 185–189 days as the late-heading bulk, and resequenced the whole genomes of the two bulks. After quality control, the filtered 2.5 billion reads per bulk were aligned to the ‘Norin 61’ reference sequence. Of the filtered reads, 2.1 billion reads from the early-heading bulk and 2.2 billion reads from the late-heading bulk could be aligned (Table 1 ). The aligned reads of each bulk covered 96.71% (21.47× average depth-of-coverage) and 96.83% (21.66× \(\:)\:\) of the reference sequences. QTL-seq analysis revealed that the moving average of ΔSNP-index in the regions of 31–37 Mbp on chromosome 2D, 490–520 Mbp on chromosome 3D, and 610–650 Mbp on chromosome 4A exceeded the 99% or 95% confidence intervals (CIs) (Fig. 3 , Supplementary Fig. S2 ). These QTL regions were designated as QDth.kupg-2D , QDth.kupg-3D , and QDth.kupg-4A . QDth.kupg-2D houses the known heading-related gene Ppd-D1 . The alignment of the reads to the ‘Norin 61’ reference sequence revealed that ‘Kitahonami’ has a 5-bp deletion in exon 7 of Ppd-D1 , resulting in a premature stop codon and a truncated protein. ‘Norin 61’ has a photoperiod-insensitive allele Ppd-D1 . Therefore, to investigate the alleles of Ppd-D1 in ‘Kitahonami’ and ‘B1-0393’, we used a reference sequence of ‘CS’, which carries the photoperiod-sensitive allele Ppd-D1b . Read alignment to the ‘CS’ reference sequence showed that ‘B1-0393’ has a 2089-bp deletion in the 5′ upstream region and confirmed the 5-bp deletion in exon 7 in ‘Kitahonami’. In addition, a 16-bp deletion in exon 8 was identified in both parents (Supplementary Fig. S3). The Ppd-D1 allele in ‘B1-0393’ was identical to that in ‘Winter-Abukumawase’, which is reportedly the photoperiod-insensitive Ppd-D1a allele (Haplotype Ppd-D1a.1 ) (Nishida et al. 2013 ). The Ppd-D1 allele in ‘Kitahonami’ was the same as that in ‘Norstar’ and ‘Chihokukomugi’, which is referred to as the photoperiod-sensitive Ppd-D1b allele (Haplotype Ppd-D1b.2 ) (Beales et al. 2007 ; Nishida et al. 2013 ). According to Nishida et al. ( 2013 ), plants with the Ppd-D1a allele head earlier than plants with the Ppd-D1b allele. Therefore, we assumed Ppd-D1 to be the causal gene in QDth.kupg-2D . In the QDth.kupg-4A region, we found 9518 variants with the SNP-index and indel-index exceeding the 95% CIs. Among them, 6 variants caused frameshift, nonsense, or stop-loss mutations, and 54 caused missense mutations (Supplementary Tables S3, S4). In the QDth.kupg-3D region, we found 1504 variants with the SNP-index and indel-index exceeding the 95% CIs. Three variants caused missense mutations (Supplementary Table S5). No known genes related to wheat heading timing were found in the QDth.kupg-3D or QDth.kupg-4A regions. Validation of QDth.kupg-2D , QDth.kupg-3D , and QDth.kupg-4A To confirm the association between the detected QTLs and heading date, we conducted a segregation analysis using 130 F 2 plants. Days to heading differed significantly depending on the genotypes of Ppd-D1 , QDth.kupg-3D , and QDth.kupg-4A (Fig. 4 a, Supplementary Table S6). Plants with the ‘B1-0393’ alleles at all three QTLs headed earlier. The phenotypic variation explained (PVE) was 10.25% by Ppd-D1 , 2% by QDth.kupg-3D , and 8.24% by QDth.kupg-4A (Table 3 ). The dominance effect of Ppd-D1 was 0.06, indicating that Ppd-D1a and Ppd-D1b have similar influence on heading, and heterozygous individuals had phenotypes intermediate between those of individuals homozygous for the parental alleles (Fig. 4 a). The additive effects were 0.51 for QDth.kupg-3D and 1.01 for QDth.kupg-4A , and dominance effects were 0.21 and 0.75, respectively (Table 3 ). The segregation analysis indicated incomplete dominance by QDth.kupg-3D and a dominant phenotype of the ‘Kitahonami’ allele of QDth.kupg-4A (Fig. 4 a). The PVE and dominant effect values indicated that Ppd-D1 had the largest and QDth.kupg-3D had the smallest effect on heading date. Table 3 Quantitative trait loci (QTLs) for days to heading detected by QTL-seq in F 2 population. Locus Chr. QTL range (Mbp) Marker position (bp) PVE c Additive effect Dominance effect QDth.kupg-2D 2D 31–37 a 34,671,330 10.25% 1.35 0.06 QDth.kupg-3D 3D 490–520 a 496,096,772 2.00% 0.51 0.21 QDth.kupg-4A 4A 610–650 b 632,986,650 8.24% 1.01 0.75 a ∆SNP-index > 95% confidence interval. b ∆SNP-index > 99% confidence interval. c Phenotypic variance explained. Since no significant interactions between these QTLs were found in ANOVA, we expected these QTLs to function in an additive manner (Supplementary Table S6). To investigate the additive effects of the three QTLs on heading date, we determined the genotypes of 598 F 3 plants. The plants carrying the three ‘B1-0393’-type alleles headed an average of 3.19 days earlier than the plants carrying three ‘Kitahonami’-type alleles (Fig. 4 b). The plants carrying the ‘B1-0393’-type allele at QDth.kupg-3D or QDth.kupg-4A with the ‘B1-0393’-type Ppd-D1a allele also headed earlier (by 2.62 and 2.92 days on average) than those with three ‘Kitahonami’-type alleles. The effect of the ‘B1-0393’-type allele was strongest for Ppd-D1 and weakest for QDth.kupg-3D . The effect on heading date of plants carrying homozygous ‘B1-0393’-type Ppd-D1a alleles and homozygous ‘Kitahonami’-type alleles at QDth.kupg-3D and QDth.kupg-4A were similar to that of the plants carrying both homozygous ‘B1-0393’-type alleles at QDth.kupg-3D and QDth.kupg-4A and homozygous ‘Kitahonami’-type alleles at Ppd-D1 (Fig. 4 b). Discussion In this study, we analyzed an early-heading line ‘B1-0393’ from a gamma-irradiated population. However, whole-genome sequencing showed that the distribution of SNP density in ‘B1-0393’ was distinct from that in ‘Kitahonami’ (Fig. 1 c), and its pattern was also considerably different from our previous results on gamma-irradiated ‘Kitahonami’ mutants (Komura et al. 2022 ). Therefore, we speculated that ‘B1-0393’ arose from contamination with a strain of unknown genotype and was not a pure mutant of ‘Kitahonami’, but continued further analyses because it has the potential to be useful for improving the heading time of Japanese wheat cultivars. We used the QTL-seq approach (Takagi et al. 2013 ; Sugihara et al. 2022 ) to identify QTLs associated with heading date using ‘Kitahonami’ and ‘B1-0393’ (Fig. 1 ). QTL-seq mapped three QTLs on chromosomes 2D, 3D, and 4A, designated as QDth.kupg-2D , QDth.kupg-3D , and QDth.kupg-4A , respectively (Fig. 3 ). Through the segregation analysis of F 2 and F 3 plants, we confirmed that ‘B1-0393’ alleles at all QTLs promote heading and that these QTLs have additive effects on heading date (Fig. 4 ). The QDth.kupg-2D region was mapped to the 31–37 Mbp interval on chromosome 2D (Fig. 3 , Table 3 ). This region coincides with the location of the known heading-controlling gene Ppd-D1 . Ppd-D1 is a main regulator of photosensitivity in wheat (Beales et al. 2007 ; Shaw et al. 2013). ‘Kitahonami’ had the same sequences as ‘CS’ except for a 5-bp deletion in exon 7, and ‘B1-0393’ had a 2089-bp deletion in the 5′-UTR region of Ppd-D1 (Supplementary Fig. S3). These mutations are identical to those in the reported photoperiod-sensitive Ppd-D1b allele and photoperiod-insensitive Ppd-D1a allele, respectively (Beales et al. 2007 ; Nishida et al. 2013 ). The deleted 2089-bp region contains several cis -elements, such as the LUX-binding site, CHE motif, and G-box motif, which act as light-response or clock gene–mediated transcription regulators (Pruneda-Paz et al. 2009 ; Helfer et al. 2011 ; Ezer et al. 2017 ). Therefore, the 2089-bp deletion in Ppd-D1 is assumed to disrupt its expression, resulting in early heading (Nishida et al. 2013 ). Our segregation analysis showed that F 2 plants with the ‘B1-0393’-type Ppd-D1a allele headed significantly earlier than those with the ‘Kitahonami’-type Ppd-D1b allele (Fig. 4 a). These results are consistent with previous studies on Ppd-D1 (Beales et al. 2007 ; Nishida et al. 2013 ). Therefore, we consider Ppd-D1 to be one of the causal genes responsible for the difference in heading dates between ‘Kitahonami’ and ‘B1-0393’. QDth.kupg-3D was located on chromosome 3D in the 490–520 Mbp interval (Fig. 3 , Table 3 ). We found eight studies that reported QTLs associated with days to heading on chromosome 3D; their estimated physical positions in the ‘Norin 61’ reference genome are shown in Supplementary Table S7 and Figure S5. Among these QTLs, TaHd066 is close to QDth.kupg-3D , with a peak marker at 533.9 Mbp (Benaouda et al. 2022 ). However, as the markers flanking TaHd066 were located in the 527–536 Mbp interval, TaHd066 would be different from QDth.kupg-3D . The other previously reported QTLs were more than 60 Mbp away from QDth.kupg-3D (Sherman et al. 2014 ; Shukla et al. 2015 ; Benaouda et al. 2022 ). Therefore, comparison with the QTL positions in the ‘Norin 61’ reference genome suggests that QDth.kupg-3D is a novel QTL for heading date. On chromosome 3D, 325 genes are annotated in the 490–520 Mbp interval where QDth.kupg-3D was detected. Among them, five are predicted as candidate heading-related genes in PlantCFG, a database of candidate flowering genes in plants (Liu et al. 2023 ) (Supplementary Tables S8). However, no amino acid substitutions were detected between ‘Kitahonami’ and ‘B1-0393’ in these five genes. On the other hand, eight genes in the QDth.kupg-3D region had missense mutations between ‘Kitahonami’ and ‘B1-0393’ (Supplemental Table S5). One of them, TraesNOR3D01G421200 ( TraesCS3D03G0867000 in CS), annotated as “Calcium sensing receptor, chloroplastic”, had the highest expression levels in shoot, leaf, and spike according to wGRN, a platform for guiding functional gene discovery using integrative gene regulatory networks in wheat (Chen et al. 2023 ) (Supplementary Fig. S6). A missense mutation occurred at position 107 in the protein encoded by TraesNOR3D01G421200 , where threonine in ‘Kitahonami’ was replaced with isoleucine in ‘B1-0393’ (Supplementary Table S5). The gene regulatory network estimated by wGRN suggests that TraesNOR3D01G421200 interacts with known heading regulating genes, such as NUCLEAR FACTOR-YB , CONSTANS-like genes, and the REVEILLE8 clock gene (Nemoto et al. 2003 ; Li et al. 2011 ; Gray et al. 2017 ). In Arabidopsis , a calcium-sensing receptor has been estimated to affect flowering time via photoperiod and the circadian clock (Bonnot et al. 2021 ; Li et al. 2022 ). Therefore, TraesNOR3D01G421200 may be a candidate gene for QDth.kupg-3D . Further studies, including fine mapping of this region and knockout of the candidate gene, will be needed to validate this hypothesis. QDth.kupg-4A was detected in the 610–650 Mbp interval on chromosome 4A (Fig. 3 , Table 3 ). We found that this region overlapped with several reported QTLs for heading date by estimating the physical positions of these QTLs in the ‘Norin 61’ reference genome (Supplementary Table S7, Fig. S5). For instance, QFlt.dms-4A was positioned in the 625–628 Mbp interval on chromosome 4A of ‘Norin 61’ (Zou et al. 2017 ). Sherman et al. ( 2014 ) identified a QTL with the flanking markers Xbarc1158 and wmc262 spanning 612–673 Mbp. A meta-QTL (Hanocq et al. 2007 ) that integrated QTLs reported in four independent studies (Börner et al. 2002 ; Charmet, personal communication; Hanocq et al. 2003 ; Kulwal et al. 2003 ) was positioned in the 598–623 Mbp interval. The barley chromosomal region containing the heading-time QTL eps4L (Laurie et al. 1995 ) is syntenic with the QDth.kupg-4A region. As multiple QTLs involved in heading time have been reported around this region, the causal gene of QDth.kupg-4A is likely to have a stable effect on heading date. On chromosome 4A, 595 genes are annotated in the QDth.kupg-4A region in the ‘Norin 61’ reference genome (Walkowiak et al. 2020 ; Shimizu et al 2021 ). Among them, missense mutations were detected in 38 genes, and 6 genes had deleterious variants, such as frameshift, nonsense, and stop-loss mutations (Supplementary Tables S3, S4). According to PlantCFG (Liu et al. 2023 ), 12 genes are predicted as candidate genes controlling heading (Supplementary Table S8). For instance, TraesNOR4A01G399300 is homologous to Arabidopsis TERMINAL FLOWER 1 , which is a key regulator of flower development (Hanano and Goto 2011 ). Genes annotated as encoding “BTB/POZ domain-containing proteins” form a cluster containing at least five paralogous genes (TraesNOR4A01G381300–1900 in Supplementary Table S8). Among them, a frameshift mutation was detected in TraesNOR4A01G381500 of ‘Kitahonami’ (Supplemental Tables S4, S8). The Arabidopsis BTB/POZ domain–encoding gene LIGHT-RESPONSE BTB1 promotes flowering via photomorphogenesis and vernalization (Christians et al. 2012 ; Hu et al. 2014 ). Although this finding suggests that TraesNOR4A01G381500 may be a candidate gene for QDth.kupg-4A , further genetic mapping and analysis of responsiveness to vernalization and photoperiod of the candidate gene are required to validate it. Ppd-D1 has been used to manipulate heading time, but in the breeding history, alleles of Ppd-D1 have been selected for local adaptation. For example, Langer et al. ( 2014 ) reported that Ppd-D1 is the main gene determining heading date in European winter wheat, with 82% of European cultivars carrying the photoperiod-sensitive allele Ppd-D1b . This biased distribution suggests that Ppd-D1a decreases yield in relatively high-latitude areas by shortening the vegetative phase. The frequency of the photoperiod-insensitive Ppd-D1a allele is low in northern Japan (Seki et al., 2011 ), although Ppd-D1a can accelerate heading in that area (Mizuno et al. 2022 ). This bias in the distribution of Ppd-D1 is probably due to the same reason as in Europe. However, there are cultivars with Ppd-D1a or other insensitive alleles in Europe and northern Japan, which suggest the existence of an unknown genetic mechanism controlling heading time to maintain yield. Although the ‘B1-0393’-type alleles of QDth.kupg-3D and QDth.kupg-4A accelerated heading, the ‘Kitahonami’-type alleles may suppress heading and prolong the vegetative phase. We guessed that the causal genes for these QTLs are the heading-related genes carrying missense and frameshift mutations in ‘Kitahonami’, which would indicate that these are natural variations related to late heading present in the modern cultivars, which may contribute to maintaining yield. Generally speaking, Ppd-1 and Vrn-1 significantly affect heading times, interacting with various genetic backgrounds and environments, making it difficult to fine-tune heading with a few-day accuracy. For example, plants with the photoperiod-insensitive Ppd-D1a allele headed around 3 days earlier than those with the photoperiod-sensitive Ppd-D1b allele in our results (Fig. 4 ). On the other hand, Nishida et al. ( 2013 ) and Mizuno et al. ( 2022 ) reported that the Ppd-D1a allele promotes heading by more than 10 days. This difference is probably due to variations in genetic background or growth environment, and is too large to optimize the heading time. Therefore, fine-tuning heading time for local adaptation is desirable by combining Ppd-1 and Vrn-1 with additional genes with minor effects on heading time (Snape et al. 2001 ; Sheehan and Bentley 2021 ). In our study, QDth.kupg-3D and QDth.kupg-4A had relatively minor effects on heading, and both had additive effects on heading dates independent of the type of Ppd-D1 (Fig. 4 ). Our results suggest that the QTLs reported here and their combination could contribute to the optimization of heading time. Identification of other heading-related genes with minor effects such as WPCL1 , TaELF3 , and WWDR1 (Mizuno et al. 2012 , 2016 , 2023 ; Alvarez et al. 2016 ; Zikhali et al. 2016 ; Hashimoto et al. 2021 ; Komura et al. 2024 ) implyies that the combination of their alleles plays an important role in local adaption that cannot be explained by the major genes. Understanding this genetic mechanism may enable more precise heading control in wheat. Declarations Competing interests The authors declare that they have no competing interests. Funding This study was supported by a grant from the MAFF-commissioned study on “Genomics-Based Technology for Agricultural Improvement [IVG1003]” to YO, HH, and FK, “Smart-breeding System for Innovative Agriculture [DIT1002]” (grant number JP007142) to FK, YO and KY, and JSPS KAKENHI Grant Number 22KJ1943 to SK. Author Contribution YO, HH, and FK generated the early-heading line. FK evaluated the phenotypes. FK generated and SK and KY analyzed genome sequencing data. SK performed genotyping. SK, FK, YI, and KY designed the experiments and genome sequencing analyses. SK, KY, and FK wrote the manuscript. All authors read and approved the final manuscript. Acknowledgement We thank Hiroki Matsuo, from the Graduate School of Agriculture, Kyoto University, for valuable discussions on the manuscript. We also thank Miyuki Oda, Megumi Araki, and Emiko Nakashima (Institute of Crop Science, NARO) for their assistance with data acquisition. Computations were partially performed on the NIG supercomputer at the ROIS National Institute of Genetics, Japan. References Alvarez MA, Tranquilli G, Lewis S, Kippes N, Dubcovsky J (2016) Genetic and physical mapping of the earliness per se locus Eps-Am 1 in Triticum monococcum identifies EARLY FLOWERING 3 ( ELF3 ) as a candidate gene. Funct Integr Genomics 16:365–382 Beales J, Turner A, Griffiths S, Snape JW, Laurie DA (2007) A Pseudo-Response Regulator is misexpressed in the photoperiod insensitive Ppd-D1a mutant of wheat ( Triticum aestivum L.). Theor Appl Genet 115:721–733 Benaouda S, Dadshani S, Koua P, Léon J, Ballvora A (2022) Identification of QTLs for wheat heading time across multiple-environments. Theor Appl Genet 135:2833–2848 Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120 Bonnot T, Blair EJ, Cordingley SJ, Nagel DH (2021) Circadian coordination of cellular processes and abiotic stress responses. Curr Opin Plant Biol 64:102133 Börner A, Schumann E, Fürste A, Cöster H, Leithold B, Röder MS, Weber WE (2002) Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat ( Triticum aestivum L.). Theor Appl Genet 105:921–936 Broman KW, Wu H, Sen Ś, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889–890 Bushnell B (2014) BBMap: a fast, accurate, splice-aware aligner. Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, CA Chen A, Dubcovsky J (2012) Wheat TILLING Mutants Show That the Vernalization Gene VRN1 Down-Regulates the Flowering Repressor VRN2 in Leaves but Is Not Essential for Flowering. PLoS Genet. doi: 10.1371/journal.pgen.1003134 Chen S, Wang J, Deng G, Chen L, Cheng X, Xu H, Zhan K (2018) Interactive effects of multiple vernalization ( Vrn-1 )- and photoperiod ( Ppd-1 )-related genes on the growth habit of bread wheat and their association with heading and flowering time. BMC Plant Biol 18:1–16 Chen Y, Guo Y, Guan P, Wang Y, Wang X, Wang Z, Qin Z, Ma S, Xin M, Hu Z, Yao Y, Ni Z, Sun Q, Guo W, Peng H (2023) A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement. Mol Plant 16:393–414 Christians MJ, Gingerich DJ, Hua Z, Lauer TD, Vierstra RD (2012) The light-response BTB1 and BTB2 proteins assemble nuclear ubiquitin ligases that modify phytochrome B and D signaling in Arabidopsis. Plant Physiol 160:118–134 Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly (Austin) 6:80–92 Díaz A, Zikhali M, Turner AS, Isaac P, Laurie DA (2012) Copy number variation affecting the photoperiod-B1 and vernalization-A1 genes is associated with altered flowering time in wheat ( Triticum aestivum ). PLoS One. doi: 10.1371/journal.pone.0033234 Ezer D, Shepherd SJK, Brestovitsky A, Dickinson P, Cortijo S, Charoensawan V, Box MS, Biswas S, Jaeger KE, Wigge PA (2017) The G-box transcriptional regulatory code in arabidopsis. Plant Physiol 175:628–640 Gray JA, Shalit-Kaneh A, Chu DN, Hsu PY, Harmer SL (2017) The REVEILLE clock genes inhibit growth of juvenile and adult plants by control of cell size. Plant Physiol 173:2308–2322 Hanano S, Goto K (2011) Arabidopsis TERMINAL FLOWER1 is involved in the regulation of flowering time and inflorescence development through transcriptional repression. Plant Cell 23:3172–3184 Hanocq E, Sayers EJ, Niarquin M, Le Gouis J, Charmet G, Gervais L, Dedryver F, Duranton N, Marty N, Dufour P, Rousset M, Worland AJ (2003) A QTL analysis for earliness under field and controlled conditions in a bread wheat doubled-haploid population. In: Proceedings of the 12th EWAC conference, UK, p 57 Hanocq E, Laperche A, Jaminon O, Lainé AL, Le Gouis J (2007) Most significant genome regions involved in the control of earliness traits in bread wheat, as revealed by QTL meta-analysis. Theor Appl Genet 114:569–584 Hartl D (2020) Essential Genetics, 7th edn. Jones & Bartlett Learning, U. S. A. Hashimoto K, Kazama Y, Ichida H, Abe T, Murai K (2021) Einkorn Wheat ( Triticum monococcum ) Mutant Extra-Early Flowering 4 , Generated by Heavy-Ion Beam Irradiation, Has a Deletion of the LIGHT-REGULATED WD1 Homolog. Cytologia (Tokyo) 86:297–302 Helfer A, Nusinow DA, Chow BY, Gehrke AR, Bulyk ML, Kay SA (2011) LUX ARRHYTHMO encodes a nighttime repressor of circadian gene expression in the Arabidopsis core clock. Curr Biol 21:126–133 Hisano H, Sakamoto K, Takagi H, Terauchi R, Sato K (2017) Exome QTL-seq maps monogenic locus and QTLs in barley. BMC Genomics 18:1–9 Hoshino T, Seko H (1996) History of wheat breeding for a half century in Japan. Euphytica 89:215–221 Hu X, Kong X, Wang C, Ma L, Zhao J, Wei J, Zhang X, Loake GJ, Zhang T, Huang J, Yang Y (2014) Proteasome-mediated degradation of FRIGIDA modulates flowering time in Arabidopsis during vernalization. Plant Cell 26:4763–4781 Hyles J, Bloomfield MT, Hunt JR, Trethowan RM, Trevaskis B (2020) Phenology and related traits for wheat adaptation. Heredity (Edinb). doi: 10.1038/s41437-020-0320-1 Illa-Berenguer E, Van Houten J, Huang Z, van der Knaap E (2015) Rapid and reliable identification of tomato fruit weight and locule number loci by QTL-seq. Theor Appl Genet 128:1329–1342 International Wheat Genome Sequencing Consortium (IWGSC) (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. Iwaki K, Haruna S, Niwa T, Kato K (2001) Adaptation and ecological differentiation in wheat with special reference to geographical variation of growth habit and Vrn genotype. Plant Breed 120:107–114 Ji G, Xu Z, Fan X, Zhou Q, Chen L, Yu Q, Liao S, Jiang C, Feng B, Wang T (2023) Identification and validation of major QTL for grain size and weight in bread wheat ( Triticum aestivum L.). Crop J 11:564–572 Jiang Y, Duan L, Guan F, Yao F, Long L, Wang Y, Zhao X, Li H, Li W, Xu Q, Jiang Q, Wang J, Wei Y, Ma J, Kang H, Qi P, Deng M, Zheng Y, Chen G (2022) Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All-Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace “Xiaohemai.” Plant Dis 106:1209–1215 Kato K, Shunji Y (1991) Varietal Variation in Photoperiodic Response, Chilling Requirement and Narrow-Sense Earliness and their Relation to Heading Time in Wheat ( Triticum aesticum L.). Japanese J Breed 17:1460–1462 Kiss T, Balla K, Veisz O, Láng L, Bedő Z, Griffiths S, Isaac P, Karsai I (2014) Allele frequencies in the VRN-A1 , VRN-B1 and VRN-D1 vernalization response and PPD-B1 and PPD-D1 photoperiod sensitivity genes, and their effects on heading in a diverse set of wheat cultivars ( Triticum aestivum L.). Mol Breed 34:297–310 Komura S, Jinno H, Sonoda T, Oono Y, Handa H, Takumi S, Yoshida K, Kobayashi F (2022) Genome sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants. BMC Genomics 23:1–21 Komura S, Yoshida K, Jinno H, Oono Y (2024) Identification of the causal mutation in early heading mutant of bread wheat ( Triticum aestivum L.) using MutMap approach. Mol Breed 44, 41 Kulwal PL, Roy JK, Balyan HS, Gupta PK (2003) QTL mapping for growth and leaf characters in bread wheat. Plant Sci 164:267–277 Langer SM, Longin CFH, Würschum T (2014) Flowering time control in European winter wheat. Front Plant Sci 5:1–11 Laurie DA, Pratchett N, Bezant JH, Snape JW (1995) RFLP mapping of five major genes and eight quantitative trait loci controlling flowering time in a winter × spring barley ( Hordeum vulgare L.) cross. Genome 38:575–585 Li B, Hou L, Song C, Wang Z, Xue Q, Li Y, Qin J, Cao N, Jia C, Zhang Y, Shi W (2022) Biological function of calcium-sensing receptor (CAS) and its coupling calcium signaling in plants. Plant Physiol Biochem 180:74–80 Li C, Distelfeld A, Comis A, Dubcovsky J (2011) Wheat flowering repressor VRN2 and promoter CO2 compete for interactions with NUCLEAR FACTOR-Y complexes. Plant J 67:763–773 Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079 Li Y, Xiong H, Guo H, Zhou C, Xie Y, Zhao L, Gu J, Zhao S (2020) Identification of the vernalization gene VRN-B1 responsible for heading date variation by QTL mapping using a RIL population in wheat. BMC Plant Biol 1–15 Liu D, Li J, Wang S, Huang T, Tao F, Lin Y, Lin W, Zhao X, Huang Y, Jia Y, Yang Z, Luo C, Zhu Q, Sung WK, Wu J, Yang QY (2023) PlantCFG: A comprehensive database with web tools for analyzing candidate flowering genes in multiple plants. Plant Commun 10–13 Ma S, Wang M, Wu J, Guo W, Chen Y, Li G, Wang Y, Shi W, Xia G, Fu D, Kang Z, Ni F (2021) WheatOmics: A platform combining multiple omics data to accelerate functional genomics studies in wheat. Mol Plant 14:1965–1968 Maeda AE, Nakamichi N (2022) Plant clock modifications for adapting flowering time to local environments. Plant Physiol 190:952–967. doi: 10.1093/plphys/kiac107 Mizuno N, Nitta M, Sato K, Nasuda S (2012) A wheat homologue of PHYTOCLOCK 1 is a candidate gene conferring the early heading phenotype to einkorn wheat. Genes Genet Syst 87:357–367 Mizuno N, Kinoshita M, Kinoshita S, Nishida H, Fujita M, Kato K, Murai K, Nasuda S (2016) Loss-of-function mutations in three homoeologous PHYTOCLOCK 1 genes in common wheat are associated with the extra-early flowering phenotype. PLoS One 11:1–19 Mizuno N, Matsunaka H, Yanaka M, Nakata M, Nakamura K, Nakamaru A, Kiribuchi-Otobe C, Ishikawa G, Chono M, Hatta K, Fujita M, Kobayashi F (2022) Allelic variations of Vrn-1 and Ppd-1 genes in Japanese wheat varieties reveal the genotype-environment interaction for heading time. Breed Sci 72:343–354 Mizuno N, Matsunaka H, Yanaka M, Ishikawa G, Kobayashi F, Nakamura K (2023) Natural variations of wheat EARLY FLOWERING 3 highlight their contributions to local adaptation through fine-tuning of heading time. Theor Appl Genet 136:1–15 Nemoto Y, Kisaka M, Fuse T, Yano M, Ogihara Y (2003) Characterization and functional analysis of three wheat genes with homology to the CONSTANS flowering time gene in transgenic rice. Plant J 36:82–93 Ni J, You C, Chen Z, Tang D, Wu H, Deng W, Wang X, Yang J, Bao R, Liu Z, Meng P, Rong T, Liu J (2023) Deploying QTL-seq rapid identification and separation of the major QTLs of tassel branch number for fine-mapping in advanced maize populations. Mol Breed 43, 88. doi: 10.1007/s11032-023-01431-y Nishida H, Yoshida T, Kawakami K, Fujita M, Long B, Akashi Y, Laurie DA, Kato K (2013) Structural variation in the 5′ upstream region of photoperiod-insensitive alleles Ppd-A1a and Ppd-B1a identified in hexaploid wheat ( Triticum aestivum L.), and their effect on heading time. Mol Breed 31:27–37 Ochagavía H, Prieto P, Zikhali M, Griffiths S, Slafer GA (2019) Earliness Per Se by Temperature Interaction on Wheat Development. Sci Rep. doi: 10.1038/s41598-019-39201-6 Pruneda-Paz JL, Breton G, Para A, Kay SA (2009) A functional genomics approach reveals CHE as a component of the Arabidopsis circadian clock. Science 323:1481–1485 Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP (2011) Integrative Genome Viewer. Nat Biotechnol 29:24–6 Seki M, Chono M, Matsunaka H, Fujita M, Oda S, Kubo K, Kiribuchi-Otobe C, Kojima H, Nishida H, Kato K (2011) Distribution of photoperiod-insensitive alleles Ppd-B1a and Ppd-D1a and their effect on heading time in Japanese wheat cultivars. Breed Sci 61:405–412 Seki M, Chono M, Nishimura T, Sato M, Yoshimura Y, Matsunaka H, Fujita M, Oda S, Kubo K, Kiribuchi-Otobe C, Kojima H, Nishida H, Kato K (2013) Distribution of photoperiod-insensitive allele Ppd-A1a and its effect on heading time in Japanese wheat cultivars. Breed Sci 63:309–316 Shaw LM, Turner AS, Laurie DA (2012) The impact of photoperiod insensitive Ppd-1a mutations on the photoperiod pathway across the three genomes of hexaploid wheat ( Triticum aestivum ). Plant J 71:71–84 Sheehan H, Bentley A (2021) Changing times: Opportunities for altering winter wheat phenology. Plants, People, Planet 3:113–123 Sherman JD, Martin JM, Blake NK, Lanning SP, Talbert LE (2014) Genetic basis of agronomic differences between a modern and a historical spring wheat cultivar. Crop Sci 54:1–13 Shimada S, Ogawa T, Kitagawa S, Suzuki T, Ikari C, Shitsukawa N, Abe T, Kawahigashi H, Kikuchi R, Handa H, Murai K (2009) A genetic network of flowering-time genes in wheat leaves, in which an APETALA1/FRUITFULL-like gene, VRN1 , is upstream of FLOWERING LOCUS T . Plant J 58:668–681 Shimizu KK, Copetti D, Okada M, Wicker T, Tameshige T, Hatakeyama M et al. (2021) De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes. Plant Cell Physiol 62:8–27 Shukla S, Singh K, Patil R V, Kadam S, Bharti S, Prasad P, Singh NK, Khanna-Chopra R (2015) Genomic regions associated with grain yield under drought stress in wheat ( Triticum aestivum L.). Euphytica 203:449–467 Singh K, Saini DK, Saripalli G, Batra R, Gautam T, Singh R, Pal S, Kumar M, Jan I, Singh S, Kumar A, Sharma H, Chaudhary J, Kumar K, Kumar S, Singh VK, Singh VP, Kumar D, Sharma S, Kumar S, Kumar R, Sharma S, Gaurav SS, Sharma PK, Balyan HS, Gupta PK (2022) WheatQTLdb V2.0: a supplement to the database for wheat QTL. Mol Breed 42:1–7 Singh V, Sinha P, Obala J, Khan AW, Chitikineni A, Saxena RK, Varshney RK (2022) QTL-seq for the identification of candidate genes for days to flowering and leaf shape in pigeonpea. Heredity (Edinb) 128:411–419 Snape J, Butterworth K, Whitechurch E, Worland AJ (2001) Waiting for Fine Times: Genetics of Flowering Time in Wheat. 67–74 Sugihara Y, Young L, Yaegashi H, Natsume S, Shea DJ, Takagi H, Booker H, Innan H, Terauchi R, Abe A (2022) High-performance pipeline for MutMap and QTL-seq. PeerJ 4–9 Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, Uemura A, Utsushi H, Tamiru M, Takuno S, Innan H, Cano LM, Kamoun S, Terauchi R (2013) QTL-seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J 74:174–183 Walkowiak S, Gao L, Monat C, Haberer G, Kassa MT, Brinton J et al. (2020) Multiple wheat genomes reveal global variation in modern breeding. Nature 588:277–283 Wang Z, Yan L, Chen Y, Wang X, Huai D, Kang Y, Jiang H, Liu K, Lei Y, Liao B (2022) Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut. Theor Appl Genet 135:1779–1795 Wickham H (2016). ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. ISBN 978-3-319–24277–4, https://ggplot2.tidyverse.org . Yan L, Loukoianov A, Tranquilli G, Helguera M, Fahima T, Dubcovsky J (2003) Positional cloning of the wheat vernalization gene VRN1 . Proc Natl Acad Sci U S A 100:6263–6268 Yan L, Helguera M, Kato K, Fukuyama S, Sherman J, Dubcovsky J (2004) Allelic variation at the VRN-1 promoter region in polyploid wheat. Theor Appl Genet 109:1677–1686 Yanagisawa A, Yoshimura Y, Amano Y, Kobayashi S, Nishimura T, Nakamichi K, Araki K, Tanifuji K, Tabiki T, Mikami K, Ikenaga M, Sato N (2007) A New Winter Wheat Variety “Kitahonami.” Rep. Hokkaido Prefect. Agric. Exp. Station. 91:1–13. Zikhali M, Wingen LU, Griffiths S (2016) Delimitation of the Earliness per se D1 ( Eps-D1 ) flowering gene to a subtelomeric chromosomal deletion in bread wheat (Triticum aestivum). J Exp Bot 67:287–299 Zou J, Semagn K, Iqbal M, Chen H, Asif M, N’Diaye A, Navabi A, Perez-Lara E, Pozniak C, Yang RC, Randhawa H, Spaner D (2017) QTLs associated with agronomic traits in the Attila × CDC Go spring wheat population evaluated under conventional management. PLoS One 12:1–20 Additional Declarations No competing interests reported. Supplementary Files Supplementaryfile1.pdf Supplementrayfile2.xlsx Cite Share Download PDF Status: Published Journal Publication published 14 Nov, 2024 Read the published version in Euphytica → Version 1 posted Editorial decision: Revision requested 26 Oct, 2024 Reviews received at journal 26 Oct, 2024 Reviewers agreed at journal 19 Oct, 2024 Reviews received at journal 16 Oct, 2024 Reviewers agreed at journal 27 Sep, 2024 Reviewers invited by journal 27 Sep, 2024 Editor assigned by journal 17 Aug, 2024 Submission checks completed at journal 17 Aug, 2024 First submitted to journal 16 Aug, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4923172","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":344735810,"identity":"dbcbd3ce-c4ca-4912-98db-6baca4d3992a","order_by":0,"name":"Shoya Komura","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABXElEQVRIie2RP0vDQBiH31C5LIfzHZHzKyQEUkul+So5CnWJOLh0KHKl0A7GzoX6IVKEglugkCzpnuLSTl0cMhUHqV7aiv2j4iiYZ7g77t6He+9+ADk5fxDysdDXE2KgCjkHWzWK+FE5NgEHskYq+JcK44I4e8oetHNb0+oNsIuqZ5EXsMwinc+nEI+YrQpEoFGBQn/nGg2PhzQOgT96sUU9qLFS3zUEJCMT40AqYRWU++2XASOXQyoQOHriWgRDaOpPrmwsHXGPXC0IoACUnnOoLMHOFPoKb9yfRLON4shblgeKlinNNii+VDQMiPsJrBpbK0r7QKHe+KHU7BLux+F1+URHZslzjZ4TX5g4DgpnvFvFe28h0d1gIhbnth61BpPnOmJFNZqmaVhmakcoSbqoMGP3x7bjOSJZOquAHGitD2RL2Oh9kxBAIc3GTaY3n/un5KvqnJycnP/DOyQqdjr+CcoMAAAAAElFTkSuQmCC","orcid":"","institution":"Kyoto University","correspondingAuthor":true,"prefix":"","firstName":"Shoya","middleName":"","lastName":"Komura","suffix":""},{"id":344735811,"identity":"ac352270-a789-4cf3-8b0f-04dc7d37078f","order_by":1,"name":"Fuminori Kobayashi","email":"","orcid":"","institution":"National Agriculture and Food Research Organization","correspondingAuthor":false,"prefix":"","firstName":"Fuminori","middleName":"","lastName":"Kobayashi","suffix":""},{"id":344735812,"identity":"ea167ba4-7c76-4ffa-a9c6-ff6aba12277c","order_by":2,"name":"Youko Oono","email":"","orcid":"","institution":"National Agriculture and Food Research Organization","correspondingAuthor":false,"prefix":"","firstName":"Youko","middleName":"","lastName":"Oono","suffix":""},{"id":344735813,"identity":"33e5f194-493c-41bb-8675-91596022331c","order_by":3,"name":"Hirokazu Handa","email":"","orcid":"","institution":"Kyoto Prefectural University","correspondingAuthor":false,"prefix":"","firstName":"Hirokazu","middleName":"","lastName":"Handa","suffix":""},{"id":344735814,"identity":"6bc4a408-71ed-4d6e-aa33-2b341f7dad50","order_by":4,"name":"Yoshihiro Inoue","email":"","orcid":"","institution":"Kyoto University","correspondingAuthor":false,"prefix":"","firstName":"Yoshihiro","middleName":"","lastName":"Inoue","suffix":""},{"id":344735815,"identity":"c96efcd9-f284-417f-b9c7-eedb10219eb5","order_by":5,"name":"Kentaro Yoshida","email":"","orcid":"","institution":"Kyoto University","correspondingAuthor":false,"prefix":"","firstName":"Kentaro","middleName":"","lastName":"Yoshida","suffix":""}],"badges":[],"createdAt":"2024-08-16 07:22:56","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4923172/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4923172/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10681-024-03441-z","type":"published","date":"2024-11-14T15:56:58+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":64630407,"identity":"52d3454c-d5b3-4fe1-8ab1-be65ad52d368","added_by":"auto","created_at":"2024-09-16 20:16:09","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":484163,"visible":true,"origin":"","legend":"\u003cp\u003eComparison of ‘Kitahonami’ and ‘B1-0393’ (a) Phenotypes. (b) Violin plots of days to heading in the 2018–19 and 2019–20 seasons. Dots indicate the data of individuals plants. (c) Distribution and density of SNPs between ‘Kitahonami’ and ‘B1-0393’ over the chromosomes. SNP density was calculated as the number of SNPs per 10 Mbp.\u003c/p\u003e","description":"","filename":"Fig1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/2b4f5a73f2bd3852cca3b9a9.jpg"},{"id":64630684,"identity":"a86e1cec-180f-4b7e-8687-4241fc279593","added_by":"auto","created_at":"2024-09-16 20:24:09","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":137021,"visible":true,"origin":"","legend":"\u003cp\u003eHistograms of days to heading of F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e populations derived from a cross between ‘Kitahonami’ and ‘B1-0393’. The ranges enclosed by the dashed lines show the heading times of the F\u003csub\u003e3\u003c/sub\u003e individuals used for the early-heading bulk (EHB) and late-heading bulk (LHB).\u003c/p\u003e","description":"","filename":"Fig2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/d82f91417d7b1731ccaa604e.jpg"},{"id":64630412,"identity":"8582120f-1b6d-4dfd-b2b4-16e099cc15eb","added_by":"auto","created_at":"2024-09-16 20:16:09","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":677751,"visible":true,"origin":"","legend":"\u003cp\u003eΔSNP-index distribution on chromosomes 2D, 3D, and 4A. Red arrows indicate the \u003cem\u003eQDth.kupg-2D\u003c/em\u003e,\u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e loci. Yellow, green, and blue dots represent SNP-index or ΔSNP-index values. Red lines indicate moving averages. The yellow and green lines in the bottom panels show the 95% and 99% confidential intervals, respectively.\u003c/p\u003e","description":"","filename":"Fig3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/5aef65b71581ef84166beb2a.jpg"},{"id":64630409,"identity":"1096eec2-ac29-4041-91f2-b0b89e34dc3c","added_by":"auto","created_at":"2024-09-16 20:16:09","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":339539,"visible":true,"origin":"","legend":"\u003cp\u003eAssociation between days to heading and the genotypes of \u003cem\u003ePpd-D1\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e. (a) Violin plots showing the results of segregation analysis of the three QTLs for days to heading using an F\u003csub\u003e2\u003c/sub\u003e population generated from a cross between ‘Kitahonami’ and ‘B1-0393.’ K, homozygous for the ‘Kitahonami’ allele; H, heterozygous; B, homozygous for the ‘B1-0393’ allele. Dots indicate days to heading of the individual plants tested. The same letters indicate no significant differences between the means by Tukey–Kramer test (\u003cem\u003eP\u003c/em\u003e \u0026lt; 0.05). (b) Comparison of days to heading among eight genotype combinations of \u003cem\u003ePpd-D1\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e,\u003cem\u003e \u003c/em\u003eand \u003cem\u003eQDth.kupg-4A\u003c/em\u003e in the F\u003csub\u003e3\u003c/sub\u003e population. K, homozygous for ‘Kitahonami’ alleles; B, homozygous for the ‘B1-0393’ alleles. Red points indicate means of days to heading. The same letters indicate no significant differences between means by Steel–Dwass test (\u003cem\u003eP\u003c/em\u003e \u0026lt; 0.05).\u003c/p\u003e","description":"","filename":"Fig4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/25b19b5bd9d9f9ed84bed470.jpg"},{"id":69274767,"identity":"aa2e0af5-07e0-482c-bdfd-b0d4882d8600","added_by":"auto","created_at":"2024-11-18 16:22:26","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2452606,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/ff03cce0-9c02-40b1-94a1-d73185c70925.pdf"},{"id":64630685,"identity":"93ea52bf-8d14-4fba-b6df-177d1be4507a","added_by":"auto","created_at":"2024-09-16 20:24:09","extension":"pdf","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":1228309,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementaryfile1.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/92b7c2916ea204186f8ed905.pdf"},{"id":64631118,"identity":"6f5bd37b-6742-4032-8a33-82b58a3f4b93","added_by":"auto","created_at":"2024-09-16 20:32:09","extension":"xlsx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":31608,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementrayfile2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-4923172/v1/aef733ebf78272317e060af7.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Identification of three QTLs that additively affect heading time in bread wheat (Triticum aestivum L.) by QTL-seq approach","fulltext":[{"header":"Introduction","content":"\u003cp\u003eThe appropriate timing of the transition from vegetative to reproductive growth and subsequent development is crucial for ensuring grain yields and quality. In bread wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.), the optimal heading date is one of the aims in breeding, as it allows growers to avoid winter cold stress while mitigating subsequent heat and drought stresses. Especially, the harvest season is close to the rainy season in Japan. Such conditions can lead to pre-harvest sprouting and \u003cem\u003eFusarium\u003c/em\u003e head blight. Thus, early heading is one of the most essential aims in the breeding program in Japan to avoid losing yield and quality (Hoshino and Seko \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e1996\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eVernalization and photoperiod are major wheat heading components (Kato and Shunji \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e1991\u003c/span\u003e; Hyles et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). In winter wheat, heading requires prolonged exposure to cold, which ensures the transition to reproductive growth in warm conditions. In contrast, plants with spring growth habit can head without exposure to cold. Vernalization is genetically controlled by \u003cem\u003eVERNALIZATION-1\u003c/em\u003e (\u003cem\u003eVrn-1\u003c/em\u003e), \u003cem\u003eVERNALIZATION-2\u003c/em\u003e (\u003cem\u003eVrn-2\u003c/em\u003e), and \u003cem\u003eWheat FLOWERING LOCUS T\u003c/em\u003e (\u003cem\u003eWFT\u003c/em\u003e) (Yan et al. \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; Shimada et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Chen and Dubcovsky \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). The spring growth habit is conferred by recessive loss-of-function alleles of \u003cem\u003eVrn-2\u003c/em\u003e or by at least one dominant allele of \u003cem\u003eVrn-1\u003c/em\u003e (Yan et al. \u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; Hyles et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe sensitivity to the duration of daylight, also known as photoperiod sensitivity, is also a crucial element in deciding the timing of heading. Photoperiod sensitivity is regulated mainly by the \u003cem\u003ePhotoperiod-1\u003c/em\u003e (\u003cem\u003ePpd-1\u003c/em\u003e) gene, which is homologous to \u003cem\u003eArabidopsis thaliana PSEUDO RESPONSE REGULATOR 7\u003c/em\u003e (\u003cem\u003ePRR7\u003c/em\u003e) (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Shaw et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). \u003cem\u003ePpd-1\u003c/em\u003e alleles are classified into photoperiod sensitive (\u003cem\u003ePpd-A1b\u003c/em\u003e, \u003cem\u003ePpd-B1b\u003c/em\u003e, and \u003cem\u003ePpd-D1b\u003c/em\u003e) and photoperiod insensitive (\u003cem\u003ePpd-A1a\u003c/em\u003e, \u003cem\u003ePpd-B1a\u003c/em\u003e, and \u003cem\u003ePpd-D1a\u003c/em\u003e) (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Shaw et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Hyles et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The photoperiod-sensitive alleles show diurnal expression patterns and promote heading under long-day conditions. In contrast, the photoperiod-insensitive alleles are hypermorphic because of deletions in their promoter regions, resulting in dysregulation of \u003cem\u003ePpd-1\u003c/em\u003e expression (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; D\u0026iacute;az et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). The misexpression of \u003cem\u003ePpd-1\u003c/em\u003e induces an increased \u003cem\u003eWFT\u003c/em\u003e expression and promotes heading regardless of daylength (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Shaw et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). In the case of \u003cem\u003ePpd-B1\u003c/em\u003e, photoperiod insensitivity is caused by an increased copy number, which results in high basal expression levels of \u003cem\u003ePpd-B1\u003c/em\u003e (D\u0026iacute;az et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Since \u003cem\u003eVrn-1\u003c/em\u003e and \u003cem\u003ePpd-1\u003c/em\u003e significantly affect heading timing, allelic variations of these genes have been used to adapt plant heading to various environments (Kiss et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Chen et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Mizuno et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The \u003cem\u003ePpd-D1a\u003c/em\u003e allele is more common in wheat cultivars grown in southern Japan than in northern regions (Seki et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2011\u003c/span\u003e), whereas the \u003cem\u003ePpd-A1a\u003c/em\u003e allele is present only in cultivars in northern Japan (Seki et al. \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Iwaki et al. (\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2001\u003c/span\u003e) found that \u003cem\u003eVrn-D1\u003c/em\u003e is more common in South Asia than other spring habitat alleles. Kiss et al. (\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) and Mizuno et al. (\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) revealed that the environment significantly modifies the effects of \u003cem\u003eVrn-1\u003c/em\u003e and \u003cem\u003ePpd-1\u003c/em\u003e, and some alleles alter heading dates only in specific environments. Therefore, an appropriate adoption of \u003cem\u003eVrn-1\u003c/em\u003e and \u003cem\u003ePpd-1\u003c/em\u003e alleles and their combinations is essential to optimize the timing of heading.\u003c/p\u003e \u003cp\u003eEven if the vernalization requirement and day length are satisfied, heading time can vary to some extent. These minor variations are caused by genes referred to as earliness \u003cem\u003eper se\u003c/em\u003e (\u003cem\u003eEps\u003c/em\u003e) genes. Since \u003cem\u003eEps\u003c/em\u003e genes alter heading dates only slightly, they contribute to local adaptation and can be used for fine-tuning (Maeda and Nakamichi \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Although numerous \u003cem\u003eEps\u003c/em\u003e quantitative trait loci (QTLs) have been identified, only a limited number of \u003cem\u003eEps\u003c/em\u003e genes have been cloned (Snape et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2001\u003c/span\u003e; Singh K. et al. \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). In wheat, circadian clock genes \u003cem\u003eWHEAT LUX ARRHYTHMO/PHYTOCLOCK 1\u003c/em\u003e (\u003cem\u003eWPCL1\u003c/em\u003e), \u003cem\u003eEARLY FLOWERING 3\u003c/em\u003e (\u003cem\u003eELF3\u003c/em\u003e), and \u003cem\u003eWHEAT WD REPEAT 1\u003c/em\u003e (\u003cem\u003eWWDR1\u003c/em\u003e) have been identified as \u003cem\u003eEps\u003c/em\u003e genes (Mizuno et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Alvarez et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Zikhali et al. \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Although loss-of-function in one or two homoeologs has a small effect, the null mutants of \u003cem\u003eWPCL1\u003c/em\u003e or \u003cem\u003eWWDR1\u003c/em\u003e show extremely early heading, which results in severe yield reduction (Mizuno et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Hashimoto et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The degree of earliness influenced by \u003cem\u003eELF3\u003c/em\u003e is altered by the prevailing temperature (Ochagav\u0026iacute;a et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Therefore, further exploration of loci that could be used to manipulate heading is important to enable flexible breeding to optimize cultivars to the target geographic region or environment.\u003c/p\u003e \u003cp\u003eThe detection of a QTL is an essential strategy for improving breeding cultivars and identifying causal genes. QTL-seq is a rapid and efficient approach to detect QTLs based on bulked segregant analysis (Takagi et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Using two parents that show opposite extreme phenotypes and pooled DNA bulks from their progeny, QTL-seq detects genomic regions linked to the phenotype of interest on the basis of ∆SNP-index distribution. This strategy has detected QTLs in rice, pigeon pea, tomato, maize, barley, and oilseed rape (Takagi et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Illa-Berenguer et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Hisano et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Singh V. et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Wang et al. \u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Ni et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). In bread wheat, because of its large genome (2\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;6\u003cem\u003ex\u003c/em\u003e\u0026thinsp;=\u0026thinsp;42, ~\u0026thinsp;17 GB), QTL-seq has been used mainly in combination with exome sequencing or RNA sequencing (Li et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Jiang et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Ji et al. \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). These approaches have detected QTLs for agronomically important traits, such as grain size, stripe rust resistance, and heading date. However, with the release of pangenome reference sequences for 10 wheat cultivars (Walkowiak et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) and the reduction in cost of next-generation sequencing, whole-genome sequencing of the bread wheat genome has become more accessible.\u003c/p\u003e \u003cp\u003eIn this study, we aimed to identify candidate genomic regions responsible for the difference in days to heading between the Japanese elite bread wheat \u0026lsquo;Kitahonami\u0026rsquo; and the breeding line \u0026lsquo;B1-0393\u0026rsquo;. Since \u0026lsquo;B1-0393\u0026rsquo; headed 5 days earlier than \u0026lsquo;Kitahonami\u0026rsquo;, the early-heading gene(s) in \u0026lsquo;B1-0393\u0026rsquo; may be used to fine-tune heading date. To identify the causal genomic regions, we performed whole-genome sequencing, conducted QTL-seq, and detected three QTLs associated with heading dates, including the \u003cem\u003ePpd-D1\u003c/em\u003e locus. Segregation tests in the F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e populations confirmed that the three QTLs had additive effects, and the F\u003csub\u003e3\u003c/sub\u003e plants harboring triple early-heading alleles had the earliest heading dates. This study provides valuable insights into QTL mapping in wheat to control heading date.\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003ePlant materials\u003c/h2\u003e \u003cp\u003e\u0026lsquo;B1-0393\u0026rsquo; was found in a gamma-irradiated wheat population and heads earlier than \u0026lsquo;Kitahonami\u0026rsquo;, whichwhi was released by the Kitami Agricultural Experiment Station, Hokkaido (Yanagisawa et al. \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). The heading dates of \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; were evaluated in the experimental field of the Institute of Crop Science, NARO (Tsukuba, Japan; lat. 36.1\u0026deg; N, long. 140.6\u0026deg; E), under natural field conditions in the 2018\u0026ndash;19 (sown 24 October 2018) and 2019\u0026ndash;20 (sown 31 October 2019) cropping seasons. \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; were crossed, generating F\u003csub\u003e2\u003c/sub\u003e (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;130) and F\u003csub\u003e3\u003c/sub\u003e (\u003cem\u003en\u003c/em\u003e\u0026thinsp;=\u0026thinsp;598) populations. To measure heading time, \u0026lsquo;Kitahonami\u0026rsquo;, \u0026lsquo;B1-0393\u0026rsquo;, and their F\u003csub\u003e2\u003c/sub\u003e progeny were sown on 27 October 2020 (2020\u0026ndash;21 season), and their F\u003csub\u003e3\u003c/sub\u003e progeny were sown on 29 October 2021 (2021\u0026ndash;22 season) in the same field under natural conditions. The first appearance of the spike tip was recorded as the heading date in both seasons.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eGenome sequencing and bioinformatics analysis\u003c/h2\u003e \u003cp\u003eTotal DNA was extracted from a leaf of \u0026lsquo;B1-0393\u0026rsquo; using a DNeasy Plant Mini Kit (Qiagen, Hilden, Germany). For QTL-seq, 20 individuals with the earliest heading dates and 20 with the latest heading dates were selected from 598 F\u003csub\u003e3\u003c/sub\u003e progeny. The total DNA of each individual was extracted as above. and bulked to produce an early heading bulk and a late heading bulk. DNA libraries of \u0026lsquo;B1-0393\u0026rsquo; and each bulk were constructed using a TruSeq DNA library Preparation Kit (Illumina, San Diego, CA, USA). Genome sequencing was performed on an Illumina NovaSeq in paired-end 150-bp mode. The sequencing reads of the wild-type \u0026lsquo;Kitahonami\u0026rsquo; were obtained from our previous study (Komura et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe quality of sequence reads was checked with FASTQC v. 0.11.7 and the reads were filtered by Trimmomatic v. 0.33 (Bolger et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) with an average minimum Phred quality score per 4 bp of \u0026lt;\u0026thinsp;20 and a length of \u0026lt;\u0026thinsp;40 bp. We used BWA v. 0.7.17 with default options to align the filtered paired-end reads to the reference genomes of bread wheat cultivars \u0026lsquo;Norin 61\u0026rsquo; and \u0026lsquo;Chinese Spring\u0026rsquo; (\u0026lsquo;CS\u0026rsquo;) (IWGSC 2018; Walkowiak et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Shimizu et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). PCR duplications were removed using SAMtools v. 1.9 (Li et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). The average number and genome coverage of aligned reads were calculated using BBMap v. 37.77 (Bushnell \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). SNPs and indels were called using BCFtools v. 1.9 (Li et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2009\u003c/span\u003e) with options \u0026ldquo;-A\u0026rdquo; and \u0026ldquo;-Q 0\u0026rdquo;. Variants differing between the parents were extracted using BCFtools v. 1.9. We removed variants that did not meet the criteria of (1) mapping quality of reads containing variants of \u0026ge;\u0026thinsp;40, (2) a homozygous variant detection site, (3)\u0026thinsp;\u0026ge;\u0026thinsp;5 aligned reads with the altered genotype at the variant detection site but no reads of the reference allele in one parent, and (4)\u0026thinsp;\u0026ge;\u0026thinsp;5 reads with the reference allele but no variant-type read counts in the other parent. The effects of SNPs and indels were predicted in SnpEff v. 4.3t (Cingolani et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). To calculate SNP density, the number of SNPs per 1 Mbp was counted in Python v. 3.9.0. The R package \u0026lsquo;ggplot2\u0026rsquo; (Wickham, \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) was used to visualize the distribution of the moving average of SNP density over the chromosomes. QTL-seq analysis (Takagi et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Sugihara et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) was performed according to the pipeline of Sugihara et al. (\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The reference genome of \u0026lsquo;Norin 61\u0026rsquo; was used for QTL-seq analysis. The ΔSNP-index of the two bulks were estimated at the positions where SNPs were detected between the \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;. To narrow down causal variants in the candidate regions, indel-index was also estimated. The SNP-index and indel-index were calculated by dividing the number of reads supporting the alternative sequence by the total number of reads aligned to the position. The ΔSNP-index was calculated as the difference between the SNP-index of the early-heading bulk and that of the late-heading bulk.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eMarker development and genotyping\u003c/h2\u003e \u003cp\u003eFor segregation analysis of QTL regions, we designed cleaved amplified polymorphic sequence (CAPS) markers on the basis of SNPs between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;. For CAPS markers, we designed subgenome-specific primers by BLAST searches against the \u0026lsquo;Norin 61\u0026rsquo; reference sequence through WheatOmics (Ma et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The markers used in this study are listed in Supplementary Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e. \u0026lsquo;Kitahonami\u0026rsquo;, \u0026lsquo;B1-0393\u0026rsquo;, and their 130 F\u003csub\u003e2\u003c/sub\u003e and 598 F\u003csub\u003e3\u003c/sub\u003e progeny were used for genotyping. Total DNA was extracted as above and PCR was performed using Quick-Taq HS DyeMix (Toyobo, Osaka, Japan) as follows: pre-denaturation at 94\u0026deg;C for 2 min; 35 cycles of denaturation at 94\u0026deg;C for 30 s, annealing at 62\u0026deg;C for 30 s, and extension at 68\u0026deg;C for 30 s; and post-extension at 68\u0026deg;C for 1 min. The PCR products were digested with the restriction enzymes listed in Supplementary Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e and electrophoresed in 2% agarose gel. To confirm the sequences of \u003cem\u003ePpd-D1\u003c/em\u003e in \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;, we checked the read alignment using Integrative Genomics Viewer (Robinson et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). \u003cem\u003ePpd-D1\u003c/em\u003e genotyping was performed using allele-specific primers for \u003cem\u003ePpd-D1a\u003c/em\u003e or \u003cem\u003ePpd-D1b\u003c/em\u003e, as described by Nishida et al. (\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eStatistical analysis\u003c/h2\u003e \u003cp\u003eThe normality of the distributions in the F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e populations was confirmed using the Shapiro\u0026ndash;Wilk test. Variance among genotypes was assessed by an analysis of variance (ANOVA) or the Kruskal\u0026ndash;Wallis test. For multiple comparisons, the Tukey\u0026ndash;Kramer test or the Steel\u0026ndash;Dwass test was used.\u003c/p\u003e \u003cp\u003eTo detect epistatic interactions between three QTLs, we used ANOVA analysis with the following general linear model (GLM):\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equa\" name=\"EquationSource\"\u003e\n$$\\:{GLM}_{iJl}=\\mu\\:+{QTL}_{1i}+{QTL}_{2j}+{QTL}_{1}{QTL}_{2ij}+{e}_{iil}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003e\u0026micro;\u003c/em\u003e is the population mean value; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{QTL}_{1i}\\:\\)\u003c/span\u003e\u003c/span\u003eand \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{QTL}_{2j}\\)\u003c/span\u003e\u003c/span\u003e indicate the effects of the \u003cem\u003ei\u003c/em\u003eth and \u003cem\u003ej\u003c/em\u003eth QTLs; and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{e}_{iil}\\)\u003c/span\u003e\u003c/span\u003e is experimental error.\u003c/p\u003e \u003cp\u003eThe additive and dominant effects and phenotypic variation explained by each locus were determined using the R package \u0026lsquo;R/qtl\u0026rsquo; (Broman et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2003\u003c/span\u003e) with the following model:\u003cdiv id=\"Equb\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equb\" name=\"EquationSource\"\u003e\n$$\\:{HD}_{ijl}=\\mu\\:+{Ppd\\text{-}D1}_{i}+{QDth.kupg\\text{-}3D}_{j}+{QDth.kupg\\text{-}4A}_{l}+{e}_{ijl}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\:{Ppd\\text{-}D1}_{i}\\)\u003c/span\u003e\u003c/span\u003e, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{QDth.kupg\\text{-}3D}_{j}\\)\u003c/span\u003e\u003c/span\u003e, and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{QDth.kupg\\text{-}4A}_{l}\\)\u003c/span\u003e\u003c/span\u003e indicate the effects of the \u003cem\u003ePpd-D1\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e loci; and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{e}_{ijl}\\)\u003c/span\u003e\u003c/span\u003e is experimental error.\u003c/p\u003e \u003cp\u003eAdditive effect (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:a\\)\u003c/span\u003e\u003c/span\u003e) and dominance effect (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:d\\)\u003c/span\u003e\u003c/span\u003e) were calculated as follows:\u003cdiv id=\"Equc\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equc\" name=\"EquationSource\"\u003e\n$$\\:a=\\frac{({A}_{1}{A}_{1}-{A}_{2}{A}_{2})}{2}$$\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equd\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equd\" name=\"EquationSource\"\u003e\n$$\\:d={A}_{1}{A}_{2}-\\frac{({A}_{1}{A}_{1}+{A}_{2}{A}_{2})}{2}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{A}_{1}{A}_{1}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{A}_{2}{A}_{2}\\)\u003c/span\u003e\u003c/span\u003e are the mean values of plants homozygous for the \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;, respectively, alleles; and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{A}_{1}{A}_{2}\\)\u003c/span\u003e\u003c/span\u003e is the mean of heterozygous plants.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eEstimation of heritability\u003c/h2\u003e \u003cp\u003ePhenotypic variance in the progeny (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{V}_{F}\\)\u003c/span\u003e\u003c/span\u003e) includes genetic and environmental variances, whereas that in the parents (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{V}_{P1}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{V}_{P2}\\)\u003c/span\u003e\u003c/span\u003e) includes only environmental variances. Therefore, the broad-sense heritability (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{H}^{2}\\)\u003c/span\u003e\u003c/span\u003e) was estimated as (Hartl \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2020\u003c/span\u003e):\u003cdiv id=\"Eque\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Eque\" name=\"EquationSource\"\u003e\n$$\\:{H}^{2}=\\frac{{V}_{F}-\\raisebox{1ex}{$1$}\\!\\left/\\:\\!\\raisebox{-1ex}{$2$}\\right.({V}_{P1}+{V}_{P2})}{{V}_{F}}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eThe narrow-sense heritability (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{h}^{2}\\)\u003c/span\u003e\u003c/span\u003e) was estimated as (Hartl \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2020\u003c/span\u003e):\u003cdiv id=\"Equf\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equf\" name=\"EquationSource\"\u003e\n$$\\:{M}_{F3\\_late}^{\u0026rdquo;}={M}_{F2}+{h}^{2}({M}_{F2\\_late}^{{\\prime\\:}}-{M}_{F2})$$\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equg\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equg\" name=\"EquationSource\"\u003e\n$$\\:{M}_{F3\\_early}^{\u0026rdquo;}={M}_{F2}-{h}^{2}({M}_{F2}-{M}_{F2\\_early}^{{\\prime\\:}})$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eBy rearranging the above equations, the following equation can be obtained:\u003cdiv id=\"Equh\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equh\" name=\"EquationSource\"\u003e\n$$\\:{h}^{2}=\\frac{({M}_{F3\\_late}^{\u0026rdquo;}-{M}_{F3\\_early}^{\u0026rdquo;})}{({M}_{F2\\_late}^{{\\prime\\:}}-{M}_{F2\\_early}^{{\\prime\\:}})}$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{M}_{F2}\\)\u003c/span\u003e\u003c/span\u003e is the mean of the F\u003csub\u003e2\u003c/sub\u003e population; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{M}_{F2\\_late}^{{\\prime\\:}}\\)\u003c/span\u003e\u003c/span\u003e is the mean of F\u003csub\u003e2\u003c/sub\u003e individuals that headed later than a threshold; \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{M}_{F2\\_early}^{{\\prime\\:}}\\)\u003c/span\u003e\u003c/span\u003e is the mean of F\u003csub\u003e2\u003c/sub\u003e individuals that headed earlier than a threshold; and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{M}_{F3\\_late}^{\u0026rdquo;}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{M}_{F3\\_early}^{\u0026rdquo;}\\)\u003c/span\u003e\u003c/span\u003e are the means of F\u003csub\u003e3\u003c/sub\u003e individuals generated from the selected F\u003csub\u003e2\u003c/sub\u003e individuals. We set the thresholds for selecting F\u003csub\u003e2\u003c/sub\u003e individuals as \u0026gt;\u0026thinsp;179 days to heading for late heading and \u0026lt;\u0026thinsp;174 days to heading for early heading.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eCharacterization of early-heading line \u0026lsquo;B1-0393\u0026rsquo;\u003c/h2\u003e \u003cp\u003e\u0026lsquo;B1-0393\u0026rsquo; headed 5 days earlier than \u0026lsquo;Kitahonami\u0026rsquo; in the 2018\u0026ndash;19 season and 3 days earlier in the 2019\u0026ndash;20 season (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). To characterize nucleotide variations between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;, we resequenced the whole genome of \u0026lsquo;B1-0393\u0026rsquo; with short reads and used the sequencing reads of \u0026lsquo;Kitahonami\u0026rsquo; from our previous study (Komura et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). After quality control, 2.7\u0026nbsp;billion reads of \u0026lsquo;Kitahonami\u0026rsquo; and 2.1\u0026nbsp;billion reads of \u0026lsquo;B1-0393\u0026rsquo; were filtered, and 2.3 and 1.8\u0026nbsp;billion reads, respectively, were aligned to the \u0026lsquo;Norin 61\u0026rsquo; reference genome (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). These reads covered 96.3% and 96.1%, respectively, of the reference genome, with an average depth-of-coverage of 15.08 and 18.56. We detected 255,369 indels and 4,618,884 SNPs between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; (Supplementary Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). To evaluate genetic divergence between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;, we calculated SNP density per 10 Mbp (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). The SNPs were unevenly distributed on the chromosomes, with interspersed regions of high and low SNP density. For example, SNP density was high in the 520\u0026ndash;640 Mbp interval of chromosome 4A, but it was low in the 0\u0026ndash;520 Mbp and 640\u0026ndash;720 Mbp intervals, suggesting that the genetic background was similar in these regions. Fewer SNPs were detected in the D genome than in the other genomes, suggesting a closer genetic background of the D genome (Supplementary Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). These results indicate that the genetic background of \u0026lsquo;B1-0393\u0026rsquo; is similar to that of \u0026lsquo;Kitahonami\u0026rsquo; in some regions but is divergent in others (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eSummary of alignments and coverages of genome sequencing of \u0026lsquo;Kitahonami\u0026rsquo; and the early-heading line \u0026lsquo;B1-0393\u0026rsquo;.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eName\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eTotal filtered reads\u0026nbsp;\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eAligned reads after removal of low-quality reads\u0026nbsp;\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eReference bases covered\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eAverage depth-of-coverage\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eKitahonami\u0026nbsp;\u003csup\u003ec\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,724,430,844\u003c/p\u003e \u003cp\u003e(93.63%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e2,334,502,744\u003c/p\u003e \u003cp\u003e(85.69%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e96.32%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e15.08\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB1-0393\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,151,447,818\u003c/p\u003e \u003cp\u003e(93.26%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e1,846,489,772\u003c/p\u003e \u003cp\u003e(85.83%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e96.10%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e18.56\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eEarly heading bulk\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,477,029,056\u003c/p\u003e \u003cp\u003e(85.58%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e2,148,834,659\u003c/p\u003e \u003cp\u003e(86.75%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e96.71%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e21.47\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLate heading bulk\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2,530,545,270\u003c/p\u003e \u003cp\u003e(92.51%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e2,246,339,991\u003c/p\u003e \u003cp\u003e(88.77%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e96.83%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e21.66\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003ea\u003c/sup\u003e Rate of filtered reads\u0026thinsp;=\u0026thinsp;Total filtered reads / Total raw reads \u0026times; 100.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003eb\u003c/sup\u003e Rate of filtered reads\u0026thinsp;=\u0026thinsp;Aligned reads after removal of low-quality reads / Total filtered reads \u0026times; 100.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"5\"\u003e\u003csup\u003ec\u003c/sup\u003e \u0026lsquo;Kitahonami\u0026rsquo; sequencing data are reproduced from Komura et al. (\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe F\u003csub\u003e2\u003c/sub\u003e population showed a Gaussian distribution of days to heading ranging from 169 to 185, while \u0026lsquo;Kitahonami\u0026rsquo; headed at 179 days and \u0026lsquo;B1-0393\u0026rsquo; at 174 days (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). In the F\u003csub\u003e3\u003c/sub\u003e population, progeny headed from 175 to 191 days, while \u0026lsquo;Kitahonami\u0026rsquo; headed at 180 days and \u0026lsquo;B1-0393\u0026rsquo; at 177 days. The distribution of days to heading was skewed towards early heading and showed transgressive segregation (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), suggesting that multiple loci determine the timing of heading. The broad-sense heritability was 0.65 in the F\u003csub\u003e2\u003c/sub\u003e population and 0.57 in the F\u003csub\u003e3\u003c/sub\u003e population, while the narrow-sense heritability was 0.29 (Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The former includes additive and dominant variances, while the latter includes only additive variances. The narrow-sense heritability was almost half of the broad-sense heritability, the difference indicating the degree of dominant variances. Therefore, the difference in heading time between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; was inferred to be a combination of dominance and additive effects.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eDays to heading in the F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e populations.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"8\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePopulation\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNo. of plants\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eMean\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eSD\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSkewness\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eKurtosis\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eBroad-sense heritability\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c8\"\u003e \u003cp\u003eNarrow-sense heritability\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF\u003csub\u003e2\u003c/sub\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e130\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e176.13\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e2.73\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.24\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.23\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e0.65\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c8\" morerows=\"1\" rowspan=\"2\"\u003e \u003cp\u003e0.29\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eF\u003csub\u003e3\u003c/sub\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e596\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e179.47\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1.22\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e3.91\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e0.57\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eDetection of QTLs for heading date by QTL-seq\u003c/h2\u003e \u003cp\u003eFrom the F\u003csub\u003e3\u003c/sub\u003e population, we selected 20 plants that headed on 176 or 177 days as the early-heading bulk and 20 plants that headed on 185\u0026ndash;189 days as the late-heading bulk, and resequenced the whole genomes of the two bulks. After quality control, the filtered 2.5\u0026nbsp;billion reads per bulk were aligned to the \u0026lsquo;Norin 61\u0026rsquo; reference sequence. Of the filtered reads, 2.1\u0026nbsp;billion reads from the early-heading bulk and 2.2\u0026nbsp;billion reads from the late-heading bulk could be aligned (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). The aligned reads of each bulk covered 96.71% (21.47\u0026times; average depth-of-coverage) and 96.83% (21.66\u0026times;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:)\\:\\)\u003c/span\u003e\u003c/span\u003eof the reference sequences. QTL-seq analysis revealed that the moving average of ΔSNP-index in the regions of 31\u0026ndash;37 Mbp on chromosome 2D, 490\u0026ndash;520 Mbp on chromosome 3D, and 610\u0026ndash;650 Mbp on chromosome 4A exceeded the 99% or 95% confidence intervals (CIs) (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Supplementary Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e). These QTL regions were designated as \u003cem\u003eQDth.kupg-2D\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eQDth.kupg-2D\u003c/em\u003e houses the known heading-related gene \u003cem\u003ePpd-D1\u003c/em\u003e. The alignment of the reads to the \u0026lsquo;Norin 61\u0026rsquo; reference sequence revealed that \u0026lsquo;Kitahonami\u0026rsquo; has a 5-bp deletion in exon 7 of \u003cem\u003ePpd-D1\u003c/em\u003e, resulting in a premature stop codon and a truncated protein. \u0026lsquo;Norin 61\u0026rsquo; has a photoperiod-insensitive allele \u003cem\u003ePpd-D1\u003c/em\u003e. Therefore, to investigate the alleles of \u003cem\u003ePpd-D1\u003c/em\u003e in \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;, we used a reference sequence of \u0026lsquo;CS\u0026rsquo;, which carries the photoperiod-sensitive allele \u003cem\u003ePpd-D1b\u003c/em\u003e. Read alignment to the \u0026lsquo;CS\u0026rsquo; reference sequence showed that \u0026lsquo;B1-0393\u0026rsquo; has a 2089-bp deletion in the 5\u0026prime; upstream region and confirmed the 5-bp deletion in exon 7 in \u0026lsquo;Kitahonami\u0026rsquo;. In addition, a 16-bp deletion in exon 8 was identified in both parents (Supplementary Fig. S3). The \u003cem\u003ePpd-D1\u003c/em\u003e allele in \u0026lsquo;B1-0393\u0026rsquo; was identical to that in \u0026lsquo;Winter-Abukumawase\u0026rsquo;, which is reportedly the photoperiod-insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e allele (Haplotype \u003cem\u003ePpd-D1a.1\u003c/em\u003e) (Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). The \u003cem\u003ePpd-D1\u003c/em\u003e allele in \u0026lsquo;Kitahonami\u0026rsquo; was the same as that in \u0026lsquo;Norstar\u0026rsquo; and \u0026lsquo;Chihokukomugi\u0026rsquo;, which is referred to as the photoperiod-sensitive \u003cem\u003ePpd-D1b\u003c/em\u003e allele (Haplotype \u003cem\u003ePpd-D1b.2\u003c/em\u003e) (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). According to Nishida et al. (\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), plants with the \u003cem\u003ePpd-D1a\u003c/em\u003e allele head earlier than plants with the \u003cem\u003ePpd-D1b\u003c/em\u003e allele. Therefore, we assumed \u003cem\u003ePpd-D1\u003c/em\u003e to be the causal gene in \u003cem\u003eQDth.kupg-2D\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eIn the \u003cem\u003eQDth.kupg-4A\u003c/em\u003e region, we found 9518 variants with the SNP-index and indel-index exceeding the 95% CIs. Among them, 6 variants caused frameshift, nonsense, or stop-loss mutations, and 54 caused missense mutations (Supplementary Tables S3, S4). In the \u003cem\u003eQDth.kupg-3D\u003c/em\u003e region, we found 1504 variants with the SNP-index and indel-index exceeding the 95% CIs. Three variants caused missense mutations (Supplementary Table S5). No known genes related to wheat heading timing were found in the \u003cem\u003eQDth.kupg-3D\u003c/em\u003e or \u003cem\u003eQDth.kupg-4A\u003c/em\u003e regions.\u003c/p\u003e \u003cp\u003e \u003cb\u003eValidation of\u003c/b\u003e \u003cb\u003eQDth.kupg-2D\u003c/b\u003e, \u003cb\u003eQDth.kupg-3D\u003c/b\u003e, \u003cb\u003eand\u003c/b\u003e \u003cb\u003eQDth.kupg-4A\u003c/b\u003e\u003c/p\u003e \u003cp\u003eTo confirm the association between the detected QTLs and heading date, we conducted a segregation analysis using 130 F\u003csub\u003e2\u003c/sub\u003e plants. Days to heading differed significantly depending on the genotypes of \u003cem\u003ePpd-D1\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea, Supplementary Table S6). Plants with the \u0026lsquo;B1-0393\u0026rsquo; alleles at all three QTLs headed earlier. The phenotypic variation explained (PVE) was 10.25% by \u003cem\u003ePpd-D1\u003c/em\u003e, 2% by \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and 8.24% by \u003cem\u003eQDth.kupg-4A\u003c/em\u003e (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). The dominance effect of \u003cem\u003ePpd-D1\u003c/em\u003e was 0.06, indicating that \u003cem\u003ePpd-D1a\u003c/em\u003e and \u003cem\u003ePpd-D1b\u003c/em\u003e have similar influence on heading, and heterozygous individuals had phenotypes intermediate between those of individuals homozygous for the parental alleles (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). The additive effects were 0.51 for \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and 1.01 for \u003cem\u003eQDth.kupg-4A\u003c/em\u003e, and dominance effects were 0.21 and 0.75, respectively (Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). The segregation analysis indicated incomplete dominance by \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and a dominant phenotype of the \u0026lsquo;Kitahonami\u0026rsquo; allele of \u003cem\u003eQDth.kupg-4A\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). The PVE and dominant effect values indicated that \u003cem\u003ePpd-D1\u003c/em\u003e had the largest and \u003cem\u003eQDth.kupg-3D\u003c/em\u003e had the smallest effect on heading date.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eQuantitative trait loci (QTLs) for days to heading detected by QTL-seq in F\u003csub\u003e2\u003c/sub\u003e population.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLocus\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eChr.\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eQTL range (Mbp)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eMarker position (bp)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePVE\u003csup\u003ec\u003c/sup\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eAdditive effect\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003eDominance effect\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eQDth.kupg-2D\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e2D\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e31\u0026ndash;37\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e34,671,330\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e10.25%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e1.35\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e0.06\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eQDth.kupg-3D\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e3D\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e490\u0026ndash;520\u003csup\u003ea\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e496,096,772\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e2.00%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.51\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e0.21\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003eQDth.kupg-4A\u003c/em\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e4A\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e610\u0026ndash;650\u003csup\u003eb\u003c/sup\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e632,986,650\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e8.24%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e1.01\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e0.75\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"7\"\u003e\u003csup\u003ea\u003c/sup\u003e ∆SNP-index\u0026thinsp;\u0026gt;\u0026thinsp;95% confidence interval.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"7\"\u003e\u003csup\u003eb\u003c/sup\u003e ∆SNP-index\u0026thinsp;\u0026gt;\u0026thinsp;99% confidence interval.\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"7\"\u003e\u003csup\u003ec\u003c/sup\u003e Phenotypic variance explained.\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eSince no significant interactions between these QTLs were found in ANOVA, we expected these QTLs to function in an additive manner (Supplementary Table S6). To investigate the additive effects of the three QTLs on heading date, we determined the genotypes of 598 F\u003csub\u003e3\u003c/sub\u003e plants. The plants carrying the three \u0026lsquo;B1-0393\u0026rsquo;-type alleles headed an average of 3.19 days earlier than the plants carrying three \u0026lsquo;Kitahonami\u0026rsquo;-type alleles (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb). The plants carrying the \u0026lsquo;B1-0393\u0026rsquo;-type allele at \u003cem\u003eQDth.kupg-3D\u003c/em\u003e or \u003cem\u003eQDth.kupg-4A\u003c/em\u003e with the \u0026lsquo;B1-0393\u0026rsquo;-type \u003cem\u003ePpd-D1a\u003c/em\u003e allele also headed earlier (by 2.62 and 2.92 days on average) than those with three \u0026lsquo;Kitahonami\u0026rsquo;-type alleles. The effect of the \u0026lsquo;B1-0393\u0026rsquo;-type allele was strongest for \u003cem\u003ePpd-D1\u003c/em\u003e and weakest for \u003cem\u003eQDth.kupg-3D\u003c/em\u003e. The effect on heading date of plants carrying homozygous \u0026lsquo;B1-0393\u0026rsquo;-type \u003cem\u003ePpd-D1a\u003c/em\u003e alleles and homozygous \u0026lsquo;Kitahonami\u0026rsquo;-type alleles at \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e were similar to that of the plants carrying both homozygous \u0026lsquo;B1-0393\u0026rsquo;-type alleles at \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e and homozygous \u0026lsquo;Kitahonami\u0026rsquo;-type alleles at \u003cem\u003ePpd-D1\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb).\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn this study, we analyzed an early-heading line \u0026lsquo;B1-0393\u0026rsquo; from a gamma-irradiated population. However, whole-genome sequencing showed that the distribution of SNP density in \u0026lsquo;B1-0393\u0026rsquo; was distinct from that in \u0026lsquo;Kitahonami\u0026rsquo; (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec), and its pattern was also considerably different from our previous results on gamma-irradiated \u0026lsquo;Kitahonami\u0026rsquo; mutants (Komura et al. \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Therefore, we speculated that \u0026lsquo;B1-0393\u0026rsquo; arose from contamination with a strain of unknown genotype and was not a pure mutant of \u0026lsquo;Kitahonami\u0026rsquo;, but continued further analyses because it has the potential to be useful for improving the heading time of Japanese wheat cultivars. We used the QTL-seq approach (Takagi et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Sugihara et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) to identify QTLs associated with heading date using \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). QTL-seq mapped three QTLs on chromosomes 2D, 3D, and 4A, designated as \u003cem\u003eQDth.kupg-2D\u003c/em\u003e, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e, respectively (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Through the segregation analysis of F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e plants, we confirmed that \u0026lsquo;B1-0393\u0026rsquo; alleles at all QTLs promote heading and that these QTLs have additive effects on heading date (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eQDth.kupg-2D\u003c/em\u003e region was mapped to the 31\u0026ndash;37 Mbp interval on chromosome 2D (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). This region coincides with the location of the known heading-controlling gene \u003cem\u003ePpd-D1\u003c/em\u003e. \u003cem\u003ePpd-D1\u003c/em\u003e is a main regulator of photosensitivity in wheat (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Shaw et al. 2013). \u0026lsquo;Kitahonami\u0026rsquo; had the same sequences as \u0026lsquo;CS\u0026rsquo; except for a 5-bp deletion in exon 7, and \u0026lsquo;B1-0393\u0026rsquo; had a 2089-bp deletion in the 5\u0026prime;-UTR region of \u003cem\u003ePpd-D1\u003c/em\u003e (Supplementary Fig. S3). These mutations are identical to those in the reported photoperiod-sensitive \u003cem\u003ePpd-D1b\u003c/em\u003e allele and photoperiod-insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e allele, respectively (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). The deleted 2089-bp region contains several \u003cem\u003ecis\u003c/em\u003e-elements, such as the LUX-binding site, CHE motif, and G-box motif, which act as light-response or clock gene\u0026ndash;mediated transcription regulators (Pruneda-Paz et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Helfer et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Ezer et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Therefore, the 2089-bp deletion in \u003cem\u003ePpd-D1\u003c/em\u003e is assumed to disrupt its expression, resulting in early heading (Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Our segregation analysis showed that F\u003csub\u003e2\u003c/sub\u003e plants with the \u0026lsquo;B1-0393\u0026rsquo;-type \u003cem\u003ePpd-D1a\u003c/em\u003e allele headed significantly earlier than those with the \u0026lsquo;Kitahonami\u0026rsquo;-type \u003cem\u003ePpd-D1b\u003c/em\u003e allele (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). These results are consistent with previous studies on \u003cem\u003ePpd-D1\u003c/em\u003e (Beales et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Nishida et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Therefore, we consider \u003cem\u003ePpd-D1\u003c/em\u003e to be one of the causal genes responsible for the difference in heading dates between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQDth.kupg-3D\u003c/em\u003e was located on chromosome 3D in the 490\u0026ndash;520 Mbp interval (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). We found eight studies that reported QTLs associated with days to heading on chromosome 3D; their estimated physical positions in the \u0026lsquo;Norin 61\u0026rsquo; reference genome are shown in Supplementary Table S7 and Figure S5. Among these QTLs, \u003cem\u003eTaHd066\u003c/em\u003e is close to \u003cem\u003eQDth.kupg-3D\u003c/em\u003e, with a peak marker at 533.9 Mbp (Benaouda et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). However, as the markers flanking \u003cem\u003eTaHd066\u003c/em\u003e were located in the 527\u0026ndash;536 Mbp interval, \u003cem\u003eTaHd066\u003c/em\u003e would be different from \u003cem\u003eQDth.kupg-3D\u003c/em\u003e. The other previously reported QTLs were more than 60 Mbp away from \u003cem\u003eQDth.kupg-3D\u003c/em\u003e (Sherman et al. \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Shukla et al. \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Benaouda et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Therefore, comparison with the QTL positions in the \u0026lsquo;Norin 61\u0026rsquo; reference genome suggests that \u003cem\u003eQDth.kupg-3D\u003c/em\u003e is a novel QTL for heading date.\u003c/p\u003e \u003cp\u003eOn chromosome 3D, 325 genes are annotated in the 490\u0026ndash;520 Mbp interval where \u003cem\u003eQDth.kupg-3D\u003c/em\u003e was detected. Among them, five are predicted as candidate heading-related genes in PlantCFG, a database of candidate flowering genes in plants (Liu et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) (Supplementary Tables S8). However, no amino acid substitutions were detected between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; in these five genes. On the other hand, eight genes in the \u003cem\u003eQDth.kupg-3D\u003c/em\u003e region had missense mutations between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo; (Supplemental Table S5). One of them, \u003cem\u003eTraesNOR3D01G421200\u003c/em\u003e (\u003cem\u003eTraesCS3D03G0867000\u003c/em\u003e in CS), annotated as \u0026ldquo;Calcium sensing receptor, chloroplastic\u0026rdquo;, had the highest expression levels in shoot, leaf, and spike according to wGRN, a platform for guiding functional gene discovery using integrative gene regulatory networks in wheat (Chen et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) (Supplementary Fig. S6). A missense mutation occurred at position 107 in the protein encoded by \u003cem\u003eTraesNOR3D01G421200\u003c/em\u003e, where threonine in \u0026lsquo;Kitahonami\u0026rsquo; was replaced with isoleucine in \u0026lsquo;B1-0393\u0026rsquo; (Supplementary Table S5). The gene regulatory network estimated by wGRN suggests that \u003cem\u003eTraesNOR3D01G421200\u003c/em\u003e interacts with known heading regulating genes, such as \u003cem\u003eNUCLEAR FACTOR-YB\u003c/em\u003e, \u003cem\u003eCONSTANS-like\u003c/em\u003e genes, and the \u003cem\u003eREVEILLE8\u003c/em\u003e clock gene (Nemoto et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; Li et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Gray et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). In \u003cem\u003eArabidopsis\u003c/em\u003e, a calcium-sensing receptor has been estimated to affect flowering time via photoperiod and the circadian clock (Bonnot et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Li et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Therefore, \u003cem\u003eTraesNOR3D01G421200\u003c/em\u003e may be a candidate gene for \u003cem\u003eQDth.kupg-3D\u003c/em\u003e. Further studies, including fine mapping of this region and knockout of the candidate gene, will be needed to validate this hypothesis.\u003c/p\u003e \u003cp\u003e \u003cem\u003eQDth.kupg-4A\u003c/em\u003e was detected in the 610\u0026ndash;650 Mbp interval on chromosome 4A (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e, Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). We found that this region overlapped with several reported QTLs for heading date by estimating the physical positions of these QTLs in the \u0026lsquo;Norin 61\u0026rsquo; reference genome (Supplementary Table S7, Fig. S5). For instance, \u003cem\u003eQFlt.dms-4A\u003c/em\u003e was positioned in the 625\u0026ndash;628 Mbp interval on chromosome 4A of \u0026lsquo;Norin 61\u0026rsquo; (Zou et al. \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Sherman et al. (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) identified a QTL with the flanking markers \u003cem\u003eXbarc1158\u003c/em\u003e and \u003cem\u003ewmc262\u003c/em\u003e spanning 612\u0026ndash;673 Mbp. A meta-QTL (Hanocq et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2007\u003c/span\u003e) that integrated QTLs reported in four independent studies (B\u0026ouml;rner et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Charmet, personal communication; Hanocq et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; Kulwal et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2003\u003c/span\u003e) was positioned in the 598\u0026ndash;623 Mbp interval. The barley chromosomal region containing the heading-time QTL \u003cem\u003eeps4L\u003c/em\u003e (Laurie et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e1995\u003c/span\u003e) is syntenic with the \u003cem\u003eQDth.kupg-4A\u003c/em\u003e region. As multiple QTLs involved in heading time have been reported around this region, the causal gene of \u003cem\u003eQDth.kupg-4A\u003c/em\u003e is likely to have a stable effect on heading date.\u003c/p\u003e \u003cp\u003eOn chromosome 4A, 595 genes are annotated in the \u003cem\u003eQDth.kupg-4A\u003c/em\u003e region in the \u0026lsquo;Norin 61\u0026rsquo; reference genome (Walkowiak et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Shimizu et al \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Among them, missense mutations were detected in 38 genes, and 6 genes had deleterious variants, such as frameshift, nonsense, and stop-loss mutations (Supplementary Tables S3, S4). According to PlantCFG (Liu et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), 12 genes are predicted as candidate genes controlling heading (Supplementary Table S8). For instance, \u003cem\u003eTraesNOR4A01G399300\u003c/em\u003e is homologous to \u003cem\u003eArabidopsis TERMINAL FLOWER 1\u003c/em\u003e, which is a key regulator of flower development (Hanano and Goto \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Genes annotated as encoding \u0026ldquo;BTB/POZ domain-containing proteins\u0026rdquo; form a cluster containing at least five paralogous genes (TraesNOR4A01G381300\u0026ndash;1900 in Supplementary Table S8). Among them, a frameshift mutation was detected in TraesNOR4A01G381500 of \u0026lsquo;Kitahonami\u0026rsquo; (Supplemental Tables S4, S8). The \u003cem\u003eArabidopsis\u003c/em\u003e BTB/POZ domain\u0026ndash;encoding gene \u003cem\u003eLIGHT-RESPONSE BTB1\u003c/em\u003e promotes flowering via photomorphogenesis and vernalization (Christians et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Hu et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Although this finding suggests that TraesNOR4A01G381500 may be a candidate gene for \u003cem\u003eQDth.kupg-4A\u003c/em\u003e, further genetic mapping and analysis of responsiveness to vernalization and photoperiod of the candidate gene are required to validate it.\u003c/p\u003e \u003cp\u003e \u003cem\u003ePpd-D1\u003c/em\u003e has been used to manipulate heading time, but in the breeding history, alleles of \u003cem\u003ePpd-D1\u003c/em\u003e have been selected for local adaptation. For example, Langer et al. (\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) reported that \u003cem\u003ePpd-D1\u003c/em\u003e is the main gene determining heading date in European winter wheat, with 82% of European cultivars carrying the photoperiod-sensitive allele \u003cem\u003ePpd-D1b\u003c/em\u003e. This biased distribution suggests that \u003cem\u003ePpd-D1a\u003c/em\u003e decreases yield in relatively high-latitude areas by shortening the vegetative phase. The frequency of the photoperiod-insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e allele is low in northern Japan (Seki et al., \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2011\u003c/span\u003e), although \u003cem\u003ePpd-D1a\u003c/em\u003e can accelerate heading in that area (Mizuno et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). This bias in the distribution of \u003cem\u003ePpd-D1\u003c/em\u003e is probably due to the same reason as in Europe. However, there are cultivars with \u003cem\u003ePpd-D1a\u003c/em\u003e or other insensitive alleles in Europe and northern Japan, which suggest the existence of an unknown genetic mechanism controlling heading time to maintain yield. Although the \u0026lsquo;B1-0393\u0026rsquo;-type alleles of \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e accelerated heading, the \u0026lsquo;Kitahonami\u0026rsquo;-type alleles may suppress heading and prolong the vegetative phase. We guessed that the causal genes for these QTLs are the heading-related genes carrying missense and frameshift mutations in \u0026lsquo;Kitahonami\u0026rsquo;, which would indicate that these are natural variations related to late heading present in the modern cultivars, which may contribute to maintaining yield.\u003c/p\u003e \u003cp\u003eGenerally speaking, \u003cem\u003ePpd-1\u003c/em\u003e and \u003cem\u003eVrn-1\u003c/em\u003e significantly affect heading times, interacting with various genetic backgrounds and environments, making it difficult to fine-tune heading with a few-day accuracy. For example, plants with the photoperiod-insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e allele headed around 3 days earlier than those with the photoperiod-sensitive \u003cem\u003ePpd-D1b\u003c/em\u003e allele in our results (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). On the other hand, Nishida et al. (\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) and Mizuno et al. (\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) reported that the \u003cem\u003ePpd-D1a\u003c/em\u003e allele promotes heading by more than 10 days. This difference is probably due to variations in genetic background or growth environment, and is too large to optimize the heading time. Therefore, fine-tuning heading time for local adaptation is desirable by combining \u003cem\u003ePpd-1\u003c/em\u003e and \u003cem\u003eVrn-1\u003c/em\u003e with additional genes with minor effects on heading time (Snape et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2001\u003c/span\u003e; Sheehan and Bentley \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). In our study, \u003cem\u003eQDth.kupg-3D\u003c/em\u003e and \u003cem\u003eQDth.kupg-4A\u003c/em\u003e had relatively minor effects on heading, and both had additive effects on heading dates independent of the type of \u003cem\u003ePpd-D1\u003c/em\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Our results suggest that the QTLs reported here and their combination could contribute to the optimization of heading time. Identification of other heading-related genes with minor effects such as \u003cem\u003eWPCL1\u003c/em\u003e, \u003cem\u003eTaELF3\u003c/em\u003e, and \u003cem\u003eWWDR1\u003c/em\u003e (Mizuno et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2012\u003c/span\u003e, \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2016\u003c/span\u003e, \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2023\u003c/span\u003e; Alvarez et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Zikhali et al. \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Hashimoto et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Komura et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2024\u003c/span\u003e) implyies that the combination of their alleles plays an important role in local adaption that cannot be explained by the major genes. Understanding this genetic mechanism may enable more precise heading control in wheat.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e \u003ch2\u003eCompeting interests\u003c/h2\u003e \u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eThis study was supported by a grant from the MAFF-commissioned study on \u0026ldquo;Genomics-Based Technology for Agricultural Improvement [IVG1003]\u0026rdquo; to YO, HH, and FK, \u0026ldquo;Smart-breeding System for Innovative Agriculture [DIT1002]\u0026rdquo; (grant number JP007142) to FK, YO and KY, and JSPS KAKENHI Grant Number 22KJ1943 to SK.\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eYO, HH, and FK generated the early-heading line. FK evaluated the phenotypes. FK generated and SK and KY analyzed genome sequencing data. SK performed genotyping. SK, FK, YI, and KY designed the experiments and genome sequencing analyses. SK, KY, and FK wrote the manuscript. All authors read and approved the final manuscript.\u003c/p\u003e\u003ch2\u003eAcknowledgement\u003c/h2\u003e\u003cp\u003eWe thank Hiroki Matsuo, from the Graduate School of Agriculture, Kyoto University, for valuable discussions on the manuscript. We also thank Miyuki Oda, Megumi Araki, and Emiko Nakashima (Institute of Crop Science, NARO) for their assistance with data acquisition. Computations were partially performed on the NIG supercomputer at the ROIS National Institute of Genetics, Japan.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eAlvarez MA, Tranquilli G, Lewis S, Kippes N, Dubcovsky J (2016) Genetic and physical mapping of the earliness \u003cem\u003eper se\u003c/em\u003e locus \u003cem\u003eEps-Am 1\u003c/em\u003e in \u003cem\u003eTriticum monococcum\u003c/em\u003e identifies \u003cem\u003eEARLY FLOWERING 3\u003c/em\u003e (\u003cem\u003eELF3\u003c/em\u003e) as a candidate gene. Funct Integr Genomics 16:365\u0026ndash;382\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBeales J, Turner A, Griffiths S, Snape JW, Laurie DA (2007) A Pseudo-Response Regulator is misexpressed in the photoperiod insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e mutant of wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). Theor Appl Genet 115:721\u0026ndash;733\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBenaouda S, Dadshani S, Koua P, L\u0026eacute;on J, Ballvora A (2022) Identification of QTLs for wheat heading time across multiple-environments. Theor Appl Genet 135:2833\u0026ndash;2848\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBolger AM, Lohse M, Usadel B (2014) Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30:2114\u0026ndash;2120\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBonnot T, Blair EJ, Cordingley SJ, Nagel DH (2021) Circadian coordination of cellular processes and abiotic stress responses. Curr Opin Plant Biol 64:102133\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eB\u0026ouml;rner A, Schumann E, F\u0026uuml;rste A, C\u0026ouml;ster H, Leithold B, R\u0026ouml;der MS, Weber WE (2002) Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). Theor Appl Genet 105:921\u0026ndash;936\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBroman KW, Wu H, Sen Ś, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889\u0026ndash;890\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBushnell B (2014) BBMap: a fast, accurate, splice-aware aligner. Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, CA\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen A, Dubcovsky J (2012) Wheat TILLING Mutants Show That the Vernalization Gene \u003cem\u003eVRN1\u003c/em\u003e Down-Regulates the Flowering Repressor \u003cem\u003eVRN2\u003c/em\u003e in Leaves but Is Not Essential for Flowering. PLoS Genet. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pgen.1003134\u003c/span\u003e\u003cspan address=\"10.1371/journal.pgen.1003134\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen S, Wang J, Deng G, Chen L, Cheng X, Xu H, Zhan K (2018) Interactive effects of multiple vernalization (\u003cem\u003eVrn-1\u003c/em\u003e)- and photoperiod (\u003cem\u003ePpd-1\u003c/em\u003e)-related genes on the growth habit of bread wheat and their association with heading and flowering time. BMC Plant Biol 18:1\u0026ndash;16\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen Y, Guo Y, Guan P, Wang Y, Wang X, Wang Z, Qin Z, Ma S, Xin M, Hu Z, Yao Y, Ni Z, Sun Q, Guo W, Peng H (2023) A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement. Mol Plant 16:393\u0026ndash;414\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChristians MJ, Gingerich DJ, Hua Z, Lauer TD, Vierstra RD (2012) The light-response \u003cem\u003eBTB1\u003c/em\u003e and \u003cem\u003eBTB2\u003c/em\u003e proteins assemble nuclear ubiquitin ligases that modify phytochrome B and D signaling in Arabidopsis. Plant Physiol 160:118\u0026ndash;134\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly (Austin) 6:80\u0026ndash;92\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eD\u0026iacute;az A, Zikhali M, Turner AS, Isaac P, Laurie DA (2012) Copy number variation affecting the \u003cem\u003ephotoperiod-B1\u003c/em\u003e and \u003cem\u003evernalization-A1\u003c/em\u003e genes is associated with altered flowering time in wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e). PLoS One. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0033234\u003c/span\u003e\u003cspan address=\"10.1371/journal.pone.0033234\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEzer D, Shepherd SJK, Brestovitsky A, Dickinson P, Cortijo S, Charoensawan V, Box MS, Biswas S, Jaeger KE, Wigge PA (2017) The G-box transcriptional regulatory code in arabidopsis. Plant Physiol 175:628\u0026ndash;640\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGray JA, Shalit-Kaneh A, Chu DN, Hsu PY, Harmer SL (2017) The \u003cem\u003eREVEILLE\u003c/em\u003e clock genes inhibit growth of juvenile and adult plants by control of cell size. Plant Physiol 173:2308\u0026ndash;2322\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHanano S, Goto K (2011) \u003cem\u003eArabidopsis\u003c/em\u003e TERMINAL FLOWER1 is involved in the regulation of flowering time and inflorescence development through transcriptional repression. Plant Cell 23:3172\u0026ndash;3184\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHanocq E, Sayers EJ, Niarquin M, Le Gouis J, Charmet G, Gervais L, Dedryver F, Duranton N, Marty N, Dufour P, Rousset M, Worland AJ (2003) A QTL analysis for earliness under field and controlled conditions in a bread wheat doubled-haploid population. In: Proceedings of the 12th EWAC conference, UK, p 57\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHanocq E, Laperche A, Jaminon O, Lain\u0026eacute; AL, Le Gouis J (2007) Most significant genome regions involved in the control of earliness traits in bread wheat, as revealed by QTL meta-analysis. Theor Appl Genet 114:569\u0026ndash;584\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHartl D (2020) Essential Genetics, 7th edn. Jones \u0026amp; Bartlett Learning, U. S. A.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHashimoto K, Kazama Y, Ichida H, Abe T, Murai K (2021) Einkorn Wheat (\u003cem\u003eTriticum monococcum\u003c/em\u003e) Mutant \u003cem\u003eExtra-Early Flowering 4\u003c/em\u003e, Generated by Heavy-Ion Beam Irradiation, Has a Deletion of the \u003cem\u003eLIGHT-REGULATED WD1\u003c/em\u003e Homolog. Cytologia (Tokyo) 86:297\u0026ndash;302\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHelfer A, Nusinow DA, Chow BY, Gehrke AR, Bulyk ML, Kay SA (2011) \u003cem\u003eLUX ARRHYTHMO\u003c/em\u003e encodes a nighttime repressor of circadian gene expression in the \u003cem\u003eArabidopsis\u003c/em\u003e core clock. Curr Biol 21:126\u0026ndash;133\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHisano H, Sakamoto K, Takagi H, Terauchi R, Sato K (2017) Exome QTL-seq maps monogenic locus and QTLs in barley. BMC Genomics 18:1\u0026ndash;9\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHoshino T, Seko H (1996) History of wheat breeding for a half century in Japan. Euphytica 89:215\u0026ndash;221\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHu X, Kong X, Wang C, Ma L, Zhao J, Wei J, Zhang X, Loake GJ, Zhang T, Huang J, Yang Y (2014) Proteasome-mediated degradation of \u003cem\u003eFRIGIDA\u003c/em\u003e modulates flowering time in \u003cem\u003eArabidopsis\u003c/em\u003e during vernalization. Plant Cell 26:4763\u0026ndash;4781\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHyles J, Bloomfield MT, Hunt JR, Trethowan RM, Trevaskis B (2020) Phenology and related traits for wheat adaptation. Heredity (Edinb). doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41437-020-0320-1\u003c/span\u003e\u003cspan address=\"10.1038/s41437-020-0320-1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIlla-Berenguer E, Van Houten J, Huang Z, van der Knaap E (2015) Rapid and reliable identification of tomato fruit weight and locule number loci by QTL-seq. Theor Appl Genet 128:1329\u0026ndash;1342\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eInternational Wheat Genome Sequencing Consortium (IWGSC) (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIwaki K, Haruna S, Niwa T, Kato K (2001) Adaptation and ecological differentiation in wheat with special reference to geographical variation of growth habit and \u003cem\u003eVrn\u003c/em\u003e genotype. Plant Breed 120:107\u0026ndash;114\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJi G, Xu Z, Fan X, Zhou Q, Chen L, Yu Q, Liao S, Jiang C, Feng B, Wang T (2023) Identification and validation of major QTL for grain size and weight in bread wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). Crop J 11:564\u0026ndash;572\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJiang Y, Duan L, Guan F, Yao F, Long L, Wang Y, Zhao X, Li H, Li W, Xu Q, Jiang Q, Wang J, Wei Y, Ma J, Kang H, Qi P, Deng M, Zheng Y, Chen G (2022) Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All-Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace \u0026ldquo;Xiaohemai.\u0026rdquo; Plant Dis 106:1209\u0026ndash;1215\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKato K, Shunji Y (1991) Varietal Variation in Photoperiodic Response, Chilling Requirement and Narrow-Sense Earliness and their Relation to Heading Time in Wheat (\u003cem\u003eTriticum aesticum\u003c/em\u003e L.). Japanese J Breed 17:1460\u0026ndash;1462\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKiss T, Balla K, Veisz O, L\u0026aacute;ng L, Bedő Z, Griffiths S, Isaac P, Karsai I (2014) Allele frequencies in the \u003cem\u003eVRN-A1\u003c/em\u003e, \u003cem\u003eVRN-B1\u003c/em\u003e and \u003cem\u003eVRN-D1\u003c/em\u003e vernalization response and \u003cem\u003ePPD-B1\u003c/em\u003e and \u003cem\u003ePPD-D1\u003c/em\u003e photoperiod sensitivity genes, and their effects on heading in a diverse set of wheat cultivars (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). Mol Breed 34:297\u0026ndash;310\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKomura S, Jinno H, Sonoda T, Oono Y, Handa H, Takumi S, Yoshida K, Kobayashi F (2022) Genome sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants. BMC Genomics 23:1\u0026ndash;21\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKomura S, Yoshida K, Jinno H, Oono Y (2024) Identification of the causal mutation in early heading mutant of bread wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.) using MutMap approach. Mol Breed 44, 41\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKulwal PL, Roy JK, Balyan HS, Gupta PK (2003) QTL mapping for growth and leaf characters in bread wheat. Plant Sci 164:267\u0026ndash;277\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLanger SM, Longin CFH, W\u0026uuml;rschum T (2014) Flowering time control in European winter wheat. Front Plant Sci 5:1\u0026ndash;11\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLaurie DA, Pratchett N, Bezant JH, Snape JW (1995) RFLP mapping of five major genes and eight quantitative trait loci controlling flowering time in a winter \u0026times; spring barley (\u003cem\u003eHordeum vulgare\u003c/em\u003e L.) cross. Genome 38:575\u0026ndash;585\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi B, Hou L, Song C, Wang Z, Xue Q, Li Y, Qin J, Cao N, Jia C, Zhang Y, Shi W (2022) Biological function of calcium-sensing receptor (CAS) and its coupling calcium signaling in plants. Plant Physiol Biochem 180:74\u0026ndash;80\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi C, Distelfeld A, Comis A, Dubcovsky J (2011) Wheat flowering repressor \u003cem\u003eVRN2\u003c/em\u003e and promoter \u003cem\u003eCO2\u003c/em\u003e compete for interactions with \u003cem\u003eNUCLEAR FACTOR-Y\u003c/em\u003e complexes. Plant J 67:763\u0026ndash;773\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078\u0026ndash;2079\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Y, Xiong H, Guo H, Zhou C, Xie Y, Zhao L, Gu J, Zhao S (2020) Identification of the vernalization gene VRN-B1 responsible for heading date variation by QTL mapping using a RIL population in wheat. BMC Plant Biol 1\u0026ndash;15\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu D, Li J, Wang S, Huang T, Tao F, Lin Y, Lin W, Zhao X, Huang Y, Jia Y, Yang Z, Luo C, Zhu Q, Sung WK, Wu J, Yang QY (2023) PlantCFG: A comprehensive database with web tools for analyzing candidate flowering genes in multiple plants. Plant Commun 10\u0026ndash;13\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa S, Wang M, Wu J, Guo W, Chen Y, Li G, Wang Y, Shi W, Xia G, Fu D, Kang Z, Ni F (2021) WheatOmics: A platform combining multiple omics data to accelerate functional genomics studies in wheat. Mol Plant 14:1965\u0026ndash;1968\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaeda AE, Nakamichi N (2022) Plant clock modifications for adapting flowering time to local environments. Plant Physiol 190:952\u0026ndash;967. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1093/plphys/kiac107\u003c/span\u003e\u003cspan address=\"10.1093/plphys/kiac107\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMizuno N, Nitta M, Sato K, Nasuda S (2012) A wheat homologue of \u003cem\u003ePHYTOCLOCK 1\u003c/em\u003e is a candidate gene conferring the early heading phenotype to einkorn wheat. Genes Genet Syst 87:357\u0026ndash;367\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMizuno N, Kinoshita M, Kinoshita S, Nishida H, Fujita M, Kato K, Murai K, Nasuda S (2016) Loss-of-function mutations in three homoeologous \u003cem\u003ePHYTOCLOCK 1\u003c/em\u003e genes in common wheat are associated with the extra-early flowering phenotype. PLoS One 11:1\u0026ndash;19\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMizuno N, Matsunaka H, Yanaka M, Nakata M, Nakamura K, Nakamaru A, Kiribuchi-Otobe C, Ishikawa G, Chono M, Hatta K, Fujita M, Kobayashi F (2022) Allelic variations of \u003cem\u003eVrn-1\u003c/em\u003e and \u003cem\u003ePpd-1\u003c/em\u003e genes in Japanese wheat varieties reveal the genotype-environment interaction for heading time. Breed Sci 72:343\u0026ndash;354\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMizuno N, Matsunaka H, Yanaka M, Ishikawa G, Kobayashi F, Nakamura K (2023) Natural variations of wheat \u003cem\u003eEARLY FLOWERING 3\u003c/em\u003e highlight their contributions to local adaptation through fine-tuning of heading time. Theor Appl Genet 136:1\u0026ndash;15\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNemoto Y, Kisaka M, Fuse T, Yano M, Ogihara Y (2003) Characterization and functional analysis of three wheat genes with homology to the \u003cem\u003eCONSTANS\u003c/em\u003e flowering time gene in transgenic rice. Plant J 36:82\u0026ndash;93\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNi J, You C, Chen Z, Tang D, Wu H, Deng W, Wang X, Yang J, Bao R, Liu Z, Meng P, Rong T, Liu J (2023) Deploying QTL-seq rapid identification and separation of the major QTLs of tassel branch number for fine-mapping in advanced maize populations. Mol Breed 43, 88. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s11032-023-01431-y\u003c/span\u003e\u003cspan address=\"10.1007/s11032-023-01431-y\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNishida H, Yoshida T, Kawakami K, Fujita M, Long B, Akashi Y, Laurie DA, Kato K (2013) Structural variation in the 5\u0026prime; upstream region of photoperiod-insensitive alleles \u003cem\u003ePpd-A1a\u003c/em\u003e and \u003cem\u003ePpd-B1a\u003c/em\u003e identified in hexaploid wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.), and their effect on heading time. Mol Breed 31:27\u0026ndash;37\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOchagav\u0026iacute;a H, Prieto P, Zikhali M, Griffiths S, Slafer GA (2019) Earliness \u003cem\u003ePer Se\u003c/em\u003e by Temperature Interaction on Wheat Development. Sci Rep. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41598-019-39201-6\u003c/span\u003e\u003cspan address=\"10.1038/s41598-019-39201-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePruneda-Paz JL, Breton G, Para A, Kay SA (2009) A functional genomics approach reveals CHE as a component of the \u003cem\u003eArabidopsis\u003c/em\u003e circadian clock. Science 323:1481\u0026ndash;1485\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRobinson JT, Thorvaldsd\u0026oacute;ttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP (2011) Integrative Genome Viewer. Nat Biotechnol 29:24\u0026ndash;6\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSeki M, Chono M, Matsunaka H, Fujita M, Oda S, Kubo K, Kiribuchi-Otobe C, Kojima H, Nishida H, Kato K (2011) Distribution of photoperiod-insensitive alleles \u003cem\u003ePpd-B1a\u003c/em\u003e and \u003cem\u003ePpd-D1a\u003c/em\u003e and their effect on heading time in Japanese wheat cultivars. Breed Sci 61:405\u0026ndash;412\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSeki M, Chono M, Nishimura T, Sato M, Yoshimura Y, Matsunaka H, Fujita M, Oda S, Kubo K, Kiribuchi-Otobe C, Kojima H, Nishida H, Kato K (2013) Distribution of photoperiod-insensitive allele \u003cem\u003ePpd-A1a\u003c/em\u003e and its effect on heading time in Japanese wheat cultivars. Breed Sci 63:309\u0026ndash;316\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShaw LM, Turner AS, Laurie DA (2012) The impact of photoperiod insensitive \u003cem\u003ePpd-1a\u003c/em\u003e mutations on the photoperiod pathway across the three genomes of hexaploid wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e). Plant J 71:71\u0026ndash;84\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSheehan H, Bentley A (2021) Changing times: Opportunities for altering winter wheat phenology. Plants, People, Planet 3:113\u0026ndash;123\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSherman JD, Martin JM, Blake NK, Lanning SP, Talbert LE (2014) Genetic basis of agronomic differences between a modern and a historical spring wheat cultivar. Crop Sci 54:1\u0026ndash;13\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShimada S, Ogawa T, Kitagawa S, Suzuki T, Ikari C, Shitsukawa N, Abe T, Kawahigashi H, Kikuchi R, Handa H, Murai K (2009) A genetic network of flowering-time genes in wheat leaves, in which an \u003cem\u003eAPETALA1/FRUITFULL-like\u003c/em\u003e gene, \u003cem\u003eVRN1\u003c/em\u003e, is upstream of \u003cem\u003eFLOWERING LOCUS T\u003c/em\u003e. Plant J 58:668\u0026ndash;681\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShimizu KK, Copetti D, Okada M, Wicker T, Tameshige T, Hatakeyama M et al. (2021) \u003cem\u003eDe Novo\u003c/em\u003e Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes. Plant Cell Physiol 62:8\u0026ndash;27\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eShukla S, Singh K, Patil R V, Kadam S, Bharti S, Prasad P, Singh NK, Khanna-Chopra R (2015) Genomic regions associated with grain yield under drought stress in wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). Euphytica 203:449\u0026ndash;467\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSingh K, Saini DK, Saripalli G, Batra R, Gautam T, Singh R, Pal S, Kumar M, Jan I, Singh S, Kumar A, Sharma H, Chaudhary J, Kumar K, Kumar S, Singh VK, Singh VP, Kumar D, Sharma S, Kumar S, Kumar R, Sharma S, Gaurav SS, Sharma PK, Balyan HS, Gupta PK (2022) WheatQTLdb V2.0: a supplement to the database for wheat QTL. Mol Breed 42:1\u0026ndash;7\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSingh V, Sinha P, Obala J, Khan AW, Chitikineni A, Saxena RK, Varshney RK (2022) QTL-seq for the identification of candidate genes for days to flowering and leaf shape in pigeonpea. Heredity (Edinb) 128:411\u0026ndash;419\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSnape J, Butterworth K, Whitechurch E, Worland AJ (2001) Waiting for Fine Times: Genetics of Flowering Time in Wheat. 67\u0026ndash;74\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSugihara Y, Young L, Yaegashi H, Natsume S, Shea DJ, Takagi H, Booker H, Innan H, Terauchi R, Abe A (2022) High-performance pipeline for MutMap and QTL-seq. PeerJ 4\u0026ndash;9\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTakagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, Uemura A, Utsushi H, Tamiru M, Takuno S, Innan H, Cano LM, Kamoun S, Terauchi R (2013) QTL-seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J 74:174\u0026ndash;183\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWalkowiak S, Gao L, Monat C, Haberer G, Kassa MT, Brinton J et al. (2020) Multiple wheat genomes reveal global variation in modern breeding. Nature 588:277\u0026ndash;283\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang Z, Yan L, Chen Y, Wang X, Huai D, Kang Y, Jiang H, Liu K, Lei Y, Liao B (2022) Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut. Theor Appl Genet 135:1779\u0026ndash;1795\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWickham H (2016). ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. ISBN 978-3-319\u0026ndash;24277\u0026ndash;4, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://ggplot2.tidyverse.org\u003c/span\u003e\u003cspan address=\"https://ggplot2.tidyverse.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYan L, Loukoianov A, Tranquilli G, Helguera M, Fahima T, Dubcovsky J (2003) Positional cloning of the wheat vernalization gene \u003cem\u003eVRN1\u003c/em\u003e. Proc Natl Acad Sci U S A 100:6263\u0026ndash;6268\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYan L, Helguera M, Kato K, Fukuyama S, Sherman J, Dubcovsky J (2004) Allelic variation at the \u003cem\u003eVRN-1\u003c/em\u003e promoter region in polyploid wheat. Theor Appl Genet 109:1677\u0026ndash;1686\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYanagisawa A, Yoshimura Y, Amano Y, Kobayashi S, Nishimura T, Nakamichi K, Araki K, Tanifuji K, Tabiki T, Mikami K, Ikenaga M, Sato N (2007) A New Winter Wheat Variety \u0026ldquo;Kitahonami.\u0026rdquo; Rep. Hokkaido Prefect. Agric. Exp. Station. 91:1\u0026ndash;13.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZikhali M, Wingen LU, Griffiths S (2016) Delimitation of the \u003cem\u003eEarliness per se D1\u003c/em\u003e (\u003cem\u003eEps-D1\u003c/em\u003e) flowering gene to a subtelomeric chromosomal deletion in bread wheat (Triticum aestivum). J Exp Bot 67:287\u0026ndash;299\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZou J, Semagn K, Iqbal M, Chen H, Asif M, N\u0026rsquo;Diaye A, Navabi A, Perez-Lara E, Pozniak C, Yang RC, Randhawa H, Spaner D (2017) QTLs associated with agronomic traits in the Attila \u0026times; CDC Go spring wheat population evaluated under conventional management. PLoS One 12:1\u0026ndash;20\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"euphytica","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"euph","sideBox":"Learn more about [Euphytica](https://www.springer.com/journal/10681)","snPcode":"10681","submissionUrl":"https://submission.springernature.com/new-submission/10681/3","title":"Euphytica","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Wheat, Heading, Quantitative trait locus, Bulked segregant analysis","lastPublishedDoi":"10.21203/rs.3.rs-4923172/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4923172/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eOptimizing the timing of heading is crucial for achieving stable production in bread wheat (\u003cem\u003eTriticum aestivum\u003c/em\u003e L.). We analyzed a breeding line, \u0026lsquo;B1-0393\u0026rsquo;, that headed 5 days earlier than the Japanese elite cultivar \u0026lsquo;Kitahonami\u0026rsquo;. To identify quantitative trait loci (QTLs) responsible for this difference, we conducted QTL-seq using F\u003csub\u003e3\u003c/sub\u003e progeny of a cross between \u0026lsquo;Kitahonami\u0026rsquo; and \u0026lsquo;B1-0393\u0026rsquo;. We detected QTLs on chromosomes 2D (\u003cem\u003eQDth.kupg-2D\u003c/em\u003e), 3D (\u003cem\u003eQDth.kupg-3D\u003c/em\u003e), and 4A (\u003cem\u003eQDth.kupg-4A\u003c/em\u003e), explaining 10.25%, 2.00%, and 8.24%, respectively, of phenotypic variance. The \u003cem\u003eQDth.kupg-2D\u003c/em\u003e locus corresponded to the major photoperiod-regulating gene \u003cem\u003ePhotoperiod-D1\u003c/em\u003e (\u003cem\u003ePpd-D1\u003c/em\u003e). \u0026lsquo;Kitahonami\u0026rsquo; had the photoperiod-sensitive \u003cem\u003ePpd-D1b\u003c/em\u003e allele and \u0026lsquo;B1-0393\u0026rsquo; had the insensitive \u003cem\u003ePpd-D1a\u003c/em\u003e allele. \u003cem\u003eQDth.kupg-4A\u003c/em\u003e overlapped with previously reported QTLs, while \u003cem\u003eQDth.kupg-3D\u003c/em\u003e was a novel QTL. Segregation analysis using F\u003csub\u003e2\u003c/sub\u003e and F\u003csub\u003e3\u003c/sub\u003e plants confirmed that alleles from \u0026lsquo;B1-0393\u0026rsquo; at all three loci accelerated heading and that the three QTLs had additive effects on days to heading. Our findings could be valuable for regulating heading time to optimize wheat yield.\u003c/p\u003e","manuscriptTitle":"Identification of three QTLs that additively affect heading time in bread wheat (Triticum aestivum L.) by QTL-seq approach","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-16 20:16:05","doi":"10.21203/rs.3.rs-4923172/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-10-27T00:51:53+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-26T14:21:45+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"162672755849647807285088762577941588539","date":"2024-10-19T11:45:49+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-17T02:31:27+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"27214523424906379631830407036128616540","date":"2024-09-27T13:03:48+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-09-27T12:09:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-08-17T07:58:21+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-08-17T07:57:54+00:00","index":"","fulltext":""},{"type":"submitted","content":"Euphytica","date":"2024-08-16T07:21:30+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"euphytica","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"euph","sideBox":"Learn more about [Euphytica](https://www.springer.com/journal/10681)","snPcode":"10681","submissionUrl":"https://submission.springernature.com/new-submission/10681/3","title":"Euphytica","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"4cff039a-d776-48ac-a5e4-d2c8512a2068","owner":[],"postedDate":"September 16th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-11-18T15:59:16+00:00","versionOfRecord":{"articleIdentity":"rs-4923172","link":"https://doi.org/10.1007/s10681-024-03441-z","journal":{"identity":"euphytica","isVorOnly":false,"title":"Euphytica"},"publishedOn":"2024-11-14 15:56:58","publishedOnDateReadable":"November 14th, 2024"},"versionCreatedAt":"2024-09-16 20:16:05","video":"","vorDoi":"10.1007/s10681-024-03441-z","vorDoiUrl":"https://doi.org/10.1007/s10681-024-03441-z","workflowStages":[]},"version":"v1","identity":"rs-4923172","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4923172","identity":"rs-4923172","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}