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 late-stage 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.

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 (Gossypium hirsutum) 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.

Summary: Global warming has reduced the productivity of many field‐grown crops, as the effects of high temperatures can lead to male sterility in such plants. Genetic regulation of the high temperature (HT) response in the major crop cotton is poorly understood. We determined the functionality and transcriptomes of the anthers of 218 cotton accessions grown under HT stress. By analyzing transcriptome divergence and implementing a genome‐wide association study (GWAS), we identified three thermal tolerance associated loci which contained 75 protein coding genes and 27 long noncoding RNAs, and provided expression quantitative trait loci (eQTLs) for 13 132 transcripts. A transcriptome‐wide association study (TWAS) confirmed six causal elements for the HT response (three genes overlapped with the GWAS results) which are involved in protein kinase activity. The most susceptible gene, GhHRK1, was confirmed to be a previously uncharacterized negative regulator of the HT response in both cotton and Arabidopsis. These functional variants provide a new understanding of the genetic basis for HT tolerance in male reproductive organs.

To understand the mechanisms of male sterility in cotton (Gossypium spp.), combined histological, biochemical and transcription analysis using RNA-Seq was carried out in the anther of the single-gene recessive genic male sterility system of male sterile line 1355A and male fertile line 1355B, which are near-isogenic lines (NILs) differing only in the fertility trait. A total of 2,446 differentially expressed genes were identified between the anthers of 1355AB lines, at three different stages of development. Cluster analysis and functional assignment of differentially expressed genes revealed differences in transcription associated with pollen wall and anther development, including the metabolism of fatty acids, glucose, pectin and cellulose. Histological and biochemical analysis revealed that a major cellular defect in the 1355A was a thicker nexine, consistent with the RNA-seq data and further gene expression studies implicated differences in fatty acids synthesis and metabolism. This study provides insight into the phenotypic characteristics and gene regulatory network of the genic male sterile line 1355A in upland cotton.