Quantitative trait loci analysis of fiber quality traits using a random-mated recombinant inbred population in Upland cotton (Gossypium hirsutum L.)
© Fang et al.; licensee BioMed Central Ltd. 2014
Received: 7 March 2014
Accepted: 19 May 2014
Published: 24 May 2014
Upland cotton (Gossypium hirsutum L.) accounts for about 95% of world cotton production. Improving Upland cotton cultivars has been the focus of world-wide cotton breeding programs. Negative correlation between yield and fiber quality is an obstacle for cotton improvement. Random-mating provides a potential methodology to break this correlation. The suite of fiber quality traits that affect the yarn quality includes the length, strength, maturity, fineness, elongation, uniformity and color. Identification of stable fiber quantitative trait loci (QTL) in Upland cotton is essential in order to improve cotton cultivars with superior quality using marker-assisted selection (MAS) strategy.
Using 11 diverse Upland cotton cultivars as parents, a random-mated recombinant inbred (RI) population consisting of 550 RI lines was developed after 6 cycles of random-mating and 6 generations of self-pollination. The 550 RILs were planted in triplicates for two years in Mississippi State, MS, USA to obtain fiber quality data. After screening 15538 simple sequence repeat (SSR) markers, 2132 were polymorphic among the 11 parents. One thousand five hundred eighty-two markers covering 83% of cotton genome were used to genotype 275 RILs (Set 1). The marker-trait associations were analyzed using the software program TASSEL. At p < 0.01, 131 fiber QTLs and 37 QTL clusters were identified. These QTLs were responsible for the combined phenotypic variance ranging from 62.3% for short fiber content to 82.8% for elongation. The other 275 RILs (Set 2) were analyzed using a subset of 270 SSR markers, and the QTLs were confirmed. Two major QTL clusters were observed on chromosomes 7 and 16. Comparison of these 131 QTLs with the previously published QTLs indicated that 77 were identified before, and 54 appeared novel.
The 11 parents used in this study represent a diverse genetic pool of the US cultivated cotton, and 10 of them were elite commercial cultivars. The fiber QTLs, especially QTL clusters reported herein can be readily implemented in a cotton breeding program to improve fiber quality via MAS strategy. The consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying fiber development.
KeywordsCotton Fiber quality traits Microsatellite markers Quantitative trait loci Random-mating Recombinant inbred lines
Cotton is the most important natural fiber crop which supports a multi-billion dollar production and processing industry [1, 2]. The cotton genus (Gossypium L.) consists of about 45 diploid species belonging to eight genome groups (A-G and K) and 5 allotetraploid (AD) species [3, 4]. Within the Gossypium genus, two of the five tetraploid species (2n = 4x = 52; G. barbadense L. and G. hirsutum L.), along with two diploid species (2n = 2x = 26; G. arboreum L. and G. herbaceum L.), were independently domesticated for cotton fiber production in the last few thousand years in the New and Old World . Of these four cultivated species, the tetraploid species G. hirsutum, also referred to as “Upland cotton”, accounts for about 95% of the global cotton production. Consequently, a great majority of world-wide cotton breeding programs have been focusing on improving Upland cotton.
With the increasing global demand for textile products, intense competition from synthetic fibers, and textile industry’s modernization by shifting to high-speed spinning technologies, the need for higher yielding Upland cotton cultivars with improved fiber quality has never been more critical . However, yield increases are often negatively associated with fiber quality within Upland cotton [6–8]. This negative association between yield and fiber quality has hampered cotton breeding efforts for the improvement of multiple traits. It has been a cotton breeder’s high desire to effectively break this negative linkage. Random mating procedures have provided an important methodology to break undesirable associations and to form new combinations in several crop plants, including tobacco , sorghum , soybean , and oats . Random mating has previously been shown to reduce correlations between traits in cotton as well [8, 13–15]. Random mating requires a considerable expenditure of time and energy. If one starts with a large diverse group of parental lines, it offers an opportunity to break up adverse linkage blocks and to form new recombinations, some of which should be superior.
Conventional cotton breeding programs, primarily relying on crossing adapted genotypes and selecting novel allele combinations based on phenotypes, have greatly contributed to the success of the cotton industry in the past century [16, 17]. However, the very complex quantitative inheritance of yield and fiber quality traits and the negative associations between them as described above require cotton breeders to develop more effective strategies in order to develop superior cultivars.
Marker-assisted selection (MAS) is one such strategy that has received great attention among plant and animal breeders during the past three decades. In cotton, many reports on mapping qualitative and quantitative traits have been published [2, 7, 17–23]. As for mapping quantitative trait loci (QTL) related to fiber quality, more than 1000 QTLs have been reported so far (recently reviewed by Said et al. ) [2, 17, 19–21, 25–32]. However, most of these QTLs were obtained based on the analysis of interspecific (mainly between G. hirsutum and G. barbadense) populations. These fiber QTLs have very limited values in Upland cotton breeding through a MAS strategy because: 1) favorable QTLs (alleles) from G. barbadense are often not present in G. hirsutum; 2) markers showing polymorphism between two species may be monomorphic within G. hirsutum, which in turn makes the markers not useful in Upland cotton breeding; 3) transfer of favorable traits from G. barbadense to G. hirsutum via inter-mating has been very difficult, and resulted in very limited success [27, 33]. In realization of these limitations, cotton scientists have been using G. hirsutum intra-specific populations to identify fiber QTLs [2, 17, 34, 35]. However, this approach faces a challenge of low intra-specific polymorphism in cotton. The average polymorphism rate between any two Upland cotton cultivars used so far in a bi-parental mapping project is about 4–8% as revealed by microsatellite markers [36–38]. Currently, there is no a high density G. hirsutum intraspecific genetic map that covers the entire genome.
To overcome this challenge, cotton scientists used a four-way cross population  or three-parent composite population  to map fiber and yield QTLs. The approach using populations involving more than two parents ensured greater genetic diversity and an increased polymorphism frequency in the mapping populations, and improved possibility of QTL analysis. In addition, research on association mapping agronomic trait QTLs using a group of varieties has been rising in cotton [2, 39, 40].
In this research, we first developed a recombinant inbred population of 550 lines that were derived from six cycles of random-mating beginning with half diallel crossing of 11 Upland cotton cultivars, followed by six generations of self-pollination. The 550 recombinant inbred lines (RILs) were planted in two years to obtain fiber quality measurements. Second, we screened 15538 simple sequence repeat (SSR) markers to identify polymorphic markers among the 11 parents. Third, we analyzed 275 RILs (Set 1) using 1582 well-distributed SSR markers. Fourth, we conducted association analysis between markers and fiber quality traits using the software package TASSEL  in order to identify fiber QTLs. Finally, for each fiber trait, we confirmed the marker-trait associations using the other 275 RILs (Set 2) by comparing the phenotypic data of the RILs that were grouped based on the genotypes of seven QTL-linked markers. The main objective of this research was to identify stable QTLs related to fiber quality traits. Since the 11 parents represent a diverse germplasm pool of the US cotton cultivars, marker-trait associations identified in this research should be very useful in Upland cotton breeding aimed at improving fiber quality.
Genetic diversity of the 11 parents
Eleven Upland cotton cultivars used for random-mated RI population development
California Planting Cotton Seed Distributors (Shafter, CA)
Coker Pedigreed Seed Co. (Hartsville, SC)
Deltapine Acala 90
Delta and Pine Land Co. (Scott, MS)
Bayer Crop Science (Lubbock, TX)
USDA-ARS (Mississippi State, MS)
Paymaster Technologies, Inc. (Plainview, TX)
Phytogen PSC 355**
Phytogen Seeds (Indianapolis, IN)
Stoneville Pedigreed Seed Co. (Stoneville, MS)
Stoneville Pedigreed Seed Co. (Stoneville, MS)
Sure-Grow Co. (Centre, AL)
Texas A&M University (College Station, TX)
Statistics and genome distribution of 1582 polymorphic SSR markers
Statistics of the 1582 polymorphic markers as revealed in 11 cotton cultivars
Genomic distribution of 1582 polymorphic SSR markers used for genotyping the 275 RILs of Set 1
Total cM covered by markers
Whole chromosome (cM)*
% of the chromosome
#Gaps larger than 20 cM (largest gap)
RIL population structure assessment
Fiber quality measurements of parents and RILs
The 11 parents and RILs were planted in triplicates in Mississippi State, MS, USA for two years (2010—2011).
Fiber quality measurements of 550 RILs and the 11 parents (2010---2011 field results)
5.27 ± 0.15
4.60 ± 0.13
7.62 ± 0.32
29.86 ± 0.75
27.94 ± 0.51
83.48 ± 0.53
5.62 ± 0.14
4.65 ± 0.10
7.48 ± 0.17
30.43 ± 0.46
28.19 ± 0.25
83.70 ± 0.31
Heritability of fiber traits based on the measurements of 2010 and 2011 field samples
Heritability (H 2 )
Identification of fiber QTLs
Fiber QTLs identified by TASSEL and their effects
#Associated marker loci
Combined total QTL effects
Positive QTL effect
Negative QTL effect
ELO: Twenty-five QTLs were identified, and located on 17 chromosomes (Additional file 7). All together, the 11 positive QTLs explained 37.4% of the total phenotypic variance with an average of 3.40% per QTL. The QTL on Chr.16 (qELO-c16-1, marker locus DC40054a) had the largest effect, and explained 5.73% of the phenotypic variance. The 14 negative QTLs explained a combined 45.4% of the phenotypic variance with an average of 3.24% per QTL. The QTL on Chr.05 (qELO-c5-1, marker locus HAU0006a) had the largest negative effect.
MIC: A total of 24 QTLs were identified, and distributed on 16 chromosomes. Chr.26 had 4 QTLs, and all were able to reduce MIC value. Thirteen QTLs had positive effects, and explained a combined 38.5% of the phenotypic variance. The QTL on Chr.20 (qMIC-c20-2, marker locus SHIN-0170a) had the largest effect, and could increase the MIC value by 3.77%. There were 11 QTLs that could reduce MIC value by a total of 32.1%. The QTL on Chr.18 (qMIC-c18-1, marker locus DPL0807a) had the largest effect, and could reduce MIC value by 4.68%.
SFC: Nineteen QTLs were identified, and located on 15 chromosomes. Nine positive QTLs explained a combined 28.3% of the phenotypic variance with an average of 3.15% per QTL. The QTL on Chr.07 (qSFC-c7, marker locus C2-0114a) had the largest effect, and could increase SFC value by 5.21%. A total of 33.9% of the phenotypic variance could be explained by the 10 negative QTLs. The QTL on Chr.16 (qSFC-c16-1, marker locus CM0066a) was responsible for 5.51% of the phenotypic variance.
STR: Twenty QTLs residing on 14 chromosomes were identified. Ten QTLs had positive effects, and explained a combined 24.4% of the phenotypic variance. The QTL on Chr.16 (qSTR-c16, marker locus CM0066a) had the largest effect, and was responsible for 4.10% of the phenotypic variance. Ten QTLs had negative effects, and could explain 40.0% of total phenotypic variance. One major QTL on Chr.07 (qSTR-c7-1, marker locus C2-0114a) had significant effect, and could decrease STR by 14.89%, the largest among all the QTLs identified.
UHM: Twenty QTLs were identified, and located on 14 chromosomes. Ten QTLs each had positive and negative effects, respectively. The largest positive effect QTL was located on Chr.18 (qUHM-c18-1, marker locus TMB1208b), and could increase UHM by 4.15%. The QTL on Chr.22 (qUHM-c22-3, marker locus HAU0086b) had the largest negative effect, and explained 6.31% of the phenotypic variance.
UI: Twenty three QTLs were found as associated with UI, and resided on 17 chromosomes. Ten QTLs could increase the UI value, and had a combined effect of 33.4% on the phenotypic variance. The QTL on Chr.16 (qUI-c16-1, marker locus CM0066a) had the largest effect, and was responsible for 4.54% of the phenotypic variance. Thirteen QTLs could reduce UI value by as much as 43.2% if combined. The QTL on Chr.07 (qUI-c7-1, marker locus C2-0114a) had the largest negative effect, and was responsible for 7.28% of the phenotypic variance.
Of the 157 marker loci associated with fiber quality traits, 37 loci were associated with more than one trait. Eighteen, thirteen, four and two were associated with 2, 3, 4 and 5 traits, respectively (Additional file 6). No single locus was associated with all 6 traits. With exception of MIC which is a complex trait (combination of fiber maturity and fineness), the effect of each QTL cluster on the traits was usually similar. In other words, a QTL cluster that was usually able to increase the values of ELO, STR, UHM or UI but reduce the SFC value is considered as a favorable QTL from a breeding point of view. An unfavorable QTL cluster usually increased SFC value but decreased the values of all other four traits. It is worth mentioning that the QTL cluster on Chr.16 (marker locus CM0066a) had the largest favorable breeding effects on these four traits (SFC, STR, UHM, UI). On the contrary, the QTL cluster on Chr.07 (marker locus C2-0114a) had the largest unfavorable breeding effects on SFC, STR, UHM and UI.
Confirmation of fiber QTLs using the 275 RILs of Set 2
Intraspecific genetic diversity and marker coverage
Although the 11 parental cultivars used in the present research were quite diverse in agronomic performance and breeding pedigrees, the genetic diversity among them as revealed by the 15538 SSR markers was low. These SSR markers were developed by many groups around the world. The sources of the markers were either from genomic DNA or EST sequences of many genotypes . The 15538 SSR primer pairs used in screening were not pre-selected. And it is reasonable to believe that these markers cover the entire genome of tetraploid cotton as reported by Wang et al. . The present research further confirms that cultivated Upland cotton has very narrow genetic diversity and close kinships possibly because of a few bottlenecks occurred during the domestication process as reported previously [3, 48–50]. In addition, there were 23 genomic regions of greater than 20 cM where no polymorphic markers were found among the 11 parents. These regions might indicate fixation (homozygosity) in the modern Upland cotton, and may contain key adaptive genes which have been fixed through breeding and selection processes. The fixed homozygous genomic blocks are a great obstacle to construct a medium to high density intraspecific G. hirsutum linkage map. However, it is also likely that these regions may not contain important agronomic genes or QTLs. Targeted sequencing of these regions may provide better insights about their biological functions. It is worth noting that nineteen of these 23 genomic blocks occurred at the telomeric regions of 17 chromosomes. The biological implications caused by this phenomenon remain unclear, as does their impact on practical breeding. Recently, Gore et al.  observed the same phenomenon when using RILs derived from a cross between TM-1 and NM24016. In our previous research to fine map the Ligon-lintless 1 (Li 1 ) gene, we revealed a similar phenomenon , and so did Cai et al. . Chr.15 had a very big gap (80.47 cM) at the end of its long arm. After examining the consensus map , we noticed that the region between 150 and 230 cM was almost exclusively mapped with RFLP markers. Likewise, the region between 0 and 20 cM of Chr.15 was also composed of RFLP markers in the consensus map. It is very likely that the overall genetic distance of Chr.15 may be not as long as 230 cM as reported by Blenda et al. . Our marker coverage may be greater than 55%. In fact, our markers covered the range between 43,364 bp and 62,598,467 bp of the G. raimondii Chr.02 (corresponding to Chr.01 or Chr.15 of G. hirsutum) (Additional file 2), which is about 99.7% of the entire G. raimondii Chr.02 [1, 44]. A similar situation also exists for Chr.06, 23, 24 and 25.
Fiber QTL numbers, chromosomal distribution and clustering
It is difficult to compare fiber QTL numbers between any two reports in the literature, because many factors affect QTL analysis results. These include experimental populations, population structure, environments, LOD threshold or significance level, and etc.. For example, Lacape et al.  identified 651 QTLs when using LOD 2. However, this number was drastically reduced to 167 when LOD 3.5 was used. In the present research, we used p < 0.01 as a threshold to declare that a marker was associated to a trait. Using such high significance threshold would without a doubt lead to some QTLs, especially those with minor effects not being identified. In fact, the combined phenotypic variance explained by the QTLs reported herein was less than 90% for all six traits (Table 6) although some of the phenotypic variance might be due to other reasons such as genotype × environment effects or epistasis. Our objective was to identify stable QTLs with moderate to major effects, and reduce the QTL number to a more manageable level from a breeder’s perspective.
The 131 fiber QTLs identified in this research were distributed on all 26 chromosomes. However, the distribution was uneven with Chr.26 harboring the greatest number of 11 QTLs. The same result was also reported by Rong et al. . Chr. 26 is known to be rich in fiber genes , and harbors the n 2 fiberless gene . Previously, Rong et al.  reported that Dt subgenome chromosomes contained 25% more fiber QTLs than At subgenome chromosomes. A recent meta-analysis of 810 fiber QTLs  indicated that 23% more fiber QTLs resided on Dt subgenome chromosomes. Our research obtained the same conclusion showing 21% more fiber QTLs residing on the Dt subgenome chromosomes. Both Rong et al.  and Said et al.  reports were largely based on interspecific populations while the present study focused on an intraspecific population. Without a doubt, Dt subgenome chromosomes play more important roles in determining fiber quality.
Congruence with previously reported fiber QTLs
As of today, more than 1000 fiber QTLs have been reported by many researchers. Although it is not easy to compare among the reported QTLs in different studies using different populations, it is possible to compare the QTLs identified in the present research to the prior studies using the common markers and the high density consensus genetic map. We considered two QTLs the same if the associated markers were mapped within 5 cM in the consensus map . Of the 131 fiber QTLs, 77 were reported in prior researches. The detail markers and references are shown in Additional file 6. Some QTLs or clusters were reported by different research groups. Notable examples include qUI-c3, QTL clusters on Chr.05, Chr.07, Chr.16, and Chr.22. The fiber strength QTL, qSTR-c24, on Chr.24 reported by Zhang’s group [2, 25] was also identified in the present study. Most of the prior QTLs that our research agreed with were originally identified in G. hirsutum cultivars , intraspecific populations [17, 55], chromosome substitution inbred population , or interspecific populations but with a few generations of backcross with G. hirsutum[56–58]. QTLs identified in an interspecific population such as those reported by Lacape et al. or Yu et al. [7, 19, 27] agreed poorly with our QTLs. This might be due to the reasons mentioned in the Background section such as QTLs from G. barbadense not being present in Upland cotton. This also presents a challenge using the QTLs identified in an interspecific population in Upland cotton breeding. Fifty-four QTLs reported in this research were not found to be close to any previously-reported fiber QTLs, and could be considered as new fiber QTLs.
Said et al.  identified 36 fiber hotspots. We compared the fiber QTLs identified in the present research with these hotspots to see whether any of these QTLs belong to these hotspots. We used the high density consensus map and common markers as the guidance. If the QTL is within 10 cM range of the hotspot, the QTL is considered as part of the hotspot. Fourteen QTLs including 2, 7, 1, 3, and 1 for ELO, MIC, STR, UHM, and UI, respectively, likely belong to the fiber hotspots described by Said et al. . The detail list of QTLs and hotspots is shown in the Additional file 6.
One of the major challenges that cotton breeders have been facing when using MAS to improve quantitative traits such as fiber quality is the difficulty to transfer a QTL from other species such as G. barbadense into Upland cotton without compromise of other traits. While QTLs identified in an intraspecific population will have fewer obstacles to implement in breeding, the low genetic diversity within Upland cotton presents another challenge for researchers to accurately identify and map stable QTLs with moderate to high effects. In this study, we used a random-mated recombinant inbred population involving 11 diverse Upland cotton cultivars to identify fiber QTLs. This approach ensured higher polymorphism which enabled us to genotype the RIL population with 1582 polymorphic SSR markers, the most as of today. Coupled with a larger population (550 RILs, the largest single population ever reported in Upland cotton), we have identified 131 fiber QTLs and 37 QTL clusters at p < 0.01. Seventy-seven QTLs were mapped at the same or similar positions as previously identified QTLs in the literature, while 54 were new. By selecting favorable alleles of only 7 QTLs (markers), the mean trait value of the selected RILs could be significantly improved (higher for ELO, STR, UHM or UI, and lower for MIC or SFC). The QTL clusters on Chr.07 and Chr.16 had major effects on 4 fiber traits. Using these two loci could simultaneously improve 4 fiber traits. We are continuing to evaluate these QTLs across a panel of cotton varieties. In addition, we are testing these RILs in multiple locations to obtain yield data, and to assess the effects of random-mating on breaking the negative linkage between yield and fiber quality traits.
Random-mated recombinant inbred population development
Eleven Upland cotton (G. hirsutum L.) lines (10 cultivars and 1 breeding line) (Table 1) that represent a diverse group of non-related cotton lines from major breeding programs across the United States were chosen as parents for random mating so that the populations developed and QTLs identified should have more applicability to the entire US Cotton Belt. Pedigrees for all except M240RNR can be found in . The 11 selected parents were crossed in a half diallel to produce 55 half-sib families in 2002 in Mississippi State, MS, USA. Using 55 F1 crosses as the founding populations, pollens from all crosses were collected, mixed, and pollinated to each founding population . Random mating of the F1 from the half diallel was designated Cycle 0 (C0). Five more cycles of random-mating were made in both Mississippi State, MS, USA and the Cotton Winter Nursery in Tecoman, Mexico using the bulked pollen method . These 55 families were kept separate during the random mating process. The detail of producing the random mated population was previously reported . The original germplasm derived from random-mating was self-pollinated one time to increase the amount of seed for distribution and was released in 2008 under the name RMUP-C5 (Random Mated Upland Population Cycle 5) . After six cycles of random mating, self-pollination was followed for six generations using single seed descent. Ten lines were randomly selected from each of these 55 founding populations and a new population including 550 RILs was created. These 550 RILs (C5S6) were used in the present research.
Field planting and fiber quality measurement
Seeds of 550 RILs along with their 11 parents were planted as three replicates in a randomized complete block on the Plant Science Research Farm at Mississippi State, MS, in 2010--2011. Each plot was 12 m long with about 120 plants. Standard field practices were applied during the plant growing seasons. Twenty-five health-looking naturally-open bolls from the central part of a plant were hand harvested from each RIL and parent in both years. Boll samples were ginned on a 10-saw laboratory gin, and fiber samples were used for fiber property analyses. Elongation (ELO,%), micronaire (MIC), short fiber content (SFC,%), bundle strength (STR, g/tex), upper-half mean fiber length (UHM, mm), and uniformity index (UI,%) were measured by Cotton Incorporated’s fiber measurement laboratory using a High Volume Instrument (HVI, USTER Technologies Inc., Charlotte, NC). Refer to Said et al.  for the description of each fiber quality attribute.
DNA extraction and SSR marker analysis
Young leaves were collected from 15 plants of each RIL, and bulked. Leaves were freeze-dried using a lyophilizer (The Virtis Company, Inc., Gardiner, NY) and crushed to powder. Total DNA was extracted from lyophilized leaf powder as previously described .
In order to identify polymorphic markers, we first screened 15538 SSR primer pairs using the DNAs of the 11 parents. These primer pairs were randomly selected, and their sequences were obtained from the Cotton Marker Database (http://www.cottonmarker.org). Two thousand one hundred thirty two (2132) markers were polymorphic among the 11 parents. Second, we identified the positions of these markers according to the high density consensus genetic map . In addition, we used the BLASTN version 2.2.26+ algorithm to align context nucleotide sequences for 2007 of the 2132 SSR markers to the D5G. raimondii Ulbr. reference genome sequence  with an E-value cutoff of 1e−20. The clone sequences for the remainder 125 markers were not available in the public databases. Third, we selected 1582 polymorphic markers that covered as much tetraploid cotton genome as possible and were uniformly distributed across the genome. The marker selection was based on both high density genetic map  and sequence alignment against D5G. raimondii reference genome sequence. We also included 13 markers that were neither mapped nor had clone sequences in the public database (Additional file 2) to see whether these non-mapped markers would be as useful as mapped-markers. And fourth, we genotyped 275 RILs of Set 1 using the selected 1582 markers. These 275 RILs of Set 1 consisted of 5 randomly selected lines from each of the 55 founding populations. The other 275 RILs that were not genotyped with full 1582 markers were referring to as Set 2.
For each marker primer pair, forward primer was fluorescently-labeled at 5’ end with 6-FAM (6-carboxyfluorescein), HEX (4,7,2′,4′,5,7-hexachloro-carboxyfluorescein) or NED (7′,8’-benzo-5-fluoro 2′,4,7,-trichloro-5-carboxyfluorescein). SSR primer oligos were purchased from Sigma Genosys (Woodlands, Texas) or Life Technologies (Foster City, CA). Triplex PCR (three pairs of primers with different dyes) was performed when conducting primer screening. Hexaplex (two markers for each dye) PCR was conducted when genotyping the RILs. PCR products were separated using an ABI Genetic Analyzer 3730 xl with ROX Gene-Scan 500® as internal size standard. The PCR amplification conditions and marker data acquisition were previously described . If a marker revealed two loci within the population, the duplicate marker loci were designated by adding a lower-case letter in alphabetical order after the primer name.
Genetic diversity statistics
A pairwise matrix of genetic similarity values was calculated in NTSYSpc2.2  using DICE coefficient . A dendrogram tree was constructed based on Unweighted Pair Group Method with Arithmatic Mean (UPGMA) using the NTSYSpc2.2 software package.
RIL population structure assessment
The efficacy of random-mating was estimated by evaluating the existence of potential population structure of the RILs using software packages STRUCTURE 2.3  and JMP Genomics 6.0 (SAS Institute, Cary, NC). The software STRUCTURE 2.3 was implemented for a model-based clustering method for inferring population structure using 1582 polymorphic SSR marker data. Ten thousand run length was applied for the analysis to generate Q matrix, which was used for association analysis as well when using software program TASSEL . When using JMP Genomics 6.0, the relatedness analysis was implemented. The relationship matrix process computes a symmetric matrix of pairwise relatedness measures for RIL entries across all SSR marker loci. Identity By Descent (IBD) was used in the study to estimate the genome-wide relatedness per JMP Genomics. The RIL population structure assessment was conducted in two ways: 1) using 275 RILs of Set 1 based on 1582 markers; 2) using 550 RILs based on 270 SSR markers.
Phenotypic data analysis and QTL association mapping
One objective of phenotypic data analysis is to obtain precise estimate of the entry’s genetic value, where the phenotypic data were generated from the field test(s) which is confounded by environmental factors. The six fiber traits were initially screened for outliers using SAS version 9.3 software package  (SAS Institute, Cary, NC) by examining the Studentized deleted residuals  obtained from mixed linear models fitted with environment, line, and replication nested with environment as random effects. For each trait, a best linear unbiased predictor (BLUP) for each line was predicted from a mixed linear model fitted across environments with ASReml version 3.0 . The statistical model is:
Where y hij is an observed value, μ is the population mean, E n is an environmental effect, G i is an genotypic effect, GE hi is a genotype-by-environment interaction effect, Bj(h) is a block effect, and e hij is a random error. In this study, we treated all effects as random except population mean and environmental effects. This model was used to estimate variance components (including heritability) and to predict genotypic value with the following equation: . The predicted genotypic values were further used to (1) estimate the Pearson correlations among these fiber traits for parents and RIL lines, respectively, and (2) conduct association mapping and results validation as detailed in the following paragraphs. Pearson correlation analysis was conducted by SAS (SAS Institute, Cary, NC).
The association analysis for fiber QTLs was conducted using the software TASSEL 3.0  by following the default protocols. We first ran a combined data analysis by a linear mixed model approach for variance component estimation and a jackknife approach for statistical test. Our results showed that GxE interaction effects were present; however, the contributions of GxE effects to these fiber traits were numerically small (ranged from 3 ~ 5%, data not shown). Thus, we used the means from the two tests for our association mapping analysis. Association analysis was first conducted using general linear model. After removing monomorphic loci and minor allele with minor allele frequency (MAF) threshold of 0.05, the Q matrix of population structure with number of estimated populations was added. Then, association analysis by using mixed linear model was implemented. In order to reduce the possibility of false QTLs, a marker locus was considered as associated with a trait only when p value was smaller than 0.01. A significance level was determined as following: highly significant (HS) p < 0.01, and extremely significant (ES) p < 0.002. When multiple marker loci that were mapped within a 5 cM interval were associated with a trait, they were considered as a single QTL. For the non-mapped markers, we used their physical alignments on the G. raimondii genome sequences to assign to a QTL. This assignment was tentative due to the difficulty to separate two homeologous chromosomes of tetraploid cotton based on only G. raimondii sequence. Nine non-mapped markers that were associated to a trait were not assigned to a QTL due to the uncertainty of their genetic or physical locations. In the present report, they are not included in the QTL counts, but will be included when more information becomes available. The QTL nomenclature was according to .
Confirmation of fiber QTLs
Although the 275 RILs of Set 2 were not genotyped with full set of 1582 SSR markers, they were planted in the field, and phenotypic data were obtained. Based on the QTL analysis results obtained using 275 RILs of Set 1 described above, we selected a total of 117 markers to genotype the 275 RILs of Set 2 for the purpose to confirm fiber QTLs. We also genotyped the Set 2 RILs using an additional 153 SSR markers in order to assess the structure of the whole population. For each trait, we selected 7 markers (QTLs) to group the Set 2 RILs into two groups: RILs with 3 or more favorable alleles, and the remaining RILs. We also did the same using all 550 RILs.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its seven Additional files. In addition, the data about the fiber QTLs and their genomic locations reported in this article are also available at CottonGen database (http://www.cottongen.org).
Marker assisted selection
Quantitative trait locus (loci)
Recombinant inbred lines
Short fiber content
Simple sequence repeat
Fiber bundle strength
Upper-half mean fiber length
This research was funded by United States Department of Agriculture-Agricultural Research Service CRIS project 6435-21000-017-00D and Cotton Incorporated project 10–747. We thank Dr. Russell Hayes for assisting with the field experiments, Dr. Gregory Thyssen for his assistance to align marker sequences against the G. raimondii reference genome. We also thank Drs. Kater Hake and Don Jones from Cotton Incorporated for their support from the onset of this project. Our appreciation goes to Mrs. Sheron Simpson and Dr. Brian Scheffler at the Genomics and Bioinformatics Research Unit at Stoneville, MS for their excellent support in SSR marker analysis. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U. S. Department of Agriculture which is an equal opportunity provider and employer.
- Paterson AH, Wendel JF, Gundlach H, Guo H, Jenkins J, Jin D, Llewellyn D, Showmaker KC, Shu S, Udall J, Yoo MJ, Byers R, Chen W, Doron-Faigenboim A, Duke MV, Gong L, Grimwood J, Grover C, Grupp K, Hu G, Lee TH, Li J, Lin L, Liu T, Marler BS, Page JT, Roberts AW, Romanel E, Sanders WS, Szadkowski E, et al: Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature. 2012, 492 (7429): 423-427. 10.1038/nature11798.PubMedView ArticleGoogle Scholar
- Zhang T, Qian N, Zhu X, Chen H, Wang S, Mei H, Zhang Y: Variations and transmission of QTL alleles for yield and fiber qualities in upland cotton cultivars developed in China. PLoS One. 2013, 8 (2): e57220-10.1371/journal.pone.0057220.PubMed CentralPubMedView ArticleGoogle Scholar
- Wendel JF, Cronn RC: Polyploidy and the evolutionary history of cotton. Adv Agron. 2003, 78: 139-186.View ArticleGoogle Scholar
- Fryxell PA: A revised taxonomic interpretation of Gossypium L. (Malvaceae). Rheedea. 1992, 2: 108-165.Google Scholar
- Brubaker CL, Bourland FM, Wendel JF: The origin and domestication of cotton. Cotton: Origin, history, technology, and production. Edited by: Smith CW, Cothren JT. 1999, New York: Wiley & Sons, 3-32.Google Scholar
- Meredith WR: Registration of MD 52ne high fiber quality cotton germplasm and recurrent parent MD 90ne. Crop Sci. 2005, 45: 806-807. 10.2135/cropsci2005.0806.View ArticleGoogle Scholar
- Yu J, Zhang K, Li S, Yu S, Zhai H, Wu M, Li X, Fan S, Song M, Yang D, Li Y, Zhang J: 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. 2013, 126 (1): 275-287. 10.1007/s00122-012-1980-x.PubMedView ArticleGoogle Scholar
- Meredith MR: Quantitative genetics. Cotton. Edited by: Koh JW, Lewis CF. 1984, Madison, WI, USA: Crop Science Sciety of America, 132-147.Google Scholar
- Humphrey AB, Matzinger DF, Cockerham CC: Effects of random intercrossing in a naturally self-fertilizing species, Nicotiana tabacum L. Crop Sci. 1969, 9: 495-497. 10.2135/cropsci1969.0011183X000900040033x.View ArticleGoogle Scholar
- Nordquist PT, Webster OJ, Gardner CO, Ross VW: Registration of three sorghum germplasm random mating populations. Crop Sci. 1973, 13: 132-10.2135/cropsci1973.0011183X001300010056x.View ArticleGoogle Scholar
- Burton JW, Brim CA: Registration of two soybean germplasm populations. Crop Sci. 1981, 21: 801-10.2135/cropsci1981.0011183X002100050052x.View ArticleGoogle Scholar
- Frey KY, Holland JB: Nine cycles of recurrent selection for increased groat-oil content in oat. Crop Sci. 1999, 39: 1636-1641. 10.2135/cropsci1999.3961636x.View ArticleGoogle Scholar
- Jenkins JN, McCarty JC, Gutierrez OA, Hayes RW, Bowman DT, Watson CE, Jones DC: Registration of RMUP-C5, a random mated population of upland cotton germplasm. J Plant Regist. 2008, 2: 239-242. 10.3198/jpr2008.02.0080crg.View ArticleGoogle Scholar
- Meredith WR, Bridge RR: Breakup of linkage blocks in cotton. Crop Sci. 1971, 11: 695-697. 10.2135/cropsci1971.0011183X001100050027x.View ArticleGoogle Scholar
- Miller PA, Rawlings JO: Breakup of linkage blocks through intermating in a cotton breeding population. Crop Sci. 1967, 7: 199-204. 10.2135/cropsci1967.0011183X000700030009x.View ArticleGoogle Scholar
- Green CC, Culp TW: Simultaneous improvement of yield, fiber quality, and yarn stregth in Upland cotton. Crop Sci. 1990, 30: 66-69. 10.2135/cropsci1990.0011183X003000010015x.View ArticleGoogle Scholar
- Zhang K, Zhang J, Ma J, Tang S, Liu D, Teng Z, Liu D, Zhang Z: 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. 2012, 29 (2): 335-348. 10.1007/s11032-011-9549-y.View ArticleGoogle Scholar
- Lacape JM, Gawrysiak G, Cao TV, Viot C, Llewellyn D, Liu S, Jacobs J, Becker D, Vianna Barroso PA, Assunção JH, Palaï O, Georges S, Jean J, Gibanda M: Mapping QTLs for traits related to phenology, morphology and yield components in an inter-specific Gossypium hirsutum × G. barbadense cotton RIL population. Field Crops Res. 2013, 144: 256-267.View ArticleGoogle Scholar
- Lacape JM, Llewellyn D, Jacobs J, Arioli T, Becker D, Calhoun S, Al-Ghazi Y, Liu S, Palai O, Georges S, Giband M, de Assuncao H, Barroso PA, Claverie M, Gawryziak G, Jean J, Vialle M, Viot C: Meta-analysis of cotton fiber quality QTLs across diverse environments in a Gossypium hirsutum x G. barbadense RIL population. BMC Plant Biol. 2010, 10: 132-10.1186/1471-2229-10-132.PubMed CentralPubMedView ArticleGoogle Scholar
- Paterson AH, Saranga Y, Menz M, Jiang CX, Wright RJ: QTL analysis of genotype x environment interactions affecting cotton fiber quality. Theor Appl Genet. 2003, 106 (3): 384-396.PubMedGoogle Scholar
- Rong J, Feltus FA, Waghmare VN, Pierce GJ, Chee PW, Draye X, Saranga Y, Wright RJ, Wilkins TA, May OL, Smith CW, Gannaway JR, Wendel JF, Paterson AH: Meta-analysis of polyploid cotton QTL shows unequal contributions of subgenomes to a complex network of genes and gene clusters implicated in lint fiber development. Genetics. 2007, 176 (4): 2577-2588. 10.1534/genetics.107.074518.PubMed CentralPubMedView ArticleGoogle Scholar
- Fang DD, Xiao J, Canci PC, Cantrell RG: A new SNP haplotype associated with blue disease resistance gene in cotton (Gossypium hirsutum L.). Theor Appl Genet. 2010, 120 (5): 943-953. 10.1007/s00122-009-1223-y.PubMedView ArticleGoogle Scholar
- Rong J, Pierce GJ, Waghmare VN, Rogers CJ, Desai A, Chee PW, May OL, Gannaway JR, Wendel JF, Wilkins TA, Paterson AH: Genetic mapping and comparative analysis of seven mutants related to seed fiber development in cotton. Theor Appl Genet. 2005, 111 (6): 1137-1146. 10.1007/s00122-005-0041-0.PubMedView ArticleGoogle Scholar
- Said JI, Lin Z, Zhang X, Song M, Zhang J: A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton. BMC Genomics. 2013, 14: 776-10.1186/1471-2164-14-776.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen H, Qian N, Guo W, Song Q, Li B, Deng F, Dong C, Zhang T: Using three overlapped RILs to dissect genetically clustered QTL for fiber strength on Chro.D8 in Upland cotton. Theor Appl Genet. 2009, 119 (4): 605-612. 10.1007/s00122-009-1070-x.PubMedView ArticleGoogle Scholar
- Kumar P, Singh R, Lubbers EL, Shen X, Paterson AH, Campbell BT, Jones DC, Chee PW: Mapping and validation of fiber strength quantitative trait loci on chromosome 24 in upland cotton. Crop Sci. 2012, 52: 1115-1122. 10.2135/cropsci2011.09.0524.View ArticleGoogle Scholar
- Lacape JM, Nguyen TB, Courtois B, Belot JL, Giband M, Gourlot JP, Gawryziak G, Roques S, Hau B: QTL analysis of cotton fiber quality using multiple G. hirsutum × G. barbadense backcross generations. Crop Sci. 2005, 45: 123-140. 10.2135/cropsci2005.0123a.View ArticleGoogle Scholar
- Mei M, Syed NH, Gao W, Thaxton PM, Smith CW, Stelly DM, Chen ZJ: Genetic mapping and QTL analysis of fiber-related traits in cotton (Gossypium). Theor Appl Genet. 2004, 108 (2): 280-291. 10.1007/s00122-003-1433-7.PubMedView ArticleGoogle Scholar
- Qin H, Guo W, Zhang YM, Zhang T: QTL mapping of yield and fiber traits based on a four-way cross population in Gossypium hirsutum L. Theor Appl Genet. 2008, 117 (6): 883-894. 10.1007/s00122-008-0828-x.PubMedView ArticleGoogle Scholar
- Wang B, Guo W, Zhu X, Wu Y, Huang N, Zhang T: QTL mapping of fiber quality in an elite hybrid derived-RIL population of Upland cotton. Euphytica. 2006, 152: 367-378. 10.1007/s10681-006-9224-2.View ArticleGoogle Scholar
- Wang B, Guo W, Zhu X, Wu Y, Huang N, Zhang T: QTL mapping of yield and yield components for elite hybrid derived-RILs in Upland cotton. J Genet Genomics. 2007, 34 (1): 35-45. 10.1016/S1673-8527(07)60005-8.PubMedView ArticleGoogle Scholar
- Wang P, Zhu Y, Song X, Cao Z, Ding Y, Liu B, Zhu X, Wang S, Guo W, Zhang T: Inheritance of long staple fiber quality traits of Gossypium barbadense in G. hirsutum background using CSILs. Theor Appl Genet. 2012, 124 (8): 1415-1428. 10.1007/s00122-012-1797-7.PubMedView ArticleGoogle Scholar
- Gore MA, Percy RG, Zhang J, Fang DD, Cantrell RG: Registration of the TM-1/NM24016 cotton recombinant inbred mapping population. J Plant Regist. 2012, 6 (1): 124-127. 10.3198/jpr2011.06.0334crmp.View ArticleGoogle Scholar
- Ulloa M, Meredith WR: Genetic linkage map and QTL analysis of agronomic and fiber quality traits in an intraspecific population. J Cotton Sci. 2000, 4: 161-170.Google Scholar
- Zeng L, Meredith WR, Gutierrez OA, Boykin DL: Identification of associations between SSR markers and fiber traits in an exotic germplasm derived from multiple crosses among Gossypium tetraploid species. Theor Appl Genet. 2009, 119 (1): 93-103. 10.1007/s00122-009-1020-7.PubMedView ArticleGoogle Scholar
- Fang DD, Hinze LL, Percy RG, Li P, Deng D, Thyssen G: A microsatellite-based genome-wide analysis of genetic diversity and linkage disequilibrium in Upland cotton (Gossypium hirsutum L.) cultivars from major cotton-growing countries. Euphytica. 2013, 191: 391-401. 10.1007/s10681-013-0886-2.View ArticleGoogle Scholar
- Lacape JM, Dessauw D, Rajab M, Noyer JL, Hau B: Microsatellite diversity in tetraploid Gossypium germplasm: assembling a highly informative genotyping set of cotton SSRs. Mol Breed. 2007, 19: 45-58.View ArticleGoogle Scholar
- Hinze LL, Dever JK, Percy RG: Molecular variation among and within improved cultivars in the US cotton germplasm collection. Crop Sci. 2012, 52: 222-230. 10.2135/cropsci2011.04.0202.View ArticleGoogle Scholar
- Abdurakhmonov IY, Saha S, Jenkins JN, Buriev ZT, Shermatov SE, Scheffler BE, Pepper AE, Yu JZ, Kohel RJ, Abdukarimov A: Linkage disequilibrium based association mapping of fiber quality traits in G. hirsutum L. variety germplasm. Genetica. 2009, 136: 401-417. 10.1007/s10709-008-9337-8.PubMedView ArticleGoogle Scholar
- Mei H, Zhu X, Zhang T: Favorable QTL alleles for yield and its components identified by association mapping in Chinese Upland cotton cultivars. PLoS One. 2013, 8: e82193-10.1371/journal.pone.0082193.PubMed CentralPubMedView ArticleGoogle Scholar
- Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007, 23 (19): 2633-2635. 10.1093/bioinformatics/btm308.PubMedView ArticleGoogle Scholar
- Dice LR: Measures of the amount of ecologic association between species. Ecology. 1945, 26: 297-302. 10.2307/1932409.View ArticleGoogle Scholar
- Blenda A, Fang DD, Rami JF, Garsmeur O, Luo F, Lacape JM: A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check. PLoS One. 2012, 7 (9): e45739-10.1371/journal.pone.0045739.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Z, Zhang D, Wang X, Tan X, Guo H, Paterson AH: A whole-genome DNA marker map for cotton based on the D-genome sequence of Gossypium raimondii L. G3 , Genes|Genomes|Genetics. 2013, 3 (10): 1759-1767.PubMed CentralPubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMed CentralPubMedGoogle Scholar
- Gore MA, Fang DD, Poland JA, Zhang J, Percy RG, Cantrell RG, Thyssen G: Linkage map construction and QTL analysis of agronomic and fiber quality traits in an introgressed recombinant inbred population of cotton (Gossypium hirsutum L.). Plant Genome. 2014, 7 (1): doi:10.3835/plantgenome2013.3807.0023Google Scholar
- Blenda A, Scheffler J, Scheffler B, Palmer M, Lacape JM, Yu JZ, Jesudurai C, Jung S, Muthukumar S, Yellambalase P, Ficklin S, Staton M, Eshelman R, Ulloa M, Saha S, Burr B, Liu SL, Zhang TZ, Fang DQ, Pepper A, Kumpatla S, Jacobs J, Tomkins J, Cantrell R, Main D: CMD: a Cotton Microsatellite Database resource for Gossypium genomics. BMC Genomics. 2006, 7: 132-10.1186/1471-2164-7-132.PubMed CentralPubMedView ArticleGoogle Scholar
- Brubaker CL, Wendel JF: Re-evaluating the origin of domesticated cotton (Gossypium hirsutum: Malvaceae) using nuclear restriction fragment length polymorphisms (RFLPs). Am J Bot. 1994, 81: 1309-1326. 10.2307/2445407.View ArticleGoogle Scholar
- Iqbal AE, Reddy OUK, El-Zik KM, Pepper AE: A genetic bottleneck in the ’evolution under domestication’ of upland cotton Gossypium hirsutum L. examined using DNA fingerprinting. Theor Appl Genet. 2001, 103: 547-554. 10.1007/PL00002908.View ArticleGoogle Scholar
- May OL, Bowman DT, Calhoun DS: Genetic diversity of US upland cotton cultivars released between 1980 and 1990. Crop Sci. 1995, 35: 1570-1574. 10.2135/cropsci1995.0011183X003500060009x.View ArticleGoogle Scholar
- Gilbert MK, Turley RB, Kim HJ, Li P, Thyssen G, Tang Y, Delhom CD, Naoumkina M, Fang DD: Transcript profiling by microarray and marker analysis of the short cotton (Gossypium hirsutum L.) fiber mutant Ligon lintless-1 (Li 1 ). BMC Genomics. 2013, 14: 403-10.1186/1471-2164-14-403.PubMed CentralPubMedView ArticleGoogle Scholar
- Cai C, Tong X, Liu F, Lv F, Wang H, Zhang T, Guo W: Discovery and identification of a novel Ligon lintless-like mutant (Lix) similar to the Ligon lintless (Li1) in allotetraploid cotton. Theor Appl Genet. 2013, 126 (4): 963-970. 10.1007/s00122-012-2029-x.PubMedView ArticleGoogle Scholar
- Xu Z, Kohel RJ, Song G, Cho J, Yu J, Yu S, Tomkins J, Yu JZ: An integrated genetic and physical map of homoeologous chromosomes 12 and 26 in Upland cotton (G. hirsutum L.). BMC Genomics. 2008, 9: 108-10.1186/1471-2164-9-108.PubMed CentralPubMedView ArticleGoogle Scholar
- An C, Jenkins JN, Wu J, Guo Y, McCartin JA: Use of fiber and fuzz mutants to detect QTL for yield components, seed and fiber traits of upland cotton. Euphytica. 2010, 172: 21-34. 10.1007/s10681-009-0009-2.View ArticleGoogle Scholar
- Shen X, Guo W, Lu Q, Zhu X, Yuan Y, Zhang T: Genetic mapping of quantitative trait loci for fiber quality and yield trait by RIL approach in Upland cotton. Euphytica. 2007, 155 (3): 371-380. 10.1007/s10681-006-9338-6.View ArticleGoogle Scholar
- Draye X, Chee P, Jiang CX, Decanini L, Delmonte TA, Bredhauer R, Smith CW, Paterson AH: Molecular dissection of interspecific variation between Gossypium hirsutum and G. barbadense (cotton) by a backcross-self approach: II. Fiber fineness. Theor Appl Genet. 2005, 111 (4): 764-771. 10.1007/s00122-005-2061-1.PubMedView ArticleGoogle Scholar
- Chee PW, Draye X, Jiang CX, Decanini L, Delmonte TA, Bredhauer R, Smith CW, Paterson AH: Molecular dissection of phenotypic variation between Gossypium hirsutum and Gossypium barbadense (cotton) by a backcross-self approach: III. Fiber length. Theor Appl Genet. 2005, 111 (4): 772-781. 10.1007/s00122-005-2062-0.PubMedView ArticleGoogle Scholar
- Chee P, Draye X, Jiang CX, Decanini L, Delmonte TA, Bredhauer R, Smith CW, Paterson AH: Molecular dissection of interspecific variation between Gossypium hirsutum and Gossypium barbadense (cotton) by a backcross-self approach: I. Fiber elongation. Theor Appl Genet. 2005, 111 (4): 757-763. 10.1007/s00122-005-2063-z.PubMedView ArticleGoogle Scholar
- Bowman DT, Gutierrez OA, Percy RG, Calhoun DS, May OL: Pedigrees of upland and pima cotton cultivars released between 1970 and 2005. 2007, Mississippi State, MS: Mississippi State Univ, Mississippi Agricultural and Forestry Experimental Station Bulletin # 1155: http://msucares.com/pubs/bulletins/b1155.pdfGoogle Scholar
- Gutiérrez OA, Bowman DT, Cole CB, Jenkins JN, McCarty JC, Wu J, Watson CE: Development of random-mated populations using bulked pollen methodology: cotton as a model. J Cotton Sci. 2006, 10: 175-179.Google Scholar
- Rohlf FJ: NTSYSpc: Numerical taxonomy and multivariate analysis system, version 2.2. 2010, Setauket, NY: Exeter SoftwareGoogle Scholar
- SAS Institute: The SAS system for Windows. Release 9.3. 2012, Cary, NC, USA: SAS InstituteGoogle Scholar
- Kutner MH, Nachtsheim CJ, Neter J, Li W: Applied linear statistical models. 2004, Boston, MA: McGraw-Hill, 4Google Scholar
- Gilmour AR, Gogel B, Cullis B, Thompson R, Butler D, Cherry M, Collins D, Dutkowski D, Harding S, Haskard K: ASReml user guide release 3.0. 2009, Hemel Hempstead, UK: VSN International LtdGoogle Scholar
- McCouch SR, Cho YG, Yano PE, Blinstrub M, Morishima H, Kinoshita T: Report on QTL nomenclature. Rice Genet Newslet. 1997, 14: 11-13.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.