Leaf senescence is an important developmental programmed degeneration process that dramatically affects crop quality and yield. The regulation of senescence is highly complex. Although senescence regulatory genes have been well characterized in model species such as Arabidopsis and rice, there is little information on the control of this process in cotton. Here, the senescence process in cotton (Gossypium hirsutum L.) leaves was investigated over a time course including young leaf, mature leaf and leaf samples from different senescence stages using RNA-Seq. Of 24,846 genes detected by mapping the tags to Gossypium genomes, 3,624 genes were identified as differentially expressed during leaf senescence. There was some overlap between the genes identified here and senescence-associated genes previously identified in other species. Most of the genes related to photosynthesis, chlorophyll metabolism and carbon fixation were downregulated; whereas those for plant hormone signal transduction were upregulated. Quantitative real-time PCR was used to evaluate the results of RNA-Seq for gene expression profiles. Furthermore, 519 differentially expressed transcription factors were identified, notably WRKY, bHLH and C3H. In addition, 960 genes involved in the metabolism and regulation of eight hormones were identified, of which many genes involved in the abscisic acid, brassinosteroid, jasmonic acid, salicylic acid and ethylene pathways were upregulated, indicating that these hormone-related genes might play crucial roles in cotton leaf development and senescence. However, most auxin, cytokinin and gibberellin pathway-related genes were downregulated, suggesting that these three hormones may act as negative regulators of senescence. This is the first high-resolution, multiple time-course, genome-wide comprehensive analysis of gene expression in cotton. These data are the most comprehensive dataset currently available for cotton leaf senescence, and will serve as a useful resource for unraveling the functions of many specific genes involved in cotton leaf development and senescence.

Summary: The transition from vegetative to reproductive growth is very important for early maturity in cotton. However, the genetic control of this highly dynamic and complex developmental process remains unclear. A high‐resolution tissue‐ and stage‐specific transcriptome profile was generated from six developmental stages using 72 samples of two early‐maturing and two late‐maturing cotton varieties. The results of histological analysis of paraffin sections showed that flower bud differentiation occurred at the third true leaf stage (3TLS) in early‐maturing varieties, but at the fifth true leaf stage (5TLS) in late‐maturing varieties. Using pairwise comparison and weighted gene co‐expression network analysis, 5312 differentially expressed genes were obtained, which were divided into 10 gene co‐expression modules. In the MElightcyan module, 46 candidate genes regulating cotton flower bud differentiation were identified and expressed at the flower bud differentiation stage. A novel key regulatory gene related to flower bud differentiation, GhCAL, was identified in the MElightcyan module. Anti‐GhCAL transgenic cotton plants exhibited late flower bud differentiation and flowering time. GhCAL formed heterodimers with GhAP1‐A04/GhAGL6‐D09 and regulated the expression of GhAP1‐A04 and GhAGL6‐D09. GhAP1‐A04‐ and GhAGL6‐D09‐silenced plants also showed significant late flowering. Finally, we propose a new flowering regulatory pathway mediated by GhCAL. This study elucidated the molecular mechanism of cotton flowering regulation and provides good genetic resources for cotton early‐maturing breeding.

Background: Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.). Results: In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers. Conclusion: We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

Background: Genetic male sterility (GMS) in cotton (Gossypium hirsutum) plays an important role in the utilization of hybrid vigor. However, the molecular mechanism of the GMS is still unclear. While numerous studies have demonstrated that microRNAs (miRNA) regulate flower and anther development, whether different small RNA regulations exist in GMS and its wild type is unclear. A deep sequencing approach was used to investigate the global expression and complexity of small RNAs during cotton anther development in this study. Results: Three small RNA libraries were constructed from the anthers of three development stages each from fertile wild type (WT) and its GMS mutant cotton, resulting in nearly 80 million sequence reads. The total number of miRNAs and short interfering RNAs in the three WT libraries was significantly greater than that in the corresponding three mutant libraries. Sixteen conserved miRNA families were identified, four of which comprised the vast majority of the expressed miRNAs during anther development. In addition, six conserved miRNA families were significantly differentially expressed during anther development between the GMS mutant and its WT. Conclusions: The present study is the first to deep sequence the small RNA population in G. hirsutum GMS mutant and its WT anthers. Our results reveal that the small RNA regulations in cotton GMS mutant anther development are distinct from those of the WT. Further results indicated that the differently expressed miRNAs regulated transcripts that were distinctly involved in anther development. Identification of a different set of miRNAs between the cotton GMS mutant and its WT will facilitate our understanding of the molecular mechanisms for male sterility.

Premature senescence can reduce the yield and quality of crops. WRKY transcription factors (TFs) play important roles during leaf senescence, but little is known about their ageing mechanisms in cotton. In this study, a group III WRKY TF, GhWRKY27, was isolated and characterized. The expression of GhWRKY27 was induced by leaf senescence and was higher in an early-ageing cotton variety than in a non-early-ageing cotton variety. Overexpression of GhWRKY27 in Arabidopsis promoted leaf senescence, as determined by reduced chlorophyll content and elevated expression of senescence-associated genes (SAGs). Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that GhWRKY27 interacted with an MYB TF, GhTT2. Putative target genes of GhWRKY27 were identified via chromatin immunoprecipitation followed by sequencing (ChIP-seq). Yeast one-hybrid (Y1H) assay and electrophoretic mobility shift assay (EMSA) revealed that GhWRKY27 binds directly to the promoters of cytochrome P450 94C1 (GhCYP94C1) and ripening-related protein 2 (GhRipen2–2). In addition, the expression patterns of GhTT2, GhCYP94C1 and GhRipen2–2 were identified during leaf senescence. Transient dual-luciferase reporter assay indicated that GhWRKY27 could activate the expression of GhCYP94C1 and GhRipen2–2. Our work lays the foundation for further study of the functional roles of WRKY genes during leaf senescence in cotton. In addition, our data provide new insights into the senescence-associated mechanisms of WRKY genes in cotton.