Exome-Wide Association Study Identifies East Asian-Specific Missense Variant MTHFR C136T Influencing Homocysteine Levels in Chinese Populations RH: ExWAS of tHCY in a Chinese Population

• Published in: Frontiers in genetics 2021-10, Vol.12, p.717621-717621
• Main Authors: Liu, Tianzi, Momin, Mohetaboer, Zhou, Huiyue, Zheng, Qiwen, Fan, Fangfang, Jia, Jia, Liu, Mengyuan, Bao, Minghui, Li, Jianping, Huo, Yong, Liu, Jialin, Zhang, Yaning, Mao, Xuemei, Han, Xiao, Hu, Zhiyuan, Zeng, Changqing, Liu, Fan, Zhang, Yan
• Format: Article
• Language: English
• Subjects: coding variants; Genetics; homocysteine; hyperhomocysteinemia; MTHFR C136T; MTHFR C677T; population heterogeneity
• ISSN: 1664-8021; 1664-8021
• DOI: 10.3389/fgene.2021.717621

Plasma total homocysteine (tHCY) is a known risk factor of a wide range of complex diseases. No genome scans for tHCY have been conducted in East Asian populations. Here, we conducted an exome-wide association study (ExWAS) for tHCY in 5,175 individuals of Chinese Han origin, followed by a replication study in 668 Chinese individuals. The ExWAS identified two loci, 1p36.22 (lead single-nucleotide polymorphism (SNP) rs1801133, MTHFR C677T) and 16q24.3 (rs1126464, DPEP1 ), showing exome-wide significant association with tHCY (p < 5E−7); and both loci have been previously associated with tHCY in non-East Asian populations. Both SNPs were replicated in the replication study ( p < 0.05). Conditioning on the genotype of C677T and rs1126464, we identified a novel East Asian-specific missense variant rs138189536 (C136T) of MTHFR ( p = 6.53E−10), which was also significant in the replication study ( p = 9.8E−3). The C136T and C677T variants affect tHCY in a compound heterozygote manner, where compound heterozygote and homozygote genotype carriers had on average 43.4% increased tHCY than had other genotypes. The frequency of the homozygote C677T genotype showed an inverse-U-shaped geospatial pattern globally with a pronounced frequency in northern China, which coincided with the high prevalence of hyperhomocysteinemia (HHCY) in northern China. A logistic regression model of HHCY status considering sex, age, and the genotypes of the three identified variants reached an area under the receiver operating characteristic curve (AUC) value of 0.74 in an independent validation cohort. These genetic observations provide new insights into the presence of multiple causal mutations at the MTHFR locus, highlight the role of genetics in HHCY epidemiology among different populations, and provide candidate loci for future functional studies.