《Clinical Pharmacology & Therapeutics》：Resolving CYP2D6 Structural Complexity with Long-Read Sequencing: Implications for Tamoxifen Precision Dosing in Thai Breast Cancer Patients
编辑推荐：
细胞色素P450 2D6 (CYP2D6) 介导的他莫昔芬生物活化为内昔芬是雌激素受体阳性乳腺癌治疗疗效的关键决定因素。在亚洲人群中，结构复杂的CYP2D6等位基因的高患病率使得常规基因分型复杂化，限制了准确的基因型-表型关联。为弥补这些不足，本研究应用牛津纳
细胞色素P450 2D6 (CYP2D6) 介导的他莫昔芬生物活化为内昔芬是雌激素受体阳性乳腺癌治疗疗效的关键决定因素。在亚洲人群中，结构复杂的CYP2D6等位基因的高患病率使得常规基因分型复杂化，限制了准确的基因型-表型关联。为弥补这些不足，本研究应用牛津纳米孔长读长测序来解析CYP2D6的结构复杂性并评估其对内昔芬代谢的影响。研究队列包括492名泰国乳腺癌患者，并对其中361名接受标准剂量他莫昔芬治疗的个体进行多变量分析，以确定CYP2D6活性评分与血浆内昔芬浓度之间的关联。长读长测序识别出82种不同的二倍型，并成功解析了复杂的结构变异，包括普遍存在的CYP2D6*36 + *10非同源重复。总体而言，46.55%的患者被归类为酶功能降低，主要为中间代谢型 (45.53%)。CYP2D6活性评分是决定内昔芬暴露的最强因素 (P < 0.001)。功能降低的基因型与较低的内昔芬水平显著相关（使用≤5.97 ng/mL的探索性基准），且与母体药物浓度无关。这些发现表明，长读长测序能够精确解析CYP2D6的结构复杂性。通过准确捕获罕见和复杂的CYP2D6等位基因，结果表明标准基因分型可能错误分类某些患者，可能使研究人群中很大一部分面临未被识别的他莫昔芬低暴露风险。
**研究背景与问题**
乳腺癌是全球女性最常见的恶性肿瘤，也是泰国重大的公共卫生负担。约70%的乳腺癌为激素受体阳性，辅助内分泌治疗是降低复发和提高生存率的关键。他莫昔芬作为选择性雌激素受体调节剂，是该治疗方案的核心药物。然而，他莫昔芬是一种前体药物，其治疗疗效依赖于在肝脏中广泛生物活化为活性代谢物，特别是4-羟基-N-去甲基他莫昔芬（内昔芬）。内昔芬的循环浓度与治疗效果密切相关，但其生物活化存在显著的个体间差异。未能达到治疗浓度的个体与治疗失败风险增加相关。他莫昔芬的生物活化主要由高度多态性的细胞色素P450 2D6 (CYP2D6) 酶催化。CYP2D6是人类基因组中遗传最复杂的基因座之一，拥有超过180个等位基因变异，导致酶活性范围广泛。尽管药代动力学研究支持CYP2D6基因型是内昔芬暴露的主要决定因素，但CYP2D6指导的他莫昔芬治疗的临床应用仍存在争议。临床研究结果不一致，越来越多地归因于常规基因分型平台的局限性。传统检测方法通常无法捕获CYP2D6变异的全部范围，导致代谢表型错误分类和基因型-表型关系模糊。这一局限性在亚洲人群中尤为突出，因为该人群中功能降低的等位基因和结构复杂的变异体发生频率很高。此外，准确的表型预测需要单倍型定相以区分变异是否位于同一等位基因（顺式）或不同染色体（反式）。短读长测序和基于PCR的方法通常产生未定相的数据，缺乏解析这些复杂结构或区分CYP2D6与CYP2D7（假基因）的能力。这些技术限制可能掩盖CYP2D6基因型与药物代谢之间的真实关系。
**研究目的与方法概述**
为应对这些挑战，本研究对泰国乳腺癌患者队列应用了靶向长读长纳米孔测序，以全面解析CYP2D6的结构复杂性。通过整合全长基因分型与血浆代谢物的定量分析，该研究旨在克服标准基因分型的局限性，确定高风险复杂等位基因的流行率，并研究准确的结构解析对预测低内昔芬暴露的影响。研究队列包括来自泰国国家癌症研究所的492名乳腺癌患者。主要技术方法包括：1）使用长距离PCR扩增全长CYP2D6基因，并通过牛津纳米孔技术进行长读长测序，以实现直接单倍型定相和结构变异检测；2）根据PharmVar联盟和临床药物遗传学实施联盟指南分配CYP2D6活性评分和代谢表型；3）使用液相色谱-串联质谱法量化他莫昔芬及其主要代谢物（包括内昔芬）的血浆浓度；4）应用多变量线性和逻辑回归模型分析CYP2D6活性评分与内昔芬浓度之间的关联，并评估低暴露风险。
**研究结果**
**1. 人口统计学和临床特征**
患者特征显示，他莫昔芬治疗亚组与非治疗亚组在年龄和身体质量指数方面相似，但非治疗组绝经后女性比例更高。两组患者主要为早期、中等级别肿瘤。他莫昔芬治疗组的雌激素受体和孕激素受体阳性率显著更高。
**2. CYP2D6结构变异分析**
通过长距离PCR在所有492个样本中筛查CYP2D6基因的结构变异。全长基因扩增在491个样本中成功。缺失分析鉴定出74个样本携带CYP2D6*5等位基因。非同源重复构型（涉及CYP2D6*36 + *10）普遍存在，在267个样本中通过特征性6.6 kb扩增子被鉴定。
**3. 纳米孔测序性能指标**
使用MinION和PromethION流动池进行测序。PromethION显示出显著更高的通量和平均每碱基覆盖深度，但两个平台都提供了下游分析和稳健基因分型所需的高深度覆盖。
**4. CYP2D6等位基因多样性和预测表型分布**
长读长测序揭示了CYP2D6基因座的高水平遗传多样性，共鉴定出82种不同的二倍型。最常见的三个等位基因分别是*36 + *10 (25.71%；功能降低)、*1 (22.15%；功能正常) 和*10 (21.54%；功能降低)。二倍型分布涵盖所有预测的活性类别。正常代谢型患者表现出广泛的等位基因组合，超快代谢型罕见。中间代谢型构成了患者的很大一部分 (45.53%)，并表现出显著的等位基因复杂性。弱代谢型罕见。约6.91%患者的表型分配不确定。
**5. 纳米孔测序解析复杂CYP2D6结构**
长读长测序生成了全长扩增子，有效解析了CYP2D6基因座的复杂基因组结构。生成长的可定相读长能够稳健地区分CYP2D6与其同源假基因，并精确识别多拷贝等位基因。单倍型定相对解析中间二倍型中的多等位基因结构至关重要。可视化分析证实了同源和非同源重复构型。除了普遍的*36 + *10重排外，测序还能检测罕见的串联构型。
**6. 新型结构变异的鉴定**
鉴定出一种包含*10核心单核苷酸多态性与*153串联的非同源重复构型，暂定名为*new_2。此外，在一个样本中，一条单倍型上存在非同源重复，而配对单倍型上存在*new_2等位基因。通过桑格测序对罕见二倍型进行了正交验证，显示出良好的一致性。
**7. 新型变异的计算机分析**
长读长测序检测到多个PharmVar中未收录的单核苷酸多态性。计算机分析（包括分子动力学模拟和他莫昔芬分子对接）表明，新发现的错义突变并未显著改变整体蛋白质稳定性或活性位点结构，临床数据显示携带者血浆他莫昔芬或代谢物水平无显著偏差，提示这些变异在表型上呈中性。
**8. CYP2D6预测表型与代谢物的关联**
对392名患者的血浆他莫昔芬及其代谢物浓度进行了量化。最终361名具有明确CYP2D6活性评分和雌激素受体状态的患者被纳入关联分析。结果显示，虽然母体他莫昔芬水平在各代谢表型组间可比，但代谢物谱显示出清晰的表型依赖性梯度。正常代谢型的内昔芬平均浓度为10.8 ng/mL。相比之下，中间代谢型的内昔芬浓度显著降低（平均5.3 ng/mL），4-羟基他莫昔芬浓度降低，同时N-去甲基他莫昔芬水平升高，表明下游转化受损。弱代谢型的内昔芬平均浓度最低（2.4 ng/mL），而超快代谢型最高（14.2 ng/mL）。不确定组的内昔芬浓度也低于正常代谢型。内昔芬/N-去甲基他莫昔芬代谢比值随CYP2D6活性增加而逐步升高，证实了该酶在内昔芬形成中的核心作用。
**9. 基因型特异性变异**
评估不同CYP2D6基因型的代谢物分布揭示了显著的表型内变异性。与参考基因型*1/*1相比，大多数正常代谢型维持了可比的内昔芬水平，但特定组合（如*1/*36 + *10）例外。在中间代谢型组内，携带特定无功能或功能降低等位基因的个体内昔芬浓度显著降低。对于弱代谢型，罕见无功能二倍型的描述性评估显示内昔芬水平极低。
**10. 内昔芬浓度和代谢比值的决定因素**
多变量线性回归分析显示，较低的CYP2D6活性评分类别与显著降低的内昔芬浓度相关，而较高的活性评分则与浓度增加相关。在协变量中，母体他莫昔芬浓度是强烈的正向预测因子，而身体质量指数则呈适度的负相关。该模型解释了结果变异的很大一部分。对内昔芬/N-去甲基他莫昔芬代谢比值的分析得出了高度一致的结果。敏感性分析证实了这些关系的稳健性。
**11. 低内昔芬暴露风险**
多变量逻辑回归分析评估了CYP2D6状态与低内昔芬暴露风险之间的关联。在基于活性评分的模型中，较高的活性评分与较低的低内昔芬风险强烈相关。在使用临床药物遗传学实施联盟衍生表型类别的替代模型中，弱代谢型和中间代谢型与正常或超快代谢型相比，低于阈值（5.97 ng/mL）的几率显著更高。在两种参数化模型中，较高的母体他莫昔芬浓度独立预测较低的低内昔芬暴露风险。当使用基于观察分布最低五分位数的替代定义时，这些关联仍然稳健。
**讨论与结论总结**
本研究是泰国人群中同类研究中规模最大的，全面描述了492名乳腺癌患者的CYP2D6遗传变异和预测代谢表型。结果揭示了具有精准医学关键意义的特定人群药物基因组学特征。共鉴定出超过80种不同的二倍型，包括几种先前未表征的变异。这种高度的遗传多样性凸显了标准基因分型平台的局限性。此外，这些发现强调，在结构复杂性高的人群中，他莫昔芬给药的精准性从根本上受限于基础基因组数据的分辨率。
一个关键发现是功能降低等位基因的高患病率，特别是*10和*36 + *10构型，它们共同构成了该队列中大部分中间代谢型。与东亚和东南亚人群的先前报告一致，完全功能等位基因的频率显著低于欧洲和非洲人群。*36 + *10单倍型由于其非同源结构排列，在药物遗传学分析中提出了独特挑战。常规方法通常缺乏分辨率来区分这种混合串联与单纯的功能降低*10等位基因单拷贝。本研究证明，长读长测序通过实现全长、可定相的读长，有效克服了这些限制。这种高分辨率方法能够精确区分CYP2D6与其高度同源的假基因，并解析短读长方法通常忽略的复杂结构变异。
线性和逻辑回归模型共同证实，CYP2D6功能是他莫昔芬代谢个体间变异的主要决定因素。较高的活性评分强烈预测较高的内昔芬浓度和代谢比值，而功能降低则显著增加低暴露风险。重要的是，用于定义低暴露的5.97 ng/mL阈值应被解释为探索性的、队列依赖的基准，而非临床验证的截断值。尽管这些模型解释了很大比例的变异，表明CYP2D6活性降低从根本上降低了内昔芬暴露，突显出泰国人群中很大一部分可能面临药物暴露降低的风险。
虽然先前关于CYP2D6基因分型临床效用的研究结果存在矛盾，但本研究表明，基因分型分辨率和等位基因覆盖范围的差异可能解释了这些分歧。当标准检测方法遗漏罕见序列变异并错误地默认分配为正常等位基因时，常常会出现这种不一致。这在亚洲人群中尤为关键，因为大量罕见但功能重要的等位基因未被常见检测 panel 覆盖。最终，复杂CYP2D6等位基因分辨率的提高揭示，依赖标准基因分型可能错误分类某些患者，掩盖其真实的代谢表型和暴露风险。
除了解析结构变异，长读长测序的无偏性促进了新型错义和内含子变异的发现。尽管计算机模拟表明这些错义变异保持稳定且未显著改变蛋白质构象结构，但它们的鉴定突显了该平台的发现潜力。相反，约7%样本的表型分配仍不确定，突显了在亚洲人群中扩展罕见等位基因功能注释框架的持续需求。
本研究存在一些局限性。首先，虽然基因型与代谢物水平之间的关联很强，但缺乏纵向临床结果数据，阻碍了在该特定队列内直接评估复发风险或无病生存期。其次，本研究仅关注CYP2D6。虽然他莫昔芬代谢是一个复杂的多基因过程，但CYP2D6是内昔芬形成的主要限速酶，也是当前临床药物遗传学实施联盟临床指南的主要焦点。因此，仅CYP2D6基因型无法捕获内昔芬暴露个体间变异的全部范围，需要未来进行全面的多基因研究。最后，尽管纳米孔测序目前比标准PCR检测需要更高的专业技术，但成本效益和自动化流程的快速改进正使其临床整合日益可行。
**研究结论**
本研究表明，长读长测序能够精确解析CYP2D6的结构复杂性。通过准确捕获罕见和复杂的CYP2D6等位基因，结果表明标准基因分型可能错误分类某些患者，可能使研究人群中很大一部分面临未被识别的他莫昔芬低暴露风险。
