Viruses and Immunity in Transplant PatientsView this Special Issue
Molecular Networks Involved in the Immune Control of BK Polyomavirus
BK polyomavirus infection is the important cause of virus-related nephropathy following kidney transplantation. BK virus reactivates in 30%–80% of kidney transplant recipients resulting in BK virus-related nephropathy in 1%–10% of cases. Currently, the molecular processes associated with asymptomatic infections in transplant patients infected with BK virus remain unclear. In this study we evaluate intrarenal molecular processes during different stages of BKV infection. The gene expression profiles of 90 target genes known to be associated with immune response were evaluated in kidney graft biopsy material using TaqMan low density array. Three patient groups were examined: control patients with no evidence of BK virus reactivation (), infected asymptomatic patients (), and patients with BK virus nephropathy (). Analysis of biopsies from asymptomatic viruria patients resulted in the identification of 5 differentially expressed genes (CD3E, CD68, CCR2, ICAM-1, and SKI) (), and functional analysis showed a significantly heightened presence of costimulatory signals (e.g., CD40/CD40L; ). Gene ontology analysis revealed several biological networks associated with BKV immune control in comparison to the control group. This study demonstrated that asymptomatic BK viruria is associated with a different intrarenal regulation of several genes implicating in antiviral immune response.
Innovations to immunosuppressive regimens have improved patient and kidney transplant survival rates; however, drug-induced immune suppression has also resulted in significant increases in complications associated with infections. BK polyomavirus (BKV) infections have emerged as an important cause of virus-related nephropathy following kidney transplantation in the era of modern immunosuppressive therapies . BKV has been shown to reactivate in 30%–80% of kidney transplant recipients but only in 1%–10% of cases resulted in the development of BKV nephropathy (BKVN) associated with subsequent kidney graft deterioration and failure [2–5].
Recently, the polymerase chain reaction (PCR) has been used for routine monitoring of BKV replication in peripheral blood . However, PCR screening has demonstrated that a majority of kidney transplant recipients are BKV positive in urine but not blood and never develop BKV nephropathy with graft function deterioration. Furthermore, it was shown that patients with asymptomatic viruria presented with significant BKV viral loads in kidney graft biopsy specimens , suggesting successful control of the BKV infection by the host immune system.
Generally, innate immunity and nonspecific effector mechanisms represent a first line of defense against infection, followed by the development of more specific, acquired immune responses. Active viral replication in tissues has been considered to be a trigger point for inflammation, and cellular immunity has been suggested to play a critical role in viral control. Recently, it was shown that patients with self-limited BKV reactivation without therapeutic intervention developed BKV-specific T cells and cleared BKV rapidly compared to patients suffering from BKV associated-nephropathy .
However, the process of successful BKV control on molecular level has not been fully described.
In this study, we investigated the intrarenal specific transcripts during a different outcome of BKV infection in kidney graft.
2. Materials and Methods
Based on the BKV monitoring study  and histopathological archive of the Institute for Clinical and Experimental Medicine, Prague, different patient groups were established as follows. (i) The asymptomatic viruria group (AV) () presented with BK viral urine loads higher than 107/mL  on the day of protocol biopsy (3 months after transplant). These patients never developed BKV-associated nephropathy and were negative for BKV replication screening in both urine and blood 12 months after transplantation. Protocol kidney graft biopsies revealed normal findings free of rejection and patients had stable graft function. (ii) Control group (NCG) () consisted of patients who were neither BKV positive in urine nor in serum at any time-point after kidney transplantation. None of the patients exhibited allograft rejection; protocol biopsy was normal and patients had stable f graft function. (iii) The BK virus nephropathy (BKVN) group consisted of patients () with BKV-associated nephropathy. Eight patients presented with BKV-associated nephropathy prior to the start of the study and archived biopsy specimens for molecular biology were used for analysis. In 2 cases, BKV-associated nephropathy developed during the course of the study. Histological confirmation of BKV-associated nephropathy was defined as detection of viral cytopathic changes with intranuclear inclusion bodies, associated renal tubular epithelial cell injury, including tubular epithelial cell necrosis, and denudation of basement membranes, as well as positive immunohistochemical staining for the SV40 T large antigen. Patients’ clinical and demographic variables are to be found in Table 1.
Three months after transplantation, protocol kidney graft biopsies were performed in all patients from group I and II and a part of the biopsy specimen fixed for molecular biology analysis at a later date. The study protocol was approved by the Ethics Committee of the Institute for Clinical and Experimental Medicine in Prague and a written informed consent was obtained from all patients.
2.2. RNA Isolation and TaqMan Low Density Array (TLDA)
Small portions (~2 mm) of the cortical or juxtamedullary zone from biopsy specimens were immediately stored in RNA later (Ambion Corporation, Austin, TX). Renal tissues were homogenized; total RNA were extracted using StrataPrep Total RNA Microprep Kit (Stratagene, La Jolla, CA, USA) and reverse transcribed into complementary DNA (cDNA) as described elsewhere .
The gene expression profile of 90 candidate gene targets known to play roles in the elicitation of immune responses (e.g., genes involved in cytokine expression, costimulatory molecules, growth factors, chemokines, immune regulation, apoptosis markers, and ischemia markers) was determined using real-time RT-PCR () with GAPDH as internal control and cDNA from a control kidney serving as the calibrator in 30 renal biopsy specimen analyses. All evaluated genes are described in Table S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2012/972102. Each immune TLDA profile contains lyophilized gene expression reagents (primers and probes (FAM labeled)) in a preconfigured 384 well format. Two samples in duplicate were analyzed per card. Each loading port was filled with a 100 μL cDNA, nuclease free water, and 2X TaqMan universal PCR master mix. Following centrifugation, cards were sealed with a TLDA sealer (Applied Biosystems, Foster City, CA) to prevent cross-contamination. RT-PCR amplification was performed using an ABI Prism 7900 H.T. Sequence Detection system (Applied Biosystems). TLDA cards were analyzed as relative quantification (RQ) and RQ manager 1.2. software for automated data analysis was used (Applied Biosystems).
2.3. Statistical and Functional Analyses
Continuous variables were presented as the mean ± SD (standard deviation). Statistical analysis of categorical characteristics was performed using the test, of continuous parametric and nonparametric variables using student’s -test and Mann-Whitney tests. The Bonferroni correction was used when appropriate. Supervised hierarchical clustering was performed using the MeV () software in order to visualize results. Gene expression data were compared using the Mann-Whitney test followed by the Bonferroni correction.
For functional analysis, large-scale data management was used to identify specific transcript patterns. In the set analysis, genes were grouped into sets determined by their annotation and then compared between defined groups. For the functional analysis, 2 types of set analyses were used. Since a limited number of genes were assessed, the fully coupled flux analysis that reflects the analysis of the pathway fragment was used as the first set analysis . Fully coupled flux represents a gene network that corresponds to a pathway in which non-zero flux for one reaction implies a nonzero flux for a second reaction, and vice versa. This flux represents the strongest qualitative connectivity that can be identified in a network. The genes coupled by their enzymatic fluxes were shown to have similar expression patterns, share transcriptional regulators, and frequently reside in the same operon. The second gene set type analysis included genes sharing a common gene ontology . For the use of functional analysis, gene expression data were log transformed and compared between groups with one-way analysis of variance (ANOVA) followed by the Bonferroni correction to appropriately account for multiple comparisons. Set-level analysis was carried out using XGENE.ORG [13, 14].
3.1. Gene Expression Profiles
We first determined whether the 90 candidate target transcripts identified differed in their levels of expression between the groups. Using these criteria, cohorts of patients with asymptomatic viruria and the negative control group were established; we also examined archived biopsies from patients with histologically proven BKV-associated nephropathy.
Using hierarchical clustering, different gene transcript profiles were identified among the study groups. Hierarchical clustering demonstrated that gene transcript profiles expressed in kidney graft tissues in patients not testing BKV positive were similar to profiles observed in patients presenting with asymptomatic BKV viruria who did not present with BKV positive blood samples or BKV-associated nephropathy. However, gene transcript profiles identified in patients with BKV-associated nephropathy formed a significantly different cluster (Figure 1).
3.2. Characterization of Intrarenal Gene Transcript Profiles in Patients with Asymptomatic BK Viruria
In order to identify the nature of intrarenal immune mechanisms associated with the control of BKV infections, intrarenal graft transcripts from patients with asymptomatic viruria or patients in the negative control group were analyzed. Five differentially expressed genes were identified, including upregulation of the T cell (CD3E) and macrophage (CD68) markers, in addition to genes encoding the receptor for the monocyte chemoattractant protein-1 (CCR2; a chemokine which specifically mediates monocyte chemotaxis) and the adhesion molecule, ICAM-1. The SKI protooncogene, involved in downregulation of TGF-β1 gene transcripts, was significantly downregulated in the asymptomatic BKV infection group () (Table 2).
Functional analysis revealed that patients with asymptomatic viruria exhibited significantly higher expression levels of the costimulatory signals CD40/CD40L () compared to the negative control group. In addition, groups of genes sharing a common ontology were analyzed, revealing that several biological networks were involved in BKV immune control. These networks were involved primarily with B cell proliferation (BCL2, CD40, CD40L, IL10), T cell proliferation (CD28, CD3E, IL12A, IL4, PTPRC), transmembrane receptor tyrosine kinase pathways (CD4, CD8A, FN1), proteolysis (ACE, ECE1, GZMB, LTA, REN), protein kinase binding (CD3E, CD4, PTPRC), antiapoptotic processes (BCL2, BCL2L, CCL2, CD40L, FAS, IL10, IL1A), and leukocyte adhesion (CD40L, ICAM1) (Table 3).
3.3. Identification of Intrarenal Gene Transcripts Associated with BKV-Associated Nephropathy
The BKVN group displayed differential regulation of 33/90 genes analyzed. The most differentially regulated genes were CCL2, CCL5, CCR2, CD4, CD68, FASL, GNLY, IL1B, IL2RA, IL8, PRF1, PTPRC, and TNF (all compared to the negative control group). These genes are primarily involved in T cell signaling, chemotaxis, activation, and cytotoxicity (Table 4).
Functional analysis revealed that patients with BKVN had significantly higher expression levels of the flux for the FASL/FAS () and CD28, CD80, CD86 (), signaling molecules associated with apoptosis, and costimulation. Moreover, a broad spectrum of molecular networks were shown to have an identical ontology, that is, genes associated with the positive regulation of NF-κB transcription and intracellular signal transduction (TNF, TGFB1), chemotaxis (CCL3, CCL5, IL18, IL1B, IL8, VEGF), activation of MAPK activity (IKBKB, TGFB1, TNF), and negative regulation of viral genome replication (CD80, IL8). A detailed description of ontology-related genes associated with BKV-associated nephropathy is listed in Table 5.
Polyoma BK virus-associated nephropathy represents the one of the most challenging infectious complications associated with kidney transplantation .
A better understanding of the molecular processes associated with the immune control of BKV infections may facilitate improvement of clinical management strategies . Evaluation of BKV in urine samples in this study was performed in a blinded fashion; therefore no changes in clinical management were carried out based on results obtained  since it would have influenced the transcript expression profiles present in biopsies collected 3 months after transplantation. To the best of our knowledge, this is the first molecular study in such patient cohort. This analysis demonstrated that effective BKV control 3 months after transplantation was associated with gene expression profiles with the potential of affecting cellular immune responses, including B and T cell signaling and anti-apoptotic gene networks whereas, in late BKV-associated nephropathy, the profound gene upregulation in networks covering T cell signaling, chemoattraction, activation, and cytotoxicity along with many other inflammatory networks were detected.
Specifically, 5 genes likely to be associated with the successful control of viral replication were identified. The heightened presence of T cells (CD3E), monocyte macrophages (CD68), and their chemoattractant chemokine CCR2 (as well as presence of adhesion molecule ICAM-1) were described. Previous studies have identified these molecules to be associated with viral infections [17, 18].
It is broadly known that T cell expansion and cytokine production are needed for the generation of effective antiviral immune responses . CD8+ cytotoxic T cells secreting interferon-gamma (IFN-γ) or/and tumor necrosis factor alpha (TNF-α) are important components in mediating host immunity against viral infections and have been shown to play critical roles in BKV clearance . In our study, however, the expression pattern of IFN-γ and TNF-α during asymptomatic viruria was just marginal () that may reflect just a limited burden of immune injury. CD8+ T lymphocyte activation is tightly regulated, especially during primary responses elicited following positive and negative costimulation following BKV infections [20, 21]. In our study, molecules associated with costimulation were consistently upregulated in kidney tissues of asymptomatic viruria patients and also in biopsies from patients with confirmed BKV-associated nephropathy. Therefore, in order to further identify additional genes associated with protection from BKV infection, data from our study was further analyzed by carrying out functional analyses. Since a low number of genes were assessed compared to the number of genes that could be analyzed following a microarray analysis, we focused our research on describing smaller interactional and functional units using a gene ontology approach and network flux that corresponds to a part of the pathway.
Flux for CD40/CD40L costimulatory signal was significantly upregulated in biopsies from patients who successfully controlled BKV infection and presented with asymptomatic viruria. It is well known that activation of CD40 on antigen presenting cells following ligation of CD40L (expressed mainly on CD4+ T lymphocytes) contributes to proinflammatory responses necessary for eradication of infections caused by certain types of pathogens . Studies focused on defining cellular immune responses with the potential of controlling BKV replication determined that the majority of the BKV-specific T cells expressed CD40L (CD154) .
Gene targets uncovered by gene ontology analysis identified genes typically associated with the elicitation of immune responses specific to viral infections where the interplay between B and T cell function and effective cellular proliferation represents a basic protective strategy.
The activation of cellular mechanisms in response to BKV infections represents an injury-repair immune response with fibrosis development as a consequence . In the current study, patients with asymptomatic viruria (who never developed BKV-associated nephropathy) presented with normal graft function at the 36-month followup. This meant that successful immune control of viral infection was likely associated with limited or transient cytokine upregulation compared to BKV-associated nephropathy where the burden of injury initiated fibrosis development.
Currently, there is limited information regarding the molecular processes associated with BKV-associated nephropathy. The upregulation of large number of genes involved in cell cycle and proliferation was shown in vitro in BKV infected primary kidney epithelial cells that suggests stimulatory nature of BKV proteins . In vivo, the analysis of genes upregulated in renal allografts affected by BKV-associated nephropathy identified proinflammatory genes (CD8 and related molecules associated with graft fibrosis) similar to the profile observed during cases of acute rejection; however, expression levels were larger in magnitude . In our study, intraparenchymal upregulation of 33 genes was observed in BKV-associated nephropathy, confirming previous results , as well as demonstrating a similar expression profile to that observed during acute rejection [25–27]. Moreover, using the functional analysis approach, another 50 biological processes were described in kidneys affected by BKV-associated nephropathy. Using hierarchical clustering, gene expression in BKV-associated nephropathy formed clearly different group.
In this study, the quantitative PCR analysis was performed. Compared to microarray-based analyses, this technique was fast and quantitative, and the results were more reliable. More specific tools for the study of specific immune responses associated with BKV infections would include ELISPOT and multiparameter flow cytometry analyses [7, 17, 28, 29].
In conclusion, this study demonstrated that asymptomatic BKV viruria reflecting successful immune system control of viral infections was associated with specific gene transcripts and immune processes, specifically transcripts associated with B lymphocyte signaling and costimulation. Furthermore, the degree of associated immune responses was much higher in patients presenting with BKV-associated nephropathy.
Conflict of Interests
The author declare that there is no conflict of interests.
Multicenter Ethics Committee of the Thomayer Hospital and Institute for Clinical and Experimental Medicine, Prague, approved the study protocol.
The authors would like to thank F. Zelezny and M. Holec for their help with data analysis using the X-gene platform. In addition, the authors are indebted to the patients and nurses for their cooperation and help. This work received a Grant from the Internal Grant Agency from the Ministry of Health (NS/9714/2008 and institutional support 00023001) and MZO 00023001.
The gene expression profile of 90 candidate gene targets known to play a role in the elicitation of immune responses (e.g., genes involved in cytokine expression, costimulatory molecules, growth factors, chemokines, immune regulation, apoptosis and ischemia markers) was determined using real-time RT-PCR.
V. Nickeleit, H. H. Hirsch, I. F. Binet et al., “Polyomavirus infection of renal allograft recipients: from latent infection to manifest disease,” Journal of the American Society of Nephrology, vol. 10, no. 5, pp. 1080–1089, 1999.View at: Google Scholar
T. Schachtner, K. Muller, M. Stein et al., “BK virus-specific immunity kinetics: a predictor of recovery from polyomavirus BK-associated nephropathy,” American Journal of Transplantation, vol. 11, no. 11, pp. 2443–2452, 2011.View at: Google Scholar
C. Platzer, S. Ode-Hakim, P. Reinke, W. D. Docke, R. Ewert, and H. D. Volk, “Quantitative PCR analysis of cytokine transcription patterns in peripheral mononuclear cells after anti-CD3 rejection therapy using two novel multispecific competitor fragments,” Transplantation, vol. 58, no. 2, pp. 264–268, 1994.View at: Google Scholar
M. Holec, F. Zelezny, J. Klema, and J. Tolar, “Integrating multiple-platform expression data through gene set features,” in Proceedings of the 5th International Conference on Bioinformatics Research and Applications (ISBRA '09), 2009.View at: Google Scholar
M. Holec, F. Zelezny, J. Klema, and J. Tolar, “Cross-genome knowledge-based expression data fusion,” in Proceedings of the International Conference on Bioinformatics, Computational Biology, Genomics and Chemoinformatics (BCBGC '09), 2009.View at: Google Scholar
G. van Kooten and J. Banchereau, “CD40-CD40 ligand,” Journal of Leukocyte Biology, vol. 67, no. 1, pp. 2–17, 2000.View at: Google Scholar
P. Pandiyan, J. K. E. Hegel, M. Krueger, D. Quandt, and M. C. Brunner-Weinzierl, “High IFN-γ production of individual CD8 T lymphocytes is controlled by CD152 (CTLA-4),” Journal of Immunology, vol. 178, no. 4, pp. 2132–2140, 2007.View at: Google Scholar