Cotton (Gossypium spp.) is one of the major fibre crops of the world. Although it is classified as salt tolerant crop, cotton growth and productivity are adversely affected by high salinity, especially at germination and seedling stages. Identification of genes and miRNAs responsible for salt tolerance in upland cotton (Gossypium hirsutum L.) would help reveal the molecular mechanisms of salt tolerance. We performed physiological experiments and transcriptome sequencing (mRNA-seq and small RNA-seq) of cotton leaves under salt stress using Illumina sequencing technology. We investigated two distinct salt stress phases—dehydration (4 h) and ionic stress (osmotic restoration; 24 h)—that were identified by physiological changes of 14-day-old seedlings of two cotton genotypes, one salt tolerant and the other salt sensitive, during a 72-h NaCl exposure. A comparative transcriptomics was used to monitor gene and miRNA differential expression at two time points (4 and 24 h) in leaves of the two cotton genotypes under salinity conditions. The expression patterns of differentially co-expressed unigenes were divided into six groups using short time-servies expression miner software. During a 24-h salt exposure, 819 transcription factor unigenes were differentially expressed in both genotypes, with 129 unigenes specifically expressed in the salt-tolerant genotype. Under salt stress, 108 conserved miRNAs from known families were differentially expressed at two time points in the salt-tolerant genotype. We further analyzed the predicted target genes of these miRNAs along with the transcriptome for each time point. Important expressed genes encoding membrane receptors, transporters, and pathways involved in biosynthesis and signal transduction of calcium-dependent protein kinase, mitogen-activated protein kinase, and hormones (abscisic acid and ethylene) were up-regulated. We also analyzed the salt stress response of some key miRNAs and their target genes and found that the expressions of five of nine target genes exhibited significant inverse correlations with their corresponding miRNAs. On the basis of these results, we constructed molecular regulatory pathways and a potential regulatory network for these salt-responsive miRNAs. Our comprehensive transcriptome analysis has provided new insights into salt-stress response of upland cotton. The results should contribute to the development of genetically modified cotton with salt tolerance.

MYB family proteins are one of the most abundant transcription factors in the cotton plant and play diverse roles in cotton growth and evolution. Previously, few studies have been conducted in upland cotton, Gossypium hirsutum. The recent release of the G. hirsutum genome sequence provides a great opportunity to identify and characterize the entire upland cotton MYB protein family. In this study, we undertook a comprehensive genome-wide characterization and expression analysis of the MYB transcription factor family during cotton fiber development. A total of 524 non-redundant cotton MYB genes, among 1986 MYB and MYB-related putative proteins, were identified and classified into four subfamilies including 1R-MYB, 2R-MYB, 3R-MYB, and 4R-MYB. Based on phylogenetic tree analysis, MYB transcription factors were divided into 16 subgroups. The results showed that the majority (69.1 %) of GhMYBs genes belong to the 2R-MYB subfamily in upland cotton. Our comparative genomics analysis has provided novel insights into the roles of MYB transcription factors in cotton fiber development. These results provide the basis for a greater understanding of MYB regulatory networks and to develop new approaches to improve cotton fiber development.

Verticillium dahliae is a phytopathogenic fungal pathogen that causes vascular wilt diseases responsible for considerable decreases in cotton yields. The complex mechanism underlying cotton resistance to Verticillium wilt remains uncharacterized. Identifying an endogenous resistance gene may be useful for controlling this disease. We cloned the ribosomal protein L18 (GaRPL18) gene, which mediates resistance to Verticillium wilt, from a wilt-resistant cotton species (Gossypium arboreum). We then characterized the function of this gene in cotton and Arabidopsis thaliana plants. GaRPL18 encodes a 60S ribosomal protein subunit important for intracellular protein biosynthesis. However, previous studies revealed that some ribosomal proteins are also inhibitory toward oncogenesis and congenital diseases in humans and play a role in plant disease defense. Here, we observed that V. dahliae infections induce GaRPL18 expression. Furthermore, we determined that the GaRPL18 expression pattern is consistent with the disease resistance level of different cotton varieties. GaRPL18 expression is upregulated by salicylic acid (SA) treatments, suggesting the involvement of GaRPL18 in the SA signal transduction pathway. Virus-induced gene silencing technology was used to determine whether the GaRPL18 expression level influences cotton disease resistance. Wilt-resistant cotton species in which GaRPL18 was silenced became more susceptible to V. dahliae than the control plants because of a significant decrease in the abundance of immune-related molecules. We also transformed A. thaliana ecotype Columbia (Col-0) plants with GaRPL18 according to the floral dip method. The plants overexpressing GaRPL18 were more resistant to V. dahliae infections than the wild-type Col-0 plants. The enhanced resistance of transgenic A. thaliana plants to V. dahliae is likely mediated by the SA pathway. Our findings provide new insights into the role of GaRPL18, indicating that it plays a crucial role in resistance to cotton “cancer”, also known as Verticillium wilt, mainly regulated by an SA-related signaling pathway mechanism.

Salinity is a major abiotic stress that limits upland cotton growth and reduces fibre production worldwide. To reveal genetic regulation via transcript and protein levels after salt stress, we comprehensively analysed the global changes in mRNA, miRNA, and protein profiles in response to salt stress in two contrasting salt-tolerant cotton genotypes. In the current study, proteomic and mRNA-seq data were combined to reveal that some genes are differentially expressed at both the proteomic and mRNA levels. However, we observed no significant change in mRNA corresponding to most of the strongly differentially abundant proteins. This finding may have resulted from global changes in alternative splicing events and miRNA levels under salt stress conditions. Evidence was provided indicating that several salt stress-responsive proteins can alter miRNAs and modulate alternative splicing events in upland cotton. The results of the stringent screening of the mRNA-seq and proteomic data between the salt-tolerant and salt-sensitive genotypes identified 63 and 85 candidate genes/proteins related to salt tolerance after 4 and 24 h of salt stress, respectively, between the tolerant and sensitive genotype. Finally, we predicted an interaction network comprising 158 genes/proteins and then discovered that two main clusters in the network were composed of ATP synthase (CotAD_74681) and cytochrome oxidase (CotAD_46197) in mitochondria. The results revealed that mitochondria, as important organelles involved in energy metabolism, play an essential role in the synthesis of resistance proteins during the process of salt exposure. We provided a plausible schematic for the systematic salt tolerance model; this schematic reveals multiple levels of gene regulation in response to salt stress in cotton and provides a list of salt tolerance-related genes/proteins. The information here will facilitate candidate gene discovery and molecular marker development for salt tolerance breeding in cotton.

Background: C2H2-zinc finger protein family is commonly found in the plant, and it is known as the key actors in the regulation of transcription and vital component of chromatin structure. A large number of the C2H2-zinc finger gene members have not been well characterized based on their functions and structure in cotton. However, in other plants, only a few C2H2-zinc finger genes have been studied.Results: In this work, we performed a comprehensive analysis and identified 386, 196 and 195 C2H2-zinc finger genes inGossypium hirsutum(upland cotton),Gossypium arboreumandGossypium raimondii, respectively. Phylogenetic tree analysis of the C2H2-zinc finger proteins encoding the C2H2-zinc finger genes were classified into seven (7) subgroups. Moreover, the C2H2-zinc finger gene members were distributed in all cotton chromosomes though with asymmetrical distribution patterns. All the orthologous genes were detected between tetraploid and the diploid cotton, with 154 orthologous genes pair detected between upland cotton andGossypium arboreumwhile 165 orthologous genes were found between upland cotton andGossypium raimondii. Synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates (Ka/Ks) analysis indicated that the cotton C2H2-zinc finger genes were highly influenced mainly by negative selection, which maintained their protein levels after the duplication events. RNA-seq data and RT-qPCR validation of the RNA seq result revealed differential expression pattern of some the C2H2-zinc finger genes at different stages of cotton fiber development, an indication that the C2H2-zinc finger genes play an important role in initiating and regulating fiber development in cotton.Conclusions: This study provides a strong foundation for future practical genome research on C2H2-zinc finger genes in upland cotton. The expression levels of C2H2-zinc finger genes family is a pointer of their involvement in various biochemical and physiological functions which are directly related to cotton fiber development during initiation and elongation stages. This work not only provides a basis for determining the nominal role of the C2H2-zinc finger genes in fiber development but also provide valuable information for characterization of potential candidate genes involved in regulation of cotton fiber development.

Over 6% of agricultural land is affected by salinity. It is becoming obligatory to use saline soils, so growing salt-tolerant plants is a priority. To gain an understanding of the genetic basis of upland cotton tolerance to salinity at seedling stage, an intra-specific cross was developed from CCRI35, tolerant to salinity, as female with Nan Dan (NH), sensitive to salinity, as the male. A genetic map of 5178 SNP markers was developed from 277 F_2:3 populations. The map spanned 4768.098 cM, with an average distance of 0.92 cM. A total of 66 QTLs for 10 traits related to salinity were detected in three environments (0, 110, and 150 mM salt treatment). Only 14 QTLs were consistent, accounting for 2.72% to 9.87% of phenotypic variation. Parental contributions were found to be in the ratio of 3:1, 10 QTLs from the sensitive and four QTLs from the resistant parent. Five QTLs were located in A_t and nine QTLs in the D_t sub-genome. Moreover, eight clusters were identified, in which 12 putative key genes were found to be related to salinity. The GBS-SNPs-based genetic map developed is the first high-density genetic map that has the potential to provide deeper insights into upland cotton salinity tolerance. The 12 key genes found in this study could be used for QTL fine mapping and cloning for further studies.