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.

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.

Improving the yield and fiber quality of upland cotton is a goal of plant breeders. However, increasing the yield and quality of cotton fibers is becoming more urgent. While the growing human population needs more cotton fiber, climate change is reducing the amount of land on which cotton can be planted, or making it difficult to ensure that water and other resources will be available in optimal quantities. The most logical means of improving yield and quality is understanding and manipulating the genes involved. Here, we used comparative transcriptomics to explore differences in gene expression between long- and short-fiber cotton lines to identify candidate genes useful for cotton improvement. Light and electron microscopy revealed that the initial fiber density was significantly greater in our short-fiber group (SFG) than in our long-fiber group (LFG). Compared with the SFG fibers, the LFG fibers were longer at all developmental stages. Comparison of the LFG and SFG transcriptomes revealed a total of 3538 differentially expressed genes (DEGs). Notably, at all three developmental stages examined, two expression patterns, consistently downregulated (profile 0) and consistently upregulated (profile 7), were identified, and both were significantly enriched in the SFG and LFG. Twenty-two DEGs known to be involved in fiber initiation were detected in profile 0, while 31 DEGs involved in fiber elongation were detected in profile 7. Functional annotation suggested that these DEGs, which included ERF1, TUA2, TUB1, and PER64, affect fiber elongation by participating in the ethylene response, microtubule synthesis, and/or the peroxidase (POD) catalytic pathway. qRT-PCR was used to confirm the RNA sequencing results for select genes. A comparison of SFG and LFG transcription profiles revealed modest but important differences in gene expression between the groups. Notably, our results confirm those of previous studies suggesting that genes involved in ethylene, tubulin, and POD pathways play important roles in fiber development. The 22 consistently downregulated DEGs involved in fiber initiation and the 31 consistently upregulated genes involved in fiber elongation are seemingly good candidate genes for improving fiber initiation and elongation in cotton.

Verticillium wilt will seriously affect cotton yield and fiber quality. BEL1-Like transcription factors are involved in the regulation of secondary cell wall (SCW) formation, especially the biosynthesis of lignin that also plays a key role in cotton disease resistance. However, there is no report on the role of BEL1-Like transcription factor in the regulation of plant biological stress. In this study, tissue expression pattern analysis showed that a BEL1-Like transcription factor GhBLH7-D06 was predominantly expressed in vascular tissues and the SCW thickening stage of fiber development, while its expression could also respond to Verticillium dahliae infection and the phytohormone MeJA treatment, which indicated that GhBLH7-D06 might be involved in the defense response of Verticillium wilt. Using virus-induced gene silencing (VIGS) technology, we found silencing the expression of GhBLH7-D06 could enhance the resistance of cotton plants to Verticillium wilt, and the acquisition of resistance might be mainly due to the significant overexpression of genes related to lignin biosynthesis and JA signaling pathway, which also proves that GhBLH7-D06 negatively regulates the resistance of cotton to Verticillium wilt. Based on the results of yeast two-hybrid (Y2H) library screening and confirmation by bimolecular fluorescence complementary (BiFC) experiment, we found an Ovate Family Protein (OFP) transcription factor GhOFP3-D13 which was also a negative regulator of cotton Verticillium wilt resistance could that interacts with GhBLH7-D06. Furthermore, the dual-luciferase reporter assay and yeast one-hybrid (Y1H) experiment indicated that GhBLH7-D06 could target binding to the promoter region of GhPAL-A06 to suppress its expression and eventually lead to the inhibition of lignin biosynthesis. In general, the GhBLH7-D06/GhOFP3-D13 complex can negatively regulate resistance to Verticillium wilt of cotton by inhibiting lignin biosynthesis and JA signaling pathway.

Many BURP domain-containing proteins, which are unique to plants, have been identified. They performed diverse functions in plant development and the stress response. To date, only a few BURP domain-containing genes have been studied, and no comprehensive analysis of the gene family in cotton has been reported. In this study, 18, 17 and 30 putative BURP genes were identified in G. raimondii (D5), G. arboreum (A2) and G. hirsutum (AD1), respectively. These BURP genes were phylogenetically classified into eight subfamilies, which were confirmed by analyses of gene structures, motifs and protein domains. The uneven distribution of BURPs in chromosomes and gene duplication analysis indicated that segmental duplication might be the main driving force of the GhBURP family expansion. Promoter regions of all GhBURPs contained at least one putative stress-related cis-elements. Analysis of transcriptomic data and qRT-PCR showed that GhBURPs showed different expression patterns in different organs, and all of them, especially the members of the RD22-like subfamily, could be induced by different stresses, such as abscisic acid (ABA) and salicylic acid (SA), which indicated that the GhBURPs may performed important functions in cotton’s responses to various abiotic stresses. Our study comprehensively analyzed BURP genes in G. hirsutum, providing insight into the functions of GhBURPs in cotton development and adaptation to stresses.