Background:Invasive fungal infection (IFI) is one of the leading causes of early death after renal transplantation. Voriconazole (VRC) is the first-line drug of IFI. Because of the large inter- and intraindividual variability in VRC plasma concentrations and the narrow therapeutic window for treating patients with IFIs, it is crucial to study the factors which could influence pharmacokinetic variability. We performed a population pharmacokinetics (PPK) study of VRC for personalized medicine.Methods:A total of 125 trough concentrations (Cmin) from 56 patients were evaluated, retrospectively. Nonlinear mixed effect model was used to describe a PPK model that was internally validated by bootstrap method. Potential covariates included demographic characteristics, physiological and pathological data, concomitant medications, and CYP2C19 genotype.Results:A 1-compartment model with first-order absorption and elimination was fit to characterize the VRC pharmacokinetics in renal transplant recipients (RTRs). Aspartate aminotransferase (AST) had a significant influence on clearance (CL) while CYP2C19 genotype had a major impact on the volume of distribution (V). The parameters of CL and V were 4.76 L/h and 22.47 L, respectively. The final model was V (L) = 22.47 × [1 + 2.21 × (EM = 1)] × [1 + 4.67 × (IM = 1)] × [1 + 3.30 × (PM = 1)] × exp (0.96); CL (L/h) = 4.76 × (AST/33)^(−0.23) × exp (0.14). VRC Cmin in intermediate metabolizers was significantly higher than in extensive metabolizers.Conclusions:Liver function and CYP2C19 polymorphism are major determinants of VRC pharmacokinetic variability in RTRs. Genotypes and clinical biomarkers can determine the initial scheme. Subsequently, therapeutic drug monitoring can optimize clinical efficacy and minimize toxicity. Hence, this is a feasible way to facilitate personalized medicine in RTRs. In addition, it is the first report about PPK of VRC in RTRs.