Next Article in Journal
Carbonized Lysine-Nanogels Protect against Infectious Bronchitis Virus
Next Article in Special Issue
Integrative Analyses of Circulating Small RNAs and Kidney Graft Transcriptome in Transplant Glomerulopathy
Previous Article in Journal
Remodeling of Arterial Tone Regulation in Postnatal Development: Focus on Smooth Muscle Cell Potassium Channels
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 11
10.3390/ijms22115414
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

An Integrated Transcriptomic Approach to Identify Molecular Markers of Calcineurin Inhibitor Nephrotoxicity in Pediatric Kidney Transplant Recipients

by
Erika T. Rhone
1,†,
Elissa Bardhi
2,†,
Sai Vineela Bontha
3,†,
Patrick D. Walker
4,
Jorge A. Almenara
5,
Catherine I. Dumur
5,
Helen Cathro
6,
Daniel Maluf
7,8 and
Valeria Mas
2,*
1
Department of Pediatrics, Eastern Virginia Medical School, Norfolk, VA 23507, USA
2
Surgical Sciences Division, Department of Surgery, School of Medicine, University of Maryland, Baltimore, MD 21201, USA
3
Cardiovascular Research Center, University of Virginia, Charlottesville, VA 22908, USA
4
Arkana Laboratories, Little Rock, AR 72211, USA
5
Aurora Diagnostics–Sonic Healthcare, Bernhardt Laboratories, Jacksonville, FL 32216, USA
6
Department of Pathology, University of Virginia, Charlottesville, VA 22908, USA
7
Division of Transplantation, Department of Surgery, School of Medicine, University of Maryland, Baltimore, MD 21201, USA
8
Program of Transplantation, School of Medicine, University of Maryland, Baltimore, MD 21201, USA
*
Author to whom correspondence should be addressed.
†
Authors contributed equally.
Int. J. Mol. Sci. 2021, 22(11), 5414; https://doi.org/10.3390/ijms22115414
Submission received: 11 April 2021 / Revised: 11 May 2021 / Accepted: 17 May 2021 / Published: 21 May 2021
(This article belongs to the Special Issue RNA Modification Dynamics in Cellular Dysregulation and Clinical Diseases)
Browse Figures
Versions Notes

Abstract

:
Calcineurin inhibitors are highly efficacious immunosuppressive agents used in pediatric kidney transplantation. However, calcineurin inhibitor nephrotoxicity (CNIT) has been associated with the development of chronic renal allograft dysfunction and decreased graft survival. This study evaluated 37 formalin-fixed paraffin-embedded biopsies from pediatric kidney transplant recipients using gene expression profiling. Normal allograft samples (n = 12) served as negative controls and were compared to biopsies exhibiting CNIT (n = 11). The remaining samples served as positive controls to validate CNIT marker specificity and were characterized by other common causes of graft failure such as acute rejection (n = 7) and interstitial fibrosis/tubular atrophy (n = 7). MiRNA profiles served as the platform for data integration. Oxidative phosphorylation and mitochondrial dysfunction were the top molecular pathways associated with overexpressed genes in CNIT samples. Decreased ATP synthesis was identified as a significant biological function in CNIT, while key toxicology pathways included NRF2-mediated oxidative stress response and increased permeability transition of mitochondria. An integrative analysis demonstrated a panel of 13 significant miRNAs and their 33 CNIT-specific gene targets involved with mitochondrial activity and function. We also identified a candidate panel of miRNAs/genes, which may serve as future molecular markers for CNIT diagnosis as well as potential therapeutic targets.

1. Introduction

Kidney transplantation (KT) remains the preferred therapy for children with end-stage renal disease [1,2,3,4]. Yet despite significant advancements in short-term outcomes [5,6,7], long-term outcomes remain suboptimal. Currently, patient survival is approximately 90% at 10 years post-transplant, but only 50–60% of allografts survive this long [7,8,9]. The mismatch between patient and graft longevity presents challenges for both pediatric patients and transplant centers, since repeat transplant not only increases the likelihood of morbidity and mortality for the child but also reduces the supply of donor organs available to other transplant candidates. Currently, 13.4% of pediatric kidney transplant waitlist candidates await a second transplant [8].
Therefore, strategies to reduce premature graft loss are of paramount importance. Late graft failure is most frequently caused by chronic renal allograft dysfunction (CRAD), which refers to the final common pathway of aggregate immunologic and non-immunologic insults resulting in the insidious loss of graft function and manifested histologically by fibrosis [10,11,12,13]. Calcineurin inhibitor nephrotoxicity (CNIT) is considered one of the primary non-immunologic factors contributing to CRAD [13,14]. In 2003, Nankivell et al. evaluated the natural history of CRAD using longitudinal protocol biopsies and found nearly universal histologic evidence of CNIT at 10 years post-KT [15]. Despite their known nephrotoxicity, calcineurin inhibitors (CNI) remain the cornerstone of maintenance immunosuppressive protocols in both adult and pediatric KT recipients [16,17].
Currently, the diagnosis of CNIT requires a kidney biopsy. However, the histologic features used to identify CNIT include permanent architectural disruption, limiting the utility of this tool to modify CNI therapy prior to the onset of irreversible damage [18,19,20,21]. Additionally, concerns regarding the specificity of CNIT histology, biopsy sampling error, and inter-observer agreement among pathologists further tarnish this gold standard and highlight the need for sensitive molecular makers of CNIT [22,23,24,25].
To date, the development of CNIT biomarkers has been hampered by an incomplete understanding of its underlying molecular mechanisms. Current reliance on late and insensitive markers such as serum creatinine precludes the identification of CNIT prior to the onset of irreversible damage. The lack of surrogate markers highlights a knowledge gap regarding the underlying molecular mechanisms of CNIT. Therefore, we sought to characterize a molecular phenotype of CNIT using an integrative approach. We hypothesized that the interplay of differentially expressed miRNAs and target mRNAs would provide the basis for a unique molecular signature for accurate pediatric CNIT diagnosis. MiRNAs, which are non-coding, endogenous single-stranded RNA that post-transcriptionally repress gene expression, have emerged as a promising class of biomarkers utilized in molecular diagnostics and have garnered interest as potential therapeutic tools [26,27]. With the utilization of miRNA and gene expression tissue profiles, the specific aim of this study is to identify (1) the critical gene pathways associated with CNIT development in pediatric KT recipients and (2) key miRNA–mRNA interactions using data integration between these two molecular layers of regulation.

2. Results

2.1. Molecular Profiling Using Microarrays

For the gene expression analysis, 1837 probe sets were differentially expressed between CNIT and Normal samples (Table S1). This corresponded to 1483 distinct mapped genes with 714 (48%) upregulated and 769 (52%) downregulated. Two additional comparisons for marker specificity were also executed: acute rejection (AR) vs. Normal yielded 678 differentially expressed genes (60% upregulated) and interstitial fibrosis/tubular atrophy (IFTA) vs. Normal demonstrated 285 genes (31% upregulated). Figure 1 summarizes the overall study design and flow.
The miRNA microarray analysis demonstrated 118 miRNAs were differentially expressed (Table S2) between Normal and CNIT samples, with a majority being upregulated (72%). A supervised hierarchical clustering was performed using miRNA expression data, demonstrating the correct grouping of CNIT and Normal cohorts (Figure 2).

2.2. Canonical Pathways and Upstream Regulators

The top gene pathways in the comparison between CNIT and Normal samples included oxidative phosphorylation (p = 3.52 × 10−34), mitochondrial dysfunction (p = 2.90 × 10−32), EIF2 signaling (p = 2.8 × 10−16), protein ubiquitin pathway (p = 2.04 × 10−7), and aryl hydrocarbon receptor signaling (p = 6.71 × 10−7) (Figure S1). Oxidative phosphorylation was the topmost affected pathway in this dataset, with an overlap of 54 (49.5%) out of 109 molecules and a predicted inhibition of NAD+, FAD+, and ATP (Figure S2). The 54 dysregulated genes are listed in Table S3.
Mitochondrial dysfunction was the next most significant canonical pathway. The genes associated with complexes I-V of the electron transport chain were downregulated, with predicted inhibition of ATP production (Figure 3). There was predicted activation of Caspases 3, 8, 9, as well as cytochrome C, in its association with activation of apoptosis. Predicted directionality of EIF2 signaling indicated strong inhibition in our data (activation z-score of −3.9).
Analysis of the mostly highly significant putative upstream regulators was based upon an activation z-score and statistically significant ‘overlap p-values’ in IPA (p < 0.01). RICTOR (RPTOR independent companion of MTOR complex 2) was identified to be the top upstream regulator in the CNIT vs. Normal dataset (p = 4.12 × 10−45), with a predicted activated state (z-score 9.301). There were 97 downstream target molecules of RICTOR present in the dataset, of which 48 (49.5%) were members of the mitochondrial dysfunction/oxidative phosphorylation pathways. Additional upstream regulator analysis demonstrated a mechanistic network based on downstream gene expression patterns, which identified a connection between the activation of RICTOR and transcription factor FOXO1 inhibition. Other activated upstream regulators with the highest scores include kinases such as MAP4K4 (z-score = 5.298, p = 3.34 × 10−6) and GCK (z-score = 3.162, p = 2.94 × 10−6), the transcription regulators KDM5A (z-score = 5.298, p = 2.37 × 10−10) and PML (z-score = 2.612, p = 1.07 × 10−3), as well as mature miR-124-3p (z-score = 2.220, p = 1.87 × 10−2). VEGF was also identified as an activated top upstream regulator (z-score = 3.358, p = 2.85 × 10−2). Top upstream regulators predicted to be inhibited were comprised mainly of transcription regulators (HNF4A, z-score= −3.573, p = 2.23 × 10−23; PPARGC1A, z-score= −4.807, p = 5.41 × 10−7; TP53, z-score= −2.783, p = 1.34 × 10−18; MYCN, z-score= −4.121, p = 2.84 × 10−15).

2.3. Biological Characterization

The top molecular and cellular functions of the CNIT vs. Normal datasets pertained to cell death and survival (p-value range: 5.72 × 10−4–1.69 × 10−28), cellular growth and proliferation (p-value range: 5.52 × 10−4–2.66 × 10−22), protein synthesis (p-value range: 4.16 × 10−4–9.10 × 10−16), and protein degradation (p-value range: 1.41× 10−4–6.15 × 10−16). Functions pertaining to protein synthesis had negative activation z-scores (specifically, protein metabolism, translation, synthesis, expression). When filtered by activation z-score, free radical scavenging (specifically, the production of reactive oxygen species) was one of the highest functions (z-score = 3.263), along with cell movement (z-score = 2.593) and the differentiation of cells (z-score = 2.628). Inhibited functions were highly overexpressed in the CNIT samples related to nucleic acid metabolism and energy production (synthesis of ATP, z-score = −2.142; metabolism of nucleoside triphosphate, z-score −2.142, and synthesis of purine ribonucleotide, z-score= −2.0147).
Functional analysis with a particular focus on nephrotoxicity yielded the following results. Among the most significant renal “tox” functions included renal necrosis/cell death (p-value range: 5.98 × 10−1–2.23× 10−4, 62 dataset molecules), renal damage (p-value range: 2.96 × 10−1–4.58× 10−3), and renal tubule injury (p-value range 2.96 × 10−1–1.12 × 10−2). Specific functions within these categories that were especially activated included ‘apoptosis of kidney cells’ (z-score = 2.052), involving upregulation of TGFB1, TLR4, CSF1, and GFPT1 genes, in addition to ‘fibrosis of kidney’ (z-score = 2.178), with associated upregulation of TGFB1, SMAD4, TLR4, and IL2RA genes.

2.4. Marker Specificity, Unique Genes Associated with CNIT, and Comparison Analyses

In a comparison of IFTA vs. Normal samples, a majority of the top canonical pathways were related to metabolism: fatty acid ß-oxidation I (p = 1.5 × 10−7), TCA Cycle II (Eukaryotic) (p = 2.05 × 10−4), Tryptophan Degradation III (Eukaryotic) (p = 2.43 × 10−4), and D-glucuronate Degradation I (p = 5.02 × 10−4), in addition to mitochondrial dysfunction (p = 8.96 × 10−10). An analysis of AR vs. Normal differentially expressed genes demonstrated a predominance of immune-related canonical pathways, including Natural Killer Cell Signaling (p = 1.16 × 10−18), Th1 and Th2 Activation Pathway (p = 1.6 × 10−18), and T Cell Receptor Signaling (p = 2.82 × 10−11).
Then, a comparison analysis evaluating CNIT vs. Normal, IFTA vs. Normal, and AR vs. Normal was assessed, with a particular emphasis on canonical pathways among these three comparisons in an effort to identify distinguishing expression patterns. When classified by statistical significance, we found mitochondrial dysfunction, fatty acid ß-oxidation I, and oxidative phosphorylation to be the top pathways. Mitochondrial dysfunction and oxidative phosphorylation were the most highly enriched among CNIT vs. Normal samples (p = 2.92 × 10−32 and p = 3.52 × 10−43, respectively), while the fatty acid B-oxidation I pathway (p = 1.5 × 10−7) was predominantly expressed among IFTA samples. By activation z-score, Fc receptor-mediated phagocytosis in macrophages and monocytes, phospholipase C signaling, integrin signaling, and EIF2 signaling were all canonical pathways that were distinctly inhibited in the CNIT vs. Normal comparison.
Then, differentially expressed genes for each of the three comparisons were analyzed for overlapping (common) genes. By eliminating overlapping genes, a list of unique genes was identified for each of the conditions. A total of 66 genes were common among all three conditions. There were 418 (61.7%) unique AR genes and 107 (37.5%) unique IFTA genes. The list of 1245 (84.0%) unique CNIT genes was analyzed separately in IPA and was used for consequent data integration analyses with miRNA microarrays.
Analysis of CNIT-specific genes demonstrated a preservation of the top canonical pathways seen previously in the original CNIT vs. Normal data, including mitochondrial dysfunction, oxidative phosphorylation, and EIF2 signaling pathways. Additional gene pathways highly expressed among these CNIT genes included regulation of eIF4 and p70S6K signaling and remodeling of epithelial adherens junctions. RICTOR also remained the top upstream regulator (p = 1.01 × 10−26), while HNF4A, TP53, MYCN, and MAPT continued to be among the list of the top ten significant upstream regulators. TGFB1 and LONP1, a gene that encodes a mitochondrial matrix peptidase, were also among this list. Evaluation of the top tox lists in IPA identified mitochondrial dysfunction, renal necrosis/cell death, NRF-mediated oxidative stress response, and aryl hydrocarbon receptor signaling as overexpressed gene sets in CNIT samples.

2.5. Integrative Analysis (miRNA–mRNA Interactions)

Data integration was performed using unique CNIT genes and differentially expressed miRNAs between CNIT and Normal groups. This generated a filtered list (experimentally observed and with appropriate expression directionality) of 13 miRNAs and their 33 mRNA targets. An integrated network of each of these 13 miRNAs and their downstream up- and down-regulated genes are represented in Figure 4.
Then, these 13 miRNAs and 33 gene targets were further analyzed for biologic significance. First, the genes were examined for their role in canonical pathways. SOD2 and APP were genes participating in mitochondrial dysfunction, while AGO4, CCND1, and MYC were identified as being involved with EIF2 signaling. The SOD2 gene is also included in superoxide radical degradation and NRF2-mediated oxidative stress response pathways. Then, specific tox lists were evaluated for both miRNA:mRNA pairs, with SOD2, APP, CCND1, miR-16-5p, and MYC being associated with the decreased transmembrane potential of mitochondria. APP was additionally identified as a gene related to increasing the permeability transition of the mitochondrial membrane. With particular attention to nephrotoxicity tox lists, APP, GFPT1, GRB10, MYC, SOD2, and TLR4 were identified as genes related to renal necrosis/cell death and apoptosis. Based on a functional analysis for nephrotoxicity, the TLR4 gene was associated with ischemia–reperfusion injury and tubular injury, while GFPT1, TLR4, GFPT1, and SMAD4 genes were identified to be involved with glomerular injury. Other functions included renal fibrosis (SMAD4), renal cell proliferation (PNN), renal cell viability (APP), and hyperplasia (MYC). Using QPCR, we validated the increase in expression of SOD2 (p = 0.0161), CCND1 (p = 0.0113), GRB10 (p = 0.0653), and TLR4 (p = 0.0577) in CNIT samples when compared to Normal (Figure 5). GRB10 and TLR4 were not statistically significant, although a trend in significance was observed. This is likely due to the relatively small sample size and supports the need of larger studies for further validation of findings.

3. Discussion

This study presents a novel integrative approach to evaluate the molecular signature of CNIT in pediatric KT recipients. This study is the first to examine tissue profiles of CNIT in an exclusively pediatric KT cohort. CNIT is of particular clinical significance given its contributions to CRAD and associations with graft longevity.
Previously, our group reported the molecular profile of CNIT in adult KT recipients [28]. There is reason to believe that age-related and developmental differences in host immune system responses, drug handling (both pharmacokinetics and pharmacodynamics), primary kidney disease, and presence of concurrent cardiovascular co-morbidities may impact gene expression [29,30,31,32,33,34,35,36]. For instance, KT recipients less than 5–6 years of age often require higher CNI doses than older patients, which were potentially related to CYP3A maturation and activity [37,38,39,40]. Other reports have demonstrated more intense immune responsiveness among pediatric recipients, suggesting that the immunological milieu for the allograft may be affected by age [36,41,42]. Furthermore, children are prioritized to receive deceased donor kidneys with lower kidney donor profile index (KDPI) scores (<35%), and thus, the quality of the graft may influence gene expression [8]. Furthermore, the need for sensitive markers of allograft injury may be greater among children given that adult-sized grafts may mask damage that is not accompanied by an increase in serum creatinine [41]. While most biomarker research initially occurs in adult populations and then is extrapolated to children, this study demonstrates the value of primary pediatric biomarker discovery in transplantation.
The focus of this study was to identify molecular pathways associated with CNIT among pediatric KT recipients. We found mitochondrial dysfunction and oxidative phosphorylation to be the two central gene pathways in kidneys with CNIT, even following the elimination of commonly differentially expressed genes among IFTA and AR. Specifically, ATP depletion, increased membrane permeability, and the production of reactive oxygen species emerged as key biological processes enriched in CNIT samples. Based on in silico modeling, inhibition of genes involved in complexes I–IV of the electron transport chain predicted the downstream effect of apoptosis activation via caspases 3, 8, and 9 and cytochrome C. EIF2 signaling was also robustly inhibited in CNIT-specific genes. Among the functions of the EIF2 pathway is to serve as an mRNA translational checkpoint as part of an integrated cellular stress response [43]. Upon exposure to environmental stress, the phosphorylation of EIF2 precludes the formation of the 43S pre-initiation complex required for translation, which is a homeostatic mechanism to down-regulate protein synthesis in cellular injury [44]. Here, the inhibited expression pattern of EIF2 suggests a global reduction in protein synthesis.
Our upstream regulator analysis identified RICTOR as the top activated regulator among CNIT-specific genes. RICTOR, as a subunit of MTORC2, has been shown to have a significant role in actin cytoskeleton organization [45,46,47]. Recently, RICTOR has been implicated in renal fibrosis development, specifically in association with TGFB1-induced fibroblast activation and epithelial–mesenchymal transition [48,49]. Previous studies have highlighted the biologic connection between RICTOR/MTORC2 and mitochondrial activity [50]. MTORC2 has been localized to both endoplasmic reticulum and mitochondria, but it has a functional relationship to mitochondria through Akt phosphorylation [50,51,52]. The RICTOR/MTORC2 complex phosphorylates Akt, a serine/threonine kinase involved in cell death and cell cycle progression [52,53,54]. Both the inhibition and activation of Akt have been associated with apoptosis and increased vulnerability to oxidative stress [55,56]. Our results demonstrate that RICTOR may be an important upstream regulator involved in the cross-talk among these pathways, with roles in both mitochondrial dysfunction and fibrosis development.
Although it was not our primary objective to characterize the molecular profiles of AR and IFTA samples, we did observe distinct gene expression patterns within these groups. Unsurprisingly, the principal gene pathways in AR were related to immunity and inflammation, including activation of Th1/Th2 alloresponse, iCOS signaling in T Helper Cells, and NK cell signaling. Predicted significant upstream regulators were cytokines (IL-2, IL-15, and interferon-gamma). IFTA samples were distinguished by gene pathways and molecular functions associated with decreased metabolism (specifically, lipid metabolism) and energy production, which is largely in concordance with previous reports [56,57,58].
There has been an emerging focus on mitochondria as key arbitrators of both acute kidney injury and chronic kidney diseases [59,60,61]. With regard to CNIT, the association between mitochondrial dysfunction and CNIs has been previously described by others [18,62,63,64,65,66,67,68]. While it remains unknown whether the observed mitochondrial dysfunction is secondary to endothelial damage associated with ischemia (i.e., CNI vasculopathy) or direct tubular toxicity, the findings of our study support the assertion that mitochondrial oxidative phosphorylation defects play a role as one of the central mechanisms of CNIT development. Given that the proximal tubular epithelium is mitochondria-rich and highly ATP-dependent, these cells may be particularly damaged as a result of CNIT [69,70,71]. However, oxidative stress has also been linked with endothelial dysfunction in CNIT, with a recent paper highlighting the role of TLR4 signaling and induction of vascular inflammation (TLR4 gene up-regulated in our data) [72].
Previous investigations of mitochondria and CNIT have utilized animal or in vitro models, making this among the first studies to demonstrate mitochondrial dysfunction in association with CNIT by evaluating gene pathways in human tissue. Our focus on canonical pathways permits the study of expression patterns of genes of interest within those pathways. Furthermore, our approach allows for the identification of important mRNA–miRNA interactions in association with these gene pathways and biological functions, incorporating yet another -omics layer to further delineate the molecular mechanisms of CNIT. Of the 13 miRNAs identified in the integrative network, several have been cited in the literature as being regulators of mitochondrial activity, including miR-16-5p, the miR-30 family, miR-26, and miR-24 [73,74,75,76].
The findings in this study are strengthened by the systems biology holistic approach. Systems biology is an essential schema for the evaluation of interrelated networks derived from multi-omic data representing multiple layers of genomic regulation [77,78]. Such a framework is necessary for understanding the KT model, given the complex and dynamic interactions between recipient host and donor allograft. This study investigated the interaction between two layers of –omic regulation (miRNAs and their mRNA targets in the same samples) to further understand the molecular underpinnings of the CNIT transcriptome. Our secondary objective was to identify a panel of candidate miRNA–mRNA targets among CNIT-specific genes to provide the basis to facilitate diagnostic biomarker discovery and prospective study. A biomarker panel comprised of miRNAs holds additional promise as a therapeutic target by utilizing antagomirs as drugs to suppress miRNA expression [79,80]. Therefore, identifying key miRNAs in CNIT may represent a path toward drug development in ameliorating its nephrotoxicity.
It is important to note that this pilot study does have the additional strength of inclusion of archival tissue samples. Formalin-fixed paraffin-embedded (FFPE) tissue samples are readily available, pose no additional risk to patients, and mitigate the challenge of obtaining adequate fresh tissue specimens for study. Additionally, archived FFPE tissue blocks represent a vast resource given its existence as part of highly annotated repositories/tissue banks. The amount of tissue used as RNA input was minimal (two to three 10-micron sections), preventing the risk of “exhausting” the sample and allowing use of small biopsy cores. Over the past several years, archived FFPE samples have emerged as a suitable tissue source for various -omics platforms as a result of advances in RNA and DNA isolation technology [81,82,83,84]. This study demonstrates the successful utilization of FFPE samples as an RNA source for both gene expression and miRNA microarrays.
The limitations of the study include its small sample size and cross-sectional study design. Nonetheless, our study is one of the first of its kind to be performed in an exclusively pediatric KT cohort, with the use of only well-characterized biopsy samples (i.e., only samples without mixed diagnoses were utilized). The reported findings are based on a cross-sectional examination of transcriptomic and epigenetic changes associated with pediatric CNIT, and they are intended to serve as a framework for future studies that will improve the knowledge of pertinent molecular pathways. An essential next step will be to validate the identified panel of markers in a large, prospective pediatric cohort with protocol biopsies and well-defined clinical and demographic characteristics.

4. Materials and Methods

4.1. Study Cohort

Thirty-eight FFPE kidney biopsies from pediatric KT recipients were available for study, 12 of which represented cases with CNIT. Normal allografts served as negative controls (n = 12), while 14 samples made up the positive control group, which functioned to determine marker/pathway specificity, including AR group (n = 7) and IFTA group (n = 7). The AR group included 5 patients with acute cellular rejection (ACR) and 2 patients with antibody-mediated rejection (AMR).
Of the 38 biopsy samples utilized for gene expression, 1 CNIT sample did not pass hybridization quality control with a %P of < 50%. The final dataset for gene expression was comprised of 37 total samples (11 CNIT; 12 Normal; 7 AR; 7 IFTA). All 15 miRNA microarrays were utilized for analysis.
No subjects were recruited as part of this study, and only de-identified archival FFPE tissue blocks were utilized. This study was verified by the Institutional Review Board at the University of Virginia (IRB-HSR # 18482) to meet the criteria for review exemption. Inclusion criteria for participation in this study include patients ≤18 years of age with biopsies taken over 6 months post-transplantation. Normal allograft protocol biopsies were taken 12 months post-transplantation.
CNIT was defined by the histological features of isometric vacuolization of the proximal convoluted tubules, nodular hyalinization of arterioles and small arteries, or striped interstitial fibrosis—in the absence of rejection, acute tubular necrosis, and/or IFTA [20,85]. ACR, AMR, and IFTA were defined as specified by the Banff criteria [85,86,87].

4.2. RNA Isolation and Microarray Hybridization

Total RNA was isolated (High Pure RNA Paraffin Kit, Roche, IN, USA) per manufacturer instructions. For each sample, three 10-micron sections (de-paraffinized in xylene) were later pooled in order to increase RNA yield. Assessment of RNA concentration and purity was performed using spectrophotometry (NanoDrop ® 1000). Only those samples demonstrating (1) enough quantity of RNA (>200 ng of RNA input) and (2) high-quality RNA (ratios of absorbance at 260 nm to 280 nm between 1.9 and 2.1) were considered suitable for downstream reactions. RNA was amplified and labeled (SensationPlus™ Kit; Affymetrix, Santa Clara, CA, USA), and consequently hybridized to miRNA and gene expression microarrays (Affymetrix™GeneChip HG-U133A 2.0 microarray, GeneChip miRNA v4.0 array, Santa Clara, CA, USA). Arrays were scanned with a GeneChip™ Scanner 3000 (GEO accession number GSE174020).

4.3. Gene Expression/miRNA Microarrays

All 38 biopsies were evaluated using gene expression microarrays. A subset of 15 biopsies was used for miRNA evaluation (n = 10 CNIT, n = 5 Normal). In order to further refine the molecular signature for CNIT, AR and IFTA samples were included for marker specificity as a “clean-up technique” to identify specific genes and pathways uniquely associated with CNIT and exclude common genes of allograft injury. The AR group was a composite of both ACR and AMR as representative of immune-mediated injury. IFTA was included as an archetype for chronic renal injury.

4.4. Gene Expression/miRNA Microarray Analysis and Quality Control

Raw intensities for each probeset were stored as electronic files (.CEL format). The robust multiarray average (RMA) method was used for background correction, which was followed by quantile normalization in the R environment [88,89,90,91]. A probeset level t-test was used to compare conditions. A p-value < 0.001 (<0.005 for miRNA microarrays) and q-value probeset-specific false discovery rate (FDR) ≤ 5% (Benjamini and Hochberg method [92,93]) were utilized to identify differentially expressed genes and miRNAs, respectively. Fold-change (genes ≥ 1.5 fold change and miRNAs ≥ 2) was used for relative quantification of gene and miRNA expression differences between groups. Only microarrays that passed quality control using intensity values of RNA spike-in controls were included.
Given the high likelihood of degraded RNA present in FFPE tissue, stringent hybridization quality control measures were applied [94]. The percentage of probesets declared present (%P) by the detection call algorithm: (cut-off > 50% ± 15%) and 3′:5′ ratios of the intensity values for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) housekeeping gene (cut-off: ≤3).

4.5. Interaction Networks, Functional/Pathway Analysis, Upstream Regulators

Ingenuity Pathway Analysis™ (IPA, http://ingenuity.com, accessed 10 May 2021) software was used to identify experimentally verified and high confidence target gene interactions to define miRNA:mRNA pairs for input into pathway enrichment analysis. MiRNA and mRNA datasets were also integrated and subjected to network analyses by IPA to identify plausible associations and potential regulatory networks relating to the progression of CNIT. Specifically, gene ontology, disease relevance, and functional and network analyses were assessed. Predicted upstream regulators with significant overlap p-value (defined as the p-value associated with the degree of overlap between observed genes present in the dataset and predicted genes of an upstream regulator) were also identified. Separate analyses were performed for CNIT vs. Normal, IFTA vs. Normal, and AR vs. Normal. Then, a comparison analysis of the three conditions was evaluated, where overlapping differentially expressed genes represented common genes, while non-overlapping genes were considered unique. Statistical significance associated with networks, pathways, and functions was calculated in IPA using a Fisher’s exact test, with p-values < 0.05 considered significant. Predictions regarding the activation state of a function, molecule, or pathway were inferred using the IPA regulation z-score algorithm (positive and negative z-scores predict activation and inhibition respectively), which is computed based upon the observed directionality gene expression changes present in the dataset.

4.6. Data Integration Analysis

A separate analysis was performed using unique genes for CNIT (those remaining following the elimination of common genes expressed by IFTA and AR vs. Normal samples). Then, an integration step was carried out between differentially expressed miRNA (CNIT vs. Normal) and this dataset of unique CNIT genes (IPA microRNA Target Filter tool, which is based on content from TarBase, TargetScan, miRecords, and published literature). This analysis generated a signature of differentially expressed miRNAs and their target mRNAs present in both datasets. The final miRNA:mRNA pairs represented a filtered list of only those experimentally observed relationships pairing in the proper directions (upregulated miRNAs and their downregulated gene targets and vice versa).

5. Conclusions

Our results have identified mitochondrial dysfunction and oxidative phosphorylation defects to be central gene pathways involved in the pathogenesis of CNIT, with renal cell death/necrosis and energy failure/ATP depletion as salient biological functions. Additionally, we discovered a panel of miRNAs and their gene targets representing a candidate molecular signature of CNIT in pediatric KT recipients. Delineation of the molecular pathways specific for CNIT offers mechanistic insights into a clinically significant problem and is a necessary first step for the development of minimally invasive biomarkers, which may allow individualized immunosuppressive therapy and improve graft survival and quality of life for children.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22115414/s1. Figure S1: Top 10 canonical pathways of CNIT. Figure S2: Expected vs. observed oxidative phosphorylation pathway. Table S1: Differentially expressed probesets between CNIT and Normal samples (provided separately as Excel Sheet). Table S2: Differentially expressed CNIT-specific genes related to oxidative phosphorylation. Table S3: Differentially expressed miRNAs between CNIT and Normal samples.

Author Contributions

Conceptualization, E.T.R. and V.M.; methodology, E.T.R., S.V.B. and V.M.; data analysis, S.V.B. and C.I.D.; histological evaluation, J.A.A., H.C., P.D.W. and D.M.; sample acquisition, H.C. and P.D.W.; data curation, E.T.R., E.B. and S.V.B.; writing—original draft preparation, E.T.R. and E.B.; writing—review and editing, V.M. and E.B.; supervision, D.M. and V.M.; project administration, D.M. and V.M.; funding acquisition, E.T.R., D.M. and V.M. All authors have read and agreed to the published version of the manuscript.

Funding

The research included in this manuscript was supported by the following funding sources: National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) training grant, T32DK072922, in addition to R21DK100678 (PI: Maluf), R01DK109581 and R01DK122682 (PI: Mas). E Rhone was also the recipient of a University of Virginia Children’s Hospital Grant-In-Aid Award.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of the University of Virginia (IRB-HSR# 18482).

Informed Consent Statement

Written informed consent has not been obtained from the patients because only de-identified tissue blocks were utilized.

Data Availability Statement

Arrays were scanned with a GeneChip™ Scanner 3000 (GSE174020).

Acknowledgments

The authors would like to thank Ellen Satter, Research Manager at Arkana Laboratories, for all her assistance in this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McDonald, S.P.; Craig, J.C. Australian and New Zealand Paediatric Nephrology Association. Long-term survival of children with end-stage renal disease. N. Engl. J. Med. 2004, 350, 2654–2662. [Google Scholar] [CrossRef]
  2. Shroff, R.; Ledermann, S. Long-term outcome of chronic dialysis in children. Pediatr. Nephrol. 2009, 24, 463–474. [Google Scholar] [CrossRef] [Green Version]
  3. Gillen, D.L.; Stehman-Breen, C.O.; Smith, J.M.; McDonald, R.A.; Warady, B.A.; Brandt, J.R.; Wong, C.S. Survival Advantage of Pediatric Recipients of a First Kidney Transplant Among Children Awaiting Kidney Transplantation. Arab. Archaeol. Epigr. 2008, 8, 2600–2606. [Google Scholar] [CrossRef]
  4. Tonelli, M.; Wiebe, N.; Knoll, G.; Bello, A.; Browne, S.; Jadhav, D.; Klarenbach, S.; Gill, J. Systematic Review: Kidney Transplantation Compared With Dialysis in Clinically Relevant Outcomes. Arab. Archaeol. Epigr. 2011, 11, 2093–2109. [Google Scholar] [CrossRef]
  5. Meier-Kriesche, H.-U.; Schold, J.D.; Srinivas, T.R.; Kaplan, B. Lack of Improvement in Renal Allograft Survival Despite a Marked Decrease in Acute Rejection Rates Over the Most Recent Era. Arab. Archaeol. Epigr. 2004, 4, 378–383. [Google Scholar] [CrossRef]
  6. Marcén, R.; Fernández-Rodriguez, A.; Rodríguez-Mendiola, N.; Ponte, B.; Galeano, C.; Villafruela, J.; Teruel, J.; Burgos, F.; Ortuño, J. Evolution of Rejection Rates and Kidney Graft Survival: A Historical Analysis. Transplant. Proc. 2009, 41, 2357–2359. [Google Scholar] [CrossRef] [Green Version]
  7. Van Arendonk, K.J.; Boyarsky, B.J.; Orandi, B.J.; James, N.T.; Smith, J.M.; Colombani, P.M.; Segev, D.L. National trends over 25 years in pediatric kidney transplant outcomes. Pediatrics 2014, 133, 594–601. [Google Scholar] [CrossRef] [Green Version]
  8. Hart, A.; Lentine, K.L.; Smith, J.M.; Miller, J.M.; Skeans, M.A.; Prentice, M.; Robinson, A.; Foutz, J.; Booker, S.E.; Israni, A.K.; et al. OPTN/SRTR 2019 Annual Data Report: Kidney. Am. J. Transplant 2021, 21 (Suppl. S2), 21–137. [Google Scholar] [CrossRef] [PubMed]
  9. Laskin, B.L.; Mitsnefes, M.M.; Dahhou, M.; Zhang, X.; Foster, B.J. The mortality risk with graft function has decreased among children receiving a first kidney transplant in the United States. Kidney Int. 2015, 87, 575–583. [Google Scholar] [CrossRef] [Green Version]
  10. Sayegh, M.H.; Carpenter, C.B. Transplantation 50 years later: Progress, challenges, and promises. N. Engl. J. Med. 2004, 351, 2761–2766. [Google Scholar] [CrossRef]
  11. Chapman, J.R.; O’Connell, P.J.; Nankivell, B.J. Chronic renal allograft dysfunction. J. Am. Soc. Nephrol. 2005, 16, 3015–3026. [Google Scholar] [CrossRef]
  12. Langewisch, E.; Mannon, R.B. Chronic Allograft Injury. Clin. J. Am. Soc. Nephrol. 2021, CJN.15590920. [Google Scholar] [CrossRef] [PubMed]
  13. Nankivell, B.J.; Kuypers, D.R. Diagnosis and prevention of chronic kidney allograft loss. Lancet 2011, 378, 1428–1437. [Google Scholar] [CrossRef]
  14. Farouk, S.S.; Rein, J.L. The Many Faces of Calcineurin Inhibitor Toxicity-What the FK? Adv. Chronic Kidney Dis. 2020, 27, 56–66. [Google Scholar] [CrossRef] [PubMed]
  15. Nankivell, B.J.; Borrows, R.J.; Fung, C.L.-S.; O’Connell, P.J.; Allen, R.D.; Chapman, J.R. The Natural History of Chronic Allograft Nephropathy. N. Engl. J. Med. 2003, 349, 2326–2333. [Google Scholar] [CrossRef] [Green Version]
  16. Myers, B.D.; Ross, J.; Newton, L.; Luetscher, J.; Perlroth, M. Cyclosporine-Associated Chronic Nephropathy. N. Engl. J. Med. 1984, 311, 699–705. [Google Scholar] [CrossRef]
  17. Calne, R.Y.; White, D.J.; Thiru, S.; Evans, D.B.; McMaster, P.; Dunn, D.C.; Pentlow, B.D.; Rolles, K. Cyclosporin A in patients receiving renal allograft from cadaveric donors. Lancet 1978, 2, 1323–1327. [Google Scholar] [CrossRef]
  18. Naesens, M.; Kuypers, D.R.J.; Sarwal, M. Calcineurin Inhibitor Nephrotoxicity. Clin. J. Am. Soc. Nephrol. 2009, 4, 481–508. [Google Scholar] [CrossRef] [Green Version]
  19. Kambham, N.; Nagarajan, S.; Shah, S.; Li, L.; Salvatierra, O.; Sarwal, M.M. A Novel, Semiquantitative, Clinically Correlated Calcineurin Inhibitor Toxicity Score for Renal Allograft Biopsies. Clin. J. Am. Soc. Nephrol. 2007, 2, 135–142. [Google Scholar] [CrossRef] [Green Version]
  20. Liptak, P.; Ivanyi, B. Primer: Histopathology of calcineurin-inhibitor toxicity in renal allografts. Nat. Clin. Pract. Nephrol. 2006, 2, 398–404. [Google Scholar] [CrossRef]
  21. Devarajan, P. Biomarkers for the early detection of acute kidney injury. Curr. Opin. Pediatr. 2011, 23, 194–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Snanoudj, R.; Royal, V.; Elie, C.; Rabant, M.; Girardin, C.; Morelon, E.; Kreis, H.; Fournet, J.-C.; Noël, L.-H.; Legendre, C. Specificity of histological markers of long-term CNI nephrotoxicity in kidney-transplant recipients under low-dose cyclosporine therapy. Am. J. Transpl. 2011, 11, 2635–2646. [Google Scholar] [CrossRef]
  23. Huang, Y.; Farkash, E. Protocol Biopsies: Utility and Limitations. Adv. Chronic Kidney Dis. 2016, 23, 326–331. [Google Scholar] [CrossRef] [PubMed]
  24. Furness, P.N.; Taub, N. International variation in the interpretation of renal transplant biopsies: Report of the CERTPAP Project. Kidney Int. 2001, 60, 1998–2012. [Google Scholar] [CrossRef] [Green Version]
  25. Serón, D.; Moreso, F.; Fulladosa, X.; Hueso, M.; Carrera, M.; Grinyó, J.M. Reliability of chronic allograft nephropathy diagnosis in sequential protocol biopsies. Kidney Int. 2002, 61, 727–733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Scian, M.J.; Maluf, D.G.; Mas, V.R. MiRNAs in kidney transplantation: Potential role as new biomarkers. Expert Rev. Mol. Diagn. 2013, 13, 93–104. [Google Scholar] [CrossRef]
  27. Fan, Y.; Chen, H.; Huang, Z.; Zheng, H.; Zhou, J. Emerging role of miRNAs in renal fibrosis. RNA Biol. 2020, 17, 1–12. [Google Scholar] [CrossRef]
  28. Maluf, D.G.; Dumur, C.I.; Suh, J.L.; Lee, J.K.; Cathro, H.P.; King, A.L.; Gallon, L.; Brayman, K.L.; Mas, V.R. Evaluation of Molecular Profiles in Calcineurin Inhibitor Toxicity Post-Kidney Transplant: Input to Chronic Allograft Dysfunction. Arab. Archaeol. Epigr. 2014, 14, 1152–1163. [Google Scholar] [CrossRef] [Green Version]
  29. Goldman, J.; Becker, M.L.; Jones, B.; Clements, M.; Leeder, J.S. Development of biomarkers to optimize pediatric patient management: What makes children different? Biomark. Med. 2011, 5, 781–794. [Google Scholar] [CrossRef] [Green Version]
  30. Burckart, G.J.; Green, D.J. The personalized medicine revolution: Getting it right for children. Pediatr. Transplant. 2012, 16, 530–532. [Google Scholar] [CrossRef]
  31. Savage, W.J.; Everett, A.D. Biomarkers in pediatrics: Children as biomarker orphans. Proteom. Clin. Appl. 2010, 4, 915–921. [Google Scholar] [CrossRef]
  32. Van Driest, S.L.; McGregor, T.L. Pharmacogenetics in clinical pediatrics: Challenges and strategies. PerMed 2013, 10, 661–671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dharnidharka, V.R.; Fiorina, P.; Harmon, W.E. Kidney transplantation in children. N. Engl. J. Med. 2014, 371, 549–558. [Google Scholar] [CrossRef]
  34. Leeder, J.S. Developmental and pediatric pharmacogenomics. Pharmacogenomics 2003, 4, 331–341. [Google Scholar] [CrossRef] [PubMed]
  35. Lipshultz, S.E.; Chandar, J.J.; Rusconi, P.G.; Fornoni, A.; Abitbol, C.L.; Burke, G.W.; Zilleruelo, G.E.; Pham, S.M.; Perez, E.E.; Karnik, R.; et al. Issues in solid-organ transplantation in children: Translational research from bench to bedside. Clinics 2014, 69, 55–72. [Google Scholar] [CrossRef]
  36. Ettenger, R.B. Age and the immune response in pediatric renal transplantation. Eur. J. Nucl. Med. Mol. Imaging 1992, 151, S7–S8. [Google Scholar] [CrossRef]
  37. Przepiorka, D.; Blamble, D.; Hilsenbeck, S.; Danielson, M.; Krance, R.; Chan, K.W. Tacrolimus clearance is age-dependent within the pediatric population. Bone Marrow Transplant. 2000, 26, 601–605. [Google Scholar] [CrossRef] [Green Version]
  38. Montini, G.; Ujka, F.; Varagnolo, C.; Ghio, L.; Ginevri, F.; Murer, L.; Thafam, B.S.; Carasi, C.; Zacchello, G.; Plebani, M. The pharmacokinetics and immunosuppressive response of tacrolimus in paediatric renal transplant recipients. Pediatr. Nephrol. 2006, 21, 719–724. [Google Scholar] [CrossRef]
  39. Naesens, M.; Salvatierra, O.; Li, L.; Kambham, N.; Concepcion, W.; Sarwal, M. Maturation of Dose-Corrected Tacrolimus Predose Trough Levels in Pediatric Kidney Allograft Recipients. Transplant 2008, 85, 1139–1145. [Google Scholar] [CrossRef]
  40. Marfo, K.; Altshuler, J.; Lu, A. Tacrolimus Pharmacokinetic and Pharmacogenomic Differences between Adults and Pediatric Solid Organ Transplant Recipients. Pharmaceutics 2010, 2, 291–299. [Google Scholar] [CrossRef] [Green Version]
  41. Birk, P.E. Surveillance biopsies in children post-kidney transplant. Pediatr. Nephrol. 2011, 27, 753–760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ettenger, R.B.; Blifeld, C.; Prince, H.; Gradus, D.B.-E.; Cho, S.; Sekiya, N.; Salusky, I.B.; Fine, R.N. The pediatric nephrologist’s dilemma: Growth after renal transplantation and its interaction with age as a possible immunologic variable. J. Pediatr. 1987, 111, 1022–1025. [Google Scholar] [CrossRef]
  43. Holcik, M.; Sonenberg, N. Translational control in stress and apoptosis. Nat. Rev. Mol. Cell Biol. 2005, 6, 318–327. [Google Scholar] [CrossRef] [PubMed]
  44. DuRose, J.B.; Scheuner, D.; Kaufman, R.J.; Rothblum, L.I.; Niwa, M. Phosphorylation of Eukaryotic Translation Initiation Factor 2α Coordinates rRNA Transcription and Translation Inhibition during Endoplasmic Reticulum Stress. Mol. Cell. Biol. 2009, 29, 4295–4307. [Google Scholar] [CrossRef] [Green Version]
  45. Sarbassov, D.D.; Ali, S.M.; Kim, D.H.; Guertin, D.A.; Latek, R.R.; Erdjument-Bromage, H.; Tempst, P.; Sabatini, D.M. Rictor, a novel binding partner of mTOR, defines a rapamycin-insensitive and raptor-independent pathway that regulates the cytoskeleton. Curr. Biol. 2004, 14, 1296–1302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Sato, T.; Ishii, J.; Ota, Y.; Sasaki, E.; Shibagaki, Y.; Hattori, S. Mammalian target of rapamycin (mTOR) complex 2 regulates filamin A-dependent focal adhesion dynamics and cell migration. Genes Cells 2016, 21, 579–593. [Google Scholar] [CrossRef] [Green Version]
  47. Yasuda, K.; Takahashi, M.; Mori, N. Mdm20 Modulates Actin Remodeling through the mTORC2 Pathway via Its Effect on Rictor Expression. PLoS ONE 2015, 10, e0142943. [Google Scholar] [CrossRef]
  48. Lamouille, S.; Connolly, E.; Smyth, J.W.; Akhurst, R.J.; Derynck, R. TGF-ß-induced activation of MTOR complex 2 drives epithelial-mesenchymal transition and cell invasion. J. Cell Sci. 2012, 125, 1259–1273. [Google Scholar] [CrossRef] [Green Version]
  49. Li, J.; Ren, J.; Liu, X.; Jiang, L.; He, W.; Yuan, W.; Yang, J.; Dai, C. RICTOR/mTORC2 signaling mediates TGFß1-induced fibroblast activation and kidney fibrosis. Kidney Int. 2015, 88, 515–527. [Google Scholar] [CrossRef] [Green Version]
  50. Betz, C.; Stracka, D.; Prescianotto-Baschong, C.; Frieden, M.; Demaurex, N.; Hall, M.N. mTOR complex 2-Akt signaling at mitochondria-associated endoplasmic reticulum membranes (MAM) regulates mitochondrial physiology. Proc. Natl. Acad. Sci. USA 2013, 110, 12526–12534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Appenzeller-Herzog, C.; Hall, M.N. Bidirectional crosstalk between endoplasmic reticulum stress and MTOR signaling. Trends Cell Biol. 2012, 22, 274–282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Arciuch, V.G.A.; Galli, S.; Franco, M.C.; Lam, P.Y.; Cadenas, E.; Carreras, M.C.; Poderoso, J.J. Akt1 intramitochondrial cycling is a crucial step in the redox modulation of cell cycle progression. PLoS ONE 2009, 4, e7523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Laplante, M.; Sabatini, D.M. mTOR signaling at a glance. J. Cell Sci. 2009, 122, 3589–3594. [Google Scholar] [CrossRef] [Green Version]
  54. Los, M.J.; Maddika, S.; Erb, B.; Schulze-Osthoff, K. Switching Akt: From survival signaling to deadly response. BioEssays 2009, 31, 492–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Nogueira, V.; Park, Y.; Chen, C.-C.; Xu, P.-Z.; Chen, M.-L.; Tonic, I.; Unterman, T.; Hay, N. Akt Determines Replicative Senescence and Oxidative or Oncogenic Premature Senescence and Sensitizes Cells to Oxidative Apoptosis. Cancer Cell 2008, 14, 458–470. [Google Scholar] [CrossRef] [Green Version]
  56. Venner, J.M.; Famulski, K.S.; Reeve, J.; Chang, J.; Halloran, P.F. Relationships among injury, fibrosis, and time in human kidney transplants. JCI Insight 2016, 1, e85323. [Google Scholar] [CrossRef] [Green Version]
  57. Kang, H.M.; Ahn, S.H.; Choi, P.; Ko, Y.-A.; Han, S.H.; Chinga, F.; Park, A.S.D.; Tao, J.; Sharma, K.; Pullman, J.; et al. Defective fatty acid oxidation in renal tubular epithelial cells has a key role in kidney fibrosis development. Nat. Med. 2015, 21, 37–46. [Google Scholar] [CrossRef]
  58. Zhao, L.; Dong, M.; Liao, S.; Du, Y.; Zhou, Q.; Zheng, H.; Chen, M.; Ji, J.; Gao, H. Identification of key metabolic changes in renal interstitial fibrosis rats using metabolomics and pharmacology. Sci. Rep. 2016, 6, 27194. [Google Scholar] [CrossRef] [Green Version]
  59. Gamboa, J.L.; Billings, F.T.; Bojanowski, M.T.; Gilliam, L.A.; Yu, C.; Roshanravan, B.; Roberts, L.J.; Himmelfarb, J.; Ikizler, T.A.; Brown, N.J. Mitochondrial dysfunction and oxidative stress in patients with chronic kidney disease. Physiol. Rep. 2016, 4, e12780. [Google Scholar] [CrossRef]
  60. Ishimoto, Y.; Inagi, R. Mitochondria: A therapeutic target in acute kidney injury. Nephrol. Dial. Transplant. 2016, 31, 1062–1069. [Google Scholar] [CrossRef]
  61. Zhan, M.; Brooks, C.; Liu, F.; Sun, L.; Dong, Z. Mitochondrial dynamics: Regulatory mechanisms and emerging role in renal pathophysiology. Kidney Int. 2013, 83, 568–581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Simon, N.; Morin, C.; Urien, S.; Tillement, J.P.; Bruguerolle, B. Tacrolimus and sirolimus decrease oxidative phosphorylation of isolated rat kidney mitochondria. Br. J. Pharmacol. 2003, 138, 369–376. [Google Scholar] [CrossRef] [Green Version]
  63. Servais, H.; Ortiz, A.; Devuyst, O.; Denamur, S.; Tulkens, P.M.; Mingeot-Leclercq, M.-P.; Arduan, A.O. Renal cell apoptosis induced by nephrotoxic drugs: Cellular and molecular mechanisms and potential approaches to modulation. Apoptosis 2007, 13, 11–32. [Google Scholar] [CrossRef] [PubMed]
  64. Palacín, M.; Coto, E.; Llobet, L.; Pacheu-Grau, D.; Montoya, J.; Ruiz-Pesini, E. FK506 affects mitochondrial protein synthesis and oxygen consumption in human cells. Cell Biol. Toxicol. 2013, 29, 407–414. [Google Scholar] [CrossRef] [PubMed]
  65. Hortelano, S.; Castilla, M.; Torres, A.M.; Tejedor, A.; Boscá, L. Potentiation by nitric oxide of cyclosporin A and FK506-induced apoptosis in renal proximal tubule cells. J. Am. Soc. Nephrol. 2000, 11, 2315–2323. [Google Scholar] [CrossRef]
  66. Henke, W.; Jung, K. Comparison of the effects of the immunosuppressive agents FK 506 and cyclosporin a on rat kidney mitochondria. Biochem. Pharmacol. 1993, 46, 829–832. [Google Scholar] [CrossRef]
  67. Ruiz-Cabello, J.; Buss, W.C.; Collier, S.W.; I Glazer, R.; Cohen, J.S. Changes in ATP after cyclosporin A treatment in a renal epithelial cell line in the rat studied by 31P-NMR spectroscopy. Res. Commun. Mol. Pathol. Pharmacol. 1994, 86, 3–13. [Google Scholar]
  68. Ortiz, A.; Tejedor, A.; Caramelo, C. Nephrotoxicity. In Drug-Induced Mitochondrial Dysfunction; Dykens, J.A., Will, Y., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2008; pp. 291–310. [Google Scholar]
  69. Forbes, J.M. Mitochondria-Power players in kidney function? Trends Endocrinol. Metab. 2016, 27, 441–442. [Google Scholar] [CrossRef]
  70. Higgins, G.C.; Coughlan, M.T. Mitochondrial dysfunction and mitophagy: The beginning and end to diabetic nephropathy? Br. J. Pharmacol. 2014, 171, 1917–1942. [Google Scholar] [CrossRef] [Green Version]
  71. Balaban, R.S.; Mandel, L.J. Metabolic substrate utilization by rabbit proximal tubule. An NADH fluorescence study. Am. J. Physiol. 1988, 254, F407–F416. [Google Scholar] [CrossRef]
  72. Rodrigues-Diez, R.; González-Guerrero, C.; Ocaña-Salceda, C.; Rodrigues-Diez, R.R.; Egido, J.; Ortiz, A.; Ruiz-Ortega, M.; Ramos, A.M. Calcineurin inhibitors cyclosporine A and tacrolimus induce vascular inflammation and endothelial activation through TLR4 signaling. Sci. Rep. 2016, 6, 27915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Li, P.; Jiao, J.; Gao, G.; Prabhakar, B.S. Control of mitochondrial activity by miRNAs. J. Cell. Biochem. 2011, 113, 1104–1110. [Google Scholar] [CrossRef] [Green Version]
  74. Sripada, L.; Tomar, D.; Singh, R. Mitochondria: One of the destinations of miRNAs. Mitochondrion 2012, 12, 593–599. [Google Scholar] [CrossRef] [PubMed]
  75. Tomasetti, M.; Neuzil, J.; Dong, L. MicroRNAs as regulators of mitochondrial function: Role in cancer suppression. Biochim. Biophys. Acta Gen. Subj. 2014, 1840, 1441–1453. [Google Scholar] [CrossRef] [PubMed]
  76. Duarte, F.V.; Palmeira, C.M.; Rolo, A.P. The Role of microRNAs in Mitochondria: Small Players Acting Wide. Genes 2014, 5, 865–886. [Google Scholar] [CrossRef] [Green Version]
  77. Bontha, S.V.; Maluf, D.G.; Mueller, T.F.; Mas, V.R. Systems Biology in Kidney Transplantation: The Application of Multi-Omics to a Complex Model. Arab. Archaeol. Epigr. 2016, 17, 11–21. [Google Scholar] [CrossRef] [Green Version]
  78. Liang, M.; Cowley, A.W., Jr.; Hessner, M.J.; Lazar, J.; Basile, D.P.; Pietrusz, J.L. Transcriptome analysis and kidney research: Toward systems biology. Kidney Int. 2005, 67, 2114–2122. [Google Scholar] [CrossRef] [Green Version]
  79. Li, Z.; Rana, T.M. Therapeutic targeting of microRNAs: Current status and future challenges. Nat. Rev. Drug Discov. 2014, 13, 622–638. [Google Scholar] [CrossRef]
  80. Mattes, J.; Yang, M.; Foster, P.S. Regulation of microRNA by antagomirs: A new class of pharmacological antagonists for the specific regulation of gene function? Am. J. Respir. Cell Mol. Biol. 2007, 36, 8–12. [Google Scholar] [CrossRef]
  81. Belder, N.; Coşkun, Ö.; Erdoğan, B.D.; Savaş, B.; Ensari, A.; Özdağ, H. Optimization of gene expression microarray protocol for formalin-fixed paraffin-embedded tissues. Genom. Data 2016, 7, 303–306. [Google Scholar] [CrossRef] [Green Version]
  82. Gehrau, R.C.; Mas, V.R.; Villamil, F.G.; Dumur, C.I.; Suh, J.L.; Maluf, D.G.; Mehta, N.K. MicroRNA Signature at the Time of Clinical HCV Recurrence Associates With Aggressive Fibrosis Progression Post-Liver Transplantation. Am. J. Transplant. 2013, 13, 729–737. [Google Scholar] [CrossRef]
  83. April, C.; Klotzle, B.; Royce, T.; Wickham-Garcia, E.; Boyaniwsky, T.; Izzo, J.; Cox, D.; Jones, W.; Rubio, R.; Holton, K.; et al. Whole-Genome Gene Expression Profiling of Formalin-Fixed, Paraffin-Embedded Tissue Samples. PLoS ONE 2009, 4, e8162. [Google Scholar] [CrossRef] [Green Version]
  84. Klopfleisch, R.; Weiss, A.T.; Gruber, A.D. Excavation of a buried treasure—DNA, mRNA, miRNA and protein analysis in formalin-fixed, paraffin-embedded tissues. Histol. Histolpathol. 2011, 26, 797–810. [Google Scholar]
  85. Solez, K.; Colvin, R.B.; Racusen, L.C.; Haas, M.; Sis, B.; Mengel, M.; Halloran, P.F.; Baldwin, W.; Banfi, G.; Collins, A.B.; et al. Banff 07 Classification of Renal Allograft Pathology: Updates and Future Directions. Arab. Archaeol. Epigr. 2008, 8, 753–760. [Google Scholar] [CrossRef]
  86. Solez, K.; Racusen, L.C. The Banff classification revisited. Kidney Int. 2013, 83, 201–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Haas, M.; Sis, B.; Racusen, L.C.; Solez, K.; Glotz, D.; Colvin, R.B.; Castro, M.C.R.; David, D.S.R.; David-Neto, E.; Bagnasco, S.M.; et al. Banff 2013 meeting report: Inclusion of c4d-negative antibody-mediated rejection and antibody-associated arterial lesions. Am. J. Transplant. 2014, 14, 272–283. [Google Scholar] [CrossRef]
  88. Barash, Y.; Dehan, E.; Krupsky, M.; Franklin, W.; Geraci, M.; Friedman, N.; Kaminski, N. Comparative analysis of algorithms for signal quantitation from oligonucleotide microarrays. Bioinformatics 2004, 20, 839–846. [Google Scholar] [CrossRef] [PubMed]
  89. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  90. Irizarry, R.A.; Hobbs, B.; Collin, F.; Beazer-Barclay, Y.D.; Antonellis, K.J.; Scherf, U.; Speed, T.P. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4, 249–264. [Google Scholar] [CrossRef] [Green Version]
  91. Carvalho, B.S.; Irizarry, R.A. A framework for oligonucleotide microarray preprocessing. Bioinformatics 2010, 26, 2363–2367. [Google Scholar] [CrossRef] [PubMed]
  92. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  93. Benjamini, Y.; Yekutieli, D. The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 2001, 29, 1165–1188. [Google Scholar] [CrossRef]
  94. Dumur, C.I.; Nasim, S.; Best, A.M.; Archer, K.J.; Ladd, A.C.; Mas, V.R.; Wilkinson, D.S.; Garrett, C.T.; Ferreira-Gonzalez, A. Evaluation of Quality-Control Criteria for Microarray Gene Expression Analysis. Clin. Chem. 2004, 50, 1994–2002. [Google Scholar] [CrossRef]
Ijms 22 05414 g001 550
Figure 1. Gene expression microarrays of 37 formalin-fixed paraffin-embedded (FFPE) kidney biopsies, including Normal, calcineurin inhibitor nephrotoxicity (CNIT), interstitial fibrosis and tubular atrophy (IFTA), and acute rejection (AR), were analyzed. Then, CNIT-specific differentially expressed genes (DEGs) were integrated with 15 miRNA microarray profiles to identify miRNA–mRNA interactions.
Figure 1. Gene expression microarrays of 37 formalin-fixed paraffin-embedded (FFPE) kidney biopsies, including Normal, calcineurin inhibitor nephrotoxicity (CNIT), interstitial fibrosis and tubular atrophy (IFTA), and acute rejection (AR), were analyzed. Then, CNIT-specific differentially expressed genes (DEGs) were integrated with 15 miRNA microarray profiles to identify miRNA–mRNA interactions.
Ijms 22 05414 g001
Ijms 22 05414 g002 550
Figure 2. Hierarchical clustering analysis showing differentially expressed microRNAs between Normal (WNL) and CNIT allografts. MiRNA signatures in tissue samples differentiate between histological conditions.
Figure 2. Hierarchical clustering analysis showing differentially expressed microRNAs between Normal (WNL) and CNIT allografts. MiRNA signatures in tissue samples differentiate between histological conditions.
Ijms 22 05414 g002
Ijms 22 05414 g003 550
Figure 3. Molecule activity predictor (MAP) pathway analysis demonstrating mitochondrial dysfunction (p = 2.90 × 10−32) as one of the top pathways affected in CNIT samples. Genes associated with complexes I-V of the electron transport chain (green) were significantly downregulated. Caspases 3, 8, 9, and Cytochrome C (yellow box) are predicted to be activated, demonstrating the in-silico activation of apoptosis.
Figure 3. Molecule activity predictor (MAP) pathway analysis demonstrating mitochondrial dysfunction (p = 2.90 × 10−32) as one of the top pathways affected in CNIT samples. Genes associated with complexes I-V of the electron transport chain (green) were significantly downregulated. Caspases 3, 8, 9, and Cytochrome C (yellow box) are predicted to be activated, demonstrating the in-silico activation of apoptosis.
Ijms 22 05414 g003
Ijms 22 05414 g004 550
Figure 4. (A) List of 13 differentially expressed miRNAs in CNIT samples. (B) Integrated network of 13 miRNAs and their downstream up- (red) and down- (green) regulated genes.
Figure 4. (A) List of 13 differentially expressed miRNAs in CNIT samples. (B) Integrated network of 13 miRNAs and their downstream up- (red) and down- (green) regulated genes.
Ijms 22 05414 g004
Ijms 22 05414 g005 550
Figure 5. QPCR validation of four candidate genes using CNIT and Normal samples. Plotted with mean values and 95% confidence intervals. Gene expression normalized using the GAPDH housekeeping gene. A decrease in ΔCt represents an increase in gene expression. SOD2 and CCND1 were significantly upregulated, while GRB10 and TLR4 showed a trend in significance.
Figure 5. QPCR validation of four candidate genes using CNIT and Normal samples. Plotted with mean values and 95% confidence intervals. Gene expression normalized using the GAPDH housekeeping gene. A decrease in ΔCt represents an increase in gene expression. SOD2 and CCND1 were significantly upregulated, while GRB10 and TLR4 showed a trend in significance.
Ijms 22 05414 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rhone, E.T.; Bardhi, E.; Bontha, S.V.; Walker, P.D.; Almenara, J.A.; Dumur, C.I.; Cathro, H.; Maluf, D.; Mas, V. An Integrated Transcriptomic Approach to Identify Molecular Markers of Calcineurin Inhibitor Nephrotoxicity in Pediatric Kidney Transplant Recipients. Int. J. Mol. Sci. 2021, 22, 5414. https://doi.org/10.3390/ijms22115414

AMA Style

Rhone ET, Bardhi E, Bontha SV, Walker PD, Almenara JA, Dumur CI, Cathro H, Maluf D, Mas V. An Integrated Transcriptomic Approach to Identify Molecular Markers of Calcineurin Inhibitor Nephrotoxicity in Pediatric Kidney Transplant Recipients. International Journal of Molecular Sciences. 2021; 22(11):5414. https://doi.org/10.3390/ijms22115414

Chicago/Turabian Style

Rhone, Erika T., Elissa Bardhi, Sai Vineela Bontha, Patrick D. Walker, Jorge A. Almenara, Catherine I. Dumur, Helen Cathro, Daniel Maluf, and Valeria Mas. 2021. "An Integrated Transcriptomic Approach to Identify Molecular Markers of Calcineurin Inhibitor Nephrotoxicity in Pediatric Kidney Transplant Recipients" International Journal of Molecular Sciences 22, no. 11: 5414. https://doi.org/10.3390/ijms22115414

APA Style

Rhone, E. T., Bardhi, E., Bontha, S. V., Walker, P. D., Almenara, J. A., Dumur, C. I., Cathro, H., Maluf, D., & Mas, V. (2021). An Integrated Transcriptomic Approach to Identify Molecular Markers of Calcineurin Inhibitor Nephrotoxicity in Pediatric Kidney Transplant Recipients. International Journal of Molecular Sciences, 22(11), 5414. https://doi.org/10.3390/ijms22115414

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