E3 ligases are key enzymes required for protein degradation. Here we identified a C3H2C3 RING domain-containing E3 ubiquitin ligase gene named GhATL68b. It is preferentially and highly expressed in developing cotton fiber cells, and is more conserved in plants than in animals or in archaea. All four orthologs copies of this gene in various diploid cottons and eight in the allotetraploid G. hirsutum were found to originate from a single common ancestor that can be traced back to C. reinhardii at about 992 million years ago (MYA). Structural variations (SVs) occurred in the promoter regions of G. hirsutum, G. herbaceum, G. arboreum and G. raimondii correlated with significantly different methylation patterns. Homozygous CRISPR-Cas9 knock-out cotton lines produced significant poor fiber quality in terms of upper half mean length, elongation at break, uniformity and mature fiber weight. GhATL68b was shown to modulate the homeostasis of 2,4-dienoyl-CoA reductase (GhDECR), a rate-limiting enzyme for β-oxidation of polyunsaturated fatty acids (PUFAs) via the ubiquitin proteasome pathway through in vitro ubiquitination and cell-free protein degradation assays. Fiber cells harvested from these knockout mutants contained significantly lower levels of PUFAs important for glycerophospholipids production and also for plasma membrane fluidity regulations. Finally, the mutant fiber-growth defective phenotype can be fully compensated by adding linolenic acid (C18:3), the most abundant type of PUFA externally in ovule culture media. To our knowledge, this is the first experimentally characterized C3H2C3 type E3 ubiquitin ligase that is involved in regulating fiber cell elongation, and it may thus provide us with a new genetic target for improved cotton lint production.

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

Mitogen-activated protein kinase (MAPK) cascade is the center of plant signal transduction system that amplify immune signals into cellular responses by phosphorylating diverse substrates. The MAPK cascade consisting of MAPK kinase kinases (MAPKKKs), MAPK kinases (MAPKKs), and MAPKs is well characterized in plants, in which Raf-like kinases are generally regarded as MAPKKKs. However, it is rarely reported that Raf-like MAPKKKs function as middle regulators to link MAPK and its downstream transcription factors in plant immunity. Verticillium wilt, caused by the soil-borne vascular fungus Verticillium dahliae, is a serious disease in many plants, including cotton. The previous studies showed that GhMPK9 (a MAPK) is involved in the response to Verticillium wilt. Here, the Raf-like kinase GhRAF39_1 is reported as helper regulates the phosphorylation of WRKY transcription factor GhWRKY40a by GhMPK9. The phosphorylated GhWRKY40a can further activate the transcription of GhERF1b to up-regulate defense-related genes while inhibit the transcription of GhABF2 to regulate the stomatal opening, thus improving the resistance to Verticillium wilt in cotton. This study reveals a new signaling module of GhMPK9-GhRAF39_1-GhWRKY40a to regulate GhERF1b- and GhABF2-mediated defense responses, which triggers plant defense against Verticillium wilt.

Partial root-zone drying irrigation (PRD) can improve water-use efficiency (WUE) without reductions in photosynthesis; however, the mechanism by which this is attained is unclear. To amend that, PRD conditions were simulated by polyethylene glycol 6000 in a root-splitting system and the effects of PRD on cotton growth were studied. Results showed that PRD decreased stomatal conductance (gs) but increased mesophyll conductance (gm). Due to the contrasting effects on gs and gm, net photosynthetic rate (AN) remained unaffected, while the enhanced gm/gs ratio facilitated a larger intrinsic WUE. Further analyses indicated that PRD-induced reduction of gs was related to decreased stomatal size and stomatal pore area in adaxial and abaxial surface which was ascribed to lower pore length and width. PRD-induced variation of gm was ascribed to the reduced liquid-phase resistance, due to increases in chloroplast area facing to intercellular airspaces and the ratio of chloroplast surface area to total mesophyll cell area exposed to intercellular airspaces, as well as to decreases in the distance between cell wall and chloroplast, and between adjacent chloroplasts. The above results demonstrate that PRD, through alterations to stomatal and mesophyll structures, decoupled gs and gm responses, which ultimately increased intrinsic WUE and maintained AN.

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Leaf shape is considered to be one of the most significant agronomic traits in crop breeding. However, the molecular basis underlying leaf morphogenesis in cotton is still largely unknown. In this study, through genetic mapping and molecular investigation using a natural cotton mutant cu with leaves curling upward, the causal gene GHCU is successfully identified as the key regulator of leaf flattening. Knockout of GHCU or its homolog in cotton and tobacco using CRISPR results in abnormal leaf shape. It is further discovered that GHCU facilitates the transport of the HD protein KNOTTED1-like (KNGH1) from the adaxial to the abaxial domain. Loss of GHCU function restricts KNGH1 to the adaxial epidermal region, leading to lower auxin response levels in the adaxial boundary compared to the abaxial. This spatial asymmetry in auxin distribution produces the upward-curled leaf phenotype of the cu mutant. By analysis of single-cell RNA sequencing and spatiotemporal transcriptomic data, auxin biosynthesis genes are confirmed to be expressed asymmetrically in the adaxial-abaxial epidermal cells. Overall, these findings suggest that GHCU plays a crucial role in the regulation of leaf flattening through facilitating cell-to-cell trafficking of KNGH1 and hence influencing the auxin response level.

DEAD-box RNA helicases, a prominent subfamily within the RNA helicase superfamily 2 (SF2), play crucial roles in the growth, development, and abiotic stress responses of plants. This study identifies 146 DEAD-box RNA helicase genes (GhDEADs) and categorizes them into four Clades (Clade A-D) through phylogenetic analysis. Promoter analysis reveals cis-acting elements linked to plant responses to light, methyl jasmonate (MeJA), abscisic acid (ABA), low temperature, and drought. RNA-seq data demonstrate that Clade C GhDEADs exhibit elevated and ubiquitous expression across different tissues, validating their connection to leaf development through real-time quantitative polymerase chain reaction (RT-qPCR) analysis. Notably, over half of GhDEADs display up-regulation in the leaves of virus-induced gene silencing (VIGS) plants of GhVIR-A/D (members of m6A methyltransferase complex, which regulate leaf morphogenesis). In conclusion, this study offers a comprehensive insight into GhDEADs, emphasizing their potential involvement in leaf development.

The reniform nematode, Rotylenchulus reniformis, is a major yield-limiting pest of upland cotton (Gossypium hirsutum) in the United States that has been steadily increasing in incidence in many states. Reniform nematode-resistant cotton cultivars have recently become commercially available for cotton producers; however, few field trials have evaluated their efficacy as a nematode management tool. The aim of this study was to evaluate reniform nematode population development, plant growth, and seed cotton yield of reniform nematode-resistant cotton cultivars in two nematode-infested fields in Louisiana. Replicated small-plot field trials were conducted in St. Joseph, LA (NERS field) and Winnsboro, LA (MRRS field) during the 2022 and 2023 growing seasons. In 2022, cultivars evaluated included: (1) DP 1646 B2XF (susceptible/tolerant), (2) DP 2141NR B3XF (resistant), (3) PHY 332 W3FE (resistant), (4) PHY 411 W3FE (resistant), and (5) PHY 443 W3FE (resistant). In 2023, an additional susceptible cotton cultivar, PHY 340 W3FE, was also included. All nematode-resistant cotton cultivars evaluated provided suppression of reniform nematode population development relative to that of the susceptible cotton cultivars, with suppression of nematode soil population densities at harvest ranging from 49 - 81% relative to DP 1646 B2XF. The resistant cultivar PHY 411 W3FE provided the most consistent suppression of reniform nematode population development, reducing reniform nematode soil population densities at harvest in both field locations and both trial years. In contrast, DP 2141NR B3XF only reduced soil population densities at harvest in the NERS field in 2023. Despite relatively consistent nematode suppression and improvements in plant vigor ratings and canopy coverage associated with the resistant cotton cultivars, a yield increase was only observed with PHY 332 W3FE and PHY 411 W3FE planted at the NERS field in 2023. Despite strong resistance to reniform nematode in the evaluated cotton cultivars, nematode soil population densities still increased during the growing season in plots planted with resistant cotton cultivars, emphasizing the need for additional management tactics to use alongside host resistance. This study indicates that new reniform nematode-resistant cotton cultivars show promising potential to reduce nematode population development during the growing season in Louisiana.

The cotton boll weevil (CBW, Anthonomus grandis) stands as one of the most significant threats to cotton crops (Gossypium hirsutum). Despite substantial efforts, the development of a commercially viable transgenic cotton event for effective open-field control of CBW has remained elusive. This study describes a detailed characterization of the insecticidal toxins Cry23Aa and Cry37Aa against CBW. Our findings reveal that CBW larvae fed on artificial diets supplemented exclusively with Cry23Aa decreased larval survival by roughly by 69%, while supplementation with Cry37Aa alone displayed no statistical difference compared to the control. However, the combined provision of both toxins in the artificial diet led to mortality rates approaching 100% among CBW larvae (LC50 equal to 0.26 PPM). Additionally, we engineered transgenic cotton plants by introducing cry23Aa and cry37Aa genes under control of the flower bud-specific pGhFS4 and pGhFS1 promoters, respectively. Seven transgenic cotton events expressing high levels of Cry23Aa and Cry37Aa toxins in flower buds were selected for greenhouse bioassays, and the mortality rate of CBW larvae feeding on their T0 and T1 generations ranged from 75% to 100%. Our in silico analyses unveiled that Cry23Aa displays all the hallmark characteristics of β-pore-forming toxins (β-PFTs) that bind to sugar moieties in glycoproteins. Intriguingly, we also discovered a distinctive zinc-binding site within Cry23Aa, which appears to be involved in protein-protein interactions. Finally, we discuss the major structural features of Cry23Aa that likely play a role in the toxin\'s mechanism of action. In view of the low LC50 for CBW larvae and the significant accumulation of these toxins in the flower buds of both T0 and T1 plants, we anticipate that through successive generations of these transgenic lines, cotton plants engineered to overexpress cry23Aa and cry37Aa hold promise for effectively managing CBW infestations in cotton crops.

Naturally-colored brown cotton (NBC) fiber is an environmentally friendly raw source of fiber for textile applications. The fiber of some NBC cultivars exhibits flame-retardant properties, which can be used in textiles that require flame resistance. Proanthocyanidins or their derivatives are responsible for the brown pigment in NBC; however, how flame retardancy is related to pigmentation in NBC is poorly understood. To gain insight into brown pigment biosynthesis, we conducted comparative transcripts and metabolites profiling analysis of developing cotton fibers between the brown (MC-BL) and white (MC-WL) cotton near-isogenic lines (NILs), genetically different only in the Lc1 locus. In this study, mass spectrometry was used to detect metabolites in BL and WL developing fibers at 8, 12, 16, 20, 24, 36, and 40 days post anthesis (DPA) and mature fibers. Transcripts analysis was performed at two critical fiber developmental points, 8 DPA (fiber elongation) and 20 DPA (secondary cell wall deposition). We found 5836 (ESI MS positive mode) and 4541 (ESI MS negative mode) metabolites significantly different accumulated between BL and WL. Among them, 142 were known non-redundant metabolites, including organic acids, amino acids, and derivatives of the phenylpropanoid pathway. Transcript analysis determined 1691 (8 DPA) and 5073 (20 DPA) differentially expressed genes (DEGs) between BL and WL, with the majority of DEGs down-regulated at 20 DPA. Organic acids of the citric acid cycle were induced, while most of the detected amino acids were reduced in the MC-BL line. Both cis- and trans-stereoisomers of flavan-3-ols were detected in developing MC-WL and MC-BL fibers; however, the gallocatechin and catechin accumulated multiple times higher. Gas chromatography-mass spectrometry (GC-MS) analysis of fatty acids determined that palmitic acid long-chain alcohols were the main constituents of waxes of mature fibers. Energy-dispersive X-ray spectrometry (EDS) analysis of mature fibers revealed that potassium accumulated three times greater in MC-BL than in MC-WL mature fibers. This study provides novel insights into the biosynthesis of pigments and its association with flame retardancy in NBC fibers.