Pleiotropic Effects of Common and Rare GCKR Exonic Mutations on Cardiometabolic Traits

Metabolic phenotyping of patients with advanced chronic liver disease for better characterization of cirrhosis regression

BACKGROUND & AIMS

Regression of cirrhosis has been observed in patients with viral and non-viral etiologies of liver disease in whom the underlying cause of liver injury was effectively suppressed. However, the understanding of the factors contributing to reversibility of fibrosis and cirrhosis is limited. Our aims were to assess clinical factors, perform genotyping of known variants, and comprehensive metabolic phenotyping to characterize the regression of fibrosis in patients with compensated advanced chronic liver disease (cACLD).

METHODS

In a case-control pilot study of 81 patients with cACLD, we compared individuals exhibiting histological or clinical evidence of cACLD regression ("regressors"; n = 44) with those showing no improvement ("non-regressors"; n = 37) after a minimum of 24 months of successful treatment of the cause of liver disease. Data were validated using an external validation cohort (n = 30).

RESULTS

Regardless of the cause of cACLD, the presence of obesity (odds ratio [OR] 0.267 95% CI 0.072-0.882; p = 0.049), high liver stiffness (OR 0.960, 95% CI 0.925-0.995; p = 0.032), and carriage of GCKR variant rs1260326 (OR 0.148, 95% CI 0.030-0.773; p = 0.019) are associated with a reduced likelihood of fibrosis regression in a subgroup of 60 patients with ACLD genotyped for known genetic variants. Using liver tissue transcriptomics, we identified metabolic pathways differentiating regressors from non-regressors, with top pathways associated with lipid metabolism - especially fatty acids, bile acids, phospholipids, triacylglycerides (biosynthesis), and the carnitine shuttle. In the entire discovery cohort, we further measured metabolites within the defined pathways, which led to the identification of 33 circulating markers differentiating regressors from non-regressors after etiological therapy. The validation cohort confirmed 14 of the differentially expressed markers.

CONCLUSIONS

We identified and validated a group of lipid biomarkers associated with regression of fibrosis that could be used as non-invasive biomarkers for detecting regression of fibrosis in cACLD.

IMPACT AND IMPLICATIONS

Regression of cirrhosis/advanced chronic liver disease (ACLD) after removal of the underlying cause of liver injury has been observed in human cirrhosis. However, detailed characterization of ACLD regression remains an unmet need. In this study, we provide a comprehensive phenotyping of individuals likely to experience ACLD regression. While obesity, carriage of GCKR variant rs1260326 and high liver stiffness were associated with lower likelihood of regression of ACLD, a signature of circulating lipid metabolites enabled differentiation of regressors from non-regressors after effective etiologic therapy. The lipid signature we discovered and externally validated could be used as non-invasive biomarker to detect regression of fibrosis in patients with compensated ACLD.

Interaction between genetic risk score and dietary fat intake on lipid-related traits in Brazilian young adults

The occurrence of dyslipidaemia, which is an established risk factor for cardiovascular diseases, has been attributed to multiple factors including genetic and environmental factors. We used a genetic risk score (GRS) to assess the interactions between genetic variants and dietary factors on lipid-related traits in a cross-sectional study of 190 Brazilians (mean age: 21 ± 2 years). Dietary intake was assessed by a trained nutritionist using three 24-h dietary recalls. The high GRS was significantly associated with increased concentration of TAG (beta = 0·10 mg/dl, 95 % CI 0·05-0·16; P < 0·001), LDL-cholesterol (beta = 0·07 mg/dl, 95 % CI 0·04, 0·11; P < 0·0001), total cholesterol (beta = 0·05 mg/dl, 95 % CI: 0·03, 0·07; P < 0·0001) and the ratio of TAG to HDL-cholesterol (beta = 0·09 mg/dl, 95 % CI: 0·03, 0·15; P = 0·002). Significant interactions were found between the high GRS and total fat intake on TAG:HDL-cholesterol ratio (Pinteraction = 0·03) and between the high GRS and SFA intake on TAG:HDL-cholesterol ratio (Pinteraction = 0·03). A high intake of total fat (>31·5 % of energy) and SFA (>8·6 % of energy) was associated with higher TAG:HDL-cholesterol ratio in individuals with the high GRS (beta = 0·14, 95 % CI: 0·06, 0·23; P < 0·001 for total fat intake; beta = 0·13, 95 % CI: 0·05, 0·22; P = 0·003 for SFA intake). Our study provides evidence that the genetic risk of high TAG:HDL-cholesterol ratio might be modulated by dietary fat intake in Brazilians, and these individuals might benefit from limiting their intake of total fat and SFA.

Hyperuricemia Subtypes Classified According to Renal Uric Acid Handling Manifesting Distinct Phenotypic and Genetic Profiles in People With Gout

OBJECTIVE

Hyperuricemia can be stratified into four subtypes according to renal uric acid handling. The aim of this study was to comprehensively describe the biologic characteristics (including genetic background) of clinically defined hyperuricemia subtypes in two large geographically independent gout cohorts.

METHODS

Hyperuricemia subtype was defined as renal uric acid overload (ROL), renal uric acid underexcretion (RUE), combined, or renal normal. Twenty single nucleotide polymorphisms (SNPs) previously identified as gout risk loci or associated with serum urate (SU) concentration in the East Asian population were genotyped. Weighted polygenic risk scores were calculated to assess the cumulative effect of genetic risks on the subtypes.

RESULTS

Of the 4,873 participants, 8.8% had an ROL subtype, 60.9% RUE subtype, 23.1% combined subtype, and 7.2% normal subtype. The ROL subtype was independently associated with older age at onset, lower SU, tophi, and diabetes mellitus; RUE was associated with lower body mass index (BMI) and non-diabetes mellitus; the combined subtype was associated with younger age at onset, higher BMI, SU, estimated glomerular filtration rate (eGFR), and smoking; and the normal subtype was independently associated with older age at onset, lower SU, and eGFR. Thirteen SNPs were associated with gout with 6 shared loci and subtype-dependent risk loci patterns. High polygenic risk scores were associated with ROL subtype (odds ratio [OR] = 9.63, 95% confidence interval [95% CI] 4.53-15.12), RUE subtype (OR = 2.18, 95% CI 1.57-3.03), and combined subtype (OR = 6.32, 95% CI 4.22-9.48) compared with low polygenic risk scores.

CONCLUSION

Hyperuricemia subtypes classified according to renal uric acid handling have subtype-specific clinical and genetic features, suggesting subtype-unique pathophysiologic mechanisms.

GlycA and CRP Are Genetically Correlated: Insight into the Genetic Architecture of Inflammageing

Inflammageing is a condition of perpetual low-grade inflammation induced by ageing. Inflammageing may be predicted by the C-reactive protein (CRP) or by a recently described biomarker which measures N-glycosylated side chains of the carbohydrate component of several acute-phase proteins known as GlycA. The objective of this study was to examine in depth the genetic relationships between CRP and GlycA as well as between each of them and other selected cytokines, which may shed light on the mechanisms of inflammageing. Using the Olink 96 Inflammation panel, data on inflammatory mediators for 1518 twins from the TwinsUK dataset were acquired. Summary statistics for genome-wide association studies for several cytokines as well as CRP and GlycA were collected from public sources. Extensive genetic correlation analyses, colocalization and genetic enrichment analyses were carried out to detect the shared genetic architecture between GlycA and CRP. Mendelian randomization was carried out to assess potential causal relationships. GlycA predicted examined cytokines with a magnitude twice as great as that of CRP. GlycA and CRP were significantly genetically correlated (Rg = 0.4397 ± 0.0854, p-value = 2.60 × 10-7). No evidence of a causal relationship between GlycA and CRP, or between these two biomarkers and the cytokines assessed was obtained. However, the aforementioned relationships were explained well by horizontal pleiotropy. Five exonic genetic variants annotated to five genes explain the shared genetic architecture observed between GlycA and CRP: IL6R, GCKR, MLXIPL, SERPINA1, and MAP1A. GlycA and CRP possess a shared genetic architecture, but the relationship between them appears to be modest, which may imply the promotion of differing inflammatory pathways. GlycA appears to be a more robust predictor of cytokines compared to CRP.

Genetic influences on circulating retinol and its relationship to human health

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

Causal links of α-thalassemia indices and cardiometabolic traits and diabetes: MR study

Our study aimed to investigate if genetic variants around 16p13.3's HBA1 locus, associated with erythrocyte indices and HbA1c levels, predict α-thalassemia-related erythrocyte indices, cardiometabolic traits, and diabetes risk in Taiwanese individuals. We analyzed Taiwan Biobank data, including whole-genome sequencing from 1,493 participants and genotyping arrays from 129,542 individuals. First, we performed regional association analysis using whole-genome sequencing data to identify genetic variants significantly associated with erythrocyte indices, confirming their linkage disequilibrium with the α0 thalassemia --SEA deletion mutation, a common cause of α-thalassemia in Southeast Asian populations. Deletion mutation sequencing further validated these variants' association with α-thalassemia. Subsequently, we analyzed genotyping array data, revealing associations between specific genetic variants and cardiometabolic traits, including lipid profiles, HbA1c levels, bilirubin levels, and diabetes risk. Using Mendelian randomization, we established causal relationships between α-thalassemia-related erythrocyte indices and cardiometabolic traits, elucidating their role in diabetes susceptibility. Our findings highlight genetic variants around the α-globin genes as surrogate markers for common α-thalassemia mutations in Taiwan, emphasizing the causal links between α-thalassemia-related erythrocyte indices, cardiometabolic traits, and heightened diabetes risk.

Prevalence of GCKR rs1260326 Variant in Subjects with Obesity Associated NAFLD and T2DM: A Case-Control Study in South Punjab, Pakistan

The glucokinase regulatory protein (GCKR) regulates glycogen metabolism and insulin secretion, and the GCKR rs1260326 is a putative single nucleotide polymorphism (SNP) associated with metabolic disorders including nonalcoholic fatty liver disease (NAFLD) and type 2 diabetes mellitus (T2DM). This study was conducted to investigate the genetic association of the GCKR rs1260326 in NAFLD and T2DM in our population. NAFLD (n = 103), T2DM (n = 100), and control (n = 100) samples were collected and genotyped for GCKR rs1260326 by tetra-arm PCR. The genetic variant GCKR rs1260326 was significantly linked with NAFLD and T2DM, while the GCKR rs1260326 was significantly associated with the progression of obesity only in NAFLD subjects. The frequency of the C allele (mutant) was higher in both NAFLD (f = 0.69) and T2DM (f = 0.66) subjects as compared to healthy controls of NAFLD (0.52) and T2DM (f = 0.32). The frequency of the C allele was also positively linked with the progression of obesity in both diseases. The frequency of the C allele was 0.66, 0.67, and 0.74 in NAFLD normal weight, overweight, and obese subjects, respectively, while the frequency of the C allele was 0.60, 0.60, and 0.74 in T2DM in normal weight, overweight, and obese subjects, respectively. Homozygous mutant (CC) was 53% in both NAFLD and T2DM subjects, while heterozygous mutant (CT) was 15.53% in NAFLD and 22% in T2DM subjects. Wild-type allele (TT) was 31.06% in NAFLD and 25% in T2DM subjects. In conclusion, the GCKR rs1260326 is a highly prevalent SNP in NAFLD and T2DM subjects, which possibly contributed to obesity, insulin resistance, and metabolic disorders in our population.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels

BACKGROUND

Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total.

METHODS

Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes.

RESULTS

Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation.

CONCLUSIONS

The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

Epidemiological-molecular profile of variants associated with type 2 diabetes mellitus in indigenous populations from the Brazilian Amazon

AIMS

While lifestyle factors are strongly associated with Type 2 diabetes (T2DM), genetic characteristics also play a role. However, much of the research on T2DM genetics focuses on European and Asian populations, leaving underrepresented groups, such as indigenous populations with high diabetes prevalence, understudied.

METHODS

We characterized the molecular profile of 10 genes involved in T2DM risk through complete exome sequencing of 64 indigenous individuals belonging to 12 different Amazonian ethnic groups.

RESULTS

The analysis revealed 157 variants, including four exclusive variants in the indigenous population located in the NOTCH2 and WFS1 genes with a modifier or moderate impact on protein effectiveness. Furthermore, a high impact variant in NOTCH2 was also found. Additionally, the frequency of 10 variants in the indigenous group showed significant differences when compared to other global populations that were evaluated.

CONCLUSION

Our study identified 4 novel variants associated with T2DM in the NOTCH2 and WFS1 genes in the Amazonian indigenous populations we studied. In addition, a variant with a high predicted impact in NOTCH2 was also observed. These findings represent a valuable starting point for conducting further association and functional studies, which could help to improve our understanding of the unique characteristics of this population.

Differential Effects of ABCG5/G8 Gene Region Variants on Lipid Profile, Blood Pressure Status, and Gallstone Disease History in Taiwan

ABCG5 and ABCG8 are two key adenosine triphosphate-binding cassette (ABC) proteins that regulate whole-body sterol trafficking. This study aimed to elucidate the association between ABCG5/G8 gene region variants and lipid profile, cardiometabolic traits, and gallstone disease history in Taiwan. A total of 1494 Taiwan Biobank participants with whole-genome sequencing data and 117,679 participants with Axiom Genome-Wide CHB Array data were enrolled for analysis. Using genotype-phenotype and stepwise linear regression analyses, we found independent associations of four Asian-specific ABCG5 variants, rs119480069, rs199984328, rs560839317, and rs748096191, with total, low-density lipoprotein (LDL), and non-high-density lipoprotein (HDL) cholesterol levels (all p ≤ 0.0002). Four other variants, which were in nearly complete linkage disequilibrium, exhibited genome-wide significant associations with gallstone disease history, and the ABCG8 rs11887534 variant showed a trend of superiority for gallstone disease history in a nested logistic regression model (p = 0.074). Through regional association analysis of various other cardiometabolic traits, two variants of the PLEKHH2, approximately 50 kb from the ABCG5/G8 region, exhibited significant associations with blood pressure status (p < 10^-6). In conclusion, differential effects of ABCG5/G8 region variants were noted for lipid profile, blood pressure status, and gallstone disease history in Taiwan. These results indicate the crucial role of individualized assessment of ABCG5/G8 variants for different cardiometabolic phenotypes.

Genome-wide association and Mendelian randomization study of fibroblast growth factor 21 reveals causal associations with hyperlipidemia and possibly NASH

BACKGROUND

Fibroblast growth factor 21 (FGF21) is a hepatokine that produces metabolic benefits, such as improvements of lipid profile. We performed a genome-wide association study (GWAS) to identify genetic variants associated with circulating FGF21 and investigated the causal effects of FGF21 on pertinent outcomes using Mendelian randomization (MR).

METHODS

We conducted a GWAS testing ∼7.8 million DNA sequence variants with circulating FGF21 in a discovery cohort of 6259 Swedish adults with replication in 4483 Swedish women. We then performed two-sample MR analyses of genetically predicted circulating FGF21 in relation to alcohol and nutrient intake, cardiovascular and metabolic biomarkers and diseases, and liver function biomarkers using publicly available GWAS summary statistics data.

RESULTS

Our GWAS identified multiple single-nucleotide polymorphisms with genome-wide significant associations (P < 5 × 10^-8) with circulating FGF21 on chromosomes 2 and 19 in or near the GCKR and FGF21 genes, respectively. The strongest signal at the FGF21 locus (rs2548957, β = 0.181, P < 2.18 × 10^-42) displayed in two-sample MR analyses robust associations with lower alcohol intake, lower circulating low-density lipoprotein cholesterol, apolipoprotein B, C-reactive protein, gamma-glutamyl transferase, and galectin-3 concentrations, and higher circulating insulin-like growth factor-I and alkaline phosphatase concentrations after correcting for multiple testing (P < 0.0018) whereas associations with fat mass, type 2 diabetes, and cardiovascular disease were largely null.

CONCLUSIONS

We identified robust associations of certain genetic variants in or near the GCKR and FGF21 genes with circulating FGF21 concentrations. Furthermore, our results support a strong causal effect of FGF21 on improved lipid profile, reduced alcohol consumption and C-reactive protein concentrations, and liver function biomarkers including fibrosis. We found largely null or weak positive associations with fat mass, diabetes, and cardiovascular disease as well as higher insulin-like growth factor-I concentrations, which could indicate a compensatory increase to regulate the above FGF21 resistant states in humans.

Pleiotropic Effects of APOB Variants on Lipid Profiles, Metabolic Syndrome, and the Risk of Diabetes Mellitus

Apolipoprotein B (ApoB) plays a crucial role in lipid and lipoprotein metabolism. The effects of APOB locus variants on lipid profiles, metabolic syndrome, and the risk of diabetes mellitus (DM) in Asian populations are unclear. We included 1478 Taiwan Biobank participants with whole-genome sequence (WGS) data and 115,088 TWB participants with Axiom genome-wide CHB array data and subjected them to genotype-phenotype analyses using APOB locus variants. Five APOB nonsynonymous mutations, including Asian-specific rs144467873 and rs13306194 variants, were selected from participants with the WGS data. Using a combination of regional association studies, a linkage disequilibrium map, and multivariate analysis, we revealed that the APOB locus variants rs144467873, rs13306194, and rs1367117 were independently associated with total, low-density lipoprotein (LDL), and non-high-density lipoprotein (non-HDL) cholesterol levels; rs1318006 was associated with HDL cholesterol levels; rs13306194 and rs35131127 were associated with serum triglyceride levels; rs144467873, rs13306194, rs56213756, and rs679899 were associated with remnant cholesterol levels; and rs144467873 and rs4665709 were associated with metabolic syndrome. Mendelian randomization (MR) analyses conducted using weighted genetic risk scores from three or two LDL-cholesterol-level-associated APOB variants revealed significant association with prevalent DM (p = 0.0029 and 8.2 × 10^-5, respectively), which became insignificant after adjustment for LDL-C levels. In conclusion, these results indicate that common and rare APOB variants are independently associated with various lipid levels and metabolic syndrome in Taiwanese individuals. MR analyses supported APOB variants associated with the risk of DM through their associations with LDL cholesterol levels.

Pleiotropic genetic architecture and novel loci for C-reactive protein levels

C-reactive protein is involved in a plethora of pathophysiological conditions. Many genetic loci associated with C-reactive protein are annotated to lipid and glucose metabolism genes supporting common biological pathways between inflammation and metabolic traits. To identify novel pleiotropic loci, we perform multi-trait analysis of genome-wide association studies on C-reactive protein levels along with cardiometabolic traits, followed by a series of in silico analyses including colocalization, phenome-wide association studies and Mendelian randomization. We find 41 novel loci and 19 gene sets associated with C-reactive protein with various pleiotropic effects. Additionally, 41 variants colocalize between C-reactive protein and cardiometabolic risk factors and 12 of them display unexpected discordant effects between the shared traits which are translated into discordant associations with clinical outcomes in subsequent phenome-wide association studies. Our findings provide insights into shared mechanisms underlying inflammation and lipid metabolism, representing potential preventive and therapeutic targets.

En masse organoid phenotyping informs metabolic-associated genetic susceptibility to NASH

Genotype-phenotype associations for common diseases are often compounded by pleiotropy and metabolic state. Here, we devised a pooled human organoid-panel of steatohepatitis to investigate the impact of metabolic status on genotype-phenotype association. En masse population-based phenotypic analysis under insulin insensitive conditions predicted key non-alcoholic steatohepatitis (NASH)-genetic factors including the glucokinase regulatory protein (GCKR)-rs1260326:C>T. Analysis of NASH clinical cohorts revealed that GCKR-rs1260326-T allele elevates disease severity only under diabetic state but protects from fibrosis under non-diabetic states. Transcriptomic, metabolomic, and pharmacological analyses indicate significant mitochondrial dysfunction incurred by GCKR-rs1260326, which was not reversed with metformin. Uncoupling oxidative mechanisms mitigated mitochondrial dysfunction and permitted adaptation to increased fatty acid supply while protecting against oxidant stress, forming a basis for future therapeutic approaches for diabetic NASH. Thus, "in-a-dish" genotype-phenotype association strategies disentangle the opposing roles of metabolic-associated gene variant functions and offer a rich mechanistic, diagnostic, and therapeutic inference toolbox toward precision hepatology. VIDEO ABSTRACT.

Genetic Variants at the APOE Locus Predict Cardiometabolic Traits and Metabolic Syndrome: A Taiwan Biobank Study

Several apolipoprotein genes are located at the APOE locus on chromosome 19q13.32. This study explored the genetic determinants of cardiometabolic traits and metabolic syndrome at the APOE locus in a Taiwanese population. A total of 81,387 Taiwan Biobank (TWB) participants were enrolled to undergo genotype−phenotype analysis using data from the Axiom Genome-Wide CHB arrays. Regional association analysis with conditional analysis revealed lead single-nucleotide variations (SNVs) at the APOE locus: APOE rs7412 and rs429358 for total, low-density lipoprotein (LDL), and high-density lipoprotein (HDL) cholesterol levels; CLPTM1 rs3786505 and rs11672748 for LDL and HDL cholesterol levels; and APOC1 rs438811 and APOE-APOC1 rs439401 for serum triglyceride levels. Genotype−phenotype association analysis revealed a significant association of rs429358 and rs438811 with metabolic syndrome and of rs7412, rs438811, and rs439401 with serum albumin levels (p < 0.0015). Stepwise regression analysis indicated that CLPTM1 variants were independently associated with LDL and HDL cholesterol levels (p = 3.10 × 10−15 for rs3786505 and p = 1.48 × 10−15 for rs11672748, respectively). APOE rs429358 and APOC1 rs438811 were also independently associated with metabolic syndrome (p = 2.29 × 10−14) and serum albumin levels (p = 3.80 × 10−6), respectively. In conclusion, in addition to APOE variants, CLPTM1 is a novel candidate locus for LDL and HDL cholesterol levels at the APOE gene region in Taiwan. Our data also indicated that APOE and APOC1 variants were independently associated with metabolic syndrome and serum albumin levels, respectively. These results revealed the crucial role of genetic variants at the APOE locus in predicting cardiometabolic traits and metabolic syndrome.