ArticlePDF Available

QTL mapping and genetic effect of chromosome segment substitution lines with excellent fiber quality from Gossypium hirsutum × Gossypium barbadense

Authors:
  • Institute of Cotton Research, Chinese Academy of Agricultural Sciences

Abstract and Figures

Chromosome segment substitution lines (CSSLs) are ideal materials for identifying genetic effects. In this study, CSSL MBI7561 with excellent fiber quality that was selected from BC4F3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC5F2 and BC5F2:3). Eighty-one polymorphic markers related to 33 chromosome introgressive segments on 18 chromosomes were finally screened using 2292 SSR markers which covered the whole tetraploid cotton genome. A total of 129 quantitative trait loci (QTL) associated with fiber quality (103) and yield-related traits (26) were detected on 17 chromosomes, explaining 0.85–30.35% of the phenotypic variation; 39 were stable (30.2%), 53 were common (41.1%), 76 were new (58.9%), and 86 had favorable effects on the related traits. More QTL were distributed in the Dt subgenome than in the At subgenome. Twenty-five stable QTL clusters (with stable or common QTL) were detected on 22 chromosome introgressed segments. Finally, the 6 important chromosome introgressed segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3, Seg-D07-3, and Seg-D06-2) were identified as candidate chromosome regions for fiber quality, which should be given more attention in future QTL fine mapping, gene cloning, and marker-assisted selection (MAS) breeding.
This content is subject to copyright. Terms and conditions apply.
Vol.:(0123456789)
1 3
Molecular Genetics and Genomics (2019) 294:1123–1136
https://doi.org/10.1007/s00438-019-01566-8
ORIGINAL ARTICLE
QTL mapping andgenetic eect ofchromosome segment substitution
lines withexcellent ber quality fromGossypium hirsutum ×
Gossypium barbadense
Shao‑qiLi1· Ai‑yingLiu1· Ling‑leiKong1· Ju‑wuGong1· Jun‑wenLi1· Wan‑kuiGong1· Quan‑weiLu2· Peng‑taoLi2·
QunGe1· Hai‑hongShang1· Xiang‑huiXiao1· Rui‑xianLiu1· QiZhang1· Yu‑zhenShi1· You‑luYuan1
Received: 6 November 2018 / Accepted: 3 April 2019 / Published online: 27 April 2019
© The Author(s) 2019
Abstract
Chromosome segment substitution lines (CSSLs) are ideal materials for identifying genetic effects. In this study, CSSL
MBI7561 with excellent fiber quality that was selected from BC4F3:5 of CCRI45 (Gossypium hirsutum) × Hai1 (Gossypium
barbadense) was used to construct 3 secondary segregating populations with 2 generations (BC5F2 and BC5F2:3). Eighty-one
polymorphic markers related to 33 chromosome introgressive segments on 18 chromosomes were finally screened using
2292 SSR markers which covered the whole tetraploid cotton genome. A total of 129 quantitative trait loci (QTL) associ-
ated with fiber quality (103) and yield-related traits (26) were detected on 17 chromosomes, explaining 0.85–30.35% of the
phenotypic variation; 39 were stable (30.2%), 53 were common (41.1%), 76 were new (58.9%), and 86 had favorable effects
on the related traits. More QTL were distributed in the Dt subgenome than in the At subgenome. Twenty-five stable QTL
clusters (with stable or common QTL) were detected on 22 chromosome introgressed segments. Finally, the 6 important
chromosome introgressed segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3, Seg-D07-3, and Seg-D06-2) were identi-
fied as candidate chromosome regions for fiber quality, which should be given more attention in future QTL fine mapping,
gene cloning, and marker-assisted selection (MAS) breeding.
Keywords CSSLs· Chromosome introgressed segments· Fiber quality· QTL· Genetic effects
Introduction
As an allopolyploid cash crop, cotton is important for
genetic research and provides the textile industry with the
most important natural fiber raw material. Among cultivated
cotton species, two tetraploid cottons, Upland cotton (Gos-
sypium hirsutum L./G.h) with higher yield and Sea-island
cotton (Gossypium barbadense L./G.b) with better fiber
quality are the most widely cultivated in agricultural pro-
duction (Ulloa etal. 2005; Zhang etal. 2009). Therefore,
it is interesting to introgress the favorable genes for fiber
quality from G.b to G.h to improve the fiber quality and yield
simultaneously. However, introgression is very difficult for
breeders to implement using conventional breeding, because
all of the related traits are quantitative traits controlled by
multiple genetic loci, and the fiber quality and yield-related
traits are usually negatively correlated (Clement etal. 2012;
Ma etal. 2014; Yu etal. 2016). Fortunately, with the rapid
development of high-precision molecular marker technology
and gene mapping, an increasing number of genetic maps
Communicated by S. Hohmann.
Shao-qi Li and Ai-ying Liu have contributed equally to this work
Electronic supplementary material The online version of this
article (https ://doi.org/10.1007/s0043 8-019-01566 -8) contains
supplementary material, which is available to authorized users.
* Yu-zhen Shi
shiyzmb@126.com
* You-lu Yuan
yuanyoulu@caas.cn
1 State Key Laboratory ofCotton Biology, Key Laboratory
ofBiologiacl andGenetic Breeding ofCotton, The Ministry
ofAgriculture, Institute ofCotton Research, Chinese
Academy ofAgricultural Science, Anyang455000, Henan,
China
2 School ofBiotechnology andFood Engineering, Anyang
Institute ofTechnology, Anyang455000, Henan, China
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1124 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
and QTL have been identified (Jia etal. 2016; Kushanov
etal. 2016).
After the first molecular genetics map of cotton was con-
structed (Reinisch etal. 1994), a wide variety of populations
were used to perform QTL mapping, which predominantly
consists of F2 (Brown etal. 2005; Wang etal. 2010; Yu
etal. 2013), double haploid (DH) (Lu etal. 2014; Cai etal.
2015), backcross (BC) (Chen etal. 2012; Shang etal. 2016;
Wang etal. 2016b), backcross inbred lines (BILs) (Pang
etal. 2012; Nie etal. 2015; Zhang etal. 2015a), recombi-
nant inbred lines (RILs) (Zhang etal. 2015c; Jamshed etal.
2016; Shang etal. 2016) and CSSLs (Liang etal. 2010; Li
etal. 2016; Shi etal. 2016; Zhai etal. 2016). The complex
genetic background in most populations makes it difficult to
estimate the positions and effects of QTL. There are a few
differences between CSSLs and recurrent parents, which is
favorable for QTL mapping, genetic effect identification and
gene cloning (Wang etal. 2016c; Lu etal. 2017). Therefore,
more attention has been paid to the development and utiliza-
tion of CSSLs in genetic research, although CSSLs are time
consuming and costly to construct (Wan etal. 2008; Zhao
etal. 2009).
The first CSSLs were constructed in tomatoes by Eshed
and Zamir (1994). Subsequently, researchers launched a
wide range of studies and applications of CSSLs in rice
(Wan etal. 2004), corn (Li etal. 2014), wheat (Liu etal.
2006) and other crops. The first set of cotton CSSLs was
developed with TM-1 as the recipient parent by Stelly etal.
(2005), and was used to analyze genetic effects, as well as
genetic relationships for fiber quality and yield component
traits with substitution of one chromosome (Stelly etal.
2005). After this study, a series of studies were carried out
on cotton CSSLs (Wang etal. 2008, 2016c Chen etal. 2012;
Su etal. 2014; Lu etal. 2017). Yang etal. (2009) detected
51 QTLs using 116 CSSLs originating from CCRI45 (G.h)
and Hai1 (G.b). Wang etal. (2012) indicated the inheritance
of long staple fiber qualities using the CSSLs developed by
TM-1 and Hai7124. Fu etal. (2013) detected 12 QTLs asso-
ciated with fiber quality and yield using the CSSLs from
TM-1 and Sub18.
Although CSSLs are effective in QTL mapping, there
is less information for detecting genetic effects from
introgressing chromosome segments from Island cotton
into Upland cotton. A series of CSSLs were constructed
through the hybridization of CCRI45 (G.h), CCRI36 (G.h)
and Hai1 (G.b) by our team(Shi etal. 2008). Subsequently,
a high-density genetic linkage map was constructed that
contained 2292 SSR markers and covered 5115.16 cM of
the cotton AD genome (Shi etal. 2015), and many QTL
were identified using populations with various generations
(BC4F1, BC4F3, BC4F3:5 and BC4F4) (Yang etal. 2009; Ge
etal. 2012; Ma 2014; Guo etal. 2015; Lan etal. 2015). The
genetic effects and heterosis of yield and yield component
traits were analyzed through hybridizing 10 CSSLs accord-
ing to North Carolina Design II (Li etal. 2016). A total
of 70 QTL and their genetic effects for fiber yield-related
traits and 29 QTL for fiber quality traits on 13 chromosomes
were detected using CSSLs (BC5F3, BC5F3:4 and BC5F3:5)
(He 2014). Twenty two QTL associated with fiber quality
and yield traits on seven chromosomes were detected in F2
and F2:3 with two CSSLs of MBI9749 and MBI9915 as par-
ents (Guo etal. 2018). A total of 24 QTLs for fiber qual-
ity and 11 QTLs for yield traits were detected in the three
segregating generations of double-crossed F1 and F2 and
F2:3, which were constructed using four CSSLs as parents
(MBI9804×MBI9855) × (MBI9752×MBI9134) (Zhai etal.
2016). Eighteen QTL for fiber quality and 6 QTL for yield-
related traits across 7 chromosomes were detected using
BC6F2 and BC6F2:3 with two parents of CCRI36 (recurrent
parent) and MBI9915(CSSL) (Song etal. 2017).
In the present study, BC5F2 and BC5F2:3 populations were
constructed by hybridization of CCRI45 (recurrent parent)
and MBI7561 (BC4F3:5) with excellent and stable fiber qual-
ity, to evaluate the genetic effects of the introgressed seg-
ments by SSR markers. This study is expected to lay the
foundation for future studies, such as genetic mechanism
exploring, QTL fine mapping, gene cloning and MAS breed-
ing applications.
Materials andmethods
Plant materials andpopulation development
MBI7561 as the female parent was selected from the CSSLs
BC4F3:5 which was constructed by advanced backcross and
selfing of combination of CCRI45 (G.h) and Hai1 (G.b).
The recurrent parent CCRI45 was a glandular cotton cultivar
widely grown with high yield and resistance to budworm
(Ma 2014), which was developed by the Institute of Cotton
Research of Chinese Academy of Agricultural Sciences (Shi
etal. 2008, 2015). The donor parent Hai1 was a dominant
glandless G. barbadense with excellent fiber quality. The
fiber quality and yield component traits of MBI7561 were
excellent and stable (Table4).
We constructed F1 (BC5F1) by backcrossing CCRI45
(male) and MBI7561 in Anyang in the summer of 2013.
BC5F1 was planted and self-crossing seeds (BC5F2) were
harvested in Hainan province in the winter of 2013. In
2014, a total of 310 BC5F2 (PopE1) individual plants were
developed and fiber and seeds (BC5F2:3) were collected
from individual plants. Both BC5F1 and BC5F2 populations
were planted in the Anyang experimental farm of the Insti-
tute of Cotton Research of Henan Province of China, with
row length of 8 m, row spacing of 0.8 m, and plant spacing
of 0.25 m. In 2015, a total of 253 BC5F2:3 (PopE2) family
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1125Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
lines were planted in one-row plot with a row length of
5 m on the Anyang experimental farm, and another 250
BC5F2:3 (PopE3) family lines were planted in two-narrow-
row plots, with row length of 3 m and plant spacing of 0.12
m on the Alar experimental farm of the Institute of Cotton
Research of Xinjiang Autonomous Region of China. The
plastic-membrane covering technique and a wide/narrow
row spacing pattern were used. The row spacing alterna-
tion was 0.2 m and 0.6 m. In addition, there were 212
common BC5F2:3 lines in the two different environments.
Investigation ofphenotypic traits
Collection ofphenotypic data
In 2014, naturally opened bolls were collected from the
F2 individual plants for phenotype evaluation, including
boll weight (BW), lint percentage (LP), fiber length (FL),
fiber strength (FS), fiber micronaire (FM), fiber uniform-
ity (FU), and fiber elongation (FE). In 2015, 30 naturally
opened bolls were harvested from the plot for F2:3 lines.
BW and LP were calculated using seed cotton weight,
fiber weight, and boll number. The fiber quality traits were
tested with HFT9000 using HVICC international calibra-
tion cotton samples in the Cotton Quality Supervision and
Testing Center of the Ministry of Agriculture of China.
Analysis ofphenotypic traits
The observed phenotypic data were analyzed using 3 soft-
wares. BW, LP, and transgressive rate over the recurrent
parent (TRORP) were calculated using EXCEL 2010. The
SAS statistical software (version 8.1, SAS Institute, Cary
NC) was used to perform the descriptive statistical analy-
sis of phenotypes. The statistical values included mean,
maximum, minimum, standard deviation (SD), skewness,
kurtosis, and coefficient of variation (CV). For ANOVA,
correlation analysis and all significance tests were per-
formed using the SPSS20.0 software (SPSS, Chicago,
IL, USA). All correlation values were Pearson Correla-
tion Coefficients. A higher FM value is not necessarily for
indicative of better fiber quality, and the fineness of the
fiber is evaluated by the FM value grading, (A level (best):
[3.7,4.2]; B level (better): [3.5, 3.6] and[4.3, 4.9]; C level
(bad): [0,3.4] and [5.0, ∞)). The FM values of CCRI45
in all environments were in the C level, thus, we selected
individuals in the A level and B level for the calculation
of transgressive rate values, and the range of FM values of
the selected individuals should be ≥ 3.5 and ≤ 4.9.
Identication ofgenotypes
DNA extraction andmarker detection
The young leaves of F2 individual plants and parents were
sampled in 2014. Genomic DNA was extracted using the
modified CTAB method(Paterson etal. 1993). The products
of PCR amplification were isolated and identified by 8%
non-denaturing vertical polyacrylamide gel electrophoresis.
The DNA segments in the gel were visualized by the silver
staining method (Zhang etal. 2000). In this study, we used
2292 SSR markers from the genetic linkage map with a total
genetic length of 5115.16 cM (Shi etal. 2015) to screen the
recurrent parent of CCRI45, donor parent of Hai1, MBI7561
and F1(MBI7561*CCRI45). The polymorphic markers were
used to identify genotypes of the F2 individual plants. The
sequences of each primer were obtained from the Cotton
Genome Database (www.cotto ngen.org) and synthesized by
Bioethics Engineering (Shanghai) Co., Ltd.
Analysis ofgenotypes
Genotyping analysis of each sample and distribution analysis
of chromosome introgressed segments were performed by
GGT2.0 software (http://www.plant breed ing.wur.nl) (Van
Berloo 2008), including the number, length and positions of
introgressed segments, and the genetic background recovery
rate of each sample. The nomenclature of segments was as
follows: Seg + the serial number of the chromosome (AD)
+ the serial number of the cluster on the chromosome.
Genetic eects analysis
QTL mapping
The QTL IciMappingV4.1 (www.isbre eding .net/softw
are/?type=detai l&id=18) software was used to perform
QTL mapping. A likelihood of odds (LOD) threshold of 2.5
was used to declare significant additive QTL. The resulting
linkage map and QTL were drawn using MapChart2.2 soft-
ware (Voorrips 2002). The QTL nomenclature was: q + the
English abbreviation of trait + the serial number of chromo-
some + the serial number of the QTL on the chromosome
associated with the same trait + (the direction of the additive
genetic effect). For example, qFL-16-2 (+) represents the
second QTL associated with the FL on chromosome 16 with
a positive additive genetic effect from the G. barbadense
introgressing segments.
QTL‑cluster analysis
The QTL cluster analysis was performed by Biomercator 4.2
software (Arcade etal. 2004). QTL were projected on the
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1126 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
genetic map and QTL cluster analysis were performed for
all traits. Four models were thus generated, and each had an
Akaike information criterion (AIC) value. The model with
the lowest AIC value was selected and used for the position
identification of QTL clusters. The nomenclature of QTL
cluster was: Clu + the English abbreviation of trait +the
serial number of chromosome + the serial number of the
cluster on the chromosome associated with the same trait.
Results
Phenotypic performance ofCSSLs populations
The descriptive statistics of phenotypic data for fiber yield
and fiber quality traits, as well as their recurrent parent
CCRI45, are presented in Table1. In the three populations,
the average performance of BW was smaller than that of
CCRI45, with a significant difference in PopE3; the aver-
age performance of FU was higher than that of CCRI45,
with a significant difference in PopE3; and the other traits
were significantly better than those of CCRI45 except for
LP in PopE2. The TRORP for BW was 17.52–32.83% and
53.62–99.26% for other traits. Overall, there were many
transgressive separations for most of traits in the popula-
tions. All traits showed a continuous normal distribution in
three populations (Table1, Fig. S1), which was consistent
with the characteristics of quantitative traits. Significant
positive correlations were found among most traits (FL and
FS/FE/FU, FS and FU/FE, FU and FE/FM/BW, and FM and
FE/LP/BW), whereas significant negative correlations were
found between FS and FM, LP, and FE / FL in most popula-
tions (Table2).
Analysis ofintrogressive segments
A total of 81 pairs of polymorphic markers were screened
between the parents and 33 introgressive segments distrib-
uted on 18 chromosomes (Fig.1). The genetic background
recovery rate of MBI7561 was 95.60%, and the proportion
of homozygous introgressive segments (138.11 cM, 2.7%)
was significantly higher than that of heterozygous seg-
ments (86.96 cM, 1.7%). Four chromosomes [Chr15(D1),
Ch25(D6), Chr07(A7), and Chr16(D7)] had more introgres-
sive segments than others.
The 81 pairs of SSR markers were used to screen the gen-
otype of the BC5F2 population, and 6 pairs of markers did
not show the polymorphism. The average rate of background
Table 1 Phenotypic performance of fiber quality and yield-related traits in three populations
Sta statistic, SE std. error, TRORP transgressive rate over the recurrent parent, Ske skewness, Kur kurtosis, FL fiber length, FS fiber strength, FM
fiber micronaire, FU fiber uniformity, FE fiber elongation, BW boll weight, LP lint percentage
Environment Trait CCRI45 Population
Mean SD Gener. Mean Min. Max. SD Ske. Kur. TRORP. (%) CV (%)
PopE1 BW (g) 5.63 1.02 BC5F25.30 2.48 9.37 0.88 0.36 0.84 32.83 16.53
LP (%) 37.54 3.22 39.90** 25.52 52.69 2.76 −0.28 1.55 66.46 6.91
FL (mm) 29.48 1.19 30.77** 25.84 35.29 1.20 −0.21 0.55 86.39 3.90
FS (cN/tex) 28.02 1.38 31.47** 23.50 39.60 2.13 −0.02 0.33 94.55 6.77
FM (unit) 4.83 0.56 4.38** 2.56 6.09 0.60 −0.37 −0.16 72.11 13.71
FU (%) 83.98 1.52 84.41 78.20 88.00 1.53 −0.67 0.69 65.34 1.81
FE (%) 6.79 0.06 6.84** 6.50 7.10 0.06 −0.22 0.39 94.60 0.94
PopE2 BW (g) 5.38 0.26 BC5F2:3 5.02 3.30 6.57 0.68 −0.01 −0.55 31.77 13.54
LP (%) 36.74 3.34 38.63 32.61 44.13 2.17 −0.14 −0.22 79.51 5.62
FL (mm) 28.57 1.65 30.86** 26.30 34.90 1.73 −0.18 −0.73 88.45 5.62
FS (cN/tex) 28.65 3.39 34.94** 28.00 44.00 2.78 0.22 −0.07 99.26 7.94
FM (unit) 5.07 0.31 4.58** 3.00 5.70 0.47 −0.32 0.24 77.70 10.35
FU (%) 84.35 1.37 85.26 82.10 87.90 1.25 −0.33 −0.39 77.23 1.47
FE (%) 6.80 0.14 6.94** 6.60 7.30 0.13 0.16 −0.11 76.72 1.82
PopE3 BW (g) 5.74 0.16 BC5F2:3 5.36** 3.54 7.23 0.43 0.06 1.45 17.52 8.08
LP (%) 38.20 1.45 40.84** 32.99 48.00 2.15 −0.25 0.21 89.40 5.28
FL (mm) 28.07 0.91 29.80** 26.40 34.10 1.25 0.08 0.06 92.59 4.18
FS (cN/tex) 26.25 1.08 28.96** 24.40 35.00 1.83 0.20 −0.06 94.18 6.31
FM (unit) 5.09 0.37 4.80** 3.40 5.80 0.43 −0.34 −0.02 60.67 8.95
FU (%) 84.30 0.90 85.24** 82.00 88.10 1.16 −0.25 −0.15 78.23 1.36
FE (%) 6.79 0.09 6.88** 6.50 7.20 0.10 −0.13 0.34 91.18 1.48
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1127Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
Table 2 Correlation of
phenotypic traits in the three
populations
FL fiber length, FS fiber strength, FM fiber micronaire, FU fiber uniformity, FE fiber elongation, BW boll
weight, LP lint percentage
*Significant correlation at the 0.05 level (2-tailed), **Significant correlation at the 0.01 level (2-tailed)
BW LP FL FS FM FU
PopE1
LP −0.19**
FL 0.07 −0.36**
FS −0.08 −0.14 0.5**
FM 0.47** 0.04 −0.05 −0.36**
FU 0.29** −0.19** 0.32** 0.20** 0.37**
FE 0.23** −0.19** 0.46** 0.35** 0.21** 0.43**
PopE2
LP 0.33**
FL 0.46** 0.20**
FS 0.01 −0.03 0.61**
FM 0.57** 0.39** 0.11 −0.22**
FU 0.26** 0.18** 0.32** 0.22** 0.24**
FE 0.50** 0.29** 0.79** 0.55** 0.29** 0.37**
PopE3
LP −0.05
FL −0.12 −0.38**
FS −0.13 −0.16* 0.66**
FM 0.33** 0.34** −0.49** −0.38**
FU −0.04 0.18* 0.13 0.25** −0.05
FE 0.03 −0.16* 0.62** 0.50** −0.13 0.24**
Fig. 1 The graphical genotypes of MBI7561. Note A: genetic background (recurrent parent CCRI45); B: homozygous introgressive segments
(donor parent Hai1); H: Heterozygous introgressive segments
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1128 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
recovery in the BC5F2 population was 97.95% and ranged
from 97.3 to 99.2%. The average rate of homozygous intro-
gressive segments was 1.11% and ranged from 0 to 2%. The
average rate of heterozygous introgressive segments was
0.94% and ranged from 0.3 to 2.4%.
QTL mapping
In total, 65 markers on 17 chromosomes were associated
with the QTL of the seven traits, and 48 of these markers
were associated with QTL in multiple populations. Nineteen
markers located on 8 chromosomes existed in the At subge-
nome, and forty-six markers were located on 9 chromosomes
in the Dt subgenome. Based on the method of inclusive com-
posite interval mapping (ICIM), a total of 129 QTL were
identified in 29 introgressive segments of 17 chromosomes
in the three populations (Table3, TableS1, Fig. S2), with
each explaining 0.85% to 30.35% of the phenotypic variation
(PVE). There were 103 and 26 QTL related to the five fiber
quality traits and two yield component traits, respectively.
Forty-one QTL were distributed in the At subgenome, and
88 were distributed in the Dt subgenome. In addition, there
were 107 QTL distributed in 7 pairs of homologous chro-
mosomes. Forty-five QTL (35% of the total number) were
detected in multiple environments, and 39 of them were
stable. Eighty-six QTL showed positive additive effects, 36
showed negative additive effects, and seven showed unstable
additive effects.
Fiber length
There were 21 QTL for FL on 11 chromosomes with the
PVE ranging from 1.67% to 11.93%; 7 and 14 of these
QTL were distributed in the At- and Dt subgenomes,
respectively. Chr16 and Chr25 were the top 2 chromo-
somes with the largest number of QTL. Sixteen QTL with
additive effects from 0.06 to 0.75 mm indicated that Hai1
alleles increased FL, and four QTL with additive effects
from −0.60 to −0.22 mm indicated that CCRI45 alleles
increased FL. The qFL-16-3 had unstable genetic effects,
which could be detected in two environments, with addi-
tive effect of −0.23 mm in PopE2 and 0.58 mm in PopE3.
Three QTL (qFL-16-5, qFL-25-2and qFL-25-3) could
be stably detected in more environments. The qFL-25-
2 linked to CGR5525 could explain 2.84%, 2.57%, and
5.51% of the observed phenotypic variations with addi-
tive effects of −0.57, −0.50, and −0.57 mm in PopE1,
PopE2, and PopE3, respectively. The qFL-16-5 linked to
NAU5408 and NAU3594 could explain 5.93% and 9.60%
of the observed phenotypic variations in PopE1 and PopE3
with the additive effect of 0.23 and 0.29 mm in two gener-
ations, respectively. The qFL-25-3 linked to GH537 could
explain 3.39% and 6.43% of the observed phenotypic vari-
ations in PopE1 and PopE3 with the additive effect of 0.22
and 0.50 mm in two generations, respectively.
Table 3 Distribution of QTL on
chromosomes
FL fiber length, FS fiber strength, FM fiber micronaire, FU fiber uniformity, FE fiber elongation, BW boll
weight, LP lint percentage
BW (g) LP (%) FL (mm) FS (cN·tex-1) FM (unit) FU (%) FE (%) Total
Chr01(A1) 0 1 0 1 1 1 1 5
Chr02(A2) 1 0 1 2 1 0 1 6
Chr04(A4) 1 0 1 0 0 0 0 3
Chr06(A6) 1 0 1 1 1 0 1 5
Chr07(A7) 1 1 2 3 1 1 2 11
Chr08(A8) 0 0 0 0 0 1 0 1
Chr09(A9) 0 1 1 1 0 1 2 6
Chr10(A10) 0 1 1 1 1 0 0 4
Chr15(D1) 1 4 2 3 1 4 5 20
Chr16(D7) 4 0 5 5 3 0 1 18
Chr17(D3) 0 0 1 2 1 2 0 6
Chr19(D5) 1 0 2 1 2 2 2 10
Chr20(D10) 0 0 0 0 1 0 0 1
Chr22(D4) 1 1 0 0 1 1 0 4
Chr23(D9) 0 1 0 1 0 0 0 2
Chr24(D8) 0 0 0 1 1 1 1 4
Chr25(D6) 3 2 4 5 3 2 4 23
Total 14 12 21 27 18 16 21 129
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1129Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
Fiber strength
There were 27 QTL for FS on 13 chromosomes with the
PVE ranging from 0.85% to 13.19%; 9 and 18 of them
were distributed in the At- and Dt subgenomes, respec-
tively. Chr16 and Chr25 were the first 2 chromosomes with
the most QTL. Seventeen QTL with additive effects from
0.03 to 3.11 cN tex−1 indicated that Hai1 alleles increased
FS, nine QTL with additive effects from −1.25 to −0.07
cN tex−1 indicated that CCRI45 alleles increased FS.
The qFS-17-1 had unstable genetic effects, which could
be detected in two environments, with additive effects of
0.35 cN tex−1 in PopE2 and −0.07 cN tex−1 in PopE3.
Ten QTL (qFS-02-2, qFS-06-1, qFS-07-2, qFS-10-1,
qFS-16-1, qFS-16-4, qFS-16-5, qFS-19-1, qFS-25-3, and
qFS-25-4) could be stably detected in more environments.
The qFS-16-1 linked to CGR6894a could explain 6.06%,
7.25%, and 7.38% of the observed phenotypic variations
with the additive effect of 0.44, 0.66, and 0.42 cN tex−1
in PopE1, PopE2, and PopE3, respectively. The qFS-16-5
linked to NAU5408 and NAU3594 could explain 6.55%,
7.42%, and 10.92% of the observed phenotypic variations
with additive effects of 0.66, 0.78, and 0.76 cN tex−1 in
PopE1, PopE2, and PopE3, respectively. The qFS-25-
3 linked to DPL0166a and SHIN-0885 could explain
13.09%, 8.16%, and 7.91% of the observed phenotypic
variations with additive effects of 0.04, 3.11, and 0.81 cN
tex−1 in PopE1, PopE2, and PopE3, respectively. The qFS-
25-4 linked to CGR5201 and SHIN1131 could explain
11.95%, 6.95%, and 4.90% of the observed phenotypic
variations with additive effects of 0.89, 0.81, and 0.54 cN
tex−1 in PopE1, PopE2, and PopE3, respectively. The qFS-
02-2 linked to DPL0450 and PGML04760 could explain
2.10% and 0.85% of the observed phenotypic variations in
PopE2 and PopE3 with additive effects of 1.10 and 0.65cN
tex−1, respectively. The qFS-06-1 linked to DC40067 and
DPL0918 could explain 3.99% and 5.92% of the observed
phenotypic variations in PopE2 and PopE3 with additive
effects of 1.48 and 1.16 cN tex−1, respectively. The qFS-
07-2 linked to NAU2002 and CGR6381 could explain
5.85% and 10.38% of the observed phenotypic variations
with additive effects of 0.58 and 0.78 cN tex−1 in PopE1
and PopE3, respectively. The qFS-10-1 linked to DPL0468
could explain 5.03% and 5.95% of the observed pheno-
typic variations with additive effects of 0.66 and 0.48 cN
tex−1 in PopE1 and PopE3, respectively. The qFS-16-4
linked to PGML02608 could explain 6.98% and 6.32% of
the observed phenotypic variations with additive effect of
0.66 and 0.55 cN tex−1 in PopE1 and PopE3, respectively.
The qFS-19-1 linked to DC40425 and BNL3089 could
explain 4.90% and 2.32% of the observed phenotypic vari-
ations with additive effects of 0.21 and 0.15 cN tex−1 in
PopE1 and PopE2, respectively.
Fiber micronaire
A total of 18 QTL for FM on 13 chromosomes with the
PVE ranging from 1.60 to 10.28% were identified; 5 and
13 of them were distributed in the At- and Dt subgenomes,
respectively. Chr16 and Chr25 were the two most prominent
chromosomes containing the largest number of QTL. Two
QTL (qFM-16-1 and qFM-17-1) with additive effects of 0.03
and 0.18 indicated that CCRI45 alleles decreased FM, and
the remaining 16 QTL with additive effects from −0.52 to
−0.06 indicated that Hai1 alleles decreased FM. Nine QTL
(qFM-02-1, qFM-10-1, qFM-15-1, qFM-16-2, qFM-19-1,
qFM-22-1, qFM-24-1, qFM-25-1, and qFM-25-2) could be
stably detected in more environments. The qFM-02-1 linked
to PGML02861 and DPL0450 could explain 2.07%, 9.24%,
and 6.25% of the observed phenotypic variations with addi-
tive effects of −0.14, −0.52, and −0.15 in PopE1, PopE2,
and PopE3, respectively. The qFM-15-1 linked to NAU3177
could explain 2.35%, 3.22%, and 4.73% of the observed phe-
notypic variations with additive effects of −0.14, −0.08 and
−0.16 in PopE1, PopE2 and PopE3, respectively. The qFM-
10-1 linked to DPL0468 could explain 4.56% and 3.56% of
the observed phenotypic variations with additive effects of
−0.17 and −0.13 in PopE1 and PopE3, respectively. The
qFM-16-2 linked to HAU1836 and BNL2634 could explain
9.88% and 1.98% of the observed phenotypic variations with
additive effect of −0.29 and −0.11 in PopE2 and PopE3,
respectively. The qFM-19-1 linked to PGML01289 could
explain 1.88% and 1.72% of the observed phenotypic vari-
ations with additive effects of −0.10 for both PopE2 and
PopE3, respectively. The qFM-22-1 linked to JESPR230 and
DPL0489 could explain 5.76% and 3.28% of the observed
phenotypic variations in PopE2 and PopE3, respectively,
with the additive effect of −0.08 and −0.15. The qFM-24-
1 linked to NAU2914 could explain 1.87% and 1.60% of
the observed phenotypic variations with additive effects of
−0.09 and −0.08 in PopE1 and PopE3, respectively. The
qFM-25-1 linked to DPL0166a and Gh537 could explain
2.05% and 6.22% of the observed phenotypic variations
with additive effects of −0.06 and −0.25 in PopE1 and
PopE3, respectively. The qFM-25-2 linked to BNL3806 and
TMB0313 could explain 10.28% and 5.91% of the observed
phenotypic variations in PopE2 and PopE3, respectively,
with additive effects of −0.18 and −0.09.
Fiber uniformity
Sixteen QTL for FU were identified on 10 chromosomes
with PVE ranging from 1.06% to 12.67%; 4 and 12 of them
were distributed in the At- and Dt subgenomes, respec-
tively. Chr15 was the most prominent chromosome with the
most QTL. Eleven QTL with additive effects from 0.01 to
0.65% indicated that Hai1 alleles increased FU. Two QTL
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1130 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
(qFU-19-1 and qFU-24-1) with additive effects of −0.49%
and −0.28% indicated that CCRI45 alleles increased FU.
Three QTL (qFU-15-1, qFU-15-3 and qFU-15-4) could
be detected in two environments, but had unstable genetic
effects. Two QTL (qFU-01-1 and qFU-07-1) could be sta-
bly detected in two environments. The qFU-01-1 linked to
BNL2921 and NAU3901 could explain 8.44% and 5.59%
of the observed phenotypic variations in PopE2 and PopE3,
respectively, with additive effects of 0.37% and 0.12%. The
qFU-07-1 linked to NAU1048 and CICR0226 could explain
1.89% and 1.82% of the observed phenotypic variations in
PopE2 and PopE3, respectively, with the additive effect of
0.38% and 0.1%.
Fiber elongation
There were 21 QTL for FE on 11 chromosomes with PVE
ranging from 1.32% to 10.06%; 8 and 13 of these QTL were
distributed in the At- and Dt subgenomes, respectively.
Chr15 and Chr25 were the top 2 chromosomes with the larg-
est number of QTL. Twelve QTL with additive effects from
0.001 to 0.05% indicated that Hai1 alleles increased FE, and
eight QTL with additive effects from −0.14 to −0.01% indi-
cated that CCRI45 alleles increased FE. The qFE-02-1 could
be detected in PopE1 and PopE2, but had unstable genetic
effects. Five QTL (qFE-06-1, qFE-16-1, qFE-19-2, qFE-
25-2, and qFE-25-4) could be stably detected in two envi-
ronments. The qFE-06-1 linked to DC40067 and DPL0918
could explain 6.87% and 3.25% of the observed phenotypic
variations in PopE2 and PopE3, respectively, with addi-
tive effects of 0.02% and 0.04%. The qFE-16-1 linked to
BNL2634 could explain 6.5% and 3.52% of the observed
phenotypic variations in PopE1 and PopE2, respectively,
with additive effects of −0.04% and −0.09%. The qFE-19-
2 linked to DC40425 and HAU1400 could explain 7.13%
and 6.31% of the observed phenotypic variations in PopE2
and PopE3, respectively, with additive effects of −0.01%
and −0.04%. The qFE-25-2 linked to CICR0701 and Gh537
could explain 7.14% and 4.64% of the observed phenotypic
variations in PopE1 and PopE3, respectively, with addi-
tive effects of 0.01% and 0.04%. The qFE-25-4 linked to
DPL0124 could explain 2.21% and 5.97% of the observed
phenotypic variations in PopE2 and PopE3, respectively,
with additive effects of −0.05% and −0.03%.
Boll weight
A total of 14 QTL for BW were identified on 9 chromosomes
with PVE ranging from 5.80 to 15.57%; 4 and 10 of them
were distributed in the At- and Dt subgenomes, respectively.
Chr16 and Chr25 were the first 2 chromosomes with the
most of QTL. Five QTL with additive effects from 0.01 to
0.28 g indicated that Hai1 alleles increased BW, and nine
QTL with additive effects from −0.57 to −0.01 g indicated
that CCRI45 alleles increased BW. Five QTL (qBW-06-1,
qBW-07-1, qBW-16-2, qBW-16-4 and qBW-25-3) could be
stably detected in more environments. The qBW-06-1 linked
to DC40067 and DPL0918 could explain 6.61%, 6.72%,
and 6.07% of the observed phenotypic variations with addi-
tive effects of −0.11, −0.22, and −0.08 g in the PopE1,
PopE2, and PopE3, respectively. The qBW-07-1 linked to
NAU2002 and NAU1085 could explain 8.61% and 6.98% of
the observed phenotypic variations with additive effects of
−0.36 and −0.20 g in PopE1 and PopE2, respectively. The
qBW-16-2 linked to HAU1836 and BNL2634 could explain
6.18% and 8.55% of the observed phenotypic variations
with additive effects of −0.31 and −0.18 gin PopE1 and
PopE2, respectively. The qBW-16-4 linked to NAU5408 and
NAU3594 could explain 7.78% and 6.92% of the observed
phenotypic variations with additive effects of −0.35 and
−0.21 gin PopE1 and PopE2, respectively. The qBW-25-3
linked to BNL3806 and SHIN1131 could explain 7.85% and
15.57% of the observed phenotypic variations in the PopE2
and PopE3, respectively, with additive effects of −0.22 and
−0.20 g.
Lint percentage
Twelve QTL for LP were identified on 8 chromosomes with
the PVE ranging from 3.39% to 30.35%, 4 and 8 of them
were distributed in At- and Dt subgenomes, respectively.
Chr15 was the most prominent chromosome with the larg-
est number of QTL. Nine QTL with the additive effect from
0.01 to 0.84% indicated that Hai1 alleles increased LP, and
four QTL with the additive effect from −2.28 to −0.52%
indicated that CCRI45 alleles increased LP. Five QTL
(qLP-09-1, qLP-15-1, qLP-15-2, qLP-15-3 and qLP-15-4)
could be stably detected in more environments. The qLP-
09-1 linked to DPL0171 could explain 30.35% and 4.92%
of the observed phenotypic variations with additive effects
of −2.28% and −0.60% in PopE1 and PopE3, respectively.
The qLP-15-1 linked to DPL0346a could explain 4.19% and
5.85% of the observed phenotypic variations with additive
effects of 0.01% and 0.52% in PopE1 and PopE3 respec-
tively. The qLP-15-2 linked to MUSS085 and SWU0280b
could explain 5.80% and 11.47% of the observed pheno-
typic variations with additive effects of 0.78% and 0.84%
in the PopE1 and PopE3, respectively. The qLP-15-3 linked
to MUCS410 and HAU0059 could explain 5.55%, 6.54%
and 5.42%of the observed phenotypic variations with addi-
tive effects of 0.01%, 0.25% and 0.72% in PopE1, PopE2
and PopE3, respectively. The qLP-15-4 linked to NAU5138
could explain 6.04%, 6.82% and 7.14% of the observed phe-
notypic variations with additive effects of 0.01%, 0.73% and
0.79% in PopE1, PopE2 and PopE3, respectively.
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1131Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
QTL cluster
QTL clusters are chromosome regions that contain multiple
QTL (≥ 3) related to various traits (Rong etal. 2007). In
the present study, 26 QTL clusters including 115 QTL were
identified on 23 introgressive segments of 14 chromosomes
(Table3, Fig. S2); 8 and 19 of them were distributed in the
At- and Dt subgenomes, respectively. The genetic length of
the clusters varied from 1 to 22 cM and was concentrated
between 1 and 5 cM. There were more QTL clusters on
Chr15, Chr16, and Chr25 than on the other chromosomes.
Twenty-three QTL clusters had stable QTL with the same
additive effect direction in more environments, and 11 of
these QTL clusters had stable QTL for FS or FL. Two QTL
clusters (Clu-16-5 and Clu-25-3) had stable QTL both for
FS and FL. The Clu-16-5 on Chr16 at 78-80 cM had 4 QTL,
the additive effects indicated that Hai1 alleles increased FL
and FS, but decreased FM and BW. The Clu-25-3 on Chr25
at 25-28 cM had five QTL, and the additive effects indicated
that Hai1 alleles increased FL, FS, FU and LP, but decreased
BW. Eight QTL clusters (Clu-02-1, Clu-06-1, Clu-07-2,
Clu-10-1, Clu-16-1, Clu-16-4, Clu-19-2, and Clu-25-4) had
stable QTL for FS. Except for Clu-16-1, and the additive
effects indicated that Hai1 alleles increased FL and FS, but
decreased FM and BW. One QTL cluster (Clu-25-1) had
stable QTL for FL, the additive effects indicated that Hai1
decreased FL and FS.
Discussion
Selection ofgenetic linkage map andimportance
ofCSSLs
In total, 23,569 pairs of SSR markers distributed in the
whole genome were used to screen for polymorphisms
between CCRI36 (G.h) and Hai1 (G.b). A genetic link-
age map with 2292 SSR loci on the 26 cotton chromo-
somes was developed from a BC1F1 population of CCRI36
× Hai1, covering the whole tetraploid cotton genome of
5115.16 cM with an average distance of 2.23 cM between
adjacent markers (Shi etal. 2015). Genetic diversity is
considered a critical issue in QTL mapping of complex
traits. CSSLs have the potential to enrich the diversity of
genetic background and uncover favorable alleles related
to important fiber yield and quality traits(Ali etal. 2010;
Wu etal. 2010; Tyagi etal. 2014). In addition, CSSLs are
ideal materials for QTL mapping, genetic effects identify-
ing and gene cloning(Yang etal. 2015), and are more con-
venient to study the minor and dominant effects of genes
(Wan etal. 2004; He etal. 2015; Li et al. 2016; Qiao
etal. 2016). In this study, the female parent MBI7561 was
selected from a CSSL constructed by G.h and G.b, which
had stable and significant advantages for fiber quality com-
pared with the recurrent parent of CCRI45, to produce
BC5F1 (Table4). The length of introgressed segments from
Hai1 for MBI7561 was 229.47 cM, of which homozygous
introgressed segments were significantly longer than het-
erozygous introgressed segments. The total proportion of
the introgressed segments for MBI7561 was small at the
whole genome (4.40%). However, the proportion of intro-
gression in this study was relatively larger than that in
the previous studies (Song etal. 2017; Guo etal. 2018).
Thus, it is necessary to purify the genetic background
using self-crossing and backcrossing. The diversity of the
MBI7561 genotypes and the dominance of the phenotype
were mutually beneficial, which showed that introgression
segments had important effects on the phenotype. Moreo-
ver, this study laid an important foundation to construct a
segregating population for identifying a large number of
QTL related to fiber yield and quality. The performance of
phenotypes (Tables1, 2, Fig. S1) and genotypes (Fig.1)
Table 4 Phenotypic
performance of fiber quality and
yield-related traits for parents
FL fiber length, FS fiber strength, FM fiber micronaire, FU fiber uniformity, FE fiber elongation, BW boll
weight, LP lint percentage
*Significant difference at the 0.05 level (2-tailed), **significant difference at the 0.01 level (2-tailed)
Year Parents BW (g) LP (%) FL (mm) FS (cN/tex) FM (unit) FU (%) FE (%)
2012 CCRI45 5.69 36.83 28.58 28.52 4.2 84.39 6.5
MBI7561 4.70** 36.25 31.83** 34.75** 3.85** 86.4** 6.5
Hai1 3.18** 32.31** 32.18** 37.7** 4.81**
2013 CCRI45 6.02 33.67 29.73 29.96 4.86 85.41 6.1
MBI7561 4.8** 37.04** 32.13** 36.65** 3.97** 86.40* 5.9**
Hai1 2.7** 30.9** 33.53** 37.7** 4.21** 83.30** 6.0**
2014 CCRI45 5.63 37.54 29.48 28.02 4.83 83.98 6.8
MBI7561 4.99** 39.56** 31.58** 33.59** 4.15** 84.86* 6.8
2015 CCRI45 5.38 36.74 28.57 28.65 5.07 84.35 6.8
MBI7561 4.72** 39.28** 31.5** 32.6** 4.40** 85.00** 7.0**
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1132 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
in the BC5F2 and BC5F2:3 populations showed a large
number of superior materials and stable QTL in multiple
environments.
Common QTL shared withprevious research
A total of 129 QTL were identified in the present study,
and 39 of them were stable (TableS1, Fig. S2). When QTL
had the same linkage markers or confidence interval overlap
between previous studies and our studies, we defined it as
a common QTL interval (Zhang etal. 2016b). Fifty three
common QTL detected in the present study were reported in
the previous researches (Zhang etal. 2008, 2015b, c, 2016a,
b; Liu 2009; Qin etal. 2009; Yang etal. 2009; Liang etal.
2010; Jia etal. 2011; Zhang 2012; Wang etal. 2013, 2016a,
b; He 2014; Ma 2014; Cao etal. 2015; Guo etal. 2015; Nie
etal. 2015; Rong etal. 2015; You 2015; Jamshed etal. 2016;
Ademe etal. 2017; Li 2017; Song etal. 2017; Guo etal.
2018), and 15 of them were stable (TableS3). The remain-
ing 76 of the 129 QTL were identified for the first time,
and 24 of the 39 stable QTL(qFL-16-5(+), qFL-25-2(−),
qFS-02-2(+), qFS-16-1(+), qFS-16-4(+), qFS-16-5(+),
qFS-19-1(+), qFM-02-1(−), qFM-10-1(−), qFM-15-1(−),
qFM-16-2(−), qFM-19-1(−), qFU-01-1(+), qFE-06-1(+),
qFE-16-1(−), qFE-19-2(−), qFE-25-4(−), qBW-07-1(−),
qBW-16-2(−), qBW-16-4(−), qBW-25-3(−), qLP-09-1(−),
qLP-15-2(+) and qLP-15-4(+)) were newly identified stable
QTL in the present study. Both stable QTL and common
QTL were suggestive of the stabilities of the genetic effects
(Zhai etal. 2016). Therefore, these 77 stable or common
QTL were very important for marker assisted breeding and
exploration of genetic mechanisms.
QTL clusters withcommon andstable QTL
QTL clusters are common phenomena in cotton (Said etal.
2015; Wang etal. 2015; Zhai etal. 2016; Song etal. 2017).
A large number of QTL were enriched in these hotspots,
which indicated that these chromosome segments might
contain pleiotropic or linked genes related to different traits
(Abdelraheem etal. 2017). These clustered QTL may belong
to the same genetic factor group contributing to the complex
network of fiber development and affecting the multiple fiber
traits(Lacape etal. 2010). QTL clusters allow cotton breed-
ers to focus their efforts on regions with pleiotropic or linked
loci. In the present study, a total of 26 QTL clusters were
identified; 23 of them had stable QTL and 22 of them had
49 common QTL (TableS2, TableS3) (Zhang etal. 2008,
2015b, c, 2016a, b; Liu 2009; Qin etal. 2009; Yang etal.
2009; Liang etal. 2010; Jia etal. 2011; Zhang 2012; Wang
etal. 2013; He 2014; Ma 2014; Cao etal. 2015; Guo etal.
2015; Nie etal. 2015; Rong etal. 2015; You 2015; Jamshed
etal. 2016; Wang etal. 2016a, b; Ademe etal. 2017; Li
2017; Song etal. 2017; Guo etal. 2018). The remaining
four clusters (Clu-02-1, Clu-04-1, Clu-16-4 and Clu-16-
5) were new. Twenty-five QTL clusters contained the 77
stable or common QTL with stable genetic effects, which
were regarded as stable QTL clusters. There were 19 stable
QTL clusters for FL or FS, 7 (Clu-07-1, Clu-07-2, Clu-16-
1, Clu-16-5, Clu-19-2, Clu-25-3 and Clu-25-4) of them for
FL and FS, 5 (Clu-15-3, Clu-16-2, Clu-16-3, Clu-17-1 and
Clu-25-1) of them mainly for FL, and 7 (Clu-02-1, Clu-06-
1, Clu-10-1, Clu-15-1, Clu-15-4, Clu-16-4 and Clu-24-1) of
them mainly for FS.
Fifteen of the 19 stable QTL clusters had common QTL
for FS or FL. Four QTL clusters (Clu-07-1, Clu-07-2, Clu-
25-3 and Clu-25-4) had common QTL for both FS and
FL, and the additive effects indicated most of Hai1 alleles
increased FL, FS and LP, but decreased FM and BW. Six
QTL clusters (Clu-15-3, Clu-16-1, Clu-16-2, Clu-16-3, Clu-
17-1 and Clu-19-2) had common QTL for FL. Except for
Clu-15-3 and Clu-16-3, the additive effects of the others
indicated Hai1 alleles increased FL and FS. Five QTL clus-
ters (Clu-06-1, Clu-10-1, Clu-15-1, Clu-15-4 and Clu-24-1)
had common QTL for FS. The additive effects of Clu-15-
1, Clu-15-4 and Clu-24-1 indicated that most of the Hai1
alleles increased FL and FS.
Three (Clu-02-1, Clu-16-4 and Clu-16-5) of the 4 new
QTL clusters were new stable QTL clusters. Their addi-
tive effects were similar to the effects of the 4 QTL clusters
which had common QTL for both FS and FL.
Genetic eects oftheimportant chromosome
introgressed segments
Analyzing the genetic effects of the chromosome segments is
necessary to develop a breeding strategy with precise direc-
tion in MAS (Zhai etal. 2016; Song etal. 2017), and to
explore genetic mechanisms. In this study, the genetic effects
of 29 chromosome introgressed segments were identified for
fiber quality and yield-related traits. The genetic effects of 23
chromosome introgressed segments were identical to those
of the QTL clusters associated with them. Four segments
(Seg-A08-2, Seg-A09-1, Seg-D09-1, and Seg-D06-3) had
stable genetic effects without any QTL clusters (TableS2).
A total of 13 chromosome introgressed segments were
important for fiber quality (10) and LP (4) improvement on
the 7 chromosomes [Chr15(D1), Chr02(A2), Chr19(D5),
Chr06(A6)–Ch25(D6), Chr07(A7)–Chr16(D7)]. The
additive effects of Seg-D06-2 with 3 stable QTL clusters
(with stable or common QTL) on Chr25(D6) at 23-33 cM
indicated that most of the introgressed Hai1 alleles stably
increased FL, FS, FU and FE, but stably decreased FM.
The additive effects of 9 segments (Seg-A07-1, Seg-A07-2,
Seg-A07-3, Seg-A02-1, Seg-A06-1, Seg-A10-2, Seg-D07-1,
Seg-D07-3, and Seg-D05-2), with single stable QTL cluster
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1133Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
(with stable or common QTL), indicated that most of the
introgressed Hai1 alleles stably increased FL and FS, but
stably decreased FM (except for Seg-D07-1). In addition,
the additive effects of 4 important chromosome introgressed
segments (Seg-A07-1, Seg-D01-2, Seg-D01-3, and Seg-D01-
4) with single-stable QTL cluster (with stable or common
QTL) indicated that the introgressed Hai1 alleles stably
increased LP.
Distribution ofimportant chromosome introgressed
segments forexcellent lines
To further verify which of the 13 important chromosome
introgressed segments had major genetic effects, we ana-
lyzed the distribution of introgressed segments for 5 excel-
lent lines with the best comprehensive phenotypes of FL,
FS, FM, and LP in multiple populations (Table5, Fig.2).
There were 6 common important chromosome introgressed
segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3,
Seg-D07-3 and Seg-D06-2) in all five lines, and the additive
effects for the related QTL indicated that the introgressed
Hai1 alleles stably improved FL, FS and FM. There were
three important chromosome introgressed segments (Seg-
D01-2, Seg-D01-3 and Seg-D01-4), with stable genetic
effects for LP, enriched in the region (169-189 cM) on
Chr15(D1). The introgressed region was detected in all 5
excellent lines. The distribution of the important chromo-
some introgressed segments (or region) in excellent lines
further confirmed that further research should be focused
on the 7 common important chromosome introgressed
segments(or region) for fine mapping, genetic mechanisms
exploring, and MAS breeding applications.
In conclusion, developing CSSLs is an effective method
for identifying genetic effects. We selected the CSSLs
Table 5 Phenotypics and introgressed segments distribution of 5 lines with excellent fiber quality
Ind. ID Phenotype value of fiber quality in three environments Chromosome introgressed seg-
ments
Fiber length (mm) Fiber strength (cN/tex) Fiber micronaire Lint percentage (%) Number Len (cM) Rate (%)
Pop
E1
Pop
E2
Pop
E3
Pop
E1
Pop
E2
Pop
E3
Pop
E1
Pop
E2
Pop
E3
Pop
E1
Pop
E2
Pop
E3
288-02 34.58 31.30 32.50 34.70 36.70 32.70 4.51 4.10 5.20 41.50 36.20 39.90 14 106.00 2.07
291-06 35.29 31.20 34.10 37.90 36.70 33.40 3.90 3.80 3.90 38.18 35.72 39.91 14 100.00 1.95
291-19 32.62 32.60 31.20 35.90 40.50 30.90 3.83 4.30 4.80 34.86 36.01 41.47 13 104.00 2.03
292-05 34.43 33.10 33.10 39.30 38.40 32.40 3.24 4.00 5.10 33.45 35.45 38.86 14 107.00 2.09
243-04 30.06 33.90 31.00 31.70 39.40 33.30 3.95 3.90 3.80 37.17 37.31 38.61 11 71.00 1.39
MBI7561 31.58 31.50 31.01 33.59 32.60 30.40 4.15 4.40 4.42 39.56 39.28 40.80 33 127.00 2.48
CCRI45 29.48 28.57 28.07 28.02 28.65 26.25 4.83 5.07 5.09 37.54 36.74 38.20 –
Fig. 2 Genotype distribution of individuals with excellent fiber quality
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1134 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
containing 33 chromosome introgression segments as a
parent to construct the three segregated populations. A
total of 129 QTL associated with fiber quality (103) and
yield-related traits (26) were detected on 17 chromosomes,
explaining 0.85–30.35% of the phenotypic variation, 39
were stable, 53 were common, 76 were new, and 86 had
favorable effects on the related traits. More QTL were dis-
tributed in the Dt subgenome than in the At subgenome.
Twenty-five stable QTL clusters (with stable or common
QTL) were detected on 22 chromosome introgressed seg-
ments. Finally, the 6 important chromosome introgressed
segments (Seg-A02-1, Seg-A06-1, Seg-A07-2, Seg-A07-3,
Seg-D07-3 and Seg-D06-2) were identified as candidate
chromosome regions for fiber quality, which should be given
more attention in future QTL fine mapping, gene cloning,
and MAS breeding.
Acknowledgements This study was funded by the National Key R
& D Program for Crop Breeding (2016YFD0100306), the National
Natural Science Foundation of China (31101188) and the Agricul-
tural Science and Technology Innovation Program for CAAS (CAAS-
ASTIP-ICRCAAS). Thanks to the Quantitative Genetics Group of
CAAS (Beijing, China) providing the software ICIMapping and help
in QTL identification.
Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of
interest.
Ethical approval This article does not contain any studies with human
participants or animals performed by any of the authors.
Informed consent Informed consent was obtained from all individual
participants included in the study.
Open Access This article is distributed under the terms of the Crea-
tive Commons Attribution 4.0 International License (http://creat iveco
mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribu-
tion, and reproduction in any medium, provided you give appropriate
credit to the original author(s) and the source, provide a link to the
Creative Commons license, and indicate if changes were made.
References
Abdelraheem A, Liu F, Song M, Zhang JF (2017) A meta-analysis of
quantitative trait loci for abiotic and biotic stress resistance in
tetraploid cotton. Mol Genet Genomics 292:1221–1235
Ademe MS, He S, Pan Z, Sun J, Wang Q, Qin H, Liu J, Liu H, Yang
J, Xu D, Yang J, Ma Z, Zhang J, Li Z, Cai Z, Zhang X, Zhang X,
Huang A, Yi X, Zhou G, Li L, Zhu H, Pang B, Wang L, Jia Y, Du
X (2017) Association mapping analysis of fiber yield and qual-
ity traits in Upland cotton (Gossypium hirsutum L.). Mol Genet
Genomics 292:1267–1280
Ali ML, Sanchez PL, Yu SB, Lorieux M, Eizenga GC (2010) Chro-
mosome segment substitution lines: a powerful tool for the
introgression of valuable genes from Oryza wild species into cul-
tivated rice (O. sativa). Rice 3:218–234
Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Char-
cosset A, Joets J (2004) BioMercator: integrating genetic maps
and QTL towards discovery of candidate genes. Bioinformatics
20:2324–2326
Brown JS, Schnell RJ, Motamayor JC, Lopes U, Kuhn DN, Borrone
JW (2005) Resistance gene mapping for witches’ broom disease
in Theobroma cacao L. in an F2 population using SSR markers
and candidate genes. J Am Soc Hortic Sci 130:366–373
Cai G, Yang Q, Yi B, Fan C, Zhang C, Edwards D, Batley J, Zhou Y
(2015) A bi-filtering method for processing single nucleotide
polymorphism array data improves the quality of genetic map
and accuracy of quantitative trait locus mapping in doubled hap-
loid populations of polyploid Brassica napus. BMC Genomics.
https ://doi.org/10.1186/s1286 4-015-1559-4
Cao Z, Zhu X, Chen H, Zhang T (2015) Fine mapping of clustered
quantitative trait loci for fiber quality on chromosome 7 using
a Gossypium barbadense introgressed line. Mol Breed. https ://
doi.org/10.1007/s1103 2-015-0393-3
Chen X, Guo W, Liu B, Zhang Y, Song X, Cheng Y, Zhang L, Zhang
T (2012) Molecular mechanisms of fiber differential develop-
ment between G. barbadense and G. hirsutum revealed by
genetical genomics. PLoS ONE 7:e30056
Clement JD, Constable GA, Stiller WN, Liu SM (2012) Negative
associations still exist between yield and fibre quality in cot-
ton breeding programs in Australia and USA. Field Crops Res
128:1–7
Eshed Y, Zamir D (1994) A genomic library of Lycopersicon pennel-
lii in L.esculentum: a tool for fine mapping of genes. Euphytica
79:175–179
Fu Y, Yuan D, Hu W, Cai C, Guo W (2013) Development of Gossyp-
ium barbadense chromosome 18 segment substitution lines in the
genetic standard line TM-1 of Gossypium hirsutum and mapping
of QTLs related to agronomic traits. Acta Agron Sin 39:21–28
Ge R, Lan M, Shi Y, Li J, Liu A, Wang T, Yuan Y (2012) Correla-
tion and path coefficient analysis of main agronomic characters
in BC4F3 and BC4F4 generations from Gossypium hirsutum L.
×Gossypium barbadense L. Chin Agric Sci Bull 28:127–130
Guo L, Shi Y, Li J, Gong J, Liu A, Shang H, Gong W, Chen T, Ge Q,
Sun J, Yuan Y (2015) Mapping QTL of fiber yield and quality
traits in F2 populations of chromosome segment substitution lines
from Gossypium hirsutum × Gossypium barbadense. Cotton Sci
27:550–560
Guo L, Shi Y, Gong J, Liu A, Tan Y, Gong W, Li J, Chen T, Shang H,
Ge Q, Lu Q, Sun J, Yuan Y (2018) Genetic analysis of the fiber
quality and yield traits in G. hirsutum background using chro-
mosome segments substitution lines (CSSLs) from Gossypium
barbadense. Euphytica 214:82. https ://doi.org/10.1007/s1068
1-018-2158-7
He R (2014) The evaluation and identifying QTL of chromosome
segment substitution lines (BC5F3, BC5F3:4, BC5F3:5) in CCRI36
background of Gossypium hirsutum L. [D], Chinese Academy of
Agricultural Sciences cnki.net (S562, in Chinese)
He Q, Yang H, Xiang S, Tian D, Wang W, Zhao T, Gai J (2015) Fine
mapping of the genetic locus L1 conferring black pods using a
chromosome segment substitution line population of soybean.
Plant Breed 134:437–445
Jamshed M, Jia F, Gong J, Palanga KK, Shi Y, Li J, Shang H, Liu A,
Chen T, Zhang Z, Cai J, Ge Q, Liu Z, Lu Q, Deng X, Tan Y, Or
Rashid H, Sarfraz Z, Hassan M, Gong W, Yuan Y (2016) Identifi-
cation of stable quantitative trait loci (QTLs) for fiber quality traits
across multiple environments in Gossypium hirsutum recombinant
inbred line population. BMC Genomics. https ://doi.org/10.1186/
s1286 4-016-2560-2
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1135Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
Jia F, Sun F, Li J, Liu A, Shi Y, Gong J, Shang H, Gong W, Wang
T, Liu Z, Yuan Y (2011) Identification of QTL for boll weight
and lint percentage of Upland cotton RIL population in multiple
environments. Mol Plant Breed 9:318–326
Jia X, Pang C, Wei H, Wang H, Ma Q, Yang J, Cheng S, Su J, Fan
S, Song M, Wusiman N, Yu S (2016) High-density linkage map
construction and QTL analysis for earliness-related traits in Gos-
sypium hirsutum L. BMC Genomics. https ://doi.org/10.1186/
s1286 4-016-3269-y
Kushanov FN, Pepper AE, Yu JZ, Buriev ZT, Shermatov SE, Saha S,
Ulloa M, Jenkins JN, Abdukarimov A, Abdurakhmonov IY (2016)
Development, genetic mapping and QTL association of cotton
PHYA, PHYB, and HY5-specific CAPS and dCAPS markers. BMC
Genet. https ://doi.org/10.1186/s1286 3-016-0448-4
Lacape JM, Llewellyn D, Jacobs J, Arioli T, Becker D, Calhoun S,
Al-Ghazi Y, Liu S, Palai O, Georges S, Giband M, de Assun-
cao H, Augusto P, Barroso V, Claverie M, Gawryziak G, Jean
J, Vialle M, Viot C (2010) Meta-analysis of cotton fiber qual-
ity QTLs across diverse environments in a Gossypium hirsutum
× G. Barbadense RIL population. BMC Plant Biol. https ://doi.
org/10.1186/1471-2229-10-132
Lan M, Shi Y, Ge R, Li J, Gong J, Liu A, Shang H, Gong W, Wang
T, Yuan Y (2015) Evaluation of fiber yield and quality traits of
chromosome segment substitution lines population(BC4F3 and
BC4F3:4) from Gossypium hirsutum × Gossypium barbadense. J
Plant Genet Resour 16:257–263
Li B (2017) Genetic analysis and mapping of major QTL on yield and
fiber quality related traits in cotton chromosome segment substitu-
tion lines. [D], Chinese Academy of Agricultural Sciences cnki.
net (S562, in Chinese)
Li F, Jia H, Liu L, Zhang C, Liu Z, Zhang Z (2014) Quantitative trait
loci mapping for kernel row number using chromosome segment
substitution lines in maize. Genet Mol Res 13:1707–1716
Li B, Shi Y, Gong J, Li J, Liu A, Shang H, Gong W, Chen T, Ge Q, Jia
C, Lei Y, Hu Y, Yuan Y (2016) Genetic effects and heterosis of
yield and yield component traits based on Gossypium barbadense
chromosome segment substitution lines in two Gossypium hirsu-
tum backgrounds. PLoS ONE 11:e0157978
Liang Y, Jia Y, Li A, Zhang B, Liu G, Li J, Shi Y, Li J, Liu A, Gong
J, Wang T, Shang H, Gong W, Yuan Y (2010) Phenotyping traits
related to yield and quality of BC5F2 substitution lines in cotton
and their QTL mapping. Mol Plant Breed 8:221–230
Liu G (2009) Development and morphology of cotton fiber and locat-
ing the elite QTL alleles of G.barbadense L. fiber. [D], Hunan
Agricultural University cnki.net (S562, in Chinese)
Liu SB, Zhou RG, Dong YC, Li P, Jia JZ (2006) Development, utiliza-
tion of introgression lines using a synthetic wheat as donor. Theor
Appl Genet 112:1360–1373
Lu L, Wei Q, Wang F, Liu G, Qin D, Guan P, Ni Z, Yao Y, Sun Q, Peng
H (2014) QTL mapping for plant height and internode lengths
in Nongda3338 × Jingdong6 DH populations of wheat (Triticum
aestivum L.). J China Agric Univ 19:1–8
Lu Q, Shi Y, Xiao X, Li P, Gong J, Gong W, Liu A, Shang H, Li J,
Ge Q, Song W, Li S, Zhang Z, Rashid MHO, Peng R, Yuan Y,
Huang J (2017) Transcriptome analysis suggests that chromosome
introgression fragments from Sea Island cotton (Gossypium bar-
badense) increase fiber strength in Upland cotton (Gossypium
hirsutum). G3 7:3469–3479
Ma L (2014) The evaluation and QTL identifying of chromosome
segment substitution lines in advanced backcross of Gossypium
hirsutum × G.barbadense. [D], Chinese Academy of Agricultural
Sciences cnki.net (S562, in Chinese)
Ma Q, Su J, Chen H, Deng F (2014) Research progress of cotton molec-
ular marker-assisted breeding. Guangdong Agric Sci 41:138–143
Nie X, Tu J, Wang B, Zhou X, Lin Z (2015) A BIL population derived
from G. hirsutum and G. barbadense provides a resource for cot-
ton genetics and breeding. PLoS ONE 10:e0141064
Pang MX, Percy RG, Hughs SE, Jones DC, Zhang JF (2012) Identifi-
cation of genes that were differentially expressed and associated
with fiber yield and quality using cDNA-AFLP and a backcross
inbred line population. Mol Breed 30:975–985
Paterson AH, Brubaker CL, Wendel JF (1993) A rapid method for
extraction of cotton (Gossypium spp.) genomic DNA suitable for
RFLP or PCR analysis. Plant Mol Biol Report 11:122–127
Qiao W, Qi L, Cheng Z, Su L, Li J, Sun Y, Ren J, Zheng X, Yang Q
(2016) Development and characterization of chromosome segment
substitution lines derived from Oryza rufipogon in the genetic
background of O. sativa spp. indica cultivar 9311. BMC Genom-
ics 17:89. https ://doi.org/10.1186/s1286 4-016-2987-5
Qin Y-s, Ye W-x, Liu R-z, Zhang T-z, Guo W-z (2009) QTL mapping
for fiber quality properties in Upland cotton (Gossypium hirsutum
L.). Sci Agric Sin 42:4145–4154
Reinisch AJ, Dong JM, Brubaker CL, Stelly DM, Wendel JF, Paterson
AH (1994) A detailed RFLP map of cotton, Gossypium hirsutum
× Gossypium barbadense: chromosome organization and evolu-
tion in a disomic polyploid genome. Genetics 138:829–847
Rong J, Feltus EA, Waghmare VN, Pierce GJ, Chee PW, Draye X,
Saranga Y, Wright RJ, Wilkins TA, May OL, Smith CW, Ganna-
way JR, Wendel JR, Paterson AH (2007) Meta-analysis of poly-
ploid cotton QTL shows unequal contributions of subgenomes to
a complex network of genes and gene clusters implicated in lint
fiber development. Genetics 176:2577–2588
Rong F, Tang L, Tang Y, Li Z, Wei Y (2015) Analysis of QTL map-
ping for fiber quality traits in introgression lines cotton. Molecular
Plant Breeding 13:1509–1516
Said JI, Song M, Wang H, Lin Z, Zhang X, Fang DD, Zhang J (2015)
A comparative meta-analysis of QTL between intraspecific Gos-
sypium hirsutum and interspecific G. hirsutum × G. barbadense
populations. Mol Genet Genomics 290:1003–1025
Shang L, Wang Y, Cai S, Wang X, Li Y, Abduweli A, Hua J (2016)
Partial dominance, overdominance, epistasis and QTL by envi-
ronment interactions contribute to heterosis in two Upland cotton
hybrids. G3 6:499–507
Shi Y, Zhang B, Li J, Liu A, Yuan Y (2008) Cotton fiber quality traits
were controlled mainly by maternal plant genotype. Hereditas
(Beijing) 30:1466–1476
Shi Y, Li W, Li A, Ge R, Zhang B, Li J, Liu G, Li J, Liu A, Shang H,
Gong J, Gong W, Yang Z, Tang F, Liu Z, Zhu W, Jiang J, Yu X,
Wang T, Wang W, Chen T, Wang K, Zhang Z, Yuan Y (2015)
Constructing a high-density linkage map for Gossypium hirsutum
× Gossypium barbadense and identifying QTLs for lint percent-
age. J Integr Plant Biol 57:450–467
Shi Y, Zhang B, Liu A, Li W, Li J, Lu Q, Zhang Z, Li S, Gong W,
Shang H, Gong J, Chen T, Ge Q, Wang T, Zhu H, Liu Z, Yuan Y
(2016) Quantitative trait loci analysis of Verticillium wilt resist-
ance in interspecific backcross populations of Gossypium hirsutum
× Gossypium barbadense. BMC Genomics 17:877. https ://doi.
org/10.1186/s1286 4-016-3128-x
Song W, Wang M, Su W, Lu Q, Xiao X, Cai J, Zhang Z, Li S, Li P,
Gong J, Gong W, Shang H, Liu A, Li J, Chen T, Ge Q, Shi Y,
Yuan Y (2017) Genetic and phenotypic effects of chromosome
segments introgressed from Gossypium barbadense into Gos-
sypium hirsutum. PloS ONE 12(9):e0184882
Stelly DM, Saha S, Raska DA, Jenkins JN, McCarty JC, Gutierrez OA
(2005) Registration of 17 Upland (Gossypium hirsutum) cotton
germplasm lines disomic for different G.barbadense chromosome
or arm substitutions. Crop Sci 45:2663–2665
Su J, Chen H, Deng F, Ma Q, Pang C, Wu M, Wang C, Yu S (2014)
Establishment and application of PCR-SSR based method to
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1136 Molecular Genetics and Genomics (2019) 294:1123–1136
1 3
quickly detect exogenous chromosome fragments in cotton. J
Agric Biotechnol 22:632–641
Tyagi P, Gore MA, Bowman DT, Campbell BT, Udall JA, Kurapar-
thy V (2014) Genetic diversity and population structure in the
US Upland cotton (Gossypium hirsutum L.). Theor Appl Genet
127:283–295
Ulloa M, Saha S, Jenkins JN, Meredith WR, McCarty JC, Stelly DM
(2005) Chromosomal assignment of RFLP linkage groups har-
boring important QTLs on an intraspecific cotton (Gossypium
hirsutum L.) joinmap. J Hered 96:132–144
Van Berloo R (2008) GGT 2.0: versatile software for visualization and
analysis of genetic data. J Hered 99:232–236
Voorrips RE (2002) MapChart: software for the graphical presentation
of linkage maps and QTLs. J Hered 93:77–78
Wan XY, Wan JM, Su CC, Wang CM, Shen WB, Li JM, Wang HL,
Jiang L, Liu SJ, Chen LM, Yasui H, Yoshimura A (2004) QTL
detection for eating quality of cooked rice in a population of chro-
mosome segment substitution lines. Theor Appl Genet 110:71–79
Wan X, Weng J, Zhai H, Wang J, Lei C, Liu X, Guo T, Jiang L, Su
N, Wan J (2008) Quantitative trait loci (QTL) analysis for rice
grain width and fine mapping of an identified QTL allele gw-5
in a recombination hotspot region on chromosome 5. Genetics
179:2239–2252
Wang P, Ding Y, Lu Q, Guo W, Zhang T (2008) Development of Gos-
sypium barbadense chromosome segment substitution lines in the
genetic standard line TM-1 of Gossypium hirsutum. Chin Sci Bull
53:1512–1517
Wang PZ, Shi LF, Su L, Hu BM (2010) Quantitative trait loci for resist-
ance against fusarium wilt based on three cotton F2 populations.
Agric Sci China 9:1799–1806
Wang P, Zhu Y, Song X, Cao Z, Ding Y, Liu B, Zhu X, Wang S, Guo
W, Zhang T (2012) Inheritance of long staple fiber quality traits
of Gossypium barbadense in G.hirsutum background using CSILs.
Theor Appl Genet 124:1415–1428
Wang F, Xu Z, Sun R, Gong Y, Liu G, Zhang J, Wang L, Zhang C,
Fan S, Zhang J (2013) Genetic dissection of the introgressive
genomic components from Gossypium barbadense L. that con-
tribute to improved fiber quality in Gossypium hirsutum L. Mol
Breed 32:547–562
Wang H, Huang C, Guo H, Li X, Zhao W, Dai B, Yan Z, Lin Z (2015)
QTL mapping for fiber and yield traits in Upland cotton under
multiple environments. PLoS ONE 10:e0130742
Wang F, Zhang C, Liu G, Chen Y, Zhang J, Qiao Q, Yuan Z, Fan S,
Zhang J (2016a) Phenotypic variation analysis and QTL mapping
for cotton (Gossypium hirsutum L.) fiber quality grown in different
cotton-producing regions. Euphytica 211:169–183
Wang H, Huang C, Zhao W, Dai B, Shen C, Zhang B, Li D, Lin Z
(2016b) Identification of QTL for fiber quality and yield traits
using two immortalized backcross populations in Upland cotton.
PloS ONE 11:e0166970
Wang Y, Wang X, Li Z, Yang X, Zhang Y, Wu L, Wu J, Zhang G,
Ma Z (2016c) Development of Pima cotton chromosome segment
substitution lines with Gossypium hirsutum background. J Plant
Genet Resour 17:114–119
Wu J, McCarty JC, Jenkins JN, Meredith WR (2010) Breeding poten-
tial of introgressions into Upland cotton: genetic effects and het-
erosis. Plant Breed 129:526–532
Yang Z, Li J, Li A, Zhang B, Liu G, Li J, Shi Y, Liu A, Jiang J, Wang
T, Yuan Y (2009) Developing chromosome segment substitution
lines(CSSLs) in cotton (Gossypium) using advanced backcross
and MAS. Mol Plant Breed 7:233–241
Yang X, Zhou X, Wang X, Li Z, Zhang Y, Liu H, Wu L, Zhang G, Yan
G, Ma Z (2015) Mapping QTL for cotton fiber quality traits using
simple sequence repeat markers, conserved intron-scanning prim-
ers, and transcript-derived fragments. Euphytica 201:215–230
You S (2015) Using chromosome segment introgresion lines to fine
mapping QTL of cotton fiber quality traits on the Chr.A2, Chr.A7,
Chr.A8 and Chr.D5. [D], Nanjing Agricultural University cnki.net
(S562, in Chinese)
Yu J, Yu S, Gore M, Wu M, Zhai H, Li X, Fan S, Song M, Zhang J
(2013) Identification of quantitative trait loci across interspecific
F2, F2:3 and testcross populations for agronomic and fiber traits in
tetraploid cotton. Euphytica 191:375–389
Yu S, Fan S, Wang H, Wei H, Pang C (2016) Progresses in research on
cotton high yield breeding in China. Sci Agric Sin 49:3465–3476
Zhai H, Gong W, Tan Y, Liu A, Song W, Li J, Deng Z, Kong L, Gong
J, Shang H, Chen T, Ge Q, Shi Y, Yuan Y (2016) Identification of
chromosome segment substitution lines of Gossypium barbadense
introgressed in G. hirsutum and quantitative trait locus mapping
for fiber quality and yield traits. PLoS ONE 11:e0159101
Zhang J, Yao ting WU, Guo WZ, Zhang TZ (2000) Fast screening of
microsatellite markers in cotton with PAGE/silver staining. Acta
Gossypii Sin 05:267–269+282
Zhang J, Wang S, Shi Y, Zhang G, Liu A, Li J, Ma Z, Yuan Y (2008)
Molecular marker and QTL for yield-related traits in transgenic
Bacillus thuringiensis (Bt) hybrid cotton. Cotton Sci 20:179–185
Zhang ZS, Hu MC, Zhang J, Liu DJ, Zheng J, Zhang K, Wang W, Wan
Q (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
Zhang J (2012) The evaluation and identifying QTL of fiber yield and
quality traits of chromosome segment substitution lines (BC5F3,
BC5F3:4, BC5F3:5) in Upland cotton. [D], Chinese Academy of
Agricultural Sciences cnki.net (S562, in Chinese)
Zhang J, Yu J, Pei W, Li X, Said J, Song M, Sanogo S (2015a) Genetic
analysis of Verticillium wilt resistance in a backcross inbred line
population and a meta-analysis of quantitative trait loci for disease
resistance in cotton. BMC Genomics. https ://doi.org/10.1186/
s1286 4-015-1682-2
Zhang M, Chen Q, Su X, Geng H, Qu Y (2015b) QTL mapping for 5
fiber quality traits in cotton hybrid populations between Sea Island
cotton and Upland cotton(Gossypium barbadense L. and Gossyp-
ium hirsutum L.). Acat Agric Boreali-Occidentalis Sin 24:64–71
Zhang Z, Li J, Muhammad J, Cai J, Jia F, Shi Y, Gong J, Shang H, Liu
A, Chen T, Ge Q, Palanga KK, Lu Q, Deng X, Tan Y, Li W, Sun
L, Gong W, Yuan Y (2015c) High resolution consensus mapping
of quantitative trait loci for fiber strength, length and micronaire
on chromosome 25 of the Upland cotton (Gossypium hirsutum L.).
PLoS ONE 10:e0135430
Zhang S, Feng L, Xing L, Yang B, Gao X, Zhu X, Zhang T, Zhou B
(2016a) New QTLs for lint percentage and boll weight mined
in introgression lines from two feral landraces into Gossypium
hirsutum acc TM-1. Plant Breed 135:90–101
Zhang SW, Zhu XF, Feng LC, Gao X, Yang B, Zhang TZ, Zhou BL
(2016b) Mapping of fiber quality QTLs reveals useful variation
and footprints of cotton domestication using introgression lines.
Sci Rep 6:31954
Zhao F, Zhang G, Zeng R, Yang Z, Zhu H, Zhong B, Ling Y, He
G (2009) Additive effects and epistasis effects of QTL for plant
height and its components using single segment substitution lines
(SSSLs) in rice. Acta Agron Sin 35:48–56
Publisher’s Note Springer Nature remains neutral with regard to
jurisdictional claims in published maps and institutional affiliations.
Content courtesy of Springer Nature, terms of use apply. Rights reserved.
1.
2.
3.
4.
5.
6.
Terms and Conditions
Springer Nature journal content, brought to you courtesy of Springer Nature Customer Service Center GmbH (“Springer Nature”).
Springer Nature supports a reasonable amount of sharing of research papers by authors, subscribers and authorised users (“Users”), for small-
scale personal, non-commercial use provided that all copyright, trade and service marks and other proprietary notices are maintained. By
accessing, sharing, receiving or otherwise using the Springer Nature journal content you agree to these terms of use (“Terms”). For these
purposes, Springer Nature considers academic use (by researchers and students) to be non-commercial.
These Terms are supplementary and will apply in addition to any applicable website terms and conditions, a relevant site licence or a personal
subscription. These Terms will prevail over any conflict or ambiguity with regards to the relevant terms, a site licence or a personal subscription
(to the extent of the conflict or ambiguity only). For Creative Commons-licensed articles, the terms of the Creative Commons license used will
apply.
We collect and use personal data to provide access to the Springer Nature journal content. We may also use these personal data internally within
ResearchGate and Springer Nature and as agreed share it, in an anonymised way, for purposes of tracking, analysis and reporting. We will not
otherwise disclose your personal data outside the ResearchGate or the Springer Nature group of companies unless we have your permission as
detailed in the Privacy Policy.
While Users may use the Springer Nature journal content for small scale, personal non-commercial use, it is important to note that Users may
not:
use such content for the purpose of providing other users with access on a regular or large scale basis or as a means to circumvent access
control;
use such content where to do so would be considered a criminal or statutory offence in any jurisdiction, or gives rise to civil liability, or is
otherwise unlawful;
falsely or misleadingly imply or suggest endorsement, approval , sponsorship, or association unless explicitly agreed to by Springer Nature in
writing;
use bots or other automated methods to access the content or redirect messages
override any security feature or exclusionary protocol; or
share the content in order to create substitute for Springer Nature products or services or a systematic database of Springer Nature journal
content.
In line with the restriction against commercial use, Springer Nature does not permit the creation of a product or service that creates revenue,
royalties, rent or income from our content or its inclusion as part of a paid for service or for other commercial gain. Springer Nature journal
content cannot be used for inter-library loans and librarians may not upload Springer Nature journal content on a large scale into their, or any
other, institutional repository.
These terms of use are reviewed regularly and may be amended at any time. Springer Nature is not obligated to publish any information or
content on this website and may remove it or features or functionality at our sole discretion, at any time with or without notice. Springer Nature
may revoke this licence to you at any time and remove access to any copies of the Springer Nature journal content which have been saved.
To the fullest extent permitted by law, Springer Nature makes no warranties, representations or guarantees to Users, either express or implied
with respect to the Springer nature journal content and all parties disclaim and waive any implied warranties or warranties imposed by law,
including merchantability or fitness for any particular purpose.
Please note that these rights do not automatically extend to content, data or other material published by Springer Nature that may be licensed
from third parties.
If you would like to use or distribute our Springer Nature journal content to a wider audience or on a regular basis or in any other manner not
expressly permitted by these Terms, please contact Springer Nature at
onlineservice@springernature.com
... At present, the important restricted factors for the further improvement of cotton varieties are the narrow genetic background of Gh [3][4][5] and the negative genetic correlations between fiber quality and yield [6,7]. In order to explore more abundant genetic variation, and develop new cotton varieties with better comprehensive performance of agronomic characters, it is an effective way to introgress the excellent fiber quality major genes of Gb into the genetic background of Gh through constructing cotton chromosome segment substitution lines (CSSLs) [4,[8][9][10][11][12][13]. CSSLs populations are powerful tools for studying the genetic basis of quantitative traits and the effects of gene pyramiding and interaction [8,14], which were first reported as introgression lines (ILs) in tomato [15]. ...
... With the wide application of cotton Simple Sequence Repeats (SSR) and Single Nucleotide Polymorphism (SNP) molecular markers, a large number of important QTL and clusters were identified by CSSLs, focusing on yield-related traits [10,[22][23][24], fiber quality traits [8,19,20,[25][26][27][28], and resistance-related traits [29,30]. However, only a few studies focused on simultaneous improvement of fiber quality and yield related traits in cotton, and most of them have shown that yield component traits were negatively genetically correlated with fiber quality [6,7,[9][10][11][12][31][32][33][34][35][36][37]. Even though some of them have addressed QTL clusters with the positive or negative additive effects for fiber quality and yield component traits but did not combine these results for further analysis of simultaneous improvement [11,31,[33][34][35][36]38]. ...
... However, only a few studies focused on simultaneous improvement of fiber quality and yield related traits in cotton, and most of them have shown that yield component traits were negatively genetically correlated with fiber quality [6,7,[9][10][11][12][31][32][33][34][35][36][37]. Even though some of them have addressed QTL clusters with the positive or negative additive effects for fiber quality and yield component traits but did not combine these results for further analysis of simultaneous improvement [11,31,[33][34][35][36]38]. Studies on the loci or introgression segments that can improve cotton yield and quality simultaneously are rarely reported. ...
Article
Full-text available
Introduction The simultaneous improvement of fiber quality and yield for cotton is strongly limited by the narrow genetic backgrounds of Gossypium hirsutum (Gh) and the negative genetic correlations among traits. An effective way to overcome the bottlenecks is to introgress the favorable alleles of Gossypium barbadense (Gb) for fiber quality into Gh with high yield. Objectives This study was to identify superior loci for the improvement of fiber quality and yield. Methods Two sets of chromosome segment substitution lines (CSSLs) were generated by crossing Hai1 (Gb, donor-parent) with cultivar CCRI36 (Gh) and CCRI45 (Gh) as genetic backgrounds, and cultivated in 6 and 8 environments, respectively. The kmer genotyping strategy was improved and applied to the population genetic analysis of 743 genomic sequencing data. A progeny segregating population was constructed to validate genetic effects of the candidate loci. Results A total of 68,912 and 83,352 genome-wide introgressed kmers were identified in the CCRI36 and CCRI45 populations, respectively. Over 90% introgressions were homologous exchanges and about 21% were reverse insertions. In total, 291 major introgressed segments were identified with stable genetic effects, of which 66(22.98%), 64(21.99%), 35(12.03%), 31(10.65%) and 18(6.19%) were beneficial for the improvement of fiber length (FL), strength (FS), micronaire, lint-percentage (LP) and boll-weight, respectively. Thirty-nine introgression segments were detected with stable favorable additive effects for simultaneous improvement of 2 or more traits in Gh genetic background, including 6 could increase FL/FS and LP. The pyramiding effects of 3 pleiotropic segments (A07:C45Clu-081, D06:C45Clu-218, D02:C45Clu-193) were further validated in the segregating population. Conclusion The combining of genome-wide introgressions and kmer genotyping strategy showed significant advantages in exploring genetic resources. Through the genome-wide comprehensive mining, a total of 11 clusters (segments) were discovered for the stable simultaneous improvement of FL/FS and LP, which should be paid more attention in the future.
... Recently, Wu et al. (2022) revealed that GhPRE1A promoted cotton fiber elongation by activating the DNAbinding bHLH factor GhPAS1, which suggested that the GhPRE1A-GhAIF3-GhPAS1 module could be located downstream of the BR signaling pathway to regulate fiber elongation. Although numerous FL QTLs and putative genes have been identified, the functional genes and molecular markers that could be utilized in MAS are still in short supply (Li et al. 2019;Feng et al. 2020;Wang et al. 2020;Song et al. 2021;Li et al. 2022;Wang et al. 2022;Chen et al. 2022;Zhang et al. 2022b). ...
Article
Full-text available
A c c e p t e d m a n u s c r i p t ACCEPTED MANUSCRIPT Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L). Author contribution statement JZ, XZ and SZ conceived and designed the experiments. SZ, HW, XL, LT, XC, and CL performed the experiments. SZ analyzed the data and wrote the paper. All authors read and approved the final version of the paper. Key message GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length.
... To date, genetic recombination between G. hirsutum and G. barbadense is commonly used to select progeny with a long FL trait, and several hundred candidate QTLs have been identified for FL (www.cottonqtldb.org) (Chandnani et al., 2018;Chen et al., 2018;He et al., 2006;Lacape et al., 2009;Li et al., 2019a;Said et al., 2015;Shi et al., 2020;Yu et al., 2013). The elongation of fibres is a subtle and complex regulatory process that occurs after initiation (À3 to 3 days post-anthesis, dpa) and lasts until 25-30 dpa (Hu et al., 2019). ...
Article
Full-text available
Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, expression new alleles (ENAs) from crossovers within genic regions and their significance in FL improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (G. hirsutum) × Hai7124 (G. barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35,420 expressed genes in developing fibres, 10,366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1,056–1,255 ENAs showed transgressive expression in fibres at 5–15 dpa (days post anthesis) of some BILs, 520 of which were located in FL‐QTLs and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL‐QTLs. Additionally, 456 extremely expressed genes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1,286, 1,089 and 1,059 expressed genes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR‐cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.
Article
Full-text available
Key message GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Abstract Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3’s expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.
Article
Full-text available
Cotton (Gossypium spp.) is the most important natural fiber source in the world. The genetic potential of cotton can be successfully and efficiently exploited by identifying and solving the complex fundamental problems of systematics, evolution, and phylogeny, based on interspecific hybridization of cotton. This study describes the results of interspecific hybridization of G. herbaceum L. (A1-genome) and G. mustelinum Miers ex Watt (AD4-genome) species, obtaining fertile hybrids through synthetic polyploidization of otherwise sterile triploid forms with colchicine (C22H25NO6) treatment. The fertile F1C hybrids were produced from five different cross combinations: (1) G. herbaceum subsp. frutescens × G. mustelinum; (2) G. herbaceum subsp. pseudoarboreum × G. mustelinum; (3) G. herbaceum subsp. pseudoarboreum f. harga × G. mustelinum; (4) G. herbaceum subsp. africanum × G. mustelinum; (5) G. herbaceum subsp. euherbaceum (variety A-833) × G. mustelinum. Cytogenetic analysis discovered normal conjugation of bivalent chromosomes in addition to univalent, open, and closed ring-shaped quadrivalent chromosomes at the stage of metaphase I in the F1C and F2C hybrids. The setting of hybrid bolls obtained as a result of these crosses ranged from 13.8–92.2%, the fertility of seeds in hybrid bolls from 9.7–16.3%, and the pollen viability rates from 36.6–63.8%. Two transgressive plants with long fiber of 35.1–37.0 mm and one plant with extra-long fiber of 39.1–41.0 mm were identified in the F2C progeny of G. herbaceum subsp. frutescens × G. mustelinum cross. Phylogenetic analysis with 72 SSR markers that detect genomic changes showed that tetraploid hybrids derived from the G. herbaceum × G. mustelinum were closer to the species G. mustelinum. The G. herbaceum subsp. frutescens was closer to the cultivated form, and its subsp. africanum was closer to the wild form. New knowledge of the interspecific hybridization and synthetic polyploidization was developed for understanding the genetic mechanisms of the evolution of tetraploid cotton during speciation. The synthetic polyploids of cotton obtained in this study would provide beneficial genes for developing new cotton varieties of the G. hirsutum species, with high-quality cotton fiber and strong tolerance to biotic or abiotic stress. In particular, the introduction of these polyploids to conventional and molecular breeding can serve as a bridge of transferring valuable genes related to high-quality fiber and stress tolerance from different cotton species to the new cultivars.
Article
A R T I C L E I N F O Keywords: Interspecific hybrid cotton Fiber heterosis MiRNA-mRNA Expression patterns in vitro ovule culture A B S T R A C T Interspecific hybridization contributes to improving cotton fiber quality, thus addressing the parental genetic diversity bottleneck problems. To elucidate the molecular mechanisms underlying fiber heterosis in interspecific hybrid cotton, an integrated analysis of phenotypes, mRNAs, and miRNAs was performed to compare the immature fibers of the parents (Gossypium hirsutum and Gossypium barbadense) with those of their hybrid at various development stages. The results showed that fiber heterosis was most obvious at 10 days post anthesis. Comprehensive analysis of mRNA and miRNA demonstrated that the main expression pattern was transgressive up-regulation (TUR) at mRNA level, while it was transgressive down-regulation (TDR) at miRNA level. The miR160-ARF and miR319-TCP pairs were mainly responsible for the changes in endogenous auxin content and secondary cell wall deposition in the fiber development stage, and they played an important role in the formation of interspecific fiber heterosis. Meanwhile, the results of indole-3-acetic acid (IAA) content determination and in vitro ovule culture indicated that auxin affected fiber heterosis. This study reveals the mechanism by which miRNA-mRNA regulates fiber heterosis, and it will promote the applications of interspecific hybrid cotton between G. hirsutum and G. barbadense.
Article
Full-text available
Key message The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Abstract Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC4F3:5). A single segment substitution line (CSSL-106) screened from BC6F2 was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC7F2 individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.
Article
Full-text available
The recombinant inbred lines of inter-specific cross, Gossypium hirsutum cv. DS-28 × G. barbadense cv. SBYF-425 was evaluated in three consecutive rainy seasons of 2017–18 (F13), 2018–19 (F14) and 2019–20 (F15) in an augmented design. The preponderance of huge continuous variability for both productivity and fiber quality traits was recorded. The principal component analysis revealed that the mapping population was well suited for mapping of productivity and fiber quality traits. On the basis the Z-scores for skewness and kurtosis, 178 RILs with normal distribution were selected for genetic linkage mapping. A high-density saturated linkage map was constructed using SNP arrays of CottonSNP63K, an Illumina’s infinium array and CottonSNP50K, CSIR-National Botanical Research Institute’s Axiom array with a total spanned length of 2402.65 cM, an average marker density of 1.54 and with map coverage of 96.99% of the reference genome. The developed genetic map of inter specific cross of Indian cotton varieties is a highly saturated in terms of coverage and highly comparable to the published maps. In QTL analysis, altogether 99 QTLs were identified for productivity and fiber quality traits. Among those, eight were stable and 38 were major QTLs. Cluster 1, 4 and 6 respectively on chromosome AD_chr.03, AD_chr.14 and AD_chr.18 were the biggest QTL clusters each with four QTLs and cluster 4 and 6 were QTL hotspots for fiber quality traits.
Article
Full-text available
Developing chromosome segments substitution lines (CSSLs) is an effective method for broadening the cotton germplasm resource, and improving the fiber quality and yield traits. In this study, the 1054 F2 individual plants and 116 F2:3 lineages were generated from the two parents of MBI9749 and MBI9915 selected from BC5F3:5 lines which originated from hybridization of CCRI36 and Hai1, and advanced backcrossing and repeated selfing. Genotypes of the parents and F2 population were analyzed. The results showed that 19 segments were introgressed for MBI9749 and 12 segments were introgressed for MBI9915, distributing on 17 linkage groups. The average background recovery rate to the recurrent parent CCRI36 was 96.70% for the two parents. An average of 16.46 segments was introgressed in F2 population. The average recovery rate of 1054 individual plants was 96.85%, and the mean length of sea island introgression segments was 157.18 cM, accounting for 3.15% of detection length. QTL mapping analysis detected 22 QTLs associated with fiber quality and yield traits in the F2 and F2:3 populations. These QTLs distributed on seven chromosomes, and the phenotypic variation was explained ranging from 1.20 to 14.61%. Four stable QTLs were detected in F2 and F2:3 populations, simultaneously. We found that eight QTLs were in agreement with the previous research. Six QTL-clusters were identified for fiber quality and yield traits, in which five QTL-clusters were on chromosome20. The results indicated that most of QTL-clusters always improve the fiber quality and have negative additive effect for yield related traits. This study demonstrated that CSSLs provide basis for fine mapping of the fiber quality and yield traits in future, and could be efficiently used for pyramiding favourable alleles to develop the new germplasms for breeding by molecular marker-assisted selection.
Article
Full-text available
MBI9915 is an introgression cotton line with excellent fiber quality. It was obtained by advanced backcrossing and continuous inbreeding from an interspecific cross between the upland cotton (Gossypium hirsutum) cultivar CCRI36 as the recurrent parent and the sea island cotton (G. barbadense) cultivar Hai1, as the donor parent. To study the genetic effects of the introgressed chromosome segments in G. hirsutum, an F2 secondary segregating population of 1537 individuals was created by crossing MBI9915 and CCRI36, and an F2:3 population was created by randomly selecting 347 individuals from the F2 generation. Quantitative trait locus (QTL) mapping and interaction for fiber length and strength were identified using IciMapping software. The genotype analysis showed that the recovery rate for MBI9915 was 97.9%, with a total 6 heterozygous segments and 13 homozygous segments. A total of 18 QTLs for fiber quality and 6 QTLs for yield related traits were detected using the two segregating generations. These QTLs were distributed across 7 chromosomes and collectively explained 0.81%–9.51% of the observed phenotypic variations. Six QTLs were consistently detected in two generations and 6 QTLs were identified in previous studies. A total of 13 pairs of interaction for fiber length and 13 pairs of interaction for fiber strength were identified in two generations. Among them, 3 pairs of interaction for fiber length and 3 pairs of interaction for fiber strength could be identified in all generations; 4 pairs of interactions affected fiber length and fiber strength simultaneously. The results clearly showed that 5 chromosome segments (Seg-5-1, Seg-7-1, Seg-8-1, Seg-20-2 and Seg-20-3) have important effects on fiber yield and quality. This study provides the useful information for gene cloning and marker-assisted breeding for excellent fiber related quality.
Article
Full-text available
As high-strength cotton fibers are critical components of high quality cotton, developing cotton cultivars with high strength fibers as well as high yield is a top priority for cotton development. Recently, chromosome segment substitution lines (CSSLs) have been developed from high-yield Upland cotton (Gossypium hirsutum) crossed with high-quality Sea Island cotton (G. barbadense). Here, we constructed a CSSL population by crossing CCRI45, a high-yield Upland cotton cultivar, with Hai1, a Sea Island cotton cultivar with superior fiber quality. We then selected two CSSLs with significantly higher fiber strength than CCRI45 (MBI7747 and MBI7561), and one CSSL with lower fiber strength than CCRI45 (MBI7285) for further analysis. We sequenced all four transcriptomes at four different time points post-anthesis, and clustered the 44,678 identified genes by function. We identified 2200 "common differentially expressed genes (DEGs)": those that were found in both high quality CSSLs (MBI7747 and MBI7561), but not in the low quality CSSL (MBI7285). Many of these genes were associated with various metabolic pathways that affect fiber strength. Upregulated DEGs were associated with polysaccharide metabolic regulation, single-organism localization, cell wall organization, and biogenesis, while the downregulated DEGs were associated with microtubule regulation, the cellular response to stress, and the cell cycle. Further analyses indicated that three genes, XLOC_036333 (mannosyl-oligosaccharide-alpha-mannosidase mns1), XLOC_029945 (FLA8), and XLOC_075372 (snakin-1) were potentially important for the regulation of cotton fiber strength. Our results suggest that these genes may be good candidates for future investigation of the molecular mechanisms of fiber strength formation, and for the improvement of cotton fiber quality through molecular breeding.
Article
Full-text available
Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p < 0.001). The genome-wide LD decayed within the genetic distance of ~30 to 32 cM at r² = 0.1, and decreased to ~1 to 2 cM at r² = 0.2, indicating the potential for association mapping. Analysis based on a mixed linear model detected 57 significant (p < 0.01) marker–trait associations, including seven associations for fiber length, ten for fiber micronaire, nine for fiber strength, eight for fiber elongation, five for fiber uniformity index, five for fiber uniformity ratio, six for boll weight and seven for lint percent, for a total of 35 SSR markers, of which 11 markers were associated with more than one trait. Among marker–trait associations, 24 associations coincided with the previously reported quantitative trait loci (QTLs), the remainder were newly identified QTLs/genes. The QTLs identified in this study will potentially facilitate improvement of fiber yield and quality in the future cotton molecular breeding programs.
Article
Full-text available
The number and location of mapped quantitative trait loci (QTL) depend on genetic populations and testing environments. The identification of consistent QTL across genetic backgrounds and environments is a pre-requisite to marker-assisted selection. This study analyzed a total of 661 abiotic and biotic stress resistance QTL based on our previous work and other publications using the meta-analysis software Biomercator. It identified chromosomal regions containing QTL clusters for different resistance traits and hotspots for a particular resistance trait in cotton from 98 QTL for drought tolerance under greenhouse (DT) and 150 QTL in field conditions (FDT), 80 QTL for salt tolerance in the greenhouse conditions (ST), 201 QTL for resistance to Verticillium wilt (VW, Verticillium dahliae), 47 QTL for resistance to Fusarium wilt (FW, Fusarium oxysporum f. sp. vasinfectum), and 85 QTL for resistance to root-knot nematodes (RKN, Meloiodogyne incognita) and reniform nematodes (RN, Rotylenchulus reniformis). The traits used in QTL mapping for abiotic stress tolerance included morphological traits—plant height and fresh and dry shoot and root weights, physiological traits—chlorophyll content, osmotic potential, carbon isotope ratio, stomatal conductance, photosynthetic rate, transpiration, canopy temperature, and leaf area index, agronomic traits—seedcotton yield, lint yield, boll weight, and lint percent, and fiber quality traits—fiber length, uniformity, strength, elongation, and micronaire. The results showed that resistance QTL are not uniformly distributed across the cotton genome; some chromosomes carried disproportionally more QTL, QTL clusters, or hotspots. Twenty-three QTL clusters were found on 15 chromosomes (c3, c4, c5, c6, c7, c11, c14, c15, c16, c19, c20, c23, c24, c25, and c26). Moreover, 28 QTL hotshots were associated with different resistance traits including one hotspot on c4 for Verticillium wilt resistance, two QTL hotspots on c24 for chlorophyll content measured under both drought and salt stress conditions, and three other hotspots on c19 for the resistance to Verticillium wilt and Fusarium wilt, and micronaire under drought stress conditions. This meta-analysis of stress tolerance QTL provides an important foundation for cotton breeding and further studies on the genetic mechanisms of abiotic and biotic stress resistance in cotton. https://link.springer.com/epdf/10.1007/s00438-017-1342-0?author_access_token=ZZFramuBFPe8LScq2b6VOPe4RwlQNchNByi7wbcMAY6z78kZ3WHMaKL_EdwvA873sMqAagBy3Cf-UuQJd5jZSnKJMpXYRzAnG8cAVNyvr34n3sWrIRwky_Lzng2H1yupqa-2K2zgglIoFrvvmNDvGQ==
Article
Full-text available
Two immortalized backcross populations (DHBCF1s and JMBCF1s) were developed using a recombinant inbred line (RIL) population crossed with the two parents DH962 and Jimian5 (as the males), respectively. The fiber quality and yield component traits of the two backcross populations were phenotyped at four environments (two locations, two years). One hundred seventy-eight quantitative trait loci (QTL) were detected including 76 for fiber qualities and 102 for yield components, explaining 4.08–17.79% of the phenotypic variation (PV). Among the 178 QTL, 22 stable QTL were detected in more than one environment or population. A stable QTL, qFL-c10-1, was detected in the previous F2 population, a RIL population in 3 environments and the current two BCF1 populations in this study, explaining 5.79–37.09% of the PV. Additionally, 117 and 110 main-effect QTL (M-QTL) and 47 and 191 digenic epistatic QTL (E-QTL) were detected in the DHBCF1s and JMBCF1s populations, respectively. The effect of digenic epistasis played a more important role on lint percentage, fiber length and fiber strength. These results obtained in the present study provided more resources to obtain stable QTL, confirming the authenticity and reliability of the QTL for molecular marker-assisted selection breeding and QTL cloning.
Article
Full-text available
Background Gossypium hirsutum L., or upland cotton, is an important renewable resource for textile fiber. To enhance understanding of the genetic basis of cotton earliness, we constructed an intra-specific recombinant inbred line population (RIL) containing 137 lines, and performed linkage map construction and quantitative trait locus (QTL) mapping. ResultsUsing restriction-site associated DNA sequencing, a genetic map composed of 6,434 loci, including 6,295 single nucleotide polymorphisms and 139 simple sequence repeat loci, was developed from RIL population. This map spanned 4,071.98 cM, with an average distance of 0.63 cM between adjacent markers. A total of 247 QTLs for six earliness-related traits were detected in 6 consecutive years. In addition, 55 QTL coincidence regions representing more than 60 % of total QTLs were found on 22 chromosomes, which indicated that several earliness-related traits might be simultaneously improved. Fine-mapping of a 2-Mb region on chromosome D3 associated with five stable QTLs between Marker25958 and Marker25963 revealed that lines containing alleles derived from CCRI36 in this region exhibited smaller phenotypes and earlier maturity. One candidate gene (EMF2) was predicted and validated by quantitative real-time PCR in early-, medium- and late-maturing cultivars from 3- to 6-leaf stages, with highest expression level in early-maturing cultivar, CCRI74, lowest expression level in late-maturing cultivar, Bomian1. Conclusions We developed an SNP-based genetic map, and this map is the first high-density genetic map for short-season cotton and has the potential to provide deeper insights into earliness. Cotton earliness-related QTLs and QTL coincidence regions will provide useful materials for QTL fine mapping, gene positional cloning and MAS. And the gene, EMF2, is promising for further study.
Article
Full-text available
Background Verticillium wilt (VW) caused by Verticillium dahliae (Kleb) is one of the most destructive diseases of cotton. The identification of highly resistant QTLs or genes in the whole cotton genome is quite important for developing a VW-resistant variety and for further molecular design breeding. Results In the present study, BC1F1, BC1S1, and BC2F1 populations derived from an interspecific backcross between the highly resistant line Hai1 (Gossypium barbadense L.) and the susceptible variety CCRI36 (G. hirsutum L.) as the recurrent parent were constructed. Quantitative trait loci (QTL) related to VW resistance were detected in the whole cotton genome using a high-density simple sequence repeat (SSR) genetic linkage map from the BC1F1 population, with 2292 loci covering 5115.16 centiMorgan (cM) of the cotton (AD) genome, and the data concerning VW resistance that were obtained from four dates of BC2F1 in the artificial disease nursery and one date of BC1S1 and BC2F1 in the field. A total of 48 QTLs for VW resistance were identified, and 37 of these QTLs had positive additive effects, which indicated that the G. barbadense alleles increased resistance to VW and decreased the disease index (DI) by about 2.2–10.7. These QTLs were located on 19 chromosomes, in which 33 in the A subgenome and 15 QTLs in the D subgenome. The 6 QTLs were found to be stable. The 6 QTLs were consistent with those identified previously, and another 42 were new, unreported QTLs, of which 31 QTLs were from G. barbadense. By meta-analysis, 17 QTL hotspot regions were identified and 10 of them were new, unreported hotspot regions. 29 QTLs in this paper were in 12 hotspot regions and were all from G. barbadense. Conclusions These stable or consensus QTL regions warrant further investigation to better understand the genetics and molecular mechanisms underlying VW resistance. This study provides useful information for further comparative analysis and marker-assisted selection in the breeding of disease-resistant cotton. It may also lay an important foundation for gene cloning and further molecular design breeding for the entire cotton genome. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3128-x) contains supplementary material, which is available to authorized users.
Article
We employ a detailed restriction fragment length polymorphism (RFLP) map to investigate chromosome organization and evolution in cotton, a disomic polyploid. About 46.2% of nuclear DNA probes detect RFLPs distinguishing Gossypium hirsutum and Gossypium barbadense; and 705 RFLP loci are assembled into 41 linkage groups and 4675 cM. The subgenomic origin (A vs. D) of most, and chromosomal identity of 14 (of 26), linkage groups is shown. The A and D subgenomes show similar recombinational length, suggesting that repetitive DNA in the physically larger A subgenome is recombinationally inert. RFLPs are somewhat more abundant in the D subgenome. Linkage among duplicated RFLPs reveals 11 pairs of homoelogous chromosomal regions-two appear homosequential, most differ by inversions, and at least one differs by a translocation. Most homoeologies involve chromosomes from different subgenomes, putatively reflecting the n = 13 to n = 26 polyploidization event of 1.1-1.9 million years ago. Several observations suggest that another, earlier, polyploidization event spawned n = 13 cottons, at least 25 million years ago. The cotton genome contains about 400-kb DNA per cM, hence map-based gene cloning is feasible. The cotton map affords new opportunities to study chromosome evolution, and to exploit Gossypium genetic resources for improvement of the world's leading natural fiber.