With radical global climate change and global warming, high temperature stress has become one of major factors exerting a major influence on crop production. In the cotton (Gossypium hirsutum L.)‐growing areas of China, especially in the Yangtze River valley, unexpected periodic episodes of extreme heat stress usually occur in July and August, the peak time of cotton flowering and boll loading, resulting in lower boll set and lint yield. Breeding programs for screening high temperature‐tolerant cotton germplasm and cultivars are urgent in order to stabilize yield in the current and future warmer weather conditions. In the present study, 14 cotton cultivars were quantified for in vitro pollen germination and pollen tube growth in response to temperatures ranging from 10 to 50 °C at 5 °C intervals. Different cotton genotypes varied in their in vitro pollen germination and pollen tube length responses to the different temperatures. Maximum pollen germination and pollen tube length ranged from 25.2% to 56.2% and from 414 to 682 μm, respectively. The average cardinal temperatures (T_min, T_opt, and T_max) also varied among the 14 cultivars and were 11.8, 27.3, and 42.7 °C for pollen germination and 11.8, 27.8, and 44.1 °C for maximum pollen tube length. Variations in boll retention and boll numbers per plant in field experiments were found for the 14 cotton cultivars and the boll retention and boll retained per plant on 20 August varied considerably in different years according to weather conditions. Boll retention on 20 August was highly correlated with maximum pollen germination (R²= 0.84) and pollen tube length (R²=0.64). A screening method based on principle component analysis of the combination of pollen characteristics in an in vitro experiment and boll retention testing in the field environment was used in the present study and, as a result, the 14 cotton cultivars could be classified as tolerant, moderately tolerant, moderately susceptible and susceptible to high temperature.(Managing editor: Li‐Hui Zhao)

MicroRNAs (miRNAs) are small molecule RNAs widely involved in responses to plant abiotic stresses. We performed small RNA sequencing of cotton anthers at four developmental stages under normal and high temperature (NT and HT, respectively) conditions to investigate the stress response characteristics of miRNA to HT. A total of 77 miRNAs, including 33 known miRNAs and 44 novel miRNAs, were identified, and 41 and 28 miRNAs were differentially expressed under NT and HT stress conditions, respectively. The sporogenous cell proliferation (SCP), meiotic phase (MP), microspore release period (MRP), and pollen maturity (PM) stages had 10 (including 12 miRNAs), four (including six miRNAs), four (including five miRNAs), and seven (including 11 miRNAs) HT stress-responsive miRNA families, respectively, which were identified after removing the changes in genotype-specific miRNAs under NT condition. Seven miRNA families (miR2949, miR167, and miR160 at the SCP stage; miR156 and miR172 at the MP stage; miR156 at the MRP stage; and miR393 and miR3476 at the PM stage), which had expression abundance of more than 10% of the total expression abundance, served as the main regulators responding to HT stress with positive or negative regulation patterns. These miRNAs orchestrated the expression of the corresponding target genes and led to different responses in the HT-tolerant and the HT-sensitive lines. The results revealed that the HT stress response of miRNAs in cotton anthers were stage-specific and differed with the development of anthers. Our study may enhance the understanding of the response of miRNAs to HT stress in cotton anthers and may clarify the mechanism of plant tolerance to HT stress.

SET (Su(var), E(z), and Trithorax) domain-containing proteins play an important role in plant development and stress responses through modifying lysine methylation status of histone. Gossypium raimondii may be the putative contributor of the D-subgenome of economical crops allotetraploid G. hirsutum and G. barbadense and therefore can potentially provide resistance genes. In this study, we identified 52 SET domain-containing genes from G. raimondii genome. Based on conserved sequences, these genes are grouped into seven classes and are predicted to catalyze the methylation of different substrates: GrKMT1 for H3K9me, GrKMT2 and GrKMT7 for H3K4me, GrKMT3 for H3K36me, GrKMT6 for H3K27me, but GrRBCMT and GrS-ET for nonhistones substrate-specific methylation. Seven pairs of GrKMT and GrRBCMT homologous genes are found to be duplicated, possibly one originating from tandem duplication and five from a large scale or whole genome duplication event. The gene structure, domain organization and expression patterns analyses suggest that these genes’ functions are diversified. A few of GrKMTs and GrRBCMTs, especially for GrKMT1A;1a, GrKMT3;3 and GrKMT6B;1 were affected by high temperature (HT) stress, demonstrating dramatically changed expression patterns. The characterization of SET domain-containing genes in G. raimondii provides useful clues for further revealing epigenetic regulation under HT and function diversification during evolution.