Single-Cell Transcriptome Atlas and Regulatory Dynamics in Developing Cotton Anthers.

Plant anthers are composed of different specialized cell types with distinct roles in plant reproduction. High temperature (HT) stress causes male sterility, resulting in crop yield reduction. However, the spatial expression atlas and regulatory dynamics during anther development and in response to HT remain largely unknown. Here, the first single-cell transcriptome atlas and chromatin accessibility survey in cotton anther are established, depicting the specific expression and epigenetic landscape of each type of cell in anthers. The reconstruction of meiotic cells, tapetal cells, and middle layer cell developmental trajectories not only identifies novel expressed genes, but also elucidates the precise degradation period of middle layer and reveals a rapid function transition of tapetal cells during the tetrad stage. By applying HT, heterogeneity in HT response is shown among cells of anthers, with tapetal cells responsible for pollen wall synthesis are most sensitive to HT. Specifically, HT shuts down the chromatin accessibility of genes specifically expressed in the tapetal cells responsible for pollen wall synthesis, such as QUARTET 3 (QRT3) and CYTOCHROME P450 703A2 (CYP703A2), resulting in a silent expression of these genes, ultimately leading to abnormal pollen wall and male sterility. Collectively, this study provides substantial information on anthers and provides clues for heat-tolerant crop creation.

Genome-Wide Identification and Analysis of the WNK Kinase Gene Family in Upland Cotton.

With-No-Lysine (WNK) kinases are a subfamily of serine/threonine protein kinases. WNKs are involved in plant abiotic stress response and circadian rhythms. However, members of the WNK subfamily and their responses to abiotic and biotic stresses in Gossypium hirsutum have not been reported. In this study, 26 GhWNKs were identified in G. hirsutum. The gene structure, conserved motifs, and upstream open reading frames (uORFs) of GhWNKs were identified. Moreover, GhWNKs regulation is predicted to be regulated by cis-acting elements, such as ABA responsive element (ABRE), MBS, and MYC. Furthermore, transcription factors including MIKC_MADS, C2H2, TALE, bZIP, Dof, MYB, bHLH, and HD-ZIP are projected to play a regulatory role in GhWNKs. The expression patterns of GhWNKs under normal conditions and biotic and abiotic stresses were evaluated, and their expression was found to vary. The expression patterns of several GhWNKs were induced by infiltration with Verticillium dahliae, suggesting that several GhWNKs may play important roles in the response of cotton to V. dahliae. Interestingly, a homoeologous expression bias within the GhWNKs was uncovered in upland cotton. Homoeologous expression bias within GhWNKs provides a framework to assist researchers and breeders in developing strategies to improve cotton traits by manipulating individual or multiple homeologs.

Robust All-Waterborne Superhydrophobic Coating with Photothermal Deicing and Passive Anti-icing Properties.

The compelling integration of superhydrophobic coatings with light-to-heat conversion capabilities has garnered substantial interest due to their dual functionality encompassing passive anti-icing and deicing attributes. However, the insufficient mechanical stability and the environmental and human health concerns stemming from the extensive use of organic solvents limit their practical application. In this study, an all-waterborne superhydrophobic photothermal coating (PCPAS) was prepared through the synergy of composite micro-nanoparticles derived from carbon nanotubes (CNT), polydopamine (PDA), and Ag particles with fluorine-containing polyacrylic emulsion (PFA). The PDA provided active sites for Ag+ reduction reaction and enhanced the interfacial interaction between CNT and Ag particles. The interfacial enhancement enabled the coating to maintain stable superhydrophobicity after 260 times sandpaper abrasion and 240 times tape peeling. Simultaneously, the composite micro-nanoparticle's light-to-heat conversion ability gave the coating excellent anti-icing/deicing capabilities. Under the condition of -20 °C, the freezing time of 30 µL of water droplets was extended to 392 s, and 2 × 2 × 2 cm ice cubes placed on the surface of the coating could completely melt after only 1142 s under simulated sunlight irradiation with a 1 kW/m2 intensity. In addition, the coating also had suitable self-cleaning properties and substrate applicability. The comprehensive attributes of this all-waterborne photothermal superhydrophobic coating render it a promising contender for anti-icing and deicing applications in challenging outdoor environments.

Genome-wide association study and transcriptome analysis reveal key genes controlling fruit branch angle in cotton.

Fruit branch angle (FBA), a pivotal component of cotton plant architecture, is vital for field and mechanical harvesting. However, the molecular mechanism of FBA formation is poorly understood in cotton. To uncover the genetic basis for FBA formation in cotton, we performed a genome-wide association study (GWAS) of 163 cotton accessions with re-sequencing data. A total of 55 SNPs and 18 candidate genes were significantly associated with FBA trait. By combining GWAS and transcriptome analysis, four genes underlying FBA were identified. An FBA-associated candidate gene Ghi_A09G08736, which is homologous to SAUR46 in Arabidopsis thaliana, was detected in our study. In addition, transcriptomic evidence was provided to show that gravity and light were implicated in the FBA formation. This study provides new insights into the genetic architecture of FBA that informs architecture breeding in cotton.

Cytochrome P450 mono-oxygenase CYP703A2 plays a central role in sporopollenin formation and ms5ms6 fertility in cotton.

The double-recessive genic male-sterile (ms) line ms5 ms6 has been used to develop cotton (Gossypium hirsutum) hybrids for many years, but its molecular-genetic basis has remained unclear. Here, we identified the Ms5 and Ms6 loci through map-based cloning and confirmed their function in male sterility through CRISPR/Cas9 gene editing. Ms5 and Ms6 are highly expressed in stages 7-9 anthers and encode the cytochrome P450 mono-oxygenases CYP703A2-A and CYP703A2-D. The ms5 mutant carries a single-nucleotide C-to-T nonsense mutation leading to premature chain termination at amino acid 312 (GhCYP703A2-A312aa ), and ms6 carries three nonsynonymous substitutions (D98E, E168K, and G198R) and a synonymous mutation (L11L). Enzyme assays showed that GhCYP703A2 proteins hydroxylate fatty acids, and the ms5 (GhCYP703A2-A312aa ) and ms6 (GhCYP703A2-DD98E,E168K,G198R ) mutant proteins have decreased enzyme activities. Biochemical and lipidomic analyses showed that in ms5 ms6 plants, C12-C18 free fatty acid and phospholipid levels are significantly elevated in stages 7-9 anthers, while stages 8-10 anthers lack sporopollenin fluorescence around the pollen, causing microspore degradation and male sterility. Overall, our characterization uncovered functions of GhCYP703A2 in sporopollenin formation and fertility, providing guidance for creating male-sterile lines to facilitate hybrid cotton production and therefore exploit heterosis for improvement of cotton.

Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton.

Flowering is a prerequisite for flowering plants to complete reproduction, and flowering time has an important effect on the high and stable yields of crops. However, there are limited reports on flowering-related genes at the genomic level in cotton. In this study, genomewide analysis of the evolutionary relationship of flowering-related genes in different cotton species shows that the numbers of flowering-related genes in the genomes of tetraploid cotton species Gossypium hirsutum and Gossypium barbadense were similar, and that these numbers were approximately twice as much as the number in diploid cotton species Gossypium arboretum. The classification of flowering-related genes shows that most of them belong to the photoperiod and circadian clock flowering pathway. The distribution of flowering-related genes on the chromosomes of the At and Dt subgenomes was similar, with no subgenomic preference detected. In addition, most of the flowering-related core genes in Arabidopsis thaliana had homologs in the cotton genome, but the copy numbers and expression patterns were disparate; moreover, flowering-related genes underwent purifying selection throughout the evolutionary and selection processes. Although the differentiation and reorganization of many key genes of the cotton flowering regulatory network occurred throughout the evolutionary and selection processes, most of them, especially those involved in the important flowering regulatory networks, have been relatively conserved and preferentially selected.

Genetic Mapping and Analysis of a Compact Plant Architecture and Precocious Mutant in Upland Cotton.

With the promotion and popularization of machine cotton-picking, more and more attention has been paid to the selection of early-maturity varieties with compact plant architecture. The type of fruit branch is one of the most important factors affecting plant architecture and early maturity of cotton. Heredity analysis of the cotton fruit branch is beneficial to the breeding of machine-picked cotton. Phenotype analysis showed that the types of fruit branches in cotton are controlled by a single recessive gene. Using an F2 population crossed with Huaxin102 (normal branch) and 04N-11 (nulliplex branch), BSA (Bulked Segregant Analysis) resequencing analysis and GhNB gene cloning in 04N-11, and allelic testing, showed that fruit branch type was controlled by the GhNB gene, located on chromosome D07. Ghnb5, a new recessive genotype of GhNB, was found in 04N-11. Through candidate gene association analysis, SNP 20_15811516_SNV was found to be associated with plant architecture and early maturity in the Xinjiang natural population. The GhNB gene, which is related to early maturity and the plant architecture of cotton, is a branch-type gene of cotton. The 20_15811516_SNV marker, obtained from the Xinjiang natural population, was used for the assisted breeding of machine-picked cotton varieties.

High temperature induces male sterility via MYB66-MYB4-Casein kinase I signaling in cotton.

High temperature (HT) causes male sterility and decreases crop yields. Our previous works have demonstrated that sugar and auxin signaling pathways, Gossypium hirsutum Casein kinase I (GhCKI), and DNA methylation are all involved in HT-induced male sterility in cotton. However, the signaling mechanisms leading to distinct GhCKI expression patterns induced by HT between HT-tolerant and HT-sensitive cotton anthers remain largely unknown. Here, we identified a GhCKI promoter (ProGhCKI) region that functions in response to HT in anthers and found the transcription factor GhMYB4 binds to this region to act as an upstream positive regulator of GhCKI. In the tapetum of early-stage cotton anthers, upregulated expression of GhMYB4 under HT and overexpressed GhMYB4 under normal temperature both led to severe male sterility phenotypes, coupled with enhanced expression of GhCKI. We also found that GhMYB4 interacts with GhMYB66 to form a heterodimer to enhance its binding to ProGhCKI. However, GhMYB66 showed an expression pattern similar to GhMYB4 under HT but did not directly bind to ProGhCKI. Furthermore, HT reduced siRNA-mediated CHH DNA methylations in the GhMYB4 promoter, which enhanced the expression of GhMYB4 in tetrad stage anthers and promoted the formation of the GhMYB4/GhMYB66 heterodimer, which in turn elevated the transcription of GhCKI in the tapetum, leading to male sterility. Overall, we shed light on the GhMYB66-GhMYB4-GhCKI regulatory pathway in response to HT in cotton anthers.

Fast anther dehiscence status recognition system established by deep learning to screen heat tolerant cotton.

BACKGROUND: From an economic perspective, cotton is one of the most important crops in the world. The fertility of male reproductive organs is a key determinant of cotton yield. Anther dehiscence or indehiscence directly determines the probability of fertilization in cotton. Thus, rapid and accurate identification of cotton anther dehiscence status is important for judging anther growth status and promoting genetic breeding research. The development of computer vision technology and the advent of big data have prompted the application of deep learning techniques to agricultural phenotype research. Therefore, two deep learning models (Faster R-CNN and YOLOv5) were proposed to detect the number and dehiscence status of anthers. RESULT: The single-stage model based on YOLOv5 has higher recognition speed and the ability to deploy to the mobile end. Breeding researchers can apply this model to terminals to achieve a more intuitive understanding of cotton anther dehiscence status. Moreover, three improvement strategies are proposed for the Faster R-CNN model, where the improved model has higher detection accuracy than the YOLOv5 model. We have made three improvements to the Faster R-CNN model and after the ensemble of the three models and original Faster R-CNN model, R2 of "open" reaches to 0.8765, R2 of "close" reaches to 0.8539, R2 of "all" reaches to 0.8481, higher than the prediction results of either model alone, which are completely able to replace the manual counting results. We can use this model to quickly extract the dehiscence rate of cotton anthers under high temperature (HT) conditions. In addition, the percentage of dehiscent anthers of 30 randomly selected cotton varieties were observed from the cotton population under normal conditions and HT conditions through the ensemble of the Faster R-CNN model and manual counting. The results show that HT decreased the percentage of dehiscent anthers in different cotton lines, consistent with the manual method. CONCLUSIONS: Deep learning technology have been applied to cotton anther dehiscence status recognition instead of manual methods for the first time to quickly screen HT-tolerant cotton varieties. Deep learning can help to explore the key genetic improvement genes in the future, promoting cotton breeding and improvement.

Degradation of de-esterified pctin/homogalacturonan by the polygalacturonase GhNSP is necessary for pollen exine formation and male fertility in cotton.

The pollen wall exine provides a protective layer for the male gametophyte and is largely composed of sporopollenin, which comprises fatty acid derivatives and phenolics. However, the biochemical nature of the external exine is poorly understood. Here, we show that the male sterile line 1355A of cotton mutated in NO SPINE POLLEN (GhNSP) leads to defective exine formation. The GhNSP locus was identified through map-based cloning and confirmed by genetic analysis (co-segregation test and allele prediction using the CRISPR/Cas9 system). In situ hybridization showed that GhNSP is highly expressed in tapetum. GhNSP encodes a polygalacturonase protein homologous to AtQRT3, which suggests a function for polygalacturonase in pollen exine formation. These results indicate that GhNSP is functionally different from AtQRT3, the latter has the function of microspore separation. Biochemical analysis showed that the percentage of de-esterified pectin was significantly increased in the 1355A anthers at developmental stage 8. Furthermore, immunofluorescence studies using antibodies to the de-esterified and esterified homogalacturonan (JIM5 and JIM7) showed that the Ghnsp mutant exhibits abundant of de-esterified homogalacturonan in the tapetum and exine, coupled with defective exine formation. The characterization of GhNSP provides new understanding of the role of polygalacturonase and de-esterified homogalacturonan in pollen exine formation.

A combination of genome-wide and transcriptome-wide association studies reveals genetic elements leading to male sterility during high temperature stress in cotton.

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.

High day and night temperatures distinctively disrupt fatty acid and jasmonic acid metabolism, inducing male sterility in cotton.

High temperature stress is an inevitable environmental factor in certain geographical regions. To study the effect of day and night high temperature stress on male reproduction, the heat-sensitive cotton line H05 was subjected to high temperature stress. High day/normal night (HN) and normal day/high night (NH) temperature treatments were compared with normal day/normal night (NN) temperature as a control. At the anther dehiscence stage, significant differences were observed, with a reduction in flower size and filament length, and sterility in pollen, seen in NH more than in HN. A total of 36 806 differentially expressed genes were screened, which were mainly associated with fatty acid and jasmonic acid (JA) metabolic pathways. Fatty acid and JA contents were reduced more in NH than HN. Under NH, ACYL-COA OXIDASE 2 (ACO2), a JA biosynthesis gene, was down-regulated. Interestingly, aco2 CRISPR-Cas9 mutants showed male sterility under the NN condition. The exogenous application of methyl jasmonate to early-stage buds of mutants rescued the sterile pollen and indehiscent anther phenotypes at the late stage. These data show that high temperature at night may affect fatty acid and JA metabolism in anthers by suppressing GhACO2 and generate male sterility more strongly than high day temperature.

Proteomic analysis reveals that sugar and fatty acid metabolisms play a central role in sterility of the male-sterile line 1355A of cotton.

Cotton (Gossypium spp.) is one of the most important economic crops and exhibits yield-improving heterosis in specific hybrid combinations. The genic male-sterility system is the main strategy used for producing heterosis in cotton. To better understand the mechanisms of male sterility in cotton, we carried out two-dimensional electrophoresis (2-DE) and label-free quantitative proteomics analysis in the anthers of two near-isogenic lines, the male-sterile line 1355A and the male-fertile line 1355B. We identified 39 and 124 proteins that were significantly differentially expressed between these two lines in the anthers at the tetrad stage (stage 7) and uninucleate pollen stage (stage 8), respectively. Gene ontology-based analysis revealed that these differentially expressed proteins were mainly associated with pyruvate, carbohydrate, and fatty acid metabolism. Biochemical analysis revealed that in the anthers of line 1355A, glycolysis was activated, which was caused by a reduction in fructose, glucose, and other soluble sugars, and that accumulation of acetyl-CoA was increased along with a significant increase in C14:0 and C18:1 free fatty acids. However, the activities of pyruvate dehydrogenase and fatty acid biosynthesis were inhibited and fatty acid ß-oxidation was activated at the translational level in 1355A. We speculate that in the 1355A anther, high rates of glucose metabolism may promote fatty acid synthesis to enable anther growth. These results provide new insights into the molecular mechanism of genic male sterility in upland cotton.

Transcriptome, cytological and biochemical analysis of cytoplasmic male sterility and maintainer line in CMS-D8 cotton.

Key message This research based on RNA-seq, biochemical, and cytological analyses sheds that ROS may serve as important signaling molecules of cytoplasmic male sterility in CMS-D8 cotton. To understand the mechanism of cytoplasmic male sterility in cotton (Gossypium hirsutum), transcriptomic, cytological, and biochemical analysis were performed between the cytoplasmic male sterility CMS-D8 line, Zhong41A, and its maintainer line Zhong41B. A total of 2335 differentially expressed genes (DEGs) were identified in the CMS line at three different stages of anther development. Bioinformatics analysis of these DEGs indicated their relationship to reactive oxygen species (ROS) homeostasis, including reduction-oxidation reactions and the metabolism of glutathione and ascorbate. At the same time, DEGs associated with tapetum development, especially the transition to secretory tapetum, were down-regulated in the CMS line. Biochemical analysis indicated that the ability of the CMS line to eliminate ROS was decreased, which led to the rapid release of H2O2. Cytological analysis revealed that the most crucial defect in the CMS line was the abnormal tapetum. All these results are consistent with the RNA sequencing data. On the basis of our findings, we propose that ROS act as signal molecules, which are released from mitochondria and transferred to the nucleus, triggering the formation of abnormal tapetum.

Promoters of Arabidopsis Casein kinase I-like 2 and 7 confer specific high-temperature response in anther.

KEY MESSAGE: (1) We systematically analyze the promoter activities of AtCKLs in various tissues; (2) AtCKL2 and AtCKL7 were expressed in early developmental anthers under high temperature (HT) conditions; (3) AtMYB24 may function as a positive regulator of AtCKL2 and AtCKL7 expression under HT. High temperature (HT) can seriously impede plant growth and development, causing severe loss of crop production. In Arabidopsis, AtCKL genes show high similarity to GhCKI, a gene reported to disrupt tapetal programmed cell death in cotton. However, most of AtCKL genes are not well characterized. Here, we systematically analyzed the expression patterns of AtCKLs in various tissues. The expression of AtCKL2 and AtCKL7 was induced in early anther development under HT, which is similar to the case of GhCKI. In silico analysis of AtCKL2 and AtCKL7 promoters indicated that four types of transcription factors (TFs) (MADS, NAC, WRKY and R2R3-MYB) might bind to AtCKL2 and AtCKL7 promoters. Furthermore, three MADS, three NAC, one WRKY, and three R2R3-MYB TFs were up-regulated in stage 1-8 anthers and three R2R3-MYB TFs were up-regulated in stage 9-10 anthers under HT, implying the important roles of R2R3-MYB genes in the response of anthers to HT. Among the R2R3-MYB genes, AtMYB24 showed the similar expression as AtCKL2 and AtCKL7 in the anthers under HT. Additionally, yeast one-hybrid and dual-luciferase reporter system assays verified that AtMYB24 could bind to AtCKL2 and AtCKL7 promoters and activate the expression of these two genes. In brief, this study provides the overall expression profiles of AtCKLs, useful information for unraveling the molecular mechanism of AtCKL2 and AtCKL7 gene expression in early anther development under HT, and important clues for elucidating the mechanism of transcriptional regulation of CKI genes in plant anther under HT, which are critical to the reduction of crop yield loss resulting from HT.

Disrupted Genome Methylation in Response to High Temperature Has Distinct Affects on Microspore Abortion and Anther Indehiscence.

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.

microRNAs involved in auxin signalling modulate male sterility under high-temperature stress in cotton (Gossypium hirsutum).

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.

Defective pollen wall contributes to male sterility in the male sterile line 1355A of cotton.

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.

Sugar and auxin signaling pathways respond to high-temperature stress during anther development as revealed by transcript profiling analysis in cotton.

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.