Go to The Journal of Clinical Investigation
  • About
  • Editors
  • Consulting Editors
  • For authors
  • Publication ethics
  • Publication alerts by email
  • Transfers
  • Advertising
  • Job board
  • Contact
  • Physician-Scientist Development
  • Current issue
  • Past issues
  • By specialty
    • COVID-19
    • Cardiology
    • Immunology
    • Metabolism
    • Nephrology
    • Oncology
    • Pulmonology
    • All ...
  • Videos
  • Collections
    • In-Press Preview
    • Resource and Technical Advances
    • Clinical Research and Public Health
    • Research Letters
    • Editorials
    • Perspectives
    • Physician-Scientist Development
    • Reviews
    • Top read articles

  • Current issue
  • Past issues
  • Specialties
  • In-Press Preview
  • Resource and Technical Advances
  • Clinical Research and Public Health
  • Research Letters
  • Editorials
  • Perspectives
  • Physician-Scientist Development
  • Reviews
  • Top read articles
  • About
  • Editors
  • Consulting Editors
  • For authors
  • Publication ethics
  • Publication alerts by email
  • Transfers
  • Advertising
  • Job board
  • Contact
Benchmarking urinary cell transcriptomes for noninvasive differentiation of BK polyomavirus–associated nephropathy from T cell–mediated rejection
Franco B. Mueller, Carol Li, Darshana M. Dadhania, Surya V. Seshan, Thalia Salinas, Vijay K. Sharma, Jenny Z. Xiang, Hans H. Hirsch, Thangamani Muthukumar, Manikkam Suthanthiran
Franco B. Mueller, Carol Li, Darshana M. Dadhania, Surya V. Seshan, Thalia Salinas, Vijay K. Sharma, Jenny Z. Xiang, Hans H. Hirsch, Thangamani Muthukumar, Manikkam Suthanthiran
View: Text | PDF
Research Article Immunology Nephrology

Benchmarking urinary cell transcriptomes for noninvasive differentiation of BK polyomavirus–associated nephropathy from T cell–mediated rejection

  • Text
  • PDF
Abstract

BK polyomavirus–associated nephropathy (BKVN) adversely impacts kidney allograft survival and often mimics acute T cell–mediated rejection (TCMR), confounding diagnosis and management. To address this conundrum, we performed unbiased RNA sequencing of urinary cells matched to biopsies classified as BKVN with intragraft inflammation (BKVN-P), BKVN without inflammation (BKVN-N), TCMR, or no rejection (NR). BKVN-N displayed dominant host DNA replication, cell cycle, and repair programs, while BKVN-P samples exhibited expansive innate immune activation, antigen presentation, chemokine upregulation, and epithelial injury. Both BKVN subtypes shared signatures of T cell exhaustion and mature and tolerogenic dendritic cell activation but differed in immune orientation — Th1 predominance in BKVN-N versus Treg and CD8 enrichment in BKVN-P. Compared with TCMR samples, BKVN-P lacked robust TCR/CD28 signaling and was enriched for viral and innate modules; BKVN-N lacked alloimmune activation. B cell exhaustion characterized BKVN-N, while BKVN-P displayed robust B cell activation with metabolic downregulation. A ratiometric urinary cell biomarker, CXCL10 mRNA/CD3E mRNA, distinguished both BKVN subtypes from TCMR with diagnostic accuracy, replicated by quantitative reverse transcription PCR for clinical translation, and confirmed in an independent cohort. These findings demonstrate the utility of urinary cell transcriptomics for resolving viral injury from alloimmunity, enabling precision diagnostics and targeted immunomodulation in kidney transplantation.

Authors

Franco B. Mueller, Carol Li, Darshana M. Dadhania, Surya V. Seshan, Thalia Salinas, Vijay K. Sharma, Jenny Z. Xiang, Hans H. Hirsch, Thangamani Muthukumar, Manikkam Suthanthiran

×

Figure 3

Differential gene expression analysis of urinary cell transcriptomes.

Options: View larger image (or click on image) Download as PowerPoint
Differential gene expression analysis of urinary cell transcriptomes.
(A...
(A) Volcano plot of urinary cell gene expression in 7 urines matched to 6 BKVN-P biopsies versus 26 urines matched to 26 NR biopsies. (B) CDF plot of gene expression ratios in 7 urines matched to 6 BKVN-P biopsies versus 26 urines matched to 26 NR biopsies. (C) Volcano plot of urinary cell gene expression in 4 urines matched to 4 BKVN-N biopsies versus 26 urines matched to 26 NR biopsies. (D) CDF plot of gene expression ratios in 4 urines matched to 4 BKVN biopsies versus 26 urines matched to 26 NR biopsies. A total of 17,756 genes were analyzed. In volcano plots, each red dot represents an overexpressed gene, and each blue dot represents an underexpressed gene, based on FDR-adjusted P value < 0.05. Gray/black dots indicate genes that are not differentially expressed. The x axis shows the log2 fold change (FC), calculated as the ratio of mRNA counts between the 2 groups being compared. The y axis represents the –log10 of the P values. In the CDF plots, the blue lines represent all 17,756 genes, the red lines represent DEGs, and the green lines represent 3,659 ubiquitously expressed housekeeping genes (24). The shift of the red curve relative to the respective green curve is statistically significant based on the Kolmogorov-Smirnov test. Supplemental Tables 2 and 3 list the BKVN-P versus NR DEGs shown in A and B. Supplemental Tables 4 and 5 list the BKVN-N versus NR DEGs shown in C and D.

Copyright © 2026 American Society for Clinical Investigation
ISSN 2379-3708

Sign up for email alerts