- 1Key Laboratory of Wheat Biology and Genetic Improvement for Low & Middle Yangtze Valley Ministry of Agriculture, Lixiahe Agricultural Institute of Jiangsu Province, Yangzhou, China
- 2National Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, China
- 3College of Agronomy, Yangzhou University, Yangzhou, China
- 4College of Agronomy, Henan Agricultural University, Zhengzhou, China
Fusarium head blight (FHB) is a destructive wheat disease present throughout the world, and host resistance is an effective and economical strategy used to control FHB. Lack of adequate resistance resource is still a main bottleneck for FHB genetics and wheat breeding research. The synthetic-derived bread wheat line C615, which does not carry the Fhb1 gene, is a promising source of FHB resistance for breeding. A population of 198 recombinant inbred lines (RILs) produced by crossing C615 with the susceptible cultivar Yangmai 13 was evaluated for FHB response using point and spray inoculations. As the disease phenotype is frequently complicated by other agronomic traits, we used both traditional and multivariate conditional QTL mapping approaches to investigate the genetic relationships (at the individual QTL level) between FHB resistance and plant height (PH), spike compactness (SC), and days to flowering (FD). A linkage map was constructed from 3,901 polymorphic SNP markers, which covered 2,549.2 cM. Traditional and conditional QTL mapping analyses found 13 and 22 QTL for FHB, respectively; 10 were identified by both methods. Among these 10, three QTL from C615 were detected in multiple years; these QTL were located on chromosomes 2AL, 2DS, and 2DL. Conditional QTL mapping analysis indicated that, at the QTL level, SC strongly influenced FHB in point inoculation; whereas PH and SC contributed more to FHB than did FD in spray inoculation. The three stable QTL (QFhbs-jaas.2AL, QFhbp-jaas.2DS, and QFhbp-jaas.2DL) for FHB were partly affected by or were independent of the three agronomic traits. The QTL detected in this study improve our understanding of the genetic relationships between FHB response and related traits at the QTL level and provide useful information for marker-assisted selection for the improvement of FHB resistance in breeding.
Introduction
Fusarium head blight (FHB), also known as head scab, is one of the most destructive fungal diseases of wheat and is found in temperate and subtropical regions around the world. The most important fungal pathogens associated with FHB are Fusarium graminearum and Fusarium culmorum (Parry et al., 1995). FHB reduces grain yield and quality due to shriveled kernels with low weight and contamination of several mycotoxins, such as deoxynivalenol (DON). Infected kernels with excessive DON pose a severe threat to human and animal health (McMullen et al., 2012). In 2015, the Codex Alimentarius Commission (CAC) enacted standards regulating DON content, with maximum permitted levels of 2 ppm and 1 ppm for unprocessed and finished wheat products, respectively, and 0.2 ppm for baby food (Codex Alimentarius Commission [CAC], 2015). FHB outbreaks have become more serious and more frequent in recent decades, possibly due to climate change and agronomic practices (especially maize–wheat rotations). Globally, FHB causes about 10–70% yield loss in epidemic years (Zhang et al., 2011). In the Middle-Lower Yangtze River Valley Region of China, FHB causes 5–10% yield loss in most years but can cause losses of nearly 100% in severe epidemic years (Cheng et al., 2012). Host resistance is currently recognized as the most effective and environmentally friendly method of FHB management (He et al., 2013), justifying the use of resistant varieties to control FHB.
The genetics of FHB resistance are complex since this trait is under multigenic control and is subject to genotype × environment interactions. Five types of FHB resistance mechanisms have been described (Mesterhazy et al., 1999): Type I (resistance to initial infection), Type II (resistance to disease spread within infected heads), Type III (resistance to DON accumulation), Type IV (resistance to kernel damage), and Type V (tolerance). More than 250 FHB QTL have been reported, covering all 21 wheat chromosomes (Buerstmayr et al., 2009; Liu et al., 2009). Among these, Fhb1 from the Chinese variety Sumai 3 is the resistance gene most often used by breeding programs to combat FHB (Cuthbert et al., 2006; Lv et al., 2014). Fhb2 from Sumai 3 (Yang et al., 2006; Lu et al., 2010) and Fhb7 (Guo et al., 2015) from Thinopyrum ponticum have also been used by resistance breeding programs. However, it is still a formidable challenge to transfer resistance genes into susceptible cultivars because major sources of resistance (such as Sumai 3 and Wangshuibai) also carry undesirable agronomic traits (Dvorjak, 2014; Li et al., 2016).
Several studies have showed that agronomic traits may be associated with FHB resistance. Plant height (PH) and days to flowering (FD) are often negatively correlated to FHB severity (Mao et al., 2010; Lu et al., 2011; Buerstmayr et al., 2012; Giancaspro et al., 2016; McCartney et al., 2016). Positive correlations between spike compactness (SC) and FHB severity have also been commonly reported (Buerstmayr et al., 2011; Giancaspro et al., 2016). Other traits associated with FHB include heading time (Emrich et al., 2008), degree of anther extrusion (Skinnes et al., 2010; Lu et al., 2013), and presence/absence of awns (Ban and Suenaga, 2000). Among these reports, researchers have generally analyzed the genetic relationship between FHB response and agronomic traits using QTL mapping approaches that first ask whether QTL are closely linked or pleiotropic. However, traditional QTL mapping analysis of genetic relationships of complex traits is often confounded with variations involving other traits (Cui et al., 2012; Li et al., 2015). Multivariable conditional analysis has been utilized in evaluating the conditional phenotypic values of a target trait by excluding the effects of related traits as well as in determining the contribution of related traits to a target trait. By comparatively analyzing both conditional and unconditional QTL, we may be able to identify the genetic relationships among different traits at the QTL level (Zhu, 1995; Wen and Zhu, 2005). Thus far, various studies on wheat using this methodology have examined PH (Cui et al., 2011; Yu et al., 2014), kernel weight (Li et al., 2015), yield (Xu et al., 2014; Fan et al., 2015), and seedling traits (Zhang et al., 2014). Such studies indicate that conditional QTL analysis is helpful for understanding the influence of one trait on another complex trait at the QTL level. However, conditional QTL mapping analyses of the genetic relationships between FHB response and agronomic traits is limited.
The identification of new FHB-resistant germplasms has broadened the resistance gene pool as well as facilitated the improvement of resistance in cultivars, as wheat breeding has thus far relied heavily on Sumai 3 and its derivatives, with Fhb1 as donor parents. Several synthetic hexaploid wheat (SHW) lines and their derivatives developed by the International Maize and Wheat Improvement Center (CIMMYT) have been used in resistance breeding programs around the world (He et al., 2013, 2016; Zhu et al., 2016). The SHW-derived CIMMYT line C615 exhibits moderate resistance to FHB in the field, but haplotype analysis revealed that it appeared to lack Fhb1. In this study, a recombinant inbred line (RIL) population was developed from a cross between C615 and a susceptible parent, Yangmai 13. The objectives of the study were to (1) dissect the FHB resistance QTL using a high-density SNP map; (2) identify SNP markers closely linked to resistance loci; and (3) analyze the genetic relationships between FHB response and morphological traits and FD at the QTL level using multivariable conditional QTL mapping analysis.
Materials and Methods
Plant Materials
C615 is a CIMMYT synthetic derived line (kindly provided by Prof. A. Mujeeb-kazi, CIMMYT, Mexico) with moderate FHB resistance and good adaptability at Yangzhou. It has the pedigree SABUF/3/BCN//CETA/AE.SQUARROSA (895), where SABUF and BCN are CIMMYT bread wheat genotypes, and CETA is a durum variety. Yangmai 13 is a FHB-susceptible Chinese soft wheat cultivar and has been widely planted in the Middle-Lower Yangtze River region. One hundred and ninety eight F7 RILs were produced by single seed descent from cross C615/Yangmai 13. For FHB evaluation, Sumai 3, Yangmai 158, and Annong 8455 were used as resistant, moderately resistant and susceptible controls, respectively.
Phenotypic Evaluation
The RILs and parents (C615 and Yangmai 13) were evaluated for FHB response and related agronomic traits at Lixiahe Agricultural Institute of Jiangsu Province, Yangzhou, during 2014–2015 (E1), 2015–2016 (E2), and 2016–2017 (E3). Field experiments were designed as randomized complete blocks with two replicates per environment. The RILs in each replication were sown in two 133 cm rows with 40 seeds per row, with a row spacing of 25 cm. The field trials were managed following local practices.
Point Inoculation
All materials were point inoculated with four F. graminearum strains (F4, F15, F34, and F0609) kindly provided by Prof. Huaigu Chen, Jiangsu Academy of Agricultural Sciences. Inoculations were performed at the late heading stage when 5 μl of macroconidial suspension (1.0 × 105 conidia/ml) was injected into a single floret in the middle of each spike; 30 spikes were inoculated per row (Lu et al., 2011). After inoculation, the disease nursery was mist-irrigated for 5 min every half-hour from 7:00 am to 6:00 pm each day to provide high humid conditions favorable for FHB infection. FHB severities were recorded 20 days after inoculation as the number of symptomatic spikelets per infected spike, and mean data were used for analysis.
Spray Inoculation
The two-replicate field nursery for spray inoculation at flowering were grown with two replicates for 2 years (E2 and E3). Central plants in each plot were inoculated by spraying a mixed macroconidial suspension (same as for point inoculation) with a backpack sprayer and was repeated 2–3 days later (He et al., 2016). The same misting system as for point-inoculated experiments was applied. FHB severity was calculated by recording the percentage of symptomatic spikelets of 30 spikes per plot (Lu et al., 2013).
The RILs and parents were also evaluated for three agronomic traits reported to be related to FHB response, i.e., PH, SC, and FD over 3 years (E1, E2, and E3). PH was measured from the ground to the top of the spikes excluding awns. SC was calculated from the spike length (SL) and spikelet number per spike (SNS) according to the equation: SC = SNS/SL (Lv et al., 2014). PH and SC were recorded as average values of twenty individual plants per line. FD was recorded when 50% of spikes a line were at anthesis.
Genotyping and Linkage Map Construction
Genomic DNA for SNP assays was extracted from young leaf tissues by the CTAB method (Stacey and Isaac, 1994). The 198 RILs and two parents were genotyped using the wheat 90K iSelect array with 81,587 SNP (Wang et al., 2014). Genotyping assays were carried out on the Illumina iScan reader and made genotypic clusters for each SNP using GenomeStudio software 1.9.4 (Illumina; http://www.illumina.com).
Prior to mapping, SNP data were evaluated following Zhai et al. (2015). SNPs with more than 20% missing values or strong segregation distortion were excluded from linkage mapping. Linkage groups were constructed using Joinmap V4.0 (Stam, 1993) with a minimum independent logarithm of odds (LOD) threshold of 10.0. The Kosambi mapping function was used to estimate genetic distances (in cM) between markers, with a maximum recombination threshold of 0.4 and a jump threshold of 5.0 (Kosambi, 1943). Linkage maps were generated with the software MapChart 2.21 (Voorrips, 2002). The long (L) and short (S) arms of each chromosome were identified from the wheat 90K consensus SNP map (Wang et al., 2014).
Data Analysis and QTL Mapping
Analysis of variance (ANOVA) and phenotypic correlation coefficients were conducted using SAS v.9.2 software (SAS Institute Inc., Cary, NC, United States). Broad-sense heritability (h2) of each trait was calculated using h2 = σ2g/[σ2g + σ2gy/r + σ2e/ry] for multiple years, where σ2g is the estimate of genetic variance, σ2gy is the estimate of genotype × year interaction variance, σ2e is the estimate of residual error variance, r is the number of replicates per line, and y is the number of years.
Both conditional and traditional QTL analyses were conducted using composite interval mapping (CIM) using QTL Cartographer 2.52 (Wang et al., 2006). The LOD value was set at 3.0 after 1,000 permutations to declare of putative QTL. The QTL (2.5 < LOD < 3.0) were also reported for other environments when these detected in at least one environment reaching the significance level. QTL intervals were estimated that a 2-LOD fell off support interval with a confidence threshold 95% (van Ooijen, 1992).
Unconditional and conditional QTL were estimated based on phenotypic and conditional phenotypic values of the traits, respectively. Conditional phenotypic values of y(FHBjFRATs) were obtained from QGAStation 2.03 as described by Zhu (1995) and were estimated through using a two-step procedure described by Cui et al. (2012); here, FHB|FRATs refers to FHB conditioned on FHB-related agronomic traits (FRATs; e.g., FHB|PH means FHB conditioned on PH). The QTL data menu settings of QGAStation 2.0 were implemented based on a method described by Cui et al. (2012).
QTL were named as follows: Qtrait-lab designation.chromosome location-X, a number distinguishing multiple linkage groups (LGs) within the same chromosome. QTL for FHB resistance identified by point inoculation were designed as QFhbp- whereas QFhbs- was used for those identified by spray inoculation. For example, QFhbp-jaas.2B-2 indicates an FHB resistance QTL on the second LGs of chromosome 2B that was identified following point inoculation.
Results
Phenotypic Evaluation and Correlations Between FHB and Related Traits
The mean phenotypic values, RIL range, parental response, and broad-sense heritability for each trait are presented in Table 1. C615 was taller than Yangmai 13, whereas the latter had higher FHB severity, SC, and FD in all environments. The broad-sense heritabilities (h2) based on RIL mean data were 0.87 for point inoculation, 0.78 for spray inoculation, 0.97 for PH, 0.77 for SC, and 0.88 for FD. Large variations among RILs were observed for all traits. Based on the mean data across environments, the frequency distributions of each trait for the RILs were continuous and strong transgressive segregation was evident (Figure 1).
TABLE 1. Phenotypic data and broad-sense heritabilities for FHB severity and correlated traits in C615/Yangmai 13 RILs and parents across environments.
FIGURE 1. Frequency distributions for FHB, Plant height, Spike compactness, and Days to Flowering in the C615/Yangmai 13 RIL population based on the mean data across environments. C: C615, Y: Yangmai 13.
Several response variables, including FHB in both the point and spray inoculated experiments, as well as PH, SC, and FD significantly differed (P < 0.001) among genotypes, G × E (genotype × year) interactions, and years (environments), except G × E in SC (Supplementary Table S1). After spray inoculation, negative correlations between FHB severity and PH (r = -0.37, P < 0.01) and FD (r = -0.18 to -0.20, P < 0.05) were found to be significant; in contrast, point-inoculated experiments found non-significant correlations between FHB severity and PH or FD. FHB severity also showed significant positive correlations with SC (r = 0.15 to 0.31, P < 0.05) in both types of inoculated experiments. Finally, FHB and related agronomic traits showed weak correlations (i.e., r < 0.5) in both the point or spray inoculated experiments (Table 2).
TABLE 2. Correlation coefficients between FHB severity and associated agronomic traits in the RIL population derived from C615/Yangmai 13.
SNP Genotyping and Linkage Map
Among the 81,587 identified SNPs from the Illumina wheat 90K SNP chip, 6,321 SNPs (7.8%) were polymorphic between the two parents. Among these polymorphic markers, those that were unanchored or linked in small linkage groups (LGs, usually smaller than 5 cM) were excluded from further analyses. A final panel of 3,901 SNPs was used for the construction of a linkage map and for QTL mapping.
Detailed information for the linkage map is provided in Supplementary Tables S2, S3, and the linkage map showed good synteny with the recently released wheat reference genome, with the exception of chromosome 4B, 5A, and 7B (Supplementary Figure S1). The entire linkage map consisted of 34 LGs representing all 21 wheat chromosomes. Chromosomes 1D, 2B, 2D, 3D, 4A, 5A, 5B, 5D, 6D, 7A, and 7B were each assembled into two LGs, and chromosome 6A consisted of three LGs. The overall linkage map covered 2,549.2 cM, with chromosome length ranging from 30.7 cM (4D) to 179.9 cM (4A) with a mean of 121.4 cM. The A genome had 1,810 SNPs (46.4%), with a total length of 1,012.6 cM, and an SNP density of 1.79 markers/cM; the B genome included 1,441 SNPs (36.9%) covering 924.3 cM, and an SNP density of 1.56 markers/cM; and the D genome had 650 SNPs (16.7%) with a total length of 612.3 cM and an SNP density of 1.06 markers/cM. The number of SNP markers on each chromosome ranged from 6 (6B) to 582 (5B) and mean overall SNP density was 1.53 markers/cM, ranging from 0.07 (6B) to 3.64 (5B).
Traditional QTL Mapping
In the point-inoculated experiments, seven traditional QTL for FHB resistance were identified on chromosomes 1A, 2D (2), 4A, 5A, 5D, and 6A (Table 3 and Figure 2). A stable QTL on chromosome 2DL (LG, 2D-1, closest marker: TA001163-0861) was detected in all 3 years, with LOD values ranging from 4.80 to 6.39 that explained up to 11.79% of the phenotypic variation. QFhbp-jaas.1AL (closest marker: RAC875_c6338_1887) was detected in 3 years and accounted for 4.83–10.83% of the phenotypic variation. Four QTL identified in 2 years accounted for 4.86–9.28% of the phenotypic variation and were located on chromosomes 2DS, 4AL, 5AL, and 6AS, respectively. The remaining QTL on chromosome 5DL were detected only in 1 year and explained 7.29% of the phenotypic variation. Alleles increasing FHB resistance from C615 were found at four loci (2DS, 2DL, 4AL, and 5DL), and resistance was also conferred by QTL on chromosomes 1AL, 1BL, 5AL, and 6AS that were derived from the susceptible parent, Yangmai 13.
TABLE 3. Composite interval mapping for FHB resistance and agronomic traits in the C615/Yangmai 13 RIL population.
FIGURE 2. Genetic linkage maps and locations of QTL for FHB resistance in the C615/Yangmai 13 population. Red and blue regions on chromosomes indicate confidence intervals of QTL. QTL names are on the right with different colors for different traits. FHB point, FHB evaluated after point inoculation; FHB spray, FHB evaluated after spray inoculation; FHB point|FRATs, FHB point without the influence of related agronomic traits (e.g., PH, SC, and FD); FHB spray|FRATs, FHB spray without the influence of three related agronomic traits.
In the spray-inoculated experiments, six traditional QTL for FHB resistance were identified on chromosomes 1A (2), 2A, 2D, 4A, and 6A (Table 3 and Figure 2). QFhbs-jaas.1AL (closest marker: wsnp_CAP12_c2438_1180601) and QFhbs-2AL (closest marker: BS00022896_51) were detected in 2 years, accounting for 6.32–7.54% and 5.23–9.34% of the phenotypic variance, respectively. The remaining four QTL were identified in 1 year and explained 5.65–10.34% of the phenotypic variance. Alleles for FHB resistance were derived from C615 for all QTL except for the one on chromosome 1AL, where Yangmai 13 contributed the resistance allele. The QTL on chromosomes 1AL, 2DL, and 4AL detected in the spray-inoculated experiments were located within the same chromosomal regions as those detected in the point-inoculated experiments.
Conditional QTL Mapping
After point inoculation, FHB severity only exhibited a significant correlation with SC (Table 2), and three of seven putative traditional QTL for FHB were mapped using conditional QTL analysis (Table 4). QFhbp|Sc-jaas.1AL showed increased contributions to those of the corresponding traditional QTL in 3 years; Moreover, QFhbp|Sc-jaas.2DS-1 and QFhbp|Sc-jaas.2DL-1 showed increased or equal contributions to those of the corresponding traditional QTL in 3 years, and explained 9.78–11.80% and 8.52–11.65% of the phenotypic variations, respectively.
TABLE 4. Genetic analysis of FHB and associated traits at the individual QTL level in the C615/Yangmai 13 RIL population.
After spray inoculation in E2 and E3, four, three, and one of the six putative traditional FHB QTL with respect to PH, SC, and FD were unmapped, respectively (Table 4). QFhbs-jaas.1AS was identified in E2 when conditioned on FD with a reduced contribution and was also detected in E3 when conditioned on PH and FD. QFhbs-jaas.1AL was detected with reduced and equal contributions in E2 or E3, when FHB was conditioned on the three related traits. QFhbs-jaas.2AL showed equal contributions by a conditional QTL mapping analysis when conditioned on FD, and while whereas a decreased contribution in E3 when conditioned on SC. Excluding the influence of SC and FD on FHB, QFhbs-jaas.2DL-1 exhibited equal contributions in E3 to those of traditional QTL. QFhbs-jaas.4AL-1 explained 5.36% of the phenotypic variation when FHB was conditioned on FD in E2 and showed a contribution equal to that of the corresponding traditional QTL.
Two and three additional QTL for FHB conditioned on SC after point and spray inoculation, were detected only by conditional QTL mapping analysis, respectively (Table 5 and Figure 2). Among the five additional conditional QTL, three were detected with negative additive values, indicating that the favorable alleles were from the FHB-resistant parent C615.
TABLE 5. Extra conditional QTL for FHB resistance with respect to the related agronomic traits in the C615/Yangmai 13 RIL population.
Discussion
Linkage Map Construction
High-density linkage maps are required for genetic studies on common wheat and its large, complex genome (Peleg et al., 2008; Zhai et al., 2016). Single nucleotide polymorphisms (SNPs) that are widely distributed throughout the genome have been used in QTL mapping (McCartney et al., 2016; Cui et al., 2017). In the current study, QTL mapping was conducted using the wheat 90K chip (Wang et al., 2014), and 6,321 polymorphic SNPs markers were allocated to 34 LGs. One reason for the high number of LGs is that the construction of this linkage map adopted a high LOD (LOD > 10). Another reason might be that some chromosome segments in bi-parents consisted of numerous monomorphic SNPs, and the original SNPs in the 90K chip were unevenly distributed among chromosomes. Comparable results were reported by Giancaspro et al. (2016) and He et al. (2016). Among these polymorphic SNPs, 650 SNPs (16.7%) were located in the D genome (especially on 4D, 5D, 6D, and 7D), much less abundant than those for the A (46.4%) and B (36.9%) genomes (Supplementary Tables S2, S3). The low coverage of the D genome was in agreement with previous reports (Wu et al., 2015; Zou et al., 2016).
Traditional and Conditional QTL for FHB Resistance
In this study, conditional QTL mapping was used to analyze FHB response conditioned on related agronomic traits, including PH, SC, and FD at the individual QTL level. The conditional QTL could be divided into four types by comparing the genetic effects of traditional QTL (Cui et al., 2011; Li et al., 2015). When our analysis of QTL for FHB resistance conditioned on SC, four different outcomes are generated: (1) a QTL was (for example, QFhbp-jaas.4AL-1 in E1 and E2) detected only in traditional QTL analysis, indicating that this QTL is completely influenced by SC; (2) a QTL (for example, QFhbs-jaas.2DL-1 in E3) was identified both in traditional and conditional QTL analysis with a very similar effect, showing that this QTL only improves FHB resistance and is not influenced by SC; (3) a QTL for FHB resistance showed either a greatly reduced or enhanced effect, implying that this QTL is partly influenced by SC. For example, QFhbp-jaas.2DS-1 exhibited enhanced effects in E1 and E2; (4) an “extra” QTL (for example, QFhbp|Sc-jaas.7BL-1 in E3) for FHB resistance was detected only in conditional QTL analysis, indicating that this QTL is completely masked by SC. Hence, the extra QTL was detected by its effect on FHB response, but has an opposite effect from SC.
Among the seven traditional QTL for FHB resistance detected following point inoculation, four, two, and one were completely, partly, and partly or not influenced by SC, respectively (Table 4). However, after spray inoculation, four, three, and one of six traditional QTL were completely influenced by PH, SC, and FD, respectively; and zero, zero, and three were independent of PH, SC, and FD (Table 4). At the QTL level, these results indicated that SC strongly influences FHB resistance after point inoculation. However, after spray inoculation, the FHB response was affected by the three agronomic traits and both PH and SC had greater contributions to FHB resistance than did FD. This was probably due to the fact that the interior of one spike was injected with a fixed amount of inoculum. Similar results were also found in correlation analysis between FHB response and agronomic traits (Table 2). To date, some hypotheses have been reported to explain these relationships (Klahr et al., 2007; Lu et al., 2011; Giancaspro et al., 2016; Malihipour et al., 2016). In general, spikes of higher lines may dry faster and be less infected by disease. The earlier- or later-flowering lines may be in the special environments (e.g., low humidity and temperature) that are not suitable for disease development. However, these traits are very undesirable in wheat breeding. Using the conditional QTL analysis to evaluate their relationships at the QTL level, we may resolve this contradiction through selecting appropriate resistance QTL. We also found that five additional conditional QTL for FHB resistance were entirely suppressed by SC (Table 5). These findings suggest that conditional QTL analysis reduce confounding QTL and other traits and thereby facilitates in the elucidation of the genetic mechanism underlying FHB resistance.
Comparison to Previous Studies
Thirteen QTL for FHB resistance were identified in this study (Table 3 and Figure 2). These QTL mainly from C615, but a few were from Yangmai 13 and accounted for low proportion of phenotypic variance (<12%). This low phenotypic variance might be affected by the quality of inoculation and the environments (Buerstmayr et al., 2009; He et al., 2016; Zhu et al., 2016). Among these, the resistance QTL on chromosome 2AL, 2DS, and 2DL were detected in both traditional and conditional QTL mapping, respectively.
The major QTL on chromosome 2DL was detected in both point and spray inoculated experiments across 3 years; this QTL was also closely linked to the SNP marker TA001163-0861. He et al. (2016) detected a QTL for FHB resistance on chromosome 2DL from CIMMYT line Soru#1, and in that case the QTL was in the marker interval Kukri_c36639_186–Xgwm539 and explained 14–20% of the phenotypic variation. SNP markers TA001163-0861 and Kukri_c36639_186 were on the 90K consensus map with a genetic distance of 5.4 cM corresponding to a physical interval of 21 Mb in the Chinese Spring RefSeq v1.0 sequence (Wang et al., 2014). In addition, Lu et al. (2013) identified a QTL on chromosome 2DL for Type I and Type II resistance from Chinese line Shanghai-3; here the QTL was tightly linked to SSR marker Xgwm539. Considering that both C615 and Soru#1 have Shanghai-3 in their pedigrees, these two QTL are likely the same. Additionally, this QTL that is responsible for different types of resistance to FHB has been detected in several other FHB-resistant germplasms, including Wuhan-1 (Somers et al., 2003), Wangshuibai (Lin et al., 2004), CJ9306 (Jiang et al., 2007), VA00W-38 (Liu et al., 2012), and SYN1 (Zhu et al., 2016). These reports collectively indicate that this QTL might be an important “true QTL”, with great potential for marker-assisted selection.
The QTL on chromosome 2AL found after spray inoculation was linked to SNP marker BS00022896_51 and explained 5.23–9.34% of the phenotypic variation. Three QTL clusters for FHB resistance were previously detected on chromosomes 2AS (2) and 2AL, and these were closely linked to SSR markers Xbarc124, Xgwm122, and Xgwm311, respectively (Gervais et al., 2003; Ma et al., 2006; Holzapfel et al., 2008; Lu et al., 2013; Zhang et al., 2014; Giancaspro et al., 2016; Petersen et al., 2016, 2017). Markers Xgwm311 (at 772 Mb) and BS00022896_51 (at 612 Mb), separated by a physical distance of 160 Mb, are in different deletion bins (Qi et al., 2004). Hence, QFhbs-jaas.2AL is likely a novel QTL. To our knowledge, only one QTL on chromosome 2DS was found in previous studies (Liu et al., 2009; Buerstmayr et al., 2011; Cai and Bai, 2014; McCartney et al., 2016) and was closely linked to the SSR marker Xgwm261 (close to Rht8) at about 20 Mb. Here, QFhbp-jaas.2DS-1 for Type II resistance was closely linked to the SNP marker Kukri_c60627_74 at 74 Mb. The SSR and SNP markers were separated by a physical distance of 44 Mb. This QTL was in a similar position to a QTL previously reported for anther extrusion in a German cultivar (He et al., 2016), which was not linked to marker Xgwm261. Thus, QFhbp-jaas.2DS-1 is likely to be new, but these two resistance QTL should be verified by future work.
Implications for Wheat Breeding
Resistance to FHB in wheat is a complex trait and marker-assisted selection is a valuable tool to improve FHB resistance. In the present study, the three stable resistance QTL from C615 on chromosomes 2AL, 2DS, and 2DL reported here belong to the second or third type of conditional QTL. These QTL were partly affected by or were independent of three related agronomic traits, and thus should be given more consideration for developing resistant varieties with good agronomic traits. Indications are that when the three resistance QTL are combined, these can reduce susceptibility by 30–40% (Supplementary Table S4). These SNP markers closely linked to the QTL detected here may be effectively used in marker-assisted selection for improving FHB resistance when C615 is used as resistance donor. The results of this study also indicate that conditional QTL mapping analysis can improve our understanding of complex traits.
Ethics Statement
We declare that these experiments comply with the ethical standards in China.
Author Contributions
YZ and SCheng conceived and designed the experiments. XY and JC performed the experiments. XY and ZJ analyzed the data. WH, TB, DG, DL, RW, YL, SChen, XC, and JL contributed reagents, materials, or analysis tools. XY wrote the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (2017YFD0100801), the Natural Science Foundation of Jiangsu Province (BK20171279), the China Agricultural Research System (CARS-03-03B), the Independent Innovation Fund for Agricultural Science and Technology in Jiangsu Province (CX (14) 5080), the Key R&D Program in Jiangsu Province (BE2015352-3).
Conflict of Interest Statement
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.
Acknowledgments
The authors are grateful to Prof. Robert A. McIntosh (Plant Breeding Institute, University of Sydney) and Prof. Tao Li (College of Agronomy, Yangzhou University), for critically reviewing this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00573/full#supplementary-material
Footnotes
- ^ http://www.earthatlas.mapchart.com
- ^ http://statgen.ncsu.edu/qtlcart/WQTLCart.htm
- ^ http://ibi.zju.edu.cn/software/qga/
References
Ban, T., and Suenaga, K. (2000). Genetic analysis of resistance to Fusarium head blight caused by Fusarium graminearum in Chinese wheat cultivar Sumai 3 and the Japanese cultivar Saikai 165. Euphytica 113, 87–99. doi: 10.1023/A:1003951509797
Buerstmayr, H., Ban, T., and Anderson, J. A. (2009). QTL mapping and marker assisted selection for Fusarium head blight resistance in wheat: a review. Plant Breed. 128, 1–26. doi: 10.1111/j.1439-0523.2008.01550.x
Buerstmayr, M., Huber, K., Heckmann, J., Steiner, B., Nelson, J. C., and Buerstmayr, H. (2012). Mapping of QTL for Fusarium head blight resistance and morphological and developmental traits in three backcross populations derived from Triticum dicoccum × Triticum durum. Theor. Appl. Genet. 125, 1751–1765. doi: 10.1007/s00122-012-1951-2
Buerstmayr, M., Lemmens, M., Steiner, B., and Buerstmayr, H. (2011). Advanced backcross QTL mapping of resistance to Fusarium head blight and plant morphological traits in a Triticum macha × T. aestivum population. Theor. Appl. Genet. 123, 293–306. doi: 10.1007/s00122-011-1584-x
Cai, J., and Bai, G. H. (2014). Quantitative trait loci for Fusarium head blight resistance in Huangcandou × ‘Jagger’ wheat population. Crop Sci. 54, 2520–2528. doi: 10.2135/cropsci2013.12.0835
Cheng, S. H., Zhang, Y., Bie, T. D., Gao, D. R., and Zhang, B. Q. (2012). Damage of wheat Fusarium head blight epidemics and genetic improvement of wheat for scab resistance in China. Jiangsu J. Agric. Sci. 28, 938–942.
Codex Alimentarius Commission [CAC] (2015). Codex Alimentarius Commission, Codex General Standard for Contaminants and Toxins in Food and Feed (Codex Stan 193–1995). 1–44. http://www.fao.org/fao-who-codexalimentarius/codex-home/en/
Cui, F., Li, J., Ding, A. M., Zhao, C. H., Wang, L., Wang, X. Q., et al. (2011). Conditional QTL mapping for plant height with respect to the length of the spike and internode in two mapping populations of wheat. Theor. Appl. Genet. 122, 1517–1536. doi: 10.1007/s00122-011-1551-6
Cui, F., Zhang, N., Fan, X. L., Zhang, W., Zhao, C. H., Yang, L. J., et al. (2017). Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number. Sci. Rep. 7:3788. doi: 10.1038/s41598-017-04028-6
Cui, F., Zhao, C. H., Li, J., Ding, A. M., Li, X. F., Bao, Y. G., et al. (2012). Kernel weight per spike: what contributes to it at the individual QTL level? Mol. Breed. 31, 265–278. doi: 10.1007/s11032-012-9786-8
Cuthbert, P. A., Somers, D. J., Thomas, J., Cloutier, S., and Brulé-Babel, A. (2006). Fine mapping Fhb1, a major gene controlling Fusarium head blight resistance in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 112, 1465–1472. doi: 10.1007/s00122-006-0249-7
Dvorjak, D. S. (2014). Fusarium Head Blight Resistance and Agronomic Performance in Soft Red Winter Wheat Populations. Ph.D. dissertation, University of Kentucky, Lexington, KY.
Emrich, K., Wilde, F., Miedaner, T., and Piepho, H. P. (2008). REML approach for adjusting the Fusarium head blight rating to a phenological date in inoculated selection experiments of wheat. Theor. Appl. Genet. 117, 65–73. doi: 10.1007/s00122-008-0753-z
Fan, X. L., Cui, F., Zhao, C. H., Zhang, W., Yang, L. J., Zhao, X. Q., et al. (2015). QTL for flag leaf size and their influence on yield-related traits in wheat (Triticum aestivum L.). Mol. Breed. 35:24. doi: 10.1007/s11032-015-0205-9
Gervais, L., Dedryver, F., Morlais, J. Y., Bodusseau, V., Negre, S., Bilous, M., et al. (2003). Mapping of quantitative trait loci for field resistance to Fusarium head blight in an European winter wheat. Theor. Appl. Genet. 106, 961–970. doi: 10.1007/s00122-002-1160-5
Giancaspro, A., Giove, S. L., Zito, D., Blanco, A., and Gadaleta, A. (2016). Mapping QTLs for Fusarium head blight resistance in an interspecific wheat population. Front. Plant Sci. 7:1381. doi: 10.3389/fpls.2016.01381
Guo, J., Zhang, X. L., Hou, Y. L., Cai, J. J., Shen, X. R., Zhou, T. T., et al. (2015). High-density mapping of the major FHB resistance gene Fhb7 derived from Thinopyrum ponticum and its pyramiding with Fhb1 by marker-assisted selection. Theor. Appl. Genet. 128, 2301–2316. doi: 10.1007/s00122-015-2586-x
He, X. Y., Lillemo, M., Shi, J. R., Wu, J. R., Bjørnstad,Å., Belova, T., et al. (2016). QTL characterization of Fusarium head blight resistance in CIMMYT bread wheat line Soru#1. PLoS One 11:e0158052. doi: 10.1371/journal.pone.0158052
He, X. Y., Singh, P. K., Duveiller, E., Dreisigacker, S., and Singh, R. P. (2013). “Development and characterization of international maize and wheat improvement center (CIMMYT) germplasm for Fusarium head blight resistance,” in Fusarium Head Blight in Latin America, eds T. Alconada Magliano and S. Chulze (Dordrecht: Springer), 241–262.
Holzapfel, J., Mohler, V., Haeberle, J., Schweizer, G., Miedaner, T., Voss, H. H., et al. (2008). “Genome distribution of QTL for Fusarium head blight resistance in European wheat germplasm,” in Proceedings of the 11th International Wheat Genetics Symposium, Brisbane, QLD, Australia, 24-29 August 2008, eds R. Appels, R. Eastwood, E. Lagudah, P. Langridge, M. Mackay, L. McIntyre, et al. (Sydney: Sydney University Press).
Jiang, G. L., Shi, J. R., and Ward, R. W. (2007). QTL analysis of resistance to Fusarium head blight in the novel wheat germplasm CJ 9306. I. Resistance to fungal spread. Theor. Appl. Genet. 116, 3–13. doi: 10.1007/s00122-007-0641-y
Klahr, A., Zimmermann, G., Wenzel, G., and Mohler, V. (2007). Effects of environment, disease progress, plant height and heading date on the detection of QTLs for resistance to Fusarium head blight in an European winter wheat cross. Euphytica 154, 17–28. doi: 10.1007/s10681-006-9264-7
Kosambi, D. D. (1943). The estimation of map distances from recombination values. Ann. Eugen. 12, 172–175. doi: 10.1111/j.1469-1809.1943.tb02321.x
Li, Q. F., Zhang, Y., Liu, T. T., Wang, F. F., Liu, K., Chen, J. S., et al. (2015). Genetic analysis of kernel weight and kernel size in wheat (Triticum aestivum L.) using unconditional and conditional QTL mapping. Mol. Breed. 35, 1–15. doi: 10.1007/s11032-015-0384-4
Li, T., Zhang, D. D., Zhou, X. L., Bai, G. H., Li, L., and Gu, S. L. (2016). Fusarium head blight resistance loci in a stratified population of wheat landraces and varieties. Euphytica 207, 551–561. doi: 10.1007/s10681-015-1539-4
Lin, F., Kong, Z. X., Zhu, H. L., Xue, S. L., Wu, J. Z., Tian, D. G., et al. (2004). Mapping QTL associated with resistance to Fusarium head blight in the Nanda2419 x Wangshuibai population. I. Type II resistance. Theor. Appl. Genet. 112, 528–535. doi: 10.1007/s00122-004-1772-z
Liu, S., Christopher, M. D., Griffey, C. A., Hall, M. D., Gundrum, P. G., and Brooks, W. S. (2012). Molecular characterization of resistance to Fusarium head blight in U.S. soft red winter wheat breeding line VA00W-38. Crop Sci. 52, 2283–2292. doi: 10.2135/cropsci2012.03.0144
Liu, S. Y., Hall, M. D., Griffey, C. A., and McKendry, A. L. (2009). Meta-analysis of QTL associated with fusarium head blight resistance in wheat. Crop Sci. 49, 1955–1968. doi: 10.2135/cropsci2009.03.0115
Lu, C. B., Cheng, S. H., Wu, R. L., Hu, Y. H., Fan, J. P., Wang, C. S., et al. (2010). Breeding the lines with resistance to Fusarium head blight of wheat cultivar Yangmai 13 by molecular marker-assisted selection. J. Triticeae Crops 30, 1058–1064.
Lu, Q. X., Lillemo, M., Skinnes, H., He, X. Y., Shi, J. R., Ji, F., et al. (2013). Anther extrusion and plant height are associated with Type I resistance to Fusarium head blight in bread wheat line ’Shanghai-3/Catbird’. Theor. Appl. Genet. 126, 317–334. doi: 10.1007/s00122-012-1981-9
Lu, Q. X., Szabo-Hever, A., Bjørnstad,Å., Lillemo, M., Semagn, K., Mesterhazy, A., et al. (2011). Two major resistance quantitative trait loci are required to counteract the increased susceptibility to Fusarium head blight of the Rht-D1b dwarfing gene in wheat. Crop Sci. 51, 2430–2438. doi: 10.2135/cropsci2010.12.0671
Lv, C., Song, Y. X., Gao, L. F., Yao, Q., Zhou, R. H., Xu, R. G., et al. (2014). Integration of QTL detection and marker assisted selection for improving resistance to Fusarium head blight and important agronomic traits in wheat. Crop J. 2, 70–78. doi: 10.1016/j.cj.2013.10.004
Ma, H. X., Zhang, K. M., Gao, L., Bai, G. H., Chen, H. G., Cai, Z. X., et al. (2006). Quantitative trait loci for resistance to Fusarium head blight and deoxynivalenol accumulation in Wangshuibai wheat under field conditions. Plant Pathol. 55, 739–745. doi: 10.1111/j.1365-3059.2006.01447.x
Malihipour, A., Gilbert, J., Fedak, G., Brûlé-Babel, A., and Cao, W. G. (2016). Characterization of agronomic traits in a population of wheat derived from Triticum timopheevii, and their association with Fusarium head blight. Eur. J. Plant Pathol. 144, 31–43. doi: 10.1007/s10658-015-0744-2
Mao, S. L., Wei, Y. M., Cao, W., Lan, X. J., Yu, M., Chen, Z. M., et al. (2010). Confirmation of the relationship between plant height and Fusarium head blight resistance in wheat (Triticum aestivum L.) by QTL meta-analysis. Euphytica 174, 343–356. doi: 10.1007/s10681-010-0128-9
McCartney, C. A., Brûlé-Babel, A. L., Fedak, G., Martin, R. A., Mccallum, B. D., Gilbert, J., et al. (2016). Fusarium head blight resistance QTL in the spring wheat cross Kenyon/86ISMN 2137. Front. Microbiol. 7:1542. doi: 10.3389/fmicb.2016.01542
McMullen, M., Bergstrom, G., Wolf, E. D., Dill-Macky, R., Hershman, D., Shaner, G., et al. (2012). A unified effort to fight an enemy of wheat and barley: Fusarium head blight. Plant Dis. 96, 1712. doi: 10.1094/PDIS-03-12-0291-FE
Mesterhazy, A., Bartok, T., Mirocha, C. G., and Komoroczy, R. (1999). Nature of wheat resistance to Fusarium head blight and the role of deoxynivalenol for breeding. Plant Breed. 118, 97–110. doi: 10.1046/j.1439-0523.1999.118002097.x
Parry, D. W., Jenkinson, P., and McLeod, L. (1995). Fusarium ear blight (scab) in small grain cereals–a review. Plant Pathol. 44, 207–238. doi: 10.1111/j.1365-3059.1995.tb02773.x
Peleg, Z., Saranga, Y., Suprunova, T., Ronin, Y., Röder, M. S., Kilian, A., et al. (2008). High-density genetic map of durum wheat × wild emmer wheat based on SSR and DArT markers. Theor. Appl. Genet. 117, 103–115. doi: 10.1007/s00122-008-0756-9
Petersen, S., Lyerly, J. H., Maloney, P. V., Brown-Guedira, G., Cowger, C., Costa, J. M., et al. (2016). Mapping of Fusarium head blight resistance quantitative trait loci in winter wheat cultivar NC-Neuse. Crop Sci. 56, 1473–1483. doi: 10.2135/cropsci2015.05.0312
Petersen, S., Lyerly, J. H., Mckendry, A. L., Islam, M. S., Brown-Guedira, G., Cowger, C., et al. (2017). Validation of fusarium head blight resistance QTL in US winter wheat. Crop Sci. 57, 1–12. doi: 10.2135/cropsci2015.07.0415
Qi, L. L., Echalier, B., Chao, S., Lazo, G. R., Butler, G. E., Anderson, O. D., et al. (2004). A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat. Genetics 168, 701–712. doi: 10.1534/genetics.104.034868
Skinnes, H., Semagn, K., Tarkegne, Y., Marøy, A. G., and Bjørnstad, Å. (2010). The inheritance of anther extrusion in hexaploid wheat and its relationship to Fusarium head blight resistance and deoxynivalenol content. Plant Breed. 129, 149–155. doi: 10.1111/j.1439-0523.2009.01731.x
Somers, D. J., Fedak, G., and Savard, M. (2003). Molecular mapping of novel genes controlling Fusarium head blight resistance and deoxynivalenol accumulation in spring wheat. Genome 46, 555–564. doi: 10.1139/g03-033
Stam, P. (1993). Construction of integrated genetic linkage maps by means of a new computer package: join Map. Plant J. 3, 739–744. doi: 10.1111/j.1365-313X.1993.00739.x
van Ooijen, J. W. (1992). Accuracy of mapping quantitative trait loci in autogamous species. Theor. Appl. Genet. 84, 803–811. doi: 10.1007/BF00227388
Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTL. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77
Wang, S., Basten, C. J., and Zeng, Z. B. (2006). Windows QTL Cartographer 2.5. Department of Statistics. Raleigh, NC: North Carolina State University.
Wang, S., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat genomic diversity using the high-density 90,000 SNP array. Plant Biotechnol. J. 12, 787–796. doi: 10.1111/pbi.12183
Wen, Y. X., and Zhu, J. (2005). Multivariable conditional analysis for complex trait and its components. Acta Genet. Sin. 32, 289–296.
Wu, Q. H., Chen, Y. X., Zhou, S. H., Fu, L., Chen, J. 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
Xu, Y. F., Wang, R. F., Tong, Y. P., Zhao, H. T., Xie, Q. G., Liu, D. C., et al. (2014). Mapping QTL for yield and nitrogen-related traits in wheat: influence of nitrogen and phosphorus fertilization on QTL expression. Theor. Appl. Genet. 127, 59–72. doi: 10.1007/s00122-013-2201-y
Yang, Z., Gilbert, J., and Procunier, J. D. (2006). Genetic diversity of resistance genes controlling fusarium head blight with simple sequence repeat markers in thirty-six wheat accessions from east asian origin. Euphytica 148, 345–352. doi: 10.1007/s10681-005-9047-6
Yu, M., Mao, S. L., Chen, G. Y., Pu, Z. E., Wei, Y. M., and Zheng, Y. L. (2014). QTL for uppermost internode and spike length in two wheat RIL populations and their affect upon plant height at an individual QTL level. Euphytica 200, 95–108. doi: 10.1007/s10681-014-1156-7
Zhai, H. J., Feng, Z. Y., Li, J., Liu, X. Y., Xiao, S. H., Ni, Z. F., et al. (2016). QTL analysis of spike morphological traits and plant height in winter wheat (Triticum aestivum L.) using a high-density SNP and SSR-based linkage map. Front. Plant Sci. 7:1617. doi: 10.3389/fpls.2016.01617
Zhai, S. N., He, Z. H., Wen, W. E., Jin, H., Liu, J. D., Zhang, Y., et al. (2015). Genome-wide linkage mapping of flour color-related traits and polyphenol oxidase activity in common wheat. Theor. Appl. Genet. 129, 377–394. doi: 10.1007/s00122-015-2634-6
Zhang, L., Luo, P. G., Ren, Z. L., and Zhang, H. Y. (2011). Controlling fusarium head blight of wheat (Triticum aestivum L.) with genetics. Adv. Biosci. Biotechnol. 2, 263–270. doi: 10.4236/abb.2011.24038
Zhang, Q., Axtman, J. E., Faris, J. D., Chao, S., Zhang, Z., Friesen, T. L., et al. (2014). Identification and molecular mapping of quantitative trait loci for Fusarium head blight resistance in emmer and durum wheat using a single nucleotide polymorphism-based linkage map. Mol. Breed. 34, 1677–1687. doi: 10.1007/s11032-014-0180-6
Zhu, J. (1995). Analysis of conditional genetic effects and variance components in developmental genetics. Genetics 141, 1633–1639.
Zhu, Z. W., Bonnett, D., Ellis, M., He, X. Y., Heslot, N., Dreisigacker, S., et al. (2016). Characterization of Fusarium head blight resistance in a CIMMYT synthetic-derived bread wheat line. Euphytica 208, 1–9. doi: 10.1007/s10681-015-1612-z
Keywords: conditional QTL analysis, Fusarium head blight, SNP marker, traditional QTL analysis, Triticum aestivum
Citation: Yi X, Cheng J, Jiang Z, Hu W, Bie T, Gao D, Li D, Wu R, Li Y, Chen S, Cheng X, Liu J, Zhang Y and Cheng S (2018) Genetic Analysis of Fusarium Head Blight Resistance in CIMMYT Bread Wheat Line C615 Using Traditional and Conditional QTL Mapping. Front. Plant Sci. 9:573. doi: 10.3389/fpls.2018.00573
Received: 11 December 2017; Accepted: 11 April 2018;
Published: 01 May 2018.
Edited by:
Chengdao Li, Murdoch University, AustraliaReviewed by:
Yongqing Jiao, Oil Crops Research Institute (CAAS), ChinaYan Long, Chinese Academy of Agricultural Sciences, China
Copyright © 2018 Yi, Cheng, Jiang, Hu, Bie, Gao, Li, Wu, Li, Chen, Cheng, Liu, Zhang and Cheng. 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 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: Yong Zhang, zy@wheat.org.cn Shunhe Cheng, yzcsh1939@126.com
†These authors have contributed equally to this work.