Systematical characterization of Rab7 gene family in Gossypium and potential functions of GhRab7B3-A gene in drought tolerance.

BACKGROUND: Cotton serves as a primary source of natural fibers crucial for the textile industry. However, environmental elements such as drought have posed challenges to cotton cultivation, resulting in adverse impacts on both production and fiber quality. Improving cotton's resilience to drought could mitigate yield losses and foster the expansion of cotton farming. Rab7 protein, widely present in organisms, controls the degradation and recycling of cargo, and has a potential role in biotic and abiotic tolerance. However, comprehensive exploration of the Rab7 gene family in Gossypium remains scarce. RESULTS: Herein, we identified a total of 10, 10, 20, and 20 Rab7 genes through genome-wide analysis in Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense, respectively. Collinearity analysis unveiled the pivotal role of whole genome or segmental duplication events in the expansion of GhRab7s. Study of gene architecture, conserved protein motifs, and domains suggested the conservation of structure and function throughout evolution. Exploration of cis-regulatory elements revealed the responsiveness of GhRab7 genes to abiotic stress, corroborated by transcriptome analysis under diverse environmental stresses. Notably, the greatly induced expression of GhRab7B3-A under drought treatment prompted us to investigate its function through virus-induced gene silencing (VIGS) assays. Silencing GhRab7B3-A led to exacerbated dehydration and wilting compared with the control. Additionally, inhibition of stomatal closure, antioxidant enzyme activities and expression patterns of genes responsive to abiotic stress were observed in GhRab7B3-A silenced plants. CONCLUSIONS: This study sheds light on Rab7 members in cotton, identifies a gene linked to drought stress, and paves the way for additional investigation of Rab7 genes associated with drought stress tolerance.

GhJUB1_3-At positively regulate drought and salt stress tolerance under control of GhHB7, GhRAP2-3 and GhRAV1 in Cotton.

Climate change severely affects crop production. Cotton is one of the primary fiber crops in the world and its production is susceptible to various environmental stresses, especially drought and salinity. Development of stress tolerant genotypes is the only way to escape from these environmental constraints. We identified sixteen homologs of the Arabidopsis JUB1 gene in cotton. Expression of GhJUB1_3-At was significantly induced in the temporal expression analysis of GhJUB1 genes in the roots of drought tolerant (H177) and susceptible (S9612) cotton genotypes under drought. The silencing of the GhJUB1_3-At gene alone and together with its paralogue GhJUB1_3-Dt reduced the drought tolerance in cotton plants. The transgenic lines exhibited tolerance to the drought and salt stress as compared to the wildtype (WT). The chlorophyll and relative water contents of wildtype decreased under drought as compared to the transgenic lines. The transgenic lines showed decreased H2O2 and increased proline levels under drought and salt stress, as compared to the WT, indicating that the transgenic lines have drought and salt stress tolerance. The expression analysis of the transgenic lines and WT revealed that GAI was upregulated in the transgenic lines in normal conditions as compared to the WT. Under drought and salt treatment, RAB18 and RD29A were strongly upregulated in the transgenic lines as compared to the WT. Conclusively, GhJUB1_3-At is not an auto activator and it is regulated by the crosstalk of GhHB7, GhRAP2-3 and GhRAV1. GhRAV1, a negative regulator of abiotic stress tolerance and positive regulator of leaf senescence, suppresses the expression of GhJUB1_3-At under severe circumstances leading to plant death.

Calcium-Dependent Protein Kinase GhCDPK16 Exerts a Positive Regulatory Role in Enhancing Drought Tolerance in Cotton.

Cotton is essential for the textile industry as a primary source of natural fibers. However, environmental factors like drought present significant challenges to its cultivation, adversely affecting both production levels and fiber quality. Enhancing cotton's drought resilience has the potential to reduce yield losses and support the growth of cotton farming. In this study, the cotton calcium-dependent protein kinase GhCDPK16 was characterized, and the transcription level of GhCDPK16 was significantly upregulated under drought and various stress-related hormone treatments. Physiological analyses revealed that the overexpression of GhCDPK16 improved drought stress resistance in Arabidopsis by enhancing osmotic adjustment capacity and boosting antioxidant enzyme activities. In contrast, silencing GhCDPK16 in cotton resulted in increased dehydration compared with the control. Furthermore, reduced antioxidant enzyme activities and downregulation of ABA-related genes were observed in GhCDPK16-silenced plants. These findings not only enhanced our understanding of the biological functions of GhCDPK16 and the mechanisms underlying drought stress resistance but also underscored the considerable potential of GhCDPK16 in improving drought resilience in cotton.

Genome-wide identification and mining elite allele variation of the Monoacylglycerol lipase (MAGL) gene family in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Monoacylglycerol lipase (MAGL) genes belong to the alpha/beta hydrolase superfamily, catalyze the terminal step of triglyceride (TAG) hydrolysis, converting monoacylglycerol (MAG) into free fatty acids and glycerol. RESULTS: In this study, 30 MAGL genes in upland cotton have been identified, which have been classified into eight subgroups. The duplication of GhMAGL genes in upland cotton was predominantly influenced by segmental duplication events, as revealed through synteny analysis. Furthermore, all GhMAGL genes were found to contain light-responsive elements. Through comprehensive association and haplotype analyses using resequencing data from 355 cotton accessions, GhMAGL3 and GhMAGL6 were detected as key genes related to lipid hydrolysis processes, suggesting a negative regulatory effect. CONCLUSIONS: In summary, MAGL has never been studied in upland cotton previously. This study provides the genetic mechanism foundation for the discover of new genes involved in lipid metabolism to improve cottonseed oil content, which will provide a strategic avenue for marker-assisted breeding aimed at incorporating desirable traits into cultivated cotton varieties.

Overexpression and knockdown of cotton GhdadD gene reveals its drought and salt stress tolerance role.

The 5'-deoxyadenosine deaminase (DADD), a member of the amidohydrolase family regulates biological purine metabolism. In this study, bioinformatic analysis, overexpression and knockdown of GhdadD gene were detected to identify its potential role in drought and salt stress tolerance. The results revealed that GhdadD was induced by ABA, Auxin, MBS and light responsive elements. In transgenic Arabidopsis, seed germination rate and root length were increased under drought or salt stress. GhdadD overexpressed seedlings resulted in higher plant height, less leaf damage and lower ion permeability. The expression of osmotic stress and ABA-responsive genes were up regulated. While in GhdadD-silenced cotton seedlings, CAT, SOD activity and soluble sugar content were reduced, MDA content was increased, and the stoma opening was depressed under drought or salt stress. Some osmics stress marker genes were also up regulated. These data indicating that GhdadD enhanced plant resistance to drought and salt stress through ABA pathways.

GhVOZ1-AVP1 module positively regulates salt tolerance in upland cotton (Gossypium hirsutum L.).

Vascular Plant One­zinc Finger (VOZ) transcription factor can respond to a variety of abiotic stresses, however its function in cotton and the molecular mechanisms of response to salt tolerance remained unclear. In this study, we found that GhVOZ1 is highly expressed in stamen and stem of cotton under normal conditions. The expression of GhVOZ1 increased significantly after 3 h of salt treatment in three-leaf staged upland cotton. Overexpressed transgenic lines of GhVOZ1 in Arabidopsis and upland cotton were treated with salt stress and we found that GhVOZ1 could respond positively to salt stress. GhVOZ1 can regulate Arabidopsis Vacuolar Proton Pump Pyrophosphatase (H+-PPase) gene (AVP1) expression through specific binding to GCGTCTAAAGTACGC site on GhAVP1 promoter, which was examined through Dual-luciferase assay and Electrophoretic mobility shift assay (EMSA). AVP1 expression was significantly increased in Arabidopsis with GhVOZ1 overexpression, while GhAVP1 expression was decreased in virus induced gene silenced (VIGS) cotton plants of GhVOZ1. Knockdown of GhAVP1 expression in cotton plants by VIGS showed decreased superoxide dismutase (SOD) and peroxidase (POD) activities, whereas an increased malondialdehyde (MDA) content and ultimately decreased salt tolerance. The GhVOZ1-AVP1 module could maintain sodium ion homeostasis through cell ion transport and positively regulate the salt tolerance in cotton, providing new ideas and insights for the study of salt tolerance.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus.

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

GhMYB30-GhMUR3 affects fiber elongation and secondary wall thickening in cotton.

Xyloglucan, an important hemicellulose, plays a crucial role in maintaining cell wall structure and cell elongation. However, the effects of xyloglucan on cotton fiber development are not well understood. GhMUR3 encodes a xyloglucan galactosyltransferase that is essential for xyloglucan synthesis and is highly expressed during fiber elongation. In this study, we report that GhMUR3 participates in cotton fiber development under the regulation of GhMYB30. Overexpression GhMUR3 affects the fiber elongation and cell wall thickening. Transcriptome showed that the expression of genes involved in secondary cell wall synthesis was prematurely activated in OE-MUR3 lines. In addition, GhMYB30 was identified as a key regulator of GhMUR3 by Y1H, Dual-Luc, and electrophoretic mobility shift assay (EMSA) assays. GhMYB30 directly bound the GhMUR3 promoter and activated GhMUR3 expression. Furthermore, DAP-seq of GhMYB30 was performed to identify its target genes in the whole genome. The results showed that many target genes were associated with fiber development, including cell wall synthesis-related genes, BR-related genes, reactive oxygen species pathway genes, and VLCFA synthesis genes. It was demonstrated that GhMYB30 may regulate fiber development through multiple pathways. Additionally, GhMYB46 was confirmed to be a target gene of GhMYB30 by EMSA, and GhMYB46 was significantly increased in GhMYB30-silenced lines, indicating that GhMYB30 inhibited GhMYB46 expression. Overall, these results revealed that GhMUR3 under the regulation of GhMYB30 and plays an essential role in cotton fiber elongation and secondary wall thickening. Additionally, GhMYB30 plays an important role in the regulation of fiber development and regulates fiber secondary wall synthesis by inhibiting the expression of GhMYB46.

Genome-wide analysis of SET domain genes and the function of GhSDG51 during salt stress in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Cotton, being extensively cultivated, holds immense economic significance as one of the most prominent crops globally. The SET (Su(var), E, and Trithorax) domain-containing protein is of significant importance in plant development, growth, and response to abiotic stress by modifying the lysine methylation status of histone. However, the comprehensive identification of SET domain genes (SDG) have not been conducted in upland cotton (Gossypium hirsutum L.). RESULTS: A total of 229 SDGs were identified in four Gossypium species, including G. arboretum, G. raimondii, G. hirsutum, and G. barbadense. These genes could distinctly be divided into eight groups. The analysis of gene structure and protein motif revealed a high degree of conservation among the SDGs within the same group. Collinearity analysis suggested that the SDGs of Gossypium species and most of the other selected plants were mainly expanded by dispersed duplication events and whole genome duplication (WGD) events. The allopolyploidization event also has a significant impact on the expansion of SDGs in tetraploid Gossypium species. Furthermore, the characteristics of these genes have been relatively conserved during the evolution. Cis-element analysis revealed that GhSDGs play a role in resistance to abiotic stresses and growth development. Furthermore, the qRT-PCR results have indicated the ability of GhSDGs to respond to salt stress. Co-expression analysis revealed that GhSDG51 might co-express with genes associated with salt stress. In addition, the silencing of GhSDG51 in cotton by the virus-induced gene silencing (VIGS) method suggested a potential positive regulatory role of GhSDG51 in salt stress. CONCLUSIONS: The results of this study comprehensively analyze the SDGs in cotton and provide a basis for understanding the biological role of SDGs in the stress resistance in upland cotton.

Identification of novel candidate loci and genes for seed vigor-related traits in upland cotton (Gossypium hirsutum L.) via GWAS.

Seed vigor (SV) is a crucial trait determining the quality of crop seeds. Currently, over 80% of China's cotton-planting area is in Xinjiang Province, where a fully mechanized planting model is adopted, accounting for more than 90% of the total fiber production. Therefore, identifying SV-related loci and genes is crucial for improving cotton yield in Xinjiang. In this study, three seed vigor-related traits, including germination potential, germination rate, and germination index, were investigated across three environments in a panel of 355 diverse accessions based on 2,261,854 high-quality single-nucleotide polymorphisms (SNPs). A total of 26 significant SNPs were detected and divided into six quantitative trait locus regions, including 121 predicted candidate genes. By combining gene expression, gene annotation, and haplotype analysis, two novel candidate genes (Ghir_A09G002730 and Ghir_D03G009280) within qGR-A09-1 and qGI/GP/GR-D03-3 were associated with vigor-related traits, and Ghir_A09G002730 was found to be involved in artificial selection during cotton breeding by population genetic analysis. Thus, understanding the genetic mechanisms underlying seed vigor-related traits in cotton could help increase the efficiency of direct seeding by molecular marker-assisted selection breeding.

A systematic analysis of the phloem protein 2 (PP2) proteins in Gossypium hirsutum reveals that GhPP2-33 regulates salt tolerance.

BACKGROUND: Phloem protein 2 (PP2) proteins play a vital role in the Phloem-based defense (PBD) and participate in many abiotic and biotic stress. However, research on PP2 proteins in cotton is still lacking. RESULTS: A total of 25, 23, 43, and 47 PP2 genes were comprehensively identified and characterized in G.arboretum, G.raimondii, G.barbadense, and G.hirsutum. The whole genome duplication (WGD) and allopolyploidization events play essential roles in the expansion of PP2 genes. The promoter regions of GhPP2 genes contain many cis-acting elements related to abiotic stress and the weighted gene co-expression network analysis (WGCNA) analysis displayed that GhPP2s could be related to salt stress. The qRT-PCR assays further confirmed that GhPP2-33 could be dramatically upregulated during the salt treatment. And the virus-induced gene silencing (VIGS) experiment proved that the silencing of GhPP2-33 could decrease salt tolerance. CONCLUSIONS: The results in this study not only offer new perspectives for understanding the evolution of PP2 genes in cotton but also further explore their function under salt stress.

GhCDPK60 positively regulates drought stress tolerance in both transgenic Arabidopsis and cotton by regulating proline content and ROS level.

Calcium-Dependent Protein Kinases (CDPKs) involved in regulating downstream components of calcium signaling pathways play a role in tolerance to abiotic stresses and seed development in plants. However, functions of only a few cotton CDPKs have been clarified at present. In this study, 80 conserved CDPKs in Gossypium hirsutum L. were identified and characterized, which was divided into four subgroups. Among them, the transcript level of GhCDPK60 was significantly upregulated under drought and several hormone treatments. And we found that the expression levels of several stress-inducible genes down-regulated in GhCDPK60-silence cotton and up-regulated in GhCDPK60-overexpressing Arabidopsis. In addition, physiological analyses demonstrated that GhCDPK60 improved drought stress tolerance by improving the osmotic adjustment ability and reducing the accumulation of reactive oxygen species (ROS) in plants. These findings broaden our understanding of the biological roles of GhCDPK60 and mechanisms underlying drought stress tolerance in cotton.

A Comprehensive Gene Co-Expression Network Analysis Reveals a Role of GhWRKY46 in Responding to Drought and Salt Stresses.

Abiotic stress, such as drought and salinity stress, seriously inhibit the growth and development of plants. Therefore, it is vital to understand the drought and salinity resistance mechanisms to enable cotton to provide more production under drought and salt conditions. In this study, we identified 8806 and 9108 differentially expressed genes (DEGs) through a comprehensive analysis of transcriptomic data related to the PEG-induced osmotic and salt stress in cotton. By performing weighted gene co-expression network analysis (WGCNA), we identified four co-expression modules in PEG treatment and five co-expression modules in salinity stress, which included 346 and 324 predicted transcription factors (TFs) in these modules, respectively. Correspondingly, whole genome duplication (WGD) events mainly contribute to the expansion of those TFs. Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) analyses revealed those different modules were associated with stress resistance, including regulating macromolecule metabolic process, peptidase activity, transporter activity, lipid metabolic process, and responses to stimulus. Quantitative RT-PCR analysis was used to confirm the expression levels of 15 hub TFs in PEG6000 and salinity treatments. We found that the hub gene GhWRKY46 could alter salt and PEG-induced drought resistance in cotton through the virus-induced gene silencing (VIGS) method. Our results provide a preliminary framework for further investigation of the cotton response to salt and drought stress, which is significant to breeding salt- and drought-tolerant cotton varieties.

Insights into the Genomic Regions and Candidate Genes of Senescence-Related Traits in Upland Cotton via GWAS.

Senescence is the last stage of plant development and is controlled by both internal and external factors. Premature senescence significantly affects the yield and quality of cotton. However, the genetic architecture underlying cotton senescence remains unclear. In this study, genome-wide association studies (GWAS) were performed based on 3,015,002 high-quality SNP markers from the resequencing data of 355 upland cotton accessions to detect genomic regions for cotton senescence. A total of 977 candidate genes within 55 senescence-related genomic regions (SGRs), SGR1-SGR55, were predicted. Gene ontology (GO) analysis of candidate genes revealed that a set of biological processes was enriched, such as salt stress, ethylene processes, and leaf senescence. Furthermore, in the leaf senescence GO term, one candidate gene was focused on: Gohir.A12G270900 (GhMKK9), located in SGR36, which encodes a protein of the MAP kinase kinase family. Quantitative real-time PCR (qRT-PCR) analysis showed that GhMKK9 was up-regulated in old cotton leaves. Overexpression of GhMKK9 in Arabidopsis accelerated natural leaf senescence. Virus-induced gene silencing (VIGS) of GhMKK9 in cotton increased drought tolerance. These results suggest that GhMKK9 is a positive regulator and might be involved in drought-induced senescence in cotton. The results provide new insights into the genetic basis of cotton senescence and will be useful for improving cotton breeding in the future.

Genomewide Identification and Characterization of the Genes Involved in the Flowering of Cotton.

Flowering is a prerequisite for flowering plants to complete reproduction, and flowering time has an important effect on the high and stable yields of crops. However, there are limited reports on flowering-related genes at the genomic level in cotton. In this study, genomewide analysis of the evolutionary relationship of flowering-related genes in different cotton species shows that the numbers of flowering-related genes in the genomes of tetraploid cotton species Gossypium hirsutum and Gossypium barbadense were similar, and that these numbers were approximately twice as much as the number in diploid cotton species Gossypium arboretum. The classification of flowering-related genes shows that most of them belong to the photoperiod and circadian clock flowering pathway. The distribution of flowering-related genes on the chromosomes of the At and Dt subgenomes was similar, with no subgenomic preference detected. In addition, most of the flowering-related core genes in Arabidopsis thaliana had homologs in the cotton genome, but the copy numbers and expression patterns were disparate; moreover, flowering-related genes underwent purifying selection throughout the evolutionary and selection processes. Although the differentiation and reorganization of many key genes of the cotton flowering regulatory network occurred throughout the evolutionary and selection processes, most of them, especially those involved in the important flowering regulatory networks, have been relatively conserved and preferentially selected.

Detection of Stable Elite Haplotypes and Potential Candidate Genes of Boll Weight Across Multiple Environments via GWAS in Upland Cotton.

Boll weight (BW) is a key determinant of yield component traits in cotton, and understanding the genetic mechanism of BW could contribute to the progress of cotton fiber yield. Although many yield-related quantitative trait loci (QTLs) responsible for BW have been determined, knowledge of the genes controlling cotton yield remains limited. Here, association mapping based on 25,169 single-nucleotide polymorphisms (SNPs) and 2,315 insertions/deletions (InDels) was conducted to identify high-quality QTLs responsible for BW in a global collection of 290 diverse accessions, and BW was measured in nine different environments. A total of 19 significant markers were detected, and 225 candidate genes within a 400 kb region (± 200 kb surrounding each locus) were predicted. Of them, two major QTLs with highly phenotypic variation explanation on chromosomes A08 and D13 were identified among multiple environments. Furthermore, we found that two novel candidate genes (Ghir_A08G009110 and Ghir_D13G023010) were associated with BW and that Ghir_D13G023010 was involved in artificial selection during cotton breeding by population genetic analysis. The transcription level analyses showed that these two genes were significantly differentially expressed between high-BW accession and low-BW accession during the ovule development stage. Thus, these results reveal valuable information for clarifying the genetic basics of the control of BW, which are useful for increasing yield by molecular marker-assisted selection (MAS) breeding in cotton.

Characterization and Functional Analysis of GhNAC82, A NAM Domain Gene, Coordinates the Leaf Senescence in Upland Cotton (Gossypium hirsutum L.).

In the process of growth and development, cotton exhibits premature senescence under various abiotic stresses, impairing yield and fiber quality. NAC (NAM, ATAF1,2, and CUC2) protein widely distributed in land plants, play the critical role in responding to abiotic stress and regulating leaf senescence. We have identified and functional analyzed a NAM domain gene GhNAC82 in upland cotton, it was located on the A11 chromosome 4,921,702 to 4,922,748 bp, only containing one exon. The spatio-temporal expression pattern analysis revealed that it was highly expressed in root, torus, ovule and fiber development stage. The results of qRT-PCR validated that GhNAC82 negatively regulated by salt stress, drought stress, H2O2 stress, IAA treatment, and ethylene treatment, positively regulated by the ABA and MeJA treatment. Moreover, heterologous overexpression of GhNAC82 results in leaf premature senescence and delays root system development in Arabidopsis thaliana. The phenotype of delayed-senescence was performed after silencing GhNAC82 by VIGS in premature cotton. Taken together, GhNAC82 was involved in different abiotic stress pathways and play important roles in negatively regulating leaf premature senescence.

A Comprehensive Identification and Function Analysis of Serine/Arginine-Rich (SR) Proteins in Cotton (Gossypium spp.).

As one of the most important factors in alternative splicing (AS) events, serine/arginine-rich (SR) proteins not only participate in the growth and development of plants but also play pivotal roles in abiotic stresses. However, the research about SR proteins in cotton is still lacking. In this study, we performed an extensive comparative analysis of SR proteins and determined their phylogeny in the plant lineage. A total of 169 SR family members were identified from four Gossypium species, and these genes could be divided into eight distinct subfamilies. The domain, motif distribution and gene structure of cotton SR proteins are conserved within each subfamily. The expansion of SR genes is mainly contributed by WGD and allopolyploidization events in cotton. The selection pressure analysis showed that all the paralogous gene pairs were under purifying selection pressure. Many cis-elements responding to abiotic stress and phytohormones were identified in the upstream sequences of the GhSR genes. Expression profiling suggested that some GhSR genes may involve in the pathways of plant resistance to abiotic stresses. The WGCNA analysis showed that GhSCL-8 co-expressed with many abiotic responding related genes in a salt-responding network. The Y2H assays showed that GhSCL-8 could interact with GhSRs in other subfamilies. The subcellular location analysis showed that GhSCL-8 is expressed in the nucleus. The further VIGS assays showed that the silencing of GhSCL-8 could decrease salt tolerance in cotton. These results expand our knowledge of the evolution of the SR gene family in plants, and they will also contribute to the elucidation of the biological functions of SR genes in the future.

Overexpression of a Cotton Aquaporin Gene GhTIP1;1-like Confers Cold Tolerance in Transgenic Arabidopsis.

Cold stress can significantly affect the development, yield, and quality of crops and restrict the geographical distribution and growing seasons of plants. Aquaporins are the main channels for water transport in plant cells. Abiotic stresses such as cold and drought dehydrate cells by changing the water potential. In this study, we cloned a gene GhTIP1;1-like encodes tonoplast aquaporin from the transcriptome database of cotton seedlings after cold stress. Expression analysis showed that GhTIP1;1-like not only responds to cold stress but was also induced by heat, drought and salt stress. Subcellular localization showed that the protein was anchored to the vacuole membrane. Promoter deletion analysis revealed that a MYC motif within the promoter region of GhTIP1;1-like were the core cis-elements in response to low temperature. Virus-induced gene silencing (VIGS) and histochemical staining indicate that GhTIP1;1-like plays a positive role in plant cold tolerance. Overexpression of GhTIP1;1-like in Arabidopsis delayed the senescence process and enhanced the cold tolerance of transgenic plants. Compared with the wild type, the soluble protein concentration and peroxidase activity of the transgenic lines under cold stress were higher, while the malondialdehyde content was lower. In addition, the expression levels of cold-responsive genes were significantly increased in transgenic plants under cold stress. Our results indicate that GhTIP1;1-like could respond to different abiotic stresses and be positively involved in regulating the cold tolerance of cotton.