Skip to main content

Replication analysis of genetic association of the NCAN-CILP2 region with plasma lipid levels and non-alcoholic fatty liver disease in Asian and Pacific ethnic groups

Abstract

Background

The Neurocan-cartilage intermediate layer protein 2 (NCAN-CILP2) region forms a tight linkage disequilibrium (LD) block and is associated with plasma lipid levels and non-alcoholic fatty liver disease (NAFLD) in individuals of European descent but not in the Malay and Japanese ethnic groups. Recent genome-wide resequence studies identified a missense single-nucleotide polymorphism (SNP) (rs58542926) of the transmembrane 6 superfamily member 2 (TM6SF2) gene in the NCAN-CILP2 region related to hepatic triglyceride content. This study aims to analyze the influences of SNPs in this region on NAFLD and plasma lipid levels in the Asian and Pacific ethnic groups and to reveal the reasons behind positive and negative genetic associations dependent on ethnicity.

Methods

Samples and characteristic data were collected from 3,013 Japanese, 119 Palauan, 947 Mongolian, 212 Thai and 401 Chinese people. Hepatic sonography data was obtained from the Japanese individuals. Genotyping data of five SNPs, rs58542926, rs735273, rs1009136, rs1858999, and rs16996148, were used to verify the effect on serum lipid levels by multiple linear regression, and the association with NAFLD in the Japanese population was examined by logistic regression analysis.

Results

rs58542926 showed significant association with the plasma triglyceride (TG) level in Japanese (P = 0.0009, effect size = 9.5 (±3.25) mg/dl/allele) and Thai (P = 0.0008, effect size = 31.6 (±11.7) mg/dl/allele) study subjects. In Mongolian individuals, there was a significant association of rs58542926 with total cholesterol level (P = 0.0003, 11.7 (±3.2) mg/dl/allele) but not with TG level. In multiple comparisons in Chinese individuals, rs58542926 was weakly (P = 0.022) associated with TG levels, although the threshold for statistical significance was not reached. In Palauan individuals, there was no significant association with the studied SNPs. rs58542926 also showed significant association with Japanese NAFLD. The minor allele (t) increased NAFLD risk (OR 1.682, 95 % CI 1.289–2.196, p value 0.00013).

Conclusion

This study confirmed the genetic association of missense SNP of TM6SF2, rs58542926, with plasma lipid levels in multiple East Asian ethnic groups and with NAFLD in Japanese individuals.

Background

Dyslipidemia, which refers to abnormally high or low lipid levels in circulating blood, is a major heritable risk factor of cardiovascular and cerebrovascular diseases. Many studies have investigated genetic factors of plasma lipid levels to explore the responsible metabolic pathways. Through genome-wide association studies (GWAS), abundant genetic variants associated with plasma lipid levels have been identified, among which the Neurocan-cartilage intermediate layer protein 2 (NCAN-CILP2) region has high statistical significance and a relatively large effect size per allele on low-density lipoprotein cholesterol (LDL) and triglyceride (TG) levels [1, 2]. The NCAN-CILP2 region spans 300 kb on chromosome 19 and forms a tight linkage disequilibrium (LD) block. Eleven genes and one miRNA are encoded in this region. This locus showed consistent and deep association with serum lipid levels in the subsequent studies for individuals of European and Chinese descent [35]. However, the tag SNP in this region, rs16996148, was not correlated with plasma lipids in Malaysian [6] and Japanese populations [7, 8].

Fig. 1
figure 1

a Pairwise linkage disequilibrium analysis of the five studied SNPs in the NCAN-CILP2 region in European (CEU), Japanese, Chinese, Mongolian, Thai and Palauan individuals. Numerals in the left-lower diagonal show |D’| values and in the right-upper diagonal half show r2 values. b Haplotype frequencies in the five studied populations and referential CEU population. The chart on the right shows the estimated haplotypes in the order of rs58542926, rs735273, rs1009136, rs1858999 and rs16996148. Capital letters indicate major alleles, and small letters indicate minor alleles

In addition to plasma lipid levels, the NCAN-CILP2 region was identified as a non-alcoholic fatty liver disease (NALFD)–associated locus by GWAS in individuals of European descent [9] but not in Japanese individuals [10, 11]. The risk allele for NAFLD in this region decreased plasma lipid levels in individuals of European descent [9]. Recent genome-wide resequence studies identified a missense single-nucleotide polymorphism (SNP) (rs58542926) in transmembrane 6 superfamily member 2 (TM6SF2), which is related to hepatic triglyceride content [12, 13]. Functional assessment of Tm6sf2 in mouse and human hepatocytes revealed that it regulated plasma TG and cholesterol levels in mice and lipid accumulation in hepatocytes [1214].

rs16996148 and rs58542926 are on opposite sides of the LD block. There was no significant difference in LD structure between European and Japanese populations in the HapMap database (http://hapmap.ncbi.nlm.nih.gov/). The present study aims to analyze the influences of SNPs in this region on NAFLD and plasma lipid levels in Asian and Pacific ethnic groups and to reveal the reasons behind positive and negative genetic correlations dependent on ethnicity.

Results

Genetic association of NCAN-CILP2 region with East Asian plasma lipid levels

Genotyping success rates were >97 % for all of the SNPs. There were no deviations from Hardy-Weinberg equilibrium for any of the SNPs in all ethnic groups (P > 0.05). Details of all 5 SNPs and minor allele frequencies (MAFs) for each population are described in Table 1. rs58542926 showed significant associations with plasma lipid levels in multiple ethnic groups, which persisted after applying the Bonferroni correction for multiple testing (5 SNPs and 4 phenotypes). In Japanese individuals, rs58542926 was strongly associated with TG level. The minor allele (t) decreased the TG level by 9.5 (±3.25) mg/dl/allele. rs58542926, rs735273, rs1009136 and rs1858999 were weakly associated with TG level in Chinese individuals, but the associations did not reach the statistical significance threshold. In Mongolian individuals, a significant association of rs58542926 was observed with total cholesterol and LDL level, but not with TG level. The minor allele decreased total cholesterol by 11.7 (±3.2) mg/dl/allele and LDL cholesterol by 8.3 (±2.9) mg/dl/allele. However, the minor allele solely decreased the TG level in Thai individuals by 31.6 (±11.7) mg/dl/allele. In Palauan individuals, there were no significant associations with the studied SNPs. In combined analysis, the association of rs58542926 with total cholesterol and TG levels becomes further significant. These results support the previous data that the studied missense SNP in TM6SF2 is the functional SNP in the NCAN-CILP2 region.

Table 1 Association study of SNPs in NCAN-CILP2 region with serum lipid levels in East Asian ethnic groups

Genetic association of the NCAN-CILP2 region with Japanese NAFLD

Among the studied SNPs, rs58542926 showed the most significant association with Japanese NAFLD. The minor allele (t) increased NAFLD risk (OR 1.682, 95 % CI 1.289–2.196, p value 0.00013) (Table 2). The allele that reduced plasma lipid levels was the risk allele for NAFLD.

Table 2 Association study of SNPs in NCAN-CILP2 region with NAFLD in Japanese individuals

Haplotype frequency and linkage disequilibrium structures

To evaluate the LD structure of NCAN-CILP2 in the Asian and Pacific populations, a pairwise linkage disequilibrium map was drawn for each ethnic group, and the maps were compared with |D’| and r 2 values of the CEU population downloaded from the HapMap database (Fig. 1a). The middle three SNPs, rs735273, rs1009136, and rs1858999, showed relatively high r 2 values in all of the ethnic groups. But, the most telomeric side SNP, rs58542926, showed weak LD with the middle LD block, even though the physical distance between rs58542926 and rs735273 was close, 5.8 kb, and they showed a high |D’| value, indicating rare recombinations between them. Conversely, the most centromeric side SNP, rs16996148, showed a moderate r 2 score with rs58542926, despite the long distance between them, 289 kb. The r 2 score between rs58542926 and rs16996148 in the CEU population was relatively high. The estimated haplotype frequency showed a higher prevalence of the haplotype encoding the NAFLD risk allele in the Palauan individuals as compared to other ethnic groups (Fig. 1b). Haplotypes carrying minor allele (t) of rs16996148 in the CEU population mostly (94.7 %) encoded as minor allele (t) of rs58542926, but in the Japanese, Mongolian, Thai and Palauan populations, substantial haplotypes (21.5–37.1 %) encoded major allele (C) of rs58542926 with (t) of rs16996148.

Discussion

This study replicated the genetic association of a missense SNP in TM6SF2 with plasma lipid levels in multiple ethnic groups, while the other four SNPs did not show significant associations. Lysine at residue 167 in TM6SF2 reduced TG but did not reduce cholesterol in the plasma of Japanese and Thai individuals. In Mongolian individuals, rs58542926 reduced total cholesterol but not TG. Despite the higher minor allele frequency in Palauan individuals, rs58542926 was not associated with lipid levels. A small population size was considered as the main cause of the lack of association, but other ethnic-specific factors could have affected the results (e.g., relatively higher prevalence of type 2 diabetes (20.2 %) and obesity (87.8 %, body mass index (BMI) > 25)) in the Palauan population [15]. Population-dependent preference to associated lipid fraction was also observed in previous studies. Norwegians showed association only with total cholesterol, [12] but European Americans showed over 400-fold smaller p value with TG than LDL cholesterol [13]. The Hispanic group did not show any significant association with lipid levels [13]. These results suggest that the regulation of cholesterol and TG in plasma by TM6SF2 is affected by ethnicity. The variable effect of rs58542926 on plasma lipid profiles may be feasibly explained by differences in dietary fat or neutraceutical intake [16] according to ethnic tradition.

This study also replicated the genetic association of TM6SF2 with NAFLD in Japanese. The lysine allele of rs58542926 increased the risk diagnosed of a bright liver in ultrasonography. Although the exact biochemical function of TM6SF2 protein has not yet been revealed, knocking down of Tm6sf2 expression in mouse liver reduced cholesterol and TG levels in mouse plasma and increased hepatic lipid droplet content [12, 13]. In human hepatocyte cell lines, disturbed expression of TM6SF2 reduced VLDL excretion and increased cytoplasmic lipid droplets [14]. Furthermore, amino acid substitution at residue 167 from glutamine to lysine diminished the protein level in the endoplasmic reticulum membrane of the transfected cells, probably via misfolded protein degradation [13]. Therefore, it was suggested that TM6SF2 protein is involved in lipid transportation through the ER membrane. Our genetic results support the theory that the SNP induced a change of the physiological function of the TM6SF2 gene.

rs16996148 was initially identified as an SNP associated with plasma cholesterol and TG in a population of European descent [2]. However, this finding was not replicated in Japanese individuals [7]. The LD score between rs16996148 and rs58542926 in the CEU population was more rigid than in East Asian populations, due to the low frequencies of Ccact and CTGct haplotypes. This distinct LD structure and haplotype frequencies in East Asian populations may have resulted in the lack of association. The haplotype frequency in our studied ethnic groups suggested that the functional SNP in TM6SF2, rs58542926, emerged in the haplotype encoding minor alleles at the remaining four SNPs. It was intriguing that the minor allele haplotype was highly prevalent in Palauan individuals as compared to other ethnic groups. Although rs58542926 is a risk allele for NAFLD, it has been reported to be beneficial for CAD [12, 17]. The NCAN-CILP2 region was also suspected to be a type 2 diabetes–associated locus [18], although the pathogenic locus has not been identified. One explanation for these observations may be that the NAFLD risk allele was dramatically increased in the ancestors of Palauan individuals by natural selection, which would be related to energy metabolism. The NAFLD risk allele may be associated with a phenotype that stores excessive energy in the liver and might be advantageous under conditions of highly restricted food availability. Further population genetic studies are required to determine if the higher prevalence of the Lys167 allele in the Palauan population is a consequence of natural selection or simple genetic drift.

Conclusion

In conclusion, our study replicated the genetic association of missense SNP of TM6SF2, rs58542926, with plasma lipid levels in East Asian multiple ethnic groups. While associated lipid, cholesterol, or TG levels were dependent on ethnicity, the effect sizes were relatively large. The SNP was also associated with Japanese NAFLD. These results support previous genetic studies and can be used to estimate the physiological function of TM6SF2, which could be a personal biomarker for medical management of dyslipidemia and NAFLD.

Methods

Study populations

The present study design was approved by the ethical committee of Jichi Medical University and Kagawa Nutrition University. All participants provided written informed consent. Samples were collected from 3,013 Japanese, 119 Palauan, 947 Mongolian, 212 Thai and 401 Chinese people. Population characteristics have been previously described [15, 19]. Participant age, sex, alcohol consumption (obtained by self-questionnaire), body mass index (BMI), and levels of plasma triglyceride (TG), total cholesterol, low-density lipoprotein cholesterol (LDL), and high-density lipoprotein cholesterol (HDL-c) (mg/dl) were collected. Hepatic sonography data were simultaneously obtained in Japanese samples for fatty liver evaluation, which was examined by an experienced clinical technologist and reviewed by clinicians.

Genotyping and association analysis

The following five SNPs of the NCAN-CILP2 LD block were chosen considering the LD structure of the Japanese population and the results of recent studies [913]: rs58542926 (most adjacent gene; TM6SF2), rs735273 (TM6SF2), rs1009136 (MAU2), rs1858999 (GATAD2A), and rs16996148 (CILP2). All five SNPs were genotyped by the TaqMan Genotyping Assay Systems (Applied Biosystems, Foster City, CA, USA). Multiple linear regression analysis was employed to calculate the effect of each minor allele on four serum lipid concentration levels: TG, total cholesterol, LDL, and HDL. The analysis was calculated separately in each ethnic group. Age, sex, and BMI were included for adjustment as covariant factors. Triglyceride levels were transformed to logarithm to reduce the overestimated effect from skewed distribution, but untransformed values of triglycerides were also assessed for the purpose of showing an effect size of SNPs. To study the association of these five SNPs and NAFLD in the Japanese population, participants who had excessive alcohol consumption (>20 g/day) history were excluded. Then, genotyping results of 2,281 of the 3,013 Japanese DNA panels were studied. Logistic regression analysis was applied with adjustments for age, sex, BMI, and type 2 diabetes mellitus (DM2). All statistical analyses were performed with SPSS software. Haplotype frequencies for each population were calculated, and LD blocks were drawn by Haploview software [20].

Abbreviations

BMI:

body mass index

DM2:

type 2 diabetes mellitus

GWAS:

genome-wide association studies

HDL-c:

high-density lipoprotein cholesterol

LD:

linkage disequilibrium

LDL:

low-density lipoprotein cholesterol

MAF:

minor allele frequency

NAFLD:

non-alcoholic fatty liver disease

NCAN-CILP2:

Neurocan-cartilage intermediate layer protein 2

SNP:

single-nucleotide polymorphism

TG:

triglyceride

TM6SF2:

transmembrane 6 superfamily member 2

References

  1. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, et al. Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008;40:161–9.

    Article  PubMed  CAS  Google Scholar 

  2. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, et al. Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008;40:189–97.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Waterworth DM, Ricketts SL, Song K, Chen L, Zhao JH, Ripatti S, et al. Genetic variants influencing circulating lipid levels and risk of coronary artery disease. Arterioscler Thromb Vasc Biol. 2010;30(11):2264–76.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Keebler ME, Deo RC, Surti A, Konieczkowski D, Guiducci C, Burtt N, et al. Fine-Mapping in African Americans of Eight Recently Discovered Genetic Loci for Plasma Lipids: The Jackson Heart Study. Circ Cardiovasc Genet. 2010;3(4):358–64.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Yan TT, Yin RX, Li Q, Huang P, Zeng XN, Huang KK, et al. Sex-specific association of rs16996148 SNP in the NCAN/CILP2/PBX4 and serum lipid levels in the Mulao and Han populations. Lipids Health Dis. 2011;10:248.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Tai ES, Sim XL, Ong TH, Wong TY, Saw SM, Aung T, et al. Polymorphisms at newly identified lipid-associated loci are associated with blood lipids and cardiovascular disease in an Asian Malay population. J Lipid Res. 2009;50(3):514–20.

    Article  PubMed  CAS  Google Scholar 

  7. Kamatani Y, Matsuda K, Okada Y, Kubo M, Hosono N, Daigo Y, et al. Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat Genet. 2010;42:210–5.

    Article  PubMed  CAS  Google Scholar 

  8. Nakayama K, Bayasgalan T, Yamanaka K, Kumada M, Gotoh T, Utsumi N, et al. Large scale replication analysis of loci associated with lipid concentrations in a Japanese population. J Med Genet. 2009;46:370–4.

    Article  PubMed  CAS  Google Scholar 

  9. Speliotes EK, Yerges-Armstrong LM, Wu J, Hernaez R, Kim LJ, Palmer CD, et al. Genome-wide association analysis identifies variants associated with nonalcoholic Fatty liver disease that have distinct effects on metabolic traits. PLoS Genet. 2011;7:e1001324.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Kawaguchi T, Sumida Y, Umemura A, Matsuo K, Takahashi M, Takamura T, et al. Genetic Polymorphisms of the Human PNPLA3 Gene Are Strongly Associated with Severity of Non-Alcoholic Fatty Liver Disease in Japanese. PLoS One. 2012;7(6):e38322.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Kitamoto T, Kitamoto A, Yoneda M, Hyogo H, Ochi H, Nakamura T, et al. Genome-wide scan revealed that polymorphisms in the PNPLA3, SAMM50, and PARVB genes are associated with development and progression of nonalcoholic fatty liver disease in Japan. Hum Genet. 2013;132:783–92.

    Article  PubMed  CAS  Google Scholar 

  12. Holmen OL, Zhang H, Fan Y, Hovelson DH, Schmidt EM, Zhou W, et al. Systematic evaluation of coding variation identifies a candidate causal variant in TM6SF2 influencing total cholesterol and myocardial infarction risk. Nat Genet. 2014;46(4):345–51.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Kozlitina J, Smagris E, Stender S, Nordestgaard BG, Zhou HH, Tybjærg-Hansen A, et al. Exome-wide association study identifies a TM6SF2 variant that confers susceptibility to nonalcoholic fatty liver disease. Nat Genet. 2014;46:352–6.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Mahdessian H, Taxiarchis A, Popov S, Silveira A, Franco-Cereceda A, Hamsten A, et al. TM6SF2 is a regulator of liver fat metabolism influencing triglyceride secretion and hepatic lipid droplet content. Proc Natl Acad Sci U S A. 2014;111:8913–8.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  15. Nakayama K, Yanagisawa Y, Ogawa A, Ishizuka Y, Munkhtulga L, Charupoonphol P, et al. High prevalence of an anti-hypertriglyceridemic variant of the MLXIPL gene in Central Asia. J Hum Genet. 2011;56(12):828–33.

    Article  PubMed  CAS  Google Scholar 

  16. Scicchitano P, Cameli M, Maiello M, Modesti PA, Muiesan ML, Novo S, et al. Nutraceuticals and dyslipidaemia: Beyond the common therapeutics. J Funct Food. 2014;6:11–32.

    Article  CAS  Google Scholar 

  17. Zhou L, Ding H, Zhang X, He M, Huang S, Xu Y, et al. Genetic Variants at Newly Identified Lipid Loci Are Associated with Coronary Heart Disease in a Chinese Han Population. PLoS One. 2011;6(11):e27481.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Saxena R, Elbers CC, Guo Y, Peter I, Gaunt TR, Mega JL, et al. Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci. Am J Hum Genet. 2012;90(3):410–25.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Nakayama K, Ogawa A, Miyashita H, Tabara Y, Igase M, Kohara K, et al. Positive natural selection of TRIB2, a novel gene that influences visceral fat accumulation, in East Asia. Hum Genet. 2013;132(2):201–17.

    Article  PubMed  CAS  Google Scholar 

  20. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

This study was supported by grant-in-aid (Nos. 25461377 and 26291096) from the Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT); a MEXT-supported program for the strategic research foundation at private universities (2013-2017); and a grant from the Takeda Science Foundation (No. 3-2313-004), Osaka, Japan. We are grateful to Kayo Nagashima and Yukiko Ohashi for their excellent technical assistance.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sadahiko Iwamoto.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SB: designed and performed experiments, analyzed data, wrote and revised manuscript; KN, YK: collected and analyzed data, wrote and revised manuscript; SM, KW, HM, ML: collected data and prepared experiments; SI: designed study, analyzed data, and critically revised manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Boonvisut, S., Nakayama, K., Makishima, S. et al. Replication analysis of genetic association of the NCAN-CILP2 region with plasma lipid levels and non-alcoholic fatty liver disease in Asian and Pacific ethnic groups. Lipids Health Dis 15, 8 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12944-016-0181-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12944-016-0181-z

Keywords