Acessibilidade / Reportar erro

Expression analysis of transcription factors in sugarcane during cold stress

Análise de expressão de fatores de transcrição em cana-de-açúcar sob estresse por frio

Abstract

Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

Keywords:
environmental stress; differentially expressed genes; SRA NCBI; molecular reaction

Resumo

Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade de transcrição de muitos genes. Portanto, este estudo fornece base molecular para melhorar a tolerância ao frio em cana-de-açúcar e outras gramíneas economicamente importantes.

Palavras-chave:
estresse ambiental; genes diferencialmente expressos; SRA NCBI; reação molecular

1. Introduction

Saccharum hybrid (Sugarcane) belongs to the Poaceae family that stores a high concentration of sugar in its stem and it is cultivated in both hemispheres about 35° N and 35° S from the equator (Cheavegatti-Gianotto et al., 2011CHEAVEGATTI-GIANOTTO, A., ABREU, H.M.C., ARRUDA, P., BESPALHOK FILHO, J.C., BURNQUIST, W.L., CRESTE, S., CIERO, L., FERRO, J.A., FIGUEIRA, A.V.O., FILGUEIRAS, T.S., GROSSI-DE-SÁ, M.F., GUZZO, E.C., HOFFMANN, H.P., LANDELL, M.G.A., MACEDO, N., MATSUOKA, S., REINACH, F.C., ROMANO, E., SILVA, W.J., SILVA FILHO, M.C. and ULIAN, E.C., 2011. Sugarcane (Saccharum X officinarum): a reference study for the regulation of genetically modified cultivars in Brazil. Tropical Plant Biology, vol. 4, no. 1, pp. 62-89. http://dx.doi.org/10.1007/s12042-011-9068-3. PMid:21614128.
http://dx.doi.org/10.1007/s12042-011-906...
; OECD, 2016ORGANIZATION FOR ECONOMIC COOPERATION AND DEVELOPMENT – OECD, 2016.Harmonisation of regulatory oversight in biotechnology safety assessment of transgenic organisms in the environment, volume 6: OECD Consensus Documents. Paris: OECD Publishing.). It has been known to mankind since pre-history due to the production of sugar and now it is an important economic crop cultivated in several countries including Pakistan (James, 2008JAMES, G., 2008. Sugarcane. 1st ed. Hoboken: John Wiely and Sons.; Figueroa-Rodriguez et al., 2019FIGUEROA-RODRÍGUEZ, K.A., HERNÁNDEZ-ROSAS, F., FIGUEROA-SANDOVAL, B., VELASCO-VELASCO, J. and AGUILAR RIVERA, N., 2019. What has been the focus of sugarcane research? A bibliometric overview. International Journal of Environmental Research and Public Health, vol. 16, no. 18, pp. 3326. http://dx.doi.org/10.3390/ijerph16183326. PMid:31509963.
http://dx.doi.org/10.3390/ijerph16183326...
). Moreover, sugarcane is one of the most important commercial crops that gives essential by products such as sugar, ethanol, molasses, and bagasse (Sipaúba-Tavares et al., 2020SIPAÚBA-TAVARES, L.H., TEDESQUE, M.G. and SCARDOELI-TRUZZI, B., 2020. Evaluation of the effects of sugarcane molasses as a carbon source for Ankistrodesmus gracilis and Haematococcus pluvialis (Chlorophyceae). Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 80, no. 3, pp. 594-600. http://dx.doi.org/10.1590/1519-6984.216820. PMid:31644653.
http://dx.doi.org/10.1590/1519-6984.2168...
). According to the UNFAO, the world’s top 10 sugarcane producing countries includes Brazil, India followed by China, Thailand, Pakistan, Mexico, Colombia, Indonesia, Philippines, and United State (Sabiha-Hanim and Halim, 2018SABIHA-HANIM, S. and HALIM, N.A.A., 2018. Sugarcane bagasse pretreatment methods for ethanol production. In T.P. BASSO and L.C. BASSO, eds. Fuel ethanol production from sugarcane. London: INTECHOPEN, pp. 63-80.). Pakistan ranked fifth with average sugarcane production of 73.40 million tons in an area of 1.2 million hectares (Khan et al., 2019KHAN, M.T., KHAN, I.A., YASMEEN, S., NIZAMANI, G.S. and AFGHAN, S., 2019. Sugarcane Biofuels and Bioenergy Production in Pakistan: Current Scenario, Potential, and Future Avenues. In: M.T. KHAN and I.A. KHAN, eds. Sugarcane biofuels. Cham: Springer, pp. 175-202.).

During the process of development and growth, plants may expose to various external environmental challenges such as cold, drought, high salinity, and pathogens, which can adversely affect the plant yield, growth, and development (Agarwal et al., 2006AGARWAL, P.K., AGARWAL, P., REDDY, M.K. and SOPORY, S.K., 2006. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports, vol. 25, no. 12, pp. 1263-1274. http://dx.doi.org/10.1007/s00299-006-0204-8. PMid:16858552.
http://dx.doi.org/10.1007/s00299-006-020...
; Hussain et al., 2004HUSSAIN, A., KHAN, Z.I., GHAFOOR, M.Y., ASHRAF, M. and RASHID, M.H., 2004. Review sugarcane, sugar metabolism and some abiotic stresses. International Journal of Agriculture And Biology, vol. 6, no. 4, pp. 732-742., Srivastava and Kumar, 2020SRIVASTAVA, S. and KUMAR, P., 2020. Abiotic stress responses and tolerance mechanisms for sustaining crop productivity in sugarcane. In: M. HASANUZZAMAN, ed. Agronomic crops. Singapore: Springer, pp. 29-47. http://dx.doi.org/10.1007/978-981-15-0025-1_3). Cold stress or low temperature is key environmental stress that influences plant development and yield (Palva et al., 2002PALVA, E.T., TAEHTIHARJU, S., TAMMINEN, I., PUHAKAINEN, T., LAITINEN, R., SVENSSON, J. and HEINO, P., 2002. Biological mechanisms of low temperature STRESS Responsse: cold acclimation and development of freezing tolerance in plants. Ibaraki: JIRCAS, pp. 9-15.; Rehman et al., 2019REHMAN, S.U., MUHAMMAD, K., QUE, Y., REHMAN, A.U., NOVAES, E. and KHAN, S., 2019. presence of seventeen genes potentially involved in cold tolerance in sugarcane and Saccharum spontaneum genotypes. International Journal of Biosciences, vol. 14, no. 1, pp. 346-355. http://dx.doi.org/10.12692/ijb/14.1.346-355.
http://dx.doi.org/10.12692/ijb/14.1.346-...
; Rehman et al., 2020REHMAN, S.U., QUE, Y., KHAN, S., NOVAES, E. and 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, vol. 30, no. 4, pp. 913-922.). The cold stress comprises of freezing stress (below 0°C) and chilling stress (ranges from 0°C to 15°C) and it can restrict the geographical distribution and cultivation of crops (Tsvetanov et al., 2000TSVETANOV, S., ATANASSOV, A. and NAKAMURA, C., 2000. Gold responsive gene/protein families and cold/freezing tolerance in cereals. Biotechnology, Biotechnological Equipment, vol. 14, no. 1, pp. 3-11. http://dx.doi.org/10.1080/13102818.2000.10819054.
http://dx.doi.org/10.1080/13102818.2000....
). Various plants in the region of tropical and subtropical are chilling sensitive and are injured by nonfreezing low temperatures, while temperate plants are chilling tolerant as these exhibit complex mechanisms to acclimatize themselves (Pareek et al., 2017PAREEK, A., KHURANA, A., SHARMA, A.K. and KUMAR, R., 2017. An overview of signaling regulons during cold stress tolerance in plants. Current Genomics, vol. 18, no. 6, pp. 498-511. http://dx.doi.org/10.2174/1389202918666170228141345. PMid:29204079.
http://dx.doi.org/10.2174/13892029186661...
; Wani and Herath, 2018WANI, S.H. and HERATH, V., 2018. Cold tolerance in plants, physiological, molecular and genetic perspectives. 1st ed. Cham: Springer. http://dx.doi.org/10.1007/978-3-030-01415-5
http://dx.doi.org/10.1007/978-3-030-0141...
).

Sugarcane is a cold-sensitive plant and cold stress causes severe damage to this plant (Li et al., 2015LI, Y., WU, J., LI, X., ZHANG, R., LIU, X. and YANG, L., 2015. Damage in sugarcane production caused by long duration of chilling and frost in Guangxi China. International Journal of Agriculture Innovation and Research, vol. 3, no. 4, pp. 1139-1144.). According to Ebrahim et al (1998)EBRAHIM, M.K., ZINGSHEIM, O., EL-SHOURBAGY, M.N., MOORE, P.H. and KOMOR, E., 1998. Growth and sugar storage in sugarcane grown at temperatures below and above optimum. Journal of Plant Physiology, vol. 153, no. 5-6, pp. 593-602. http://dx.doi.org/10.1016/S0176-1617(98)80209-5.
http://dx.doi.org/10.1016/S0176-1617(98)...
, the chilling causes less growth and development of sugarcane, and thus sugarcane develops fewer shorter internodes and leaves, moreover, 37°C was considered the optimum temperature for better growth and development. According to field observation studies, sugarcane sensitivity to cold varies among varieties, moreover, subtropical hybrid plant species are reasonably cold tolerant in contrast to tropical species (Du et al., 1999DU, Y.C., NOSE, A. and WASANO, K., 1999. Thermal characteristics of C4 photosynthetic enzymes from leaves of three sugarcane species differing in cold sensitivity. Plant & Cell Physiology, vol. 40, no. 3, pp. 298-304. http://dx.doi.org/10.1093/oxfordjournals.pcp.a029541.
http://dx.doi.org/10.1093/oxfordjournals...
). Recently, sugarcane has attained the interest of researchers for improving its resilience against cold, especially the molecular mechanism in response to cold stress has become an important issue for its improvement (Yang et al., 2017YANG, Y., ZHANG, X., SU, Y., ZOU, J., WANG, Z., XU, L. and QUE, Y., 2017. miRNA alteration is an important mechanism in sugarcane response to low-temperature environment. BMC Genomics, vol. 18, no. 1, pp. 833. http://dx.doi.org/10.1186/s12864-017-4231-3. PMid:29084515.
http://dx.doi.org/10.1186/s12864-017-423...
; Su et al., 2020SU, W., REN, Y., WANG, D., SU, Y., FENG, J., ZHANG, C., TANG, H., XU, L., MUHAMMAD, K. and QUE, Y., 2020. The alcohol dehydrogenase gene family in sugarcane and its involvement in cold stress regulation. BMC Genomics, vol. 21, no. 1, pp. 521. http://dx.doi.org/10.1186/s12864-020-06929-9. PMid:32727370.
http://dx.doi.org/10.1186/s12864-020-069...
; Rehman et al., 2020REHMAN, S.U., QUE, Y., KHAN, S., NOVAES, E. and 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, vol. 30, no. 4, pp. 913-922.).

In plants, cryoprotective proteins encoded by various (COR) cold-responsive genes, protect plants cells from cold damage, which is one of the vital mechanisms of cold acclimation (Thomashow, 1999THOMASHOW, M.F., 1999. Plant cold acclimation: freezing tolerance genes and regulatory mechanism. Annual Review of Plant Physiology and Plant Molecular Biology, vol. 50, no. 1, pp. 571-599. http://dx.doi.org/10.1146/annurev.arplant.50.1.571. PMid:15012220.
http://dx.doi.org/10.1146/annurev.arplan...
). Presently, ICE1_CBF_COR transcriptional cascade is an obvious and well-known cold acclimation pathway, and in this pathway, C-repeat binding factors (CBF) also called dehydration responsive element binding factors (DREB) are induced rapidly during cold stress and bind to the cold-responsive genes promoter region and activate the transcription of downstream genes (Thomashow, 1999THOMASHOW, M.F., 1999. Plant cold acclimation: freezing tolerance genes and regulatory mechanism. Annual Review of Plant Physiology and Plant Molecular Biology, vol. 50, no. 1, pp. 571-599. http://dx.doi.org/10.1146/annurev.arplant.50.1.571. PMid:15012220.
http://dx.doi.org/10.1146/annurev.arplan...
; Chinnusamy et al., 2006CHINNUSAMY, V., ZHU, J. and ZHU, J.K., 2006. Gene regulation during cold acclimation in plants. Physiologia Plantarum, vol. 126, no. 1, pp. 52-61. http://dx.doi.org/10.1111/j.1399-3054.2006.00596.x.
http://dx.doi.org/10.1111/j.1399-3054.20...
). Low temperature inducible, CBF/DREB elements are the transcription factors binding to C-repeats binding site that get activated under drought and cold stresses are direct evidence of ABA-dependent pathway. Ppy CBFs (Pyrus pyrifolia CBF) studies of the Asian pear suggested that all PpyCBFs showed expression during salinity, drought, abscisic acid stress, and cold or low temperature. Moreover, PpyCBFs showed expression both in ABA-mediated and independent signaling pathways (Ahmad et al., 2019AHMAD, M., LI, J., YANG, Q., JAMIL, W., TENG, Y. and BAI, S., 2019. Phylogenetic, molecular, and functional characterization of PpyCBF proteins in Asian Pears (Pyrus pyrifolia). International Journal of Molecular Sciences, vol. 20, no. 9, pp. 2074. http://dx.doi.org/10.3390/ijms20092074. PMid:31035490.
http://dx.doi.org/10.3390/ijms20092074...
). This CBF/DRE transcription factors family includes CBF1, CBF2, and CBF3 that contain (AP2) DNA-binding domain (Akhtar et al., 2012AKHTAR, M., JAISWAL, A., TAJ, G., JAISWAL, J.P., QURESHI, M.I. and SINGH, N.K., 2012. DREB1/CBF transcription factors: their structure, function and role in abiotic stress tolerance in plants. Journal of genetics, vol. 91, no. 3, pp. 385-395. PMid:23271026.). The CBF transcription factors are recognized and bound to the C-repeat/dehydration response element motif (CRT/DRE) having conserved sequence (CCGAC) present in many cold-regulated genes promoters that mediate cold tolerance in plants (Shi et al., 2018SHI, Y., DING, Y. and YANG, S., 2018. Molecular regulation of CBF signaling in cold acclimation. Trends in Plant Science, vol. 23, no. 7, pp. 623-637. http://dx.doi.org/10.1016/j.tplants.2018.04.002. PMid:29735429.
http://dx.doi.org/10.1016/j.tplants.2018...
).

Recently, exploration of transcription factors involved in abiotic and biotic stresses to understand functional divergence in sugarcane has become the main interest of various researchers (Pereira-Santana et al., 2017PEREIRA-SANTANA, A., ALVARADO-ROBLEDO, E.J., ZAMORA-BRISEÑO, J.A., AYALA-SUMUANO, J.T., GONZALEZ-MENDOZA, V.M., ESPADAS-GIL, F., ALCARAZ, L.D., CASTAÑO, E., KEB-LLANES, M.A., SANCHEZ-TEYER, F. and RODRIGUEZ-ZAPATA, L.C., 2017. Transcriptional profiling of sugarcane leaves and roots under progressive osmotic stress reveals a regulated coordination of gene expression in a spatiotemporal manner. PLoS One, vol. 12, no. 12, pp. e0189271. http://dx.doi.org/10.1371/journal.pone.0189271. PMid:29228055.
http://dx.doi.org/10.1371/journal.pone.0...
; Wang et al., 2020WANG, D., WANG, L., SU, W., REN, Y., YOU, C., ZHANG, C., QUE, Y. and SU, Y., 2020. A class III WRKY transcription factor in sugarcane was involved in biotic and abiotic stress responses. Scientific Reports, vol. 10, no. 1, pp. 20964. http://dx.doi.org/10.1038/s41598-020-78007-9. PMid:33262418.
http://dx.doi.org/10.1038/s41598-020-780...
). Relatively, few cold-responsive genes have been reported and RNA-sequencing technology can provide useful transcriptional data for a deeper understanding of the molecular process in response to specific stress. Experiments performed on sugarcane (Saccharum spp) under cold stress have identified 34 EST and out of which 20 are novel COR genes and further one of them have shown similarity with CB (Miura and Furumoto, 2013MIURA, K. and FURUMOTO, T., 2013. Cold signaling and cold response in plants. International Journal of Molecular Sciences, vol. 14, no. 3, pp. 5312-5337. http://dx.doi.org/10.3390/ijms14035312. PMid:23466881.
http://dx.doi.org/10.3390/ijms14035312...
).

The growth period of sugarcane crop in subtropics regions is restricted to eight to nine months due to winter’s low temperature and it may usually occur from November to February every year (Dharshini et al., 2016DHARSHINI, S., CHAKRAVARTHI, M., NARAYAN, A.J., MANOJ, V.M., NAVEENARANI, M., KUMAR, R., MEENA, M., RAM, B. and APPUNU, C., 2016. De novo sequencing and transcriptome analysis of a low temperature tolerant Saccharum spontaneum clone IND 00-1037. Journal of Biotechnology, vol. 231, pp. 280-294. http://dx.doi.org/10.1016/j.jbiotec.2016.05.036. PMid:27269250.
http://dx.doi.org/10.1016/j.jbiotec.2016...
). In sugarcane breeding programs, the tolerance to cold and development of cold-tolerant cultivars is important and requires comprehensive knowledge on molecular mechanisms naturally adapted by tolerant genotypes during cold stress (Dharshini et al., 2016DHARSHINI, S., CHAKRAVARTHI, M., NARAYAN, A.J., MANOJ, V.M., NAVEENARANI, M., KUMAR, R., MEENA, M., RAM, B. and APPUNU, C., 2016. De novo sequencing and transcriptome analysis of a low temperature tolerant Saccharum spontaneum clone IND 00-1037. Journal of Biotechnology, vol. 231, pp. 280-294. http://dx.doi.org/10.1016/j.jbiotec.2016.05.036. PMid:27269250.
http://dx.doi.org/10.1016/j.jbiotec.2016...
). Very limited studies have been conducted in transcription factors activity during cold stress in sugarcane. Therefore, the present study aims to investigate the expression pattern of transcription factors in sugarcane during cold stress. Further, metabolic pathways were also identified and studied.

2. Methods and Materials

2.1. Experimental materials and cold stress treatment

Two sugarcane cultivars namely ROC22 and GT08-1108 were selected (Tang et al., 2018TANG, S., YANG, L. and LI, Y., 2018. Comparative analysis on transcriptome among different sugarcane cultivars under low temperature stress. Shengwu Jishu Tongbao, vol. 34, no. 12, pp. 116-124.). Cold stress experiment was performed at the Guangxi Sugarcane Research Institute, Guangxi Academy of Agricultural Sciences, China. RNA sequencing data (57.41 GB) of the selected sugarcane cultivars were downloaded from the SRA, NCBI database, where the samples were arranged with three biological replicates with a cold stress treatment time point (4°C) and control (28°C) (Table 1). Total RNA was extracted from the selected sugarcane cultivars (cold treated as well as normal sugarcane plants) by using an RNA extraction kit (Tiangen). The extracted RNA was further confirmed by agarose gel electrophoresis and for RNA degradation and contamination detection, Nanodrop was used, while Qubit Fluorometer was used for RNA concentration. Furthermore, Agilent 210 was used to detect RNA integrity. For mRNA enrichment, magnetic beads with oligo (dT) were used, and fragmentation solution was added to interrupt mRNA. cDNA with a six-base random primer was synthesized from the mRNA as a template. Buffer, dNTPS, and DNA polymerase were used to synthesize a double-stranded cDNA library. After cDNA library construction, the samples were sent for the Illumina High throughput sequencing to Beijing Nuohe Zhiyuan Biological Company.

Table 1
Samples of sugarcane analyzed.

2.2. Reading and splicing of transcriptome raw data

The transcriptomic raw data of the selected sugarcane cultivars were analyzed for quality control with FASTQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and low-quality sequences, as well as Illumina adaptors, were removed with Trimmomatic (Bolger et al., 2014BOLGER, A.M., LOHSE, M. and USADEL, B., 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, vol. 30, no. 15, pp. 2114-2120. http://dx.doi.org/10.1093/bioinformatics/btu170. PMid:24695404.
http://dx.doi.org/10.1093/bioinformatics...
). The sugarcane reference genome was used to align the reads by using STAR software and HTSeq-count was used for counting the reads aligned in each gene (Zhang et al., 2018ZHANG, J., ZHANG, X., TANG, H., ZHANG, Q., HUA, X., MA, X., ZHU, F., JONES, T., ZHU, X., BOWERS, J., WAI, C.M., ZHENG, C., SHI, Y., CHEN, S., XU, X., YUE, J., NELSON, D.R., HUANG, L., LI, Z., XU, H., ZHOU, D., WANG, Y., HU, W., LIN, J., DENG, Y., PANDEY, N., MANCINI, M., ZERPA, D., NGUYEN, J.K., WANG, L., YU, L., XIN, Y., GE, L., ARRO, J., HAN, J.O., CHAKRABARTY, S., PUSHKO, M., ZHANG, W., MA, Y., MA, P., LV, M., CHEN, F., ZHENG, G., XU, J., YANG, Z., DENG, F., CHEN, X., LIAO, Z., ZHANG, X., LIN, Z., LIN, H., YAN, H., KUANG, Z., ZHONG, W., LIANG, P., WANG, G., YUAN, Y., SHI, J., HOU, J., LIN, J., JIN, J., CAO, P., SHEN, Q., JIANG, Q., ZHOU, P., MA, Y., ZHANG, X., XU, R., LIU, J., ZHOU, Y., JIA, H., MA, Q., QI, R., ZHANG, Z., FANG, J., FANG, H., SONG, J., WANG, M., DONG, G., WANG, G., CHEN, Z., MA, T., LIU, H., DHUNGANA, S.R., HUSS, S.E., YANG, X., SHARMA, A., TRUJILLO, J.H., MARTINEZ, M.C., HUDSON, M., RIASCOS, J.J., SCHULER, M., CHEN, L.Q., BRAUN, D.M., LI, L., YU, Q., WANG, J., WANG, K., SCHATZ, M.C., HECKERMAN, D., VAN SLUYS, M.A., SOUZA, G.M., MOORE, P.H., SANKOFF, D., VANBUREN, R., PATERSON, A.H., NAGAI, C. and MING, R., 2018. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nature Genetics, vol. 50, no. 11, pp. 1565-1573. http://dx.doi.org/10.1038/s41588-018-0237-2. PMid:30297971.
http://dx.doi.org/10.1038/s41588-018-023...
). The number of reads from each library aligned in each gene was used as the measurement of gene expression levels.

2.3. Differential expression analysis and normalization of data and transcription factor

For differential expression analysis between treatments and normalization of RNA sequence data, R library edgeR was used (Robinson et al., 2010ROBINSON, M.D., MCCARTHY, D.J. and SMYTH, G.K., 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, vol. 26, no. 1, pp. 139-140. http://dx.doi.org/10.1093/bioinformatics/btp616. PMid:19910308.
http://dx.doi.org/10.1093/bioinformatics...
). Initially, the samples were filtered for normalization, and the expression of genes with 1 count per million (CPM) for at least two samples were excluded. TMM (Trimmed Mean) (Robinson and Oshlack, 2010ROBINSON, M.D. and OSHLACK, A., 2010. A scaling normalization method for differential expression analysis of RNA-Seq Data. Genome Biology, vol. 11, no. 3, pp. R25. http://dx.doi.org/10.1186/gb-2010-11-3-r25. PMid:20196867.
http://dx.doi.org/10.1186/gb-2010-11-3-r...
) method was used for estimating normalization factors of each sample, and to equalize and minimize the volume and high expression of genes and differences in sample size.

During differential expression analysis, the package edgeR assumes the negative binomial (NBD) distribution of the sequences count digital discrete data related to each sample. Therefore, Conditional Maximum Likelihood (CML) estimation was used for sequential analysis of gene to gene (tag-wise) dispersions, and Bayes empirical strategy assumed the shrinkage in tag-wise (gene-wise) distribution of residual variation. Differential gene expression analysis (DGE) and over dispersed data were estimated by Fisher Exact test (Robinson and Smyth, 2008ROBINSON, M.D. and SMYTH, G.K., 2008. Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics, vol. 9, no. 2, pp. 321-332. http://dx.doi.org/10.1093/biostatistics/kxm030. PMid:17728317.
http://dx.doi.org/10.1093/biostatistics/...
). The likelihood ratio test (LRT) was carried out for both low temperature and control experimental conditions., Moreover, False Discovery Rate (FDR) (P-value of) gave us the significance of the differential expression with a threshold of an FDR<0.05. Log fold change ratio (LogFC) measured the magnitude of differential expression (Quackenbush, 2002QUACKENBUSH, J., 2002. Microarray data normalization and transformation. Nature Genetics, suppl. 32, pp. 496-501. http://dx.doi.org/10.1038/ng1032. PMid:12454644.
http://dx.doi.org/10.1038/ng1032...
). The upregulated and downregulated genes under cold stress were defined as the genes expressed with LogFC ≥ 2 and LogFC value ≤-2 value respectively. Further, the DEGs were grouped with cluster analysis, annotated, and visualized within biochemical pathways.

2.4. Gene ontology and gene ontology enrichment analysis

Agrigo (Tian et al., 2017TIAN, T., LIU, Y., YAN, H., YOU, Q., YI, X., DU, Z., XU, W. and SU, Z., 2017. Agrigo v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Research, vol. 45, no. W1, pp. W122-W129. http://dx.doi.org/10.1093/nar/gkx382. PMid:28472432.
http://dx.doi.org/10.1093/nar/gkx382...
) gene enrichment analysis was performed for the differentially expressed transcription factors in sugarcane cultivars during cold stress. Besides, Agrigo is web-based online database and kit of tools to provide deep support to the agriculture community for gene enrichment analysis (Zhou et al., 2010ZHOU, D., ZHOU, X., LING, Y., ZHANG, Z. and SU, Z., 2010. Agrigo: a GO analysis toolkit for the agricultural community. Nucleic Acids Research, vol. 38, suppl. 2, pp. W64-W70. PMid:20435677.; Tian et al., 2017TIAN, T., LIU, Y., YAN, H., YOU, Q., YI, X., DU, Z., XU, W. and SU, Z., 2017. Agrigo v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Research, vol. 45, no. W1, pp. W122-W129. http://dx.doi.org/10.1093/nar/gkx382. PMid:28472432.
http://dx.doi.org/10.1093/nar/gkx382...
). The Gene Ontology and enrichment analysis were categorized into three that are: 1) Biological process, 2) Molecular function, 3) Cellular components.

3. Results

3.1. RNA sequence data analysis

The RNA sequence data of both cultivars (ROC22 and GT081108) were downloaded from SRA, NCBI database. The adapter and unwanted reads were removed, and the reads were mapped to a recently published reference genome (Zhang et al., 2018ZHANG, J., ZHANG, X., TANG, H., ZHANG, Q., HUA, X., MA, X., ZHU, F., JONES, T., ZHU, X., BOWERS, J., WAI, C.M., ZHENG, C., SHI, Y., CHEN, S., XU, X., YUE, J., NELSON, D.R., HUANG, L., LI, Z., XU, H., ZHOU, D., WANG, Y., HU, W., LIN, J., DENG, Y., PANDEY, N., MANCINI, M., ZERPA, D., NGUYEN, J.K., WANG, L., YU, L., XIN, Y., GE, L., ARRO, J., HAN, J.O., CHAKRABARTY, S., PUSHKO, M., ZHANG, W., MA, Y., MA, P., LV, M., CHEN, F., ZHENG, G., XU, J., YANG, Z., DENG, F., CHEN, X., LIAO, Z., ZHANG, X., LIN, Z., LIN, H., YAN, H., KUANG, Z., ZHONG, W., LIANG, P., WANG, G., YUAN, Y., SHI, J., HOU, J., LIN, J., JIN, J., CAO, P., SHEN, Q., JIANG, Q., ZHOU, P., MA, Y., ZHANG, X., XU, R., LIU, J., ZHOU, Y., JIA, H., MA, Q., QI, R., ZHANG, Z., FANG, J., FANG, H., SONG, J., WANG, M., DONG, G., WANG, G., CHEN, Z., MA, T., LIU, H., DHUNGANA, S.R., HUSS, S.E., YANG, X., SHARMA, A., TRUJILLO, J.H., MARTINEZ, M.C., HUDSON, M., RIASCOS, J.J., SCHULER, M., CHEN, L.Q., BRAUN, D.M., LI, L., YU, Q., WANG, J., WANG, K., SCHATZ, M.C., HECKERMAN, D., VAN SLUYS, M.A., SOUZA, G.M., MOORE, P.H., SANKOFF, D., VANBUREN, R., PATERSON, A.H., NAGAI, C. and MING, R., 2018. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nature Genetics, vol. 50, no. 11, pp. 1565-1573. http://dx.doi.org/10.1038/s41588-018-0237-2. PMid:30297971.
http://dx.doi.org/10.1038/s41588-018-023...
). The edgeR package was used to perform expression analysis. Our targeted genes (Transcription factor activity genes) were selected, and their expression analysis was studied. In our data analysis, the transcription factor was found significantly upregulated in both cultivars (ROC22 and GT08-1108). In ROC22, 963 transcription factor genes were differentially upregulated in 5558 genes, while 283 genes were downregulated in 3252 genes (Table 2). While in GT08-1108, 973 transcription factor genes were upregulated in 5649 genes, and 293 genes were down-regulated in 3289 genes (Table 2).

Table 2
Differentially expressed upregulated and downregulated transcription factors during cold stress in sugarcane cultivars ROC22 and GT08-1108.

3.2. Transcription factors upregulated in sugarcane cultivars during cold stress

The transcription factor activities were found upregulated in protein binding, transcription factor binding, the activity of DNA-sequence specific binding transcription factor, activity in transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, transcription factor binding of RNA polymerase II to DNA sequence-specific regulatory region, the activity of TF of RNA polymerase II, transcription factor cofactor activity, transcription factor activity at plastid promoter, transcription factor activity from RNA polymerase I promoter, transcription factor activity from the promoter of RNA polymerase II and RNA polymerase III, in both cultivars (Table 3 and 4). Mostly, the transcription factor activity was found upregulated in sequence-specific DNA binding up to 57% in total DGE (Figure 1).

Table 3
Transcription factor activity of upregulated genes during cold stress in ROC22 cultivar.
Table 4
Transcription factor activity of upregulated genes during cold stress in GT08-1108 cultivar.
Figure 1
Differentially upregulated transcription factors in selected sugarcane cultivars ROC22 and GT08-1108.

3.3. Transcription factors downregulated in sugarcane cultivars during cold stress

In the present data analysis, the transcription factor activities were found highly upregulated. Although transcription factor activity in our analysis was also found downregulated in both cultivars. Genes (TF activity) are differentially downregulated in protein binding, transcription binding, the activity of sequence-specific DNA binding transcription factor, transcription factor complex, transcription factor activity of RNA polymerase II, the activity of nucleic acid binding transcription factor (Table 5 and 6). In downregulated transcription factor (TF) activity, the activity of transcription factor in protein binding was found 31% among all downregulated genes (Figure 2).

Table 5
Transcription factor activity of downregulated genes during cold stress in sugarcane cultivar ROC22.
Table 6
Transcription factor activity of downregulated genes during cold stress in sugarcane cultivar GT08-1108.
Figure 2
Differential downregulated transcription factors in selected sugarcane cultivars ROC22 and GT08-1108.

3.4. Metabolic pathways of the transcription factor in sugarcane during cold stress

The gene ontology and gene enrichment analysis were performed using Agrigo gene enrichment analysis. The transcription activity was enriched in molecular function and cellular components while in biological processes no transcription factor activity was found. In both cultivars, the transcription factor activity was found highly upregulated in molecular function and cellular components (Figure 3 and 4). Further, the transcription factor activity was also found downregulated in molecular function and cellular components in both cultivars (Figure 5 and 6).

Figure 3
Schematic view of the metabolic pathway for transcription factors upregulated in sugarcane cultivar GT08-1108 during cold stress.
Figure 4
Schematic view of the metabolic pathway for transcription factors upregulated in sugarcane cultivar ROC22 during cold stress.
Figure 5
Schematic view of the metabolic pathway for transcription factors downregulated in sugarcane cultivar GT08-1108 during cold stress.
Figure 6
Schematic view of the metabolic pathway for transcription factors downregulated in sugarcane cultivar ROC22 during cold stress.

4. Discussion

Cold stress is one of the major abiotic stresses which reduces the quality, yield, and growth of many economically important crops worldwide (Zhou et al., 2011ZHOU, M.Q., SHEN, C., WU, L.H., TANG, K.X. and LIN, J., 2011. CBF-dependent signaling pathway: a key responder to low temperature stress in plants. Critical Reviews in Biotechnology, vol. 31, no. 2, pp. 186-92. http://dx.doi.org/10.3109/07388551.2010.505910. PMid:20919819.
http://dx.doi.org/10.3109/07388551.2010....
; Knight and Knight, 2012KNIGHT, M.R. and KNIGHT, H., 2012. Low-temperature perception leading to gene expression and cold tolerance in higher plants. The New Phytologist, vol. 195, no. 4, pp. 737-751. http://dx.doi.org/10.1111/j.1469-8137.2012.04239.x. PMid:22816520.
http://dx.doi.org/10.1111/j.1469-8137.20...
; Rehman et al., 2020REHMAN, S.U., QUE, Y., KHAN, S., NOVAES, E. and 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, vol. 30, no. 4, pp. 913-922.). During cold stress, plants adjust or adjust their biochemical and physiological activities by reprogramming the expression of many cold-responsive genes (Bohnert et al., 1995BOHNERT, H.J., NELSON, D.E. and JENSEN, R.G., 1995. Adaptation to environmental stresses. The Plant Cell, vol. 7, no. 7, pp. 1099-1111. http://dx.doi.org/10.2307/3870060. PMid:12242400.
http://dx.doi.org/10.2307/3870060...
; Chinnusamy et al., 2007CHINNUSAMY, V., ZHU, J. and ZHU, J.K., 2007. Cold stress regulation of gene expression in plants. Trends in Plant Science, vol. 12, no. 10, pp. 444-451. http://dx.doi.org/10.1016/j.tplants.2007.07.002. PMid:17855156.
http://dx.doi.org/10.1016/j.tplants.2007...
). In subtropical regions, chilling stress is one of the major constraints for sugarcane cultivation, which hampers yield and sugar production (Mathur et al., 2021MATHUR, S., SUNOJ, V.S.J., ELSHEERY, N.I., REDDY, V.R., JAJOO, A. and CAO, K.F., 2021. Regulation of Photosystem II heterogeneity and Photochemistry in two cultivars of C4 crop Sugarcane under Chilling stress. Frontiers in Plant Science, vol. 12, pp. 627012. http://dx.doi.org/10.3389/fpls.2021.627012. PMid:33643354.
http://dx.doi.org/10.3389/fpls.2021.6270...
). In plants, mainly the physiological responses come from (COR) cold-responsive genes. The expression of these genes is controlled by transcription factors (Fu et al., 2018FU, J., YU, J., JIANG, C. and CHENG, B., 2018. g‐C3N4‐Based heterostructured photocatalysts. Advanced Energy Materials, vol. 8, no. 3, pp. 1701503. http://dx.doi.org/10.1002/aenm.201701503.
http://dx.doi.org/10.1002/aenm.201701503...
).

Currently, the most well-defined cold signaling pathway is the C-repeats binding transcriptional factor (CBF) regulatory cascade (Thomashow, 1999THOMASHOW, M.F., 1999. Plant cold acclimation: freezing tolerance genes and regulatory mechanism. Annual Review of Plant Physiology and Plant Molecular Biology, vol. 50, no. 1, pp. 571-599. http://dx.doi.org/10.1146/annurev.arplant.50.1.571. PMid:15012220.
http://dx.doi.org/10.1146/annurev.arplan...
; Chinnusamy et al., 2003CHINNUSAMY, V., OHTA, M., KANRAR, S., LEE, B., HONG, X., AGARWAL, A. and ZHU, J.K., 2003. ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes & Development, vol. 17, no. 8, pp. 1043-1054. http://dx.doi.org/10.1101/gad.1077503. PMid:12672693.
http://dx.doi.org/10.1101/gad.1077503...
, 2007CHINNUSAMY, V., ZHU, J. and ZHU, J.K., 2007. Cold stress regulation of gene expression in plants. Trends in Plant Science, vol. 12, no. 10, pp. 444-451. http://dx.doi.org/10.1016/j.tplants.2007.07.002. PMid:17855156.
http://dx.doi.org/10.1016/j.tplants.2007...
; Fu et al., 2018FU, J., YU, J., JIANG, C. and CHENG, B., 2018. g‐C3N4‐Based heterostructured photocatalysts. Advanced Energy Materials, vol. 8, no. 3, pp. 1701503. http://dx.doi.org/10.1002/aenm.201701503.
http://dx.doi.org/10.1002/aenm.201701503...
; Knight and Knight, 2012KNIGHT, M.R. and KNIGHT, H., 2012. Low-temperature perception leading to gene expression and cold tolerance in higher plants. The New Phytologist, vol. 195, no. 4, pp. 737-751. http://dx.doi.org/10.1111/j.1469-8137.2012.04239.x. PMid:22816520.
http://dx.doi.org/10.1111/j.1469-8137.20...
). In our study, the transcription factor activity was found highly upregulated in sugarcane cultivar ROC22, 963 TF genes were significantly upregulated among a total of 5,649 upregulated genes, while 283 TF genes were found downregulated among the total of 3,289 downregulated genes. Similarly, in cultivar GT08-1108, about 974 TF genes were identified among 5,649 upregulated genes, while 293 TF genes were found among 3,289 downregulated genes. Transcription factors play a significant role under diverse stress conditions in plants by either down-regulation or up-regulation of genes at the transcriptional level (Javed et al., 2020JAVED, T., SHABBIR, R., ALI, A., AFZAL, I., ZAHEER, U. and GAO, S.J., 2020. Transcription factors in plant stress responses: challenges and potential for sugarcane improvement. Plants, vol. 9, no. 4, pp. 491. http://dx.doi.org/10.3390/plants9040491. PMid:32290272.
http://dx.doi.org/10.3390/plants9040491...
). Almost in the vascular plant genome 7% coding capacity is attributed to TFs for regulating the gene expression (Rushton et al., 2008RUSHTON, P.J., BOKOWIEC, M.T., LAUDEMAN, T.W., BRANNOCK, J.F., CHEN, X. and TIMKO, M.P., 2008. TOBFAC: the database of tobacco transcription factors. BMC Bioinformatics, vol. 9, pp. 53. http://dx.doi.org/10.1186/1471-2105-9-53. PMid:18221524.
http://dx.doi.org/10.1186/1471-2105-9-53...
). In plants, thousands of TFs families have been identified, while the major TFs families including WRKY, NAC, MYB, AP2/ERF, etc. mediated through signal transduction pathways have been used to cope with biotic and abiotic stresses in diverse crops (Javed et al., 2020JAVED, T., SHABBIR, R., ALI, A., AFZAL, I., ZAHEER, U. and GAO, S.J., 2020. Transcription factors in plant stress responses: challenges and potential for sugarcane improvement. Plants, vol. 9, no. 4, pp. 491. http://dx.doi.org/10.3390/plants9040491. PMid:32290272.
http://dx.doi.org/10.3390/plants9040491...
).

Although genome-wide RNA expression analysis has become a routine tool in biological research, extracting biological insight from such information remains a major challenge a powerful analytical method called Gene Enrichment Analysis (GEA) for interpreting gene expression data. The method derives its power by focusing on gene sets, that is, groups of genes that share common biological functions, chromosomal locations, or regulation (Subramanian et al., 2005SUBRAMANIAN, A., TAMAYO, P., MOOTHA, V.K., MUKHERJEE, S., EBERT, B.L., GILLETTE, M.A., PAULOVICH, A., POMEROY, S.L., GOLUB, T.R., LANDER, E.S. and MESIROV, J.P., 2005. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545-15550. http://dx.doi.org/10.1073/pnas.0506580102. PMid:16199517.
http://dx.doi.org/10.1073/pnas.050658010...
). In the GO enrichment data analysis, the transcription factor activity was annotated with GO categories related to protein binding, transcription factor binding, the activity of sequence-specific DNA binding transcription factor, activity in transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, transcription factor binding of RNA polymerase II to the regulatory region of sequence-specific DNA, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, transcription factor activity from the promoter of RNA polymerase II and RNA polymerase III.

Transcription factors (CBFs) are also known as DRE/CRT (Dehydration responsive elements) in the promoter region of cold-responsive genes, which regulate the expression as well as function (Stockinger et al., 1997STOCKINGER, E.J., GILMOUR, S.J. and THOMASHOW, M.F., 1997. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proceedings of the National Academy of Sciences, vol. 94, no. 3, pp. 1035-1040. https://doi.org/10.1073/pnas.94.3.1035.
https://doi.org/10.1073/pnas.94.3.1035...
; Liu et al., 1998LIU, Q., KASUGA, M., SAKUMA, Y., ABE, H., MIURA, S., YAMAGUCHI-SHINOZAKI, K. and SHINOZAKI, K., 1998. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. The Plant Cell, vol. 10, no. 8, pp. 1391-1406. http://dx.doi.org/10.1105/tpc.10.8.1391. PMid:9707537.
http://dx.doi.org/10.1105/tpc.10.8.1391...
; Vazquez-Hernandez et al., 2017VAZQUEZ-HERNANDEZ, M., ROMERO, I., ESCRIBANO, M.I., MERODIO, C. and SANCHEZ-BALLESTA, M.T., 2017. Deciphering the role of CBF/DREB transcription factors and dehydrins in maintaining the quality of table grapes cv. autumn royal treated with high CO2 levels and stored at 0 C. Frontiers of Plant Science, vol. 8, pp. 1591. http://dx.doi.org/10.3389/fpls.2017.01591. PMid:28970842.
http://dx.doi.org/10.3389/fpls.2017.0159...
). ICE1 (Inducer of CBF expression 1), an MYC-like bHLH transcriptional activator, acts as a positive regulator of transcription factors during cold signaling by binding to MYC recognition elements in the CBF promoter (Chinnusamy et al., 2003CHINNUSAMY, V., OHTA, M., KANRAR, S., LEE, B., HONG, X., AGARWAL, A. and ZHU, J.K., 2003. ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes & Development, vol. 17, no. 8, pp. 1043-1054. http://dx.doi.org/10.1101/gad.1077503. PMid:12672693.
http://dx.doi.org/10.1101/gad.1077503...
, 2010CHINNUSAMY, V., ZHU, J.K. and SUNKAR, R., 2010. Gene regulation during cold stress acclimation in plants. Methods in Molecular Biology, vol. 639, pp. 39-55. http://dx.doi.org/10.1007/978-1-60761-702-0_3. PMid:20387039.
http://dx.doi.org/10.1007/978-1-60761-70...
; Fursova et al., 2009FURSOVA, O.V., POGORELKO, G.V. and TARASOV, V.A., 2009. Identification of ICE2, a gene involved in cold acclimation which determines freezing tolerance in Arabidopsis thaliana. Gene, vol. 429, no. 1-2, pp. 98-103. http://dx.doi.org/10.1016/j.gene.2008.10.016. PMid:19026725.
http://dx.doi.org/10.1016/j.gene.2008.10...
). In our analysis, the transcription factor complex was highly upregulated by 4.9 logFc2. In GT08-1108, 630 upregulated genes for transcription factor activity, 43 for transcription factor activity in protein binding, 39 genes in RNAII polymerase transcription factor activity, 40 genes in transcription factor binding, 270 genes for transcription complex were recorded. the numbers of expressed upregulated transcription factors in ROC22 are 610 genes were upregulated for transcription factor activity, 43 for transcription factor activity in protein binding, 38 genes in RNAII polymerase transcription factor activity, 40 genes in transcription factor binding, 270 genes for transcription complex.

5. Conclusion

In Saccharum hybrid (Sugarcane), the lack of genome and transcriptome data has limited molecular analysis. In this research work, we provided the investigation about transcription factor activity during cold stress in two sugarcane cultivars namely GT08-1108 and ROC22. The results of this study showed that transcription factor genes play an important role in the expression of cold-responsive genes during cold stress in sugarcane. We further suggested to investigate the function of these TFS in drought, salinity, climate change, and extreme temperature in economically important grass crops like sugarcane, wheat, and rice. Finally, the development of cold and other environmental stress-tolerant sugarcane cultivars will be substantially beneficial for the sugar industry.

Acknowledgments

The authors are thankful to National Key Research and Development Program of China (2018YFD1000503), and the Higher Education Commission, Pakistan for supporting and funding this project.

References

  • AGARWAL, P.K., AGARWAL, P., REDDY, M.K. and SOPORY, S.K., 2006. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports, vol. 25, no. 12, pp. 1263-1274. http://dx.doi.org/10.1007/s00299-006-0204-8 PMid:16858552.
    » http://dx.doi.org/10.1007/s00299-006-0204-8
  • AHMAD, M., LI, J., YANG, Q., JAMIL, W., TENG, Y. and BAI, S., 2019. Phylogenetic, molecular, and functional characterization of PpyCBF proteins in Asian Pears (Pyrus pyrifolia). International Journal of Molecular Sciences, vol. 20, no. 9, pp. 2074. http://dx.doi.org/10.3390/ijms20092074 PMid:31035490.
    » http://dx.doi.org/10.3390/ijms20092074
  • AKHTAR, M., JAISWAL, A., TAJ, G., JAISWAL, J.P., QURESHI, M.I. and SINGH, N.K., 2012. DREB1/CBF transcription factors: their structure, function and role in abiotic stress tolerance in plants. Journal of genetics, vol. 91, no. 3, pp. 385-395. PMid:23271026.
  • BOHNERT, H.J., NELSON, D.E. and JENSEN, R.G., 1995. Adaptation to environmental stresses. The Plant Cell, vol. 7, no. 7, pp. 1099-1111. http://dx.doi.org/10.2307/3870060 PMid:12242400.
    » http://dx.doi.org/10.2307/3870060
  • BOLGER, A.M., LOHSE, M. and USADEL, B., 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, vol. 30, no. 15, pp. 2114-2120. http://dx.doi.org/10.1093/bioinformatics/btu170 PMid:24695404.
    » http://dx.doi.org/10.1093/bioinformatics/btu170
  • CHEAVEGATTI-GIANOTTO, A., ABREU, H.M.C., ARRUDA, P., BESPALHOK FILHO, J.C., BURNQUIST, W.L., CRESTE, S., CIERO, L., FERRO, J.A., FIGUEIRA, A.V.O., FILGUEIRAS, T.S., GROSSI-DE-SÁ, M.F., GUZZO, E.C., HOFFMANN, H.P., LANDELL, M.G.A., MACEDO, N., MATSUOKA, S., REINACH, F.C., ROMANO, E., SILVA, W.J., SILVA FILHO, M.C. and ULIAN, E.C., 2011. Sugarcane (Saccharum X officinarum): a reference study for the regulation of genetically modified cultivars in Brazil. Tropical Plant Biology, vol. 4, no. 1, pp. 62-89. http://dx.doi.org/10.1007/s12042-011-9068-3 PMid:21614128.
    » http://dx.doi.org/10.1007/s12042-011-9068-3
  • CHINNUSAMY, V., OHTA, M., KANRAR, S., LEE, B., HONG, X., AGARWAL, A. and ZHU, J.K., 2003. ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes & Development, vol. 17, no. 8, pp. 1043-1054. http://dx.doi.org/10.1101/gad.1077503 PMid:12672693.
    » http://dx.doi.org/10.1101/gad.1077503
  • CHINNUSAMY, V., ZHU, J. and ZHU, J.K., 2006. Gene regulation during cold acclimation in plants. Physiologia Plantarum, vol. 126, no. 1, pp. 52-61. http://dx.doi.org/10.1111/j.1399-3054.2006.00596.x
    » http://dx.doi.org/10.1111/j.1399-3054.2006.00596.x
  • CHINNUSAMY, V., ZHU, J. and ZHU, J.K., 2007. Cold stress regulation of gene expression in plants. Trends in Plant Science, vol. 12, no. 10, pp. 444-451. http://dx.doi.org/10.1016/j.tplants.2007.07.002 PMid:17855156.
    » http://dx.doi.org/10.1016/j.tplants.2007.07.002
  • CHINNUSAMY, V., ZHU, J.K. and SUNKAR, R., 2010. Gene regulation during cold stress acclimation in plants. Methods in Molecular Biology, vol. 639, pp. 39-55. http://dx.doi.org/10.1007/978-1-60761-702-0_3 PMid:20387039.
    » http://dx.doi.org/10.1007/978-1-60761-702-0_3
  • DHARSHINI, S., CHAKRAVARTHI, M., NARAYAN, A.J., MANOJ, V.M., NAVEENARANI, M., KUMAR, R., MEENA, M., RAM, B. and APPUNU, C., 2016. De novo sequencing and transcriptome analysis of a low temperature tolerant Saccharum spontaneum clone IND 00-1037. Journal of Biotechnology, vol. 231, pp. 280-294. http://dx.doi.org/10.1016/j.jbiotec.2016.05.036 PMid:27269250.
    » http://dx.doi.org/10.1016/j.jbiotec.2016.05.036
  • DU, Y.C., NOSE, A. and WASANO, K., 1999. Thermal characteristics of C4 photosynthetic enzymes from leaves of three sugarcane species differing in cold sensitivity. Plant & Cell Physiology, vol. 40, no. 3, pp. 298-304. http://dx.doi.org/10.1093/oxfordjournals.pcp.a029541
    » http://dx.doi.org/10.1093/oxfordjournals.pcp.a029541
  • EBRAHIM, M.K., ZINGSHEIM, O., EL-SHOURBAGY, M.N., MOORE, P.H. and KOMOR, E., 1998. Growth and sugar storage in sugarcane grown at temperatures below and above optimum. Journal of Plant Physiology, vol. 153, no. 5-6, pp. 593-602. http://dx.doi.org/10.1016/S0176-1617(98)80209-5
    » http://dx.doi.org/10.1016/S0176-1617(98)80209-5
  • FIGUEROA-RODRÍGUEZ, K.A., HERNÁNDEZ-ROSAS, F., FIGUEROA-SANDOVAL, B., VELASCO-VELASCO, J. and AGUILAR RIVERA, N., 2019. What has been the focus of sugarcane research? A bibliometric overview. International Journal of Environmental Research and Public Health, vol. 16, no. 18, pp. 3326. http://dx.doi.org/10.3390/ijerph16183326 PMid:31509963.
    » http://dx.doi.org/10.3390/ijerph16183326
  • FU, J., YU, J., JIANG, C. and CHENG, B., 2018. g‐C3N4‐Based heterostructured photocatalysts. Advanced Energy Materials, vol. 8, no. 3, pp. 1701503. http://dx.doi.org/10.1002/aenm.201701503
    » http://dx.doi.org/10.1002/aenm.201701503
  • FURSOVA, O.V., POGORELKO, G.V. and TARASOV, V.A., 2009. Identification of ICE2, a gene involved in cold acclimation which determines freezing tolerance in Arabidopsis thaliana. Gene, vol. 429, no. 1-2, pp. 98-103. http://dx.doi.org/10.1016/j.gene.2008.10.016 PMid:19026725.
    » http://dx.doi.org/10.1016/j.gene.2008.10.016
  • HUSSAIN, A., KHAN, Z.I., GHAFOOR, M.Y., ASHRAF, M. and RASHID, M.H., 2004. Review sugarcane, sugar metabolism and some abiotic stresses. International Journal of Agriculture And Biology, vol. 6, no. 4, pp. 732-742.
  • JAMES, G., 2008. Sugarcane 1st ed. Hoboken: John Wiely and Sons.
  • JAVED, T., SHABBIR, R., ALI, A., AFZAL, I., ZAHEER, U. and GAO, S.J., 2020. Transcription factors in plant stress responses: challenges and potential for sugarcane improvement. Plants, vol. 9, no. 4, pp. 491. http://dx.doi.org/10.3390/plants9040491 PMid:32290272.
    » http://dx.doi.org/10.3390/plants9040491
  • KHAN, M.T., KHAN, I.A., YASMEEN, S., NIZAMANI, G.S. and AFGHAN, S., 2019. Sugarcane Biofuels and Bioenergy Production in Pakistan: Current Scenario, Potential, and Future Avenues. In: M.T. KHAN and I.A. KHAN, eds. Sugarcane biofuels Cham: Springer, pp. 175-202.
  • KNIGHT, M.R. and KNIGHT, H., 2012. Low-temperature perception leading to gene expression and cold tolerance in higher plants. The New Phytologist, vol. 195, no. 4, pp. 737-751. http://dx.doi.org/10.1111/j.1469-8137.2012.04239.x PMid:22816520.
    » http://dx.doi.org/10.1111/j.1469-8137.2012.04239.x
  • LI, Y., WU, J., LI, X., ZHANG, R., LIU, X. and YANG, L., 2015. Damage in sugarcane production caused by long duration of chilling and frost in Guangxi China. International Journal of Agriculture Innovation and Research, vol. 3, no. 4, pp. 1139-1144.
  • LIU, Q., KASUGA, M., SAKUMA, Y., ABE, H., MIURA, S., YAMAGUCHI-SHINOZAKI, K. and SHINOZAKI, K., 1998. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. The Plant Cell, vol. 10, no. 8, pp. 1391-1406. http://dx.doi.org/10.1105/tpc.10.8.1391 PMid:9707537.
    » http://dx.doi.org/10.1105/tpc.10.8.1391
  • MATHUR, S., SUNOJ, V.S.J., ELSHEERY, N.I., REDDY, V.R., JAJOO, A. and CAO, K.F., 2021. Regulation of Photosystem II heterogeneity and Photochemistry in two cultivars of C4 crop Sugarcane under Chilling stress. Frontiers in Plant Science, vol. 12, pp. 627012. http://dx.doi.org/10.3389/fpls.2021.627012 PMid:33643354.
    » http://dx.doi.org/10.3389/fpls.2021.627012
  • MIURA, K. and FURUMOTO, T., 2013. Cold signaling and cold response in plants. International Journal of Molecular Sciences, vol. 14, no. 3, pp. 5312-5337. http://dx.doi.org/10.3390/ijms14035312 PMid:23466881.
    » http://dx.doi.org/10.3390/ijms14035312
  • ORGANIZATION FOR ECONOMIC COOPERATION AND DEVELOPMENT – OECD, 2016.Harmonisation of regulatory oversight in biotechnology safety assessment of transgenic organisms in the environment, volume 6: OECD Consensus Documents Paris: OECD Publishing.
  • PALVA, E.T., TAEHTIHARJU, S., TAMMINEN, I., PUHAKAINEN, T., LAITINEN, R., SVENSSON, J. and HEINO, P., 2002. Biological mechanisms of low temperature STRESS Responsse: cold acclimation and development of freezing tolerance in plants. Ibaraki: JIRCAS, pp. 9-15.
  • PAREEK, A., KHURANA, A., SHARMA, A.K. and KUMAR, R., 2017. An overview of signaling regulons during cold stress tolerance in plants. Current Genomics, vol. 18, no. 6, pp. 498-511. http://dx.doi.org/10.2174/1389202918666170228141345 PMid:29204079.
    » http://dx.doi.org/10.2174/1389202918666170228141345
  • PEREIRA-SANTANA, A., ALVARADO-ROBLEDO, E.J., ZAMORA-BRISEÑO, J.A., AYALA-SUMUANO, J.T., GONZALEZ-MENDOZA, V.M., ESPADAS-GIL, F., ALCARAZ, L.D., CASTAÑO, E., KEB-LLANES, M.A., SANCHEZ-TEYER, F. and RODRIGUEZ-ZAPATA, L.C., 2017. Transcriptional profiling of sugarcane leaves and roots under progressive osmotic stress reveals a regulated coordination of gene expression in a spatiotemporal manner. PLoS One, vol. 12, no. 12, pp. e0189271. http://dx.doi.org/10.1371/journal.pone.0189271 PMid:29228055.
    » http://dx.doi.org/10.1371/journal.pone.0189271
  • QUACKENBUSH, J., 2002. Microarray data normalization and transformation. Nature Genetics, suppl. 32, pp. 496-501. http://dx.doi.org/10.1038/ng1032 PMid:12454644.
    » http://dx.doi.org/10.1038/ng1032
  • REHMAN, S.U., MUHAMMAD, K., QUE, Y., REHMAN, A.U., NOVAES, E. and KHAN, S., 2019. presence of seventeen genes potentially involved in cold tolerance in sugarcane and Saccharum spontaneum genotypes. International Journal of Biosciences, vol. 14, no. 1, pp. 346-355. http://dx.doi.org/10.12692/ijb/14.1.346-355
    » http://dx.doi.org/10.12692/ijb/14.1.346-355
  • REHMAN, S.U., QUE, Y., KHAN, S., NOVAES, E. and 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, vol. 30, no. 4, pp. 913-922.
  • ROBINSON, M.D. and OSHLACK, A., 2010. A scaling normalization method for differential expression analysis of RNA-Seq Data. Genome Biology, vol. 11, no. 3, pp. R25. http://dx.doi.org/10.1186/gb-2010-11-3-r25 PMid:20196867.
    » http://dx.doi.org/10.1186/gb-2010-11-3-r25
  • ROBINSON, M.D. and SMYTH, G.K., 2008. Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics, vol. 9, no. 2, pp. 321-332. http://dx.doi.org/10.1093/biostatistics/kxm030 PMid:17728317.
    » http://dx.doi.org/10.1093/biostatistics/kxm030
  • ROBINSON, M.D., MCCARTHY, D.J. and SMYTH, G.K., 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, vol. 26, no. 1, pp. 139-140. http://dx.doi.org/10.1093/bioinformatics/btp616 PMid:19910308.
    » http://dx.doi.org/10.1093/bioinformatics/btp616
  • RUSHTON, P.J., BOKOWIEC, M.T., LAUDEMAN, T.W., BRANNOCK, J.F., CHEN, X. and TIMKO, M.P., 2008. TOBFAC: the database of tobacco transcription factors. BMC Bioinformatics, vol. 9, pp. 53. http://dx.doi.org/10.1186/1471-2105-9-53 PMid:18221524.
    » http://dx.doi.org/10.1186/1471-2105-9-53
  • SABIHA-HANIM, S. and HALIM, N.A.A., 2018. Sugarcane bagasse pretreatment methods for ethanol production. In T.P. BASSO and L.C. BASSO, eds. Fuel ethanol production from sugarcane London: INTECHOPEN, pp. 63-80.
  • SHI, Y., DING, Y. and YANG, S., 2018. Molecular regulation of CBF signaling in cold acclimation. Trends in Plant Science, vol. 23, no. 7, pp. 623-637. http://dx.doi.org/10.1016/j.tplants.2018.04.002 PMid:29735429.
    » http://dx.doi.org/10.1016/j.tplants.2018.04.002
  • SIPAÚBA-TAVARES, L.H., TEDESQUE, M.G. and SCARDOELI-TRUZZI, B., 2020. Evaluation of the effects of sugarcane molasses as a carbon source for Ankistrodesmus gracilis and Haematococcus pluvialis (Chlorophyceae). Brazilian Journal of Biology = Revista Brasileira de Biologia, vol. 80, no. 3, pp. 594-600. http://dx.doi.org/10.1590/1519-6984.216820 PMid:31644653.
    » http://dx.doi.org/10.1590/1519-6984.216820
  • SRIVASTAVA, S. and KUMAR, P., 2020. Abiotic stress responses and tolerance mechanisms for sustaining crop productivity in sugarcane. In: M. HASANUZZAMAN, ed. Agronomic crops Singapore: Springer, pp. 29-47. http://dx.doi.org/10.1007/978-981-15-0025-1_3
  • STOCKINGER, E.J., GILMOUR, S.J. and THOMASHOW, M.F., 1997. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proceedings of the National Academy of Sciences, vol. 94, no. 3, pp. 1035-1040. https://doi.org/10.1073/pnas.94.3.1035
    » https://doi.org/10.1073/pnas.94.3.1035
  • SU, W., REN, Y., WANG, D., SU, Y., FENG, J., ZHANG, C., TANG, H., XU, L., MUHAMMAD, K. and QUE, Y., 2020. The alcohol dehydrogenase gene family in sugarcane and its involvement in cold stress regulation. BMC Genomics, vol. 21, no. 1, pp. 521. http://dx.doi.org/10.1186/s12864-020-06929-9 PMid:32727370.
    » http://dx.doi.org/10.1186/s12864-020-06929-9
  • SUBRAMANIAN, A., TAMAYO, P., MOOTHA, V.K., MUKHERJEE, S., EBERT, B.L., GILLETTE, M.A., PAULOVICH, A., POMEROY, S.L., GOLUB, T.R., LANDER, E.S. and MESIROV, J.P., 2005. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545-15550. http://dx.doi.org/10.1073/pnas.0506580102 PMid:16199517.
    » http://dx.doi.org/10.1073/pnas.0506580102
  • TANG, S., YANG, L. and LI, Y., 2018. Comparative analysis on transcriptome among different sugarcane cultivars under low temperature stress. Shengwu Jishu Tongbao, vol. 34, no. 12, pp. 116-124.
  • THOMASHOW, M.F., 1999. Plant cold acclimation: freezing tolerance genes and regulatory mechanism. Annual Review of Plant Physiology and Plant Molecular Biology, vol. 50, no. 1, pp. 571-599. http://dx.doi.org/10.1146/annurev.arplant.50.1.571 PMid:15012220.
    » http://dx.doi.org/10.1146/annurev.arplant.50.1.571
  • TIAN, T., LIU, Y., YAN, H., YOU, Q., YI, X., DU, Z., XU, W. and SU, Z., 2017. Agrigo v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Research, vol. 45, no. W1, pp. W122-W129. http://dx.doi.org/10.1093/nar/gkx382 PMid:28472432.
    » http://dx.doi.org/10.1093/nar/gkx382
  • TSVETANOV, S., ATANASSOV, A. and NAKAMURA, C., 2000. Gold responsive gene/protein families and cold/freezing tolerance in cereals. Biotechnology, Biotechnological Equipment, vol. 14, no. 1, pp. 3-11. http://dx.doi.org/10.1080/13102818.2000.10819054
    » http://dx.doi.org/10.1080/13102818.2000.10819054
  • VAZQUEZ-HERNANDEZ, M., ROMERO, I., ESCRIBANO, M.I., MERODIO, C. and SANCHEZ-BALLESTA, M.T., 2017. Deciphering the role of CBF/DREB transcription factors and dehydrins in maintaining the quality of table grapes cv. autumn royal treated with high CO2 levels and stored at 0 C. Frontiers of Plant Science, vol. 8, pp. 1591. http://dx.doi.org/10.3389/fpls.2017.01591 PMid:28970842.
    » http://dx.doi.org/10.3389/fpls.2017.01591
  • WANG, D., WANG, L., SU, W., REN, Y., YOU, C., ZHANG, C., QUE, Y. and SU, Y., 2020. A class III WRKY transcription factor in sugarcane was involved in biotic and abiotic stress responses. Scientific Reports, vol. 10, no. 1, pp. 20964. http://dx.doi.org/10.1038/s41598-020-78007-9 PMid:33262418.
    » http://dx.doi.org/10.1038/s41598-020-78007-9
  • WANI, S.H. and HERATH, V., 2018. Cold tolerance in plants, physiological, molecular and genetic perspectives 1st ed. Cham: Springer. http://dx.doi.org/10.1007/978-3-030-01415-5
    » http://dx.doi.org/10.1007/978-3-030-01415-5
  • YANG, Y., ZHANG, X., SU, Y., ZOU, J., WANG, Z., XU, L. and QUE, Y., 2017. miRNA alteration is an important mechanism in sugarcane response to low-temperature environment. BMC Genomics, vol. 18, no. 1, pp. 833. http://dx.doi.org/10.1186/s12864-017-4231-3 PMid:29084515.
    » http://dx.doi.org/10.1186/s12864-017-4231-3
  • ZHANG, J., ZHANG, X., TANG, H., ZHANG, Q., HUA, X., MA, X., ZHU, F., JONES, T., ZHU, X., BOWERS, J., WAI, C.M., ZHENG, C., SHI, Y., CHEN, S., XU, X., YUE, J., NELSON, D.R., HUANG, L., LI, Z., XU, H., ZHOU, D., WANG, Y., HU, W., LIN, J., DENG, Y., PANDEY, N., MANCINI, M., ZERPA, D., NGUYEN, J.K., WANG, L., YU, L., XIN, Y., GE, L., ARRO, J., HAN, J.O., CHAKRABARTY, S., PUSHKO, M., ZHANG, W., MA, Y., MA, P., LV, M., CHEN, F., ZHENG, G., XU, J., YANG, Z., DENG, F., CHEN, X., LIAO, Z., ZHANG, X., LIN, Z., LIN, H., YAN, H., KUANG, Z., ZHONG, W., LIANG, P., WANG, G., YUAN, Y., SHI, J., HOU, J., LIN, J., JIN, J., CAO, P., SHEN, Q., JIANG, Q., ZHOU, P., MA, Y., ZHANG, X., XU, R., LIU, J., ZHOU, Y., JIA, H., MA, Q., QI, R., ZHANG, Z., FANG, J., FANG, H., SONG, J., WANG, M., DONG, G., WANG, G., CHEN, Z., MA, T., LIU, H., DHUNGANA, S.R., HUSS, S.E., YANG, X., SHARMA, A., TRUJILLO, J.H., MARTINEZ, M.C., HUDSON, M., RIASCOS, J.J., SCHULER, M., CHEN, L.Q., BRAUN, D.M., LI, L., YU, Q., WANG, J., WANG, K., SCHATZ, M.C., HECKERMAN, D., VAN SLUYS, M.A., SOUZA, G.M., MOORE, P.H., SANKOFF, D., VANBUREN, R., PATERSON, A.H., NAGAI, C. and MING, R., 2018. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nature Genetics, vol. 50, no. 11, pp. 1565-1573. http://dx.doi.org/10.1038/s41588-018-0237-2 PMid:30297971.
    » http://dx.doi.org/10.1038/s41588-018-0237-2
  • ZHOU, D., ZHOU, X., LING, Y., ZHANG, Z. and SU, Z., 2010. Agrigo: a GO analysis toolkit for the agricultural community. Nucleic Acids Research, vol. 38, suppl. 2, pp. W64-W70. PMid:20435677.
  • ZHOU, M.Q., SHEN, C., WU, L.H., TANG, K.X. and LIN, J., 2011. CBF-dependent signaling pathway: a key responder to low temperature stress in plants. Critical Reviews in Biotechnology, vol. 31, no. 2, pp. 186-92. http://dx.doi.org/10.3109/07388551.2010.505910 PMid:20919819.
    » http://dx.doi.org/10.3109/07388551.2010.505910

Publication Dates

  • Publication in this collection
    20 Dec 2021
  • Date of issue
    2023

History

  • Received
    28 Aug 2020
  • Accepted
    03 May 2021
Instituto Internacional de Ecologia R. Bento Carlos, 750, 13560-660 São Carlos SP - Brasil, Tel. e Fax: (55 16) 3362-5400 - São Carlos - SP - Brazil
E-mail: bjb@bjb.com.br