Tangtanatakul Pattarin, Thumarat Chisanu, Satproedprai Nusara, Kunhapan Punna, Chaiyasung Tassamonwan, Klinchanhom Siriwan, Wang Yong-Fei, Wei Wei, Wongshinsri Jeerapat, Chiewchengchol Direkrit, Rodsaward Pongsawat, Ngamjanyaporn Pintip, Suangtamai Thanitta, Mahasirimongkol Surakameth, Pisitkun Prapaporn, Hirankarn Nattiya
Department of Transfusion Sciences and Clinical Microbiology, Faculty of Allied Health Sciences, Chulalongkorn University, Bangkok, Thailand.
Section of Translational Medicine, Faculty of Medicine, Mahidol University, Ramathibodi Hospital, Bangkok, Thailand.
Arthritis Res Ther. 2020 Aug 8;22(1):185. doi: 10.1186/s13075-020-02276-y.
Differences in the expression of variants across ethnic groups in the systemic lupus erythematosus (SLE) patients have been well documented. However, the genetic architecture in the Thai population has not been thoroughly examined. In this study, we carried out genome-wide association study (GWAS) in the Thai population.
Two GWAS cohorts were independently collected and genotyped: discovery dataset (487 SLE cases and 1606 healthy controls) and replication dataset (405 SLE cases and 1590 unrelated disease controls). Data were imputed to the density of the 1000 Genomes Project Phase 3. Association studies were performed based on different genetic models, and pathway enrichment analysis was further examined. In addition, the performance of disease risk estimation for individuals in Thai GWAS was assessed based on the polygenic risk score (PRS) model trained by other Asian populations.
Previous findings on SLE susceptible alleles were well replicated in the two GWAS. The SNPs on HLA class II (rs9270970, A>G, OR = 1.82, p value = 3.61E-26), STAT4 (rs7582694, C>G, OR = 1.57, p value = 8.21E-16), GTF2I (rs73366469, A>G, OR = 1.73, p value = 2.42E-11), and FAM167A-BLK allele (rs13277113, A>G, OR = 0.68, p value = 1.58E-09) were significantly associated with SLE in Thai population. Meta-analysis of the two GWAS identified a novel locus at the FBN2 that was specifically associated with SLE in the Thai population (rs74989671, A>G, OR = 1.54, p value = 1.61E-08). Functional analysis showed that rs74989671 resided in a peak of H3K36me3 derived from CD14+ monocytes and H3K4me1 from T lymphocytes. In addition, we showed that the PRS model trained from the Chinese population could be applied in individuals of Thai ancestry, with the area under the receiver-operator curve (AUC) achieving 0.76 for this predictor.
We demonstrated the genetic architecture of SLE in the Thai population and identified a novel locus associated with SLE. Also, our study suggested a potential use of the PRS model from the Chinese population to estimate the disease risk for individuals of Thai ancestry.
系统性红斑狼疮(SLE)患者中不同种族群体间变异表达的差异已有充分记录。然而,泰国人群的遗传结构尚未得到全面研究。在本研究中,我们对泰国人群进行了全基因组关联研究(GWAS)。
独立收集两个GWAS队列并进行基因分型:发现数据集(487例SLE患者和1606名健康对照)和复制数据集(405例SLE患者和1590名无关疾病对照)。数据被推算至千人基因组计划第3阶段的密度。基于不同遗传模型进行关联研究,并进一步进行通路富集分析。此外,基于其他亚洲人群训练的多基因风险评分(PRS)模型,评估泰国GWAS中个体的疾病风险估计性能。
先前关于SLE易感等位基因的发现在两个GWAS中得到了很好的重复。HLA II类基因上的单核苷酸多态性(SNPs)(rs9270970,A>G,比值比[OR]=1.82,p值=3.61E-26)、信号转导和转录激活因子4(STAT4)(rs7582694,C>G,OR=1.57,p值=8.21E-16)、通用转录因子2I(GTF2I)(rs73366469,A>G,OR=1.73,p值=2.42E-11)以及FAM167A-BLK等位基因(rs13277113,A>G,OR=0.68,p值=1.58E-09)在泰国人群中与SLE显著相关。对两个GWAS的荟萃分析在原纤蛋白2(FBN2)基因座发现了一个新位点,该位点在泰国人群中与SLE特异性相关(rs74989671,A>G,OR=1.54,p值=1.61E-08)。功能分析表明,rs74989671位于来自CD14+单核细胞的H3K36me3峰值和来自T淋巴细胞的H3K4me1峰值处。此外,我们表明,从中国人群训练得到的PRS模型可应用于泰国血统个体,该预测模型的受试者工作特征曲线下面积(AUC)达到0.76。
我们展示了泰国人群中SLE的遗传结构,并鉴定出一个与SLE相关的新位点。此外,我们的研究表明,中国人群的PRS模型有可能用于估计泰国血统个体的疾病风险。