Cytomegalovirus (CMV) viremia continues to cause significant morbidity and mortality in kidney transplant patients with clinical complications including organ rejection and death. Whole blood gene expression dynamics in CMV viremic patients from onset of DNAemia through convalescence has not been well studied to date in humans. To evaluate how CMV infection impacts whole blood leukocyte gene expression over time, we evaluated a matched cohort of 62 kidney transplant recipients with and without CMV DNAemia using blood samples collected at multiple time points during the 12-month period after transplant. While transcriptomic differences were minimal at baseline between DNAemic and non-DNAemic patients, hundreds of genes were differentially expressed at the long-term timepoint, including genes enriching for pathways important for macrophages, interferon, and IL-8 signaling. Amongst patients with CMV DNAemia, the greatest amount of transcriptomic change occurred between baseline and 1-week post-DNAemia, with increase in pathways for interferon signaling and cytotoxic T cell function. Time-course gene set analysis of these differentially expressed genes revealed that most of the enriched pathways had a significant time-trend. While many pathways that were significantly down- or upregulated at 1 week returned to baseline-like levels, we noted that several pathways important in adaptive and innate cell function remained upregulated at the long-term timepoint after resolution of CMV DNAemia. Differential expression analysis and time-course gene set analysis revealed the dynamics of genes and pathways involved in the immune response to CMV DNAemia in kidney transplant patients. Understanding transcriptional changes caused by CMV DNAemia may identify the mechanism behind patient vulnerability to CMV reactivation and increased risk of rejection in transplant recipients and suggest protective strategies to counter the negative immunologic impact of CMV. These findings provide a framework to identify immune correlates for risk assessment and guiding need for extending antiviral prophylaxis.

Diabetic kidney disease (DKD) is a key microvascular complication of diabetes, with few therapies for targeting renal disease pathogenesis and progression. We performed transcriptional and protein studies on 103 unique blood and kidney tissue samples from patients with and without diabetes to understand the pathophysiology of DKD injury and its progression. The study was based on the use of 3 unique patient cohorts: peripheral blood mononuclear cell (PBMC) transcriptional studies were conducted on 30 patients with DKD with advancing kidney injury; Gene Expression Omnibus (GEO) data was downloaded, containing transcriptional measures from 51 microdissected glomerulous from patients with DKD. Additionally, 12 independent kidney tissue sections from patients with or without DKD were used for validation of target genes in diabetic kidney injury by kidney tissue immunohistochemistry and immunofluorescence. PBMC DKD transcriptional analysis, identified 853 genes (p < 0.05) with increasing expression with progression of albuminuria and kidney injury in patients with diabetes. GEO data was downloaded, normalized, and analyzed for significantly changed genes. Of the 325 significantly up regulated genes in DKD glomerulous (p < 0.05), 28 overlapped in PBMC and diabetic kidney, with perturbed FcER1 signaling as a significantly enriched canonical pathway. FcER1 was validated to be significantly increased in advanced DKD, where it was also seen to be specifically co-expressed in the kidney biopsy with tissue mast cells. In conclusion, we demonstrate how leveraging public and private human transcriptional datasets can discover and validate innate immunity and inflammation as key mechanistic pathways in DKD progression, and uncover FcER1 as a putative new DKD target for rational drug design.

Kidney fibrosis constitutes the shared final pathway of nearly all chronic nephropathies, but biomarkers for the non-invasive assessment of kidney fibrosis are currently not available. To address this, we characterize five candidate biomarkers of kidney fibrosis: Cadherin-11 (CDH11), Sparc-related modular calcium binding protein-2 (SMOC2), Pigment epithelium-derived factor (PEDF), Matrix-Gla protein, and Thrombospondin-2. Gene expression profiles in single-cell and single-nucleus RNA-sequencing (sc/snRNA-seq) datasets from rodent models of fibrosis and human chronic kidney disease (CKD) were explored, and Luminex-based assays for each biomarker were developed. Plasma and urine biomarker levels were measured using independent prospective cohorts of CKD: the Boston Kidney Biopsy Cohort, a cohort of individuals with biopsyconfirmed semiquantitative assessment of kidney fibrosis, and the Seattle Kidney Study, a cohort of patients with common forms of CKD. Ordinal logistic regression and Cox proportional hazards regression models were used to test associations of biomarkers with interstitial fibrosis and tubular atrophy and progression to end-stage kidney disease and death, respectively. Sc/snRNA-seq data confirmed cell-specific expression of biomarker genes in fibroblasts. After multivariable adjustment, higher levels of plasma CDH11, SMOC2, and PEDF and urinary CDH11 and PEDF were significantly associated with increasing severity of interstitial fibrosis and tubular atrophy in the Boston Kidney Biopsy Cohort. In both cohorts, higher levels of plasma and urinary SMOC2 and urinary CDH11 were independently associated with progression to end-stage kidney disease. Higher levels of urinary PEDF associated with end-stage kidney disease in the Seattle Kidney Study, with a similar signal in the Boston Kidney Biopsy Cohort, although the latter narrowly missed statistical significance. Thus, we identified CDH11, SMOC2, and PEDF as promising non-invasive biomarkers of kidney fibrosis.

Maintenance of systemic homeostasis by kidney requires the coordinated response of diverse cell types. The use of single-cell RNA sequencing (scRNAseq) for patient tissue samples remains fraught with difficulties with cell isolation, purity, and experimental bias. The ability to characterize immune and parenchymal cells during transplant rejection will be invaluable in defining transplant pathology where tissue availability is restricted to needle biopsy fragments. Herein, we present feasibility data for multiplexing approach for droplet scRNAseq (Mux-Seq). Mux-Seq has the potential to minimize experimental batch bias and variation even with very small sample input. In this first proof-of-concept study for this approach, explant tissues from six normal and two transplant recipients after multiple early post-transplant rejection episodes leading to nephrectomy due to aggressive antibody mediated rejection, were pooled for Mux-Seq. A computational tool, Demuxlet was applied for demultiplexing the individual cells from the pooled experiment. Each sample was also applied individually in a single microfluidic run (singleplex) to correlate results with the pooled data from the same sample. Our applied protocol demonstrated that data from Mux-Seq correlated highly with singleplex (Pearson coefficient 0.982) sequencing results, with the ability to identify many known and novel kidney cell types including different infiltrating immune cells. Trajectory analysis of proximal tubule and endothelial cells demonstrated separation between healthy and injured kidney from transplant explant suggesting evolving stages of cell-specific differentiation in alloimmune injury. This study provides the technical groundwork for understanding the pathogenesis of alloimmune injury and host tissue response in transplant rejection and normal human kidney and provides a protocol for optimized processing precious and low input human kidney biopsy tissue for larger scale studies.

Rejection remains the main cause of premature graft loss after kidney transplantation, despite the use of potent immunosuppression. This highlights the need to better understand the composition and the cell-to-cell interactions of the alloreactive inflammatory infiltrate. Here, we performed droplet-based single-cell RNA sequencing of 35,152 transcriptomes from 16 kidney transplant biopsies with varying phenotypes and severities of rejection and without rejection, and identified cell-type specific gene expression signatures for deconvolution of bulk tissue. A specific association was identified between recipient-derived FCGR3A+ monocytes, FCGR3A+ NK cells and the severity of intragraft inflammation. Activated FCGR3A+ monocytes overexpressed CD47 and LILR genes and increased paracrine signaling pathways promoting T cell infiltration. FCGR3A+ NK cells overexpressed FCRL3, suggesting that antibody-dependent cytotoxicity is a central mechanism of NK-cell mediated graft injury. Multiplexed immunofluorescence using 38 markers on 18 independent biopsy slides confirmed this role of FcγRIII+ NK and FcγRIII+ nonclassical monocytes in antibody-mediated rejection, with specificity to the glomerular area. These results highlight the central involvement of innate immune cells in the pathogenesis of allograft rejection and identify several potential therapeutic targets that might improve allograft longevity.