In this study, we aimed to improve upon a published population pharmacokinetic (PK) model for venlafaxine (VEN) in the treatment of depression in older adults, then investigate whether CYP2D6 metabolizer status affected modelâestimated PK parameters of VEN and its active metabolite Oâdesmethylvenlafaxine. The model included 325 participants from a clinical trial in which older adults with depression were treated with openâlabel VEN (maximum 300âmg/day) for 12âweeks and plasma levels of VEN and Oâdesmethylvenlafaxine were assessed at weeks 4 and 12. We fitted a nonlinear mixedâeffect PK model using NONMEM to estimate PK parameters for VEN and Oâdesmethylvenlafaxine adjusted for CYP2D6 metabolizer status and age. At both lower doses (up to 150âmg/day) and higher doses (up to 300âmg/day), CYP2D6 metabolizers impacted PK modelâestimated VEN clearance, VEN exposure, and active moiety (VEN + Oâdesmethylvenlafaxine) exposure. Specifically, compared with CYP2D6 normal metabolizers, (i) CYP2D6 ultraârapid metabolizers had higher VEN clearance; (ii) CYP2D6 intermediate metabolizers had lower VEN clearance; (iii) CYP2D6 poor metabolizers had lower VEN clearance, higher VEN exposure, and higher active moiety exposure. Overall, our study showed that including a pharmacogenetic factor in a population PK model could increase model fit, and this improved model demonstrated how CYP2D6 metabolizer status affected VENârelated PK parameters, highlighting the importance of genetic factors in personalized medicine.