Gossypium hirsutum has proven difficult to sequence owing to its complex allotetraploid (A_tD_t) genome. Here we produce a draft genome using 181-fold paired-end sequences assisted by fivefold BAC-to-BAC sequences and a high-resolution genetic map. In our assembly 88.5% of the 2,173-Mb scaffolds, which cover 89.6% ∼ 96.7% of the A_tD_t genome, are anchored and oriented to 26 pseudochromosomes. Comparison of this G. hirsutum A_tD_t genome with the already sequenced diploid Gossypium arboreum (AA) and Gossypium raimondii (DD) genomes revealed conserved gene order. Repeated sequences account for 67.2% of the A_tD_t genome, and transposable elements (TEs) originating from D_t seem more active than from A_t. Reduction in the A_tD_t genome size occurred after allopolyploidization. The A or A_t genome may have undergone positive selection for fiber traits. Concerted evolution of different regulatory mechanisms for Cellulose synthase (CesA) and 1-Aminocyclopropane-1-carboxylic acid oxidase1 and 3 (ACO1,3) may be important for enhanced fiber production in G. hirsutum.

Upon assembling the firstGossypium herbaceum(A₁) genome and substantially improving the existingGossypium arboreum(A₂) andGossypium hirsutum((AD)₁) genomes, we showed that all existing A-genomes may have originated from a common ancestor, referred to here as A₀, which was more phylogenetically related to A₁than A₂. Further, allotetraploid formation was shown to have preceded the speciation of A₁and A₂. Both A-genomes evolved independently, with no ancestor–progeny relationship. Gaussian probability density function analysis indicates that several long-terminal-repeat bursts that occurred from 5.7 million years ago to less than 0.61 million years ago contributed compellingly to A-genome size expansion, speciation and evolution. Abundant species-specific structural variations in genic regions changed the expression of many important genes, which may have led to fiber cell improvement in (AD)₁. Our findings resolve existing controversial concepts surrounding A-genome origins and provide valuable genomic resources for cotton genetic improvement.

Genetic linkage maps play fundamental roles in understanding genome structure, explaining genome formation events during evolution, and discovering the genetic bases of important traits. A high-density cotton (Gossypium spp.) genetic map was developed using representative sets of simple sequence repeat (SSR) and the first public set of single nucleotide polymorphism (SNP) markers to genotype 186 recombinant inbred lines (RILs) derived from an interspecific cross between Gossypium hirsutum L. (TM-1) and G. barbadense L. (3-79). The genetic map comprised 2072 loci (1825 SSRs and 247 SNPs) and covered 3380 centiMorgan (cM) of the cotton genome (AD) with an average marker interval of 1.63 cM. The allotetraploid cotton genome produced equivalent recombination frequencies in its two subgenomes (At and Dt). Of the 2072 loci, 1138 (54.9%) were mapped to 13 At-subgenome chromosomes, covering 1726.8 cM (51.1%), and 934 (45.1%) mapped to 13 Dt-subgenome chromosomes, covering 1653.1 cM (48.9%). The genetically smallest homeologous chromosome pair was Chr. 04 (A04) and 22 (D04), and the largest was Chr. 05 (A05) and 19 (D05). Duplicate loci between and within homeologous chromosomes were identified that facilitate investigations of chromosome translocations. The map augments evidence of reciprocal rearrangement between ancestral forms of Chr. 02 and 03 versus segmental homeologs 14 and 17 as centromeric regions show homeologous between Chr. 02 (A02) and 17 (D02), as well as between Chr. 03 (A03) and 14 (D03). This research represents an important foundation for studies on polyploid cottons, including germplasm characterization, gene discovery, and genome sequence assembly.

Cultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) × 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits.

Alkaline stress is one of the abiotic stresses limiting cotton production. Though RNA-Seq analyses, have been conducted to investigate genome-wide gene expression in response to alkaline stress in plants, the response of sodium bicarbonate (NaHCO3) stress-related genes in cotton has not been reported. To explore the mechanisms of cotton response to this alkaline stress, we used next-generation sequencing (NGS) technology to study transcriptional changes of cotton under NaHCO3 alkaline stress. A total of 18,230 and 11,177 differentially expressed genes (DEGs) were identified in cotton roots and leaves, respectively. Gene ontology (GO) analysis indicated the enrichment of DEGs involved in various stimuli or stress responses. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs associated with plant hormone signal transduction, amino acid biosynthesis, and biosynthesis of secondary metabolites were regulated in response to the NaHCO3 stress. We further analyzed genes enriched in secondary metabolic pathways and found that secondary metabolites were regulated to eliminate the reactive oxygen species (ROS) and improve the cotton tolerance to the NaHCO3 stress. In this study, we learned that the toxic effect of NaHCO3 was more profound than that of NaOH at the same pH. Thus, Na+, HCO3- and pH had a great impact on the growth of cotton plant. The novel biological pathways and candidate genes for the cotton tolerance to NaHCO3 stress identified from the study would be useful in the genetic improvement of the alkaline tolerance in cotton.

High concentrations of heavy metals in the soil should be removed for environmental safety. Cadmium (Cd) is a heavy metal that pollutes the soil when its concentration exceeds 3.4 mg/kg. Although the potential use of cotton to remediate heavy Cd-polluted soils is known, little is understood about the molecular mechanisms of Cd tolerance. In this study, transcriptome analysis was used to identify Cd tolerance genes and their potential mechanisms in cotton. We exposed cotton plants to excess Cd and identified 4627 differentially expressed genes (DEGs) in the root, 3022 DEGs in the stem and 3854 DEGs in the leaves through RNA-Seq analysis. Among these genes were heavy metal transporter coding genes (ABC, CDF, HMA, etc.), annexin genes and heat shock genes (HSP), amongst others. Gene ontology (GO) analysis showed that the DEGs were mainly involved in the oxidation–reduction process and metal ion binding. The DEGs were mainly enriched in two pathways, the influenza A and pyruvate pathway. GhHMAD5, a protein containing a heavy-metal binding domain, was identified in the pathway to transport or to detoxify heavy metal ions. We constructed a GhHMAD5 overexpression system in Arabidopsis thaliana that showed longer roots compared to control plants. GhHMAD5-silenced cotton plants showed more sensitivity to Cd stress. The results indicate that GhHMAD5 is involved in Cd tolerance, which gives a preliminary understanding of the Cd tolerance mechanism in upland cotton. Overall, this study provides valuable information for the use of cotton to remediate soils polluted with Cd and potentially other heavy metals.