Search

Drought stress significantly affects crop productivity. Carotenoids are essential photosynthetic pigment for plants, bacteria, and algae, with signaling and antioxidant functions. Lutein is a crucial branch product in the carotenoid synthesis pathway, which effectively improves the stress tolerance of higher plants. lycopene cyclase, a central enzyme for lutein synthesis, holds great significance in regulating lutein production. This research establishes a correlation between lutein content and stress resistance by measuring the drought resistance and lutein content of various cotton materials. To identify which crucial genes are associated with lutein, the lycopene cyclase family (LCYs) was analyzed. The research found that LCYs form a highly conserved family divided into two subfamilies, LCY-ε (lycopene ε-cyclase) and LCY-β (lycopene β-cyclase). Most members of the LCY family contain photoresponsive elements and abscisic acid elements. qRT-PCR demonstrates showed that most genes responded positively to drought stress, and GhLCYε-3 was expressed significantly differently under drought stress. Virus-induced gene silencing (VIGS) assay showed that the content of GhLCYε-3 was significantly increased with MDA and PRO, and the contents of chlorophyll and lutein were significantly decreased in pYL156 plants. The decrease in GhLCYε-3 expression is speculated to lead to reduced lutein content in vivo, resulting in the accumulation of reactive oxygen species (ROS) and decreased drought tolerance. This research enriched the understanding of LCY gene family and lutein function, and provided a new reference for cotton planting in arid areas. Synopsis: Lycopene cyclase plays an important role in enhancing the ability of scavenging ROS and drought resistance of plants.

Abstract Vitamin C, also referred to as ascorbic acid (AsA), is recognized for its capacity to cure and avert scurvy, and it is crucial for regular human growth and development. In various crops, AsA participates in stress response mechanisms mediated by abscisic acid and has been discovered to have a crucial function in the morphogenesis, growth, development, and production of male gametes in plants. GDP-D-mannose 3′,5′-epimerase (GME) is essential in the synthesis of vitamin C. Our research identified 91, 83, 51, and 46 genes, respectively, found in G. barbadense (GbGMEs), G. hirsutum (GhGMEs), G. arboretum (GaGMEs), and G. raimondii (GrGMEs). Plants resulting from VIGS infection with GhGME31D clearly showed yellowing, water loss and wilting of leaves and black spots on stems. Measurement of MDA and AsA levels indicated that the plants were more damaged. This indicates that AsA has a substantial impact on plant growth and development.

Abstract The VIM (belonged to E3 ubiquitin ligase) gene family is crucial for plant growth, development, and stress responses, yet their role in salt stress remains unclear. We analyzed phylogenetic relationships, chromosomal localization, conserved motifs, gene structure, cis-acting elements, and gene expression patterns of the VIM gene family in four cotton varieties. Our findings reveal 29, 29, 17, and 14 members in Gossypium hirsutum (G.hirsutum), Gossypium barbadense (G.barbadense), Gossypium arboreum (G.arboreum), and Gossypium raimondii (G. raimondii), respectively, indicating the maturity and evolution of this gene family. motifs among GhVIMs genes were observed, along with the presence of stress-responsive, hormone-responsive, and growth-related elements in their promoter regions. Gene expression analysis showed varying patterns and tissue specificity of GhVIMs genes under abiotic stress. Silencing GhVIM28 via virus-induced gene silencing revealed its role as a salt-tolerant negative regulator. This work reveals a mechanism by which the VIM gene family in response to salt stress in cotton, identifying a potential negative regulator, GhVIM28, which could be targeted for enhancing salt tolerance in cotton. The objective of this study was to explore the evolutionary relationship of the VIM gene family and its potential function in salt stress tolerance, and provide important genetic resources for salt tolerance breeding of cotton.

Abstract Background DNA methylation is an important epigenetic mode of genomic DNA modification and plays a vital role in maintaining epigenetic content and regulating gene expression. Cytosine-5 DNA methyltransferase (C5-MTase) are the key enzymes in the process of DNA methylation. However, there is no systematic analysis of the C5-MTase in cotton so far, and the function of DNMT2 genes has not been studied. Methods In this study, the whole genome of cotton C5-MTase coding genes was identified and analyzed using a bioinformatics method based on information from the cotton genome, and the function of GhDMT6 was further validated by VIGS experiments and subcellular localization analysis. Results 33 C5-MTases were identified from three cotton genomes, and were divided into four subfamilies by systematic evolutionary analysis. After the protein domain alignment of C5-MTases in cotton, 6 highly conserved motifs were found in the C-terminus of 33 proteins involved in methylation modification, which indicated that C5-MTases had a basic catalytic methylation function. These proteins were divided into four classes based on the N-terminal difference, of which DNMT2 lacks the N-terminal regulatory domain. The expression of C5-MTases in different parts of cotton was different under different stress treatments, which indicated the functional diversity of cotton C5-MTase gene family. Among the C5-MTases, the GhDMT6 had a obvious up-regulated expression. After silencing GhDMT6 with VIGS, the phenotype of cotton seedlings under different stress treatments showed a significant difference. Compared with cotton seedlings that did not silence GhDMT6, cotton seedlings silencing GhDMT6 showed significant stress resistance. Conclusion The results show that C5-MTases plays an important role in cotton stress response, which is beneficial to further explore the function of DNMT2 subfamily genes.

Abstract Cytokinin oxidase/dehydrogenase (CKX) is a key regulatory enzyme for the irreversible degradation of the plant hormone cytokinin (CK), which is important in growth and development and response to abiotic stresses in cotton. In this study, 27, 28, 14 and 14 CKXs were screened by FAD structural domain and cytokinin binding structural domain in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii, respectively. Their phylogenetic relationships and expression patterns were analyzed, and most GhCKXs were found to be tissue-specific and responsive to various abiotic stresses such as cold, heat, salt and PEG. GhCKX6b-Dt was selected for gene silencing in evolutionary branch II for salt stress, because its expression increased after salt stress in cotton plants. An increase in PRO and MDA content and a decrease in SOD activity due to this gene were found after inducing salt stress, contributing to oxidative damage and decreased salt tolerance. In this study, CKXs were analyzed to reveal the possible role of GhCKXs against abiotic stresses in cotton, which provides a basis for further understanding of the biological functions of CK in plants such as growth and development and stress resistance.

Abstract Background Inositol monophosphates (IMP) are key enzymes in the ascorbic acid (AsA) synthesis pathways, which play vital roles in regulating plant growth and development and stresses tolerance. To date, no comprehensive analysis of the expression profile of IMP genes and their functions under abiotic stress in cotton has been reported. Results In this study, the genetic characteristics, phylogenetic evolution, cis-acting elements and expression patterns of IMP gene family in cotton were systematically analyzed. A total of 28, 27, 13 and 13 IMP genes were identified in Gossypium hirsutum (G. hirsutum), Gossypium barbadense (G. barbadense), Gossypium arboreum (G. arboreum), and Gossypium raimondii (G. raimondii), respectively. Phylogenetic analysis showed that IMP family genes could cluster into 3 clades. Structure analysis of genes showed that GhIMP genes from the same subgroup had similar genetic structure and exon number. And most GhIMP family members contained hormone-related elements (abscisic acid response element, MeJA response element, gibberellin response element) and stress-related elements (low temperature response element, defense and stress response element, wound response element). After exogenous application of abscisic acid (ABA), some GhIMP genes containing ABA response elements positively responded to alkaline stress, indicating that ABA response elements played an important role in response to alkaline stress. qRT-PCR showed that most of GhIMP genes responded positively to alkaline stress, and GhIMP10D significantly upregulated under alkaline stress, with the highest up-regulated expression level. Virus-induced gene silencing (VIGS) experiment showed that compared with 156 plants, MDA content of pYL156:GhIMP10D plants increased significantly, while POD, SOD, chlorophyII and AsA content decreased significantly. Conclusions This study provides a thorough overview of the IMP gene family and presents a new perspective on the evolution of this gene family. In particular, some IMP genes may be involved in alkaline stress tolerance regulation, and GhIMP10D showed high expression levels in leaves, stems and roots under alkaline stress, and preliminary functional verification of GhIMP10D gene suggested that it may regulate tolerance to alkaline stress by regulating the activity of antioxidant enzymes and the content of AsA. This study contributes to the subsequent broader discussion of the structure and alkaline resistance of IMP genes in cotton.

Abstract Background Flooding is among the most severe abiotic stresses in plant growth and development. The mechanism of submergence tolerance of cotton in response to submergence stress is unknown. Results The transcriptome results showed that a total of 6,893 differentially expressed genes (DEGs) were discovered under submergence stress. Gene Ontology (GO) enrichment analysis showed that DEGs were involved in various stress or stimulus responses. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that DEGs related to plant hormone signal transduction, starch and sucrose metabolism, glycolysis and the biosynthesis of secondary metabolites were regulated by submergence stress. Eight DEGs related to ethylene signaling and 3 ethylene synthesis genes were identified in the hormone signal transduction. For respiratory metabolism, alcohol dehydrogenase (ADH, GH_A02G0728) and pyruvate decarboxylase (PDC, GH_D09G1778) were significantly upregulated but 6-phosphofructokinase (PFK, GH_D05G0280), phosphoglycerate kinase (PGK, GH_A01G0945 and GH_D01G0967) and sucrose synthase genes (SUS, GH_A06G0873 and GH_D06G0851) were significantly downregulated in the submergence treatment. Terpene biosynthetic pathway-related genes in the secondary metabolites were regulated in submergence stress. Conclusions Regulation of terpene biosynthesis by respiratory metabolism may play a role in enhancing the tolerance of cotton to submergence under flooding. Our findings showed that the mevalonate pathway, which occurs in the cytoplasm of the terpenoid backbone biosynthesis pathway (ko00900), may be the main response to submergence stress.

Abstract Background Cotton is an important industrial crop and a pioneer crop for saline-alkali land restoration. However, the molecular mechanism underlying the cotton response to salt is not completely understood. Methods Here, we used metabolome data and transcriptome data to analyze the salt tolerance regulatory network of cotton and metabolic biomarkers. Results In this study, cotton was stressed at 400 m M NaCl for 0 h, 3 h, 24 h and 48 h. NaCl interfered with cotton gene expression, altered metabolite contents and affected plant growth. Metabolome analysis showed that NaCl stress increased the contents of amino acids, sugars and ABA, decreased the amount of vitamin and terpenoids. K-means cluster analysis of differentially expressed genes showed that the continuously up-regulated genes were mainly enriched in metabolic pathways such as flavonoid biosynthesis and amino acid biosynthesis. Conclusion The four metabolites of cysteine (Cys), ABA(Abscisic acid), turanose, and isopentenyladenine-7-N-glucoside (IP7G) were consistently up-regulated under salt stress, which may indicate that they are potential candidates for cotton under salt stress biomarkers. Combined transcriptome and metabolome analysis revealed accumulation of cysteine, ABA, isopentenyladenine-7-N-glucoside and turanose were important for salt tolerance in cotton mechanism. These results will provide some metabolic insights and key metabolite biomarkers for salt stress tolerance, which may help to understanding of the metabolite response to salt stress in cotton and develop a foundation for cotton to grow better in saline soil.

Abstract Background Flavonols play important roles in antioxidation and anticancer activities, longevity, and cardiovascular protection. Flavonol synthase (FLS) is a key enzyme for flavonol synthesis. Result Phenotypic, transcriptional and metabolic data were analyzed, which showed that there was a close relationship between salt stress and flavonoids, and flavonols were significantly upregulated under salt stress. Nine, seven, four, and four FLS genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum, and Gossypium raimondii, respectively. The results of subcellular localization showed that FLS existed in the nucleus and cytoplasmic. Through phylogenetic analysis, 24 FLS genes were divided into three subfamilies. The results of the RNA sequencing showed that the expression of GhFLS genes was mainly induced by salt, drought, low temperature, and heat stress. GhFLS promoter mainly comprised plant hormone response elements and abiotic stress elements, indicating that the GhFLS gene may play a key role in abiotic stress response. The proline contents of pYL156:GhFLS1 was reduced significantly compared to pYL156 under salt stress, thereby reducing the resistance of cotton to salt stress. Conclusion This study will lay a foundation for further study on the antioxidant regulation mechanism of the FLS gene under abiotic stress.

A-galactosidases (AGALs), the oligosaccharide (RFO) catabolic genes of the raffinose family, play crucial roles in plant growth and development and in adversity stress. They can break down the non-reducing terminal galactose residues of glycolipids and sugar chains. In this study, the whole genome of AGALs was analyzed. Bioinformatics analysis was conducted to analyze members of the AGAL family in Gossypium hirsutum, Gossypium arboreum, Gossypium barbadense, and Gossypium raimondii. Meanwhile, RT-qPCR was carried out to analyze the expression patterns of AGAL family members in different tissues of terrestrial cotton. It was found that a series of environmental factors stimulated the expression of the GhAGAL3 gene. The function of GhAGAL3 was verified through virus-induced gene silencing (VIGS). As a result, GhAGAL3 gene silencing resulted in milder wilting of seedlings than the controls, and a significant increase in the raffinose content in cotton, indicating that GhAGAL3 responded to NaCl stress. The increase in raffinose content improved the tolerance of cotton. Findings in this study lay an important foundation for further research on the role of the GhAGAL3 gene family in the molecular mechanism of abiotic stress resistance in cotton.

Abstract Background 2-oxoglutarate-dependent dioxygenase (2ODD) is the second largest family of oxidases involved in various oxygenation/hydroxylation reactions in plants. Many members in the family regulate gene transcription, nucleic acid modification/repair and secondary metabolic synthesis. The 2ODD family genes also function in the formation of abundant flavonoids during anthocyanin synthesis, thereby modulating plant development and response to diverse stresses. Results Totally, 379, 336, 205, and 204 2ODD genes were identified in G. barbadense (Gb), G. hirsutum (Gh), G. arboreum (Ga), and G. raimondii (Gb), respectively. The 336 2ODDs in G. hirsutum were divided into 15 subfamilies according to their putative functions. The structural features and functions of the 2ODD members in the same subfamily were similar and evolutionarily conserved. Tandem duplications and segmental duplications served essential roles in the large-scale expansion of the cotton 2ODD family. Ka/Ks values for most of the gene pairs were less than 1, indicating that 2ODD genes undergo strong purifying selection during evolution. Gh2ODDs might act in cotton responses to different abiotic stresses. GhLDOX3 and GhLDOX7, two members of the GhLDOX subfamily from Gh2ODDs, were significantly down-regulated in transcription under alkaline stress. Moreover, the expression of GhLDOX3 in leaves was significantly higher than that in other tissues. These results will provide valuable information for further understanding the evolution mechanisms and functions of the cotton 2ODD genes in the future. Conclusions Genome-wide identification, structure, and evolution and expression analysis of 2ODD genes in Gossypium were carried out. The 2ODDs were highly conserved during evolutionary. Most Gh2ODDs were involved in the regulation of cotton responses to multiple abiotic stresses including salt, drought, hot, cold and alkali.

Abstract Background Abscisic acid (ABA) is an important stress hormone, the changes of abscisic acid content can alter plant tolerance to stress, abscisic acid is crucial for studying plant responses to abiotic stress. The abscisic acid aldehyde oxidase (AAO) plays a vital role in the final step in the synthesis of abscisic acid, therefore, understanding the function of AAO gene family is of great significance for plants to response to abiotic stresses. Result In this study, 6, 8, 4 and 4 AAO genes were identified in four cotton species. According to the structural characteristics of genes and the traits of phylogenetic tree, we divided the AAO gene family into 4 clades. Gene structure analysis showed that the AAO gene family was relatively conservative. The analysis of cis-elements showed that most AAO genes contained cis-elements related to light response and plant hormones. Tissue specificity analysis under NaHCO3 stress showed that GhAAO2 gene was differentially expressed in both roots and leaves. After GhAAO2 gene silencing, the degree of wilting of seedlings was lighter than that of the control group, indicating that GhAAO2 could respond to NaHCO3 stress. Conclusions In this study, the AAO gene family was analyzed by bioinformatics, the response of GhAAO gene to various abiotic stresses was preliminarily verified, and the function of the specifically expressed gene GhAAO2 was further verified. These findings provide valuable information for the study of potential candidate genes related to plant growth and stress.

Abstract Background Phosphorylation regulated by the two-component system (TCS) is a very important approach signal transduction in most of living organisms. Histidine phosphotransfer (HP) is one of the important members of the TCS system. Members of the HP gene family have implications in plant stresses tolerance and have been deeply studied in several crops. However, upland cotton is still lacking with complete systematic examination of the HP gene family. Results A total of 103 HP gene family members were identified. Multiple sequence alignment and phylogeny of HPs distributed them into 7 clades that contain the highly conserved amino acid residue “XHQXKGSSXS”, similar to the Arabidopsis HP protein. Gene duplication relationship showed the expansion of HP gene family being subjected with whole-genome duplication (WGD) in cotton. Varying expression profiles of HPs illustrates their multiple roles under altering environments particularly the abiotic stresses. Analysis is of transcriptome data signifies the important roles played by HP genes against abiotic stresses. Moreover, protein regulatory network analysis and VIGS mediated functional approaches of two HP genes (GhHP23 and GhHP27) supports their predictor roles in salt and drought stress tolerance. Conclusions This study provides new bases for systematic examination of HP genes in upland cotton, which formulated the genetic makeup for their future survey and examination of their potential use in cotton production.

Abstract Background Under abiotic stress conditions, cotton growth is inhibited and yield losses are severe. Identification of calnexin family members and function analysis under abiotic stress laid the foundation for the screening of stress-related candidate genes. Results A total of 60 CNX family members have been identified in Gossypium hirsutum, G. barbadense, G. arboreum, and G. raimondii, and they were divided into two categories: CNX and CRT genes. Through the construction of a phylogenetic tree, they were subdivided into three classes. Further analysis of chromosome localization, conserved promoters, gene structure and selection under pressure showed that the family members were highly conserved in the evolution process. Analysis of cis-acting elements in the promoter regions showed that CNX family genes contain regulatory elements for growth and development, anaerobic, drought, defense and stress response, and plant hormones. Using RNA-seq data to study the expression pattern of GhCNX genes under cold, hot, salt stress and Polyethylene glycol, it was observed that the gene expression levels changed by different degrees under different stress conditions, indicating that GhCNX members were involved in the regulation of multiple biological stresses. Conclusion This study provides an insight into the members of cotton CNX genes. The results of this study suggested that CNX family members play a role in defense against adversity and provide a foundation for the discovery of stress-related genes.

Cysteine, an early sulfur-containing compound in plants, is of significant importance in sulfur metabolism. CYS encodes cysteine synthetase that further catalyzes cysteine synthesis. In this investigation, CYS genes, identified from genome-wide analysis of Gossypium hirsutum bioinformatically, led to the discovery of GhCYS2 as the pivotal gene responsible for Cd2+ response. The silencing of GhCYS2 through virus-induced gene silencing (VIGS) rendered plants highly susceptible to Cd2+ stress. Silencing GhCYS2 in plants resulted in diminished levels of cysteine and glutathione while leading to the accumulation of MDA and ROS within cells, thereby impeding the regular process of photosynthesis. Consequently, the stomatal aperture of leaves decreased, epidermal cells underwent distortion and deformation, intercellular connections are dramatically disrupted, and fissures manifested between cells. Ultimately, these detrimental effected culminating in plant wilting and a substantial reduction in biomass. The association established between Cd2+ and cysteine in this investigation offered a valuable reference point for further inquiry into the functional and regulatory mechanisms of cysteine synthesis genes.

Abstract Background The cationic amino acid transporters (CAT) play indispensable roles in maintaining metabolic functions, such as synthesis of proteins and nitric oxide (NO), biosynthesis of polyamine, and flow of amino acids, by mediating the bidirectional transport of cationic amino acids in plant cells. Results In this study, we performed a genome-wide and comprehensive study of 79 CAT genes in four species of cotton. Localization of genes revealed that CAT genes reside on the plasma membrane. Seventy-nine CAT genes were grouped into 7 subfamilies by phylogenetic analysis. Structure analysis of genes showed that CAT genes from the same subgroup have similar genetic structure and exon number. RNA-seq and real-time PCR indicated that the expression of most GhCAT genes were induced by salt, drought, cold and heat stresses. Cis-elements analysis of GhCAT promoters showed that the GhCAT genes promoters mainly contained plant hormones responsive elements and abiotic stress elements, which indicated that GhCAT genes may play key roles in response to abiotic stress. Moreover, we also conducted gene interaction network of the GhCAT proteins. Silencing GhCAT10D expression decreased the resistance of cotton to salt stress because of a decrease in the accumulation of NO and proline. Conclusion Our results indicated that CAT genes might be related with salt tolerance in cotton and lay a foundation for further study on the regulation mechanism of CAT genes in cationic amino acids transporting and distribution responsing to abiotic stress.

Abstract Background Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS). Method From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed. Result A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H2O2 and other reactive oxygen removal systems. Conclusion This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Identification, evolution, and expression patterns of BSK (BR signaling kinase) family genes revealed that BSKs participated in the response of cotton to abiotic stress and maintained the growth of cotton in extreme environment. The steroidal hormone brassinosteroids (BR) play important roles in different plant biological processes. This study focused on BSK which were downstream regulatory element of BR, in order to help to decipher the functions of BSKs genes from cotton on growth development and responses to abiotic stresses and lean the evolutionary relationship of cotton BSKs. BSKs are a class of plant-specific receptor-like cytoplasmic kinases involved in BR signal transduction. In this study, bioinformatics methods were used to identify the cotton BSKs gene family at the cotton genome level, and the gene structure, promoter elements, protein structure and properties, gene expression patterns and candidate interacting proteins were analyzed. In the present study, a total of 152 BSKs were identified by a genome-wide search in four cotton species and other 11 plant species, and phylogenetic analysis revealed three evolutionary clades. It was identified that BSKs contain typical PKc and TPR domains, the N-terminus is composed of extended chains and helical structures. Cotton BSKs genes show different expression patterns in different tissues and organs. The gene promoter contains numerous cis-acting elements induced by hormones and abiotic stress, the hormone ABA and Cold-inducing related elements have the highest count, indicating that cotton BSK genes may be regulated by various hormones at different growth stages and involved in the response regulation of cotton to various stresses. The expression analysis of BSKs in cotton showed that the expression levels of GhBSK06, GhBSK10, GhBSK21 and GhBSK24 were significantly increased with salt-inducing. This study is helpful to analyze the function of cotton BSKs genes in growth and development and in response to stress.

Abstract Background Carboxylesterase (CXE) is a type of hydrolase with α/β sheet hydrolase activity widely found in animals, plants and microorganisms, which plays an important role in plant growth, development and resistance to stress. Results A total of 72, 74, 39, 38 CXE genes were identified in Gossypium barbadense, Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively. The gene structure and expression pattern were analyzed. The GBCXE genes were divided into 6 subgroups, and the chromosome distribution of members of the family were mapped. Analysis of promoter cis-acting elements showed that most GBCXE genes contain cis-elements related to plant hormones (GA, IAA) or abiotic stress. These 6 genes we screened out were expressed in the root, stem and leaf tissues. Combined with the heat map, GBCXE49 gene was selected for subcellular locate and confirmed that the protein was expressed in the cytoplasm. Conclusions The collinearity analysis of the CXE genes of the four cotton species in this family indicated that tandem replication played an indispensable role in the evolution of the CXE gene family. The expression patterns of GBCXE gene under different stress treatments indicated that GBCXE gene may significantly participate in the response to salt and alkaline stress through different mechanisms. Through the virus-induced gene silencing technology (VIGS), it was speculated that GBCXE49 gene was involved in the response to alkaline stress in G. barbadense.

Abstract Background Nucleobase–ascorbate transporters (NAT), synonymously called nucleobase–cation symporter 2 (NCS2) proteins, were earlier reported to be involved in plant growth, development and resistance to stress. Previous studies concluded that s a polymorphic SNP associated with NAT12 was significant different between salt-tolerant and salt-sensitive materials of upland cotton. In current study, a comprehensive analysis of NAT family genes was conducted for the first time in cotton. Results In this study, we discovered 32, 32, 18, and 16 NAT genes in Gossypium hirsutum, G. barbadense, G. raimondii and G. arboreum, respectively, which were classified into four groups (groups I–IV) based on the multiple sequence analysis. These GhNAT genes were unevenly distributed on At and Dt sub-genome in G. hirsutum. Most GhNAT members in the same group had similar gene structure characteristics and motif composition. The collinearity analysis revealed segmental duplication as well as tandem duplication contributing to the expansion of the GhNATs. The analysis of cis-acting regulatory elements of GhNATs showed that the function of GhNAT genes in cotton might be related to plant hormone and stress response. Under different conditions, the expression levels further suggested the GhNAT family genes were associated with plant response to various abiotic stresses. GhNAT12 was detected in the plasma membrane. And it was validated that the GhNAT12 gene played an important role in regulating cotton resistance to salt and drought stress through the virus-induced gene silencing (VIGS) analysis. Conclusions A comprehensive analysis of NAT gene family was performed in cotton, including phylogenetic analysis, chromosomal location, collinearity analysis, motifs, gene structure and so on. Our results will further broaden the insight into the evolution and potential functions of NAT genes in cotton. Current findings could make significant contribution towards screening more candidate genes related to biotic and abiotic resistance for the improvement in cotton.

Abstract Background The internal NAD(P)H dehydrogenase (NDA) gene family was a member of the NAD(P)H dehydrogenase (ND) gene family, mainly involved in the non-phosphorylated respiratory pathways in mitochondria and played crucial roles in response to abiotic stress. Methods The whole genome identification, structure analysis and expression pattern of NDA gene family were conducted to analyze the NDA gene family. Results There were 51, 52, 26, and 24 NDA genes identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. According to the structural characteristics of genes and traits of phylogenetic tree, we divided the NDA gene family into 8 clades. Gene structure analysis showed that the NDA gene family was relatively conservative. The four Gossypium species had good collinearity, and segmental duplication played an important role in the evolution of the NDA gene family. Analysis of cis-elements showed that most GhNDA genes contained cis-elements related to light response and plant hormones (ABA, MeJA and GA). The analysis of the expression patterns of GhNDA genes under different alkaline stress showed that GhNDA genes were actively involved in the response to alkaline stress, possibly through different molecular mechanisms. By analyzing the existing RNA-Seq data after alkaline stress, it was found that an NDA family gene GhNDA32 was expressed, and then theGhNDA32 was silenced by virus-induced gene silencing (VIGS). By observing the phenotype, we found that the wilting degree of silenced plants was much higher than that of the control plant after alkaline treatment, suggesting that GhNDA32 gene was involved in the response to alkaline stress. Conclusions In this study, GhNDAs participated in response to alkaline stress, especially NaHCO3 stress. It was of great significance for the future research on the molecular mechanism of NDA gene family in responding to abiotic stresses.

Abstract Background Melatonin 2-hydroxylase (M2H) is the first enzyme in the catabolism pathway of melatonin, which catalyzes the production of 2-hydroxymelatonin (2-OHM) from melatonin. The content of 2-hydroxymelatonin in plants is much higher than that of melatonin. So M2H may be a key enzyme in the metabolic pathway of melatonin. Method We conducted a systematic analysis of the M2H gene family in Gossypium hirsutum based on the whole genome sequence by integrating the structural characteristics, phylogenetic relationships, expression profile, and biological stress of the members of the Gossypium hirsutum M2H gene family. Result We identified 265 M2H genes in the whole genome of Gossypium hirsutum, which were divided into 7 clades (clades I-VII) according to phylogenetic analysis. Most M2H members in each group had similar motif composition and gene structure characteristics. More than half of GhM2H members contain ABA-responsive elements and MeJA-responsive elements. Under different stress conditions, the expression levels of the gene changed, indicating that GhM2H members were involved in the regulation of abiotic stress. Some genes in the GhM2H family were involved in regulating melatonin levels in cotton under salt stress, and some genes were regulated by exogenous melatonin. Conclusion This study is helpful to explore the function of GhM2H, the downstream metabolism gene of melatonin in cotton, and lay the foundation for better exploring the molecular mechanism of melatonin improving cotton's response to abiotic stress.

Drought stress is one of the abiotic stresses that limits crop production and greatly affects crop yield. Enhancement of plant stress resistance by NaCl pretreatment has been reported, but the mechanism by which NaCl pretreatment activates cotton stress resistance remains unclear. In this study, upland cotton (Gossypium hirsutum cv H177) was used as the material to conducted the treatments with three replications: 0 Mm NaCl + 0% PEG6000 (Polyethylene glycol), 0 mM NaCl + 15% PEG6000, 50 mM NaCl + 15% PEG6000 to explore the molecular mechanism by which NaCl improves the drought tolerance of cotton. The results showed that pretreatment with 50 mM NaCl could alleviate the adverse effects of PEG on cotton seeds while promoting the elongation of root length. RNA-seq showed that NaCl specifically induced the expression of carotenoid-related genes. By silencing the upstream gene GHLUT2 of lutein synthesis, it was found that the chlorophyll of silenced plants decreased, and leaf wilting was more sensitive to drought. We found that NaCl enhanced the drought resistance of cotton by regulating genes related to the carotenoid and abscisic acid downstream synthesis pathways. This study provides a new reference for the study of drought resistance in cotton and a theoretical basis for the molecular breeding of cotton.

Abstract ZAT (Zinc Finger of Arabidopsis thaliana) proteins are composed of a plant-specific transcription factor family, which play an important role in plant growth, development, and stress resistance. To study the potential function of ZAT family in cotton, the whole genome identification, expression, and structure analysis of ZAT gene family were carried out. In this study, our analysis revealed the presence of 115, 56, 59, and 115 ZAT genes in Gossypium hirsutum, G. raimondii, G. arboreum and G. barbadense, respectively. According to the number of domains and phylogenetic characteristics, we divided ZAT genes of four Gossypium species into 4 different clades, and further divided them into 11 subfamilies. The results of collinearity analysis showed that segmental duplication was the main method to amplify the cotton ZAT genes family. Analysis of cis-elements of promoters indicated that most GhZAT genes contained cis-elements related to plant hormones and abiotic stress. According to heatmap analysis, the expression patterns of GhZAT genes under different stresses indicated that GhZAT genes were significantly involved in the response to cold, heat, salt, and PEG stress, possibly through different mechanisms. Among the highly expressed genes, we cloned a G. hirsutum gene GhZAT67. Through virus-induced gene silencing (VIGS), we found that its expression level decreased significantly after being silenced. Under alkaline treatment, the wilting degree of silenced plants was even greater than the wild type, which proved that GhZAT67 gene was involved in the response to alkaline stress.

Abstract Background The abiotic stress such as soil salinization and heavy metal toxicity has posed a major threat to sustainable crop production worldwide. Previous studies revealed that halophytes were supposed to tolerate other stress including heavy metal toxicity. Though HMAD (heavy-metal-associated domain) was reported to play various important functions in Arabidopsis, little is known in Gossypium. Results A total of 169 G. hirsutum genes were identified belonging to the HMAD gene family with the number of amino acids ranged from 56 to 1011. Additionally, 84, 76 and 159 HMAD genes were identified in each G. arboreum, G. raimondii and G. barbadense, respectively. The phylogenetic tree analysis showed that the HMAD gene family were divided into five classes, and 87 orthologs of HMAD genes were identified in four Gossypium species, such as genes Gh_D08G1950 and Gh_A08G2387 of G. hirsutum are orthologs of the Gorai.004G210800.1 and Cotton_A_25987 gene in G. raimondii and G. arboreum, respectively. In addition, 15 genes were lost during evolution. Furthermore, conserved sequence analysis found the conserved catalytic center containing an anion binding (CXXC) box. The HMAD gene family showed a differential expression levels among different tissues and developmental stages in G. hirsutum with the different cis-elements for abiotic stress. Conclusions Current study provided important information about HMAD family genes under salt-stress in Gossypium genome, which would be useful to understand its putative functions in different species of cotton.

Glycosyltransferases mainly catalyse the glycosylation reaction in living organisms and widely exists in plants. UGTs have been identified from G. raimondii, G. arboreum and G. hirsutum. However, Genome-wide systematic analysis of UGTs superfamily have not been studied in G. barbadense. 752 UGTs were identified from four cotton species and grouped into 18 clades, of which R was newly discovered clades. Most UGTs were clustered at both ends of the chromosome and showed a heterogeneous distribution. UGT proteins were widely distributed in cells, with the highest distribution in chloroplasts. UGTs of the same clade shared similar intron/exon structural features. During evolution, the gene family has undergone strong selection for purification. UGTs were significantly enriched in “transcriptional activity (GO:0016758)” and “metabolic processes (GO:0008152)”. Genes from the same clade differed in function under various abiotic stresses. The analysis of cis-acting element and qRT–PCR may indicate that GHUGTs play important roles in plant growth, development and abiotic stress. We further found that GHUGT74-2 plays an important role under submergence. The study broadens the understanding of UGTs in terms of gene characteristics, evolutionary processes, and gene function in cotton and provides a new way to systematically and globally understand the structure–function relationship of multigene families in the evolutionary process.

Glutamate decarboxylase (GAD) mainly regulated the biosynthesis of γ-aminobutyric acid (GABA) and played an important role in plant growth and stress resistance. To explore the potential function of GAD in cotton growth, the genome-wide identification, structure, and expression analysis of GAD genes were performed in this study. There were 10, 9, 5, and 5 GAD genes identified in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. GAD was divided into four clades according to the protein motif composition, gene structure, and phylogenetic relationship. The segmental duplication was the main way of the GAD gene family evolution. Most GhGADs respond to abiotic stress. Clade Ⅲ GAD was induced by Cd2+ stress, especially GhGAD6, and silencing GhGAD6 would lead to more serious Cd2+ poisoning in cotton. The oxidative damage caused by Cd2+ stress was relieved by increasing the GABA content. It was speculated that the decreased expression of GhGAD6 reduced the content of GABA in vivo and caused the accumulation of ROS. This study will further expand our understanding of the relationship between the evolution and function of the GhGAD gene family and provide new genetic resources for cotton breeding under environmental stress and phytoremediation.

Abstract Background The Raffinose synthetase (RAFS) genes superfamily is critical for the synthesis of raffinose, which accumulates in plant leaves under abiotic stress. However, it remains unclear whether RAFS contributes to resistance to abiotic stress in plants, specifically in the Gossypium species. Results In this study, we identified 74 RAFS genes from G. hirsutum, G. barbadense, G. arboreum and G. raimondii by using a series of bioinformatic methods. Phylogenetic analysis showed that the RAFS gene family in the four Gossypium species could be divided into four major clades; the relatively uniform distribution of the gene number in each species ranged from 12 to 25 based on species ploidy, most likely resulting from an ancient whole-genome polyploidization. Gene motif analysis showed that the RAFS gene structure was relatively conservative. Promoter analysis for cis-regulatory elements showed that some RAFS genes might be regulated by gibberellins and abscisic acid, which might influence their expression levels. Moreover, we further examined the functions of RAFS under cold, heat, salt and drought stress conditions, based on the expression profile and co-expression network of RAFS genes in Gossypium species. Transcriptome analysis suggested that RAFS genes in clade III are highly expressed in organs such as seed, root, cotyledon, ovule and fiber, and under abiotic stress in particular, indicating the involvement of genes belonging to clade III in resistance to abiotic stress. Gene co-expressed network analysis showed that GhRFS2A-GhRFS6A, GhRFS6D, GhRFS7D and GhRFS8A-GhRFS11A were key genes, with high expression levels under salt, drought, cold and heat stress. Conclusion The findings may provide insights into the evolutionary relationships and expression patterns of RAFS genes in Gossypium species and a theoretical basis for the identification of stress resistance materials in cotton.

Abstract Salinity is an extensive and adverse environmental stress to crop plants across the globe, and a major abiotic constraint responsible for limited crop production threatening the crop security. Soil salinization is a widespread problem across the globe, threatening the crop production and food security. Salinity impairs plant growth and development via reduction in osmotic potential, cytotoxicity due to excessive uptake of ions such as sodium (Na+) and chloride (Cl−), and nutritional imbalance. Cotton, being the most cultivated crop on saline-alkaline soils, it is of great importance to elucidate the mechanisms involved in Na2SO4 tolerance which is still lacking in upland cotton. Zhong 9835, a Na2SO4 resistant cultivar was screened for transcriptomic studies through various levels of Na2SO4 treatments, which results into identification of 3329 differentially expressed genes (DEGs) in roots, stems and leave at 300 mM Na2SO4 stress till 12 h in compared to control. According to gene functional annotation analysis, genes involved in reactive oxygen species (ROS) scavenging system including osmotic stress and ion toxicity were significantly up-regulated, especially GST (glutathione transferase). In addition, analysis for sulfur metabolism, results in to identification of two rate limiting enzymes [APR (Gh_D05G1637) and OASTL (Gh_A13G0863)] during synthesis of GSH from SO4 2−. Furthermore, we also observed a crosstalk of the hormones and TFs (transcription factors) enriched in hormone signal transduction pathway. Genes related to IAA exceeds the rest of hormones followed by ubiquitin related genes which are greater than TFs. The analysis of the expression profiles of diverse tissues under Na2SO4 stress and identification of relevant key hub genes in a network crosstalk will provide a strong foundation and valuable clues for genetic improvements of cotton in response to various salt stresses.

Abstract Background Cotton (Gossypium hirsutum) is considered a fairly salt tolerant crop however, salinity can still cause significant economic losses by affecting the yield and deteriorating the fiber quality. We studied a salt-tolerant upland cotton cultivar under temporal salt stress to unfold the salt tolerance molecular mechanisms. Biochemical response to salt stress (400 mM) was measured at 0 h, 3 h, 12 h, 24 h and 48 h post stress intervals and single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) combined with the unique molecular identifiers approach was used to identify differentially expressed genes (DEG). Results Antioxidant enzymes including, catalase (CAT), peroxidase (POD), superoxide dismutase (SOD) were found significantly induced under temporal salt stress, suggesting that reactive oxygen species scavenging antioxidant machinery is an essential component of salt tolerance mechanism in cotton. We identified a wealth of novel transcripts based on the PacBio long reads sequencing approach. Prolonged salt stress duration induces high number of DEGs. Significant numbers of DEGs were found under key terms related to stress pathways such as “response to oxidative stress”, “response to salt stress”, “response to water deprivation”, “cation transport”, “metal ion transport”, “superoxide dismutase”, and “reductase”. Key DEGs related to hormone (abscisic acid, ethylene and jasmonic acid) biosynthesis, ion homeostasis (CBL-interacting serine/threonine-protein kinase genes, calcium-binding proteins, potassium transporter genes, potassium channel genes, sodium/hydrogen exchanger or antiporter genes), antioxidant activity (POD, SOD, CAT, glutathione reductase), transcription factors (myeloblastosis, WRKY, Apetala 2) and cell wall modification were found highly active in response to salt stress in cotton. Expression fold change of these DEGs showed both positive and negative responses, highlighting the complex nature of salt stress tolerance mechanisms in cotton. Conclusion Collectively, this study provides a good insight into the regulatory mechanism under salt stress in cotton and lays the foundation for further improvement of salt stress tolerance.

Abstract Background The U-box gene encodes a ubiquitin ligase that contain U-box domain. The plant U-box gene (PUB) plays an important role in the response to stresses, but few reports about PUBs in cotton were available. Therefore research on PUBs is of great importance and a necessity when studying the mechanisms of stress- tolerance in cotton. Results In this study, we identified 93, 96, 185 and 208 PUBs from four sequenced cotton species G. raimondii (D5), G. arboreum (A2), G. hirsutum (AD1) and G. barbadense (AD2), respectively. Prediction analysis of subcellular localization showed that the PUBs in cotton were widely localized in cells, but primarily in the nucleus. The PUBs in cotton were classified into six subfamilies (A-F) on the basis of phylogenetic analysis, which was testified by the analysis of conserved motifs and exon-intron structures. Chromosomal localization analysis showed that cotton PUBs were unevenly anchored on all chromosomes, varying from 1 to 14 per chromosome. Through multiple sequence alignment analysis, 3 tandem duplications and 28 segmental duplications in cotton genome D5, 2 tandem duplications and 25 segmental duplications in A2, and 143 homologous gene pairs in A2 and D5 were found; however no tandem duplications in A2 or D5 were found. Additionally, 105, 14 and 17 homologous gene pairs were found in the intra-subgenome of At and Dt, At sub-genome and Dt sub-genome of G. hirsutum, respectively. Functional analysis of GhPUB85A and GhPUB45D showed that these genes positively responded to abiotic stresses, but the expression patterns were different. In addition, although the expression levels of these two homologous genes were similar, their contributions were different when responding to stresses, specifically showing different responses to abiotic stresses and functional differences between the two subgenomes of G. hirsutum. Conclusions This study reported the genome-wide identification, structure, evolution and expression analysis of PUBs in cotton, and the results showed that the PUBs were highly conserved throughout the evolutionary history of cotton. All PUB genes were involved in the response to abiotic stresses (including salt, drought, hot and cold) to varying degrees.

Serotonin N-acetyltransferase (SNAT) is a key enzyme in the biosynthesis of melatonin, and plays an important role in the regulation of melatonin synthesis. The study of SNAT is of great significance to understand the function of melatonin. In this study, we analyzed the structural characteristics, phylogenetic relationship, gene structure, expression pattern, evolutionary relationship and stress response of the members of the SNAT gene family in upland cotton through bioinformatics. A putative Serotonin n-acetyltransferase gene GhSNAT3D was identified, and preliminarily function of GhSNAT3D was verified by virus-induced gene silencing. We identified a total of 52 SNAT genes in the whole genome of G. hirsutum, and part of the GhSNATs were regulated by exogenous melatonin. The content of melatonin, antioxidant enzyme activity and Ca2+ content of GhSNAT3D gene silenced plants decreased, and the salt tolerance of GhSNAT3D gene silenced plants was reduced. Exogenous melatonin supplementation restored the salt tolerance of GhSNAT3D gene silenced plants. GhSNAT3D may interact with GhSNAT25D and ASMT to regulate melatonin synthesis. This study provided an important basis for further study on the regulation of melatonin in cotton against abiotic stress.

Plant 3-ketoacyl-CoA synthase (KCS) gene family catalyzed a β ketoacyl-CoA synthase, which was the rate-limiting enzyme for the synthesis of very long chain fatty acids (VLCFAs). Gossypium barbadense was well-known not only for high-quality fiber, which was perceived as a cultivated species of Gossypium. In this study, a total of 131 KCS genes were identified in four cotton species, there were 38, 44, 26, 23 KCS genes in the G. barbadense, the G. hirsutum, the G. arboreum and G. raimondii, respectively. The gene structure and expression pattern were analyzed. GBKCS genes were divided into six subgroups, the chromosome distribution of members of the family were mapped. The prediction of cis-acting elements of the GBKCS gene promoters suggested that the GBKCS genes may be involved in hormone signaling, defense and the stress response. Collinearity analysis on the KCS genes of the four cotton species were formulated. Tandem duplication played an indispensable role in the evolution of the KCS gene family. Specific expression analysis of 20 GBKCS genes indicated that GBKCS gene were widely expressed in the first 25 days of fiber development. Among them, GBKCS3, GBKCS8, GBKCS20, GBKCS34 were expressed at a high level in the initial long-term level of the G. barbadense fiber. This study established a foundation to further understanding of the evolution of KCS genes and analyze the function of GBKCS genes.

Gossypium barbadense is a cultivated cotton not only known for producing superior fiber but also for its salt and alkaline resistance. Here, we used Whole Genome Bisulfite Sequencing (WGBS) technology to map the cytosine methylation of the whole genome of the G. barbadense hypocotyl at single base resolution. The methylation sequencing results showed that the mapping rates of the three samples were 75.32, 77.54, and 77.94%, respectively. In addition, the Bisulfite Sequence (BS) conversion rate was 99.78%. Approximately 71.03, 53.87, and 6.26% of the cytosine were methylated at CG, CHG, and CHH sequence contexts, respectively. A comprehensive analysis of DNA methylation and transcriptome data showed that the methylation level of the promoter region was a positive correlation in the CHH context. Saline-alkaline stress was related to the methylation changes of many genes, transcription factors (TFs) and transposable elements (TEs), respectively. We explored the regulatory mechanism of DNA methylation in response to salt and alkaline stress during cotton hypocotyl elongation. Our data shed light into the relationship of methylation regulation at the germination stage of G. barbadense hypocotyl cell elongation and salt-alkali treatment. The results of this research help understand the early growth regulation mechanism of G. barbadense in response to abiotic stress.

Abstract Background DNA methylation is an important epigenetic factor that maintains and regulates gene expression. The mode and level of DNA methylation depend on the roles of DNA methyltransferase and demethylase, while DNA demethylase plays a key role in the process of DNA demethylation. The results showed that the plant’s DNA demethylase all contained conserved DNA glycosidase domain. This study identified the cotton DNA demethylase gene family and analyzed it using bioinformatics methods to lay the foundation for further study of cotton demethylase gene function. Results This study used genomic information from diploid Gossypium raimondii JGI (D), Gossypium arboreum L. CRI (A), Gossypium hirsutum L. JGI (AD1) and Gossypium barbadebse L. NAU (AD2) to Arabidopsis thaliana. Using DNA demethylase genes sequence of Arabidopsis as reference, 25 DNA demethylase genes were identified in cotton by BLAST analysis. There are 4 genes in the genome D, 5 genes in the genome A, 10 genes in the genome AD1, and 6 genes in the genome AD2. The gene structure and evolution were analyzed by bioinformatics, and the expression patterns of DNA demethylase gene family in Gossypium hirsutum L. were analyzed. From the phylogenetic tree analysis, the DNA demethylase gene family of cotton can be divided into four subfamilies: REPRESSOR of SILENCING 1 (ROS1), DEMETER (DME), DEMETER-LIKE 2 (DML2), and DEMETER-LIKE3 (DML3). The sequence similarity of DNA demethylase genes in the same species was higher, and the genetic relationship was also relatively close. Analysis of the gene structure revealed that the DNA demethylase gene family members of the four subfamilies varied greatly. Among them, the number of introns of ROS1 and DME subfamily was larger, and the gene structure was more complex. For the analysis of the conserved domain, it was known that the DNA demethylase family gene member has an endonuclease III (ENDO3c) domain. Conclusion The genes of the DNA demethylase family are distributed differently in different cotton species, and the gene structure is very different. High expression of ROS1 genes in cotton were under abiotic stress. The expression levels of ROS1 genes were higher during the formation of cotton ovule. The transcription levels of ROS1 family genes were higher during cotton fiber development.

As one of the cash crops, cotton is facing the threat of abiotic stress during its growth and development. It has been reported that melatonin is involved in plant defense against salt stress, but whether melatonin can improve cotton salt tolerance and its molecular mechanism remain unclear. We investigated the role of melatonin in cotton salt tolerance by silencing melatonin synthesis gene and exogenous melatonin application in upland cotton. In this study, applicating of melatonin can improve salt tolerance of cotton seedlings. The content of endogenous melatonin was different in cotton varieties with different salt tolerance. The inhibition of melatonin biosynthesis related genes and endogenous melatonin content in cotton resulted in the decrease of antioxidant enzyme activity, Ca2+ content and salt tolerance of cotton. To explore the protective mechanism of exogenous melatonin against salt stress by RNA-seq analysis. Melatonin played an important role in the resistance of cotton to salt stress, improved the salt tolerance of cotton by regulating antioxidant enzymes, transcription factors, plant hormones, signal molecules and Ca2+ signal transduction. This study proposed a regulatory network for melatonin to regulate cotton’s response to salt stress, which provided a theoretical basis for improving cotton’s salt tolerance.

The CLAVATA3/endosperm surrounding region-related (CLE) small peptides are a group of C-terminally encoded and post-translationally modified signal molecules involved in regulating the growth and development of various plants. However, the function and evolution of these peptides have so far remained elusive in cotton. In this study, 55, 56, and 86 CLE genes were identified in the Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum genomes, respectively, and all members were divided into seven groups. These groups were distinctly different in their protein characteristics, gene structures, conserved motifs, and multiple sequence alignment. Whole genome or segmental duplications played a significant role in the expansion of the CLE family in cotton, and experienced purifying selection during the long evolutionary process in cotton. Cis-acting regulatory elements and transcript profiling revealed that the CLE genes of cotton exist in different tissues, developmental stages, and respond to abiotic stresses. Protein properties, structure prediction, protein interaction network prediction of GhCLE2, GhCLE33.2, and GhCLE28.1 peptides were, respectively, analyzed. In addition, the overexpression of GhCLE2, GhCLE33.2, or GhCLE28.1 in Arabidopsis, respectively, resulted in a distinctive shrub-like dwarf plant, slightly purple leaves, large rosettes with large malformed leaves, and lack of reproductive growth. This study provides important insights into the evolution of cotton CLEs and delineates the functional conservatism and divergence of CLE genes in the growth and development of cotton.

Abstract Background DNA methylation, with a cryptic role in genome stability, gene transcription and expression, is involved in the drought response process in plants, but the complex regulatory mechanism is still largely unknown. Results Here, we performed whole-genome bisulfite sequencing (WGBS) and identified long non-coding RNAs on cotton leaves under drought stress and re-watering treatments. We obtained 31,223 and 30,997 differentially methylated regions (representing 2.48% of the genome) after drought stress and re-watering treatments, respectively. Our data also showed that three sequence contexts, including mCpG, mCHG, mCHH, all presented a hyper-methylation pattern under drought stress and were nearly restored to normal levels after the re-watering treatment. Among all the methylation variations, asymmetric CHH methylation was the most consistent with external environments, suggesting that methylation/demethylation in a CHH context may constitute a novel epigenetic modification in response to drought stress. Combined with the targets of long non-coding RNAs, we found that long non-coding RNAs may mediate variations in methylation patterns by splicing into microRNAs. Furthermore, the many hormone-related genes with methylation variations suggested that plant hormones might be a potential mechanism in the drought response. Conclusions Future crop-improvement strategies may benefit by taking into account not only the DNA genetic variations in cotton varieties but also the epigenetic modifications of the genome.

Upland cotton (Gossypium hirsutum L.), an important source of natural fiber, can tolerate relatively high salinity and drought stresses. In the present study, a plasma membrane Na+/H+ antiporter gene, GhSOS1, was cloned from a salt-tolerant genotype of G. hirsutum, Zhong 9807. The expression level of GhSOS1 in cotton roots was significantly upregulated in the presence of high concentrations of NaCl (200 mM), while its transcript abundance was increased when exposed to low temperature and drought stresses. Localization analysis using onion epidermal cells showed that the GhSOS1 protein was localized to the plasma membrane. The overexpression of GhSOS1 in Arabidopsis enhanced tolerance to salt stress, as indicated by a lower MDA content and decreased Na+/K+ ratio in transgenic plants. Moreover, the transcript levels of stress-related genes were significantly higher in GhSOS1 overexpression lines than in wild-type plants under salt treatment. Hence, GhSOS1 may be a potential target gene for enhancing salt tolerance in transgenic plants.

Recent researches on long noncoding RNAs (lncRNAs) have expanded our horizon of gene regulation and the cellular complexity. However, the number, characteristics and expression patterns of lncRNAs remain poorly characterized and how these lncRNAs biogenesis are regulated in response to drought stress in cotton are still largely unclear. In the study, using a reproducibility-based RNA-sequencing and bioinformatics strategy to analyze the lncRNAs of 9 samples under three different environment stresses (control, drought stress and re-watering, three replications), we totally identified 10,820 lncRNAs of high-confidence through five strict steps filtration, of which 9,989 were lincRNAs, 153 were inronic lncRNAs, 678 were anti-sense lncRNAs. Coding function analysis showed 6,470 lncRNAs may have the ability to code proteins. Small RNAs precursor analysis revealed that 196 lncRNAs may be the precursors to small RNAs, most of which (35.7%, 70) were miRNAs. Expression patterns analysis showed that most of lncRNAs were expressed at a low level and most inronic lncRNAs (75.95%) had a consistent expression pattern with their adjacent protein-coding genes. Further analysis of transcriptome data uncovered that lncRNAs XLOC_063105 and XLOC_115463 probably function in regulating two adjacent coding genes CotAD_37096 and CotAD_12502, respectively. Investigations of the content of plant hormones and proteomics analysis under drought stress also complemented the prediction. We analyzed the characteristics and the expression patterns of lncRNAs under drought stress and re-watering treatment, and found lncRNAs may be likely to involve in regulating plant hormones pathway in response to drought stress.

Salinity stress is a major abiotic factor that affects crop output, and as a pioneer crop in saline and alkaline land, salt tolerance study of cotton is particularly important. In our experiment, four salt-tolerance varieties with different salt tolerance indexes including CRI35 (65.04%), Kanghuanwei164 (56.19%), Zhong9807 (55.20%) and CRI44 (50.50%), as well as four salt-sensitive cotton varieties including Hengmian3 (48.21%), GK50 (40.20%), Xinyan96-48 (34.90%), ZhongS9612 (24.80%) were used as the materials. These materials were divided into salt-tolerant group (ST) and salt-sensitive group (SS). Illumina Cotton SNP 70K Chip was used to detect SNP in different cotton varieties. SNPv (SNP variation of the same seedling pre- and after- salt stress) in different varieties were screened; polymorphic SNP and SNPr (SNP related to salt tolerance) were obtained. Annotation and analysis of these SNPs showed that (1) the induction efficiency of salinity stress on SNPv of cotton materials with different salt tolerance index was different, in which the induction efficiency on salt-sensitive materials was significantly higher than that on salt-tolerant materials. The induction of salt stress on SNPv was obviously biased. (2) SNPv induced by salt stress may be related to the methylation changes under salt stress. (3) SNPr may influence salt tolerance of plants by affecting the expression of salt-tolerance related genes.