Skip to main content
Download PDF

Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection

BMC Plant Biology volume 23, Article number: 52 (2023) Cite this article

Abstract

Background

Diseases are the major factor affecting the quality and yield of sugarcane during its growth and development. However, our knowledge about the factors regulating disease responses remain limited. The present study focuses on identifying genes regulating transcriptional mechanisms responsible for resistance to leaf scald caused by Xanthomonas albilineans in S. spontaneum and S. officinarum.

Results

After inoculation of the two sugarcane varieties SES208 (S. spontaneum) and LA Purple (S. officinarum) with Xanthomonas albilineans, SES208 exhibited significantly greater resistance to leaf scald caused by X. albilineans than did LA Purple. Using transcriptome analysis, we identified a total of 4323 and 1755 differentially expressed genes (DEGs) in inoculated samples of SES208 and LA Purple, respectively. Significantly, 262 DEGs were specifically identified in SES208 that were enriched for KEGG pathway terms such as plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction. Furthermore, we built a transcriptional regulatory co-expression network that specifically identified 16 and 25 hub genes in SES208 that were enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction. All of these essential genes might be significantly involved in resistance-regulating responses in SES208 after X. albilineans inoculation. In addition, we found allele-specific expression in SES208 that was associated with the resistance phenotype of SES208 when infected by X. albilineans. After infection with X. albilineans, a great number of DEGs associated with the KEGG pathways ‘phenylpropanoid biosynthesis’ and ‘flavonoid biosynthesis’ exhibited significant expression changes in SES208 compared to LA Purple that might contribute to superior leaf scald resistance in SES208.

Conclusions

We provided the first systematical transcriptome map that the higher resistance of SES208 is associated with and elicited by the rapid activation of multiple clusters of defense response genes after infection by X. albilineans and not merely due to changes in the expression of genes generically associated with stress resistance. These results will serve as the foundation for further understanding of the molecular mechanisms of resistance against X. albilineans in S. spontaneum.

Peer Review reports

Very high yields of sucrose make Saccharum officinarum and other sugarcane species among the most important crops as the largest sources of sugar and biofuel feedstocks from warm climates around the world [1, 2]. In the past few decades, sugarcane has been used not only to generate sugars for food but also to develop potentially sustainable processes to shift the global energy mix toward cleaner and renewable resources [3]. Current cultivated sugarcane hybrids were derived from interspecific crosses between several parental species including S. barberi, S. sinense, S. officinarum, and wild accessions of S. spontaneum and S. robustum that resulted in complex autopolyploids with 2n = 100–130 chromosomes [4, 5]. The growth and productivity of sugarcane can be very adversely impacted by biotic stress caused by microorganisms and also by abiotic stresses [6]. Researchers have taken several decades to breed cultivated varieties for high sucrose yield and disease resistance. In particular, the high sucrose content of modern sugarcane cultivars was derived from S. officinarum, which exhibits poor disease resistance, while their disease resistance was derived from its relative S. spontaneum, which exhibits low sugar content. S. spontaneum is a key contributor of outstanding traits such as tolerance to environmental stress and ratoon capacity and was therefore used during a program of backcrossing to S. officinarum, referred to as nobilization [7]. Thus, considering the genetic background of sugarcane, analysis of the mechanisms of its parental responses to pathogenic bacteria has become necessary. We therefore consider the S. officinarum cultivar LA Purple and the S. spontaneum cultivar SES208 ideal materials for genetic research into high sugar content and disease resistance traits.

As one of the main diseases of sugarcane (Saccharum spp. hybrids), leaf scald can cause the death of the entire sugarcane plant. Due to the potential for leaf scald to cause severe economic losses, this disease has severely restricted the cultivation of susceptible sugarcane varieties [8]. Leaf scald was first recognized as a disease caused by an organism in the Proteobacteria, Xanthomonas albilineans. Symptoms of infection with this pathogen in sugarcane have allowed the progression of this disease to be classified into three different stages including incubation period, chronic period, and acute period [9]. X. albilineans colonizes the plant vascular system and parenchyma cells [10] and leads to the inhibition of chloroplast development and disruption of photosynthesis [11], followed by “scalding” of stalks and stalk death. A hypersensitive response (HR) is induced early during infection with X. albilineans. This pathogen-induced HR is thought to reduce infiltration and spread of the pathogen through local cell death at the site of infection [12, 13].

During their evolution, plants have established an immune system against pathogen infection [14]. Plants live in environments that also host some microorganisms that, as pathogens, can adversely affect their growth and development. To avoid infection by pathogens in their environments, plants respond via an innate immune system [15]. Plants first recognize molecules originating from both pathogenic and non-pathogenic microbes and then respond to virulence factors from pathogens. The primary immune response in plants, known as PAMP-triggered immunity (PTI), involves the recognition of conserved microbial features. PTI is initiated by the perception of pathogens at the cell surface and pathogens that cross the host plasma membrane, which are recognized according to their pathogen-associated molecular patterns (PAMPs) [16]. The hypersensitive response (HR) is the hallmark response of effector-triggered immunity, or ETI, which is initiated directly or indirectly when effectors from pathogens are recognized by NLRs, which are receptors containing nucleotide-binding domains and leucine-rich repeats [15]. Crosstalk between PTI and ETI involves mitogen-activated protein kinase (MAPK) signaling [17] and results in massive accumulation of active oxygen within a short time [18]. A substantial influx of Ca2+ [19] and the biosynthesis of plant hormones, including salicylic acid (SA), abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA) also occurs [20]. The beta-hydroxy phenolic compound SA is known as a defense-related hormone that mediates systemic acquired resistance (SAR) [13]. In 1979, aspirin (acetyl-SA) was first reported to increase the resistance of a susceptible tobacco line to infection with TMV (tobacco mosaic virus) [21]. Further, JA- and SA-induced systemic resistance (ISR) can reprogram plant induction mechanisms to alleviate stresses [22]. Plant secondary metabolites, which can be expressed constitutively or in response to pathogens, mediate interactions between plants and microbes in particular [23]. Studies conducted on the activities of these metabolites during disease resistance in plants have shown that many metabolites including phenylpropanoids, terpenoid quinones, and flavonoids can potentially combat plant pathogens [23, 24].

Several studies of plants infected with X.albilineans have been carried out in species such as rice [25, 26], elephant grass [27], and wheat [28], and shown that the expression of some genes involved in pathways controlling the biosynthesis and activity of ethylene, jasmonic acid, or secondary metabolites respond to inoculation with this pathogen. Although transcripts of more than 5000 genes with altered expression were previously revealed in an analysis comparing the transcriptomes of sugarcane cultivars that are resistant or susceptible to leaf scald after infection with X. albilineans [29], the cellular mechanisms of resistance to this pathogen in sugarcane are still unknown. The first global sugarcane proteome dataset revealed that seven candidate genes associated with “photosynthesis”, “glycolytic process”, “glycosylation process”, “plant innate immune system”, “plant cytochrome P450”, and “non-specific lipid transfer proteins” are expressed in response to X. albilineans inoculation in SES208 [30]. However, although some infected plants do not display external symptoms, especially under conditions unsuitable for the pathogen, the reasons for such latent progression through infection phases are not known [9]. Therefore, a qualitative PCR assay was established to detect X. albilineans infection, exclude the presence of other pathogens, and monitor epidemics and evaluate resistance in the field [31]. Despite such progress, the mechanisms of the molecular responses of modern sugarcane cultivars to infection with X.albilineans are still unclear.

To analyze the origins of X. albilineans resistance in modern sugarcane varieties from a cross between S. spontaneum exhibiting disease resistance and S. officinarum exhibiting high sugar and biomass yields, the leaf scald resistance of these two parental species was compared. In the present study, we used RNA-Seq and weighted gene correlation network analysis (WGCNA) to compare the transcriptomes of the representative sugarcane cultivars SES208 (Saccharum spontaneum, 2n = 64) and LA Purple (Saccharum officinarum, 2n = 80) during their genetic responses to infection with X. albilineans. Our analysis identified transcripts of DEGs associated with X. albilineans resistance that will be of great significance for improving our understanding of X. albilineans resistance in sugarcane.

Results

Differences in responses of sugarcane cultivars LA purple and SES208 to inoculation with X. albilineans

To investigate the differences in resistance to leaf scald caused by X. albilineans between the sugarcane cultivars SES208 (S. spontaneum), a disease-resistant variety, and LA Purple (S. officinarum), a disease-susceptible variety, we experimentally inoculated plants with X. albilineans. We observed the phenotypic differences between inoculated and uninoculated plants of these two species and collected leaf tissue samples at 24 hpi (hpi: hours post inoculation), 48 hpi, and 72 hpi for subsequent RNA extraction. Plants of SES208 and LA Purple exhibited different responses to inoculation with X. albilineans. Although there were no leaf scald symptoms on SES208 plants after infection, there were thin streaks on the leaves. LA Purple plants exhibited symptoms of leaf scald at 72 hpi infection (Fig. 1a), thus confirming the differences between LA Purple and SES208 in resistance to leaf scald caused by X. albilineans.

Fig. 1
figure 1

a Leaf scald symptoms observed in sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. Pencil line streak (red arrow). b Box diagram of gene expression (log10(FPKM+ 1)). c The numbers of differentially expressed genes (DEGs) with increased or decreased transcript abundance in SES208 and LA Purple (|FDR| > 2,p-value< 0.005). d Venn diagrams of DEGs in SES208 and LA Purple, respectively. SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

Identification of DEGs and core conserved DEGs in resistant and susceptible sugarcane cultivars in response to inoculation with X. albilineans

Global gene expression profiles in leaves of SES208 (S. spontaneum) and LA Purple (S. officinarum) inoculated with X. albilineans and uninoculated control plants were analyzed at 24 hpi, 48 hpi, and 72 hpi. In general, global gene expression in SES208 was higher than that in LA Purple (Fig. 1b). Transcriptome comparisons at each time point were made within each of the cultivars. Compared to time point 0, at 24 hpi, transcripts of 1480 DEGs were more abundant and 939 DEGs were less abundant in leaves of SES208 while transcripts of 104 DEGs were more abundant and 182 DEGs were less abundant in leaves of LA-Purple (Fig. 1c). Trends in the number of genes with increased transcript abundance and those with decreased transcript abundance were similar between 0 and 72 hpi in LA Purple but differed during this time period in SES208 (Fig. 1c). Between 0 and 48 hpi, SES208 showed 735 genes with increased transcript abundance and 1538 genes with decreased transcript abundance (Fig. 1c). The number of genes with increased transcript abundance in SES208 was significantly greater than that in LA Purple (Fig. 1c). The highest number of genes with increased transcript abundance in SES208 occurred at 24 hpi, but in LA Purple occurred at 72 hpi (Fig. 1c). The highest number of DEGs with decreased transcript abundance occurred at 48 hpi in leaves of SES208 and at 72 hpi in leaves of LA Purple (Fig. 1c). A total of 4323 DEGs were expressed in leaves of SES208 in response to X. albilineans, nearly twice the number of DEGs in LA Purple in response to this pathogen (Fig. 1d). The expression of 975, 646, and 714 DEGs specifically changed at 24, 48, and 72 hpi, respectively, in SES208, while the expression of 70, 158, and 1053 DEGs specifically changed in LA Purple at 24, 48, and 72 hpi, respectively (Fig. 1d). These results indicate that the transcript expression of more genes in SES208 than in LA Purple changed rapidly during infection with X. albilineans.

Further, we annotated the genes encoding TFs among the DEGs responding to inoculation with X. albilineans in LA Purple and SES208. We identified genes representing a total of 23 TF families in both cultivars, and found that more TF families exhibited changed transcript expression in SES208 than in LA Purple (Fig. 2a). It is worth noting that gene families encoding WRKY, NAC, MYB, ERF, and bHLH were significantly overrepresented among DEGs in SES208 inoculated with X. albilineans. All DEGs encoding TFs showed higher expression in SES208 than in LA Purple except for three genes encoding ERF, MYB, and WRKY, respectively (Fig. 2b). The transcript abundances of 15 genes encoding MYBs, 17 encoding NACs, 8 encoding ERFs, and 8 encoding WRKYs significantly increased, while the transcript abundances of 12 genes encoding bHLHs and 12 encoding ERFs significantly decreased in SES208 compared to LA Purple after inoculation with X. albilineans (Fig. 2b), suggesting that members of these TF families might be involved in the regulation of responses to this pathogen in sugarcane.

Fig. 2
figure 2

Transcription factors (TFs) encoded by DEGs in sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. a The number of TFs encoded by DEGs. b Heat map of TF expression after inoculation with X. albilineans. Expression values were calculated as the relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

GO term and KEGG pathway enrichment analysis of DEGs responding to inoculation with X. albilineans in SES208 and LA purple

To further analyze the possible functions of DEGs during X. albilineans infection in LA Purple and SES208 sugarcane, KEGG and GO analyses were performed. A total of 6098 DEGs in SES208 and LA Purple were significantly enriched in 119 KEGG pathways (Fig. 3a). By cross-comparison of all DEGs between SES208 and LA Purple, we found that these two cultivars shared 715 DEGs across the comparison, while 3804 showed specifically changed expression in SES208 and 1579 DEGs showed specifically changed expression in LA Purple (Fig. 3a). We also found that the overall expression of DEGs responding to X. albilineans was higher in SES208 than in LA Purple (Fig. 3b). KEGG analysis was performed on the core DEGs and cultivar-specific DEGs in SES208 and LA Purple. In SES208, we found that 3804 cultivar-specific DEGs were enriched in 110 pathways, particularly in plant-pathogen interaction (Fig. 3c). Notably, 58, 124, 80, and 28 genes were significantly enriched for pathways including MAPK signaling, plant-pathogen interactions, plant hormone-mediated signal transduction, and the biosynthesis of ubiquinone and other terpenoid quinones, respectively. Another 95, 71, 39, and 37 DEGs were mainly concentrated in functions such as the biosynthesis of phenylpropanoids, metabolism of starch and sucrose, metabolism of phenylalanine, and the biosynthesis of flavonoids, respectively (Fig. 3c). In LA Purple, 1579 cultivar-specific DEGs were also significantly enriched for terms in metabolic pathways and the biosynthesis of secondary metabolites (Fig. 3d). A total of 715 DEGs in common between SES208 and LA Purple were enriched in 76 pathways and 298 DEGs were enriched for terms related to metabolic pathways and secondary metabolite biosynthesis (Fig. 3e). Our GO term enrichment analysis showed that, compared to LA Purple, DEGs in SES208 were significantly enriched for terms in the biological processes domain including those related to responses to stresses or plant hormones (Additional file 1).

Fig. 3
figure 3

Analysis of differential gene expression in sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. a Venn diagrams of DEGs between SES208 with LA Purple. b Heat map of conserved core DEGs expressed in common between SES208 with LA Purple after inoculation with Xanthomonas albilineans. c-e Numbers of DEGs in major KEGG pathways after inoculation with Xanthomonas albilineans. c Specific DEGs in SES208. The red font represents the pathway specifically enriched in SES208; The blue font indicates that the number of DEGs in the enriched pathway in SES208 is greater than that in LA Purple. d Specific DEGs in LA Purple. e DEGs expressed in common between SES208 and LA Purple. Expression values were calculated as relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

WGCNA analysis of coexpressed gene modules responding to inoculation with X. albilineans in SES208 and LA purple

To identify the key regulatory genes related to X. albilineans resistance in SES208 and LA Purple, a WGCNA was performed. A total of 6098 DEGs in SES208 and LA Purple were chosen for constructing the network, and of the 18 modules with correlated expression detected by WGCNA, seven had correlation coefficients higher than 0.7 (p < 0.001) (Fig. 4a). DEGs in SES208 enriched in the MEblack module (r = 0.83) and the MEdarkolivegreen module (r = 0.87) were involved in plant hormone signaling, MAPK signaling, plant-pathogen interactions, biosynthesis of terpenoid quinones, metabolism of starch and sucrose, and phenylpropanoid biosynthesis (Additional file 2). These genes were then chosen to construct coexpression networks to identify candidate regulatory genes that respond to X. albilineans in SES208. Our co-expression network indicated 16 and 25 hub genes in the MEblack module and the MEdarkolivegreen module, respectively (Fig. 4b, c). These hub genes are enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction, all of which might be involved in regulating the responses of sugarcane to infection with X. albilineans (Fig. 4b, c).

Fig. 4
figure 4

Weighted gene co-expression network analysis (WGCNA) of sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. a Heat map of the correlation between WGCNA modules in SES208 and LA Purple. SES: SES208; LA: LA Purple. The left lane indicates 18 modules. The right lane indicates the module-trait correlation from − 1 to 1. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively. b The co-expression networks in SES208 in MEblack module. c The co-expression networks in SES208 in MEdarkolivegreen module. Light blue: plant pathogen interaction; Purple: MAPK signaling pathway; Red: plant hormone signal transduction; Orange: terpenoid quinone biosynthesis; Green: phenylpropanoid biosynthesis; Pink:starch and sucrose metabolism. The blue font indicates hub genes

Full size image

Allelic analysis of DEGs that respond to inoculation with X. albilineans in SES208 and LA purple

Allelic diversity including ASE (allele-specific expression) may be assessed in diploid or polyploid plants. ASE can lead to dramatic changes in gene expression profiles [32]. To identify ASE of defense response genes in sugarcane, we analyzed 345 DEGs from SES208 and LA Purple sugarcane that were identified in response to X. albilineans by KEGG pathway enrichment analysis (Fig. 5). According to the allele table https://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Saccharum%20spont%20alleleTable-Jan%202019.csv [4], all of these 345 DEGs have one to four alleles and include dispersed duplicated genes and tandem duplicated genes. These 345 DEGs include a total of 281 differentially expressed alleles, 57 dispersed duplicated genes, and 7 tandem duplicated genes (Fig. 5). Many of the DEGs that might be involved in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction exhibited decreased transcript abundance after X. albilineans infection, indicating that the expression of these DEGs was likely negatively regulated in response to pathogen invasion (Fig. 5a). However, the transcript abundances of most DEGs in metabolic pathways increased, indicating that the expression of these DEGs was likely positively regulated in response to pathogen invasion (Fig. 5b).

Fig. 5
figure 5

Alleles analysis of DEGs involved in major pathways with function annotation information and expression pattern between sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. a Alleles analysis of DEGs involved in MAPK signaling pathway, plant-pathogen interaction and plant hormone signal transduction. b Alleles analysis of DEGs in phenylpropanoid biosynthesis, phenylalanine metabolism, flavonoid biosynthesis and terpenoid quinone biosynthesis. The blue font indicates hub genes in Fig. 4. Genes shown in red font exhibited expression of one specific allele. The expression value was calculated by relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

Of these 345 DEGs, 21 exhibited expression of one specific allele (Fig. 5). We also examined the expression of differentially expressed alleles of DEGs with two, three, or four alleles in response to X. albilineans, and found that 64.86, 49.02, and 30.61% of the haplotypes were differentially expressed among these multi-allelic DEGs, respectively, indicating an allelic dominance effect in these two cultivars responding to X. albilineans treatment (Fig. 5; Additional file 3). After analyzing their transcript expression in SES208, we found that 85 out of 281 DEGs exhibited ASE, 42 out of these 85 DEGs expressed alleles with increased transcript abundance, and 43 out of these 85 DEGs exhibited decreased transcript abundance for one differentially expressed allele at 24 hpi (Additional file 4). It is noteworthy that we identified 62 genes with two differentially expressed alleles, 93 genes with three differentially expressed alleles, and 44 genes with four differentially expressed alleles (Addition file 3). The transcript abundances of alleles of about 37.2% of the genes with two to four alleles were strictly up-regulated at 24 hpi in response to X.albilineans in SES208 (Additional file 4). Thus, our results indicate that ASE is associated with resistance to infection by X. albilineans in SES208. Further, we discovered that the genes exhibiting ASE showed similar expression patterns in response to X. albilineans in SES208 (Additional file 4). Further, 19 DEGs showed increased transcript abundances and 38 DEGs showed decreased transcript abundances at 24 hpi in dispersed duplicated DEGs (Additional file 5). The transcript abundances of five tandem duplicated DEGs showed decreased transcript abundances at 24 hpi and two tandem duplicated DEG showed increased transcript abundances at 48 hpi in SES208 during X. albilineans infection (Additional file 5). These results indicated the expression of these alleles, dispersed duplicated genes, and tandem duplicated genes showed significant and consistent responses to infection with X. albilineans.

Potential functions of major metabolic pathways to X. albilineans infection in sugarcane

Our KEGG analysis of DEGs responding to inoculation with X. albilineans infection identified a total of 165 DEGs between SES208 and LA Purple that might participate in the biosynthesis of terpenoid quinones, phenylpropanoids, and flavonoids, or the metabolism of starch and sucrose or phenylalanine, among which 122 of 165 DEGs exhibited high transcript abundances and 43 out of 165 DEGs exhibited low transcript abundances in SES208 compared to LA Purple at 24 hpi in response to X. albilineans infection (Fig. 5b).

The transcript abundances of DEGs encoding the key enzymes involved in flavonoid biosynthesis, phenylalanine metabolism, and terpenoid quinone biosynthesis were mainly increased in SES208 compared to LA Purple in response to X. albilineans infection and the alleles of these DEGs showed the same expression patterns (Fig. 5b). Of the DEGs putatively involved in the phenylpropanoid biosynthesis pathway, the transcript abundances of 81% were increased at 24 hpi or 48 hpi in SES208 but were less highly expressed in LA Purple (Fig. 5b). The transcripts of all alleles of the beta-glucosidase(BglX) and Bright Trichomes 1 (BRT1) genes, which are both essential in phenylpropanoid synthesis, increased significantly after infection with X. albilineans in SES208 but not in LA Purple (Fig. 5B; Additional file 6). The CYP450 family is involved in the synthesis and metabolism of phenylpropanoids [33] and one of its members, CYP73, also showed increased transcript expression after 24 hpi in SES208 (Fig. 5B; Additional file 7). The CYP75 gene family plays an important role in flavonoid biosynthesis in plants [34]. Members of the CYP 98 family catalyze the 3-hydroxylation step in the phenylpropanoid pathway (add reference). Interestingly, all alleles of these two genes showed increased transcript abundance after 24 hpi in SES208 (Fig. 5B), indicating that their expression responds positively to infection with X. albilineans in SES208. However, the expression of these two genes showed little change in LA Purple (Fig. 5B), which could be consistent with the stronger resistance of SES208 than LA Purple to infection with X.albilineans.

DEGs associated with starch and sucrose metabolism pathways showed differential expression trends in SES208 and LA Purple after X. albilineans infection. For example, glucan endo-1,3-beta-glucosidase 4 (GN4) and glucose-1-phosphate adenylyltransferase (glgC) exhibited increased transcript expression in SES208 at 24 hpi (Additional file 7). However, glycogenin (GYG1) exhibited decreased transcript expression at 24 hpi (Additional file 6). The transcript abundances of otsB, which encodes trehalose 6-phosphate phosphatase, sharply decreased for five of eight alleles at 24 hpi in response to X. albilineans infection. The dynamic changes in the expression of the above DEGs in the starch and sucrose metabolism pathway might provide evidence of a direct relationship between higher X. albilineans resistance and lower biomass in SES208.

Potential functions of the MAPK signaling pathway and plant-pathogen interaction pathway during X. albilineans infection

Based on the above results, we identified candidate genes from the plant-pathogen interaction and MAPK signaling pathways whose expression responds to X. albilineans infection. In PTI, respiratory burst oxidase homologues (RBOHs) are a major source of reactive oxygen species (ROS) during plant-microbe interactions [35, 36]. We found that the Rboh gene exhibited significantly increased expression in response to inoculation with X. albilineans in SES208 (Fig. 6a). Further, transcript expression of the catalase1/ascorbate peroxidase1 (CAT1/APX1) genes related to cell death and H2O2 production increased, which might induce MAPK signaling to protect plant cells from oxidative stress and death. In the early stages of X. albilineans infection of sugarcane, the transcript expression of SsMEKK, SsMKK, and SsMPK decreased (Fig. 6b). The transcript expression of SsWRKY22 and SsWRKY33, TF-encoding core genes with two and three alleles, respectively, decreased sharply in response to infection with X. albilineans, suggesting that the decreased expression of these genes could be an important aspect of the sugarcane response to X. albilineans (Fig. 5a; Fig. 6b). The expression of genes encoding cyclic nucleotide gated channels (CNGCs) that are normally involved in the regulation of Ca2+ channels and stomata increased in SES208, while the expression of CNGCs showed no response in LA Purple (Fig. 6a). The expression of the majority of genes encoding CDPK and CML, except for Sspon.03G0026290-3P in SES208, dramatically decreased at 24 hpi (Fig. 6a).

Fig. 6
figure 6

Models for regulation of gene expression in plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction in sugarcane varieties SES208 and LA Purple after inoculation with Xanthomonas albilineans. a Plant-pathogen interaction. b MAPK signaling pathway. c Plant hormone signal transduction. Expression values were calculated by relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

The transcript abundance of the ethylene receptor ripening inhibitor (RIN) gene decreased in SES208 but not in LA Purple (Fig. 6a). The disease resistance genes RPM1 (RESISTANCE TO P. SYRINGAE PV MACULICOLA 1) and RPS2 (RIBOSOMAL PROTEIN S2), which exhibit ASE, decrease in transcript expression after X. albilineans infection in SES208. In contrast, downstream resistance genes involved in R gene-mediated plant disease resistance responses showed significantly decreased expression (Fig. 6a). Further, the 3-ketoacyl-CoA synthase (KCS) genes involved in the synthesis of ultra-long-chain fatty acids increased sharply at 24 hpi (Fig. 6a). These results indicate that differential expression of genes involved in the PTI and ETI signaling pathways is likely essential for X. albilineans response in SES208.

Potential functions of plant hormone signal transduction in resistance to X. albilineans infection

The plant hormones SA, JA, ET, and ABA act in pathways that protect plants against infection with pathogens [37]. NPR1-LIKE PROTEIN 3 (NPR3) is an important regulator of plant SAR signaling, and we found that the transcript expression of the NPR3 genes decreased sharply in SES208 and LA Purple at 24 hpi (Fig. 6c; Additional file 7). In SES208, the transglutaminase (TGA) gene downstream of the SA pathway exhibited significantly increased transcript expression that peaked at 48 hpi in SES208, but showed a slight decrease at 24 hpi in LA Purple (Fig. 6c; Additional file 7). The transcript abundance of the downstream gene PR1 (pathogenesis-related protein 1), a core R gene with ASE, increased in both SES208 and LA Purple in response to X. albilineans infection (Fig. 6c). The JAZ (jasmonate ZIM domain) gene functions during the regulation of responses involving the plant hormone JA [38]. At 48 hpi and 72 hpi, the transcript abundances of JAZ genes were significantly increased in SES208 (Fig. 6c). All four alleles of MYC2 showed decreased expression at 24 hpi in SES208 (Fig. 5a; Fig. 6c; Additional file 7).

In the ethylene signal transduction pathway, response to antagonist 1 (RAN1) helps ethylene receptor (ETR) to interact with ethylene [39]. After X. albilineans infection, the expression of ETR decreased in LA Purple, but increased in SES208 (Fig. 6c). Transcript expression of the EBF1 (EIN3-BINDING F BOX PROTEIN 1) gene decreased significantly in SES208, but there was no obvious change in the transcript abundance of this gene in LA Purple (Fig. 6c; Fig.8). The expression of ETHYLENE-INSENSITIVE3 (EIN3) was likely inhibited by EBF1 in SES208 while the expression of ETHYLENE RESPONSE FACTOR (ERF) genes increased in both SES208 and LA Purple (Fig. 6c). Transcripts encoding PYL (pyr1-like), an ABA receptor that is involved in the sensing and transduction of ABA signals [40] were expressed at lower levels after infection in SES208. However, genes encoding type 2C protein phosphatases (PP2C) were all expressed at higher levels at 24 hpi, and the downstream genes encoding SnRK2 (class III SNF-1-related protein kinase 2) exhibited increased expression of two genes and decreased expression of two other genes after inoculation of SES208 with X. albilineans (Additional file 7). Further, the expression of genes encoding abscisic acid-responsive TFs (ABF) increased and their transcripts exhibited increased abundance at 24 hpi in SES208 (Fig. 6c).

Discussion

Analysis of sugarcane cultivars transcriptomes during X. albilineans infection

We constructed 24 cDNA libraries to elucidate genetic responses that are triggered by X. albilineans infection in the sugarcane species S. spontaneum and S. officinarum. Comparative transcriptome and functional enrichment analyses of SES208 and LA Purple revealed sharp X. albilineans infection-induced gene expression changes in SES208 compared to LA Purple. These changes involved genes in key plant defense pathways such as plant-pathogen interaction, MAPK signaling, and plant hormone function pathways in addition to metabolic pathways for flavonoid, phenylpropanoid, and terpenoid quinone biosynthesis as well as for phenylalanine, starch, and sucrose metabolism.

Polyploidy results in changes in cellular structure, function [41, 42], and allelic diversity [43] that result in dramatic changes in gene expression profiles [44, 45], such as ASE [32]. Here, we identified a total of 345 DEGs involved in X. albilineans resistance in SES208 (Fig. 5). The hub DEGs identified were candidates that might be crucial for responses to X. albilineans infection in sugarcane. Among these 345 DEGs, 85 exhibited ASE while 196 showed balanced allelic expression (Additional file 4). We suggest that ASE in SES208 is largely a genotype-specific phenomenon and might be common in some individual wild species. Notably, the expression patterns of genes with multiple alleles were similar, which indicates that these alleles might act synergistically in resistance to X. albilineans. We predict that the genes showing ASE in SES208 reflect the breeding history of sugarcane, and might have undergone greater selection pressure. Furthermore, the occurrence of ASE among these defense genes might be associated with traits inherited from wild genotypes as discussed previously [46]. The higher frequency of ASE in defense-responsive genes is consistent with observations in wheat, in which many such genes exhibit homologous gene expression bias in plants attacked by pathogenic fungi [47]. Sugarcane hybridization and breeding have long focused mainly on selection for disease resistance and cultivars thus vary widely in their response to pathogen infection.

Transcriptional regulation is an important aspect of gene expression regulation in plants that controls many crucial biological processes via gene regulatory networks [48, 49]. Transcripts encoding five types of TFs including WRKY, NAC, MYB, ERF, and bHLH were enriched in SES208 in response to X. albilineans infection, which suggests that these TFs could be directly or indirectly involved in resistance to X. albilineans. Previous studies have indicated that WRKY TFs can interact with receptors for PTI and become involved in MAPK signaling to activate or regulate gene expression [50]. The decreased transcript expression of SsWRKYs in SES208 indicates that these genes might be involved in some regulatory responses during defense to X. albilineans (Fig. 2; Fig.8). NAC TFs ANAC019 and ANAC055 might be involved in plant defense responses mediated by JA and in regulation of the expression of the JA-induced defense genes VSP1 and LOX2 [51]. R2R3-MYB TFs such as AtMYB96 [50] and OsMYB4 [52] could increase P. syringae resistance and increase the expression of PR protein-encoding genes in Arabidopsis and rice, respectively. The transcript expression of about 90% of the SsNACs and all of the SsMYBs in SES208 increased, suggesting that these TFs might activate gene expression during resistance to X. albilineans (Fig. 2). ERF19 plays a role with NINJA in ERF19-mediated regulation of Arabidopsis innate immunity [53]. SsERF genes show ASE and all increase in transcript abundance in SES208, suggesting that they might activate gene expression to enhance resistance to X. albilineans (Fig. 2). The bHLH transcription factor GmPIB1 increases soybean resistance to P. sojae infection by inhibiting the expression of GmSPOD1 [54]. Further, the bHLH TF HBI1 negatively regulates some aspects of disease defense responses in Arabidopsis [55]. However, the transcript abundances of about 85% of SsbHLHs actually decreased in SES208, suggesting that the expression of the bHLH TFs themselves might be regulated negatively during responses to X. albilineans in sugarcane. This might also indicate that increased or decreased transcript expression of genes encoding different classes of TFs affect various resistance responses to X. albilineans in SES208 (Fig. 2).

Plant-pathogen interaction pathway responsive X. albilineans infection in sugarcane

Plant immunity requires the action of multiple receptors and ligands to provide adequate protection against pathogen attack. In the early stage of oxidative burst when plants are infected by pathogenic microorganisms, they sense PAMPs, abiotic stress, and extracellular signals through surface-localized PRRs [56]. Flagellin receptor FLS2 can recognize bacterial proteins [57] and elevate ROS production [58]. APX1 is a key factor in removing ROS produced in response to stresses and disease progression in plants [59]. In addition, CAT can also decompose H2O2 to generate O2, and trigger benzoic acid synthesis for generation of SA, leading to the elaboration of a systemic defense response (Fig. 7). CAT activity in transgenic tobacco can result in serious HR and effectively control bacterial infection [60, 61]. The expression of the SsCAT1 gene decreased continuously at 24 hpi upon X. albilineans infection (Fig. 8), and the transcript abundance of two SsAPX1 genes increased at 24 hpi in SES208 (Fig. 6a). These results demonstrate that SsWRKY genes might negatively regulate the expression of SsPR1 and plants might mainly rely on SsAPX1 to generate H2O2 in early stages of infection and then rely on SsCAT1 activated by the MAPK signaling pathway during later stages. CEBiP and CERK act as LysM-RP and LysM-RK proteins, respectively, and are also regarded as PRRs [62]. Interestingly, the response of CEBiP to the pathogen in SES208 and LA Purple was similar to that of OsLYP4 and OsLYP6 in response to X. oryzae pv. oryzae (PGNXoo), X. oryzae pv. oryzaecola (PGNXoc), and Pseudomonas syringae pv. tomato (Pto) DC3000 (PGNPto) in rice. OsLYP4 and OsLYP6 are normally complexed but can dissociate and when complexed with OsCERK can sense bacterial peptidoglycan (PGN). The susceptibility of transgenic plants was positively correlated with the transcript abundances of OsLYP4 and OsLYP6 genes [63]. At 24 hpi, the transcript abundance of CEBiP is lower in SES208 and higher in LA Purple (Additional file 7).

Fig. 7
figure 7

Hypothetical model of the molecular mechanism of the response to X. albilineans in sugarcane

Full size image
Fig. 8
figure 8

Verification by qRT-PCR of the expression patterns of four DEGs involved in different pathways during infection of sugarcane with Xanthomonas albilineans. The right y-axis indicates FPKM value; the left y-axis indicates relative expression; the x-axis indicates different time periods after infection. 0 h: uninfected plants; 24 h, 48 h, and 72 h represent 24, 48, and 72 hours post inoculation of sugarcane with X. albilineans, respectively

Full size image

As a cation channel protein, CNGC can participate in the regulation of Ca2+ influx to prevent pathogen invasion [64]. In SES208, the expression of CNGS transcripts showed an increasing trend. If CNGS participates in X. albilineans resistance in SES208, it might act to open calcium channels, promote Ca2+ influx, and activate intracellular disease resistance genes (Fig. 6a; Additional file 7). Calcium-dependent protein kinases (CDPKs) can be activated by Ca2+ [65] and overexpression of the CPK4 gene in tobacco could enhance the transcription of SAR-related genes, lead to HR, and improve the broad-spectrum disease resistance of plants [66]. In addition, the increase in calcium can directly or through CaM/CML, result in the phosphorylation or dephosphorylation of TFs [67]. In the early stages of infection in sugarcane, metabolites from pathogenic bacteria might inhibit the binding of SsCML and SsCDPK to calcium ions, resulting in decreased expression of SsCML and SsCDPK transcripts (Additional file 7). Expression of SsRboh transcripts increases, which might enhance production of reactive oxygen species (H2O2) (Fig. 6a). The processes above constitute our proposed model for a pathogen-associated molecular PTI response that might occur during the response of SES208 against pathogens.

When a pathogen suppresses or acts on PTI by secreting effector proteins into the host, R proteins can sense or recognize the effector protein, thereby triggering ETI [15]. During the ETI response, effectors (Avirulence, Avr) activate the expression of numerous genes encoding disease-resistance proteins that specifically recognize the effector and allow a resistance response [15]. The RPM1 protein can recognize the RIN4 protein that becomes phosphorylated due to the presence of the effectors AvrRPM1 and AvrB, which in turn results in the HR [68]. The RPS2 gene that functions in disease resistance in Arabidopsis was used to identify a virulence gene avrRpt2 and to select the mutants with altered resistance to Pseudomonas syringae, a bacterial plant pathogen [69]. In SES208, the transcript abundances of SsRIN4, SsRPM1, SsRPS2, and genes encoding R proteins significantly decreased within 0-24 hpi (Fig. 6a). It is possible that when plants are strongly attacked by pathogenic organisms, these genes cannot be expressed quickly enough to resist invasion by pathogens. The expression of these genes responds relatively weakly to pathogen invasion and a slight HR occurs. The protein encoded by XopD from X. campestris pv. campestris strain B100 targets the MYB30 TF in Arabidopsis that activates plant defense and cell death induced by bacteria [70]. In our analysis of TFs, we observed increased transcript expression for 15 genes encoding MYB TFs (Fig. 2; Fig. 6a). Transcript abundances of MYB genes and KSC both increased in SES208 (Fig. 6a), which suggest that MYBs might act as positive regulators of KCS to activate SES208 defenses against X. albilineans infection.

Possible roles of plant hormone-mediated signal transduction genes responsive to X. albilineans infection in sugarcane

Both PTI and ETI can cause increased biosynthesis of SA and defense phytohormones that mediate SAR (Fig. 7) [71]. The SA receptor NPR3 is known to promote the degradation of the transcriptional repressor JAZ, thereby activating JA synthesis and inducing JA-responsive gene expression to allow plants to resist pathogens [72]. In SES208, the expression of Sspon.01G0024890-2C (SsNPR3) decreased, which suggests that defense responses in SES208 might not depend on this signaling pathway (Fig. 6c; Additional file 7). MYC2 is a positive regulator in the JA signal pathway [73]. At 0 hpi, the abundances of MYC2 and JAZ transcripts were significantly higher in SES208 than in LA Purple (Fig. 8), but the expression of both these genes decreased in response to X. albilineans in SES208 (Fig. 6c; Additional file 7), suggesting that JAZ might interact with MYC to inhibit the transcription of JA-responsive genes. During infection of plants with viruses or other pathogens, HR responses are triggered by Avr/R protein interactions that can result in changes in the levels of disease defense-related hormones including ET, JA, NO, SA, and in the accumulation of ROS via MAPK signaling (Fig. 7) [74].

Some studies have shown that ABA can also have positive or negative regulatory effects during plant immunity [75, 76]. In our study, we regarded the SsPYL and SsPP2C genes as putative negative and positive regulators of responses to X. albilineans infection (Fig. 6c; Additional file 7). Therefore, we hypothesize that several members of the MAPK family might be activated directly or indirectly by the ABA core signaling pathway during response to X. albilineans infection (Fig. 7) [77]. The above results suggest that ETI and ABA might play necessary roles in response to X. albilineans infection in SES208.

Major metabolic pathways responsive to X. albilineans infection in sugarcane

Studies of different plant metabolites have shown the potential defensive activity of these compounds against plant pathogens [24]. The phenylpropanoid pathway is involved in the production of flavonoids, lignans, phytoalexins, phenolic acids, coumarins, and stilbene, which are directly and indirectly involved in plant development and disease responses [78]. Thus, altered levels of JA, SA, and secondary metabolites might influence HR and other disease resistance-related processes (Fig. 7) [79]. We found that the transcript expression of essential genes BglX and BRT1 encoding important enzymes in the phenylpropanoid biosynthesis pathway increased in SES208 during X. albilineans infection, which might have contributed to the defense of plants against this pathogen (Additional file 7). As important phenylpropanoid compounds, flavonoids and their derivatives have been shown to act as antimicrobial and signaling molecules during plant defense responses [80]. Enhanced levels of flavonoids in cotton suppress plant invasion by the fungal pathogen Verticillium dahliae [81]. The expression of flavonoid biosynthesis-related enzymes such as chalcone synthase (CHS), 5-O-(4-coumaroyl)-D-quinate 3′-monooxygenase (CYP98A/C3’H), flavonoid 3′-monooxygenase (CYP75B1) encoded by SsCHS, SsCYP98A, and SsCYP75B1, respectively, increased in SES208 after X. albilineans infection (Fig. 5b; Additional file 7), suggesting that increased flavonoid synthesis might be part of the defense response against X. albilineans infection in this sugarcane cultivar.

In maize, terpenoids represent another class of metabolites with disease defense functions in plants. Therefore, the higher expression of important genes in the terpenoid biosynthetic pathway and the increased accumulation of terpenoids have both been correlated with disease resistance [82, 83]. The transcript expression of four genes encoding 4-hydroxyphenylpyruvate dioxygenase (HPD) increased in SES208 after X. albilineans infection (Fig. 5b; Additional file 7), suggesting that the synthesis of terpenoids might aid defense against X. albilineans infection in SES208. Interestingly, many plant hormones, including cytokinins, ABA, and gibberellins are derived from terpenoids (Fig. 7) [84, 85].

During various metabolic processes, TFs can regulate the transcript abundances of genes that encode the regulatory proteins and structural enzymes that participate in the biosynthesis of secondary metabolites. Six major TF families including AP2/ERF, bHLH, MYB, NAC, WRKY, and bZIP are involved in the metabolism of terpenoids in plants [86], and the MYB-bHLH-WD40 complex regulates flavonoids bio-synthesis in plants [87, 88](Fig.7). Overall, the expression of DEGs encoding TFs and enzymes involved in primary and secondary metabolism and hormone signaling was more rapidly and strongly induced in SES208 than in LA Purple after inoculation with X. albilineans. Increased expression of transcripts related to the accumulation of metabolites also significantly increased in SES208, which is consistent with some previous findings that altered metabolite abundances can be associated with greater disease resistance in plants [89, 90]. These results suggest avenues for future research related to the effects of metabolic changes on the resistance of sugarcane to X. albilineans.

Conclusions

Infection with Xanthomonas albilineans can have destructive consequences for the vascular system of sugarcane plants. So far, the mechanisms by which sugarcane plants defend themselves against pathogens post inoculation have been unclear. Therefore, research into and cloning of resistance genes and breeding resistant varieties should provide effective methods for preventing infection with X. albilineans. The present study compared transcriptome profiles in the sugarcane cultivars SES208 and LA Purple before and after X. albilineans infection, revealing the apparent activation of an immune response to X. albilineans infection in SES208 but not in LA Purple. Genes showing changes in expression during X. albilineans infection were annotated with functions including plant-pathogen interactions, MAPK signaling, plant hormone signal transduction, and flavonoid and phenylpropanoid biosynthesis. The expression of many genes often associated with activated disease resistance decreased sharply after infection with X. albilineans in SES208, suggesting that its defense mechanisms are worthy of further study. Further, the involvement of putative core resistance genes that exhibit ASE was observed less frequently than that of multi-allelic genes with balanced allele expression, suggesting that the expression of genes with balanced allelic expression might allow synergistic maintenance of this resistance phenotype, especially as the existence of genes exhibiting ASE in S. saccharum appears to be a largely genotype-specific phenomenon.

In conclusion, the results of the present study provide useful knowledge for the improvement of resistance to X. albilineans in sugarcane and provide new insights into sugarcane responses against X. albilineans infection.

Methods

Sample collection and X. albilineans inoculation

For the present study, seedings of the sugarcane cultivars SES208 and LA Purple were grown in a standard greenhouse at the Center for Genomics and Biotechnology in Fujian Province with a 16-h light (28 °C) and 8-h dark (28 °C) cycle. A single colony of X. albilineans strain Xa-FJ1 was suspended in 1 mL of XAS liquid medium [91] and cultured at 28 °C with constant shaking at 200 rpm for 48 h. Then we added 1 μL of the latter culture to 40 mL of fresh XAS liquid medium and incubated at 28 °C for approximately 10 h. The bacterial suspension was diluted to 108 CFU/mL for plant inoculation. The leaves of twenty-day-old seedlings of SES208 and LA Purple with 3–5 leaves were inoculated with X. albilineans by cutting the leaf blades at mid-length with sterile scissors dipped in the bacterial suspension. Five control plants were treated in a similar way, but the bacterial suspension was replaced by sterile liquid medium XAS [29]. Leaf samples were collected from the two cultivars at 0 h post inoculation (hpi) hpi (named AP0 hpi and LA0 hpi), 24 hpi (named AP24 hpi and LA24 hpi), 48 hpi (named AP48 hpi and LA48 hpi) and 72 hpi (named AP72 hpi and LA72 hpi). Three biological and three technical replicates were collected for all treatments, immediately frozen in liquid nitrogen, and then stored at − 80 °C prior to RNA extraction.

RNA extraction and library construction

Total RNA was extracted using the TRIzol® kit (Invitrogen, Carlsbad, CA, USA). The integrity and quality of RNA were assessed using 1.0% agar gel electrophoresis and the NanoDrop 2000 (Thermo Fisher Scientific, USA). An NEBNext® Ultra™ RNA Library Prep Kit for Illumina® was used to construct the RNA-Seq libraries. The sample library concentrations were evaluated using the Qubit® 2.0 fluorimeter with the Qubit dsDNA HS Assay Kit before sending to Novogene (City, Country) for sequencing on the Illumina NovaSeq 6000 System.

Differential gene expression analysis, weighted correlation network analysis (WGCNA), and function enrichment analysis

The assembly and annotation of the sugarcane AP85-441 genome was downloaded from Ray Ming’s laboratory at (http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/) [92]. We conducted quantitative analysis of RNA-Seq data using Trinity Transcript Quantification. For detailed methods, please refer to the website: https://github.com/trinityrnaseq/trinityranseq/wiki/Trinity-Transcript-Quantification. The default comparison software Bowtie2 [93] and then RSEM [94] were used to calculate TPM (Transcripts Per Kilobase Million) values, and genes with an adjusted p-value < 0.001 (|log2(fold change)| ≥ 2) found using edgeR [95] were considered as DEGs.

Differentially expressed TFs were predicted using the PlantTFDB website (http://planttfdb.cbi.pku.edu.cn/prediction.php). Then WGCNA analysis was performed using the WGCNA (v1.70-3, https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/) package to find clusters (modules) of genes with highly correlated expression [96].

GO (Gene ontology) enrichment analysis [97] and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis [98] of the DEGs was performed using the EggNOG (v5.0) Database (http://eggnog5.embl.de/#/app/home) and OmicShare online tools (https://www.omicshare.com/).

Quantitative real-time PCR analysis

Each cDNA sample was synthesized using a Prime Script RT Reagent Kit (TaKaRa) and diluted with nuclease-free water to a concentration of 200 ng/μl. Primers to amplify genes chosen for quantitative real-time PCR (qRT-PCR) were designed by IDT (https://sg.idtdna.com/pages) (Additional file 8). Each qRT-PCR was conducted in a 12.5 μl reaction system with the following components: 6.25 μL SYBR Premix Ex TaqII (2 ×), 0.5 μL 0.4 μM forward primer, 0.5 μL 0.4 μM reverse primer, 0.5 μL diluted cDNA, and 4.75 μL ddH2O. Each qRT-PCR amplification was conducted on a Bio-Rad CFX96 Real Time PCR System (Bio-Rad, USA) with the following thermal cycling procedure: 95 °C for 3 min, 40 cycles of 95 °C for 10 s, 55 °C for 30 s. The temperature was increased by 0.5 °C every 5 s and the melting curve was generated from 65 °C to 95 °C. The transcript expression of each gene we analyzed was detected in three biological and three technical replications. The eukaryotic elongation factor 1a (eEF-1a) gene was used as an internal control for normalizing the expression data [99]. Relative gene expression was calculated using the 2−∆∆CT method [100].

Availability of data and materials

The raw sequencing reads of transcriptome data have been deposited to the NCBI Sequence Read Archive (SRA) as Bioproject PRJNA893415 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA893415).

Abbreviations

DEGs:

differentially expressed genes

TF:

transcription factor

ASE:

Allele-specific expression

LA:

Louisiana

HR:

hypersensitive response

PTI:

PAMP-triggered immunity

PAMPs:

pathogen-associated molecular patterns

ETI:

effector-triggered immunity

MAPK:

mitogen-activated protein kinase

SA:

salicylic acid

ABA:

abscisic acid

JA:

jasmonic acid

ET:

ethylene

SAR:

systemic acquired resistance

TMV:

tobacco mosaic virus

ISR:

induced systemic resistance

WGCNA:

weighted gene correlation network analysis

hpi:

hours post inoculation

PAL:

phenylalanine ammonia lyase

4CL:

4-coumarate:CoA ligase

F5H:

ferulate-5-hydroxylase

GN4:

glucan endo-1,3-beta-glucosidase 4

AMY:

alpha-amylase

INV:

beta-fructofuranosidase

ISA:

isoamylase

glgA:

starch synthase

glgC:

glucose-1-phosphate adenylyltransferase

GYG1:

glycogenin

E2.4.1.34:

1,3-beta-glucan synthase

PYG:

glycogen phosphorylase

HK:

hexokinase

PRRs:

pattern recognition receptors

FLS2:

FLAGELLIN SENSING2

RBOHs:

respiratory burst oxidase homologues

CAT1/APX1:

catalase1/ascorbate peroxidase1

CNGCs:

cyclic nucleotide gated channels

RIN:

ethylene receptor ripening inhibitor

RPM1:

RESISTANCE TO P. SYRINGAE PV MACULICOLA 1

RPS2:

RIBOSOMAL PROTEIN S2

KCS:

3-ketoacyl-CoA synthase

NPR1:

NONEXPRESSER OF PR GENES 1

NPR3:

NPR1-LIKE PROTEIN 3

TGA:

transglutaminase

PR1:

pathogenesis-related protein 1

JAZ:

jasmonate ZIM domain

RAN1:

response to antagonist 1

ETR:

ethylene receptor

CTR1:

CONSTITUTIVE TRIPLE RESPONSE 1

EBF1:

EIN3-BINDING F BOX PROTEIN 1

EIN3:

ETHYLENE-INSENSITIVE3

ERF:

ETHYLENE RESPONSE FACTOR

PP2C:

type 2C protein phosphatases

SnRK2:

class III SNF-1-related protein kinase 2

ABF:

abscisic acid-responsive transcription factor

CDPKs:

Calcium-dependent protein kinases

LRR:

leucine-rich repeat

CHS:

chalcone synthase

ANR:

anthocyanidin reductase

CYP98A/C3’H:

5-O-(4-coumaroyl)-D-quinate 3′-monooxygenase

CYP75B1:

flavonoid 3′-monooxygenase

HPD:

4-hydroxyphenylpyruvate dioxygenase

HPT:

homogentisate phytyltransferase

TAT:

tyrosine aminotransferase

TPS6:

TERPENE SYNTHASE 6

TPM:

Transcripts Per Kilobase Million

GO:

Gene Ontology.

References

  1. Botha FC, Moore PH. biomass and bioenergy. In: In: Sugarcane: Physiology, Biochemistry, and Functional Biology: John Wiley & Sons Ltd; 2013. p. 521–40.

    Chapter  Google Scholar 

  2. Chandra A. Physio-biochemical and molecular approaches to enhance sucrose content in sugarcane: Indian initiatives. Sugar Tech. 2011;13(4):315–21.

    Article  CAS  Google Scholar 

  3. Aguiar A, Milessi TS, Mulinari DR, Lopes MS, da Costa SM, Candido RG. Sugarcane straw as a potential second generation feedstock for biorefinery and white biotechnology applications. Biomass Bioenergy. 2021;144:105896.

    Article  CAS  Google Scholar 

  4. Zhang J, Zhang X, Tang H, Zhang Q, Hua X, Ma X, et al. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nat Genet. 2018;50(11):1565–73.

    Article  CAS  Google Scholar 

  5. D'Hont A, Grivet L, Feldmann P, Rao S, Berding N, Glaszmann JC. Characterisation of the double genome structure of modern sugarcane cultivars (Saccharum spp.) by molecular cytogenetics. Mol Gen Genet. 1996;250(4):405–13.

    Article  CAS  Google Scholar 

  6. O'Hara I, Mundree S. Sugarcane-based biofuels and bioproducts: John Wiley & Sons; 2016.

    Book  Google Scholar 

  7. Dillon SL, Shapter FM, Henry RJ, Cordeiro G, Izquierdo L, Lee LS. Domestication to crop improvement: genetic resources for Sorghum and saccharum (Andropogoneae). Ann Bot-London. 2007;100(5):975–89.

    Article  Google Scholar 

  8. Brandes EW, B. Sartoris G: sugarcane: its origin and improvement. Yearb US Dep Agric. 1936:561–624.

  9. Rott P, Bailey RA, Comstock JC, Croft BJ, Girard J. Saumtally AS: a guide to sugarcane diseases. D-CAS 1.2: a guide to sugarcane diseases. D-CAS. 2000;1:2.

    Google Scholar 

  10. Mensi I, Vernerey MS, Gargani D, Nicole M, Rott P. Breaking dogmas: the plant vascular pathogen Xanthomonas albilineans is able to invade non-vascular tissues despite its reduced genome. Open Biol. 2014;4(2):130116.

  11. Hashimi SM, Wall MK, Smith AB, Maxwell A, Birch RG. The phytotoxin albicidin is a novel inhibitor of DNA gyrase. Antimicrob Agents Ch. 2007;51(1):181–7.

    Article  CAS  Google Scholar 

  12. Johansson ON, Nilsson AK, Gustavsson MB, Backhaus T, Andersson MX, Ellerstrom M. A quick and robust method for quantification of the hypersensitive response in plants. PeerJ. 2015;3:e1469.

    Article  Google Scholar 

  13. Ding P, Ding Y. Stories of salicylic acid: a plant defense hormone. Trends Plant Sci. 2020;25(6):549–65.

    Article  CAS  Google Scholar 

  14. Han GZ. Origin and evolution of the plant immune system. New Phytol. 2019;222(1):70–83.

    Article  Google Scholar 

  15. Jones JD, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9.

    Article  CAS  Google Scholar 

  16. Chisholm ST, Coaker G, Day B, Staskawicz BJ. Host-microbe interactions: shaping the evolution of the plant immune response. Cell. 2006;124(4):803–14.

    Article  CAS  Google Scholar 

  17. Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, et al. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002;415(6875):977–83.

    Article  CAS  Google Scholar 

  18. Torres MA, Dangl JL, Jones JD. Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. Proc Natl Acad Sci U S A. 2002;99(1):517–22.

    Article  CAS  Google Scholar 

  19. Yu X, Feng BM, He P, Shan LB. From Chaos to harmony: responses and signaling upon microbial pattern recognition. Annu Rev Phytopathol. 2017;55:109–37.

    Article  CAS  Google Scholar 

  20. Yuan M, Ngou BPM, Ding P, Xin XF. PTI-ETI crosstalk: an integrative view of plant immunity. Curr Opin Plant Biol. 2021;62:102030.

    Article  CAS  Google Scholar 

  21. White RF. Acetylsalicylic acid (aspirin) induces resistance to tobacco mosaic virus in tobacco. Virology. 1979;99(2):410–2.

    Article  CAS  Google Scholar 

  22. Choudhary DK, Prakash A, Johri BN. Induced systemic resistance (ISR) in plants: mechanism of action. Indian J Microbiol. 2007;47(4):289–97.

    Article  CAS  Google Scholar 

  23. Piasecka A, Jedrzejczak-Rey N, Bednarek P. Secondary metabolites in plant innate immunity: conserved function of divergent chemicals. New Phytol. 2015;206(3):948–64.

    Article  Google Scholar 

  24. Zaynab M, Fatima M, Abbas S, Sharif Y, Umair M, Zafar MH, et al. Role of secondary metabolites in plant defense against pathogens. Microb Pathog. 2018;124:198–202.

    Article  CAS  Google Scholar 

  25. Zhang F, Huang LY, Zhang F, Ali J, Cruz CV, Zhuo DL, et al. Comparative transcriptome profiling of a rice line carrying Xa39 and its parents triggered by Xanthomonas oryzae pv. Oryzae provides novel insights into the broad-spectrum hypersensitive response. BMC Genomics. 2015;16(1):111.

  26. Tariq R, Wang C, Qin T, Xu F, Tang Y, Gao Y, et al. Comparative transcriptome profiling of Rice near-isogenic line carrying Xa23 under infection of Xanthomonas oryzae pv. Oryzae. Int J Mol Sci. 2018;19(3):717.

  27. Meng JY, Ntambo MS, Luo LM, Huang MT, Fu HY, Gao SJ. Molecular identification of Xanthomonas albilineans infecting elephant grass (Pennisetum purpureum) in China. Crop Prot. 2019;124:353.

  28. Garcia-Seco D, Chiapello M, Bracale M, Pesce C, Bagnaresi P, Dubois E, et al. Transcriptome and proteome analysis reveal new insight into proximal and distal responses of wheat to foliar infection by Xanthomonas translucens. Sci Rep-Uk. 2017;7(1):10157.

  29. Ntambo MS, Meng JY, Rott PC, Henry RJ, Zhang HL, Gao SJ. Comparative transcriptome profiling of resistant and susceptible sugarcane cultivars in response to infection by Xanthomonas albilineans. Int J Mol Sci. 2019;20(24):6138.

    Article  CAS  Google Scholar 

  30. Meng JY, Ntambo MS, Rott PC, Fu HY, Huang MT, Zhang HL, et al. Identification of differentially expressed proteins in sugarcane in response to infection by Xanthomonas albilineans using iTRAQ quantitative proteomics. Microorganisms. 2020;8(1):76.

    Article  CAS  Google Scholar 

  31. Wang HB, Xiao NY, Wang YJ, Guo JL, Zhang JS. Establishment of a qualitative PCR assay for the detection of Xanthomonas albilineans (Ashby) Dowson in sugarcane. Crop Prot. 2020;130:105053.

    Article  CAS  Google Scholar 

  32. Gaur U, Li K, Mei S, Liu G. Research progress in allele-specific expression and its regulatory mechanisms. J Appl Genet. 2013;54(3):271–83.

    Article  CAS  Google Scholar 

  33. Werck-Reichhart D. Cytochromes P450 in phenylpropanoid metabolism. Drug Metabol Drug Interact. 1995;12(3-4):221–43.

    Article  CAS  Google Scholar 

  34. Nelson D, Werck-Reichhart D. A P450-centric view of plant evolution. Plant J. 2011;66(1):194–211.

    Article  CAS  Google Scholar 

  35. Torres MA, Onouchi H, Hamada S, Machida C, Hammond-Kosack KE, Jones JD. Six Arabidopsis thaliana homologues of the human respiratory burst oxidase (gp91phox). Plant J. 1998;14(3):365–70.

    Article  CAS  Google Scholar 

  36. Noirot E, Der C, Lherminier J, Robert F, Moricova P, Kieu K, et al. Dynamic changes in the subcellular distribution of the tobacco ROS-producing enzyme RBOHD in response to the oomycete elicitor cryptogein. J Exp Bot. 2014;65(17):5011–22.

    Article  CAS  Google Scholar 

  37. Alazem M, Lin NS. Roles of plant hormones in the regulation of host-virus interactions. Mol Plant Pathol. 2015;16(5):529–40.

    Article  CAS  Google Scholar 

  38. Pauwels L, Goossens A. The JAZ proteins: a crucial interface in the jasmonate signaling cascade. Plant Cell. 2011;23(9):3089–100.

    Article  CAS  Google Scholar 

  39. Binder BM. Ethylene signaling in plants. J Biol Chem. 2020;295(22):7710–25.

    Article  CAS  Google Scholar 

  40. Dittrich M, Mueller HM, Bauer H, Peirats-Llobet M, Rodriguez PL, Geilfus CM, et al. The role of Arabidopsis ABA receptors from the PYR/PYL/RCAR family in stomatal acclimation and closure signal integration. Nat Plants. 2019;5(9):1002–11.

    Article  CAS  Google Scholar 

  41. Comai L. Genetic and epigenetic interactions in allopolyploid plants. Plant Mol Biol. 2000;43(2-3):387–99.

    Article  CAS  Google Scholar 

  42. Renny-Byfield S, Wendel JF. Doubling down on genomes: polyploidy and crop plants. Am J Bot. 2014;101(10):1711–25.

    Article  Google Scholar 

  43. Yant L, Bomblies K. Genome management and mismanagement--cell-level opportunities and challenges of whole-genome duplication. Genes Dev. 2015;29(23):2405–19.

    Article  CAS  Google Scholar 

  44. Feldman M, Levy AA. Genome evolution in allopolyploid wheat--a revolutionary reprogramming followed by gradual changes. J Genet Genomics. 2009;36(9):511–8.

    Article  CAS  Google Scholar 

  45. Chen ZJ, Ni Z. Mechanisms of genomic rearrangements and gene expression changes in plant polyploids. Bioessays. 2006;28(3):240–52.

    Article  Google Scholar 

  46. Correr FH, Furtado A, Garcia AAF, Henry RJ, Margarido GRA. Allele expression biases in mixed-ploid sugarcane accessions. Sci Rep. 2022;12(1):8778.

    Article  CAS  Google Scholar 

  47. Powell JJ, Fitzgerald TL, Stiller J, Berkman PJ, Gardiner DM, Manners JM, et al. The defence-associated transcriptome of hexaploid wheat displays homoeolog expression and induction bias. Plant Biotechnol J. 2017;15(4):533–43.

    Article  CAS  Google Scholar 

  48. Valliyodan B, Nguyen HT. Understanding regulatory networks and engineering for enhanced drought tolerance in plants. Curr Opin Plant Biol. 2006;9(2):189–95.

    Article  CAS  Google Scholar 

  49. Hobert O. Gene regulation by transcription factors and microRNAs. Science. 2008;319(5871):1785–6.

    Article  CAS  Google Scholar 

  50. Liu YL, Schiff M, Dinesh-Kumar SP. Involvement of MEK1 MAPKK, NTF6 MAPK, WRKY/MYB transcription factors, COI1 and CTR1 in N-mediated resistance to tobacco mosaic virus. Plant J. 2004;38(5):800–9.

    Article  CAS  Google Scholar 

  51. Bu Q, Jiang H, Li CB, Zhai Q, Zhang J, Wu X, et al. Role of the Arabidopsis thaliana NAC transcription factors ANAC019 and ANAC055 in regulating jasmonic acid-signaled defense responses. Cell Res. 2008;18(7):756–67.

    Article  CAS  Google Scholar 

  52. Vannini C, Iriti M, Bracale M, Locatelli F, Faoro F, Croce P, et al. The ectopic expression of the rice Osmyb4 gene in Arabidopsis increases tolerance to abiotic, environmental and biotic stresses. Physiol Mol Plant P. 2006;69(1-3):26–42.

    Article  CAS  Google Scholar 

  53. Huang PY, Zhang JS, Jiang BE, Chan C, Yu JH, Lu YP, et al. NINJA-associated ERF19 negatively regulates Arabidopsis pattern-triggered immunity. J Exp Bot. 2019;70(3):1033–47.

    Article  CAS  Google Scholar 

  54. Cheng Q, Dong LD, Gao TJ, Liu TF, Li NH, Wang L, et al. The bHLH transcription factor GmPIB1 facilitates resistance to Phytophthora sojae in Glycine max. J Exp Bot. 2018;69(10):2527–41.

    Article  CAS  Google Scholar 

  55. Fan M, Bai MY, Kim JG, Wang TN, Oh E, Chen L, et al. The bHLH transcription factor HBI1 mediates the trade-off between growth and pathogen-associated molecular pattern-triggered immunity in Arabidopsis. Plant Cell. 2014;26(2):828–41.

    Article  CAS  Google Scholar 

  56. Nejat N, Mantri N. Plant immune system: crosstalk between responses to biotic and abiotic stresses the missing link in understanding plant Defence. Curr Issues Mol Biol. 2017;23:1–15.

    Article  Google Scholar 

  57. Wang D, Liang X, Bao Y, Yang S, Zhang X, Yu H, et al. A malectin-like receptor kinase regulates cell death and pattern-triggered immunity in soybean. EMBO Rep. 2020;21(11):e50442.

    Article  CAS  Google Scholar 

  58. Zou Y, Wang S, Zhou Y, Bai J, Huang G, Liu X, et al. Transcriptional regulation of the immune receptor FLS2 controls the ontogeny of plant innate immunity. Plant Cell. 2018;30(11):2779–94.

    Article  CAS  Google Scholar 

  59. Davletova S, Rizhsky L, Liang HJ, Zhong SQ, Oliver DJ, Coutu J, et al. Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. Plant Cell. 2005;17(1):268–81.

    Article  CAS  Google Scholar 

  60. Gaffney T, Friedrich L, Vernooij B, Negrotto D, Nye G, Uknes S, et al. Requirement of salicylic acid for the induction of systemic acquired resistance. Science. 1993;261(5122):754–6.

    Article  CAS  Google Scholar 

  61. Chen Z, Silva H, Klessig DF. Active oxygen species in the induction of plant systemic acquired resistance by salicylic acid. Science. 1993;262(5141):1883–6.

    Article  CAS  Google Scholar 

  62. Tang DZ, Wang GX, Zhou JM. Receptor kinases in plant-pathogen interactions: more than pattern recognition. Plant Cell. 2017;29(4):618–37.

    Article  CAS  Google Scholar 

  63. Liu B, Li JF, Ao Y, Qu J, Li Z, Su J, et al. Lysin motif-containing proteins LYP4 and LYP6 play dual roles in peptidoglycan and chitin perception in rice innate immunity. Plant Cell. 2012;24(8):3406–19.

    Article  CAS  Google Scholar 

  64. Hetherington AM, Brownlee C. The generation of Ca2+ signals in plants. Annu Rev Plant Biol. 2004;55:401–27.

    Article  CAS  Google Scholar 

  65. Xu W, Huang W. Calcium-dependent protein kinases in Phytohormone signaling pathways. Int J Mol Sci. 2017;18(11):2436.

  66. Heo WD, Lee SH, Kim MC, Kim JC, Chung WS, Chun HJ, et al. Involvement of specific calmodulin isoforms in salicylic acid-independent activation of plant disease resistance responses. Proc Natl Acad Sci U S A. 1999;96(2):766–71.

    Article  CAS  Google Scholar 

  67. Reddy ASN, Ali GS, Celesnik H, Day IS. Coping with stresses: roles of calcium- and calcium/calmodulin-regulated gene expression. Plant Cell. 2011;23(6):2010–32.

    Article  CAS  Google Scholar 

  68. Mackey D, Holt BF 3rd, Wiig A, Dangl JL. RIN4 interacts with pseudomonas syringae type III effector molecules and is required for RPM1-mediated resistance in Arabidopsis. Cell. 2002;108(6):743–54.

    Article  CAS  Google Scholar 

  69. Kunkel BN, Bent AF, Dahlbeck D, Innes RW, Staskawicz BJ. RPS2, an Arabidopsis disease resistance locus specifying recognition of pseudomonas syringae strains expressing the avirulence gene avrRpt2. Plant Cell. 1993;5(8):865–75.

    CAS  Google Scholar 

  70. Canonne J, Marino D, Jauneau A, Pouzet C, Briere C, Roby D, et al. The Xanthomonas type III effector XopD targets the Arabidopsis transcription factor MYB30 to suppress plant defense. Plant Cell. 2011;23(9):3498–511.

    Article  CAS  Google Scholar 

  71. Johnson C, Boden E, Arias J. Salicylic acid and NPR1 induce the recruitment of trans-activating TGA factors to a defense gene promoter in Arabidopsis. Plant Cell. 2003;15(8):1846–58.

    Article  CAS  Google Scholar 

  72. Liu LJ, Sonbol FM, Huot B, Gu YN, Withers J, Mwimba M, et al. Salicylic acid receptors activate jasmonic acid signalling through a non-canonical pathway to promote effector-triggered immunity. Nat Commun. 2016;7:13099.

  73. Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007;19(7):2225–45.

    Article  CAS  Google Scholar 

  74. Mandadi KK, Scholthof KB. Plant immune responses against viruses: how does a virus cause disease? Plant Cell. 2013;25(5):1489–505.

    Article  CAS  Google Scholar 

  75. Adie BAT, Perez-Perez J, Perez-Perez MM, Godoy M, Sanchez-Serrano JJ, Schmelz EA, et al. ABA is an essential signal for plant resistance to pathogens affecting JA biosynthesis and the activation of defenses in Arabidopsis. Plant Cell. 2007;19(5):1665–81.

    Article  CAS  Google Scholar 

  76. Song W, Ma X, Tan H, Zhou J. Abscisic acid enhances resistance to Alternaria solani in tomato seedlings. Plant Physiol Biochem. 2011;49(7):693–700.

    Article  CAS  Google Scholar 

  77. de Zelicourt A, Colcombet J, Hirt H. The role of MAPK modules and ABA during abiotic stress signaling. Trends Plant Sci. 2016;21(8):677–85.

    Article  Google Scholar 

  78. Bellincampi D, Cervone F, Lionetti V. Plant cell wall dynamics and wall-related susceptibility in plant-pathogen interactions. Front Plant Sci. 2014;5:228.

    Article  Google Scholar 

  79. Li J, Chen M, Fan T, Mu X, Gao J, Wang Y, et al. Underlying mechanism of accelerated cell death and multiple disease resistance in a maize lethal leaf spot 1 allele. J Exp Bot. 2022;73(12):3991–4007.

    Article  CAS  Google Scholar 

  80. Ververidis F, Trantas E, Douglas C, Vollmer G, Kretzschmar G, Panopoulos N. Biotechnology of flavonoids and other phenylpropanoid-derived natural products. Part I: chemical diversity, impacts on plant biology and human health. Biotechnol J. 2007;2(10):1214–34.

    Article  CAS  Google Scholar 

  81. Long L, Liu J, Gao Y, Xu FC, Zhao JR, Li B, et al. Flavonoid accumulation in spontaneous cotton mutant results in red coloration and enhanced disease resistance. Plant Physiol Biochem. 2019;143:40–9.

    Article  CAS  Google Scholar 

  82. Christensen SA, Huffaker A, Sims J, Hunter CT, Block A, Vaughan MM, et al. Fungal and herbivore elicitation of the novel maize sesquiterpenoid, zealexin A4, is attenuated by elevated CO2. Planta. 2018;247(4):863–73.

    Article  CAS  Google Scholar 

  83. Huffaker A, Kaplan F, Vaughan MM, Dafoe NJ, Ni X, Rocca JR, et al. Novel acidic sesquiterpenoids constitute a dominant class of pathogen-induced phytoalexins in maize. Plant Physiol. 2011;156(4):2082–97.

    Article  CAS  Google Scholar 

  84. Chen F, Ludwiczuk A, Wei G, Chen XL, Crandall-Stotler B, Bowman JL. Terpenoid secondary metabolites in bryophytes: chemical diversity, biosynthesis and biological functions. Crit Rev Plant Sci. 2018;37(2-3):210–31.

    Article  Google Scholar 

  85. Amini H, Naghavi MR, Shen T, Wang Y, Nasiri J, Khan IA, et al. Tissue-Specific Transcriptome Analysis Reveals Candidate Genes for Terpenoid and Phenylpropanoid Metabolism in the Medicinal Plant Ferula assafoetida. G3 (Bethesda). 2019;9(3):807–16.

    Article  CAS  Google Scholar 

  86. Xu Y, Zhu C, Xu C, Sun J, Grierson D, Zhang B, et al. Integration of metabolite profiling and transcriptome analysis reveals genes related to volatile Terpenoid metabolism in finger citron (C. medica var. sarcodactylis). Molecules. 2019;24(14):2564.

  87. Zhao J. Flavonoid transport mechanisms: how to go, and with whom. Trends Plant Sci. 2015;20(9):576–85.

    Article  CAS  Google Scholar 

  88. Xie Y, Tan H, Ma Z, Huang J. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressors of the MYB/bHLH/WD40 complex in Arabidopsis thaliana. Mol Plant. 2016;9(5):711–21.

    Article  CAS  Google Scholar 

  89. Mu X, Li J, Dai Z, Xu L, Fan T, Jing T, et al. Commonly and specifically activated defense responses in maize disease lesion mimic mutants revealed by integrated transcriptomics and metabolomics analysis. Front Plant Sci. 2021;12:638792.

    Article  Google Scholar 

  90. Konig S, Feussner K, Kaever A, Landesfeind M, Thurow C, Karlovsky P, et al. Soluble phenylpropanoids are involved in the defense response of Arabidopsis against Verticillium longisporum. New Phytol. 2014;202(3):823–37.

    Article  Google Scholar 

  91. Lin LH, Ntambo MS, Rott PC, Wang QN, Lin YH, Fu HY, et al. Molecular detection and prevalence of Xanthomonas albilineans, the causal agent of sugarcane leaf scald, in China. Crop Prot. 2018;109:17–23.

    Article  CAS  Google Scholar 

  92. Zhang JS, Zhang XT, Tang HB, Zhang Q, Hua XT, Ma XK, et al. Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nat Genet. 2018;50(12):1565–73.

    Article  CAS  Google Scholar 

  93. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  CAS  Google Scholar 

  94. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12(1):1–6.

    Article  Google Scholar 

  95. Chen Y, Lun AT, Smyth GK. From reads to genes to pathways: differential expression analysis of RNA-Seq experiments using Rsubread and the edgeR quasi-likelihood pipeline. F1000Res. 2016;5:1438.

    Google Scholar 

  96. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC bioinformatics. 2008;9:559.

    Article  Google Scholar 

  97. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.

    Article  Google Scholar 

  98. Kanehisa M. KEGG bioinformatics resource for plant genomics and metabolomics. Methods Mol Biol. 2016;1374:55–70.

    Article  CAS  Google Scholar 

  99. Ling H, Wu Q, Guo J, Xu L, Que Y. Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative rt-PCR. PLoS One. 2014;9(5):e97469.

    Article  Google Scholar 

  100. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

Xanthomonas albilineans was provided by Professor Gao Sanji of Fujian Agriculture and Forestry University.

Funding

This project was supported by a startup fund from Fujian Agriculture and Forestry University to R.M.

Author information

Author notes

  1. Yaying Ma and Hongying Yu contributed equally to this work.

Authors and Affiliations

  1. Center for Genomics and Biotechnology, Fujian Provincial Key Laboratory of Haixia Applied Plant Systems Biology, Key Laboratory of Genetics, Breeding and Multiple Utilization of Crops, Ministry of Education, Fujian Agriculture and Forestry University, Fuzhou, 350002, Fujian, China

    Yaying Ma & Yijing Lu

  2. College of Agriculture, Fujian Agriculture and Forestry University, Fuzhou, 350002, China

    Yaying Ma, Yijing Lu, Mahpara Fatima, Ray Ming & Jingjing Yue

  3. Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, 518120, China

    Hongying Yu

  4. National Engineering Research Center for Sugarcane, Fujian Agriculture and Forestry University, Fuzhou, 350002, China

    Sanji Gao

Authors
  1. Yaying Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Hongying Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yijing Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Sanji Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mahpara Fatima
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ray Ming
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jingjing Yue
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YYM, HYY, RM, and JJY designed the experiments and performed the analyses. YYM, HYY, and YJL performed the experiments. YYM designed the figures and tables. YYM, HYY, and JJY wrote the manuscript. SJG provided the X. albilineans strain. JJY and RM supervised the experiments and revised the manuscript. MF revised the manuscript. All authors have approved the final article.

Corresponding authors

Correspondence to Ray Ming or Jingjing Yue.

Ethics declarations

Ethics approval and consent to participate

This project uses plant materials and complies with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. 

Gene Ontology (GO) analysis of DEGs in SES208 and LA Purple infected by X. albilineans.

Additional file 2. 

Number of DEGs in major KEGG pathways in MEblack and MEdarkolivegreen module after inoculation with Xanthomonas albilineans.

Additional file 3.

Proportion of alleles in differential genes.

Additional file 4. 

Heatmaps of differentially expressed alleles involved in MAPK signaling pathway, plant-pathogen interaction, plant hormone signal transduction, phenylpropanoid biosynthesis, phenylalanine metabolism, starch and sucrose metabolism, flavonoid biosynthesis and terpenoid-quinone biosynthesis. The expression value was calculated by relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h,and 72 h represent 24, 48 and 72 hours post inoculation of sugarcane with X. albilineans, respectively.

Additional file 5. 

Heatmaps of dispersely and tandem duplicated DEGs involved in MAPK signaling pathway, plant-pathogen interaction, plant hormone signal transduction, phenylpropanoid biosynthesis, starch and sucrose metabolism, flavonoid biosynthesis and terpenoid-quinone biosynthesis. The expression value was calculated by relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h,and 72 h represent 24, 48 and 72 hours post inoculation of sugarcane with X. albilineans, respectively.

Additional file 6. 

Heatmaps of DEGs involved in phenylpropanoid biosynthesis and starch and sucrose metabolism. The expression value was calculated by relative expression of genes (log2 Fold Change from − 2 to + 2). SES: SES208; LA: LA Purple. 0 h: uninfected plants; 24 h, 48 h,and 72 h represent 24, 48 and 72 hours post inoculation of sugarcane with X. albilineans, respectively.

Additional file 7. 

Significant differentially expressed genes (DEGs) information involved in sugarcane response to X. albilineans infection.

Additional file 8.

Primes sequences used in qRT-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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, Y., Yu, H., Lu, Y. et al. Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection. BMC Plant Biol 23, 52 (2023). https://doi.org/10.1186/s12870-023-04073-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04073-6

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us