As calcium signal sensors, calcium-dependent protein kinases (CPKs) play vital roles in stimulating the production of secondary metabolites to participate in plant development and response to environmental stress. However, investigations of the Glycyrrhiza uralensis CPK family genes and their multiple functions are rarely reported. In this study, a total of 23 GuCPK genes in G. uralensis were identified, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, and promoter cis-acting elements were analyzed. Ten GuCPKs showed root-specific preferential expressions, and GuCPKs indicated different expression patterns under treatments of CaCl₂ and NaCl. In addition, under 2.5 mM of CaCl₂ and 30 mM of NaCl treatments, the diverse, induced expression of GuCPKs and significant accumulations of glycyrrhizic acid and flavonoids suggested the possible important function of GuCPKs in regulating the production of glycyrrhizic acid and flavonoids. Our results provide a genome-wide characterization of CPK family genes in G. uralensis, and serve as a foundation for understanding the potential function and regulatory mechanism of GuCPKs in promoting the biosynthesis of glycyrrhizic acid and flavonoids under salt stress.

Monodehydroasorbate reductase (MDHAR) (EC1.6.5.4), a key enzyme in ascorbate-glutathione recycling, plays important roles in cell growth, plant development and physiological response to environmental stress via control of ascorbic acid (AsA)-mediated reduction/oxidation (redox) regulation. Until now, information regarding MDHAR function and regulatory mechanism in Gossypium have been limited. Herein, a genome-wide identification and comprehensive bioinformatic analysis of 36 MDHAR family genes in four Gossypium species, Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, were performed, indicating their close evolutionary relationship. Expression analysis of GhMDHARs in different cotton tissues and under abiotic stress and phytohormone treatment revealed diverse expression features. Fiber-specific expression analysis showed that GhMDHAR1A/D, 3A/D and 4A/D were preferentially expressed in fiber fast elongating stages to reach peak values in 15-DPA fibers, with corresponding coincident observances of MDHAR enzyme activity, AsA content and ascorbic acid/dehydroascorbic acid (AsA/DHA) ratio. Meanwhile, there was a close positive correlation between the increase of AsA content and AsA/DHA ratio catalyzed by MDHAR and fiber elongation development in different fiber-length cotton cultivars, suggesting the potential important function of MDHAR for fiber growth. Following H₂O₂ stimulation, GhMDHAR demonstrated immediate responses at the levels of mRNA, enzyme, the product of AsA and corresponding AsA/DHA value, and antioxidative activity. These results for the first time provide a comprehensive systemic analysis of the MDHAR gene family in plants and the four cotton species and demonstrate the contribution of MDHAR to fiber elongation development by controlling AsA-recycling-mediated cellular redox homeostasis.

RNA sequencing (RNA-Seq)-based gene expression analysis is applicable to a wide range of biological purposes in various species. Reverse transcription quantitative PCR (RT-qPCR) is also used to assess target gene expression utilizing stably expressed reference genes as internal control under a given set of conditions. However, investigations of the reference genes for RT-qPCR normalization in the process of somatic embryogenesis (SE) initial dedifferentiation in Gossypium hirsutum are rarely reported. In this study, on the basis of our previous transcriptome data of three different induction stages during SE initial dedifferentiation process in four G. hirsutum cultivars that have different SE capability, 15 candidate genes were selected during SE initial dedifferentiation process, and their expression stability was evaluated by geNorm, NormFinder, and BestKeeper. The results indicated that the two genes of endonuclease 4 (ENDO4) and 18S ribosomal RNA (18S rRNA) showed stable expression in the four different G. hirsutum cultivars, endowing them to be appropriate reference genes during three induction stages in the four cotton cultivars. In addition, the stability and reliability of the two reference genes of ENDO4 and 18S rRNA were further verified by comparing the expressions of auxin-responsive protein 22 (AUX22) and ethylene-responsive transcription factor 17 (ERF17) between RT-qPCR results and the RNA-seq data, which showed strong positive correlation coefficient (R² = 0.8396–0.9984), validating again the steady expression of ENDO4 and 18S rRNA as the reliable reference genes. Our results provide effective reference genes for RT-qPCR normalization during SE process in different G. hirsutum cultivars.

Ascorbate oxidase (AO) and skewed5 (SKU5)-similar (SKS) proteins belong to the multicopper oxidase (MCO) family and play important roles in plants in response to environmental stress via modulation of oxidoreduction homeostasis. Currently, reports on the response of Gossypium barbadense MCO to Verticillium wilt (VW) caused by Verticillium dahliae are still limited. Herein, RNA sequencing of two G. barbadense cultivars of VW-resistant XH21 and VW-susceptible XH7 under V. dahliae treatment, combined with physiological and genetic analysis, was performed to analyze the function and mechanism of multicopper oxidases GbAO and GbSKS involved in V. dahliae resistance. The identified differentially expressed genes are mainly involved in the regulation of oxidoreduction reaction, and extracellular components and signaling. Interestingly, ascorbate oxidase family members were discovered as the most significantly upregulated genes after V. dahliae treatment, including GbAO3A/D, GbSKS3A/D, and GbSKS16A/D. H₂O₂ and Asc contents, especially reductive Asc in both XH21 and XH7, were shown to be increased. Silenced expression of respective GbAO3A/D, GbSKS3A/D, and GbSKS16A/D in virus-induced gene silencing (VIGS) cotton plants significantly decreased the resistance to V. dahliae, coupled with the reduced contents of pectin and lignin. Our results indicate that AO might be involved in cotton VW resistance via the regulation of cell wall components.

Chilling stress generates significant inhibition of normal growth and development of cotton plants and lead to severe reduction of fiber quality and yield. Currently, little is known for the molecular mechanism of brown-fiber cotton (BFC) to respond to chilling stress. Herein, RNA-sequencing (RNA-seq)-based comparative analysis of leaves under 4°C treatment in two different-tolerant BFC cultivars, chilling-sensitive (CS) XC20 and chilling-tolerant (CT) Z1612, was performed to investigate the response mechanism. A total of 72650 unigenes were identified with eight commonly used databases. Venn diagram analysis identified 1194 differentially expressed genes (DEGs) with significant up-regulation in all comparison groups. Furthermore, enrichment analyses of COG and KEGG, as well as qRT-PCR validation, indicated that 279 genes were discovered as up-regulated DEGs (UDEGs) with constant significant increased expression in CT cultivar Z1612 groups at the dimensions of both each comparison group and treatment time, locating in the enriched pathways of signal transduction, protein and carbohydrate metabolism, and cell component. Moreover, the comprehensive analyses of gene expression, physiological index and intracellular metabolite detections, and ascorbate antioxidative metabolism measurement validated the functional contributions of these identified candidate genes and pathways to chilling stress. Together, this study for the first time report the candidate key genes and metabolic pathways responding to chilling stress in BFC, and provide the effective reference for understanding the regulatory mechanism of low temperature adaptation in cotton.

Cotton is a crucial economic crop that supplies natural fibers for the textile industry, with fiber quality being greatly impacted by abiotic stress throughout its growth stages. The Golden2-Like (GLK) gene family plays a key role in plant development and adaptation to abiotic stress. However, the specific functions and regulatory mechanisms of GLK members in cotton remain largely unexplored. In this study, a thorough analysis of GLK in four cotton species (Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense) was conducted. A total of 198 GLK genes were identified in cotton. Conserved sequence analysis revealed that most GLK proteins contain two highly conserved domains: a MYB DNA-binding domain and a C-terminal (GCT) box. Promoter element analysis results show that the GLK gene family contains many stress response-related elements. Expression analysis demonstrated that GhGLK2, GhGLK11, GhGLK16, and GhGLK30 responded significantly to drought, salt, and temperature stresses. And GhGLK2, GhGLK13, GhGLK38, GhGLK42, and GhGLK46 responded significantly to cotton development. Yeast one-hybrid, yeast two-hybrid, and dual-luciferase assay results indicate that GhGLK2 interacts with GhGUN5, GhPIL6, GhNAC6, GhTPX2, and GhERF10. These findings suggest that these GhGLKs may play crucial roles in regulating the response to abiotic stress. Overall, this study provides a solid theoretical foundation for understanding the role of the GLK gene family in cotton’s response to abiotic stress.