Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo–assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Trichomes originate from epidermal cells and can be classified as either glandular or non-glandular. Gossypium species are characterized by the presence of small and darkly pigmented lysigenous glands that contain large amounts of gossypol. Here, using a dominant glandless mutant, we characterize GoPGF, which encodes a basic helix-loop-helix domain-containing transcription factor, that we propose is a positive regulator of gland formation. Silencing GoPGF leads to a completely glandless phenotype. A single nucleotide insertion in GoPGF, introducing a premature stop codon is found in the duplicate recessive glandless mutant (gl₂gl₃). The characterization of GoPGF helps to unravel the regulatory network of glandular structure biogenesis, and has implications for understanding the production of secondary metabolites in glands. It also provides a potential molecular basis to generate glandless seed and glanded cotton to not only supply fibre and oil but also provide a source of protein for human consumption.