Characterization of ALTO-encoding circular RNAs expressed by Merkel cell polyomavirus and trichodysplasia spinulosa polyomavirus.

Circular RNAs (circRNAs) are a conserved class of RNAs with diverse functions, including serving as messenger RNAs that are translated into peptides. Here we describe circular RNAs generated by human polyomaviruses (HPyVs), some of which encode variants of the previously described alternative large T antigen open reading frame (ALTO) protein. Circular ALTO RNAs (circALTOs) can be detected in virus positive Merkel cell carcinoma (VP-MCC) cell lines and tumor samples. CircALTOs are stable, predominantly located in the cytoplasm, and N6-methyladenosine (m6A) modified. The translation of MCPyV circALTOs into ALTO protein is negatively regulated by MCPyV-generated miRNAs in cultured cells. MCPyV ALTO expression increases transcription from some recombinant promoters in vitro and upregulates the expression of multiple genes previously implicated in MCPyV pathogenesis. MCPyV circALTOs are enriched in exosomes derived from VP-MCC lines and circALTO-transfected 293T cells, and purified exosomes can mediate ALTO expression and transcriptional activation in MCPyV-negative cells. The related trichodysplasia spinulosa polyomavirus (TSPyV) also expresses a circALTO that can be detected in infected tissues and produces ALTO protein in cultured cells. Thus, human polyomavirus circRNAs are expressed in human tumors and infected tissues and express proteins that have the potential to modulate the infectious and tumorigenic properties of these viruses.

Humoral response to SARS-CoV-2 infection among liver transplant recipients.

OBJECTIVE: Immunosuppressive agents are known to interfere with T and/or B lymphocytes, which are required to mount an adequate serologic response. Therefore, we aim to investigate the antibody response to SARS-CoV-2 in liver transplant (LT) recipients after COVID-19. DESIGN: Prospective multicentre case-control study, analysing antibodies against the nucleocapsid protein, spike (S) protein of SARS-CoV-2 and their neutralising activity in LT recipients with confirmed SARS-CoV-2 infection (COVID-19-LT) compared with immunocompetent patients (COVID-19-immunocompetent) and LT recipients without COVID-19 symptoms (non-COVID-19-LT). RESULTS: Overall, 35 LT recipients were included in the COVID-19-LT cohort. 35 and 70 subjects fulfilling the matching criteria were assigned to the COVID-19-immunocompetent and non-COVID-19-LT cohorts, respectively. We showed that LT recipients, despite immunosuppression and less symptoms, mounted a detectable antinucleocapsid antibody titre in 80% of the cases, although significantly lower compared with the COVID-19-immunocompetent cohort (3.73 vs 7.36 index level, p<0.001). When analysing anti-S antibody response, no difference in positivity rate was found between the COVID-19-LT and COVID-19-immunocompetent cohorts (97.1% vs 100%, p=0.314). Functional antibody testing showed neutralising activity in 82.9% of LT recipients (vs 100% in COVID-19-immunocompetent cohort, p=0.024). CONCLUSIONS: Our findings suggest that the humoral response of LT recipients is only slightly lower than expected, compared with COVID-19 immunocompetent controls. Testing for anti-S antibodies alone can lead to an overestimation of the neutralising ability in LT recipients. Altogether, routine antibody testing against separate SARS-CoV-2 antigens and functional testing show that the far majority of LT patients are capable of mounting an adequate antibody response with neutralising ability.

Rotavirus-Induced Neonatal Epileptic Encephalopathy-A Disease Spectrum Illustrated by Monochorionic Twins.

Rotavirus has been associated with neonatal seizures and specific white matter magnetic resonance imaging (MRI) abnormalities. We describe monochorionic twins who not only tested positive for rotavirus with these white matter MRI abnormalities but who also showed an electroencephalogram (EEG) pattern characteristic of early infantile epileptic encephalopathy (EIEE), which has so far solely been described in epileptic encephalopathies with a poor prognosis. This report suggests that rotavirus infection must be added to the list of causes of EIEE EEG, and that the outcome then is likely more favorable. As MRI and EEG signs of rotavirus encephalopathy were present in one twin with only subtle neurologic symptoms, rotavirus may well cause insidious central nervous system complications more often. We suggest considering rotavirus infection in neonates presenting with seizures, and to add rotavirus infection to the differential diagnosis of EIEE.

Impact of HPyV9 and TSPyV coinfection on the development of BK polyomavirus viremia and associated nephropathy after kidney transplantation.

BACKGROUND: BK polyomavirus (BKPyV) persistently infects the urinary tract and causes viremia and nephropathy in kidney transplantation (KTx), recipients. In a previous study, we observed an increased incidence and load of BKPyV viremia in KTx patients coinfected with human polyomavirus 9 (HPyV9). Here we sought confirmation of this observation and explored whether novel HPyVs that have been detected in urine (HPyV9 and trichodysplasia spinulosa polyomavirus [TSPyV]) potentially aggravate BKPyV infection. METHODS: A well-characterized cohort of 209 KTx donor-recipient pairs was serologically and molecularly analyzed for HPyV9 and TSPyV coinfection. These data were correlated with the occurrence of BKPyV viremia and BKPyVAN in the recipients within a year after KTx. RESULTS: Seropositivity for HPyV9 (19%) and TSPyV (89%) was comparable between donors and recipients and did not correlate with BKPyV viremia and BKPyVAN that developed in 25% and 3% of the recipients, respectively. Two recipients developed TSPyV viremia and none HPyV9 viremia. Modification of the predictive effect of donor BKPyV seroreactivity on recipient BKPyV viremia by HPyV9 and TSPyV was not observed. CONCLUSIONS: Our data provide no evidence for a promoting effect of HPyV9 and TSPyV on BKPyV infection and BKPyVAN in renal allograft patients. Therefore, we do not recommend including HPyV9 and TSPyV screening in KTx patients.

Primary trichodysplasia spinulosa polyomavirus infection in a kidney transplant child displaying virus-infected decoy cells in the urine.

We report a case of primary trichodysplasia spinulosa (TS) infection in a kidney transplant child and describe for the first time the presence of degenerated TS-associated polyomavirus (TSPyV)-infected cells in a TS patient's urine that are morphologically different from BK or JC polyomavirus-infected decoy cells.

Human papillomavirus and posttransplantation cutaneous squamous cell carcinoma: A multicenter, prospective cohort study.

Organ transplant recipients (OTRs) have a 100-fold increased risk of cutaneous squamous cell carcinoma (cSCC). We prospectively evaluated the association between ß genus human papillomaviruses (ßPV) and keratinocyte carcinoma in OTRs. Two OTR cohorts without cSCC were assembled: cohort 1 was transplanted in 2003-2006 (n = 274) and cohort 2 was transplanted in 1986-2002 (n = 352). Participants were followed until death or cessation of follow-up in 2016. ßPV infection was assessed in eyebrow hair by using polymerase chain reaction-based methods. ßPV IgG seroresponses were determined with multiplex serology. A competing risk model with delayed entry was used to estimate cumulative incidence of histologically proven cSCC and the effect of ßPV by using a multivariable Cox regression model. Results are reported as adjusted hazard ratios (HRs). OTRs with 5 or more different ßPV types in eyebrow hair had 1.7 times the risk of cSCC vs OTRs with 0 to 4 different types (HR 1.7, 95% confidence interval 1.1-2.6). A similar risk was seen with high ßPV loads (HR 1.8, 95% confidence interval 1.2-2.8). No significant associations were seen between serum antibodies and cSCC or between ßPV and basal cell carcinoma. The diversity and load of ßPV types in eyebrow hair are associated with cSCC risk in OTRs, providing evidence that ßPV is associated with cSCC carcinogenesis and may present a target for future preventive strategies.

Development and Evaluation of a Broad Bead-Based Multiplex Immunoassay To Measure IgG Seroreactivity against Human Polyomaviruses.

The family of polyomaviruses, which cause severe disease in immunocompromised hosts, has expanded substantially in recent years. To accommodate measurement of IgG seroresponses against all currently known human polyomaviruses (HPyVs), including the Lyon IARC polyomavirus (LIPyV), we extended our custom multiplex bead-based HPyV immunoassay and evaluated the performance of this pan-HPyV immunoassay. The VP1 proteins of 15 HPyVs belonging to 13 Polyomavirus species were expressed as recombinant glutathione S-transferase (GST) fusion proteins and coupled to fluorescent Luminex beads. Sera from healthy blood donors and immunocompromised kidney transplant recipients were used to analyze seroreactivity against the different HPyVs. For BK polyomavirus (BKPyV), the GST-VP1 fusion protein-directed seroresponses were compared to those obtained against BKPyV VP1 virus-like particles (VLP). Seroreactivity against most HPyVs was common and generally high in both test populations. Low seroreactivity against HPyV9, HPyV12, New Jersey PyV, and LIPyV was observed. The assay was reproducible (Pearson's r2 > 0.84, P < 0.001) and specific. Weak but consistent cross-reactivity between the related viruses HPyV6 and HPyV7 was observed. The seroresponses measured by the GST-VP1-based immunoassay and a VP1 VLP-based enzyme-linked immunosorbent assay were highly correlated (Spearman's ρ = 0.823, P < 0.001). The bead-based pan-HPyV multiplex immunoassay is a reliable tool to determine HPyV-specific seroresponses with high reproducibility and specificity and is suitable for use in seroepidemiological studies.

Clinically Relevant Reactivation of Polyomavirus BK (BKPyV) in HLA-A02-Positive Renal Transplant Recipients Is Associated with Impaired Effector-Memory Differentiation of BKPyV-Specific CD8+ T Cells.

Polyomavirus BK (BKPyV) frequently reactivates in immunosuppressed renal transplant recipients (RTRs) and may lead to graft loss due to BKPyV-induced interstitial nephritis (BKVN). Little is known on the differentiation of CD8+ T cells targeting BKPyV in RTRs. Here we investigated whether BKPyV-specific CD8+ T cell differentiation differs in RTRs with varying degrees of BKPyV reactivation and/or BKVN. Using combinatorial encoding with tetramers carrying BKPyV major capsid protein (VP1) and large T antigen protein (LTAG) epitopes, we investigated CD8+ T cell responses to BKPyV in longitudinally obtained PBMC samples from 46 HLA-A02-positive RTRs and 20 healthy adults. We were also able to isolate BKPyV-specific CD8+ T cells from five renal allografts, two of which were affected by BKVN. Before transplantation, BKPyV-specific CD8+ T cells targeting VP1 and LTAG epitopes appeared predominantly as central-memory and CD27+/CD28+ effector-memory (TEM), and naïve-like PD-1-expressing cells, respectively. After viral reactivation, BKPyV-specific CD8+ T cells assumed CD28- TEM and TEMRA states in patients who were able to control BKPyV, whereas differentiation lagged behind in patients with severe viral reactivation or BKVN. Furthermore, VP1-specific CD69+/CD103+ tissue-resident memory (TRM) cells accumulated in BKVN-affected allografts but lacked signs of effector differentiation. In contrast, granzyme B-expressing effector cells were detected in allografts not affected by BKVN. In conclusion, effector-memory differentiation of BKPyV-specific CD8+ T cells in patients with high viral load or BKVN is impaired. Further characterization of the specific mechanisms behind this altered cellular differentiation is necessary to develop therapies that can prevent the emergence of BKVN.

Clinical implications of rapid ePlex® Respiratory Pathogen Panel testing compared to laboratory-developed real-time PCR.

Rapid diagnosis of respiratory infections is of great importance for adequate isolation and treatment. Due to the batch-wise testing, laboratory-developed real-time polymerase chain reaction (PCR) assays (LDT) often result in a time to result of one day. Here, LDT was compared with rapid ePlex® Respiratory Pathogen (RP) Panel testing of GenMark Diagnostics (Carlsbad, CA, USA) with regard to time to result, installed isolation precautions, and antibacterial/antiviral treatment. Between January and March 2017, 68 specimens of 64 patients suspected of an acute respiratory infection were tested with LDT and the ePlex® RP panel. The time to result was calculated as the time between sample reception and result reporting. Information regarding isolation and antibacterial/antiviral treatment was obtained from the patient records. Thirty specimens tested LDT positive (47%) and 29 ePlex® RP panel positive (45%). The median time to result was 27.1 h (range 6.5-96.6) for LDT versus 3.4 h (range 1.5-23.6) for the RP panel, p-value < 0.001. In 14 out of 30 patients, isolation was discontinued based on the ePlex® RP panel results, saving 21 isolation days. ePlex® RP panel test results were available approximately one day ahead of the LDT results in the 19 patients receiving antiviral/antibacterial treatment. In addition, two bacterial pathogens, not requested by the physician, were detected using the RP panel. Analysis of respiratory infections with the ePlex® RP panel resulted in a significant decrease in time to result, enabling a reduction in isolation days in half of the patients. Furthermore, syndromic RP panel testing increased the identification of causative pathogens.

Interpreting cancer genomes using systematic host network perturbations by tumour virus proteins.

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype-phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations. Genome sequencing efforts have identified numerous germline mutations, and large numbers of somatic genomic alterations, associated with a predisposition to cancer. However, it remains difficult to distinguish background, or 'passenger', cancer mutations from causal, or 'driver', mutations in these data sets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. Here we test the hypothesis that genomic variations and tumour viruses may cause cancer through related mechanisms, by systematically examining host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways, such as Notch signalling and apoptosis, that go awry in cancer. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on a par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches increase the specificity of cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate the prioritization of cancer-causing driver genes to advance the understanding of the genetic basis of human cancer.

Primary Polyomavirus Infection, Not Reactivation, as the Cause of Trichodysplasia Spinulosa in Immunocompromised Patients.

Classic human polyomaviruses (JC and BK viruses) become pathogenic when reactivating from latency. For the rare skin disease trichodysplasia spinulosa, we show that manifestations of the causative polyomavirus (TSPyV) occur during primary infection of the immunosuppressed host. High TSPyV loads in blood and cerebrospinal fluid, sometimes coinciding with cerebral lesions and neuroendocrine symptoms, marked the acute phase of trichodysplasia spinulosa, whereas initiation and maturation of TSPyV seroresponses occurred in the convalescent phase. TSPyV genomes lacked the rearrangements typical for reactivating polyomaviruses. These findings demonstrate the clinical importance of primary infection with this rapidly expanding group of human viruses and explain the rarity of some novel polyomavirus-associated diseases.

ICTV Virus Taxonomy Profile: Polyomaviridae.

The Polyomaviridae is a family of small, non-enveloped viruses with circular dsDNA genomes of approximately 5 kbp. The family includes four genera whose members have restricted host range, infecting mammals and birds. Polyomavirus genomes have also been detected recently in fish. Merkel cell polyomavirus and raccoon polyomavirus are associated with cancer in their host; other members are human and veterinary pathogens. Clinical manifestations are obvious in immunocompromised patients but not in healthy individuals. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the taxonomy of the Polyomaviridae, which is available at www.ictv.global/report/polyomaviridae.

Trichodysplasia spinulosa-Associated Polyomavirus Uses a Displaced Binding Site on VP1 to Engage Sialylated Glycolipids.

Trichodysplasia spinulosa-associated Polyomavirus (TSPyV) was isolated from a patient suffering from trichodysplasia spinulosa, a skin disease that can appear in severely immunocompromised patients. While TSPyV is one of the five members of the polyomavirus family that are directly linked to a human disease, details about molecular recognition events, the viral entry pathway, and intracellular trafficking events during TSPyV infection remain unknown. Here we have used a structure-function approach to shed light on the first steps of TSPyV infection. We established by cell binding and pseudovirus infection studies that TSPyV interacts with sialic acids during attachment and/or entry. Subsequently, we solved high-resolution X-ray structures of the major capsid protein VP1 of TSPyV in complex with three different glycans, the branched GM1 glycan, and the linear trisaccharides α2,3- and α2,6-sialyllactose. The terminal sialic acid of all three glycans is engaged in a unique binding site on TSPyV VP1, which is positioned about 18 Å from established sialic acid binding sites of other polyomaviruses. Structure-based mutagenesis of sialic acid-binding residues leads to reduction in cell attachment and pseudovirus infection, demonstrating the physiological relevance of the TSPyV VP1-glycan interaction. Furthermore, treatments of cells with inhibitors of N-, O-linked glycosylation, and glycosphingolipid synthesis suggest that glycolipids play an important role during TSPyV infection. Our findings elucidate the first molecular recognition events of cellular infection with TSPyV and demonstrate that receptor recognition by polyomaviruses is highly variable not only in interactions with sialic acid itself, but also in the location of the binding site.

Interspecific adaptation by binary choice at de novo polyomavirus T antigen site through accelerated codon-constrained Val-Ala toggling within an intrinsically disordered region.

It is common knowledge that conserved residues evolve slowly. We challenge generality of this central tenet of molecular biology by describing the fast evolution of a conserved nucleotide position that is located in the overlap of two open reading frames (ORFs) of polyomaviruses. The de novo ORF is expressed through either the ALTO protein or the Middle T antigen (MT/ALTO), while the ancestral ORF encodes the N-terminal domain of helicase-containing Large T (LT) antigen. In the latter domain the conserved Cys codon of the LXCXE pRB-binding motif constrains codon evolution in the overlapping MT/ALTO ORF to a binary choice between Val and Ala codons, termed here as codon-constrained Val-Ala (COCO-VA) toggling. We found the rate of COCO-VA toggling to approach the speciation rate and to be significantly accelerated compared to the baseline rate of chance substitution in a large monophyletic lineage including all viruses encoding MT/ALTO and three others. Importantly, the COCO-VA site is located in a short linear motif (SLiM) of an intrinsically disordered region, a typical characteristic of adaptive responders. These findings provide evidence that the COCO-VA toggling is under positive selection in many polyomaviruses, implying its critical role in interspecific adaptation, which is unprecedented for conserved residues.

Characterization of T Antigens, Including Middle T and Alternative T, Expressed by the Human Polyomavirus Associated with Trichodysplasia Spinulosa.

UNLABELLED: The polyomavirus tumor (T) antigens play crucial roles in viral replication, transcription, and cellular transformation. They are encoded by partially overlapping open reading frames (ORFs) located in the early region through alternative mRNA splicing. The T expression pattern of the trichodysplasia spinulosa-associated polyomavirus (TSPyV) has not been established yet, hampering further study of its pathogenic mechanisms and taxonomic relationship. Here, we characterized TSPyV T antigen expression in human cell lines transfected with the TSPyV early region. Sequencing of T antigen-encoded reverse transcription-PCR (RT-PCR) products revealed three splice donor and acceptor sites creating six mRNA splice products that potentially encode the antigens small T (ST), middle T (MT), large T (LT), tiny T, 21kT, and alternative T (ALTO). Except for 21kT, these splice products were also detected in skin of TSPyV-infected patients. At least three splice products were confirmed by Northern blotting, likely encoding LT, MT, ST, 21kT, and ALTO. Protein expression was demonstrated for LT, ALTO, and possibly MT, with LT detected in the nucleus and ALTO in the cytoplasm of transfected cells. Splice site and start codon mutations indicated that ALTO is encoded by the same splice product that encodes LT and uses internal start codons for initiation. The genuineness of ALTO was indicated by the identification of acetylated N-terminal ALTO peptides by mass spectrometry. Summarizing, TSPyV exhibits an expression pattern characterized by both MT and ALTO expression, combining features of rodent and human polyomaviruses. This unique expression pattern provides important leads for further study of polyomavirus-related disease and for an understanding of polyomavirus evolution. IMPORTANCE: The human trichodysplasia spinulosa-associated polyomavirus (TSPyV) is distinguished among polyomaviruses for combining productive infection with cell-transforming properties. In the research presented here, we further substantiate this unique position by indicating expression of both middle T antigen (MT) and alternative T antigen (ALTO) in TSPyV. So far, none of the human polyomaviruses was shown to express MT, which is considered the most important viral oncoprotein of rodent polyomaviruses. Coexpression of ALTO and MT, which involves a conserved, recently recognized overlapping ORF subject to positive selection, has not been observed before for any polyomavirus. As a result of our findings, this study provides valuable new insights into polyomavirus T gene use and expression. Obviously, these insights will be instrumental in further study and gaining an understanding of TSPyV pathogenicity. More importantly, however, they provide important leads with regard to the interrelationship, functionality, and evolution of polyomaviruses as a whole, indicating that TSPyV is a suitable model virus to study these entities further.

A taxonomy update for the family Polyomaviridae.

Many distinct polyomaviruses infecting a variety of vertebrate hosts have recently been discovered, and their complete genome sequence could often be determined. To accommodate this fast-growing diversity, the International Committee on Taxonomy of Viruses (ICTV) Polyomaviridae Study Group designed a host- and sequence-based rationale for an updated taxonomy of the family Polyomaviridae. Applying this resulted in numerous recommendations of taxonomical revisions, which were accepted by the Executive Committee of the ICTV in December 2015. New criteria for definition and creation of polyomavirus species were established that were based on the observed distance between large T antigen coding sequences. Four genera (Alpha-, Beta, Gamma- and Deltapolyomavirus) were delineated that together include 73 species. Species naming was made as systematic as possible - most species names now consist of the binomial name of the host species followed by polyomavirus and a number reflecting the order of discovery. It is hoped that this important update of the family taxonomy will serve as a stable basis for future taxonomical developments.

Human Beta-papillomavirus infection and keratinocyte carcinomas.

Although the role of oncogenic human Alpha-papillomaviruses (HPVs) in the development of mucosal carcinomas at different body sites (eg cervix, anus, oropharynx) is fully recognized, a role for HPV in keratinocyte carcinomas (KCs; basal and squamous cell carcinomas) of the skin is not yet clear. KCs are the most common cancers in Caucasians, with the major risk factor being ultraviolet (UV) light exposure. A possible role for Beta-HPV types (BetaPV) in the development of KC was suggested several decades ago, supported by a number of epidemiological studies. Our current review summarizes the recent molecular and histopathological evidence in support of a causal association between BetaPV and the development of KC, and outlines the suspected synergistic effect of viral gene expression with UV radiation and immune suppression. Further insights into the molecular pathways and protein interactions used by BetaPV and the host cell is likely to extend our understanding of the role of BetaPV in KC.

Human polyomavirus 9 infection in kidney transplant patients.

Several human polyomaviruses of unknown prevalence and pathogenicity have been identified, including human polyomavirus 9 (HPyV9). To determine rates of HPyV9 infection among immunosuppressed patients, we screened serum samples from 101 kidney transplant patients in the Netherlands for HPyV9 DNA and seroreactivity. A total of 21 patients had positive results for HPyV9 DNA; positivity rates peaked at 3 months after transplantation, but the highest viral loads were measured just after transplantation. During 18 months of follow-up, HPyV9 seroprevalence increased from 33% to 46% among transplant patients; seroprevalence remained stable at ≈30% in a control group of healthy blood donors in whom no HPyV9 DNA was detected. Further analysis revealed an association between detection of HPyV9 and detection of BK polyomavirus but not of cytomegalovirus. Our data indicate that HPyV9 infection is frequent in kidney transplant patients, but the nature of infection-endogenous or donor-derived-and pathogenic potential of this virus remain unknown.

Immune responsiveness in stable kidney transplantation patients: Complete inhibition of T-cell proliferation but residual T-cell activity during maintenance immunosuppressive treatment.

The recommended immunosuppressive treatment after kidney transplantation consists of tacrolimus, mycophenolate mofetil, and low-dose corticosteroids. Drug concentrations are monitored using therapeutic drug monitoring (TDM), which does not necessarily correlate with pharmacodynamic activity. To find the balance between optimal efficacy and minimal toxicity, it might be more informative to monitor patients' immunological status rather than drug concentrations. We selected a panel of T-cell-based immune assays, which were used for immunomonitoring of 14 stable kidney transplantation patients. Whole blood was incubated with a T-cell stimulus, after which T-cell proliferation, T-cell activation marker expression and cytokine production were measured to study residual immune activity in vitro (before drug intake; drug added to the incubation) and ex vivo (after drug intake). T-cell proliferation was completely suppressed in all patients over the full day, while IL-2, IFN-γ, CD71, and CD154 showed fluctuations over the day with a strong inhibition (75%-25%) at 2 h post-dose. The level of inhibition was variable between patients and could not be related to pharmacokinetic parameters or the presence of regulatory or senescence immune cells. Moreover, the level of inhibition did not correlate with the in vitro tacrolimus drug effect as studied by incubating pre-dose blood samples with additional tacrolimus. Overall, IL-2, IFN-γ, CD71, and CD154 seem to be good markers to monitor residual immune activity of transplantation patients. To evaluate the correlation between these pharmacodynamic biomarkers and clinical outcome, prospective observational studies are needed.