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.

ATP-binding cassette transporter G (ABCG) has been shown to be engaged in export of broad-spectrum compounds with structural differences, but little is known concerning its role in cutin formation of cotton (Gossypium spp.). In this study, we conduct a genome-wide survey and detected 69, 71, 124 and 131 ABCG genes within G. arboretum, G. raimondii, G. hirsutum and G. barbadense, separately. The above ABCGs could be divided into four groups (Ia, Ib, Ic, II). Some ABCG genes such as GhABCG15, whose homologous gene transports cuticular lipid in Arabidopsis, was preferentially expressed in the development of fiber. A weighted gene co-expression network analysis (WGCNA) demonstrated that GhABCG expression was significantly associated with the amount of 16-Hydroxypalmitate (a main component of cutin precursor) in cotton fibers. Further, silencing of GhABCG15 by virus-induced gene silencing (VIGS) in cotton generated brightened and crinkled leaves as well as reduced thickness of cuticle and increased permeability. Chemical composition analysis showed the cutin content in GhABCG15-silenced leaves had decreased while the wax content had increased. Our results provide an insight for better understanding of the role of the Gossypium ABCG family and revealed the essential role of GhABCGs in cotton cutin formation.