Transcriptome Analysis Reveals Drought Stress response Mediated by Biological Processes and Key Pathways in Gossypium darwinii
Abstract Background: The allotetraploid wild cotton species, Gossypium darwinii, possess important traits of stress tolerance and good genetic stability that can be used for cotton improvement. In this study, we analyzed the RNA-seq transcriptomes from leaves of G. darwinii seedlings with and without drought stress. Results: A total of 86.7 million valid reads with an average length of 95.79 bp were generated from the two samples and 58,960 transcripts with a length of more than 500 bp were assembled. We searched the known proteins on the strength of sequence similarity; these transcripts were annotated with COG, KEGG, and GO functional categories. According to gene expression abundance RPKM value, we carried out RT-qPCR analysis to determine the expression pattern of the obtained transcription factors. A total of 58,961 genes was differentially expressed (DEG), with 32,693 and 25,919 genes found to be upregulated and downregulated, respectively. Through gene ontology and KEGG pathways, the upregulated genes were found to associate with all the GO terms, molecular functions (MF), biological process (BP) and cellular components (CC), which are highly linked to enhancing drought stress tolerance. The results obtained provide valuable information and a large number of transcripts, which can be exploited by cotton breeders in developing a more drought stress tolerant cotton germplasm.