Search

Cotton seeds, as the main by-product of cotton, are not only an important raw material for edible oil and feed but also a source of biofuel. The quality of cotton seeds directly affects cotton planting and is closely related to the yield and fiber quality. However, the molecular mechanism governing cotton seed size remains largely unexplored. This study investigates the regulatory mechanisms of cotton seed size by focusing on two cotton genotypes, N10 and N12, which exhibit notable phenotypic variations across multiple environments. Developing seeds were sampled at various stages (5, 20, 30, and 35 DPA) and subjected to RNA-seq. Temporal pattern clustering and WGCNA on differentially expressed genes identified 413 candidate genes, including these related to sugar metabolism that were significantly enriched in transcriptional regulation. A genetic transformation experiment indicated that the overexpression of the GhUXS5 gene encoding UDP-glucuronate decarboxylase 5 significantly increased seed size, suggesting an important role of GhUXS5 in regulating cotton seed size. This discovery provides crucial insights into the molecular mechanisms controlling cotton seed size, helping to unravel the complex regulatory network and offering new strategies and targets for cotton breeding to enhance the economic value of cotton seeds and overall cotton yield.

Drought stress significantly affects the growth, development, and yield of cotton, triggering the response of multiple genes. Among them, ascorbate peroxidase (APX) is one of the important antioxidant enzymes in the metabolism of reactive oxygen species in plants, and APX enhances the ability of plants to resist oxidation, thus increasing plant stress tolerance. Therefore, enhancing the activity of APX in cells is crucial to improving plant stress resistance. Previous studies have isolated differentially expressed proteins under drought stress (GhAPX7) in drought-resistant (KK1543) and drought-sensitive (XLZ26) plants. Thus, this study analyzed the expression patterns of GhAPX7 in different cotton tissues to verify the drought resistance function of GhAPX7 and explore its regulatory pathways. GhAPX7 had the highest expression in cotton leaves, which significantly increased under drought stress, suggesting that GhAPX7 is essential for improving antioxidant capacity and enzyme activities in cotton. GhAPX7 silencing indirectly affects pronounced leaf yellowing and wilting in drought-resistant and drought-sensitive plants under drought stress. Malondialdehyde (MDA) content was significantly increased and chlorophyll and proline content and APX enzyme activity were generally decreased in silenced plants compared to the control. This result indicates that GhAPX7 may improve drought resistance by influencing the contents of MDA, chlorophyll, proline, and APX enzyme activity through increased expression levels. Transcriptome analysis revealed that the drought-related differentially expressed genes between the control and treated groups enriched plant hormone signal transduction, MAPK signaling, and plant–pathogen interaction pathways. Therefore, the decreased expression of GhAPX7 significantly affects the expression levels of genes in these three pathways, reducing drought resistance in plants. This study provides insights into the molecular mechanisms of GhAPX7 and its role in drought resistance and lays a foundation for further research on the molecular mechanisms of response to drought stress in cotton.

DUSPs, a diverse group of protein phosphatases, play a pivotal role in orchestrating cellular growth and development through intricate signaling pathways. Notably, they actively participate in the MAPK pathway, which governs crucial aspects of plant physiology, including growth regulation, disease resistance, pest resistance, and stress response. DUSP is a key enzyme, and it is the enzyme that limits the rate of cell metabolism. At present, complete understanding of the DUSP gene family in cotton and its specific roles in resistance to Verticillium wilt (VW) remains elusive. To address this knowledge gap, we conducted a comprehensive identification and analysis of four key cotton species: Gossypium arboreum, Gossypium barbadense, Gossypium hirsutum, and Gossypium raimondii. The results revealed the identification of a total of 120 DUSP genes in the four cotton varieties, which were categorized into six subgroups and randomly distributed at both ends of 26 chromosomes, predominantly localized within the nucleus. Our analysis demonstrated that closely related DUSP genes exhibited similarities in terms of the conserved motif composition and gene structure. A promoter analysis performed on the GhDUSP gene promoter revealed the presence of several cis-acting elements, which are associated with abiotic and biotic stress responses, as well as hormone signaling. A tissue expression pattern analysis demonstrated significant variations in GhDUSP gene expression under different stress conditions, with roots exhibiting the highest levels, followed by stems and leaves. In terms of tissue-specific detection, petals, leaves, stems, stamens, and receptacles exhibited higher expression levels of the GhDUSP gene. The gene expression analysis results for GhDUSPs under stress suggest that DUSP genes may have a crucial role in the cotton response to stress in cotton. Through Virus-Induced Gene Silencing (VIGS) experiments, the silencing of the target gene significantly reduced the resistance efficiency of disease-resistant varieties against Verticillium wilt (VW). Consequently, we conclude that GH_A11G3500-mediated bispecific phosphorylated genes may serve as key regulators in the resistance of G. hirsutum to Verticillium wilt (VW). This study presents a comprehensive structure designed to provide an in-depth understanding of the potential biological functions of cotton, providing a strong foundation for further research into molecular breeding and resistance to plant pathogens.

Cotton is a critical crop with massive economic implications worldwide. Verticillium wilt is a soil-borne ailment caused by Verticillium dahliae, which harms the growth and development of cotton. Therefore, investigating the genes associated with resistance to verticillium wilt is of particular significance. In this study, we identified the GhIQD1 gene through transcriptome analysis and experimentally characterized the role of the GhIQD1 gene in cotton against V. dahliae. The findings indicated that GhIQD1 acts as a calmodulin-binding protein. The expression of GhIQD1 was the highest in stems, and the expression level increased significantly following infection with V. dahliae. The expression in resistant cotton varieties was higher than in susceptible cotton varieties. Through overexpression of the GhIQD1 gene in tobacco, these transgenic plants exhibited improved resistance to V. dahliae. In contrast, by silencing the GhIQD1 gene in cotton through VIGS, the resistance to V. dahliae was reduced. Following inoculation, the leaves yellowed, and the disease index was higher. Transcriptome analysis of transgenic tobacco 72 h after inoculation indicated that overexpression of GhIQD1 increased the enrichment of the calmodulin pathway and stimulated the production of plant hormones alongside secondary metabolites. Consequently, we investigated the relationship between the GhIQD1 gene and plant disease-resistant hormones SA, JA, and ABA. In summary, this study uncovered the mechanism by which GhIQD1 conferred resistance to V. dahliae in cotton through positive regulation of JA and ABA, providing crucial information for further research on the adaptation of plants to pathogen invasion.

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

Fiber quality traits, especially fiber strength, length, and micronaire (FS, FL, and FM), have been recognized as critical fiber attributes in the textile industry, while the lint percentage (LP) was an important indicator to evaluate the cotton lint yield. So far, the genetic mechanism behind the formation of these traits is still unclear. Quantitative trait loci (QTL) identification and candidate gene validation provide an effective methodology to uncover the genetic and molecular basis of FL, FS, FM, and LP. A previous study identified three important QTL/QTL cluster loci, harboring at least one of the above traits on chromosomes A01, A07, and D12 via a recombinant inbred line (RIL) population derived from a cross of Lumianyan28 (L28) × Xinluzao24 (X24). A secondary segregating population (F2) was developed from a cross between L28 and an RIL, RIL40 (L28 × RIL40). Based on the population, genetic linkage maps of the previous QTL cluster intervals on A01 (6.70–10.15 Mb), A07 (85.48–93.43 Mb), and D12 (0.40–1.43 Mb) were constructed, which span 12.25, 15.90, and 5.56 cM, with 2, 14, and 4 simple sequence repeat (SSR) and insertion/deletion (Indel) markers, respectively. QTLs of FL, FS, FM, and LP on these three intervals were verified by composite interval mapping (CIM) using WinQTL Cartographer 2.5 software via phenotyping of F2 and its derived F2:3 populations. The results validated the previous primary QTL identification of FL, FS, FM, and LP. Analysis of the RNA-seq data of the developing fibers of L28 and RIL40 at 10, 20, and 30 days post anthesis (DPA) identified seven differentially expressed genes (DEGs) as potential candidate genes. qRT-PCR verified that five of them were consistent with the RNA-seq result. These genes may be involved in regulating fiber development, leading to the formation of FL, FS, FM, and LP. This study provides an experimental foundation for further exploration of these functional genes to dissect the genetic mechanism of cotton fiber development.

Sulfotransferases (SOTs) (EC 2.8.2.-) are sulfate regulatory proteins in a variety of organisms that have been previously shown to be involved in regulating a variety of physiological and biological processes, such as growth, development, adaptation to land, stomatal closure, drought tolerance, and response to pathogen infection. However, there is a lack of comprehensive identification and systematic analysis of SOT in cotton, especially in G. barbadense. In this study, we used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structure, expression patterns, evolutionary relationships, selection pressure and stress response of SOT gene family members in G. barbadense. In this study, a total of 241 SOT genes were identified in four cotton species, among which 74 SOT gene members were found in G. barbadense. According to the phylogenetic tree, 241 SOT protein sequences were divided into five distinct subfamilies. We also mapped the physical locations of these genes on chromosomes and visualized the structural information of SOT genes in G. barbadense. We also predicted the cis-acting elements of the SOT gene in G. barbadense, and we identified the repetitive types and collinearity analysis of SOT genes in four cotton species. We calculated the Ka/Ks ratio between homologous gene pairs to elucidate the selective pressure between SOT genes. Transcriptome data were used to explore the expression patterns of SOT genes, and then qRT-PCR was used to detect the expression patterns of GBSOT4, GBSOT17 and GBSOT33 under FOV stress. WGCNA (weighted gene co-expression network analysis) showed that GB_A01G0479 (GBSOT4) belonged to the MEblue module, which may regulate the resistance mechanism of G. barbadense to FOV through plant hormones, signal transduction and glutathione metabolism. In addition, we conducted a VIGS (virus-induced gene silencing) experiment on GBSOT4, and the results showed that after FOV inoculation, the plants with a silenced target gene had more serious leaf wilting, drying and cracking than the control group, and the disease index of the plants with the silenced target gene was significantly higher than that of the control group. This suggests that GBSOT4 may be involved in protecting the production of G. barbadense from FOV infection. Subsequent metabolomics analysis showed that some flavonoid metabolites, such as Eupatorin-5-methylether (3′-hydroxy-5,6,7,4′-tetramethoxyflavone, were accumulated in cotton plants in response to FOV infection.

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

The increasing water scarcity associated with environmental change brings significant negative impacts to the growth of cotton plants, whereby it is urgent to enhance plant tolerance to drought. Here, we overexpressed the com58276 gene isolated from the desert plant Caragana korshinskii in cotton plants. We obtained three OE plants and demonstrated that com58276 confers drought tolerance in cotton after subjecting transgenic seeds and plants to drought. RNA-seq revealed the mechanisms of the possible anti-stress response, and that the overexpression of com58276 does not affect growth and fiber content in OE cotton plants. The function of com58276 is conserved across species, improving the tolerance of cotton to salt and low temperature, and demonstrating its applicability to improve plant resistance to environmental change.

Protein palmitoylation is an ability of the frame of the cell marker protein is one of the most notable reversible changes after translation. However, studies on protein palmitoylation in cotton have not yet been performed. In our current research, the PAT gene family was systematically identified and bioinformatically analyzed in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, and 211 PAT genes were authenticated and classified into six subfamilies. Sixty-nine PAT genes were identified in upland cotton, mainly at the ends of its the 26 chromosomes of upland cotton. The majority of these genes are located in the nucleus of the plant. Gene structure analysis revealed that each member encodes a protein that which contains at least one DHHC structural domain. Cis-acting element analysis indicated that GhPATs genes are mainly involved in hormone production, light response and stress response. Gene expression pattern analysis indicated that most GhPATs genes were differentially expressed upon induction by pathogenic bacteria, drought, salt, hot and cold stresses, and some GhPATs could be induced by multiple abiotic stresses simultaneously. GhPATs genes showed different expression patterns in tissue-specific assays and were found to be preferentially expressed in roots, followed by expression in stems and leaves. Virus-induced gene silencing (VIGS) experiments showed that cotton was significantly less resistant to Verticillium dahliae when GhPAT27 was silenced. We conclude that the GhPAT27 gene, which mediates S-palmitoylation acetylation, may be involved in the regulation of upland cotton resistance to Verticillium wilt (VW). Overall, this work has provided a fundamental framework for understanding the latent capabilities of GhPATs and a solid foundation for molecular breeding and plant pathogen resistance in cotton.

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

Ethylene plays a pivotal role in plant stress resistance and 1-aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis. Upland cotton (Gossypium hirsutum L.) is the most important natural fiber crop, but the function of ACS in response to abiotic stress has rarely been reported in this plant. We identified 18 GaACS, 18 GrACS, and 35 GhACS genes in Gossypiumarboreum, Gossypium raimondii and Gossypiumhirsutum, respectively, that were classified as types I, II, III, or IV. Collinearity analysis showed that the GhACS genes were expanded from diploid cotton by the whole-genome-duplication. Multiple alignments showed that the C-terminals of the GhACS proteins were conserved, whereas the N-terminals of GhACS10 and GhACS12 were different from the N-terminals of AtACS10 and AtACS12, probably diverging during evolution. Most type II ACS genes were hardly expressed, whereas GhACS10/GhACS12 were expressed in many tissues and in response to abiotic stress; for example, they were highly and hardly expressed at the early stages of cold and heat exposure, respectively. The GhACS genes showed different expression profiles in response to cold, heat, drought, and salt stress by quantitative PCR analysis, which indicate the potential roles of them when encountering the various adverse conditions, and provide insights into GhACS functions in cotton’s adaptation to abiotic stress.

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant (‘Xinluzao36′) and heat-sensitive (‘Che61-72′) cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

Cotton (Gossypium hirsutum L.) is one of the most important cash crops worldwide. Cytoplasmic male sterility (CMS) is an excellent breeding system for exploitation of heterosis, which has great potential to increase crop yields. To understand the molecular mechanism of CMS in cotton, we compared transcriptome, cytomorphological, physiological and bioinformatics data between the CMS line C2P5A and its maintainer line C2P5B. By using high-throughput sequencing technology, 178,166 transcripts were assembled and 2013 differentially expression genes (DEGs) were identified at three different stages of C2P5A anther development. In this study, we identified DEGs associated with reactive oxygen species (ROS), peroxisomes, aldehyde dehydrogenases (ALDH), cytochrome oxidase subunit VI, and cytochrome P450, and DEGs associated with tapetum development, Jojoba acyl-CoA reductase-related male sterility protein, basic helix-loop-helix (bHLH) and MYB transcription factors. The abnormal expression of one of these genes may be responsible for the CMS C2P5A line. In gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, DEGs were mainly related to carbohydrate metabolism, amino acid metabolism, transport and catabolism, and signal transduction. Carbohydrate metabolism provides energy for anther development, starch and sucrose metabolism, fatty acid biosynthesis and metabolism and ascorbate and aldarate metabolism. These results showed that numerous genes and multiple complex metabolic pathways regulate cotton anther development. Weighted correlation network analysis (WGCNA) indicated that three modules, ‘turquoise,’ ‘blue,’ and ‘green,’ were specific for the CMS C2P5A line. The ‘turquoise’ and ‘blue’ modules were mainly related to carbohydrate metabolism, amino acid metabolism, energy metabolism, peroxisomes, pyruvate metabolism as well as fatty acid degradation. The ‘green’ module was mainly related to energy metabolism, carbon metabolism, translation, and lipid metabolism. RNA-sequencing and WGCNA polymerization modules were screened for key genes and pathways related to CMS in cotton. This study presents a new perspective for further research into the metabolic pathways of pollen abortion in the CMS C2P5A line and also provides a theoretical basis for its breeding and production.