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.

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.

Dwarf cottons are more resistant to damage from wind and rain and associated with stable, increased yields, and also desirable source for breeding the machine harvest varieties. In an effort to uncover the transcripts and miRNA networks involved in plant height, the transcriptome and small RNA sequencing were performed based on dwarf mutant Ari1327 (A1), tall-culm mutant Ari3697 (A3) and wild type Ari971 (A9) in Gossypium hirsutum. The stem apexes of wild-type upland cotton (Ari971) and its dwarf mutant (Ari1327) and tall-culm mutant (Ari3697) at the fifth true leaf stage were extracted for RNA, respectively. Transcriptome and small RNA libraries were constructed and subjected to next generation sequencing. The transcriptome sequencing analysis showed that the enriched pathways of top 3 differentially expressed genes (DEGs) were categorized as carotenoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction in both A1–A9 and A3–A9. The ABA and IAA related factors were differentially expressed in the mutants. Importantly, we found the lower expressed SAUR and elevated expressed GH3, and ABA related genes such as NCED and PP2C maybe relate to reduced growth of the plant height in Ari1327 which was consistent with the higher auxin and ABA content in this mutant. Furthermore, miRNA160 targeted to the auxin response factor (ARF) and miRNA166 (gma-miR166u and gma-miR166h-3p) targeted to ABA responsive element binding factor were related to the mutation in cotton. We have noticed that the cell growth related factors (smg7 targeted by gra-miR482 and 6 novel miRNAs and pectate-lyases targeted by osa-miR159f), the redox reactions related factors (Cytochrome P450 targeted by miR172) and MYB genes targeted by miR828, miR858 and miR159 were also involved in plant height of the cotton mutants. A total of 226 conserved miRNAs representing 32 known miRNA families were obtained, and 38 novel miRNAs corresponding to 23 unique RNA sequences were identified. Total 531 targets for 211 conserved miRNAs were obtained. Using PAREsnip, 27 and 29 miRNA/target conserved interactions were validated in A1–A9 and A3–A9, respectively. Furthermore, miRNA160, miRNA858 and miRNA172 were validated to be up-regulated in A1–A9 but down-regulated in A3–A9, whereas miRNA159 showed the opposite regulation. This comprehensive interaction of the transcriptome and miRNA at tall-culm and dwarf mutant led to the discovery of regulatory mechanisms in plant height. It also provides the basis for in depth analyses of dwarf mutant genes for further breeding of dwarf cotton.