A view of transcriptome during cold stress in sugarcane using Saccharum spontaneum genome


  • Shafee UR REHMAN Hazara University, Department of Biotechnology and Genetic Engineering, Mansehra, KP (PK)
  • Khushi MUHAMMAD Hazara University, Department of Biotechnology and Genetic Engineering, Mansehra, KP; Imperial College London, Department of Life Sciences, Sir Alex Fleming Building, South Kensington Campus, SW7 2AZ, London (PK)
  • Antonio C. MOTA PORTO Universidade Federal de Lavras, Department of Biology, Institute of Natural Science, MG (BR)
  • Youxiong QUE Fujian Agriculture and Forestry University Fuzhou, Key Lab of Sugarcane Biology and Genetic Breeding, Ministry of Agriculture, 350002-Fujian Province (CN)
  • Evandro NOVAES Universidade Federal de Lavras, Department of Biology, Institute of Natural Science, MG (BR)




bioinformatics, pipeline, RNA-Seq


Transcriptomic data of two sugarcane cultivars ‘ROC22’ and ‘GT08-1108’ were investigated for the expression analysis of cold responsive genes. The raw RNA Seq data of the sugarcane cultivars were downloaded from the SRA NCBI database and were reanalyzed and mapped by using Saccharum spontaneum genome.  In the Saccharum spontaneum reference genome, 83826 unigenes were annotated and, among these, 46,159 (55%) were functionally annotated with Gene Ontology (GO) categories. In the transcriptome-based analysis, 183,515 unigenes were assembled and, among these, 110,021 (60%) were functionally annotated with GO categories. For the cultivar GT08-1108, using the reference genome pipeline, 11,652 (13.9%) unigenes were differentially expressed (7,238 upregulated; 4,414 downregulated), while 16,145 (8.8%) were differentially expressed (8,965 upregulated; 7,180 downregulated) using transcriptome-based pipeline. In the cultivar ROC22, 11,516 (13.7%) genes were differentially expressed (7,174 upregulated; 4,342 downregulated) and 20,317 (11.1%) (10,898 upregulated; 9,419 downregulated) for the genome and transcriptome-based analysis, respectively. In the genome analyses, among downregulated genes, 3,248 were coincident between the two cultivars, the remaining 1,166 differentially expressed only in ‘GT-1180’ and 1,094 only in ‘ROC22’. With the transcriptome assembly, 13,113 genes were deferentially expressed in both cultivars, the remaining 3,032 unique to ‘GT08-1108’ and 7,204 in ‘ROC22’. We concluded that sugarcane in response to cold stress expresses many genes, although the transcriptome assembly overestimated the number of unigenes and, consequently, a higher number of differentially expressed genes. This may be due to difficulties in separating homeologues from paralogue genes. When a reference genome is available, we recommend its use since genes predicted on a reference genome tend to be more accurate.


Metrics Loading ...


Agarwal PK, Agarwal P, Reddy MK, Sopory SK (2006). Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Report 25:1263-1274. https://doi.org/10.1007/s00299-006-0204-8

Bolger AM, Lohse M, Usadel B (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30(15):2114-2120. https://doi.org/10.1093/bioinformatics/btu170

Cheavegatti-Gianotto A, De-Abreu HMC, Arruda P, Bespalhok Filho JC, Burnquist WL, Creste DE, Fatima Grossi-De-Sa SM (2011). Sugarcane (Saccharum X officinarum): a reference study for the regulation of genetically modified cultivars in Brazil. Tropical Plant Biology 4:62-89. https://doi.org/10.1007/s12042-011-9068-3

Chu Y, Corey DR (2012). RNA sequencing: platform selection, experimental design, and data interpretation. Nucleic Acid Therapeutics 22:271-274. https://doi.org/10.1089/nat.2012.0367

Conde A, Chaves MM, Geros H (2011). Membrane transport, sensing and signaling in plant adaptation to environmental stress. Plant and Cell Physiology 52:1583-1602. https://doi.org/10.1093/pcp/pcr107

Degand H, Faber AM, Dauchot N, Mingeot D, Watillon B, Cutsem PV, Morsomme P (2009). Proteomic analysis of chicory root identifies proteins typically involved in cold acclimation. Proteomics 9:2903-2907. https://doi.org/10.1002/pmic.200800744

Figueroa-Rodríguez KA, Hernández-Rosas F, Figueroa-Sandoval B, Velasco-Velasco J, Aguilar Rivera N (2019). What has been the focus of sugarcane research? A bibliometric overview. International Journal of Environmental Research and Public Health 16:3326. https://doi.org/10.3390/ijerph16183326

Fryer MJ, Oxborough K, Martin B, Ort DR, Baker NR (1995). Factors associated with depression of photosynthetic quantum efficiency in maize at low growth temperature. Plant Physiology 108:761-767. https://doi.org/10.1104/pp.108.2.761

Fowler S, Thomashow MF (2002). Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. The Plant Cell 14:1675-1690. https://doi.org/10.1105/tpc.003483

Gill SS, Tuteja N (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiology and Biochemistry 48:909-930. https://doi.org/10.1016/j.plaphy.2010.08.016

Grivet L, Arruda P (2002). Sugarcane genomics: depicting the complex genome of an important tropical crop. Current Opinion in Plant Biology 5:122-127. https://doi.org/10.1016/S1369-5266(02)00234-0

Gubis J, Vankova R, Cervena V, Dragunova M, Hudcovicova M, Lichtnerova H, Dokupil T Jurekova Z (2007). Transformed tobacco plants with increased tolerance to drought. South African Journal of Botany 73:505-511. https://doi.org/10.1016/j.sajb.2007.03.011

Hotta CT, Lembke CG, Domingues DS, Ochoa EA, Cruz GM, Melotto Passarin DM, … Crivellari AC (2010). The biotechnology roadmap for sugarcane improvement. Tropical Plant Biology 3:75-87. https://doi.org/10.1007/s12042-010-9050-5

Hoang NV, Furtado A, Mason PJ, Marquardt A, Kasirajan L, Thirugnanasambandam PP, ... Henry R J (2017). A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing. BMC Genomics 18(1):1-22. https://doi.org/10.1186/s12864-017-3757-8

Hoang NV, Furtado A, Mcqualter RB, Henry RJ (2015). Next generation sequencing of total DNA from sugarcane provides no evidence for chloroplast heteroplasmy. New Negatives in Plant Science 1(2):33-45. https://doi.org/10.1016/j.neps. 2015.10.001

Hussain A, Khan ZI, Ghafoor MY, Ashraf M (2004). Review sugarcane, sugar metabolism and some abiotic stresses. https://doi.org/10.3920/QAS2015.0824

James G (2008). Sugarcane.1st ed. John Wiely and Sons.

Khan MT, Khan IA, Yasmeen S, Nizamani GS, Afghan S (2019). Sugarcane Biofuels and Bioenergy Production in Pakistan. In: Current Scenario, Potential, and Future Avenues Sugarcane Biofuels. Springer, Cham pp 175-202.

Kollipara KP, Saab IN, Wych RD, Lauer MJ, Singletary GW (2002). Expression profiling of reciprocal maize hybrids divergent for cold germination and desiccation tolerance. Plant Physiology 129:974-992. https://doi.org/10.1104/pp.000729

Kovtun Y, Chiu WL, Tena G, Sheen J (2000). Functional analysis of oxidative stress-activated mitogen-activated protein kinase cascade in plants. Proceedings of the National Academy of Sciences 97:2940-2945. https://doi.org/10.1073/pnas.97.6.2940

Luo HL, Weili X Jianzhong L (2002). Study on the acceptation of Acacia mangium to low temperature stress. Journal of South China Agriculture University 23:51-53.

Lyzenga WJ, Stone SL (2012). Abiotic stress tolerance mediated by protein ubiquitination. Journal of Experimental Botany 63:599-616. https://doi.org/10.1093/jxb/err310

Lyzenga WJ, Stone SL (2012). Abiotic stress tolerance mediated by protein ubiquitination. Journal of Experimental Botany 63(2):599-616. https://doi.org/10.1093/jxb/err310

Machado DFSP, Ribeiro RV, Silveira JAG, da Filho JRM, Machado EC (2013). Rootstocks induce contrasting photosynthetic responses of orange plants to low night temperature without affecting the antioxidant metabolism. Theoretical and Experimental Plant Physiology 25:26-35.

Moore PH (1987). Breeding for stress resistance. In: Heinz DJ (Ed). Sugarcane Improvement through Breeding. pp 503-542.

Morin RM, Bainbridge A, Fejes M, Hirst M, Krzywinski M, Pugh T, McDonald H, Varhol R, Jones S, Marra M (2008). Profiling the HeLa S3 transcriptome using randomly primed cDNA and massively parallel short-read sequencing. Biotechniques 45(1):81-94. https://doi.org/10.2144/000112900

Moore JP, Nguema-Ona EE, Vicre-Gibouin M, Sorensen I, Willats WGT, Driouich A, Farrant JM (2013). Arabinose-rich polymers as an evolutionary strategy to plasticize resurrection plant cell walls against desiccation. Planta 237:739-754. https://doi.org/10.1007/s00425-012-1785-9

Manimekalai R, Suresh G, Singaravelu B (2022). Sugarcane transcriptomics in response to abiotic and biotic stresses: a review. Sugar Tech 15:1-24. https://doi.org/10.1007/s12355-021-01098-9

Nie GY, Robertson EJ, Fryer MJ, Leech RM, Baker NR (1995). Response of the photosynthetic apparatus in maize leaves grown at low temperature on transfer to normal growth temperature. Plant, Cell and Environment 18:1-12. https://doi.org/10.1111/j.1365-3040.1995.tb00538.x

Nogueira FT, De Rosa VE, Menossi M, Ulian EC, Arruda P (2003). RNA expression profiles and data mining of sugarcane response to low temperature. Plant Physiology 132:1811-1824. https://doi.org/10.1104/pp.102.017483

OECD (2016). Safety Assessment of Transgenic Organisms in the environment. Vol: 6. OECD Consensus Documents, Harmonisation of Regulatory Oversight in Biotechnology, OECD Publishing, Paris. https://www.oecd.org/chemicalsafety/biotrack/safety-assessment-of-transgenic-organisms-in-the-environment-volume-6-9789264253421-en.htm

Pang T, Ye CY, Xia X, Yin W (2013). De novo sequencing and transcriptome analysis of the desert shrub, Ammopiptanthus mongolicus, during cold acclimation using Illumina/Solexa. BMC Genomics 14:488. https://doi.org/10.1186/1471-2164-14-488

Quackenbush J (2002). Microarray data normalization and transformation. Nature Genetics 32:496-501. https://doi.org/10.1038/ng1032

Que Y, Su Y, Guo J, Wu Q, Xu L (2014). A global view of transcriptome dynamics during Sporisorium scitamineum Challenge in sugarcane by RNAseq. PLoS One Journal pone 9:106476. http://dx.doi.org/10.1186/1471-2164-15-996.

Rehman SU, Que Y, Khan S, Novaes E, Muhammad K (2020). Molecular study of CBF and WCOR14 genes in selected sugarcane cultivars and its wild relative Saccharum spontaneum L. Journal of Animal and Plant Sciences 30:913-922. https://doi.org/10.36899/JAPS.2020.4.0107

Reddy VS, Reddy AS (2004). Proteomics of calcium-signaling components in plants. Photochemistry 65:1745-1776. https://doi.org/10.1016/j.phytochem.2004.04.033

Robinson MD, Mccarthy DJ, Smyth GK (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139-40. https://doi.org/10.1093/bioinformatics/btp616

Robinson MD, Oshlack A (2010). A scaling normalization method for differential expression analysis of RNA-Seq data. Genome Biology 11:1-9. https://doi.org/10.1186/gb-2010-11-3-r25

Robinson MD, Smyth GK (2008). Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics 9:321-332. https://doi.org/10.1093/biostatistics/kxm030

Renaut J, Lutts S, Hofmann L, Hausman JF (2004). Responses of poplar to chilling temperatures: proteomic and physiological aspects. Plant Biology 7:81-90. https://doi.org/10.1055/s-2004-815733

Sabiha-Hanim S, Halim NAA (2018). Sugarcane bagasse pretreatment methods for ethanol production. In: Fuel Ethanol Production from Sugarcane. Intechopen. http://dx.doi.org/10.5772/intechopen.81656

Sales CRG, Ribeiro RV, Machado DFSP, Machado RS, Dovis VL, Lagôa AMMA (2012). Gas exchange and carbohydrate balance in sugarcane plants under root stressful conditions. Bragantia 71:319-327.

Serba DD, Meng X, Schnable J, Bashir E, Michaud JP, Vara Prasad PV, Perumal R (2021). Comparative transcriptome analysis reveals genetic mechanisms of sugarcane aphid resistance in grain Sorghum. International Journal of Molecular Sciences 22(13):7129. https://doi.org/10.3390/ijms22137129

Shi-Yun T, Li-Tao Y, Yang-Rui L (2018). Comparative analysis on transcriptome among different sugarcane cultivars under low temperature stress. Biotechnology Bulletin 34:116-124. http://dx.doi.org/10.13560/j.cnki.biotech.bull.1985.2018-0522

Souza GM, Berges H, Bocs S, Casu R, D’Hont A, Ferreira JE, Henry R, Ming R, Potier B, Sluys MV, Vincentz M, Paterson AH (2011). The sugarcane genome challenge: strategies for sequencing a highly complex genome. Tropical Plant Biology 4:145. https://doi.org/10.1007/s12042-011-9079-0

Strauss G, Hauser H (1986). Stabilization of lipid bilayer vesicles by sucrose during freezing. Proceedings of the National Academy of Sciences 83:2422-2426. https://doi.org/10.1073/pnas.83.8.2422

Sze H, Liang F, Hwang I, Curran AC, Harper JF (2000). Diversity and regulation of plant Ca2+ pumps: insights from expression in yeast. Annual Review of Plant Biology 51:433-462. https://doi.org/10.1146/annurev.arplant.51.1.433

Sato Y (2001). Heat shock-mediated APX gene expression and protection against chilling injury in rice seedlings. Journal of Experimental Botany 52:145-151. https://doi.org/10.1093/jexbot/52.354.145

Strauss G, Hauser H (1986). Stabilization of lipid bilayer vesicles by sucrose during freezing. Proceedings of the National Academy of Sciences 83(8):2422-2426. https://doi.org/10.1073/pnas.83.8.2422

Szekely G, Abraham E, Cseplo A, Rigo G, Zsigmond L, Csiszar J, … Koncz C (2008). Duplicated P5CS genes of Arabidopsis play distinct roles in stress regulation and developmental control of proline biosynthesis. Plant Journal 53:11-28.

Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, Xu W, Su Z (2017). Agrigo v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Research 45:122-129. https://doi.org/10.1093/nar/gkx382

Timperio AM, Egidi MG, Zolla L (2008). Proteomics applied on plant abiotic stresses: role of heat shock proteins (HSP). Journal of Proteomics 71:391-411. https://doi.org/10.1016/j.jprot.2008.07.005

Vettore AL, da Silva FR, Kemper EL, Souza GM, da Silva AM, Ferro MI, … Nobrega MP (2003). Analysis and functional annotation of an expressed sequence tag collection for tropical crop sugarcane. Genome Research 13(12):2725-2735. http://www.genome.org/cgi/doi/10.1101/gr.1532103

Wang Z, Gerstein M, Snyder M (2009). RNA-Seq: A Revolutionary Tool for Transcriptomics. Nature Review Genetics 10:57-63. https://doi.org/10.1038/nrg2484

Xiong L, Schumaker KS, Zhu JK (2002). Cell signaling during cold, drought, and salt stress. The Plant Cell 14:165-183. https://doi.org/10.1105/tpc.000596

Xu S, Wang J, Shang H, Huang Y, Yao W, Chen B, Zhang M (2018). Transcriptomic characterization and potential marker development of contrasting sugarcane cultivars. Scientific Reports 8(1):1-1. https://doi.org/10.1038/s41598-018-19832-x

Yan J, He C, Wang J, Mao Z, Holaday SA, Allen RD, Zhang H (2004). Overexpression of the Arabidopsis 14-3-3 protein GF14 lambda in cotton leads to a “stay-green” phenotype and improves stress tolerance under moderate drought conditions. Plant Cell Physiology 45:1007-1014.

You J, Chan Z (2015). ROS regulation during abiotic stress responses in crop plants. Frontiers in Plant Science 6:1092. https://doi.org/10.3389/fpls.2015.01092

Zhang J, Zhang X, Tang H, Zhang Q, Hua X, Ma X, ... Ming R (2018). Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nature Genetics 50(11):1565-1573. https://doi.org/10.1038/s41588-018-0237-2

Zhou D, Zhou X, Ling Y, Zhang Z, Su Z (2010). Agrigo: a GO analysis toolkit for the agricultural community. Nucleic Acids Research. https://doi.org/10.1093/nar/gkq31.



How to Cite

UR REHMAN, S., MUHAMMAD, K., MOTA PORTO, A. C., QUE, Y., & NOVAES, E. (2022). A view of transcriptome during cold stress in sugarcane using Saccharum spontaneum genome. Notulae Botanicae Horti Agrobotanici Cluj-Napoca, 50(3), 12765. https://doi.org/10.15835/nbha50312765



Research Articles
DOI: 10.15835/nbha50312765