Portal de Periódicos da CAPES

microRNA-4331 Promotes Transmissible Gastroenteritis Virus (TGEV)-induced Mitochondrial Damage Via Targeting RB1, Upregulating Interleukin-1 Receptor Accessory Protein (IL1RAP), and Activating p38 MAPK Pathway In Vitro

2017; Elsevier BV; Volume: 17; Issue: 2 Linguagem: Inglês

10.1074/mcp.ra117.000432

1535-9484

Xiaomin Zhao, Xiaoyuan Bai, Lijuan Guan, Juejun Li, Xiangjun Song, Xuelian Ma, Jianxiong Guo, Zhichao Zhang, Qian Du, Yong Huang, Dewen Tong,

Virus-based gene therapy research

Transmissible gastroenteritis virus (TGEV), a member of the coronaviridae family, could cause fatal diarrhea of piglets and result in numerous economic losses. Previous studies demonstrated that TGEV infection could lead to mitochondrial damage and upregulate miR-4331 level. So miR-4331 may play an important regulatory role in the control of mitochondrial function. To explore the potential role of miR-4331 in mitochondrial damage, we adopted a strategy consisting of quantitative proteomic analysis of porcine kidney (PK-15) cells in response to miR-4331 and TGEV infection. Eventually, 69 differentially expressed proteins were gained. The target of miR-4331 was identified. The effects of miR-4331 and its target RB1 on mitochondrial Ca2+ level, mitochondrial membrane potential (MMP), interleukin-1 receptor accessory protein (IL1RAP), p38 MAPK signaling pathway were investigated. The results showed that miR-4331 elevated mitochondrial Ca2+ level, reduced MMP, targets Retinoblastoma 1 (RB1), upregulated IL1RAP, and induced activation of p38 MAPK pathway during TGEV infection. RB1 was identified as the direct targets of miR-4331 and downregulated IL1RAP, suppressed the activation of p38 MPAK, and attenuated TGEV-induced mitochondrial damage. In addition, IL1RAP played a positive role in activating p38 MAPK signaling and negative role in TGEV-induced mitochondrial damage. The data indicate that miR-4331 aggravates TGEV-induced mitochondrial damage by repressing expression of RB1, promoting IL1RAP, and activating p38 MAPK pathway. Transmissible gastroenteritis virus (TGEV), a member of the coronaviridae family, could cause fatal diarrhea of piglets and result in numerous economic losses. Previous studies demonstrated that TGEV infection could lead to mitochondrial damage and upregulate miR-4331 level. So miR-4331 may play an important regulatory role in the control of mitochondrial function. To explore the potential role of miR-4331 in mitochondrial damage, we adopted a strategy consisting of quantitative proteomic analysis of porcine kidney (PK-15) cells in response to miR-4331 and TGEV infection. Eventually, 69 differentially expressed proteins were gained. The target of miR-4331 was identified. The effects of miR-4331 and its target RB1 on mitochondrial Ca2+ level, mitochondrial membrane potential (MMP), interleukin-1 receptor accessory protein (IL1RAP), p38 MAPK signaling pathway were investigated. The results showed that miR-4331 elevated mitochondrial Ca2+ level, reduced MMP, targets Retinoblastoma 1 (RB1), upregulated IL1RAP, and induced activation of p38 MAPK pathway during TGEV infection. RB1 was identified as the direct targets of miR-4331 and downregulated IL1RAP, suppressed the activation of p38 MPAK, and attenuated TGEV-induced mitochondrial damage. In addition, IL1RAP played a positive role in activating p38 MAPK signaling and negative role in TGEV-induced mitochondrial damage. The data indicate that miR-4331 aggravates TGEV-induced mitochondrial damage by repressing expression of RB1, promoting IL1RAP, and activating p38 MAPK pathway. Transmissible gastroenteritis virus (TGEV) 1The abbreviations used are: TGEV, transmissible gastroenteritis virus; PK-15 cells, porcine kidney (PK) 15 cells; miRNA, microRNA; MMP, mitochondrial membrane potential; UTR, untranslated region; GO, gene ontology; KEGG, kyoto encyclopedia of genes and genomes; MOI, multiplicity of infection; TMT, tandem mass tags; hpi, hours post-infection; hpt, hours post-transfection; HPLC, high performance liquid chromatography; MS, mass spectrometry. 1The abbreviations used are: TGEV, transmissible gastroenteritis virus; PK-15 cells, porcine kidney (PK) 15 cells; miRNA, microRNA; MMP, mitochondrial membrane potential; UTR, untranslated region; GO, gene ontology; KEGG, kyoto encyclopedia of genes and genomes; MOI, multiplicity of infection; TMT, tandem mass tags; hpi, hours post-infection; hpt, hours post-transfection; HPLC, high performance liquid chromatography; MS, mass spectrometry.is an enveloped enteropathogenic coronavirus with positive-sense single-stranded RNA genome (1.Weiss S.R. Leibowitz J.L. Coronavirus Pathogenesis. 2011; 81: 85-164Google Scholar). TGEV infects pigs and especially causes piglets up to 14 days of age high mortality, which can reach to 100% (2.Hu W.W. Zhu L.Q. Yang X. Lin J. Yang Q. The epidermal growth factor receptor regulates cofilin activity and promotes transmissible gastroenteritis virus entry into intestinal epithelial cells.Oncotarget. 2016; 7: 12206-12221Crossref PubMed Scopus (38) Google Scholar, 3.Wang X. Hu W. Zhu L. Yang Q. Bacillus subtilis and surfactin inhibit the transmissible gastroenteritis virus from entering the intestinal epithelial cells.Biosci. Rep. 2017; 37Crossref Scopus (34) Google Scholar) and result in numerous economic losses. Our previous work demonstrated that TGEV infection can reduce MMP, leading to mitochondrial damage (4.Ding L. Xu X. Huang Y. Li Z. Zhang K. Chen G. Yu G. Wang Z. Li W. Tong D. Transmissible gastroenteritis virus infection induces apoptosis through FasL- and mitochondria-mediated pathways.Vet. Microbiol. 2012; 158: 12-22Crossref PubMed Scopus (50) Google Scholar, 5.Ding L. Zhao X. Huang Y. Du Q. Dong F. Zhang H. Song X. Zhang W. Tong D. Regulation of ROS in transmissible gastroenteritis virus-activated apoptotic signaling.Biochem. Biophys. Res. Commun. 2013; 442: 33-37Crossref PubMed Scopus (27) Google Scholar). Mitochondrion, an organelle of eukaryocyte, not only acts as the energy metabolism factory, but also is involved in many key biological processes, such as apoptosis, pathogenic infection (6.Torraco A. Peralta S. Iommarini L. Diaz F. Mitochondrial Diseases Part I: mouse models of OXPHOS deficiencies caused by defects in respiratory complex subunits or assembly factors.Mitochondrion. 2015; 21: 76-91Crossref PubMed Scopus (34) Google Scholar, 7.Su W.C. Lin Y.F. Yu X.P. Wang Y.X. Lin X.D. Su Q.Z. Shen D.Y. Chen Q.X. Mitochondria-associated apoptosis in human melanoma cells induced by cardanol monoene from cashew nut shell liquid.J. Agric. Food Chem. 2017; 65: 5620-5631Crossref PubMed Scopus (16) Google Scholar). Mitochondrion is a very sensitive organelle to microenvironmental changes in cells and much easier becomes dysfunction than other organelles (8.Gomez-Cabrera M.C. Sanchis-Gomar F. Garcia-Valles R. Pareja-Galeano H. Gambini J. Borras C. Vina J. Mitochondria as sources and targets of damage in cellular aging.Clin. Chem. Lab. Med. 2012; 50: 1287-1295Crossref PubMed Scopus (56) Google Scholar). Mitochondria, the major source of energy, in the form of ATP, are closely interconnected with cell death. Ca2+ is essential for maintain mitochondrial function. However, Ca2+ accumulation in mitochondria can impair mitochondrial function, resulting in a transient depolarization of MMP and reducing ATP production (9.Duchen M.R. Mitochondria and calcium: from cell signalling to cell death.J. Physiol. 2000; 529: 57-68Crossref PubMed Scopus (932) Google Scholar, 10.Santulli G. Xie W. Reiken S.R. Marks A.R. Mitochondrial calcium overload is a key determinant in heart failure.Proc. Natl. Acad. Sci. U.S.A. 2015; 112: 11389-11394Crossref PubMed Scopus (320) Google Scholar, 11.Peng T.I. Jou M.J. Oxidative stress caused by mitochondrial calcium overload.Ann. N.Y. Acad. Sci. 2010; 1201: 183-188Crossref PubMed Scopus (352) Google Scholar). Both Ca2+ accumulation and MMP depolarization can cause mitochondrial damage, so the alterations of Ca2+ level and MMP are used to assess mitochondrial function (12.Perez M.J. Quintanilla R.A. Development or disease: duality of the mitochondrial permeability transition pore.Dev. Biol. 2017; 426: 1-7Crossref PubMed Scopus (88) Google Scholar). miRNAs are small noncoding RNA species containing about 22 nucleotides and contribute to regulating many cellular processes including apoptosis, development, differentiation, and cell cycle (13.Wilson R.C. Doudna J.A. Molecular mechanisms of RNA interference.Annu. Rev. Biophys. 2013; 42: 217-239Crossref PubMed Scopus (681) Google Scholar). We previously reported that miR-4331 level was upregulated during TGEV-induced mitochondrial dysfunction (14.Song X. Zhao X. Huang Y. Xiang H. Zhang W. Tong D. Transmissible gastroenteritis virus (TGEV) infection alters the expression of cellular microRNA species that affect transcription of TGEV gene 7.Int. J. Biol. Sci. 2015; 11: 913-922Crossref PubMed Scopus (19) Google Scholar). Whereas, whether miR-4331 participate in regulating TGEV-induced mitochondrial damage in PK-15 cells is unclear. According to GO enrichment and KEGG pathway analysis of miR-4331 targets, miR-4331 targets can be enriched in mitochondria and mitochondria-related pathways. Therefore, we postulate miR-4331 can participate in regulating mitochondrial dysfunction in PK-15 cells. To investigate the regulatory effects of miR-4331 on TGEV-induced mitochondrial damage, miR-4331 mimics or mimics control was transfected into PK-15 cells and the proteomic analysis was performed using LC-MS/MS labeled by Tandem Mass Tags (TMT). Then we identified the targets of miR-4331 and detected the effects of miR-4331 target on mitochondrial damage and p38 MAPK pathway. We demonstrated that miR-4331 increased mitochondrial Ca2+ level and decreased MMP via activating p38 MAPK pathway and upregulating expression of IL1RAP. Moreover, miR-4331 activated p38 MAPK pathway via targeting RB1 and upregulating IL1RAP. Our data suggest that miR-4331 could facilitate TGEV-induced mitochondrial damage in PK-15 cells via targeting RB1, promoting IL1RAP, and activating p38 pathway. Anti-RB1 was purchased from Santa Cruz Biotechnology. Horseradish peroxidase (HRP)-conjugated secondary antibody was purchased from Pierce (US). Anti-p-p38 and anti-β-actin antibody were purchased from Cell Signaling Technology (US). Anti-IL1RAP and anti-GAPDH primary antibody was purchased from GeneTex (US). The psiCHECK-2 plasmid and dual-luciferase reporter assay system were purchased from Promega (US). The primers and siRNAs were synthesized by Ribo Biotech (RiboBio, China). The TGEV Shaanxi strain was isolated from TGEV-infected piglets (15.Ding L. Chen G.D. Xu X.G. Tong D.W. Isolation and identification of porcine transmissible gastroenteritis virus Shaanxi strain and sequence analysis of its N gene.Chin. J. Vet. Med. 2011; 47 (09–11)Google Scholar). PK-15 cells were obtained from ATCC (CCL-33) and grown in Dulbecco's Minimal Essential Medium (DMEM) supplemented with 10% fetal bovine serum (Hyclone, US), 100 IU of penicillin, and 100 mg/ml streptomycin, at 37 °C in 5% CO2 atmosphere incubator. In order to investigating the effect of miR-4331 on mitochondrial damage in response to TGEV infection, PK-15 cells were transfected with 100 nm miR-4331 mimics or negative control using Lipofectamine 3000 and infected with TGEV at an MOI of 1.0 at 24 h post-transfection (hpt). The cells were collected for quantitative proteomic analysis at 24 h post-infection (hpi). Two biological replicates preparations labeled with TMT were analyzed. To evaluate the transfection efficiency, the miR-4331 level was measured by real-time PCR. Reverse transcription reactions and real-time PCR were performed as described previously (14.Song X. Zhao X. Huang Y. Xiang H. Zhang W. Tong D. Transmissible gastroenteritis virus (TGEV) infection alters the expression of cellular microRNA species that affect transcription of TGEV gene 7.Int. J. Biol. Sci. 2015; 11: 913-922Crossref PubMed Scopus (19) Google Scholar). The relative quantification of miRNAs was normalized to U6 using the two-ddCt method (16.Livak K.J. Schmittgen T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method.Methods. 2001; 25: 402-408Crossref PubMed Scopus (124148) Google Scholar). The PK-15 cells were sonicated on ice in lysis buffer (8 m urea, 1% Triton-100, 65 mm DTT and 0.1% protease inhibitor mixture) and centrifuged at 12,000 × g for 20 min at 4 °C. The clarified supernatant was collected. Finally, the protein was precipitated with cold 15% TCA for 2 h at −20 °C and centrifuged at 4 °C for 10 min. The precipitate was redissolved in buffer (8 m urea, 100 mm TEAB, pH 8.0). One hundred micrograms protein of each sample was digested with trypsin. The two control samples and two treatment samples were respectively labeled with TMT (126, 129, 127, and 130). The samples were fractionated by high pH reverse-phase High Performance Liquid Chromatography (HPLC). Briefly, peptides were dried by vacuum centrifugation. Then the Peptides were dissolved in 0.1% formic acid and directly loaded onto a reversed-phase precolumn (Acclaim PepMap 100, Thermo Scientific). Peptides were separated using a reversed-phase analytical column (Acclaim PepMap RSLC, Thermo Scientific) and analyzed by Q ExactiveTM hybrid quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific). The peptides were subjected to NSI source followed by tandem mass spectrometry (MS/MS) in Q ExactiveTM (Thermo Fisher Scientific) coupled online to the HPLC. Intact peptides were detected in the Orbitrap at a resolution of 70,000. Peptides were selected for MS/MS using NCE setting as 32. Ion fragments were detected in the Orbitrap at a resolution of 17,500. A data-dependent procedure that alternated between one MS scan followed by 20 MS/MS scans was applied for the top 20 precursor ions above a threshold ion count of 2E4 in the MS survey scan with 30.0 s dynamic exclusion. The electrospray voltage was 2.0 kV. Automatic gain control (AGC) was used to prevent overfilling of the ion trap. 5E4 ions were accumulated for generation of MS/MS spectra. For MS scans, the m/z scan range was 350 to 1800. Fixed first mass was set as 100 m/z. MS/MS data were analyzed by Mascot (version 2.3.0) according to a previously described protocol (17.Cox J. Matic I. Hilger M. Nagaraj N. Selbach M. Olsen J.V. Mann M. A practical guide to the MaxQuant computational platform for SILAC-based quantitative proteomics.Nat. Protoc. 2009; 4: 698-705Crossref PubMed Scopus (624) Google Scholar). Briefly, raw MS data files were processed using the LC/MS software Proteome Discoverer (version 1.3.0.339) (Thermo scientific, US) and converted into the Mascot generic format (mgf) files. Mascot software performed peak generation, precursor mass recalibration, extraction of TMT 6-plex reporter ions intensity, and calculation of TMT 6plex reporter ions intensity ratios. For each MS/MS spectra, top 10 most intense peaks in every 100 Da window were extracted for database search. Then Mascot search was performed by searching tandem mass spectra against Uniprot Sus scrofa database (UniProt release 2015_04, 26054 entries) (http://www.uniprot.org) concatenated with reversed decoy database and protein sequences of common contaminants. Cleavage enzyme was specified as trypsin/P, maximum number of missing cleavages was set to 2. Mass tolerance for precursor ion was set to 10 ppm and mass tolerance for MS/MS ions was set to 0.02 Da. Carbamyiodomethylation on Cys and TMT6plex (N-term) and TMT6plex (K) were specified as fixed modification and oxidation on Met, acetylation on protein N-terminal were specified as variable modifications. False discovery rate (FDR) thresholds for protein, peptide and modification site were specified at 1%. Mascot software assembled the peptide/protein groups, calculated false discovery rates, filter the identifications. Protein quantification was set to use only unique peptides bearing any modification. The median ratio of TMT 6-plex reporter intensity of unique peptides was set as protein relative abundance changes. In this study, we prepared two duplicate samples. For each protein, we set the average ratio calculated from two duplicate samples as the final quantitation of protein. Student's test was explored to calculate differential significance degree of protein relative abundance changes. t test p value less than 0.05 was considered significantly differentially. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD008174 (http://www.ebi.ac.uk/pride/archive/) (18.Vizcaino J.A. Csordas A. Del-Toro N. Dianes J.A. Griss J. Lavidas I. Mayer G. Perez-Riverol Y. Reisinger F. Ternent T. Xu Q.W. Wang R. Hermjakob H. 2016 update of the PRIDE database and its related tools.Nucleic Acids Res. 2016; 44: 11033Crossref PubMed Scopus (21) Google Scholar). Gene Ontology (GO) annotation of proteome was derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA). Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate proteins pathways. The Wolf Psort database (https://wolfpsort.hgc.jp/) was used to predict subcellular localization. The differentially expressed protein name identifiers (Uniprot accession) were searched against the STRING database (version 10.0) for protein-protein interactions (19.Szklarczyk D. Franceschini A. Kuhn M. Simonovic M. Roth A. Minguez P. Doerks T. Stark M. Muller J. Bork P. Jensen L.J. von Mering C. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored.Nucleic Acids Res. 2011; 39: D561-D568Crossref PubMed Scopus (2716) Google Scholar). To overexpress RB1 and IL1RAP, the full-length sequences of RB1 and IL1RAP were obtained by PCR from PK-15 cells and cloned into plasmid pCI-neo (Promega, USA). The constructions were respectively named pCI-neo-RB1 and pCI-neo-IL1RAP. The primer sequences are shown in supplemental Information S1. PK-15 cells were treated with RIPA lysis buffer containing phenylmethyl sulfonylfluoride (PMSF). The proteins were separated using 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred onto the polyvinylidene difluoride (PVDF) membrane. The membrane was blocked with 5% BSA for 2 h at room temperature and then incubated with primary antibodies overnight at 4 °C. The HRP-conjugated secondary antibodies were used to incubate the membrane for 1 h at room temperature. Total RNA was extracted with TRIzol (Invitrogen, US). 2 μg of total RNA was treated with DNase I (Fermentas, Germany) for 30 min at 37 °C and reversely transcribed using the First-strand cDNA synthesis kit (Invitrogen). The cDNA was used as the template for real-time PCR. The real-time PCR was performed on Bio-Rad iQ5 Real-Time PCR System (Bio-Rad). The reverse transcription primers and real-time PCR primers are shown in supplemental Information S1. 3′ UTRs of candidate targets containing the miR-4331 binding sites were respectively amplified by PCR and were cloned into the vector psiCHECK-2 (Promega). The binding sites of miR-4331 in RB1-wt 3′ UTR were mutated following the mutagenesis protocol (20.Heckman K.L. Pease L.R. Gene splicing and mutagenesis by PCR-driven overlap extension.Nat. Protoc. 2007; 2: 924-932Crossref PubMed Scopus (857) Google Scholar) to generate RB1-mut. The PCR primers are shown in supplemental Information S1. PK-15 cells were co-transfected with plasmid RB1-wt (or RB1-mut) and miR-4331 mimics (or miR-4331 inhibitors) using Lipofectamine 3000 (Invitrogen, US). The luciferase activities were detected at 48 hpt using a Dual-Glo Luciferase Assay System (Promega). The siRNAs of RB1 and IL1RAP were synthesized by GenePharma (GenePharma, China) (supplemental Information S1). PK-15 cells were transfected with 100 nm siRNA or irrelevant siRNA using Lipofectamine 3000 (Invitrogen, US). MMP of PK-15 cells was detected using JC-1 as an indicator (21.Joiner M.L.A. Koval O.M. Li J.D. He B.J. Allamargot C. Gao Z. Luczak E.D. Hall D.D. Fink B.D. Chen B.Y. Yang J.Y. Moore S.A. Scholz T.D. Strack S. Mohler P.J. Sivitz W.I. Song L.S. Anderson M.E. CaMKII determines mitochondrial stress responses in heart.Nature. 2012; 491: 269-273Crossref PubMed Scopus (303) Google Scholar). PK-15 cells were transfected with miR-4331 mimics or miR-4331 inhibitors and were subsequently infected with TGEV at 1 MOI. PK-15 cells were treated with PBS containing JC-1 and incubated at 37 °C for 20 min. The absorbance was measured at 550 ex/590 em. The mitochondrial calcium level was detected using Rhod-2 kit (GENMED, China) following the manufacturer's protocol. The absorbance was measured at 490 ex/590 em. The relative fluorescence unit (RFU) was calculated. PK-15 cells were infected with TGEV at 1 MOI for 24 h. The mitochondrial Ca2+ level and MMP of PK-15 cells were evaluated. The results showed that TGEV infection increased mitochondrial Ca2+ level (Fig. 1A) and decreased mitochondrial MMP (Fig. 1B). miR-4331 targets were predicted using TargetScan and miRanda, by which 618 targets were obtained (supplemental Information S2). The 618 targets of miR4331 were searched against Gene Ontology (GO) database to provide enrichment information on biological processes, molecular functions, and cellular components. Among total 618 predicted targets, 26 targets were related to mitochondria, 4.2% of total targets. GO enrichment of the 618 targets of miR-4331 showed that 89.6% of targets were enriched in metabolic process and 31.6% of targets were linked to mitochondria in cellular component (Fig. 2). We previously reported that TGEV infection caused the reduction of MMP through inducing ROS accumulation and increased miR-4331 level (5.Ding L. Zhao X. Huang Y. Du Q. Dong F. Zhang H. Song X. Zhang W. Tong D. Regulation of ROS in transmissible gastroenteritis virus-activated apoptotic signaling.Biochem. Biophys. Res. Commun. 2013; 442: 33-37Crossref PubMed Scopus (27) Google Scholar, 14.Song X. Zhao X. Huang Y. Xiang H. Zhang W. Tong D. Transmissible gastroenteritis virus (TGEV) infection alters the expression of cellular microRNA species that affect transcription of TGEV gene 7.Int. J. Biol. Sci. 2015; 11: 913-922Crossref PubMed Scopus (19) Google Scholar). To investigate the effects of miR-4331 on mitochondria damage during TGEV infection, the PK-15 cells were transfected with miR-4331 mimics, or miRNA mimics control, miR-4331 inhibitors, inhibitors, and subsequently infected with TGEV at 1 MOI for 24 h. The miR-4331 level remarkably increased by miR-4331 mimics (Fig. 3A) and suppressed by miR-4331 inhibitors (Fig. 3B). miR-4331 mimics led to a decrease of mitochondrial Ca2+ level (Fig. 3C) and an increase of MMP (Fig. 3E). miR-4331 inhibitors reduced mitochondrial Ca2+ level (Fig. 3D) and increased MMP (Fig. 3F). A total of 4209 proteins were identified by high-resolution LC-MS/MS analysis, among which 4165 proteins were quantified (supplemental Information S3). When setting ≥1.2-fold as the upregulated threshold and ≤0.83-fold as the downregulated threshold, 69 differentially expressed proteins were obtained, including 33 upregulated proteins and 36 downregulated proteins (Table I). Raw data of MS are available via ProteomeXchange with identifier PXD008174.Table IDifferentially expressed proteins of PK-15 cells induced by miR-4331 during TGEV infectionaMore information is available in supplemental Information S4.No.ProteinFold changep valueNo.ProteinFold changep value1GBP11.6362.86082E-0836LEMD30.8270.0338922332SDC41.60750.01093548537TPP10.8270.0081982123Uncharacterized protein1.46650.00341164938PHLDB10.8210.0463236424MRPL111.44650.02336227639L3MBTL30.81850.0017537695HK21.4171.38772E-0740CDH60.8180.0045080346LOC1005228871.4110.00401410741CAV10.8188.42356E-067LAMA31.3820.01463565942MAGED20.8170.0156306598HELLS1.3220.00953039743PTPN120.8170.002782659LOC1005124201.31850.02947360344PDLIM50.8170.00968389410LOC1005162931.31750.01102695645DDX540.8130.02263713411SLC7A21.3120.00594739346SNX50.8120.00179892112HMGCS11.3110.00056065747AHNAK0.8052.69598E-1013P4HA11.3030.00203645148G3BP10.7990.00353594614NDUFA131.2940.01811083749LOC1005225540.7950.00142585815IL1RAP1.280.00018053150LOC1006213790.790.00376904816TACSTD21.2750.00401177251CPNE30.7860.02939187317E-cadherin1.2720.00432070252CAMK2D0.7810.0183876118PLOD11.2720.02220411853ERAP10.7770.00016180319ITM2B1.26650.01187829254DAB20.7710.01701997220COX5A1.2620.01795237255MX10.7677.97178E-0921HSPG21.2560.00191595256SNX20.7630.00548833522F11R1.24650.01873858657PSAT10.7550.03840071523TOP11.2436.08874E-0858LOC1005229260.7550.0076631224NDUFB51.2380.04815895959BRI3BP0.74750.04154183325ANXA11.2270.03394556960PCK20.7390.0296672126NDUFS11.2243.01379E-0661SRP720.72250.00033077827MRPL471.2220.00682747762BCAT10.72050.0029141428RELA1.2190.00250815363Jup0.6960.01561187929CD471.21850.04218522564LOC1005241440.61950.00026073230LOC1005137561.2140.0294067965RAPH10.6140.01830628631MRPL241.2140.02314237166ADFP0.580.003031532CYP51A11.2130.00227282367CWF19L10.5610.00189293533FN11.2042.01222E-0668DSG20.5040.00025213134STAT5a0.83150.00520743169PGCP0.3880.04549130435TAP20.830.005291123a More information is available in supplemental Information S4. Open table in a new tab The 69 differentially expressed proteins were searched against Wolf Psort database (https://wolfpsort.hgc.jp/) for prediction of subcellular localization. The subcellular localization analysis revealed that the distribution of differentially expressed proteins was distributed in cytosol (43.48%), nuclear (15.49%), mitochondria (13.04%), extracellular (11.59%), plasma membrane (7.25%), cytosol and nuclear (2.9%), endoplasmic reticulum (2.9%), and peroxisome (2.9%) (Fig. 4A). The 69 differentially expressed proteins were annotated by UniProt-GOA database or InterProScan. The proteins were enriched by GO annotation based on three categories: biological process, cellular component, and molecular function. The differentially expressed proteins are primarily enriched in cellular process (66.67%), metabolic process (55.07%), and single-organism process (53.62%). The molecular functions of the differentially expressed proteins are mainly involved in binding (52.17%) and catalytic activity (42.03%). In addition, the differentially expressed proteins are the component of cell part (86.96%), cell (86.96%), and organelle (65.22%) (Fig. 4B). Kyoto Encyclopedia of Genes and Genomes (KEGG) online service tools KAAS were used to annotate protein's KEGG database description. Then the annotations were mapped on KEGG pathway database using KEGG mapper. KEGG enrichment analysis revealed the differentially expressed proteins were enriched in Glucagon signaling pathway, Inflammatory mediator regulation of TRP channels, Aldosterone synthesis and secretion, ECM-receptor interaction, cAMP signaling pathway, Oxidative phosphorylation, Ca2+ signaling pathway, MAPK signaling, and NF-κB signaling pathway (Fig. 5). The differentially expressed proteins were introduced into the web-tool STRING to generate protein-protein interaction networks (Fig. 6).Fig. 6Interactions of differentially expressed proteins. Specific interaction networks analysis of differentially expressed proteins. Network nodes represent proteins. Splice isoforms or post-translational modifications are collapsed, i.e. each node represents the protein produced by a single, protein-coding gene locus. Colored nodes represent query proteins and the first shell of interactors. White nodes represent the second shell of interactions. Small nodes are the proteins with unknown 3D structure. Large nodes are the proteins with known or predicted 3D structure. Lines represent protein-protein interactions. Known Interactions: represents from curated databases; represents experimentally determined. Predicted Interactions: represents gene neighborhood; represents gene fusions; gene co-occurrence. Others: represents textmining; represents co-expression; represents protein homology.View Large Image Figure ViewerDownload Hi-res image Download (PPT) The transcription and expression levels of some differentially expressed proteins were verified using Western blot and real-time PCR. The Western blot analysis revealed that the protein levels of IL1RAP and RELA in PK-15 cells were dramatically increased by miR-4331 overexpression (Fig. 7A). The mRNA levels of 9 mitochondria-related differentially expressed proteins, including IL1RAP, COX5A, NDUFB5, MARPL24, MARPL1, RELA, CAMKA2D, PCK2, and BCAT1, were confirmed by real-time PCR (Fig. 7B). The Western blot and real-time PCR results were consistent with the proteomic data. According proteomic analysis, overexpression of miR-4331 may lead to an increase of IL1RAP protein level. It is reported that IL1RAP is a component of IL-1 protein complex (22.Boraschi D. Tagliabue A. The interleukin-1 receptor family.Semin. Immunol. 2013; 25: 394-407Crossref PubMed Scopus (172) Google Scholar) and IL-1 is a activator of p38 MAPK (23.Wang X. Wu H. Miller A.H. Interleukin 1alpha (IL-1alpha) induced activation of p38 mitogen-activated protein kinase inhibits glucocorticoid receptor function.Mol. Psychiatry. 2004; 9: 65-75Crossref PubMed Scopus (119) Google Scholar). We previously found that TGEV infection reduced MMP and activated p38 MAPK pathway (5.Ding L. Zhao X. Huang Y. Du Q. Dong F. Zhang H. Song X. Zhang W. Tong D. Regulation of ROS in transmissible gastroenteritis virus-activated apoptotic signaling.Biochem. Biophys. Res. Commun. 2013; 442: 33-37Crossref PubMed Scopus (27) Google Scholar). Therefore, we presume that IL1RAP might play a role in TGEV-induced decrease of MMP and activation of p38 MAPK. To explore the effects of IL1RAP on mitochondria and p38 MAPK pathway, we constructed pCI-neo-IL1RAP to overexpress IL1RAP and synthesized siRNA of IL1RAP, siIL1RAP, to silence IL1RAP. The results showed that pCI-neo-IL1RAP led to an increase IL1RAP at mRNA level (Fig. 8A) and protein level (Fig. 8B) and siIL1RAP reduced mRNA level (Fig. 8A) and protein level (Fig. 8B). Mitochondrial Ca2+ level was upregulated by pCI-neo-IL1RAP and downregulated by siIL1RAP (Fig. 8C). In addition, MMP level was suppressed by pCI-neo-IL1RAP and increased by siIL1RAP (Fig. 8D). As expected, phosphorylation level of p38 was enhanced by pCI-neo-IL1RAP and attenuated by siIL1RAP (Fig. 8E). To investigate whether the activation of p38 MAPK promotes TGEV-induced mitochondrial damage, PK-15 cells were treated with SB203580, the specific inhibitor of p38 MAPK, and infected with TGEV. As predicted, TGEV-induced phosphorylation of p38 was attenuated by SB203580 (Fig. 9A). Inhibition of p38 phosphorylation led to a decrease of mitochondrial Ca2+ level (Fig. 9B) and an increase of MMP during TGEV infection (Fig. 9C), indicating that p38 promotes TGEV-induced mitochondrial damage. We demonstrated that miR-4331 could upregulated IL1RAP expression via quantitative proteomic analysis and IL1RAP contributed to activation of p38 MAPK pathway via phosphorylating p38. However, whether miR-4331 can activate p38 pathway is unknown. To investigate the effect of miR-4331 on p38 MAPK, miR-4331 mimics or inhibitors was introduced into PK-15 cells. The phosphorylation level of p38 was detected. The result showed that TGEV-induced phosphorylation of p38 was enhanced by miR-4331 mimics and attenuated by miR-4331 inhibitors (Fig. 10), suggesting miR-4331 caused activation of p38 MAPK pathway. To identify whether miR-4331 directly binds to the 3′ UTRs of miR-4331 targets, the miR-4331 binding sites of target mRNA were mutated with a 4-bp substitution (Fig. 11A). The wild type sequences of 3′ UTRs containing miR-4331 binding sites and mutated 3′ UTRs sequences of miR-4331 targets, in which miR-4331 seed sequence was mutated, were cloned into the 3′ UTR of Renilla luciferase in dual-luciferase reporter plasmid psiCHECK-2 to generate constructions, RB1-wt and RB1-mut (Fig. 11B). The constructs were co-transfected into PK-15 cells with either miR-4331 mimics/mimics control or miR-4331 inhibitors/inhibitors control. The Renilla luciferase activities of RB1-wt and RB1-mut were respectively reduced and improved by miR-4331 mimics (normalized to Firefly luciferase activity) (Fig. 11C). In contrast, the Renilla luciferase activities of RB1-wt and RB1-mut were respectively raised and decreased by miR-4331 inhibitors (normalized to Firefly luciferase activity) (Fig. 11D). The results revealed that miR-4331 directly binds to 3′ UTR of RB1 mRNA. To determine whether miR-4331 inhibits RB1 expression, either miR-4331 mimics or miRNA inhibitors control were transfected into PK-15 cells. RB1 protein level was tested using Western blot, showing that expression of RB1 was suppressed by miR-4331 mimics and improved by miR-4331 inhibitors (Fig. 11E). We demonstrated that miR-4331 facilitated mitochondrial damage and directly targeted RB1. However, whether miR-4331 affects mitochondria through targeting RB1 is unclear. To investigate the effect of RB1 on TGEV-induced mitochondrial damage, three siRNAs of RB1, siRB1–1, siRB1–2, and siRB1–3, were synthesized and respectively transfected into PK-15 cells to silence RB1. The silencing efficiency of RB1 siRNAs were evaluated, showing siRB1–3 was the most effective siRNA for RB1 silencing (Fig. 12A and 12B). Moreover, RB1 gene was cloned into eukaryotic expression plasmid pCI-neo, named pCI-neo-RB1, to overexpress RB1. pCI-neo-RB1 and pCI-neo were respectively transfected into PK-15 cells, followed by evaluation of mRNA and protein level. Expectedly, both the mRNA and protein level were upregulated (Fig. 12C and 12D). Then, PK-15 cells were transfected with either siRB1–3 or pCI-neo-RB1 and subsequently infected with TGEV for 24 h. The mitochondrial Ca2+ level and MMP were measured. The results showed that mitochondrial Ca2+ level was decreased by pCI-neo-RB1 and increased by siRB1–3 (Fig. 12E). In addition, MMP was upregulated by pCI-neo-RB1 and decreased by siRB1–3 (Fig. 12F). The results reveal that miR-4331 promotes TGEV-induced mitochondrial damage via targeting RB1. We showed that miR-4331 could aggravate TGEV-induced mitochondrial damage via targeting RB1 and activating p38 MAPK pathway. Therefore, we presumed that miR-4331 is likely to regulate IL1RAP and p38 MAPK signaling through its target RB1. We silenced RB1 using siRB1–3 and overexpressed RB1 using pCI-neo-RB1. IL1RAP, p38, and p-p38 were analyzed by Western blotting. It was found that RB1 was suppressed by siRB1–3 (Fig 13). Silencing RB1 caused an upregulation of IL1RAP and p-p38 (Fig. 13). Reversely, overexpression of RB1 led to a suppress of IL1RAP and p-p38 (Fig. 13). In this study, we investigated the effects of miR-4331 on mitochondrial damage induced by TGEV infection via targeting RB1. The results showed that miR-4331 promoted TGEV-induced mitochondrial damage via targeting RB1, upregulating IL1RAP, and enhancing TGEV-induced activation of p38 pathway. Our findings reveal a novel regulatory effect of miR-4331 on TGEV-induced mitochondrial damage. It is reported that TGEV infection caused the alterations of proteomes in PK-15 cells, including upregulated STAT1 and many ISGs (ISG15, IFIT1, IFIT2, IFIT3, IFIT5, OAS1, OAS2, and Mx1) (24.An K. Fang L. Luo R. Wang D. Xie L. Yang J. Chen H. Xiao S. Quantitative proteomic analysis reveals that transmissible gastroenteritis virus activates the JAK-STAT1 signaling pathway.J. Proteome Res. 2014; 13: 5376-5390Crossref PubMed Scopus (42) Google Scholar). In this study, Mx1 was differentially decreased by miR-4331 in contrast to upregulation in TGEV-infected PK-15 cells. It indicates that miR-4331 plays a negative role in TGEV-induced Mx1 expression. In addition, although STAT1 and many ISGs, including ISG15, IFIT1, IFIT2, IFIT3, OAS1, OAS2, Mx1, were detected using TMT approach, only Mx1 was altered by miR-4331 in TGEV-infected PK-15 cells. The likely reason is that miR-4331 downregulated TGEV-induced STAT1 and these ISGs, including ISG15, IFIT1, IFIT2, IFIT3, OAS1, OAS2, to normal level. In contrast to upregulation in TGEV-infected PK-15 cells, Mx1 was differentially decreased by miR-4331 in TGEV-infected PK-15 cells. It indicates that miR-4331 plays a negative role in TGEV-induced Mx1 expression. In this study, RB1 is identified as the target of miR-4331 and affects mitochondrial function. RB1 is a suppressor of cell death via binding and repressing transcription factor E2Fs. The hypophosphorylated RB1, the active form of RB1, binds to and sequesters transcription factor E2F1 to operate as an inhibitor of E2Fs (25.Gordon G.M. Du W. Conserved RB functions in development and tumor suppression.Protein Cell. 2011; 2: 864-878Crossref PubMed Scopus (40) Google Scholar). C-terminal region of RB1 has a E2F1-specific binding site that is sufficient to repress E2F1-induced apoptosis, so the C-terminal interaction site on RB1 acts as a potent inhibitor of E2F1 (26.Dick F.A. Dyson N. pRB contains an E2F1-specific binding domain that allows E2F1-induced apoptosis to be regulated separately from other E2F activities.Mol. Cell. 2003; 12: 639-649Abstract Full Text Full Text PDF PubMed Scopus (114) Google Scholar). A pool of evidence shows RB1 localizes in mitochondria (27.Hilgendorf K.I. Leshchiner E.S. Nedelcu S. Maynard M.A. Calo E. Ianari A. Walensky L.D. Lees J.A. The retinoblastoma protein induces apoptosis directly at the mitochondria.Genes Dev. 2013; 27: 1003-1015Crossref PubMed Scopus (88) Google Scholar). Therefore, a set of reports highlights a series of links between RB1 and mitochondrial function. RB1 promotes mitochondrial biogenesis for erythropoiesis (28.Sankaran V.G. Orkin S.H. Walkley C.R. Rb intrinsically promotes erythropoiesis by coupling cell cycle exit with mitochondrial biogenesis.Genes Dev. 2008; 22: 463-475Crossref PubMed Scopus (113) Google Scholar). Loss of RB1 caused a decrease of mitochondrial mass, downregulated mitochondrial function, oxidative phosphorylation, MMP, and accumulation of hypopolarized mitochondria (29.Nicolay B.N. Danielian P.S. Kottakis F. Lapek Jr, J.D. Sanidas I. Miles W.O. Dehnad M. Tschop K. Gierut J.J. Manning A.L. Morris R. Haigis K. Bardeesy N. Lees J.A. Haas W. Dyson N.J. Proteomic analysis of pRb loss highlights a signature of decreased mitochondrial oxidative phosphorylation.Genes Dev. 2015; 29: 1875-1889Crossref PubMed Scopus (65) Google Scholar), implying that suppression of RB1 leads to mitochondrial damage. These results further support our conclusion that RB1 weakens TGEV-induced mitochondrial damage. Here, RB1 was identified as the target of miR-4331 and inhibited by miR-4331. Whereas, RB1 was not identified as the differentially expressed protein using TMT, it may be because of low expression level of RB1 and the insufficient sensitivity of this technology. Cytokines of the IL-1 family are considered critical regulators of intestinal mucositis (30.de Mooij C.E.M. Netea M.G. van der Velden W. Blijlevens N.M.A. Targeting the interleukin-1 pathway in patients with hematological disorders.Blood. 2017; 129: 3155-3164Crossref PubMed Scopus (44) Google Scholar). IL-1 is a key proinflammatory cytokine that can initiate a series of signaling resulting in activation of NF-κB signaling pathway. IL-1 interacts with L1RAP and IL-1 receptor type I (IL-1RI) to form IL-1 complex (31.Acuner Ozbabacan S.E. Gursoy A. Nussinov R. Keskin O. The structural pathway of interleukin 1 (IL-1) initiated signaling reveals mechanisms of oncogenic mutations and SNPs in inflammation and cancer.PLoS Comput. Biol. 2014; 10: e1003470Crossref PubMed Scopus (52) Google Scholar), which mediates the activation of p38 MAPK signaling in presence of IL-1β treatment (32.Chen R. Li M. Zhang Y. Zhou Q. Shu H.B. The E3 ubiquitin ligase MARCH8 negatively regulates IL-1beta-induced NF-kappaB activation by targeting the IL1RAP coreceptor for ubiquitination and degradation.Proc. Natl. Acad. Sci. U.S.A. 2012; 109: 14128-14133Crossref PubMed Scopus (75) Google Scholar). It is in consistent with our finding that IL1RAP can lead to the activation of p38 MAPK pathway. Moreover, we found miR-4331 caused an upregulation of RELA, which is the p65 subunit of NF-κB. NF-κB is a key activator of inflammation and is activated during TGEV infection (33.Ding Z. An K. Xie L. Wu W. Zhang R. Wang D. Fang Y. Chen H. Xiao S. Fang L. Transmissible gastroenteritis virus infection induces NF-kappaB activation through RLR-mediated signaling.Virology. 2017; 507: 170-178Crossref PubMed Scopus (36) Google Scholar, 34.Guo J. Li F. Qian S. Bi D. He Q. Jin H. Luo R. Li S. Meng X. Li Z. TGEV infection upregulates FcRn expression via activation of NF-kappaB signaling.Sci Rep. 2016; 6: 32154Crossref PubMed Scopus (21) Google Scholar). Therefore, we speculate miR-4331 may function as a regulator of TGEV-induced NF-κB pathway through regulating IL1RAP and RELA. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD008174 (http://www.ebi.ac.uk/pride/archive/) (18). Download .zip (4.22 MB) Help with zip files