Fiber Development Genes Research Articles

Transcriptome analysis combined with metabolome analysis reveals the significant functions of CesA genes in cotton (Gossypium hirsutum) fiber length development

Cotton is widely distributed worldwide, and improving the quality of its fiber is one of the most important tasks in cotton breeding. Cotton fibers are primarily composed of cellulose, which is synthesized by CesA complexes (CSCs). However, the functions of CesA genes in cotton fiber development have not been comprehensively analysed. In this study, the cotton transcriptome and metabolome were used to investigate the function of CesA genes in fiber development. Finally, 321 metabolites were obtained, 84 of which were associated with the corresponding genes. Interestingly, a target gene named Gh_A08G144300, one of the CesA gene family members, was closely correlated with the development of cotton fibers. The target CesA gene Gh_A08G144300 was analysed to determine its specific function in cotton fiber development. High-level gene expression of Gh_A08G144300 was found at different fiber development stages by RNA-seq analysis, and the silencing of Gh_A08G144300 visibly inhibited the growth of cotton fibers, showing that it is critical for their growth. This study provides an important reference for research on the gene function of Gh_A08G144300 and the regulatory mechanism of fiber development in cotton.

Genome-wide identification, characterization, and expression analysis of aluminum-activated malate transporter genes (ALMTs) in Gossypium hirsutum L.

Aluminum-activated malate transporters (ALMT) are widely involved in plant growth and metabolic processes, including adaptation to acid soils, guard cell regulation, anion homeostasis, and seed development. Although ALMT genes have been identified in Arabidopsis, wheat, barley, and Lotus japonicus, little is known about its presence in Gossypium hirsutum L. In this study, ALMT gene recognition in diploid and tetraploid cotton were done using bioinformatics analysis that examined correlation between homology and evolution. Differentially regulated ALMT genetic profile in G. hirsutum was examined, using RNA sequencing and qRT-PCR, during six fiber developmental time-points, namely 5 d, 7 d, 10 d, 15 d, 20 d, and 25 d. We detected 36 ALMT genes in G. hirsutum, which were subsequently annotated and divided into seven sub-categories. Among these ALMT genes, 34 had uneven distribution across 14/26 chromosomes. Conserved domains and gene structure analysis indicated that ALMT genes were highly conserved and composed of exons and introns. The GhALMT gene expression profile at different DPA (days post anthesis) in different varieties of G. hirsutum is indicative of a crucial role of ALMT genes in fiber development in G. hirsutum. This study provides basis for advancements in the cloning and functional enhancements of ALMT genes in enhancing fiber development and augmenting high quality crop production.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes

BackgroundDirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton.ResultsIn this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls.ConclusionsThese findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

The Ligon lintless -2 Short Fiber Mutation Is Located within a Terminal Deletion of Chromosome 18 in Cotton.

Extreme elongation distinguishes about one-fourth of cotton (Gossypium sp.) seed epidermal cells as "lint" fibers, useful for the textile industry, from "fuzz" fibers (<5 mm). Ligon lintless-2 (Li 2 ), a dominant mutation that results in no lint fiber but normal fuzz fiber, offers insight into pathways and mechanisms that differentiate spinnable cotton from its progenitors. A genetic map developed using 1,545 F2 plants showed that marker CISP15 was 0.4 cM from Li 2 , and "dominant" simple sequence repeat (SSR) markers (i.e. with null alleles in the Li 2 genotype) SSR7 and SSR18 showed complete linkage with Li 2 Nonrandom distribution of markers with null alleles suggests that the Li 2 phenotype results from a 176- to 221-kb deletion of the terminal region of chromosome 18 that may have been masked in prior pooled-sample mapping strategies. The deletion includes 10 genes with putative roles in fiber development. Two Glycosyltransferase Family 1 genes showed striking expression differences during elongation of wild-type versus Li 2 fiber, and virus-induced silencing of these genes in the wild type induced Li 2 -like phenotypes. Further, at least 7 of the 10 putative fiber development genes in the deletion region showed higher expression in the wild type than in Li 2 mutants during fiber development stages, suggesting coordinated regulation of processes in cell wall development and cell elongation, consistent with the hypothesis that some fiber-related quantitative trait loci comprise closely spaced groups of functionally diverse but coordinately regulated genes.

Genome-wide analysis of cotton C2H2-zinc finger transcription factor family and their expression analysis during fiber development

BackgroundC2H2-zinc finger protein family is commonly found in the plant, and it is known as the key actors in the regulation of transcription and vital component of chromatin structure. A large number of the C2H2-zinc finger gene members have not been well characterized based on their functions and structure in cotton. However, in other plants, only a few C2H2-zinc finger genes have been studied.ResultsIn this work, we performed a comprehensive analysis and identified 386, 196 and 195 C2H2-zinc finger genes in Gossypium hirsutum (upland cotton), Gossypium arboreum and Gossypium raimondii, respectively. Phylogenetic tree analysis of the C2H2-zinc finger proteins encoding the C2H2-zinc finger genes were classified into seven (7) subgroups. Moreover, the C2H2-zinc finger gene members were distributed in all cotton chromosomes though with asymmetrical distribution patterns. All the orthologous genes were detected between tetraploid and the diploid cotton, with 154 orthologous genes pair detected between upland cotton and Gossypium arboreum while 165 orthologous genes were found between upland cotton and Gossypium raimondii. Synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates (Ka/Ks) analysis indicated that the cotton C2H2-zinc finger genes were highly influenced mainly by negative selection, which maintained their protein levels after the duplication events. RNA-seq data and RT-qPCR validation of the RNA seq result revealed differential expression pattern of some the C2H2-zinc finger genes at different stages of cotton fiber development, an indication that the C2H2-zinc finger genes play an important role in initiating and regulating fiber development in cotton.ConclusionsThis study provides a strong foundation for future practical genome research on C2H2-zinc finger genes in upland cotton. The expression levels of C2H2-zinc finger genes family is a pointer of their involvement in various biochemical and physiological functions which are directly related to cotton fiber development during initiation and elongation stages. This work not only provides a basis for determining the nominal role of the C2H2-zinc finger genes in fiber development but also provide valuable information for characterization of potential candidate genes involved in regulation of cotton fiber development.

Genomic analysis of marker-associated fiber development genes in upland cotton (Gossypium hirsutum L)

Understanding genetic control of fiber development is essential for cotton improvement. This study analyzed several thousand fiber Unigenes in Upland cotton (Gossypium hirsutum). They included 2423 Unigenes for fiber initiation (IN), 1841 for fiber elongation (EL) and 86 for secondary cell wall deposition (SCWD). These Unigenes were located on recombination hotspots between At and Dt subgenomes with nearly equal frequency (51% vs. 49%). Sequence analysis with Upland cotton TM-1 genome assembly showed that 477 IN, 464 EL, 20 SCWD and 367 non-fiber Unigenes intersect with predicted coding sequences. There were 140 IN, 118 EL and one SCWD Unigenes that do not overlap between individual stages of fiber development. Annotation of top scoring Unigenes with Kyoto encyclopedia of genes and genomes indicated mitochondrial biogenesis for IN Unigene #18205, conversion of PRPP to histidine for IN #25329, and spliceosome pathways for IN #24380. EL Unigene #11669 appeared to be responsible for septum site-determining protein, EL #11916 for vacuolar protein, and EL #17251 for protein trichome birefringence-like 33. One SCWD Unigene #190099 was responsible for legumin A-like protein. In addition, Unigene #03070 overlapped three fiber development stages and involved in fatty acid biosynthesis. Unigenes that shared same genomic positions for all fiber stages are of particular interest as they may have possible pleiotropic function. On the other hand, Unigenes from one fiber development stage that did not share with any other stage may have stage-specific function. This information would be useful for plant breeders to exploit fiber genes in cotton breeding programs.

Synchronous Transcription of Cytoskeleton-Associated Genes is Critical to Cotton Fiber Elongation

Cytoskeleton-associated genes’ expression in the modern cotton fibers has selectively been evolved under domestication. Since the magnitude of transcriptional dynamicity of cytoskeletal genes in fibers has not yet been thoroughly explored, comprehensive expression analyses of the domestication-driven profilins (GhPRFs) were performed by assaying 12 tissues from various developmental stages of elite cotton genotypes. Results from 3240 successful profilin (GhPRF1–GhPRF6) expression assays revealed enormous spatial transcriptional biases with their predominant expression in the elongating fibers. Moreover, fiber developmental stages have more pronounced impacts on differential gene expression than the genetic backgrounds. Transcriptional comparisons of GhPRFs-associated key cytoskeletal genes GhAnnex and GhACT show perfect (r = 1.0) to strong positive (r > 0.95) correlation in the developing fibers. RNAi-mediated constitutive suppression of domestication-driven GhPRF1 gene resulted in significant expression modulation of GhAnnex and GhACT genes in floral/fiber tissues and consequently sustained the cell wall development. These results highlight that temporal expression of GhPRF–GhAnnex–GhACT triad is exceptionally coordinated which illustrates fascinating crosstalk in diverse cellular contexts and cell types. Temporal fold-expression variations of other cytoskeletal-‘structural’ genes Tub–Exp–FLA triad were transversely synchronous to the GhPRF–GhAnnex–GhACT genic triad. However, the CMW-‘modifying’ genes GalT1, ABP1, and CesA were unconditionally divergent from the synchronous temporal dynamics of genic triad investigated beforehand. To our knowledge, this study represents the first in-depth description of synchronous transcriptional dynamicity of cytoskeletal-‘structural’ genes in the developing fibers. These results present a novel viewpoint on the temporal expression evolution of cytoskeletal-‘structural’ genes, augmenting our understanding of the evolutionary genetics of cytoskeletal genes in fiber development.

Genome-wide SSR-based association mapping for fiber quality in nation-wide upland cotton inbreed cultivars in China.

BackgroundSince upland cotton was introduced into China during the 1920s–1950s, hundreds of inbreed cultivars have been developed. To explore the molecular diversity, population structure and elite alleles, 503 inbred cultivars developed in China and some foreign cultivars from the United States and the Soviet Union were collected and analyzed by 494 genome-wide SSRs (Simple Sequence Repeats).MethodsFour hundred and ninety-four pairs of SSRs with high polymorphism and uniform distribution on 26 chromosomes were used to scan polymorphisms in 503 nation-wide upland cottons. The programming language R was used to make boxplots for the phenotypic traits in different environments. Molecular marker data and 6 fiber quality traits were analyzed by the method of MLM (mixed linear model) (P + G + Q + K) in the TASSEL software package on the basis of the population structure and linkage disequilibrium analysis. The loci of elite allelic variation and typical materials carrying elite alleles were identified based on phenotypic effect values.ResultsA total of 179 markers were polymorphic and generated 426 allele loci; the population based on molecular diversity was classified into seven subpopulations corresponding to pedigree origin, ecological and geographical distribution. The attenuation distance of linkage disequilibrium dropped significantly up to 0–5 cM. Association mapping for fiber quality showed that 216 marker loci were associated with fiber quality traits (P < 0.05) explaining 0.58 % ~ 5.12 % of the phenotypic variation, with an average of 2.70 %. Thirteen marker loci were coincident with other studies, and three were detected for the same trait. Seven quantitative trait loci were related to known genes in fiber development. Based on phenotypic effects, 48 typical materials that contained the elite allele loci related to fiber quality traits were identified and are widely used in practical breeding.ConclusionsThe molecular diversity and population structure of 503 nation-wide upland cottons in China were evaluated by 494 genome-wide SSRs, and association mapping for fiber quality revealed known and novel elite alleles. The molecular diversity provides a guide for parental mating in cotton breeding, and the association mapping results will aid in the fine-mapping genes related to fiber quality traits and facilitate further studies on candidate genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2662-x) contains supplementary material, which is available to authorized users.

Developing fiber specific promoter-reporter transgenic lines to study the effect of abiotic stresses on fiber development in cotton.

Cotton is one of the most important cash crops in US agricultural industry. Environmental stresses, such as drought, high temperature and combination of both, not only reduce the overall growth of cotton plants, but also greatly decrease cotton lint yield and fiber quality. The impact of environmental stresses on fiber development is poorly understood due to technical difficulties associated with the study of developing fiber tissues and lack of genetic materials to study fiber development. To address this important question and provide the need for scientific community, we have generated transgenic cotton lines harboring cotton fiber specific promoter (CFSP)-reporter constructs from six cotton fiber specific genes (Expansin, E6, Rac13, CelA1, LTP, and Fb late), representing genes that are expressed at different stages of fiber development. Individual CFSP::GUS or CFSP::GFP construct was introduced into Coker 312 via Agrobacterium mediated transformation. Transgenic cotton lines were evaluated phenotypically and screened for the presence of selectable marker, reporter gene expression, and insertion numbers. Quantitative analysis showed that the patterns of GUS reporter gene activity during fiber development in transgenic cotton lines were similar to those of the native genes. Greenhouse drought and heat stress study showed a correlation between the decrease in promoter activities and decrease in fiber length, increase in micronaire and changes in other fiber quality traits in transgenic lines grown under stressed condition. These newly developed materials provide new molecular tools for studying the effects of abiotic stresses on fiber development and may be used in study of cotton fiber development genes and eventually in the genetic manipulation of fiber quality.

Distribution and evolution of cotton fiber development genes in the fibreless Gossypium raimondii genome

Cotton fiber represents the largest single cell in plants and they serve as models to study cell development. This study investigated the distribution and evolution of fiber Unigenes anchored to recombination hotspots between tetraploid cotton (Gossypium hirsutum) At and Dt subgenomes, and within a parental diploid cotton (Gossypium raimondii) D genome. Comparative analysis of At vs D and Dt vs D showed that 1) the D genome provides many fiber genes after its merger with another parental diploid cotton (Gossypium arboreum) A genome although the D genome itself does not produce any spinnable fiber; 2) similarity of fiber genes is higher between At vs D than between Dt vs D genomic hotspots. This is the first report that fiber genes have higher similarity between At and D than between Dt and D. The finding provides new insights into cotton genomic regions that would facilitate genetic improvement of natural fiber properties.

Comparison of the transcriptome between two cotton lines of different fiber color and quality.

To understand the mechanism of fiber development and pigmentation formation, the mRNAs of two cotton lines were sequenced: line Z128 (light brown fiber) was a selected mutant from line Z263 (dark brown fiber). The primary walls of the fiber cell in both Z263 and Z128 contain pigments; more pigments were laid in the lumen of the fiber cell in Z263 compared with that in Z128. However, Z263 contained less cellulose than Z128. A total of 71,895 unigenes were generated: 13,278 (20.26%) unigenes were defined as differentially expressed genes (DEGs) by comparing the library of Z128 with that of Z263; 5,345 (8.16%) unigenes were up-regulated and 7,933 (12.10%) unigenes were down-regulated. qRT-PCR and comparative transcriptional analysis demonstrated that the pigmentation formation in brown cotton fiber was possibly the consequence of an interaction between oxidized tannins and glycosylated anthocyanins. Furthermore, our results showed the pigmentation related genes not only regulated the fiber color but also influenced the fiber quality at the fiber elongation stage (10 DPA). The highly expressed flavonoid gene in the fiber elongation stage could be related to the fiber quality. DEGs analyses also revealed that transcript levels of some fiber development genes (Ca2+/CaM, reactive oxygen, ethylene and sucrose phosphate synthase) varied dramatically between these two cotton lines.

Structure, expression differentiation and evolution of duplicated fiber developmental genes in Gossypium barbadense and G. hirsutum.

BackgroundBoth Gossypium hirsutum and G. barbadense probably originated from a common ancestor, but they have very different agronomic and fiber quality characters. Here we selected 17 fiber development-related genes to study their structures, tree topologies, chromosomal location and expression patterns to better understand the interspecific divergence of fiber development genes in the two cultivated tetraploid species.ResultsThe sequence and structure of 70.59% genes were conserved with the same exon length and numbers in different species, while 29.41% genes showed diversity. There were 15 genes showing independent evolution between the A- and D-subgenomes after polyploid formation, while two evolved via different degrees of colonization. Chromosomal location showed that 22 duplicate genes were located in which at least one fiber quality QTL was detected. The molecular evolutionary rates suggested that the D-subgenome of the allotetraploid underwent rapid evolutionary differentiation, and selection had acted at the tetraploid level. Expression profiles at fiber initiation and early elongation showed that the transcripts levels of most genes were higher in Hai7124 than in TM-1. During the primary-secondary transition period, expression of most genes peaked earlier in TM-1 than in Hai7124. Homeolog expression profile showed that A-subgenome, or the combination of A- and D-subgenomes, played critical roles in fiber quality divergence of G. hirsutum and G. barbadense. However, the expression of D-subgenome alone also played an important role.ConclusionIntegrating analysis of the structure and expression to fiber development genes, suggests selective breeding for certain desirable fiber qualities played an important role in divergence of G. hirsutum and G. barbadense.

Polyploidization Altered Gene Functions in Cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

Characterization of 19 novel cotton FLA genes and their expression profiling in fiber development and in response to phytohormones and salt stress

Fasciclin-like arabinogalactan proteins (FLAs), a subclass of arabinogalactan proteins (AGPs), are usually involved in cell development in plants. To investigate the expression profiling as well as the role of FLA genes in fiber development, 19 GhFLA genes (cDNAs) were isolated from cotton (Gossypium hirsutum). Among them, 15 are predicted to be glycosylphosphatidylinositol anchored to the plasma membranes. The isolated cotton FLAs could be divided into four groups. Real-time quantitative reverse transcriptase polymerase chain reaction results indicated that the GhFLA genes are differentially expressed in cotton tissues. Three genes (GhFLA1/2/4) were specifically or predominantly expressed in 10 days post-anthesis fibers, and the transcripts of the other four genes (GhFLA6/14/15/18) were accumulated at relatively high levels in cotton fibers. Furthermore, expressions of the GhFLA genes are regulated in fiber development and in response to phytohormones and NaCl. The identification of cotton FLAs will facilitate the study of their roles in cotton fiber development and cell wall biogenesis.

Cotton (Gossypium spp.) R2R3-MYB transcription factors SNP identification, phylogenomic characterization, chromosome localization, and linkage mapping

R2R3-MYB transcription factors of plants are involved in the regulation of trichome length and density. Several of them are differentially expressed during initiation and elongation of cotton fibers. We report sequence phylogenomic characterization of the six MYB genes, their chromosomal localization, and linkage mapping via SNP marker in AD-genome cotton (2n = 52). Phylogenetic grouping and comparison to At- and Dt-genome putative ancestral diploid species of allotetraploid cotton facilitated differentiation between genome-specific polymorphisms (GSPs) and marker-suitable locus-specific polymorphisms (LSPs). The SNP frequency averaged one per 77 bases overall, and one per 106 and 30 bases in coding and non-coding regions, respectively. SNP-based multivariate relationships conformed to independent evolution of the six MYB homoeologous loci in the four tetraploid species. Nucleotide diversity analysis indicated that the six MYB loci evolved more quickly in the Dt- than At-genome. The greater variation in the Dt-D genome comparisons than that in At-A genome comparisons showed no significant bias among synonymous substitution, non-synonymous substitution, and nucleotide change in non-coding regions. SNPs were concordantly mapped by deletion analysis and linkage mapping, which confirmed their value as candidate gene markers and indicated the reliability of the SNP discovery strategy in tetraploid cotton species. We consider that these SNPs may be useful for genetic dissection of economically important fiber and yield traits because of the role of these genes in fiber development.

Molecular cloning and characterization of the actin-depolymerizing factor gene in Gossypium barbadense

The Sea Island cotton (Gossypium barbadense L.) has been highly valued in verticillium wilt resistance and many fiber qualities including fiber length, strength and fineness. To identify whether it had some special genes in fiber development in comparison with the Upland cotton (G. hirsutum L.), an actin-depolymerizing factor (ADF) gene was cloned and characterized in this research. A 420 bp open reading frame of the cloned gene, named GbADF1, encoded a protein of 139 amino acids, including 39.57% nonpolar amino acids, 17.27% acidic amino acids, 15.83% basic amino acids and 31.92% hydrophobic amino acids. Its molecular weight was about 15 kDa, and pI 5.04. GbADF1 contained two conserved domains, 6-Ser and the PIP2/actin binding site. Its amino acid sequence was similar to the ADF/cofilin family of other plants. Compared with cDNA sequence, the GbADF1 gene contained one intron near the 3' end in genomic sequence. Semi-quantitative RT-PCR result showed that GbADF1 was a constitutive expression gene in cotton, and higher expression level was detected in fibers than in trophic tissues. The GbADF1 was successfully expressed as a fusion protein in Escherichia coli BL21 (DE3). The molecular weight was firstly calculated by SDS-PAGE. Western blotting analysis confirmed the existence of a protein corresponding to GbADF1. The structure of GbADF1 was different from that of other ADF genes in higher plant, although the coding sequences of all cloned ADFs were highly conserved.

Microsatellite markers associated with lint percentage trait in cotton, Gossypium hirsutum

Molecular markers associated with fiber development traits have the potential to play a key role in understanding of cotton fiber development. Seventeen SSRs out of 304 markers tested from MGHES (EST-SSR), JESPR (genomic SSR), and TMB (BAC-derived SSR) collections showed significant linkage associations (using a Kurskal-Wallis non-parametric test) with lint percentage QTL in a set of recombinant inbred cotton lines (RILs) segregating for lint percentage. The permutation test of these potential markers associated with lint percentage QTL(s) determined that 12 SSR markers have stable estimates, exceeding empirically chosen threshold significance values at or above α = 0.01. Interval mapping demonstrated that 9 SSRs with stable critical LOD threshold values at α = 0.01 have significant QTL effect. Multiple QTL-mapping (MQM) revealed that at least, two highly significant fiber development QTLs exist around regions TMB0471 and MGHES–31 (explained about 23–59% of the phenotypic variation of lint percentage) and around markers MGHES–31 and TMB0366 (accounted for 5.4–12.5% phenotypic variation of lint percentage). These markers, in particular fiber-specific EST-SSRs, might be the possible ‘candidate’ loci contributing for fiber development in cotton. BAC-derived SSRs associated with fiber trait are the possible markers that are useful for the identification of physical genomic contigs that contain fiber development genes. Several lint percentage trait associated SSR markers have been located to chromosomes 12, 18, 23, and 26 using deletion analysis in aneuploid chromosome substitution lines. Outcomes of the work may prove useful in understanding and revealing the molecular basis of the fiber development, and the utilization of these markers for development of superior cotton cultivars through marker-assisted selection (MAS) programs.