Background RNA-seq is poised to play a major role in the management of kidney transplant patients. Rigorous definition of housekeeping genes (HKG) is essential for further progress in this field. Using single genes or a limited set HKG is inherently problematic since their expression might be altered by specific diseases in the patients being studied. Methods To generate a HKG set specific for kidney transplantation, we performed RNA-sequencing from renal allograft biopsies collected in a variety of clinical settings. Various normalization methods were applied to identify transcripts that had a coefficient of variation of expression that was below the 2nd percentile across all samples, and the corresponding genes were designated as housekeeping genes. Comparison with transcriptomic data from the Gene Expression Omnibus (GEO) database, pathway analysis and molecular biological functions were utilized to validate the housekeeping genes set. Results We have developed a bioinformatics solution to this problem by using nine different normalization methods to derive large HKG gene sets from a RNA-seq data set of 47,611 transcripts derived from 30 biopsies. These biopsies were collected in a variety of clinical settings, including normal function, acute rejection, interstitial nephritis, interstitial fibrosis/tubular atrophy and polyomavirus nephropathy. Transcripts with coefficient of variation below the 2nd percentile were designated as HKG, and validated by showing their virtual absence in diseased allograft derived transcriptomic data sets available in the GEO. Pathway analysis indicated a role for these genes in maintenance of cell morphology, pyrimidine metabolism, and intracellular protein signaling. Conclusions Utilization of these objectively defined HKG data sets will guard against errors resulting from focusing on individual genes like 18S RNA, actin & tubulin, which do not maintain constant expression across the known spectrum of renal allograft pathology.

This study interrogates the antigen-specificity of inflammatory infiltrates in renal biopsies with BK polyomavirus (BKPyV) viremia (BKPyVM) with or without allograft nephropathy (BKPyVN). PBMC from 5 healthy HLA-A0101 subjects were stimulated by peptides derived from the BKPYV proteome or polymorphic regions of HLA. Next generation sequencing (NGS) of the T-cell receptor (TCR) cDNA was performed on peptide stimulated PBMC and 23 biopsies with T-cell mediated rejection (TCMR) or BKPyVN. Biopsies from patients with BKPyVM or BKVPyVN contained 7.7732 times more alloreactive than virus reactive clones. Biopsies with TCMR also contained BKPyV-specific clones, presumably a manifestation of heterologous immunity. The mean cumulative T-cell clonal frequency was 0.1378 for alloreactive clones and 0.0375 for BKPyV reactive clones. Samples with BKPyVN and TCMR clustered separately in dendrograms of V-family and J-gene utilization patterns. Dendrograms also revealed that V-gene, J-gene, and D-gene usage patterns were a function of HLA type. In conclusion, biopsies with BKPyVN contain abundant allospecific clones that exceed the number of virus reactive clones. The T-cell component of tissue injury in viral nephropathy appears to be mediated primarily by an ‘innocent bystander’ mechanism in which the principal element is secondary T-cell influx triggered by both anti-viral and anti-HLA immunity.

BackgroundRecent work using DNA microarrays has suggested that genes related to DNA replication, RNA polymerase assembly, and pathogen recognition receptors can serve as surrogate tissue biomarkers for polyomavirus BK nephropathy (BKPyVN).MethodsWe have examined this premise by looking for differential regulation of these genes using a different technology platform (RNA-seq) and an independent set 25 biopsies covering a wide spectrum of diagnoses.ResultsRNA-seq could discriminate T cell–mediated rejection from other common lesions seen in formalin fixed biopsy material. However, overlapping RNA-seq signatures were found among all disease processes investigated. Specifically, genes previously reported as being specific for the diagnosis of BKPyVN were found to be significantly upregulated in T cell–mediated rejection, inflamed areas of fibrosis/tubular atrophy, as well as acute tubular injury.ConclusionsIn conclusion, the search for virus specific molecular signatures is confounded by substantial overlap in pathogenetic mechanisms between BKPyVN and nonviral forms of allograft injury. Clinical heterogeneity, overlapping exposures, and different morphologic patterns and stage of disease are a source of substantial variability in “Omics” experiments. These variables should be better controlled in future biomarker studies on BKPyVN, T cell–mediated rejection, and other forms of allograft injury, before widespread implementation of these tests in the transplant clinic.