Analysis of ESTs from multiple Gossypium hirsutum tissues and identification of SSRs

Comparative Transcriptomic Analysis of Biological Process and Key Pathway in Three Cotton (Gossypium spp.) Species Under Drought Stress

Drought is one of the most important abiotic stresses that seriously affects cotton growth, development, and production worldwide. However, the molecular mechanism, key pathway, and responsible genes for drought tolerance incotton have not been stated clearly. In this research, high-throughput next generation sequencing technique was utilized to investigate gene expression profiles of three cotton species (Gossypium hirsutum, Gossypium arboreum, and Gossypium barbadense L.) under drought stress. A total of 6968 differentially expressed genes (DEGs) were identified, where 2053, 742, and 4173 genes were tested as statistically significant; 648, 320, and 1998 genes were up-regulated, and 1405, 422, and 2175 were down-regulated in TM-1, Zhongmian-16, and Pima4-S, respectively. Total DEGs were annotated and classified into functional groups under gene ontology analysis. The biological process was present only in tolerant species(TM-1), indicating drought tolerance condition. The Kyoto encyclopedia of genes and genomes showed the involvement of plant hormone signal transduction and metabolic pathways enrichment under drought stress. Several transcription factors associated with ethylene-responsive genes (ICE1, MYB44, FAMA, etc.) were identified as playing key roles in acclimatizing to drought stress. Drought also caused significant changes in the expression of certain functional genes linked to abscisic acid (ABA) responses (NCED, PYL, PP2C, and SRK2E), reactive oxygen species (ROS) related in small heat shock protein and 18.1 kDa I heat shock protein, YLS3, and ODORANT1 genes. These results will provide deeper insights into the molecular mechanisms of drought stress adaptation in cotton.

Re-evaluation of the inheritance for root-knot nematode resistance in the Upland cotton germplasm line M-120 RNR revealed two epistatic QTLs conferring resistance

We report a second major QTL for root-knot nematode resistance in the highly resistant Upland cotton line M-120RNR and show epistasis between two resistant QTLs with different mechanisms conferring resistance. In an earlier study, we identified a major QTL on Chromosome 11 associated with resistance to root-knot nematode in the M-120 RNR Upland cotton line (Gossypium hirsutum L.) of the Auburn 623 RNR source. Herein, we re-evaluated the genetics of the resistance to root-knot nematode in the M-120 RNR × Pima S-6 population by linkage mapping using recently published SSR markers. The QTL analysis detected two regions significantly associated with the resistance phenotype. In addition to the QTL previously identified on Chromosome 11 (qMi-C11), a major QTL was identified on Chromosome 14 (qMi-C14). The resistance locus on qMi-C11 originated from the Clevewilt parent, while the qMi-C14 locus originated from the other resistant parent, Mexico Wild Jack Jones. The qMi-C14 locus had logarithms of odds score of 17 and accounted for 45 % of the total phenotype variation in egg production. It was also associated with galling index, but the percent variation explained was only 6 %, suggesting that the qMi-C11 locus had a much stronger effect on root gall suppression than egg production, while the qMi-C14 locus had a stronger effect on egg production than galling. The results also suggest that the transgressive segregation observed in the development of Auburn 623 RNR was due to the pyramiding of at least two main effect QTLs as well as an additive-by-additive epistatic effects between the two resistant loci. The SSRs markers tightly linked to the qMi-C11 and qMi-C14 loci will greatly facilitate the improvement of RKN resistance in cotton via marker-assisted breeding.

Genome-wide functional analysis of cotton (Gossypium hirsutum) in response to drought

Cotton is one of the most important crops for its natural textile fibers in the world. However, it often suffered from drought stress during its growth and development, resulting in a drastic reduction in cotton productivity. Therefore, study on molecular mechanism of cotton drought-tolerance is very important for increasing cotton production. To investigate molecular mechanism of cotton drought-resistance, we employed RNA-Seq technology to identify differentially expressed genes in the leaves of two different cultivars (drought-resistant cultivar J-13 and drought-sensitive cultivar Lu-6) of cotton. The results indicated that there are about 13.38% to 18.75% of all the unigenes differentially expressed in drought-resistant sample and drought-sensitive control, and the number of differentially expressed genes was increased along with prolonged drought treatment. DEG (differentially expression gene) analysis showed that the normal biophysical profiles of cotton (cultivar J-13) were affected by drought stress, and some cellular metabolic processes (including photosynthesis) were inhibited in cotton under drought conditions. Furthermore, the experimental data revealed that there were significant differences in expression levels of the genes related to abscisic acid signaling, ethylene signaling and jasmonic acid signaling pathways between drought-resistant cultivar J-13 and drought-sensitive cultivar Lu-6, implying that these signaling pathways may participate in cotton response and tolerance to drought stress.

Generation and analysis of a large-scale expressed sequence Tag database from a full-length enriched cDNA library of developing leaves of Gossypium hirsutum L

Background

Cotton (Gossypium hirsutum L.) is one of the world’s most economically-important crops. However, its entire genome has not been sequenced, and limited resources are available in GenBank for understanding the molecular mechanisms underlying leaf development and senescence.

Methodology/Principal Findings

In this study, 9,874 high-quality ESTs were generated from a normalized, full-length cDNA library derived from pooled RNA isolated from throughout leaf development during the plant blooming stage. After clustering and assembly of these ESTs, 5,191 unique sequences, representative 1,652 contigs and 3,539 singletons, were obtained. The average unique sequence length was 682 bp. Annotation of these unique sequences revealed that 84.4% showed significant homology to sequences in the NCBI non-redundant protein database, and 57.3% had significant hits to known proteins in the Swiss-Prot database. Comparative analysis indicated that our library added 2,400 ESTs and 991 unique sequences to those known for cotton. The unigenes were functionally characterized by gene ontology annotation. We identified 1,339 and 200 unigenes as potential leaf senescence-related genes and transcription factors, respectively. Moreover, nine genes related to leaf senescence and eleven MYB transcription factors were randomly selected for quantitative real-time PCR (qRT-PCR), which revealed that these genes were regulated differentially during senescence. The qRT-PCR for three GhYLSs revealed that these genes express express preferentially in senescent leaves.

Conclusions/Significance

These EST resources will provide valuable sequence information for gene expression profiling analyses and functional genomics studies to elucidate their roles, as well as for studying the mechanisms of leaf development and senescence in cotton and discovering candidate genes related to important agronomic traits of cotton. These data will also facilitate future whole-genome sequence assembly and annotation in G. hirsutum and comparative genomics among Gossypium species.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Genome-wide identification of differentially expressed genes under water deficit stress in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cotton is the world's primary fiber crop and is a major agricultural commodity in over 30 countries. Like many other global commodities, sustainable cotton production is challenged by restricted natural resources. In response to the anticipated increase of agricultural water demand, a major research direction involves developing crops that use less water or that use water more efficiently. In this study, our objective was to identify differentially expressed genes in response to water deficit stress in cotton. A global expression analysis using cDNA-Amplified Fragment Length Polymorphism was conducted to compare root and leaf gene expression profiles from a putative drought resistant cotton cultivar grown under water deficit stressed and well watered field conditions.

RESULTS

We identified a total of 519 differentially expressed transcript derived fragments. Of these, 147 transcript derived fragment sequences were functionally annotated according to their gene ontology. Nearly 70 percent of transcript derived fragments belonged to four major categories: 1) unclassified, 2) stress/defense, 3) metabolism, and 4) gene regulation. We found heat shock protein-related and reactive oxygen species-related transcript derived fragments to be among the major parts of functional pathways induced by water deficit stress. Also, twelve novel transcripts were identified as both water deficit responsive and cotton specific. A subset of differentially expressed transcript derived fragments was verified using reverse transcription-polymerase chain reaction. Differential expression analysis also identified five pairs of duplicated transcript derived fragments in which four pairs responded differentially between each of their two homologues under water deficit stress.

CONCLUSIONS

In this study, we detected differentially expressed transcript derived fragments from water deficit stressed root and leaf tissues in tetraploid cotton and provided their gene ontology, functional/biological distribution, and possible roles of gene duplication. This discovery demonstrates complex mechanisms involved with polyploid cotton's transcriptome response to naturally occurring field water deficit stress. The genes identified in this study will provide candidate targets to manipulate the water use characteristics of cotton at the molecular level.

Increasing cotton genome coverage with polymorphic SSRs as revealed by SSCP

Simple sequence repeat (SSR) markers are widely used in plant genetics and breeding. However, there are many SSR markers that do not reveal polymorphism in cotton. Traditional SSR genotyping methods only provide information on product sizes. This leaves many marker polymorphism undetected, thus, lowering the utility of SSRs. In the present study, monomorphic SSRs between two mapping parents, 'Emian22' and 3-79, were subjected to single-strand conformation polymorphism (SSCP) analysis to reveal polymorphism. Of the 4194 monomorphic SSR primer pairs, 158 pairs (3.77%) showed polymorphism and revealed 174 polymorphic loci. Sequence analysis showed that the differences in PCR products between the mapping parents were solely due to base transition or transversion, which was in agreement with SSCP principles. SSCP also revealed SSRs with motifs of AT/TA and GAA/CTT were more polymorphic in dinucleotides and trinucleotides, respectively. Genetic mapping integrated 160 loci into our interspecific BC(1) linkage map, 5 of which associated with QTLs related to cotton fiber quality. The technique discussed in the present study enables us to detect polymorphism of monomorphic SSRs, and increase the utilization efficiency of the existing SSR primers.

Generation of ESTs for flowering gene discovery and SSR marker development in upland cotton

Background

Upland cotton, Gossypium hirsutum L., is one of the world's most important economic crops. In the absence of the entire genomic sequence, a large number of expressed sequence tag (EST) resources of upland cotton have been generated and used in several studies. However, information about the flower development of this species is rare.

Methodology/Principal Findings

To clarify the molecular mechanism of flower development in upland cotton, 22,915 high-quality ESTs were generated and assembled into 14,373 unique sequences consisting of 4,563 contigs and 9,810 singletons from a normalized and full-length cDNA library constructed from pooled RNA isolated from shoot apexes, squares, and flowers. Comparative analysis indicated that 5,352 unique sequences had no high-degree matches to the cotton public database. Functional annotation showed that several upland cotton homologs with flowering-related genes were identified in our library. The majority of these genes were specifically expressed in flowering-related tissues. Three GhSEP (G. hirsutum L. SEPALLATA) genes determining floral organ development were cloned, and quantitative real-time PCR (qRT-PCR) revealed that these genes were expressed preferentially in squares or flowers. Furthermore, 670 new putative microsatellites with flanking sequences sufficient for primer design were identified from the 645 unigenes. Twenty-five EST–simple sequence repeats were randomly selected for validation and transferability testing in 17 Gossypium species. Of these, 23 were identified as true-to-type simple sequence repeat loci and were highly transferable among Gossypium species.

Conclusions/Significance

A high-quality, normalized, full-length cDNA library with a total of 14,373 unique ESTs was generated to provide sequence information for gene discovery and marker development related to upland cotton flower development. These EST resources form a valuable foundation for gene expression profiling analysis, functional analysis of newly discovered genes, genetic linkage, and quantitative trait loci analysis.

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 an ADP-glucose pyrophosphorylase small subunit gene expressed in developing cotton (Gossypium hirsutum) fibers

ADP-glucose pyrophosphorylase (ADPGp, EC 2.7.7.27) is a tetrameric protein composed of two small and two large subunits that catalyzes the biosynthesis of ADP-glucose from glucose-phosphate which is used to provide the glucose subunits for starch biosynthesis. A second cotton gene encoding an ADPGp small subunit has been cloned and characterized. The gene contains eight introns similar to previously reported potato and cotton ADPGp small subunit genes. The deduced translation of the gene contained a poorly conserved transit peptide and well conserved catalytic and regulatory elements typical of other plant ADPGps. The 5' end of the mRNA was cloned and sequenced to identify the transcriptional start site (TSS). The promoter region upstream of the TSS did not contain the core promoter sequence in the typical positions indicating this gene may not use a standard core promoter. Other sequence motifs associated with tissue specific expression and phytohormone response were present. Reverse transcription (RT)-PCR with gene specific primers identified the sites of expression of this gene. Expression was most abundant in the meristem region, and immature stem and relatively lower in starch accumulating roots demonstrating that this gene has a different pattern of expression than the previously reported cotton ADPGp small subunit gene. Additionally this gene was differentially expressed in cotton fibers. The presence of starch was confirmed in developing cotton fibers suggesting that starch metabolism plays a role in cotton fiber development.

Characteristics of 11 polymorphic microsatellite markers in the red imported fire ant, Solenopsis invicta Buren

We have characterized 11 polymorphic microsatellite loci in the invasive ant Solenopsis invicta. Primer pairs were evaluated on fire ants collected from monogyne mounds in Lauderdale County, Mississippi. The observed and effective number of alleles ranged from two to six and from 1.31 to 2.64, respectively. The observed and expected heterozygosity values ranged from 0.1613 to 0.7826 and from 0.1491 to 0.6242, respectively. The polymorphism information content of the microsatellites ranged from 0.1482 to 0.6208. Probability tests indicated significant deviations from the Hardy-Weinberg equilibrium at three loci. Pairwise tests did not detect linkage disequilibrium between any pair of loci.

Isolation and characterization of a cotton cdh-like gene

Cotton fiber cells elongate without dividing to form economically valuable spinnable fiber. Reports of the ploidy level of fiber cells are variable. Early reports indicated an increase in nuclear DNA content in young fibers; however, subsequent reports failed to observe such a significant increase in ploidy level. Evaluation and analysis of genes involved in regulation of DNA synthesis and other aspects of cell cycle regulation identified relevant genes that were present in fiber cells though usually at low levels. We report the isolation and characterization of another gene likely to be involved in cell cycle/DNA synthesis control. This gene was similar to a gene from Medicago species that controls entry into anaphase by regulating the activity of the anaphase promoting complex ability to ubiquinate selected proteins. The cotton gene was composed of nine exons and the deduced translational sequences have motifs similar to a Medicago gene expressed in highly polyploid cells. Based on this similarity the cotton gene was designated Ghcdh. Ghcdh is highly expressed in meristems and leaves but is present at much lower levels in fiber cells. These data are consistent with the lower levels of polyploidy reported for cotton fiber. A simple sequence repeat was identified in the gene that may be exploited as a marker to map this gene and associate it with important traits in cotton.

Expression of genes associated with carbohydrate metabolism in cotton stems and roots

BACKGROUND

Cotton (Gossypium hirsutum L) is an important crop worldwide that provides fiber for the textile industry. Cotton is a perennial plant that stores starch in stems and roots to provide carbohydrates for growth in subsequent seasons. Domesticated cotton makes these reserves available to developing seeds which impacts seed yield. The goals of these analyses were to identify genes and physiological pathways that establish cotton stems and roots as physiological sinks and investigate the role these pathways play in cotton development during seed set.

RESULTS

Analysis of field-grown cotton plants indicated that starch levels peaked about the time of first anthesis and then declined similar to reports in greenhouse-grown cotton plants. Starch accumulated along the length of the stem and the shape and size of the starch grains from stems were easily distinguished from transient starch. Microarray analyses compared gene expression in tissues containing low levels of starch with tissues rapidly accumulating starch. Statistical analysis of differentially expressed genes indicated increased expression among genes associated with starch synthesis, starch degradation, hexose metabolism, raffinose synthesis and trehalose synthesis. The anticipated changes in these sugars were largely confirmed by measuring soluble sugars in selected tissues.

CONCLUSION

In domesticated cotton starch stored prior to flowering was available to support seed production. Starch accumulation observed in young field-grown plants was not observed in greenhouse grown plants. A suite of genes associated with starch biosynthesis was identified. The pathway for starch utilization after flowering was associated with an increase in expression of a glucan water dikinase gene as has been implicated in utilization of transient starch. Changes in raffinose levels and levels of expression of genes controlling trehalose and raffinose biosynthesis were also observed in vegetative cotton tissues as plants age.

Characteristics and analysis of simple sequence repeats in the cotton genome based on a linkage map constructed from a BC1 population between Gossypium hirsutum and G. barbadense

In the past decade, several molecular maps of cotton have been constructed using diverse DNA molecular markers and mapping populations. In this study, an interspecific linkage map of allotetraploid cotton was developed using a BC1 population ((Gossypium hirsutum x G. barbadense) x G. hirsutum). This map was genome-wide and was based entirely on simple sequence repeat (SSR) markers. Forty-four linkage groups were assigned to 26 chromosomes, with 917 loci spanning 5452.2 cM of the genome. The average distance between loci was 5.9 cM, providing uniform coverage of the A subgenome and D subgenome. Characteristics of this map were analyzed in detail, including the distributions of genomic SSRs, expressed sequence tag (EST)-SSRs, and distorted markers. Furthermore, the relationships between motif characteristics (size, type, length) and the level of polymorphism in EST-SSRs were also surveyed. The results showed that tetranucleotide and dinucleotide repeats had similar levels of polymorphism, and ACAT, AC, and ACT repeats had the highest polymorphism rates. Loci with lengths of 27 bp, 33 bp, and 24 bp were more likely to be polymorphic. This work will provide information to assist in designing future EST-SSRs.

Cotton benzoquinone reductase: up-regulation during early fiber development and heterologous expression and characterization in Pichia pastoris

Benzoquinone reductase (BR; EC 1.6.5.7) is an enzyme which catalyzes the bivalent redox reactions of quinones without the production of free radical intermediates. Using 2D-PAGE comparisons, two proteins were found to be up-regulated in wild-type cotton ovules during the fiber initiation stage but not in the fiberless line SL 1-7-1. These proteins were excised from the gel, partially sequenced and identified to be BR isoforms. PCR was used to amplify both full length coding regions of 609bp and once cloned, the restriction enzyme HindIII was used to distinguish the clones encoding the BR1 (one site) and BR2 (two sites) isoforms. Both deduced protein sequences had 203 residues which differed at 14 residues. The molecular mass and pIs were similar between the measured protein (2D-PAGE) and the theoretical protein (deduced). Heterologous proteins BR1 and BR2 were produced for further study by ligating the BR1 and BR2 clones in frame into the alpha-factor secretion sequence in pPICZalphaA vector and expressed with Pichia pastoris. Both BR1 and BR2 were approximately 26.5kDa and did enzymatically reduce 2,6-dimethoxybenzoquinone similar to the fungal BR.

Analysis of gene expression in cotton fiber initials

BACKGROUND

Cotton (Gossypium hirsutum L.) fibers are trichomes that initiate from the ovule epidermis. Little is known about the developmental pathway causing fiber to differentiate from ovular epidermal cells even though limits on the number of cells that differentiate into fiber will limit yield.

RESULTS

A method was developed to isolate RNA from fiber initials 1 day post anthesis (dpa). Complementary DNA libraries representing 1 dpa fibers and other cotton tissues were sequenced and analyzed. Assembly of G. hirsutum Expressed Sequenced Tags (ESTs) identified over 11,000 sequences not previously represented in GenBank. New genes identified among these ESTs were represented on microarrays. The microarrays were used to identify genes enriched in fiber initials (1 dpa fibers) and elongating fibers. Analyses of Gene Ontologies (GO) of differentially expressed genes determined that terms associated with the "membranes" were statistically over represented among genes increased in expression in fiber initials and 10 dpa fibers. Staining ovules with a fluorescent dye confirmed an increase in Endoplasmic Reticulum (ER) occurred in fiber initials on the day of anthesis, persisted through 3 dpa and was absent in a fiberless mutant. Two genes similar to the CAPRICE/TRIPTYCHON (CPC) gene that inhibits differentiation of leaf trichomes in Arabidopsis were also characterized. Genes associated with novel regulation of brassinosterols, GTP mediated signal transduction and cell cycle control and components of a Ca+2 mediated signaling pathway were identified. Staining of cellular Ca+2 indicated that fiber initials had more Ca+2 than other ovule cells supporting a role for Ca+2 in fiber development.

CONCLUSION

Analysis of genes expressed in fiber initials identified a unique stage in fiber development characterized by an increase in ER and Ca+2 levels that occurred between 0 and 1 dpa. The gene similar to CPC has a MYB domain but appears to lack a transcription activating domain similar to the Arabisopsis gene. The method used to stain the ER also can be used to count fiber initials and showed fiber cells develop from adjacent cells unlike leaf trichomes.

A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in Gossypium

The mapping of functional genes plays an important role in studies of genome structure, function, and evolution, as well as allowing gene cloning and marker-assisted selection to improve agriculturally important traits. Simple sequence repeats (SSRs) developed from expressed sequence tags (ESTs), EST-SSR (eSSR), can be employed as putative functional marker loci to easily tag corresponding functional genes. In this paper, 2218 eSSRs, 1554 from G. raimondii-derived and 754 from G. hirsutum-derived ESTs, were developed and used to screen polymorphisms to enhance our backbone genetic map in allotetraploid cotton. Of the 1554 G. raimondii-derived eSSRs, 744 eSSRs were able to successfully amplify polymorphisms between our two mapping parents, TM-1 and Hai7124, presenting a polymorphic rate of 47.9%. However, only a 23.9% (159/754) polymorphic rate was produced from G. hirsutum-derived eSSRs. No relationship was observed between the level of polymorphism, motif type, and tissue origin, but the polymorphism appeared to be correlated with repeat type. After integrating these new eSSRs, our enhanced genetic map consists of 1790 loci in 26 linkage groups and covers 3425.8 cM with an average intermarker distance of 1.91 cM. This microsatellite-based, gene-rich linkage map contains 71.96% functional marker loci, of which 87.11% are eSSR loci. There were 132 duplicated loci bridging 13 homeologous At/Dt chromosome pairs. Two reciprocal translocations after polyploidization between A2 and A3, and between A4 and A5, chromosomes were further confirmed. A functional analysis of 975 ESTs producing 1122 eSSR loci tagged in the map revealed that 60% had clear BLASTX hits (<1e(-10)) to the Uniprot database and that 475 were associated mainly with genes belonging to the three major gene ontology categories of biological process, cellular component, and molecular function; many of the ESTs were associated with two or more category functions. The results presented here will provide new insights for future investigations of functional and evolutionary genomics, especially those associated with cotton fiber improvement.