- Open Access
Comparative physiological and transcriptome analysis in cultivated and wild sugarcane species in response to hydrogen peroxide-induced oxidative stress
BMC Genomics volume 24, Article number: 155 (2023)
Sugarcane is an important energy crop grown worldwide,supplementing various renewable energy sources. Cultivated and wild sugarcane species respond differently to biotic and abiotic stresses. Generally, wild species are tolerant to various abiotic stresses. In the present study, the physiological and molecular responses of cultivated and wild sugarcane species to oxidative stress at the transcriptional levels were compared. Transcriptional responses were determined using RNAseq. The representative RNA-seq transcript values were validated by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and confirmed through physiological responses.
Oxidative stress causes leaf-rolling and -tip drying in cultivated sugarcane, but the wild species are tolerant. Higher chlorophyll fluorescence was observed in the wild species than that in the cultivated varieties under stress. Wild species can maintain a higher chlorophyll stability index than the cultivated species, which was confirmed by the lower transcripts of the chlorophyllase gene in the wild species than that in the cultivated variety. Transcription factor genes (NAC, MYB, and WRKY) were markedly expressed in response to oxidative stress, revealing their involvement in stress tolerance. The analysis revealed synchronized expression of acetyl-transferase, histone2A, cellulose synthase, and secondary cell wall biosynthetic genes in the wild species. The validation of selected genes and 15 NAC transcription factors using RT-qPCR revealed that their expression profiles were strongly correlated with RNA-seq. To the best of our knowledge, this is the first report on the oxidative stress response in cultivated and wild sugarcane species.
Physiological and biochemical changes in response to oxidative stress markedly differ between cultivated and wild sugarcane species. The differentially expressed stress-responsive genes are grouped intothe response to oxidative stress, heme-binding, peroxidase activity, and metal ion binding categories. Chlorophyll maintenance is a stress tolerance response enhanced by the differential regulation of the chlorophyllase gene.There is a considerable difference in the chlorophyll stability index between wild and cultivated varieties. We observed a substantial regulation of secondary wall biosynthesis genes in the wild species compared with that in the cultivated variety, suggesting differences in stress tolerance mechanisms.
Sugarcane is an important energy crop cultivated in tropical and subtropical regions globally. Sugarcane is a significant component of bioenergy crops in many tropical and subtropical economies . Approximately 100 countries produce sugarcane covering22 million ha (FAOSTAT, 2008). Increasing sugarcane productivity is essential for producing renewable energy sources and food sustainability.
Cultivated sugarcane is a complex hybrid of several species and subspecies. The majority of the commercial sugarcane varieties in India have approximately70–80%, 10–20%, and 8–13% of their chromosomes derived from Saccharum officinarum, Saccharum spontaneum, and interspecific recombination, respectively. A draft genome sequence assembly of the commercial sugarcane cultivar SP80-3280  and a reference sequence corresponding to the gene-rich regions of the primary (monoploid) sugarcane genome (R 570) of 382 Mb  are currently available. Sorghum (Sorghum bicolor)is an excellent alternative reference genome for the comparative analysis and annotation of genes. Comparative mapping between sugarcane and sorghum revealed a 95.2% sequence identity in the coding region.
Oxidative stress is an important factor in various abiotic stresses such as high temperature, drought, salinity, light, flooding, and heavy metals . Oxidative stress occurs due to an imbalance between the production of reactive oxygen species (ROS) and the antioxidant defense system in plant cells, leading to physiological and metabolic changes. Plants have evolved various molecular mechanisms to detect rapid changesand adapt. Depending on the dosage, H2O2 induces multiple responses in the cell. A high dosage leads to hypersensitive cell death [5,6,7,8], whereas a low dosage blocks cell cycle progression  which functions as a developmental signal for differentiation in the secondary cell wall .
Wild sugarcane species, such as S.spontaneum and Erianthus sp. tolerate several abiotic stresses and are sources of the tolerance genes in cultivated sugarcane through introgressive breeding. However, repeated backcrossing and selection for high yield have resulted in a narrow gene pool in the cultivated varieties; hence, the need to identify the resistant genes in the wild species foruse in breeding programs to develop stress-tolerant cultivated varieties. Several reports explained the use of hydrogen peroxide simulating oxidative stress in plants [11, 12]. Vijayalakshmi  used varying concentrations of H2O2 (0.05, 0.1, 0.15, 0.2 mM) to induce oxidative stress in rice.
Gene expression profiling through next-generation sequencing (NGS) has become an important tool for investigating how organisms respond to environmental changes. Understanding these reprogramming events is essential for gene discovery and introgression breeding programs aiming to enhance tolerance to multiple abiotic stresses in cultivated sugarcane. Several reports are available on the global gene expression profiles of sugarcane in response to drought , salinity and osmotic  and cold stresses . In our previous study , we reported the expression profiles of several genes involved in oxidative stress tolerance in sugarcane. The present study focused on comparative oxidative stress responses in cultivated and stress-tolerant sugarcane wild species (Erianthus sp.and S.spontaneum) to identify oxidative stress-responsive genes using transcriptomics. The present study provides an understanding of the molecular pathways of global gene expressions involved in abiotic stress tolerance, enabling the selection of sugarcane varieties adapted to marginal environments.
Morphological and physiological response of plants under oxidative stress
The cultivated and wild sugarcane species exhibited varied responses to oxidative stress
Erianthus sp. had an average of 27.1% chlorophyll fluorescence (CF) and S. spontaneum had a 23.7% increase in chlorophyll stability index (CSI) over the cultivated variety under stress conditions (500 ppm). (Fig. 1a). There were substantial differences in CSI and CF between different genotypes and H2O2 concentrations (Table 1). S. spontaneum exhibited considerable increases in H2O2 and superoxide dismutase (SOD) activity under stress conditions (Fig. 1b and c). Both wild species had low peroxidation values (Fig. 1d). There were considerable differences in SOD, peroxidase (POX), and lipid peroxidation (LP) concentrations between the different genotypes and H2O2 concentrations. S. spontaneum exhibited high proline concentration under stress conditions. A 22.5% average increase in proline concentration was observed in wild species (500 ppm 72 h treatment) compared with that in the cultivated variety (Fig. 1e) (Table 1). Overall, the wild species (IMP-564 and SES-90) and Co 86,032 cultivated variety exhibited considerable differences in physiological and biochemical parameters. There were no visible morphological changes between Erianthus sp.and S. spontaneum compared with those of the control, at 500 and 1,000 ppm, but leaf-lamina drying and leaf-rolling were observed at 500 and 1,000 ppm treatments (Fig. 2).
De novo assembly of transcripts and annotation of expressed genes
RNA samples from Erianthus sp., S. spontaneum, and the cultivated variety were sequenced using the Illumina NextSeq500 platform and generated ~ 121, ~ 120, and ~ 117 million reads (150 bases) with GC bases ranging from 55 to 57%. The N50 was 1,111 bp for each sample. Contigs were joined and assembled into unigenes (Table 2). The total number of assembled unigenes 3,07,358 out of which 2,23,896 were annotated based on sorghum genome (Sb) (TableS1) and sugarcane (varietyR570 / SP80-3280) as the reference (Table S2). Gene Ontology (GO) categorization of Erianthus sp. (Fig. 3a and b), S. spontaneum (Fig. 3c and d), and cultivated variety (Fig. 3e and f) exhibited that themajority of gene transcriptswere in the “molecular function” (46, 45, and 45%, respectively) and “cellular component” (41, 42, and 41%, respectively) and the least number of transcripts were in the “biological process” (12.9, 13.4, and 13.7%, respectively). There were moregene transcripts under the molecular function and biological process categories for cultivated varieties than those for the wild species. The GO for Saccharum sp., Erianthus sp. and cultivated variety is provided in Table S3.
Differentially expressed gene transcripts
We used the DESeq tool to identify differentially expressed genes between the control and oxidative stress conditions and filtered them based on their fold changes (> 2) and p-values (p < 0.05). We identified differentially expressed gene transcripts between species under oxidative stress conditions. There were 25, 51, and 30 upregulated and 105, 45, and 54 down regulated (farnesol dehydrogenase) gene transcripts in the Erianthus sp.; S. spontaneum; and cultivated variety, respectively. S. spontaneum exhibited relatively high expression values for the endo glucanasegene (FPKM [18,646]), glycan biosynthesis, lipid metabolism, and starch and sucrose metabolism. Few of the upregulated genes are ferredoxin nitrate reductase, cellulose synthase and peroxisome related genes. The complete lists of annotated differential gene expressions (DGEs) and their FPKM are provided in Table S4. The gene transcripts related to photosynthetic carbon fixation exhibited lower expression levels in the cultivated than that in wild species. Many differentially expressed transcripts were annotated as hypothetical genes; 94, 31, and 30 upregulated and 291, 139, and 297 down-regulated transcripts in Erianthus sp.; S. spontaneum; and cultivated variety, respectively, may be novel in the Saccharum complex (Table S5).
Differentially expressed transcription factor genes
In the present study, we explored important biotic and abiotic stress-responsive transcription factor families in sugarcane and other species. An average of 736 C3H (Zn finger domain) and 650 far-red responsive 1 (FAR1) transcripts were identified. FAR are positive regulators of chlorophyll biosynthesis via activating the expression of the HEMB1 gene . FAR1, C3H, MADS, bHLH, and NAC had high percentages in the control Erianthus sp. and low percentages in stressed S. Spontaneum (Table 3). Similarly, high and low percentages of orphans were observed in non-stressed Erianthus sp. and the cultivated variety, respectively. A high and low C2H2 TF percentage was observed in the stressed cultivated variety and stressed S.spontaneum, respectively. Other TFs, such as MYB-related and AP2- ethylene-responsive gene (EREBP), were high in then on-stressed cultivated variety and low in the stressed Erianthus sp. and non-stressed S.spontaneum. Among the other large TF families, NAC, MYB, and WRKY putative TFs were discussed considering FPKM. An average of 446 MYB, 434 WRKY, and 433 NAC transcripts were differentially expressed. The distribution of differentially regulated transcription factors (TFs) is illustrated in Fig. 4, and the details of the differentially expressed TF genes are provided in Table S6.The expression levels of three upregulated and nine downregulated (> twofold change in up- or down-regulated under oxidative stress conditions) NAC transcripts were validated using RT-qPCR. Table 4]). The comparative expression profiles of NAC genes (Fig. 5) exhibited ScNAC59 and ScNAC80 down regulation in cultivated varieties and up-regulation in wild species. Contrastingly, the ScNAC65 and ScNAC46 transcripts were down-regulated in wild species and upregulated in cultivated varieties. The comparison between the NGS values and qPCR-validated profiles for the cultivated and wild sugarcane species is provided in Table S7.
Gene expression validation through RT-qPCR
The expression values of four upregulated genes (dehydration responsive protein, cellulose synthase, calmodulin-binding protein, and cytochrome P450) and four down-regulated unigenes (ATP-dependent 6-phosphofructokinase (ATP-PFK), POX (EC 22.214.171.124), defensin, and ammonium transporter) were validated using RT-qPCR (Table S8). The list of genes selected for RT-qPCR validation ispresented in (Table 5), and their expression profiles are presented. The NGS data were consistent with those of the RT-qPCR for the transcripts corresponding to the dehydration-responsive protein,calmodulin-binding protein, and cellulose synthase; hence, their values are valid. Cytochrome P450 exhibited a change in expression based onthe NGS data. The POX and defensin expressions were down-regulated, confirming the reliability of the sequencing data. However, the ammonium transporter gene and ATP-PFK were upregulated in Erianthus sp. and the cultivated variety, respectively.
The present study reveals the physiological and molecular responsesof different species of sugarcane under oxidative stress. The entire plant molecular response was captured by transcriptome sequencing using the Illumina NextSeq500 platform. We selected two wild species (Erianthus arundinaceus and S. spontaneum) and hybrid cultivated varieties.
In our previous studies on oxidative , heat  and drought stresses , higher CF were observed in the wild species than that in the cultivated sugarcane under oxidative stress. The prolinecontent was 69 and 16.6% higher in S.Spontaneum and Erianthus sp., respectively, compared with those in their respective controls. LP is directly proportional to the extentof the plant damage. High H2O2 levels,accompanied by an increase in LP have been observed in young sugarcane plants during the initial growth phase under severe water stress . In the present study,the wild species exhibited low peroxidation values. For example, Erianthus sp. had a 23% decrease in peroxidation, and S.spontaneum exhibited 36% lower peroxidation than that of the cultivated variety at 72 h after treatment. Several similar observations have been reported in other plants concerning abiotic stress. A lower levels oflipid peroxidation was observed in heat stress-tolerant sugarcane varieties [18, 21] and oxidative stress-tolerant sugarcane varieties . Plants increase POX activity as a protective response during environmental stresses . A previous study reported increased POX activity during moisture stress in sugarcane . In the present study, cultivated and wild sugarcane species exhibited increased POX. In our previous study, we observed a decrease in POX at high H2O2 concentration (1,000 ppm) in cultivated varieties except S. spontaneum. . The prolonged exposure to oxidative stress, wild species maintained higher enzyme activities under stress conditionsthan that in the cultivated variety. We observed that of the wild species, S.spontaneum had a better physiological and biochemical support for increased stress tolerance.
In the differential expressed transcripts under oxidative stress, there were reduced transcripts in the “cellular component” category associated with chloroplast in all species. Most transcripts were associated with chloroplast thylakoids; however, some chloroplast genes were upregulated under stress conditions in S. spontaneum compared with those in the cultivated variety. Under the molecular function, transcripts belonging to ATP binding, transferase activity, and heme binding had a higher expression, but the transcripts belonging to the chlorophyllase gene, alpha-amylase, and calcium ion binding had a lower expression in the cultivated variety than that in Erianthus sp. . Based on CF and CSI physiological studies, we observed significant differences among the three genotypes (p < 0.001). However, E. arundinaceus had a higher CSI (27.1% increase in CF over that of the control) than that of the cultivated variety. Through RNAseq, it was observed that chlorophyllase gene transcripts were lower in wild species and higher in cultivated varieties, which may account for the degradation of chlorophyll. In the present physiological study, the CSI in the wild species was higher than that in the cultivated variety. CSI is an indicator of the maintenance of photosynthetic pigments. A previous study reported higher CSI in heat-tolerant sugarcane varieties . A previous study on rice  indicated that chlorophyll content was reduced in susceptible varieties compared with that in tolerant varieties under salt stress.
From the most significant differentially expressed genes (DEGs), we identified gene transcripts associated with stress-related pathways. A similar study reported the upregulation of DEGs related to ethylene, jasmonic acid, oxidative burst, NBS-LRR, cell wall modification, systematic acquired resistance, and pathogen-related proteins in sugarcane plantlets . Among the upregulated DEGs in Erianthus sp., GNAT (N-acetyl transferase) exhibited a 5-fold change. High expression of acetyl transferases increased the expression of cellulose synthase genesand conferred salt stress tolerance in plants . We observed an accumulation of cellulose synthase transcripts, which increased FPKM from 239 to 6,890 in S.spontaneum. S.spontaneum exhibited upregulation of the histone H2A gene (fivefold increase) under stress conditions. The histone acetylation-associated upregulation of cell wall-related genes during salt stress was observed in maize . The oxylipin biosynthetic process (GO:0,031,408) is essential for the metabolism of fatty acids in signal transduction for stomatal closure . We identified upregulation of oxylipin-related genes in Erianthus sp. and S.spontaneum, but not in cultivated varieties.
Wild species and sugarcane exhibited upregulation of gene transcripts related nitrate reductase (NR) is involved in the production of NO which act as signal molecular to induce the activities of antioxidant enzymes ). In S. spontaneum, we observed the upregulation of genes related to 2-oxocarboxylic acid metabolism, andthe downregulation of genes involved in phenylpropanoid biosynthesis and lipid metabolism. In the cultivated variety, upregulated genes were related to fatty acid degradation and phenylpropanoid biosynthesis, and downregulated genes were related to peroxisomes. The upregulation of genes related to the phenylpropanoid pathway increases lignin content .The regulation of phenylpropanoid biosynthetic genes was downregulation in wild species and upregulation in the cultivated variety. Although the wild species have higher lignin content than that in the cultivated variety , stress regulates the lignin synthesis genes to favor energy metabolism over lignin synthesis in wild species. Peroxisome metabolism improves plant stress tolerance and is involvedin many important metabolic reactions, particularly abiotic stress responses . Manipulation of peroxisomal scavenging systems for ROS may enhance plant fitness under environmental stress conditions .
Genes related to fatty acid metabolism were differentially expressed under oxidative stress. Previous reports have stated that cell membrane fatty acids are modulators of many signal transduction pathways activated by environmental stimuli . The high EREBP expression in wild species may promotes early floral meristem identity  and the transition from an inflorescence meristem to a floral meristem .Other genes (AP2_ARATH) were upregulated (AIL5) or downregulated (AP2) in the wild species. S.spontaneum revealed a high regulation of MYB-related genes (MY1R1_SOLTU) involved in the regulation of drought and salt tolerance and enhanced stomatal closure in response to abscisic acid (ABA ) and the MYB4 gene involved in cold stress .
The stress-responsive WRKY gene expression was similar in both wild species. Seven WRKY genes were upregulated in the wild species but downregulated in cultivated varieties. Differential regulation ofgenes between wild species and cultivated variety (upregulation of WRKY 35 and down regulation of WRKY 46) signals the conversion of sucrose (storage form) to reducing sugars in the cultivated variety . Stress-responsive NACs [40,41,42,43] show differential expression between wild and cultivated sugarcane species. The upregulation of ScNAC81 and ScNAC46 revealed the involvement of ABA-inducible leaf senescence signaling  and leaf senescence  in wild species. Over-expression of the NAC81 (GmNAC81) gene in soybean is involved in stress-induced leaf senescence . Down regulation of ScNAC89 may affect plant cell division .The expression of ScNAC78 revealed the induction of genes related to biosynthesis and was required for the accumulation of anthocyanins in response to high light stress . ScNAC55 was upregulated under stress in Erianthus sp., and similar results were reported in Brassica napus L . as NAC55 modulated ROS levels and cell death. Differential regulation of ScNAC69 (salt stress response , ScNAC59 (a senescence-associated gene involvedin salt and H2O2-dependent signaling pathways , ScNAC7 (involved in secondary wall biosynthesis and programmed cell death [52, 53] and ScNAC23 (a low-temperature responsive gene ) revealed their involvement in stress response in sugarcane. Putative stress-associated NAC genes have previously been reported in sugarcane . We observed varying expressions of NAC genes in response to oxidative stress, such as NAC78 . NAC7delays senescence and increases yield in maize . In our previous study, NAC expression was correlated with oxidative stress tolerance .
RT-qPCR analysis was conducted to verify the RNA-Seq data for 8 representative and 15 NAC genes. The expression profiles of the validated genes were strongly correlated with the RNAseq expression values, confirming the reliability of the transcriptome data. Dehydration-responsive proteins, upregulated in both Erianthus sp. and cultivated varieties, play a significant role in plant stress tolerance . The calmodulin-binding protein plays a crucial role in salt stress tolerance  in rice. Cytochrome P450 is upregulatedin all species involved in ROS metabolism . ATP-PFK and POX expression was downregulated under oxidative stress. Down regulation of the ATP-PFK gene inhibits glycolysis, promotes the production of nicotinamide adenine dinucleotide phosphate, thereby protecting against oxidative stress . The differences in expression between NGS and RT-qPCR may be due to the high sensitivity of RT-qPCR amplification.
The present study presents the differential responses of cultivated sugarcane and wild relatives to oxidative stress through physiological and transcriptomic analyses. The present study revealed that oxidative stress induced morphological changes in the sugarcane cultivated variety and was damaged at > 500 ppm H2O2. Molecular analyses of the DEGs indicated that most genes were under theresponse to oxidative stress, heme binding, POX activity, and metal ion binding categories.The differential regulation of chlorophyllase transcripts suggests substantial differences in the CSI for stress tolerance. At the cellular level, H2O2 induces secondary wall biosynthesis in wild species compared to cultivated plants, suggesting different stress tolerance mechanisms between cultivated and wild sugarcane species. The present study discovered highly expressed stress-responsive genes in wild types. Providing a basis for the incorporation of these genes in sugarcane to develop climate-resilient crops in the future.
Materials and methods
Induction of oxidative stress
A pot culture experiment was designed to measure the physiological response to oxidative stress in cultivated (Co 86032) and wild sugarcane types (S. spontaneum, [SES-90] and Erianthus sp. [IK -76–91]). The plants were allowed to grow for 65 daysunder greenhouse conditions and treated with two different concentrations of 30% H2O2 (500 ppm and 1,000 ppm). Pots were sprayed with H2O2 (approximately 1.5 L /pot) consecutively for 3 days (48 h and 72 h) at 8:00 am Indian Standard Time, and the experiment was replicated thrice.
Physiological response of plants under oxidative stress
Physiological changes in plants under oxidative stress were evaluated as previously described . Membrane damage due to peroxidation of the lipid bilayer was evaluated by measuring the malondialdehyde content . Proline,which may act as a direct scavenger of ROS or for osmolyte maintenance,was measured using the sulfosalicylic acid-acid ninhydrin method proposed by . The CSI,an indicator of photosynthetic pigment maintenance, was calculated as previously described . The evaluation was replicated thrice. The data on physiological parameters, such as CF, LP, and proline content and enzyme activities, such as POX and SOD, were obtained at 48 h and 72 h of H2O2 treatment. Significant differences between the treatments were tested statistically using the ANOVA (one-way/ two-way) method using JMP 9.0 software (JMP Statistical Discovery, Cary, NC, USA).
Sample collection and total RNA extraction
Tissues (young leaves and meristem tissues) were collected in triplicate (from 65 days old plants) from plants 48 h and 72 h after spraying with H2O2, and total RNA was extracted from all the samples using TRIzol reagent (Invitrogen, Waltham, MA, USA). The isolated RNAs were treated with DNase (Promega, Madison, WI, USA) to remove residual DNA contamination. RNA samples were separatedusing 1.5% (w/v) agarose gel electrophoresis (Takara Bio, Shiga, Japan). The quality and quantity of RNA were determined using the Nanodrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and Qubit fluorometer (Thermo Fisher Scientific), and the A260/A280 ratio was 1.9, indicating the purity of RNA (devoid of DNA and protein contamination). Biological replicate samples of the isolated RNAs were pooled before RNA-Seq library preparation.
Library preparation and sequencing
The Illumina compact transcriptome library was constructed using the NEXTflex Rapid Directional RNA-Seq library protocol (Cat # 5138–08 [PerkinElmer Applied Genomics, Waltham, MA, USA]) according to the manufacturer’s instructions. The prepared library was quantified using Qubit (Thermo Fisher Scientific) and validated for quality by running an aliquot on TapeStation system (High Sensitivity D1000 [Agilent Technologies, Santa Clara, CA, USA]).The observed RNA integrity value for all libraries was > 9 for all RNA samples, and the transcriptome libraries were sequenced (150 bp chemistry) using an Illumina NextSeq500 instrument (Illumina, San Diego, CA, USA) at Genotypic Technologies, Bangalore, India.
De-novo assembly of sequenced transcriptome
The initial quality of paired-end raw reads obtained from the Illumina sequencer was confirmed using the FASTQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) tool (Illumina). Unwanted regions in the reads (adapters, low-quality reads, and ambiguous bases ‘N) were trimmed, andhigh-quality trimmed reads were obtained for further analysis. The reads from each sample were normalized and assembled de novo separately using Trinity  (K-mer25[GitHub, San Francisco, CA, USA]). Trinity-generated assemblies were clustered based on sequence similarity. Transcripts were clustered using CD-HIT (cluster database at high identity with tolerance [GitHub]) at 95% identity and query coverage to reduce the redundancy without exclusion of sequence diversity. Clustered transcripts were used for further annotation.
Functional annotation of transcripts
Assembled transcripts were similarly searched against the non-redundant National Center for Biotechnology Information (NCBI), Clusters of Orthologous Group  and Uniprot  databases (European Bioinformatics Institute, Cambridge, UK) using the Basic Local Alignment Search Tool (BLAST; e-value:0.00001 [PubMed, Bethesda, MD, USA]).GO terms associated with the transcripts were identified using the UniProt database. The metabolic pathway enzymes expressed in the transcriptome were identified using the Kyoto encyclopaedia of genes and genomes (KEGG) database (  https://www.genome.jp/kegg/pathway.html) and the KEGG-Automatic Annotation Server tool.In the pathway analysis, a few reference organisms were selected, such as Arabidopsis thaliana, Citrus sinensis, Fragaria vesca, Glycine max, Oryza sativa japonica, Theobroma cacao, and Vitis vinifera.
Differential gene expression
The differential expression of RNA-Seq (Bioconductor, Boston, MA, USA) was used for DGE analysis between samples . Once the DGE was calculated, the genes with the expression log2-fold change greater than or less than twofold (up/downregulated) and with a p-value support of less than 0.01 were considered significant results.We used sugarcane (variety R570 / SP80-3280) and S. bicolor genomes as references for the annotation of DEGs. TF genes were classified into families, and orthologous relationships were identified using the TF genes (WRKY, MYB-related, AP2-EREBP, and others) from the PlantTFDB 5.0 (CBI, New Delhi, India ) comprising 320,370 sequences from 165 species of green plants and NAC TF sequences from the Grassius database . All these data were used as a backend database, and TFs were identified from the transcriptome using the NCBI BLAST command and ortholog identification. The pathways were mapped for all significant transcripts in all the DGE comparison sets.
RT-qPCR validation of DEGs
Differentially expressed transcripts between the three samples identified by RNA-Seq were confirmed by RT-qPCR to validate gene expression. Fifteen NAC genes (Table 4) and other representative DEGs (Table 5) were used. cDNA was prepared using a cDNA reverse transcription kit (Thermo Fisher). Forward and reverse primers were designed for these genes, and the elongation factor 2 gene was used as the reference gene. Three biological replicates were used for each sample. Among the identified NAC TFs, 15 were validated by RT-qPCR. The PrimeScript™ 1st strand cDNA synthesis kit (Takara Bio) was used for cDNA conversion. The 25S rRNA gene was used as a reference for the expression analysis of NAC genes. All RT-qPCR experiments were performed using the Rotor-GeneQ Qiagen Real-Time PCR system (Qiagen, Hilden, Germany) and QuantiNova™ SYBR Green PCR Kit(Qiagen). Relative quantification was performed using thedelta-delta Ct normalization method  and log2fold change (Log2FC).
Availability of data and materials
All data from the presented study are included in this published article and supplementary information files. The transcriptome data generated in the present study were deposited in the NCBI Sequence Read Archive database under accession number PRJNA606674. https://www.ncbi.nlm.nih.gov/sra.
Reactive oxygen species
National Center for Biotechnology Information
Differential gene expression
Basic local alignment search tool
- RT qPCR:
Reverse transcriptase quantitative polymerase chain reaction
Chlorophyll stability index
Kyoto encyclopedia of genes and genomes
Indian Council of Agricultural Research, Sugarcane Breeding Institute
Waclawovsky AJ, Sato PM, Lembke CG, Moore PH, Glaucia M. Sugarcane for bioenergy production : an assessment of yield and regulation of sucrose content. Plant Biotechnol J. 2010;263–76.
Riano-Pachon DM, Mattiello L. Draft genome sequencing of the sugarcane hybrid SP80-3280. F1000Res. 2017;6:861.
Garsmeur O, Droc G, Antonise R, Grimwood J, Potier B, Aitken K, et al. A mosaic monoploid reference sequence for the highly complex genome of sugarcane. Nat Commun. 2018;9:2638.
Kido ÉA, Ferreira Neto JRC, Silva RLDO, Pandolfi V, Guimares ACR, Veiga DT, et al. New insights in the sugarcane transcriptome responding to drought stress as revealed by supersage. Sci World J. 2012;2012:821062
Lamb C, Dixon RA. The oxidative burst in plant disease resistance. Annu Rev Plant Biol. 1997;48:251–75.
Jabs T, Dietrich RA, Dangl JL. Initiation of runaway cell death in an Arabidopsis mutant by extracellular superoxide. Science. 1996;273:1853–6 (80- ). (80- ).
Wang H, Li J, Bostock RM, Gilchrist DG. Apoptosis: A functional paradigm for programmed plant cell death induced by a host-selective phytotoxin and invoked during development. Plant Cell. 1996;8:375–91.
Alvarez ME, Pennell RI, Meijer PJ, Ishikawa A, Dixon RA, Lamb C. Reactive oxygen intermediates mediate a systemic signal network in the establishment of plant immunity. Cell. 1998;92:773–84.
Reichheld JP, Vernoux T, Lardon F, Van Montagu M, Inzé D. Specific checkpoints regulate plant cell cycle progression in response to oxidative stress. Plant J. 1999;17:647–56.
Potikha TS, Collins CC, Johnson DI, Delmer DP, Levine A. The involvement of hydrogen peroxide in the differentiation secondary walls in cotton fibers. Plant Physiol. 1999;119:849–58.
Rather BA, Masood A, Sehar Z, Majid A, Anjum NA, Khan NA. Mechanisms and Role of Nitric Oxide in Phytotoxicity-Mitigation of Copper. Front Plant Sci. 2020;11:675.
Vijayalakshmi D, Srividhya S, Muthulakshmi S, Satishraj R. Induction of Oxidative Stress by Hydrogen Peroxide Treatment in Rice Genotypes to Study the Osmolyte Accumulation Pattern and Antioxidant Capacity. J Stress Physiol Biochem. 2014;10:37–46.
Selvi A, Devi K, Manimekalai R, Prathima PT. Comparative analysis of drought-responsive transcriptomes of sugarcane genotypes with differential tolerance to drought. 3 Biotech. 2020;10:1–14.
Pereira-Santana A, Alvarado-Robledo EJ, Zamora-Briseño JA, Ayala-Sumuano JT, Gonzalez-Mendoza VM, Espadas-Gil F, et al. Transcriptional profiling of sugarcane leaves and roots under progressive osmotic stress reveals a regulated coordination of gene expression in a spatiotemporal manner. PLoS One. 2017;12:e0189271
Selvarajan D, Mohan C, Dhandapani V, Nerkar G, Jayanarayanan AN, Vadakkancherry Mohanan M, et al. Differential gene expression profiling through transcriptome approach of Saccharum spontaneum L. under low temperature stress reveals genes potentially involved in cold acclimation. 3 Biotech. 2018;8.
Manimekalai R, Narayanan J, Gokul M, Selvi A, Gomathi R, Arun kumar R. Biochemical and physiological response to oxidative stress in cultivated sugarcane and wild genera. Indian J Plant Physiol. 2018;23:261–70.
Lin R, Ding L, Casola C, Ripoll DR, Feschotte C, Wang H. Transposase-derived transcription factors regulate light signaling in Arabidopsis. Science. 2007;318:1302–5. (80- ).
Kohila S, Gomathi R. Adaptive physiological and biochemical response of sugarcane genotypes to high-temperature stress. Indian J Plant Physiol. 2018;23:1–16.
Sanghera GS, Kashyap L, Kumar R. Genetic Improvement of Sugarcane Through Non- Conventional Approaches. 2016. p. 1–25.
Boaretto LF, Carvalho G, Borgo L, Creste S, Landell MGA, Mazzafera P, et al. Water stress reveals differential antioxidant responses of tolerant and non-tolerant sugarcane genotypes. Plant Physiol Biochem PPB. 2014;74:165–75.
Abbas S, Ahmed S, Sabir S, Shah A. Detection of drought tolerant sugarcane genotypes (Saccharum officinarum) using lipid peroxidation, antioxidant activity, glycine-betaine and proline contents. J. Soil Sci. Plant Nutr. vol.14 no.1 Temuco mar. 2014 Epub 19-Ene-2014. http://dx.doi.org/10.4067/S0718-95162014005000019.
Narayanan J, Manimekalai R, Selvi A, et al. Physiological, Biochemical and Molecular Responses to Oxidative Stress in Saccharum spontaneum. Sugar Tech. 2023;25:282–93. https://doi.org/10.1007/s12355-022-01189-1.
Pujari DS, Chanda SV. Effect of salinity stress on growth, peroxidase and IAA oxidase activities in vigna seedlings. Acta Physiol Plant. 2002;24:435–9.
Vasantha S, Alarmelu S, Hemaprabha G, Shanthi RM. Evaluation of Promising Sugarcane Genotypes for Drought. Sugar Tech. 2005;7:82–3.
Liu X, Li L, Li M, Su L, Lian S, Zhang B, et al. AhGLK1 affects chlorophyll biosynthesis and photosynthesis in peanut leaves during recovery from drought. Sci Rep. 2018;8:2250.
Raju G, Shanmugam K, Kasirajan L. High-throughput sequencing reveals genes associated with high-temperature stress tolerance in sugarcane. 3 Biotech. 2020;10:1–13.
Mm M. Chlorophyll Stability Index (CSI): its impact on salt tolerance in rice. Int Rice Res Notes. 2000;25:38–9.
Santa Brigida AB, Rojas CA, Grativol C, de Armas EM, Entenza JOP, Thiebaut F, et al. Sugarcane transcriptome analysis in response to infection caused by Acidovorax avenae subsp. avenae. PLoS One. 2016;11:e0166473.
Zheng M, Liu X, Lin J, Liu X, Wang Z, Xin M, et al. Histone acetyltransferase GCN5 contributes to cell wall integrity and salt stress tolerance by altering the expression of cellulose synthesis genes. Plant J. 2019;97:587–602.
Savchenko T, Kolla VA, Wang CQ, Nasafi Z, Hicks DR, Phadungchob B, et al. Functional convergence of oxylipin and abscisic acid pathways controls stomatal closure in response to drought. Plant Physiol. 2014;164:1151–60.
Jardim-Messeder D, da Franca ST, Fonseca JP, Junior JN, Barzilai L, Felix-Cordeiro T, et al. Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes. Mol Genet Genomics. 2020;295:717–39.
Llerena JPP, Figueiredo R, Brito M dos S, Kiyota E, Mayer JLS, Araujo P, et al. Deposition of lignin in four species of Saccharum. Sci Rep. 2019. https://doi.org/10.1038/s41598-019-42350-3.
Kessel-Vigelius SK, Wiese J, Schroers MG, Wrobel TJ, Hahn F, Linka N. An engineered plant peroxisome and its application in biotechnology. Plant Sci. 2013;210:232–40.
Nielsen J. Metabolic engineering: Techniques for analysis of targets for genetic manipulations. Biotechnol Bioeng. 1998;58:125–32.
Jofuku KD, Den Boer BGW, Van Montagu M, Okamuro JK. Control of arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994;6:1211–25.
Drews GN, Bowman JL, Meyerowitz EM. Negative regulation of the Arabidopsis homeotic gene AGAMOUS by the APETALA2 product. Cell. 1991. https://doi.org/10.1016/0092-8674(91)90551-9.
Shin D, Moon SJ, Han S, Kim BG, Park SR, Lee SK, et al. Expression of StMYB1R-1, a novel potato single MYB-like domain transcription factor, increases drought tolerance. Plant Physiol. 2011;155:421–32.
Vannini C, Locatelli F, Bracale M, Magnani E, Marsoni M, Osnato M, et al. Overexpression of the rice Osmyb4 gene increases chilling and freezing tolerance of Arabidopsis thaliana plants. Plant J. 2004. https://doi.org/10.1046/j.1365-313X.2003.01938.x.
Sun C, Palmqvist S, Olsson H, Borén M, Ahlandsberg S, Jansson C. A novel WRKY transcription factor, SUSIBA2, participates in sugar signaling in barley by binding to the sugar-responsive elements of the iso1 promoter. Plant Cell. 2003;15:2076–92.
Shao H, Wang H, Tang X. NAC transcription factors in plant multiple abiotic stress responses: Progress and prospects. Front Plant Sci. 2015;6:902
Nuruzzaman M, Manimekalai R, Sharoni AM, Satoh K, Kondoh H, Ooka H, et al. Genome-wide analysis of NAC transcription factor family in rice. Gene. 2010;465:30–44.
Fujita M, Fujita Y, Maruyama K, Seki M, Hiratsu K, Ohme-Takagi M, et al. A dehydration-induced NAC protein, RD26, is involved in a novel ABA-dependent stress-signaling pathway. Plant J. 2004;39:863–76.
Kadier YZuY, yi, Dai Q min, Song G, Lin S wen, Sun Q peng, et al. Genome-wide identification, classification and expression analysis of NAC family of genes in sorghum [Sorghum bicolor (L.) Moench]. Plant Growth Regul. 2017;83:301–12.
Takasaki H, Maruyama K, Takahashi F, Fujita M, Yoshida T, Nakashima K, et al. SNAC-As, stress-responsive NAC transcription factors, mediate ABA-inducible leaf senescence. Plant J. 2015;84:1114–23.
Oda-Yamamizo C, Mitsuda N, Sakamoto S, Ogawa D, Ohme-Takagi M, Ohmiya A. The NAC transcription factor ANAC046 is a positive regulator of chlorophyll degradation and senescence in Arabidopsis leaves. Sci Rep. 2016;6:23609.
Pimenta MR, Silva PA, Mendes GC, Alves JR, Caetano HDN, Machado JPB, et al. The Stress-Induced Soybean NAC Transcription Factor GmNAC81 Plays a Positive Role in Developmentally Programmed Leaf Senescence. Plant Cell Physiol. 2016;57:1098–114.
Smyczynski C, Roudier F, Gissot L, Vaillant E, Grandjean O, Morin H, et al. The C terminus of the immunophilin PASTICCINO1 is required for plant development and for interaction with a NAC-like transcription factor. J Biol Chem. 2006;281:25475–84.
Morishita T, Kojima Y, Maruta T, Nishizawa-Yokoi A, Yabuta Y, Shigeoka S. Arabidopsis NAC transcription factor, ANAC078, regulates flavonoid biosynthesis under high-light. Plant Cell Physiol. 2009;50:2210–22.
Niu F, Wang C, Yan J, Guo X, Wu F, Yang B, et al. Functional characterization of NAC55 transcription factor from oilseed rape (Brassica napus L.) as a novel transcriptional activator modulating reactive oxygen species accumulation and cell death. Plant Mol Biol. 2016;92:89–104.
Park J, Kim YS, Kim SG, Jung JH, Woo JC, Park CM. Integration of auxin and salt signals by the NAC transcription factor NTM2 during seed germination in arabidopsis. Plant Physiol. 2011;156:537–49.
Balazadeh S, Kwasniewski M, Caldana C, Mehrnia M, Zanor MI, Xue GP, et al. ORS1, an H2O2-responsive NAC transcription factor, controls senescence in arabidopsis thaliana. Mol Plant. 2011;4:346–60.
Kubo M, Udagawa M, Nishikubo N, Horiguchi G, Yamaguchi M, Ito J, et al. Transcription switches for protoxylem and metaxylem vessel formation. Genes Dev. 2005. https://doi.org/10.1101/gad.1331305.
Zhou J, Zhong R, Ye ZH. Arabidopsis NAC domain proteins, VND1 to VND5, are transcriptional regulators of secondary wall biosynthesis in vessels. PLoS One. 2014;9:e105726
Ditt RF, Gentile A, Tavares RG, Camargo SR, Fernandez JH, da Silva MJ, et al. Analysis of the stress-inducible transcription factor SsNAC23 in sugarcane plants. Sci Agric. 2011;68:454–61.
Ramaswamy M, Narayanan J, Manickavachagam G, Athiappan S, Arun M, Gomathi R, et al. Genome wide analysis of NAC gene family ‘sequences’ in sugarcane and its comparative phylogenetic relationship with rice, sorghum, maize and Arabidopsis for prediction of stress associated NAC genes. Agri Gene. 2017;3:1–11.
Gladman NP, Marshall RS, Lee KH, Vierstra RD. The proteasome stress regulon is controlled by a pair of NAC transcription factors in arabidopsis. Plant Cell. 2016;28:1279–96.
Zhang J, Fengler KA, Van Hemert JL, Gupta R, Mongar N, Sun J, et al. Identification and characterization of a novel stay-green QTL that increases yield in maize. Plant Biotechnol J. 2019;17:2272–85.
Manimekalai R, Narayanan J, Ranjini R, Gokul M, Selvi A, Kumar P, et al. Hydrogen Peroxide-Induced Oxidative Stress in Sugarcane and Response Expression Pattern of Stress-Responsive Genes Through Quantitative RT-PCR. Sugar Tech. 2018;20:681–91.
Lata C, Prasad M. Role of DREBs in regulation of abiotic stress responses in plants. J Exp Bot. 2011;62:4731–48.
Yuenyong W, Chinpongpanich A, Comai L, Chadchawan S, Buaboocha T. Downstream components of the calmodulin signaling pathway in the rice salt stress response revealed by transcriptome profiling and target identification. BMC Plant Biol. 2018;18:335.
Veith A, Moorthy B. Role of cytochrome P450s In the generation and metabolism of reactive oxygen species. Current Opinion in Toxicology. 2018;7:44–51.
Mullarky E, Cantley LC. Diverting Glycolysis to Combat Oxidative Stress. In: Innovative Medicine. 2015. p. 3–23.
Dhindsa RS, Plumb-dhindsa P, Thorpe TA. Leaf senescence: Correlated with increased levels of membrane permeability and lipid peroxidation, and decreased levels of superoxide dismutase and catalase. J Exp Bot. 1981;32:93–101.
Bates LS, Waldren RP, Teare IDPS. 39:205–7. Rapid determination of free proline for water-stress studies. Plant Soil. 1973;39:205–7.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.
Galperin MY, Makarova KS, Wolf YI, Koonin E V. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43 Database issue:D261–9.
The UniProt Consortium. UniProt: a worldwide hub of protein knowledge Nucleic Acids Res. 2019;47:D506–515.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28:27–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Jin J, Tian F, Yang DC, Meng YQ, Kong L, Luo J, et al. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017;45:D1040-5.
Yilmaz A, Nishiyama MY, Fuentes BG, Souza GM, Janies D, Gray J, et al. GRASSIUS: A platform for comparative regulatory genomics across the grasses. Plant Physiol. 2009;149:171–80.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001. https://doi.org/10.1006/meth.2001.1262.
The authors would like to thank the Director of theICAR-SBI, and Coimbatore for permissionto carry out the presentstudy.
The research was funded by the ICAR-Sugarcane Breeding institute.
Ethics approval and consent to participate
The research followed the Ethics guidelines, and the authors consented to publish the manuscript. The study protocols complied with the relevant institutional, national, and international guidelines and legislation. Plant materials were obtained from the Indian Council of Agricultural Research, Sugarcane Breeding Institute (ICAR-SBI) germplasm collections. These descriptors are available in the repository of the ICAR-SBI. Permissions to use the wild type (S. spontaneum, SES-90) for the research were obtained from the Director of the ICAR-SBI, Coimbatore, India.
Consent for publication
The authors declare no financial and non-financial competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Supplementary file. I:
Summary of Annotation for DEGs of Erianthus arundinaceus, Saccharum spontaneum and cultivated variety.
Additional file 2. Supplementary file. II:
Annotation of unigenes of Erianthus arundinaceus, Saccharum spontaneum and cultivated variety based on Saccharum sp as reference.
Additional file 3. Supplementary file. III:
GO terms for DEGs of Erianthus arundinaceus, Saccharum spontaneum and cultivated variety.
Additional file 4. Supplementary file. IV:
List of annotated DGEs with their FPKM for Erianthus arundinaceus, Saccharum spontaneum and cultivated variety.
Additional file 5. Supplementary file. V:
List of novel DEG of Erianthus arundinaceus, Saccharum spontaneum and cultivated variety.
Additional file 6. Supplementary file. VI:
Annotated differentially expressed AP2-ER-EBP, MYB related, WRKY and NAC transcription factor genes with FPKM values.
Additional file 7. Supplementary File. VII:
Validation of expression of NAC genes through RT-PCR.
Additional file 8. Supplementary File. VIII:
Validation of expression of representative genes through RT-PCR.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Manimekalai, R., Selvi, A., Narayanan, J. et al. Comparative physiological and transcriptome analysis in cultivated and wild sugarcane species in response to hydrogen peroxide-induced oxidative stress. BMC Genomics 24, 155 (2023). https://doi.org/10.1186/s12864-023-09218-3
- Erianthus sp
- Saccharum sp
- sugarcane transcriptome
- NAC genes