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 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 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 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.

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.

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.

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 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.

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 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.

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 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.

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.

Abstract Background Members of the AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) family are involved in various plant biological processes via protein-DNA and protein-protein interaction. However, no the systematic identification and analysis of AHL gene family have been reported in cotton. Results To investigate the potential functions of AHLs in cotton, genome-wide identification, expressions and structure analysis of the AHL gene family were performed in this study. 48, 51 and 99 AHL genes were identified from the G.raimondii, G.arboreum and G.hirsutum genome, respectively. Phylogenetic analysis revealed that the AHLs in cotton evolved into 2 clades, Clade-A with 4–5 introns and Clade-B with intronless (excluding AHL20–2). Based on the composition of the AT-hook motif(s) and PPC/DUF 296 domain, AHL proteins were classified into three types (Type-I/−II/−III), with Type-I AHLs forming Clade-B, and the other two types together diversifying in Clade-A. The detection of synteny and collinearity showed that the AHLs expanded with the specific WGD in cotton, and the sequence structure of AHL20–2 showed the tendency of increasing intron in three different Gossypium spp. The ratios of non-synonymous (Ka) and synonymous (Ks) substitution rates of orthologous gene pairs revealed that the AHL genes of G.hirsutum had undergone through various selection pressures, purifying selection mainly in A-subgenome and positive selection mainly in D-subgenome. Examination of their expression patterns showed most of AHLs of Clade-B expressed predominantly in stem, while those of Clade-A in ovules, suggesting that the AHLs within each clade shared similar expression patterns with each other. qRT-PCR analysis further confirmed that some GhAHLs higher expression in stems and ovules. Conclusion In this study, 48, 51 and 99 AHL genes were identified from three cotton genomes respectively. AHLs in cotton were classified into two clades by phylogenetic relationship and three types based on the composition of motif and domain. The AHLs expanded with segmental duplication, not tandem duplication. The expression profiles of GhAHLs revealed abundant differences in expression levels in various tissues and at different stages of ovules development. Our study provided significant insights into the potential functions of AHLs in regulating the growth and development in cotton.

Drought-induced 19 (Di19) protein is a Cys2/His2 (C2H2) type zinc-finger protein, which plays a crucial role in plant development and in response to abiotic stress. This study systematically investigated the characteristics of the GhDi19 gene family, including the member number, gene structure, chromosomal distribution, promoter cis-elements, and expression profiles. Transcriptomic analysis indicated that some GhDi19s were up-regulated under heat and salt stress. Particularly, two nuclear localized proteins, GhDi19-3 and GhDi19-4, were identified as being in potential salt stress responsive roles. GhDi19-3 and GhDi19-4 decreased sensitivity under salt stress through virus-induced gene silencing (VIGS), and showed significantly lower levels of H2O2, malondialdehyde (MDA), and peroxidase (POD) as well as significantly increased superoxide dismutase (SOD) activity. This suggested that their abilities were improved to effectively reduce the reactive oxygen species (ROS) damage. Furthermore, certain calcium signaling and abscisic acid (ABA)-responsive gene expression levels showed up- and down-regulation changes in target gene-silenced plants, suggesting that GhDi19-3 and GhDi19-4 were involved in calcium signaling and ABA signaling pathways in response to salt stress. In conclusion, GhDi19-3 and GhDi19-4, two negative transcription factors, were found to be responsive to salt stress through calcium signaling and ABA signaling pathways.

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.

Cotton (Gossypium) seed fibers can be divided into lint (long) or fuzz (very short). Using fiberless (fuzzless-lintless) mutants, the lint initiation gene Li3 was identified by map-based cloning. The gene is an R2R3-MYB transcription factor located on chromosome D12 (GhMML3_D12). Sequence analysis revealed that li3 is a loss-of-function allele containing a retrotransposon insertion in the second exon that completely blocks the gene’s expression. The genetic loci n2 and n3 underlying the recessive fuzzless phenotype in Gossypium hirsutum were also mapped. The genomic location of n3 overlapped with that of the dominant fuzzless locus N1, and n3 appeared to be a loss-of-function allele caused by a single nucleotide polymorphism (SNP) mutation in the coding region of GhMML3_A12. The n2 allele was found to be co-located with li3 and originated from G. babardense. n2 and li3 are possibly the multiple alleles of the GhMML3_D12 gene. Genetic analysis showed that Li3 and N3 are a pair of homologs with additive effects for the initiation of fibers (fuzz or lint). In addition, the presence of another locus was speculated, and it appeared to show an inhibitory effect on the expression of GhMML3. These findings provide new information about the genetic factors affecting the initiation of fibers in cotton.

The open-bud (ob) mutants in cotton display abnormal flower buds with the stigma and upper anthers exposed before blooming. This characteristic is potentially useful for the efficient production of hybrid seeds. The recessive inheritance pattern of the ob phenotype in allotetraploid cotton is determined by duplicated recessive loci (ob1ob1ob2ob2). In this study, ob1, which is a MIXTA-like MYB gene on chromosome D13 (MML10_Dt), was identified by map-based cloning. In Gossypium barbadense (Gb) acc. 3–79, a single nucleotide polymorphism (SNP) (G/A) at the splice site of the first intron and an 8-bp deletion in the third exon of MML10_Dt were found, which are the causative mutations at the ob1 loci. A 1783-bp deletion that leads to the loss of the third exon and accounts for the causal variation at the ob2 loci was found in MML10_At of Gossypium hirsutum (Gh) acc. TM-1. The ob phenotype results from the combination of these two loss-of-function loci. Genotyping assays showed that the ob1 and ob2 loci appeared after the formation of allotetraploid cotton and were specific for Gb and Gh, respectively. All Gb lines and most Gh cultivars carry the single corresponding mutant alleles. Genome-wide transcriptome analysis showed that some of the MYB genes and genes related to cell wall biogenesis, trichome differentiation, cytokinin signal transduction, and cell division were repressed in the ob mutants, which may lead to suppression of petal growth. These findings should be of value for breeding superior ob lines in cotton.