Advanced Backcross QTL Analysis of Fiber Strength and Fineness in a Cross between Gossypium hirsutum and G. mustelinum

The molecular genetic basis of cotton fiber strength and fineness in crosses between Gossypium mustelinum and Gossypium hirsutum (Upland cotton) was dissected using 21 BC₃F₂ and 12 corresponding BC₃F_2:3 and BC₃F_2:4 families. The BC₃F₂ families were genotyped with simple sequence repeat markers from a G. hirsutum by G. mustelinum linkage map, and the three generations of BC₃-derived families were phenotyped for fiber strength (STR) and fineness (Micronaire, MIC). A total of 42 quantitative trait loci (QTLs) were identified through one-way analysis of variance, including 15 QTLs for STR and 27 for MIC, with the percentage of variance explained by individual loci averaging 13.86 and 14.06%, respectively. Eighteen of the 42 QTLs were detected at least twice near the same markers in different generations/families or near linked markers in the same family, and 28 of the 42 QTLs were identified in both mixed model-based composite interval mapping and one-way variance analyses. Alleles from G. mustelinum increased STR for eight of 15 and reduced MIC for 15 of 27 QTLs. Significant among-family genotypic effects (P < 0.001) were detected in 13 and 10 loci for STR and MIC respectively, and five loci showed significant (P < 0.001) genotype × family interaction for MIC. These results support the hypothesis that fiber quality improvement for Upland cotton could be realized by introgressing G. mustelinum alleles although complexities due to the different effects of genetic background on introgressed chromatin might be faced. Building on prior work with G. barbadense, G. tomentosum, and G. darwinii, QTL mapping involving introgression of G. mustelinum alleles offers new allelic variation to Upland cotton germplasm.

Introduction

This is the fourth report describing the interspecific G. mustelinum by G. hirsutum (Upland cotton) genetic map and mapping and introgressing quantitative trait loci (QTLs) of fiber quality traits from G. mustelinum into Upland cotton. In three previous papers, we reported the interspecific G. hirsutum by G. mustelinum genetic linkage map (Wang et al., 2016b), and described 24 QTLs for fiber elongation (Wang et al., 2016a,b) and 65 QTLs for fiber length traits (Wang et al., 2017). Some alleles from G. mustelinum improved fiber properties, demonstrating their potential value for improving fiber quality in Upland cotton breeding.

Here interspecific QTL mapping and introgression were reported for two important fiber quality traits, fiber strength (STR) and fineness measured in Micronaire (MIC). STR is determined as the necessary force breaking a fiber “beard.” Following the measurement of length by the Uster High Volume Instrument (HVI), the measurement of strength is performed on the same fiber beard with two sets of jaws clamping the beard by using a gauge length of 3.175 mm. The breaking force is measured directly and normalized with an assessment of the fiber mass from the optical sensor (in combination with the MIC value) to give the strength in cN/tex (Naylor et al., 2014). Fiber tenacity usually affects yarn tenacity more than any other properties of fiber at optimum yarn twist, and when fiber strength increases one cN/tex, yarn strength will increase about 0.5 cN/tex or even more. Generally bundle tenacity higher than 30 cN/tex is considered desirable (Estur and Knappe, 2007). Usually STR should be increased in cotton breeding programs.

MIC is a measure of the air permeability of compressed cotton fibers, which reflects both maturity in the development degree of cell walls and fiber fineness measured in linear density. A constant cotton fiber mass is compressed into a space with known volume. This compressed sample is used to measure air permeability, which is converted to appropriate numbers denoting MIC values. Generally MIC readings of 3.7–4.2 are premium, 3.5–3.6 or 4.3–4.9 are base, and 3.4-and-under or 5.0-and-higher are substandard and may result in a discounted price to the producer (http://www.cottoninc.com/fiber/quality/Classification-Of-Cotton/Classing-booklet.pdf). When the measurement comes in too low, the cotton is more susceptible to entangling around debris, which means too much of the good fiber will also be lost. When it is too high, it also causes problems since a coarser fiber negatively affects the spinning process, as well as overall quality is undesirable from the aspect of yarn evenness and spinning (Montalvo, 2005). The fineness determines how many fibers are present in the cross-section of a yarn of given thickness. More fibers available in the cross section of yarn generally results in stronger yarn, usually produced with finer fibers (Estur and Knappe, 2007). Additional fibers in the cross-section provide not only additional strength but also better evenness in the yarn. Many Upland cotton varieties have the problem of high MIC values currently, so it is necessary to breed varieties with comparatively low MIC or finer fiber.

This research aims to map QTLs for STR and MIC in a set of advanced-backcross G. hirsutum × G. mustelinum population. The QTLs mapped in this research will enhance our understanding of the molecular genetic basis of cotton fiber quality. They will also benefit cotton molecular breeding to improve STR and MIC with G. mustelinum alleles and to ascertain the specific genetic basis of these important traits.

Materials and Methods

Population Development and Field Evaluation

Three generations of interspecific advanced- backcross populations, namely 21 BC₃F₂ families and 12 BC₃F_2:3 and BC₃F_2:4 families were developed as follows: a G. hirsutum acc., PD94042 was crossed with G. mustelinum (AD4-8), and then F₁ plants were independently backcrossed to the G. hirsutum parent for three cycles; a total of 21 lineages that produced BC₃F₁ plants were self-pollinated and generated 21 BC₃F₂ families with size ranging from 127 to 160 plants (totally 3,203 BC₃F₂ progenies; Wang et al., 2016a); in addition, 12 BC₃F_2:3 and BC₃F_2:4 families with size ranging from 130 to 160 lines (totally 1,826 lines) were developed for 12 of the 21 BC₃F₂ families with enough seeds (Table S1). The 21 BC₃F₂ families were planted in 2006; completely randomized designs were used for the 12 BC₃F_2:3 and BC₃F_2:4 families with two replicate plots in two years (2008 and 2009) in Tifton, Georgia. All cultural practices were performed as described in Wang et al. (2016a). For the BC₃F₂ generation, seed-cotton of all bolls was hand-harvested for each plant; for the BC₃F_2:3 and BC₃F_2:4 generations, seed-cotton was hand-picked from two random replicate plots as two samples. Seed-cottons were ginned on a saw gin, and then STR and MIC were tested by the Cotton Incorporated Textile Services Laboratory (Cotton Incorporated, Cary, N.C.) by using the Uster High Volume Instrument (HVI; Zellweger-Uster, Knoxville, Tenn.).

Genotyping and Data Analysis

The BC₃F₁ plants were genotyped with 218 SSR markers selected for even representation of our interspecific G. mustelinum by G. hirsutum map which comprised 1,055 loci (Wang et al., 2016b), constructed with an F₂ population of the same two parents. The markers with introgression from G. mustelinum in the BC₃F₁s were then screened in the corresponding BC₃F₂ families for genotyping (Wang et al., 2016a). An average of 58 markers were used to genotype each BC₃F₂ family, with marker numbers ranging from 47 to 81, and the size of BC₃F₂ families ranged from 127 to 160. This set of genotype data based on each BC₃F₂ individual also constitutes the genotype of the corresponding BC₃F_2:3 and BC₃F_2:4 lines.

For every marker locus that segregated within BC₃F₂ families, one-way variance analyses were used to test associations between phenotypes and marker genotypes for statistical significance. The GLM procedure in the SAS ver.8 software package (SAS Institute, 1999) was used to carry out the analyses, and the significance threshold was set at P < 0.001 for F-test. The gene action modes (additive, a and dominant, d) for individual QTLs were evaluated with their significance levels estimated as described by Paterson et al. (1990). The dominance/additive (d/a) ratio of 3 was used as the threshold to determine whether the QTLs were over- or under-dominant (Chee et al., 2005a).

QTLs were also detected with the software QTLNetwork V2.1 (Yang et al., 2008) for STR and MIC in each of the BC₃-derived families, so as to map epistatic QTLs and also help confirm the reliability of the QTLs identified by one-way variance analyses. The critical F value calculated based on 1,000 permutation tests was used in the mixed model-based composite interval mapping (MCIM) method, with walk speed and window size set at 1 and 10 cM, respectively. A putative main-effect or epistatic QTL was claimed with the significance threshold set at P = 0.001. Considering environmental effects, QTLs were also mapped in joint analysis for the 12 families that were grown in three different environments/generations (BC₃F₂, BC₃F_2:3, and BC₃F_2:4) with the software QTLNetwork V2.1 (Yang et al., 2008). QTLs sharing a common marker between the two methods (QTLNetwork 2.1 and one-way variance analyses) were regarded as the same QTLs. QTL Nomenclature of QTLs was performed as previously described by McCouch et al. (1997); the QTL name began with a “q,” representing a QTL, followed by an abbreviation of the trait name, the chromosome name, and consecutive numbers indicating the QTL number of the same trait on the same chromosome (Wang et al., 2017).

For the loci that segregated in two or more families, two-way mixed model variance analyses were also applied, using the MIXED procedure of the SAS ver. 8 package, including genotype (G) as a fixed factor and family (F) and genotype × family (G × F) interaction as random factors, and the residual maximum likelihood (REML) method was used to estimate model parameters. Marker-trait association (genotype factor) was evaluated with an F statistic using a general Satterthwaite approximation for the denominator degrees of freedom (SAS Institute, 1999). A likelihood-ratio (ChiSq, χ²) test was carried out for the G × F interaction (Self and Liang, 1987; Chee et al., 2005a). P < 0.001 was set as the significance level for both G and G × F effects.

Results

Phenotypic Distribution and Correlations

The distributions of STR and MIC of the BC₃ progenies are shown in Figure 1 and Table 1. Both traits expressed significant segregation in the three generations, and their distribution ranges were wider in BC₃F₂ than BC₃F_2:3 and BC₃F_2:4 (Table 1). Although G. mustelinum, the wild parent cannot produce spinnable fiber, many BC₃ progenies show higher STR and lower MIC (usually lower MIC is preferred) than the cultivated parent (G. hirsutum, Table 1). For STR, five of 21 BC₃F₂, seven of 12 BC₃F_2:3, and nine of 12 BC₃F_2:4 families showed higher mean values than the G. hirsutum parent; four families, namely POP20, POP27, POP34, and POP35 had higher mean values than the G. hirsutum parent in all three generations (Figure 1). For MIC, 16 of 21 BC₃F₂, seven of 12 BC₃F_2:3, and also seven of 12 BC₃F_2:4 families showed mean values lower than that of the G. hirsutum parent; seven families, namely POP11, POP16, POP17, POP20, POP27, POP32, and POP35 showed lower mean values than the G. hirsutum parent in all three generations (Figure 1). Correlations were calculated to evaluate relationships between STR and MIC. The correlation coefficients were −0.005, −0.129, and −0.206 in BC₃F₂, BC₃F_2:3, and BC₃F_2:4 respectively, and the correlations reached significant level of P < 0.01 in BC₂F_2:3 and BC₃F_2:4.

FIGURE 1

Figure 1. Phenotypic distribution of fiber strength (STR) and fineness (Micronaire, MIC) demonstrated by family mean deviation from recurrent Gossypium hirsutum (GH) parent.

TABLE 1

Table 1. Summary statistics of fiber strength and fineness measured on the G. hirsutum parent and BC₃ progenies.

Main-Effect QTLs Detected for Each Trait

By assuming that each block of linked markers showing significant marker-trait association (P < 0.001) within a family denoted a single QTL, a total of 42 non-overlapping QTLs were identified in the three generations of BC₃-derived families (Figure 2 and Table 2). These QTLs were mapped to 20 chromosomes, with 22 located on 12 A-subgenome chromosomes and 20 located on eight D-subgenome chromosomes. One or more QTLs for STR and MIC were identified in 13 of the 21 families, with a maximum of five QTLs in each of three families (POP12, POP15, and POP32). Eighteen of the 42 QTLs could be identified at least twice near the same markers in different generations/families or near linked markers in the same family. Twenty-eight of the 42 QTLs were also identified by QTLNetwork. The detailed QTL information is listed in Figure 2 and Table 2.

FIGURE 2

Figure 2. QTLs of fiber strength (STR) and fineness (Micronaire, MIC). The solid black bars indicated QTLs of STR, and the green bars filled with slashes indicated QTLs of MIC.

TABLE 2

Table 2. Biometrical parameters of QTLs affecting fiber strength and fineness.

QTLs for STR

A total of 15 non-overlapping QTLs were identified on 12 chromosomes for STR (Figure 2 and Table 2), with seven located in the A-subgenome and eight located in the D-subgenome. Seven QTLs were found at least twice near the same markers in different generations/families or near linked markers in the same family, namely qSTR-2-1, qSTR-11-1, qSTR-19-1, qSTR-19-2, qSTR-22-1, qSTR-22-2, and qSTR-25-1. The percentage of variance explained (PVE) by individual QTLs ranged from 8.74% (qSTR-23-1) to 29.22% (qSTR-11-1), with an average of 13.86%. G. mustelinum alleles increased STR for eight of the 15 QTLs (Table 2). Ten of the 15 QTLs could also be identified by the MCIM method of QTLNetwork (Table 2).

QTLs for MIC

A total of 27 non-overlapping QTLs for MIC were identified on 15 chromosomes with 15 located on nine A-subgenome chromosomes and 12 located on six D-subgenome chromosomes (Figure 2 and Table 2). Eleven QTLs were found at least twice near the same markers in different generations/families or near linked markers in the same family. The PVE of each individual locus ranged from 7.92% (qMIC-5-1) to 22.95% (qMIC-25-1), with an average of 14.06%. For 15 of the 27 QTLs, G. mustelinum alleles reduced MIC and therefore contributed to finer fiber (Table 2). Eighteen of the 27 QTLs could also be identified by the MCIM method of QTLNetwork (Table 2).

Consistency of QTLs across Families

A total of 211 of the 218 SSR marker loci were segregating in two or more families; consequently, two-way ANOVA was performed to identify marker-trait associations, so as to evaluate their consistency among different families. Significant (P < 0.001) among-family G effects were identified at 13 and 10 loci for STR and MIC respectively (Table S2). For STR, the 13 loci appear to represent eight non-overlapping genomic regions, and QTLs were identified in four regions in within-family analysis (qSTR-5-1, qSTR-6-1, qSTR-19-2, and qSTR-23-1). For MIC, the 10 loci represent three non-overlapping genomic regions, for which QTLs were detected in three regions in within-family analysis (qMIC-5-1, qMIC-7-1, and qMIC-16-1; Table 2, Table S2).

For MIC, a total of five loci were significant (P < 0.001) for G × F interactions (Table S3). QTLs were detected at all five loci (qMIC-4-1, qMIC-4-2, qMIC-10-1, qMIC-15-2, qMIC-19-3, and qMIC-25-2; Table 2, Table S3) in different segregating families. No locus was significant for G × F interactions for STR.

Epistatic QTLs and Their Interactions with Environments

A total of 13 epistatic QTLs were identified with significant additive × additive (AA) effects (P < 0.001), most of which (10/13) involved loci that were not linked to any main-effect QTLs (Table 3). For STR, the interaction between a region on Chr5 and another region on Chr19 was identified in both BC₃F_2:3 and joint analysis in POP34 simultaneously. G. mustelinum alleles increased STR (with negative AA effects) for three epistatic QTLs, and alleles from G. mustelinum decreased MIC (with positive AA effects) for five epistatic QTLs (Table 3). The interactions between epistatic QTLs with environment were listed in Table 3.

TABLE 3

Table 3. Estimated epistasis and epistasis × environment interaction effects of QTLs for fiber strength and fineness.

Discussion

Cultivated Gossypium hirsutum has a narrow gene pool, having experienced genetic bottlenecks during polyploid formation and divergence from its sister polyploid species, as well as during domestication, dispersal by humans, and scientific breeding. The domestication of a comparatively small subgroup of wild species and, in more recent years, over-exploitation of only a few genetic backgrounds in breeding programs of modern Upland cotton by crossing a few closely-related genotypes repeatedly to develop new cultivars has resulted in a genetically depauperate cotton germplasm. Slow genetic progress in improving fiber productivity and quality is indicative that many favorable alleles have reached fixation in the elite gene pool. The deficiency of genetic variation in current germplasm has enhanced the difficulty for breeders to provide low-cost intrinsic genetic solutions to cotton fiber production such as resistance to abiotic and biotic hazards or new needs in fiber quality or agronomic traits.

Because the narrowing of the cotton genetic base, new sources of genetic variation need to be introduced into the cotton gene pool to ensure future success in breeding new cotton cultivars. As a wild tetraploid cotton species diverged far from G. hirsutum (Wendel and Cronn, 2003) but sexually compatible with cultivated cotton, G. mustelinum may harbor elite alleles beneficial to the improvement of fiber quality traits in Upland cotton (Alves et al., 2013).

In this research, advanced backcross QTL (AB-QTL) analysis was carried out based on SSR markers and phenotypic data collected from three generations of BC₃-derived families with introgression from G. mustelinum. Phenotypic assessment of the advanced backcross populations indicated significant segregation for STR and MIC, which indicated that both negative and positive alleles existed for each trait in both parents (Table 1, Figure 1).

Since the wild parent, G. mustelinum cannot produce spinnable fiber, it is not a surprise to find out that for STR, many BC₃-derived families (16 of 21 in BC₃F₂, five of 12 in BC₃F_2:3, and three of 12 in BC₃F_2:4) showed mean STR values lower than the recurrent G. hirsutum parent, PD94042; for MIC, there were five BC₃-derived families in each generation having higher mean MIC values or coarser fiber than the recurrent parent. This “negative” transgression, yielding a poorer phenotype than that of the recurrent parent, suggests that interspecific hybridization formed many undesirable new gene combinations. Nonetheless, for STR, five of 21 BC₃F₂, seven of 12 BC₃F_2:3, and nine of 12 BC₃F_2:4 families showed mean values higher than that of the G. hirsutum parent, and four families (POP20, POP27, POP34, and POP35) had higher mean values than the G. hirsutum parent in all three generations; for MIC, most BC₃-derived families, namely 16 of 21 BC₃F₂, seven of 12 BC₃F_2:3, and also seven of 12 BC₃F_2:4 families outperformed the recurrent parent with lower MIC (finer fiber: Figure 1), and seven families (POP11, POP16, POP17, POP20, POP27, POP32, and POP35) showed mean values lower than the G. hirsutum parent in the three generations, showing good stability across environments. Many individual plants/lines in these families have better STR and MIC traits than those of the recurrent parent (Figure 1). This “positive” transgression, yielding a superior phenotype than that of the recurrent parent, suggests that interspecific hybridization formed some desirable new gene combinations, which was also found in previous reports on introgression of G. barbadense (Chee et al., 2005a,b; Draye et al., 2005), G. tomentosum (Zhang et al., 2011), and G. darwinii (Wang et al., 2012). As suggested by previous research (Jiang et al., 2000), gene transfer between gene pools is a significant consequence of interspecific hybridization, which will increase selectable genetic variation and introduce genes for adaptive traits. These results support the hypothesis that fiber quality improvement for Upland cotton may be accomplished by introgressing elite genes from G. mustelinum and other tetraploid cotton species.

A total of 15 and 27 non-overlapping QTLs were mapped in one-way analysis of variance with PVE of 13.86 and 14.06% on average for STR and MIC respectively. The effects of many QTLs showed good reproducibility, with 18 of the 42 QTLs detected at least twice near the same markers in different generations/families or near linked markers in the same family. In addition, 28 of the 42 QTLs were also identified by QTLNetwork (Table 2). The detection of QTLs near different markers at corresponding chromosomal locations or across various generations/families with different methods further supports the likelihood that these QTLs are real.

Alleles from the wild G. mustelinum increased STR for 53% (eight) of 15 QTLs and decreased MIC (conferred finer fiber) for 56% (15) of 27 QTLs. These QTLs are of great importance to be deeply exploited to transfer elite genes from G. mustelinum into Upland cotton. Efforts to improve STR and MIC are now in progress by constructing near-isogenic lines allowing these alleles from G. mustelinum to be more easily manageable in cotton breeding programs.

In addition to main-effect QTLs, epistatic QTLs were identified for STR and MIC using data collected over three different generations (Table 3). The results indicated that both epistatic QTLs and main-effect play key genetic roles in STR and MIC (Tables 2, 3). Epistatic QTLs are more complicated to manage compared to utilization of main-effect alleles, naturally at least twice the difficulty will be encountered in testing and introgressing of epistatic QTLs. Furthermore, it is interesting that the majority of epistatic interactions (10/13) were identified between genetic loci that were not linked to any QTL (Table 3), similar to our previous results derived from fiber length, where 14 of 17 epistatic QTLs involved loci not linked to any main-effect QTLs (Wang et al., 2017), suggesting high complexity of fiber quality inheritance. In joint analysis, effects of epistasis × environment were detected for both STR and MIC (Table 3). Although, generally the PVE of epistasis × environment was smaller than that of epistasis, epistasis × environment interactions may add to the difficulty of breeding.

Since more than one family was often segregating for the same chromosomal segment, it was possible to explore genetic background effects on introgressed chromatin. For STR and MIC, significant (P < 0.001) among-family G effects were identified at 13 and 10 loci (Table S2), with four and three loci revealing QTLs in within-family analysis for STR (qSTR-5-1, qSTR-6-1, qSTR-19-2, and qSTR-23-1) and MIC (qMIC-5-1, qMIC-7-1, and qMIC-16-1). Some among-family G effects demonstrated good reproducibility, with three of 13 for STR, and two of 10 for MIC detected in different generations (Table S2). A total of five loci were significant (P < 0.001) for G × F interactions for MIC (Table S3), with QTLs detected in all the five loci (qMIC-4-1, qMIC-4-2, qMIC-10-1, qMIC-15-2, qMIC-19-3, and qMIC-25-2; Table 2, Table S3). The most extreme case of G × F interaction was detected at the locus MUSB1050 on Chr4 (Table 2, Table S3). In family POP32, alleles from G. mustelinum at this locus conferred a decrease of 0.41 for additive effect in MIC that accounted for 21.45% of PVE. Interestingly, this same locus conferred an additive increase of 0.09 for MIC that accounted for 10.78% of PVE in family POP35; this locus also segregated in family POP10, POP11, and POP27 but showed no significant association with MIC.

This study also adds to prior information on the significant influence of the tetraploid D-subgenome on fiber quality traits, although the D-subgenome was derived from a diploid ancestor without the ability to produce spinnable fiber (Jiang et al., 1998; Chee et al., 2005b; Zhang et al., 2008, 2011). In this current study, among the 42 QTLs affecting STR and MIC traits, the D-subgenome (20) had slightly fewer QTLs than the A-subgenome (22). Considering the QTLs for fiber elongation and fiber length detected in our previous reports (Wang et al., 2016a,b, 2017), 73 QTLs were identified in the D-subgenome, more than the A-subgenome (58), which collectively supports the finding that for fiber quality traits, more QTLs occurred on the D-subgenome than the A-subgenome(Jiang et al., 1998; Paterson et al., 2003).

With this fourth report on QTLs for fiber traits from G. mustelinum, this series of papers collectively describes six fiber quality traits investigated in 21 BC₃F₂, 12 BC₃F_2:3, and 12 BC₃F_2:4 families, namely fiber elongation (EL), fiber uniformity index (UI), upper-half mean length (UHM), short fiber content (SFC), STR, and MIC. Mean values of some families outperformed the recurrent G. hirsutum parent, PD94042, for each trait in each generation; likewise, many individual plants/lines showed superior fiber quality performance than the G. hirsutum parent, showing promise that our goal of introgressing alleles from G. mustelinum to improve Upland cotton may work.

Two loci (JESPR224 and BNL3264) may be of special interest in this series of reports. Two QTLs were detected near the locus JESPR224 (qMIC-25-1 in POP01 in BC₃F₂, and qSTR-25-2 in POP02 in BC₃F₂), where G. mustelinum alleles decreased MIC and increased STR. Four QTLs (qUHM-25-1, qSTR-25-1, qMIC-25-1, and EL25.1) were detected near the loci BNL3264 on Chr25 in POP17, for which G. mustelinum alleles increased EL, UHM, and STR and decreased MIC (Wang et al., 2016a,b, 2017). These QTLs hold promise for improving cotton fiber elongation, length, strength and fineness simultaneously through marker-assisted selection.

More co-locations of QTLs with both desirable and undesirable effects on different traits were observed. For instance, two QTLs (qELO-19-1 and qSTR-19-2) were detected in POP15 near the locus BNL3811 on Chr19, for which G. mustelinum alleles increased STR but decreased EL; three QTLs (qMIC-12-1, qSFC-12-2, and qUI-12-3) were detected near the locus DPL0866 on Chr12 in POP11, for which G. mustelinum alleles decreased MIC but increased SFC and decreased UI; four QTLs (qUI-10-1, qSFC-10-2, qMIC-10-1, and qELO-10-1) were detected near the locus JESPR6 on Chr10, for which G. mustelinum alleles increased EL but they also increased MIC and SFC and decreased UI (Wang et al., 2016a,b, 2017). The co-location of QTLs with opposite effects on different traits indicates the difficulty of improving diverse traits simultaneously in breeding programs, and may also account for the challenges that have been faced with using exotic germplasm in the absence of DNA marker information.

Building on prior work on G. barbadense, G. tomentosum, and G. darwinii, QTL mapping involving introgression of G. mustelinum alleles offers new allelic variation to the Upland cotton gene pool; in addition, the new germplasm created here offers an opportunity for the cotton community to explore G. mustelinum alleles in an elite cultivated background, and also provided materials potentially useful in cotton breeding programs. Wild cotton species represent a repository of divergent and in some cases favorable alleles for a variety of traits including fiber quality. The sixth and seventh tetraploid cotton species found and confirmed, G. ekmanianum Wittmack (endemic to the Dominican Republic; Krapovickas and Seijo, 2008; Wendel and Grover, 2015) and Gossypium sp. nov. (found from two islands, Wake and Peale in the Wake Atoll in the Pacific Ocean; Wendel and Grover, 2015) offer additional scope for exploration of gene introgression.

Author Contributions

AP conceived and designed the experiments; PC and OM oversaw crossing, phenotyping, and genotyping, which were performed by BW, ZhiZ, ZheZ and EL; DJ oversaw fiber analysis; BW, XD, LS, and TS analyzed the data; BW, AP, and PC wrote the paper.

Funding

This work was supported by National Science Foundation (PFI award IIP 0917856; AIR award IIP 1127755), Cotton Incorporated, Key Research and Development Project of Jiangsu Province, China (Modern Agriculture, BE2015353), the Applying Basic Research Fund of Nantong City-Agricultural and Social Development (MS12016043).

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

We thank Jennifer McCurdy for help in field work and Cornelia Lemke for technical assistance in the lab.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.01848/full#supplementary-material

References

Alves, M., Barroso, P., Ciampi, A., Hoffmann, L., Azevedo, V., and Cavalcante, U. (2013). Diversity and genetic structure among subpopulations of Gossypium mustelinum (Malvaceae). Genet. Mol. Res. 12, 597–609. doi: 10.4238/2013.February.27.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chee, P., Draye, X., Jiang, C. X., Decanini, L., Delmonte, T. A., Bredhauer, R., et al. (2005a). Molecular dissection of interspecific variation between Gossypium hirsutum and Gossypium barbadense (cotton) by a backcross-self approach: I. Fiber elongation. Theor. Appl. Genet. 111, 757–763. doi: 10.1007/s00122-005-2063-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Chee, P. W., Draye, X., Jiang, C. X., Decanini, L., Delmonte, T. A., Bredhauer, R., et al. (2005b). Molecular dissection of phenotypic variation between Gossypium hirsutum and Gossypium barbadense (cotton) by a backcross-self approach: III. Fiber length. Theor. Appl. Genet. 111, 772–781. doi: 10.1007/s00122-005-2062-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Draye, X., Chee, P., Jiang, C. X., Decanini, L., Delmonte, T. A., Bredhauer, R., et al. (2005). Molecular dissection of interspecific variation between Gossypium hirsutum and G. barbadense (cotton) by a backcross-self approach: II. Fiber fineness. Theor. Appl. Genet. 111, 764–771. doi: 10.1007/s00122-005-2061-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Estur, G., and Knappe, M. (2007). Cotton Exporter's Guide. Geneva: The International Trade Centre (ITC), Palais des Nations.

Jiang, C. X., Chee, P. W., Draye, X., Morrell, P. L., Smith, C. W., and Paterson, A. H. (2000). Multilocus interactions restrict gene introgression in interspecific populations of polyploid Gossypium (cotton). Evolution 54, 798–814. doi: 10.1111/j.0014-3820.2000.tb00081.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C. X., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Krapovickas, A., and Seijo, G. (2008). Gossypium ekmanianum (Malvaceae), algodon silvestre de la Republica Dominicana. Bonplandia 17, 55–63. Available online at: http://ibone.unne.edu.ar/objetos/up/documentos/bonplandia/public/17_1/55_63.pdf

Google Scholar

McCouch, S., Cho, Y., Yano, M., Paul, E., Blinstrub, M., Morishima, H., et al. (1997). Report on QTL nomenclature. Rice Genet. Newsl. 14:11.

Google Scholar

Montalvo, J. Jr. (2005). Relationships between micronaire, fineness, and maturity. I. Fundamentals. J. Cotton Sci. 9, 81–88. Available online at: http://www.leadership.cotton.org/journal/2005-09/2/upload/jcs09-081.pdf

Google Scholar

Naylor, G. R., Delhom, C. D., Cui, X., Gourlot, J.-P., and Rodgers, J. (2014). Understanding the influence of fiber length on the High Volume Instrument™ measurement of cotton fiber strength. Textile Res. J. 84, 979–988. doi: 10.1177/0040517513515318

CrossRef Full Text | Google Scholar

Paterson, A. H., Deverna, J. W., Lanini, B., and Tanksley, S. D. (1990). Fine mapping of quantitative trait loci using selected qverlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics 124, 735–742.

PubMed Abstract | Google Scholar

Paterson, A. H., Saranga, Y., Menz, M., Jiang, C. X., and Wright, R. J. (2003). QTL analysis of genotype x environment interactions affecting cotton fiber quality. Theor. Appl. Genet. 106, 384–396. doi: 10.1007/s00122-002-1025-y

PubMed Abstract | CrossRef Full Text | Google Scholar

SAS Institute (1999). SAS/STAT User's Guide, Version 8. Cary, NC: SAS Institute.

Self, S. G., and Liang, K.-Y. (1987). Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. J. Am. Stat. Assoc. 82, 605–610. doi: 10.1080/01621459.1987.10478472

CrossRef Full Text | Google Scholar

Wang, B., Draye, X., Zhang, Z., Zhuang, Z., May, O. L., Paterson, A. H., et al. (2016a). Advanced backcross QTL analysis of fiber elongation in a cross between Gossypium hirsutum and G. mustelinum. Crop Sci. 56, 1760–1768. doi: 10.2135/cropsci2015.12.0753

CrossRef Full Text

Wang, B., Draye, X., Zhuang, Z., Zhang, Z., Liu, M., Lubbers, E. L., et al. (2017). QTL analysis of cotton fiber length in advanced backcross populations derived from a cross between Gossypium hirsutum and G. mustelinum. Theor. Appl. Genet. 6, 1297–1308. doi: 10.1007/s00122-017-2889-1

CrossRef Full Text | Google Scholar

Wang, B., Liu, L., Zhang, D., Zhuang, Z., Guo, H., Qiao, X., et al. (2016b). A genetic map between Gossypium hirsutum and the Brazilian endemic G. mustelinum and its application to QTL mapping. G3 6, 1673–1685. doi: 10.1534/g3.116.029116

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Nie, Y. C., Lin, Z. X., Zhang, X. L., Liu, J. J., and Bai, J. (2012). Molecular diversity, genomic constitution, and QTL mapping of fiber quality by mapped SSRs in introgression lines derived from Gossypium hirsutum x G. darwinii Watt. Theor. Appl. Genet. 125, 1263–1274. doi: 10.1007/s00122-012-1911-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wendel, J. F., and Cronn, R. C. (2003). Polyploidy and the evolutionary history of cotton. Adv. Agron. 78, 139–186. doi: 10.1016/S0065-2113(02)78004-8

CrossRef Full Text | Google Scholar

Wendel, J. F., and Grover, C. E. (2015). “Taxonomy and evolution of the cotton genus, Gossypium,” in Cotton, 2nd Edn, eds D. Fang and R. Percy (Madison, WI: American Society of Agronomy, Inc.; Crop Science Society of America, Inc.; and Soil Science Society of America, Inc).

Google Scholar

Yang, J., Hu, C. C., Hu, H., Yu, R. D., Xia, Z., Ye, X. Z., et al. (2008). QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics 24, 721–723. doi: 10.1093/bioinformatics/btm494

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H. B., Li, Y., Wang, B., and Chee, P. W. (2008). Recent advances in cotton genomics. Intern. J. Plant Genomics 2008:742304. doi: 10.1155/2008/742304

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z. S., Rong, J. K., 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

PubMed Abstract | CrossRef Full Text | Google Scholar