Literature DB >> 26249758

RNA Profiling in Human and Murine Transplanted Hearts: Identification and Validation of Therapeutic Targets for Acute Cardiac and Renal Allograft Rejection.

L N L Van Aelst1, G Summer2, S Li3, S K Gupta4, W Heggermont1, K De Vusser5, P Carai1,2, M Naesens5, J Van Cleemput1, F Van de Werf1, J Vanhaecke1, T Thum4, M Waer3, A-P Papageorgiou1,2, B Schroen2, S Heymans1,2,6.   

Abstract

Acute cellular rejection (ACR) is the adverse response of the recipient's immune system against the allogeneic graft. Using human surveillance endomyocardial biopsies (EMBs) manifesting ACR and murine allogeneic grafts, we profiled implicated microRNAs (miRs) and mRNAs. MiR profiling showed that miR-21, -142-3p, -142-5p, -146a, -146b, -155, -222, -223, and -494 increased during ACR in humans and mice, whereas miR-149-5p decreased. mRNA profiling revealed 70 common differentially regulated transcripts, all involved in immune signaling and immune-related diseases. Interestingly, 33 of 70 transcripts function downstream of IL-6 and its transcription factor spleen focus forming virus proviral integration oncogene (SPI1), an established target of miR-155, the most upregulated miR in human EMBs manifesting rejection. In a mouse model of cardiac transplantation, miR-155 absence and pharmacological inhibition attenuated ACR, demonstrating the causal involvement and therapeutic potential of miRs. Finally, we corroborated our miR signature in acute cellular renal allograft rejection, suggesting a nonorgan specific signature of acute rejection. We concluded that miR and mRNA profiling in human and murine ACR revealed the shared significant dysregulation of immune genes. Inflammatory miRs, for example miR-155, and transcripts, in particular those related to the IL-6 pathway, are promising therapeutic targets to prevent acute allograft rejection. © Copyright 2015 The American Society of Transplantation and the American Society of Transplant Surgeons.

Entities:  

Keywords:  basic (laboratory) research/science; heart (allograft) function/dysfunction; heart transplantation/cardiology; microarray/gene array; molecular biology: mRNA/mRNA expression; molecular biology: micro RNA; rejection: acute

Mesh:

Substances:

Year:  2015        PMID: 26249758      PMCID: PMC5054886          DOI: 10.1111/ajt.13421

Source DB:  PubMed          Journal:  Am J Transplant        ISSN: 1600-6135            Impact factor:   8.086


acute cellular rejection endomyocardial biopsy guanylate‐binding proteins heart failure heart transplantation in situ hybridization International Society of Heart and Lung Transplantation microRNA postoperative day suppressor of cytokine signaling 1 spleen focus forming virus proviral integration oncogene transporter associated with antigen processing TCRα subunit constant gene 3′ untranslated region

Introduction

Heart transplantation (HTX) is the only curative treatment available for end‐stage heart failure (HF). Despite contemporary immunosuppressive regimens, ACR occurs in approximately one quarter of recipients during the first year and accounts for 10% of deaths between 1 month and 1 year following transplantation 1. Furthermore, ACR is a risk factor for 5‐year posttransplant mortality and cardiac allograft vasculopathy, the latter itself a major impediment to survival beyond 3 years after transplantation 1. Surveillance endomyocardial biopsies (EMBs) are the standard of care to diagnose acute rejection before graft function declines based on histopathological evidence of inflammation and cardiomyocyte necrosis; the International Society of Heart and Lung Transplantation (ISHLT) issues guidelines to standardize pathology reports 2. Despite decreasing rates of ACR due to improved immunosuppressants 1, new insights into the cellular changes associated with ACR remain important to diagnose ACR earlier, to identify patients at risk for ACR and to develop novel therapeutics without side effects 3. In the past decade, miRs emerged as a popular gene family to serve as a target in a wide range of pathologies including cardiovascular diseases. MiRs are phylogenetically conserved, noncoding RNAs which inhibit mRNA translation to protein by interacting with its 3′ untranslated region (UTR) 4. They are essential for cardiac development and homeostasis; miRs are also crucially involved in cardiac fibrosis, hypertrophy, and electrical remodeling 5. Several reports demonstrated altered serological and tissular miR expression following rejection of cardiac 6, 7, hepatic 8, renal 9, 10, 11, and intestinal 12, 13 allografts, yet these studied only one species, restricted the analysis to miR profiling without mRNA profiling and did not include a therapeutic intervention to demonstrate a causal role for specific miRs. One very recent report implicated miR‐155 in murine cardiac allograft rejection through its interaction with suppressor of cytokine signaling 1 (SOCS1) in dendritic cells. Apart from not including a translation to human medicine, this study lacked unbiased miR and mRNA expression analysis 14. In our study, we first determined common differentially expressed miRs and mRNAs in human ACR and a mouse model of ACR. Next, we unraveled mRNA pathways involved in ACR despite optimal therapy. Notably, the IL‐6 pathway was markedly induced. MiR‐155 expression was strongly upregulated in both humans and mice and its targeted genomic deletion or pharmacological inhibition attenuated inflammation and improved graft survival in a murine model. We concluded that RNA targeting, in particular miR‐155 and the IL‐6 functional pathway more specifically, represent an attractive novel treatment option for ACR.

Materials and Methods

Animals and surgery

Experiments were approved by the ethical committee of the KULeuven (Leuven, Belgium) and performed according to the Belgian legislation on laboratory animals. Eight to ten‐week old BALB/cJ (Janvier), C57Bl/6J (Janvier), bicm2/m2 (miR155−/−) and WT littermates (miR155+/+) were subjected to cervical heterotopic HTX and sham operations 15, 16.

Human endomyocardial biopsies

Human material was obtained during sampling for clinical purposes and available for research according to the Declaration of Helsinki and the ethical committees at UZ Leuven (Leuven, Belgium) and Maastricht University Medical Center (Maastricht, the Netherlands). Surveillance EMBs included in our analyses were preferentially taken later than 6 weeks after transplantation to avoid the perioperative phase. Control biopsies consisted of age‐matched patients with unexplained ventricular tachyarrhythmias yet with a normal ejection fraction, cardiac morphology, and no systemic or cardiac inflammation or viral persistence. Biopsies were snap‐frozen for RNA analysis and formalin or Bouin fixed for histology.

Histological and morphometrical analysis

Mice were sacrificed 3, 5, or 7 days after HTX. Organs were removed, rinsed, blotted dry, weighed, and snap frozen. Before weighing, atria were separated from the ventricles. Ventricles were divided in basal, midventricular, and apical parts. The apical part was split and snap frozen. The midventricular part was used for in situ hybridization (ISH). The basal part was fixed in 4% paraformaldehyde (PFA; Sigma, Diegem, Belgium) in phosphate‐buffered saline (PBS) overnight, post‐fixed in ethanol, dehydrated, and embedded in paraffin. Four micrometer thick sections were made for histological (hematoxylineosin) and immunohistochemical stainings 17. Immunohistochemical stainings were performed with antibodies directed against CD45 (BD Biosciences, Erembodegem, Belgium), CD3 (AbD Serotec, Kidlington, UK), Mac3 (BD Biosciences). Omission of the primary antibody served as negative control.

In situ hybridization

The midventricular part of the harvested tissue was fixed in 10% neutral buffered formalin for 48 h, postfixed in ethanol, dehydrated, and embedded in paraffin. ISH was performed as described previously 18. Briefly, sections were dewaxed, dehydrated, air‐dried, and incubated in 1.3 mg/mL Pepsin/0.1 mM HCl solution for 5 min at 37°C. Following a wash in MQ, sections were incubated at 42.5°C overnight in 1× hybridization buffer (ENZO Life Sciences, Antwerpen, Belgium) with 100 nM mmu‐miR‐155 double DIG‐labeled probe or scrambled control (Exiqon). Afterwards, sections were washed in 5xSSC and in 0.2×SSC at hybridization temperature. Subsequently, sections were blocked in DIG blocking buffer from the DIG Wash and Block Buffer Set (Roche, Vilvoorde, Belgium) for 20 min followed by incubation in 1:500 anti‐DIG‐AP, Fab fragments (Roche) for 45 min. Sections were washed 3 × 5 min with TBS and incubated 3 min in DIG detection buffer at room temperature. AP signal was detected using NBT/BCIP tablets (Roche) for 4 h at 30° in darkness. Finally, sections were counterstained with fast red (Sigma), washed in MQ, dehydrated through a graded ethanol series and mounted with entellan (EMS).

Western blotting

Tissue was homogenized in RIPA buffer; phosphatase (Sigma) and protease (Complete®, Roche) inhibitors were added. Lysates were resolved by SDS–PAGE and immunoblotted for SPI1/PU.1 (Cell Signaling Technology, Leiden, the Netherlands) overnight at 4°C. GAPDH (Fitzgerald Industries, Huissen, the Netherlands) was used as a loading control. Imaging was performed on a ChemiDoc XRS+ System (BioRad, Temse, Belgium). Quantification was done using Image J software (NIH).

RNA isolation and expression

RNA was isolated with the miRVana kit (Ambion, Gent, Belgium). We controlled its integrity and concentration using Nanodrop (ThermoScientific, Erembodegem, Belgium) and BioAnalyzer (Agilent, Diegem, Belgium). For murine samples, 1 μg of total RNA was reverse transcribed with the miScript cDNA synthesis kit (Qiagen, Venlo, the Netherlands); for human samples, 0.5 μg of total RNA was reverse transcribed. Real‐time quantitative PCR was performed with SYBR green mix (Applied Biosystems, Gent, Belgium) on an ABI Prism 7500 (Applied Biosystems). LNA primers (Exiqon, Vedbaek, Denmark) were used to detection mature miRs. U6 was used as an endogenous control for miRs; GAPDH was used as a housekeeping gene for mRNAs.

MiR and mRNA expression profiling

Total RNA from human biopsies was subjected to miR expression profiling on NanoStrings' nCounter platform by the Nucleomics Core of the Flemish Institute of Biotechnology (Merelbeke, Belgium). The gene target list consisted of 654 human miRs, 80 human‐associated viral miRs, 6 positive controls, 8 negative controls, and 5 housekeeping genes. Further identification of differentially regulated miRs was performed using the R pipeline provided by the manufacturer. Total RNA from mouse hearts was hybridized to mouse miRNA microarray G4472B (Agilent) based on miRBase release 12.0, containing 627 mouse miRs, 39 mouse viral miRs, 6 negative controls, 6 positive controls, and 5 housekeeping genes. Further identification of differentially regulated miRs was performed using software provided by the manufacturer at Hannover Medical School (MHH, Hannover, Germany). Total RNA from human biopsies and mouse hearts was subjected to mRNA expression profiling on Affymetrix Human PrimeView and Mouse MoGene 1.0 ST v1 microarrays, respectively. The expression data was normalized using the Affymetrix R package and appropriate CustomCDFs from the Microarray Lab (University of Michigan) 19. Differential expression was assessed in R software(20) using the limma package 21. MiRs were considered differentially expressed with a p‐value of <0.05; mRNAs were considered differentially expressed with a p‐value of <0.05 and a fold change of >2.0.

Design, synthesis, and application of antagomiR‐155 and scram‐antagomiR

AntagomiRs were synthesized by Fidelity Systems (Gaithersburg, MD) 22. The antagomiR‐155 (anti‐155) oligonucleotide 5′‐cscsccuaucacaauuagcaususasas‐Chol‐3′ is complementary to nucleotides 4‐25 in mmu‐miR‐155. The scrambled antagomiR (anti‐scram) has the following sequence: 5′‐uscsugacccuaaacauucaascsusas‐Chol‐3′. Letters represent 2′‐O‐Me‐modified nucleotides; subscript ‘s’ represent a phosphorothioate linkage; ‘Chol’ represents cholesterol linked through a hydroxyprolinol linkage.

Statistical analysis

Data are expressed as mean ± SEM. Normal distribution of continuous variables was tested with the Kolmogorov–Smirnov method. An unpaired Student's t test was used when groups passed the normality test, otherwise, a Mann–Whitney test was used. ANOVA with Bonferroni posthoc analysis was performed when more than two groups were compared. Graft outcome data were compared by the log rank test. A two‐sided p‐value of <0.05 was considered statistically significant.

Results

MiR expression profile in murine and human ACR following HTX

We compared miR expression in human biopsies from control hearts (n = 7) and surveillance EMBs manifesting severe rejection (ISHLT 3; n = 7), from patients during the first year following HTX (Figure 1A).
Figure 1

ACR of human cardiac allografts. (A) Hematoxylin‐eosin staining of human EMBs manifesting different grades of rejection (ISHLT 0 = no rejection; ISHLT 1 = mild rejection; ISHLT 3 = moderate to severe rejection). (B) Heat map of differentially expressed miRs in control hearts and ISHLT 3 grafts. (C) Volcano plot showing fold change (log2 values) and probability (log10 values) for individual miRs, comparing control and ISHLT 3 hearts. Common dysregulated miRs in humans and mice are highlighted. ACR, acute cellular rejection; EMB, endomyocardial biopsy; ISHLT, International Society of Heart and Lung Transplantation.

ACR of human cardiac allografts. (A) Hematoxylineosin staining of human EMBs manifesting different grades of rejection (ISHLT 0 = no rejection; ISHLT 1 = mild rejection; ISHLT 3 = moderate to severe rejection). (B) Heat map of differentially expressed miRs in control hearts and ISHLT 3 grafts. (C) Volcano plot showing fold change (log2 values) and probability (log10 values) for individual miRs, comparing control and ISHLT 3 hearts. Common dysregulated miRs in humans and mice are highlighted. ACR, acute cellular rejection; EMB, endomyocardial biopsy; ISHLT, International Society of Heart and Lung Transplantation. Microarray analysis identified 47 significantly differentially regulated miRs in EMBs with severe rejection compared with control biopsies, of which 20 were upregulated and 27 downregulated; unsupervised clustering of samples based on miR expression results in perfect discrimination of control and ISHLT 3 samples (Figure 1B and C; Table S1). Pathway analysis revealed that the 47 dysregulated human miRs have previously been implicated in cancer (both in situ growth and metastasis), inflammation, and autoimmune disorders (Table S2). In parallel, we studied differential miR expression in murine ACR using a full MHC haplotype‐mismatched model of HTX involving male BALB/cJ hearts (MHC haplotype H‐2d) to male C57Bl/6J mice (H‐2b)(Figure 2A). In the absence of immunosuppression, graft failure, defined as the cessation of beating of the graft, occurred on the seventh postoperative day (Allo 7d) while inflammation was already present at postoperative day 3 (Allo3d; Figure 2B–E). Cardiac miR expression profiles were analyzed in Allo3d and Allo7d grafts. Here, 173 miRs were differentially regulated in Allo7d versus untransplanted BALB/cJ hearts, of which 86 were upregulated and 87 were downregulated (Figure 2F and G; Table S3). These dysregulated miRs have previously been implicated in the pathogenesis of infectious diseases, diabetes mellitus, and nonalcoholic fatty liver disease (Table S4).
Figure 2

ACR of murine cardiac allografts and comparison of human and murine array results. (A) Schema of included groups to study the murine miR expression profile in ACR following HTX. (B) Representative CD45 staining of corresponding regions in control hearts and cardiac allografts, 3d (Allo3d) and 7d (Allo7d) after HTX, with quantification in (C). (D and E) depict CD4 and Mac3 staining, respectively. (F) Heat map of differentially expressed miRs in murine control hearts and failed grafts. (G) Volcano plot showing fold change (log2 values) and probability (log10 values) for individual miRs, comparing control hearts and Allo7d. (H) ISH of miR‐155 during murine ACR showed miR‐155 staining in inflammatory cells (black arrowheads), whereas adjacent ventricular myocardium didn't stain (white arrowheads). (I) Venn diagrams showing the numbers and overlap of differentially regulated miRs between humans and mice. (J) Fold changes for all miRs in humans and mice during ACR; common dysregulated miRs in humans and mice during ACR are highlighted. (K) Following HTX, splenic SPI1 protein levels markedly decreased, whereas splenic RNA levels significantly increased (L) (*p < 0.05). ACR, acute cellular rejection; HTX, heart transplantation.

ACR of murine cardiac allografts and comparison of human and murine array results. (A) Schema of included groups to study the murine miR expression profile in ACR following HTX. (B) Representative CD45 staining of corresponding regions in control hearts and cardiac allografts, 3d (Allo3d) and 7d (Allo7d) after HTX, with quantification in (C). (D and E) depict CD4 and Mac3 staining, respectively. (F) Heat map of differentially expressed miRs in murine control hearts and failed grafts. (G) Volcano plot showing fold change (log2 values) and probability (log10 values) for individual miRs, comparing control hearts and Allo7d. (H) ISH of miR‐155 during murine ACR showed miR‐155 staining in inflammatory cells (black arrowheads), whereas adjacent ventricular myocardium didn't stain (white arrowheads). (I) Venn diagrams showing the numbers and overlap of differentially regulated miRs between humans and mice. (J) Fold changes for all miRs in humans and mice during ACR; common dysregulated miRs in humans and mice during ACR are highlighted. (K) Following HTX, splenic SPI1 protein levels markedly decreased, whereas splenic RNA levels significantly increased (L) (*p < 0.05). ACR, acute cellular rejection; HTX, heart transplantation. Subsequently, we compared murine and human miR expression profiles of ACR. Our analysis revealed 10 consistently differentially regulated miRs in both species: miR‐21, miR‐146a, miR‐146b, miR‐155, miR‐142‐3p, miR‐142‐5p, miR‐222, miR‐223, and miR‐494 (consistently upregulated), and miR‐149‐5p (consistently downregulated) (Figure 2I and J; Table 1). We validated their differential expression by qPCR in mice and humans. In mice, we compared control hearts, Allo3d and Allo7d grafts; qPCR results were in line with the microarray expression profiles and proved the time‐dependent change in the expression of the 10 miRs, paralleling the severity of inflammation (Figure S1). In humans, we compared EMBs manifesting ISHLT 0, 1, and 3 rejection for validation qPCRs (n = 19, 20, and 10, respectively; Table S5 summarizes patient characteristics). Validation qPCRs revealed that miRs previously linked to inflammation paralleled histopathological ACR severity (miR‐142‐3p, miR‐142‐5p, miR‐146a, and miR‐155) or were able to separate EMBs with rejection from those without ACR (miR‐223) (Figure S2).
Table 1

Ten differentially regulated miRs in mice and humans during ACR

Homo sapiens Mus musculus
Fold change ISHLT 3 versus controlp‐valueFold change Allo 7d versus controlp‐value
miR‐15516.959.38 10−5 5.571.77 10−5
miR‐146a6.919.71 10−4 3.421.96 10−4
miR‐219.671.39 10−3 10.207.72 10−6
miR‐142‐5p4.801.48 10−3 13.542.99 10−5
miR‐2238.066.00 10−3 5.192.14 10−4
miR‐142‐3p9.516.65 10−3 11.572.55 10−5
miR‐146b7.101.99 10−2 255.717.83 10−12
miR‐2223.812.47 10−2 230.261.83 10−3
miR‐4941.932.66 10−2 9.585.94 10−4
miR‐1490.204.57 10−2 0.101.25 10−5

ACR, acute cellular rejection; ISHLT, International Society of Heart and Lung Transplantation; miR, microRNA.

Ten differentially regulated miRs in mice and humans during ACR ACR, acute cellular rejection; ISHLT, International Society of Heart and Lung Transplantation; miR, microRNA. In conclusion, cardiac ACR induced distinct changes in miR expression in murine and human grafts. MiR‐155 is upregulated in human and murine ACR, and the most upregulated miR in human grafts manifesting rejection.

mRNA expression profile in murine and human ACR following HTX

Apart from miR profiling, we performed comparative mRNA arrays in both human and murine control hearts and rejected cardiac allografts. Only transcripts with a significant twofold up or downregulation in rejected allografts were considered differentially expressed. This analysis revealed 70 commonly differentially regulated mRNAs. Remarkably, several highly upregulated common mRNAs in Allo7d grafts were already significantly upregulated in Allo3d grafts. For example, TCRα subunit constant gene (Trac), the most differentially expressed gene in Allo7d grafts, was already 3.6‐fold upregulated in Allo3d grafts. Other examples include Cxcl9 and ‐10, the Guanylate‐binding proteins (Gbp) 1 to −5, transporter associated with antigen processing (Tap) 1, CD74, and Stat1, all of which are involved in immune signaling. Upstream regulator and mechanistic network analysis revealed that 33/70 common transcripts are involved in downstream signaling of IL‐6 and its transcription factor SPI1 (Figures 3, 4; Table S6–8). Interestingly, SPI1 has a conserved binding site for miR‐155 in its 3′ UTR 18.
Figure 3

ACR of human and murine cardiac allografts: mRNA array results. (A) Volcano plot showing fold changes (log2 values) and probability (log10 values) for individual human mRNAs, comparing control and ISHLT3 EMBs. (B) Volcano plot showing fold changes (log2 values) and probability (log10 values) for individual murine mRNAs, comparing control hearts and Allo7d. (C) Volcano plot showing fold changes (log2 values) and probability values (log10) for individual murine mRNAs, comparing control hearts and Allo 3d. ACR, acute cellular rejection; ISHLT, International Society of Heart and Lung Transplantation; mRNA, messenger RNA.

Figure 4

Differentially regulated human and mouse mRNA transcripts from the micro array studies. (A) Pathway enrichment analysis of human mRNA arrays (ISHLT3 vs. control) showed significant enrichment of predominantly immune‐related pathways (Table S7). The heat map, showing in rows all genes represented in the enriched pathways, confirm that the ISHLT3 cardiac samples are highly regulated except for one. In addition, gene expression in ISHLT0 cardiac samples does not significantly differ from control gene expression, as becomes apparent from the mixed unsupervised clustering of these groups. (B) Pathway enrichment analysis of murine mRNA arrays (Allo7d vs. control) showed significant enrichment of predominantly metabolic and immune‐related pathways (Table S8). The heat maps, showing in rows all genes represented in these pathways, confirm that the Allo7d samples are highly regulated and that the Allo3d samples are in a transition phase in which the respective genes are mildly regulated. Moreover, the unsupervised clustering of the samples based on the gene expression of these pathways results in the perfect discrimination of control, Allo3d, and Allo7d samples. ISHLT, International Society of Heart and Lung Transplantation; mRNA, messenger RNA.

ACR of human and murine cardiac allografts: mRNA array results. (A) Volcano plot showing fold changes (log2 values) and probability (log10 values) for individual human mRNAs, comparing control and ISHLT3 EMBs. (B) Volcano plot showing fold changes (log2 values) and probability (log10 values) for individual murine mRNAs, comparing control hearts and Allo7d. (C) Volcano plot showing fold changes (log2 values) and probability values (log10) for individual murine mRNAs, comparing control hearts and Allo 3d. ACR, acute cellular rejection; ISHLT, International Society of Heart and Lung Transplantation; mRNA, messenger RNA. Differentially regulated human and mouse mRNA transcripts from the micro array studies. (A) Pathway enrichment analysis of human mRNA arrays (ISHLT3 vs. control) showed significant enrichment of predominantly immune‐related pathways (Table S7). The heat map, showing in rows all genes represented in the enriched pathways, confirm that the ISHLT3 cardiac samples are highly regulated except for one. In addition, gene expression in ISHLT0 cardiac samples does not significantly differ from control gene expression, as becomes apparent from the mixed unsupervised clustering of these groups. (B) Pathway enrichment analysis of murine mRNA arrays (Allo7d vs. control) showed significant enrichment of predominantly metabolic and immune‐related pathways (Table S8). The heat maps, showing in rows all genes represented in these pathways, confirm that the Allo7d samples are highly regulated and that the Allo3d samples are in a transition phase in which the respective genes are mildly regulated. Moreover, the unsupervised clustering of the samples based on the gene expression of these pathways results in the perfect discrimination of control, Allo3d, and Allo7d samples. ISHLT, International Society of Heart and Lung Transplantation; mRNA, messenger RNA. In conclusion, ACR induced transcripts involved in immune signaling and immune‐related diseases. Specifically, components of the IL‐6 signaling pathway including the miR‐155 target SPI1 were similarly dysregulated in human and murine ACR.

Genomic absence of miR‐155 in cardiac inflammatory cells attenuates ACR

We performed in situ hybridization (ISH) on sections of murine inflamed Allo7d grafts to identify which cells produce miR‐155 during ACR. We observed perinuclear miR‐155 staining in pericardial inflammatory patches and infiltrating cells, whereas we did not detect any signal in adjacent cardiomyocytes (Figure 2H). These findings are in line with previous studies showing miR‐155 upregulation in leukocytes following activation 16, 18, 23. To support the pathogenic role of miR‐155 in ACR, we transplanted BALB/cJ hearts into C57Bl/6J mice with a targeted deletion of miR‐155 (miR155−/−) and into their WT littermates (miR155+/+) (Figure 5A). Median graft survival time increased significantly in miR155−/− versus miR155+/+ mice (miR‐155−/−: 10 days; miR‐155+/+: 7.5 days; p = 0.004) (Figure 5B). Despite longer exposure to the recipient's immune system and ensuing tissue edema, grafts from miR155−/− mice weighed significantly less upon rejection (miR‐155−/−: 10.1 mg/g ± 1.2; miR‐155+/+: 12.7 mg/g ± 1.7; p = 0.03) (Figure 5C). However, the amount of infiltrating inflammatory cells in these failed grafts, harvested at different postoperative time points, showed no significant between genotype difference (p = 0.6) (Figure S3A).
Figure 5

ACR in miR‐155 +/+ and miR‐155 −/− mice following HTX. (A) Experimental setup of miR‐155 +/+ and miR‐155 −/− mice in spontaneous survival experiments. (B) Graft failure occurred later when BALB/cJ hearts were transplanted into miR‐155 −/− versus miR‐155 +/+ mice. (C) Despite longer exposure to the host's immune system, BALB/cJ grafts weighed significantly less when transplanted into miR‐155 −/− mice. (D) Experimental setup of miR‐155 +/+ and miR‐155 −/− mice sacrificed on a fixed time point (5d) following HTX. (E) Whereas the number of infiltrating CD45 and Mac3‐positive cells was significantly less in BALB/cJ hearts transplanted into miR155 −/− versus miR‐155 +/+ mice, the number of CD3‐positive cells did not significantly differ between genotypes (*p < 0.05). ACR, acute cellular rejection; MST, median survival time; HTX, heart transplantation; miR, microRNA.

ACR in miR‐155 +/+ and miR‐155 −/− mice following HTX. (A) Experimental setup of miR‐155 +/+ and miR‐155 −/− mice in spontaneous survival experiments. (B) Graft failure occurred later when BALB/cJ hearts were transplanted into miR‐155 −/− versus miR‐155 +/+ mice. (C) Despite longer exposure to the host's immune system, BALB/cJ grafts weighed significantly less when transplanted into miR‐155 −/− mice. (D) Experimental setup of miR‐155 +/+ and miR‐155 −/− mice sacrificed on a fixed time point (5d) following HTX. (E) Whereas the number of infiltrating CD45 and Mac3‐positive cells was significantly less in BALB/cJ hearts transplanted into miR155 −/− versus miR‐155 +/+ mice, the number of CD3‐positive cells did not significantly differ between genotypes (*p < 0.05). ACR, acute cellular rejection; MST, median survival time; HTX, heart transplantation; miR, microRNA. In another cohort, we sacrificed all animals on postoperative day (POD) 5 (i.e. before occurrence of graft failure) (Figure 5D). Gravimetric analysis showed a numerically lower normalized graft and spleen weight in miR155−/− recipients, yet these differences did not reach statistical significance (grafts, miR155−/−: 5.3 mg/g ± 0.4 vs. miR155+/+: 6.2 mg/g ± 0.6, p = 0.2; spleens, miR‐155−/−: 4.6 mg/g ± 0.6 vs. miR155+/+: 5.6 mg/g ± 0.5; p = 0.3) (Figure S3B and C). More importantly, influx of CD45‐positive leukocytes and Mac3‐positive macrophages was significantly impaired in BALB/c hearts transplanted into miR155−/− recipients as opposed to miR155+/+ recipients (p = 0.03 and 0.04, respectively). CD3‐positive area (T cells) was numerically lower in grafts transplanted into miR155−/− recipients, yet this didn't reach statistical significance (p = 0.3) (Figure 5E). Taken together, these results indicate that higher leukocytic expression levels of miR‐155 are associated with increasing numbers of infiltrating leukocytes in general and macrophages more specifically.

AntagomiR‐155 treatment inhibits acute cardiac allograft rejection

We inhibited miR‐155 in vivo using three injections of antagomiR‐155 (anti‐155) targeting mature mmu‐miR‐155 and compared its beneficial effects with antagomiR‐scrambled oligonucleotide (anti‐scram) injected recipients. Intravenous injections were administered 1 day before scheduled transplantation and on POD 1 and 3. These time points were arbitrarily chosen to simulate perioperative treatment as could be applied in patients. Mice were sacrificed on POD 5 (Figure 6A).
Figure 6

ACR in anti‐155 and anti‐scram injected mice following HTX. (A) Experimental setup of anti‐155 and anti‐scram injected mice sacrificed on a fixed time point (5d) following HTX. (B) Normalized expression levels of miR‐155 in (un) transplanted hearts following anti‐155 and anti–scram injection. (C) Normalized splenic expression levels of miR‐155 in (un) transplanted mice following injection with anti‐155 and anti‐scram. (D) BALB/cJ grafts weighed significantly less when transplanted into anti‐155 injected mice. The number of infiltrating CD45 (E) and Mac3‐positive (F) cells was significantly less in BALB/cJ hearts transplanted into anti‐155 versus anti‐scram injected mice. (G) The number of CD3‐positive cells was not significantly different between treatments. (H) Following HTX, splenic SPI1 decreased in anti‐scram mice, whereas it remained elevated in anti‐155 mice. Splenic SPI1 was significantly higher in transplanted anti‐155 versus anti‐scram injected mice (*p < 0.05; **p < 0.01). ACR, acute cellular rejection; HTX, heart transplantation.

ACR in anti‐155 and anti‐scram injected mice following HTX. (A) Experimental setup of anti‐155 and anti‐scram injected mice sacrificed on a fixed time point (5d) following HTX. (B) Normalized expression levels of miR‐155 in (un) transplanted hearts following anti‐155 and anti–scram injection. (C) Normalized splenic expression levels of miR‐155 in (un) transplanted mice following injection with anti‐155 and anti‐scram. (D) BALB/cJ grafts weighed significantly less when transplanted into anti‐155 injected mice. The number of infiltrating CD45 (E) and Mac3‐positive (F) cells was significantly less in BALB/cJ hearts transplanted into anti‐155 versus anti‐scram injected mice. (G) The number of CD3‐positive cells was not significantly different between treatments. (H) Following HTX, splenic SPI1 decreased in anti‐scram mice, whereas it remained elevated in anti‐155 mice. Splenic SPI1 was significantly higher in transplanted anti‐155 versus anti‐scram injected mice (*p < 0.05; **p < 0.01). ACR, acute cellular rejection; HTX, heart transplantation. In unoperated mice, cardiac miR‐155 levels, known to be very low in healthy hearts, did not decrease significantly following anti‐155 treatment; however, splenic levels decreased significantly (Figure 5C). In transplanted animals, miR‐155 levels decreased significantly both in cardiac allografts and spleens upon anti‐155 treatment (Figure 6B and C). Biometrical analysis revealed significantly lower graft weights in anti‐155 treated mice on POD 5 (anti‐scram: 4.8 mg/g ± 0.1; anti‐155: 4.0 mg/g ± 0.2; p = 0.02) (Figure 6D). More importantly, CD45 staining area decreased significantly following anti‐155 treatment, as did the number of Mac3‐positive cells infiltrating the graft (Figure 6E and F). Though numerically lower, CD3‐positive area in the grafts was not significantly lower in anti‐155 treated animals (Figure 5G). In conclusion, pharmacological miR‐155 inhibition was feasible and significantly decreased the number of infiltrating leukocytes and macrophages, which is in line with our results obtained in miR‐155−/− mice.

MiRs as a nonorgan specific signature of ACR?

The 10 miRs identified as a common signature of ACR following HTX were also quantified in human renal allograft biopsies. We compared rejecting and nonrejecting samples. Comparable with our observations in cardiac allografts, the subset of five miRs previously implicated in inflammatory processes, miR‐142‐3p, miR‐142‐5p, miR‐146b, miR‐155, miR‐223, was able to distinguish between rejecting and nonrejecting samples (Figure 7).
Figure 7

Quantitative RT‐PCR validation of miR array results in human rejecting and nonrejecting renal allografts (p < 0.05)

Quantitative RT‐PCR validation of miR array results in human rejecting and nonrejecting renal allografts (p < 0.05) In conclusion, we corroborated our miR hallmark of acute cellular cardiac allograft rejection in acute cellular renal allograft rejection, hinting at a nonspecies and nonorgan specific miR signature of ACR.

Protein levels of miR‐155 targets involved in innate immunity decrease during ACR

We quantified the expression of SPI 1 early (Allo3d) and late (Allo7d) after transplantation. SPI 1 is a distinct leukocytic miR‐155 target implicated in the repression of IL‐6 signaling 24 and acting as an inhibitor of dendritic cell pathogen binding and antigen presentation of dendritic cells to T cells 25. Splenic expression levels of SPI1 decreased in a time‐dependent way (Figure 2K). Importantly, RNA levels of SPI1 increased, providing evidence of miR‐155's repressive activity with regard to the translation of SPI1 mRNA to protein (Figure 2L). Therefore, during ACR, miR‐155 may promote antigen presentation by dendritic cells to T cells and stimulate IL‐6 signaling by repressing SPI 1, a pivotal inhibitor of these pathways. To corroborate the decreased influx of inflammatory cells into the allograft following repeated intravenous anti‐155 administration and elaborate on the implication of miR‐155 in ACR, we compared SPI1 expression in anti‐scram and −155 injected animals. Without transplantation, splenic SPI1 levels were similar in both treatment groups (p = 0.4). Following cardiac transplantation, splenic SPI1 decreased in anti‐scram injected animals, whereas in anti‐155 injected mice, SPI1 protein levels remained elevated and significantly higher compared with anti‐scram injected mice (p = 0.03; Figure 6H). These data provide evidence for a role of SPI1 and miR‐155 during ACR and at the same time demonstrate the feasibility of interfering with miR expression to modulate critical signaling pathways involved in allograft rejection.

Discussion

This study is the first to unravel both miR and mRNA expression profiles in normal hearts and ACR in humans and mice, leading to a miR signature that discriminates rejected from nonrejected grafts. We identify miR‐155 as a key contributor to ACR: in humans, miR‐155 is the most upregulated miR following rejection and, using an experimental model, we showed that reduced levels of miR‐155 attenuated intragraft inflammation and delayed rejection. Apart from establishing miR‐155 as a candidate target for novel therapeutics, we provide evidence for a new paradigm to combat rejection, namely by antimiR 18 or antagomiR‐based therapeutics. MiR‐155 has an established role in inflammation and immunity 26, 27. Macrophage activation and dendritic cell maturation require its upregulation 16, 25, 28. Furthermore, miR‐155 functions at the interface of innate and adaptive immunity as miR‐155 deficient dendritic cells are unable to mount appropriate B and T cell responses 14, 16, 23. These immunological mechanisms all have a pivotal role in ACR. In our study, genomic absence and pharmacological inhibition of miR‐155 delayed graft failure by reducing influx of leukocytes (CD45‐positive cells) in general and macrophages (Mac3‐staining cells) more specifically. In another recent report, antagomiR‐155 attenuated allograft rejection through the inhibition of dendritic cell maturation, yet this study did not address miR expression in human transplant samples 14. Taken together, these two studies nonetheless provide evidence for a causal relationship between miR‐155 and ACR and the potential of miR‐155 targeting to prevent organ rejection. Apart from miR‐155, our miR expression arrays and validation qPCRs revealed additional miRs with comparable changes in human and murine inflamed cardiac grafts: miR‐149‐5p (downregulated), miR‐21, miR‐146a, miR‐146b, miR‐142‐3p, miR‐142‐5p, miR‐222, miR‐223, and miR494 (upregulated). MiR‐142‐3p, miR‐142‐5p, miR‐146a, miR‐146b, and miR‐223 are notorious for their implication in inflammation and the innate immune system 27, 29, 30. More specifically, upregulation of miR‐142‐3p, a hematopoietic cell‐specific miR, may modulate inflammatory responses to promote graft tolerance, whereas upregulation of miR‐142‐5p was previously linked to graft rejection 29, 31. Increased miR‐21 and miR‐222 expression have been implicated in fibrosis 32, 33, providing a hypothetical link between ACR and myocardial interstitial fibrosis, a specific manifestation of chronic rejection 1, 34. The implication of miR‐494 and miR‐149‐5p in ACR deserves further study. Importantly, TNFα, which increases upon cellular rejection 34, upregulates miR‐494 35 and downregulates miR‐149‐5p 36, in line with our findings. The upregulation of miR‐494 might be beneficial as it prevents ischemia‐reperfusion injury 37. The downregulation of miR‐149‐5p may contribute to the induction of genes downstream of the proinflammatory cytokine IL‐6 38, which we also observed in our mRNA arrays. Interestingly, others have shown identical miRs as in our study to be changed during ACR of transplanted organs 7, 9, 10, 11, 12, 13. We, therefore, hypothesize that a generic change in miR expression occurs during ACR, irrespective of the organ transplanted. This generic change, occurring despite contemporary immunosuppressive regimens, merits further study as it hints at common pro‐inflammatory pathways in these diverse manifestations of ACR. Admittedly, we restricted our analysis to tissular expression of mRNAs and miRs, whereas others studied their serological expression 7, 39. Furthermore, we only compared rejecting and nonrejecting human samples collected during the first year following HTX, the typical “vulnerable phase” for ACR, whereas others addressed miR expression around 1 year after transplantation 7. These two caveats may explain why miR expression profiles do not completely overlap between studies. Our mRNA arrays also revealed interesting new targets to prevent and treat ACR. The comparison of significantly dysregulated transcripts in transplanted mice and patients manifesting rejection yielded 70 common dysregulated mRNAs, with 33 of these 70 coding for proteins involved in the IL‐6 pathway. These findings are in line with previous reports on the role of IL‐6 in cardiac rejection 40, 41 and urge for further efforts to determine the true potential of IL‐6 targeting to prevent ACR. Importantly, the pro‐inflammatory effects of miR‐155 during ACR implicate that it targets mRNAs coding for anti‐inflammatory proteins. Intriguingly, our mRNA arrays showed that the 33 dysregulated transcripts involved in the IL‐6 pathway also function downstream of SPI1 18, 25 an established target of miR‐155. Whereas the anti‐inflammatory transcription factor SPI1 decreased during ACR, miR‐155 inhibition de‐repressed SPI1. We therefore postulate that miR‐155 inhibition attenuates allograft rejection partially by de‐repression of SPI1, apart from its recently unraveled inhibitory effect on dendritic cell maturation through SOCS1 upregulation 14. Our study has some limitations. In our derivation cohort, we decided to use untransplanted BALB/cJ hearts as controls rather than BALB/cJ our C57Bl/6J syngeneic grafts. Earlier research indicated that during normal development 42 and following physiological 43 and pathophysiological 44 stimuli applied to the heart, strain‐specific differences are undeniable. In addition, we hypothesize that our comparison with human samples intentionally collected beyond 6 weeks after transplant surgery cancels out the inflammatory signature associated with the surgical procedure itself. In retrospect, grafting a BALB/cJ heart into C57Bl/6J Rag −/− common gamma‐chain −/− mice would make the optimal control, as this would cancel out the nonimmune inflammatory signature. Furthermore, although we used a sample set of 49 EMBs for validation qPCRs, the scarcity of EMBs manifesting severe rejection in our human cohort obliged us to reuse samples of the derivation cohort in the validation cohort. Future research should validate our markers in an independent set of samples. In conclusion, ACR induced marked miR and mRNA expression changes in human and murine grafts. MiR‐155 was the most upregulated miR during ACR in humans and contributes to adverse cardiac inflammation, as evidenced by prolonged graft survival and diminished leukocyte influx in mice with a targeted genomic deletion of miR‐155 and following pharmacological inhibition of miR‐155 using antagomiRs. Additional comparison of dysregulated mRNAs in human and murine grafts during ACR‐revealed marked induction of genes downstream of the proinflammatory cytokine IL‐6 and SPI1, a conserved target of miR‐155. As such, this study identifies miRs in general, miR‐155 in particular, as well as the IL‐6 signaling pathway, as promising targets for the treatment of ACR.

Disclosure

The authors of this manuscript have no conflicts of interest to disclose as described by the American Journal of Transplantation. Figure S1: Quantitative RT‐PCR validation of miR array results in murine cardiac allografts. Click here for additional data file. Figure S2: Quantitative RT‐PCR validation of miR array results in human cardiac allografts. Click here for additional data file. Figure S3: Acute cellular rejection in miR‐155 +/+ and miR‐155 −/− mice following HTX. Click here for additional data file. Table S1: Differentially regulated human miRs in ISHLT3 versus Control. Click here for additional data file. Table S2: Ingenuity Pathway Analysis of differentially regulated human miRs in ISHLT3 versus Control EMBs. Click here for additional data file. Table S3: Differentially regulated mouse miRs in Allo7d versus Control. Click here for additional data file. Table S4: Ingenuity Pathway Analysis of differentially regulated mouse miRs in Allo7d versus Control. Click here for additional data file. Table S5: Patient characteristics in validation cohort. Click here for additional data file. Table S6: Differentially regulated human and mouse mRNA transcripts from the micro array studies. Click here for additional data file. Table S7: Pathway analysis of differentially regulated human mRNAs in ISHLT 3 versus Control EMBs. Click here for additional data file. Table S8: Pathway analysis of differentially regulated mouse mRNAs in ACR. Click here for additional data file.
  43 in total

1.  Upregulation of miR-142-3p in peripheral blood mononuclear cells of operationally tolerant patients with a renal transplant.

Authors:  Richard Danger; Annaïck Pallier; Magali Giral; Marc Martínez-Llordella; Juan José Lozano; Nicolas Degauque; Alberto Sanchez-Fueyo; Jean-Paul Soulillou; Sophie Brouard
Journal:  J Am Soc Nephrol       Date:  2012-01-26       Impact factor: 10.121

2.  Revision of the 1990 working formulation for the standardization of nomenclature in the diagnosis of heart rejection.

Authors:  Susan Stewart; Gayle L Winters; Michael C Fishbein; Henry D Tazelaar; Jon Kobashigawa; Jacki Abrams; Claus B Andersen; Annalisa Angelini; Gerald J Berry; Margaret M Burke; Anthony J Demetris; Elizabeth Hammond; Silviu Itescu; Charles C Marboe; Bruce McManus; Elaine F Reed; Nancy L Reinsmoen; E Rene Rodriguez; Alan G Rose; Marlene Rose; Nicole Suciu-Focia; Adriana Zeevi; Margaret E Billingham
Journal:  J Heart Lung Transplant       Date:  2005-06-20       Impact factor: 10.247

3.  A technique of cervical heterotopic heart transplantation in mice.

Authors:  Z H Chen
Journal:  Transplantation       Date:  1991-12       Impact factor: 4.939

4.  The Registry of the International Society for Heart and Lung Transplantation: 29th official adult heart transplant report--2012.

Authors:  Josef Stehlik; Leah B Edwards; Anna Y Kucheryavaya; Christian Benden; Jason D Christie; Anne I Dipchand; Fabienne Dobbels; Richard Kirk; Axel O Rahmel; Marshall I Hertz
Journal:  J Heart Lung Transplant       Date:  2012-10       Impact factor: 10.247

5.  MicroRNA-155 modulates the pathogen binding ability of dendritic cells (DCs) by down-regulation of DC-specific intercellular adhesion molecule-3 grabbing non-integrin (DC-SIGN).

Authors:  Rocio T Martinez-Nunez; Fethi Louafi; Peter S Friedmann; Tilman Sanchez-Elsner
Journal:  J Biol Chem       Date:  2009-04-22       Impact factor: 5.157

6.  miRNA profiling discriminates types of rejection and injury in human renal allografts.

Authors:  Julia Wilflingseder; Heinz Regele; Paul Perco; Alexander Kainz; Afschin Soleiman; Ferdinand Mühlbacher; Bernd Mayer; Rainer Oberbauer
Journal:  Transplantation       Date:  2013-03-27       Impact factor: 4.939

Review 7.  Non-invasive approaches for the diagnosis of acute cardiac allograft rejection.

Authors:  Christopher A Miller; James E Fildes; Simon G Ray; Helen Doran; Nizar Yonan; Simon G Williams; Matthias Schmitt
Journal:  Heart       Date:  2012-12-20       Impact factor: 5.994

Review 8.  The emerging role of MIR-146A in the control of hematopoiesis, immune function and cancer.

Authors:  Catherine Labbaye; Ugo Testa
Journal:  J Hematol Oncol       Date:  2012-03-27       Impact factor: 17.388

9.  Distinct microRNA profiles are associated with the severity of hepatitis C virus recurrence and acute cellular rejection after liver transplantation.

Authors:  Deepak Joshi; Siamak Salehi; Helen Brereton; Matthew Arno; Alberto Quaglia; Nigel Heaton; John O'Grady; Kosh Agarwal; Varuna Aluvihare
Journal:  Liver Transpl       Date:  2013-03-22       Impact factor: 5.799

10.  Expression of miR-142-5p in peripheral blood mononuclear cells from renal transplant patients with chronic antibody-mediated rejection.

Authors:  Richard Danger; Chloé Paul; Magali Giral; Amélie Lavault; Yohann Foucher; Nicolas Degauque; Annaïck Pallier; Maxim Durand; Stéphanie Castagnet; Jean-Paul Duong Van Huyen; Michel Delahousse; Karine Renaudin; Jean-Paul Soulillou; Sophie Brouard
Journal:  PLoS One       Date:  2013-04-05       Impact factor: 3.240

View more
  14 in total

1.  High Dimensional Renal Profiling: Towards a Better Understanding or Renal Transplant Immune Suppression.

Authors:  Cyd M Castro-Rojas; Rita R Alloway; E Steve Woodle; David A Hildeman
Journal:  Curr Transplant Rep       Date:  2019-01-14

2.  Regulation of Endothelial-to-Mesenchymal Transition by MicroRNAs in Chronic Allograft Dysfunction.

Authors:  Emily K Glover; Nina Jordan; Neil S Sheerin; Simi Ali
Journal:  Transplantation       Date:  2019-04       Impact factor: 4.939

Review 3.  MicroRNAs in kidney injury and disease.

Authors:  Nassim Mahtal; Olivia Lenoir; Claire Tinel; Dany Anglicheau; Pierre-Louis Tharaux
Journal:  Nat Rev Nephrol       Date:  2022-08-16       Impact factor: 42.439

Review 4.  Opportunities for microRNAs in the Crowded Field of Cardiovascular Biomarkers.

Authors:  Perry V Halushka; Andrew J Goodwin; Marc K Halushka
Journal:  Annu Rev Pathol       Date:  2018-10-17       Impact factor: 23.472

Review 5.  MicroRNAs in Heart Failure, Cardiac Transplantation, and Myocardial Recovery: Biomarkers with Therapeutic Potential.

Authors:  Palak Shah; Michael R Bristow; J David Port
Journal:  Curr Heart Fail Rep       Date:  2017-12

6.  Association of Serum MiR-142-3p and MiR-101-3p Levels with Acute Cellular Rejection after Heart Transplantation.

Authors:  Ihdina Sukma Dewi; Zsuzsanna Hollander; Karen K Lam; Janet-Wilson McManus; Scott J Tebbutt; Raymond T Ng; Paul A Keown; Robert W McMaster; Bruce M McManus; Olof Gidlöf; Jenny Öhman
Journal:  PLoS One       Date:  2017-01-26       Impact factor: 3.240

7.  Acute Rejection After Kidney Transplantation Associates With Circulating MicroRNAs and Vascular Injury.

Authors:  Roel Bijkerk; Barend W Florijn; Meriem Khairoun; Jacques M G J Duijs; Gurbey Ocak; Aiko P J de Vries; Alexander F Schaapherder; Marko J K Mallat; Johan W de Fijter; Ton J Rabelink; Anton Jan van Zonneveld; Marlies E J Reinders
Journal:  Transplant Direct       Date:  2017-06-19

8.  Functional Screening Identifies MicroRNAs as Multi-Cellular Regulators of Heart Failure.

Authors:  Robin Verjans; Wouter J A Derks; Kerstin Korn; Birte Sönnichsen; Rick E W van Leeuwen; Blanche Schroen; Marc van Bilsen; Stephane Heymans
Journal:  Sci Rep       Date:  2019-04-15       Impact factor: 4.379

9.  miR-149-5p Inhibits Vascular Smooth Muscle Cells Proliferation, Invasion, and Migration by Targeting Histone Deacetylase 4 (HDAC4).

Authors:  Boya Zhang; Yang Dong; Ming Liu; Lei Yang; Zhuo Zhao
Journal:  Med Sci Monit       Date:  2019-10-09

10.  A Combined microRNA and Chemokine Profile in Urine to Identify Rejection After Kidney Transplantation.

Authors:  Els M Gielis; Jacqueline D H Anholts; Els van Beelen; Geert W Haasnoot; Hans W De Fijter; Ingeborg Bajema; Sebastiaan Heidt; Mathijs van de Vrie; Luuk B Hilbrands; Marko J K Mallat; Kristien J Ledeganck; Frans H J Claas; Michael Eikmans
Journal:  Transplant Direct       Date:  2021-06-10
View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.