Exploring cotton SFR2's conundrum in response to cold stress.

Cotton is an important agricultural crop to many regions across the globe but is sensitive to low-temperature exposure. The activity of the enzyme SENSITIVE TO FREEZING 2 (SFR2) improves cold tolerance of plants and produces trigalactosylsyldiacylglycerol (TGDG), but its role in cold sensitive plants, such as cotton remains unknown. Recently, it was reported that cotton SFR2 produced very little TGDG under normal and cold conditions. Here, we investigate cotton SFR2 activation and TGDG production. Using multiple approaches in the native system and transformation into Arabidopsis thaliana, as well as heterologous yeast expression, we provide evidence that cotton SFR2 activates differently than previously found among other plant species. We conclude with the hypothesis that SFR2 in cotton is not activated in a similar manner regarding acidification or freezing like Arabidopsis and that other regions of SFR2 protein are critical for activation of the enzyme than previously reported.

Positive roles of the Ca2+ sensors GbCML45 and GbCML50 in improving cotton Verticillium wilt resistance.

As a universal second messenger, cytosolic calcium (Ca2+) functions in multifaceted intracellular processes, including growth, development and responses to biotic/abiotic stresses in plant. The plant-specific Ca2+ sensors, calmodulin and calmodulin-like (CML) proteins, function as members of the second-messenger system to transfer Ca2+ signal into downstream responses. However, the functions of CMLs in the responses of cotton (Gossypium spp.) after Verticillium dahliae infection, which causes the serious vascular disease Verticillium wilt, remain elusive. Here, we discovered that the expression level of GbCML45 was promoted after V. dahliae infection in roots of cotton, suggesting its potential role in Verticillium wilt resistance. We found that knockdown of GbCML45 in cotton plants decreased resistance while overexpression of GbCML45 in Arabidopsis thaliana plants enhanced resistance to V. dahliae infection. Furthermore, there was physiological interaction between GbCML45 and its close homologue GbCML50 by using yeast two-hybrid and bimolecular fluorescence assays, and both proteins enhanced cotton resistance to V. dahliae infection in a Ca2+-dependent way in a knockdown study. Detailed investigations indicated that several defence-related pathways, including salicylic acid, ethylene, reactive oxygen species and nitric oxide signalling pathways, as well as accumulations of lignin and callose, are responsible for GbCML45- and GbCML50-modulated V. dahliae resistance in cotton. These results collectively indicated that GbCML45 and GbCML50 act as positive regulators to improve cotton Verticillium wilt resistance, providing potential targets for exploitation of improved Verticillium wilt-tolerant cotton cultivars by genetic engineering and molecular breeding.

Salivary gland transcriptomics of the cotton aphid Aphis gossypii and comparative analysis with other sap-sucking insects.

Aphids are sap-sucking insects responsible for crop losses and a severe threat to crop production. Proteins in the aphid saliva are integral in establishing an interaction between aphids and plants and are responsible for host plant adaptation. The cotton aphid, Aphis gossypii (Hemiptera: Aphididae) is a major pest of Gossypium hirsutum. Despite extensive studies of the salivary proteins of various aphid species, the components of A. gossypii salivary glands are unknown. In this study, we identified 123,008 transcripts from the salivary gland of A. gossypii. Among those, 2933 proteins have signal peptides with no transmembrane domain known to be secreted from the cell upon feeding. The transcriptome includes proteins with more comprehensive functions such as digestion, detoxification, regulating host defenses, regulation of salivary glands, and a large set of uncharacterized proteins. Comparative analysis of salivary proteins of different aphids and other insects with A. gossypii revealed that 183 and 88 orthologous clusters were common in the Aphididae and non-Aphididae groups, respectively. The structure prediction for highly expressed salivary proteins indicated that most possess an intrinsically disordered region. These results provide valuable reference data for exploring novel functions of salivary proteins in A. gossypii with their host interactions. The identified proteins may help develop a sustainable way to manage aphid pests.

Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development.

Sulfate transporter (SULTR) proteins are in charge of the transport and absorption on sulfate substances, and have been reported to play vital roles in the biological processes of plant growth and stress response. However, there were few reports of genome-wide identification and expression-pattern analysis of SULTRs in Hibiscus mutabilis. Gossypium genus is a ideal model for studying the allopolyploidy, therefore two diploid species (G. raimondii and G. arboreum) and two tetraploid species (G. hirsutum and G. barbadense) were chosen in this study to perform bioinformatic analyses, identifying 18, 18, 35, and 35 SULTR members, respectively. All the 106 cotton SULTR genes were utilized to construct the phylogenetic tree together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 8 Zea mays ones, which was divided into Group1-Group4. The clustering analyses of gene structures and 10 conserved motifs among the cotton SULTR genes showed the consistent evolutionary relationship with the phylogenetic tree, and the results of gene-duplication identification among the four representative Gossypium species indicated that genome-wide or segment duplication might make main contributions to the expansion of SULTR gene family in cotton. Having conducted the cis-regulatory element analysis in promoter region, we noticed that the existing salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) elements could have influences with expression levels of cotton SULTR genes. The expression patterns of GhSULTR genes were also investigated on the 7 different tissues or organs and the developing ovules and fibers, most of which were highly expressed in root, stem, sepal, receptacel, ovule at 10 DPA, and fiber at 20 and 25 DPA. In addition, more active regulatory were observed in GhSULTR genes responding to multiple abiotic stresses, and 12 highly expressed genes showed the similar expression patterns in the quantitative Real-time PCR experiments under cold, heat, salt, and drought treatments. These findings broaden our insight into the evolutionary relationships and expression patterns of the SULTR gene family in cotton, and provide the valuable information for further screening the vital candidate genes on trait improvement.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton.

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

Investigation of salt tolerance in cotton germplasm by analyzing agro-physiological traits and ERF genes expression.

The development of genotypes that can tolerate high levels of salt is crucial for the efficient use of salt-affected land and for enhancing crop productivity worldwide. Therefore, incorporating salinity tolerance is a critical trait that crops must possess. Salt resistance is a complex character, controlled by multiple genes both physiologically and genetically. To examine the genetic foundation of salt tolerance, we assessed 16 F1 hybrids and their eight parental lines under normal and salt stress (15 dS/m) conditions. Under salt stress conditions significant reduction was observed for plant height (PH), bolls/plant (NBP), boll weight (BW), seed cotton yield (SCY), lint% (LP), fiber length (FL), fiber strength (FS), potassium to sodium ratio (K+/Na+), potassium contents (K+), total soluble proteins (TSP), carotenoids (Car) and chlorophyll contents. Furthermore, the mean values for hydrogen peroxide (H2O2), sodium contents (Na+), catalase (CAT), superoxide dismutase (SOD), peroxidase (POD), and fiber fineness (FF) were increased under salt stress. Moderate to high heritability and genetic advancement was observed for NBP, BW, LP, SCY, K+/Na+, SOD, CAT, POD, Car, TSP, FL, and FS. Mean performance and multivariate analysis of 24 cotton genotypes based on various agro-physiological and biochemical parameters suggested that the genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, along with crosses FBS-FALCON × JSQ-White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 achieved the maximum values for K+/Na+, K+, TSP, POD, Chlb, CAT, Car, LP, FS, FL, PH, NBP, BW, and SCY under salt stress and declared as salt resistant genotypes. The above-mentioned genotypes also showed relatively higher expression levels of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 at 15 dS/m and proved the role of these ERF genes in salt tolerance in cotton. These findings suggest that these genotypes have the potential for the development of salt-tolerant cotton varieties with desirable fiber quality traits.

GhMYB52 Like: A Key Factor That Enhances Lint Yield by Negatively Regulating the Lignin Biosynthesis Pathway in Fibers of Upland Cotton (Gossypium hirsutum L.).

In the context of sustainable agriculture and biomaterial development, understanding and enhancing plant secondary cell wall formation are crucial for improving crop fiber quality and biomass conversion efficiency. This is especially critical for economically important crops like upland cotton (Gossypium hirsutum L.), for which fiber quality and its processing properties are essential. Through comprehensive genome-wide screening and analysis of expression patterns, we identified a particularly high expression of an R2R3 MYB transcription factor, GhMYB52 Like, in the development of the secondary cell wall in cotton fiber cells. Utilizing gene-editing technology to generate a loss-of-function mutant to clarify the role of GhMYB52 Like, we revealed that GhMYB52 Like does not directly contribute to cellulose synthesis in cotton fibers but instead represses a subset of lignin biosynthesis genes, establishing it as a lignin biosynthesis inhibitor. Concurrently, a substantial decrease in the lint index, a critical measure of cotton yield, was noted in parallel with an elevation in lignin levels. This study not only deepens our understanding of the molecular mechanisms underlying cotton fiber development but also offers new perspectives for the molecular improvement of other economically important crops and the enhancement of biomass energy utilization.

GhCLCc-1, a Chloride Channel Gene from Upland Cotton, Positively Regulates Salt Tolerance by Modulating the Accumulation of Chloride Ions.

The ionic toxicity induced by salinization has adverse effects on the growth and development of crops. However, researches on ionic toxicity and salt tolerance in plants have focused primarily on cations such as sodium ions (Na+), with very limited studies on chloride ions (Cl-). Here, we cloned the homologous genes of Arabidopsis thaliana AtCLCc, GhCLCc-1A/D, from upland cotton (Gossypium hirsutum), which were significantly induced by NaCl or KCl treatments. Subcellular localization showed that GhCLCc-1A/D were both localized to the tonoplast. Complementation of Arabidopsis atclcc mutant with GhCLCc-1 rescued its salt-sensitive phenotype. In addition, the silencing of the GhCLCc-1 gene led to an increased accumulation of Cl- in the roots, stems, and leaves of cotton seedlings under salt treatments, resulting in compromised salt tolerance. And ectopic expression of the GhCLCc-1 gene in Arabidopsis reduced the accumulation of Cl- in transgenic lines under salt treatments, thereby enhancing salt tolerance. These findings elucidate that GhCLCc-1 positively regulates salt tolerance by modulating Cl- accumulation and could be a potential target gene for improving salt tolerance in plants.

GhHDZ76, a cotton HD-Zip transcription factor, involved in regulating the initiation and early elongation of cotton fiber development in G. hirsutum.

In this study, the whole HD-Zip family members of G. hirsutum were identified, and GhHDZ76 was classified into the HD-Zip IV subgroup. GhHDZ76 was predominantly expressed in the 0-5 DPA of fiber development stage and localized in the nucleus. Overexpression of GhHDZ76 significantly increased the length and density of trichomes in Arabidopsis thaliana. The fiber length of GhHDZ76 knockout lines by CRISPR/Cas9 was significantly shorter than WT at the early elongation and mature stage, indicating that GhHDZ76 positively regulate the fiber elongation. Scanning electron microscopy showed that the number of ovule surface protrusion of 0 DPA of GhHDZ76 knockout lines was significantly lower than WT, suggesting that GhHDZ76 can also promote the initiation of fiber development. The transcript level of GhWRKY16, GhRDL1, GhEXPA1 and GhMYB25 genes related to fiber initiation and elongation in GhHDZ76 knockout lines were significantly decreased. Yeast two-hybrid and Luciferase complementation imaging (LCI) assays showed that GhHDZ76 can interact with GhWRKY16 directly. As a transcription factor, GhHDZ76 has transcriptional activation activity, which could bind to L1-box elements of the promoters of GhRDL1 and GhEXPA1. Double luciferase reporter assay showed that the GhWRKY16 could enhance the transcriptional activity of GhHDZ76 to pGhRDL1, but it did not promote the transcriptional activity of GhHDZ76 to pGhEXPA1. GhHDZ76 protein may also promote the transcriptional activity of GhWRKY16 to the downstream target gene GhMYB25. Our results provided a new gene resource for fiber development and a theoretical basis for the genetic improvement of cotton fiber quality.

Characterization and expression analysis of GLABRA3 (GL3) genes in cotton: insights into trichome development and hormonal regulation.

BACKGROUND: GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) genes encode a typical helix-loop-helix (bHLH) transcription factors that primarily regulate trichome branching and root hair development, DNA endoreduplication, trichoblast size, and stomatal formation. The functions of GL3 genes in cotton crop have been poorly characterized. In this study, we performed comprehensive genome-wide scans for GL3 and EGL3 homologs to enhance our comprehension of their potential roles in trichome and fiber development in cotton crop. METHODS AND RESULTS: Our findings paraded that Gossypium hirsutum and G. barbadense have 6 GL3s each, unevenly distributed on 4 chromosomes whereas, G. arboreum, and G. raimondii have 3 GL3s each, unevenly distributed on 2 chromosomes. Gh_A08G2088 and Gb_A09G2187, despite having the same bHLH domain as the other GL3 genes, were excluded due to remarkable short sequences and limited number of motifs, indicating a lack of potential functional activity. The phylogenetic analysis categorized remaining 16 GL3s into three subfamilies (Group I-III) closely related to A. thaliana. The 16 GL3s have complete bHLH domain, encompassing 590-631 amino acids, with molecular weights (MWs) ranging from 65.92 to 71.36 kDa. Within each subfamily GL3s depicted shared similar gene structures and motifs, indicating conserved characteristics within respective groups. Promoter region analysis revealed 27 cis-acting elements, these elements were responsive to salicylic acid, abscisic acid (ABA), methyl jasmonate (MeJA), and gibberellin. The expression of GL3 genes was analyzed across 12 tissues in both G. barbadense and G. hirsutum using the publicly available RNA-seq data. Among GL3s, Gb_D11G0219, Gb_D11G0214, and Gb_D08G2182, were identified as relatively highly expressed across different tissues, consequently selected for hormone treatment and expression validation in G. barbadense. RT-qPCR results demonstrated significant alterations in the expression levels of Gb_D11G0219 and Gb_D11G0214 following MeJA, GA, and ABA treatment. Subcellular localization prediction revealed that most GL3 proteins were predominantly expressed in the nucleus, while a few were localized in the cytoplasm and chloroplasts. CONCLUSIONS: In summary, this study lays the foundation for subsequent functional validation of GL3 genes by identifying hormonal regulation patterns and probable sites of action in cotton trichome formation and fiber development. The results stipulate a rationale to elucidate the roles and regulatory mechanisms of GL3 genes in the intricate process of cotton fibre and trichome development.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton.

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

A Genome-Wide Analysis of the CEP Gene Family in Cotton and a Functional Study of GhCEP46-D05 in Plant Development.

C-TERMINALLY ENCODED PEPTIDEs (CEPs) are a class of peptide hormones that have been shown in previous studies to play an important role in regulating the development and response to abiotic stress in model plants. However, their role in cotton is not well understood. In this study, we identified 54, 59, 34, and 35 CEP genes from Gossypium hirsutum (2n = 4x = 52, AD1), G. barbadense (AD2), G. arboreum (2n = 2X = 26, A2), and G. raimondii (2n = 2X = 26, D5), respectively. Sequence alignment and phylogenetic analyses indicate that cotton CEP proteins can be categorized into two subgroups based on the differentiation of their CEP domain. Chromosomal distribution and collinearity analyses show that most of the cotton CEP genes are situated in gene clusters, suggesting that segmental duplication may be a critical factor in CEP gene expansion. Expression pattern analyses showed that cotton CEP genes are widely expressed throughout the plant, with some genes exhibiting specific expression patterns. Ectopic expression of GhCEP46-D05 in Arabidopsis led to a significant reduction in both root length and seed size, resulting in a dwarf phenotype. Similarly, overexpression of GhCEP46-D05 in cotton resulted in reduced internode length and plant height. These findings provide a foundation for further investigation into the function of cotton CEP genes and their potential role in cotton breeding.

Genome-Wide Identification of Calmodulin-Binding Protein 60 Gene Family and the Function of GhCBP60B in Cotton Growth and Development and Abiotic Stress Response.

The calmodulin-binding protein 60 (CBP60) family is a gene family unique to plants, and its members play a crucial role in plant defense responses to pathogens and growth and development. Considering that cotton is the primary source of natural cotton textile fiber, the functional study of its CBP60 gene family members is critical. In this research, we successfully identified 162 CBP60 members from the genomes of 21 species. Of these, 72 members were found in four cotton species, divided into four clades. To understand the function of GhCBP60B in cotton in depth, we conducted a detailed analysis of its sequence, structure, cis-acting elements, and expression patterns. Research results show that GhCBP60B is located in the nucleus and plays a crucial role in cotton growth and development and response to salt and drought stress. After using VIGS (virus-induced gene silencing) technology to conduct gene silencing experiments, we found that the plants silenced by GhCBP60B showed dwarf plants and shortened stem nodes, and the expression of related immune genes also changed. In further abiotic stress treatment experiments, we found that GhCBP60B-silenced plants were more sensitive to drought and salt stress, and their POD (peroxidase) activity was also significantly reduced. These results imply the vital role of GhCBP60B in cotton, especially in regulating plant responses to drought and salt stress. This study systematically analyzed CBP60 gene family members through bioinformatics methods and explored in depth the biological function of GhCBP60B in cotton. These research results lay a solid foundation for the future use of the GhCBP60B gene to improve cotton plant type and its drought and salt resistance.

Cotton sphingosine kinase GhLCBK1 participates in fiber cell elongation by affecting sphingosine-1-phophate and auxin synthesis.

Sphingolipids serve as essential components of biomembrane and possess significant bioactive properties. Sphingosine-1-phophate (S1P) plays a key role in plant resistance to stress, but its specific impact on plant growth and development remains to be fully elucidated. Cotton fiber cells are an ideal material for investigating the growth and maturation of plant cells. In this study, we examined the content and composition of sphingosine (Sph) and S1P throughout the progression of fiber cell development. The content of S1P elevated gradually during fiber elongation but declined during the transition stage. Exogenous application of S1P promoted fiber elongation while using of FTY720 (an antagonist of S1P), and DMS (an inhibitor of LCBK) hindered fiber elongation. Cotton Long Chain Base Kinase 1 (GhLCBK1) was notably expressed during the fiber elongation stage, containing all conserved domains of LCBK protein and localized in the endoplasmic reticulum. Overexpression GhLCBK1 increased the S1P content and promoted fiber elongation while retarded secondary cell wall (SCW) deposition. Conversely, downregulation of GhLCBK1 reduced the S1P levels, and suppressed fiber elongation, and accelerated SCW deposition. Transcriptome analysis revealed that upregulating GhLCBK1 or applying S1P induced the expression of GhEXPANSIN and auxin related genes. Furthermore, the levels of IAA were elevated and reduced in the fibers when up-regulating or down-regulating GhLCBK1, respectively. Our investigation demonstrated that GhLCBK1 and its product S1P facilitated the elongation of fiber cells by affecting auxin biosynthesis. This study contributes novel insights into the intricate regulatory pathways involved in fiber cell elongation, identifying GhLCBK1 as a potential target gene and laying the groundwork for enhancing fiber quality via genetic manipulation.

Systematic characterization of Gossypium GLN family genes reveals a potential function of GhGLN1.1a regulates nitrogen use efficiency in cotton.

The enzyme glutamine synthetase (GLN) is mainly responsible for the assimilation and reassimilation of nitrogen (N) in higher plants. Although the GLN gene has been identified in various plants, there is little information about the GLN family in cotton (Gossypium spp.). To elucidate the roles of GLN genes in cotton, we systematically investigated and characterized the GLN gene family across four cotton species (G. raimondii, G. arboreum, G. hirsutum, and G. barbadense). Our analysis encompassed analysis of members, gene structure, cis-element, intragenomic duplication, and exploration of collinear relationships. Gene duplication analysis indicated that segmental duplication was the primary driving force for the expansion of the GhGLN gene family. Transcriptomic and quantitative real-time reverse-transcription PCR (qRT-PCR) analyses indicated that the GhGLN1.1a gene is responsive to N induction treatment and several abiotic stresses. The results of virus-induced gene silencing revealed that the accumulation and N use efficiency (NUE) of cotton were affected by the inactivation of GhGLN1.1a. This study comprehensively analyzed the GhGLN genes in Gossypium spp., and provides a new perspective on the functional roles of GhGLN1.1a in regulating NUE in cotton.

The effect of farming techniques on degradation of DDT in historical cotton farms.

DDT was used in the mid 20th century for crop and livestock production. After use, DDT and its degradates DDE and DDD (collectively DDX) remain in the environment for decades. A few studies have reported that the rate of degradation of DDT into its metabolites is affected by various farming techniques like tillage, irrigation, and use of fertilizers. However, most of these studies did not evaluate active farms, and none of them focused on the Southeast US or historical cotton farms. Therefore, in this study, we aimed to determine if different farming techniques affect the decomposition of DDT in Walton County, Georgia, where farms historically grew cotton. Five Walton County farms were sampled for soil, and churches were sampled as control sites. The extensive land history of the farms was recorded, and the soil levels of p,p'-DDT, p,p'-DDE, p,p'-DDD, o,p'-DDT, and o,p'-DDE were measured using gas chromatography-tandem mass spectrometry. All farm sites had detectable levels of p,p'-DDT, p,p'-DDE, and p,p'-DDD, while few sites had detectable levels of o,p'-DDT and o,p'-DDE. Tillage was found to speed up p,p'-DDE degradation, but there was no effect on p,p'-DDT degradation. Plowing was associated with an increase in decomposition of p,p'-DDT, but p,p'-DDE and p,p'-DDD were not significantly increased. The largest difference in the degradation of DDT was based on the fertilizer type. Natural fertilizer sped up degradation of p,p'-DDT and p,p'-DDE; synthetic fertilizer increased p,p'-DDE degradation, but not p,p'-DDT degradation.

Identification of Crucial Modules and Genes Associated with Bt Gene Expression in Cotton.

The expression of Bacillus thuringiensis (Bt) toxins in transgenic cotton confers resistance to insect pests. However, it has been demonstrated that its effectiveness varies among cotton cultivars and different tissues. In this study, we evaluated the expression of Bt protein in 28 cotton cultivars and selected 7 cultivars that differed in Bt protein expression for transcriptome analysis. Based on their Bt protein expression levels, the selected cultivars were categorized into three groups: H (high Bt protein expression), M (moderate expression), and L (low expression). In total, 342, 318, and 965 differentially expressed genes were detected in the H vs. L, M vs. L, and H vs. M comparison groups, respectively. And three modules significantly associated with Bt protein expression were identified by weighted gene co-expression network analysis. Three hub genes were selected to verify their relationships with Bt protein expression using virus-induced gene silencing (VIGS). Silencing GhM_D11G1176, encoding an MYC transcription factor, was confirmed to significantly decrease the expression of Bt protein. The present findings contribute to an improved understanding of the mechanisms that influence Bt protein expression in transgenic cotton.

Effect of Arbuscular Mycorrhizal Fungi (AMF) on photosynthetic characteristics of cotton seedlings under saline-alkali stress.

The study aimed to find the best Arbuscular Mycorrhizal Fungi (AMF) strain for cotton growth in Xinjiang's salinity and alkali conditions. Cotton (Xinluzao 45) was treated with Funneliformis mosseae (GM), Rhizophagus irregularis (GI), and Claroideoglomus etunicatum (GE) as treatments, while untreated cotton served as the control (CK). Salinity stress was applied post-3-leaf stage in cotton. The study analyzed cotton's reactions to diverse saline-alkali stresses, focusing on nutrient processes and metabolism. By analyzing the growth and photosynthetic characteristics of plants inoculated with Funneliformis mosseae to evaluate its salt tolerance. Saline-alkali stress reduced chlorophyll and hindered photosynthesis, hampering cotton growth. However, AMF inoculation mitigated these effects, enhancing photosynthetic rates, CO2 concentration, transpiration, energy use efficiency, and overall cotton growth under similar stress levels. GM and GE treatments yielded similar positive effects. AMF inoculation enhanced cotton plant height and biomass. In GM treatment, cotton exhibited notably higher root length than other treatments, showing superior growth under various conditions. In summary, GM-treated cotton had the highest infection rate, followed by GE-treated cotton, with GI-treated cotton having the lowest rate (GM averaging 0.95). Cotton inoculated with Funneliformis mosseae, Rhizophagus irregularis, and Claroideoglomus etunicatum juvenile showed enhanced chlorophyll and photosynthetic levels, reducing salinity effects. Funneliformis mosseae had the most significant positive impact.

Comprehensive analysis of MAPK gene family in upland cotton (Gossypium hirsutum) and functional characterization of GhMPK31 in regulating defense response to insect infestation.

KEY MESSAGE: The transcriptomic, phenotypic and metabolomic analysis of transgenic plants overexpressing GhMPK31 in upland cotton revealed the regulation of H2O2 burst and the synthesis of defensive metabolites by GhMPK31. Mitogen-activated protein kinases (MAPKs) are a crucial class of protein kinases, which play an essential role in various biological processes in plants. Upland cotton (G. hirsutum) is the most widely cultivated cotton species with high economic value. To gain a better understanding of the role of the MAPK gene family, we conducted a comprehensive analysis of the MAPK gene family in cotton. In this study, a total of 55 GhMPK genes were identified from the whole genome of G. hirsutum. Through an investigation of the expression patterns under diverse stress conditions, we discovered that the majority of GhMPK family members demonstrated robust responses to abiotic stress, pathogen stress and pest stress. Furthermore, the overexpression of GhMPK31 in cotton leaves led to a hypersensitive response (HR)-like cell death phenotype and impaired the defense capability of cotton against herbivorous insects. Transcriptome and metabolomics data analysis showed that overexpression of GhMPK31 enhanced the expression of H2O2-related genes and reduced the accumulation of defensive related metabolites. The direct evidence of GhMPK31 interacting with GhRBOHB (H2O2-generating protein) were found by Y2H, BiFC, and LCI. Therefore, we propose that the increase of H2O2 content caused by overexpression of GhMPK31 resulted in HR-like cell death in cotton leaves while reducing the accumulation of defensive metabolites, ultimately leading to a decrease in the defense ability of cotton against herbivorous insects. This study provides valuable insights into the function of MAPK genes in plant resistance to herbivorous insects.

Gh4CL20/20A involved in flavonoid biosynthesis is essential for male fertility in cotton (Gossypium hirsutum L.).

Flavonoids have been shown to play an essential role in plant growth and fertility. 4-Coumarate CoA ligase (4CL) is one of the indispensable enzymes involved in the biosynthesis of flavonoids. However, the role of 4CL and flavonoids in impact on cotton fertility is still unknown. In this study, on the basis of identification of an additional Gh4CL gene, Gh4CL20A, by using an updated G. hirsutum genome, we found that Gh4CL20A and its homologous Gh4CL20 were preferentially expressed in petals and stamens. The petals of the loss-of-function Gh4CL20/Gh4CL20A mutant generated by CRISPR/Cas9 gene editing remained white until wilting. Notably, the mutant showed indehiscent anthers, reduced number of pollen grains and pollen viability, leading to male sterility. Histological analysis revealed that abnormal degradation of anther tapetum at the tetrad stage and abnormal pollen grain development at the mature stage caused male sterility of the gene editing mutant. Analysis of the anther transcriptome identified a total of 10574 and 11962 genes up- and down-regulated in the mutant, respectively, compared to the wild-type. GO, KEGG, and WGCNA analyses linked the abnormality of the mutant anthers to the defective flavonoid biosynthetic pathway, leading to decreased activity of 4CL and chalcone isomerase (CHI) and reduced accumulation of flavonoids in the mutant. These results imply a role of Gh4CL20/Gh4CL20A in assuring proper development of cotton anthers by regulating flavonoid metabolism. This study elucidates a molecular mechanism underlying cotton anther development and provides candidate genes for creating cotton male sterile germplasm that has the potential to be used in production of hybrid seeds.