|Home | About | Journals | Submit | Contact Us | Français|
The pathogenesis of chronic rejection, Bronchiolitis Obliterans Syndrome (BOS) following lung transplantation (LT) is poorly understood. We hypothesized that development of antibodies to HLA (DSA) is associated with dysregulation of microRNA (miRNA) that predisposes BOS. Towards this, miRNA profiling of mononuclear cells from 10 stable LT (DSA−BOS−), 10 LT with DSA+BOS− (DSA group) and 10 LT with DSA+BOS+ (BOS group) were performed. Prediction by mirPath indicated that differential miRNAs in DSA+BOS− compared to stable are significantly up-regulated (relative fold >2, p < 0.05) for TGF-β and B cell receptor signal pathways. A total of seventy-four miRNAs were up-regulated and six miRNAs were down regulated in LT with DSA+BOS+ when compared to stable (relative fold >2, p < 0.05). There was also significant enrichment of cell cycle and gap junction pathways. An inverse correlation between expression of two key miRNAs and their target genes were observed: miR-369-5p and miR-548d were down regulated in DSA+ LT while their gene targets in TGF-β signal pathways were up-regulated. In addition, miR-628-5p and miR-134 were down regulated and their target genes (B cell development) were up-regulated. Therefore, we conclude that alloimmunity induced changes in miRNAs affecting the TGF-β and B cell receptor signal pathways play important roles in BOS development.
Lung transplantation (LT) is a treatment option for end-stage lung diseases including idiopathic pulmonary fibrosis, chronic obstructive pulmonary disease, and cystic fibrosis (1,2). However, long-term survival remains a challenge due to development of Bronchiolitis Obliterans Syndrome (BOS) (3). Approximately 50% of LT recipients (LTx) develop BOS by five years after transplantation (3). Many etiologies for the pathogenesis of BOS have been proposed and it is accepted that development of BOS is due to immunological insults. Our study demonstrated a correlation between the development of antibodies (Abs) to donor mismatched HLA (DSA) during post-transplantation and the risk for BOS (4). Furthermore, Abs against HLA has been eluted from lungs that underwent chronic rejection (5). Further support for the possible role for anti-HLA in the immunopathogenesis of obliterative airway disease (OAD) has been provided by studies in which administration of anti-HLA class I intra-bronchially into produced OAD in the native lungs (6). However, the mechanisms leading to BOS following de novo development of anti-HLA still remains unknown.
MiRNAs are a class of non-coding small RNAs, 22–25 nucleotides in length, which bind to the 3′UTR of target genes and thereby repress translation and/or induce degradation of target gene mRNA (7). Aberrant expression of miRNAs is associated with initiation and progression of pathological processes including immune-system disorders like asthma and rheumatoid arthritis (8,9). Furthermore, miRNAs have been shown to influence toll-like receptor signaling, inflammatory gene expression, T cell and B cell differentiation, lineage specificity and development of fibrosis, cirrhosis, and graft loss (10). Studies also demonstrated that miRNA dysregulation may result in chronic antibody-mediated rejection following kidney transplantation (11–13). Therefore, to test the hypothesis that immune responses to alloantigens as evidenced by development of Abs to donor HLA will dysregulate miRNAs involved in regulating immune responses, we determined miRNA expression profiles following human LT, including the LTx who developed DSA to mismatched MHC and developed BOS.
Thirty LTx at Barnes Jewish Hospital/Washington University School of Medicine were enrolled between January 2009 and March 2013. This study was approved by the institutional review board and written informed consent was obtained. Among the 30 patients included, 10 developed DSA but had normal FEV1 (Forced Expiratory Volume in 1 s) and had no evidence of BOS two years after LT (DSA+BOS− group, 20 ± 10 months from LT); 10 developed DSA and had a clinical diagnosis of BOS (DSA+BOS+ group, 25 ± 9 months from LT), and 10 didn't develop DSA and did not have any evidence of BOS within a two year period (stable group, 22 ± 8 months from LT) (Table 1). Stable group and DSA+BOS− group were available for follow-up. There were no significant differences of mean fluorescent intensity (MFI) values of DSA and antibodies to MHC class I or class II in the patients between DSA+BOS− and DSA+BOS+ groups. There were also no significant differences (p > 0.05) in the development of antibodies to self-antigens (Ka1 T and Col V) among three groups. BOS was defined by the criteria set forth by the International Society for Heart & Lung Transplantation (ISHLT) (14). LTx from these three groups were similar in age, gender, race, underlying diseases, and had bilateral LT. The standard immunotherapy protocol consisted of tacrolimus, mycophenolate mofetil, and prednisone.
PBMCs were separated from blood using the Ficoll-hypaque density separation (Amersham, Uppsala, Sweden) (15). BAL fluid from bronchoscopy was centrifuged at 300 g and 4°C for 10 min to separate cells. There were no statistically significant differences in total blood mononuclear cell (PBMC) counts and PBMC profiles among stable, DSA+BOS−, and DSA+BOS+ patient groups. PBMCs and BAL cells were suspended in TRIzol (Agilent Technologies, Morrisville, NC) and stored at −80°C until RNA extraction. Total RNA was extracted using the mirVana miRNA isolation kit according to the manufacturer's specification (Ambion, Austin, TX).
The TaqMan Low-Density Array Human miRNA Panel v1.0 (Applied Biosystems, Foster City, CA) was utilized for global miRNA profiling. The panel includes two 384-well microfluidic cards (human miRNA pool A and pool B) that contain primers and probes for 746 different human miRNAs in addition to six small nucleolar RNAs that function as endogenous controls for data normalization. RNA extracted from PBMC among stable, DSA+BOS− and DSA+BOS+ patients were pooled separately and used for miRNA profiling. Equal quantity of RNA (30 ng) from PBMCs (n = 10) from each of the groups were pooled and reverse transcribed for cDNA synthesis using the TaqMan Multiplex RT set (Applied Biosystems) for TaqMan Array Human miRNA Panels. Each RT reaction was diluted 62.5-fold with water and 55 μL of each diluted product was combined with 55 μL of TaqMan 2× Universal PCR Master Mix, No AmpErase UNG (Applied Biosystems). One-hundred microliters of the sample/master mix for each multiplex pool were loaded into fill reservoirs on the microfluidic card. The array was then centrifuged and mechanically sealed with the Applied Biosystems sealer device. Quantitative PCR was performed on an Applied BioSystems 7900HT thermocycler (Applied Biosystems) using the manufacturer's recommendations.
Independent sets of PBMC and BAL were used for qPCR confirmation. Total RNA was isolated, and six individual miRNAs were further quantitated by TaqMan miRNA assays (Applied Biosystems) and RNU44 was used as controls to normalize data.
To determine the expression of target genes, cDNA from independent sets of PBMC was synthesized using Super-Script II reverse transcriptase with random hexamers (Life Technologies, Carlsbad, CA). Real-time PCR was performed using iQ SYBR Green Supermix (Bio-Rad, Hercules, CA). GAPDH was detected as endogenous control.
DIANA-mirPath (16) was employed to perform the enrichment analysis of predicted target genes by one or more miRNAs in biological pathways. The software performs analysis of multiple miRNA target genes to all known KEGG pathways (Kyoto Encyclopedia of Genes and Genomes) based on two algorithms, microT-CDS (17,18) and miRTarBase (19). The input dataset enrichment was performed by a Pearson's Chi-squared test, and FDR correction was used for the number of pathways. The graphical output of the program provides an overview of the parts of the pathway modulated by selected miRNAs facilitating the interpretation and presentation of the results. The statistical significance associated with the identified biological pathways was calculated by the mirPath software (http://microrna.gr/mirpath).
The selection of 10 samples within each group was based upon sample availability and statistical test power. The efficiency of pooling total RNA in microarray experiments was based upon published methods (20,21) and the RNA pooling enabled us to reduce number of arrays without a loss of precision. Data analysis on mi RNA expression levels from TaqMan Low-Density Array was performed by the SDS software version 2.2.2 (Applied Biosystems). Assays that had Ct values >35 were removed from the analysis. The delta Ct values were calculated by using RNU44 as the endogenous control. We initially designed the experiments to minimize the impact of covariates by matching the samples for key confounding factors. Therefore, we used a univariate test to screen for differentially expressed miRNAs, and validated the result with a separate set of samples. A Student's t-test was performed to detect differentially expressed miRNAs. Bonferroni procedure was used to calculate adjusted p values to control false discovery rate (FDR) among stable, DSA+BOS− and DSAþBOS+ patients groups, respectively (22). Heatmap was performed by R software. The Mann–Whitney test was performed using GraphPad Prism 5 to determine the significance of miRNA levels of PBMC and BAL for independent validation. The association between expression levels of mRNAs and miRNAs was analyzed using the Pearson's correlation. All data were expressed as mean ± SE unless otherwise specified. p values less than 0.05 were considered significant.
To identify differentially expressed miRNAs in LT patients, we profiled the expression of 746 miRNAs using the TaqMan miRNA array. The study design is illustrated in Figure 1. The total RNA obtained from PBMCs from 10 individuals in each group (stable, DSA+BOS− and DSA+BOS+) was pooled with equal quantity (30 ng). The selection of pooled samples (10 per group) are approximately equivalent to 4–7 independent samples in each group, if samples were not pooled, and also assumes a wide range of values for the ratio of biological variations among the samples over the technical variations in the assays (21) (Figure S1). Global miRNA expression profiling was performed on the pooled RNA using TaqMan Low-density Human miRNA Array in combination with Megaplex RT and Megaplex pre-amplification techniques. Preliminary evaluation of candidate differentially expressed miRNAs was performed by Heatmap, potential target prediction and pathway enrichment and finally independent validation of differentially expressed miRNAs using the TaqMan miRNA assay.
To investigate the relative abundance of detected miRNAs, they were normalized in each pool (stable, DSA+BOS− and DSA+BOS+) to RNU44 in the corresponding pool. The results identified 349 miRNAs (~46.8%) (Ct <35 in at least one pool were classed as detectable), indicating that global miRNA expression profiling is a feasible method to identify aberrantly expressed miRNAs. Subsequently we investigated the miRNA expression profile after de novo development of DSA. Comparison of the miRNA expression profile in the DSA+BOS− to the stable patients demonstrated that 31 miRNAs were up regulated and 14 miRNAs were down regulated in the DSA+BOS− group (relative fold >2, p < 0.05) (Figure 2A and Table S1), suggesting that these miRNAs might be important in DSA mediated immune responses after LT. The heatmap, which represents the mean fold change of differential miRNA signature in DSA+BOS− patients, was generated by 45 differential miRNAs (31 over-expressed miRNAs and 14 under-expressed miRNAs in the DSA+BOS− group, Figure 2B). In order to examine the biologic pathways modulated by the differentially expressed miRNAs in DSA+BOS− patients, we applied DIANA-mirPath on the DSA-related deregulated miRNA signature. Thirty-two KEGG biological processes were significantly modulated by the differentially expressed miRNAs in DSA+BOS− patients (p < 0.05, FDR corrected). Among them, the TGF beta signal pathway (p < 0.0001) was the most prominent pathway enriched in quantities with DSA-related differentially expressed miRNA patterns (Table 2 and Figure S2). Furthermore, B cell receptor signal pathway (p = 0.0013, Figure S4), T cell receptor signal pathway (p < 0.0001), chemokine signaling pathway (p < 0.0001), cytokine-cytokine receptor interaction (p = 0.0075) were significantly modulated by the differentially expressed miRNAs in DSA+BOS− patients (Table 2). These results strongly suggest that immune response associated biologic pathways were affected following de novo development of DSA.
In order to identify the differentially expressed miRNAs in the DSA+BOS+ group, we compared miRNA expression levels between the DSA+BOS+ and the stable group. There were 74 miRNAs up regulated and 6 miRNAs down regulated in the DSA+BOS+ group, compared to those of stable (relative fold >2, p < 0.05) (Figure 2A and Table S2). In addition, 26 miRNAs were up regulated and 3 miRNAs were down regulated in both DSA+BOS+ versus stable and DSA+BOS− versus stable comparisons (relative fold >2, p < 0.05, FDR corrected) (Figure 3A), indicating that these differentially expressed miRNAs are associated with de novo development of DSA in LTx with BOS.
We next performed comparative pathway analysis using DIANA-mirPath on differentially expressed miRNAs in DSA+BOS+ patients. There were 35 biologic pathways enriched in differentially expressed miRNAs in DSA+BOS+ patients (p < 0.05, FDR corrected). Compared with the pathways enriched by differential miRNAs in DSA+BOS− patients, there are 25 dysregulated biologic pathways enriched in DSA- and BOS-related dysregulated miRNA signatures in DSA+BOS− and DSA+BOS+ patients (Table 2). In addition, the TGF beta signal pathway (DSA+BOS−: p < 0.0001, Figure S2 and DSA+BOS+: p < 0.0001, Figure S3) was also enriched prominently (Table 2). Furthermore, several other pathways involved in the dysregulation of immune-response including B cell receptor (Figure S5) and T cell receptor signal pathways, chemokine signaling pathway and cytokine-cytokine receptor interaction pathway were also affected (Table 2). These results demonstrate that the dysregulation of TGF beta signaling and immune response pathways persisted in LTx who had persistent DSA leading to BOS.
Finally, we found 10 additional signal pathways that were dysregulated only in the DSA+BOS+ patients. They were inositol phosphate metabolism (p < 0.0001), gap junction (p < 0.0001), cell cycle (p < 0.0001), pertussis (p = 0.0002), phosphatidylinositol signaling system (p = 0.0072), basal transcription factors (p = 0.0074), pentose phosphate pathway (p = 0.037), N-glycan biosynthesis (p = 0.037), Toxoplasmosis (p = 0.037), and epithelial cell signaling in Helicobacter pylori infection (p = 0.049) (Table 3).
By comparing the miRNA profiles between DSA+BOS+ and DSA+BOS− groups, we identified 52 candidate miRNAs up regulated and 6 miRNAs down regulated (relative fold >2, p < 0.05, Figure 2A and Table S3). We then compared two differentially expressed miRNA patterns from DSA+BOS− versus stable and DSA+BOS+ versus DSA+BOS−. Twelve miRNAs were up regulated in both comparisons (Figure 3B), indicating that these differentially expressed miRNAs may play important roles in DSA mediated BOS development.
Comparative pathway analysis by DIANA-mirPath demonstrated that 39 dysregulated biologic pathways were enriched (p < 0.05, FDR corrected, Table 4) in differentially expressed miRNAs in DSA+BOS+ LTx compared to DSA+BOS− LTx. Among them, Wnt signaling (p < 0.0001), TGF-beta signal (p < 0.0001), and B cell receptor signal pathways (p < 0.0001) were ranked as the most enriched by comparing miRNA profiles between the DSA+BOS+ and DSA+BOS− groups. In addition, inositol phosphate metabolism (p < 0.0001), gap junction (p < 0.0001), and cell cycle (p < 0.0001) were also enriched when comparing the DSA+BOS+ and DSA+BOS− groups (Table 4).
We employed TaqMan RT-PCR to validate the expression level of differentially expressed miRNAs. Six miRNAs (miR-518f, miR-98, miR-10a, miR-369-5p, miR-548d, and miR-133b) were selected for validation based on published association with T and B cell differentiation, lineage specificity and development of post-transplant fibrosis, cirrhosis, and graft loss (10). They were validated by TaqMan RT-qPCR using individual PBMC and BAL samples from three independent groups of 9 stable, 7 DSA+BOS−, and 10 DSA+BOS+ LTx. There were no significant differences in age, gender, race, underlying diseases, and medical treatment among the three groups analyzed (Table S4). In agreement with the miRNA array data, miR-518f, miR-98 and miR-10a were significantly increased in PBMCs from DSA+BOS− and DSA+BOS+ groups when compared to the stable group. MiR-369-5p, miR-548d, and miR-133b were significantly under expressed in the DSA+BOS− compared to the stable (Figure 4). Furthermore, four of them (miR-518f, miR-10a, miR-369-5p, and miR-133b) were found to be differentially expressed in BAL cells. Expression of miR-10a and miR-518f was up regulated in DSA+BOS− and DSA+BOS+ patients compared to in the stable. Expression of miR-369-5p and miR-133b was under expressed in DSA+BOS− and DSA+BOS+ patients relative to the stable (Figure 5). Taken together, independent validation by TaqMan RT-qPCR confirmed the validity of previous finding of differential expression of miRNAs using TaqMan Low-Density Array. In addition, the results demonstrate that miR-518f and miR-98 were up regulated in the DSA+BOS− group and continue to have higher levels in DSA+BOS+ group (Figures 4 and and55).
We further determined the potential target genes for differentially expressed miRNAs based on pathway analysis via the mirPath algorithm. We focused on the TGF beta signaling pathway for further investigation since up regulation of the TGF-beta signaling pathway has been known to be a risk factor for BOS after LT (23). Three down-regulated miRNAs (miR-133b, miR-369-3p, and miR-548d) were associated with up regulation of the TGF beta signaling pathway by regulating their gene targets (Figure 6A). Further RT-qPCR confirmed that miR-369-3p was down-regulated while its gene target, latent TGF-beta binding proteins (LTBP1) (encoding latent TGF-beta binding proteins), mRNA expression was up regulated in DSA+BOS+ versus stable patients (Pearson's correlation p-value: p < 0.001, R2 = 0.981, Figure 6B); Similarly, there was an inverse correlation between expression of miR-369-3p and LTBP1 mRNA in DSA+BOS− versus stable patients (p < 0.001, R2 = 0.969, Figure 6B). We further observed that miR-548d was down regulated, consistent with up regulation of its potential target, TGFRB1 (encoding transforming growth factor, beta receptor 1) mRNA expression in DSA+BOS+ versus stable patients (p < 0.001, R2 = 0.431, Figure 6C). An inverse correlation between expression of miR-548d and TGFBR1 mRNA were also found in DSA+BOS− patients compared to those in the stable group (p < 0.001, R2 = 0.672, Figure 6C). These results clearly demonstrate that differentially expressed miRNAs regulate the TGF beta signaling pathway which in turn leads to deregulation of their target genes, de novo DSA development and BOS.
B cell receptor signaling pathways play a crucial role in humoral responses and have been implicated in chronic rejection (24, 25). We therefore selected two differentially expressed miRNAs (miR-628–5p and miR-134) to determine their roles in B cell receptor signaling and development of DSA and BOS (Figure 7A). MiR-628–5p was down-regulated while the mRNA for its potential target gene, BTK (Bruton agammaglobulinemia tyrosine kinase), was up regulated in DSA+BOS− patients compared to stable (p < 0.001, R2 = 0.623, Figure 7B). There was an inverse correlation between expression of miR-628–5p and BTK mRNA in DSA+BOS+ versus stable patients (p < 0.001, R2 = 0.752, Figure 7B). Similarly, miR-134 which was reduced and its target gene, BLNK (B cell linker) mRNA expression was increased in DSA+BOS- subjects relative to that of the stable (p < 0.001, R2 = 0.992, Figure 7C). We also noted an inverse correlation between expression of miR-134 and BLNK mRNA found in DSA+BOS+ patients compared to those in the stable (p < 0.001, R2 = 0.863, Figure 7C). These data indicate that differentially expressed miRNAs regulating B cell receptor signal pathway can dysregulate their target genes in LTx with de novo development of DSA leading to BOS.
The goal of the present study was to investigate whether miRNAs are aberrantly expressed in LT recipients who develop DSA leading to BOS. In order to determine differential expression of miRNAs, we generated three pools of PBMCs from stable, DSA+BOS− and DSA+BOS+ LT patients. Using Taqman Low-Density miRNA Array method we evaluated the miRNA expression profiles and a relative change of at least twofold was deemed significant. Comprehensive pathway analyses indicated that dysregulated biologic processes including B cell, Toll-like, and T cell receptor pathways, TGF-β and chemokine signal pathways, and cytokine-cytokine receptor interactions were significantly altered in DSA+BOS− and DSA+BOS+ LT patients. Chronic allograft rejection has been attributed to immune responses to donor HLA (26,27). A strong correlation between de novo development of DSA and BOS following human LT has been reported (28,29). Results from the current study provide evidence that development of DSA after LT, can induce differential expression of miRNA that influences the overall immunity. We further demonstrate a differential expression of four miRNAs (miR-10a, miR-518f, miR-369-5p, and miR-133b) at the local site, that is, in the cells from BAL fluid. MiR-10a is known to be involved in the commitment and differentiation of hematopoietic cells (30,31), and has a crucial role in the stability and functioning of regulatory T cells (Tregs) (31). In this study, miR-10a expression was significantly increased in PBMC from DSA+BOS− and DSA+BOS+ patients. The miR-133b located in the IL17A/F locus is known to be co-regulated with IL-17 production in αβ and γδ T cells (32). We found that miR-133b was significantly under expressed in the DSA+BOS− group compared to that of stable. Therefore, dysregulation of these two miRNAs: miR-10a and miR-133b in LT patients possibly lead to Treg cells dysfunction and aberrant IL-17 production, which have been proposed as the key mediators in the development of BOS. Analysis of Treg and IL-17 dependent pathways following DSA development after LT and its correlation to BOS need further investigation.
Of particular importance is our finding of significant up regulation of the TGF beta signal pathway both in the DSA+BOS− as well as DSA+BOS+ groups. Persistent dysregulation in LT patients who developed BOS is of great interest as it suggests that development of alloimmunity induces this change early and persists during BOS development. TGF-β expression has been reported to be increased in patients with BOS even before the onset of airflow obstruction (33,34). TGF-β itself was not up regulated in the sample set of miRNAs chosen, but downstream mediators including the TGF beta receptors and LTBP1, were markedly up regulated (Figure 6B and Figure 6C). Additionally, the Wnt signaling pathway is also of interest since Wnts have diverse roles in governing cell fate, proliferation, migration, polarity, and death (35). The significance of the Wnt pathway in allograft rejection and in BOS following human LT is currently unknown. Moreover, DSA+BOS+ patients also demonstrated dysregulation of miRNAs, which are involved in the cell cycle and gap junction signaling pathways. It is of interest that the development of Abs to lung associated self-antigen K-alpha 1 tubulin, an epithelial gap junction, significantly correlates with the development of BOS (36).
There are limitations in the present study. First, using miRNA array methods to analyze miRNA expression patterns is an evolving area. The standard methods of statistical analysis and analysis of significance have not yet been firmly established. In particular, choosing the appropriate threshold for determining statistical significance is, in part, dependent on the purpose of the analysis. We evaluated the miRNA expression profiles with greater than twofold increase because of the optimization of numbers of candidate miRNAs and high specificity of results afforded by this level of statistical stringency. Secondly, we failed to validate the elevated expression of miR-875-5p, miR-548c-5p, and miR-548d-5p in the DSA+BOS+ patients. A probable reason for this discrepancy is that although the TaqMan miRNA array is promising for simultaneously measuring miRNA expression in a high-throughput manner (37,38) and the introduction of an miRNA cDNA pre-amplification step significantly reduces the amount of RNA needed (39). It has a potential for pre-amplification bias that may be responsible for variations in the Ct values. This may diminish the reproducibility with miRNAs of low abundance. Future studies should include stratification analysis of a larger cohort to eliminate sample selection biases. Finally, a larger sample size may be needed to more accurately determine how anti-HLA can dysregulate the transcriptional activation and regulation of miRNA expression patterns in LT patients.
In summary, our study, for the first time, demonstrated that the alloimmune response as evidenced by development of DSA following human LT leads to changes in miRNA expression affecting primarily the TGF-beta and B cell receptor signal pathways. This is noted early following LT when there is no evidence of BOS. These miRNAs changes persist and expand to other cell cycle signal pathways in LT patients during development of BOS. The clinical significance of these findings include potential use of miRNAs as early markers for identifying risk factors in development of BOS which needs further investigation.
Figure S1: Method employed for estimating the number of patients by pooled and non-pooled designs. The efficiency of pooling ten patients into one sample was determined based on biological and technical variations within the assay as described earlier (21). The plot shows equivalent number of patients when ten samples were not pooled for independent microRNA arrays. The varying lambda (λ = σb2/σe2) is the ratio of underlying biological variations among the ten samples over the technical variations among the assays.
Figure S2: The KEGG pathway “TGF-beta signaling pathway” was significantly altered between DSAþBOS-LTx patients and Stable patients subjects at baseline, with 17 miRNAs (hsa-let-7c, hsa-miR-98, hsa-miR-125a-5p, hsa-miR-137, hsa-miR-203, hsa-miR-299-5p, hsa-miR-338-3p, hsa-miR-381, hsa-miR-548c-5p, hsa-miR-874, hsa-miR-875-5p, hsa-miR-133b, hsa-miR-134, hsa-miR-433, hsa-miR-494, hsa-miR-503, hsa-miR-628-5p) targeting 46 genes. Targeted genes are indicated by yellow (one miRNA) and brown (multiple miRNAs) colors in the pathway map.
Figure S3: The KEGG pathway “TGF-beta signaling pathway” was significantly altered between DSA+BOS+LTx patients and Stable patients subjects at baseline, with 29 miRNAs (hsa-let-7c, hsa-miR-1, hsa-miR-98, hsa-miR-125a-5p, hsa-miR-137, hsa-miR-140-3p, hsa-miR-152, hsa-miR-193a-5p, hsa-miR-203, hsa-miR-329, hsa-miR-331-5p, hsa-miR-338-3p, hsa-miR-372, hsa-miR-381, hsa-miR-410, hsa-miR-422a, hsa-miR-450b-5p, hsa-miR-542-3p, hsa-miR-548c-5p, hsa-miR-548d-5p, hsa-miR-582-3p, hsa-miR-627, hsa-miR-671-3p, hsa-miR-219-1-3p, hsa-miR-376c, hsa-miR-511, hsa-miR-941, hsa-miR-875-5p, hsa-miR-542-5p) targeting 126 genes. Targeted genes are indicated by yellow (one miRNA) and brown (multiple miRNAs) colors in the pathway map.
Figure S4: The KEGG pathway “B cell receptor signaling pathway” was significantly altered between DSA+BOS-LTx patients and Stable patients subjects at baseline, with 20 miRNAs (hsa-let-7c, hsa-miR-98, hsa-miR-125a-5p, hsa-miR-137, hsa-miR-184, hsa-miR-203, hsa-miR-299-5p, hsa-miR-338-3p, hsa-miR-381, hsa-miR-548c-5p, hsa-miR-627, hsa-miR-874, hsa-miR-133b, hsa-miR-134, hsa-miR-433, hsa-miR-489, hsa-miR-494, hsa-miR-503, hsa-miR-542-5p, hsa-miR-628-5p) targeting 34 genes. Targeted genes are indicated by yellow (one miRNA) and brown (multiple miRNAs) colors in the pathway map.
Figure S5: The KEGG pathway “B cell receptor signaling pathway” was significantly altered between DSA+BOS+LTx patients and Stable patients subjects at baseline, with 24 miRNAs (hsa-let-7c, hsa-miR-1, hsa-miR-98, hsa-miR-125a-5p, hsa-miR-137, hsa-miR-140-3p, hsa-miR-152, hsa-miR-203, hsa-miR-329, hsa-miR-338-3p, hsa-miR-372, hsa-miR-381, hsa-miR-410, hsa-miR-450b-5p, hsa-miR-542-3p, hsa-miR-548c-5p, hsa-miR-548d-5p, hsa-miR-582-3p, hsa-miR-627, hsa-miR-219-1-3p, hsa-miR-376c, hsa-miR-511, hsa-miR-941, hsa-miR-542-5p) targeting 34 genes. Targeted genes are indicated by yellow (one miRNA) and brown (multiple miRNAs) colors in the pathway map.
Table S1: Differential miRNAs when DSA+BOS- versus stable.
Table S2: Differential miRNAs when DSA+BOS+ versus stable.
Table S3: Differential microRNAs when DSA+BOS+ versus DSA+BOS-.
Table S4: Demographic data on the lung transplant recipients in the replicative independent cohort.
We thank Billie Glasscock for her assistance in preparing this manuscript. This work is supported by NIH RO1 HL056643-13 and the BJC foundation (TM).