Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade, and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a “relay race”-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

Male reproduction in flowering plants is highly sensitive to high temperature (HT). To investigate molecular mechanisms of the response of cotton (Gossypium hirsutum) anthers to HT, a relatively complete comparative transcriptome analysis was performed during anther development of cotton lines 84021 and H05 under normal temperature and HT conditions. In total, 4,599 differentially expressed genes were screened; the differentially expressed genes were mainly related to epigenetic modifications, carbohydrate metabolism, and plant hormone signaling. Detailed studies showed that the deficiency in S-ADENOSYL-L-HOMOCYSTEINE HYDROLASE1 and the inhibition of methyltransferases contributed to genome-wide hypomethylation in H05, and the increased expression of histone constitution genes contributed to DNA stability in 84021. Furthermore, HT induced the expression of CASEIN KINASEI (GhCKI) in H05, coupled with the suppression of starch synthase activity, decreases in glucose level during anther development, and increases in indole-3-acetic acid (IAA) level in late-stage anthers. The same changes also were observed in Arabidopsis (Arabidopsis thaliana) GhCKI overexpression lines. These results suggest that GhCKI, sugar, and auxin may be key regulators of the anther response to HT stress. Moreover, PHYTOCHROME-INTERACTING FACTOR genes (PIFs), which are involved in linking sugar and auxin and are regulated by sugar, might positively regulate IAA biosynthesis in the cotton anther response to HT. Additionally, exogenous IAA application revealed that high background IAA may be a disadvantage for latestage cotton anthers during HT stress. Overall, the linking of HT, sugar, PIFs, and IAA, together with our previously reported data on GhCKI, may provide dynamic coordination of plant anther responses to HT stress.

Male sterility caused by long-term high-temperature (HT) stress occurs widely in crops. MicroRNAs (miRNAs), a class of endogenous non-coding small RNAs, play an important role in the plant response to various abiotic stresses. To dissect the working principle of miRNAs in male sterility under HT stress in cotton, a total of 112 known miRNAs, 270 novel miRNAs and 347 target genes were identified from anthers of HT-insensitive (84021) and HT-sensitive (H05) cotton cultivars under normal-temperature and HT conditions through small RNA and degradome sequencing. Quantitative reverse transcriptase-polymerase chain reaction and 5'-RNA ligase-mediated rapid amplification of cDNA ends experiments were used to validate the sequencing data. The results show that miR156 was suppressed by HT stress in both 84021 and H05; miR160 was suppressed in 84021 but induced in H05. Correspondingly, SPLs (target genes of miR156) were induced both in 84021 and H05; ARF10 and ARF17 (target genes of miR160) were induced in 84021 but suppressed in H05. Overexpressing miR160 increased cotton sensitivity to HT stress seen as anther indehiscence, associated with the suppression of ARF10 and ARF17 expression, thereby activating the auxin response that leads to anther indehiscence. Supporting this role for auxin, exogenous Indole-3-acetic acid (IAA) leads to a stronger male sterility phenotype both in 84021 and H05 under HT stress. Cotton plants overexpressing miR157 suppressed the auxin signal, and also showed enhanced sensitivity to HT stress, with microspore abortion and anther indehiscence. Thus, we propose that the auxin signal, mediated by miRNAs, is essential for cotton anther fertility under HT stress.

High-temperature (HT) stress induces male sterility, leading to yield reductions in crops. DNA methylation regulates a range of processes involved in plant development and stress responses, but its role in male sterility under HT remains unknown. Here, we investigated DNA methylation levels in cotton () anthers under HT and normal temperature (NT) conditions by performing whole-genome bisulfite sequencing to investigate the regulatory roles of DNA methylation in male fertility under HT. Global disruption of DNA methylation, especially CHH methylation (where H = A, C, or T), was detected in an HT-sensitive line. Changes in the levels of 24-nucleotide small-interfering RNAs were significantly associated with DNA methylation levels. Experimental suppression of DNA methylation led to pollen sterility in the HT-sensitive line under NT conditions but did not affect the normal dehiscence of anther walls. Further transcriptome analysis showed that the expression of genes in sugar and reactive oxygen species (ROS) metabolic pathways were significantly modulated in anthers under HT, but auxin biosynthesis and signaling pathways were only slightly altered, indicating that HT disturbs sugar and ROS metabolism via disrupting DNA methylation, leading to microspore sterility. This study opens up a pathway for creating HT-tolerant cultivars using epigenetic techniques.

SummaryCotton fibre is an important natural fibre for the textile industry. The number of fibre initials determines the lint percentage, which is an important factor for cotton fibre yield. Although fibre development has been described by transcriptomic analysis, the mechanism by which the long noncoding RNA manipulates the initiation of lint and fuzz fibres remains unknown. In this study, three lines with different lint percentages were developed by crossing Xu142 with its fibreless mutant Xu142 fl. We collected the epidermal cells from the ovules with attached fibres at 0 and 5 days post anthesis (DPA) from Xu142, the fibreless mutant Xu142 fl and the three lint percent diversified lines for deep transcriptome sequencing. A total of 2641 novel genes, 35 802 long noncoding RNAs (lncRNAs) and 2262 circular RNAs (circRNAs) were identified, of which 645 lncRNAs were preferentially expressed in the fibreless mutant Xu142 fl and 651 lncRNAs were preferentially expressed in the fibre‐attached lines. We demonstrated the functional roles of the three lncRNAs in fibre development via a virus‐induced gene silencing (VIGS) system. Our results showed that silencing XLOC_545639 and XLOC_039050 in Xu142 fl increased the number of fibre initials on the ovules, but silencing XLOC_079089 in Xu142 resulted in a short fibre phenotype. This study established the transcriptomic repertoires in cotton fibre initiation and provided evidence for the potential functions of lncRNAs in fibre development.