Abstract The endoplasmic reticulum (ER) is the entry portal of the conventional secretory pathway where the newly synthesized polypeptides fold, modify, and assemble. The ER responses to the unfolded proteins in its lumen (ER stress) by triggering intracellular signal transduction pathways include the ER-associated degradation (ERAD) pathway and the unfolded protein response (UPR) pathway. In yeast and mammals, the ubiquitin ligase Hrd1 is indispensable for the ERAD pathway, and also Hrd1-mediated ERAD pathway plays a crucial role in maintaining homeostasis and metabolism of human beings. However, the underlying physiological roles and regulatory mechanism of the Hrd1-involved ERAD pathway in the plant pathogenic fungi are still unclear. Here, we identified the Hrd1 orthologous proteins from 9 different fungi and noticed that these Hrd1 orthologs are conserved. Through identification of MoHrd1 putative interacting proteins by co-immunoprecipitation assays and enrichment analysis, we found that MoHrd1 is involved in the secretory pathway, energy synthesis, and metabolism. Taken together, our results suggest that MoHrd1 is conserved among fungi and play an important role in cellular metabolism and infection-related development. Our finding helps uncover the mechanism of Hrd1-involved ERAD pathway in fungi and sheds a new light to understand the pathogenic mechanism of Magnaporthe oryzae. Keywords: endoplasmic reticulum–associated degradation (ERAD), the ubiquitin ligase Hrd1, mass spectrometry, Magnaporthe oryzae Introduction Approximately one-third of the cellular proteomes is destined for the plasma membrane, extracellular space, or secretory compartments via targeting to the entry portal of the secretory pathway-endoplasmic reticulum (ER).^[45]1 The newly synthesized polypeptides enter into the ER through the Sec61 translocon complex for folding and N-linked glycosylation. Then, a diverse set of chaperones and accessory proteins assist immature polypeptides in folding to the native conformation, whereas chaperone interactions keep these proteins soluble to prevent aggregation in the ER lumen.^[46]2–[47]5 The ER is a checkpoint of the secretory pathway, with hosting a stringent and specialized protein quality control (PQC) system, which delivers the correctly folded proteins to the Golgi complex while retaining the abnormal folding proteins in the ER.^[48]1,[49]6,[50]7 The ER PQC machinery regulates both the ER-associated degradation (ERAD) pathway for eliminating the ultimately unfolded proteins and unfolded protein response (UPR) pathway, a signaling cascade pathway for reducing protein biosynthesis and facilitating ER folding and degradation ability.^[51]7–[52]9 ERAD pathway removed the aberrant proteins via a series of tightly coupled stages consisting of substrate recognition, retrotranslocation, and ubiquitin proteasome system (UPS)-dependent degradation^[53]4,[54]8,[55]10,[56]11 In the budding yeast, ER luminal substrates are recognized by the molecular chaperones such as Kar2, Scj1, Jem1, and Pdi1, retrotranslocated, and ubiquitinated by the Hrd1 complex and are finally delivered to the UPS for degradation.^[57]4,[58]12–[59]16 The multitransmembrane protein Hrd1 is the E3 ubiquitin ligase for the substrate ubiquitination. Moreover, cryo-electron microscopic structure suggested that the retrotranslocation channel for the misfolded proteins crossing the ER membrane is formed by the membrane protein Hrd1.^[60]17–[61]20 In the mammals, Hrd1 is indispensable for the misfolded proteins which retrotranslocated from the ER lumen to the cytosolic UPS for degradation.^[62]6,[63]21–[64]23 Moreover, Hrd1-mediated ERAD pathway is essential for maintaining metabolic homeostasis and provides a breathtaking new challenge of disease control such as Alzheimer disease and Parkinson disease.^[65]24,[66]25 However, the physiological roles of the Hrd1-involved ERAD pathway in the plant pathogenic fungi are still unclear. The hemibiotrophic filamentous fungus Magnaporthe oryzae (Pyricularia oryzae) is the causal agent of rice blast disease which is one of the most destructive diseases of cultivated rice worldwide and threats rice production and global food security seriously. Magnaporthe oryzae has been used as a primary model organism for studying plant-fungi interaction.^[67]26–[68]29 Similar to other fungi and oomycetes, rice blast fungus secretes some protein termed effectors during infection development for suppressing plant immunity and establishing plant disease. Sequentially, M. oryzae has evolved two distinct secretion systems: on one hand, apoplastic effectors such as Bas4 and Slp1 are delivered into the space between the fungal cell wall and extra-invasive hyphal membrane (EIHM) by conventional ER-Golgi secretion pathway; on the other hand, cytoplasmic effectors (such as Pwl2, AvrPiz-t, and AvrPi9) are secreted from invasive hyphae into the extracellular compartment and the biotrophic interfacial complex (BIC).^[69]30–[70]36 The ER-associated secretory pathway is important for rice blast disease development. Previous studies elaborated that MoHac1, MoLhs1, and MoKar2 involved in UPR pathway are essential for asexual development and infection-related morphogenesis,^[71]37,[72]38 but little is known about another cascade pathway of the ER PQC system in M. oryzae. Protein-protein interactions play a vital role in many cellular processes, such as translation and signal transduction.^[73]39–[74]41 Immunoprecipitation (IP) is a proven technique for the enrichment and isolation of interacting proteins via an antigen binding to a specific antibody.^[75]42 Tandem mass spectrometry (MS) could identify the proteins which directly or indirectly bind each other.^[76]39,[77]43 Thus, IP-MS is a feasible and complementary approach to screen interaction proteins and investigate the potential biological function of the target protein. For example, Li et al. screened out MoArk1 interaction proteins MoCapA and MoCapB using IP-MS. Further study showed that MoArk1 regulated MoCapA and MoCapB to involve in control of growth, conidiation, and pathogenesis in M. oryzae.^[78]44 In this study, we identified the Hrd1 orthologous proteins in different fungi and noticed that these proteins are conserved. Moreover, we demonstrated that MoHrd1 could be crucial for the secretory pathway, energy synthesis, and metabolism in M. oryzae both by bioinformatics analysis and by the co-immunoprecipitation (Co-IP) assays. Materials and Methods Strains and culture conditions The M. oryzae Guy11 strain was used as wild type (WT) in this study. All strains were cultured on complete medium (CM) agar plates at 28°C (CM: 10 g D-glucose, 2 g peptone, 1 g yeast extract, 1 g casamino acid, 50 mL 20× nitrate salts, 1 mL trace elements, 1 mL vitamin solution, 15 g agar, and 1 L distilled water). Liquid CM medium was used to prepare the mycelia for DNA and protein extraction. Construction of MoHRD1::GFP To establish MoHRD1::GFP construct, full-length genomic DNA of MoHRD1 with the native promoter (1.5 kb upstream fragment) was cloned from Guy11 using primer pair HRD1 Com-F and HRD1 Com-R which was then co-introduced with Xho I (Takara, Shiga, Japan)-digested pYF11 into the XK125 yeast competent cell. Plasmids of MoHRD1::GFP were recovered from generating Trp^+ yeast transformants and verified by performing polymerase chain reaction (PCR) assays with primer pair HRD1 Con-F and GFP-R. IP assays The MoHRD1-GFP fusion construct was introduced by transformation into protoplasts of the Mohrd1 null mutant. Mycelia were ground into fine powder in liquid nitrogen and resuspended in 20 mL lysis buffer (10 mM Tris-HCl [pH 7.5], 150 mM NaCl, 0.5 mM EDTA, 0.5% NP-40) with freshly added 1 mM phenylmethylsulfonyl fluoride (PMSF) and 10 μL of protease inhibitor cocktail (Sigma Aldrich , Shanghai, China). Then, total proteins were incubated with anti-GFP (green fluorescent protein) agarose beads (ChromoTek, Planegg-Martinsried, Germany). Proteins bound to GFP agarose beads were eluted after a series of washing steps according to the manufacturer’s protocol. Analysis of bound proteins was performed by the Beijing Genomics Institute using MS. Pathway enrichment analysis and Gene Ontology enrichment analysis of MoHrd1 putative interaction proteins KEGG (Kyoto Encyclopedia of Genes and Genome; [79]http://www.genome.jp/kegg/) is a primary biological database of the pathway and an efficient instrument for characterization of metabolism and metabolic network.^[80]45 The most important metabolic pathways and signal transduction pathways that putative interacting proteins are involved could be identified via pathway enrichment analysis. Providing functional information of gene product, GO (Gene Ontology; [81]http://www.geneontology.org/) is a published bioinformatics database that classifies functions along molecular function, biological process, and cellular component in 3 aspects. The enrichment analysis of MoHrd1 putative interaction proteins was performed using ClueGO^[82]46 which is a Cytoscape plugin to improve the biological interpretation of gene lists and a functional organization network was constructed. Results Identification of Hrd1p orthologous proteins in different fungi Using the amino acid sequence of Hrd1p from Saccharomyces cerevisiae as a query, we performed a homology search in 9 fungal species including M. oryzae, Ustilago maydis, Fusarium graminearum, Botrytis cinerea, Neurospora crassa, Aspergillus fumigatus, Aspergillus nidulans, Aspergillus oryzae, and Aspergillus niger in the NCBI (National Center for Biotechnology Information) database ([83]http://blast.ncbi.nlm.nih.gov/). Nine potential candidates including the orthologous protein of M. oryzae encoded by MGG_09205 were found. Based on the sequences of these 9 Hrd1 orthologous proteins, a phylogenetic tree was constructed using the ClustalX and the MEGA version 6 software^[84]47,[85]48 ([86]Figure 1). The phylogenetic dendrogram shows that the 4 potential Hrd1 orthologs in the Aspergillus species are most closely related, whereas the Hrd1 orthologous protein in M. oryzae is closely related to orthologs in N. crassa, F. graminearum, and B. cinerea. Figure 1. [87]Figure 1. [88]Open in a new tab Phylogenetic tree of the Hrd1 orthologous proteins in different fungi. The Hrd1 protein sequences from different fungi were aligned using the ClustalX software. Phylogenetic analysis was performed using MEGA version 6. Hrd1 is highly conserved in different fungi To investigate and compare the function of these Hrd1 orthologous proteins, functional domains were identified using the SMART software^[89]49,[90]50 ([91]http://smart.embl-heidelberg.de/; [92]Figure 2). Identification of domains from protein sequences demonstrated that Hrd1 proteins are highly conserved among selected fungi and all contain a RING domain and several transmembrane regions. Moreover, 4 to 6 transmembrane regions and a RING domain localize to the similar loci from N-terminus, whereas most of them hold a coiled-coil domain between the last 2 transmembrane regions. These observations suggest that the function of these Hrd1 orthologous proteins could be highly conserved. Figure 2. [93]Figure 2. [94]Open in a new tab Domain architecture of Hrd1 orthologous proteins in different fungi. The protein functional domains were identified using the SMART software ([95]http://smart.embl-heidelberg.de/). We further performed a sequence alignment using the ClustalX software.^[96]47 Amino acid alignments of Hrd1 orthologs display high conservation. Almost one-half of the amino acids are highly conserved among these 10 fungi, especially in the N-terminus transmembrane regions and RING region ([97]Figure 3). In the RING region, almost one-third of the amino acids are identical among all selected fungi; meanwhile, more than one-half of the amino acids are uniform among F. graminearum, B. cinerea, N. crassa, M. oryzae, A. fumigatus, A. nidulans, A. oryzae, and A. niger ([98]Figure 4A). For visually displaying the conservation of the RING domain, we built a sequence logo diagram using the WebLogo software ([99]http://weblogo.berkeley.edu/logo.cgi).^[100]51 Simultaneously, the amino acids logo diagram showed that Hrd1 RING domain is highly conserved in these selected fungi ([101]Figure 4B). Figure 3. [102]Figure 3. [103]Open in a new tab Hrd1 is conserved in different fungi. Sequence alignment of Hrd1 in different fungi. Sequence identity is shaded in black; shading in dark gray shows that amino acid sequence are conserved over 75% species, whereas shading in light gray shows only over 50% species. The alignment was performed using ClustalX software. Figure 4. [104]Figure 4. [105]Open in a new tab Hrd1 RING domain is conserved in different fungi. (A) Sequence alignment of the Hrd1 RING domains in different fungi. Red shows sequence identity, blue shows sequence is conserved over 75% species, whereas purple displays only over a half, and the consensus exhibited the conserved amino acids. (B) Sequence logo diagram shows Hrd1 RING domain ortholog in different fungi. The alignment was performed using the ClustalX software. The logo was built using the WebLogo software ([106]http://weblogo.berkeley.edu/logo.cgi). In the human pathogenic fungi A. fumigatus, loss of HrdA (the Hrd1 ortholog in A. fumigatus) exhibited defects in degradation of a folding-defective ERAD substrate, as well as activation of the UPR.^[107]52 In the filamentous fungus A. oryzae, HrdA (the Hrd1 ortholog in A. oryzae) mediated ERAD pathway is required for the degradation of the moderate-level MsdS mutant.^[108]53 We found that approximately 80% of the MoHrd1 RING domain amino acid sequences are identical to those in A. fumigatus and A. oryzae ([109]Figure 4A). Meanwhile, Hrd1 orthologs in M. oryzae, A. fumigatus, and A. oryzae contain similar functional domain combination ([110]Figure 2). In brief, these results suggest that MoHrd1 may play an important role in the ER PQC system of M oryzae. Putative interaction proteins of MoHrd1 were identified by IP-MS To identify MoHrd1 putative interaction proteins, Co-IP assays were performed. A GFP-encoding gene was fused at the C-terminus of the MoHrd1-coding sequence using the native promoter and was transformed into the ΔMohrd1 mutant (unpublished data). Wild type expressing GFP gene driven by a constitutive promoter was used as a negative control. Western blots were used to detect Co-IP assays with GFP-tag proteins extracted from MoHrd1::GFP and the GFP control strain of the WT. As shown in [111]Figure 5, the anti-GFP antibody detected the 117-kDa MoHrd1::GFP-fusing protein and the GFP, 26-kDa bands, respectively, from total proteins isolated from the corresponding strains and proteins eluted from anti-GFP beads, indicating that GFP-tag proteins were successfully co-immunoprecipitated in the elution buffer. Figure 5. Figure 5. [112]Open in a new tab Western blots of total proteins (total) and proteins eluted (elution) from the anti-GFP agarose beads from the wild type constitutively expressing GFP and transformant expressing MoHrd1::GFP were probed with the anti-GFP (GFP-Tag 7G9 Mouse mAb; Abmart, Shanghai, China) antibody. GFP indicates green fluorescent protein. Meanwhile, tandem MS was employed to identify MoHrd1 putative interacting proteins in the immunoprecipitates. The candidate interaction proteins present in the different replicates but not present in the negative control were regarded as potential interaction proteins; in total, 161 putative interaction proteins were identified ([113]Supplementary Table S2). In S. cerevisiae, 60 Hrd1p interaction proteins were identified using matrix-assisted laser desorption/ionization-time of flight MS and liquid chromatography-tandem MS.^[114]54 Using S cerevisiae Hrd1 interaction proteins as a query, we further performed a homologous search in M. oryzae; in total, 60 homologous proteins were found ([115]Supplementary Table S3). Moreover, 5 of the MoHrd1 putative interaction proteins consisting of MGG_05193, MGG_13508, MGG_01790, MGG_04400, and MGG_02297 were present in these 60 homologous proteins. MGG_13508, MGG_02297, and MGG_05193 were the homologous proteins of S. cerevisiae Hrd3, Yos9, and Cdc48, respectively. In yeast, Hrd3 was required for ERAD pathway and was indispensable for Hrd1 stability^[116]55; Yos9 was essential for ERAD substrate recognition and took part in delivering substrates into Hrd1 complex^[117]56–[118]58; Cdc48 formed a heterotrimeric complex with Npl4 and Ufd1 to facilitate substrate retrotranslocation.^[119]59 In short, S. cerevisiae Hrd3, Yos9, and Cdc48 were required for Hrd1-involved ERAD pathway. Taken together, this may suggest that MoHrd1 synergizes with MGG_13508, MGG_02297, and MGG_05193 to play a similar role in the ERAD pathway of M. oryzae. Pathway enrichment analysis and GO enrichment analysis of the MoHrd1 putative interaction proteins Pathway enrichment analysis (KEGG) and functional enrichment analysis (GO) were applied to elucidating the biological functions of the putative interaction proteins related to MoHrd1. Enriched results were subjected to multiple testing adjustment with a threshold value (adjusted P value with Bonferroni correction) less than 0.05. In total, 5 significantly enriched pathways were annotated, and their descriptions suggest that they are highly correlated with energy metabolisms such as tricarboxylic acid (TCA) cycle, ribosome biogenesis, and secretory pathways such as protein export and protein processing in ER in eukaryotes ([120]Figure 6A). To better exhibit functional consequence, only the top 10 significant enriched GO terms are shown in [121]Figure 6B. Consistent with pathway enrichment results, several GO terms such as gene expression, translation, ribosome biogenesis, metabolic process, and biosynthetic process indicate that these protein data sets belong to the secretory pathway, energy synthesis, and metabolism. Figure 6. [122]Figure 6. [123]Open in a new tab Pathway enrichment analysis and GO (Gene Ontology) enrichment analysis of MoHrd1 putative interaction proteins. (A) Advanced bubble chart shows enrichment of MoHrd1 putative interaction proteins in different signaling pathways. Y-axis label represents pathway, and X-axis label represents rich factor (rich factor = amount of putative interaction proteins enriched in the pathway/amount of all proteins in background proteins set). Size and color of the bubble represent the amount of putative interaction proteins enriched in pathway and enrichment significance, respectively. (B) Column diagram shows enrichment of MoHrd1 putative interaction proteins in different metabolic pathways. Y-axis label represents pathway, and X-axis label represents the ratio of putative interaction proteins. Discussion In this study, we identified the Hrd1 orthologous proteins in 9 different fungi and found that these orthologs are conserved. Further enrichment analysis suggested that the Hrd1 ortholog in M. oryzae should be essential for the secretory pathway, energy synthesis, and metabolism. Previous studies on the ubiquitin ligase Hrd1 revealed that the multispanning membrane Hrd1 is required for the ER luminal substrates retrotranslocation and ubiquitination, and its autoubiquitination triggered protein retrotranslocation. Moreover, this membrane-anchored protein is involved in the constitution of the retrotranslocation channel for the movement of the abnormal polypeptides through the ER membrane.^[124]4,[125]17,[126]18,[127]20,[128]60 And, the Hrd1 RING domain is indispensable for the function of Hrd1 in the ERAD pathway; the site-specific mutant (C399S) is not functional.^[129]61 More than that recent studies revealed that the Hrd1 RING domain plays an important role in the activation of its cognate E2 ubiquitin-conjugating enzyme.^[130]62 Our results showed that all 10 Hrd1 orthologs contain similar functional domain combination ([131]Figure 2) and the RING domain is highly conserved ([132]Figure 4). This may suggest that the Hrd1 orthologs including MoHrd1 play a similar role in the ERAD pathway among different fungi. Recent studies have shown that M. oryzae secreted effectors to suppress plant immunity and support rice blast development.^[133]30,[134]33 For instance, M. oryzae secreted Slp1 (secreted LysM protein 1) to overcome PTI (pathogen-associated molecular patterns–triggered immunity) via suppressing chitin-induced plant immune responses.^[135]63 Hence, correct secretion of effectors is necessary for the M oryzae infection–related development and abnormal secretory proteins should be eliminated in time. Consequently, MoHrd1-mediated ERAD pathway may be important for the infection-related morphogenesis. Similarly, our results showed that MoHrd1 is related to the conventional secretory pathway (translation, protein processing in ER, protein export; [136]Figure 6). This may suggest that MoHrd1 is required for the pathogenesis of M. oryzae via mediating secretion of effectors. Misfolded proteins undergo 2 sequential steps after retrotranslocation and ubiquitination in the cytosol, N-linked glycans are removed from misfolded proteins, and remaining proteins delivered to the proteasome for degradation into small peptides; meanwhile, free oligosaccharides are generated with the Png1 catalyzation.^[137]4 Our results showed that MoHrd1 putative interaction proteins mainly enriched the energy synthesis and metabolism ([138]Figure 6). This may reveal that the product of the ERAD pathway such as small peptides and free oligosaccharides is recycled into the cellular metabolism such as the TCA cycle in M. oryzae. Hrd1 (HMG-CoA reductase degradation) was discovered in the genetic analysis of Hmg2 degradation. Hmg2 is one of the HMG-CoA (3-hydroxy-3-methyl-glutaryl-CoA) reductase isozymes which is considered as a rate-controlling enzyme of sterol synthesis (mevalonate pathway) in eukaryotes. Hrd1 contributes to the regulation of sterol synthesis via degradation of Hmg2.^[139]64,[140]65 Similar to these studies, KEGG enrichment analysis and GO enrichment analysis results supported that Hrd1 contributes to the energy synthesis and metabolism ([141]Figure 6). In summary, we found that MoHrd1 is conserved with other fungi and may play an important role in cellular metabolism and infection-related development. This study helps unveil the physiological roles and regulatory mechanism of Hrd1-involved ERAD pathway in fungi and provides a new idea to understand the pathogenic mechanism of M. oryzae. Supplementary Material Supplementary material [142]Table_S1.pdf^ (160.9KB, pdf) Supplementary material [143]Table_S2.xlsx^ (68.6KB, xlsx) Supplementary material [144]Table_S3.xlsx^ (22.8KB, xlsx) Acknowledgments