Summary: Plants endure challenging environments in which they are constantly threatened by diverse pathogens. The soil‐borne fungus Verticillium dahliae is a devastating pathogen affecting many plant species including cotton, in which it significantly reduces crop yield and fiber quality. Melatonin involvement in plant immunity to pathogens has been reported, but the mechanisms of melatonin‐induced plant resistance are unclear. In this study, the role of melatonin in enhancing cotton resistance to V. dahliae was investigated. At the transcriptome level, exogenous melatonin increased the expression of genes in phenylpropanoid, mevalonate (MVA), and gossypol pathways after V. dahliae inoculation. As a result, lignin and gossypol, the products of these metabolic pathways, significantly increased. Silencing the serotonin N‐acetyltransferase 1 (GhSNAT1) and caffeic acid O‐methyltransferase (GhCOMT) melatonin biosynthesis genes compromised cotton resistance, with reduced lignin and gossypol levels after V. dahliae inoculation. Exogenous melatonin pre‐treatment prior to V. dahliae inoculation restored the level of cotton resistance reduced by the above gene silencing effects. Melatonin levels were higher in resistant cotton cultivars than in susceptible cultivars after V. dahliae inoculation. The findings indicate that melatonin affects lignin and gossypol synthesis genes in phenylpropanoid, MVA, and gossypol pathways, thereby enhancing cotton resistance to V. dahliae.

Glutathione S-transferases (GSTs) play versatile functions in multiple aspects of plant growth and development. A comprehensive genome-wide survey of this gene family in the genomes of G. raimondii and G. arboreum was carried out in this study. Based on phylogenetic analyses, the GST gene family of both two diploid cotton species could be divided into eight classes, and approximately all the GST genes within the same subfamily shared similar gene structure. Additionally, the gene structures between the orthologs were highly conserved. The chromosomal localization analyses revealed that GST genes were unevenly distributed across the genome in both G. raimondii and G. arboreum. Tandem duplication could be the major driver for the expansion of GST gene families. Meanwhile, the expression analysis for the selected 40 GST genes showed that they exhibited tissue-specific expression patterns and their expression were induced or repressed by salt stress. Those findings shed lights on the function and evolution of the GST gene family in Gossypium species.

Calcium-dependent protein kinases (CDPKs) are one of the largest protein kinases in plants and participate in different physiological processes through regulating downstream components of calcium signaling pathways. In this study, 41 CDPK genes, from GrCPK1 to GrCPK41, were identified in the genome of the diploid cotton, Gossypium raimondii. The phylogenetic analysis indicated that all these genes were divided into four subgroups and members within the same subgroup shared conserved exon-intron structures. The expansion of GrCPKs family in G. raimondii was due to the segmental duplication events, and the analysis of Ka/Ks ratios implied that the duplicated GrCPKs had mainly undergone strong purifying selection pressure with limited functional divergence. The cold-responsive elements in promoter regions were detected in the majority of GrCPKs. The expression analysis of 11 selected genes showed that GrCPKs exhibited tissue-specific expression patterns and the expression of GrCPKs were induced or repressed by cold treatment. These observations would lay an important foundation for functional and evolutionary analysis of CDPK gene family in Gossypium species.

Soil salinization poses a threat to the sustainability of agricultural production and has become a global issue. Cotton is an important cash crop and plays an important role in economic development. Salt stress has been harming the yield and quality of many crops, including cotton, for many years. In recent years, soil salinization has been increasing. It is crucial to study the mechanism of cotton salt tolerance and explore diversified materials and methods to alleviate the salt stress of cotton for the development of the cotton industry. Nanoparticles (NPs) are an effective means to alleviate salt stress. In this study, zinc oxide NPs (ZnO NPs) were sprayed on cotton leaves with the aim of investigating the intrinsic mechanism of NPs to alleviate salt stress in cotton. The results show that the foliar spraying of ZnO NPs significantly alleviated the negative effects of salt stress on hydroponic cotton seedlings, including the improvement of above-ground and root dry and fresh weight, leaf area, seedling height, and stem diameter. In addition, ZnO NPs can significantly improve the salt-induced oxidative stress by reducing the levels of MDA, H₂O₂, and O₂⁻ and increasing the activities of major antioxidant enzymes, such as superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT). Furthermore, RNA-seq showed that the foliar spraying of ZnO NPs could induce the expressions of CNGC, NHX2, AHA3, HAK17, and other genes, and reduce the expression of SKOR, combined with the CBL-CIPK pathway, which alleviated the toxic effect of excessive Na⁺ and reduced the loss of excessive K⁺ so that the Na⁺/K⁺ ratio was stabilized. In summary, our results indicate that the foliar application of ZnO NPs can alleviate high salt stress in cotton by adjusting the Na⁺/K⁺ ratio and regulating antioxidative ability. This provides a new strategy for alleviating the salt stress of cotton and other crops, which is conducive to the development of agriculture.

To produce unsaturated fatty acids, membrane-bound fatty acid desaturases (FADs) can be exploited to introduce double bonds into the acyl chains of fatty acids. In this study, 19 membrane-bound FAD genes were identified in Gossypium raimondii through database searches and were classified into four different subfamilies based on phylogenetic analysis. All 19 membrane-bound FAD proteins shared three highly conserved histidine boxes, except for GrFAD2.1, which lost the third histidine box in the C-terminal region. In the G. raimondii genome, tandem duplication might have led to the increasing size of the FAD2 cluster in the Omega Desaturase subfamily, whereas segmental duplication appeared to be the dominant mechanism for the expansion of the Sphingolipid and Front-end Desaturase subfamilies. Gene expression analysis showed that seven membrane-bound FAD genes were significantly up-regulated and that five genes were greatly suppressed in G. raimondii leaves exposed to low temperature conditions.

Main conclusion: A comprehensive account of theLBDgene family ofGossypiumwas provided in this work. Expression analysis and functional characterization revealed thatLBDgenes might play different roles inG. hirsutumandG. barbadense. The Lateral Organ Boundaries Domain (LBD) proteins comprise a plant-specific transcription factor family, which plays crucial roles in physiological processes of plant growth, development, and stress tolerance. In the present work, a systematical analysis ofLBDgene family from two allotetraploid cotton species,G. hirsutumandG. barbadense, together with their genomic donor species,G. arboreumandG. raimondii, was conducted. There were 131, 128, 62, and 68LBDsidentified inG. hirsutum,G. barbadense,G. arboreumandG. raimondii, respectively. The LBD proteins could be classified into two main classes, class I and class II, based on the structure of their lateral organ boundaries domain and traits of phylogenetic tree, and class I was further divided into five subgroups. The gene structure and motif composition analyses conducted in bothG. hirsutumandG. barbadenserevealed thatLBDgenes kept relatively conserved within the subfamilies. Synteny analysis suggested that segmental duplication acted as an important mechanism in expansion of the cottonLBDgene family. Cis-element analysis predicated the possible functions ofLBDgenes. Public RNA-seq data were investigated to analyze the expression patterns of cottonLBDgenes in various tissues as well as gene expression under abiotic stress treatments. Furthermore, RT-qPCR results found thatGhLBDshad various expression regulation under MeJA treatments. Expression analysis indicated the differential functions of cottonLBDgenes in response to abiotic stress and hormones.