Next Article in Journal
The Effect of Ca2+, Lobe-Specificity, and CaMKII on CaM Binding to NaV1.1
Previous Article in Journal
Evaluation of the Genetic Diversity and Differentiation of Black Locust (Robinia pseudoacacia L.) Based on Genomic and Expressed Sequence Tag-Simple Sequence Repeats
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 9
10.3390/ijms19092493
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exosomal miRNA Signatures for Late-Onset Acute Graft-Versus-Host Disease in Allogenic Hematopoietic Stem Cell Transplantation

by
Seiichiro Yoshizawa
1,*,
Tomohiro Umezu
1,2,
Yuu Saitoh
1,
Moritaka Gotoh
1,
Daigo Akahane
1,
Chiaki Kobayashi
1,
Junko H. Ohyashiki
2 and
Kazuma Ohyashiki
1,2
1
Department of Hematology, Tokyo Medical University, Tokyo 160-0023, Japan
2
Department of Advanced Cellular Therapy, Tokyo Medical University, Tokyo 160-0023, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(9), 2493; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092493
Submission received: 27 July 2018 / Revised: 15 August 2018 / Accepted: 21 August 2018 / Published: 23 August 2018
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)
Browse Figures
Versions Notes

Abstract

:
Recent studies have demonstrated that exosomal microRNAs (miRNAs) have the potential of facilitating molecular diagnosis. Currently, little is known about the underlying mechanism behind late-onset acute graft-versus-host disease (LA GVHD). Identifying differentially expressed miRNAs in exosomes should be useful for understanding the role of miRNAs in this disease. This study was established to investigate the relevance of miRNAs in exosomes derived from patients developing LA GVHD after allogeneic hematopoietic stem cell transplantation (HSCT). Plasma samples were collected from patients with LA GVHD (n = 5), non-GVHD (n = 5), and controls (n = 8) for exosomal miRNA expression profiling using a TaqMan low-density array; the results were validated by quantitative reverse transcription polymerase chain reaction (RT-PCR). We analyzed exosomal miRNAs differentially expressed among these three groups. MirTarBase was employed to predict potential target genes of the miRNAs specific for LA GVHD. We detected 55 miRNAs that were differentially expressed with a significant change >2.0-fold between LA GVHD and non-GVHD. Of these, we selected the 10 miRNAs (miR-423-5p, miR-19a, miR-142-3p, miR-128, miR-193b, miR-30c, miR-193a, miR-191, miR-125b, and miR-574-3p) with the most significant differential expression. Using quantitative RT-PCR, we further identified that miR-128 was significantly upregulated at the onset of LA GVHD compared with that in normal controls and is a promising diagnostic marker of LA GVHD, with an area under the curve (AUC) value of 0.975. MirTarBase analysis revealed that the predicted target genes of miR-128 are involved in the immune system and inflammation. Increased expression of miR-128 may serve as a novel, noninvasive biomarker for early LA GVHD diagnosis.

Graphical Abstract

1. Introduction

The occurrence of severe graft-versus-host disease (GVHD) after allogeneic hematopoietic stem cell transplantation (HSCT) leads to a poor prognosis. Recently, with this expansion in the utilization of allogeneic HSCT, such as the introduction of reduced-intensity conditioning and donor lymphocyte infusions, late-onset acute and chronic GVHD has been reclassified based on the exact clinical manifestation instead of the time of onset [1]. According to the new classification using National Institutes of Health (NIH) consensus criteria established in 2005, patients with de novo acute GVHD occurring beyond day 100 post-transplantation are considered to have late-onset acute GVHD (LA GVHD) [1]. In accordance with the NIH definition, 20–40% of chronic GVHD cases are reconsidered to be LA GVHD. A few reports have shown that LA GVHD is also associated with greater nonrelapse mortality after allogeneic HSCT than non-GVHD [2,3]. However, there is still an insufficient understanding of the association between LA GVHD and a lower incidence of relapse caused by the induction of a beneficial graft-versus-tumor (GVT) effect. To improve long-term survival, intensive efforts have been made to prevent GVHD without losing the GVT effect. Thus, several molecular biomarkers for diagnosis and prediction of the severity and prognosis of GVHD have been explored [4,5,6,7].
MicroRNAs (miRNAs) are small noncoding RNA molecules that regulate post-transcriptional gene expression by degrading or suppressing their target mRNA [8]. Recently, it has been widely recognized that some circulating miRNAs are included in the content of endosome-derived extracellular vesicles including exosomes. Exosomes are secreted by many types of cells into various biological fluids, such as serum and plasma [9,10]. The demonstration that exosomes contain not only proteins, RNAs, and lipids, but also functional miRNAs indicated that miRNAs could participate in intercellular communication and elicit immune responses by their transport between cells via exosomes [11,12]. Therefore, analysis of the alterations of exosomal miRNAs in the context of various diseases has the potential to facilitate molecular diagnosis [13,14].
In recent studies, miRNA profiling identified a GVHD-specific miRNA signature that can control the regulation of classic acute GVHD in allogeneic immunity. In particular, miR-155, known as an immune-related microRNA, and miR-15a/16 and miR-17-92a were shown to be closely involved in the pathophysiology behind the development of classic acute GVHD in a xenogeneic GVHD murine model [15,16,17]. Moreover, recent studies have demonstrated that specific plasma miRNA signatures, such as miR-423, miR-199a-3p, miR-93, miR-377, and miR-586, could serve as independent biomarkers for the prediction, diagnosis, and prognosis of hematologic diseases and classic acute GVHD [18,19]. Despite the fact that the roles of miRNAs in classic acute GVHD after allogeneic HSCT have been extensively studied over the last few years, the relevance and expression profile of miRNAs in LA GVHD have not yet been elucidated. Because a number of previous studies investigated the role of miRNAs as potential predictive biomarkers in acute GVHD, we aimed to determine the clinically relevant exosomal miRNA profile in patients developing LA GVHD after allogeneic HSCT.

2. Results

2.1. Patient Characteristics

A summary of patients’ clinical characteristics is presented in Table 1. The median age of all patients was 58.5 years (range, 26–70 years) and the median onset of LA GVHD was 165 days (range, 117–238 days) after allogeneic HSCT. GVHD prophylaxis included mainly calcineurin inhibitor with or without short-term methotrexate. Overall, four of the five of patients with LA GVHD developed severe GVHD (grades III–IV) and no patients achieved remission despite second-line therapy. Four patients died from GVHD deterioration, while the remaining one patient died from post-transplantation lymphoproliferative disorder.

2.2. miRNA Expression Profiling at LA GVHD Diagnosis

The procedure for selecting candidate miRNAs is illustrated in Figure 1. To identify candidate miRNAs with altered expression in LA GVHD patients, miRNA expression profiling was performed on exosome fractions taken from plasma samples of 10 allogeneic HSCT patients (5 LA GVHD and 5 non-GVHD) and 8 normal subjects. We then compared exosomal miRNA expression in LA GVHD patients with that in non-GVHD patients using a TaqMan miRNA Array (NCBI, Gene Expression Omnibus (GEO); GSE110435). The heat map in Figure 2 demonstrates that the miRNA signature in LA GVHD group tended to show high expression. Among these three groups, a total of 163 differentially expressed miRNAs were identified using GeneSpring software. We confined our analysis to miRNAs for which the expression level exhibited a fold change (FC) of >2.0 in exosomes of the LA GVHD group compared with that in the non-GVHD group. Interestingly, for all of the 55 selected miRNAs, the expression was upregulated in the LA GVHD group. We then evaluated the 10 miRNAs with the highest FC (miR-423-5p, miR-19a, miR-142-3p, miR-128, miR-193b, miR-30c, miR-193a, miR-191, miR-125b, and miR-574-3p), as shown in Table 2. Of these miRNAs, miR-423-5p, miR-19a, and miR-142, which were most strikingly differentially expressed, were previously reported to play important roles in acute GVHD [18,20,21].

2.3. Comparison of Exosomal miRNAs Differentially Expressed among LA GVHD, Non-GVHD, and Healthy Volunteers

To validate the findings on the LA GVHD-specific exosomal miRNAs, we assessed the 10 exosomal miRNAs that were most highly differentially expressed among the three groups. We adopted the screening method for selecting LA GVHD-specific exosomal miRNAs as follows: (1) significant difference in the expression of miRNAs is observed between LA GVHD and healthy volunteers (normal control); or (2), no significant difference in the expression of miRNAs is observed between the non-GVHD and healthy volunteers (normal control). We firstly compared samples obtained from patients at the diagnosis of LA GVHD with those from healthy volunteers (normal control) and found that three miRNAs (miR-423-5p, miR-128, and miR-125b) were significantly overexpressed in the exosomes of those with LA GVHD (p < 0.05) when compared with the levels in healthy subjects (Figure 3A). Meanwhile, none of the remaining seven miRNAs showed significant differences in expression between LA GVHD and healthy volunteers (Figure 3B). Subsequently, we compared samples obtained from non-GVHD patients with those from healthy volunteers in Figure 3A. No significant difference in the expression of miR-128 and miR-125b was observed between the non-GVHD and healthy volunteers, while miR-423-5p showed downregulation in the non-GVHD group when compared with healthy volunteers. Therefore, two miRNAs, miR-128 and miR-125b, were chosen as candidate miRNAs specific for LA GVHD; For further analysis, we focused on the clinical implications of these two dysregulated miRNAs in LA GVHD.

2.4. Quantification of Individual miRNA by Real-Time Quantitive RT-PCR

To confirm the specific deregulated miRNAs change obtained from the TaqMan low-density array (TLDA), we quantified miR-128 and miR-125b expression by quantitative real-time PCR (qRT-PCR). There were no significant differences in miR128 (Figure 4A) and miR-125b (Figure 4B) expression between the normal control and non-GVHD group (p = 0.3842 and 0.8336, respectively). In contrast, miR-128 and miR-125b were significantly increased in LA GVHD patients (LA GVHD vs. non-GVHD, p = 0.0027 and 0.0289; LA GVHD vs. NC, p = 0.0001 and 0.0104, respectively) (Figure 4A,B). To estimate cut-off values for the two miRNAs (miR-128, miR-125b) in exosomes for distinguishing between LA GVHD and non-GVHD groups, we performed receiver-operating characteristic (ROC) curve analysis. The ROC curve for miR-128 showed extremely high diagnostic accuracy, with an area under the ROC curve of 0.975 and diagnostic sensitivity of 100% [95% confidence interval (CI), 47.82–100%] and specificity of 87.5% (95% CI, 47.35–99.68%) (Figure 5A). Similarly, the ROC curve for miR-125b showed an AUC of 0.975 and diagnostic sensitivity and specificity of 80% (95% CI, 28.36–99.49%) and 87.5% (95% CI, 47.35–99.68%), respectively (Figure 5B). The cut-off values were thus determined to be 3.183 for miR-128 and 2.29 for miR-125b, which were used in the following analysis. Our data suggest that altered expression of miRNAs may serve as a diagnostic tool for LA GVHD in HSCT practice.

2.5. Sequential Analysis of Exosomal miRNAs Post-Allogenic HSCT

To further explore whether exosomal miRNAs can be used as predictive biomarkers of LA GVHD, we next measured miR-125b and miR128 expressions sequentially in exosomes of available serum samples taken from two transplanted patients (patient number PN 3 and PN 4 in Table 1). As illustrated in Figure 6, we observed that the level of miR-128 was significantly upregulated approximately a week or two weeks before the onset of LA GVHD in PN 3 and PN 4 (Figure 6A,B). When we used the cut-off value (3.183) obtained from the ROC analysis, the levels of miR-128 expression on days +25, +46, and +70 of PN 3 were under this value but exceeded the cut-off value at day +105 in the prediagnosis of LA GVHD (Figure 6A). In addition, the levels of miR-128 expression on days +11 and +81 of PN 4 were under that value but exceeded the cut-off value at day +111 in the prediagnosis of LA GVHD (Figure 6B). On the other hand, the miR-125b expression level was upregulated immediately after HSCT, and elevated miR-125b expression was unable to predict LA GVHD (Figure 6C,D). This suggests that increased miR-128 expression in patients after HSCT may be a predictive marker prior to the onset of LA GVHD. Finally, the validity of miR-128 as a biomarker was confirmed by using it to discriminate LA GVHD patients from non-GVHD cases, although we should consider the bias of stem cell source.

2.6. Prediction of the Target Genes of Candidate miRNAs

This study indicated that upregulated miR-128 is crucial as an LA GVHD-specific miRNA. The clinical relevance of miRNAs involved in the development of LA GVHD was not evident. Thus, MirTarBase was used to find predicted target genes of miR-128 in order to determine the biological significance of this miRNA. More than 100 target genes were extracted by MirTarBase, and the target genes for which strong evidence existed in the results of the search are summarized in Table 3. Accordingly, 12 genes, including B-cell specific Moloney murine leukemia virus integration region 1 (BMI1), F-box, and tryptophan-aspartic acid (WD) repeat domain containing 7 (FBXW7), were experimentally validated by a luciferase reporter assay in the literature [22,23,24,25,26,27].

3. Discussion

Recently, the definition of GVHD was revised in the NIH consensus criteria, with the new category of “late-onset” acute (LA) GVHD being established. Given the short history of this category, the pathophysiology of LA GVHD is still poorly understood. Although the diagnosis of GVHD is judged based on clinical manifestations and pathological findings of the tissue sections from the intestines or liver, it is difficult to provide an accurate pathological diagnosis of GVHD. For this reason, pathological findings have greatly affected various factors of radiation damage, drug injury, and viral infection, presenting a diagnostic dilemma. Therefore, there is an urgent need of exploring reliable, noninvasive biomarkers of GVHD. In the context of GVHD after allogeneic HSCT, biomarker research has progressed through a variety of techniques, such as immunocompetent cells (CD4 helper T cells, regulatory T cells), plasma (soluble IL-2R), and inflammatory cytokines (IFN-γ, TNF-α). For example, Paczesny et al. identified that a panel of four biomarkers (IL-2Rα, TNFR1, IL-8, HGF) in serum was useful for the diagnosis of acute GVHD [4]. This group also found that the suppression of tumorigenicity 2 (ST2) could distinguish treatment-resistant GVHD patients, and that ST2 levels were associated with nonrelapse mortality [6]. Moreover, a specific plasma miRNA signature for classic acute GVHD was recently determined by Xiao et al. [18]. To date, the identification of altered miRNAs has been thought to contribute to the HSCT-related pathogenesis of acute GVHD; however, the mechanism underlying LA GVHD remains to be elucidated. The above-mentioned results provided us with a potential tool as a new liquid biopsy technique.
Because a number of previous studies investigated the role of exosomes and exosomal miRNAs as potential predictive biomarkers in acute GVHD [20,28], we aimed to elucidate the possible role of miRNAs specifically in LA GVHD patients. Exosomes containing miRNAs reflect the behavior of certain types of cells, so a survey of changes of miRNAs in exosomes should promote our understanding of clinical conditions. In this context, technical advances in exosome separation and characterization have also recently been achieved, enabling exosomal miRNAs to be easily detected by high-throughput techniques [29]. Additionally, Kordelas et al. reported an attractive treatment using mesenchymal stem cell (MSC)-derived exosomes for therapy-refractory GVHD [30]. We also previously reported that many miRNAs have been identified as noninvasive, diagnostic, prognostic, and predictive markers for hematologic malignancy [31,32,33]. In the current study, we explored the LA GVHD-specific exosomal miRNA signature for LA GVHD compared with that in non-GVHD. We also performed bioinformatic analysis to investigate the potential molecular mechanisms behind LA GVHD and identify molecular target genes. Our results of validation by qRT-PCR showed that miRNA-128 significantly increased in LA GVHD patients. The miRNA miR-128 is well known as a brain-enriched miRNA, which plays a crucial role in the development of the nervous system [34]. The miRNA miR-128 has also been previously demonstrated to be dysregulated in tissues and blood samples in several malignant tumor patients and to be correlated with tumor progression [35,36,37]. In the current study, it was impossible to conduct a functional analysis of miR-128. Hence, using MirTarBase, we instead searched for potential target genes to investigate the mechanism of miR-128 (Table 3).
We perfomed ROC analysis to assess the diagnostic capacity of exosomal miRNAs. Exosomal miR-128 exhibited its potential in discriminating LA GVHD from non-GVHD groups. When the expression level of candidate miRNAs was analyzed sequentially in a subset of LA GVHD patients during allogeneic HSCT, miR-128 was shown to increase over the cut-off value earlier than the development of the disease. However, our results must be interpreted with caution, because this study was limited by the relatively small number of samples from LA GVHD patients and absence of a validation cohort of any sort. After all, molecular mechanisms on the dysregulation of exosomal miRNAs in patients with LA GVHD could not be fully provided. Also, candidate miRNAs were then examined longitudinally in only two patients over time. Thus, these results need to be further studied with large cohort.
The results of the top 10 miRNAs obtained by TLDA included several characterized miRNAs (miR-423-5p, miR-19a, and miR-142) known to be dysregulated in acute GVHD [18,20,21]. In particular, miR-423-5p, which is involved in the immune response of acute GVHD incidence, was elevated up to 49-fold in exosomes of LA GVHD patients in our study. Crossland et al. reported differential expression of miR-423 in exosomes at day 14, after HSCT had shown potential as a predictive biomarker for the occurrence of acute GVHD [20]. Therefore, the pathology of LA GVHD may resemble that of acute GVHD caused by donor T-cell alloreactivity in light of miRNAs. However, the biology of GVHD is complex. In other words, miRNAs also regulate various molecular targets, including normal development, differentiation, and maturation of hematopoietic cells in the immune system [38]. The difference of miRNA expression may influence the process of immune reconstitution after allogeneic HSCT. Moreover, the changes of immune-related miRNAs may be due to the dosage of immunosuppression. The dysregulated miRNAs when LA GVHD occurs may be repressed by immunosuppressive regimen. However, it is difficult to evaluate whether the level of this miRNA returns to baseline upon therapeutic intervention because of the uncontrollable nature of LA GVHD. It is a challenge for future study to clarify the association between miRNA expression and response to therapy. Although it is difficult to define miRNAs specific to LA GVHD in very complex medical cases, and in the small number of transplanted patients in the current study, exosomal miR-128 may be related to the immunoregulation of LA GVHD.
In summary, this study demonstrated that the levels of several exosomal miRNAs change at the time of LA GVHD onset after HSCT. Although a longitudinal series of blood samples was taken from 2 patients only, upregulation of miR-128 expression level preceded of LA GVHD incidence. In particular, our results suggest that biological alterations of miR-128 in exosomes may serve as potent predictive biomarkers for close monitoring of the onset of LA GVHD. We consider that this work could provide a foundation for revealing the importance of miRNA in the pathophysiology of LA GVHD. GVHD remains a major cause of transplant-associated complications, which can lead not only to lethal organ damage, but also to a decreased quality of life. Inhibition of certain exosomal miRNAs, such as miR-128, might be considered as a novel potential treatment of LA GVHD to target the immune system and inflammatory signaling, to control GVHD. However, further studies are needed to elucidate the molecular mechanisms behind the development of LA GVHD.

4. Materials and Methods

4.1. Ethics

The use of patient samples was approved by the Institutional Review Board of Tokyo Medical University (IRB no. 1979, approved on 28 March 2011). Written informed consent was obtained from all of the participants before the collection of specimens, in accordance with the Declaration of Helsinki.

4.2. Diagnosis of Late-Onset Acute GVHD

The diagnosis of LA GVHD here was based on clinical symptoms or was histologically proven by biopsy in the target organs, as described previously [1]. Patients with typical manifestations of classic acute GVHD occurring beyond 100 days after transplantation and persistent, recurrent acute GVHD were considered as having LA GVHD, in accordance with the new classification based on NIH consensus criteria [1].

4.3. Isolation of Exosome Fractions from Plasma Sample of Patients

Ten patients with hematological malignancies who underwent allogeneic HSCT between February 2012 and November 2013 were included in this study, which consisted of five patients with LA GVHD (gut, n = 2; liver, n = 2; skin and gut, n = 1) and five stable patients without GVHD symptoms (non-GVHD) under immunosuppression. All samples were collected sequentially after HSCT. For patients with LA GVHD, we evaluated blood samples obtained at the onset of this disease. No treatment for LA GVHD had been given at the time of sampling. For comparison, blood samples from HSCT patients without GVHD symptoms and other complications were used. In addition, eight healthy controls were also randomly selected in our hospital. We extracted exosome fractions from 200 µL of plasma using Total Exosome Isolation Reagent (Invitrogen, Carlsbad, CA, USA) and then assessed exosomal miRNA expression in a subset of samples, as described previously [11].

4.4. Exosomal miRNA Profiling

Total RNA was isolated from exosomes fractions of five patients with LA GVHD and five patients without GVHD (non-GVHD) using the miRNeasy Mini kit (Qiagen, Hilden, Germany). Exosomes were diluted with 700 μL of QiaZol. After 5 min of incubation, 1 μL of 1 nM ath-miR-159 (Hokkaido System Science, Hokkaido, Japan) was added to the aliquot followed by vortexing for 30 s and incubation for 5 min. Subsequent phenol extraction and centrifuge filtration were performed, in accordance with the manufacturer’s instructions (Qiagen, Hilden, Germany). The RT reaction and pre-amplification step were set up in accordance with the manufacturer’s instructions (Thermo Fisher Science, Waltham, MA, USA). miRNAs were reverse-transcribed with Megaplex Prime Pools (Human Pools A v2:1; Thermo Fisher Sciences, Waltham, MA, USA).

4.5. TaqMan Loq-Density Array Screening

The miRNA expression profiles were determined with a TaqMan miRNA Array Human Card A (Thermo Fisher Sciences, Waltham, MA, USA). Quantitative RT-PCR was performed on an Applied Biosystem 7900 HT thermocycler, in accordance with the manufacturer’s recommended program [31]. Using SDS2.2 software and Data Assist (version 3.01, Thermo Fisher Sciences, Waltham, MA, USA), the expression of exosomal miRNAs was calculated based on cycle threshold (Ct) values normalized by those of ath-miR-159, which was spiked in each exosomal sample. Data analysis was performed using GeneSpring® software (Version 12.1, Agilent Technologies, Palo Alto, CA, USA) and R software (https://www.r-project.org, accessed on 27 Arpil 2017). The Benjamini–Hochberg algorithm was used for the estimation of false discovery rates, as we reported previously [31].

4.6. Real-Time Quantitative RT-PCR for Candidate miRNAs

We determined the levels of the individual miRNAs by real-time quantitative RT-PCR with a TaqMan MicroRNA Assay (Thermo Fisher Sciences, Waltham, MA, USA) and the following miRNA-specific stem-loop primers: hsa-miR-128 (assay ID: 002216, Ambion, Austin, TX, USA) and hsa-miR-125b (assay ID: 000449, Ambion, Austin, TX, USA). Subsequently, quantitative real-time PCR was performed with an ABI Prism 7900 sequence detection system (Thermo Fisher Sciences, Waltham, MA, USA). The reaction was initiated by incubation at 95 °C for 2 min, followed by 50 cycles of 95 °C for 15 s and then 60 °C for 1 min. All reactions were run in duplicate. Mean (Ct) values for all miRNAs were quantified with sequence detection system software (SDS version 1.02; Thermo Fisher Sciences, Waltham, MA, USA) [31]. The expression of all miRNAs was normalized to miR-159, which was stably detected in all samples, yielding a −ΔΔCt value that was calculated by subtracting the −ΔΔCt value of the normal samples from the respective −ΔCt values of patient samples. The expression of all miRNAs was normalized using the 2−ΔΔCt method.
Following the identification of differentially expressed miRNAs, the predicted target genes of these altered miRNAs were investigated using the experimentally validated miRNA–target interactions database MirTarBase (http://mrtarbase.mbc.nctu.edu.tw/, accessed on 11 September 2017).

4.7. Statistical Analysis

Data are expressed as mean ± standard deviation (SD). Multiple groups were compared by one-way analysis of variance (ANOVA) with post hoc Tukey’s test. Statistical analysis was performed using GeneSpring software, and miRNAs with both a ΔCt of either >1.0 or <−1.0 and a P value of less than 0.05 were deemed to be differentially expressed. The area under the curve (AUC) of the receiver-operating characteristic (ROC) curve was calculated. GraphPad Prism software (version 5c for Macintosh; GraphPad Software Inc., La Jolla, CA, USA) was used for statistical analysis.

Author Contributions

S.Y., K.O. and J.H.O. designed the research and wrote the paper. T.U. and C.K. performed the research and analyzed the data. Y.S., M.G. and D.A. provided patient samples.

Funding

This study was supported in part by the Private University Strategic Research-Based Support Project (S1311016, S1511011) from the Ministry of Education, Culture, Sports, Science and Technology (Tokyo, Japan). We thank Edanz Group (www.edanzediting.com/ac) for editing a draft of this manuscript.

Acknowledgments

Kazuma Ohyashiki received research support from Celgene KK, Chugai Pharmaceutical KK, and Janssen Pharma KK.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Filipovich, A.; Weisdorf, D.; Pavletic, S.; Socie, G.; Wingard, J.R.; Lee, S.J.; Martin, P.; Chien, J.; Przepiorka, D.; Couriel, D.; et al. National Institutes of Health consensus development project on criteria for clinical trials in chronic graft-versus-host disease: I. Diagnosis and staging working group report. Biol. Blood Marrow Transplant. 2005, 11, 945–956. [Google Scholar] [CrossRef] [PubMed]
  2. Omer, A.K.; Weisdorf, D.J.; Lazaryan, A.; Shanley, R.; Blazar, B.R.; MacMillan, M.L.; Brunstein, C.; Bejanyan, N.; Arora, M. Late Acute Graft-versus-Host Disease after Allogeneic Hematopoietic Stem Cell Transplantation. Biol. Blood Marrow Transplant. 2016, 22, 879–883. [Google Scholar] [CrossRef] [PubMed]
  3. Cho, B.S.; Lee, S.E.; Song, H.H.; Lee, J.H.; Yahng, S.A.; Eom, K.S.; Kim, Y.-J.; Kim, H.-J.; Lee, S.; Min, C.-K.; et al. Graft-versus-tumor according to type of graft-versus-host disease defined by National Institute of Health consensus criteria and associated outcome. Biol. Blood Marrow Transplant. 2012, 18, 1136–1143. [Google Scholar] [CrossRef] [PubMed]
  4. Paczesny, S.; Krijanovski, O.I.; Braun, T.M.; Choi, S.W.; Clouthier, S.G.; Kuick, R.; Misek, D.E.; Cooke, K.R.; Kicko, C.L.; Weyand, A.; et al. A biomarker panel for acute graft-versus-host disease. Blood 2009, 113, 273–278. [Google Scholar] [CrossRef] [PubMed]
  5. Nelson, P.R., Jr.; Khawaja, M.R.; Perkins, S.M.; Mumaw, C.L.; Orschell, C.; Paczesny, S. Prognostic biomarkers for acute graft-versus-host disease risk after cyclophosphamide-fludarabine nonmyeloablative allotransplantation. Biol Blood Marrow Transplant 2014, 20, 1861–1864. [Google Scholar] [CrossRef] [PubMed]
  6. Vader Lugt, M.T.; Braun, T.M.; Hanash, S.; Ritz, J.; Ho, V.T.; Antin, J.H.; Zhang, Q.; Wong, C.-H.; Wang, H.; Chin, A.; et al. ST2 as a marker for risk of therapy-resistant graft-versus-host disease and death. N. Engl. J. Med. 2013, 369, 529–539. [Google Scholar] [CrossRef] [PubMed]
  7. McDonald, G.B.; Tabellini, L.; Storer, B.E.; Lawler, R.L.; Martin, P.J.; Hansen, J.A. Plasma biomarkers of acute GVHD and nonrelapse mortality: Predictive value of measurements before GVHD onset and treatment. Blood 2015, 126, 113–120. [Google Scholar] [CrossRef] [PubMed]
  8. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  9. Mathivanan, S.; Ji, H.; Simpson, R.J. Exosome: Extracellular organelles important in intracellular communication. J. Proteom. 2010, 73, 1907–1920. [Google Scholar] [CrossRef] [PubMed]
  10. Yuan, A.; Farber, E.L.; Rapoport, A.L.; Tejada, D.; Deniskin, R.; Akhmedov, N.B.; Farber, D.B. Transfer of microRNAs by embryonic stem cell microvesicles. PLoS ONE 2009, 4, e4722. [Google Scholar] [CrossRef] [PubMed]
  11. Umezu, T.; Tadokoro, H.; Azuma, K.; Yoshizawa, S.; Ohyashiki, K.; Ohyashiki, J.H. Exosomal miR-135b shed from hypoxic multiple myeloma cells enhances angiogenesis by targeting factor-inhibiting HIF-1. Blood 2014, 124, 3748–3757. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Valadi, H.; Ekstrom, K.; Bossios, A.; Sjostrand, M.; Lee, J.J.; Lotvall, J.O. Exosome-mediated transfer of mRNAa and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol. 2007, 9, 654–659. [Google Scholar] [CrossRef] [PubMed]
  13. Qu, Z.; Wu, J.; Wu, J.; Ji, A.; Qiang, G.; Jiang, Y.; Jiang, C.; Ding, Y. Exosomal miR-665 as a novel minimally invasive biomarker for hepatocellular carcinoma diagnosis and prognosis. Oncotarget 2017, 8, 80666–80678. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sueta, A.; Yamamoto, Y.; Tomiguchi, M.; Takeshita, T.; Yamamoto-Ibusuki, M.; Iwase, H. Differential expression of exosomal miRNAs between breast cancer patients with and without recurrence. Oncotarget 2017, 8, 69934–69944. [Google Scholar] [CrossRef] [PubMed]
  15. Ranganathan, P.; Heaphy, C.E.; Costinean, S.; Stauffer, N.; Na, C.; Hamadani, M.; Santhanam, R.; Mao, C.; Taylor, P.A.; Sandhu, S.; et al. Regulation of acute graft-versus-host disease by microRNA-155. Blood 2012, 119, 4786–4797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Liu, X.; Robinson, S.N.; Setoyama, T.; Tung, S.S.; D’Abundo, L.; Shah, M.Y.; Yang, H.; Yvon, E.; Shah, N.; Yang, H.; et al. FOXP3 is a direct target of miR15a/16 in umbilical cord blood regulatory T cells. Bone Marrow Transplant 2014, 49, 793–799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Wu, Y.; Heinrichs, J.; Bastian, D.; Fu, J.; Nguyen, H.; Schutt, S.; Liu, Y.; Jin, J.; Liu, C.; Li, Q.-J.; et al. MicroRNA-17-92 controls T-cell responses in graft-versus-host disease and leukemia relapse in mice. Blood 2015, 126, 1314–1323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Xiao, B.; Wang, Y.; Li, W.; Baker, M.; Guo, J.; Corbet, K.; Tsalik, E.L.; Li, Q.-J.; Palmer, S.M.; Wood, C.W.; et al. Plasma microRNA signature as a noninvasive biomarker for acute graft-versus-host disease. Blood 2013, 122, 3365–3375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wang, Y.; Zhao, X.; Ye, X.; Luo, H.; Zhao, T.; Diao, Y.; Zhang, H.; Lv, M.; Zhang, W.; Huang, X.; et al. Plasma microRNA-586 is a new biomarker for acute graft-versus-host disease. Ann. Hematol. 2015, 94, 1505–1514. [Google Scholar] [CrossRef] [PubMed]
  20. Crossland, R.E.; Norden, J.; Juric, M.K.; Green, K.; Pearce, K.F.; Lendrem, C.; Greinix, H.T.; Dickinson, A.M. Expression of Serum microRNAs is Altered During Acute Graft-versus-Host Disease. Front. Immunol. 2017, 8, 308. [Google Scholar] [CrossRef] [PubMed]
  21. Sun, Y.; Oravecz-Wilson, K.; Mathewson, N.; Wang, Y.; McEachin, R.; Liu, C.; Toubai, T.; Wu, J.; Rossi, C.; Braun, T.; et al. Mature T cell responses are controlled by microRNA-142. J. Clin. Investig. 2015, 125, 2825–2840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Godlewski, J.; Nowicki, M.O.; Bronisz, A.; Williams, S.; Otsuki, A.; Nuovo, G.; Raychaudhury, A.; Newton, H.B.; Chiocca, E.A.; Lawler, S. Targeting of the Bmi-1 oncogene/stem cell renewal factor by microRNA-128 inhibits glioma proliferation and self-renewal. Cancer Res. 2008, 68, 9125–9130. [Google Scholar] [CrossRef] [PubMed]
  23. Cui, J.G.; Zhao, Y.; Sethi, P.; Li, Y.Y.; Mahta, A.; Culicchia, F.; Lukiw, W.J. Micro-RNA-128 (miRNA-128) down-regulation in glioblastoma targets ARP5 (ANGPTL6), Bmi-1 and E2F-3a, key regulators of brain cell proliferation. J. Neurooncol. 2010, 98, 297–304. [Google Scholar] [CrossRef] [PubMed]
  24. Jin, M.; Zhang, T.; Liu, C.; Badeaux, M.A.; Liu, B.; Liu, R.F.; Jeter, C.; Chen, X.; Vlassov, A.V.; Tang, D.G. miRNA-128 suppresses prostate cancer by inhibiting BMI-1 to inhibit tumor-initiating cells. Cancer Res. 2014, 74, 4183–4195. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, D.; He, X.; Zou, J.; Guo, P.; Jiang, S.; Lv, N.; Alexseyev, Y.; Luo, L.; Luo, Z. Negative regulation of Bmi-1 by AMPK and implication in cancer progression. Oncotarget 2016, 7, 6188–6200. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, W.X.; Wilfred, B.R.; Hu, Y.; Stromberg, A.J.; Nelson, P.T. Anti-Argonaute RIP-Chip shows that miRNA transfections alter global patterns of mRNA recruitment to microribonucleoprotein complexes. RNA 2010, 16, 394–404. [Google Scholar] [CrossRef] [PubMed]
  27. Hafner, M.; Landthaler, M.; Burger, L.; Khorshid, M.; Hausser, J.; Berninger, P.; Rothballer, A.; Ascano, M., Jr.; Jungkamp, A.-C.; Munschauer, M.; et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell 2010, 141, 129–141. [Google Scholar] [CrossRef] [PubMed]
  28. Lia, G.; Brunello, L.; Bruno, S.; Carpanetto, A.; Omede, P.; Festuccia, M.; Tosti, L.; Maffini, E.; Giaccone, L.; Arpinati, M.; et al. Extracellular vesicles as potentiall biomarkers of acute graft-vs-host disease. Leukemia 2018, 32, 765–773. [Google Scholar] [CrossRef] [PubMed]
  29. Ohyashiki, K.; Umezu, T.; Katagiri, S.; Kobayashi, C.; Azuma, K.; Tauchi, T.; Okabe, S.; Fukuoka, Y.; Oyashika, J.H. Down-regulation of plasma miR-125 in chronic myeloid leukemia patients with successful discontinuation of imatinib. Int. J. Mol. Sci. 2016, 17, 570. [Google Scholar] [CrossRef] [PubMed]
  30. Kordelas, L.; Rebmann, V.; Ludwig, A.K.; Radtke, S.; Ruesing, J.; Doeppner, T.R.; Epple, M.; Horn, P.A.; Beelen, D.W.; Giebel, B. MSC-derived exosomes: a novel tool to treat therapy-refractory graft-versus-host disease. Leukemia 2014, 28, 970–973. [Google Scholar] [CrossRef] [PubMed]
  31. Ohyashiki, J.H.; Ohtsuki, K.; Mizoguchi, I.; Yoshimoto, T.; Katagiri, S.; Umezu, T.; Ohyashiki, K. Downregulated microRNA-148b in circulating PBMCs in chronic myeloid leukemia patients with undetectable minimal residual disease: A possible biomarker to discontinue imatinib safely. Drug Des. Dev. Ther. 2014, 8, 1151–1159. [Google Scholar] [CrossRef] [PubMed]
  32. Yoshizawa, S.; Ohyashiki, J.H.; Ohyashiki, M.; Umezu, T.; Suzuki, K.; Inagaki, A.; Lida, S.; Ohyashiki, K. Downregulated plasma miR-92a levels have clinical impact on multiple myeloma and related disorders. Blood Cancer J. 2012, 2, e53. [Google Scholar] [CrossRef] [PubMed]
  33. Ohyashiki, K.; Umezu, T.; Yoshizawa, S.; Ito, Y.; Ohyashiki, M.; Kawashima, H.; Tanaka, M.; Kuroda, M.; Ohyashiki, J.H. Clinical impact of down-regulated plasma miR-92a levels in non-Hodgkin’s lymphoma. PLoS ONE 2011, 6, e16408. [Google Scholar] [CrossRef] [PubMed]
  34. Smirnova, L.; Gräfe, A.; Seiler, A.; Schumacher, S.; Nitsch, R.; Wulczyn, F.G. Regulation of miRNA expression during neural cell specification. Eur. J. Neurosci. 2005, 21, 1469–1477. [Google Scholar] [CrossRef] [PubMed]
  35. Khan, A.P.; Poisson, L.M.; Bhat, V.B.; Fermin, D.; Zhao, R.; Kalyana-Sundaram, S.; Michailidis, G.; Nesvizhskii, A.I.; Omenn, G.S.; Chinnaiyan, A.M.; et al. Quantitative proteomic profiling of prostate cancer reveals a role for miR-128 in prostate cancer. Mol. Cell Proteom. 2010, 9, 298–312. [Google Scholar] [CrossRef] [PubMed]
  36. Woo, H.H.; László, C.F.; Greco, S.; Chambers, S.K. Regulation of colony stimulating factor-1 expression and ovarian cancer cell behavior in vitro by miR-128 and miR-152. Mol. Cancer 2012, 11, 58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Roth, P.; Wischhsen, J.; Happold, C.; Chandran, P.A.; Hofer, S.; Eisele, G.; Weller, M.; Keller, A. A specific miRNA signature in the peripheral blood of glioblastoma patients. J. Neurochem. 2011, 118, 449–457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Razo, A.; Dontje, B.; Vellenga, E.; de Haan, G.; Schuringa, J.J. Long-term maintenance of human hematopoietic stem/progenitor cells by expression of BMI1. Blood 2008, 111, 2621–2630. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 02493 g001 550
Figure 1. Flow diagram for the identification of late-onset acute GVHD-specific miRNAs. The exosomes fractions were isolated from 18 plasma samples of HSCT patients and healthy volunteers, and the expression profile of 381 exosomal miRNAs were analyzed. To validate candidate miRNAs, we performed TLDA and further narrowed down the miRNAs to those significantly differentially expressed by at least two-fold.
Figure 1. Flow diagram for the identification of late-onset acute GVHD-specific miRNAs. The exosomes fractions were isolated from 18 plasma samples of HSCT patients and healthy volunteers, and the expression profile of 381 exosomal miRNAs were analyzed. To validate candidate miRNAs, we performed TLDA and further narrowed down the miRNAs to those significantly differentially expressed by at least two-fold.
Ijms 19 02493 g001
Ijms 19 02493 g002 550
Figure 2. miRNA expression profiling using TaqMan miRNA Array (Thermo Fisher Science, Waltham, MA, USA). Exosome fraction taken from plasma samples of 10 patients who underwent allogeneic HSCT (5 late-onset acute GVHD and 5 non-GVHD), with 8 normal controls being included. A differential expression pattern was found between late-onset acute GVHD patients and non-GVHD patients. Using SDS (version 2.4) and DataAssist software (version 3.01, Thermo Fisher, version, Waltham, MA, USA), we calculated the expression of miRNAs based on their Ct value normalized by the Ct value of ath-miR-159. Data were analyzed with GeneSpring software.
Figure 2. miRNA expression profiling using TaqMan miRNA Array (Thermo Fisher Science, Waltham, MA, USA). Exosome fraction taken from plasma samples of 10 patients who underwent allogeneic HSCT (5 late-onset acute GVHD and 5 non-GVHD), with 8 normal controls being included. A differential expression pattern was found between late-onset acute GVHD patients and non-GVHD patients. Using SDS (version 2.4) and DataAssist software (version 3.01, Thermo Fisher, version, Waltham, MA, USA), we calculated the expression of miRNAs based on their Ct value normalized by the Ct value of ath-miR-159. Data were analyzed with GeneSpring software.
Ijms 19 02493 g002
Ijms 19 02493 g003a 550Ijms 19 02493 g003b 550
Figure 3. This microRNA expression according to late-onset acute GVHD (LA GVHD) vs. non-GVHD vs. normal controls at the onset of late-onset acute GVHD. Box plots depict the 75th percentile, median, and 25th percentile values; and whiskers represent minimum to maximum expression. p values were calculated using the independent two-sample t-test. (A) The expression levels of three miRNAs (miR-128, miR-125b, and miR-423-5p) were significantly higher in the exosomes of LA GVHD than in normal controls. The expression level of miR-423-5p showed downregulation in the non-GVHD groups when compared with normal controls. Accordingly, two miRNAs (miR-128, miR-125b) were chosen as candidate miRNAs specific for LA GVHD. (B) No significant differences of expression of seven miRNAs (miR-19a, miR-142-3p, miR-193b, miR-30c, miR-193a, miR-191, and miR-574-3p) were observed between LA GVHD and normal controls.
Figure 3. This microRNA expression according to late-onset acute GVHD (LA GVHD) vs. non-GVHD vs. normal controls at the onset of late-onset acute GVHD. Box plots depict the 75th percentile, median, and 25th percentile values; and whiskers represent minimum to maximum expression. p values were calculated using the independent two-sample t-test. (A) The expression levels of three miRNAs (miR-128, miR-125b, and miR-423-5p) were significantly higher in the exosomes of LA GVHD than in normal controls. The expression level of miR-423-5p showed downregulation in the non-GVHD groups when compared with normal controls. Accordingly, two miRNAs (miR-128, miR-125b) were chosen as candidate miRNAs specific for LA GVHD. (B) No significant differences of expression of seven miRNAs (miR-19a, miR-142-3p, miR-193b, miR-30c, miR-193a, miR-191, and miR-574-3p) were observed between LA GVHD and normal controls.
Ijms 19 02493 g003aIjms 19 02493 g003b
Ijms 19 02493 g004 550
Figure 4. (A) The expression level of miR-128 was significantly upregulated in the exosomes of late-onset acute GVHD when compared with that in non-GVHD (p = 0.0027) and in healthy subjects (p < 0.005). (B) The expression level of mir-125b was significantly upregulated in the exosomes of late-onset acute GVHD when compared with that in non-GVHD (p = 0.0289) and in healthy subjects (p = 0.0104).
Figure 4. (A) The expression level of miR-128 was significantly upregulated in the exosomes of late-onset acute GVHD when compared with that in non-GVHD (p = 0.0027) and in healthy subjects (p < 0.005). (B) The expression level of mir-125b was significantly upregulated in the exosomes of late-onset acute GVHD when compared with that in non-GVHD (p = 0.0289) and in healthy subjects (p = 0.0104).
Ijms 19 02493 g004
Ijms 19 02493 g005 550
Figure 5. Association between exosomal miRNA and late-onset acute GVHD. (A) ROC curve for the identification of patients with late-onset acute GVHD versus non-GVHD using miR-128. (B) ROC curve for the identification of patients with late-onset acute GVHD versus non-GVHD using miR-125b.
Figure 5. Association between exosomal miRNA and late-onset acute GVHD. (A) ROC curve for the identification of patients with late-onset acute GVHD versus non-GVHD using miR-128. (B) ROC curve for the identification of patients with late-onset acute GVHD versus non-GVHD using miR-125b.
Ijms 19 02493 g005
Ijms 19 02493 g006 550
Figure 6. Bar graphs show relative expression levels of each miRNA at each time point after HSCT. (A,B) The dynamic changes in exosomal miR-128 prior to the late-onset acute GVHD diagnosis. The information of PN 3 was day +7, day +25, day +46, day +70, day +105, and onset. The level of miR-128 was upregulated over the cut-off level at day +105 (A). The information of PN 4 was day +11, day +81, day +111, and onset. The level of miR-128 was upregulated over the cut-off level at day +111 (B). (C,D) The changes of exosomal miR-125b varied widely after HSCT. * p < 0.005, ns: Not significant.
Figure 6. Bar graphs show relative expression levels of each miRNA at each time point after HSCT. (A,B) The dynamic changes in exosomal miR-128 prior to the late-onset acute GVHD diagnosis. The information of PN 3 was day +7, day +25, day +46, day +70, day +105, and onset. The level of miR-128 was upregulated over the cut-off level at day +105 (A). The information of PN 4 was day +11, day +81, day +111, and onset. The level of miR-128 was upregulated over the cut-off level at day +111 (B). (C,D) The changes of exosomal miR-125b varied widely after HSCT. * p < 0.005, ns: Not significant.
Ijms 19 02493 g006
Table
Table 1. Characteristics of hematopoietic stem cell transplantation (HSCT) patients.
Table 1. Characteristics of hematopoietic stem cell transplantation (HSCT) patients.
DiagnosisAge (y)GenderDonorHLA Sero CompatibilityStem Cell SourceConditioningGVHD ProphylaxisDays of LA GVHD OnsetDays of SamplingTarget OrganGrade
Late-Onset Acute GVHD Group
PN 1AML33MMatched unrelated6/6BMMyeloablativeCsA + st-MTX+238+248LiverIV
PN 2NHL38MMatched related6/6PBSCReduced intensityCsA + st-MTX+180+221LiverIV
PN 3AML57MMismatched5/6BMMyeloablativeTac + st-MTX+117+136GutIII
PN 4ALL56FMismatched4/6Umbilical cord bloodMyeloablativeTac+119+132Skin + gutIV
PN 5ALL60FMismatched4/6Umbilical cord bloodMyeloablativeTac+165+174gutII
Non-GVHD Group
PN 6ALL62MMatched unrelated6/6BMReduced intensityCsA + st-MTX +42
PN 7ALL26FMatched unrelated6/6BMMyeloablativeCsA + st-MTX +41
PN 8AML60MMismatched4/6Umbilical cord bloodMyeloablativeTac +49
PN 9AML62MMismatched5/6PBSCReduced intensityTac + st-MTX +41
PN 10AML70FMismatched4/6Umbilical cord bloodReduced intensityTac +45
AML, acute myeloid leukemia; NHL, non-Hodgkin lymphoma; ALL, acute lymphoblastic leukemia; BM, bone marrow; PBSC, peripheral blood stem cell; CsA, cyclosporine A; Tac, tacrolimus; MTX, methotrexate.
Table
Table 2. Top 10 exosomal microRNAs differentially expressed between late-onset acute GVHD and non-GVHD.
Table 2. Top 10 exosomal microRNAs differentially expressed between late-onset acute GVHD and non-GVHD.
Exosomal miRNAsFold Changep Value
hsa-miR-423-5p49.070.000198
hsa-miR-19a37.330.000013
hsa-miR-142-3p28.350.000198
hsa-miR-12824.790.000111
hsa-miR-193b16.120.000209
hsa-miR-30c14.330.003157
hsa-miR-193a12.560.000111
hsa-miR-19112.30.000081
hsa-miR-125b10.70.001064
hsa-miR-574-3p6.130.000198
Table
Table 3. Predicted target genes for miR-128 identified by miRTarBase.
Table 3. Predicted target genes for miR-128 identified by miRTarBase.
miRNAExperimentally Validated Targets by Reporter Assay (miRTarBase)
hsa-miR-128BMI1, FBXW7, DCX, RELN, WEE1, TGFBR1, NEK2, SREBF1, SREBF2, ABCA1, ABCG1, RXRA

Share and Cite

MDPI and ACS Style

Yoshizawa, S.; Umezu, T.; Saitoh, Y.; Gotoh, M.; Akahane, D.; Kobayashi, C.; Ohyashiki, J.H.; Ohyashiki, K. Exosomal miRNA Signatures for Late-Onset Acute Graft-Versus-Host Disease in Allogenic Hematopoietic Stem Cell Transplantation. Int. J. Mol. Sci. 2018, 19, 2493. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092493

AMA Style

Yoshizawa S, Umezu T, Saitoh Y, Gotoh M, Akahane D, Kobayashi C, Ohyashiki JH, Ohyashiki K. Exosomal miRNA Signatures for Late-Onset Acute Graft-Versus-Host Disease in Allogenic Hematopoietic Stem Cell Transplantation. International Journal of Molecular Sciences. 2018; 19(9):2493. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092493

Chicago/Turabian Style

Yoshizawa, Seiichiro, Tomohiro Umezu, Yuu Saitoh, Moritaka Gotoh, Daigo Akahane, Chiaki Kobayashi, Junko H. Ohyashiki, and Kazuma Ohyashiki. 2018. "Exosomal miRNA Signatures for Late-Onset Acute Graft-Versus-Host Disease in Allogenic Hematopoietic Stem Cell Transplantation" International Journal of Molecular Sciences 19, no. 9: 2493. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms19092493

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop