Results of Two Cases of Pig-to-Human Kidney Xenotransplantation.

BACKGROUND: Xenografts from genetically modified pigs have become one of the most promising solutions to the dearth of human organs available for transplantation. The challenge in this model has been hyperacute rejection. To avoid this, pigs have been bred with a knockout of the alpha-1,3-galactosyltransferase gene and with subcapsular autologous thymic tissue. METHODS: We transplanted kidneys from these genetically modified pigs into two brain-dead human recipients whose circulatory and respiratory activity was maintained on ventilators for the duration of the study. We performed serial biopsies and monitored the urine output and kinetic estimated glomerular filtration rate (eGFR) to assess renal function and xenograft rejection. RESULTS: The xenograft in both recipients began to make urine within moments after reperfusion. Over the 54-hour study, the kinetic eGFR increased from 23 ml per minute per 1.73 m2 of body-surface area before transplantation to 62 ml per minute per 1.73 m2 after transplantation in Recipient 1 and from 55 to 109 ml per minute per 1.73 m2 in Recipient 2. In both recipients, the creatinine level, which had been at a steady state, decreased after implantation of the xenograft, from 1.97 to 0.82 mg per deciliter in Recipient 1 and from 1.10 to 0.57 mg per deciliter in Recipient 2. The transplanted kidneys remained pink and well-perfused, continuing to make urine throughout the study. Biopsies that were performed at 6, 24, 48, and 54 hours revealed no signs of hyperacute or antibody-mediated rejection. Hourly urine output with the xenograft was more than double the output with the native kidneys. CONCLUSIONS: Genetically modified kidney xenografts from pigs remained viable and functioning in brain-dead human recipients for 54 hours, without signs of hyperacute rejection. (Funded by Lung Biotechnology.).

Donor and recipient genetics: Implications for the development of posttransplant diabetes mellitus.

Posttransplant diabetes mellitus (PTDM) is a prevalent complication of liver transplantation and is associated with cardiometabolic complications. We studied the consequences of genetic effects of liver donors and recipients on PTDM outcomes, focusing on the diverse genetic pathways related to insulin that play a role in the development of PTDM. One thousand one hundred fifteen liver transplant recipients without a pretransplant diagnosis of type 2 diabetes mellitus (T2D) and their paired donors recruited from 2 transplant centers had polygenic risk scores (PRS) for T2D, insulin secretion, and insulin sensitivity calculated. Among recipients in the highest T2D-PRS quintile, donor T2D-PRS did not contribute significantly to PTDM. However, in recipients with the lowest T2D genetic risk, donor livers with the highest T2D-PRS contributed to the development of PTDM (OR [95% CI] = 3.79 [1.10-13.1], P = .035). Recipient risk was linked to factors associated with insulin secretion (OR [95% CI] = 0.85 [0.74-0.98], P = .02), while donor livers contributed to PTDM via gene pathways involved in insulin sensitivity (OR [95% CI] = 0.86 [0.75-0.99], P = .03). Recipient and donor PRS independently and collectively serve as predictors of PTDM onset. The genetically influenced biological pathways in recipients primarily pertain to insulin secretion, whereas the genetic makeup of donors exerts an influence on insulin sensitivity.

Genetic versus self-reported African ancestry of the recipient and neighborhood predictors of kidney transplantation outcomes in 2 multiethnic urban cohorts.

African American (AA) kidney recipients have a higher risk of allograft rejection and failure compared to non-AAs, but to what extent these outcomes are due to genetic versus environmental effects is currently unknown. Herein, we tested the effects of recipient self-reported race versus genetic proportion of African ancestry (pAFR), and neighborhood socioeconomic status (SES) on kidney allograft outcomes in multiethnic kidney transplant recipients from Columbia University (N = 1083) and the University of Pennsylvania (N = 738). All participants were genotyped with SNP arrays to estimate genetic admixture proportions. US census tract variables were used to analyze the effect of neighborhood factors. In both cohorts, self-reported recipient AA race and pAFR were individually associated with increased risk of rejection and failure after adjustment for known clinical risk factors and neighborhood SES factors. Joint analysis confirmed that self-reported recipient AA race and pAFR were both associated with a higher risk of allograft rejection (AA: HR 1.61 (1.31-1.96), P = 4.05E-06; pAFR: HR 1.90 (1.46-2.48), P = 2.40E-06) and allograft failure (AA: HR 1.52 (1.18-1.97), P = .001; pAFR: HR 1.70 (1.22-2.35), P = .002). Further research is needed to disentangle the role of genetics versus environmental, social, and structural factors contributing to poor transplantation outcomes in kidney recipients of African ancestry.

Immune response after pig-to-human kidney xenotransplantation: a multimodal phenotyping study.

BACKGROUND: Cross-species immunological incompatibilities have hampered pig-to-human xenotransplantation, but porcine genome engineering recently enabled the first successful experiments. However, little is known about the immune response after the transplantation of pig kidneys to human recipients. We aimed to precisely characterise the early immune responses to the xenotransplantation using a multimodal deep phenotyping approach. METHODS: We did a complete phenotyping of two pig kidney xenografts transplanted to decedent humans. We used a multimodal strategy combining morphological evaluation, immunophenotyping (IgM, IgG, C4d, CD68, CD15, NKp46, CD3, CD20, and von Willebrand factor), gene expression profiling, and whole-transcriptome digital spatial profiling and cell deconvolution. Xenografts before implantation, wild-type pig kidney autografts, as well as wild-type, non-transplanted pig kidneys with and without ischaemia-reperfusion were used as controls. FINDINGS: The data collected from xenografts suggested early signs of antibody-mediated rejection, characterised by microvascular inflammation with immune deposits, endothelial cell activation, and positive xenoreactive crossmatches. Capillary inflammation was mainly composed of intravascular CD68+ and CD15+ innate immune cells, as well as NKp46+ cells. Both xenografts showed increased expression of genes biologically related to a humoral response, including monocyte and macrophage activation, natural killer cell burden, endothelial activation, complement activation, and T-cell development. Whole-transcriptome digital spatial profiling showed that antibody-mediated injury was mainly located in the glomeruli of the xenografts, with significant enrichment of transcripts associated with monocytes, macrophages, neutrophils, and natural killer cells. This phenotype was not observed in control pig kidney autografts or in ischaemia-reperfusion models. INTERPRETATION: Despite favourable short-term outcomes and absence of hyperacute injuries, our findings suggest that antibody-mediated rejection in pig-to-human kidney xenografts might be occurring. Our results suggest specific therapeutic targets towards the humoral arm of rejection to improve xenotransplantation results. FUNDING: OrganX and MSD Avenir.

Whole-exome sequencing analyses in a Saudi Ischemic Stroke Cohort reveal association signals, and shows polygenic risk scores are related to Modified Rankin Scale Risk.

Ischemic stroke represents a significant societal burden across the globe. Rare high penetrant monogenic variants and less pathogenic common single nucleotide polymorphisms (SNPs) have been described as being associated with risk of diseases. Genetic studies in Saudi Arabian patients offer a greater opportunity to detect rare high penetrant mutations enriched in these consanguineous populations. We performed whole exome sequencing on 387 ischemic stroke subjects from Saudi Arabian hospital networks with up to 20,230 controls from the Saudi Human Genome Project and performed gene burden analyses of variants in 177 a priori loci derived from knowledge-driven curation of monogenic and genome-wide association studies of stroke. Using gene-burden analyses, we observed significant associations in numerous loci under autosomal dominant and/or recessive modelling. Stroke subjects with modified Rankin Scale (mRSs) above 3 were found to carry greater cumulative polygenic risk score (PRS) from rare variants in stroke genes (standardized PRS mean > 0) compared to the population average (standardized PRS mean = 0). However, patients with mRS of 3 or lower had lower cumulative genetic risk from rare variants in stroke genes (OR (95%CI) = 1.79 (1.29-2.49), p = 0.0005), with the means of standardized PRS at or lower than 0. In conclusion, gene burden testing in Saudi stroke populations reveals a number of statistically significant signals under different disease inheritance models. However, interestingly, stroke subjects with mRS of 3 or lower had lower cumulative genetic risk from rare variants in stroke genes and therefore, determining the potential mRS cutoffs to use for clinical significance may allow risk stratification of this population.

Whole-exome sequencing of a Saudi epilepsy cohort reveals association signals in known and potentially novel loci.

BACKGROUND: Epilepsy, a serious chronic neurological condition effecting up to 100 million people globally, has clear genetic underpinnings including common and rare variants. In Saudi Arabia, the prevalence of epilepsy is high and caused mainly by perinatal and genetic factors. No whole-exome sequencing (WES) studies have been performed to date in Saudi Arabian epilepsy cohorts. This offers a unique opportunity for the discovery of rare genetic variants impacting this disease as there is a high rate of consanguinity among large tribal pedigrees. RESULTS: We performed WES on 144 individuals diagnosed with epilepsy, to interrogate known epilepsy-related genes for known and functional novel variants. We also used an American College of Medical Genetics (ACMG) guideline-based variant prioritization approach in an attempt to discover putative causative variants. We identified 32 potentially causative pathogenic variants across 30 different genes in 44/144 (30%) of these Saudi epilepsy individuals. We also identified 232 variants of unknown significance (VUS) across 101 different genes in 133/144 (92%) subjects. Strong enrichment of variants of likely pathogenicity was observed in previously described epilepsy-associated loci, and a number of putative pathogenic variants in novel loci are also observed. CONCLUSION: Several putative pathogenic variants in known epilepsy-related loci were identified for the first time in our population, in addition to several potential new loci which may be prioritized for further investigation.

Unsupervised mRNA-seq classification of heart transplant endomyocardial biopsies.

BACKGROUND: Endomyocardial biopsy (EMB) is currently considered the gold standard for diagnosing cardiac allograft rejection. However, significant limitations related to histological interpretation variability are well-recognized. We sought to develop a methodology to evaluate EMB solely based on gene expression, without relying on histology interpretation. METHODS: Sixty-four EMBs were obtained from 47 post-heart transplant recipients, who were evaluated for allograft rejection. EMBs were subjected to mRNA sequencing, in which an unsupervised classification algorithm was used to identify the molecular signatures that best classified the EMBs. Cytokine and natriuretic peptide peripheral blood profiling was also performed. Subsequently, we performed gene network analysis to identify the gene modules and gene ontology to understand their biological relevance. We correlated our findings with the unsupervised and histological classifications. RESULTS: Our algorithm classifies EMBs into three categories based solely on clusters of gene expression: unsupervised classes 1, 2, and 3. Unsupervised and histological classifications were closely related, with stronger gene module-phenotype correlations for the unsupervised classes. Gene ontology enrichment analysis revealed processes impacting on the regulation of cardiac and mitochondrial function, immune response, and tissue injury response. Significant levels of cytokines and natriuretic peptides were detected following the unsupervised classification. CONCLUSION: We have developed an unsupervised algorithm that classifies EMBs into three distinct categories, without relying on histology interpretation. These categories were highly correlated with mitochondrial, immune, and tissue injury response. Significant cytokine and natriuretic peptide levels were detected within the unsupervised classification. If further validated, the unsupervised classification could offer a more objective EMB evaluation.

Gut microbiota analyses of inflammatory bowel diseases from a representative Saudi population.

BACKGROUND: Crohn's diseases and ulcerative colitis, both of which are chronic immune-mediated disorders of the gastrointestinal tract are major contributors to the overarching Inflammatory bowel diseases. It has become increasingly evident that the pathological processes of IBDs results from interactions between genetic and environmental factors, which can skew immune responses against normal intestinal flora. METHODS: The aim of this study is to assess and analyze the taxa diversity and relative abundances in CD and UC in the Saudi population. We utilized a sequencing strategy that targets all variable regions in the 16 S rRNA gene using the Swift Amplicon 16 S rRNA Panel on Illumina NovaSeq 6000. RESULTS: The composition of stool 16 S rRNA was analyzed from 219 patients with inflammatory bowel disease and from 124 healthy controls. We quantified the abundance of microbial communities to examine any significant differences between subpopulations of samples. At the genus level, two genera in particular, Veillonella and Lachnoclostridium showed significant association with CD versus controls. There were significant differences between subjects with CD versus UC, with the top differential genera spanning Akkermansia, Harryflintia, Maegamonas and Phascolarctobacterium. Furthermore, statistically significant taxa diversity in microbiome composition was observed within the UC and CD groups. CONCLUSIONS: In conclusion we have shown that there are significant differences in gut microbiota between UC, CD and controls in a Saudi Arabian inflammatory bowel disease cohort. This reinforces the need for further studies in large populations that are ethnically and geographically diverse. In addition, our results show the potential to develop classifiers that may have add additional richness of context to clinical diagnosis of UC and CD with larger inflammatory bowel disease cohorts.

Oral microbiota analyses of paediatric Saudi population reveals signatures of dental caries.

BACKGROUND: Oral microbiome sequencing has revealed key links between microbiome dysfunction and dental caries. However, these efforts have largely focused on Western populations, with few studies on the Middle Eastern communities. The current study aimed to identify the composition and abundance of the oral microbiota in saliva samples of children with different caries levels using machine learning approaches. METHODS: Oral microbiota composition and abundance were identified in 250 Saudi participants with high dental caries and 150 with low dental caries using 16 S rRNA sequencing on a NextSeq 2000 SP flow cell (Illumina, CA) using 250 bp paired-end reads, and attempted to build a classifier using random forest models to assist in the early detection of caries. RESULTS: The ADONIS test results indicate that there was no significant association between sex and Bray-Curtis dissimilarity (p ~ 0.93), but there was a significant association with dental caries status (p ~ 0.001). Using an alpha level of 0.05, five differentially abundant operational taxonomic units (OTUs) were identified between males and females as the main effect along with four differentially abundant OTUs between high and low dental caries. The mean metrics for the optimal hyperparameter combination using the model with only differentially abundant OTUs were: Accuracy (0.701); Matthew's correlation coefficient (0.0509); AUC (0.517) and F1 score (0.821) while the mean metrics for random forest model using all OTUs were:0.675; 0.054; 0.611 and 0.796 respectively. CONCLUSION: The assessment of oral microbiota samples in a representative Saudi Arabian population for high and low metrics of dental caries yields signatures of abundances and diversity.

Gut microbiota analyses of Saudi populations for type 2 diabetes-related phenotypes reveals significant association.

BACKGROUND: Large-scale gut microbiome sequencing has revealed key links between microbiome dysfunction and metabolic diseases such as type 2 diabetes (T2D). To date, these efforts have largely focused on Western populations, with few studies assessing T2D microbiota associations in Middle Eastern communities where T2D prevalence is now over 20%. We analyzed the composition of stool 16S rRNA from 461 T2D and 119 non-T2D participants from the Eastern Province of Saudi Arabia. We quantified the abundance of microbial communities to examine any significant differences between subpopulations of samples based on diabetes status and glucose level. RESULTS: In this study we performed the largest microbiome study ever conducted in Saudi Arabia, as well as the first-ever characterization of gut microbiota T2D versus non-T2D in this population. We observed overall positive enrichment within diabetics compared to healthy individuals and amongst diabetic participants; those with high glucose levels exhibited slightly more positive enrichment compared to those at lower risk of fasting hyperglycemia. In particular, the genus Firmicutes was upregulated in diabetic individuals compared to non-diabetic individuals, and T2D was associated with an elevated Firmicutes/Bacteroidetes ratio, consistent with previous findings. CONCLUSION: Based on diabetes status and glucose levels of Saudi participants, relatively stable differences in stool composition were perceived by differential abundance and alpha diversity measures. However, community level differences are evident in the Saudi population between T2D and non-T2D individuals, and diversity patterns appear to vary from well-characterized microbiota from Western cohorts. Comparing overlapping and varying patterns in gut microbiota with other studies is critical to assessing novel treatment options in light of a rapidly growing T2D health epidemic in the region. As a rapidly emerging chronic condition in Saudi Arabia and the Middle East, T2D burdens have grown more quickly and affect larger proportions of the population than any other global region, making a regional reference T2D-microbiome dataset critical to understanding the nuances of disease development on a global scale.

Joint testing of donor and recipient genetic matching scores and recipient genotype has robust power for finding genes associated with transplant outcomes.

Genetic matching between transplant donor and recipient pairs has traditionally focused on the human leukocyte antigen (HLA) regions of the genome, but recent studies suggest that matching for non-HLA regions may be important as well. We assess four genetic matching scores for use in association analyses of transplant outcomes. These scores describe genetic ancestry distance using identity-by-state, or genetic incompatibility or mismatch of the two genomes and therefore may reflect different underlying biological mechanisms for donor and recipient genes to influence transplant outcomes. Our simulation studies show that jointly testing these scores with the recipient genotype is a powerful method for preliminary screening and discovery of transplant outcome related single nucleotide polymorphisms (SNPs) and gene regions. Following these joint tests with marginal testing of the recipient genotype and matching score separately can lead to further understanding of the biological mechanisms behind transplant outcomes. In addition, we present results of a liver transplant data analysis that shows joint testing can detect SNPs significantly associated with acute rejection in liver transplant.

Early detection of SARS-CoV-2 and other infections in solid organ transplant recipients and household members using wearable devices.

The increasing global prevalence of SARS-CoV-2 and the resulting COVID-19 disease pandemic pose significant concerns for clinical management of solid organ transplant recipients (SOTR). Wearable devices that can measure physiologic changes in biometrics including heart rate, heart rate variability, body temperature, respiratory, activity (such as steps taken per day) and sleep patterns, and blood oxygen saturation show utility for the early detection of infection before clinical presentation of symptoms. Recent algorithms developed using preliminary wearable datasets show that SARS-CoV-2 is detectable before clinical symptoms in >80% of adults. Early detection of SARS-CoV-2, influenza, and other pathogens in SOTR, and their household members, could facilitate early interventions such as self-isolation and early clinical management of relevant infection(s). Ongoing studies testing the utility of wearable devices such as smartwatches for early detection of SARS-CoV-2 and other infections in the general population are reviewed here, along with the practical challenges to implementing these processes at scale in pediatric and adult SOTR, and their household members. The resources and logistics, including transplant-specific analyses pipelines to account for confounders such as polypharmacy and comorbidities, required in studies of pediatric and adult SOTR for the robust early detection of SARS-CoV-2, and other infections are also reviewed.

Genome-wide non-HLA donor-recipient genetic differences influence renal allograft survival via early allograft fibrosis.

Donor-recipient (D-R) differences at human leukocyte antigen (HLA) loci are currently incorporated into organ sharing, allocation and immunosuppression decisions. However, while acute rejection episodes have substantially diminished, progressive histologic damage occurs in allografts and improved long-term survival remains an unrealized goal among kidney recipients. Here we tested the hypothesis that non-HLA dependent, genome-wide D-R genetic differences could contribute to unchecked alloimmunity with histologic and functional consequences, culminating in long-term allograft failure. Genome-wide single nucleotide polymorphism (SNP) array data, excluding the HLA region, was utilized from 385 transplants to study the role of D-R differences upon serial histology and allograft survival. ADMIXTURE analysis was performed to quantitatively estimate ancestry in each D-R pair and PLINK was used to estimate the proportion of genome-shared identity-by-descent (pIBD) between D-R pairs. Subsequently, quantitative measures of recipient ancestry based on non-HLA SNPs was associated with death-censored allograft survival in adjusted Cox models. In D-R pairs of similar ancestry, pIBD was significantly associated with allograft survival independent of HLA mismatches in 224 transplants. Surprisingly, pIBD and recipient ancestry were not associated with clinical or subclinical rejection at any time post-transplant. Significantly, in multivariable analysis, pIBD inversely correlated with vascular intimal fibrosis in 160 biopsies obtained less than one year which in turn was significantly associated with allograft survival. Thus, our novel data show that non-HLA D-R differences associate with early vascular intimal fibrosis and allograft survival.

Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis.

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

Contribution of non-HLA incompatibility between donor and recipient to kidney allograft survival: genome-wide analysis in a prospective cohort.

BACKGROUND: The introduction of HLA matching of donors and recipients was a breakthrough in kidney transplantation. However, half of all transplanted kidneys still fail within 15 years after transplantation. Epidemiological data suggest a fundamental role of non-HLA alloimmunity. METHODS: We genotyped 477 pairs of deceased donors and first kidney transplant recipients with stable graft function at three months that were transplanted between Dec 1, 2005, and April 30, 2015. Genome-wide genetic mismatches in non-synonymous single nucleotide polymorphisms (nsSNPs) were calculated to identify incompatibilities in transmembrane and secreted proteins. We estimated the association between nsSNP mismatch and graft loss in a Cox proportional hazard model, adjusting for HLA mismatch and clinical covariates. Customised peptide arrays were generated to screen for antibodies against genotype-derived mismatched epitopes in 25 patients with biopsy-confirmed chronic antibody-mediated rejection. FINDINGS: 59 268 nsSNPs affecting a transmembrane or secreted protein were analysed. The median number of nsSNP mismatches in immune-accessible transmembrane and secreted proteins between donors and recipients was 1892 (IQR 1850-1936). The degree of nsSNP mismatch was independently associated with graft loss in a multivariable model adjusted for HLA eplet mismatch (HLA-A, HLA-B, HLA-C, HLA-DP, HLA-DQ, and HLA-DR). Each increase by a unit of one IQR had an HR of 1·68 (95% CI 1·17-2·41, p=0·005). 5-year death censored graft survival was 98% in the quartile with the lowest mismatch, 91% in the second quartile, 89% in the third quartile, and 82% in the highest quartile (p=0·003, log-rank test). Customised peptide arrays verified a donor-specific alloimmune response to genetically predicted mismatched epitopes. INTERPRETATION: Genetic mismatch of non-HLA haplotypes coding for transmembrane or secreted proteins is associated with an increased risk of functional graft loss independently of HLA incompatibility. As in HLA alloimmunity, donor-specific alloantibodies can be identified against genotype derived non-HLA epitopes. FUNDING: Austrian Science Fund, WWTF (Vienna Science and Technology Fund), and Ministry of Health of the Czech Republic.

Genomics and Liver Transplantation: Genomic Biomarkers for the Diagnosis of Acute Cellular Rejection.

Acute cellular rejection (ACR) is a common complication in liver transplantation recipients (LTRs), especially within the first 12 months, and it is associated with increased morbidity and mortality. Although abnormalities in standard liver biochemistries may raise the clinical suspicion for ACR, it lacks specificity, and invasive liver biopsies, which are associated with numerous risks, are required for definitive diagnoses. Biomarker discovery for minimally invasive tools for diagnosis and prognostication of ACR after liver transplantation (LT) has become a rapidly evolving field of research with a recent shift in focus to omics-based biomarker discovery. Although none are yet ready to replace the standard of care, there are several promising minimally invasive, blood-derived biomarkers that are under intensive research for the diagnosis of ACR in LTRs. These omics-based biomarkers, encompassing DNA, RNA, proteins, and metabolites, hold tremendous potential. Some are likely to become integrated into ACR diagnostic algorithms to assist clinical decision making with a high degree of accuracy that is cost-effective and reduces or even obviates the need for an invasive liver biopsy.

Polygenic risk score of non-melanoma skin cancer predicts post-transplant skin cancer across multiple organ types.

Polygenic risk scores (PRSs) calculated from genome-wide association studies (GWASs) of non-melanoma skin cancer (NMSC) in a general, non-transplant setting have recently been shown to predict risk of and time to post-renal transplant skin cancer. In this study, we set out to test these findings in a cohort of heart, lung, and liver transplant patients to see whether these scores could be applied across different organ transplant types. Using the PRS from Stapleton et al (2018), PRS was calculated for each sample across a European ancestry heart, lung, and liver transplant cohorts (n = 523) and tested as predictor of time to NMSC post-transplant. The top PRS, squamous cell carcinoma (SCC) pT1 x 10-5 , (n SNPs = 1953), SCC pT1 x 10-6 , and SCC pT1 x 10-6 (n SNPs = 1061) were significantly predictive in the time to NMSC, SCC, and basal cell carcinoma (BCC) analysis across organ (P = .006, .02, and .02, respectively). We observed here a similar direction of effect and effect size [NMSC HR = 1.31(1.08-1.59)] to that in the original discovery study with increased polygenic burden leading to a faster time to developing NMSC. In summary, we found that PRS of NMSC calculated from GWAS of NMSC in non-transplant populations independently replicated in this cohort of heart, lung, and liver transplant.

Pharmacogenomics in kidney transplant recipients and potential for integration into practice.

WHAT IS KNOWN AND OBJECTIVE: Pharmacogenomic biomarkers are now used in many clinical care settings and represent one of the successes of precision medicine. Genetic variants are associated with pharmacokinetic and pharmacodynamic changes leading to medication adverse effects and changes in clinical response. Actionable pharmacogenomic variants are common in transplant recipients and have implications for medications used in transplant, but yet are not broadly incorporated into practice. METHODS: From the Clinical Pharmacogenetics Implementation Consortium and Dutch Pharmacogenetics Working Group guidelines, and PharmGKB databases, 12 pharmacogenomic genes with 30 variants were selected and used to create diplotypes and actionable pharmacogenomic phenotypes. A total of 853 kidney allograft recipients who had genomic information available from a genome-wide association study were included. RESULTS: Each recipient had at least one actionable pharmacogenomic diplotype/phenotype, whereas the majority (58%) had three or four actionable diplotypes/phenotypes and 17.4% had five or more among the 12 genes. The participants carried actionable diplotypes/phenotypes for multiple medications, including tacrolimus, azathioprine, clopidogrel, warfarin, simvastatin, voriconazole, antidepressants and proton-pump inhibitors. WHAT IS NEW AND CONCLUSION: Pharmacogenomic variants are common in transplant recipients, and transplant recipients receive medications that have actionable variants. CLINICAL TRIAL: Genomics of Transplantation, clinicaltrials.gov (NCT01714440).

Polygenic risk score as a determinant of risk of non-melanoma skin cancer in a European-descent renal transplant cohort.

Renal transplant recipients have an increased risk of non-melanoma skin cancer (NMSC) compared to in the general population. Here, we show polygenic risk scores (PRS) calculated from genome-wide association studies (GWAS) of NMSC in a general, nontransplant setting, can predict risk of, and time to posttransplant skin cancer. Genetic variants, reaching predefined P-value thresholds were chosen from published squamous cell carcinoma (SCC) and basal cell carcinoma (BCC) nontransplant GWAS. Using these GWAS, BCC and SCC PRS were calculated for each sample across three European ancestry renal transplant cohorts (n = 889) and tested as predictors of case:control status and time to NMSC posttransplant. BCC PRS calculated at P-value threshold 1 × 10-5 was the most significant predictor of case:control status of NMSC posttransplant (OR = 1.61; adjusted P = .0022; AUC [full model adjusted for clinical predictors and PRS] = 0.81). SCC PRS at P-value threshold 1 × 10-5 was the most significant predictor of time to posttransplant NMSC (adjusted P = 9.39 × 10-7 ; HR = 1.41, concordance [full model] = 0.74). PRS of nontransplant NMSC is predictive of case:control status and time to NMSC posttransplant. These results are relevant to how genomics can risk stratify patients to help develop personalized treatment regimens.

Tacrolimus troughs and genetic determinants of metabolism in kidney transplant recipients: A comparison of four ancestry groups.

Tacrolimus trough and dose requirements vary dramatically between individuals of European and African American ancestry. These differences are less well described in other populations. We conducted an observational, prospective, multicenter study from which 2595 kidney transplant recipients of European, African, Native American, and Asian ancestry were studied for tacrolimus trough, doses, and genetic determinants of metabolism. We studied the well-known variants and conducted a CYP3A4/5 gene-wide analysis to identify new variants. Daily doses, and dose-normalized troughs were significantly different between the four groups (P < .001). CYP3A5*3 (rs776746) was associated with higher dose-normalized tacrolimus troughs in all groups but occurred at different allele frequencies and had differing effect sizes. The CYP3A5*6 (rs10264272) and *7 (rs413003343) variants were only present in African Americans. CYP3A4*22 (rs35599367) was not found in any of the Asian ancestry samples. We identified seven suggestive variants in the CYP3A4/5 genes associated with dose-normalized troughs in Native Americans (P = 1.1 × 10-5 -8.8 × 10-6 ) and one suggestive variant in Asian Americans (P = 5.6 × 10-6 ). Tacrolimus daily doses and dose-normalized troughs vary significantly among different ancestry groups. We identified potential new variants important in Asians and Native Americans. Studies with larger populations should be conducted to assess the importance of the identified suggestive variants.