Skip to main content
Download PDF

A newly identified glycosyltransferase AsRCOM provides resistance to purple curl leaf disease in agave

BMC Genomics volume 24, Article number: 669 (2023) Cite this article

Abstract

Background

Purple curl leaf disease brings a significant threat to the development of agave industry, the underlying mechanism of disease-resistant Agave sisalana. hybrid 11648 (A. H11648R) is still unknown.

Results

To excavate the crucial disease-resistant genes against purple curl leaf disease, we performed an RNA-seq analysis for A.H11648R and A.H11648 during different stages of purple curl leaf disease. The DEGs (differentially expressed genes) were mainly enriched in linolenic acid metabolism, starch and sucrose mechanism, phenylpropanoid biosynthesis, hypersensitive response (HR) and systemic acquired resistance. Further analysis suggested that eight candidate genes (4’OMT2, ACLY, NCS1, GTE10, SMO2, FLS2, SQE1 and RCOM) identified by WGCNA (weighted gene co-expression network analysis) may mediate the resistance to agave purple curl disease by participating the biosynthesis of benzylisoquinoline alkaloids, steroid, sterols and flavonoids, and the regulation of plant innate immunity and systemic acquired resistance. After qPCR verification, we found that AsRCOM, coding a glycosyltransferase and relevant to the regulation of plant innate immunity and systemic acquired resistance, may be the most critical disease-resistant gene. Finally, the overexpression of AsRCOM gene in agave could significantly enhance the resistance to purple curl disease with abundant reactive oxygen species (ROS) accumulations.

Conclusions

Integrative RNA-seq analysis found that HR may be an important pathway affecting the resistance to purple curl leaf disease in agave, and identified glycosyltransferase AsRCOM as the crucial gene that could significantly enhance the resistance to purple curl leaf disease in agave, with obvious ROS accumulations.

Peer Review reports

Background

Growing primarily in semi-arid regions, Agave sisalana. hybrid 11648 (A. H11648) is a globally cultivated plant known for its tough natural fibers [1]. Besides its use in producing natural fibers, A. H11648 could also be utilized for alcoholic beverages [2], nutraceuticals [3], and nanocellulose [4], making it a highly valuable crop. Despite its widely valuable usage, A. H11648 also faces a number of challenges due to its susceptibility to various biotic and abiotic stresses, such as cold temperatures [5], root rot disease [6], high temperatures [7], low rainfall [8], low aboveground biomass coverage [9] and agave purple curl leaf disease [10]. Among these various stresses, agave purple curl leaf disease poses a significant threat to the agave industry as it severely damages the plant’s leaves which are the primary source of its natural fibers and will cause a sharp reduction in leaves production by more than 30% in severe cases [11].

Agave purple curl leaf disease is spread through Dysmicoccus neobrevipes and has become a significant constraint on the healthy development of agave industry. When D. neobrevipes infects the leaves with its piercing-sucking mouthparts, the phytoplasma that it carries enters into the plant, leading to agave purple curl leaf disease and severe damage to plant leaves [10]. Once agave purple curl leaf disease was found on a few plants in one plantation, it almost spread to the whole plantation within approximately two months [12]. Although pesticide application is currently the most effective solution for preventing the disease, it is not sustainable or advisable due to environmental and economic price. In our previous researches, the disease-resistant (A. H11648R) plants for agave purple curl leaf disease were found in A. H11648 gardens that had suffered from this disease for many consecutive years. Then, planting A. H11648R agave is the more prevalent solution. However, the ability of A. H11648R agave against purple curl leaf disease weakens over time. Therefore, it is very urgent and critical to elucidate the underlying mechanism of A. H11648R agave against purple curl leaf disease and excavate the pivotal disease-resistant genes to cultivate stable disease-resistant cultivars.

Currently, research on plant defense mechanisms against phytoplasma primarily focuses on three areas. The first area involves the external application of plant resistance elicitors, such as benzothiadiazole (BTH) which could significantly reduce the infection of phytoplasma by inducing systemic acquired resistance [13]. The second area is the study of phytoplasma effector proteins, which are secreted by phytoplasma and could affect the growth and development of plants [14, 15]. Lastly, it is research on disease-resistant genes, which enable plants to inhibit the symptoms resulting from phytoplasma infection [16]. Among these three areas, transcriptome sequencing (RNA-seq) has been proved as an effective data mining method and performed to identify a series of crucial phytoplasma effector proteins and resistance genes in a variety of plants. For instance, combined transcriptome and metabolome analysis of Nerium indicum L. found that the genes enriched in MAPK-signaling (plant), plant-pathogen interaction, plant-hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, linoleic acid and α-linoleic acid metabolism pathways may play an important role in defensing against the phytoplasma infection [17]. Comparative transcriptome analysis showed that Calmodulin-like (CML) genes and cinnamoyl-CoA reductase-like SNL6 gene may be the key genes against phytoplasma infection in Indian jujube ‘Cuimi’ [18]. The Mu-GsSRK gene identified by transcriptome and DNA methylome could enhance transgenic plant resistance to the phytoplasma in mulberry [19]. The major latex protein-like 329 gene identified by transcriptome and proteome analysis demonstrated enhanced plant resistance to phytoplasma through altering flavonoid content in mulberry [20]. Until now, the relevant studies on agave purple curl leaf disease are very rare and its resistant mechanism to phytoplasma is still unknown. Therefore, RNA-seq could be an excellent strategy to decipher the underlying mechanism of disease-resistant agave against purple curl leaf disease and helpful for breeding stable and disease-resistant agave cultivars.

In this study, we conducted a transcriptome profile analysis for the disease-susceptible A. H11648 and its disease-resistant mutant A. H11648R during different processes of purple curl leaf disease. Our analysis revealed that the systemic acquired resistance accompanied with hypersensitive response (HR) may be the primary contributor to increase the resistance to purple curl leaf disease in agave. Moreover, we identified AsRCOM as the potential resistant gene by weighted gene co-expression network analysis (WGCNA), and its overexpression in agave could significantly increase the resistance to purple curl leaf disease. This study provides a comprehensive understanding of the resistance to purple curl leaf disease in agave and could be helpful for the development of new disease-resistant and stable agave germplasm.

Results

RNA-seq profiles of A. H11648 and A. H11648R plants during different processes of agave purple curl leaf disease

In this study, RNA-seq analysis was conducted to better understand the mechanism of agave purple curl leaf disease infection during different stages of the disease. Two agave materials were selected for experimental treatment: A. H11648R and A. H11648 plants for purple curl leaf disease (Fig. 1A). A. H11648 showed yellow blade tips at 60 d, and blackened and died blade tips at 90 d, while A. H11648R showed no obvious disease symptoms throughout the infection stages (Fig. 1A). The leaves from three stages in purple curl leaf disease were selected for RNA-seq analysis, which were early (0 day, without inoculation, control), mild (60 day post inoculation, where the blade tips of the A. H11648 agave showed obvious yellow spots) and severe (90 day post inoculation, where the blade tips of the A. H11648 agave turned black and died, and the neighboring leaf was sampled). Three biological replicates were conducted, and the A. H11648R agave showed no symptoms of disease throughout all the D. neobrevipes infection stages.

A total of 18 RNA-seq libraries were constructed, comprising of 6 samples which each included 3 replications, resulting in a total of 43,825,614 to 55,003,390 raw reads and 41,657,756 to 52,397,778 high-quality clean reads per library (Table 1). The coverage of mapped reads ranged from 79.22 to 81.64%. The Q30 values, which are an important measure of RNA-seq quality, were all above 95.57% in all 18 libraries. Furthermore, the GC content of clean reads ranged from 48.85 to 51.1%. In summary, our RNA-seq data was high quality and reliable.

Table 1 RNA-seq data and quality statistics
Full size table

Furthermore, the principal components analysis (PCA) was performed based on the gene expression data obtained from fragments per kilobase of exon per million fragments mapped reads (FPKM). The results revealed that the 18 samples were divided into five subgroups. CK-1 and T-1 were classified into one subgroup, while CK-2, CK-3, T-2, and T-3 were classified into a distinct subgroup, respectively (Fig. 1B). These findings were consistent with the sample characteristics and indicated that A. H11648R and A. H11648 agave exhibit diverse transcriptional expression patterns during different stages of agave purple curl leaf disease.

To obtain the accurate differentially expressed genes (DEGs) among different samples (Suppl. Table 1), the cutoff of adjusted p-value < 0.05 was used in this research. A total of 547 (95 up-regulated and 451 down-regulated) DEGs were identified in CK-1-vs-T-1; 4,137 (1,782 up-regulated and 2,355 down-regulated) DEGs were identified in CK-2-vs-T-2; and 22,712 (14,786 up-regulated and 7,926 down-regulated) DEGs were identified in CK-3-vs-T-3 (Fig. 1C). The pairwise group CK-3-vs-T-3 had the greatest number of DEGs out of all three pairwise groups, which also was the time point possessing the most apparent phenotype of purple curl leaf disease in all three time points. In general, a total of 175, 1,182, and 19,657 DEGs were identified in CK-1-vs-T-1, CK-2-vs-T-2, and CK-3-vs-T-3, respectively. Moreover, 210 DEGs were simultaneously identified in all three pairwise groups, and 241, 341 and 2,924 DEGs were identified between CK-1-vs-T-1 and CK-2-vs-T-2, CK-1-vs-T-1 and CK-3-vs-T-3, CK-2-vs-T-2 and CK-3-vs-T-3, respectively (Fig. 1D). The above results suggested that the number of DEGs in CK-vs-T pairwise groups gradually increased along with the extended infection time, which indicated that the underlying molecular immune reactions between A. H11648R and A. H11648 plants may emerge a significant divergence over time.

Fig. 1
figure 1

Phenotypes of the A. H11648 and A. H11648R under the infection of purple curl leaf disease in agave and its RNA-seq results during different processes of purple curl leaf disease in agave. (A) The field phenotypes of A. H11648 and A. H11648R agave. The blade tips of A. H11648 and A. H11648R were sampled at 0 d (without D. neobrevipes, control), 60 d and 90 d after agave purple curl leaf disease infection. The sampled position is marked with red circles. (B) The principal components analysis of 18 samples based on the fragments per kilobase of exon per million mapped reads (FPKM). Group ‘CK-1’ indicates three biological replicates (‘F-0d-1’, ‘F-0d-2’, and ‘F-0d-3’) sampled on day 0 in A. H11648 agave. Group ‘CK-2’ indicates three biological replicates (‘F-60d-1’, ‘F-60d-2’, and ‘F-60d-3’) sampled on day 60 in A. H11648 agave. Group ‘CK-3’ indicates three biological replicates (‘F-90d-1’, ‘F-90d-2’, and ‘F-90d-3’) sampled on day 90 in A. H11648 agave. Likewise, group ‘T-1’ (‘R-0d-1’, ‘R-0d-2’, ‘R-0d-3’), group ‘T-2’ (‘R-60d-1’, ‘R-60d-2’, ‘R-60d-3’), and group ‘T-3’ (‘R-90d-1’, ‘R-90d-2’, and ‘R-90d-3’) represent the three biological replicates of A. H11648R plants sampled on day 0, day 60, and day 90, respectively. (C) The differentially expressed genes (DEGs) statistics among different samples. “CK” represents A. H11648 plants and “T” represents A. H11648R plants. “1”, “2” and “3” represent samples at 0 d, 60 d and 90 d after agave purple curl leaf disease infection, respectively. (D) The venn diagram of DEGs among different comparable groups

Full size image

Pathway enrichment analysis of DEGs among different pairwise groups

To gain a better understanding of the mechanism of the resistance to purple curl leaf disease in agave, we further analyzed the differentially expressed genes (DEGs) among different pairwise groups based on the KEGG and Gene Ontology (GO) databases, expecting to identify the pathways that are significantly relevant to agave disease resistance. In CK1 vs. T1 pairwise group, we found that the DEGs were mainly enriched into the pathways of limonene and pinene degradation, alpha-linolenic acid metabolism, starch and sucrose mechanism, and phenylpropanoid biosynthesis (Fig. 2A). In the CK2 vs. T2 pairwise group, the DEGs were mainly enriched into the pathways of plant hormone signal transduction, alpha-linolenic acid metabolism, starch and sucrose mechanism, flavonoid biosynthesis, and phenylpropanoid biosynthesis (Fig. 2B). In the CK3 vs. T3 pairwise group, the DEGs were mainly enriched into the pathways of alpha-linolenic acid metabolism, flavonoid biosynthesis, phenylpropanoid biosynthesis, and starch and sucrose metabolism (Fig. 2C). Furthermore, GO enrichment analysis for the common DEGs among these three comparable groups were performed. The results showed that the processes of systemic acquired resistance, pectin catabolic process, polysaccharide catabolic process, cell wall modification, carbohydrate catabolic process, galacturonan metabolic process, innate immune response, and defense response to bacterium were significantly enriched (Fig. 2D). Because pectin, polysaccharide, carbohydrate, and galacturonan are all important components of the plant cell wall, we speculated that plant cell wall modification and systemic acquired resistance may play an important role in regulating the resistance to agave purple curl leaf disease.

Fig. 2
figure 2

KEGG and GO enrichment analysis for the RNA-seq data during different stages of purple curl leaf disease in agave. The top 20 KEGG enriched pathways are presented for each pairwise group analysis: (A) CK-1-vs-T-1, (B) CK-2-vs-T-2, and (C) CK-3-vs-T-3. (D) The GO enriched pathways for the common DEGs among these three comparable groups. “CK” represents A. H11648 plants, and “T” represents A. H11648R plants. The numbers “1,” “2,” and “3” represent samples taken at 0 d, 60 d and 90 d after agave purple curl leaf disease infection, respectively. The “RichFactor” represents the ratio of the number of DEGs to the total number of agave genes in the pathway. Thus, a higher “RichFactor” indicates a higher degree of enrichment. Qvalue represents the Pvalue corrected by FDR (False Discovery Rate). Enriched GO pathway for the DEGs was defined as significantly enriched pathway when its Qvalue ≤ 0.05. The smaller the Qvalue, the more significant the pathway

Full size image

To strengthen our understanding for the above KEGG and GO enrichment analysis, a comparison between the agave and other species were conducted for the catalogues enriched by KEGG and GO enrichment analysis on phytoplasma resistance in plants. The results suggested that the catalogues enriched in both agave and other species include pentose and glucuronate interconversions, phenylpropanoid biosynthesis, sphingolipid metabolism, starch and sucrose metabolism, flavonoid biosynthesis, alpha-linolenic acid metabolism and other 163 catalogues (Suppl. Table 2). Considering that the number of DEGs at 60 d and 90 d after purple curl leaf disease infection was very high, we chosen 55 representative DEGs at both 60 d and 90 d after purple curl leaf disease infection, which all occupied the top sixty at the absolute value of FPKM (Fragments per Kilobase Million). The heatmap of these representative genes suggested that the down-regulated genes in A. H11648R plants were involved in the oxidation of phenols, the biosynthesis of fructan, anthocyanin and phenylpropanoid, the modification of cell walls, hypersensitive reaction, ROS metabolism, abiotic and biotic stresses, the degradation of damaged protein under stress conditions, the establishment of auxin gradients, the reinforcement of the plant cell wall and wounding or pathogen challenge, cell size determination and promoting both cambium activity and phloem specification (Fig. 3). The up-regulated genes in A. H11648R plants were involved in heavy-betal-binding protein, starch breakdown, the control of intracellular Na+ and K+ homeostasis, microtubule-associated protein, RNA silencing pathway, phytochrome B pathway, the biosynthesis of coenzyme A, GA and the maturation of mRNAs, the oxidative degradation of abscisic acid, microtubule-associated protein, development, pre-mRNA splicing, regulatory subunit of pyrophosphate-frutose, inorganic phosphate transportion, cytokinesis, nucleotide excision repair of damaged DNA and hydrolyzing glycosides and monolignol glucosides (Fig. 3). The above results provided a series of candidate genes and its relevant functions which may be involved in the regulation of agave resistance against purple curl leaf disease.

Fig. 3
figure 3

Heatmap for the crucial DEGs enriched by KEGG and GO analysis of all three stages in purple curl leaf disease. “F” represents A. H11648 plants, “R” represents A. H11648R plants. “60d” and “90d” represent samples taken at 60 d and 90 d after purple curl leaf disease infection in agave, respectively

Full size image

The important function of hypersensitive response (HR) in regulating the resistance to agave purple curl disease

To further clarify the underlying mechanism of the resistance to purple curl leaf disease in agave, we conducted a KEGG enrichment analysis on the common DEGs in the three pairwise groups (CK-1-vs-T-1, CK-2-vs-T-2 and CK-3-vs-T-3). The results of the enrichment analysis were similar to the previous enrichment results (Fig. 4A). Interestingly, the plant-pathogen interaction pathway was significantly enriched, in which pathogen/microbe-associated molecular patterns (PAMP/MAMP) immunity were triggered (Fig. 4B). Besides, Ca2+ signal and reactive oxygen species (ROS) burst are triggered, leading to hypersensitive response (HR) and cell wall reinforcement, and bacterial secretion system is also triggered by affecting RIN4 and then results in HR (Fig. 4B). Consistent with these findings, the heatmap analysis of HR-related genes showed that most of these genes exhibited significant expression between A. H11648 and A. H11648R agave, such as LOX (lipoxygenase), PR4 (wound-induced basic protein) and NPR1 (pathogenesis-related protein 1) genes (Fig. 4C). Taken together, we speculated that the resistance of A. H11648R against purple curl leaf disease may be regulated by HR in agave.

Fig. 4
figure 4

KEGG enrichment analysis of the common DEGs between the A. H11648R and A. H11648 plants during different stages of purple curl leaf disease. (A) The top twenty KEGG enrichment analysis of the common DEGs among the three pairwise groups (CK-1-vs-T-1, CK-2-vs-T-2, and CK-3-vs-T-3) during different stages of purple curl leaf disease. “CK” represents A. H11648 agave, and “T” represents A. H11648 agave. “1”, “2”, and “3” represent samples taken at 0 d, 60 d, and 90 d after purple curl leaf disease infection, respectively. The “RichFactor” indicates the ratio of DEGs to the total number of agave genes in the pathway. The higher the “RichFactor”, the greater the enrichment degree. The Pvalue represents the enriched degree of the pathway for the DEGs. The smaller the Pvalue, the more significant the pathway. (B) The plant-pathogen interaction pathway significantly enriched among the common DEGs between the A. H11648R and A. H11648 plants during different stages of purple curl leaf disease. (C) Heatmap for DEGs related to HR. “F” represents A. H11648 plants, “R” represents A. H11648R plants. “0d” “60d” and “90d” represent samples taken at 0 d, 60 d and 90 d after purple curl leaf disease infection in agave, respectively

Full size image

Weighted gene co-expression network analysis (WGCNA) of DEGs and qPCR analysis

Although we identified 55 representative DEGs in previous section, the vital disease-resistant genes are still unknown. To identify the crucial genes against purple curl disease in A. H11648, a weighted gene co-expression network analysis (WGCNA) was conducted. The results showed that the DEGs between A. H11648R and A. H11648 plants during different processes of purple curl leaf disease were divided into 14 modules based on their expression patterns (Suppl. Table 3). The bisque4 and lavenderblush2 modules were selected to build network analysis maps with Cytoscape software to identify hub genes against purple curl leaf disease, from which the DEGs were mostly up-regulated in A. H11648R compared to A. H11648 plants. Then, eight hub genes, five in the bisque4 module and three in the lavenderblush2 module, were found to be significantly related to the resistance to purple curl leaf disease in A. H11648. In the bisque4 module (Fig. 5A), the hub genes were as follow: 4’OMT2 ((S)-coclaurine N-methyl transferase isoformX1, 4’OMT2) is involved in the biosynthesis of benzylisoquinoline alkaloids and possesses insecticidal and bacteriostatic functions [21, 22]; ACLY (ATP-citrate synthase alpha chain protein, ACLY) is involved in the biosynthesis of flavonoids [23]; NCS1 (S-norcoclaurine synthase1-like, NCS1) is involved in the biosynthesis of the common precursor of all benzylisoquinoline alkaloids [24]; GTE10 (transcription factor GTE10-like) is a negative regulator of the environmental stress response of the ABA pathway in plants [25]; and SMO2 (methyl sterol monooxygenase 2-2-like) is involved in sterols biosynthesis and related the metabolism of alkaloids, flavonoids, and sterols [26]. In the lavenderblush2 module (Fig. 5B), the hub genes were as follow: FLS2 (LRR receptor-like serine/threonine-protein kinase FLS2) is a potent elicitor of the defense response to pathogen-associated molecular patterns (PAMPs) [27]; SQE1 (squalene epoxidase) is involved in steroid biosynthesis [28]; and RCOM (glycosyltransferase family 92 protein RCOM-like) may be responsible for the regulation of plant innate immunity and systemic acquired resistance, the phenylpropanoid metabolism, regulating the resistance to disease in plants [29,30,31,32].

To further investigate the function of these hub genes in regulating the resistance to purple curl leaf disease in A. H11648, qPCR was performed with the samples from leaves of A. H11648R and A. H11648 at 0 d, 30 d, 60 d, and 90 d after the infection of purple curl leaf disease. The results revealed that AsRCOM gene possess a very high gene expression level along with the whole disease processes in A. H11648R compared to A. H11648 plants (Fig. 5C). Considering that RCOM may play an important role in the regulation of plant innate immunity and systemic acquired resistance [29,30,31,32] and the process of HR may participate the resistance to purple curl leaf disease in agave (Fig. 4), we speculated that AsRCOM gene may be the most critical hub gene in regulating the resistance to purple curl leaf disease in agave.

Fig. 5
figure 5

Weighted gene co-expression network analysis of differentially expressed genes (DEGs) and qPCR analysis for the hub genes identified by WGCNA analysis. (A) The hub genes identified in the bisque4 module by WGCNA analysis. (B) The hub genes identified in the lavenderblush2 module by WGCNA analysis. (C) The relative expression levels of hub genes identified by WGCNA at different stages of agave purple curl leaf disease. “R” represents the A. H11648R plants, “S” represents the A. H11648 plants. Each experiment was performed with three biological replicates. Vertical bars represent ± the standard error of the mean (n = 3, n represents the biological replicates). *P < 0.05, **P < 0.01 and ***P < 0.001 (One-way ANOVA analysis of variance with Dunnett’s multiple comparison test), NS represents no significance

Full size image

Functional analysis of AsRCOM in agave

To confirm our hypothesis that AsRCOM is the key gene responsible for the resistance to purple curl leaf disease in agave, we created the AsRCOM overexpression (AsRCOM-OE) plants in A. H11648 mediated by Agrobacterium tumefaciens EHA105. The results of qPCR verification for the AsRCOM-OE plants found that the expression level of AsRCOM genes in AsRCOM-OE #1 and AsRCOM-OE #2 transgenic lines were significantly higher than that in the recipient A.H11648 (Supple. Figure 1), which were selected for the subsequent functional verification. Then, the resistance of AsRCOM-OE plants against agave purple curl disease was detected. The results revealed that the resistance of AsRCOM-OE plants against purple curl leaf disease was significantly higher than that in A. H11648 (Fig. 6A).

Furthermore, physiological index analyses were conducted between AsRCOM-OE plants and A. H11648 at various stages of purple curl leaf disease in agave (Fig. 6B). The results showed that the H2O2 contents in AsRCOM-OE plants increased by 33.4% and 55.8% at 60 d and 90 d after infected with purple curl leaf disease compared to that in A. H11648, respectively. Additionally, the superoxide anion contents in AsRCOM-OE plants increased by 18.1% and 25.6% at 60 d and 90 d after infected with purple curl leaf disease compared to that in A. H11648, respectively. The rates of ROS production in AsRCOM-OE plants also increased by 37% and 36.7% at 60 d and 90 d after infected with purple curl leaf disease compared to that in A. H11648, respectively. In contrast, the proline (PRO) contents in AsRCOM-OE plants decreased by 28.8% and 32.7% at 60 d and 90 d after infected with purple curl leaf disease compared to that in A. H11648, respectively. No significant differences were observed in these indexes at 0d (control) between AsRCOM-OE plants and A. H11648. These results suggested that the AsRCOM gene could significantly enhance the resistance to purple curl leaf disease in agave and the ROS contents in AsRCOM-OE plants were obviously higher than that in A. H11648 during the different stages of purple curl leaf disease in agave.

Fig. 6
figure 6

The phenotype and physiological indexes analysis of AsRCOM-OE plants exposed to purple curl leaf disease in agave. (A) The phenotype of AsRCOM-OE plants exposed to purple curl leaf disease in agave. Because purple curl leaf disease in agave was caused by D. neobrevipes, the same number (80) of D. neobrevipes were used to trigger the onset of purple curl leaf disease in the AsRCOM-OE plants and A. H11648. The recipient A. H11648 was set as control. (B) The physiological indexes analyses were conducted between AsRCOM-OE plants and A. H11648 during different stages of agave purple curl leaf disease. The leaves of AsRCOM-OE plants and A. H11648 were sampled to measure H2O2 contents, superoxide anion contents, rates of ROS production, and proline (PRO) contents at 0 d (control), 60 d, and 90 d after D. neobrevipes infection. Three biological replicates were performed for each experiment. The vertical bars represent ± the standard error of the mean (n = 3, where n represents the biological replicates). *P < 0.05, **P < 0.01, and ***P < 0.001 (One-way ANOVA analysis of variance with Dunnett’s multiple comparison test), and NS represents no significance

Full size image

Discussion

Recently, agave purple curl leaf disease has emerged as a major threat to the healthy development of agave industry [11]. Our research has identified a natural mutant from agave H.11,648, which exhibits remarkable resistance to the purple curl leaf disease. However, the underlying mechanism of this natural mutant in increased resistance is still unknown. Elucidating this mechanism and excavating the key resistant genes will be helpful for cultivating stable disease-resistant cultivar in agave. As an effective data mining method, RNA-seq analysis has been used to identify lots of powerfully stress-resistant genes [33, 34].

Therefore, to better understand the underlying mechanism of A. H11648R agave in enhancing the resistance to purple curl leaf disease, we conducted an RNA-seq analysis between A. H11648R and A. H11648 during different stages of purple curl leaf disease (Fig. 1). According to the KEGG and GO enrichment analysis results (Fig. 2), a series of common processes were shared with the RNA-seq results in other species, such as systemic acquired resistance, linoleic acid metabolism, pectin catabolic process and so on. Among these processes, we found that there were twenty representative candidate genes (Fig. 3) were also reported in regulating the resistance against phytoplasma in grapevine [35], mulberry [19], Paulownia fortunei [36], tomato [37], apple [38], coconut [39], and Mexican lime tree [40]. And there were thirty-five agave-specific genes that may play an important role in regulating the resistance against purple curl leaf disease in agave, including fructan 1-fructosyltransferase, protein indeterminate-domain 9, flavonoid 3’,5’-methyltransferase-like isoform X1, protein argonaute 16-like, abscisic acid 8’-hydroxylase 3, DNA repair protein UVH3 and so on (Fig. 3). These agave-specific genes are mainly involved in nucleotide excision repair of damaged DNA [41], the oxidative degradation of abscisic acid [42], the RNA silencing pathway [43], the reinforcement of the plant cell wall and wound or pathogen challenge [44], the establishment of auxin gradients [45] and the biosynthesis of fructan [46], which suggested that the processes of DNA repair, ABA and auxin signal pathway, the RNA silencing pathway and the reinforcement of plant cell wall may occupy a vital position in response to phytoplasma infection in agave. Besides, there were thirteen catalogues that were enriched only in agave compared to that in other species, which included non-homologous end-joining, caffeine metabolism, proteasome, glycosphingolipid biosynthesis-ganglio series, glycosphingolipid biosynthesis-globo series, thiamine metabolism, lipoic acid metabolism, glycosylphosphatidylinositol (GPI) -anchor biosynthesis, base excision repair, sulfur relay system, N-Glycan biosynthesis, detoxification, and transcription factor activity, protein binding (Suppl. Table 2). These catalogues exhibited the specificity of agave compared to other species during responding to phytoplasma infection.

Furthermore, as an important tool for excavating the central hub genes in multiple metabolism processes [47,48,49], WGCNA were conducted in order to decipher the pivotal genes relevant to the resistance to purple curl disease in agave. A total of eight hub genes were identified (Fig. 5), most of which are essential in increasing the resistance to insects and diseases in plants [18,19,20,21,22,23,24,25,26,27,28,29]. Except for the GTE10 and FLS2 genes, almost all of these genes are relevant to the metabolism of alkaloids, a kind of antibacterial compounds [50, 51], which suggest that alkaloids may be an effective external admixture for the defense against agave purple curl leaf disease. The results of qPCR analysis found that the AsRCOM gene may be the most crucial gene among all the hub genes identified by WGCNA analysis (Fig. 5C). We also found that the overexpression of AsRCOM gene in A. H11648 could significantly increase the resistance to agave purple curl leaf disease with large amounts of ROS accumulation compared to the recipient A. H11648 (Fig. 6). This finding is consistent with the previous researches that the glycosyltransferase could obviously enhance the plant immunity against biotic stresses [26,27,28,29] and ROS accumulation was an important defense responses in plants [52, 53]. Besides, reports showed that enhanced ROS production has been implicated in the hypersensitive response and plant defense signals in plants [54, 55], which further confirmed that HR may play an important role in regulating the resistance to purple curl leaf disease in agave. Similarly, the results of RNA-seq analysis also showed that HR may play a vital role in the interaction between plants and pathogenic bacteria in sugarcane, Brassica napus, Saccharum, rice, grapevine, Populus and sand pear [56,57,58,59,60,61,62]. Together, our research showed that the AsRCOM gene may regulate the resistance to purple curl leaf disease through HR pathway.

The present study has some limitations. The in-depth molecular mechanisms of AsRCOM-OE plants in enhancing the resistance to purple curl leaf disease in agave need to be excavated, including the signal regulation network, the effect of AsRCOM gene on the cellular and organic structures, the potential interactive protein or transcription factor of AsRCOM. Subsequently, the problem that whether the AsRCOM gene is related to HR should be verified. In addition, the mutants of AsRCOM gene mediated by CRISPR/Cas9 system in agave need to be created to supplement the function of AsRCOM gene in agave.

Conclusion

By RNA-seq analysis, we discovered that systemic acquired resistance accompanied with HR may have a crucial function in controlling resistance to purple curl leaf disease in agave. Additionally, we found the AsRCOM protein of glycosyltransferase family identified by WGCNA analysis could substantially enhance the resistance to purple curl leaf disease in agave. Together, our researches provide a profound understanding for the mechanisms of A. H11648R in enhancing the resistance to purple curl leaf disease and a series of important genetic resources for breeding disease-resistant cultivar in agave.

Methods

Plant materials and samples

Disease-susceptible A. H11648 and its natural mutant disease-resistant A. H11648R screened from disease-susceptible A. H11648 continuously infected with agave purple curl disease for many years were used as experimental materials. Agave purple curl leaf disease was induced by cultured D. neobrevipes. During the process of infection with D. neobrevipes, D. neobrevipes were artificial cultured with fresh agave leaves (Updated every three days) as food sources under the condition of 26 °C, 16 h/8 h light/dark photoperiod. Then, 80 live D. neobrevipes were separated from agave leaves by flicking the leaves slowly, and initially inoculated to A.H11648 and A.H11648R plants all with no D. neobrevipes infection, respectively. Subsequently, insect nets were used to exclude A.H11648 and A.H11648R plants from wild D. neobrevipes or other insects until the end of the experiment. Besides, A.H11648 and A.H11648R plants with no inoculation of D. neobrevipes were set as control to confirm that the disease performances were just caused by purple curl leaf disease. Agave blade tips were sampled at 0 d, 60 d, and 90 d post-infection (Fig. 1A), in which the sample at 0 d (On the day of D. neobrevipes inoculation, but without inoculation) was chosen as control. Blade tips were sampled for RNA-seq analysis with three biological replicates at each timepoint. All samples were collected and immediately stored at -80 °C for RNA extraction and RNA-seq analysis.

RNA extraction and Illumina sequencing

Total RNA was extracted from the sampled blade tips using the TRIzol reagent (Aidlab, Beijing, China) according to the manufacturer’s protocol. The concentration of total RNA in each sample was measured using a NanoDrop microvolume spectrophotometer (Thermo Scientific NanoDrop Products, Waltham, MA, USA). RNA samples were treated to remove genomic DNA contamination and cDNA was synthesized using a one-step reverse transcription kit (Transgen, Beijing, China). cDNA library construction and sequencing were completed by Guangzhou Genedenovo, Guangzhou, China. The Illumina HiSeq2500 platform (Illumina, San Diego, CA, America) was used for RNA-seq with 6 GB clean, paired-end data.

RNA-seq analysis

The quality of the raw RNA-seq data was evaluated using FastQC v0.11.2 [63]. Reads with low quality (N content > 5%, low-quality bases > 30% (quality < 29)) were removed, and higher quality reads were trimmed to eliminate adapters with Trimmomatic (v0.36.5) [64]. The high-quality reads were then assembled using Trinity v2.4.0 (https://github.com/trinityrnaseq/trinityrnaseq/wiki) to generate a unigene library. BlastX (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastx&PAGETYPE=BlastSearch&LINKLOC=blasthome) was used with multiple online databases, including Eggnog (http://eggnogdb.embl.de/#/app/emapper), Uniprot (https://www.uniprot.org/), GO (http://wegolgenomics.org.cn/cgi-bin/wego/index.pl), Pfam (http://pfam.xfam.org/), Clusters of Orthologous Groups (COGs, https://www.ncbi.nlm.nih.gov/research/cog), TMhmm, (https://services.healthtech.dtu.dk/service.php?TMHMM-2.0), and KEGG (http://www.genome.jp/kegg/kegg2.html) [65,66,67], to annotate the assembled unigene transcripts with an E-value cut-off of 10− 5. DEGs were identified using the following criteria: FDR ≤ 0.01 and |log2 ratio| ≥ 1 (FDR indicating the false discovery rate and log2 ratio representing log2 of fold change for the RPKM of up/down-regulation). Heatmap analysis of the DEGs in the RNA-seq data was performed using TBtools [68], based on the RPKM (reads per kilobase per million mapped reads). GO enrichment analysis with the identified DEGs was conducted in R-Studio (v3.3.0) with the ggplot2 package (v3.3.5, https://ggplot2.tidyverse.org/). The GO term categories, such as molecular function, were confirmed with a p-value cut-off of < 0.05, hypergeometric distribution, and Bonferroni correction.

Weighted gene co-expression network analysis (WGCNA) of DEGs

The WGCNA R package (v1.70-3) [69] was utilized to perform a weighted gene co-expression network analysis (WGCNA) using normalized expression values of the DEGs. A range of powers from 1 to 20 were tested for network threshold, and the resulting network was visualized using Cytoscape (v3.9.1, https://cytoscape.org/).

qPCR verification for the hub genes and its function analysis in agave

In the RNA extraction and real-time quantitative RT-PCR (qRT-PCR) experiments, the TransGen RNA extraction kit (Transgen, Beijing, China) was utilized to extract total RNA from agave leaves following the manufacturer’s instructions. The first cDNA strand was synthesized using the TransGen RT SuperMix (Transgen, Beijing, China) kit. The primers for qRT-PCR experiments and subsequent transgenic detection primers were stored in Suppl. Tables 4, with actin gene serving as an internal reference. The qRT-PCR was conducted using the TransGen SYBR® Green PCR kit (Transgen, Beijing, China) in accordance with the manufacturer’s instructions. Three biological replicates were carried out for each experiment.

The A.H11648 was selected as recipient to transform with the pCAMBIA1302-AsRCOM plasmid and pCAMBIA1302 empty vector using Agrobacterium tumefaciens strain EHA105. Shoot apex explants were chosen for callus culture in agave. Then the callus were transformed with the relative vectors. After selective cultivation with hygromycin, PCR detection and transplanting to the field, positive transgenic agave seedlings and its recipient A.H11648 that all were three months old were planted in the same area covered with insect nets to exclude them from wild D. neobrevipes or other insects until the end of the experiment. Then, 80 live D. neobrevipes were separated from agave leaves by flicking the leaves slowly, and inoculated to AsRCOM overexpression and A.H11648 plants all with no D. neobrevipes infection, respectively. The phenotype of these plants were observed and recorded to evaluate their resistance against purple curl leaf disease every seven days in agave.

Physiological and biochemical indexes analysis

The physiological and biochemical indexes analysis was performed as described previously [70]. The H2O2 contents, ROS production rates, superoxide anion contents, and PRO contents were measured in accordance with the kit instructions (COMIN, Suzhou, China). The leaves of the AsRCOM-OE and the recipient agave during different stages of purple curl leaf disease were sampled and performed for the above physiological and biochemical indexes analysis. Three biological replicates were conducted for each experiments.

Statistical analysis

All data were analyzed using GraphPad Prism 7 software. The statistical comparisons were performed using One-way ANOVA analysis of variance with Dunnett’s multiple comparison test. Differences were considered significant when the P-value was less than 0.05. All experiments were conducted with three biological replicates.

Data Availability

All data and materials are conserved in our laboratory. RNA-seq raw reads have been deposited in the SRA metadata available at https://dataview.ncbi.nlm.nih.gov/object/PRJNA852945?reviewer=fuoka4k1e6knf5434mhljh1uit (BioProject: PRJNA852945), and the RNA-seq mRNA sequences have been deposited in GenBank BankIt submission (Submission ID: 2,723,837; 2,723,846; 2,723,857; 2,723,860).

References

  1. Debnath M, Pandey M, Sharma R, Thakur GS, Lal P. Biotechnological intervention of Agave sisalana: a unique fiber yielding plant with medicinal property. J Med Plant Res. 2010;4:177–87.

    CAS  Google Scholar 

  2. Dalton R. Saving the agave. Nature. 2005;438:1070–71.

    Article  CAS  PubMed  Google Scholar 

  3. Quiñones-Muñoz TA, Villanueva-Rodríguez SJ, Torruco-Uco JG. Nutraceutical properties of Medicago sativa L., Agave spp., Zea mays L. and Avena sativa L.: a review of metabolites and mechanisms. Metabolites. 2022; 12: 806.

  4. Evdokimova OL, Alves CS, Whiffen RMK, Ortega Z, Tomás H, Rodrigues J. Cytocompatible cellulose nanofibers from invasive plant species Agave Americana L. and Ricinus communis L.: a renewable green source of highly crystalline nanocellulose. J Zhejiang Univ Sci B. 2021;22:450–61.

    Article  CAS  Google Scholar 

  5. Niechayev NA, Jones AM, Rosenthal DM, Davis SC. A model of environmental limitations on production of Agave Americana L. grown as a biofuel crop in semi-arid regions. J Exp Bot. 2019;70:6549–59.

    Article  CAS  PubMed  Google Scholar 

  6. Vega-Ramos KL, Uvalle-Bueno JX, Gomez-Leyva JF. Molecular variability among isolates of Fusarium oxysporum associated with root rot disease of Agave tequilana. Biochem Genet. 2013;51:243–55.

    Article  CAS  PubMed  Google Scholar 

  7. Luján R, Lledías F, Martínez LM, Barreto R, Cassab GI, Nieto-Sotelo J. Small heat-shock proteins and leaf cooling capacity account for the unusual heat tolerance of the central spike leaves in Agave tequilana var. Weber. Plant Cell Environ. 2009;32:1791–803.

    Article  PubMed  Google Scholar 

  8. Torres-García I, León-Jacinto A, Vega E, Moreno-Calles AI, Casas A. Integral projection models and sustainable forest management of Agave inaequidens in western Mexico. Front Plant Sci. 2020;11:1224.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Escobar-Flores JG, Sandoval S, Gámiz-Romero E. Unmanned aerial vehicle images in the machine learning for agave detection. Environ Sci Pollut Res Int. 2022;29:61662–73.

    Article  PubMed  Google Scholar 

  10. Peng-Peng L, Wei-Huai W, Jin-Long Z, Gui-Hua W, Chun-Ping H, Pei-Qun L, Xing H, Yan-Qiong L, Ke-Xian Y. Establishment and optimization of single-tube nested PCR detection technique for phytoplasma related to sisal purple leafroll disease. J Agri Biotec. 2021;29:1426–34.

    Google Scholar 

  11. Biao H, Rong Y, Li-hong X, Jia-liu Z. Etiology and pathogen identification of purple leaf roll virus and application of resistant seedlings. J Anhui Agri Sci. 2015;43:177–9.

    Google Scholar 

  12. Long ZY, Feng LJ, Sheng YQ, Biao HY, Feng L, Bo WJ, Ping LG, Jie Y. Research progress and prospects of three main diseases of sisal. Chin J TRO AGRI. 2020;40:72–82.

    Google Scholar 

  13. Bressan A, Purcell AH. Effect of benzothiadiazole on transmission of X-Disease phytoplasma by the vector colladonus montanus to Arabidopsis thaliana, a new experimental host plant. Plant Dis. 2005;89:1121–4.

    Article  CAS  PubMed  Google Scholar 

  14. Chen P, Zhang Y, Li Y, Yang Q, Li Q, Chen L, Chen Y, Ye X, Tan B, Zheng X, Cheng J, Wang W, Li J, Feng J. Jujube witches’ broom phytoplasma effector Zaofeng3, a homologous effector of SAP54, induces abnormal floral organ development and shoot proliferation. Phytopathology. 2023. https://doi.org/10.1094/PHYTO-10-21-0448-R

    Article  PubMed  Google Scholar 

  15. Chen P, Chen L, Ye X, Tan B, Zheng X, Cheng J, Wang W, Yang Q, Zhang Y, Li J, Feng J. Phytoplasma effector Zaofeng6 induces shoot proliferation by decreasing the expression of ZjTCP7 in Ziziphus jujuba. Hortic Res. 2022;9:uhab032.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Yang Q, Li Q, Gu L, Chen P, Zhang Y, Li Y, Chen Y, Ye X, Tan B, Zheng X, et al. The jujube TCP transcription factor ZjTCP16 regulates plant growth and cell size by affecting the expression of genes involved in plant morphogenesis. Forests. 2022;13:723.

    Article  Google Scholar 

  17. Wang S, Wang S, Li M, Su Y, Sun Z, Ma H. Combined transcriptome and metabolome analysis of Nerium indicum L. elaborates the key pathways that are activated in response to witches’ broom disease. BMC Plant Biol. 2022;22:291.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Xu Y, Wang C, Kong D, Cao M, Zhang Q, Tahir M, Yang Y, Yang S, Bo W, Pang X. Identification of high tolerance to jujube witches’ broom in Indian jujube (Ziziphus mauritiana Lam.) And mining differentially expressed genes related to the tolerance through transcriptome analysis. Plants. 2023;12:2082.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Liu C, Dong X, Xu Y, Dong Q, Wang Y, Gai Y, Ji X. Transcriptome and DNA methylome reveal insights into phytoplasma infection responses in mulberry (Morus multicaulis Perr). Front Plant Sci. 2021;12:697702.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Gai YP, Yuan SS, Liu ZY, Zhao HN, Liu Q, Qin RL, Fang LJ, Ji XL. Integrated phloem sap mRNA and protein expression analysis reveals phytoplasma-infection responses in mulberry. Mol Cell Proteomics. 2018;17:1702–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Chen JZ, Lu X, Hu YQ, Guo HH, Ma XL, Guo X, Jiang ZB, Wang F. Research progress on chemical constituents and pharmacological studies on root bark of Lycium barbarum. Zhongguo Zhong Yao Za Zhi. 2021;46:3066–75.

    PubMed  Google Scholar 

  22. Quiroz-Carreño S, Pastene-Navarrete E, Espinoza-Pinochet C, Muñoz-Núñez E, Devotto-Moreno L, Céspedes-Acuña CL, Alarcón-Enos J. Assessment of insecticidal activity of benzylisoquinoline alkaloids from chilean Rhamnaceae plants against fruit-fly Drosophila melanogaster and the Lepidopteran Crop Pest Cydia pomonella. Molecules. 2020;25:5094.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Fatland BL, Ke J, Anderson MD, Mentzen WI, Cui LW, Allred CC, Johnston JL, Nikolau BJ, Wurtele ES. Molecular characterization of a heteromeric ATP-citrate lyase that generates cytosolic acetyl-coenzyme A in Arabidopsis. Plant Physiol. 2002;130:740–56.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Witz S, Panwar P, Schober M, Deppe J, Pasha FA, Lemieux MJ, Möhlmann T. Structure-function relationship of a plant NCS1 member–homology modeling and mutagenesis identified residues critical for substrate specificity of PLUTO, a nucleobase transporter from Arabidopsis. PLoS ONE 9: e91343.

  25. Kim MJ, Shin R, Schachtman DP. A nuclear factor regulates abscisic acid responses in Arabidopsis. Plant Physiol. 2009;151:1433–45.

    Article  CAS  PubMed  Google Scholar 

  26. Darnet S, Bard M, Rahier A. Functional identification of sterol-4alpha-methyl oxidase cDNAs from Arabidopsis thaliana by complementation of a yeast erg25 mutant lacking sterol-4alpha-methyl oxidation. FEBS Lett. 2001;508:39–43.

    Article  CAS  PubMed  Google Scholar 

  27. Yuan M, Jiang Z, Bi G, Nomura K, Liu M, Wang Y, Cai B, Zhou JM, He SY, Xin XF. Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature. 2021;592:105–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Posé D, Castanedo I, Borsani O, Nieto B, Rosado A, Taconnat L, Ferrer A, Dolan L, Valpuesta V, Botella MA. Identification of the Arabidopsis dry2/sqe1-5 mutant reveals a central role for sterols in drought tolerance and regulation of reactive oxygen species. Plant J. 2009;59:63–76.

    Article  PubMed  Google Scholar 

  29. Ebert B, Birdseye D, Liwanag AJM, Laursen T, Rennie EA, Guo X, Catena M, Rautengarten C, Stonebloom SH, Gluza P, Pidatala VR, Andersen MCF, Cheetamun R, Mortimer JC, Heazlewood JL, Bacic A, Clausen MH, Willats WGT, Scheller HV. The three members of the Arabidopsis glycosyltransferase family 92 are functional β-1,4-galactan synthases. Plant Cell Physiol. 2018;59:2624–36.

    Article  CAS  PubMed  Google Scholar 

  30. Mohnike L, Rekhter D, Huang W, Feussner K, Tian H, Herrfurth C, Zhang Y, Feussner I. The glycosyltransferase UGT76B1 modulates N-hydroxy-pipecolic acid homeostasis and plant immunity. Plant Cell. 2021;33:735–49.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Huang XX, Wang Y, Lin JS, Chen L, Li YJ, Liu Q, Wang GF, Xu F, Liu L, Hou BK. The novel pathogen-responsive glycosyltransferase UGT73C7 mediates the redirection of phenylpropanoid metabolism and promotes SNC1-dependent Arabidopsis immunity. Plant J. 2021;107:149–65.

    Article  CAS  PubMed  Google Scholar 

  32. Chang B, Zhao L, Feng Z, Wei F, Zhang Y, Zhang Y, Huo P, Cheng Y, Zhou J, Feng H. Galactosyltransferase GhRFS6 interacting with GhOPR9 involved in defense against Verticillium wilt in cotton. Plant Sci. 2023;328:111582.

    Article  CAS  PubMed  Google Scholar 

  33. Ludwig-Müller J. What can we learn from -omics approaches to understand clubroot disease? Int J Mol Sci. 2022;23:6293.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Campos MD, Félix MDR, Patanita M, Materatski P, Varanda C. High throughput sequencing unravels tomato-pathogen interactions towards a sustainable plant breeding. Hortic Res. 2021;8:171.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Hren M, Nikolić P, Rotter A, Blejec A, Terrier N, Ravnikar M, Dermastia M, Gruden K. Bois noir’ phytoplasma induces significant reprogramming of the leaf transcriptome in the field grown grapevine. BMC Genomics. 2009;10:460.

    Article  PubMed  Google Scholar 

  36. Fan G, Dong Y, Deng M, Zhao Z, Niu S, Xu E. Plant-pathogen interaction, circadian rhythm, and hormone-related gene expression provide indicators of phytoplasma infection in Paulownia fortunei. Int J Mol Sci. 2014;15:23141–62.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Buoso S, Pagliari L, Musetti R, Martini M, Marroni F, Schmidt W, Santi S. Candidatus Phytoplasma solani’ interferes with the distribution and uptake of iron in tomato. BMC Genomics. 2019;20:703.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Weil T, Ometto L, Esteve-Codina A, Gómez-Garrido J, Oppedisano T, Lotti C, Dabad M, Alioto T, Vrhovsek U, Hogenhout S, Anfora G. Linking omics and ecology to dissect interactions between the apple proliferation phytoplasma and its psyllid vector Cacopsylla melanoneura. Insect Biochem Mol Biol. 2020;127:103474.

    Article  CAS  PubMed  Google Scholar 

  39. Nejat N, Cahill DM, Vadamalai G, Ziemann M, Rookes J, Naderali N. Transcriptomics-based analysis using RNA-Seq of the coconut (Cocos nucifera) leaf in response to yellow decline phytoplasma infection. Mol Genet Genomics. 2015;290:1899–910.

    Article  CAS  PubMed  Google Scholar 

  40. Monavarfeshani A, Mirzaei M, Sarhadi E, Amirkhani A, Khayam Nekouei M, Haynes PA, Mardi M, Salekdeh GH. Shotgun proteomic analysis of the Mexican lime tree infected with Candidatus Phytoplasma aurantifolia. J Proteome Res. 2013;12:785–95.

    Article  CAS  PubMed  Google Scholar 

  41. Liu Z, Hall JD, Mount DW. Arabidopsis UVH3 gene is a homolog of the Saccharomyces cerevisiae RAD2 and human XPG DNA repair genes. Plant J. 2001;26:329–38.

    Article  CAS  PubMed  Google Scholar 

  42. Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, Hirai N, Koshiba T, Kamiya Y, Nambara E. The Arabidopsis cytochrome P450 CYP707A encodes ABA 8’-hydroxylases: key enzymes in ABA catabolism. EMBO J. 2004;23:1647–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Rowley MJ, Avrutsky MI, Sifuentes CJ, Pereira L, Wierzbicki AT. Independent chromatin binding of ARGONAUTE4 and SPT5L/KTF1 mediates transcriptional gene silencing. PLoS Genet. 2011;7:e1002120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Martin JL, McMillan FM. SAM (dependent) I AM: the S-adenosylmethionine-dependent methyltransferase fold. Curr Opin Struct Biol. 2002;12:783–93.

    Article  CAS  Google Scholar 

  45. Cui D, Zhao J, Jing Y, Fan M, Liu J, Wang Z, Xin W, Hu Y. The Arabidopsis IDD14, IDD15, and IDD16 cooperatively regulate lateral organ morphogenesis and gravitropism by promoting auxin biosynthesis and transport. PLoS Genet. 2013;9:e1003759.

    Article  CAS  PubMed  Google Scholar 

  46. Lüscher M, Hochstrasser U, Vogel G, Aeschbacher R, Galati V, Nelson CJ, Boller T, Wiemken A. Cloning and functional analysis of sucrose: sucrose 1-fructosyltransferase from tall fescue. Plant Physiol. 2000;124:1217–28.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Jin Y, Li J, Zhu Q, Du X, Liu F, Li Y, Ahmar S, Zhang X, Sun J, Xue F. GhAPC8 regulates leaf blade angle by modulating multiple hormones in cotton (Gossypium hirsutum L). Int J Biol Macromol. 2022;195:217–28.

    Article  CAS  PubMed  Google Scholar 

  48. Lian J, Liu W, Sun Y, Men S, Wu J, Zeb A, Yang T, Ma LQ, Zhou Q. Nanotoxicological effects and transcriptome mechanisms of wheat (Triticum aestivum L.) under stress of polystyrene nanoplastics. J Hazard Mater. 2022;423:127241.

    Article  CAS  PubMed  Google Scholar 

  49. Cai H, Yu N, Liu Y, Wei X, Guo C. Meta-analysis of fungal plant pathogen Fusarium oxysporum infection-related gene profiles using transcriptome datasets. Front Microbiol. 2022;13:970477.

    Article  PubMed Central  Google Scholar 

  50. Chakraborty M, Mahmud NU, Ullah C, Rahman M, Islam T. Biological and biorational management of blast diseases in cereals caused by Magnaporthe oryzae. Crit Rev Biotechnol. 2021;41:994–1022.

    Article  CAS  PubMed  Google Scholar 

  51. Liu H, Yang S, Li T, Ma S, Wang P, Wang G, Su S, Ding Y, Yang L, Zhou X, Yang S. Design, synthesis and bioactivity evaluation of novel 2-(pyrazol-4-yl)-1,3,4-oxadiazoles containing an imidazole fragment as antibacterial agents. Molecules. 2023;28:2442.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Gao M, He Y, Yin X, Zhong X, Yan B, Wu Y, Chen J, Li X, Zhai K, Huang Y, Gong X, Chang H, Xie S, Liu J, Yue J, Xu J, Zhang G, Deng Y, Wang E, Tharreau D, Wang GL, Yang W, He Z. Ca2+ sensor-mediated ROS scavenging suppresses rice immunity and is exploited by a fungal effector. Cell. 2021;184:5391–404.

    Article  CAS  PubMed  Google Scholar 

  53. Meng X, Zhang S. MAPK cascades in plant disease resistance signaling. Annu Rev Phytopathol. 2013;51:245–66.

    Article  CAS  PubMed  Google Scholar 

  54. Liang X, Zhang L, Natarajan SK, Becker DF. Proline mechanisms of stress survival. Antioxid Redox Signal. 2013;19:998–1011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Dutta A, Chan SH, Pauli NT, Raina R. Hypersensitive response-like lesions 1 codes for AtPPT1 and regulates accumulation of ROS and defense against bacterial pathogen Pseudomonas syringae in Arabidopsis thaliana. Antioxid Redox Signal. 2015;22:785–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Wu Q, Chen Y, Zou W, Pan YB, Lin P, Xu L, Grisham MP, Ding Q, Su Y, Que Y. Genome-wide characterization of sugarcane catalase gene family identifies a ScCAT1 gene associated disease resistance. Int J Biol Macromol. 2023;232:123398.

    Article  CAS  PubMed  Google Scholar 

  57. Zhou T, Xu W, Hirani AH, Liu Z, Tuan PA, Ayele BT, Daayf F, McVetty PBE, Duncan RW, Li G. Transcriptional insight into Brassica napus resistance genes LepR3 and Rlm2-mediated defense response against the Leptosphaeria maculans infection. Front Plant Sci. 2019;10:823.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Zhao Z, Zhang R, Wang D, Zhang J, Zang S, Zou W, Feng A, You C, Su Y, Wu Q, Que Y. Dissecting the features of TGA gene family in Saccharum and the functions of under biotic stresses. Plant Physiol Biochem. 2023;200:107760.

    Article  CAS  PubMed  Google Scholar 

  59. Matić S, Bagnaresi P, Biselli C, Orru’ L, Amaral Carneiro G, Siciliano I, Valé G, Gullino ML, Spadaro D. Comparative transcriptome profiling of resistant and susceptible rice genotypes in response to the seedborne pathogen Fusarium fujikuroi. BMC Genomics. 2016;17:608.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Weng K, Li ZQ, Liu RQ, Wang L, Wang YJ, Xu Y. Transcriptome of Erysiphe necator-infected Vitis pseudoreticulata leaves provides insight into grapevine resistance to powdery mildew. Hortic Res. 2014;1:14049.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Liang H, Staton M, Xu Y, Xu T, Leboldus J. Comparative expression analysis of resistant and susceptible Populus clones inoculated with Septoria musiva. Plant Sci. 2014;223:69–78.

    Article  CAS  PubMed  Google Scholar 

  62. Wang H, Lin J, Chang Y, Jiang CZ. Comparative transcriptomic analysis reveals that ethylene/H2O2-mediated hypersensitive response and programmed cell death determine the compatible interaction of sand pear and Alternaria alternata. Front Plant Sci. 2017;8:195.

    PubMed  PubMed Central  Google Scholar 

  63. Brown J, Pirrung M, McCue LA. FQC Dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics. 2017;33:3137–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Kanehisa M, Goto S, KEGG. Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed Central  Google Scholar 

  67. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–D92.

    Article  CAS  PubMed  Google Scholar 

  68. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13:1194–202.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  70. Liu S, Liu S, Wang M, Wei T, Meng C, Wang M, Xia G. A wheat SIMILAR TO RCD-ONE gene enhances seedling growth and abiotic stress resistance by modulating redox homeostasis and maintaining genomic integrity. Plant Cell. 2014;26:164–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Thanks to all the authors for their contributions to plant material treatment, data analysis and constructive suggestions in manuscript preparation.

Funding

This research was financially supported by the earmarked fund for China Agriculture Research System (NO.CARS-16), Guangdong Provincial Team of Technical System Innovation for Sugarcane Sisal Industry (2022KJ104-03), and supported by Hainan Provincial Natural Science Foundation of China (323QN295). These funding bodies were not involved in the experimental design of the study, data collection, analysis and interpretation as well as in writing the manuscript.

Author information

Author notes

  1. Zhiwei Lu and Xiaowan Hou contributed equally to the present work.

Authors and Affiliations

  1. Zhanjiang City Key Laboratory for Tropical Crops Genetic Improvement, South Subtropical Crops Institute, Chinese Academy of Tropical Agricultural Sciences, Zhanjiang, Guangdong, 524091, China

    Zhiwei Lu, Zhi Ke, Yanmei Zhang, ZiPing Yang & Wenzhao Zhou

  2. State Key Laboratory of Plant Physiology and Biochemistry, College of Biological Sciences, China Agricultural University, Beijing, 100193, China

    Zhiwei Lu

  3. Institute of Crop Science, Key Laboratory of Biology and Genetic Improvement of Triticeae Crops, Chinese Academy of Agricultural Sciences (CAAS), National Key Facility for Crop Gene Resources and Genetic Improvement, Ministry of Agriculture, Beijing, 100081, China

    Zhiwei Lu

  4. Key Laboratory for Postharvest Physiology and Technology of Tropical Horticultural Products of Hainan Province, South Subtropical Crops Research Institute, Chinese Academy of Tropical Agricultural Sciences, Zhanjiang, Guangdong, 524091, China

    Xiaowan Hou

Authors
  1. Zhiwei Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaowan Hou
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zhi Ke
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yanmei Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. ZiPing Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Wenzhao Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Z.W.L and X.W.H conceived and designed experiments. Z.K, Y.M.Z and Z.P.Y conducted RNA-seq analyses. Z.W.L prepared the figures and tables. Z.W.L and X.W.H drafted the manuscript. All the authors critically reviewed the manuscript and provided intellectual inputs. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Wenzhao Zhou.

Ethics declarations

Ethics approval and consent to participate

All methods were performed in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Supplementary Table S1

. The RPKM of DEGs between CK-1 and T-1 group

Supplementary Material 2: Supplementary Table S2

. The catalogues enriched by KEGG or GO enrichment analysis on phytoplasma resistance in agave and other species and the relevant genes in agave

Supplementary Material 3: Supplementary Table S3

. The modules of RNA-seq data divided by WGCNA

Supplementary Material 4: Supplementary Table S4

. The primer sequences of relevant genes in qRT-PCR and PCR experiments

Supplementary Material 5: Supplementary Fig. 1

. The expression level analysis of AsRCOM gene in AsRCOM-OE transgenic lines by qPCR. A total of three transgenic lines of AsRCOM-OE were obtained. The expression level of AsRCOM gene in the recipient A.H11648 was set as 1. Three biological replicates were performed for each experiment. The vertical bars represent ± the standard error of the mean (n = 3, where n represents the biological replicates). *P < 0.05, **P < 0.01, and ***P < 0.001 (One-way ANOVA analysis of variance with Dunnett’s multiple comparison test), and NS represents no significance

Supplementary Material 6: Supplement file 1

. The amino acid sequences of RCOM proteins in multiple species. AsRCOM represents the amino acids in Agave sisalana. hybrid 11648

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

Lu, Z., Hou, X., Ke, Z. et al. A newly identified glycosyltransferase AsRCOM provides resistance to purple curl leaf disease in agave. BMC Genomics 24, 669 (2023). https://doi.org/10.1186/s12864-023-09700-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09700-y

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us