Sugarcane (Saccharum spp.) is one of the most important crops in the world. Drought stress is a major abiotic stress factor that significantly reduces sugarcane yields. However the gene network that mediates plant responses to water stress remains largely unknown in several crop species. Although several microRNAs that mediate post-transcriptional regulation during water stress have been described in other species, the role of the sugarcane microRNAs during drought stress has not been studied. The objective of this work was to identify sugarcane miRNAs that are differentially expressed under drought stress and to correlate this expression with the behavior of two sugarcane cultivars with different drought tolerances. The sugarcane cultivars RB867515 (higher drought tolerance) and RB855536 (lower drought tolerance) were cultivated in a greenhouse for three months and then subjected to drought for 2, 4, 6 or 8 days. By deep sequencing of small RNAs, we were able to identify 18 miRNA families. Among all of the miRNAs thus identified, seven were differentially expressed during drought. Six of these miRNAs were differentially expressed at two days of stress, and five miRNAs were differentially expressed at four days. The expression levels of five miRNAs (ssp-miR164, ssp-miR394, ssp-miR397, ssp-miR399-seq 1 and miR528) were validated by RT-qPCR (quantitative reverse transcriptase PCR). Six precursors and the targets of the differentially expressed miRNA were predicted using an in silico approach and validated by RT-qPCR; many of these targets may play important roles in drought tolerance. These findings constitute a significant increase in the number of identified miRNAs in sugarcane and contribute to the elucidation of the complex regulatory network that is activated by drought stress.
Chapman EJ, Carrington JC (2007) Specialization and evolution of endogenous small RNA pathways. Nat Rev Gen 8: 884–896.
[5]
Lee Y, Kim M, Han JJ, Yeom KH, Lee S, et al. (2004) MicroRNA genes are transcribed by RNA polymerase II. Embo J 23: 4051–4060.
[6]
Lee RC, Feinbaum RL, Ambros V (1993) The C-Elegans Heterochronic Gene Lin-4 Encodes Small Rnas with Antisense Complementarity to Lin-14. Cell 75: 843–854.
[7]
Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP (2002) MicroRNAs in plants. Genes Dev 16: 1616–1626.
[8]
Sunkar R, Girke T, Jain PK, Zhu JK (2005) Cloning and characterization of microRNAs from rice. Plant Cell 17: 1397–1411.
[9]
Kantar M, Lucas SJ, Budak H (2011) miRNA expression patterns of Triticum dicoccoides in response to shock drought stress. Planta 233: 471–484.
[10]
Liu HH, Tian X, Li YJ, Wu CA, Zheng CC (2008) Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA 14: 836–843.
[11]
Sunkar R, Zhu JK (2004) Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell 16: 2001–2019.
[12]
Zhang ZX, Wei LY, Zou XL, Tao YS, Liu ZJ, et al. (2008) Submergence-responsive microRNAs are potentially involved in the regulation of morphological and metabolic adaptations in maize root cells. Ann Bot 102: 509–519.
[13]
Zhang XH, Zou Z, Gong PJ, Zhang JH, Ziaf K, et al. (2011) Over-expression of microRNA169 confers enhanced drought tolerance to tomato. Biotechnol Lett 33: 403–409.
[14]
Zhao BT, Liang RQ, Ge LF, Li W, Xiao HS, et al. (2007) Identification of drought-induced microRNAs in rice. Biochem Biophys Res Commun 354: 585–590.
[15]
Kulcheski FR, de Oliveira LF, Molina LG, Almerao MP, Rodrigues FA, et al. (2011) Identification of novel soybean microRNAs involved in abiotic and biotic stresses. BMC Genomics 12: 307.
[16]
Zhou L, Liu Y, Liu Z, Kong D, Duan M, et al. (2010) Genome-wide identification and analysis of drought-responsive microRNAs in Oryza sativa. J Exp Bot 61: 4157–4168.
[17]
Inman-Bamber NG, Smith DM (2005) Water relations in sugarcane and response to water deficits. Field Crop Res 92: 185–202.
[18]
Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, et al. (2007) High-Throughput Sequencing of Arabidopsis microRNAs: Evidence for Frequent Birth and Death of MIRNA Genes. Plos One 2: e219.
[19]
Song C, Wang C, Zhang C, Korir NK, Yu H, et al. (2010) Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics 11: 431.
[20]
Zhu QH, Spriggs A, Matthew L, Fan LJ, Kennedy G, et al. (2008) A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res 18: 1456–1465.
[21]
Yang XZ, Li L (2011) miRDeep-P: a computational tool for analyzing the microRNA transcriptome in plants. Bioinformatics 27: 2614–2615.
[22]
Unver T, Budak H (2009) Conserved microRNAs and their targets in model grass species Brachypodium distachyon. Planta 230: 659–669.
[23]
Rocha FR, Papini-Terzi FS, Nishiyama MY, Vencio RZN, Vicentini R, et al. (2007) Signal transduction-related responses to phytohormones and environmental challenges in sugarcane. BMC Genomics 8.
[24]
Islam MS, Wang MH (2010) Expression Patterns of an Abiotic Stress-inducible Dehydrin Gene, LeDhnT14, in Tomato. Hortic Environ Biotechnol 51: 556–561.
[25]
Ingram J, Bartels D (1996) The Molecular Basis of Dehydration Tolerance in Plants. Annu Rev Plant Physiol Plant Mol Biol 47: 377–403.
[26]
Vaseva I, Akiscan Y, Demirevska K, Anders I, Feller U (2011) Drought stress tolerance of red and white clover-comparative analysis of some chaperonins and dehydrins. Scientia Hortic 130: 653–659.
[27]
Giordani T, Buti M, Natali L, Pugliesi C, Cattonaro F, et al. (2011) An analysis of sequence variability in eight genes putatively involved in drought response in sunflower (Helianthus annuus L.). Theor Appl Genet 122: 1039–1049.
[28]
Wang L, Liu H, Li D, Chen H (2011) Identification and characterization of maize microRNAs involved in the very early stage of seed germination. BMC Genomics 12: 154.
[29]
Calvino M, Bruggmann R, Messing J (2011) Characterization of the small RNA component of the transcriptome from grain and sweet sorghum stems. BMC Genomics 12: 356.
[30]
Chen CF, Ridzon DA, Broomer AJ, Zhou ZH, Lee DH, et al. (2005) Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acid Res 33.
[31]
Varkonyi-Gasic E, Wu RM, Wood M, Walton EF, Hellens RP (2007) Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 3.
[32]
Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, et al. (2009) The Sorghum bicolor genome and the diversification of grasses. Nature 457: 551–556.
[33]
Jannoo N, Grivet L, Chantret N, Garsmeur O, Glaszmann JC, et al. (2007) Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome. Plant J 50: 574–585.
[34]
Dillon SL, Shapter FM, Henry RJ, Cordeiro G, Izquierdo L, et al. (2007) Domestication to crop improvement: genetic resources for Sorghum and Saccharum (Andropogoneae). Ann Bot 100: 975–989.
[35]
Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, et al. (2008) Criteria for Annotation of Plant MicroRNAs. Plant Cell 20: 3186–3190.
[36]
Kim VN (2005) MicroRNA biogenesis: Coordinated cropping and dicing. Nature Reviews Molecular Cell Biol 6: 376–385.
[37]
Unver T, Budak H (2009) Conserved microRNAs and their targets in model grass species Brachypodium distachyon. Planta 230: 659–669.
[38]
Zanca AS, Vicentini R, Ortiz-Morea FA, Del Bem LEV, da Silva MJ, et al. (2010) Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane. BMC Plant Biol 10.
[39]
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP (2006) A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 20: 3407–3425.
[40]
Martinez G, Forment J, Llave C, Pallas V, Gomez G (2011) High-throughput sequencing, characterization and detection of new and conserved cucumber miRNAs. Plos One 6: e19523.
[41]
Zhang BH, Pan XP, Stellwag EJ (2008) Identification of soybean microRNAs and their targets. Planta 229: 161–182.
[42]
Mi SJ, Cai T, Hu YG, Chen Y, Hodges E, et al. (2008) Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell 133: 116–127.
[43]
Voinnet O (2009) Origin, Biogenesis, and Activity of Plant MicroRNAs. Cell 136: 669–687.
[44]
Raman S, Greb T, Peaucelle A, Blein T, Laufs P, et al. (2008) Interplay of miR164, CUP-SHAPED COTYLEDON genes and LATERAL SUPPRESSOR controls axillary meristem formation in Arabidopsis thaliana. Plant J 55: 65–76.
[45]
Guo HS, Xie Q, Fei JF, Chua NH (2005) MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell 17: 1376–1386.
[46]
Ooka H, Satoh K, Doi K, Nagata T, Otomo Y, et al. (2003) Comprehensive Analysis of NAC Family Genes in Oryza sativa and Arabidopsis thaliana. DNA Res 10: 239–247.
[47]
Kikuchi K, Ueguchi-Tanaka M, Yoshida KT, Nagato Y, Matsusoka M, et al. (2000) Molecular analysis of the NAC gene family in rice. Mol Gen Genet 262: 1047–1051.
[48]
Collinge M, Boller T (2001) Differential induction of two potato genes, Stprx2 and StNAC, in response to infection by Phytophthora infestans and to wounding. Plant Mol Biol 46: 521–529.
[49]
Hegedus D, Yu M, Baldwin D, Gruber M, Sharpe A, et al. (2003) Molecular characterization of Brassicanapus NAC domain transcriptional activators induced in response to biotic and abiotic stress. Plant Mol Biol 53: 383–397.
[50]
Ditt RF, Gentile A, Tavares RG, Camargo SR, Fernandez JH, et al. (2011) Analysis of the stress-inducible transcription factor SsNAC23 in sugarcane plants. Sci Agric 68: 454–461.
[51]
Liu QL, Xu KD, Zhao LJ, Pan YZ, Jiang BB, et al. (2011) Overexpression of a novel chrysanthemum NAC transcription factor gene enhances salt tolerance in tobacco. Biotechnol Lett 33: 2073–2082.
[52]
Berger Y, Harpaz-Saad S, Brand A, Melnik H, Sirding N, et al. (2009) The NAC-domain transcription factor GOBLET specifies leaflet boundaries in compound tomato leaves. Development 136: 823–832.
[53]
Wu G, Lewis DR, Spalding EP (2007) Mutations in Arabidopsis multidrug resistance-like ABC transporters separate the roles of acropetal and basipetal auxin transport in lateral root development. Plant Cell 19: 1826–1837.
[54]
Sunkar R (2010) MicroRNAs with macro-effects on plant stress responses. Semin Cell Dev Biol 21: 805–811.
[55]
Dharmasiri N, Dharmasiri S, Estelle M (2005) The F-box protein TIR1 is an auxin receptor. Nature 435: 441–445.
[56]
Perez-Torres CA, Lopez-Bucio J, Cruz-Ramirez A, Ibarra-Laclette E, Dharmasiri S, et al. (2008) Phosphate availability alters lateral root development in Arabidopsis by modulating auxin sensitivity via a mechanism involving the TIR1 auxin receptor. Plant Cell 20: 3258–3272.
[57]
Jeong MJ, Park SC, Kwon HB, Byun MO (2000) Isolation and characterization of the gene encoding glyceraldehyde-3-phosphate dehydrogenase. Biochem Bioph Res Co 278: 192–196.
[58]
Arumugam Pillai MAP, Lihuang ZL, Akiyama TA (2002) Molecular cloning, characterization, expression and chromosomal location of OsGAPDH, a submergence responsive gene in rice (Oryza sativa L.). Theor Appl Genet 105: 34–42.
[59]
Abdel-Ghany SE, Pilon M (2008) MicroRNA-mediated systemic down-regulation of copper protein expression in response to low copper availability in Arabidopsis. J Biol Chem 283: 15932–15945.
[60]
Ranocha P, Chabannes M, Chamayou S, Danoun Sd, Jauneau A, et al. (2002) Laccase down-regulation causes alterations in phenolic metabolism and cell wall structure in poplar. Plant Physiol 129: 145–155.
[61]
Hoopes JT, Dean JFD (2004) Ferroxidase activity in a laccase-like multicopper oxidase from Liriodendron tulipifera. Plant Physiol Biochem 42: 27–33.
[62]
Neumann PM, Azaizeh H, Leon D (1994) Hardening of root cell walls: a growth inhibitory response to salinity stress. Plant Cell Environ 17: 303–309.
[63]
Cachorro P, Ortiz A, Barcelo AR, Cerda A (1993) Lignin deposition in vascular tissues of Phaseolus vulgaris roots in response to salt stress and Ca2+ ions. Phyton-Ann Rei Bot A 33: 33–40.
[64]
Liang M, Haroldsen V, Cai X, Wu Y (2006) Expression of a putative laccase gene, ZmLAC1, in maize primary roots under stress. Plant Cell Environ 29: 746–753.
[65]
Wei JZ, Tirajoh A, Effendy J, Plant AL (2000) Characterization of salt-induced changes in gene expression in tomato (Lycopersicon esculentum) roots and the role played by abscisic acid. Plant Sci 159: 135–148.
[66]
Wang T, Chen L, Zhao M, Tian Q, Zhang W-H (2011) Identification of drought-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. BMC Genomics 12: 367.
[67]
Gaxiola RA, Li J, Undurraga S, Dang LM, Allen GJ, et al. (2001) Drought- and salt-tolerant plants result from overexpression of the AVP1 H+-pump. Proc Natl Acad Sci U S A 98: 11444–11449.
[68]
Park S, Li JS, Pittman JK, Berkowitz GA, Yang HB, et al. (2005) Up-regulation of a H+-pyrophosphatase (H+-PPase) as a strategy to engineer drought-resistant crop plants. Proc Natl Acad Sci USA 102: 18830–18835.
[69]
Zhao FY, Zhang XJ, Li PH, Zhao YX, Zhang H (2006) Co-expression of the Suaeda salsa SsNHX1 and Arabidopsis AVP1 confer greater salt tolerance to transgenic rice than the single SsNHX1. Mol Breeding 17: 341–353.
[70]
Zhang H, Shen G, Kuppu S, Gaxiola R, Payton P (2011) Creating drought- and salt-tolerant cotton by overexpressing a vacuolar pyrophosphatase gene. Plant Signal Behav 6: 861–863.
[71]
Park S, Rancour DM, Bednarek SY (2007) Protein domain-domain interactions and requirements for the negative regulation of Arabidopsis CDC48/p97 by the plant ubiquitin regulatory X (UBX) domain-containing protein, PUX1. J Biol Chem 282: 5217–5224.
[72]
Rancour DM, Park S, Knight SD, Bednarek SY (2004) Plant UBX domain-containing protein 1, PUX1, regulates the oligomeric structure and activity of Arabidopsis CDC48. J Biol Chem 279: 54264–54274.
[73]
Golldack D, Luking I, Yang O (2011) Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep 30: 1383–1391.
[74]
Yang O, Popova OV, Suthoff U, Luking I, Dietz KJ, et al. (2009) The Arabidopsis basic leucine zipper transcription factor AtbZIP24 regulates complex transcriptional networks involved in abiotic stress resistance. Gene 436: 45–55.
[75]
Yoshida T, Fujita Y, Sayama H, Kidokoro S, Maruyama K, et al. (2010) AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. Plant J 61: 672–685.
[76]
Yang XZ, Zhang HY, Li L (2011) Global analysis of gene-level microRNA expression in Arabidopsis using deep sequencing data. Genomics 98: 40–46.
[77]
Hofacker IL (2003) Vienna RNA secondary structure server. Nucleic Acids Res 31: 3429–3431.
[78]
Audic S, Claverie JM (1997) The significance of digital gene expression profiles. Genome Res 7: 986–995.
[79]
Dai XB, Zhao PX (2011) psRNATarget: a plant small RNA target analysis server. Nucleic Acid Res 39: W155–W159.
[80]
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acid Res 25: 3389–3402.
[81]
Papini-Terzi FS, Rocha FR, Vencio RZN, Oliveira KC, Felix JD, et al. (2005) Transcription profiling of signal transduction-related genes in sugarcane tissues. DNA Res 12: 27–38.
[82]
Pabinger S, Thallinger GG, Snajder R, Eichhorn H, Rader R, et al. (2009) QPCR: Application for real-time PCR data management and analysis. BMC Bioinformatics 10: 268.
