Jiajia Cao1†
Yaoyao Shang1†Dongmei Xu1Kangle Xu1
Xinran Cheng1Xu Pan1Xue Liu1Mingli Liu1Chang Gao1Shengnan Yan1Hui Yao1Wei Gao1
Jie Lu1Haiping Zhang1*
Cheng Chang1*
Xianchun Xia2Shihe Xiao2Chuanxi Ma1- 1KeyLaboratory of Wheat Biology and Genetic Improvement on Southern Yellow and Huai River Valley, Ministry of Agriculture and Rural Affairs, College of Agronomy, Anhui Agricultural University, Hefei, China
- 2Institute of Crop Sciences, National Wheat Improvement Center, Chinese Academy of Agricultural Sciences, Beijing, China
Improvement of grain weight and size is an important objective for high-yield wheat breeding. In this study, 174 recombinant inbred lines (RILs) derived from the cross between Jing 411 and Hongmangchun 21 were used to construct a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq). Three mapping methods, including inclusive composite interval mapping (ICIM), genome-wide composite interval mapping (GCIM), and a mixed linear model performed with forward–backward stepwise (NWIM), were used to identify QTLs for thousand grain weight (TGW), grain width (GW), and grain length (GL). In total, we identified 30, 15, and 18 putative QTLs for TGW, GW, and GL that explain 1.1–33.9%, 3.1%–34.2%, and 1.7%–22.8% of the phenotypic variances, respectively. Among these, 19 (63.3%) QTLs for TGW, 10 (66.7%) for GW, and 7 (38.9%) for GL were consistent with those identified by genome-wide association analysis in 192 wheat varieties. Five new stable QTLs, including 3 for TGW (Qtgw.ahau-1B.1, Qtgw.ahau-4B.1, and Qtgw.ahau-4B.2) and 2 for GL (Qgl.ahau-2A.1 and Qgl.ahau-7A.2), were detected by the three aforementioned mapping methods across environments. Subsequently, five cleaved amplified polymorphic sequence (CAPS) markers corresponding to these QTLs were developed and validated in 180 Chinese mini-core wheat accessions. In addition, 19 potential candidate genes for Qtgw.ahau-4B.2 in a 0.31-Mb physical interval were further annotated, of which TraesCS4B02G376400 and TraesCS4B02G376800 encode a plasma membrane H+-ATPase and a serine/threonine-protein kinase, respectively. These new QTLs and CAPS markers will be useful for further marker-assisted selection and map-based cloning of target genes.
Introduction
Wheat (Triticum aestivum L.) is one of the most important cereal crops and provides approximately 20% of the dietary calories for humans worldwide. The demand for wheat is continually increasing due to rapid population growth (Organization of Food and Agriculture, 2009). Grain weight and size, including thousand grain weight (TGW), grain length (GL), and grain width (GW), are positively associated with wheat grain yield (Börner et al., 2002; Dholakia et al., 2003; Breseghello and Sorrells, 2007). Therefore, dissecting the genetic basis of these grain-related traits contributes to the improvement of wheat yield.
Numerous quantitative trait loci (QTLs) underlying grain weight and size were detected on almost all wheat chromosomes by linkage mapping and genome-wide association studies (GWAS) (Cui et al., 2014; Simmonds et al., 2014; Wu et al., 2015; Thorwarth et al., 2019). Several candidate genes for TGW, GW, and GL have been isolated in wheat. TaCKX6-D1 encodes a cytokinin oxidase/dehydrogenase (CKX2) associated with grain number and grain weight (Ashikari et al., 2005; Li et al., 2011, 2013). TaGW2-6A encodes a RING-type E3 ubiquitin ligase and functions as a negative regulator of GW (Su et al., 2012; Yang et al., 2012; Vandana et al., 2015). TaCwi-A1 encodes a cell wall invertase for carbon partitioning during early grain filling (Wang et al., 2008; Ma et al., 2012). TaSAP-A1 belongs to a gene family of stress-associated proteins significantly related to TGW, number of grains per spike, and spike length (Chang et al., 2013). TaGS1a encodes a cytosolic glutamine synthetase underlying grain size (Ying et al., 2013). TaGS-D1 is associated with TGW and GL (Zhang et al., 2014). TaSus1 and TaSus2 encode two isoforms of sucrose synthase (Jiang et al., 2011; Hou et al., 2014). TaGS5-3A regulates grain size (Ma L. et al., 2015). 6-SFT-A2 is involved in fructan biosynthesis and significantly associated with TGW under rainfed conditions (Yue et al., 2015). TaSnRK2.3, TaSnRK2.9, and TaSnRK2.10, three members of the SnRK2 family, are associated with TGW (Miao et al., 2017; Zhang Z. et al., 2017; Rehman et al., 2019). More recently, Cheng et al. (2020) isolated the Tasg-D1 gene, which encodes serine/threonine protein kinase glycogen synthase kinase 3 (STKc_GSK3), leading to formation of round grains in Indian dwarf wheat (Triticum sphaerococcum Perc.). TaSDIR1-4A, which encodes a RING-type E3 ubiquitin ligase is associated with TGW in well-watered and heat-stress environments (Wang J. et al., 2020). The physical interaction of the proteins encoded by TaDA1-A and TaGW2-6B was significantly associated with grain size and weight (Liu H. et al., 2019; Liu J. et al., 2019). TaGS3-7A, a homologous gene of OsGS3, was shown to be associated with grain weight (Yang et al., 2019b). In addition, TaTGW6 (Hu et al., 2016a), Tabas1 (Zhu et al., 2016), TaCKX4 (Chang et al., 2015), TaFlo2 (Sajjad et al., 2017), TaTGW-7A (Hu et al., 2016b), TaCYP78A3 (Ma M. et al., 2015), and TaGL3 (Yang et al., 2019a) are also significantly associated with grain weight or size. All of the genes listed above are useful in genetic improvement of wheat yield. However, wheat grain weight and size are complex traits controlled by multiple genes. Identification and validation of more QTLs for grain-related traits will not only promote our understanding of the genetic basis of the two traits, but also accelerate the process of pyramiding favorable alleles in high-yield wheat breeding.
In the past decades, QTL mapping and GWAS were performed based on randomly amplified polymorphic DNA (RAPDs) (Qu and Hancock, 1997), restriction fragment length polymorphisms (RFLPs) (Keim et al., 1990), amplified fragment length polymorphisms (AFLPs) (Vos et al., 1995; Chagne et al., 2002), simple sequence repeats (SSRs) (Torada et al., 2005), and sequence-tagged site (STS) markers (Talbert et al., 1994). However, the number of these markers is limited, and this affects the accuracy and efficiency of QTL mapping and GWAS.
In recent years, several approaches for high-throughput molecular markers supported by next-generation sequencing technologies and SNP chips have developed rapidly, including DNA sequencing (RAD-seq) (Baird et al., 2008), genotyping-by-sequencing (GBS-seq) (Elshire et al., 2011), restriction site-associated specific locus amplified fragment sequencing (SLAF-seq) (Sun et al., 2013), and wheat 90K (Wang et al., 2014), 55K, and 660K SNP arrays (designed by the Chinese Academy of Agricultural Sciences) (Cui et al., 2017; Huang et al., 2019). These high-throughput genotyping methods obviously accelerate the identification of novel loci and candidate genes underlying wheat yield–related traits. Cui et al. (2017) constructed a genetic map with 4959 bin markers using a wheat 660K SNP array in a recombinant inbred line (RIL) population derived from Kenong 9204 × Jing 411 and detected a major QTL for kernel number per spike on chromosome 4A. Fiedler et al. (2017) identified a QTL for test weight on chromosome 1B by GBS-seq of 1184 lines collected from the North Dakota State University durum wheat breeding program. Wang F. et al. (2020) constructed a high-density genetic map (HDGM), including 3556 SLAF-based SNP and SSR markers, and detected 37 QTLs related to yield traits in an RIL population derived from cotton cultivars LMY 22 (high-yield) and LY 343 (superior fiber quality). Notably, SLAF-seq has been widely applied to construct a HDGM for QTL mapping as well as narrowing target regions in several crops, such as maize (Liu et al., 2016; Li et al., 2018), rice (Yang et al., 2017; Quan et al., 2018), soybeans (Li et al., 2017; Han et al., 2019), and peppers (Guo et al., 2017; Zhang et al., 2019). However, few HDGMs have been constructed by SLAF-seq for common wheat.
The aims of this study were to (1) identify a large number of SNP markers to construct a HDGM by SLAF-seq in an RIL population derived from a cross between high-TGW cultivar Jing 411 and low-TGW landrace Hongmangchun 21, (2) identify QTLs for grain weight and size using SNP markers, and (3) design cleaved amplified polymorphism sequences (CAPS) markers closely linked to new stable QTLs for grain weight and size and validate them in 192 wheat varieties and 180 Chinese mini-core wheat accessions.
Materials and Methods
Plant Materials
A RIL population (174 lines) derived from a cross between Jing 411 (average TGW: 41.94 g) and Hongmangchun 21 (average TGW: 19.98 g) (designated JH-RILs) was used for QTL mapping; it was planted during the 2005–2006, 2006–2007, 2007–2008, 2009–2010, 2010–2011, 2011–2012, 2013–2014, 2014–2015, 2015–2016, 2016–2017, and 2017–2018 cropping seasons. One hundred ninety-two wheat varieties (WVs) grown in the 2015–2016, 2016–2017, and 2017–2018 cropping seasons, containing 159 cultivars, 11 advanced breeding lines, and 22 landraces, were chosen for a GWAS panel to validate QTLs identified by linkage analysis. The origin of these varieties has been described in detail by Zhu et al. (2019). One hundred eighty Chinese mini-core collection accessions (CMCCs) planted in the 2014–2015, 2015–2016, and 2016–2017 cropping seasons comprising 74 cultivars, 11 advanced breeding lines, and 95 landraces were used to validate new stable loci for grain weight and size.
The above three populations were planted at the experimental station of Anhui Agricultural University in Hefei (31°58′N, 117°240′E), Anhui Province, China. Field trials were conducted in 2-m-long plots, each consisting of double rows spaced 0.25 m apart in randomized complete blocks with two replications.
TGW, GW, and GL Tests
Sixty spikes per plot were harvested and threshed at physiological maturity, which is characterized by yellow color on the whole plant including leaves, stems, and spikes (Kulwal et al., 2005). Three hundred threshed seeds were used to measure TGW, GL, and GW in triplicate using the SC-G wheat grain appearance quality image analysis system (Hangzhou WSeen Detection Technology Co., Ltd, Hangzhou, China) (Yin et al., 2015).
For JH-RILs, TGW was measured during the 2005–2006, 2006–2007, 2007–2008, 2009–2010, 2010–2011, 2011–2012, 2013–2014, 2014–2015, 2015–2016, 2016–2017, and 2017–2018 cropping seasons, and these measurements were designated 2006TGW-JH, 2007TGW-JH, 2008TGW-JH, 2010TGW-JH, 2011TGW-JH, 2012TGW-JH, 2014TGW-JH, 2015TGW-JH, 2016TGW-JH, 2017TGW-JH, and 2018TGW-JH, respectively. GW and GL were measured in part seasons, and these measurements were designated 2014GW-JH, 2015GW-JH, 2016GW-JH, 2017GW-JH, 2018GW-JH, 2014GL-JH, 2015GL-JH, 2016GL-JH, 2017GL-JH, and 2018GL-JH, respectively.
For 192 WVs, TGW was evaluated during the 2015–2016, 2016–2017, and 2017–2018 cropping seasons, designated 2016TGW-192, 2017TGW-192, and 2018TGW-192, respectively; GW and GL were also measured, designated 2016GW-192, 2017GW-192, 2018GW-192, 2016GL-192, 2017GL-192, and 2018GL-192, respectively.
For 180 CMCCs, TGW was evaluated during the 2014–2015, 2015–2016, and 2016–2017 cropping seasons, designated 2015TGW-CMCC, 2016TGW-CMCC, and 2017TGW-CMCC, respectively; GW and GL were measured and designated 2015GW-CMCC, 2016GW-CMCC, 2017GW-CMCC, 2015GL-CMCC, 2016GL-CMCC, and 2017GL-CMCC, respectively.
Statistical Analysis
Pearson’s correlation analysis and Mann–Whitney U-tests were performed by the software SPSS 20.0 (IBM Corporation, Armonk, NY, United States).
Genomic DNA Isolation
Dry seeds were collected to extract genomic DNA using the SDS method (Kang et al., 1998). DNA quality was tested by ND5000 spectrophotometer (NanoDrop, Wilmington, DE, United States) and in 1% agarose gels.
Construction of SLAF Library
The JH-RILs were genotyped by a modified SLAF-seq strategy to develop genome-wide SNP markers (Sun et al., 2013). The Chinese Spring reference genome IWGSC RefSeq v1.01 was utilized to simulate in silico the number of markers digested by different restriction enzymes and design a pilot SLAF experiment (Sun et al., 2013). A restriction enzyme (RsaI, New England Biolabs, United States) was used to digest the genomic DNA. The fragments ranging from 464 to 484 bp were selected and purified after agarose gel electrophoresis for 100-bp paired-end sequencing. To assess the accuracy of the experiment for SLAF library construction, Oryza sativa L. japonica was chosen as a control in comparison to mapping populations with the same experimental treatment and sequencing.
SLAF-Seq Data Analysis and Genotyping
In order to ensure the accuracy of data analysis, quality control for the raw sequencing data were performed according to the following criteria: (i) reads containing adaptor sequences were filtered out, (ii) reads with unknown bases exceeding 10% in length were filtered out, (iii) reads unaligned to the wheat reference genome IWGSC RefSeq v1.0 (see footnote 1) were filtered out, and (iv) reads with a single end aligned to the wheat reference genome IWGSC RefSeq v1.0 were also filtered out because of the inaccurate physical position.
SLAF paired-end mapped reads were clustered based on sequence similarity identified by alignments to the wheat reference genome IWGSC RefSeq v1.0 with BWA software (Chong et al., 2003). SNPs in all SLAF loci were identified between parents by the Genome Analysis Toolkit (McKenna et al., 2010)2. SLAFs with 8 or more SNPs, which are considered as high-frequency variable regions of wheat, were filtered out, and this influenced the accuracy of following steps. Minor allele frequency evaluation was used to define alleles in each SLAF. SLAF markers were generated following Sun et al. (2013). Polymorphic markers were selected by genotyping the parents and then classified into eight segregation types (ab × cd, ef × eg, hk × hk, lm × ll, nn × np, aa × bb, ab × cc, and cc × ab). For the RIL population, SLAF markers in the segregation pattern of aa × bb were selected. The quality of SLAF markers used to build genetic maps was controlled by the following criteria: (i) average sequence depths must be >6-fold in parents, (ii) markers with more than 25% missing data were filtered out, (iii) markers with significant segregation distortion (P < 0.001) were rejected, and (iv) markers were required to have a logarithm of odds (MLOD) exceeding 3.
Bulked Segregant Analysis of JH-RILs by 660K SNP Arrays
Two high-TGW bulks (HB1 and HB2) and two low-TGW bulks (LB1 and LB2) were made by the same amount of DNA from extremely high and low TGW lines (5 lines per bulk) of JH-RILs based on phenotypic values. The above four bulks and two parents were genotyped using the wheat 660K SNP array by the CapitalBio Corporation (Beijing, China). SNP genotyping and clustering were performed by Genome Studio Polyploid Clustering v1.03. Consistent SNPs without missing data in HB1, HB2, and Jing 411 were classified as group I, and consistent SNPs without missing data in LB1, LB2, and Hongmangchun 21 were classified as group II. Different SNPs between groups I and II were selected for subsequent analysis and development of PCR-based molecular markers. Physical positions of these SNPs were searched by sequence BLAST with IWGSC RefSeq v1.0 (see footnote 1).
Development and Genotyping of PCR-Based Molecular Markers
Seventy-four SNPs with polymorphisms between groups I and II were converted to CAPS markers by Primer Premier 5.04 and used to genotype JH-RILs to construct genetic maps (Supplementary Table 1). Representative CAPS markers closely linked with new stable QTLs were used to genotype 180 CMCCs. A PCR mixture with a total volume of 10 μL comprised 1.0 μL of 10 × PCR buffer, 200 μM dNTPs, 4 pmol of each primer, 0.5 U Taq DNA polymerase, and 50–100 ng template DNA. The reaction was performed in a C1000 Thermal Cycler (Bio-Rad, United States) with the following program: denaturation at 94°C for 5 min followed by 40 cycles at 94°C for 30 s, touchdown starting at 62°C for 30 s (decreasing 0.3°C per cycle), and 72°C for 30 s that were followed by a final extension at 72°C for 8 min. The PCR products were digested with the corresponding restriction enzymes (Supplementary Table 15) for 5 h and separated on 2.5% agarose gels.
Eleven functional markers specific to TaSdr-2A (Zhang Y. et al., 2017), TaSus-2A (Hou et al., 2014), Tamyb10-3A (Himi et al., 2011), TaVp1B3-3B (Yang Y. et al., 2007), TaVp1-b2 (3B) (Chang et al., 2010), Tamyb10-3D (Himi et al., 2011), TaTGW6-4A (Hu et al., 2016a), TaMKK3-4A (Torada et al., 2016), TaGW2-6A (Su et al., 2012), TaGASR7-A1 (7A) (Dong et al., 2014), and TaPTF1-7B (Zhu, 2016) were used for QTL analysis in JH-RILs in combination with SLAF and CAPS markers. PCR conditions and gel electrophoresis for functional markers were performed following the above reported studies (Supplementary Table 2).
Genetic Linkage Map Construction and QTL Mapping
SLAF markers distributed on 21 wheat chromosomes were selected by mapping using IWGSC RefSeq v1.0 (see footnote 1). To ensure the accuracy of genetic linkage maps, the MLOD scores between markers were estimated, and markers with MLOD scores less than 5 were removed. Then, SLAF, CAPS, and gene-specific markers were collectively used to construct a high-density genetic bin map (designated SLAF-map) without redundant markers by IciMapping v4.1 using the Kosambi mapping function in JH-RILs.
The inclusive composite interval mapping (ICIM) program of QTL IciMapping v4.1 (Li et al., 2007)6 was used to detect QTLs for TGW, GL, and GW in JH-RILs with a walk speed of 1 cM and a window size of 10 cM. The genome-wide composite interval mapping (GCIM) of QTL.gCIM.GUI v1.0 (Wen et al., 2019)7 in R version 3.6.0 software was applied to identify QTLs with a random model and a walk speed of 1 cM. For these two methods, an LOD score of 2.5 was used for claiming the presence of QTLs. The QTLnetwork v2.0 software (Yang J. et al., 2007; Yang et al., 2008), which is based on a mixed linear model (MLM), was used for forward–backward stepwise (NWIM) analysis with a threshold of P = 0.05 to select cofactors, multiple linear regression with a 1 cM walk speed, and a window size set at 10 cM.
Adjacent QTLs with the same sign of additive effects satisfying at least one of the following criteria were defined as the same QTL: (1) positions of QTL peaks within 10 cM (Kumar et al., 2016) and (2) QTLs with overlapped confidence intervals (Cui et al., 2017). QTLs for TGW detected in more than five environments and those for GW or GL in more than three environments were defined as stable loci (Cui et al., 2014).
Genotyping and Association Analysis of 192 WVs Using 90K SNP Array
The DNA samples of 192 WVs were genotyped by the Illumina iSelect 90K Infinium SNP array, including 81,587 SNPs (Wang et al., 2014), by the Beijing Compass Biotechnology Co., Ltd. Genotypic clusters for each SNP were determined by the Genome Studio version 2011.1 software (Illumina, see footnote 3). The physical positions of SNPs were obtained from IWGSC RefSeq v1.0 (see footnote 1). The SNPs with less than 20% missing data and a minor allele frequency exceeding 5% were used for association analysis.
Linkage disequilibrium and population structure were analyzed for 192 WVs following methods in our previous study (Zhu et al., 2019). The same K and Q matrix data were used in the present study to identify significant marker–trait associations (MTAs) for grain weight and size by the MLM. The MLM was performed by TASSEL 5.0 to detect MTAs at a significance level of P < 0.001 (Bradbury et al., 2007; Chen et al., 2016).
Gene Annotation
Databases from the IWGSC gene annotation8 were used to perform gene function annotations by BLAST.
Cloning of TraesCS4B02G376400
Primer pair PMA-1 was designed to amplify the partial sequence of TraesCS4B02G376400 based on the putative sequence from IWGSC RefSeq v1.0 (see footnote 1) (Supplementary Table 3). PCR amplifications were performed in 10-μL volumes that included 0.25 μM of each primer, 0.25 mM dNTPs, 100 ng genomic DNA, 0.5 unit Fastpfu polymerase, and 1 μL of 10 × Fastpfu PCR buffer (Beijing TransGen Biotech, Beijing, China). The amplification program consisted of an initial denaturation at 94°C for 5 min followed by 36 cycles of denaturation at 94°C for 45 s, annealing at 52°C for 50 s, and extension at 72°C for 1 min that were followed by a final extension at 72°C for 12 min. Amplified PCR fragments were separated on 1.5% agarose gels. Target fragments were recovered and cloned into a Blunt-zero vector and transformed into T1 competent cells (Beijing TransGen Biotech, Beijing, China).
Sequencing was carried out on an ABI 3500 genetic analyzer (Applied Biosystems, Shanghai, China). Sequence alignments were performed using DNAMAN 6.09.
Development of CAPS Marker for TraesCS4B02G376400
The SNP distinguished by restriction enzyme StyI found in TraesCS4B02G376400 between Jing 411 and Hongmangchun 21 was converted to a CAPS marker PMA-2 by Primer Premier 5.0 (see footnote 4) and used to genotype 180 CMCCs (Supplementary Table 3). The PCR conditions for the marker PMA-2 were consistent with those for 74 CAPS markers mentioned above.
Results
Statistical Analysis of Phenotypic Data
Phenotypic data for TGW, GW, and GL exhibited extensive variations in JH-RILs and 192 WVs across environments. In JH-RILs, the ranges of TGW, GW, and GL were 13.71–68.28 g (mean 35.97 g), 2.25–4.38 mm (mean 3.01 mm), and 3.99–8.55 mm (mean 6.36 mm), respectively. In the 192 WVs, TGW, GW, and GL ranged from 19.58 to 54.21 g (mean 38.29 g), 1.21 to 4.68 mm (mean 3.01 mm), and 4.37 to 7.50 mm (mean 6.23 mm), respectively (Supplementary Table 4).
Highly significant correlations were observed among TGW, GW, and GL values in different environments. TGW showed positive correlations with GW and GL with correlation coefficients ranging from 0.39–0.97 (average 0.67) and 0.43–0.94 (average 0.76) in JH-RILs and from 0.28–0.84 (average 0.54) and 0.19–0.72 (average 0.47) in the 192 WVs (Supplementary Tables 5, 6).
Construction of Genetic Linkage Maps
After SLAF library construction and high-throughput sequencing, 576.7 M clean reads were obtained with each read being 200 bp. The number of SLAF tags was 816,183 for Jing 411; 798,085 for Hongmangchun 21; and 462,529 for progenies. The average depths of the SLAF markers were 21.52 in Jing 411, 19.64 in Hongmangchun 21, and 4.16 in progenies. Based on the filtering criteria described above, 1529 SLAF bin markers containing 4820 high-quality polymorphic SLAF markers in combination with 74 CAPS and 11 functional markers were used for QTL mapping (Supplementary Figure 1). Finally, a high-density genetic bin map, including 1529 SLAF bin markers, 74 CAPS markers, and 11 gene-specific functional markers spanned 2014.43 cM in length with an average marker interval of 1.22 cM. The genetic distances of the A, B, and D genomes were 864.2 cM (608 markers), 900.0 cM (764 markers), and 250.2 cM (242 markers), respectively (Supplementary Figure 1 and Supplementary Table 7).
QTLs for TGW, GW, and GL Identified by Linkage Analysis
A total of 30 QTLs for TGW were identified using three mapping methods (ICIM, GCIM, and NWIM) in JH-RILs and were distributed on chromosomes 1B (3), 2A (3), 2B (2), 2D, 3A (3), 3B, 4A (2), 4B (2), 4D, 5A (3), 6A, 6B (2), 7A (2), 7B (2), and 7D (2). These QTLs explained 1.1–33.9% of the phenotypic variance, especially five QTLs (Qtgw.ahau-2A.3, Qtgw.ahau-5A.1, Qtgw.ahau-7A.2, Qtgw.ahau-7B.1, and Qtgw.ahau-7B.2) that had relatively high phenotypic variance explained (PVE) ranging from 10.3 to 16.5%. In addition, 13 QTLs for TGW were detected by two methods across environments and were distributed on chromosomes 2A (2), 3A, 3B, 4A, 4D, 5A (2), 7A (2), 7B (2), and 7D. Five QTLs, including Qtgw.ahau-1B.1, Qtgw.ahau-2D, Qtgw.ahau-4B.1, Qtgw.ahau-4B.2, and Qtgw.ahau-6B.1, were collectively detected by all three methods across environments. Particularly, Qtgw.ahau-4B.2 was identified in 10 environments, explaining 11.7% of the phenotypic variances on average and, thus, was considered a stable and major QTL (Figure 1, Table 1, and Supplementary Table 8).
Figure 1. New stable QTLs for TGW and GL in the Jing 411/Hongmangchun 21 population. Uniform million bp (Mb) scale is shown on the left. Green, orange, and red blocks indicate QTLs detected by ICIM performed by QTL IciMapping v4.1, GCIM using QTL.gCIM.GUI v1.0, and composite interval mapping based on mixed linear model using forward–backward stepwise (NWIM) in QTLNetwork v2.0, respectively.
Table 1. QTLs for grain weight and size identified by all three mapping methods in the Jing 411/Hongmangchun 21 RIL population.
Fifteen QTLs associated with GW were mapped on chromosomes 1B, 2D (3), 3B, 4A, 4B (3), 4D, 5A, 5B, 7A (2), and 7B using the aforementioned three mapping methods. Among them, 6 QTLs for GW on chromosomes 1B, 2D, 3B, 4A, 4B, and 4D were detected by two methods across environments. Notably, two QTLs (Qgw.ahau-2D.3 and Qgw.ahau-4B.3), which explained 18.5 and 20.0% of the phenotypic variances on average, respectively, were jointly detected by all three methods (Table 1 and Supplementary Table 8).
Eighteen QTLs for GL were identified on chromosomes 1B (3), 2A (3), 2B, 2D, 3A, 4B (4), 6B, 7A (3), and 7B with three mapping methods. Of these, two QTLs, Qgl.ahau-1B.1 and Qgl.ahau-4B.3, were identified by two methods in more than one environment. Three QTLs, including Qgl.ahau-2A.1, Qgl.ahau-6B, and Qgl.ahau-7A.2, were consistently detected by all three methods in several environments. Of these, Qgl.ahau-2A.1 and Qgl.ahau-7A.2 explained 12.1% and 11.9% of the phenotypic variances, respectively, and were considered as stable and major QTLs (Figure 1, Table 1, and Supplementary Table 8).
Validation of QTLs for TGW, GW, and GL by GWAS
Association analysis was performed to validate QTLs for TGW, GW, and GL detected in JH-RILs using 192 WVs in three environments based on an MLM (Supplementary Figure 2). The results indicated that 19 (63.3%) of 30 QTLs for TGW, 10 (66.7%) of 15 for GW, and 7 (38.9%) of 18 for GL detected in JH-RILs were consistent with those identified by GWAS in 192 varieties. Among them, two reported QTLs for GW (Qgw.ahau-2D.1, 272.5–420.5 Mb; Qgw.ahau-4D, 226.4–409.4 Mb) identified in JH-RILs were closely linked with two stable MTAs (D_GB5Y7FA01EHPZX_186, 353.1 Mb on chromosome 2D; IACX65, 352.4 Mb on chromosome 4D) as detected by GWAS in 192 WVs across all environments. Four reported intervals for TGW harboring Qtgw.ahau-5A.3 (559.3–596.6 Mb), Qtgw.ahau-6B.2 (679.6–698.3 Mb), Qtgw.ahau-7A.1 (99.7–200.5 Mb), Qtgw.ahau-7D.1 with positive addictive effect (256.1–383.1 Mb)/Qtgw.ahau-7D.2 with negative addictive effect (270.6–391.9 Mb) identified in JH-RILs were near to four MTAs (Excalibur_c56952_411, 569.8 Mb on chromosome 5A; RAC875_c17347_312, 694.2 Mb on chromosome 6B; RAC875_c5986_3670, 200.1 Mb on chromosome 7A; RAC875_c45846_454, 324.5 Mb on chromosome 7D) detected in 192 WVs, explained 12.6, 11.0, 10.7, and 12.5% of the phenotypic variances on average, respectively. Notably, three novel and major QTLs, including Qtgw.ahau-4B.2 (660.1–660.4 Mb), Qgl.ahau-2A.1 (44.2–54.4 Mb), and Qgl.ahau-7A.2 (446.4–496.9 Mb), were adjacent to three MTAs (RAC875_c48025_483, 666.1 Mb on chromosome 4B; Tdurum_contig55610_784, 42.5 Mb on chromosome 2A; Kukri_c5101_2636, 498.9 Mb on chromosome 7A) based on their physical positions in 192 wheat varieties, respectively (Table 1 and Supplementary Table 9).
Validation of New Stable QTLs in 180 CMCCs
Four CAPS markers (1B-JHMfeI for Qtgw.ahau-1B.1, 2A-CAPSmin for Qgl.ahau-2A.1, 4B-8621 for Qtgw.ahau-4B.1, and 7A-3738 for Qgl.ahau-7A.2) were further developed to validate the associations of four new stable QTLs with grain weight and size in 180 CMCCs (Table 1 and Supplementary Tables 1, 8). Allelic variations of the CAPS markers were detected, designated as 1B-JHMfeI-C/1B-JHMfeI-T, 2A-CAPSmin-T/2A-CAPSmin-C, 4B-8621-C/4B-8621-G, and 7A-3738-G/7A-3738-C, respectively (Supplementary Figure 3). Compared with the Jing 411 genotypes (1B-JHMfeI-C, 2A-CAPSmin-T, 4B-8621-C, and 7A-3738-G), the Hongmangchun 21 genotypes (1B-JHMfeI-T, 2A-CAPSmin-C, 4B-8621-G, and 7A-3738-C) were significantly associated with lower TGW, GW, and GL values across environments (P < 0.01 or P < 0.05) (Table 2).
Table 2. Association of allelic variations of four CAPS markers for four novel QTLs with TGW, GW, and GL in 180 CMCC.
Candidate Gene Prediction for Qtgw.ahau-4B.2
The new stable and major QTL Qtgw.ahau-4B.2 was flanked by two CAPS markers, 4B-2133 and 4B-8416, in a physical interval of 0.31 Mb (Table 1). According to IWGSC annotation databases, the target region of Qtgw.ahau-4B.2 contains 19 genes. Among them, TraesCS4B02G376400 and TraesCS4B02G376800, which encode a plasma membrane H+-ATPase and a serine/threonine-protein kinase, respectively, were considered as potential candidate genes (Supplementary Table 10).
Cloning of TraesCS4B02G376400 and Validation in 180 CMCCs
Sequence analysis indicated that TraesCS4B02G376400 comprised 20 exons and 19 introns, with a complete sequence of 6024 bp (Figure 2A). After partial sequencing of TraesCS4B02G376400 (1919–5084 bp), a SNP (A-G) was identified in the 12th intron between Jing 411 and Hongmangchun 21 (Figure 2B). Then, a gene-specific marker, PMA-2, was developed based on this SNP and used to detect allelic variations of TraesCS4B02G376400 among 180 CMCCs (Figure 2C). Compared with the Jing 411-type (PMA-2-A), the Hongmangchun 21-type (PMA-2-G) was significantly associated with lower TGW, GW, and GL values across environments (P < 0.01) (Table 3).
Figure 2. CAPS marker development based on a SNP in TraesCS4B01G376400. (A) Diagram of the gene structure of TraesCS4B01G376400. (B) SNP (A–G) found in the partial sequence (1919–5084 bp) of TraesCS4B01G376400. (C) Electrophoresis patterns of the CAPS marker PMA-2 digested by StyI between two parents Jing 411 and Hongmangchun 21. M indicates marker. “J” and “H” indicate electrophoresis patterns of Jing 411 genotype (undigested) and Hongmangchun 21 genotype (digested).
Table 3. Association of allelic variations of the CAPS marker PMA-2 for Qtgw.ahau-4B.2 with TGW, GW, and GL in 180 CMCC.
Discussion
Reliability of QTL Mapping
The reliability of QTL mapping is crucial for subsequent fine mapping and map-based cloning of target genes and is affected by the density and genotyping of markers used to construct genetic linkage maps (Cui et al., 2017; Han et al., 2019). Recently, RAD-seq (Baird et al., 2008), GBS-seq (Elshire et al., 2011), SLAF-seq (Sun et al., 2013), and SNP arrays (Wang et al., 2014)10 were performed as low-cost and efficient strategies to develop abundant SNPs to increase the density of markers. In the present study, SLAF tags (4820), CAPS (74), and gene-specific markers (11) were collectively used to build a high-density genetic linkage map harboring 1614 bin markers spanning 2014.43 cM in length with an average marker interval of 1.22 cM (Supplementary Figure 1).
It is well known that different mapping methods are based on specific genetic models as well as their corresponding statistical hypothesis, from which the result is in fact a probability statement. Therefore, it is better to use multiple mapping methods to reduce the risks of identifying ghost QTLs or missing real QTLs (Su et al., 2010). Zhang et al. (2016) identified 63 and 16 additive QTLs for fiber strength of cotton using a composite interval mapping (CIM) method with WinQTLCartographer2.5 and an ICIM model by QTL IciMapping4.1, respectively, and four QTLs detected by both CIM and ICIM were considered as stable loci for the fiber strength of cotton. In our previous study, based on both single- and multi-locus MLMs, 23 and 39 MTAs for preharvest sprouting resistance were detected by GWAS, respectively, and 6 loci were jointly identified by the two models and were, thus, considered stable loci (Zhu et al., 2019).
In addition, a combination of QTL analysis and GWAS is also an important strategy to detect reliable loci. Wang et al. (2018) identified 41 QTLs for maize tocopherol content by linkage mapping in 6 RIL populations and 32 significant loci by GWAS in a diverse panel of 508 inbred lines, and a major QTL co-localized in both linkage analysis and GWAS was finely mapped and characterized as a non-tocopherol pathway gene involved in the modulation of natural tocopherol variation.
In this study, QTL analysis with three mapping methods (ICIM, GCIM, and NWIM) and GWAS with MLM were collectively performed to identify loci for grain weight and size. Notably, five QTLs for TGW, two for GW, and four for GL were identified by all three mapping methods in JH-RILs. Six of them, including Qtgw.ahau-4B.1, Qtgw.ahau-4B.2, Qgw.ahau-2D, Qgw.ahau-4B.3, Qgl.ahau-2A.1, and Qgl.ahau-7A.2, were validated by GWAS in 192 WVs. Thereby, these loci were considered to be reliable (Table 1, Supplementary Table 9, and Supplementary Figure 2).
Comparisons With Previous Studies
All QTLs for grain weight and size detected in the present study were compared with previously reported loci based on the wheat reference genome sequence (IWGSC RefSeq 1.0; IWGSC, 2018). In total, 23 (76.7%) of 30 QTLs for TGW, 11 (73.3%) of 15 for GW, and 13 (72.2%) of 18 for GL were completely or partially coincident with previous results (Supplementary Table 8). Qtgw.ahau-2D and Qgw.ahau-2D.3 (543–561 Mb) were detected close to a major QTL for TGW (Q.TKW.ui-2D-1, which is flanked by cfd73 and IWB1093, 553 Mb) that was identified by Liu et al. (2018). Qgl.ahau-4B.1 (31–38 Mb), identified in JH-RILs, was close to a major and stable locus for TGW (TKW-AX_110713957, 42 Mb) (Li et al., 2019). Moreover, Qtgw.ahau-6B.1 and Qgl.ahau-6B (676–679 Mb) were detected near a stable QTL for TGW (TKW-AX_109917592, 675–677 Mb) that was identified by Li et al. (2019) (Table 1).
It is noteworthy that six QTL regions detected in this study contained eight reported genes for grain weight and size, including TaFlo2-A1 (2A) (Sajjad et al., 2017), TaSnRK2.9-5A (Rehman et al., 2019), TaGL3-5A (Yang et al., 2019a), TaGW2-6A (Su et al., 2012; Yang et al., 2012; Hong et al., 2014; Qin et al., 2014; Vandana et al., 2015; Wang et al., 2017), TaSus1-7A (Hou et al., 2014), TaGASR7-7A (Dong et al., 2014), and TaTGW-7A (Hu et al., 2016b). This supports the reliability of QTL mapping performed in the current study. However, most of these reported genes had minor effects on grain weight and size traits in JH-RILs except for Qtgw.ahau-7A.1 (co-located with Qgw.ahau-7A.1 and Qgl.ahau-7A.1), indicating that control of grain weight and size is multigenic (Supplementary Table 8).
In particular, several putative, novel QTLs were identified in the present study compared to previous studies. For instance, Qtgw.ahau-4B.1, Qtgw.ahau-4B.2, and Qgl.ahau-2A.2, which explained 12.4, 11.7, and 12.1% of the phenotypic variances, on average, respectively, have not been reported by previous studies and are, thus, regarded as novel loci (Table 1 and Supplementary Table 8). Interestingly, two loci for TGW and GL detected by all three mapping methods (Qtgw.ahau-1B.1 and Qgl.ahau-7A.2) in the present study overlapped with the reported QTLs associated with kernel number per spike (KNS) (Cui et al., 2014; Wang et al., 2017), indicating a pleiotropy or close linkage between them (Table 1 and Supplementary Table 8).
Pleiotropic or Co-localized QTLs
Thousand grain weight shows significant, positive correlations with GW and GL as shown previously and in the current study (Sun et al., 2009; Drikvand et al., 2013; Wang J. et al., 2020) (Supplementary Tables 5, 6), suggesting that grain weight and size were simultaneously selected in high-yield wheat breeding. In this study, several QTL intervals were identified for multiple yield-related traits. For example, two pairs of QTLs related to TGW and GW were detected in same intervals on chromosomes 2D and 4A; four QTL intervals related to both TGW and GL were identified on chromosomes 1B, 2A, 3A, and 6B; two QTL intervals for both GW and GL were detected on chromosomes 2D and 4B. Moreover, six QTL intervals associated with TGW, GW, and GL were identified on chromosomes 1B, 4B (2), 7A (2), and 7B, respectively. These include four novel QTL intervals, 648–652 Mb of chromosome 1B (Qtgw.ahau-1B.1, Qgw.ahau-1B.1, and Qgl.ahau-1B.1), 540–542 Mb of chromosome 4B (Qtgw.ahau-4B.1, Qgw.ahau-4B.2, and Qgl.ahau-4B.2), 660.11–660.42 Mb of chromosome 4B (Qtgw.ahau-4B.2, Qgw.ahau-4B.3, and Qgl.ahau-4B.4), and 446–496 Mb of chromosome 7A (Qtgw.ahau-7A.2, Qgw.ahau-7A.2, and Qgl.ahau-7A.2). Interestingly, four novel QTL intervals for grain weight and size identified in this study were co-located with QTLs for KNS detected in previous research, including Qtgw.ahau-1B.1 (co-located with Qgl.ahau-1B.1) (Cui et al., 2014), Qtgw.ahau-3A.2 (Li et al., 2019), Qtgw.ahau-4D (Cui et al., 2014), and Qtgw.ahau-7A.2 (co-located with Qgw.ahau-7A.2 and Qgl.ahau-7A.2) (Wang et al., 2017). In brief, these findings indicate that the above QTL intervals are either pleiotropic or tightly linked regions controlling multiple yield-related traits.
Candidate Gene Prediction of Qtgw.ahau-4B.2
The new stable locus Qtgw.ahau-4B.2 spanning a physical interval of 0.31 Mb was detected in 10 of 11 environments in which 19 genes were annotated (Table S10). Among them, TraesCS4B02G376400 and TraesCS4B02G376800 encode a plasma membrane H+-ATPase and a serine/threonine-protein kinase, respectively. The plasma membrane (PM) H+-ATPase is an important ion pump in the plant cell membrane. By extruding protons from the cell and generating a membrane potential, this pump energizes the PM, which is a prerequisite for plant growth. Thus, the PM H+-ATPase is regarded as a driver of growth (Falhof et al., 2016). Rober-Kleber et al. (2003) indicated that the PM H+-ATPase is involved in auxin-mediated cell elongation during wheat embryo development. Auxin activates the proton pump, resulting in apoplastic acidification that contributes to cell wall loosening and elongation of the scutellum. Therefore, the PM H+-ATPase is a component of the auxin-signaling cascade that may direct pattern formation in embryos. Moreover, several reported genes related to TGW belong to the serine/threonine-protein kinase family, such as TaSnRK2.3 (Miao et al., 2017), TaSnRK2.9 (Rehman et al., 2019), TaSnRK2.10 (Zhang Z. et al., 2017), and Tasg-D1 (Cheng et al., 2020), implying that the serine/threonine-protein kinase proteins play important roles in regulation of wheat grain development. Taken together, TraesCS4B02G376400 and TraesCS4B02G376800 are likely to be candidate genes of Qtgw.ahau-4B.2.
Conclusion
A total of 30, 15, and 18 putative additive QTLs for TGW, GW, and GL, respectively, were identified by SLAF-map in JH-RILs using three mapping methods. Particularly, five novel QTLs with stable and significant effects, including Qtgw.ahau-1B.1, Qgl.ahau-2A.1, Qtgw.ahau-4B.1, Qtgw.ahau-4B.2, and Qgl.ahau-7A.2 were identified by all three mapping methods and further validated in a natural population. In addition, TraesCS4B02G376400 and TraesCS4B02G376800 were considered as potential candidate genes underlying Qtgw.ahau-4B.2. The novel QTLs and CAPS markers developed will be helpful for map-based cloning of the target regions and gene pyramiding in breeding for wheat PHS resistance.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
JC and YS conceived the study, put into effect the main linkage analyses, and drafted the manuscript. DX, KX, and XC took part in the experiments and drafting of the manuscript. XP, XL, ML, CG, SY, HY, and WG processed the experimental data and helped to draft the manuscript. JL, HZ, CC, XX, SX, and CM conceived and guided the experiments, and helped in coordinating the project and drafting the manuscript. All the authors read and accepted the final manuscript.
Funding
Our research was supported by grants from the National Key Research and Development plan “Breeding new wheat varieties with high-yielding, high-quality and water-saving in the south of Huang-Huai River winter wheat area”—the breeding of new wheat germplasm and varieties with resistance to adversity (2017YFD0100703); The breeding and industrialization of new wheat varieties with high yield and multi-resistance (18030701182), The China Agriculture Research System (CARS-03); The National Key Research and Development Program of China (2016YFD0101802, 2017YFD0100804); Wheat Genetics and Breeding Research Platform Innovation Team of Anhui’s University (2015-), Jiangsu Collaborative Innovation Center for Modern Crop Production (JCIC-MCP); The Agriculture Research System of Anhui Province (AHCYTX-02), and The University Synergy Innovation Program of Anhui Province (GXXT-2019-033).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling Editor declared a past collaboration with one of the authors SX.
Acknowledgments
We thank Prof. Jizeng Jia for kindly providing Chinese micro-core wheat collections. We would like to thank TopEdit (www.topeditsci.com) for its linguistic assistance during the preparation of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.584859/full#supplementary-material
Supplementary Figure 1 | Genetic map containing 1614 bin markers for linkage mapping in the Jing 411/Hongmangchun 21 RILs population.
Supplementary Figure 2 | Manhattan plots and Quantile–quantile (Q-Q) plots for TGW, GW, and GL of 192 WVs (lines) by the MLM in Tassel v5.0.
Supplementary Figure 3 | Electrophoresis patterns of four CAPS markers 1B-JHMfeI, 2A-CAPSmin, 4B-8621 and 7A-3738 between two parents Jing 411 and Hongmangchun 21.
Supplementary Table 1 | Information of 74 CAPS markers developed in this study.
Supplementary Table 2 | Information on 11 gene-specific markers.
Supplementary Table 3 | Primers used for cloning and marker development of TraesCS4B02G376400.
Supplementary Table 4 | Phenotypic data of grain weight and size in the Jing 411/Hongmangchun 21 RIL population and 192 WVs.
Supplementary Table 5 | Pearson’s correlations between TGW, GW, and GL in the Jing 411/Hongmangchun 21 RIL population.
Supplementary Table 6 | Pearson’s correlations between TGW, GW, and GL in 192 WVs.
Supplementary Table 7 | Summary of the high-density genetic linkage bin maps constructed by SLAF-seq, CAPS, and gene-specific markers in the Jing 411/Hongmangchun 21 population.
Supplementary Table 8 | QTLs for grain weight and size identified in the Jing 411/Hongmangchun 21 RIL population.
Supplementary Table 9 | Thirty-one markers for grain weight and size in 192 WVs identified by GWAS.
Supplementary Table 10 | Gene annotation of the Qtgw.ahau-4B.2 interval.
Footnotes
- ^ https://urgi.versailles.inra.fr/blast_iwgsc/blast.php
- ^ https://gatk.broadinstitute.org
- ^ http://www.illumina.com
- ^ https://www.PremierBiosoft.com
- ^ https://www.neb-china.com
- ^ http://www.isbreeding.net
- ^ https://cran.r-project.org/web/packages/QTL.gCIMapping.GUI/index.html
- ^ https://wheat-urgi.versailles.inra.fr/Tools/JBrowse
- ^ http://www.lynnon.com/
- ^ https://www.cerealsdb.uk.net/cerealgenomics/CerealsDB/indexNEW.php
References
Ashikari, M., Sakakibara, H., Lin, S. Y., Yamamoto, T., Takashi, T., Nishimura, A., et al. (2005). Cytokinin oxidase regulates rice grain production. Science 309, 741–745. doi: 10.1126/science.1113373
Baird, N. A., Etter, P. D., Atwood, T. S., Currey, M., Shiver, A. L., Lewis, Z. A., et al. (2008). Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One 3:e03376. doi: 10.1371/journal.pone.0003376
Börner, A., Schumann, E., Fürste, A., Cöster, H., Leithold, B., Röder, M., et al. (2002). Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet. 105, 921–936. doi: 10.1007/s00122-002-0994-1
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Breseghello, F., and Sorrells, M. E. (2007). QTL analysis of kernel size and shape in two hexaploid wheat mapping populations. Field Crops Res. 101, 172–179. doi: 10.1016/j.fcr.2006.11.008
Cao, S., Xu, D., Hanif, M., Xia, X., and He, Z. (2020). Genetic architecture underpinning yield component traits in wheat. Theor. Appl. Genet. 133, 1811–1823. doi: 10.1007/s00122-020-03562-8
Chagne, D., Lalanne, C., Madur, D., Kumar, S., Frigerio, J., Krier, C., et al. (2002). A high density genetic map of maritime pine based on AFLPs. Ann. For. Sci. 59, 627–636. doi: 10.1051/forest:2002048
Chang, C., Feng, J. M., Si, H. Q., Yin, B., Zhang, H. P., and Ma, C. X. (2010). Validating a novel allele of viviparous-1 (Vp-1Bf) associated with high seed dormancy of Chinese wheat landrace, Wanxianbaimaizi. Mol. Breeding 25, 517–525. doi: 10.1007/s11032-009-9350-3
Chang, C., Lu, J., Zhang, H., Ma, C., and Sun, G. (2015). Copy number variation of cytokinin oxidase gene Tackx4 associated with grain weight and chlorophyll content of flag leaf in common wheat. PLoS One 10:e0145970. doi: 10.1371/journal.pone.0145970
Chang, J., Zhang, J., Mao, X., Li, A., Jia, J., and Jing, R. (2013). Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat. Planta 237, 1495–1508. doi: 10.1007/s00425-013-1860-x
Chen, G., Zhang, H., Deng, Z., Wu, R., and Tian, J. (2016). Genome-wide association study for kernel weight-related traits using SNPs in a Chinese winter wheat population. Euphytica 212, 1–13. doi: 10.1007/s10681-016-1750-y
Cheng, X., Xin, M., Xu, R., Chen, Z., Cai, W., Chai, L., et al. (2020). A single amino acid substitution in STKc_GSK3 kinase conferring semispherical grains and its implications for the origin of Triticum sphaerococcum Perc. Plant Cell 32, 923–934. doi: 10.1105/tpc.19.00580
Chong, C. F., Li, Y. C., Wang, T. L., and Hang, C. (2003). Stratification of adverse outcomes by preoperative risk factors in coronary artery bypass graft patients: an artificial neural network prediction model. Amia Annu. Symp. Proc. 2003, 160–164.
Cui, F., Zhang, N., Fan, X. L., Zhang, W., Zhao, C. H., Yang, L. J., et al. (2017). Utilization of a wheat 660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel numbe. Sci. Rep. 7:3788. doi: 10.1038/s41598-017-04028-6
Cui, F., Zhao, C., Ding, A., Li, J., Wang, L., Li, X., et al. (2014). Construction of an integrative linkage map and QTL mapping of grain yield-related traits using three related wheat RIL populations. Theor. Appl. Genet. 127, 659–675. doi: 10.1007/s00122-013-2249-8
Cuthbert, J. L., Somers, D. J., Brulebabel, A., Brown, P. D., and Crow, G. H. (2008). Molecular mapping of quantitative trait loci for yield and yield components in spring wheat (Triticum aestivum L.). Theor. Appl. Genet. 117, 595–608. doi: 10.1007/s00122-008-0804-5
Dholakia, B. B., Ammiraju, J. S. S., Singh, H., Lagu, M. D., Röder, M. S., Rao, V. S., et al. (2003). Molecular marker analysis of kernel size and shape in bread wheat. Plant Breeding 122, 392–395. doi: 10.1046/j.1439-0523.2003.00896.x
Dong, L., Wang, F., Liu, T., Dong, Z., Li, A., Jing, R., et al. (2014). Natural variation of TaGASR7-A1 affects grain length in common wheat under multiple cultivation conditions. Mol. Breeding 34, 937–947. doi: 10.1007/s11032-014-0087-2
Drikvand, R., Bihamta, M. R., Najafian, G., and Ebrahimi, A. (2013). Kernel quality association and path analysis in bread wheat. Int. J. Biol. 5, 73–79. doi: 10.5539/ijb.v5n3p73
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6:e19379. doi: 10.1371/journal.pone.0019379
Falhof, J., Pedersen, J. T., Fuglsang, A. T., and Palmgren, M. (2016). Plasma membrane H+-ATPase regulation in the center of plant physiology. Mol. Plant 9, 323–337. doi: 10.1016/j.molp.2015.11.002
Fiedler, J. D., Evan, S., Yuan, L., Monika, M. D. J., Hegstad, J. B., Chen, B., et al. (2017). Genome-wide association and prediction of grain and semolina quality traits in durum wheat breeding populations. Plant Genome 10, 1–12. doi: 10.3835/plantgenome2017.05.0038
Gao, F., Wen, W., Liu, J., Rasheed, A., Yin, G., Xia, X., et al. (2015). Genome-wide linkage mapping of QTL for yield components, plant height and yield-related physiological traits in the Chinese wheat cross Zhou 8425B/Chinese spring. Front. Plant Sci. 6:1099. doi: 10.1007/s10681-016-1682-6
Garcia, M., Eckermann, P., Haefele, S., Satija, S., Sznajder, B., Timmins, A., et al. (2019). Genome-wide association mapping of grain yield in a diverse collection of spring wheat (Triticum aestivum L.) evaluated in southern Australia. PLoS One 14:e0211730. doi: 10.1371/journal.pone.0211730
Guo, G., Wang, S., Liu, J., Pan, B., Diao, W., Ge, W., et al. (2017). Rapid identification of QTLs underlying resistance to cucumber mosaic virusin pepper (Capsicum frutescens). Theor. Appl. Genet. 130, 41–52. doi: 10.1007/s00122-016-2790-3
Han, J., Han, D., Guo, Y., Yan, H., and Qiu, L. (2019). QTL mapping pod dehiscence resistance in soybean (Glycine max L. Merr.) using specific-locus amplified fragment sequencing. Theor. Appl. Genet. 132, 2253–2272. doi: 10.1007/s00122-019-03352-x
Himi, E., Maekawa, M., Miura, H., and Noda, K. (2011). Development of PCR markers for Tamyb10 related to R-1, red grain color gene in wheat. Theor. Appl. Genet. 122, 1561–1576. doi: 10.1007/s00122-011-1555-2
Hong, Y., Chen, L., Du, L. P., Su, Z., Wang, J., Ye, X., et al. (2014). Transcript suppression of TaGW2 increased grain width and weight in bread wheat. Funct. Integr. Genomics 14, 341–349. doi: 10.1007/s10142-014-0380-5
Hou, J., Jiang, Q., Hao, C., Wang, Y., Zhang, H., and Zhang, X. (2014). Global selection on sucrose synthase haplotypes during a century of wheat breeding. Plant Physiol. 164, 1918–1929. doi: 10.1104/pp.113.232454
Hu, M. J., Zhang, H. P., Cao, J. J., Zhu, X. F., Wang, S. X., Jiang, H., et al. (2016a). Characterization of an IAA-glucose hydrolase gene TaTGW6 associated with grain weight in common wheat (Triticum aestivum L.). Mol. Breeding 36:25. doi: 10.1007/s11032-016-0449-z
Hu, M. J., Zhang, H. P., Liu, K., Cao, J. J., Wang, S. X., Jiang, H., et al. (2016b). Cloning and characterization of TaTGW-7A gene associated with grain weight in wheat via SLAF-seq-BSA. Front. Plant Sci. 7:1902. doi: 10.3389/fpls.2016.01902
Huang, S., Wu, J., Wang, X., Mu, J., Xu, Z., Zeng, Q., et al. (2019). Utilization of the genomewide wheat 55K SNP array for genetic analysis of stripe rust resistance in common wheat line P9936. Phytopathology 109, 819–827. doi: 10.1094/phyto-10-18-0388-r
Huang, X. Q., Cloutier, S., Lycar, L., Radovanovic, N., Humphreys, D. G., Noll, J. S., et al. (2006). Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (Triticum aestivum L.). Theor. Appl. Genet. 113, 753–766. doi: 10.1007/s00122-006-0346-7
Jiang, Q., Hou, J., Hao, C., Wang, L., Ge, H., Dong, Y., et al. (2011). The wheat (T. aestivum) sucrose synthase 2 gene (TaSus2) active in endosperm development is associated with yield traits. Funct. Integr. Genomics 11, 49–61. doi: 10.1007/s10142-010-0188-x
Juliana, P., Poland, J., Huerta-Espino, J., Shrestha, S., and Singh, R. P. (2019). Improving grain yield, stress resilience and quality of bread wheat using large-scale genomics. Nat. Genet. 51, 1–10. doi: 10.1038/s41588-019-0496-6
Kang, H. W., Cho, Y. G., Yoon, U. H., and Eun, M. Y. (1998). A rapid DNA extraction method for RFLP and PCR analysis from a single dry seed. Plant Mol. Biol. Report 16:90. doi: 10.1023/A:1007418606098
Keim, P., Diers, B. W., Olson, T. C., and Shoemaker, R. C. (1990). RFLP mapping in soybean: association between marker loci and variation in quantitative traits. Genetics 126, 735–742. doi: 10.1007/BF00056365
Kulwal, P. L., Kumar, N., Gaur, A., Khurana, P., Khurana, J. P., Tyagi, A. K., et al. (2005). Mapping of a major QTL for pre-harvest sprouting tolerance on chromosome 3A in bread wheat. Theor. Appl. Genet. 111, 1052–1059. doi: 10.1007/s00122-005-0021-4
Kumar, A., Mantovani, E. E., Seetan, R., Soltani, A., Echeverry-Solarte, M., Jain, S., et al. (2016). Dissection of genetic factors underlying wheat kernel shape and size in an Elite × Nonadapted cross using a high density SNP linkage map. Plant Genome 9, 1–22. doi: 10.3835/plantgenome2015.09.0081
Li, B., Fan, S., Yu, F., Chen, Y., Zhang, S., Han, F., et al. (2017). High-resolution mapping of QTL for fatty acid composition in soybean using specific-locus amplified fragment sequencing. Theor. Appl. Genet. 130, 1467–1479. doi: 10.1007/s00122-017-2902-8
Li, F., Wen, W., Liu, J., Zhang, Y., and Cao, S. (2019). Genetic architecture of grain yield in bread wheat based on genome-wide association studies. BMC Plant Biol. 19:168. doi: 10.1186/s12870-019-1781-3
Li, H., Ye, G., and Wang, J. (2007). A modified algorithm for the improvement of composite interval mapping. Genetics 175, 361–374. doi: 10.1534/genetics.106.066811
Li, M., Tang, D., Wang, K., Wu, X., Lu, L., Yu, H., et al. (2011). Mutations in the F-box gene LARGER PANICLE improve the panicle architecture and enhance the grain yield in rice. Plant Biotech. J. 9, 1002–1013. doi: 10.1111/j.1467-7652.2011.00610.x
Li, R., Song, W., Wang, B., Wang, J., Zhang, D., Zhang, Q., et al. (2018). Identification of a locus conferring dominant resistance to maize rough dwarf disease in maize. Sci. Rep. 8:3248. doi: 10.1038/s41598-018-21677-3
Li, S., Zhao, B., Yuan, D., Duan, M., Qian, Q., Tang, L., et al. (2013). Rice zinc finger protein DST enhances grain production through controlling Gn1a/OsCKX2 expression. Proc. Natl. Acad. Sci. U.S.A. 110, 3167–3172. doi: 10.1073/pnas.1300359110
Liu, C., Zhou, Q., Dong, L., Wang, H., Liu, F., Weng, J., et al. (2016). Genetic architecture of the maize kernel row number revealed by combining QTL mapping using a high-density genetic map and bulked segregant RNA sequencing. BMC Genomics 17:915. doi: 10.1186/s12864-016-3240-y
Liu, H., Li, H., Hao, C., Wang, K., Wang, Y., Qin, L., et al. (2019). TaDA1, a conserved negative regulator of kernel size, has an additive effect with TaGW2 in common wheat (Triticum aestivum L.). Plant Biotech. J. 18, 1330–1342. doi: 10.1111/pbi.13298
Liu, J., Wu, B., Singh, R. P., and Velu, G. (2019). QTL mapping for mineral nutrient concentrations and agronomic-related traits in a hexaploid wheat mapping population. J. Cereal Sci. 88, 57–64. doi: 10.1016/j.jcs.2019.05.008
Liu, Y., Wang, R., Hu, Y. G., and Chen, J. (2018). Genome-wide linkage mapping of quantitative trait loci for late-season physiological and agronomic traits in spring wheat under irrigated conditions. Agronomy 8:60. doi: 10.3390/agronomy8050060
Ma, D., Yan, J., He, Z., Ling, W., and Xia, X. (2012). Characterization of a cell wall invertase gene TaCwi-A1 on common wheat chromosome 2A and development of functional markers. Mol. Breeding 29, 43–52. doi: 10.1007/s11032-010-9524-z
Ma, L., Li, T., Hao, C., Wang, Y., and Zhang, X. (2015). TaGS5-3A, a grain size gene selected during wheat improvement for larger kernel and yield. Plant Biotechnol. J. 14, 1269–1280. doi: 10.1111/pbi.12492
Ma, M., Wang, Q., Li, Z., Cheng, H., and Zhao, H. (2015). Expression of TaCYP78A3, a gene encoding cytochrome P450 CYP78A3 protein in wheat (Triticum aestivum L.), affects seed size. Plant J. 83, 312–325. doi: 10.1111/tpj.12896
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: a mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1016/S0002-9343(99)00454-4
Miao, L., Mao, X., Wang, J., Liu, Z., Zhang, B., Li, W., et al. (2017). Elite haplotypes of a protein kinase gene TaSnRK2.3 associated with important agronomic traits in common wheat. Front. Plant Sci. 8:368. doi: 10.3389/fpls.2017.00368
Mir, R. R., Kumar, N., Jaiswal, V., Girdharwal, N., Prasad, M., Balyan, H. S., et al. (2012). Genetic dissection of grain weight in bread wheat through quantitative trait locus interval and association mapping. Mol. Breeding 29, 963–972. doi: 10.1007/s11032-011-9693-4
Mohler, V., Albrecht, T., Castell, A., Diethelm, M., Schweizer, G., and Hartl, L. (2016). Considering causal genes in the genetic dissection of kernel traits in common wheat. J. Appl. Genet. 57, 467–476. doi: 10.1007/s13353-016-0349-2
Organization of Food and Agriculture (2009). The state of food and agriculture. Livestock in the balance. State Food Agric. 79:410. doi: 10.1080/00369225508735601
Qin, L., Hao, C., Hou, J., Wang, Y., Li, T., Wang, L., et al. (2014). Homologous haplotypes, expression, genetic effects and geographic distribution of the wheat yield gene TaGW2. BMC Plant Biol. 14:107. doi: 10.1186/1471-2229-14-107
Qu, L., and Hancock, J. F. (1997). Randomly amplified polymorphic DNA- (RAPD-) based genetic linkage map of blueberry derived from an interspecific cross between diploid Vaccinium darrowi and tetraploid V. corymbosum. J. Am. Soc. Hortic. Sci. 122, 69–73. doi: 10.1023/A:1008607226805
Quan, R., Wang, J., Hui, J., Bai, H., Lyu, X., Zhu, Y., et al. (2018). Improvement of salt tolerance using wild rice genes. Front. Plant Sci. 8:2269. doi: 10.3389/fpls.2017.02269
Quarrie, S. A., Steed, A., Calestani, C., Semikhodskii, A., Lebreton, C., Chinoy, C., et al. (2005). A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross Chinese Spring×SQ1 and its use to compare QTLs for grain yield across a range of environments. Theor. Appl. Genet. 110, 865–880. doi: 10.1007/s00122-004-1902-7
Rehman, S. U., Wang, J., Chang, X., Zhang, X., Mao, X., and Jing, R. (2019). A wheat protein kinase gene TaSnRK2.9-5A associated with yield contributing traits. Theor. Appl. Genet. 132, 907–919. doi: 10.1007/s00122-018-3247-7
Rober-Kleber, N., Albrechtova, J. T. P., Fleig, S., Huck, N., Michalke, W., Wagner, E., et al. (2003). Plasma membrane H+-ATPase is involved in auxin-mediated cell elongation during wheat embryo development. Plant Physiol. 131, 1302–1312. doi: 10.1104/pp.013466
Sajjad, M., Ma, X., Khan, S. H., Shoaib, M., Song, Y., Yang, W., et al. (2017). TaFlo2-A1, an ortholog of rice Flo2, is associated with thousand grain weight in bread wheat (Triticum aestivum L.). BMC Plant Biol. 17:164. doi: 10.1186/s12870-017-1114-3
Simmonds, J., Scott, P., Leveringtonwaite, M., Turner, A. S., Brinton, J., Korzun, V., et al. (2014). Identification and independent validation of a stable yield and thousand grain weight QTL on chromosome 6A of hexaploid wheat (Triticum aestivum L.). BMC Plant Biol. 14:191. doi: 10.1186/s12870-014-0191-9
Su, C., Lu, W., Zhao, T., and Gai, J. (2010). Verification and fine-mapping of QTLs conferring days to flowering in soybean using residual heterozygous lines. Sci. Bull. 55, 499–508. doi: 10.1007/s11434-010-0032-7
Su, Z., Hao, C., Wang, L., Dong, Y., and Zhang, X. (2012). Identification and development of a functional marker of TaGW2 associated with grain weight in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 122, 211–223. doi: 10.1007/s00122-010-1437-z
Sukumaran, S., Lopes, M., Dreisigacker, S., and Reynolds, M. (2017). Genetic analysis of multi-environmental spring wheat trials identifies genomic regions for locus-specific trade-offs for grain weight and grain number. Theor. Appl. Genet. 131:999. doi: 10.1007/s00122-018-3066-x
Sun, X., Liu, D., Zhang, X., Li, W., Hui, L., Hong, W., et al. (2013). SLAF-seq: an efficient method of large-scale De novo SNP discovery and genotyping using high-throughput sequencing. PLoS One 8:e58700. doi: 10.1371/journal.pone.0058700
Sun, X. Y., Wu, K., Zhao, Y., Kong, F. M., Han, G. Z., Jiang, H. M., et al. (2009). QTL analysis of kernel shape and weight using recombinant inbred lines in wheat. Euphytica 165, 615–624. doi: 10.1007/s10681-008-9794-2
Talbert, L. E., Blake, N. K., Chee, P. W., Blake, T. K., and Magyar, G. M. (1994). Evaluation of “sequence-tagged-site” PCR products as molecular markers in wheat. Theor. Appl. Genet. 87, 789–794. doi: 10.1007/BF00221130
Thorwarth, P., Liu, G., Ebmeyer, E., Schacht, J., Schachschneider, R., Kazman, E., et al. (2019). Dissecting the genetics underlying the relationship between protein content and grain yield in a large hybrid wheat population. Theor. Appl. Genet. 132, 489–500. doi: 10.1007/s00122-018-3236-x
Torada, A., Ikeguchi, S., and Koike, M. (2005). Mapping and validation of PCR-based markers associated with a major QTL for seed dormancy in wheat. Euphytica 143, 251–255. doi: 10.1007/s10681-005-7872-2
Torada, A., Koike, M., Ogawa, T., Takenouchi, Y., Tadamura, K., Wu, J., et al. (2016). A causal gene for seed dormancy on wheat chromosome 4A encodes a MAP kinase kinase. Curr. Biol. 26, 782–787. doi: 10.1016/j.cub.2016.01.063
Vandana, J., Vijay, G., Saloni, M., Priyanka, A., Kumar, K. M., Paul, K. J., et al. (2015). Identification of novel SNP in promoter sequence of TaGW2-6A associated with grain weight and other agronomic traits in wheat (Triticum aestivum L.). PLoS One 10:e0129400. doi: 10.1371/journal.pone.0129400
Vos, P., Hogers, R., Bleeker, M., Reijans, M., and Kuiper, M. (1995). AFLPs a new technique for DNA fingerprinting. Nucleic Acids Res. 23, 4407–4414. doi: 10.1093/nar/23.21.4407
Wang, E., Wang, J., Zhu, X., Hao, W., Wang, L., Li, Q., et al. (2008). Control of rice grain-filling and yield by a gene with a potential signature of domestication. Nat. Genet. 40, 1370–1374. doi: 10.1038/ng.220
Wang, F., Zhang, J., Chen, Y., Zhang, C., Gong, J., Song, Z., et al. (2020). Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G.barbadense introgressions. Plant Biotech. J. 18, 707–720. doi: 10.1111/pbi.13237
Wang, H., Xu, S., Fan, Y., Liu, N., Zhan, W., Liu, H., et al. (2018). Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses. Plant Biotech. J. 16, 1464–1475. doi: 10.1111/pbi.12889
Wang, J., Wang, R., Mao, X., Zhang, J., Liu, Y., Xie, Q., et al. (2020). RING finger ubiquitin E3 ligase gene TaSDIR1-4A contributes to grain size in common wheat. J. Exp. Bot. 71:eraa271. doi: 10.1093/jxb/eraa271
Wang, L., Ge, H., Hao, C., Dong, Y., Zhang, X., and Hudson, M. E. (2012). Identifying loci influencing 1,000-kernel weight in wheat by microsatellite screening for evidence of selection during breeding. PLoS One 7:e29432. doi: 10.1371/journal.pone.0029432
Wang, S., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotech. J. 12, 787–796.
Wang, S. X., Zhu, Y. L., Zhang, D. X., Shao, H., and Ma, C. X. (2017). Genome-wide association study for grain yield and related traits in elite wheat varieties and advanced lines using SNP markers. PLoS One 12:e0188662. doi: 10.1371/journal.pone.0188662
Wen, Y., Zhang, Y., Zhang, J., Feng, J., Dunwell, J. M., and Zhang, Y. (2019). An efficient multi-locus mixed model framework for the detection of small and linked QTLs in F2. Brief. Bioinform. 20, 1913–1924. doi: 10.1093/bib/bby058
Wu, Q., Chen, Y., Zhou, S., Fu, L., Chen, J., Xiao, Y., et al. (2015). High-density genetic linkage map construction and QTL mapping of grain shape and size in the wheat population yanda1817 × beinong6. PLoS One 10:e0118144. doi: 10.1371/journal.pone.0118144
Xin, F., Zhu, T., Wei, S., Han, Y., Zhao, Y., Zhang, D., et al. (2020). QTL mapping of kernel traits and validation of a major QTL for kernel length-width ratio using SNP and bulked segregant analysis in wheat. Sci. Rep. 10, 1–12. doi: 10.1038/s41598-019-56979-7
Xu, Y. F., Li, S. S., Li, L. H., Ma, F. F., Fu, X. Y., Shi, Z. L., et al. (2017). QTL mapping for yield and photosynthetic related traits under different water regimes in wheat. Mol. Breeding 37:34. doi: 10.1007/s11032-016-0583-7
Yang, J., Hu, C., Hu, H., Yu, R., Xia, Z., Ye, X., et al. (2008). QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics 24, 721–723. doi: 10.1093/bioinformatics/btm494
Yang, J., Zhou, Y., Wu, Q., Chen, Y., Zhang, P., Zhang, Y., et al. (2019a). Molecular characterization of a novel TaGL3-5A allele and its association with grain length in wheat (Triticum aestivum L.). Theor. Appl. Genet. 132, 1799–1814. doi: 10.1007/s00122-019-03316-1
Yang, J., Zhou, Y., Zhang, Y., Hu, W., Wu, Q., Chen, Y., et al. (2019b). Cloning, characterization of TaGS3 and identification of allelic variation associated with kernel traits in wheat (Triticum aestivum L.). BMC Genetics 20:1–11. doi: 10.1186/s12863-019-0800-6
Yang, J., Zhu, J., and Williams, R. W. (2007). Mapping the genetic architecture of complex traits in experimental populations. Bioinformatics 23, 1527–1536. doi: 10.1093/bioinformatics/btm143
Yang, X., Nong, B., Xia, X., Zhang, Z., Zeng, Y., Liu, K., et al. (2017). Rapid identification of a new gene influencing low amylose content in rice landraces (Oryza sativa L.) using genome-wide association study with specific-locus amplified fragment sequencing. Genome 60, 465–472. doi: 10.1139/gen-2016-0104
Yang, Y., Zhao, X. L., Xia, L. Q., Chen, X. M., Xia, X. C., Yu, Z., et al. (2007). Development and validation of a Viviparous-1 STS marker for pre-harvest sprouting tolerance in Chinese wheats. Theor. Appl. Genet. 115, 971–980. doi: 10.1007/s00122-007-0624-z
Yang, Z., Bai, Z., Li, X., Wang, P., Wu, Q., Yang, L., et al. (2012). SNP identification and allelic-specific PCR markers development for TaGW2, a gene linked to wheat kernel weight. Theor. Appl. Genet. 125, 1057–1068. doi: 10.1007/s00122-012-1895-6
Yin, C., Li, H., Li, S., Xu, L., Zhao, Z., and Wang, J. (2015). Genetic dissection on rice grain shape by the two-dimensional image analysis in one japonica × indica population consisting of recombinant inbred lines. Theor. Appl. Genet. 128, 1969–1986. doi: 10.1007/s00122-015-2560-7
Ying, G., Sun, J., Zhang, G., Wang, Y., and Li, S. (2013). Haplotype, molecular marker and phenotype effects associated with mineral nutrient and grain size traits of TaGS1a in wheat. Field Crops Res. 154, 119–125. doi: 10.1016/j.fcr.2013.07.012
Yue, A., Li, A., Mao, X., Chang, X., Li, R., and Jing, R. (2015). Identification and development of a functional marker from 6-SFT-A2 associated with grain weight in wheat. Mol. Breeding 35:63. doi: 10.1007/s11032-015-0266-9
Zhang, J., Gizaw, S. A., Bossolini, E., Hegarty, J., Howell, T., Carter, A. H., et al. (2018). Identification and validation of QTL for grain yield and plant water status under contrasting water treatments in fall-sown spring wheats. Theor. Appl. Genet. 131, 1741–1759. doi: 10.1007/s00122-018-3111-9
Zhang, X., Wang, G., Dong, T., Chen, B., Du, H., Li, C., et al. (2019). High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.). BMC Plant Biol. 19:167. doi: 10.1186/s12870-019-1753-7
Zhang, Y., Liu, J., Xia, X., and He, Z. (2014). TaGS-D1, an ortholog of rice OsGS3, is associated with grain weight and grain length in common wheat. Mol. Breeding 34, 1097–1107. doi: 10.1007/s11032-014-0102-7
Zhang, Y., Xia, X., and He, Z. (2017). The seed dormancy allele TaSdr-A1a associated with pre-harvest sprouting tolerance is mainly present in Chinese wheat landraces. Theor. Appl. Genet. 130, 81–89. doi: 10.1007/s00122-016-2793-0
Zhang, Z., Lv, G., Li, B., Wang, J., Zhao, Y., Kong, F., et al. (2017). Isolation and characterization of the TaSnRK2.10 gene and its association with agronomic traits in wheat (Triticum aestivum L.). PLoS One 12:e0174425. doi: 10.1371/journal.pone.0174425
Zhang, Z., Shang, H., Shi, Y., Huang, L., Li, J., Ge, Q., et al. (2016). Construction of a high-density genetic map by specific locus amplified fragment sequencing (SLAF-seq) and its application to Quantitative Trait Loci (QTL) analysis for boll weight in upland cotton (Gossypium hirsutum.). BMC Plant Biol. 16:79. doi: 10.1186/s12870-016-0741-4
Zhu, X. F. (2016). Cloning and expression analysis of TaPTF1 and Tabas1 genes associated with thousand grain weight in common wheat (Triticum aestivum L.) [D]. in Chinese with English abstract. Anhui Agric. Univ.
Zhu, X. F., Zhang, H. P., Hu, M. J., Wu, Z. Y., Jiang, H., Cao, J. J., et al. (2016). Cloning and characterization of Tabas1-B1 gene associated with flag leaf chlorophyll content and thousand-grain weight and development of a gene-specific marker in wheat. Mol. Breeding 36:142. doi: 10.1007/s11032-016-0563-y
Keywords: common wheat, grain weight/size, QTL, SLAF-seq, GWAS, SNP
Citation: Cao J, Shang Y, Xu D, Xu K, Cheng X, Pan X, Liu X, Liu M, Gao C, Yan S, Yao H, Gao W, Lu J, Zhang H, Chang C, Xia X, Xiao S and Ma C (2020) Identification and Validation of New Stable QTLs for Grain Weight and Size by Multiple Mapping Models in Common Wheat. Front. Genet. 11:584859. doi: 10.3389/fgene.2020.584859
Received: 18 July 2020; Accepted: 21 September 2020;
Published: 11 November 2020.
Edited by:
Deepmala Sehgal, International Maize and Wheat Improvement Center, MexicoReviewed by:
Liu Jin Dong, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, ChinaHimanshu Sharma, National Agri-Food Biotechnology Institute, India
Copyright © 2020 Cao, Shang, Xu, Xu, Cheng, Pan, Liu, Liu, Gao, Yan, Yao, Gao, Lu, Zhang, Chang, Xia, Xiao and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Haiping Zhang, zhhp20@163.com; Cheng Chang, changtgw@126.com
†These authors have contributed equally to this work