
The scab disease, caused by Elsinoe perseae, poses a significant risk to avocado (Persea americana) production in countries with warm and humid climates. Although the genome has been published, the precise virulence factors accountable for the pathogenicity of E. perseae have not yet been determined. The current study employed an in silico approach to identify and functionally characterize the secretory proteins of E. perseae. A total of 654 potential secretory proteins were identified, of which 190 were classified as carbohydrate-active enzymes (CAZymes), 49 as proteases, and 155 as potential effectors. A comparison to six other closely related species identified 40 species-specific putative effectors in E. perseae, indicating their specific involvement in the pathogenicity of E. perseae on avocado. The data presented in this study might be valuable for further research focused on understanding the molecular mechanisms that contribute to the pathogenicity of E. perseae on avocado.
Citation: Biju Vadakkemukadiyil Chellappan. Comparative secretome analysis unveils species-specific virulence factors in Elsinoe perseae, the causative agent of the scab disease of avocado (Persea americana)[J]. AIMS Microbiology, 2024, 10(4): 894-916. doi: 10.3934/microbiol.2024039
[1] | Tongmeng Jiang, Pan Jin, Guoxiu Huang, Shi-Cheng Li . The function of guanylate binding protein 3 (GBP3) in human cancers by pan-cancer bioinformatics. Mathematical Biosciences and Engineering, 2023, 20(5): 9511-9529. doi: 10.3934/mbe.2023418 |
[2] | Prathibha Ambegoda, Hsiu-Chuan Wei, Sophia R-J Jang . The role of immune cells in resistance to oncolytic viral therapy. Mathematical Biosciences and Engineering, 2024, 21(5): 5900-5946. doi: 10.3934/mbe.2024261 |
[3] | Yujia Xi, Liying Song, Shuang Wang, Haonan Zhou, Jieying Ren, Ran Zhang, Feifan Fu, Qian Yang, Guosheng Duan, Jingqi Wang . Identification of basement membrane-related prognostic signature for predicting prognosis, immune response and potential drug prediction in papillary renal cell carcinoma. Mathematical Biosciences and Engineering, 2023, 20(6): 10694-10724. doi: 10.3934/mbe.2023474 |
[4] | Sufang Wu, Hua He, Jingjing Huang, Shiyao Jiang, Xiyun Deng, Jun Huang, Yuanbing Chen, Yiqun Jiang . FMR1 is identified as an immune-related novel prognostic biomarker for renal clear cell carcinoma: A bioinformatics analysis of TAZ/YAP. Mathematical Biosciences and Engineering, 2022, 19(9): 9295-9320. doi: 10.3934/mbe.2022432 |
[5] | Samantha Erwin, Stanca M. Ciupe . Germinal center dynamics during acute and chronic infection. Mathematical Biosciences and Engineering, 2017, 14(3): 655-671. doi: 10.3934/mbe.2017037 |
[6] | Jun Wang, Mingzhi Gong, Zhenggang Xiong, Yangyang Zhao, Deguo Xing . Immune-related prognostic genes signatures in the tumor microenvironment of sarcoma. Mathematical Biosciences and Engineering, 2021, 18(3): 2243-2257. doi: 10.3934/mbe.2021113 |
[7] | Jie Chen, Jinggui Chen, Bo Sun, Jianghong Wu, Chunyan Du . Integrative analysis of immune microenvironment-related CeRNA regulatory axis in gastric cancer. Mathematical Biosciences and Engineering, 2020, 17(4): 3953-3971. doi: 10.3934/mbe.2020219 |
[8] | Caiyun Wu, Cong Ma, Jing Yuan, Pei Zhou . Exploration of potential therapeutic and prognostic value of CXC chemokines in cervical squamous cell carcinoma and endocervical adenocarcinoma based on bioinformatics analysis. Mathematical Biosciences and Engineering, 2021, 18(6): 8201-8222. doi: 10.3934/mbe.2021407 |
[9] | Xin Lin, Xingyuan Li, Binqiang Ma, Lihua Hang . Identification of novel immunomodulators in lung squamous cell carcinoma based on transcriptomic data. Mathematical Biosciences and Engineering, 2022, 19(2): 1843-1860. doi: 10.3934/mbe.2022086 |
[10] | Yong Luo, Qingbiao Chen, Jingbo Lin . Identification and validation of a tumor mutation burden-related signature combined with immune microenvironment infiltration in adrenocortical carcinoma. Mathematical Biosciences and Engineering, 2022, 19(7): 7055-7075. doi: 10.3934/mbe.2022333 |
The scab disease, caused by Elsinoe perseae, poses a significant risk to avocado (Persea americana) production in countries with warm and humid climates. Although the genome has been published, the precise virulence factors accountable for the pathogenicity of E. perseae have not yet been determined. The current study employed an in silico approach to identify and functionally characterize the secretory proteins of E. perseae. A total of 654 potential secretory proteins were identified, of which 190 were classified as carbohydrate-active enzymes (CAZymes), 49 as proteases, and 155 as potential effectors. A comparison to six other closely related species identified 40 species-specific putative effectors in E. perseae, indicating their specific involvement in the pathogenicity of E. perseae on avocado. The data presented in this study might be valuable for further research focused on understanding the molecular mechanisms that contribute to the pathogenicity of E. perseae on avocado.
Globally, renal cell carcinoma (RCC) is among the most common malignant tumors. It is the 6th and 10th most frequently new-diagnosed tumor in men and women, respectively [1,2]. According to the classification recommendations of the World Health Organization, renal cell carcinoma, with highly significant histological heterogeneities, is composed of up to 16 subgroups. Among them, clear cell renal cell carcinoma accounts for about 80% of all RCC types [3,4,5]. For patients presenting with early-stage ccRCC, as revealed by computed tomography, surgical excision is always the optimal treatment option with the maximal probability of being cured [6]. However, due to its asymptomatic features during the early stages and the lack of annual health examination, nearly 30% of these cases are diagnosed as metastatic ccRCC during the first consultation, implying that these patients have lost the chance of receiving radical resection [7,8]. In general, due to insensitivities, radiotherapy and chemotherapy are not optimal for patients with terminal ccRCC [9]. In recent years, various advances, particularly the emergence and upgrade of tyrosine kinase inhibitors together with immune checkpoint inhibitors, have improved the treatment of advanced unresectable ccRCC. These drugs have been shown to significantly prolong the overall survival (OS) or disease-free survival (DFS) outcomes in a subset of these patients [10,11,12]. However, due to intertumor and intratumor heterogeneity, some patients are not fully responsive or exhibit resistance to these regimens, resulting in dismal prognostic outcomes [13,14]. Hence, it is important to develop new remedies for treatment of patients with late ccRCC.
Cancer vaccine, which is aimed at activating patients' immunity against malignant tumor cells and improving the suppression status of the tumor immune micro-environment (TIME), has drawn plenty of attention from researchers worldwide [15,16]. Different types of vaccines, mainly comprising nucleic acids (DNA or mRNA), peptides, tumor cells, dendritic cells, among others, constitute the landscape of this promising remedy [17]. Messenger RNA (mRNA), encoding at least one tumor antigen, can be transported into the cytoplasm of APCs, translated into proteins, and be presented on the surface of these cells, thereby improving anti-tumor immunity [18]. The rationality of the mRNA cancer vaccine is malignant cells expressing tumor specific and tumor associated antigens. Both antigens can provide pertinent information on development of mRNA vaccines, which can generate targeted proteins with immunogenicity to induce immunoreactions, particularly in the local tumor microenvironment [19]. Tumor associated antigens are expressed at a low level in non-malignant cells but are predominantly over-expressed in malignant cells due to some certain somatic mutations, resulting in abundant presentations on the surface of cancer cells and recognition by T cells to initiate immune responses against cancer. Therefore, they are promising options for cancer vaccine development [20].
Due to technological advances, the mRNA vaccine has been more feasible as the lower innate immunogenicity, enhancements in stability and delivery system [21]. In the era of the COVID-19 pandemic, mRNA-based vaccines were confirmed to be safe, tolerable, and effective [22]. Compared to other vaccine types, mRNA-based vaccines have some unique advantages, including simultaneous activation of humoral and cell immunity, encoding of more epitopes to be presented by APCs and without the possibility of integration into the host cell genome [23]. Some completed clinical trials involving several malignant solid tumors, such as prostate cancer, breast cancer, and melanoma, have revealed the favorable application prospects of mRNA cancer vaccine monotherapy or in combination with other immunotherapies [24,25,26]. Rittig et al. reported that mRNA-based vaccines are feasible and beneficial for patients with metastatic RCC [27]. In addition, as with other immune-based treatments, the efficacies of mRNA vaccines are also, to an extent, dependent on TIME status, in particular, the infiltration levels and functions of CD8+ T cells [28].
In the present study, we investigated the potential of using ccRCC-related tumor antigens to develop mRNA-based vaccines. Elevated expressions of a tumor antigen that was positively associated with infiltration scores of APCs, were correlated with better survival outcomes. Additionally, patients were assigned into different groups in accordance with tumor immune subtypes, which can be used to guide the appropriate patient selection for receiving the mRNA vaccine. In brief, we elucidate on the development of mRNA vaccines and personalized medicines for ccRCC patients.
Raw sequencing data in counts form derived from ccRCC patients involved in the TCGA project along with their corresponding phenotypes and survival data were downloaded from the comprehensive Xena Platform [29] (https://xena.ucsc.edu/). In addition, open data in maf format, which is only available for partial ccRCC patients in the TCGA cohort (n = 336) and containing mutation-relevant information, were accessed at the Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/). Samples that lacked survival data or had survival time less than 30 days were excluded. Genes that were expressed in more than half of the samples were retained for subsequent analyses. The R package "maftools" was used to determine the mutated genes, mutation frequency and tumor mutation burden (TMB) among the ccRCC subjects. Additionally, the "somaticinteractions" function in this package was used to evaluate expression correlations among the mutated genes of interest.
The cBioPortal (v4.1.4, http://www.cbioportal.org/) for cancer genomics is an online tool incorporating multiple data from, but not limited to, the TCGA, Cancer Cell Line Encyclopedia (CCLE) and International Cancer Genome Consortium (ICGC) [30]. We used cBioPortal to explore and visualize the ccRCC genome alteration status based on 451 samples from TCGA for identification of potential tumor antigens.
The Gene Expression Profiling Interactive Analysis (GEPIA), as a form of online web server (http://gepia2.cancer-pku.cn/), integrates the raw sequencing data from samples in the TCGA database as well as Genotype-Tissue Expression project (GTEx) and has been widely used for gene expression and survival analyses of 33 cancer types [31]. Over-expressed genes in ccRCC tumor tissues, compared to normal kidney tissues, were obtained by ANOVA with two cutoff values (Log2 fold change > 1 together with q < 0.05) and were then mapped on chromosomes for visualization. All patients with survival information were divided into high-and low-expression groups based on median expression values of the genes. The Kaplan-Meier method in combination with log rank test were performed to determine whether there were significant differences in OS or DFS outcomes between groups.
In 2017, Li et al. developed the Tumor IMmune Estimation Resource (TIMER) for comprehensive analysis of tumor-infiltrating immune cells in various cancer types [32]. TIMER, a versatile web server (https://cistrome.shinyapps.io/timer/), was used to estimate the abundances of several immune cell types in cancer tissues and has multiple functional modules, including the gene module where users can determine if the expression levels of a given gene are correlated with the infiltration degrees of some immune cells. We used this module to assess the association between potential tumor antigens and abundances of three kinds of immune infiltrates, including B cells, macrophages and dendritic cells. After eliminating the impact of tumor purity, the partial Spearman's rho value was calculated, with p ≤ 0.05 denoting statistical significance.
The single cell RNA sequencing data of one ccRCC sample was downloaded from the gene expression omnibus (GEO), accession number GSM4630028 [33], and used to determine whether the aforementioned antigens were highly expressed in tumor cells. A panel of technical functions in the "Seurat" package were successively used for confirmation at single-cell resolution. Briefly, cells with 200–5000 genes and mitochondrial gene percentage < 10% were eligible for subsequent analyses. After the count data was normalized, highly varied genes were identified and subjected to principal-component analysis for dimensionality reduction. The appropriate number of principal components was acquired by the "JackStraw" function, and these components were further used to cluster cells into different types at a resolution of 0.6. Then, the "FindAllMarkers" function was used to search for marker genes in each cell type, while the ccRCC-specific marker genes (CA9 and NDUFA4L2) [34] were selected for annotation of cancer cells.
The immune-related gene list was acquired from the ImmPort database (https://www.immport.org/shared/genelists). After the removal of duplicated genes, a total of 1,793 genes remained. Expression levels of immune relevant genes in TCGA ccRCC samples were extracted to constitute a new matrix, after which the "ConsensusClusterPlus" R package [35] was used to identify the underlying immune subtypes among these patients. The partition around medoids cluster algorithm accompanied with the "1-Pearson correlation" distance was used for clustering. The total number of subsampling was set as 500, and four in five samples were simultaneously selected for each resampling. The upper limit of the resulting clusters was set as 6, and the optimal partition was visually evaluated by virtue of the consensus matrix combined with the consensus cumulative distribution function. Then, principal component analysis was performed to assess if samples from different subtypes could be well distinguished. Inter-group differences in overall survival time of patients were determined using the "survival" R package. Correlations between immune subtypes and clinical features, molecular subtypes, tumor mutation burden were explored to describe these intrinsic features among the defined immune subtypes.
Stromal and immune scores for each sample were calculated using the "Estimate" algorithm [36]. Then, on the basis of a gene list containing twenty-eight immune cell types, the single sample gene set enrichment analysis (ssGSEA) method in the "GSVA" package was used to determine the abundance of these cell types [37]. Lastly, based on gene expression characteristics of twenty-two kinds of immune cells, the "CIBERSORT" algorithm was used to calculate the proportion of these cells in samples [38].
Two individual gene sets, encompassing immunogenic cell death (ICD) modulators and immune checkpoints (ICP), were retrieved from previously published articles [39,40]. Then, expressions of these genes in diverse immune subgroups were extracted and compared by performing the Wilcoxon rank sum test.
The co-expression modules were identified and correlated with the immune subtypes using the WGCNA package in R software [41]. Non- or low-varying genes, evaluated using the median absolute deviation (MAD), were filtered out. Thereafter, a sample tree was plotted to detect whether the samples had outliers. Then, the optimal soft threshold was calculated per the inputted gene expression profile and was served as an actual argument of the "blockwiseModules" function which was a one-step method to construct the co-expression matrix. Subsequently, the genes were clustered into different color modules based on the topological overlap matrix (TOM). After that, the correlations between modules and the immune subtypes were explored. Lastly, gene ontology (GO) analysis was carried out using the "clusterProfiler" package [42] to annotate the cellular location and biological functions of the module genes that were strongly correlated with the immune subtypes.
Genomic alterations are known to influence one's response to anticancer drugs. In this study, the Genomics of Drug Sensitivity in Cancer (GDSC) database was used to provide the gene expression profile of numerous cancer cell lines and half maximal inhibitory concentration (IC50) of hundreds of common compounds against these cell strains [43]. Using datasets from GDSC as a training set, the R package oncoPredict [44] can predict sensitivity to anticancer agents via constructing a ridge regression model based on the gene expression matrix obtained from clinical samples of tumor tissues. Drug-sensitivity analysis was investigated in patients with different immune subtypes.
To identify potential tumor antigens for ccRCC, 1,638 over-expressed genes with the potential to encode tumor associated antigens (TAA) were determined via differential gene-expression profiling between normal and malignant renal parenchyma tissues (Figure 1A). Variations in tumor mutation counts and altered genome fractions in ccRCC cases were analyzed to screen the TAA-encoding mutated genes, resulting in 11,686 mutation genes (Figure 1B, C). Mutational analysis revealed that GAS8 antisense RNA 1 (GAS8-AS1), von Hippel-Lindau tumor suppressor (VHL), mucin 16 (MUC16), ankyrin repeat, and PH domain 3 (ARAP3), heme binding protein 1 (HEBP1), phosphatase and actin regulator 1 (PHACTR1), SSX family member 3 (SSX3), teashirt zinc finger homeobox 3 (TSHZ3), ATP binding cassette subfamily A member 6 (ABCA6), as well as CD4 molecule (CD4), were the top 10 highly mutated genes in the mutation count group (Figure 1D). The top 10 highly altered genes in the fraction genome altered group included von Hippel-Lindau tumor suppressor (VHL), inhibitor of nuclear factor kappa B kinase subunit beta (IKBKB), transforming acidic coiled-coil containing protein 1 (TACC1), neuregulin 1 (NRG1), protein O-mannose kinase (POMK), pleckstrin and Sec7 domain containing 3 (PSD3), polybromo 1 (PBRM1), adrenoceptor beta 3 (ADRB3), adaptor related protein complex 3 subunit mu 2 (AP3M2) and cholinergic receptor nicotinic beta 3 subunit (CHRNB3) (Figure 1E). Taken together, 629 overexpressed and mutated genes were identified for further analysis.
Based on the afore-mentioned intersecting result, genes with prognostic values were screened for development of an mRNA vaccine. Six and seven genes were found to be respectively correlated with OS and DFS of ccRCC patients, and finally, two common genes, LDL receptor related protein 2 (LRP2) and dedicator of cytokinesis 8 (DOCK8) were identified (Figure 2A). Compared to the low expression group, patients in the high LRP2 expression group exhibited a superior survival time (Figure 2B, C), while patients with high expressions of DOCK8 were also associated with better survival outcomes (Figure 2D, E). These findings suggest that the two genes may have immune-stimulatory functions. Expressions of LRP2 and DOCK8 were positively correlated with the abundance of macrophages, dendritic cells and B cells (Figure 3A, B). Then, expressions of LRP2 and DOCK8 in different stages were determined (Figure 3C, D) using the UALCAN database [45]. It was established that with progressing stages, expressions of the two genes gradually decreased, implying that mRNA vaccines, based on the two potential tumor antigens, may yield better immunogenicity in early-stage ccRCC patients. The two genes (LRP2 and DOCK8) with potential immune-stimulating effects are potential candidates for ccRCC-mRNA vaccine development.
After quality control, 21,549 genes and 5666 cells remained, and 20 principal components were used to cluster these cells, leading to various cell types (Figure 4A). Cell clusters (0, 1, 12) with simultaneously higher expressions of CA9 and NDUFA4L2 were annotated as ccRCC cells (Figure 4B). In Figure 4C, D, LRP2 was mainly expressed in tumor cells, but DOCK8 was not, hence only LRP2 could be regarded as an authentic tumor antigen for development of an mRNA vaccine.
Immune-typing can be used to reflect, to a great extent, the immune status in heterogeneous tumor micro-environments, which informs the selection of patients to be administered with the mRNA vaccine. Therefore, consensus clustering, according to expression levels of immune-related genes, was performed in TCGA ccRCC samples. Considering the results from delta area, consensus CDF and consensus matrix, the optimal cluster number (n = 2) was determined, and accordingly, two robust immune subtypes (IS1 and IS2) were identified (Figure 5A–C). Survival analysis revealed significant differences between the two groups, with patients in the IS2 subtype having better OS outcomes than those in the IS1 subtype (Figure 5D). In addition, principal components analysis (PCA) was performed to assess subtype assignments, which revealed that most of the samples in the IS1 group were distinct from samples in the IS2 group (Figure 5E). The distribution ratio of immune subtypes identified herein in six well-known pan-cancer immune subtypes [46], encompassing C1 (Would Healing), C2 (IFN-γ Dominant), C3 (Inflammatory), C4 (Lymphocyte Depleted), C5 (Immunologically Quiet) as well as C6 (TGF-β), were explored, showing a high abundance of IS2 in C4 and C5 subtypes (Figure 5F). The distribution ratios of IS1 and IS2 in different grades and stages of ccRCC patients were investigated (Figure 5G, H). Lower grades and stages were associated with IS2 while higher grades and stages were associated with IS1. Expressions of LRP2 were significantly high in the IS2 group (Figure 5I). Taken together, immunophenotyping is a potential robust indicator for assessment of the clinical outcomes of ccRCC patients.
Tumor mutation burden and somatic mutation rates play a vital role in anti-cancer immunity, and generally, higher TMB levels are associated with stronger immunogenicity, leading to increased immunotherapeutic efficacies [47]. Thus, TMB and mutations of the two subtypes were determined and compared. The IS1 exhibited higher TMB and number of mutated genes but without statistical significance (Figure 6A, B). Further, highly mutated genes, such as VHL (47%) and PBRM1 (40%), across the two immune subtypes in ccRCC patients were revealed by somatic mutation analysis (Figure 6C). Then, expression correlations of the frequently mutated genes were explored (Figure 6D). There were no marked differences in immune subtypes with regard to the magnitude of mutation burden.
During ccRCC progression, a dysregulated TIME plays a pivotal role and largely influences the outcomes of mRNA vaccine therapy. First, based on the ESTIMATE method, immune and stromal scores were directly obtained from each included sample, and tumor purity was calculated as: Tumor purity = cos (0.6049872018 + 0.0001467884 × ESTIMATE score) [36]. In Figure 7A–C, compared with the IS2 group, the IS1 group had significantly high immune and stromal scores, which implied a lower tumor purity. Then, ssGSEA was used to assess the abundance of 28 immune cell types in the two immune groups. Differences in abundance of immune cells between the two groups were significant (Figure 7D). Compared to IS2, IS1 exhibited a significantly high abundance of activated B cells, type 1 T helper cells, myeloid-derived suppressor cells (MDSC), activated CD4 T cells, macrophages, activated CD8 T cells and regulatory T cells among others. CIBERSORT analysis (Figure 7E) revealed that compared with IS2, IS1 samples had significant high proportions of CD8 T cells, regulatory T cells, plasma cells and M0 macrophages with low proportions of M2 macrophages, CD4 memory resting T cells, resting mast cells and monocytes.
Both ICP and ICD modulators have important effects on modulation of host anti-tumor immunity, thereby affecting the immunotherapeutic efficacy, including that of the mRNA vaccine. Therefore, differential expressions of the above immune-related molecules were determined in these immune subtypes. A total of 47 ICP were detected in the TCGA ccRCC cohort and a large proportion of these genes (more than 30), including PDCD1, CD274 and CTLA4, were significantly highly expressed in the IS1 subtype with only two genes (TNFRSF14, TNFRSF4) being over-expressed in the IS2 subtype (Figure 8A). Moreover, 13 ICD modulator-related genes, including CXCL10, FRP1, IFNAR2, IFNE, CALR, LRP1, PANX1, P2RY2, ANXA1, P2RX7, HGF, MET and TLR4, were significantly elevated in IS1 tumors, while EIF2A, IFNAR1 and HMGB1 were up-regulated in IS2 tumors (Figure 8B). Therefore, immunotyping is a potential indicator of the efficacy the mRNA vaccine as the subtypes enable comprehensive assessment of the expressions of ICP and ICD modulators. Moreover, effectiveness of the mRNA vaccine can be damaged in the context of over-expressed ICP and under-expressed ICD modulators.
The classification of high-varying genes, using gene co-expression modules, was achieved and correlated with different immune subtypes. Based on the results of mean connectivity coupled with scale-free fit index, five was set as the soft-thresholding power for subsequent analyses (Figure 9A). Then, 10 modules (distinguished by distinct colors) were identified via average linkage hierarchical clustering (Figure 9B) and the number of genes in each module established (Figure 9C). Next, correlation analysis of module to trait (that is, in this situation, immune subtypes) was performed. According to the correlation coefficient and p-value (Figure 9D), the three mostly positively correlated modules for IS1 were the black (MEblack: rho: 0.74, p < 0.01), yellow (MEyellow: rho: 0.53, p < 0.01) and grey (MEgrey: rho: 0.39, p < 0.01) modules, while the blue (MEblue: rho: -0.18, p < 0.01) and green (MEgreen: rho: -0.28, p < 0.01) modules were negatively correlated with IS1. Module membership vs. gene significance scores of genes in the MEblack were visualized in the scatter plot (Figure 9E), which confirmed the highly positive relationship between the black module and IS1 (r = 0.89, p < 0.01). GO enrichment analysis of genes extracted from the black module suggested that these genes were mainly involved in the B cell receptor signaling pathway, phagocytosis and complement activation, among others (Figure 9F).
To determine whether patients with the immune subtype that was less suitable for mRNA vaccination could be more sensitive to other clinically available drugs, anti-cancer drug sensitivity for each sample in the TCGA-ccRCC cohort was predicted with the magnitude of the IC50 value (Figure 10A–F). Surprisingly, most of the drugs, including sunitinib, cabozantinib, temsirolimus, and tivozanib, were fitter for patients in the IS1 cluster, where they exhibited a relatively low IC50. There was a slight diversity with regard to the sensitivity of sorafenib and pazopanib between the two subtypes, with IS2 being more sensitive to sorafenib than IS1. Taken together, although IS1 patients may be less suitable for receiving the mRNA vaccine, they are more suitable for targeted therapy involving drugs that target the vascular endothelial growth factor receptor (VEGFR) and mammalian target of rapamycin (mTOR). The diversity in treatment sensitivity between the two subtypes also confirmed the heterogeneity among patients, which reflects the urgency of individualized treatments.
Among the different renal cancer types, ccRCC is the most common and is highly aggressive. For patients with early stage diagnosis, the 5-year survival rate after reasonable therapies is more than 90% and may reduce to approximately 12% due to concomitant metastasis [48]. Elucidation of the genesis and progression of ccRCC have enabled the development of molecularly targeted agents for applications in advanced patients. For instance, sunitinib, everolimus and nivolumab have been proven to improve the survival rates in partial patients [49]. However, the efficacy of these drugs is always limited or temporary and many patients ultimately die of metastases. Therefore, there is a need to develop new alternative or combinational therapeutics.
Tumor antigen supplemented with or without an adjuvant is the core of a representative cancer vaccine, being responsible for supporting immune cells to recognize and eliminate cancer cells [50]. Traditional chemotherapy, molecular targeted therapy and immunotherapy are associated with more development time and cost, recurrent drug resistance or side-effects and limited therapeutic potency; in contrast, the cancer vaccine is more promising for its' maximal specific effects, lower toxicity and long-lasting immunological memory [51]. Given these merits, three types of cancer vaccines against ccRCC embodying autologous tumor cells, dendritic cells and peptides have been developed and evaluated in clinical trials. Jocham et al. reported that patients in the autologous renal tumor cell vaccine group (n = 177) had a significantly high 5-year progression-free survival rate than the control group (no vaccine, n = 202), and only 12 treatment-associated adverse events were observed [52]. IMA901, which contains multiple tumor-related peptides, is a therapeutic vaccine that can induce T cell responses in patients with metastatic RCC, resulting in better disease control and longer overall survival outcomes [53]. AGS-003 is prepared using tumor RNA and dendritic cells derived from autologous monocytes of metastatic ccRCC patients. Compared with sunitinib monotherapy, the combination of AGS-003 and sunitinib was associated with prolonged survival outcomes [54]. Despite these results, none of them are currently available for routine use in clinics, suggesting the need for development of other kinds of vaccines. The mRNA-based vaccine platform has more novelty with the potential to be a major breakthrough for several reasons. First, it can be designed to synchronously deliver multiple tumor-associated antigens and induce antibody-mediated humoral as well as cell-mediated immune responses, decreasing the possibility of vaccine resistance [55]. Second, being different from the peptide vaccine, it can encode whole-length tumor antigens, which enables APCs to simultaneously present more epitopes, increasing the likelihood of a wide T cell response [56]. Last, it is non-infectious and free of proteins or contaminations in the production process, implying better tolerability and safety [57]. Generally, the open reading frame encoding vaccine antigens and structural components, including the 5'-cap, 3'-tail and two untranslated regions form the core of a typical mRNA vaccine. By lowering the immunogenicity and enhancing the stability as well as delivery efficiency of mRNA vaccines, their feasibility can be further increased. For instance, the immunogenicity of mRNA can be reduced by chemical modifications of nucleotides, while its stability can be improved via the addition of AU- or GC-rich sequences [58,59]. Studies have reported tumor regression after vaccination alone or in combination with other therapies, such as immune checkpoint blockers or chemotherapeutics [60]. However, regarding ccRCC, mRNA vaccination is still in its infancy and further investigations should aim at mining tumor antigens with superior immunogenicity.
In this study, various preliminary targeting antigens were revealed by analyses of somatic mutations and over-accumulated genes of ccRCC. Given that these antigens were predicted by altered gene profiles, and may not have significant functions in ccRCC, their prognostic value and correlations with APCs infiltrations were further investigated. Then, expressions of these potential antigens were further investigated at the single-cell level. Finally, one tumor antigen (LRP2), which was positively correlated with both survival time and richness of APCs, was identified to be a promising candidate for vaccine development. Based on the above analyses, the mRNA sequence of LRP2 can be chemically modified and synthesized using in vitro transcription (IVT) systems. After encapsulation with lipid nanoparticles, the optimized LRP2-encoding mRNA vaccine can be injected into ccRCC patients to induce immune responses against cancer cells. The potentiality of this tumor antigen for development of an anti-ccRCC mRNA vaccine has been reported. For example, LRP2 was highly expressed in malignant melanoma tissues, compared to benign nevi, and siRNA-mediated suppression of LRP2 expression significantly inhibited melanoma cell proliferation and survival [61]. Metformin decreases the proliferative abilities of human papillary thyroid cancer cells by inducing the apoptosis of TPC-1 cells via the downregulation of LRP2 [62]. In addition, LRP2 expressions gradually increased from premalignant oral leukoplakia to oral squamous cell carcinoma and were associated with metallothioneins, indicating the onco-driving potential of this interaction [63].
Only a small number of patients that have adopted cancer vaccine-based regimens have exhibited increased objective response rates and prolonged survival time; therefore, ccRCC patients were stratified according to the expression status of immune related genes to achieve appropriate uses of cancer vaccine therapy. Two immune subtypes were finally identified, with distinct clinical and molecular characteristics. Compared to patients in the IS2 cluster, patients within IS1 exhibited significantly low overall survival outcomes, indicating that immunotype is a robust prognostic indicator for ccRCC. Apart from prognosis, the abundance of the two immune subtypes was compared by grade and stage, revealing that patients with advanced ccRCC accounted for a greater percentage of IS1. A higher IS2 proportion was found in the immunological "cold" subtypes (C4 and C5), and administration of an mRNA vaccine in IS2 patients may induce host immune responses towards a "hot" state. Furthermore, compared with IS1, IS2 exhibited higher expressions of LRP2, suggesting that mRNA vaccines may confer more benefits to patients with IS2. Intriguingly, the two subtypes manifested significantly different sensitivities to several drugs that were routinely used in clinical settings, providing a reference for patients who were not suitable for vaccination and emphasizing the necessity of personalized treatment.
Various immune cells and multiple immune-associated molecules, such as ICP and ICD modulators, are vital components of the TIME that can greatly influence the outcomes of the mRNA vaccine. IS1 had higher immune scores than IS2, indicating the immune "hot" status of IS1, and was further validated by widely elevated infiltrations of different immune cells. Notably, expressions of cells with immune-suppressive abilities, especially regulatory T cells and myeloid-derived suppressor cells, were enhanced in IS1, thereby inhibiting the functions of immune effector cells by various mechanisms [64]. The majority of ICP, including PD-1, PD-L1 and CTLA4, were highly expressed in IS1, suggesting that IS1 was dominated by an immune-suppressive microenvironment, which can inhibit the performance of ICD modulators [65]. Taken together, IS1 was characterized by elevated expressions of multiple immune-suppressive cells and molecules that can compromise the effectiveness of mRNA vaccines, therefore, combination or sequential therapies with immune checkpoint blockers should be considered.
For cancer, the magnitude of the immunogenicity (in other words, the capacity for eliciting adaptive immune responses) is dependent on tumor antigens, infiltration degree of immune cells and the presence of immunomodulatory molecules in the tumor microenvironment [66]. Although it has a relatively lower tumor mutational burden, ccRCC has been described as an immunogenic cancer since survival benefits were observed, more than two decades ago, in patients treated with interleukin-2 and interferon-α. The subsequent success of immune checkpoint inhibitors further consolidated this concept [67]. Hence, mRNA-type cancer vaccine, as an alternative or assistant immunotherapy, may have broad application prospects in ccRCC.
In conclusion, LRP2 is a potential target for the development of a ccRCC mRNA vaccine, and patients in the IS2 cluster may highly benefit from this therapy. Additionally, this study provides some insights for the development of mRNA vaccines against ccRCC and the selection of suitable vaccine-receiving patients.
This study was supported by the Key Research and Development Program of Guangxi (2021AB13014), Major Project of Guangxi Innovation Driven (AA18118016), and Guangxi Key Laboratory for Genomic and Personalized Medicine [grant number 16-380-54, 17-259-45, 19-050-22, 19-185-33, 20-065-33].
The authors declare that they have no conflicts of interest.
[1] |
Dean R, Van Kan JAL, Pretorius ZA, et al. (2012) The Top 10 fungal pathogens in molecular plant pathology. Mol Plant Pathol 13: 414-430. https://doi.org/10.1111/j.1364-3703.2011.00783.x ![]() |
[2] | Biju VC, Fokkens L, Houterman PM, et al. (2017) Multiple evolutionary trajectories have led to the emergence of races in Fusarium oxysporum f. sp. lycopersici. Appl Environ Microbiol 83. https://doi.org/10.1128/AEM.02548-16 |
[3] |
Singh RP, Hodson DP, Jin Y, et al. (2015) Emergence and spread of new races of wheat stem rust fungus: continued threat to food security and prospects of genetic control. Phytopathology 105: 872-884. https://doi.org/10.1094/PHYTO-01-15-0030-FI ![]() |
[4] |
Pontes JG de M, Fernandes LS, dos Santos R Vander, et al. (2020) Virulence factors in the phytopathogen–host interactions: an overview. J Agric Food Chem 68: 7555-7570. https://doi.org/10.1021/acs.jafc.0c02389 ![]() |
[5] |
Jaswal R, Kiran K, Rajarammohan S, et al. (2020) Effector biology of biotrophic plant fungal pathogens: current advances and future prospects. Microbiol Res 241: 126567. https://doi.org/10.1016/j.micres.2020.126567 ![]() |
[6] |
Stergiopoulos I, de Wit PJGM (2009) Fungal effector proteins. Annu Rev Phytopathol 47: 233-263. https://doi.org/10.1146/annurev.phyto.112408.132637 ![]() |
[7] |
Pradhan A, Ghosh S, Sahoo D, et al. (2021) Fungal effectors, the double edge sword of phytopathogens. Curr Genet 67: 27-40. https://doi.org/10.1007/s00294-020-01118-3 ![]() |
[8] |
Jeffress S, Arun-Chinnappa K, Stodart B, et al. (2020) Genome mining of the citrus pathogen Elsinoë fawcettii; prediction and prioritisation of candidate effectors, cell wall degrading enzymes and secondary metabolite gene clusters. PLoS One 15: e0227396. https://doi.org/10.1371/journal.pone.0227396 ![]() |
[9] | Tan KLS, Mohamad SB (2024) Fungal pathogen in digital age: review on current state and trend of comparative genomics studies of pathogenic fungi. Adv Microbiol 63: 23-31. https://doi.org/10.2478/am-2024-0003 |
[10] | Azeez SO, Adeboye SE (2024) Advances in understanding plant-pathogen interactions: insights from tomato as a model system. VirusDis . https://doi.org/10.1007/s13337-024-00889-4 |
[11] |
Soanes DM, Richards TA, Talbot NJ (2007) Insights from sequencing fungal and oomycete genomes: what can we learn about plant disease and the evolution of pathogenicity?. Plant Cell 19: 3318-3326. https://doi.org/10.1105/tpc.107.056663 ![]() |
[12] |
Wang Y, Wu J, Yan J, et al. (2022) Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies. BMC Genomics 23: 34. https://doi.org/10.1186/s12864-021-08165-1 ![]() |
[13] | Li J, Ai M, Hou J, et al. (2024) Plant–pathogen interaction with root rot of Panax notoginseng as a model: Insight into pathogen pathogenesis, plant defence response and biological control. Mol Plant Pathol 25. https://doi.org/10.1111/mpp.13427 |
[14] | Pernezny K, Marlatt RB Diseases of Avocado in Florida, EDIS Publication PP21 (2007). |
[15] | Belizaire CM, Gañán-Betancur L, Gazis R (2023) Avocado scab caused by Elsinoe perseae: A diagnostic guide. Plant Health Prog . https://doi.org/10.1094/PHP-10-23-0084-DG |
[16] | Jenkins AE (1934) Sphaceloma Perseae, the cause of avocado scab. J Agric Res 49: 859-869. |
[17] | Jenkins AE (1934) A species of Sphaceloma on Avocado. Phytopathology 24. |
[18] |
Everett KR, Rees-George J, Pushparajah IPS, et al. (2011) Molecular identification of sphaceloma perseae (avocado scab) and its absence in New Zealand. J Phytopathol 159: 106-113. https://doi.org/10.1111/j.1439-0434.2010.01739.x ![]() |
[19] |
Belizaire CM, Gañán-Betancur L, Gazis R (2024) Avocado scab caused by Elsinoe perseae: a diagnostic guide. Plant Health Prog 25: 218-225. https://doi.org/10.1094/PHP-10-23-0084-DG ![]() |
[20] |
Everett KR, Stevens PS, Cutting JGM (1999) Postharvest fruit rots of avocado are reduced by benomyl applications during flowering. Proceedings of the New Zealand Plant Protection Conference 52: 153-156. https://doi.org/10.30843/nzpp.1999.52.11598 ![]() |
[21] | Everett KR, Owen SG, Cutting JGM (2005) Testing efficacy of fungicides against postharvest pathogens of avocado (Persea americana cv Hass). N Z Plant Prot 58: 89-95. https://doi.org/10.30843/nzpp.2005.58.4260 |
[22] |
Sumida CH, Fantin LH, Braga K, et al. (2020) Control of root rot (Phytophthora cinnamomi) in avocado (Persea Americana) with bioagents. Summa Phytopathol 46: 205-211. https://doi.org/10.1590/0100-5405/192195 ![]() |
[23] |
Guarnizo N, Álvarez A, Oliveros D, et al. (2022) Elicitor activity of curdlan and its potential application in protection of hass avocado plants against Phytophthora cinnamomi rands. Horticulturae 8: 646. https://doi.org/10.3390/horticulturae8070646 ![]() |
[24] |
Gañán-Betancur L, Gazis R (2023) Genome sequence resource of the avocado scab pathogen Elsinoe perseae. Microbiol Resour Announc 12. https://doi.org/10.1128/mra.00190-23 ![]() |
[25] |
Humann JL, Lee T, Ficklin S, et al. (2019) Structural and functional annotation of eukaryotic genomes with GenSAS. Methods Mol Biol : 29-51. https://doi.org/10.1007/978-1-4939-9173-0_3 ![]() |
[26] |
Li Z, Wang Y, Fan Y, et al. (2021) Transcriptome analysis of the grape- Elsinoë ampelina pathosystem reveals novel effectors and a robust defense response. Mol Plant-Microbe Interact 34: 110-121. https://doi.org/10.1094/MPMI-08-20-0227-R ![]() |
[27] |
Yun HK, Louime C, Lu J (2007) First report of anthracnose caused by Elsinoë ampelina on muscadine grapes (Vitis rotundifolia) in Northern Florida. Plant Dis 91: 905-905. https://doi.org/10.1094/PDIS-91-7-0905B ![]() |
[28] |
Stanke M, Keller O, Gunduz I, et al. (2006) AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res 34: W435-W439. https://doi.org/10.1093/nar/gkl200 ![]() |
[29] | Borodovsky M, Lomsadze A (2011) Eukaryotic Gene Prediction Using GeneMark.hmm-E and GeneMark-ES. Curr Protoc Bioinformatics 35. https://doi.org/10.1002/0471250953.bi0406s35 |
[30] | Haas BJ, Salzberg SL, Zhu W, et al. (2008) Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol 9. https://doi.org/10.1186/gb-2008-9-1-r7 |
[31] |
Li Z, Fan Y, Chang P, et al. (2020) Genome sequence resource for Elsinoë ampelina , the causal organism of grapevine anthracnose. Mol Plant-Microbe Interact 33: 576-579. https://doi.org/10.1094/MPMI-12-19-0337-A ![]() |
[32] |
Shanmugam G, Jeon J, Hyun J-W (2020) Draft genome sequences of Elsinoë fawcettii and Elsinoë australis causing scab diseases on citrus. Mol Plant-Microbe Interact 33: 135-137. https://doi.org/10.1094/MPMI-06-19-0169-A ![]() |
[33] |
Su J, Liu J, Hu Y, et al. (2022) High-quality genome sequence resource of Elsinoë arachidis strain LY-HS-1, causing scab disease of peanut. Plant Dis 106: 1506-1509. https://doi.org/10.1094/PDIS-11-21-2549-A ![]() |
[34] |
Zhang X, Zou H, Yang Y, et al. (2022) Genome resource for Elsinoë batatas, the causal agent of stem and foliage scab disease of sweet potato. Phytopathology 112: 973-975. https://doi.org/10.1094/PHYTO-08-21-0344-A ![]() |
[35] |
Wingfield BD, Berger DK, Coetzee MPA, et al. (2022) IMA genome-F17. IMA Fungus 13: 19. https://doi.org/10.1186/s43008-022-00104-3 ![]() |
[36] |
Chellappan BV, El-Ganainy SM, Alrajeh HS, et al. (2023) In silico characterization of the secretome of the fungal pathogen thielaviopsis punctulata, the causal agent of date palm black scorch disease. J Fungi 9: 303. https://doi.org/10.3390/jof9030303 ![]() |
[37] |
Teufel F, Almagro Armenteros JJ, Johansen AR, et al. (2022) SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat Biotechnol 40: 1023-1025. https://doi.org/10.1038/s41587-021-01156-3 ![]() |
[38] |
Käll L, Krogh A, Sonnhammer ELL (2004) A combined transmembrane topology and signal peptide prediction method. J Mol Biol 338: 1027-1036. https://doi.org/10.1016/j.jmb.2004.03.016 ![]() |
[39] | Hallgren J, Tsirigos KD, Pedersen MD, et al. (2022) DeepTMHMM predicts alpha and beta transmembrane proteins using deep neural networks. bioRxiv . https://doi.org/10.1101/2022.04.08.487609 |
[40] |
de Castro E, Sigrist CJ a, Gattiker A, et al. (2006) ScanProsite: Detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res 34: 362-365. https://doi.org/10.1093/nar/gkl124 ![]() |
[41] |
Almagro Armenteros JJ, Salvatore M, Emanuelsson O, et al. (2019) Detecting sequence signals in targeting peptides using deep learning. Life Sci Alliance 2: e201900429. https://doi.org/10.26508/lsa.201900429 ![]() |
[42] |
Horton P, Park K-J, Obayashi T, et al. (2007) WoLF PSORT: protein localization predictor. Nucleic Acids Res 35: W585-W587. https://doi.org/10.1093/nar/gkm259 ![]() |
[43] |
Gíslason MH, Nielsen H, Almagro Armenteros JJ, et al. (2021) Prediction of GPI-anchored proteins with pointer neural networks. Curr Res Biotechnol 3: 6-13. https://doi.org/10.1016/j.crbiot.2021.01.001 ![]() |
[44] | Finn RD, Bateman A, Clements J, et al. (2014) Pfam: The protein families database. Nucleic Acids Res 42. https://doi.org/10.1093/nar/gkt1223 |
[45] |
Zdobnov EM, Apweiler R (2001) InterProScan-An integration platform for the signature-recognition methods in InterPro. Bioinformatics 17: 847-848. https://doi.org/10.1093/bioinformatics/17.9.847 ![]() |
[46] |
Cantarel BL, Coutinho PM, Rancurel C, et al. (2009) The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res 37: D233-D238. https://doi.org/10.1093/nar/gkn663 ![]() |
[47] |
Yin Y, Mao X, Yang J, et al. (2012) dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res 40: W445-W451. https://doi.org/10.1093/nar/gks479 ![]() |
[48] | Sperschneider J, Gardiner DM, Dodds PN, et al. (2016) EffectorP: Predicting fungal effector proteins from secretomes using machine learning. New Phytologist . https://doi.org/10.1111/nph.13794 |
[49] |
Winnenburg R (2006) PHI-base: a new database for pathogen host interactions. Nucleic Acids Res 34: D459-D464. https://doi.org/10.1093/nar/gkj047 ![]() |
[50] |
Rawlings ND, Barrett AJ, Thomas PD, et al. (2018) The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res 46: D624-D632. https://doi.org/10.1093/nar/gkx1134 ![]() |
[51] |
Sun J, Lu F, Luo Y, et al. (2023) OrthoVenn3: an integrated platform for exploring and visualizing orthologous data across genomes. Nucleic Acids Res 51: W397-W403. https://doi.org/10.1093/nar/gkad313 ![]() |
[52] | Morales-García JL, López-Cornejo CI, Pedraza-Santos ME, et al. (2023) Morpho-molecular identification of the causal agent of avocado scab in Michoacán. Mexican J Phytopathol 41. https://doi.org/10.18781/R.MEX.FIT.2302-4 |
[53] | Gavande PV, Goyal A, Fontes CMGA (2023) Carbohydrates and Carbohydrate-Active enZymes (CAZyme): An overview. Glycoside Hydrolases.Elsevier 1-23. https://doi.org/10.1016/B978-0-323-91805-3.00012-5 |
[54] |
Benini S (2020) Carbohydrate-active enzymes: structure, activity, and reaction products. Int J Mol Sci 21: 2727. https://doi.org/10.3390/ijms21082727 ![]() |
[55] |
Defilippi BG, Ejsmentewicz T, Covarrubias MP, et al. (2018) Changes in cell wall pectins and their relation to postharvest mesocarp softening of “Hass” avocados (Persea americana Mill.). Plant Physiol Biochem 128: 142-151. https://doi.org/10.1016/j.plaphy.2018.05.018 ![]() |
[56] |
Platt-Aloia KA, Thomson WW, Young RE (1980) Ultrastructural changes in the walls of ripening avocados: transmission, scanning, and freeze fracture microscopy. Botanical Gazette 141: 366-373. https://doi.org/10.1086/337169 ![]() |
[57] |
Kubicek CP, Starr TL, Glass NL (2014) Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol 52: 427-451. https://doi.org/10.1146/annurev-phyto-102313-045831 ![]() |
[58] |
Lairson LL, Henrissat B, Davies GJ, et al. (2008) Glycosyltransferases: structures, functions, and mechanisms. Annu Rev Biochem 77: 521-555. https://doi.org/10.1146/annurev.biochem.76.061005.092322 ![]() |
[59] |
Horn SJ, Vaaje-Kolstad G, Westereng B, et al. (2012) Novel enzymes for the degradation of cellulose. Biotechnol Biofuels 5: 45. https://doi.org/10.1186/1754-6834-5-45 ![]() |
[60] |
Phillips CM, Beeson WT, Cate JH, et al. (2011) Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa. ACS Chem Biol 6: 1399-1406. https://doi.org/10.1021/cb200351y ![]() |
[61] |
Villares A, Moreau C, Bennati-Granier C, et al. (2017) Lytic polysaccharide monooxygenases disrupt the cellulose fibers structure. Sci Rep 7: 40262. https://doi.org/10.1038/srep40262 ![]() |
[62] |
Scheller HV, Ulvskov P (2010) Hemicelluloses. Annu Rev Plant Biol 61: 263-289. https://doi.org/10.1146/annurev-arplant-042809-112315 ![]() |
[63] |
van den Brink J, de Vries RP (2011) Fungal enzyme sets for plant polysaccharide degradation. Appl Microbiol Biotechnol 91: 1477-1492. https://doi.org/10.1007/s00253-011-3473-2 ![]() |
[64] |
Kubicek CP, Starr TL, Glass NL (2014) Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol 52: 427-451. https://doi.org/10.1146/annurev-phyto-102313-045831 ![]() |
[65] |
Wilkens C, Andersen S, Dumon C, et al. (2017) GH62 arabinofuranosidases: Structure, function and applications. Biotechnol Adv 35: 792-804. https://doi.org/10.1016/j.biotechadv.2017.06.005 ![]() |
[66] |
Chauhan PS, Gupta N (2017) Insight into microbial mannosidases: a review. Crit Rev Biotechnol 37: 190-201. https://doi.org/10.3109/07388551.2015.1128878 ![]() |
[67] |
Garron M-L, Cygler M (2014) Uronic polysaccharide degrading enzymes. Curr Opin Struct Biol 28: 87-95. https://doi.org/10.1016/j.sbi.2014.07.012 ![]() |
[68] |
Husain Q (2010) β Galactosidases and their potential applications: a review. Crit Rev Biotechnol 30: 41-62. https://doi.org/10.3109/07388550903330497 ![]() |
[69] |
Prade RA (1996) Xylanases: from Biology to BioTechnology. Biotechnol Genet Eng Rev 13: 101-132. https://doi.org/10.1080/02648725.1996.10647925 ![]() |
[70] |
Flutto L (2003) PECTIN|Properties and determination. Encyclopedia of Food Sciences and Nutrition.Elsevier 4440-4449. https://doi.org/10.1016/B0-12-227055-X/00901-9 ![]() |
[71] |
Willats WG, McCartney L, Mackie W, et al. (2001) Pectin: cell biology and prospects for functional analysis. Plant Mol Biol 47: 9-27. https://doi.org/10.1023/A:1010662911148 ![]() |
[72] |
Li J, Peng C, Mao A, et al. (2024) An overview of microbial enzymatic approaches for pectin degradation. Int J Biol Macromol 254: 127804. https://doi.org/10.1016/j.ijbiomac.2023.127804 ![]() |
[73] |
Arya GC, Cohen H (2022) The multifaceted roles of fungal cutinases during infection. J Fungi 8: 199. https://doi.org/10.3390/jof8020199 ![]() |
[74] |
Jashni MK, Mehrabi R, Collemare J, et al. (2015) The battle in the apoplast: further insights into the roles of proteases and their inhibitors in plant–pathogen interactions. Front Plant Sci 6. https://doi.org/10.3389/fpls.2015.00584 ![]() |
[75] |
Todd JNA, Carreón-Anguiano KG, Islas-Flores I, et al. (2022) Fungal effectoromics: a world in constant evolution. Int J Mol Sci 23: 13433. https://doi.org/10.3390/ijms232113433 ![]() |
[76] |
Liu W, Liu J, Ning Y, et al. (2013) Recent progress in understanding pamp- and effector-triggered immunity against the rice blast fungus magnaporthe oryzae. Mol Plant 6: 605-620. https://doi.org/10.1093/mp/sst015 ![]() |
[77] |
Moreno-Sánchez I, Pejenaute-Ochoa MD, Navarrete B, et al. (2021) Ustilago maydis secreted endo-xylanases are involved in fungal filamentation and proliferation on and inside plants. J Fungi 7: 1081. https://doi.org/10.3390/jof7121081 ![]() |
[78] |
Chen S, Songkumarn P, Venu RC, et al. (2013) Identification and characterization of in planta–expressed secreted effector proteins from Magnaporthe oryzae that induce cell death in rice. Mol Plant-Microbe Interact 26: 191-202. https://doi.org/10.1094/MPMI-05-12-0117-R ![]() |
[79] |
Polonio Á, Fernández-Ortuño D, de Vicente A, et al. (2021) A haustorial-expressed lytic polysaccharide monooxygenase from the cucurbit powdery mildew pathogen Podosphaera xanthii contributes to the suppression of chitin-triggered immunity. Mol Plant Pathol 22: 580-601. https://doi.org/10.1111/mpp.13045 ![]() |
[80] |
Kombrink A, Thomma BPHJ (2013) LysM effectors: secreted proteins supporting fungal life. PLoS Pathog 9: e1003769. https://doi.org/10.1371/journal.ppat.1003769 ![]() |
[81] |
Miya A, Albert P, Shinya T, et al. (2007) CERK1, a LysM receptor kinase, is essential for chitin elicitor signaling in Arabidopsis. Proceedings of the National Academy of Sciences 104: 19613-19618. https://doi.org/10.1073/pnas.0705147104 ![]() |
[82] |
Pazzagli L, Seidl-Seiboth V, Barsottini M, et al. (2014) Cerato-platanins: Elicitors and effectors. Plant Sci 228: 79-87. https://doi.org/10.1016/j.plantsci.2014.02.009 ![]() |
[83] |
Wösten HAB (2001) Hydrophobins: multipurpose proteins. Annu Rev Microbiol 55: 625-646. https://doi.org/10.1146/annurev.micro.55.1.625 ![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |