<title>Background</title> <p>Cotton is the world's most important natural textile fiber and a significant oilseed crop. The cotton genus (<italic>Gossypium </italic>L.) includes approximately 45 diploid species (2<italic>n </italic>= 2<italic>x </italic>= 26) differentiated cytogenetically into eight genome groups (A-G &amp; K), and five allotetraploid species (2<italic>n </italic>= 4<italic>x </italic>= 52) [##UREF##0##1##]. Diploid <italic>Gossypium </italic>species differentiated approximately 5–10 million years ago (Mya), however, polyploidization is estimated to have occurred more recently 1–2 Mya [##UREF##1##2##]. All allotetraploids were formed from interspecific hybridization events between an A-genome-like ancestral African species and a D-genome-like North American species. The closest extant relative of the original tetraploid progenitors is the A-genome species <italic>G. herbaceum </italic>L. (A1) and the D-genome species <italic>G. raimondii </italic>(D5) Ulbrich. Of these, four cotton species, including two tetraploids <italic>G. hirsutum </italic>L. (AD)1 and <italic>G. barbadense </italic>L (AD)2, and two diploids <italic>G. herbaceum </italic>L. (A1) and <italic>G. arboreum </italic>L. (A2) were independently domesticated for fiber.</p> <p>Upland cotton has the highest yield, and based on the importance of fiber, over 95% of the annual worldwide cotton crop is derived from <italic>G. hirsutum </italic>L., upland cotton, and the extra-long staple (ELS) or Pima cotton (<italic>G. barbadense </italic>L.) accounts for less than 2% [##UREF##2##3##]. Two diploid species <italic>G. herbaceum </italic>L. (A1) and <italic>G. arboreum </italic>L. (A2) are planted less often. In cultivated tetraploid cotton species, the D-subgenome plays an important role in genome structure, function and evolution. For example, many quantitative trait loci (QTL) for fiber-related traits have been detected in the D-subgenome of tetraploid cotton [##REF##9539752##4##, ####UREF##3##5##, ##UREF##4##6##, ##UREF##5##7##, ##UREF##6##8##, ##REF##15618303##9####15618303##9##]. D-genome species do not produce spinnable fiber [##REF##11256432##10##]; however important genes or regulators for fiber morphogenesis and fiber properties have been detected in this genome. Based on the above analyses, understanding the contribution of the A- and D-subgenomes to gene expression in the allotetraploids may greatly facilitate fiber trait improvement [##UREF##7##11##,##REF##16889650##12##]. To attain this goal, decoding cotton genomes will be a foundation to enhance our understanding of the functional and agronomic significance of polyploidy and genome size variation within <italic>Gossypium </italic>[##REF##18056866##13##].</p> <p>Genome size differences are evident in the tetraploids and their diploid progenitors. The haploid genome size is estimated to be ~980-Mb for <italic>G. raimondii </italic>Ulbrich, ~1.86-Gb for <italic>G. arboreum </italic>L., and ~2.83 Gb for <italic>G. hirsutum </italic>L. [##REF##15256507##14##]. Diploid species variation in DNA content reflects increases and decreases in copy numbers of various repeat families [##REF##9582192##15##], especially retrotransposon-like elements [##REF##16954538##16##]. The method most appropriate for elucidating whole-genome sequence information in cotton is either BAC-by-BAC sequencing or gene-enrichment approaches. A pilot study by the U.S. Department of Energy Joint Genome Institutes [##UREF##8##17##] has been initiated to generate the whole-genome shotgun sequence of <italic>G. raimondii</italic>. Meanwhile, gene-enrichment techniques such as methylation filtration and C<sub>o</sub>t-based cloning have also been used to compare <italic>G. raimondii, G. arboreum, G. hirsutum</italic>, and <italic>G. barbadense </italic>(B. Scheffler, Workshop communication).</p> <p>The whole-genome sequence analysis of <italic>G. hirsutum </italic>will undoubtedly provide the ultimate reference and resource for structural, functional, and evolutionary studies of the species that accounts for &gt; 95% of world cotton production. Prior to large-scale sequencing of tetraploid <italic>G. hirsutum </italic>genomes, a microcolinearity analysis of a few pairs of homoeologous BACs was completed, and indicated that sequence conservation of homoeologous BACs was high in both intergenic and genic regions [##REF##15256507##14##]. In addition, Grover et al. (2007)[##REF##17461788##18##] suggested size differences between homoeologous BACs was attributed to differential accumulation of retroelements.</p> <p>The GeneTrek approach has been proposed as an efficient way to evaluate the general properties of any genome [##REF##15491923##19##,##REF##16357197##20##] and has been successfully applied to predictions regarding components of the maize genome [##REF##17615239##21##]. To better understand the general composition and structure of the tetraploid cotton genome, in the present paper, we also employed GeneTrek and BAC tagging information approaches to analyze. This methodology facilitated our evaluation of the structure and composition of the allotetraploid genome based on 142 <italic>G. hirsutum </italic>cv. Maxxa BAC clones downloaded from the National Center for Biotechnology Information (NCBI) [##UREF##9##22##]. The study provided us the first glimpse at cotton genome complexity, and the results indicated that the gene distribution in cotton genome is uneven with gene-rich and gene-free regions, and rich in repetitive elements. Introns might have no contribution to different subgenome size in <italic>Gossypium</italic>, and a two-fold genome difference between A- and D-subgenomes, which might largely be attributed to large amplifications of transposable elements in low-density gene or gene-free regions.</p>
<title>Methods</title> <title>Cotton BACs</title> <p>One hundred forty-five cotton BAC sequences were downloaded from the National Center for Biotechnology Information (NCBI) [##UREF##9##22##] on June 2, 2007. As part of the maize sequencing project by the Genome Sequencing Center, Washington University School of Medicine, the BACs were initially submitted as <italic>Zea mays</italic>. However, further analysis, determined the clones were from <italic>Gossypium hirsutum </italic>cv. Maxxa. The sequence data used in this paper were the product of collaborative efforts by The Maize Sequencing Consortium, including the University of Arizona, Cold Spring Harbor Laboratory, Iowa State University, and the Genome Sequencing Center at Washington University School of Medicine in St. Louis. We selected 142 from 145 BACs with sizes &gt; 20 kb for gene annotation. A 32,101 bp length gap region from AC189045 was later excluded because the predicted genes were phage related and it was decided the sequence data were contaminated. Finally, 142 BACs spanning nearly 14.2 Mb (0.5%) of the cotton genome were used for the analysis.</p> <title>Genetic mapping of BAC clones based on the simple sequence repeats (SSRs)</title> <p>Each BAC was searched for SSRs with the online software SSRIT [##REF##11483586##36##]. SSRIT, written in Perl script, is a microsatellite search tool available at the USDA-ARS Center for Bioinformatics and at Comparative Genomics at Cornell University. Dinucleotide, trinucleotide, tetranucleotide, pentanucleotide and hexanucleotide SSRs were detected with SSRIT. The search standards for different repeat motifs were as described in Wang et al. (2006) [##UREF##13##37##]. Primer pairs flanking the SSRs were designed using the program Primer3.0 [##UREF##14##38##] and tested against our mapping parents, <italic>G. barbadense </italic>cv. Hai7124 and <italic>G. hirsutum </italic>acc. TM-1, standard lines for genetic and genomic research. Furthermore, the polymorphic SSRs were integrated into our backbone genetic map of allotetraploid cultivated cotton [##REF##17409069##25##] using Joinmap 3.0 software with a minimum log-of-odds (LOD) score of 6.0. The structure of known BACs was further identified using mapping results and molecular size comparisons among <italic>G. hirsutum </italic>cv. Maxxa, <italic>G. hirsutum </italic>acc.TM-1, and <italic>G. barbadense </italic>cv. Hai7124, with diploid <italic>G. herbaceum </italic>var. <italic>africanum </italic>and <italic>G. raimondii </italic>as controls.</p> <title>Annotation of LTR retrotranspons and other mobile elements</title> <p>Repetitive element prediction was accomplished through Repeatmasker [##UREF##15##39##], CENSOR [##REF##8867843##40##], and BLAST identity to characterize elements in REPBASE (version 8.5) [##REF##10973072##41##]. Compared with the results of repetitive element prediction, LTR retrotranspons were further identified by LTR_finder software [##UREF##16##42##], and manually verified by structural features such as LTR and TSD pairs, a primer binding site and a polypurine tract.</p> <title>Sequence analysis and gene annotation</title> <p>BAC sequences were subject to three <italic>ab initio </italic>gene prediction programs, FGENESH (Softberry) [##UREF##17##43##], GENSCAN+[##REF##9149143##44##] and GENEMARK.HMM [##REF##9461475##45##]. Gene models provided query sequences to search the National Center for Biotechnology Information (NCBI) non-redundant protein database and the <italic>Arabidopsis thaliana </italic>protein database [##UREF##18##46##]. All BLASTP hits were manually evaluated to determine if a gene model was likely to be a real gene or not based on <italic>e</italic>-value, query alignment and hit annotation. The integral parts of known repetitive elements were removed from the above gene models for further analysis. In addition, the sequences with gene models but no annotations were subjected to BLASTN queries against the cotton EST database released in the NCBI website [##REF##16478941##47##]. Gene Ontology (GO) of tandem duplication genes was obtained from UniProt Gene Ontology [##UREF##19##48##]. The GO values for the best homologous hits were used to determine the ontology of molecular function, cellular components and biological processes for these sequences.</p>
<title>Results</title> <title>Confirmation of 142 BACs origin</title> <p>Due to the fact that 142 BACs were result from a mistake first submitted as part of the maize sequencing project by the Genome Sequencing Center, Washington University School of Medicine and further corrected as <italic>G. hirsutum </italic>cv. Maxxa BAC clones, we downloaded these BACs from the National Center for Biotechnology Information (NCBI) [##UREF##9##22##] and confirmed their origin by developing BAC-SSR markers from 142 BAC sequences.</p> <p>Each BAC was scanned for dinucleiotide to hexanucleiotide repeats of at least 18 bp in length. A total of 694 microsatellite sequences were detected. Among them, 208 SSRs were dinucleotides, 118 trinucleotides, 69 tetranucleotides, 80 pentanucleotides and 219 hexanucleotides. In addition, 578 SSR primer pairs were developed and used to detect the amplification ability in <italic>G. hirsutum </italic>cv. Maxxa, and our two mapping parents, <italic>G. hirsutum </italic>acc. TM-1 and <italic>G. barbadense </italic>cv. Hai7124. Among them, all 578 primer pairs amplified expected fragment sizes in <italic>G. hirsutum </italic>cv. Maxxa, and 161 primer pairs from 79 BACs amplified polymorphisms between TM-1 and Hai7124, yielding a 27.85% polymorphic rate. Both the high-level transferability among <italic>G. hirsutum </italic>cv. Maxxa, acc. TM-1, and <italic>G. barbadense </italic>cv. Hai7124 and the high-level polymorphism between TM-1 and Hai7124 indicated that these 142 BAC sequences must be from Maxxa genome. Further, these genomic SSR markers also have potential for use in future cotton genomics and molecular breeding. The newly developed SSR primer sequences, Genbank accession numbers, repeat motifs and numbers, expected product size, and polymorphic data between TM-1 and Hai7124 are presented in additional file ##SUPPL##0##1##.</p> <title>Global analysis of genome structure and composition of tetraploid cotton</title> <title>Gene annotations</title> <p>Using the sequence information of the 142 BACs spanning 14.2 Mb of the cotton genome, genome structure and composition of tetraploid cotton were analyzed. Comprehensively analyzing the gene prediction results from three <italic>ab initio </italic>gene prediction programs FGENESH, GENEMARK.HMM and GENSCAN+, 3,440 gene models were predicted. Of them, 1,329 (38.6%) were identified repeat components (mostly LTR retrotransposons), which were further analyzed with mobile elements; 1,653 (47.9%) lacked homology to other NCBI protein database species. Furthermore, the putative protein-encoding gene was subjected to BLASTN queries against the cotton EST database released in NCBI and two hundred eight were partially confirmed by EST evidence; and 458 showed homology to other species in the NCBI protein database. Based on significant homology to other species (<italic>e </italic>≤ 10<sup>-10</sup>), 412 gene models were classified as verified gene candidates, with an average gene density of one gene per 34.5 kb. Forty-six gene models were classified as hypothetical proteins. If these gene numbers are extrapolated to the entire tetraploid cotton genome, with an estimated size of 2,500 Mb, tetraploid cotton contains more than 70,000 (verified) genes (Table ##TAB##0##1##). Details on the annotation of each predicted gene can be found in additional file ##SUPPL##1##2##.</p> <title>Local gene density and distribution</title> <p>Among the 142 analyzed BACs, 30 (21%) did not contain either a verified or hypothetical gene. Furthermore, gene density was estimated as the number of genes on a BAC divided by BAC length. The results showed that on different BACs, gene density varied from 0 to 33.2 per 100 kb (Figure ##FIG##0##1##, Additional file ##SUPPL##2##3##). This indicated uneven cotton gene distribution, and a higher gene density in some regions than others. In AC188398 and AC189045 BACs, 30 and 19 gene models were predicted, respectively, lacking repetitive elements.</p> <title>Gene islands</title> <p>The number of genes per one gene island can be determined by gene distribution in gene-rich regions. The number of genes on one gene island is counted according to the following criterion: the identifiable repetitive sequences in the intergenic region between two neighboring genes must be less than 5 kb [##REF##17615239##21##]. Furthermore, genes at either end of a BAC or gap within one BAC are discarded from the analysis because one boundary of the gene island is not defined. Based on this criterion, 309 gene islands from one to ten genes (both verified and hypothetical genes) were resolved (Figure ##FIG##1##2##). Two hundred twenty-four of 309 (72.5%) gene islands contained only one gene with an average of 1.5 genes per gene island.</p> <title>Tandem duplication of genes</title> <p>Thirty gene islands contained more than two genes, and in those islands, several types of tandem duplication genes encoding the same function were identified (see Additional file ##SUPPL##3##4##). According to the molecular function classification of these duplication genes, most were related to binding, such as sar1 GTP-binding secretory factor, ire kinase, RNA-binding protein 10, swi2 snf2-like protein, succinate dehydrogenase flavoprotein alpha subunit, adenylate kinase, and sll2 protein. Other genes functioned in catalytic activities, including genes coding ornithine carbamoyltransferase, glucose-methanol-cholineoxidoreductase family protein, adenylosuccinate lyase, protein phosphatase-5, protein kinase family protein, methylmalonate-semialdehyde dehydrogenase, calcineurin-like phosphoesterase family proteins and serine carboxypeptidase ii. Additional genes were determined to serve in transporter activities such as plasma membrane intrinsic proteins, structural molecule activity such as 50s ribosomal protein l15, and unknown molecular function, such as growth-regulating factor 1, among others. Several disease-resistant gene clusters resided in AC187066, AC190836 and AC202830 BACs. These specific gene clusters presumably accumulated more mutations in both coding and upstream promoter regions to favor a broader response to pathogen attack [##REF##12524363##23##]. Several QTLs related to <italic>Verticillium</italic>-resistance [##UREF##10##24##] were also found in these regions, but warrants further investigation.</p> <title>Mobile elements analysis</title> <p>RepeatMasker and CENSOR program was first applied to search for repetitive elements. In all, 1,951 mobile elements with a total length of 1,468,873 bp were predicted (Table ##TAB##1##2##). Retro-elements were a major component and accounted for 93.9% of the predicted elements with LTR/gypsy comprising 61.1%, LTR/copia 31.2% and LINE elements 1.5%. Four types of DNA transposons were identified, including hobo-Activator (0.3%), En-Spm (3.6%), MUDR-IS905 (0.6%) and Tourist/Harbinger (0.2%). DNA transposons accounted for 4.7% of all predicted elements. A scan of the 142 cotton BACs predicted 0.2% RC/Helitron and 1.1% unclassified mobile elements. Three hundred forty-four intact LTR retrotransposons were predicted and identified by LTR_finder software, leading to an additional 656,779 bp of LTR repetitive sequences. Including the 3,526,152 bp sequence length repeats identified via gene models, mobile DNA accounted for at least 5.7 Mb or 40% of the BAC sequences (Table ##TAB##0##1##).</p> <title>Sequence analysis of gene-free BACs</title> <p>Among the 142 analyzed BACs, 30 showed the complete absence of genes. To further investigate the content of such genomic regions, seven of the 30 BACs with complete sequence assembly were selected. These regions were largely comprised of LTR retrotransposons and were all the primary components of all seven BACs. In retrotransposon types, a number of fragmented <italic>gypsy</italic>-like elements were found in a large "<italic>gypsy</italic>-landing pad", indicating <italic>gypsy</italic>-like retroelements were substantial components of gene-free BACs. All retroelements identified in AC187849, AC194319, AC189931 and AC190814 were <italic>gypsy </italic>elements. The second component of retroelements was a <italic>copia</italic>-like element. The seven gene-free BACs showed the absence of any LINE element (Table ##TAB##2##3##). Most LTR retrotransposons in the seven BACs were truncated, and only 18 had two intact LTRs and target site duplications (TSD). Most of these LTRs were organized in nested structures.</p> <title>Comparative analysis of genome structure and composition between A- and D- subgenome chromosomes</title> <title>Temporal mapping of 70 BACs based on SSRs</title> <p>To compare the A- and D-subgenome chromosome structures and compositions, the BACs must be anchored into their corresponding subgenome or chromosomes. We firstly based the present research on developed SSR markers and our mapping population. By polymorphism analysis, one hundred sixty-one primer pairs developed from 79 BACs could produce 183 polymorphic loci in the two mapping parents, TM-1 and Hai7124. Because TM-1 was used as the recurrent parent in the backcrossed population, 17 dominant TM-1 loci, (of the 183 polymorphic loci) could not be used to anchor the related-BACs to their corresponding chromosomes. The remaining 166 polymorphic loci, amplified from 144 polymorphic SSR primer pairs were integrated into our previously published map containing 1,790 loci and spanning 3425.8 cM [##REF##17409069##25##]. Subsequently, a new updated genetic map composed of 2,247 loci in 26 linkage groups covering 3540.4 cM with an average inter-marker distance of 1.58 cM was produced (Figure ##FIG##2##3##, ##FIG##3##4##, ##FIG##4##5##, ##FIG##5##6##, ##FIG##6##7##, ##FIG##7##8##, ##FIG##8##9##). Based on the new integrated genetic map, 166 polymorphic loci developed from 70 BAC clones were anchored to their corresponding chromosomes. Further analysis of 70 BACs with tagged results revealed 18 BACs, one possessing more than two polymorphic loci, which were then mapped within 0.5 cM in one linkage group. For example, our results found six polymorphic loci produced by SSR primer pairs NAU6520, NAU6530, NAU6593, NAU6658, NAU6675, and NAU6697 from BAC AC202830 all tagged in chromosome D11 (Figure ##FIG##7##8##). We also found 13 BACs in which more than two polymorphic loci from one BAC were mapped in their homoeologous chromosome pairs. For example, among six polymorphic loci amplified by SSR primers developed from the same BAC (AC188035), two loci amplified by NAU6615 and NAU6667 were tagged in chromosome A10, however, four loci by NAU6215, NAU6476, NAU6515 and NAU6667 tagged in chromosome D10 (Figure ##FIG##6##7##). Additionally, 11 BACs with more than two polymorphic loci produced by SSR primer pairs from the same BAC were tagged in non-homoeologous chromosomes (for example, two polymorphic loci amplified by NAU6389 and NAU6626 developed from AC190805 were anchored in A4 and D8, respectively (Figure ##FIG##4##5##, ##FIG##5##6##)) and 28 BACs in which only one polymorphic locus was tagged (see Additional file ##SUPPL##0##1##). Based on these SSR mapping results, however, we could not definitively establish if the mapped polymorphic locus is at exactly the BAC clone position, since these polymorphic loci can be amplified either from the A- or D-subgenome and even other chromosomes in allotetraploid cotton.</p> <title>Identified tagging of 37 BACs based on amplified product analysis</title> <p>We are left to question if we can exactly anchor these BACs into their corresponding chromosomes by their-derived SSR markers. It is well known that allotetraploid cotton contains two distinct genomes, which resemble the extant A-genome of <italic>G. herbaceum </italic>(n = 13) and D-genome of G. <italic>raimondii </italic>Ulbrich (n = 13). The A- and D-genome species diverged from a common ancestor approximately 6–11 million years ago. Therefore, most SSR primers should easily amplify two loci, one from At- and the second from Dt-subgenome chromosomes. So, if one polymorphic locus was detected and mapped between Hai7124 and TM-1, and the amplified product was the same as that in Maxxa (from which BAC clones were isolated and SSRs developed), we concluded that this SSR derived-BAC should be anchored at its SSR tagging position. For example, SSR primer NAU6627, derived from AC187848, generated two loci from <italic>G. hirsutum </italic>acc. TM-1, and <italic>G. barbadense </italic>cv. Hai7124 (Figure ##FIG##9##10A##). One locus produced two polymorphic alleles, NAU6627<sub>-250 </sub>in Hai7124 and NAU6627<sub>-247 </sub>in TM-1, and mapped on chromosome D11 (Figure ##FIG##7##8##, Additional file ##SUPPL##0##1##). NAU6627 was designed based on sequence information from AC187848 in Maxxa and its expected product size was 247 bp in Maxxa, therefore we concluded that the BAC was anchored into chromosome D11 in the D-subgenome. Following further analysis of amplified products from primer NAU6627 in diploid species <italic>G. herbaceum </italic>var. <italic>africanum </italic>and <italic>G. raimondii</italic>, (the two closest extant relatives of the original tetraploid progenitors), we could still anchor this polymorphic locus into the D-subgenome, given <italic>G. raimondii </italic>produced almost the same fragment at the expected 247 bp product size. Furthermore, two loci were amplified from each TM-1 and Hai7124, with one monomorphic locus producing the expected product size from the BAC clone in Maxxa and the other a polymorphic locus tagged in the genetic map. Therefore, we associated this BAC clone with the subgenome chromosome by tagging comparisons of the polymorphic locus and the amplified products from <italic>G. herbaceum </italic>var. <italic>africanum </italic>and <italic>G. raimondii</italic>. The SSR primer NAU6202 derived from AC190263 generated two loci in both TM-1 and Hai7124 (Figure ##FIG##9##10B##), validating the former results. One SSR locus was monomorphic, the alleles approximately 350 bp in size, and the other locus produced two polymorphic alleles, NAU6202<sub>-380 </sub>in Hai7124 and NAU6202<sub>-400 </sub>in TM-1, which subsequently mapped on chromosome A7 (Figure ##FIG##5##6##, Additional file ##SUPPL##0##1##). Since NAU6202 was designed based on sequence information from AC190263 in Maxxa and its expected product size was 349 bp, close to the monomorphic allele size in <italic>G. raimondii</italic>, we concluded that the BAC was anchored into the D-subgenome. Finally, two loci were amplified by SSR primer pairs, which exhibited the same molecular size between TM-1 (Maxxa) and Hai7124 for each locus. In addition, the two loci were individually amplified in their diploid progenitors <italic>G. herbaceum </italic>and <italic>G. raimondii </italic>with almost the same fragment size as in the corresponding tetraploid subgenome. The sub-genome composition of the amplified product was confirmed by comparing its product with diploid <italic>G. herbaceum </italic>and <italic>G. raimondii </italic>(Figure ##FIG##9##10C##). For example, NAU6465 SSR primer pairs amplified two monomorphic loci with product sizes of 170 bp and 150 bp, respectively in TM-1 and Hai7124. Two fragments of 150 bp and 170 bp respectively from diploid <italic>G. herbaceum </italic>and <italic>G. raimondii </italic>were also generated. NAU6465 SSR primer pairs were developed based on sequence information from AC188140 in Maxxa with an expected size of 171 bp (Additional file ##SUPPL##0##1##). <italic>G. raimondii </italic>produced a 170 bp fragment very near the expected 171 bp. Therefore, we determined that the AC188140 clone was anchored in the D-subgenome.</p> <p>According to the criteria described above, 37 clones were truly anchored to the A- and D-subgenomes. Twelve BACs belonged to the A-subgenome and 25 to the D-subgenome (Table ##TAB##3##4##). The remaining 33 BACs could not be anchored to their subgenome because two distinct amplified loci corresponding to A- or D- subgenome were not resolved in tetraploid cotton. Therefore, the need for further experimentation, such as BAC-FISH analysis was recommended. Interestingly, among 25 BACs belonging to the D-subgenome, 19 BAC taggings coincided with SSR mapping. However, among 12 BACs within the A-subgenome, four (AC187470; AC187836; AC194383 and AC202821) were not consistent with SSR mapping. The four BACs were all anchored on the A-subgenome, but more than two polymorphic SSR loci from the same BAC were mapped on D-subgenome chromosomes (Table ##TAB##3##4##, Additional file ##SUPPL##0##1##). These results indicated that different evolutionary pressures acted on the A-subgenomes and D-subgenomes among different tetraploid cotton species in their corresponding homoeologous loci, and that the D-subgenomes exhibited more rapid evolutionary rates with increased nucleotide and allelic diversity than the A-subgenomes.</p> <title>Comparative sequence analysis between the A- and D-subgenome chromosomes</title> <p>Thirty-seven BACs with verified origins were identified in this study, 12 BACs belonging to the A-subgenome with a total length of 1,200,814 bp length including 69 gaps (average 5.75 gaps/BACs); and 25 BACs within the D-subgenome covering 2,374,313 bp length with 37 gaps (average 1.48 gaps/BACs). These results indicated that A-subgenome BACs possessed regions more difficult to sequence than those from the D-subgenome. Furthermore, the genes predicted from the 37 BACs were evaluated for possible intron size contributions that correlated with genome size between the A- and D-subgenome chromosomes. In the 12 BACs belonging to the A-subgenome, 67 genes were predicted with an average of 937 bp exons and 920 bp introns for each gene; however, in the 25 BACs belonging to the D-subgenome, 104 genes were predicted with an average of 1297 bp exons and 1414 bp introns for each gene. Therefore, introns might have no contribution to different subgenome size in <italic>Gossypium</italic>.</p>
<title>Discussion</title> <title>Characteristics of genome structure in allotetraploid cotton</title> <p>Cotton is the world's most important natural textile fiber and a significant oilseed crop. Cotton fiber is also an outstanding single-cell model to study plant cell elongation, and cell wall and cellulose biosynthesis [##REF##11743074##26##]. Of all 50 cotton species, <italic>Gossypium hirsutum </italic>provides over 95% of the annual cotton crop worldwide. Elucidating the tetraploid cotton genome composition and structure, especially upland cotton, will vastly expand opportunities in cotton research and agronomic improvements worldwide. However, cotton possesses a complex genome so whole genome sequencing of tetraploid cotton represents a substantial challenge [##REF##18056866##13##]. The GeneTrek approach has been proposed as an efficient means to evaluate the general properties of any genome by annotating a small set of randomly selected BACs [##REF##15491923##19##,##REF##16357197##20##]. In maize, sequence analysis of 100 randomly selected BACs led to the prediction of 42,000–56,000 genes with at least 66% repetitive DNA [##REF##16339807##27##]. In addition, sequence analysis of 74 randomly selected BACs showed that the maize nuclear genome contains about 37,000 candidate genes and 5,500 truncated and probable pseudogenes. However, the distribution of genes and repetitive elements is uneven [##REF##17615239##21##]. In the present study, properties of the upland cotton genome, such as total gene number, amount and distribution of repetitive DNA, and gene distribution, were first predicted based on the annotation of 142 randomly sequenced BACs. Compared with a density of one gene every 7.5 kb in the <italic>CesA </italic>region of homoeologous BACs [##REF##15256507##14##], the <italic>AdhA </italic>region of homoeologous BACs exhibits one gene per 20 kb for the A-subgenome and one gene every 13 kb for the D-subgenome [##REF##17461788##18##]. These data led to the prediction of more than 70,000 genes with one gene per 34.5 kb in upland cotton. Because upland cotton is an allotetraploid and has duplicated copies of genes in homoeologous regions of the A- and D-subgenomes, approximately 35,000 genes were predicted in each subgenome. In tetraploid cotton, the distribution of genes is uneven, with gene-rich and gene-free regions. We also found 21% of BACs lacked genes and 72.5% of the gene islands contained only one gene. These results indicated that selecting only gene-rich BACs for cotton genome sequencing is not adequate to cover the entire genome, owing to the fact that more than one fifth of BACs exhibit an absence of genes.</p> <p>In this study, 1,653 predicted gene models lacked homology to other species in the NCBI protein database. In addition, we verified 208 ESTs by BLASTN queries against the cotton EST database. However, we could not confirm if these transcripts were related to mobile elements, gene candidates, or special products in cotton. Therefore, we have not used the information to predict the structure and composition of the upland cotton genome. However, the functions and properties of these transcripts warrant further study to enhance the understanding of the complex upland cotton genome.</p> <title>Structure difference between A- and D-subgenome chromosomes</title> <p>In plants, the following factors have been summarized as the main mechanisms for genome size expansion: (1) long terminal repeat (LTR) retrotransposable element amplification and insertion such as that in maize [##UREF##11##28##]; (2) variation in intron size [##REF##10454621##29##]; (3) expansion of tandemly repetitive DNA sequences [##REF##11799393##30##]; (4) segmental duplications [##REF##10688139##31##]; (5) accumulation of pseudogenes [##UREF##12##32##]; and (6) transfer of organellar DNA to the nucleus [##REF##14615181##33##]. The cultivated cotton species <italic>Gossypium hirsutum </italic>has long been known as an allotetraploid possessing a nuclear A- and D-subgenome. A- and D-genome species diverged from a common ancestor approximately 5–10 Mya and acquired genomes that differ nearly twofold in size [##UREF##1##2##]. Based on the putative mechanisms of genome size expansion described above, it is uncertain which of the mechanism(s) played an important role in the composition and structure of the tetraploid cotton genomes. To explore this question, several studies have been initiated through comparative sequence analysis of specific genomic regions or by application of more global approaches [##REF##15256507##14##,##REF##16954538##16##,##REF##17461788##18##]. Grover et al. (2004)[##REF##15256507##14##] investigated A- and D-genome size evolution from tetraploid cotton in a 104 kb contiguous sequence surrounding the <italic>CesA1 </italic>gene, and demonstrated no evidence of genome size variation between the A- and D-subgenome genic regions. In a similar study, Grover et al. (2007)[##REF##17461788##18##] obtained the aligned length surrounding the <italic>AdhA </italic>gene with 101.7 kb in the A-subgenome, 49 kb in the D-subgenome, 112.3 kb from the diploid A-genome and 55 kb from the diploid D-genome. The results revealed the aligned length size variation was mainly attributed to differential accumulation of retroelements. Hawkins et al. (2006)[##REF##16954538##16##] compared diploid A- and D-genome size differences by utilizing the whole genome shotgun (WGS) method and concluded that 40%–65% of each genome is composed of transposable elements, with <italic>Copia</italic>-like sequences accumulated in smaller genomes and <italic>Gypsy</italic>-like sequences in larger genomes.</p> <p>Based on the sequence analysis of 37 subgemone-known BACs, we found no relationship between introns and different subgenome size in <italic>Gossypium</italic>. However, an average of 5.75 gaps/BAC indicated an increased number of gaps, lending difficulty to BAC assembly in the A-subgenome. The D-subgenome had an average of 1.48 gaps/BAC, demonstrating that BACs from the A-subgenome are more difficult for sequence assembly than those from the D-subgenome. This and previous studies revealed the presence of homeolog sequence and structure conservation in gene-rich regions, suggesting large amplification of transposable elements may not be in gene-rich regions, but may reside in low-density gene or gene-free regions. In future studies, the structure and function of DNA sequences in these gap regions can be confirmed by whole BAC sequence assembly analysis; and A-specific and D-specific regions related with transposable elements can be located using combined BAC-FISH technology.</p> <title>The D-subgenome has a more rapid evolutionary rate in different tetraploid cotton species</title> <p>Sequence and marker analyses from several previous studies indicated that varied evolutionary pressures might act on the D-subgenomes from different tetraploid cotton species. In both <italic>G. hirsutum </italic>and <italic>G. barbadense</italic>, the D-subgenome maintained greater nucleotide and allelic diversity than did the A-subgenome, results supported by duplicated paralogous <italic>Adh </italic>loci comparisons [##REF##11961094##34##,##REF##10331275##35##]. In addition, <italic>G. raimondii</italic>-derived EST-SSR markers had high polymorphic frequencies between <italic>G. hirsutum </italic>and <italic>G. barbadense </italic>[##REF##17409069##25##]. In this paper, we investigated whether BACs were characterized by an A- or D-subgenome. SSR marker BACs were largely tagged in the D-subgenome determined by integration of polymorphic marker loci with our tetraploid cotton backbone linkage groups. Our results further confirmed previous studies where sequence and structure conservation of homeologs between the A- and D- subgenomes was high. These data are consistent with the evolutionary history of tetraploid cotton progenitors, where diploid A- and D-genome species were derived from the same ancestor approximately 5–10 Mya. Alternatively, relaxed selection acted on the D-subgenomes from different tetraploid cotton species, evidenced by greater DNA sequence diversity among D-subgenomes than A-subgenomes in different tetraploid cotton species.</p>
<title>Conclusion</title> <p>The study provided us the first glimpse at cotton genome complexity, and the results indicated that the gene distribution in cotton genome is uneven with gene-rich and gene-free regions, and rich in repetitive elements. This study will serve as a foundation for tetraploid cotton whole genome sequencing in the future.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Upland cotton has the highest yield, and accounts for &gt; 95% of world cotton production. Decoding upland cotton genomes will undoubtedly provide the ultimate reference and resource for structural, functional, and evolutionary studies of the species. Here, we employed GeneTrek and BAC tagging information approaches to predict the general composition and structure of the allotetraploid cotton genome.</p> <title>Results</title> <p>142 BAC sequences from <italic>Gossypium hirsutum </italic>cv. Maxxa were downloaded <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov"/> and confirmed. These BAC sequence analysis revealed that the tetraploid cotton genome contains over 70,000 candidate genes with duplicated gene copies in homoeologous A- and D-subgenome regions. Gene distribution is uneven, with gene-rich and gene-free regions of the genome. Twenty-one percent of the 142 BACs lacked genes. BAC gene density ranged from 0 to 33.2 per 100 kb, whereas most gene islands contained only one gene with an average of 1.5 genes per island. Retro-elements were found to be a major component, first an enriched LTR/gypsy and second LTR/copia. Most LTR retrotransposons were truncated and in nested structures. In addition, 166 polymorphic loci amplified with SSRs developed from 70 BAC clones were tagged on our backbone genetic map. Seventy-five percent (125/166) of the polymorphic loci were tagged on the D-subgenome. By comprehensively analyzing the molecular size of amplified products among tetraploid <italic>G. hirsutum </italic>cv. Maxxa, acc. TM-1, and <italic>G. barbadense </italic>cv. Hai7124, and diploid <italic>G. herbaceum </italic>var. <italic>africanum </italic>and <italic>G. raimondii</italic>, 37 BACs, 12 from the A- and 25 from the D-subgenome, were further anchored to their corresponding subgenome chromosomes. After a large amount of genes sequence comparison from different subgenome BACs, the result showed that introns might have no contribution to different subgenome size in <italic>Gossypium</italic>.</p> <title>Conclusion</title> <p>This study provides us with the first glimpse of cotton genome complexity and serves as a foundation for tetraploid cotton whole genomesequencing in the future.</p>
<title>Abbreviations</title> <p>BAC: Bacterial Artificial Chromosome; FISH: Fluorescence In Situ Hybridization; QTL: Quantitative Trait Loci; LTR: Long Tandem Repeat; TSD: Target Site Duplications; SSR: Simple Sequence Repeat; WGS: Whole Genome Shotgun.</p> <title>Authors' contributions</title> <p>Experiments were designed by WZG with suggestions from TZZ. Experiments were performed by WZG and CPC. Bioinformatics analyses were performed by CBxW. LZ and LW helped with data analyses. WZG drafted the manuscript and TZZ revised the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We gratefully acknowledge the Maize Sequencing Consortium in USA, including the University of Arizona, Cold Spring Harbor Laboratory, Iowa State University, and the Genome Sequencing Center at Washington University School of Medicine in St. Louis for their free release of cotton BACs sequence data. This program was financially supported in part by National Science Foundation in China (30671324, 30730067), the Program for New Century Excellent Talents in University (NCET-04-0500), Jiangsu Natural Science Foundation, China (BK2007719), and the Program for 111 project (B08025).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Gene density variation among BACs.</bold> The BACs are sorted by overall gene density.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Gene distribution in gene islands.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>A updated genetic map of A1/D1 and A6/D6 homoeologous pairs</bold>. Note: Genetic map was constructed using a BC<sub>1 </sub>population obtained from the interspecific cross: <italic>G. hirsutum </italic>L. acc. TM-1 × <italic>G. barbadense </italic>L. cv. Hai7124. Chromosomes and linkage groups are arranged by homoeologous pairs and their corresponding conventional chromosome numbers denoted in bracket. Positions of loci are given in centiMorgans. Fragment sizes from Hai7124 allele (in base pair) are given next to the marker name, and the marker loci unmarked fragment size indicate tagged genes unpublished. Deviated loci are underlined. All SSR-derived BAC clones markers are indicated in bold. Homoeologous loci identified are connected by a bar.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>A updated genetic map of A2/D2 and A3/D3 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>A updated genetic map of A4/D4 and A5/D5 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>A updated genetic map of A7/D7 and A8/D8 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##. Deviated interval in A7 and D7 is boxed.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>A updated genetic map of A9/D9 and A10/D10 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>A updated genetic map of A11/D11 and A12/D12 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p><bold>A updated genetic map of A13/D13 homoeologous pairs</bold>. All legends are same as described for Figure ##FIG##2##3##.</p></caption></fig> <fig position="float" id="F10"><label>Figure 10</label><caption><p><bold>Identification of three BACs (AC187848, AC190263 and AC188140) belongings by amplified size analysis</bold>. Note: 1. <italic>G. hirsutum </italic>cv. Maxxa; 2. <italic>G. hirsutum </italic>acc. TM-1; 3. <italic>G. barbadense </italic>cv. Hai7124; 4. <italic>G. herbaceum </italic>var. <italic>africanum</italic>; 5. <italic>G. raimondii </italic>Arrow: amplified product with expected size from Maxxa. A: NAU6627 for AC187848 with expected size 247 bp. SSR tagging position was consistent with BAC clone belongings. B: NAU6202 for AC190263 with expected size 349 bp. BAC clone belongings is in homoeologous chromosome of SSR tagging position. C: NAU6465 for AC188140 with expected size 171 bp. BAC clone belongings is confirmed by comparing amplified product size in tetraploid with that in their diploid progenitors <italic>G. herbaceum </italic>and <italic>G. raimondii</italic>.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Summary of annotation results for 142 randomly selected cotton BACs</p></caption><table frame="hsides" rules="groups"><tbody><tr><td align="left">Total number of BACs analyzed</td><td align="left">142</td></tr><tr><td align="left">Combined BAC lengths</td><td align="left">14.2 Mb</td></tr><tr><td align="left">Amount of identified repetitive DNA (percentage)</td><td align="left">5.7 Mb (40.1%)</td></tr><tr><td align="left">Amount of unidentified DNA with predicted ORFs structure (Nos, percentage)</td><td align="left">1.4 Mb (1653, 9.9%)</td></tr><tr><td align="left">Unidentified ORFs that show collinearity with cotton EST database (percentage)</td><td align="left">208 (12.6%)</td></tr><tr><td align="left">Number of genes with similarity or collinearity support</td><td align="left">412</td></tr><tr><td align="left">Number of hypothetical genes with low similarity or collinearity support</td><td align="left">46</td></tr><tr><td align="left">Overall gene density</td><td align="left">One gene per 34.5 kb</td></tr><tr><td align="left">Number of estimated total cotton genes</td><td align="left">More than 70,000</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Types of Transposable elements in cotton genome</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left" colspan="2">Retroelements</td><td align="right">No. elements</td><td align="right">Length occupied (bp)</td><td align="right">Percentage</td></tr></thead><tbody><tr><td/><td align="left">Lines(L1/CIN4)</td><td align="right">30</td><td align="right">7242</td><td align="right">1.50%</td></tr><tr><td/><td align="left">LTR elements(Ty1/Copia)</td><td align="right">608</td><td align="right">384327</td><td align="right">31.20%</td></tr><tr><td/><td align="left">LTR elements(Gypsy/DIRS1)</td><td align="right">1192</td><td align="right">1052315</td><td align="right">61.10%</td></tr><tr><td/><td align="left">LTR elements(unclasiffied)</td><td align="right">2</td><td align="right">396</td><td align="right">0.10%</td></tr><tr><td/><td align="left">Total</td><td align="right">1832</td><td align="right">1444280</td><td align="right">93.90%</td></tr><tr><td align="left" colspan="2">DNA transposons</td><td/><td/><td/></tr><tr><td/><td align="left">hobo-Activator</td><td align="right">6</td><td align="right">1525</td><td align="right">0.30%</td></tr><tr><td/><td align="left">En-Spm</td><td align="right">70</td><td align="right">16371</td><td align="right">3.60%</td></tr><tr><td/><td align="left">MUDR-IS905</td><td align="right">12</td><td align="right">961</td><td align="right">0.60%</td></tr><tr><td/><td align="left">Tourist/Harbinger</td><td align="right">5</td><td align="right">1504</td><td align="right">0.20%</td></tr><tr><td/><td align="left">Total</td><td align="right">93</td><td align="right">20361</td><td align="right">4.70%</td></tr><tr><td align="left">RC/Helitron</td><td/><td align="right">4</td><td align="right">238</td><td align="right">0.20%</td></tr><tr><td align="left">Unclassified</td><td/><td align="right">22</td><td align="right">3455</td><td align="right">1.10%</td></tr><tr><td align="left">Total</td><td/><td align="right">1951</td><td/><td align="right">100%</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Repetitive elements in 7 gene-free BACs</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">BACs</td><td align="left">Length (bp)</td><td align="left">No. mobile elements</td><td align="left">Gypsy-like</td><td align="left">Copia-like</td><td align="left">RC/helitron</td><td align="left">No. intact LTRs</td></tr></thead><tbody><tr><td align="left">AC187206</td><td align="left">87776</td><td align="left">22</td><td align="left">16</td><td align="left">6</td><td/><td align="left">2</td></tr><tr><td align="left">AC187396</td><td align="left">95093</td><td align="left">16</td><td align="left">8</td><td align="left">7</td><td align="left">1</td><td align="left">0</td></tr><tr><td align="left">AC187849</td><td align="left">103016</td><td align="left">55</td><td align="left">55</td><td/><td/><td align="left">5</td></tr><tr><td align="left">AC189743</td><td align="left">110888</td><td align="left">17</td><td align="left">11</td><td align="left">6</td><td/><td align="left">2</td></tr><tr><td align="left">AC189931</td><td align="left">103756</td><td align="left">28</td><td align="left">28</td><td/><td/><td align="left">2</td></tr><tr><td align="left">AC190814</td><td align="left">103822</td><td align="left">18</td><td align="left">18</td><td/><td/><td align="left">1</td></tr><tr><td align="left">AC194319</td><td align="left">77717</td><td align="left">40</td><td align="left">40</td><td/><td/><td align="left">6</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>The subgenome belongings of BAC clones</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">BACs</td><td align="left">Chro. Tagged*</td><td align="left">Subgenome**</td><td align="left">BACs</td><td align="left">Chro. Tagged*</td><td align="left">Subgenome**</td></tr></thead><tbody><tr><td align="left">AC187225</td><td align="left">A2/D2</td><td align="left">A</td><td align="left">AC187200</td><td align="left">D13</td><td align="left">D</td></tr><tr><td align="left">AC187470</td><td align="left">D11</td><td align="left">A</td><td align="left">AC187214</td><td align="left">D11</td><td align="left">D</td></tr><tr><td align="left">AC187478</td><td align="left">A10/D10</td><td align="left">A</td><td align="left">AC187471</td><td align="left">D6</td><td align="left">D</td></tr><tr><td align="left">AC187836</td><td align="left">D5</td><td align="left">A</td><td align="left">AC187578</td><td align="left">A7/D7</td><td align="left">D</td></tr><tr><td align="left">AC188017</td><td align="left">A12/D12</td><td align="left">A</td><td align="left">AC187810</td><td align="left">D11</td><td align="left">D</td></tr><tr><td align="left">AC188842</td><td align="left">A2/D2</td><td align="left">A</td><td align="left">AC187794</td><td align="left">D11</td><td align="left">D</td></tr><tr><td align="left">AC193383</td><td align="left">A11</td><td align="left">A</td><td align="left">AC187848</td><td align="left">D11</td><td align="left">D</td></tr><tr><td align="left">AC194364</td><td align="left">A13/D13</td><td align="left">A</td><td align="left">AC188026</td><td align="left">D9</td><td align="left">D</td></tr><tr><td align="left">AC194383</td><td align="left">D10</td><td align="left">A</td><td align="left">AC188035</td><td align="left">D10</td><td align="left">D</td></tr><tr><td align="left">AC202821</td><td align="left">D11</td><td align="left">A</td><td align="left">AC188140</td><td align="left">D10</td><td align="left">D</td></tr><tr><td align="left">AC202822</td><td align="left">A9</td><td align="left">A</td><td align="left">AC188200</td><td align="left">D10</td><td align="left">D</td></tr><tr><td align="left">AC193751</td><td align="left">unknown</td><td align="left">A</td><td align="left">AC188401</td><td align="left">D8</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC188760</td><td align="left">D2</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC189748</td><td align="left">D3</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC190263</td><td align="left">A7/D7</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC190279</td><td align="left">D9</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC193505</td><td align="left">A5/D4</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC193514</td><td align="left">D1</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC193940</td><td align="left">D9</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC193994</td><td align="left">D6</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC197182</td><td align="left">A5/D5</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC202830</td><td align="left">D11</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC202831</td><td align="left">D5</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC187545</td><td align="left">unknown</td><td align="left">D</td></tr><tr><td/><td/><td/><td align="left">AC187548</td><td align="left">unknown</td><td align="left">D</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Summary of SSR primers developed from BACs generated from <italic>Gossypium hirsutum </italic>cv. Maxxa. The data provided all the markers' information developed from 142 BACs generated from <italic>Gossypium hirsutum </italic>cv. Maxxa.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Gene features predicted in 142 cotton BACs. The data provided all gene features predicted in 142 cotton BACs.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>142 BACs clones information. The data provided 142 BACs clones information in detail.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Tandem duplication genes in 30 gene islands. The data provided the ontology analysis of tandem duplication genes in 30 gene islands. F: molecular function; C: cellular component; P: biological process.</p></caption></supplementary-material>
<table-wrap-foot><p>*The chromosome tagging position of polymorphic SSR loci that SSR primer pairs were developed from corresponding BAC clone.</p><p>** The subgenome belongings of corresponding BAC clone.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-314-1"/> <graphic xlink:href="1471-2164-9-314-2"/> <graphic xlink:href="1471-2164-9-314-3"/> <graphic xlink:href="1471-2164-9-314-4"/> <graphic xlink:href="1471-2164-9-314-5"/> <graphic xlink:href="1471-2164-9-314-6"/> <graphic xlink:href="1471-2164-9-314-7"/> <graphic xlink:href="1471-2164-9-314-8"/> <graphic xlink:href="1471-2164-9-314-9"/> <graphic xlink:href="1471-2164-9-314-10"/>
<media xlink:href="1471-2164-9-314-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-314-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-314-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-314-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Fryxell"], "given-names": ["PA"], "article-title": ["A revised taxonomic interpretation of "], "italic": ["Gossypium "], "source": ["Rheedea"], "year": ["1992"], "volume": ["2"], "fpage": ["108"], "lpage": ["165"]}, {"surname": ["Wendel", "Cronn"], "given-names": ["JF", "RC"], "article-title": ["Polyploidy and the evolutionary history of cotton"], "source": ["Adv Agron"], "year": ["2003"], "volume": ["78"], "fpage": ["139"], "lpage": ["186"]}, {"article-title": ["National Cotton Council, USA"]}, {"surname": ["Kohel", "Yu", "Park", "Lazo"], "given-names": ["RJ", "J", "Y-H", "GR"], "article-title": ["Molecular mapping and characterization of traits controlling fiber quality in cotton"], "source": ["Euphytica"], "year": ["2001"], "volume": ["121"], "fpage": ["163"], "lpage": ["172"]}, {"surname": ["Park", "Alabady", "Ulloa", "Sickler", "Wilkins", "Yu", "Stelly", "Kohel", "el-Shihy", "Cantrell"], "given-names": ["YH", "MS", "M", "B", "TA", "J", "DM", "RJ", "OM", "RG"], "article-title": ["Genetic mapping of new cotton fiber loci using EST-derived microsatellites in an interspecific recombinant inbred line cotton population"], "source": ["Mol Gen Genomics"], "year": ["2005"], "volume": ["274"], "fpage": ["428"], "lpage": ["441"]}, {"surname": ["Paterson", "Saranga", "Menz", "Jiang", "Wright"], "given-names": ["AH", "Y", "M", "CX", "RJ"], "article-title": ["QTL analysis of genotype \u00d7 environmental interactions affecting cotton fiber quality"], "source": ["Theor ApplGenet"], "year": ["2003"], "volume": ["106"], "fpage": ["384"], "lpage": ["396"]}, {"surname": ["Shen", "Guo", "Zhu", "Yuan", "Zhang"], "given-names": ["XL", "WZ", "XF", "YL", "TZ"], "article-title": ["Molecularmapping of QTLs for qualities in three diverse lines in Upland cottonusing SSR markers"], "source": ["Mol Breed"], "year": ["2005"], "volume": ["15"], "fpage": ["169"], "lpage": ["181"]}, {"surname": ["Saha", "Raska", "Stelly"], "given-names": ["S", "DA", "DM"], "article-title": ["Upland cotton ("], "italic": ["Gossypium hirsutum ", "G. tomentosum "], "sub": ["1 "], "source": ["J Cotton Sci"], "year": ["2006"], "volume": ["10"], "fpage": ["146"], "lpage": ["154"]}, {"article-title": ["U.S. Department of Energy Joint Genome Institutes"]}, {"article-title": ["NCBI"]}, {"surname": ["Yang", "Guo", "Zhang"], "given-names": ["C", "WZ", "TZ"], "article-title": ["Molecular mapping of QTL for "], "italic": ["Verticillium ", "G. barbadense "], "source": ["Plant Sci"], "year": ["2008"], "volume": ["174"], "fpage": ["290"], "lpage": ["298"]}, {"surname": ["SanMiguel", "Bennetzen"], "given-names": ["P", "JL"], "article-title": ["Evidence that a recent increase in maize genome size was caused by the massive amplification of intergene retrotransposons"], "source": ["Ann Bot (Lond)"], "year": ["1998"], "volume": ["82"], "fpage": ["37"], "lpage": ["44"]}, {"surname": ["Zhang"], "given-names": ["J"], "article-title": ["Evolution by gene duplication: An update"], "source": ["Trends Ecol Evol"], "year": ["2003"], "volume": ["18"], "fpage": ["292"], "lpage": ["298"]}, {"surname": ["Wang", "Guo", "Cai", "Zhang"], "given-names": ["CB", "WZ", "CP", "TZ"], "article-title": ["Characterization, development and exploitation of EST-derived microsatellites in "], "italic": ["Gossypium raimondii "], "source": ["Chin Sci Bull"], "year": ["2006"], "volume": ["51"], "fpage": ["557"], "lpage": ["561"]}, {"article-title": ["Primer3.0"]}, {"article-title": ["Repeatmasker"]}, {"article-title": ["LTR_finder software"]}, {"article-title": ["FGENESH"]}, {"article-title": ["The "], "italic": ["Arabidopsis thaliana "]}, {"article-title": ["UniProt Gene Ontology"]}]
BMC Genomics. 2008 Jul 1; 9:314