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.

Cotton (Gossypium spp.) is an important economic crop and there is obvious heterosis in cotton, fertility has played an important role in this heterosis. However, the genes that exhibit critical roles in anther development and fertility are not well understood. Here, we report an acyl-CoA N-acyltransferase (EC2.3; GhACNAT) that plays a key role in anther development and fertility. Suppression of GhACNAT by virus-induced gene silencing in transgenic cotton (G. hirsutum L. cv. C312) resulted in indehiscent anthers that were full of pollen, diminished filaments and stamens and plant sterility. We found GhACNAT was involved in lipid metabolism and jasmonic acid (JA) biosynthesis. The genes differentially expressed in GhACNAT-silenced plants and C312 were mainly involved in catalytic activity and transcription regulator activity in lipid metabolism. In GhACNAT-silenced plants, the expression levels of genes involved in lipid metabolism and jasmonic acid biosynthesis were significantly changed, the amount of JA in leaves and reproductive organs was significantly decreased compared with the amounts in C312. Treatments with exogenous methyl jasmonate rescued anther dehiscence and pollen release in GhACNAT-silenced plants and caused self-fertility. The GhACNAT gene may play an important role in controlling cotton fertility by regulating the pathways of lipid synthesis and JA biogenesis.

Here, cytosine methylation at single-base resolution across the whole genome of cotton (Gossypium hirsutum L.) anthers was mapped using the whole-genome bisulfite sequencing technique, and the methylome changes associated with high-temperature (HT) stress were analysed in two cotton lines of the CMS system with contrasting HT stress tolerance. The cotton anther genome was found to display approximately 31.6%, 68.7%, 61.8%, and 21.8% methylation across all sequenced C sites and in the CG, CHG and CHH sequence contexts, respectively. In an integrated global methylome and transcriptome analysis, only promoter-unmethylated genes showed higher expression levels than promoter-methylated genes, whereas gene body methylation presented an obvious positive correlation with gene expression. The methylation profiles of transposable elements in cotton anthers were characterized, and more differentially methylated transposable elements were demethylated under HT stress. HT-induced promoter methylation changes caused upregulated expression of the mitochondrial respiratory chain enzyme-associated genes GhNDUS7, GhCOX6A, GhCX5B2, and GhATPBM, ultimately promoting a series of redox processes to form ATP for normal anther development under HT stress. In vitro application of the common DNA methylation inhibitor 5-azacytidine and accelerator methyl trifluoromethanesulfonate demonstrated that DNA demethylation promoted anther development, while increased methylation only partially inhibited anther development under HT stress.