aDepartment of Genetics, University of Delhi South Campus bNational Institute of Immunology cDepartment of Rheumatology, All India Institute of Medical Sciences, New Delhi dDepartment of Physics and Astrophysics, University of Delhi, Delhi eJawaharlal Nehru Center for Advanced Scientific Research, Bangalore, India fSanta Fe Institute, Santa Fe, Mexico.
Pharmacogenet Genomics. 2014 Apr;24(4):211-9. doi: 10.1097/FPC.0000000000000036.
Methotrexate (MTX) is the drug of first choice for the treatment of rheumatoid arthritis (RA), but is effective only in around 60% of the patients. Identification of genetic markers to predict response is essential for effective treatment within a critical window period of 6 months after diagnosis, but have been hitherto elusive. In this study, we used genome-wide genotype data to identify the potential risk variants associated with MTX (poor)response in a north Indian RA cohort.
Genome-wide genotyping data for a total of 457 RA patients [297 good (DAS28-3≤3.2) and 160 poor (DAS28-3≥5.1) responders] on MTX monotherapy were tested for association using an additive model. Support vector machine and genome-wide pathway analysis were used to identify additional risk variants and pathways. All risk loci were imputed to fine-map the association signals and identify causal variant(s) of therapeutic/diagnostic relevance.
Seven novel suggestive loci from genome-wide (P≤5×10(-5)) and three from support vector machine analysis were associated with MTX (poor)response. The associations of published candidate genes namely DHFR (P=0.014), FPGS (P=0.035), and TYMS (P=0.005) and purine and nucleotide metabolism pathways were reconfirmed. Imputation, followed by bioinformatic analysis indicated possible interaction between two reversely oriented overlapping genes namely ENOSF1 and TYMS at the post-transcriptional level.
In this first ever genome-wide analysis on MTX treatment response in RA patients, 10 new risk loci were identified. These preliminary findings warrant replication in independent studies. Further, TYMS expression at the post-transcriptional level seems to be probably regulated through an antisense-RNA involving the 6-bp ins/del marker in the overlapping segment at 3'UTR of TYMS-ENOSF1, a finding with impending pharmacogenetic applications.
甲氨蝶呤(MTX)是治疗类风湿关节炎(RA)的首选药物,但仅对约 60%的患者有效。在诊断后 6 个月的关键治疗窗口期内,确定预测反应的遗传标志物对于有效治疗至关重要,但迄今为止尚未确定。在这项研究中,我们使用全基因组基因型数据来鉴定与北印度 RA 队列中 MTX(差)反应相关的潜在风险变异。
对总共 457 名接受 MTX 单药治疗的 RA 患者[297 名(DAS28-3≤3.2)和 160 名(DAS28-3≥5.1)差反应者]的全基因组基因分型数据进行了关联测试,使用加性模型。支持向量机和全基因组途径分析用于识别额外的风险变异和途径。所有风险位点均被推断以精细映射关联信号并确定治疗/诊断相关的因果变异。
从全基因组(P≤5×10(-5))和支持向量机分析中发现了七个新的提示性位点与 MTX(差)反应相关。已发表的候选基因 DHFR(P=0.014)、FPGS(P=0.035)和 TYMS(P=0.005)以及嘌呤和核苷酸代谢途径的关联得到了重新确认。推断后,生物信息学分析表明,在转录后水平上,两个反向重叠基因 ENOSF1 和 TYMS 之间可能存在相互作用。
在这项 RA 患者 MTX 治疗反应的全基因组首次分析中,鉴定出了 10 个新的风险位点。这些初步发现需要在独立研究中进行复制。此外,TYMS 在转录后水平的表达似乎可能通过涉及 TYMS-ENOSF1 3'UTR 中 6-bpins/del 标记的反义 RNA 进行调节,这一发现具有即将出现的药物遗传学应用。