Multi-omics network characterization reveals novel microRNA biomarkers and mechanisms for diagnosis and subtyping of kidney transplant rejection
Journal of Translational Medicine volume 19, Article number: 346 (2021)
Kidney transplantation is an optimal method for treatment of end-stage kidney failure. However, kidney transplant rejection (KTR) is commonly observed to have negative effects on allograft function. MicroRNAs (miRNAs) are small non-coding RNAs with regulatory role in KTR genesis, the identification of miRNA biomarkers for accurate diagnosis and subtyping of KTR is therefore of clinical significance for active intervention and personalized therapy.
In this study, an integrative bioinformatics model was developed based on multi-omics network characterization for miRNA biomarker discovery in KTR. Compared with existed methods, the topological importance of miRNA targets was prioritized based on cross-level miRNA-mRNA and protein–protein interaction network analyses. The biomarker potential of identified miRNAs was computationally validated and explored by receiver-operating characteristic (ROC) evaluation and integrated “miRNA-gene-pathway” pathogenic survey.
Three miRNAs, i.e., miR-145-5p, miR-155-5p, and miR-23b-3p, were screened as putative biomarkers for KTR monitoring. Among them, miR-155-5p was a previously reported signature in KTR, whereas the remaining two were novel candidates both for KTR diagnosis and subtyping. The ROC analysis convinced the power of identified miRNAs as single and combined biomarkers for KTR prediction in kidney tissue and blood samples. Functional analyses, including the latent crosstalk among HLA-related genes, immune signaling pathways and identified miRNAs, provided new insights of these miRNAs in KTR pathogenesis.
A network-based bioinformatics approach was proposed and applied to identify candidate miRNA biomarkers for KTR study. Biological and clinical validations are further needed for translational applications of the findings.
Kidney transplantation is the most effective procedure for clinical treatment of patients with end-stage kidney failure. Based on the stable work of transplanted kidney, electrolyte and acid–base balance could be maintained and patients no longer need the support of dialysis to eliminate the metabolic waste. However, due to the immune response of human body, kidney transplant rejection (KTR), including antibody-medicated rejection (AMR) and T cell-mediated rejection (CMR), is observed to have irreversible effects or even results in the loss of renal allograft function . Although a series of clinical symptoms, e.g., increased serum creatinine, decreased urine output, distending pain in the transplantation area, and fever, are measured for indicating KTR occurrence, the specificity is still limited since the infection and drug response could cause similar phenotypes. Thus, screening sensitive biomarkers to predict and monitor the development of rejection after allogeneic kidney transplantation is of great significance.
MicroRNAs (miRNAs) are small single-stranded non-coding RNAs with the potential as biomarkers for a variety of pathologies. Recent studies found that miRNAs are functional regulators and indicators in the immune processes related to kidney transplantation. For example, Liu et al. reported that downregulated miR-10b could mediate rejection of renal allografts through inhibiting BCL2L11 expression . Jin et al. decoded a novel mechanism between miR-650 and BCL11B for preventing rejection after kidney transplantation . Based on the meta-analysis and F344-Lewis rat kidney transplantation model, Liang et al. explored the diagnostic role and dynamic change of miR-155 in acute rejection . Although the experimental evidences are accumulated, few of these studies focused on the multi-regulatory patterns between miRNAs and mRNAs, and miRNA-mRNA associations are identified without thoroughly weighting miRNA interactions as structural knowledge for biomarker prioritization .
It is widely acknowledged that the interplay between miRNAs and mRNAs forms large-scale miRNA-mRNA regulatory networks and dynamic changes in such network systems provide functional insights in disease prediction . In the era of big biomedical data and artificial intelligence, computer-aided biomarker discovery has become a new frontier to model and decipher complex miRNA-mRNA regulations [7, 8]. Many network topology-based theories are proposed and characterized to extract key miRNA signatures from the noisy background for translational researches. For example, miRNAs with more targets in the network are often found to be important in regulation, and the Hub property is commonly selected for biomarker identification . Compared with this approach, in our previous work miRNAs with significantly strong single-line regulatory power were measured for their possibility as biomarkers . These findings increase the understanding in miRNA biology, however, most of them focus solely on cancer applications, and the framework needs to be expanded by integrating feature parameters from multi-omics network systems such as the downstream protein–protein interaction (PPI) network (PPIN) for cross-level functional survey.
In this study, we proposed a novel bioinformatics approach by characterizing and integrating topological features from both miRNA-mRNA and PPI networks to screen miRNA biomarkers and mechanisms for predicting rejection after allogeneic kidney transplantation. As shown in Fig. 1, in the method a KTR-specific miRNA-mRNA network was first developed to evaluate the binary regulatory modes between miRNAs and mRNAs. Then KTR-specific mRNAs were extracted to investigate their interactions based on PPI data, and a cross-level miRNA-mRNA-PPI network (miR-PPIN) was finally constructed to identify miRNAs and infer key miRNA-mRNA regulations in KTR. Compared with traditional approaches equally treating the contribution of different miRNA targets in the network for model training , in this pipeline the structural importance of target genes was mainly prioritized based on their topological activities in PPIN, and the targets in the centre of PPIN, i.e., high degree, closeness, and betweenness, were classified as CORE factors for miRNA regulatory power measurement.
Materials and methods
Data collection and processing
The miRNA and mRNA datasets were collected and downloaded from Gene Expression Omnibus (GEO) . As illustrated in Table 1, GSE30282 contains a total of 65 post-transplant biopsy samples profiled on Affymetrix Multi-species miRNA-1 Array. In this study AMR, CMR and the normal control samples were selected . In GSE129166, peripheral blood samples and kidney allograft biopsies were chosen for genome-wide gene expression analysis using Affymetrix Human Genome U133 Plus 2.0 Array . To ensure the consistency of sample sources for model training, kidney allograft biopsy samples with clear pathologic characteristics, i.e., AMR or CMR (borderline samples were removed), were analyzed. In addition, GSE115816 with 6 AMR, 4 CMR and 6 stable graft function samples screened by Illumina HiSeq 2500 was used as an independent validation dataset to evaluate and compare the biomarker power of identified miRNAs for KTR prediction and subtyping in blood samples.
The differentially expression analysis were performed on the normalized data to extract differentially expressed (DE) miRNAs and mRNAs between KTR and normal groups. According to the comparison of methods for microarray-based DE-gene identification , the empirical bayes (eBayes) method was recommended for raw p-value calculation  and the Benjamini–Hochberg false discovery rate (FDR) was calculated for multiple testing and p-value adjustment. For the gene related to multiple probes, the probe with the most significant variation was assigned. The cutoff for DE-miRNA and DE-mRNA extraction was set as the adjusted p-value (adj.p-value) < 0.05.
The miRNA-mRNA regulatory data were retrieved from our previous work, where a human global miRNA-mRNA network was constructed based on the integration of miRNA-mRNA pairs from both experimentally validated and computationally predicted databases . For human PPI data, the online tool STRING (v11.0) was applied and PPIs with the score ≥ 0.7 (high confidence) and active interaction sources except text-mining were chosen for network modeling and analysis .
Multi-omics network characterization and miRNA biomarker identification
As described in Fig. 1, an integrated KTR-specific network was constructed by merging interactive data at two omics levels to quantify the power of miRNAs for gene regulation:
First, the DE-miRNAs and DE-mRNAs were mapped onto the human global miRNA-mRNA network, and a KTR-specific miRNA-mRNA network was extracted to measure miRNA regulation at the post-transcriptional level.
Second, it should be acknowledged that genes in the network could influence the downstream biological processes and eventually cause disease phenotypes via interacting with each other. To better prioritize the interactive activity of genes regulated by miRNAs, a KTR-specific PPI network was further developed by integrating KTR-specific mRNAs with human PPI data using STRING online . Here three topological features, i.e., degree centrality (DC), closeness centrality (CC), and betweenness centrality (BC), were calculated to indicate the structural significance of genes in PPIN :
where k is the number of edges that are linked to a given node i. Nodes with high DC are critical in the network.
where N is the total number of nodes in the network, and d(i, j) represents the distance between i and j. This index represents how close the given node is to all the other nodes in the network. In general, the node with high CC holds the optimum position for evaluating the information flow.
where nxy is the number of the shortest paths connecting node x and y, and mixy is the number of the shortest paths connecting node x and y that contain the given node i. Since the interaction of two non-adjacent nodes can be influenced by nodes that lie between them, this index indicates the significance of a given node based on the number of penetrated shortest paths.
In this study, genes whose DC, CC and BC are greater than their average values are defined as CORE set in the network (see Additional file 1), and the KTR-specific miRNA-mRNA network were upgraded by merging PPIN to a cross-level miRNA-mRNA-PPI network (miR-PPIN). In this novel framework, the structural importance of target genes is prioritized according to their contributions in the network, and miRNAs with more CORE targets have the theoretical priority for biomarker discovery.
Based on the above definitions, two feature parameters, i.e., Number of ALL Targeted Genes (NTGA) and Number of Targeted CORE Genes (NTGC), were respectively defined, where the former is numerically equal to the number of genes regulated by miRNAs and the latter is the number of CORE genes targeted by a given miRNA. As hub genes are reported to be functionally important in the network system , miRNAs with high NTGA values have strong power in regulation. Moreover, NTGC index strengthens such regulatory pattern by reasonably weighting the contribution of target genes across multi-omics networks. Hence candidate miRNA biomarkers were computationally identified with the following two steps:
Step 1: miRNAs with significantly high NTGA values (p-value < 0.05, Wilcoxon signed-rank test) were extracted from the KTR-specific network.
Step 2: miRNAs with significantly high NTGC values (p-value < 0.05, Wilcoxon signed-rank test) were selected from those screened in Step 1 as candidate biomarkers for expression and functional validation.
ROC analysis and performance validation
The receiver-operating characteristic (ROC) analysis was performed and compared using the statistical tool MedCalc (v20.009) at two levels, i.e., the identified miRNAs as single predictors and a combined signature that incorporated all identified miRNAs. Here the Logistic regression was applied to predict the value of combined signature based on the expression data of each miRNA variables. The area under ROC curve (AUC) was calculated to quantify the predictive power of biomarkers for KTR diagnosis (KTR vs. nKTR) and subtyping (AMR vs. CMR).
The potential of miRNAs as biomarkers was further validated by literature-based pathogenic survey, and keywords including “kidney/renal transplantation”, “rejection”, “immune response”, etc. were used to investigate the associations of identified miRNAs with KTR through text-mining in PubMed and Web of Science.
Knowledge-guided key regulatory module extraction and functional exploration
The pathogenesis of identified miRNAs was explored based on a three dimensional “miRNA-gene-pathway” paradigm. First, functional targets of biomarker miRNAs were retrieved from KTR-specific miR-PPIN network for Gene Ontology (GO) enrichment at the biological process (BP) level using the Database for Annotation, Visualization and Integrated Discovery (DAVID, v6.8) online . The criterion for significant term selection was set as p-value < 0.05. Then PPIs in KTR-specific miR-PPIN were clustered by the plug-in Molecular Complex Detection (MCODE) of Cytoscape for molecular complex analysis with default thresholds . Finally, significant clusters were identified and key regulatory modules were extracted by integrating BP knowledge and miRNA regulations for translational “genotype-phenotype” etiologic understanding.
Global features of KTR-specific miR-PPIN
In this study, a KTR-specific miR-PPIN network was constructed from the integration of DE-miRNA/mRNA, miRNA-mRNA and PPI data for systems-level characterization of miRNA-mediated interactions in KTR (see Additional file 2). As shown in Fig. 2a, the network contains a total of 949 miRNA-mRNA pairs and 1,447 PPIs among 31 miRNAs and 642 genes that were abnormally expressed in KTR. Based on STRING analysis  as described in Fig. 2b, the genes in the network were significantly enriched in T cell-related BP terms such as NK T cell proliferation, T cell antigen processing and presentation, positive regulation of memory T cell activation, et al., which indicated the specificity of identified network in KTR development [22, 23].
To characterize the global features of miRNAs and genes in the network, the degree distribution was analyzed at miRNA-mRNA and PPI level, respectively. As illustrated in Fig. 2c and d, miRNAs with more targets in miRNA-mRNA network were fewer, and most of the targets tended to locate at border sites in PPIN as the degree of nodes followed the power-law distribution with the slope of -1.34. Based on these findings, the structural importance of target genes was then prioritized according to their degree, closeness and betweenness properties in PPIN, and the highly prioritized targets were classified as CORE gene set to strengthen miRNA regulation (see Materials and methods). Finally, features in these two network systems were collaboratively measured, and the number of CORE targets by miRNAs was calculated in miR-PPIN. As shown in Fig. 2e, there were very few miRNAs that regulated more CORE genes and it could be a potential clue for biomarker discovery.
Biomarker miRNAs for KTR diagnosis and subtyping
In this study, a total of three miRNAs, i.e., miR-145-5p, miR-155-5p, and miR-23b-3p, were identified as candidate biomarkers for KTR prediction based on multi-omics network characterization with the newly defined topological parameters and screening criteria, i.e., significantly high NTGA and NTGC. As illustrated in Table 2 and Fig. 3 respectively, miR-145-5p and miR-23b-3p were downregulated in KTR group, whereas miR-155-5p was overexpressed in KTR especially in CMR.
The ROC analyses were performed both on the training and an independent validation dataset to evaluate the potential of identified miRNAs and the miRNA combination for KTR diagnosis and subtyping. At single miRNA level, as drawn in Fig. 4a, b and e, all the three miRNAs showed strong power (AUC > 0.8, average AUC = 0.898) in differentiating kidney allograft biopsy samples of KTR and non-KTR (nKTR, normal control), and they had comparable ability for KTR subtyping, i.e., AMR and CMR (AUC > 0.75, average AUC = 0.829). In the validation set, the miRNAs were measured using blood samples since blood can be obtained noninvasively compared with kidney tissues. As shown in Fig. 4c, d and e, miR-155-5p and miR-23b-3p tended to outperform on KTR prediction (AUC > 0.75), and miR-145-5p reached the AUC of 0.917 for AMR and CMR subtyping. Compared with single biomarkers, the combination of three miRNAs improved the overall performance in different groups, which indicated the potential of the three miRNAs as a combined signature for KTR management.
It should be noticed that the AUC values of identified miRNAs were highly heterogeneous between tissue and blood groups. In addition to the differences in sample source, the number and size of dataset is another important consideration. Therefore large-sample-based consistency analysis should be conducted for further functional validation.
Literature-based functional validation and comparison
The literature-based survey was conducted to validate the performance of the bioinformatics model and investigate the function of identified miRNAs. Among them, miR-155-5p is a well-studied star related to dysfunctional allograft status in kidney transplantation. For example, Anglicheau et al. demonstrated that the intragraft expression of miRNAs was strongly associated with mRNAs and the disorder in such regulation led to acute rejection and affected kidney allograft function. They found that miR-155-5p (namely miR-155) was overexpressed in KTR biopsy samples, which could serve as a diagnostic biomarker for indication of kidney allograft status . Similar to this result, the overexpression of miR-155-5p was detected in T-cell mediated rejection from a study by Soltaninejad et al., and the diagnostic role of this miRNA was suggested for CMR prediction . To validate the dynamic change of miR-155-5p expression in KTR, Liang et al. established a rat kidney transplantation model and the result showed that the increased level of miR-155-5p in plasma was highly consistent with the degree of rejection development, thus this miRNA could be a useful biomarker for monitoring the functional status of allograft after kidney transplantation . In addition to tissue and blood samples, the upregulation of miR-155-5p was also captured in urine from patients who experienced rejection. Based on qRT-PCR validation, the level of urinary miR-155-5p was significantly increased in patients before and during acute rejection, and this miRNA held the power to be a biomarker for early diagnosis of rejection after kidney transplantation .
The remaining two miRNAs, i.e., miR-145-5p and miR-23b-3p, were also reported to be functional in KTR progression. According to the RT-PCR analysis by Matz et al., the abnormal expression of miR-145-5p distinguished rejection from the normal controls and it was downregulated in patient cohorts with interstitial fibrosis/tubular atrophy (IFTA), which indicated its diagnostic role as an IFTA-specific biomarker in blood for fibrosis-caused etiology study . Meanwhile, this miRNA, as well as miR-23b-3p, could distinguish acute rejection and acute pyelonephritis, which provided underlying clues for rejection identification from inflammatory conditions .
To summarize, miR-155-5p is a reported biomarker for KTR diagnosis and the remaining candidates also have close associations with KTR especially with kidney acute rejection processes, supporting the overall predictive performance of the proposed model. Compared with previous approaches and findings, the expression pattern of identified miRNAs in KTR group was highly consistent. More importantly, in this study a network-based computational strategy was applied to identify miRNAs with biomarker potential both in KTR diagnosis (KTR vs. nKTR) and subtyping (AMR vs. CMR), which provided novel insights in miRNA-medicated KTR pathogenesis.
Key regulatory modules and mechanisms in KTR pathogenesis
It is widely known that immune response is a leading factor affecting the function of allograft after kidney transplantation [29, 30]. In this study immune-related BP terms, e.g., interferon-gamma-mediated signaling pathway, type I interferon signaling pathway, immune response, antigen processing and presentation, etc. were significantly enriched by biomarker miRNA targets including HLA-A, HLA-C, HLA-G, SAMHD1, ICAM1, IRF1/2, etc. in KTR-specific miR-PPIN (see Additional file 3). To explore genetic interactions associated with identified BP terms, the plug-in MCODE of Cytoscape was applied to cluster functional modules for pathogenetic understanding. As illustrated in Additional file 4, a total of 18 clusters were extracted from KTR-specific PPIN based on MCODE analysis and the Cluster 2 with identified biomarker targets was selected for etiologic study.
As shown in Fig. 5, in KTR-specific Cluster 2 the identified targets were interacted and formed as a key module in regulating immune and interferon-related pathways. Herein HLA-A, HLA-C, HLA-G, IRF1/2 and SP100 were regulated by miR-23b-3p and the remaining three were targets of miR-155-5p, respectively. Among them, HLA (human leukocyte antigen) is known to be the expression product of human major histocompatibility complex (MHC). The mismatch of HLA is a risk factor to affect kidney allograft function and to increase the chance of rejection [31,32,33]. For example, Ko et al. investigated the clinical outcome of kidney transplantation with ABO and HLA incompatibilities. The result found that ABO and HLA incompatibilities could increase the probability of infection-medicated rejection and influence the overall survival both for patient and graft . Moreover, HLA donor-specific antibodies could induce inflammation, vessel injury and AMR by binding to vascular endothelial cells of the allograft . As three important members in HLA class I heavy chain paralogues, HLA-A, HLA-C and HLA-G were found to be potentially regulated by miR-23b-3p in this study. Here HLA class I molecules have central functions in the immune system, and they can be detected by cytotoxic T cells via presenting peptides originated from the endoplasmic reticulum lumen. In particular, Janssen et al. found that polymorphisms in donor derived HLA-G had significant effects on acute rejection after kidney transplantation, and 14-bp ins/ins and the + 3142GG genotypes may be protective factors to decrease KTR . Hence the signaling of miR-23b-3p/HLA/immune axis offers novel insights in KTR genesis. In addition to HLA groups, ICAM1 and VCAM1 are well investigated in the function of transplanted kidney. For example, the immunological response after kidney transplantation would increase the endothelial expression of ICAM1 and VCAM1, and the rs5498 ICAM1 polymorphism and the C allele of the rs1041163 in VCAM1 highlighted the rejection risk and long-term allograft failure [36, 37]. Moreover, as an inflammatory protein, ICAM1 plays active roles in kidney chronic rejection, and the ICAM1-based network biomarker accelerated the precision classification of different KTR phenotypes . Last but not least, IRF1 is a transcription factor and it mediated MHC induction and regulated resistance to necrosis in KTR, which would be crucial for the interventions against rejection .
Based on these findings, the putative “miRNA-gene-pathway” axis was inferred at the computational level and more insightful mechanisms in KTR pathogenesis need to be experimentally validated in the future work.
Allograft rejection is a common complication after kidney transplantation. The early detection and precise subtyping of rejection signatures are of clinical significance for personalized therapy. As a class of non-coding RNAs with regulatory role in multiple biological processes, miRNAs are well reported to be potential biomarkers for predicting rejection response in kidney transplantation [40, 41]. However, most of the studies are experimental, which lacks deep insights in the global interactions between miRNAs as well as target genes for systems-level pathogenesis understanding.
In this study, a novel bioinformatics framework was proposed based on integration of microarray data and multi-omics network knowledge for cross-level deciphering of RNA interplay during KTR development. In methodology, a KTR-specific miR-PPIN was constructed and topological characterization was conducted to uncover theoretical hypothesis for miRNA biomarker discovery. Compared with traditional methods that equally rank the contribution of miRNAs in miRNA-mRNA network [42, 43], the network system in this study was expanded by reasonably weighting the interaction among targets of miRNAs based on downstream PPI analysis. Hence in this model, the topological importance of genes targeted by miRNAs was prioritized according to their structural characteristics in PPIN. Two feature parameters, i.e., NTGA and NTGC, were defined for biomarker identification. Finally, a total of three miRNAs, i.e., miR-145-5p, miR-155-5p, and miR-23b-3p, were screened as candidate biomarkers for monitoring the rejection after kidney transplantation.
Among the three candidates, miR-155-5p was significantly overexpressed in rejection groups, whereas the remaining two were downregulated in KTR process. The ROC analysis showed that all the three miRNAs distinguished KTR from the normal controls with the AUC > 0.8 in kidney allograft biopsy samples, and they achieved comparable performance on AMR and CMR subtyping (AUC > 0.75) as well. Meanwhile, miR-155-5p/miR-23b-3p (AUC > 0.75) and miR-145-5p (AUC > 0.9) were respectively found to be powerful for KTR prediction and subtyping in blood, indicating the noninvasive diagnostic potential of identified miRNAs. In addition to single biomarkers, the three miRNAs could serve as a combined signature to improve the overall performance for KTR management.
According to literature reports, miR-155-5p had found to be a diagnostic biomarker and a prognostic factor for predicting allograft status and rejection [4, 26, 44]. The remaining miRNAs were also found to be dysfunctional during KTR development. Compared with existing findings, in this study the miRNAs showed the predictive power both in KTR diagnosis and subtyping. Moreover, miR-23b-3p and miR-155-5p were inferred to regulate immune and interferon-medicated pathways via targeting the module of HLA-related and ICAM1/VCAM1 genes respectively, which highlighted the novel putative mechanisms for KTR understanding.
The limitations in this study should be carefully considered although meaningful results were obtained and analyzed. First, the importance and contribution of miRNAs and genes in the network were prioritized based on their structural characteristics, the functional significance of genes in immune processes related to rejection responses needs to be investigated and quantified to improve the sensitivity of miRNAs in KTR management. Second, in addition to KTR and the normal controls, samples from other conditional phenotypes such as infection-mediated or drug toxicity-induced delated graft function (DGF) should be included and compared to test the biomarker specificity of identified miRNAs in rejection since similar clinical symptoms, e.g., decreased amount of urine and increased level of serum creatinine, often occurred in the early stage both of KTR and DGF [45, 46]. Finally, it is potentially important that the biomarkers originally identified from tissue biopsies were able to be validated in blood samples as blood could be easily obtained for clinical translation. However, the ROC result from tissue and blood group tended to be heterogeneous in this study. Although many studies reported that the concordance rate for detection of gene variant and expression may be affected by the type and source of sample data [47,48,49], it should be admitted that the blood sample set used for comparison was relatively small. Hence the biomarker significance of identified miRNAs in blood needs to be further explored with the accumulation of enough public data to decrease the risk of model overfitting. On the other hand, large-sample-based multi-center experimental validation and pathogenetic survey are expected to be performed for future translational applications.
In this study, a total of three miRNAs, i.e., miR-145-5p, miR-155-5p, and miR-23b-3p, were screened as candidate biomarkers for diagnosis and subtyping of allograft rejection after kidney transplantation based on an integrated bioinformatics model with multi-omics network characterization. The potential regulatory mechanisms of HLA-related genes, immune signaling pathways and identified miRNAs were computationally discovered for pathogenic understanding. In the future work, molecular experiments using human samples will be performed for further clinical validation.
Availability of data and materials
The data generated or analyzed during this study are available from the corresponding authors upon reasonable request.
Kidney transplant rejection
Protein–protein interaction network
Gene Expression Omnibus
False discovery rate
- NTGA :
Number of ALL Targeted Genes
- NTGC :
Number of Targeted CORE Genes
Area under ROC curve
The Database for Annotation, Visualization and Integrated Discovery
Molecular Complex Detection
Interstitial fibrosis/tubular atrophy
Human leukocyte antigen
Major histocompatibility complex
Delated graft function
Cooper JE. Evaluation and treatment of acute rejection in kidney allografts. Clin J Am Soc Nephrol. 2020;15:430–8.
Liu X, Dong C, Jiang Z, Wu WK, Chan MT, Zhang J, Li H, Qin K, Sun X. MicroRNA-10b downregulation mediates acute rejection of renal allografts by derepressing BCL2L11. Exp Cell Res. 2015;333:155–63.
Jin P, Chen H, Xie J, Zhou C, Zhu X. Essential role of microRNA-650 in the regulation of B-cell CLL/lymphoma 11B gene expression following transplantation: a novel mechanism behind the acute rejection of renal allografts. Int J Mol Med. 2017;40:1840–50.
Liang J, Tang Y, Liu Z, Wang X, Tang L, Zou Z, Zhou C, Wu K, Zhang F, Zhang F, Lu Y. Increased expression of miR-155 correlates with abnormal allograft status in solid organ transplant patients and rat kidney transplantation model. Life Sci. 2019;227:51–7.
Lin Y, Qian F, Shen L, Chen F, Chen J, Shen B. Computer-aided biomarker discovery for precision medicine: data resources, models and applications. Brief Bioinform. 2019;20:952–75.
Hu B, Ma X, Fu P, Sun Q, Tang W, Sun H, Yang Z, Yu M, Zhou J, Fan J, Xu Y. miRNA-mRNA regulatory network and factors associated with prediction of prognosis in hepatocellular carcinoma. Genomics Proteomics Bioinformatics 2021.
Shen L, Bai J, Wang J, Shen B. The fourth scientific discovery paradigm for precision medicine and healthcare: challenges ahead. Precision Clinical Medicine 2021.
Shen B, Lin Y, Bi C, Zhou S, Bai Z, Zheng G, Zhou J. Translational informatics for Parkinson’s disease: from big biomedical data to small actionable alterations. Genomics Proteomics Bioinformatics. 2019;17:415–29.
Chen B, Hua Z, Qin X, Li Z. Integrated microarray to identify the hub miRNAs and constructed miRNA-mRNA network in neuroblastoma via bioinformatics analysis. Neurochem Res. 2021;46:197–212.
Lin Y, Wu W, Sun Z, Shen L, Shen B. MiRNA-BD: an evidence-based bioinformatics model and software tool for microRNA biomarker discovery. RNA Biol. 2018;15:1093–105.
Lin Y, Miao Z, Zhang X, Wei X, Hou J, Huang Y, Shen B. Identification of key MicroRNAs and mechanisms in prostate cancer evolution based on biomarker prioritization model and carcinogenic survey. Front Genet. 2020;11:596826.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.
Wilflingseder J, Regele H, Perco P, Kainz A, Soleiman A, Muhlbacher F, Mayer B, Oberbauer R. miRNA profiling discriminates types of rejection and injury in human renal allografts. Transplantation. 2013;95:835–41.
Van Loon E, Gazut S, Yazdani S, Lerut E, de Loor H, Coemans M, Noel LH, Thorrez L, Van Lommel L, Schuit F, et al. Development and validation of a peripheral blood mRNA assay for the assessment of antibody-mediated kidney allograft rejection: a multicentre, prospective study. EBioMedicine. 2019;46:463–72.
Jeffery IB, Higgins DG, Culhane AC. Comparison and evaluation of methods for generating differentially expressed gene lists from microarray data. BMC Bioinformatics. 2006;7:359.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13.
Lin Y, Yuan X, Shen B. Network-based biomedical data analysis. Adv Exp Med Biol. 2016;939:309–32.
Zhang W, Landback P, Gschwend AR, Shen B, Long M. New genes drive the evolution of gene interaction networks in the human and mouse genomes. Genome Biol. 2015;16:202.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.
Morris AB, Pinelli DF, Liu D, Wagener M, Ford ML. Memory T cell-mediated rejection is mitigated by FcgammaRIIB expression on CD8(+) T cells. Am J Transplant. 2020;20:2206–15.
Hughes AD, Zhao D, Dai H, Abou-Daya KI, Tieu R, Rammal R, Williams AL, Landsittel DP, Shlomchik WD, Morelli AE, et al. Cross-dressed dendritic cells sustain effector T cell responses in islet and kidney allografts. J Clin Invest. 2020;130:287–94.
Anglicheau D, Sharma VK, Ding R, Hummel A, Snopkowski C, Dadhania D, Seshan SV, Suthanthiran M. MicroRNA expression profiles predictive of human renal allograft status. Proc Natl Acad Sci U S A. 2009;106:5330–5.
Soltaninejad E, Nicknam MH, Nafar M, Ahmadpoor P, Pourrezagholi F, Sharbafi MH, Hosseinzadeh M, Foroughi F, Yekaninejad MS, Bahrami T, et al. Differential expression of microRNAs in renal transplant patients with acute T-cell mediated rejection. Transpl Immunol. 2015;33:1–6.
Millan O, Budde K, Sommerer C, Aliart I, Rissling O, Bardaji B, Matz M, Zeier M, Silva I, Guirado L, Brunet M. Urinary miR-155-5p and CXCL10 as prognostic and predictive biomarkers of rejection, graft outcome and treatment response in kidney transplantation. Br J Clin Pharmacol. 2017;83:2636–50.
Matz M, Heinrich F, Lorkowski C, Wu K, Klotsche J, Zhang Q, Lachmann N, Durek P, Budde K, Mashreghi MF. MicroRNA regulation in blood cells of renal transplanted patients with interstitial fibrosis/tubular atrophy and antibody-mediated rejection. PLoS ONE. 2018;13:e0201925.
Oghumu S, Bracewell A, Nori U, Maclean KH, Balada-Lasat JM, Brodsky S, Pelletier R, Henry M, Satoskar AR, Nadasdy T, Satoskar AA. Acute pyelonephritis in renal allografts: a new role for microRNAs? Transplantation. 2014;97:559–68.
Peng B, Gong H, Tian H, Zhuang Q, Li J, Cheng K, Ming Y. The study of the association between immune monitoring and pneumonia in kidney transplant recipients through machine learning models. J Transl Med. 2020;18:370.
Jouve T, Laheurte C, Noble J, Weinhard J, Daligault M, Renaudin A, Naciri Bennani H, Masson D, Gravelin E, Bugnazet M, et al. Immune responses following tocilizumab therapy to desensitize HLA-sensitized kidney transplant candidates. Am J Transplant. 2021. https://doi.org/10.1111/ajt.16709.
Ko Y, Kim JY, Kim SH, Kim DH, Lim SJ, Shin S, Kim YH, Jung JH, Park SK, Kwon H, Han DJ. Acute rejection and infectious complications in ABO- and HLA-incompatible kidney transplantations. Ann Transplant. 2020;25:e927420.
Meneghini M, Crespo E, Niemann M, Torija A, Lloberas N, Pernin V, Fontova P, Melilli E, Fava A, Montero N, et al. Donor/Recipient HLA molecular mismatch scores predict primary humoral and cellular alloimmunity in kidney transplantation. Front Immunol. 2020;11:623276.
Abuhelaiqa E, Friedlander R, Aull M, Putheti P, Sharma V, Suthanthiran M, Dadhania D. Acute rejection, kidney allograft function, and graft survival in patients with circulating pre-transplant IgG antibodies directed against donor HLA-A, -B, or -C locus determined antigens. Clin Transpl. 2016;32:83–91.
Wei X, Valenzuela NM, Rossetti M, Sosa RA, Nevarez-Mejia J, Fishbein GA, Mulder A, Dhar J, Keslar KS, Baldwin WM 3rd, et al. Antibody-induced vascular inflammation skews infiltrating macrophages to a novel remodeling phenotype in a model of transplant rejection. Am J Transplant. 2020;20:2686–702.
Janssen M, Thaiss F, Nashan B, Koch M, Thude H. Donor derived HLA-G polymorphisms have a significant impact on acute rejection in kidney transplantation. Hum Immunol. 2019;80:176–83.
Kloda K, Domanski L, Pawlik A, Kurzawski M, Safranow K, Ciechanowski K. Effect of the ICAM1 and VCAM1 gene polymorphisms on delayed graft function and acute kidney allograft rejection. Ann Transplant. 2010;15:15–20.
Kloda K, Domanski L, Pawlik A, Safranow K, Ciechanowski K. The impact of ICAM1 and VCAM1 gene polymorphisms on long-term renal transplant function and recipient outcomes. Ann Transplant. 2013;18:231–7.
Zhu D, Liu Z, Pan Z, Qian M, Wang L, Zhu T, Xue Y, Wu D. A new method for classifying different phenotypes of kidney transplantation. Cell Biol Toxicol. 2016;32:323–32.
Afrouzian M, Ramassar V, Urmson J, Zhu LF, Halloran PF. Transcription factor IRF-1 in kidney transplants mediates resistance to graft necrosis during rejection. J Am Soc Nephrol. 2002;13:1199–209.
Bijkerk R, Florijn BW, Khairoun M, Duijs J, Ocak G, de Vries APJ, Schaapherder AF, Mallat MJK, de Fijter JW, Rabelink TJ, et al. Acute rejection after kidney transplantation associates with circulating MicroRNAs and vascular injury. Transplant Direct. 2017;3:e174.
Matz M, Lorkowski C, Fabritius K, Durek P, Wu K, Rudolph B, Neumayer HH, Mashreghi MF, Budde K. Free microRNA levels in plasma distinguish T-cell mediated rejection from stable graft function after kidney transplantation. Transpl Immunol. 2016;39:52–9.
Ali A, Hadlich F, Abbas MW, Iqbal MA, Tesfaye D, Bouma GJ, Winger QA, Ponsuksili S. MicroRNA-mRNA networks in pregnancy complications: a comprehensive downstream analysis of potential biomarkers. Int J Mol Sci. 2021;22:2313.
Qi X, Yu C, Wang Y, Lin Y, Shen B. Network vulnerability-based and knowledge-guided identification of microRNA biomarkers indicating platinum resistance in high-grade serous ovarian cancer. Clin Transl Med. 2019;8:28.
Van Aelst LN, Summer G, Li S, Gupta SK, Heggermont W, De Vusser K, Carai P, Naesens M, Van Cleemput J, Van de Werf F, et al. RNA profiling in human and murine transplanted hearts: identification and validation of therapeutic targets for acute cardiac and renal allograft rejection. Am J Transplant. 2016;16:99–110.
Lim WH, Johnson DW, Teixeira-Pinto A, Wong G. Association between duration of delayed graft function, acute rejection, and allograft outcome after deceased donor kidney transplantation. Transplantation. 2019;103:412–9.
Wu WK, Famure O, Li Y, Kim SJ. Delayed graft function and the risk of acute rejection in the modern era of kidney transplantation. Kidney Int. 2015;88:851–8.
Xie F, Zhang Y, Mao X, Zheng X, Han-Zhang H, Ye J, Zhao R, Zhang X, Sun J. Comparison of genetic profiles among primary lung tumor, metastatic lymph nodes and circulating tumor DNA in treatment-naive advanced non-squamous non-small cell lung cancer patients. Lung Cancer. 2018;121:54–60.
Lim SB, Yeo T, Lee WD, Bhagat AAS, Tan SJ, Tan DSW, Lim WT, Lim CT. Addressing cellular heterogeneity in tumor and circulation for refined prognostication. Proc Natl Acad Sci U S A. 2019;116:17957–62.
Chen N, Wen S, Sun X, Fang Q, Huang L, Liu S, Li W, Qiu M. Elevated mitochondrial DNA copy number in peripheral blood and tissue predict the opposite outcome of cancer: a meta-analysis. Sci Rep. 2016;6:37404.
The authors gratefully thank the editors and reviewers for their constructive suggestions to improve this manuscript. They also thank Prof. Jiajia Chen and Dr. Xin Qi at School of Chemistry and Life Sciences, Suzhou University of Science and Technology, China, for their kind help and support in manuscript revision.
This research was funded by the National Natural Science Foundation of China (Grant Number 32070671), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant Number 20KJB180010), and the regional innovation cooperation between Sichuan and Guangxi Provinces (Grant Number 2020YFQ0019).
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lin, Y., Wang, L., Ge, W. et al. Multi-omics network characterization reveals novel microRNA biomarkers and mechanisms for diagnosis and subtyping of kidney transplant rejection. J Transl Med 19, 346 (2021). https://doi.org/10.1186/s12967-021-03025-8