Genetic dissection of tetraploid cotton resistant to Verticillium wilt using interspecific chromosome segment introgression lines

Unraveling key genes and pathways involved in Verticillium wilt resistance by integrative GWAS and transcriptomic approaches in Upland cotton

Verticillium dahliae Kleb, the cause of Verticillium wilt, is a particularly destructive soil-borne vascular disease that affects cotton, resulting in serious decline in fiber quality and causing significant losses in cotton production worldwide. However, the progress in identification of wilt-resistance loci or genes in cotton has been limited, most probably due to the highly complex genetic nature of the trait. Nevertheless, the molecular mechanism behind the Verticillium wilt resistance remains poorly understood. In the present study, we investigated the phenotypic variations in Verticillium tolerance and conducted a genome wide association study (GWAS) among a natural population containing 383 accessions of upland cotton germplasm and performed transcriptomic analysis of cotton genotypes with differential responses to Verticillium wilt. GWAS detected 70 significant SNPs and 116 genes associated with resistance loci in two peak signals on D02 and D11 in E1. The transcriptome analysis identified a total of 2689 and 13289 differentially expressed genes (DEGs) among the Verticillium wilt-tolerant (J46) and wilt-susceptible (J11) genotypes, respectively. The DEGs were predominantly enriched in metabolism, plant hormone signal transduction, phenylpropanoid pathway, MAPK cascade pathway and plant-pathogen interaction pathway in GO and KEGG analyses. The identified DEGs were found to comprise several transcription factor (TF) gene families, primarily including AP2/ERF, ZF, WRKY, NAC and MYB, in addition to pentatricopeptide repeat (PPR) proteins and Resistance (R) genes. Finally, by integrating the two results, 34 candidate genes were found to overlap between GWAS and RNA-seq analyses, associated with Verticillium-wilt resistance, including WRKY, MYB, CYP and RGA. This work contributes to our knowledge of the molecular processes underlying cotton responses to Verticillium wilt, offering crucial insights for additional research into the genes and pathways implicated in these responses and paving the way for developing Verticillium wilt-resistant cotton varieties through accelerated breeding by providing a plethora of candidate genes.

Genomic confirmation of Gossypium barbadense introgression into G. hirsutum and a subsequent MAGIC population

Introgression of superior fiber traits from Pima cotton (Gossypium barbadense, GB) into high yield Upland cotton (G. hirsutum) has been a breeding objective for many years in a few breeding programs in the world. However, progress has been very slow due to introgression barriers resulting from whole genome hybridization between the two species. To minimize such barriers, chromosome substitution lines (CS-B) from Pima cotton 3-79 in an Upland cotton cultivar TM-1 were developed. A multiparent advanced generation inter-cross (MAGIC) population consisting of 180 recombinant inbred lines (RILs) was subsequently made using the 18 CS-B lines and three Upland cotton cultivars as parents. In this research, we sequenced the whole genomes of the 21 parents and 180 RILs to examine the G. barbadense introgression. Of the 18 CS-B lines, 11 contained the target GB chromosome or chromosome segment, two contained more than two GB chromosomes, and five did not have the expected introgression. Residual introgression in non-target chromosomes was prevalent in all CS-B lines. A clear structure existed in the MAGIC population and the 180 RILs were distributed into three groups, i.e., high, moderate, and low GB introgression. Large blocks of GB chromosome introgression were still present in some RILs after five cycles of random-mating, an indication of recombination suppression or other unknown reasons present in the population. Identity by descent analysis revealed that the MAGIC RILs contained less introgression than expected. This research presents an insight on understanding the complex problems of introgression between cotton species.

Genome-wide quantitative trait loci mapping on Verticillium wilt resistance in 300 chromosome segment substitution lines from Gossypium hirsutum × Gossypium barbadense

Cotton Verticillium wilt (VW) is a devastating disease seriously affecting fiber yield and quality, and the most effective and economical prevention measure at present is selection and extension of Gossypium varieties harboring high resistance to VW. However, multiple attempts to improve the VW resistance of the most widely cultivated upland cottons have made little significant progress. The introduction of chromosome segment substitution lines (CSSLs) provide the practical solutions for merging the superior genes related with high yield and wide adaptation from Gossypium hirsutum and VW resistance and the excellent fiber quality from Gossypium barbadense. In this study, 300 CSSLs were chosen from the developed BC5F3:5 CSSLs constructed from CCRI36 (G. hirsutum) and Hai1 (G. barbadense) to conduct quantitative trait locus (QTL) mapping of VW resistance, and a total of 40 QTL relevant to VW disease index (DI) were identified. Phenotypic data were obtained from a 2-year investigation in two fields with two replications per year. All the QTL were distributed on 21 chromosomes, with phenotypic variation of 1.05%-10.52%, and 21 stable QTL were consistent in at least two environments. Based on a meta-analysis, 34 novel QTL were identified, while 6 loci were consistent with previously identified QTL. Meanwhile, 70 QTL hotspot regions were detected, including 44 novel regions. This study concentrates on QTL identification and screening for hotspot regions related with VW in the 300 CSSLs, and the results lay a solid foundation not only for revealing the genetic and molecular mechanisms of VW resistance but also for further fine mapping, gene cloning and molecular designing in breeding programs for resistant cotton varieties.

Transcriptomic and biochemical analysis of upland cotton (Gossypium hirsutum) and a chromosome segment substitution line from G. hirsutum × G. barbadense in response to Verticillium dahliae infection

BACKGROUND

Verticillium wilt (VW), also known as "cotton cancer," is one of the most destructive diseases in global cotton production that seriously impacts fiber yield and quality. Despite numerous attempts, little significant progress has been made in improving the VW resistance of upland cotton. The development of chromosome segment substitution lines (CSSLs) from Gossypium hirsutum × G. barbadense has emerged as a means of simultaneously developing new cotton varieties with high-yield, superior fiber, and resistance to VW.

RESULTS

In this study, VW-resistant investigations were first conducted in an artificial greenhouse, a natural field, and diseased nursery conditions, resulting in the identification of one stably VW-resistant CSSL, MBI8255, and one VW-susceptible G. hirsutum, CCRI36, which were subsequently subjected to biochemical tests and transcriptome sequencing during V991 infection (0, 1, and 2 days after inoculation). Eighteen root samples with three replications were collected to perform multiple comparisons of enzyme activity and biochemical substance contents. The findings indicated that VW resistance was positively correlated with peroxidase and polyphenol oxidase activity, but negatively correlated with malondialdehyde content. Additionally, RNA sequencing was used for the same root samples, resulting in a total of 77,412 genes, of which 23,180 differentially expressed genes were identified from multiple comparisons between samples. After Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the expression profiles identified using Short Time-series Expression Miner, we found that the metabolic process in the biological process, as well as the pathways of phenylpropanoid biosynthesis and plant hormone signal transduction, participated significantly in the response to VW. Gene functional annotation and expression quantity analysis indicated the important roles of the phenylpropanoid metabolic pathway and oxidation-reduction process in response to VW, which also provided plenty of candidate genes related to plant resistance.

CONCLUSIONS

This study concentrates on the preliminary response to V991 infection by comparing the VW-resistant CSSL and its VW-susceptible recurrent parent. Not only do our findings facilitate the culturing of new resistant varieties with high yield and superior performance, but they also broaden our understanding of the mechanisms of cotton resistance to VW.

Transcriptomic and biochemical analysis of upland cotton (Gossypium hirsutum) and a chromosome segment substitution line from G. hirsutum × G. barbadense in response to Verticillium dahliae infection

BACKGROUND

Verticillium wilt (VW), also known as "cotton cancer," is one of the most destructive diseases in global cotton production that seriously impacts fiber yield and quality. Despite numerous attempts, little significant progress has been made in improving the VW resistance of upland cotton. The development of chromosome segment substitution lines (CSSLs) from Gossypium hirsutum × G. barbadense has emerged as a means of simultaneously developing new cotton varieties with high-yield, superior fiber, and resistance to VW.

RESULTS

In this study, VW-resistant investigations were first conducted in an artificial greenhouse, a natural field, and diseased nursery conditions, resulting in the identification of one stably VW-resistant CSSL, MBI8255, and one VW-susceptible G. hirsutum, CCRI36, which were subsequently subjected to biochemical tests and transcriptome sequencing during V991 infection (0, 1, and 2 days after inoculation). Eighteen root samples with three replications were collected to perform multiple comparisons of enzyme activity and biochemical substance contents. The findings indicated that VW resistance was positively correlated with peroxidase and polyphenol oxidase activity, but negatively correlated with malondialdehyde content. Additionally, RNA sequencing was used for the same root samples, resulting in a total of 77,412 genes, of which 23,180 differentially expressed genes were identified from multiple comparisons between samples. After Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the expression profiles identified using Short Time-series Expression Miner, we found that the metabolic process in the biological process, as well as the pathways of phenylpropanoid biosynthesis and plant hormone signal transduction, participated significantly in the response to VW. Gene functional annotation and expression quantity analysis indicated the important roles of the phenylpropanoid metabolic pathway and oxidation-reduction process in response to VW, which also provided plenty of candidate genes related to plant resistance.

CONCLUSIONS

This study concentrates on the preliminary response to V991 infection by comparing the VW-resistant CSSL and its VW-susceptible recurrent parent. Not only do our findings facilitate the culturing of new resistant varieties with high yield and superior performance, but they also broaden our understanding of the mechanisms of cotton resistance to VW.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Gbvdr6, a Gene Encoding a Receptor-Like Protein of Cotton (Gossypium barbadense), Confers Resistance to Verticillium Wilt in Arabidopsis and Upland Cotton

Verticillium wilt is a soil-borne disease that can cause devastating losses in cotton production. Because there is no effective chemical means to combat the disease, the only effective way to control Verticillium wilt is through genetic improvement. Therefore, the identification of additional disease-resistance genes will benefit efforts toward the genetic improvement of cotton resistance to Verticillium wilt. Based on screening of a BAC library with a partial Ve homologous fragment and expression analysis, a V. dahliae-induced gene, Gbvdr6, was isolated and cloned from the Verticillium wilt-resistant cotton G. barbadense cultivar Hai7124. The gene was located in the gene cluster containing Gbve1 and Gbvdr5 and adjacent to the Verticillium wilt-resistance QTL hotspot. Gbvdr6 was induced by Verticillium dahliae Kleb and by the plant hormones salicylic acid (SA), methyl jasmonate (MeJA) and ethephon (ETH) but not by abscisic acid (ABA). Gbvdr6 was localized to the plasma membrane. Overexpression of Gbvdr6 in Arabidopsis and cotton enhanced resistance to V. dahliae. Moreover, the JA/ET signaling pathway-related genes PR3, PDF 1.2, ERF1 and the SA-related genes PR1 and PR2 were constitutively expressed in transgenic plants. Gbvdr6-overexpressing Arabidopsis was less sensitive than the wild-type plant to MeJA. Furthermore, the accumulation of reactive oxygen species and callose was triggered at early time points after V. dahliae infection. These results suggest that Gbvdr6 confers resistance to V. dahliae through regulation of the JA/ET and SA signaling pathways.

Quantitative Trait Locus Mapping and Candidate Gene Analysis for Verticillium Wilt Resistance Using Gossypium barbadense Chromosomal Segment Introgressed Line

Verticillium wilt (VW) is a soil-borne fungal disease that is caused by Verticillium dahliae Kleb and seriously damages cotton production annually in China. To date, many efforts have been made to improve the resistance of upland cotton against VW, but little progress has been achieved because of a lack of resistant upland cotton to VW. G. barbadense is known to carry high resistance to VW; however, it is difficult to transfer the resistance trait from G. barbadense to upland cotton because of linkage drag and distortion in the interspecific hybrid. In this study, a chromosomal segment introgression line (CSIL), SuVR043, containing a single and homozygous chromosome segment of G. barbadense cv. H7124 D04 (Chr 22), was created and used to construct an F₂ population for mapping of VW resistance quantitative trait loci (QTLs) in the greenhouse. Two major resistance QTLs against nondefoliating V. dahliae isolate Bp2, called qVW-Bp2-1 and qVW-Bp2-2, which were flanked by the markers cgr6409-ZHX37 and ZHX57-ZHX70 and explained an average of 16.38 and 22.36% of the observed phenotypic variation, respectively, were detected in three independent replicate experiments. The genetic distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 were 2.4 and 0.8 cM, respectively. By analyzing the genome sequence of the qVW-Bp2-1 and qVW-Bp2-2 regions, we determined that the accurate physical distances from cgr6409 to ZHX37 and from ZHX57 to ZHX70 in the G. barbadense genome are 254 and 140 kb, and that those spans 36 and 20 putative genes, respectively. The results of the expression analysis showed significant differences in the expression profiles of GbCYP450, GbTMEM214, and GbRLK among G. barbadense cv. H7124, CSIL SuVR043 and G. hirsutum acc. Sumian 8 at different times after inoculation with V. dahliae isolate Bp2. Virus-induced gene silencing (VIGS) analysis showed that silencing of GbCYP450 and GbTMEM214 decreased H7124 and CSIL SuVR043 resistance to VW. These results form a solid foundation for fine mapping and cloning of resistance genes in the substituted segment and will provide valuable assistance in future efforts to breed for VW resistance.

Regional association analysis-based fine mapping of three clustered QTL for verticillium wilt resistance in cotton (G. hirsutum. L)

Background

Verticillium wilt is one of the most destructive diseases affecting global cotton production. The most effective way to control wilt disease has been the development of new cotton varieties that are resistant to VW. VW-resistant Upland cotton cultivars have been created in both the USA and China by Gossypium barbadense introgression. More than 100 VW resistance quantitative trait loci have been detected.

Results

Three clustered VW resistance-related QTL were detected in a 120-line association population and assigned to a genome region of 14,653,469–55,190,112 bp in Dt_chr9. A regional association analysis-based fine-mapping strategy was developed to narrow down the confidence intervals of the above QTL. The estimated LD decay of the genome region of interest was much faster than those of the Dt_chr9 chromosome and the whole genome, suggesting the existence of a recombination hotspot. Thirty-seven haplotype blocks were detected. The confidence intervals of the three clustered QTL were narrowed down to a region of 937,906 bp involving QTL-i23734Gh and a region of 1,389,417 bp involving QTL- i10740Gh, respectively. Each region contained the strongest association signal. Comparative analysis redefined the confidence intervals of the other three QTLs, qDL52T2-c19, QTL-BNL4069, and QTL-JESPR0001. The broad-spectrum VW resistance QTL qVW-D9–1 was demonstrated to be closely linked with the three redefined QTL, QTL-i23734Gh, QTL- i10740Gh and QTL-JESPR0001. Twelve functional genes were detected to be located within the redefined confidence intervals of VW resistance QTL. The mRNA CotAD_60243, encoding E3 ubiquitin-protein ligase UPL2-like, responsible for plant innate immunity and broad-spectrum disease resistance, was found to be overlapped with the strongest association signal i10740Gh. Six mRNAs encoding putative disease-resistance proteins were within the redefined confidence interval of QTL-JESPR0001, suggesting a tandem arrangement of R genes.

Conclusions

Our results proved that the VW resistance effect related to three clustered VW resistance-related QTL was actually controled by two redefined major QTL and severlal minor loci. The broad-spectrum VW resistance QTL qVW-D9–1 may be closely linked with the two redefined major QTLs. The tandem arrangement of R genes were detected in the redefined confidence interval of QTL-JESPR0001. The candidate genes obtained should be helpful in identifying and characterizing defense genes related to VW resistance QTL.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4074-y) contains supplementary material, which is available to authorized users.

Quantitative Trait Locus Mapping for Verticillium wilt Resistance in an Upland Cotton Recombinant Inbred Line Using SNP-Based High Density Genetic Map

Verticillium wilt (VW) caused by Verticillium dahlia Kleb is one of the most destructive diseases of cotton. Numerous efforts have been made to improve the resistance of upland cotton against VW, with little progress achieved due to the paucity of upland cotton breeding germplasms with high level of resistance to VW. Gossypium barbadense was regarded as more resistant compared to upland cotton; however, it is difficult to apply the resistance from G. barbadense to upland cotton improvement because of the hybrid breakdown and the difficulty to fix resistant phenotype in their interspecific filial. Here we reported QTLs related to VW resistance identified in upland cotton based on 1 year experiment in greenhouse with six replications and 4 years investigations in field with two replications each year. In total, 119 QTLs of disease index (DI) and of disease incidence (DInc) were identified on 25 chromosome of cotton genome except chromosome 13 (c13). For DI, 62 QTLs explaining 3.7-12.2% of the observed phenotypic variations were detected on 24 chromosomes except c11 and c13. For DInc, 59 QTLs explaining 2.3-21.30% of the observed PV were identified on 19 chromosomes except c5, c8, c12-c13, c18-c19, and c26. Seven DI QTLs were detected to be stable in at least environments, among which six have sGK9708 alleles, while 28 DInc QTLs were detected to be stable in at least environments. Eighteen QTL clusters containing 40 QTLs were identified on 13 chromosomes (c1-c4, c6-c7, c10, c14, c17 c20-c22, and c24-c25). Most of the stable QTLs aggregated into these clusters. These QTLs and clusters identification can be an important step toward Verticillium wilt resistant gene cloning in upland cotton and provide useful information to understand the complex genetic bases of Verticillium wilt resistance.

Quantitative trait loci analysis of Verticillium wilt resistance in interspecific backcross populations of Gossypium hirsutum × Gossypium barbadense

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, BC₁F₁, BC₁S₁, and BC₂F₁ 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 BC₁F₁ 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 BC₂F₁ in the artificial disease nursery and one date of BC₁S₁ and BC₂F₁ 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.

A Ve homologous gene from Gossypium barbadense, Gbvdr3, enhances the defense response against Verticillium dahliae

The tomato Ve1 gene and several Ve1 homologues are involved in the resistance to Verticillium dahliae. Here, we report on another Ve homologous gene, Gbvdr3, from a Verticillium wilt-resistant cotton cultivar, Gossypium barbadense Hai7124, which has a 3207-bp region that encodes a predicted receptor-like protein. Transient expression analyses indicated that Gbvdr3 is localized in the plasma membrane, and virus-induced gene silencing of Gbvdr3 compromised the resistance of Hai7124 cotton to a defoliating strain of V. dahliae, V991, but not to a non-defoliating strain, BP2. This resistance pattern was further confirmed by over-expression of Gbvdr3 in transgenic Arabidopsis, which significantly elevated the expression of the ethylene-regulated gene GST2, the ethylene- and jasmonic acid-regulated defense-related genes PR3 and PDF1.2, and the salicylic acid-regulated genes PR1 and PR5, but not the PR2 gene. It also triggered the accumulation of hydrogen peroxide and callose at early time points during infection by the V991 defoliating strain. In contrast, elevated accumulation of hydrogen peroxide or callose in Gbvdr3-expressed Arabidopsis leaves was not apparent under infection by the non-defoliating strain, BP2. These results suggested that Gbvdr3 is involved in the resistance to a unique spectrum of defoliating V. dahliae strains.

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

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1682-2) contains supplementary material, which is available to authorized users.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.