Genome-Wide Association Study of Tacrolimus Pharmacokinetics Identifies Novel Single Nucleotide Polymorphisms in the Convalescence and Stabilization Periods of Post-transplant Liver Function

Dynamic immune recovery process after liver transplantation revealed by single-cell multi-omics analysis

Elucidating the temporal process of immune remodeling under immunosuppressive treatment after liver transplantation (LT) is critical for precise clinical management strategies. Here, we performed a single-cell multi-omics analysis of peripheral blood mononuclear cells (PBMCs) collected from LT patients (with and without acute cellular rejection [ACR]) at 13 time points. Validation was performed in two independent cohorts with additional LT patients and healthy controls. Our study revealed a four-phase recovery process after LT and delineated changes in immune cell composition, expression programs, and interactions along this process. The intensity of the immune response differs between the ACR and non-ACR patients. Notably, the newly identified inflamed NK cells, CD14+RNASE2+ monocytes, and FOS-expressing monocytes emerged as predictive indicators of ACR. This study illuminates the longitudinal evolution of the immune cell landscape under tacrolimus-based immunosuppressive treatment during LT recovery, providing a four-phase framework that aids the clinical management of LT patients.

Impact of genetic polymorphisms on tacrolimus trough blood concentration in Chinese liver transplant recipients

Purpose: The aim of this study was to analyze the effects of various genetic polymorphisms and clinical factors on tacrolimus (TAC) concentration in the convalescence period (CP) and stabilization period (SP) post-liver transplantation. Patients & methods: A total of 13 SNPs were genotyped in 97 Chinese liver transplant recipients. Associations between SNPs and TAC trough blood concentration/dose ratio (C₀/D) were analyzed using different genetic models in both CP and SP. Results: Only five SNPs were significantly associated with TAC log (C₀/D) in the CP, and none showed a significant association in the SP. We identified rs15524 (CYP3A5), rs9200 (C6), albumin and creatinine as independent predictors of TAC C₀/D in the CP. Furthermore, a final model in the CP explained a total of 30.5% TAC variation. Conclusion: Our study results suggest that in the early stages post-transplantation surgery, recipients' genetic and clinical factors exert a short-term impact on TAC metabolism that gradually decreases with time.

Leveraging drug perturbation to reveal genetic regulators of hepatic gene expression in African Americans

Expression quantitative locus (eQTL) studies have paved the way in identifying genetic variation impacting gene expression levels. African Americans (AAs) are disproportionately underrepresented in eQTL studies, resulting in a lack of power to identify population-specific regulatory variants especially related to drug response. Specific drugs are known to affect the biosynthesis of drug metabolism enzymes as well as other genes. We used drug perturbation in cultured primary hepatocytes derived from AAs to determine the effect of drug treatment on eQTL mapping and to identify the drug response eQTLs (reQTLs) that show altered effect size following drug treatment. Whole-genome genotyping (Illumina MEGA array) and RNA sequencing were performed on 60 primary hepatocyte cultures after treatment with six drugs (Rifampin, Phenytoin, Carbamazepine, Dexamethasone, Phenobarbital, and Omeprazole) and at baseline (no treatment). eQTLs were mapped by treatment and jointly with Meta-Tissue. We found varying transcriptional changes across different drug treatments and identified Nrf2 as a potential general transcriptional regulator. We jointly mapped eQTLs with gene expression data across all drug treatments and baseline, which increased our power to detect eQTLs by 2.7-fold. We also identified 2,988 reQTLs (eQTLs with altered effect size after drug treatment). reQTLs were more likely to overlap transcription factor binding sites, and we uncovered reQTLs for drug metabolizing genes such as CYP3A5. Our results provide insights into the genetic regulation of gene expression in hepatocytes through drug perturbation and provide insight into SNPs that effect the liver's ability to respond to transcription upregulation.

Genotype-guided model significantly improves accuracy of tacrolimus initial dosing after liver transplantation

BACKGROUND

The initial dose of tacrolimus after liver transplantation (LT) is critical for rapidly achieving the steady state of the drug concentration, minimizing the potential adverse reactions and warranting long-term patient prognosis. We aimed to develop and validate a genotype-guided model for determining personalized initial dose of tacrolimus.

METHODS

By combining pharmacokinetic modeling, pharmacogenomic analysis and multiple statistical methods, we developed a genotype-guided model to predict individualized tacrolimus initial dose after LT in the discovery (n = 150) and validation cohorts (n = 97) respectively. This model was further validated in a prospective, randomized and single-blind clinical trial from August, 2021 to February, 2022 (n = 40, ChiCTR2100050288).

FINDINGS

Our model included donor's and recipient's genotypes, recipient's weight and total bilirubin, which achieved an area under the curve of receiver operating characteristic curve (AUC of ROC) of 0.88 and 0.79 in the discovery and validation cohorts, respectively. We found that patients who were given tacrolimus within the recommended concentration range (RCR) (4-10 ng/mL), the new-onset metabolic syndromes are lower, especially for new-onset diabetes (p = 0.043). In the clinical trial, compared to those in experience-based (EB) group, patients in the model-based (MB) group were more likely to achieving the RCR (75% vs 40%, p = 0.025) with a more variable individualized dose (0.023-0.096 mg/kg/day vs 0.045-0.057 mg/kg/day). Moreover, significantly fewer medication adjustments were required for the MB group than the EB group (2.75 ± 2.01 vs 6.05 ± 3.35, p = 0.001).

INTERPRETATION

Our genotype-based model significantly improved the initial dosing accuracy of tacrolimus and reduced the number of medication adjustments, which are critical for improving the prognosis of LT patients.

FUNDING

National Natural Science Foundation of China, Shanghai three-year action plan, National Science and Technology Major Project of China.

Genomewide association study identifies a novel variant associated with tacrolimus trough concentration in Chinese renal transplant recipients

Tacrolimus (TAC) is an immunosuppressant widely used in kidney transplantation. TAC displays considerable interindividual variability in pharmacokinetics (PKs). Genetic and clinical factors play important roles in TAC PKs. We enrolled a total of 251 Chinese renal transplant recipients and conducted a genomewide association study (GWAS), linkage disequilibrium (LD), and one-way analysis of variance (ANOVA) to find genetic variants affecting log-transformed TAC trough blood concentration/dose ratio (log[C₀ /D]). In addition, we performed dual luciferase reporter gene assays and multivariate regression models to evaluate the effect of the genetic variants. The GWAS results showed that all 23 genomewide significant single-nucleotide polymorphisms (p < 5 × 10^-8 ) were located on chromosome 7, including CYP3A5*3. LD, conditional association analysis, and one-way ANOVA showed that rs75125371 T > C independently influenced TAC log(C₀ /D). Dual luciferase reporter gene assays indicated that rs75125371 minor allele (C) was significantly associated with increased normalized luciferase activity than the major allele (T) in the Huh7 cells (p = 1.2 × 10^-5 ) and HepaRG cells (p = 0.0097). A model inclusive of age, sex, hematocrit, CYP3A5*3, and rs75125371 explained 37.34% variance in TAC C₀ . These results suggest that rs75125371 T > C is a functional and population-specific variant affecting TAC C₀ in Chinese renal transplant recipients.

The Impact of Donor and Recipient Genetic Variation on Outcomes After Solid Organ Transplantation: A Scoping Review and Future Perspectives

At the outset of solid organ transplantation, genetic variation between donors and recipients was recognized as a major player in mechanisms such as allograft tolerance and rejection. Genome-wide association studies have been very successful in identifying novel variant-trait associations, but have been difficult to perform in the field of solid organ transplantation due to complex covariates, era effects, and poor statistical power for detecting donor-recipient interactions. To overcome a lack of statistical power, consortia such as the International Genetics and Translational Research in Transplantation Network have been established. Studies have focused on the consequences of genetic dissimilarities between donors and recipients and have reported associations between polymorphisms in candidate genes or their regulatory regions with transplantation outcomes. However, knowledge on the exact influence of genetic variation is limited due to a lack of comprehensive characterization and harmonization of recipients' or donors' phenotypes and validation using an experimental approach. Causal research in genetics has evolved from agnostic discovery in genome-wide association studies to functional annotation and clarification of underlying molecular mechanisms in translational studies. In this overview, we summarize how the recent advances and progresses in the field of genetics and genomics have improved the understanding of outcomes after solid organ transplantation.

Trans-eQTLs of the CYP3A4 and CYP3A5 associated with tacrolimus trough blood concentration in Chinese renal transplant patients

This study aimed to systematically investigate trans-eQTLs of CYP3A4 and CYP3A5 affecting tacrolimus trough blood concentrations in Chinese renal transplant patients. We used Plink v1.90 to perform data quality control and linear regression analysis on GTEx v8 data. SNPs with p-value < 0.05 were selected and the GTEx eQTL Calculator was used to further prioritize the eQTLs of CYP3A4 and CYP3A5 in the liver and small intestine. The eQTLs with a p-value < 5 × 10^-5 and MAF≥ 0.05 in the CHB population were selected as candidate eQTLs. The genotyping of candidate eQTLs was performed using high-resolution melting (HRM) assays and Sanger DNA sequencing. This study included 845 Chinese renal transplant patients who received tacrolimus as an immunosuppressive agent. Association between 103 candidate eQTLs and log-transformed tacrolimus concentration/dose ratio (log (C₀/D)) in this cohort was conducted using the SNPassoc package of R software. In the end, a total of 75,632 liver eQTLs of CYP3A4, 69,558 liver eQTLs of CYP3A5, 48,596 small intestine eQTLs of CYP3A4 and 28,616 small intestine eQTLs of CYP3A5 were obtained using the GTEx v8 eQTL Calculator. Of the 103 candidate eQTLs, rs75727207, rs181294422 and rs28522676 were significantly associated with tacrolimus log(C₀/D) in different genetic models. We discovered a substantial number of novel eQTLs of CYP3A4 and CYP3A5 in liver and small intestine, also found that rs75727207, rs181294422 and rs28522676 may affect tacrolimus trough blood concentrations in Chinese renal transplant patients.

The influence of recipient SLCO1B1 rs2291075 polymorphism on tacrolimus dose-corrected trough concentration in the early period after liver transplantation

Purpose

To explore the relationship between rs2291075 polymorphism in SLCO1B1 gene, which encodes an influx transmembrane protein transporter, and tacrolimus dose–corrected trough concentration (C/D, ng ml⁻¹ mg⁻¹ kg⁻¹) in the early period after liver transplantation.

Methods

CYP3A5 rs776746 and SLCO1B1 rs2291075 polymorphisms of 210 liver transplantation patients and their corresponding donor livers were assessed by PCR amplification and DNA sequencing. The influence of gene polymorphisms on C/D values of tacrolimus was analyzed. The early postoperative period after liver transplantation was divided into the convalescence phase (1–14 days) and stationary phase (15–28 days) according to the change of liver function and tacrolimus C/D values.

Results

The combined analysis of donor and recipient CYP3A5 rs776746 could distinguish the metabolic phenotype of tacrolimus into three groups: fast elimination (FE), intermediate elimination (IE), and slow elimination (SE), which was entitled the FIS classification system. Tacrolimus C/D ratios of recipient SLCO1B1 rs2291075 CT and TT carriers were very close and were significantly lower than those of recipient SLCO1B1 rs2291075 CC genotype carriers in convalescence phase (p = 0.0195) and in stationary phase (p = 0.0152). There were no statistically significant differences between tacrolimus C/D ratios of patients carried with SLCO1B1 rs2291075 CT, TT genotype donors, and those carried with SLCO1B1 rs2291075 CC genotype donors. A model consisting of tacrolimus daily dose, total bilirubin, FIS classification, and recipient SLCO1B1 rs2291075 could predict tacrolimus C/D ratios in the convalescence phase by multivariate analysis. However, recipient SLCO1B1 rs2291075 genotype failed to enter forecast model for C/D ratios in stationary phase. Recipient SLCO1B1 rs2291075 genotype had significant effect on tacrolimus C/D ratios in convalescence phase (p = 0.0300) and stationary phase (p = 0.0400) in subgroup, which excluded the interference come from donor and recipient CYP3A5 rs776746.

Conclusion

SLCO1B1 rs2291075 could be a novel genetic locus associated with tacrolimus metabolism. The combined analysis of donor and recipient CYP3A5 rs776746, recipient SLCO1B1 rs2291075 genotypes, could be helpful to guide the personalized administration of tacrolimus in early period after liver transplantation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00228-020-03058-w.

Importance of genetic polymorphisms in liver transplantation outcomes

Although, liver transplantation serves as the only curative treatment for patients with end-stage liver diseases, it is burdened with complications, which affect survival rates. In addition to clinical risk factors, contribution of recipient and donor genetic prognostic markers has been extensively studied in order to reduce the burden and improve the outcomes. Determination of single nucleotide polymorphisms (SNPs) is one of the most important tools in development of personalized transplant approach. To provide a better insight in recent developments, we review the studies published in the last three years that investigated an association of recipient or donor SNPs with most common issues in liver transplantation: Acute cellular rejection, development of new-onset diabetes mellitus and non-alcoholic fatty liver disease, hepatocellular carcinoma recurrence, and tacrolimus concentration variability. Reviewed studies confirmed previously established SNP prognostic factors, such as PNPLA3 rs738409 for non-alcoholic fatty liver disease development, or the role of CYP3A5 rs776746 in tacrolimus concentration variability. They also identified several novel SNPs, with a reasonably strong association, which have the potential to become useful predictors of post-transplant complications. However, as the studies were typically conducted in one center on relatively low-to-moderate number of patients, verification of the results in other centers is warranted to resolve these limitations. Furthermore, of 29 reviewed studies, 28 used gene candidate approach and only one implemented a genome wide association approach. Genome wide association multicentric studies are needed to facilitate the development of personalized transplant medicine.