Extreme phenotype sampling and next generation sequencing to identify genetic variants associated with tacrolimus in African American kidney transplant recipients.

African American (AA) kidney transplant recipients (KTRs) have poor outcomes, which may in-part be due to tacrolimus (TAC) sub-optimal immunosuppression. We previously determined the common genetic regulators of TAC pharmacokinetics in AAs which were CYP3A5 *3, *6, and *7. To identify low-frequency variants that impact TAC pharmacokinetics, we used extreme phenotype sampling and compared individuals with extreme high (n = 58) and low (n = 60) TAC troughs (N = 515 AA KTRs). Targeted next generation sequencing was conducted in these two groups. Median TAC troughs in the high group were 7.7 ng/ml compared with 6.3 ng/ml in the low group, despite lower daily doses of 5 versus 12 mg, respectively. Of 34,542 identified variants across 99 genes, 1406 variants were suggestively associated with TAC troughs in univariate models (p-value < 0.05), however none were significant after multiple testing correction. We suggest future studies investigate additional sources of TAC pharmacokinetic variability such as drug-drug-gene interactions and pharmacomicrobiome.

Extreme Phenotype Sampling and Next Generation Sequencing to Identify Genetic Variants Associated with Tacrolimus in African American Kidney Transplant Recipients.

African American (AA) kidney transplant recipients (KTRs) have poor outcomes, which may in-part be due to tacrolimus (TAC) sub-optimal immunosuppression. We previously determined the common genetic regulators of TAC pharmacokinetics in AAs which were CYP3A5 *3, *6, and *7. To identify low-frequency variants that impact TAC pharmacokinetics, we used extreme phenotype sampling and compared individuals with extreme high (n=58) and low (n=60) TAC troughs (N=515 AA KTRs). Targeted next generation sequencing was conducted in these two groups. Median TAC troughs in the high group were 7.7 ng/ml compared with 6.3 ng/ml in the low group, despite lower daily doses of 5 versus 12mg, respectively. Of 34,542 identified variants across 99 genes, 1,406 variants were suggestively associated with TAC troughs in univariate models (p-value <0.05), however none were significant after multiple testing correction. We suggest future studies investigate additional sources of TAC pharmacokinetic variability such as drug-drug-gene interactions and pharmacomicrobiome.

Identification of genetic variants associated with tacrolimus metabolism in kidney transplant recipients by extreme phenotype sampling and next generation sequencing.

An extreme phenotype sampling (EPS) model with targeted next-generation sequencing (NGS) identified genetic variants associated with tacrolimus (Tac) metabolism in subjects from the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort which included 1,442 European Americans (EA) and 345 African Americans (AA). This study included 48 subjects separated into 4 groups of 12 (AA high, AA low, EA high, EA low). Groups were selected by the extreme phenotype of dose-normalized Tac trough concentrations after adjusting for common genetic variants and clinical factors. NGS spanned > 3 Mb of 28 genes and identified 18,661 genetic variants (3961 previously unknown). A group of 125 deleterious variants, by SIFT analysis, were associated with Tac troughs in EAs (burden test, p = 0.008), CYB5R2 was associated with Tac troughs in AAs (SKAT, p = 0.00079). In CYB5R2, rs61733057 (increased allele frequency in AAs) was predicted to disrupt protein function by SIFT and PolyPhen2 analysis. The variants merit further validation.

Urinary microbiome associated with chronic allograft dysfunction in kidney transplant recipients.

BACKGROUND: We performed a study to identify differences in the urinary microbiome associated with chronic allograft dysfunction (CAD) and compared the urinary microbiome of male and female transplant recipients with CAD. METHODS: This case-control study enrolled 67 patients within the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort at two transplant centers. CAD was defined as a greater than 25% rise in serum creatinine relative to a 3 month post-transplant baseline. Urine samples from patients with and without CAD were analyzed using 16S V4 bacterial ribosomal DNA sequences. RESULTS: Corynebacterium was more prevalent in female and male patients with CAD compared to non-CAD female patients (P = 0.0005). A total 21 distinct Operational Taxonomic Unit (OTUs) were identified as significantly different when comparing CAD and non-CAD patients using Kruskal-Wallis (P < 0.01). A subset analysis of female patients with CAD compared to non-CAD females identified similar differentially abundant OTUs, including the genera Corynebacterium and Staphylococcus (Kruskal-Wallis; P = 0.01; P = 0.004, respectively). Male CAD vs female CAD analysis showed greater abundance of phylum Proteobacteria in males. CONCLUSION: There were differences in the urinary microbiome when comparing female and male CAD patients with their female non-CAD counterparts and these differences persisted in the subset analysis limited to female patients only.

CRISPR/Cas9 Genetic Modification of CYP3A5 *3 in HuH-7 Human Hepatocyte Cell Line Leads to Cell Lines with Increased Midazolam and Tacrolimus Metabolism.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 engineering of the CYP3A5 *3 locus (rs776746) in human liver cell line HuH-7 (CYP3A5 *3/*3) has led to three CYP3A5 *1 cell lines by deletion of the exon 3B splice junction or point mutation. Cell lines CYP3A5 *1/*3 sd (single deletion), CYP3A5 *1/*1 dd (double deletion), or CYP3A5 *1/*3 pm (point mutation) expressed the CYP3A5 *1 mRNA and had elevated CYP3A5 mRNA (P < 0.0005 for all engineered cell lines) and protein expression compared with HuH-7. In metabolism assays, HuH-7 had less tacrolimus (all P < 0.05) or midazolam (MDZ) (all P < 0.005) disappearance than all engineered cell lines. HuH-7 had less 1-OH MDZ (all P < 0.0005) or 4-OH (all P < 0.005) production in metabolism assays than all bioengineered cell lines. We confirmed CYP3A5 metabolic activity with the CYP3A4 selective inhibitor CYP3CIDE. This is the first report of genomic CYP3A5 bioengineering in human cell lines with drug metabolism analysis.

Deceased-Donor Apolipoprotein L1 Renal-Risk Variants Have Minimal Effects on Liver Transplant Outcomes.

BACKGROUND: Apolipoprotein L1 gene (APOL1) G1 and G2 renal-risk variants, common in populations with recent African ancestry, are strongly associated with non-diabetic nephropathy, end-stage kidney disease, and shorter allograft survival in deceased-donor kidneys (autosomal recessive inheritance). Circulating APOL1 protein is synthesized primarily in the liver and hydrodynamic gene delivery of APOL1 G1 and G2 risk variants has caused hepatic necrosis in a murine model. METHODS: To evaluate the impact of these variants in liver transplantation, this multicenter study investigated the association of APOL1 G1 and G2 alleles in deceased African American liver donors with allograft survival. Transplant recipients were followed for liver allograft survival using data from the Scientific Registry of Transplant Recipients. RESULTS: Of the 639 liver donors evaluated, 247 had no APOL1 risk allele, 300 had 1 risk allele, and 92 had 2 risk alleles. Graft failure assessed at 15 days, 6 months, 1 year and total was not significantly associated with donor APOL1 genotype (p-values = 0.25, 0.19, 0.67 and 0.89, respectively). CONCLUSIONS: In contrast to kidney transplantation, deceased-donor APOL1 G1 and G2 risk variants do not significantly impact outcomes in liver transplantation.

Differentially expressed gene transcripts using RNA sequencing from the blood of immunosuppressed kidney allograft recipients.

We performed RNA sequencing (RNAseq) on peripheral blood mononuclear cells (PBMCs) to identify differentially expressed gene transcripts (DEGs) after kidney transplantation and after the start of immunosuppressive drugs. RNAseq is superior to microarray to determine DEGs because it's not limited to available probes, has increased sensitivity, and detects alternative and previously unknown transcripts. DEGs were determined in 32 adult kidney recipients, without clinical acute rejection (AR), treated with antibody induction, calcineurin inhibitor, mycophenolate, with and without steroids. Blood was obtained pre-transplant (baseline), week 1, months 3 and 6 post-transplant. PBMCs were isolated, RNA extracted and gene expression measured using RNAseq. Principal components (PCs) were computed using a surrogate variable approach. DEGs post-transplant were identified by controlling false discovery rate (FDR) at < 0.01 with at least a 2 fold change in expression from pre-transplant. The top 5 DEGs with higher levels of transcripts in blood at week 1 were TOMM40L, TMEM205, OLFM4, MMP8, and OSBPL9 compared to baseline. The top 5 DEGs with lower levels at week 1 post-transplant were IL7R, KLRC3, CD3E, CD3D, and KLRC2 (Striking Image) compared to baseline. The top pathways from genes with lower levels at 1 week post-transplant compared to baseline, were T cell receptor signaling and iCOS-iCOSL signaling while the top pathways from genes with higher levels than baseline were axonal guidance signaling and LXR/RXR activation. Gene expression signatures at month 3 were similar to week 1. DEGs at 6 months post-transplant create a different gene signature than week 1 or month 3 post-transplant. RNAseq analysis identified more DEGs with lower than higher levels in blood compared to baseline at week 1 and month 3. The number of DEGs decreased with time post-transplant. Further investigations to determine the specific lymphocyte(s) responsible for differential gene expression may be important in selecting and personalizing immune suppressant drugs and may lead to targeted therapies.