Yu Chen1Guodong Liu1Hehuan Ma1,2Zhangqiang Song1Chuanyun Zhang1Jingxia Zhang1Junhao Zhang3Furong Wang1,2*
Jun Zhang1,2*- 1Key Laboratory of Cotton Breeding and Cultivation in Huang-Huai-Hai Plain, Ministry of Agriculture, Cotton Research Center of Shandong Academy of Agricultural Sciences, Jinan, China
- 2College of Life Sciences, Shandong Normal University, Jinan, China
- 3College of Agriculture, Nanjing Agricultural University, Nanjing, China
The improvement of fiber quality is an essential goal in cotton breeding. In our previous studies, several quantitative trait loci (QTLs) contributing to improved fiber quality were identified in different introgressed chromosomal regions from Sea Island cotton (Gossypium barbadense L.) in a primary introgression population (Pop. A) of upland cotton (G. hirsutum L.). In the present study, to finely map introgressed major QTLs and accurately dissect the genetic contribution of the target introgressed chromosomal segments, we backcrossed two selected recombinant inbred lines (RILs) that presented desirable high fiber quality with their high lint-yielding recurrent parent to ultimately develop two secondary mapping populations (Pop. B and Pop. C). Totals of 20 and 27 QTLs for fiber quality were detected in Pop. B and Pop. C, respectively, including four and five for fiber length, four and eight for fiber micronaire, two and four for fiber uniformity, five and four for fiber elongation, and six and four for fiber strength, respectively. Two QTLs for lint percentage were detected only in Pop. C. In addition, seven stable QTLs were identified, including two for both fiber length and fiber strength and three for fiber elongation. Five QTL clusters for fiber quality were identified in the introgressed chromosomal regions, and negative effects of these chromosomal regions on lint percentage (a major lint yield parameter) were not observed. Candidate genes with a QTL-cluster associated with fiber strength and fiber length in the introgressed region of Chr.7 were further identified. The results may be helpful for revealing the genetic basis of superior fiber quality contributed by introgressed alleles from G. barbadense. Possible strategies involving marker-assisted selection (MAS) for simultaneously improving upland cotton fiber quality and lint yield in breeding programs was also discussed.
Introduction
Cotton (Gossypium spp.) is the leading worldwide textile fiber crop and provides the majority of the raw materials for the textile industry. As both the textile industry rapidly develops and the demand for textile products continues to diversify, the high requirements for fiber quality continue to increase; therefore, the improvement of cotton fiber quality is very important (Kohel et al., 2001). The cotton genus (Gossypium) comprises 45 diploid (2n = 26) and six tetraploid (2n = 52) species (Grover et al., 2015). Among the four cultivated Gossypium species worldwide, the American allotetraploid species (G. hirsutum L. and G. barbadense L.) dominate worldwide cotton production and have nearly replaced the Old World diploid cultivars (G. arboreum L. and G. herbaceum L.). Owing to its high lint yield and broad adaptability, the upland cotton allotetraploid species, G. hirsutum L. (AD)1, is the most widely cultivated cotton species and is of primary economic importance in modern worldwide cotton production (Mehboob-ur-Rahman et al., 2012; Zhang et al., 2015).
The genetic diversity of cultivated upland cotton is limited (May et al., 1995; Mehboob-ur-Rahman et al., 2012; Tyagi et al., 2014). Natural variation among wild and cultivated relatives of crops is an under-exploited resource in plant breeding. The introgression of such genetic components or genes into the background of elite crop lines provides plant breeders an important opportunity to improve the agricultural performance of modern crop varieties. Conventional breeding efforts via interspecific hybridization have resulted in the transfer of specific genes and, consequently, useful traits such as favorable fiber quality and biotic or abiotic tolerances and so on into the upland cotton gene pool, such as G. harknessii (Meyer, 1973), G. thurberi (Culp and Harrell, 1973; Culp et al., 1979), G. stockii (Nazeer et al., 2014), G. longicalyx (Nacoulima et al., 2012) and G. arboreum (Ahmad et al., 2011; Chandra and Sreenivasan, 2011; Iqbal et al., 2015). In particular, germplasm with improved fiber quality conferred by chromosomal introgressions from other species within the genus Gossypium has been developed; one of those species is another cultivated allotetraploid cotton, Sea Island cotton (G. barbadense L.), whose lint yield is much lower than that of G. hirsutum L. but whose fiber quality is markedly superior. The developed germplasms have been used in breeding to improve the fiber quality of upland cotton (Mergeai, 2003; Zhou et al., 2003; Smith et al., 2008; Cao et al., 2014). However, due to the negative genetic correlation between fiber quality and lint yield (Zhou et al., 2003; Shen et al., 2007; Shang et al., 2015), the use of introgressed germplasm to improve fiber quality and lint yield simultaneously in upland cotton breeding is difficult (Meredith, 2005; Wang et al., 2011, 2013).
Mapping quantitative trait loci (QTLs) for fiber quality and yield have facilitated molecular marker-assisted selection (MAS) to simultaneously improve cotton yield and fiber quality (Duran et al., 2009); as such, the identification of stable QTLs is necessary to use MAS to improve cotton fiber quality. More than 1000 QTLs involved in cotton fiber quality have been identified (Said et al., 2013; Yu et al., 2014). However, most research is based on a mapping population developed by the interspecific hybridization of G. hirsutum × G. barbadense (Jiang et al., 1998; Zhang et al., 2002; Lacape et al., 2003, 2009; Nguyen et al., 2004; Guo et al., 2007; Yu et al., 2007; Zhao et al., 2012); therefore, the use of these QTLs for MAS in upland cotton breeding is extremely limited (Liang et al., 2013). Some introgressed alleles for fiber quality improvement from wild cotton species such as G. tomentosum (Zhang et al., 2011) and G. klotzschianum (Xu et al., 2012) in upland cotton have recently been identified. However, most of the studies on QTL mapping of fibre quality in which an intraspecific population of G. hirsutum introgressed with chromosomal segments from other species of Gossypium were used (Shen et al., 2007; Zhang et al., 2009, 2012; Sun et al., 2012; Li et al., 2016; Ma et al., 2017; Tan et al., 2018) did not pay close attention to the contributions of the introgressed genetic components (or alleles). Genetic dissection of the introgressed alleles for fiber quality improvement from other Gossypium species, especially from cultivated allotetraploid species with superior fiber quality, i.e., G. barbadense L, would be beneficial to MAS in upland cotton breeding with respect to the simultaneous improvement of fiber quality and lint yield.
In our previous study, using an introgression germplasm with superior fiber quality, we identified components introgressed from G. barbadense into upland cotton at a genome-wide level, and found that the majority of favorable alleles for fiber quality traits were derived from the introgression genomic components (Wang et al., 2011). By using F2 and F2:3 mapping populations, we further detected QTLs related to fiber quality and lint percentage and found that most of the favorable alleles were located on the introgressed chromosomal segments on Chr. 2, Chr. 3, Chr. 7, Chr. 16, Chr. 19, Chr. 23, and Chr. 25. In addition, major QTLs (qFS-C7-1, qFL-C16-1, qFL-C25-1, etc.) were subsequently detected (Wang et al., 2013). We had also constructed an RIL population for further searching stable QTLs for fiber quality improvement by MAS (Wang et al., 2016).
The objectives of this study were as follows: (1) Develop secondary introgression mapping population for anchoring the target introgressed chromosomal regions by backcrossing two selected RILs (Wang et al., 2016) that present desirable higher fiber quality with their high lint-yielding recurrent parent; (2) Accurately dissect the genetic contribution of the target introgressed chromosomal regions to improved fiber quality and consequently, investigate the genetic basis of the formation of excellent fiber quality traits in upland cotton; and (3) Finely map the introgressed major or stable fiber quality QTLs, and further evaluate the degree and direction (positive or negative) of the effects of those QTLs on the main lint yield traits, such as lint percentage. In addition, by using introgressed germplasm with superior fiber quality, we aimed to explore molecular breeding strategies for the simultaneous improvement of yield and fiber quality in upland cotton.
Materials and Methods
Plant Materials
Lumianyan22 (LMY22) was an upland Bt-cotton cultivar with high yield (high lint percentage) and medium staple developed by our research group. Luyuan 343 (LY343) was an upland cotton germplasm of high-fiber quality introgressed with chromosomal segments from G. barbadense cv. Ashimouni (Su et al., 2000; Wang et al., 2011). A recombinant inbred lines population was constructed using a cross of LMY22 × LY343 (Wang et al., 2016). We selected two high-fiber quality lines, R472 and R497, which include the target introgression chromosomal region (Supplementary Table S1), from the RIL population; and crossed with the recurrent parent, LMY22, for constructing mapping populations.
Construction of the Secondary Mapping Populations, Phenotypic Valuation and Identification of Putative Introgression Genomic Components
The cross of R472 × LMY22 and R497 × LMY22 were performed at Experimental Station of Shandong Cotton Research Center (ESSCRC), Linqing County, Shandong province, China, in 2011. And the F1 plants were self-pollinated at Hainan Island in winter season of the year. The F2 populations comprising 514 (for Pop. B) and 544 (for Pop. C) individuals, were grown at ESSCRC and in routine cotton trial field in summer season of the year 2012. We named the two populations as Pop. B and Pop. C to be distinguished from the primary introgression populations (defined as “Pop. A” hereinafter) constructed by LY343 × LMY22. All F2 individuals were self-pollinated and the seed cotton samples were harvested. The seeds of each F2 individuals were homogenized and divided randomly into two similar populations of F2:3 family lines. The two secondary populations were planted separately at different environments in summer season of the year 2013. One was planted at Changsha (28°12′ N, 112°59′E), Hunan province of China, and the other was planted at ESSCRC (36°81′N, 115.71°13′E), Shandong province of China. The fiber quality and agronomical traits evaluations were performed on every individual (line) of all the three environments, and three datasets (E1: 2012, ESSCRC, F2; E2: 2013, ESSCRC, F2:3; E3: 2013, Hunan, F2:3) for the followed QTL mapping.
All fiber quality parameters of fiber length (FL, mm), fiber strength (FS, cN/Tex), fiber length uniformity (FU), fiber elongation (FE) and Micronaire value (FM, an integrated fiber quality parameter of fineness and maturity) were tested by the Supervision, Inspection and Test Center of Cotton Quality (SITCCQ), Ministry of Agriculture of China (MAC) using a High-Volume Precision Instrument (HVI; Zellweger-Uster, Knoxville, TN, United States). All other agronomical traits including lint percentage were evaluated by conventional methods of cotton breeding.
The introgression marker loci and putative introgression genomic components derived from G. barbadense L cv. Ashimouni were identified followed our previous report (Wang et al., 2011, 2013).
Molecular Markers and Assays
Genomic DNA of all plant materials were extracted followed by the method reported by Paterson et al. (1993). The SSR (Simple Sequence Repeat) markers which have been exploited for QTL mapping in tetraploid cotton (Guo et al., 2007; Shen et al., 2007; Wang et al., 2008, 2013; Lacape et al., 2009; Zhang et al., 2009, 2011, 2012; Sun et al., 2012; Yu et al., 2014; Shi et al., 2015) were used in this study. All sequences of the SSR primers from CMD1 and CottonDB1 were synthesized by Sangon Biotech (Shanghai, China). The SSR markers which were polymorphic between R472 and LMY22, R497 and LMY22 were screened out as affirmative markers that used for mapping. The affirmative polymorphic primers were used to genotype F2 individual plants. The SSR assay were performed according to the procedure described by Zhang et al. (2002).
Genetic Map Construction and QTL Analysis
All phenotypic data (including five parameters of fiber quality, and the main lint yield trait, i.e., lint percentage) were analyzed using SPSS17.0. The linkage map was constructed using Joinmap4.0 with a LODs core of 6.5. QTLs were detected by WinQTL Cartographer 2.5 using composite interval mapping (Wang et al., 2007). A stringent LOD threshold value were estimated by 1000 permutation test for all traits, and was used to declare the significant QTLs. The graphic representation of the linkage group and QTL marked were created by Map Chart 2.2 (Voorrips, 2002). QTL nomenclature was employed as before (Wang et al., 2011, 2013). QTLs detected in more than two data sets were considered as common or stable QTLs.
Candidate Gene Analysis of the Stable QTL for Fiber Quality
The data of whole-genome sequences of G. hirsutum acc.TM-1 were downloaded from public databases2 (Zhang et al., 2015). The sequences between the adjacent markers for the stable QTL to genomic sequences were performed using SegHunter Software (Ye et al., 2010). Genes within the QTL mapped interval were defined as candidates responsible for the fiber quality trait. The candidate genes were identified by extracting the sequences for between the adjacent markers according to the position of the markers on corresponding chromosome, and comparing to the CDS sequences of the G. hirsutum acc.TM-1 to analyze the identity, alignment length and position. A bioinformatics analysis of annotated genes in the target regions was then performed.
Quantitative RT-PCR
Total RNA was isolated from the developing cotton fibers [10-, 15-, 20-, and 25-day post anthesis (DPA)] of LMY22 and LY343 using a quick RNA Isolation Kit (Huayueyang, Beijing) according to the manufacturer’s protocol. The RNA qualities were detected by agarose gel electrophoresis, and RNA concentrations were measured using NanoPhotometer (ND2000). One μg RNA of each samples were prepared to perform cDNAs synthesis using PrimeScript RT Reagent Kit with gDNA Eraser (TaKaRa, Japan). In order to verify the temporal expression patterns of target genes in developing fibers, qRT-PCR was performed on the grounds of the protocol of SYBG® Premix Ex TaqTM II (TaKaRa, Japan) on the LightCycler® 96 Real-Time PCR System (Roche, Switzerland). The constitutive expression gene β-Actin was used as internal reference to normalize the relative expression levels (forward primer: 5′-ATCCTCCGTCTTGACCTTG-3′, and reverse primer: 5′-TGTCCGTCAGGCAACTCAT-3′). Three biological and two technical replicates were performed to make the results credible. The relative gene expression level was calculated by the 2-ΔΔCt method.
Results
Phenotypic Variation of Fiber Quality
To deeply dissect the genetic contribution of the introgressed G. barbadense alleles to improved fiber quality in G. hirsutum, we backcrossed two higher fiber quality lines to their recurrent parent to develop two secondary mapping populations. Compared with the female parent LMY22, the two parents from the RIL population, R472 and R497, presented significantly better fiber length, strength, and micronaire, but for fiber uniformity and fiber elongation, no significant difference was detected (Supplementary Table S2).
The phenotypic variation for fiber quality in the two populations is shown in Supplementary Table S3. The average fiber length, strength and micronaire in Pop. C were superior to those in Pop. B. The skewness and kurtosis values of the mapping populations showed that the fiber quality traits in three environments were normally distributed in Pop. B. In addition, the fiber quality traits in the three environments (except that the fiber uniformity in the environment of F2 LQ) were normally distributed in Pop. C.
Construction of Genetic Map
In this study, 101 and 112 SSR markers (screened from 366 SSR primer pairs that were polymorphic in our primary mapping population; Wang et al., 2013) revealed polymorphism between the two selected RIL parents (R472 and R497) and the recurrent parent (LMY22). Herein, for Pop. B (LMY22 × R472), ninety-two marker loci were mapped to eighteen linkage groups, which were assigned to fifteen chromosomes, and nine loci were not linked to any group. The genetic map spanned 516.2 cM and accounted for 11.6% of the total allotetraploid cotton genome; the average distance between adjacent markers was approximately 5.61 cM. The minimum distance was 0.1 cM (Supplementary Table S4 and Figure 1). For Pop. C (LMY22 × R497), 105 marker loci were mapped to nineteen linkage groups on fifteen chromosomes, and seven loci were not linked to any group. The genetic map spanned 822.4 cM and accounted for 18.48% of the total cotton genome; the average distance was approximately 7.83 cM between adjacent markers. The minimum distance was 0.1 cM (Supplementary Table S4 and Figure 1).
FIGURE 1. Quantitative trait loci (QTLs) for fiber quality detected in three environments by secondary mapping populations. (A): QTLs for fiber quality detected in three environments by Pop. B. (B): QTLs for fiber quality detected in three environments by Pop. C. The linkage group bars filled with red color indicated the putative introgression chromosomal regions or alleles. Deviated loci were shown with italics in linkage groups. The rhombus, square, pentagram, circular, and triangle in the right of the corresponding linkage groups represented for fiber length (FL), fiber strength (FS), micronaire (FM), fiber elongation (FE) and fiber uniformity (FU). The different color of purple, red and green represented the QTL detected in E1, E2, and E3, respectively.
Putative Introgressed Genomic Components
A total of 21 out of 92 mapped SSR loci were identified as putative introgressed G. barbadense loci in Pop. B. The introgressed loci were distributed on 11 introgressed segments assigned to seven chromosomes, i.e., Chr. 3, Chr. 7, Chr. 13, Chr. 15, Chr. 16, Chr. 18, and Chr. 19. The introgressions accounted for approximately 21.50, 1.48, 46.04, 19.83, 32.59, 41.965, and 7.57% of each chromosome, respectively (Supplementary Table S5). In addition, twenty of one hundred and five mapped SSR loci were identified as putative introgressed loci in Pop. C. These introgressed loci were distributed on 12 introgressed segments assigned to eight chromosomes, i.e., Chr. 3, Chr. 5, Chr. 7, Chr. 18, Chr. 19, Chr. 21, Chr. 23, and Chr. 25. The introgressions accounted for approximately 13.47, 27.91, 7.72, 15.62, 1.21, 12.69, 52.32, and 88.64% of each chromosome, respectively (Supplementary Table S6).
QTLs for Fiber Quality Traits and Lint Percentage
A total of 43 QTLs for fiber quality and lint percentage traits were detected on 13 chromosomal regions; these QTLs explained 2.69–18.8% of the total phenotypic variation. Sixteen and twenty-three QTLs were detected in Pop. B and in Pop. C, respectively. Four of the QTLs were common to both Pop. B and C, and 12 QTLs were defined as stable QTLs (Table 1 and Figure 1). The favorable alleles of 25 QTLs (60.98%) were derived from the secondary introgression parents (R472 or R497), whose superior fiber quality properties were derived from our primary introgression parent, LY343.
TABLE 1. The QTLs for fiber-related traits detected in F2 and F2:3 of the secondary introgression mapping population.
Fiber Length
Eight QTLs explaining 3.43–10.77% of the phenotypic variation were detected. Three and four QTLs were detected in Pop. B and Pop. C, respectively. One QTL was found in both populations. Four QTLs (qFL-C3-1, qFL-C7-1, qFL-C16-2, and qFL-C19-1) were identified as being stable, of which qFL-C7-1 was detected in three data sets. This QTL explained 4.00–6.98% of the phenotypic variation and was tightly linked to marker DC40182 (Table 1 and Figure 1), while the other QTLs, qFL-C3-1, qFL-C16-2, and qFL-C19-1, were detected in two data sets.
Fiber Uniformity
Six QTLs for fiber uniformity were detected and explaining 2.91–5.68% of the phenotypic variation. All the QTLs were detected in only one data set. The favorable alleles of the QTLs were derived from superior-fiber quality parents except that of qFU-C22-1 (Table 1 and Figure 1).
Fiber Micronaire
Eleven QTLs were detected, and four (qFM-C7-1, qFM-C7-2, qFM-C21-1, and qFM-C25-1) of QTLs were identified as stable QTLs. The qFM-C7-1 was detected in four data sets, explaining the phenotypic variation from 10.61 to 18.80% and tightly linked to marker DC40182. Regarding the other three stable QTLs, qFM-C7-2 was detected in two data sets and explained 9.64–12.27% of the phenotypic variation, qFM-C21-1 was detected in two data sets and explained 2.91–3.13% of the phenotypic variation, and qFM-C25-1 was detected in three data sets and explained 3.01–8.52% of the phenotypic variation. With the exception of those of the QTL qFM-C7-2, the favorable alleles were derived from the high-quality secondary introgression parents. The other QTLs were detected in only one data set (Table 1 and Figure 1).
Fiber Strength
A total of nine QTLs for fiber strength were detected in the two secondary introgression mapping populations. Five and three QTLs were detected in Pop. B and Pop. C, respectively. One QTL was common to both populations. Three of the QTLs (qFS-C3-1, qFS-C7-1, and qFS-C15-1) were identified as stable QTLs, and their favorable alleles were derived from the secondary introgression parents. The stable QTL qFS-C7-1 was detected in four data sets and explained 3.00–18.47% of the phenotypic variation; this QTL was tightly linked to marker DC40182. The other stable QTLs, qFS-C3-1 and qFS-C15-1, were detected in two data sets and explained 3.47–5.21% and 3.74–4.12% of the phenotypic variation, respectively. The remaining six QTLs were detected in only one data set, and with the exception of those of the QTL qFS-C16-1, the favorable alleles were derived from the superior-fiber quality parents (Table 1 and Figure 1).
Fiber Elongation
A total of seven QTLs for fiber elongation were detected in the two mapping populations. These QTLs explained 3.27–9.06% of the phenotypic variation. One QTL, qFE-C7-1, was detected in two data sets and defined as a stable QTL; this QTL was also tightly linked to marker DC40182. The other six QTLs were detected in only one data set (Table 1 and Figure 1).
Lint Percentage
Two QTLs for lint percentage were detected in Pop. C and explained 5.57–8.87% of the phenotypic variation. Both QTLs were detected in only one data set, and the favorable alleles were derived from the high lint-percentage parent LMY22 (Table 1 and Figure 1).
Common QTLs Across Multiple Populations
The QTLs for a defined trait detected in two different populations for which the confidence intervals overlap are considered “common” QTLs (Kenis et al., 2008). In the present study, we comprehensively compared the QTLs that were detected in the two secondary introgression populations (Pop. B and Pop. C) with the previously detected QTLs in the primary population (Pop. A) (Wang et al., 2011) and found that five common QTLs were shared between different populations, including two for fiber length, one for fiber strength, and two for fiber elongation (Table 2 and Figure 2). The qFE-C3-1 QTL for fiber elongation in the introgressed region of Chr. 3 was common to both Pop. A and Pop. C, and the qFL-C16-1 QTL for fiber length on Chr. 16 was common to both Pop. A and Pop. B. In the introgressed region of Chr. 7, three common QTLs, including qFL-C7-1, qFS-C7-1, and qFE-C7-1, were common to Pop. A, Pop. B and Pop. C. These three QTLs were also detected in the RIL population developed from the same parents (Wang et al., 2016). These common QTLs are considered stable QTLs and can be applied via MAS as favorable introgressed alleles to improve fiber quality in cotton breeding.
TABLE 2. Common QTLs detected in three different populations.
FIGURE 2. The common QTLs detected among different populations. (A–C) was representing the common QTLs on Chr. 3, Chr. 7, and Chr. 16, respectively. The rhombus, square, pentagram, circular, and triangle in the side of the corresponding linkage groups represented for fiber length, fiber strength, micronaire, fiber elongation, and fiber uniformity. The different color of purple, red, and green represented that the QTLs were detected in E1, E2, and E3, respectively. The black box represented for the region of the QTLs in Pop. A.
QTL Clusters Associated With Fiber-Related Traits
A QTL cluster is defined as a chromosomal region that contains multiple QTLs for more than three associated traits. In this study, we identified five QTL clusters which comprised of 27 QTLs and were distributed on Chr. 3, Chr. 7, Chr. 16, and Chr. 25 (Table 3). Chromosome 3 contained two QTL clusters, both of which harbored three QTLs and were detected in Pop. B and Pop. C. The stable QTL qFE-C3-2 detected in Pop. C and Pop. A resided in the region of c3-cluster-2 (Table 3 and Figure 2). An important QTL cluster, the c7-cluster, was located on Chr. 7 and was detected in Pop. B and Pop. C. This QTL cluster was also detected in almost all data sets of Pop. A and of the RIL population developed from the same parents. Two major QTLs, qFS-C7-1 and qFM-C7-1, which explained more than 10% of the phenotypic variation for fiber strength and micronaire, were co-located in this cluster region, suggesting that key genes controlling fiber development may reside in this region. This QTL cluster region was flanked by markers DPL0757 and DC40182, and the approximate physical distance was 185 kb (Table 3 and Figure 2). The QTL cluster c16-cluster, contained five fiber quality QTLs, and a stable QTL, qFL-C16-1, resided within this cluster region (Table 3 and Figure 2). The c25-cluster detected in Pop. C contained four QTLs, and the approximate physical distance was 209.4 kb (Table 3).
TABLE 3. Distribution of QTL clusters in Pop. B and Pop. C.
Prediction and Comparative Analysis of Candidate Genes in the c7-Cluster Region
The 0.6 cM region between DPL0757 and DC40182 corresponds to a 185-kb genomic region on chromosome A07 of the G. hirsutum acc. TM-1 genome. According to the reference genome annotation of G. hirsutum (Zhang et al., 2015), the 185-kb genomic region contains five predicted genes: those coding for the RING/U-box superfamily protein (RHA2B, Gh_A07G1752), CONSTANS-like 9 (COL9, Gh_A07G1753), and a plant protein of unknown function (DUF868) (Gh_A07G1755), as well as other two genes with no known function (Gh_A07G1754 and Gh_A07G1756) (Supplementary Table S7). Coding sequence comparisons between the chromosomes of G. hirsutum and G. barbadense revealed some asymmetry, and the four genes in the G. hirsutum genome have homoeologous counterparts in the G. barbadense reference genome sequence, i.e., Gh_A07G1752 is homoeologous to GOBAR_AA08023, Gh_A07G1753 is homeologous to GOBAR_AA08022, Gh_A07G1755 is homeologous to GOBAR_AA08021, and Gh_A07G1756 is homeologous to GOBAR_AA08020 (Figure 3).
FIGURE 3. Annotated genes within the 185-kb in c7-cluster region between the G. hirsutum and G. barbadense. Gh_A07G1752 (RING/U-box superfamily protein, RHA2B); Gh_A07G1753 (CONSTANS-like 9, COL9); Gh_A07G1755 (Plant protein of unknown function, DUF868); Gh_A07G1754 and Gh_A07G1756 (without any function). GOBAR_AA08020 (uncharacterized membrane protein); GOBAR_AA08021 (conserved hypothetical protein); GOBAR_AA08022 (B-box type zinc finger protein with CCT domain isoform 1); GOBAR_AA08023 (RING/U-box superfamily protein, putative).
The expression levels of the five annotated genes within this QTL region were tested by RT-qPCR using fiber tissue from four developmental stages (10, 15, 20, and 25 DPA) in LMY22 and LY343. Via the sequence alignment of G. hirsutum acc. TM-1 (Zhang et al., 2015), single locus-specific primers of the five candidate genes were designed for the divergent coding region (Supplementary Table S8). The relative expression levels in both parents (LMY22 and LY343) are shown in Figure 4. The expression level of Gh_A07G1752 were gradually increased from 10 DPA to 25 DPA, especially in LY343; compared with those in LMY22, the expression levels of Gh_A07G1752 in LY343 exhibited an approximately six-fold increase at the 20 DPA stage. The gene Gh_A07G1753 was significantly differentially expressed between LMY22 and LY343, especially at 15 and 20 DPA. The expression levels of Gh_A07G1754 did not significantly differ between LMY22 and LY343 during the four developmental stages. Moreover, we did not detect any expression for Gh_A07G1755 or Gh_A07G1756 during the four developmental stages.
FIGURE 4. RT-qPCR expression assay of genes in c7-cluster region during fiber development. (A) RING/U-box superfamily protein (RHA2B, Gh_A07G1752), (B) CONSTANS-like 9 (COL9, Gh_A07G1753), and (C) Not any function (Gh_A07G1754). The white bars are L22, and the black bars are LY343. Error bars indicate the standard deviation of three biological replicates. ∗∗Denote significantly differential expression at p < 0.01.
Discussion
Genetic Effects of Stable Introgressed QTLs and QTL Clusters Detected in Secondary Mapping Populations
To identify introgressed exogenous chromosomal regions, we identified 158 out of 322 (158/322) SSR marker loci in a genetic map constructed using a primary introgression mapping population developed from a cross of LMY22 × LY343; these markers represented putative introgressed G. barbadense alleles (Wang et al., 2013). The 322 polymorphic SSR markers were screened from a total of 18,467 SSRs previously exploited in interspecific mapping (Nguyen et al., 2004; Guo et al., 2007; Lacape et al., 2009; Zhang et al., 2009, 2011). The total introgressions covered 646.91 cM, which accounted for approximately 12.86% of the tetraploid cotton genome, and were located on 41 introgressed chromosomal segments assigned to 20 chromosomes (Wang et al., 2013). QTL mapping demonstrated that most of the favorable alleles for fiber quality QTLs were located in the introgressed chromosomal regions derived from the primary introgression parent, i.e., LY343. Stable QTLs and/or QTL clusters were identified in only a few of the introgressed chromosomal regions, including the introgressed regions on Chr. 7 and Chr. 25. These introgressed chromosomal regions might contain important genes for superior fiber quality (Wang et al., 2013, 2016).
In the present study, by backcrossing secondary introgression upland cotton lines that present both desirable higher fiber quality and relatively better lint yield performance with their high lint-yielding recurrent parent, we developed a secondary introgression mapping population to anchor the target introgressed chromosomal regions. The mapping parents were individual lines from the RIL population developed from the primary introgression parent (Wang et al., 2016). The secondary introgression parent for developing Pop. B, i.e., R472, contains the introgressed regions on Chr. 7, Chr. 13, Chr. 15, Chr. 16, Chr. 18, Chr. 19, and Chr. 20, while the parent for developing Pop. C, i.e., R497, contains the introgressed regions on Chr. 3, Chr. 5, Chr. 7, Chr. 8, Chr. 18, Chr. 19, Chr. 21, Chr. 22, Chr. 23, and Chr. 25 (Supplementary Table S1). These parents contain the key introgressed chromosomal regions that have favorable alleles for fiber quality QTLs or QTL clusters previously detected by us (Wang et al., 2011, 2013, 2016).
Using these two secondary introgression mapping populations (Pop. B and Pop. C), we detected 41 QTLs for fiber quality, and these QTLs were distributed on all the above mentioned chromosomal regions (Table 1). As expected, most of the QTLs for fiber quality were mapped to the introgressed chromosomal regions, and most of the favorable alleles of these QTLs were contributed by the introgressed genomic components (Figure 1) We detected introgressed fiber quality QTLs or QTL clusters on Chr. 3, Chr. 7, Chr. 15, Chr. 16, Chr. 18, and Chr. 19 in Pop. A, and we detected introgressed QTLs or QTL clusters on Chr. 3, Chr. 7, Chr. 19, Chr. 21, and Chr. 25 in Pop. B. In contrast with the mapping results for Pop. A previously reported by us (Wang et al., 2011, 2013), the present results showed five groupings of common QTLs in the introgressed region of Chr. 7, i.e., qFL-C7-1-Pop. A, qFL-C7-1-Pop. B and qFL-C7-1-Pop. C; qFS-C7-1-Pop. A, qFS-C7-1-Pop. B and qFS-C7-1-Pop. C; qFS-C7-2-Pop. A and qFS-C7-3-Pop. C; qFE-C7-1-Pop. A, qFE-C7-1-Pop. B and qFE-C7-1-Pop. C; and qFE-C7-1-Pop. A and qFE-C7-2-Pop. C, among which the three QTLs for fiber length, strength, and elongation were detected in all the mapping populations (Figure 2). As we reported recently, these three QTLs were also detected in RIL populations in multiple environments (Wang et al., 2016) and were reasonably stable introgressed QTL clusters for fiber quality. In addition, mapping via linkage analysis (Zhang et al., 2009; Sun et al., 2012; Xu et al., 2012; Cao et al., 2015; Shi et al., 2015; Jamshed et al., 2016; Fang X. et al., 2017) and genome-wide association studies (GWAS) (Islam et al., 2016; Su et al., 2016; Fang L. et al., 2017; Sun et al., 2017; Tan et al., 2018) have both been used to detect the major QTLs for fiber quality in this chromosomal region.
We detected one common QTL on Chr. 3, namely, qFE-C3-1-Pop. A and qFE-C3-2-Pop. C, for fiber elongation (Figure 2). We also detected one pair of common QTLs on Chr. 16, namely, qFL-C16-1-Pop. A and qFL-C16-1-Pop. B, for fiber length (Figure 2). Together with the QTL pairs detected on Chr. 7, without exception, those common QTLs were located in the introgressed chromosomal regions. In addition to the major or stable QTLs and QTL clusters for fiber quality located in the introgressed region of Chr. 25 (Figure 1), these introgression chromosomal regions might contain favorable alleles (or genes) for the desirable superior fiber quality. Detailed and systematic studies of these regions (or alleles) would be helpful for understanding the genetic and/or molecular mechanisms controlling the formation of superior fiber quality in the introgression germplasm.
Candidate Genes for Superior Fiber Quality in the Putative Introgressed Chromosomal Region of Chr. 7
The fiber quality of LY343 was much better than those of LMY22 (Wang et al., 2011). We selected two high fiber-quality secondary introgression lines, i.e., R472 and R497, which contained overlapping introgressed fragments on some chromosomes (Chr. 7, Chr. 18, and Chr. 19), from the RIL population for backcrossing with LMY22 to construct secondary introgression mapping populations designated Pop. B and Pop. C. Using the secondary mapping populations, we mapped the major and stable fiber quality QTLs in the putative introgressed region of Chr. 7 within a 0.5 cM interval between markers DPL0757 and DC40182. According to the reference genome of G. hirsutum L. acc. TM-1 (Zhang et al., 2015), we found that five genes are present in the stable QTL region between DPL0757 and DC40182, namely, Gh_A07G1752, Gh_A07G1753, Gh_A07G1754, Gh_A07G1755, and Gh_A07G1756 (Supplementary Table S7). Gh_A07G1752, a small but important gene, encodes a RING/U-box superfamily protein that plays regulatory roles both in the development of a variety of organisms and in the dramatic increase in cellular abscisic acid (ABA) levels and resulting drought tolerance (Ko et al., 2006; Xia et al., 2013). In the present study, the expression levels of Gh_A07G1752 gradually increased in the four fiber developmental stages (Figure 4). Gh_A07G1753 encodes the zinc finger protein CONSTANS-LIKE 9 (COL9), a member of the CONSTANS-LIKE gene family, and by repressing the expression of CO, COL9 is involved in the regulation of flowering time, concomitantly reducing the expression of FT and delaying floral transition (Cheng and Wang, 2005). Furthermore, based on the Gene Ontology results, Gh_A07G1753 was recently predicted to encode a transcription factor involved in cotton fiber development (Islam et al., 2016). Gh_A07G1755, which encodes a plant protein of unknown function (DUF868), may be a new a gene that needs to be analyzed in the future. Gh_A07G1754 and Gh_A07G1756 have no known function according to the annotation for G. hirsutum L. acc. TM-1 (Zhang et al., 2015). Cotton fiber development is a complex process, and fiber development include four overlapping stages: fiber initiation, cell elongation, secondary wall deposition, and maturation (Lee et al., 2007). Many transcription factors which were involved in metabolic pathways such as phytohormone signaling, energy metabolism of cell, fatty acid metabolism, secondary metabolism and other signaling pathways are related in this process (Hande et al., 2017). So more detailed research is needed to identify the introgressed gene(s) that is (are) involved in cotton fiber development and contribute to superior fibre quality for this study.
QTL Clusters and QTL Hotspots for Fiber Quality and Their Relationships With Introgressed Genomic Components
The existence of QTL clusters has been reported in cotton (Qin et al., 2008; Jiang et al., 2009; Said et al., 2013, 2015; Li et al., 2016; Zhang et al., 2016; Si et al., 2017) and other organisms (Xie et al., 2008; Mallard et al., 2013). Analyses of the reported 1,223 QTLs for yield and fiber quality, disease resistance and other traits across the Gossypium genome revealedthe presence of QTL clusters and/or specific trait QTL hotspots on almost all of the tetraploid chromosomes except Chr.9, Chr.13, Chr.20, and Chr.22 (Said et al., 2013). This finding indicates that genes related to specific traits are more concentrated within certain regions of the genome than others.
In our study, several QTLs that influence fiber quality were detected in the same chromosomal regions, indicating that QTL clusters and hotspots were also present. We detected five clusters: c3-cluster-1, c3-cluster-2, c7-cluster, c16-cluster, and c25-cluster. As mentioned above, the c7-cluster contained five annotated or predicted genes (Supplementary Table S7). There are many reports about QTL clusters of fiber quality on Chr. 7 (Zhang et al., 2009; Sun et al., 2012; Xu et al., 2012; Cao et al., 2015; Shi et al., 2015). However, no reports mentioned QTL clusters within the DPL0757-DC40182 interval. In the DPL0757-DC40182 interval located on Chr. 7 in Pop. B and Pop. C, qFL-C7-1, qFS-C7-1, qFM-C7-1 and qFE-C7-1 influenced fiber quality, including fiber length, strength, micronaire and elongation. This finding indicated that the DPL0757-DC40182 interval located on Chr. 7 is a QTL cluster and a hotspot interval for fiber quality. In addition to our previously reported mapping results in Pop. A (Wang et al., 2011, 2013), the QTL qFS-C7-1 was detected in seven data sets (F2:3 in LQ and HuN in both Pop. B and Pop. C; F2 in LQ, F2:3 in LQ and HaiN in Pop. A), and except in one data set, was tightly linked to marker DC40182 (Pop. A and the F2:3 generation in LQ) (Wang et al., 2013). This QTL explained 3–26.21% of the phenotypic variation. As such, the QTL qFS-C7-1 is a stable major QTL. It provides an underlying basis for cloning QTLs for both understanding the genetic mechanism of the superior fiber formation and deep mining the introgressed alleles (genes) associated with fiber development to further design and precisely breed (via MAS) for improved fiber quality in upland cotton.
Evaluation of the Genetic Effects of the Introgressed Fiber Quality QTLs and/or QTL Cluster on Lint Yield Traits and Possible MAS-Based Strategies for Simultaneously Improving Fiber Quality and Lint Yield in Cotton Breeding
Cotton lint yield is mainly determined by the number of bolls in per area (NB), boll weight (BW) and lint percentage (or gin turnout), among these lint percentage is an important factor for determining cotton fiber production. Several studies have shown that fiber quality has a negatively genetically correlated with lint yield (Scholl and Miller, 1976; Zhou et al., 2003; Shen et al., 2007; Yu et al., 2013). Therefore, this phenomenon is a major obstacle with respect to the simultaneous improvement of lint yield (mainly lint percentage) and fiber quality in upland cotton breeding programs (Wang et al., 2013, 2016).
A total of 41 QTLs for fiber quality and two QTLs for lint percentage were detected in this study; most of those QTLs for fiber quality were related to the introgressed chromosomal regions or introgressed alleles derived from G. barbadense, and the favorable alleles of 25 QTLs are introgressions (Table 1 and Figure 1). Interestingly, the favorable alleles of fourteen of seventeen (14/17) fiber length- and strength-related QTLs were introgressions. Only the favorable alleles of three QTLs on Chr. 16 originate from LMY22, i.e., they were not introgressions. This result further indicated that the introgressed alleles derived from G. barbadense play a key role in the formation of superior fiber quality, especially fiber length and strength.
Two QTLs for lint percentage (qLP-C13-1 and qLP-C25-1) were identified in one of the secondary mapping populations (Pop. C). The favorable alleles of both QTLs are derived from LMY22. One of the QTLs, qLP-C25-1 was clustered with qFL-C25-1, and they are both tightly linked to marker HAU2049 (Table 1). The favorable alleles of qLP-C25-1 and qFL-C25-1 originate from LMY22 and the secondary introgression parent, respectively. This finding indicated that these introgressed alleles are favorable for improving fiber length but it is not desirable for improving lint percentage. Therefore, regardless of whether the negative effect is caused by linkage drag or pleiotropism, this QTL and its tightly linked marker are not suitable for the simultaneous improvement of fiber quality and lint yield in cotton breeding.
We previously detected five QTLs for lint percentage (Wang et al., 2013); these QTLs were consistent with the lint percentage QTLs reported by Shi et al. (2015). Recently, using a RIL population in multiple environments, we also detected nine QTLs for lint percentage (Wang et al., 2016). One of the QTLs, namely, qLP16.1, located on Chr. 16 (D07) was detected in four environments, and the favorable allele originated from the introgression parent; unfortunately, we detected negative correlations between lint percentage and fiber quality traits in this introgressed chromosomal region. Chr. 16 is homeologous to Chr. 7 (A07), which contain two QTL clusters for fiber quality (Figure 1 and Table 3). This finding suggests that polyploid lineages may determine the variations in subgenome contributions to QTLs (Fonceka et al., 2012). However, we detect no negative effects on lint percentage from the fiber quality-related QTLs and/or QTL clusters located on the introgressed chromosomal region of Chr. 7. Therefore, we believe that these fiber quality QTLs and/or QTL clusters are favorable for MAS to improve fiber quality in cotton breeding.
The simultaneous improvement of fiber quality and lint yield by conventional breeding methods is slow and difficult (Zhang et al., 2012). Molecular markers constitute a foundation for constructing a high-density linkage map to exploit effective markers (new SSRs, SNPs, and InDels) (Wang et al., 2015) and lay a foundation both for dissecting additional QTLs for lint percentage, and fiber quality and for stacking favorable alleles for the fiber quality and lint percentage mapping in non-homologous regions, enhancing the possibility of using MAS to improve both fiber quality and lint yield (Wang et al., 2013). We detected seven stable QTLs for fiber quality in this study (Table 3), and most of the QTLs have no negative effects on lint percentage. Notably, the fiber quality QTL cluster in the introgressed chromosomal region of Chr. 7 had significant genetic effects on fiber strength, length, micronaire. In addition, the genetic distance was 0.6 cM between the adjacent molecular markers, which is inaccordance with requirements.
Importantly, the QTLs for lint percentage, whose favorable alleles originate from non-introgression upland cotton parent, were mainly detected on Chr. 5, Chr. 15, Chr. 17, Chr. 18, and Chr. 22 using the genetic populations involved in introgressions from G. barbadense in previous studies (Wang et al., 2011, 2013, 2016) and other reports (Jenkins et al., 2006; Yu et al., 2013; Shi et al., 2015). Notably, these QTLs for lint percentage were distributed on the chromosomes that are not harboring the major QTLs for superior fiber quality. In addition, we detected no mutual negative effects on fiber quality and lint percentage in the majority of the fiber quality and lint percentage QTLs and/or QTL clusters, which is extremely important. Accordingly, we concluded that favorable elite alleles or QTLs can be pyramided via MAS for the simultaneous improvement on lint yield and quality in upland cotton breeding.
Author Contributions
JZ and FW designed the experiments. GL, YC, JhZ, CZ, FW, and JZ performed field trials, phenotypic evaluation, and data collection. YC and JxZ contributed to the preparation of cotton RNA samples. HM performed the analysis of QTL mapping data, with the help of FW for the analysis of raw data. ZS performed the qRT-PCR. YC and FW drafted the manuscript. JZ revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was financially supported in part by grants from the National Key R&D Project in China (2016YFD0100203), the Youth Science Foundation in Shandong Academy of Agricultural Sciences (2016YQN06), the National Science Foundation in China (31671742 and 31601345), National Project of Modern Agricultural Industry Technology System in China (CARS-15-05), Taishan Scholars Program of Shandong Province (No. ts201511070) and the Innovation Project in Shandong Academy of Agricultural Sciences (CXGC2016A01).
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.
The reviewer BZ declared a shared affiliation, though no other collaboration, with one of the authors JhZ to the handling Editor.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01023/full#supplementary-material
Footnotes
References
Ahmad, S., Mahmood, K., Hanif, M., Nazeer, W., Malik, W., Qayyum, A., et al. (2011). Introgression of cotton leaf curl virus-resistant genes from Asiatic cotton (Gossypium arboreum) into upland cotton (G. hirsutum). Genet. Mol. Res. 10, 2404–2414. doi: 10.4238/2011.October.7.2
Cao, Z., Wang, P., Zhu, X., Chen, H., and Zhang, T. (2014). SSR marker-assisted improvement of fiber qualities in Gossypium hirsutum using G. barbadense introgression lines. Theor. Appl. Genet. 127, 587–594. doi: 10.1007/s00122-013-2241-3
Cao, Z., Zhu, X., Chen, H., and Zhang, T. (2015). Fine mapping of clustered quantitative trait loci for fiber quality on chromosome 7 using a Gossypium barbadense introgressed line. Mol. Breed. 35:215. doi: 10.1007/s11032-015-0393-3
Chandra, M., and Sreenivasan, S. (2011). Studies on improved Gossypium arboreum cotton: part I— Fibre quality parameters. IJFTR 36, 24–34.
Cheng, X. F., and Wang, Z. Y. (2005). Overexpression of COL9, a CONSTANS-LIKE gene, delays flowering by reducing expression of CO and FT in Arabidopsis thaliana. Plant J. 43, 758–768. doi: 10.1111/j.1365-313X.2005.02491.x
Culp, T. W., and Harrell, D. C. (1973). Breeding methods for improving yield and fiber quality of upland cotton (Gossypium hirsutum L.)1. Crop Sci. 13, 686–689. doi: 10.2135/cropsci1973.0011183X001300060030x
Culp, T. W., Harrell, D. C., and Kerr, T. (1979). Some genetic implications in the transfer of high fiber strength genes to upland cotton 1. Crop Sci. 19, 481–484. doi: 10.2135/cropsci1979.0011183X001900040013x
Duran, C., Edwards, D., and Batley, J. (2009). Genetic maps and the use of synteny. Methods Mol. Biol. 513, 41–55. doi: 10.1007/978-1-59745-427-8_3
Fang, L., Wang, Q., Hu, Y., Jia, Y., Chen, J., Liu, B., et al. (2017). Genomic analyses in cotton identify signatures of selection and loci associated with fiber quality and yield traits. Nat. Genet. 49, 1089–1098. doi: 10.1038/ng.3887
Fang, X., Liu, X., Wang, X., Wang, W., Liu, D., Zhang, J., et al. (2017). Fine-mapping qFS07.1 controlling fiber strength in upland cotton (Gossypium hirsutum L.). Theor. Appl. Genet. 130, 795–806. doi: 10.1007/s00122-017-2852-1
Fonceka, D., Tossim, H. A., Rivallan, R., Vignes, H., Faye, I., Ndoye, O., et al. (2012). Fostered and left behind alleles in peanut: interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding. BMC Plant Biol. 12:26. doi: 10.1186/1471-2229-12-26
Grover, C. E., Zhu, X., Grupp, K. K., Jareczek, J. J., Gallagher, J. P., Szadkowski, E., et al. (2015). Molecular confirmation of species status for the allopolyploid cotton species, Gossypium ekmanianum Wittmack. Genet. Resour. Crop Evol. 62, 103–114. doi: 10.1007/s10722-014-0138-x
Guo, W., Cai, C., Wang, C., Han, Z., Song, X., Wang, K., et al. (2007). A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in Gossypium. Genetics 176, 527–541. doi: 10.1534/genetics.107.070375
Hande, A. S., Katageri, I. S., Jadhav, M. P., Adiger, S., Gamanagatti, S., Padmalatha, K. V., et al. (2017). Transcript profiling of genes expressed during fibre development in diploid cotton (Gossypium arboreum L.). BMC Genomics 18:675. doi: 10.1186/s12864-017-4066-y
Iqbal, M. A., Abbas, A., Zafar, Y., and Rahman, M.-U. (2015). Characterization of indigenous Gossypium Arboreum L. Genotypes for various fiber quality traits. Pak. J. Bot. 47, 2347–2354.
Islam, M. S., Thyssen, G. N., Jenkins, J. N., Zeng, L., Delhom, C. D., Mccarty, J. C., et al. (2016). A MAGIC population-based genome-wide association study reveals functional association of GhRBB1_A07 gene with superior fiber quality in cotton. BMC Genomics 17:903. doi: 10.1186/s12864-016-3249-2
Jamshed, M., Jia, F., Gong, J., Palanga, K. K., Shi, Y., Li, J., et al. (2016). Identification of stable quantitative trait loci (QTLs) for fiber quality traits across multiple environments in Gossypium hirsutum recombinant inbred line population. BMC Genomics 17:197. doi: 10.1186/s12864-016-2560-2
Jenkins, J. N., Wu, J., Mccarty, J. C., and Saha, S. (2006). Genetic effects of thirteen Gossypium barbadense L. chromosome substitution lines in topcrosses with Upland cotton cultivars: I. Yield and yield components. Crop Sci. 46, 1169–1178. doi: 10.2135/cropsci2005.08-0269
Jiang, C., Wright, R. J., El-Zik, K. M., and Paterson, A. H. (1998). Polyploid formation created unique avenues for response to selection in Gossypium (cotton). Proc. Natl. Acad. Sci. U.S.A. 95, 4419–4424. doi: 10.1073/pnas.95.8.4419
Jiang, F., Zhao, J., Zhou, L., Guo, W., and Zhang, T. (2009). Molecular mapping of Verticillium wilt resistance QTL clustered on chromosomes D7 and D9 in upland cotton. Sci. China C Life Sci. 52, 872–884. doi: 10.1007/s11427-009-0110-8
Kenis, K., Keulemans, J., and Davey, M. W. (2008). Identification and stability of QTLs for fruit quality traits in apple. Tree Genet. Genomes 4, 647–661. doi: 10.1007/s11295-008-0140-6
Ko, J. H., Yang, S. H., and Han, K. H. (2006). Upregulation of an Arabidopsis RING-H2 gene, XERICO, confers drought tolerance through increased abscisic acid biosynthesis. Plant J. 47, 343–355. doi: 10.1111/j.1365-313X.2006.02782.x
Kohel, R. J., Yu, J., Park, Y.-H., and Lazo, G. R. (2001). Molecular mapping and characterization of traits controlling fiber quality in cotton. Euphytica 121, 163–172. doi: 10.1023/A:1012263413418
Lacape, J. M., Jacobs, J., Arioli, T., Derijcker, R., Forestierchiron, N., Llewellyn, D., et al. (2009). A new interspecific, Gossypium hirsutum × G. barbadense, RIL population: towards a unified consensus linkage map of tetraploid cotton. Theor. Appl. Genet. 119, 281–292. doi: 10.1007/s00122-009-1037-y
Lacape, J. M., Nguyen, T. B., Thibivilliers, S., Bojinov, B., Courtois, B., Cantrell, R. G., et al. (2003). A combined RFLP-SSR-AFLP map of tetraploid cotton based on a Gossypium hirsutum × Gossypium barbadense backcross population. Genome 46, 612–626. doi: 10.1139/g03-050
Lee, J. J., Woodward, A. W., and Chen, Z. J. (2007). Gene expression changes and early events in cotton fibre development. Ann. Bot. 100, 1391–1401. doi: 10.1093/aob/mcm232
Li, C., Dong, Y., Zhao, T., Li, L., Li, C., Yu, E., et al. (2016). Genome-wide SNP linkage mapping and QTL analysis for fiber quality and yield traits in the upland cotton recombinant inbred lines population. Front. Plant Sci. 7:1356. doi: 10.3389/fpls.2016.01356
Liang, Q. Z., Cheng, H. U., Hua, H., Zhaohu, L. I., and Hua, J. P. (2013). Construction of a linkage map and QTL mapping for fiber quality traits in upland cotton (Gossypium hirsutum L.). Sci. Bull. 58, 3233–3243. doi: 10.1371/journal.pone.0182918
Ma, L., Zhao, Y., Wang, Y., Shang, L., and Hua, J. (2017). QTLs analysis and validation for fiber quality traits using maternal backcross population in upland cotton. Front. Plant Sci. 8:2168. doi: 10.3389/fpls.2017.02168
Mallard, S., Cantet, M., Massire, A., Bachellez, A., Ewert, S., and Lefebvre, V. (2013). A key QTL cluster is conserved among accessions and exhibits broad-spectrum resistance to Phytophthora capsici: a valuable locus for pepper breeding. Mol. Breed. 32, 349–364. doi: 10.1007/s11032-013-9875-3
May, O. L., Bowman, D. T., and Calhoun, D. S. (1995). Genetic diversity of U.S. Upland cotton cultivars released between 1980 and 1990. Crop Sci. 35, 1570–1574. doi: 10.2135/cropsci1995.0011183X003500060009x
Mehboob-ur-Rahman, S. T., Tabbasam, N., Iqbal, M. A., Ashraf, M., Zafar, Y., and Paterson, A. H. (2012). Cotton genetic resources. A review. Agron. Sustain. Dev. 32, 419–432. doi: 10.1007/s13593-011-0051-z
Meredith, W. R. (2005). Minimum number of genes controlling cotton fiber strength in a backcross population. Crop Sci. 45, 1114–1119. doi: 10.2135/cropsci2003.0425
Mergeai, G. (2003). “Forty years of genetic improvement of cotton through interspecific hybridization at gembloux agricultural university: achievement and prospects,” in Proceedings of the World Cotton Research Conference-3, in Cotton Production for the New Millennium: World Cotton Research Conference–113 Proceedings, ed. A. Swanepoel (Pretoria: Agricultural Research Council), 120–133.
Meyer, V. G. (1973). A study of reciprocal hybrids between upland cotton (Gossypium hirsutum L.) and experimental lines with cytoplasms from seven other species 1. Crop Sci. 13, 439–444. doi: 10.2135/cropsci1973.0011183X001300040015x
Nacoulima, N., Baudoin, J. P., and Mergeai, G. (2012). Introgression of improved fiber fineness trait in G. hirsutum L. from G. longicalyx Hutch. & Lee. Commun. Agric. Appl. Biol. Sci. 77, 207–211.
Nazeer, W., Ahmad, S., Mahmood, K., Tipu, A. L., Mahmood, A., and Zhou, B. (2014). Introgression of genes for cotton leaf curl virus resistance and increased fiber strength from Gossypium stocksii into upland cotton (G. hirsutum). Genet. Mol. Res. 13, 1133–1143. doi: 10.4238/2014.February.21.2
Nguyen, T. B., Giband, M., Brottier, P., Risterucci, A. M., and Lacape, J. M. (2004). Wide coverage of the tetraploid cotton genome using newly developed microsatellite markers. Theor. Appl. Genet. 109, 167–175. doi: 10.1007/s00122-004-1612-1
Paterson, A. H., Brubaker, C. L., and Wendel, J. F. (1993). A rapid method for extraction of cotton (Gossypium spp.) genomic DNA suitable for RFLP or PCR analysis. Plant Mol. Biol. Rep. 11, 122–127. doi: 10.1007/BF02670470
Qin, H., Guo, W., Zhang, Y. M., and Zhang, T. (2008). QTL mapping of yield and fiber traits based on a four-way cross population in Gossypium hirsutum L. Theor. Appl. Genet. 117, 883–894. doi: 10.1007/s00122-008-0828-x
Said, J. I., Lin, Z., Zhang, X., Song, M., and Zhang, J. (2013). A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton. BMC Genomics 14:776. doi: 10.1186/1471-2164-14-776
Said, J. I., Song, M., Wang, H., Lin, Z., Zhang, X., Fang, D. D., et al. (2015). A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum x G. barbadense populations. Mol. Genet. Genomics 290, 1003–1025. doi: 10.1007/s00438-014-0963-9
Scholl, R. L., and Miller, P. A. (1976). Genetic association between yield and fiber strength in upland cotton. Crop Sci. 16, 780–783. doi: 10.2135/cropsci1976.0011183X001600060010x
Shang, L., Liang, Q., Wang, Y., Wang, X., Wang, K., Abduweli, A., et al. (2015). Identification of stable QTLs controlling fiber traits properties in multi-environment using recombinant inbred lines in Upland cotton (Gossypium hirsutum L.). Euphytica 205, 877–888. doi: 10.1007/s10681-015-1434-z
Shen, X., Guo, W., Lu, Q., Zhu, X., Yuan, Y., and Zhang, T. (2007). Genetic mapping of quantitative trait loci for fiber quality and yield trait by RIL approach in Upland cotton. Euphytica 155, 371–380. doi: 10.1007/s10681-006-9338-6
Shi, Y., Li, W., Li, A., Ge, R., Zhang, B., Li, J., et al. (2015). Constructing a high-density linkage map for Gossypium hirsutum x Gossypium barbadense and identifying QTLs for lint percentage. J. Integr. Plant Biol. 57, 450–467. doi: 10.1111/jipb.12288
Si, Z., Chen, H., Zhu, X., Cao, Z., and Zhang, T. (2017). Genetic dissection of lint yield and fiber quality traits of G. hirsutum in G. barbadense background. Mol. Breed. 37:9. doi: 10.1007/s11032-016-0607-3
Smith, C. W., Hague, S., Hequet, E., Thaxton, P. S., and Brown, I. N. (2008). Development of extra-long staple upland cotton. Crop Sci. 48, 1823–1831. doi: 10.2135/cropsci2008.01.0052
Su, J., Fan, S., Li, L., Wei, H., Wang, C., Wang, H., et al. (2016). Detection of favorable QTL alleles and candidate genes for lint percentage by GWAS in Chinese upland cotton. Front. Plant Sci. 7:1576. doi: 10.3389/fpls.2016.01576
Su, X., Gao, G., Shi, X., and Zhu, D. (2000). Breeding of a terrestrial long-staple cotton—Luyuan 343. Acta Agric. Nucleatae Sin. 14, 180–183.
Sun, F. D., Zhang, J. H., Wang, S. F., Gong, W. K., Shi, Y. Z., Liu, A. Y., et al. (2012). QTL mapping for fiber quality traits across multiple generations and environments in upland cotton. Mol. Breed. 30, 569–582. doi: 10.1007/s11032-011-9645-z
Sun, Z., Wang, X., Liu, Z., Gu, Q., Zhang, Y., Li, Z., et al. (2017). Genome-wide association study discovered genetic variation and candidate genes of fibre quality traits in Gossypium hirsutum L. Plant Biotechnol. J. 15, 982–996. doi: 10.1111/pbi.12693
Tan, Z., Zhang, Z., Sun, X., Li, Q., Sun, Y., Yang, P., et al. (2018). Genetic map construction and fiber quality QTL mapping using the CottonSNP80k array in upland cotton. Front. Plant Sci. 9:225. doi: 10.3389/fpls.2018.00225
Tyagi, P., Gore, M. A., Bowman, D. T., Campbell, B. T., Udall, J. A., and Kuraparthy, V. (2014). Genetic diversity and population structure in the US Upland cotton (Gossypium hirsutum L.). Theor. Appl. Genet. 127, 283–295. doi: 10.1007/s00122-013-2217-3
Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77
Wang, F., Gong, Y., Zhang, C., Liu, G., Wang, L., Xu, Z., et al. (2011). Genetic effects of introgression genomic components from Sea Island cotton (Gossypium barbadense L.) on fiber related traits in upland cotton (G. hirsutum L.). Euphytica 181, 41–53. doi: 10.1007/s10681-011-0378-1
Wang, F., Xu, Z., Sun, R., Gong, Y., Liu, G., Zhang, J., et al. (2013). Genetic dissection of the introgressive genomic components from Gossypium barbadense L. that contribute to improved fiber quality in Gossypium hirsutum L. Mol. Breed. 32, 547–562. doi: 10.1007/s11032-013-9888-y
Wang, F., Zhang, C., Liu, G., Chen, Y., Zhang, J., Qiao, Q., et al. (2016). Phenotypic variation analysis and QTL mapping for cotton (Gossypium hirsutum L.) fiber quality grown in different cotton-producing regions. Euphytica 211, 169–183. doi: 10.1007/s10681-016-1728-9
Wang, P., Ding, Y., Lu, Q., Guo, W., and Zhang, T. (2008). Development of Gossypium barbadense chromosome segment substitution lines in the genetic standard line TM-1 of Gossypium hirsutum. Sci. Bull. 53, 1512–1517. doi: 10.1007/s11434-008-0220-x
Wang, Q., Fang, L., Chen, J., Hu, Y., Si, Z., Wang, S., et al. (2015). Genome-wide mining, characterization, and development of microsatellite markers in Gossypium species. Sci. Rep. 5:10638. doi: 10.1038/srep10638
Wang, S., Basten, C., Gaffney, P., and Zeng, Z. (2007). Windows QTL Cartographer 2.0. Raleigh: Bioinformatics Research Center.
Xia, Z., Su, X., Liu, J., and Wang, M. (2013). The RING-H2 finger gene 1 (RHF1) encodes an E3 ubiquitin ligase and participates in drought stress response in Nicotiana tabacum. Genetica 141, 11–21. doi: 10.1007/s10709-013-9702-0
Xie, X., Jin, F., Song, M.-H., Suh, J.-P., Hwang, H.-G., Kim, Y.-G., et al. (2008). Fine mapping of a yield-enhancing QTL cluster associated with transgressive variation in an Oryza sativa × O. rufipogon cross. Theor. Appl. Genet. 116, 613–622. doi: 10.1007/s00122-007-0695-x
Xu, P., Zhu, J., Zhang, X., Ni, W., Xu, Y., and Shen, X. (2012). Molecular mapping and identification of QTLs for fiber micronaire on chromosome 7 from Gossypium klotzschianum. Acta Agron. Sin. 38, 447–453. doi: 10.3724/SP.J.1006.2012.00447
Ye, W., Wang, Y., and Dou, D. (2010). SeqHunter: a bioinformatics toolbox for local Blast and sequence analysis. China J. Bioinform. 8, 364–367.
Yu, J., Jung, S., Cheng, C. H., Ficklin, S. P., Lee, T., Zheng, P., et al. (2014). CottonGen: a genomics, genetics and breeding database for cotton research. Nucleic Acids Res. 42, D1229–D1236. doi: 10.1093/nar/gkt1064
Yu, J., Yu, S., Lu, C., Wang, W., Fan, S., Song, M., et al. (2007). High-density linkage map of cultivated allotetraploid cotton based on SSR, TRAP, SRAP and AFLP markers. J. Integr. Plant Biol. 49, 716–724. doi: 10.1016/j.gene.2015.08.022
Yu, J., Zhang, K., Li, S., Yu, S., Zhai, H., Wu, M., et al. (2013). Mapping quantitative trait loci for lint yield and fiber quality across environments in a Gossypium hirsutum x Gossypium barbadense backcross inbred line population. Theor. Appl. Genet. 126, 275–287. doi: 10.1007/s00122-012-1980-x
Zhang, J., Guo, W., and Zhang, T. (2002). Molecular linkage map of allotetraploid cotton (Gossypium hirsutum L. × Gossypium barbadense L.) with a haploid population. Theor. Appl. Genet. 105, 1166–1174. doi: 10.1007/s00122-002-1100-4
Zhang, K., Zhang, J., Ma, J., and Tang, S. (2012). Genetic mapping and quantitative trait locus analysis of fiber quality traits using a three-parent composite population in upland cotton (Gossypium hirsutum L.). Mol. Breed. 29, 335–348. doi: 10.1007/s11032-011-9549-y
Zhang, S. W., Zhu, X. F., Feng, L. C., Xiang, G., Yang, B., Zhang, T. Z., et al. (2016). Mapping of fiber quality QTLs reveals useful variation and footprints of cotton domestication using introgression lines. Sci. Rep. 6:31954. doi: 10.1038/srep31954
Zhang, T., Hu, Y., Jiang, W., Fang, L., Guan, X., Chen, J., et al. (2015). Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 33, 531–537. doi: 10.1038/nbt.3207
Zhang, Z., Rong, J., Waghmare, V. N., Chee, P. W., May, O. L., Wright, R. J., et al. (2011). QTL alleles for improved fiber quality from a wild Hawaiian cotton, Gossypium tomentosum. Theor. Appl. Genet. 123, 1075–1088. doi: 10.1007/s00122-011-1649-x
Zhang, Z. S., Hu, M. C., Zhang, J., Liu, D. J., Zheng, J., Zhang, K., et al. (2009). Construction of a comprehensive PCR-based marker linkage map and QTL mapping for fiber quality traits in upland cotton (Gossypium hirsutum L.). Mol. Breed. 24, 49–61. doi: 10.1007/s11032-009-9271-1
Zhao, L., Lv, Y., Cai, C., Tong, X., Chen, X., Zhang, W., et al. (2012). Toward allotetraploid cotton genome assembly: integration of a high-density molecular genetic linkage map with DNA sequence information. BMC Genomics 13:539. doi: 10.1186/1471-2164-13-539
Keywords: fiber quality, introgression alleles, secondary mapping populations, quantitative trait loci (QTL), QTL cluster
Citation: Chen Y, Liu G, Ma H, Song Z, Zhang C, Zhang J, Zhang J, Wang F and Zhang J (2018) Identification of Introgressed Alleles Conferring High Fiber Quality Derived From Gossypium barbadense L. in Secondary Mapping Populations of G. hirsutum L. Front. Plant Sci. 9:1023. doi: 10.3389/fpls.2018.01023
Received: 16 January 2018; Accepted: 25 June 2018;
Published: 18 July 2018.
Edited by:
Petr Smýkal, Palacký University Olomouc, CzechiaReviewed by:
Mehboob-ur- Rahman, National Institute for Biotechnology and Genetic Engineering, PakistanXinlian Shen, Jiangsu Academy of Agricultural Sciences (JAAS), China
Baoliang Zhou, Nanjing Agricultural University, China
Copyright © 2018 Chen, Liu, Ma, Song, Zhang, Zhang, Zhang, Wang and Zhang. 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: Furong Wang, wfr1125@126.com Jun Zhang, scrczj@saas.ac.cn; zj0928@126.com