Construction and analysis of cotton (Gossypium arboreumL.) drought-related cDNA library
© Li et al; licensee BioMed Central Ltd. 2009
Received: 07 January 2009
Accepted: 02 July 2009
Published: 02 July 2009
Drought is one of the most important environmental factors causing water stress for cotton, and it greatly limits cotton growth and crop productivity. So far only a few drought-tolerance genes have been functionally characterized in details, and most efforts on this topic have been made in model organisms. Therefore, to identify more drought-related genes in cotton plays a crucial role in elucidating the underlying mechanisms of drought tolerance as well as utilizing bioengineering techniques to improve the tolerance in this organism.
Here we constructed a subtractive drought-tolerance cDNA library using suppressive subtractive hybridization (SSH). Through differential screening and bioinformatics analysis, we identified 392 positive clones with differential expression, corresponding 265 unique genes. By BLAST search against Genbank, we found that more than half of these EST sequences were homologous to those previously known drought-related genes and that there were 57 sequences with unknown functions, suggesting that many more genes are involved in this complex trait. Moreover, using RT-PCR, we examined the expression of nine representative candidate genes and confirmed that their expression levels were increased at different levels under drought stress.
Our results show that drought tolerance is a complex trait in cotton, which involves the coordination of many genes and multiple metabolism pathways. The candidate EST sequences we identified here would facilitate further functional studies of drought-related genes and provide important insights into the molecular mechanisms of drought-stress tolerance and genetic breeding in cotton.
Drought stress is a crucial limiting factor for cotton production. Hence, enhancing drought tolerance has been one of the key issues in the practice of cotton planting. Breeding has been used to improve the drought tolerance of cotton, but so far the progress with this approach has been slow and limited . Genetic engineering is another approach that could be used. However, with this approach, information about genes involved in cotton drought stress is required in advance. For this purpose, up-regulating key genes under drought stress may enhance drought tolerance.
Many drought-related genes have been reported in other plants, including maize, rice, and Arabidopsis. These genes can be mainly classified into two groups. One group contains proteins whose function is directly involved in stress tolerance, such as the enzymes required for photosynthesis enzymes [2, 3], LEA proteins , mRNA binding proteins, protein enzymes, proline-rich proteins  and various proteases . Proteins encoded by the other group appear to play regulatory roles, such as transcription factor MYB , zinc finger proteins [8, 9], heat shock proteins (HSP) [10, 11] and so on. Importantly, the wide range of these drought-stressed genes suggests that the responses to drought stress are rather complicated in plants.
Several techniques can be used to identify the genes expressed in response to drought stress, including DDRT-PCR , cDNA-AFLP , and suppression subtractive hybridization (SSH) . Among these techniques, SSH appears to produce fewer false positives . Therefore, we used SSH to construct a subtractive cDNA library of drought-stressed cotton.
Construction of suppression-subtracted cDNA library
Differential screening of SSH library and analysis of SSH cDNA sequences
To eliminate the false positive clones and quantify the relative expression level of the cloned cDNAs more accurately, we performed further cDNA differential screening. The cDNA clones of differentially expressed genes were identified by successive screenings with the subtracted tester and driver as probes, respectively. Finally, we identified 392 positive clones whose tester expression level was significantly higher (> 3.0-fold) than that with driver.
Bioinformatics sequence analysis
We sequenced 392 differentially expressed clones from different screening. The EST cluster analysis indicated that these sequences represented 265 unique ESTs. All the unique ESTs were submitted to the EST database of GenBank http://www.ncbi.nlm.nih.gov/dbEST. These 265 unique ESTs included 41 contigs and 224 singletons. Based on homology search of BLASTX and BLASTN, among the 265 non-redundant sequences, 208 clones (78.5%) are homologous to known genes; and 57 clones (21.5%) are homologous to genes with unknown function or without matches in the database.
In this study, we constructed a subtractive cDNA library containing water-stress induced transcripts in cotton seedlings. Based on our bioinformatics analysis, this SSH cDNA library contained many genes related to drought-stress tolerance. After differential screening, we identified 392 differentially expressed genes. With the homology search against Genbank, 265 of the genes had significant homologs in other plants, of which many are related to the drought-stress genes previously reported in Arabidopsis and maize.
The RT-PCR analysis on the EST expression showed that nine ESTs were up-regulated at different levels in drought-stress cotton seedlings. Among them, photosystem I psaH protein, and H+-ATPase-related gene  belong to the genes in cellular structure of leaves. These genes absorb and utilize water through adjusting the photosynthesis process. Under drought stress, the plant may reduce evaporation by closing the holes on the leaf surface, thereby fully utilizing water. Therefore, these two genes are highly induced.
The expression of glyceraldehyde-3-phosphate dehydogenase , alcohol dehydrogenase and drought-induced cysteine protease [18, 19] were also obviously increased under drought stress. These three genes have been reported to be related to water stress in other plants. Our study confirms that these genes are also involved in the drought tolerance of cotton, suggesting that the response to drought is a very complex physiological and biochemical process and involves multiple metabolism pathways and many enzymes.
Moreover, some genes induced under drought stress were found to be associated with other environmental stresses [20, 21], such as cold, salt, aluminium stress and so on. For example, our RT-PCR analysis confirmed one cold-related protein whose expression was slightly increased under drought stress. However, it should be noted that the number of such genes is very limited, which suggests that cotton is robust to drought stress and the damage is relatively moderate under drought stress.
Some studies show that under drought stress, the plant can improve the drought tolerance through adjusting osmoticum such as proline , trehalose  and glycinebetaine . In our study, the mRNA expression of proline-related drought-induced protein (DIP) was significantly increased under drought stress, suggesting that this gene plays a key role in the drought-stress defense. Heat-shock proteins (e.g., HSP 90) are known to be important in protecting plants against stresses. They can bind to denatured proteins and maintain the soluble state, then facilitate to establish proper protein configuration and prevent unwanted protein from aggregating with the existence of Mg2+ and ATP.
Compared with previously established SSH libraries in maize  and rice , our library contained many EST sequences annotated with drought tolerance. These ESTs appear to be involved in multiple metabolism pathways in the plant physiological and biochemical processes. In addition to known drought-induced genes, some differentially expressed genes are unknown, whose functional roles remain unclear and require further investigation in future.
Our cDNA collection contains a broad repertoire of drought-related genes encoding proteins involved in both initial and physiological responses during normal and water stress stages of cotton seedlings. Our study would contribute to a better understanding the molecular mechanism of water-stress tolerance and facilitate the genetic manipulation of stress tolerance in cotton.
Total RNA was obtained from the 3-to-6-leaf stage grown seedlings of cotton (Gossypium arboreum L.). Seeds were imbibed for one day in water at 30°C, and then sown in sterilized soil in plates for germination in Light-Emitted Feeding Box at 28°C. After 3–4 days, when appropriate, seedlings were transferred to soil and grown to the 3-to-6-leaf stage. In this study, 17% PEG6000 was used to induce "drought stress". Tissues were collected from seedlings maintained in water and in 17% PEG6000 at 1.5 hour after the start of the treatment and then hourly until 6.5 hour of the treatment. The tissues at different time points were pooled for RNA extraction.
Extraction and purification of RNA
A modified CTAB method was used for RNA extraction. The absorption ratio of OD260/280 was used to verify the quality of RNA. For PCR-select DNA subtraction, mRNA was purified with the Oligotex™ mRNA Mini Kit (Qiagen, Germany).
Construction of subtractive library
Suppression subtractive library was constructed according to the user manual of the PCR-Select™ cDNA Subtraction Kit (Clontech, USA). RNA from the seedlings maintained in water was used as the driver and RNA from the seedlings treated with 17% PEG 6000 was used as the tester. Double-stranded cDNA was prepared from mRNA. The cDNA was digested with RsaI for about 3–4 h and then ligated to adapters 1 and 2R provided in the kit. Two rounds of hybridization and PCR amplification were processed to normalize and enrich differentially expressed cDNA. The subtractive products were inserted into the pMD18-T Vector (Takara, China) and transformed into DH5α cells. The positive white clones were then selected and cultured on LB containing 100 μg/ml ampicillin in 96-well plates at 37°C for 7 h. Then, 30% glycerol was added, and the culture was kept at -80°C.
Amplification of cDNA inserts
The SSH products were amplified by PCR with nested primers provided in the kit. The reaction mixture contained 14.8 μl sterilized double distilled H2O (ddH2O), 2 μl 10 × ExTaq buffer (TaKaRa, China), 0.2 μl nested adapter 1, 2 R primers (10 μM each), 1.6 μl dNTP (0.5 mM), 0.2 μl ExTaq polymerase (TaKaRa, China), and 1 μl bacterial culture. PCR was performed as follows: 94°C for 5 min; 30 cycles at (95°C for 30 s; 68°C, 45 s; 72°C, 45 s; and 72°C, 5 min). The PCR products were electrophoresed on 1% agarose gel to confirm the amplification quality and quantity.
Differential screening of subtracted cDNA library
After the subtracted cDNA library is obtained, it is important to confirm that individual clones indeed represent differentially expressed genes. For this purpose, differential screening of the subtracted library helps to eliminate false positives and to save time and effort. Such a screening was performed with the Clontech™ PCR-Select Differential Screening Kit (Clontech, USA). The cDNA inserts of the positive clones were amplified by PCR. The amplified products (1 μl) were then spotted on Hybond-N nylon membrane (Amersham, UK). After air drying, the membranes were denatured with 0.6 N NaOH, neutralized with 0.5 M Tris-HCl (pH 7.5), rinsed with sterilized water for 30 s, and then baked for 2 h. Sterilized water was used as the negative control, and a housekeeping gene, Histone 3 was used as the positive control. Probes were prepared with DIG High Prime DNA Labeling and Detection Starter Kit (Roche, Switzerland), and the following probes were used: 1) forward subtracted, 2) reverse subtracted, 3) unsubtracted tester, and 4) unsubtracted driver. We selected the fragments that only hybridized with the labelled tester cDNA or with at least three-fold higher signals on these membranes compared with those hybridized with the labelled driver cDNA for further sequencing.
Bioinformatics sequence analysis
The positive clones identified by differential screening were sequenced at Beijing Genomics Institute. All the sequences were searched against the NCBI database with BLASTN and BLASTX http://www.ncbi.nlm.nih.gov/BLAST/. Functional classification of the ESTs was performed with the GO http://www.ebi.ac.uk/GO/ and COG http://www.ncbi.nlm.nih.gov/COG/ tools.
Some cDNAs were selected for RT-PCR to confirm the results of differential screening. The housekeeping gene, Histone 3, was used as an internal standard. Primers were designed with the Primer 5.0 software. PCR was conducted as follows: 95°C, 5 min; 30 cycles (94°C, 30 s; 54°C, 1 min; and 72°C, 1 min); 72°C, 10 min. To perform a quantitative analysis on RT-PCR bands, the Chemigenius2 Bio Imaging System (Syngene, USA) was used.
This study was supported by the National 863 Project of China (Grant no. 2006AA100105). We would like to thank Dr. Chee-Kok Chin for his valuable suggestions and encouragements during revision of this manuscript.
- He ZP, Zhou QQ, Zhu WJ: Higher lint percent and drought tolerant cotton line selected by radiation breeding. Shi Yan Sheng Wu Xue Bao. 2001, 34: 51-54.PubMedGoogle Scholar
- Norma ME, Eduardo CM, Anna T: Photosynthetic limitation by CO2 diffusion in drought stressed orange leaves on three rootstocks. Photosynth Res. 2008, 96: 163-172.View ArticleGoogle Scholar
- Kitao M, Lei TT: Circumvention of over-excitation of PSII by maintaining electron transport rate in leaves of four cotton genotypes developed under long-term drought. Plant Biol. 2007, 9: 69-76.View ArticlePubMedGoogle Scholar
- Xiao B, Huang Y, Tang N: Over-expression of a LEA gene in rice improves drought resistance under the field conditions. Theor Appl Genet. 2007, 115: 35-46.View ArticlePubMedGoogle Scholar
- Yu LX, Chamberland H, Tabaeizadeh Z: Negative regulation of gene expression of a novel proline-, threonine-, and glycine-rich protein by water stress in Lycopersicon chilense. Genome. 1996, 39: 1185-1193.View ArticlePubMedGoogle Scholar
- Zhang X, Liu S, Takano T: Two cysteine proteinase inhibitors from Arabidopsis thaliana, AtCYSa and AtCYSb, increasing the salt, drought, oxidation and cold tolerance. Plant Mol Biol. 2008, 68 (1-2): 131-143.View ArticlePubMedGoogle Scholar
- Cominelli E, Galbiati M, Tonelli C: A guard-cell-specific MYB transcription factor regulates stomatal movements and plant drought tolerance. Curr Biol. 2005, 15: 1196-1200.View ArticlePubMedGoogle Scholar
- Kim YO, Pan S, Kang H: A zinc finger-containing glycine-rich RNA-binding protein, atRZ-1a, has a negative impact on seed germination and seedling growth of Arabidopsis thaliana under salt or drought stress conditions. Plant Cell Physiol. 2007, 48: 1170-1181.View ArticlePubMedGoogle Scholar
- Yang X, Sun C, Lin Z: Molecular cloning and characterization of a gene encoding RING zinc finger ankyrin protein from drought-tolerant Artemisia desertorum. J Biosci. 2008, 33: 103-112.View ArticlePubMedGoogle Scholar
- Swindell WR, Huebner M, Weber AP: Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics. 2007, 8: 125-PubMed CentralView ArticlePubMedGoogle Scholar
- Ristic Z, Gifford DJ, Cass DD: Heat Shock Proteins in Two Lines of Zea mays L. That Differ in Drought and Heat Resistance. Plant Physiol. 1991, 97: 1430-1434.PubMed CentralView ArticlePubMedGoogle Scholar
- Maqbool A, Zahur M, Riazuddin S: Identification and expression of six drought responsive transcripts through differential display in desi cottion (Gossypium aroreum). Mol Biol (Mosk). 2008, 42: 559-565.Google Scholar
- Nguyen TT, Klueva N, Nguyen HT: Saturation mapping of QTL regions and identification of putative candidate genes for drought tolerance in rice. Mol Genet Genomics. 2004, 272: 35-46.View ArticlePubMedGoogle Scholar
- Diatchenko L, Lau YF, Siebert PD: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci USA. 1996, 93: 6025-6030.PubMed CentralView ArticlePubMedGoogle Scholar
- Diatchenko L, Lukyanov S, Siebert PD: Suppression subtractive hybridization: a versatile method for identifying differentially expressed genes. Methods Enzymol. 1999, 303: 349-380.View ArticlePubMedGoogle Scholar
- Hanitzsch M, Schnitzer D, Seidel T: Transcript level regulation of the vacuolar H (+)-ATPase subunit isoforms VHA-a, VHA-E and VHA-G in Arabidopsis thaliana. Mol Membr Biol. 2007, 24: 507-518.View ArticlePubMedGoogle Scholar
- Cherniad'ev II, Monakhova OF: Activity of NADP-dependent glyceraldehyde-phosphate dehydrogenase and phosphoenolpyruvate carboxylase in wheat leaves under water stress. Prikl Biokhim Mikrobiol. 2006, 42: 353-361.PubMedGoogle Scholar
- Koizumi M, Yamaguchi-Shinozaki K, Shinozaki K: Structure and expression of two genes that encode distinct drought-inducible cysteine proteinases in Arabidopsis thaliana. Gene. 1993, 129: 175-182.View ArticlePubMedGoogle Scholar
- Harrak H, Azelmat S, Tabaeizadeh Z: Isolation and characterization of a gene encoding a drought-induced cysteine protease in tomato (Lycopersicon esculentum). Genome. 2001, 44: 368-374.View ArticlePubMedGoogle Scholar
- Cuming AC, Cho SH, Quatrano RS: Microarray analysis of transcriptional responses to abscisic acid and osmotic, salt, and drought stress in the moss. New Phytol. 2007, 176: 275-287.View ArticlePubMedGoogle Scholar
- Vinocur B, Altman A: Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol. 2005, 16: 123-132.View ArticlePubMedGoogle Scholar
- He CY, Zhang JS, Chen SY: A soybean gene encoding a proline-rich protein is regulated by salicylic acid, an endogenous circadian rhythm and by various stresses. Theor Appl Genet. 2002, 104: 1125-1131.View ArticlePubMedGoogle Scholar
- Stiller I, Dulai S, Bánfalvi Z: Effects of drought on water content and photosynthetic parameters in potato plants expressing the trehalose-6-phosphate synthase gene of Saccharomyces cerevisiae. Planta. 2008, 227: 299-308.View ArticlePubMedGoogle Scholar
- Ben Hassine A, Ghanem ME, Lutts S: An inland and a coastal population of the Mediterranean xero-halophyte species Atriplex halimus L. differ in their ability to accumulate proline and glycinebetaine in response to salinity and water stress. J Exp Bot. 2008, 59: 1315-1326.View ArticlePubMedGoogle Scholar
- Zheng J, Liu YJ, Wang GY: Isolation and analysis of water stress induced genes in maize seedlings by subtractive PCR and cDNA macroarray. Plant Molecular Biology. 2004, 55: 807-823.View ArticlePubMedGoogle Scholar
- Wang Q, Guan Y, Chu C: Overexpression of a rice OsDREB1F gene increases salt, drought, and low temperature tolerance in both Arabidopsis and rice. Plant Molecular Biology. 2008, 67: 589-602.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.