Search

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.

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.

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.

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.

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.

Background Nuclear transcription factor Y subunit A (NFYA) plays an important role in plant growth, development, and response to abiotic stress. Results This study systematically analyzed the NFYA gene family. Chromosome location analysis found that some NFYA genes in Gossypium hirsutum may have been lost during evolution. Collinearity analysis and selection pressure analysis indicated that the GhNFYA gene family underwent fragment duplication and whole genome duplication during evolution. At the same time, promoter cis-element analysis and gene interaction network analysis predicted that the expression of GhNFYA gene may be regulated by plant hormones and stress. To further explore the function of the gene, Gossypium hirsutum seedlings were treated with 4 °C, 37 °C, salt and PEG stress, respectively, found that the expression of NFYA is stimulated by multiple environments. By constructing a co-expression network, interactions between genes were found to defend against salt stress. Through virus-induced gene silencing experiments, it was found that plants that silenced the GhNFYA16 gene were significantly more sensitive to salt stress. Conclusions This study found the relationship between the structure and function of NFYA gene family, provided a basis for the biological identification and functional verification of NFYA family members, and provided clues to clarify the specific roles of different types of NFYA proteins under different abiotic stress.

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.

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.

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.

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.

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.

In order to explore the molecular mechanism of salt tolerance and heterosis in cotton, a methylation-sensitive amplification polymorphism method based on capillary electrophoresis was used to analyze DNA methylation level in a cotton hybrid CCRI 29 and its two parents. The major results and conclusions are as follows: Firstly, salt tolerance test showed that CCRI 29 had higher salt-tolerance level than its both parents. The global DNA methylation level in CCRI 29 under salt treatment significantly increased, whereas the two parents did not change significantly between salt treatment and control. All kinds of variation of DNA methylation happened in cotton under salt treatment, and hypermethylation happened at a significantly higher rate than that of hypomethylation in CCRI 29 but not in its two parents. The results suggested that the increase of global DNA methylation level in cotton genome and also different methylation types played an important role in tolerance to salt treatment in cotton. Secondly, both hypomethylation and hypermethylation happened to different genes and some genes maintained the same methylation level after salt stress, which meant that complex gene expression alterations occurred when responding to salt stress in cotton, indicating the complicated characteristics of roles that specific genes played in salt tolerance. Thirdly, although most cytosine methylation sites in hybrid CCRI 29 shared the same status as that of at least one of the parents, the site number of hypomethylation is significantly higher than that of hypermethylation in CCRI 29 compared to parents under both control and salt stress, indicating that demethylation could be the mechanism to explain heterosis in cotton hybrid.

A methylation-sensitive amplification polymorphism method based on capillary electrophoresis was used to analyze DNA methylation levels in three cotton accessions, two salt-tolerant accessions CCRI 35 and Zhong 07 and one salt-sensitive accession CCRI 12. Many categories of DNA methylation happened in the three cotton accessions under salt treatment, including hypermethylation, hypomethylation, and other patterns. Hypermethylation happened at a significantly higher rate than that of hypomethylation in salt-tolerant accessions CCRI 35 and Zhong 07. On the contrary, in salt-sensitive accession CCRI 12, hypomethylation happened at a significantly higher rate than that of hypermethylation. In general, the global DNA methylation level significantly increased under salt stress in both salt-tolerant accessions CCRI 35 and Zhong 07, whereas there was no significant difference in the salt-sensitive accessions CCRI 12. Our results suggested that salt-tolerant cotton might have a mechanism of increasing the methylation level when responding to salt stress; the increase of the global level of DNA methylation and also different methylation patterns might play important roles in tolerance to salt stress in cotton. Some interesting genes were found through cloning and analysis of differently methylated DNA sequences, which might contribute to salt tolerance in cotton.

The CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 system has been widely used for genome editing in various plants because of its simplicity, high efficiency and design flexibility. However, to our knowledge, there is no report on the application of CRISPR/Cas9-mediated targeted mutagenesis in cotton. Here, we report the genome editing and targeted mutagenesis in upland cotton (Gossypium hirsutum L., hereafter cotton) using the CRISPR/Cas9 system. We designed two guide RNAs to target distinct sites of the cotton Cloroplastos alterados 1 (GhCLA1) and vacuolar H+-pyrophosphatase (GhVP) genes. Mutations in these two genes were detected in cotton protoplasts. Most of the mutations were nucleotide substitutions, with one nucleotide insertion and one substitution found in GhCLA1 and one deletion found in GhVP in cotton protoplasts. Subsequently, the two vectors were transformed into cotton shoot apexes through Agrobacterium-mediated transformation, resulting in efficient target gene editing. Most of the mutations were nucleotide deletions, and the mutation efficiencies were 47.6–81.8% in transgenic cotton plants. Evaluation using restriction-enzyme-PCR assay and sequence analysis detected no off-target mutations. Our results indicated that the CRISPR/Cas9 system was an efficient and specific tool for targeted mutagenesis of the cotton genome.

WRKY transcription factors play important roles in various stress responses in diverse plant species. In cotton, this family has not been well studied, especially in relation to fiber development. Here, the genomes and transcriptomes of Gossypium raimondii and Gossypium arboreum were investigated to identify fiber development related WRKY genes. This represents the first comprehensive comparative study of WRKY transcription factors in both diploid A and D cotton species. In total, 112 G. raimondii and 109 G. arboreum WRKY genes were identified. No significant gene structure or domain alterations were detected between the two species, but many SNPs distributed unequally in exon and intron regions. Physical mapping revealed that the WRKY genes in G. arboreum were not located in the corresponding chromosomes of G. raimondii, suggesting great chromosome rearrangement in the diploid cotton genomes. The cotton WRKY genes, especially subgroups I and II, have expanded through multiple whole genome duplications and tandem duplications compared with other plant species. Sequence comparison showed many functionally divergent sites between WRKY subgroups, while the genes within each group are under strong purifying selection. Transcriptome analysis suggested that many WRKY genes participate in specific fiber development processes such as fiber initiation, elongation and maturation with different expression patterns between species. Complex WRKY gene expression such as differential Dt and At allelic gene expression in G. hirsutum and alternative splicing events were also observed in both diploid and tetraploid cottons during fiber development process. In conclusion, this study provides important information on the evolution and function of WRKY gene family in cotton species.

The complex allotetraploid nature of the cotton genome (AADD; 2n = 52) makes genetic, genomic and functional analyses extremely challenging. Here we sequenced and assembled the Gossypium arboreum (AA; 2n = 26) genome, a putative contributor of the A subgenome. A total of 193.6 Gb of clean sequence covering the genome by 112.6-fold was obtained by paired-end sequencing. We further anchored and oriented 90.4% of the assembly on 13 pseudochromosomes and found that 68.5% of the genome is occupied by repetitive DNA sequences. We predicted 41,330 protein-coding genes in G. arboreum. Two whole-genome duplications were shared by G. arboreum and Gossypium raimondii before speciation. Insertions of long terminal repeats in the past 5 million years are responsible for the twofold difference in the sizes of these genomes. Comparative transcriptome studies showed the key role of the nucleotide binding site (NBS)-encoding gene family in resistance to Verticillium dahliae and the involvement of ethylene in the development of cotton fiber cells.

Background Trehalose (a-D-glucopyranosyl a-D-glucopyranoside) is a nonreducing disaccharide and is widely distributed in bacteria, fungi, algae, plants and invertebrates. In the study, the identification of trehalose-6-phosphate synthase (TPS) genes stress-related in cotton, and the genetic structure analysis and molecular evolution analysis of TPSs were conducted with bioinformatics methods, which could lay a foundation for further research of TPS functions in cotton. Results The genome information of Gossypium raimondii (group D), G. arboreum L. (group A), and G. hirsutum L. (group AD) was used in the study. Fifty-three TPSs were identified comprising 15 genes in group D, 14 in group A, and 24 in group AD. Bioinformatics methods were used to analyze the genetic structure and molecular evolution of TPSs. Real-time PCR analysis was performed to investigate the expression patterns of gene family members. All TPS family members in cotton can be divided into two subfamilies: Class I and Class II. The similarity of the TPS sequence is high within the same species and close within their family relatives. The genetic structures of two TPS subfamily members are different, with more introns and a more complicated gene structure in Class I. There is a TPS domain(Glyco transf_20) at the N-terminal in all TPS family members and a TPP domain(Trehalose_PPase) at the C-terminal in all except GrTPS6, GhTPS4, and GhTPS9. All Class II members contain a UDP-forming domain. The responses to environmental stresses showed that stresses could induce the expression of TPSs but the expression patterns vary with different stresses. Conclusions The distribution of TPSs varies with different species but is relatively uniform on chromosomes. Genetic structure varies with different gene members, and expression levels vary with different stresses and exhibit tissue specificity. The upregulated genes in upland cotton TM-1 is significantly more than that in G. raimondii and G. arboreum L. Shixiya 1. Electronic supplementary material The online version of this article (doi:10.1186/s12863-016-0360-y) contains supplementary material, which is available to authorized users.