ORIGINAL ARTICLE
Single nucleotide variant counts computed from RNA sequencing and cellular traffic into human kidney allografts

https://doi.org/10.1111/ajt.14870Get rights and content
Under a Creative Commons license
open archive

Advances in bioinformatics allow identification of single nucleotide polymorphisms (variants) from RNA sequence data. In an allograft biopsy, 2 genomes contribute to the RNA pool, 1 from the donor organ and the other from the infiltrating recipient’s cells. We hypothesize that imbalances in genetic variants of RNA sequence data of kidney allograft biopsies provide an objective measure of cellular infiltration of the allograft. We performed mRNA sequencing of 40 kidney allograft biopsies, selected to represent a comprehensive range of diagnostic categories. We analyzed the sequencing reads of these biopsies and of 462 lymphoblastoid cell lines from the 1000 Genomes Project, for RNA variants. The ratio of heterozygous to nonreference genome homozygous variants (Het/Hom ratio) on all autosomes was determined for each sample, and the estimation of stromal and immune cells in malignant tumors using expression data (ESTIMATE) score was computed as a complementary estimate of the degree of cellular infiltration into biopsies. The Het/Hom ratios (P = .02) and the ESTIMATE scores (P < .001) were associated with the biopsy diagnosis. Both measures correlated significantly (r = .67, P < .0001), even though the Het/Hom ratio is based on mRNA sequence variation, while the ESTIMATE score uses mRNA expression. Het/Hom ratio and the ESTIMATE score may offer unbiased and quantitative parameters for characterizing cellular traffic into human kidney allografts.

KEYWORDS

genomics
kidney transplantation/nephrology
molecular biology: mRNA/mRNA expression
monitoring: immune
translational research/science

Abbreviations

Active ABMR
active antibody-mediated rejection
acute TCMR
acute T cell–mediated rejection
BAM
binary alignment/map
CD
cluster of differentiation
cDNA
complementary deoxyribonucleic acid
chronic active ABMR
chronic active antibody-mediated rejection
dbSNP 138
single nucleotide polymorphism database build 138
DNA
deoxyribonucleic acid
ESTIMATE
estimation of stromal and immune cells in malignant tumors using expression data
FS
Fisher Strand value
GATK
genome analysis tool kit
gEUVADIS
Genetic European Variation in Health and Disease
GRCh37
Genome Reference Consortium human genome build 37
Het/Hom ratio
heterozygosity to nonreference genome homozygosity ratio
IFTA
interstitial fibrosis and tubular atrophy
LCL
lymphoblastoid cell line
mRNA
messenger ribonucleic acid
RNA
ribonucleic acid
SNP
single nucleotide polymorphism
STAR aligner
Spliced Transcripts Alignment to a Reference aligner

Cited by (0)