Fig. 2. Transcript-eGFR association and EGF expression in the kidney in patients with CKD. (A) Correlation between observed and predicted eGFR (continuous line) on the basis of intrarenal transcripts integrated by ridge regression analysis (top 1 to 30 transcripts shown), providing eGFR prediction in the validation cohort (n = 55). Gray area represents 95% confidence interval (CI). Black arrow indicates the highest correlation provided by a six-marker set, which is depicted in (B). (B) Maximal correlation (r = 0.77, P < 0.0001) was demonstrated between observed eGFR and that predicted by gene expression values of six transcripts (n = 55). Dotted line represents 95% CI. (C) Intrarenal EGF mRNA showed a significant correlation with patients’ eGFR in discovery (n = 164) and validation cohorts (n = 55 and 42). TLDA, TaqMan Low Density Array. (D) EGF mRNA expression in major human organs/tissues (a selection out of a panel of 84 human organs/tissues and cell lines), derived from BioGPS (http://biogps.org), indicating a highly kidney-specific expression pattern. For full data set, see fig. S2. (E) EGF mRNA expression pattern in adult kidney (glomeruli, inner and outer cortex, inner and outer medulla, papillary tips, and renal pelvis); data were extracted from Higgins et al. (10). (F) In situ hybridization demonstrates tubule-specific EGF mRNA in both cortex (I and III) and medulla (II and IV), with reduced expression in CKD (III and IV) compared to healthy controls (I and II). Black arrows indicate positive staining in pink. Scale bars, 50 μM.
Fig. 3. Correlation of uEGF/Cr with intrarenal EGF mRNA and eGFR. (A and B) uEGF/Cr is correlated with intrarenal EGF mRNA in patients with matching urine enzyme-linked immunosorbent assay (ELISA) data and tissue mRNA expression data in C-PROBE (A) (n = 34) and NEPTUNE (B) (n = 85) cohorts. uEGF/Cr is correlated significantly (P < 0.0001) with eGFR at the time of biopsy in patients from C-PROBE [(C) n = 349], NEPTUNE [(D) n = 141], and PKU-IgAN [(E) n = 452].
Fig. 4. Association of uEGF/Cr with tubulointerstitial damage. (A and B) Tubulointerstitial damage is reflected by the percentage of cortex affected by IF/TA. Scoring of IF/TA was based on evaluation of silver, periodic acid–Schiff (PAS), and trichrome-stained kidney sections (n = 102) by six readers who were blinded to the uEGF/Cr value. For each section, the average scores of the six readers for IF and TA were calculated as indicators of tubulointerstitial damage. uEGF/Cr is significantly correlated with IF (A) and TA (B) scores (Spearman correlation, P < 0.001).
Fig. 5. Association of EGF with eGFR slope. (A) Intrarenal EGF RNA expression correlated significantly (P < 0.001) with eGFR slope in C-PROBE patients (n = 29). (B and C) Correlation of the observed eGFR slope of CKD patients in C-PROBE (n = 344) with slope predicted by uEGF/Cr (B) or ACR (C) using a regression model (adjusted for age and gender). eGFR slope predicted by uEGF/Cr (B) exhibited a higher correlation with the observed value than slope predicted by ACR (C).
Fig. 6. Multivariable-adjusted HRs for predicting the composite end point on the basis of uEGF/Cr. HRs were adjusted by age, gender, eGFR, and ACR. Adjusted HRs and 95% CIs were obtained by separate Cox regression models in each study cohort. A one-unit decrease of uEGF/Cr (in log scale) was associated with an increased risk of CKD progression of 3.73 (1.85 to 7.69)–fold, 3.43 (1.72 to 6.67)–fold, and 1.96 (1.45 to 2.70)–fold in these three cohorts. The unadjusted HRs for EGF were 0.33 (0.21 to 0.51) (C-PROBE), 0.33 (0.21 to 0.52) (NEPTUNE), and 0.57 (0.46 to 0.70) (PKU-IgAN).
- Table 1 Demographic characteristics of CKD patients in the gene expression discovery and validation cohorts.
Age and eGFR are presented as means ± SD. SLE, lupus nephropathy; RPGN, rapidly progressive glomerulonephritis; HTN, hypertensive nephropathy; MGN, membranous glomerulonephritis; DN, diabetic nephropathy; FSGS, focal segmental glomerulosclerosis; MCD, minimal change disease; TMD, thin basement membrane disease.
Disease type CKD patients eGFR (ml/min per 1.73 m2) Age (years) Gender (% female) Discovery cohort: ERCB SLE 30 63.7 ± 29.4 34.7 ± 13.3 76.7 IgAN 24 75.9 ± 37.9 36.4 ± 14.6 25 RPGN 21 46.6 ± 31.5 58.5 ± 14.1 42.9 HTN 20 43.9 ± 25.1 57.2 ± 12.1 25 MGN 18 88.9 ± 41.4 53.4 ± 19.3 44.4 DN 17 44.3 ± 24.9 58.3 ± 10.7 29.4 FSGS 16 73.4 ± 38.4 46.2 ± 17.6 56.2 MCD 12 100.7 ± 33.9 35.8 ± 16.8 33.3 TMD 6 93.4 ± 29.4 46.0 ± 14.5 33.3 Total 164 66.4 ± 37.2 46.9 ± 17.6 43.3 Validation cohort 1: ERCB SLE 10 60.1 ± 31.5 37.1 ± 13.8 60 IgAN 17 50.8 ± 34.1 46.0 ± 17.5 29.4 RPGN 5 20.9 ± 15.0 50.4 ± 21.7 20 MGN 4 54.4 ± 32.4 58.7 ± 23.6 50 DN 1 97.3 44.8 100 FSGS 8 55.9 ± 36.0 43.6 ± 18.6 25 MCD 4 86.0 ± 44.5 44.4 ± 19.3 0 Other 6 68.4 ± 60.7 45.2 ± 16.7 33.3 Total 55 55.7 ± 38.3 45.1 ± 17.6 32.7 Validation cohort 2: C-PROBE SLE 16 106.4 ± 60.1 31.3 ± 9.9 87.5 IgAN 1 110.4 23 100 HTN 1 146.5 29 0 MGN 2 51.9 ± 21.9 53.0 ± 5.7 50 DN 1 50.0 80.8 100 FSGS 11 51.2 ± 21.0 41.9 ± 14.0 27.2 MCD 1 124.7 22 100 Other 9 41.3 ± 23.1 52.4 ± 18.9 55.6 Total 42 76.3 ± 50.8 39.6 ± 15.7 61.9 - Table 2 Baseline characteristics of patients whose urine samples were analyzed in urine studies.
Baseline characteristics C-PROBE (n = 349) NEPTUNE (n = 141) PKU-IgAN (n = 452) Age 54.8 ±16.1 45.7 ± 17.3 36.0 ± 12.0 Gender (% female) 59.3 35.5 51.3 Race (% black) 47.9 23.4 0 eGFR (ml/min per 1.73 m2) 55.4 ± 32.2 67.7 ± 37.1 84.8 ± 34.1 ACR (mg/mg) 0.8 ± 1.2 2.3 ± 2.1 0.6 ± 0.6 Log2 uEGF/Cr (ng/mg) 2.5 ± 1.1 3.0 ± 1.3 3.5 ± 1.0 - Table 3 The association of uEGF with time to composite end point.
The association was evaluated by a Cox proportional hazards regression analysis in C-PROBE, NEPTUNE, and PKU-IgAN cohorts. Renal event is defined as the presence of composite end point of ESKD or 40% reduction in baseline eGFR. The follow-up lengths (in years) for this analysis for the three cohorts were 1.8 ± 0.8, 2.0 ± 0.8, and 3.6 ± 2.2, respectively. AIC, Akaike information criterion; LR, likelihood ratio test.
Cohort Model AIC C-statistics P value (LR test) C-PROBE (25 events/189 patients) M1 eGFR + ACR* 206.3 0.75 (0.58–0.91) M2 eGFR + ACR* + EGF 192.1 0.87 (0.77–0.97) <0.0001 NEPTUNE (26 events/110 patients) M1 eGFR + ACR* 213.98 0.74 (0.61–0.86) M2 eGFR + ACR* + EGF 204.14 0.80 (0.68–0.92) 0.0006 PKU-IgAN (68 events/428 patients) M1 eGFR + ACR* 660.18 0.71 (0.59–0.83) M2 eGFR + ACR* + EGF 647.94 0.75 (0.64–0.86) 0.00016 *Adjusted for age and gender.
Supplementary Materials
www.sciencetranslationalmedicine.org/cgi/content/full/7/316/316ra193/DC1
Materials and Methods
Fig. S1. Baseline eGFR prediction by a three-marker panel in C-PROBE.
Fig. S2. EGF and two other candidate mRNAs’ expression patterns in a panel of 84 human organs, tissues, and cell lines.
Fig. S3. Correlation of EGF mRNA and uEGF/Cr with eGFR and eGFR slope in DN patients and CKD patients with diabetes.
Fig. S4. Correlation of uEGF/Cr with ACR.
Fig. S5. EGF as the top upstream regulator of genes correlated with eGFR slope.
Fig. S6. Correlation of uEGF/Cr with eGFR slope.
Fig. S7. ROC curve and corresponding AUC statistics for models with and without uEGF.
Fig. S8. mRNA localization by in situ hybridization: Negative and positive control images.
Table S1. Significantly enriched canonical pathways in intrarenal marker set.
Table S2. qRT-PCR assays used to validate expression of the intrarenal transcripts.
Table S3. Correlations of identified intrarenal transcripts with log2 eGFR of patients from the discovery and two validation cohorts.
Table S4. Demographic characteristics of NEPTUNE patients with intrarenal EGF expression data available.
Table S5. Top 10 upstream regulators of transcripts correlated with eGFR change over time (eGFR slope).
References (45–53)
Additional Files
- Supplementary Material for:
Tissue transcriptome-driven identification of epidermal growth factor as a chronic kidney disease biomarker
Wenjun Ju,* Viji Nair, Shahaan Smith, Li Zhu, Kerby Shedden, Peter X. K. Song, Laura H. Mariani, Felix H. Eichinger, Celine C. Berthier, Ann Randolph, Jennifer Yi-Chun Lai, Yan Zhou, Jennifer J. Hawkins, Markus Bitzer, Matthew G. Sampson, Martina Thier, Corinne Solier, Gonzalo C. Duran-Pacheco, Guillemette Duchateau-Nguyen, Laurent Essioux, Brigitte Schott, Ivan Formentini, Maria C. Magnone, Maria Bobadilla, Clemens D. Cohen, Serena M. Bagnasco, Laura Barisoni, Jicheng Lv, Hong Zhang, Hai-Yan Wang, Frank C. Brosius, Crystal A. Gadegbeku, Matthias Kretzler,* for the ERCB, C-PROBE, NEPTUNE, and PKU-IgAN Consortium
*Corresponding author. E-mail: wenjunj{at}med.umich.edu (W.J.); kretzler{at}umich.edu (M.K.)
Published 2 December 2015, Sci. Transl. Med. 7, 316ra193 (2015)
DOI: 10.1126/scitranslmed.aac7071This PDF file includes:
- Materials and Methods
- Fig. S1. Baseline eGFR prediction by a three-marker panel in C-PROBE.
- Fig. S2. EGF and two other candidate mRNAs’ expression patterns in a panel of 84 human organs, tissues, and cell lines.
- Fig. S3. Correlation of EGF mRNA and uEGF/Cr with eGFR and eGFR slope in DN patients and CKD patients with diabetes.
- Fig. S4. Correlation of uEGF/Cr with ACR.
- Fig. S5. EGF as the top upstream regulator of genes correlated with eGFR slope.
- Fig. S6. Correlation of uEGF/Cr with eGFR slope.
- Fig. S7. ROC curve and corresponding AUC statistics for models with and without uEGF.
- Fig. S8. mRNA localization by in situ hybridization: Negative and positive control images.
- Table S1. Significantly enriched canonical pathways in intrarenal marker set.
- Table S2. qRT-PCR assays used to validate expression of the intrarenal transcripts.
- Table S3. Correlations of identified intrarenal transcripts with log2 eGFR of patients from the discovery and two validation cohorts.
- Table S4. Demographic characteristics of NEPTUNE patients with intrarenal EGF expression data available.
- Table S5. Top 10 upstream regulators of transcripts correlated with eGFR change over time (eGFR slope).
- References (45–53)