Genetic risk score and adiposity interact to influence triglyceride levels in a cohort of Filipino women

Polygenic risk for triglyceride levels in the presence of a high impact rare variant

BACKGROUND

Elevated triglyceride (TG) levels are a heritable and modifiable risk factor for cardiovascular disease and have well-established associations with common genetic variation captured in a polygenic risk score (PRS). In young adulthood, the 22q11.2 microdeletion conveys a 2-fold increased risk for mild-moderate hypertriglyceridemia. This study aimed to assess the role of the TG-PRS in individuals with this elevated baseline risk for mild-moderate hypertriglyceridemia.

METHODS

We studied a deeply phenotyped cohort of adults (n = 157, median age 34 years) with a 22q11.2 microdeletion and available genome sequencing, lipid level, and other clinical data. The association between a previously developed TG-PRS and TG levels was assessed using a multivariable regression model adjusting for effects of sex, BMI, and other covariates. We also constructed receiver operating characteristic (ROC) curves using logistic regression models to assess the ability of TG-PRS and significant clinical variables to predict mild-moderate hypertriglyceridemia status.

RESULTS

The TG-PRS was a significant predictor of TG-levels (p = 1.52E-04), along with male sex and BMI, in a multivariable model (pmodel = 7.26E-05). The effect of TG-PRS appeared to be slightly stronger in individuals with obesity (BMI ≥ 30) (beta = 0.4617) than without (beta = 0.1778), in a model unadjusted for other covariates (p-interaction = 0.045). Among ROC curves constructed, the inclusion of TG-PRS, sex, and BMI as predictor variables produced the greatest area under the curve (0.749) for classifying those with mild-moderate hypertriglyceridemia, achieving an optimal sensitivity and specificity of 0.746 and 0.707, respectively.

CONCLUSIONS

These results demonstrate that in addition to significant effects of sex and BMI, genome-wide common variation captured in a PRS also contributes to the variable expression of the 22q11.2 microdeletion with respect to elevated TG levels.

Polygenic Risk of Hypertriglyceridemia Is Modified by BMI

Background: Genetic risk scores (GRSs) have partially improved the understanding of the etiology of moderate hypertriglyceridemia (HTG), which until recently was mainly assessed by secondary predisposing causes. The main objective of this study was to assess whether this variability is due to the interaction between clinical variables and GRS. Methods: We analyzed 276 patients with suspected polygenic HTG. An unweighted GRS was developed with the following variants: c.724C > G (ZPR1 gene), c.56C > G (APOA5 gene), c.1337T > C (GCKR gene), g.19986711A > G (LPL gene), c.107 + 1647T > C (BAZ1B gene) and g.125478730A > T (TRIB gene). Interactions between the GRS and clinical variables (body mass index (BMI), diabetes mellitus, diet, physical activity, alcohol consumption, age and gender) were evaluated. Results: The GRS was associated with triglyceride (TG) concentrations. There was a significant interaction between BMI and GRS, with the intensity of the relationship between the number of alleles and the TG concentration being greater in individuals with a higher BMI. Conclusions: GRS is associated with plasma TG concentrations and is markedly influenced by BMI. This finding could improve the stratification of patients with a high genetic risk for HTG who could benefit from more intensive healthcare interventions.

Gene-environment interactions due to quantile-specific heritability of triglyceride and VLDL concentrations

"Quantile-dependent expressivity" is a dependence of genetic effects on whether the phenotype (e.g., triglycerides) is high or low relative to its distribution in the population. Quantile-specific offspring-parent regression slopes (β_OP) were estimated by quantile regression for 6227 offspring-parent pairs. Quantile-specific heritability (h²), estimated by 2β_OP/(1 + r_spouse), decreased 0.0047 ± 0.0007 (P = 2.9 × 10^-14) for each one-percent decrement in fasting triglyceride concentrations, i.e., h² ± SE were: 0.428 ± 0.059, 0.230 ± 0.030, 0.111 ± 0.015, 0.050 ± 0.016, and 0.033 ± 0.010 at the 90th, 75th, 50th, 25th, and 10th percentiles of the triglyceride distribution, respectively. Consistent with quantile-dependent expressivity, 11 drug studies report smaller genotype differences at lower (post-treatment) than higher (pre-treatment) triglyceride concentrations. This meant genotype-specific triglyceride changes could not move in parallel when triglycerides were decreased pharmacologically, so that subtracting pre-treatment from post-treatment triglyceride levels necessarily created a greater triglyceride decrease for the genotype with a higher pre-treatment value (purported precision-medicine genetic markers). In addition, sixty-five purported gene-environment interactions were found to be potentially attributable to triglyceride's quantile-dependent expressivity, including gene-adiposity (APOA5, APOB, APOE, GCKR, IRS-1, LPL, MTHFR, PCSK9, PNPLA3, PPARγ2), gene-exercise (APOA1, APOA2, LPL), gene-diet (APOA5, APOE, INSIG2, LPL, MYB, NXPH1, PER2, TNFA), gene-alcohol (ALDH2, APOA5, APOC3, CETP, LPL), gene-smoking (APOC3, CYBA, LPL, USF1), gene-pregnancy (LPL), and gene-insulin resistance interactions (APOE, LPL).

Genetics of Hypertriglyceridemia

Hypertriglyceridemia, a commonly encountered phenotype in cardiovascular and metabolic clinics, is surprisingly complex. A range of genetic variants, from single-nucleotide variants to large-scale copy number variants, can lead to either the severe or mild-to-moderate forms of the disease. At the genetic level, severely elevated triglyceride levels resulting from familial chylomicronemia syndrome (FCS) are caused by homozygous or biallelic loss-of-function variants in LPL, APOC2, APOA5, LMF1, and GPIHBP1 genes. In contrast, susceptibility to multifactorial chylomicronemia (MCM), which has an estimated prevalence of ~1 in 600 and is at least 50-100-times more common than FCS, results from two different types of genetic variants: (1) rare heterozygous variants (minor allele frequency <1%) with variable penetrance in the five causal genes for FCS; and (2) common variants (minor allele frequency >5%) whose individually small phenotypic effects are quantified using a polygenic score. There is indirect evidence of similar complex genetic predisposition in other clinical phenotypes that have a component of hypertriglyceridemia, such as combined hyperlipidemia and dysbetalipoproteinemia. Future considerations include: (1) evaluation of whether the specific type of genetic predisposition to hypertriglyceridemia affects medical decisions or long-term outcomes; and (2) searching for other genetic contributors, including the role of genome-wide polygenic scores, novel genes, non-linear gene-gene or gene-environment interactions, and non-genomic mechanisms including epigenetics and mitochondrial DNA.

The evolution of genetic-based risk scores for lipids and cardiovascular disease

PURPOSE OF REVIEW

With improved next-generation sequencing technology, open-access genetic databases and increased awareness of complex trait genetics, we are entering a new era of risk assessment in which genetic-based risk scores (GRSs) will play a clinical role. We review the concepts underlying polygenic models of disease susceptibility and challenges in clinical implementation.

RECENT FINDINGS

Polygenic risk scores are currently used in genetic research on dyslipidemias and cardiovascular disease (CVD). Although the underlying principles for constructing polygenic scores for lipids are established, the lack of consensus on which score to use is indicated by the large number - about 50 - that have been published. Recently, large-scale polygenic scores for CVD appear to afford superior risk prediction compared to small-scale scores. Despite the potential benefits of GRSs, certain biases towards ethnicity and sex need to be worked through.

SUMMARY

We are on the verge of clinical application of GRSs to provide incremental information on dyslipidemia and CVD risk above and beyond traditional clinical variables. Additional work is required to develop a consensus of how such scores will be constructed and measured in a validated manner, as well as clinical indications for their use.

Genetic Risk Score Predictive of the Plasma Triglyceride Response to an Omega-3 Fatty Acid Supplementation in a Mexican Population

Our group built a genetic risk score (GRS) of the plasma triglyceride (TG) response to an omega-3 (n-3) fatty acid (FA) supplementation in Caucasian Canadians that explained 21.53% of the TG variance. The objective was to refine the GRS by fine mapping and to test its association with the TG response in young Mexican adults. A total of 191 participants underwent a 6-week n-3 FA supplementation providing 2.7g/day of docosahexaenoic and eicosapentaenoic acids. Using quantitative polymerase chain reaction (PCR), 103 single-nucleotide polymorphisms (SNPs) were genotyped. A stepwise regression adjusted for age, sex, and body mass index (BMI) was used to select the strongest SNPs to include in the genetic risk model. A GRS was calculated from the sum of at-risk alleles. The contribution of the GRS to the TG response was assessed by ANCOVA with age, sex, and BMI included in the model. Several differences in allele frequency were observed between Canadians and Mexicans. Five lead SNPs were included in the genetic risk model, in which the GRS accounted for 11.01% of the variance of the TG response (p < 0.0001). These findings highlight the important contribution of genetic factors to the heterogeneity of the TG response to an n-3 FA supplementation among Mexicans.

Polygenic influences on dyslipidemias

PURPOSE OF REVIEW

Rare large-effect genetic variants underlie monogenic dyslipidemias, whereas common small-effect genetic variants - single nucleotide polymorphisms (SNPs) - have modest influences on lipid traits. Over the past decade, these small-effect SNPs have been shown to cumulatively exert consistent effects on lipid phenotypes under a polygenic framework, which is the focus of this review.

RECENT FINDINGS

Several groups have reported polygenic risk scores assembled from lipid-associated SNPs, and have applied them to their respective phenotypes. For lipid traits in the normal population distribution, polygenic effects quantified by a score that integrates several common polymorphisms account for about 20-30% of genetic variation. Among individuals at the extremes of the distribution, that is, those with clinical dyslipidemia, the polygenic component includes both rare variants with large effects and common polymorphisms: depending on the trait, 20-50% of susceptibility can be accounted for by this assortment of genetic variants.

SUMMARY

Accounting for polygenic effects increases the numbers of dyslipidemic individuals who can be explained genetically, but a substantial proportion of susceptibility remains unexplained. Whether documenting the polygenic basis of dyslipidemia will affect outcomes in clinical trials or prospective observational studies remains to be determined.

Systolic Blood Pressure Response in SPRINT (Systolic Blood Pressure Intervention Trial) and ACCORD (Action to Control Cardiovascular Risk in Diabetes): A Possible Explanation for Discordant Trial Results

Background

SPRINT (Systolic Blood Pressure Intervention Trial) and the ACCORD (Action to Control Cardiovascular Risk in Diabetes) blood pressure trial used similar interventions but produced discordant results. We investigated whether differences in systolic blood pressure (SBP) response contributed to the discordant trial results.

Methods and Results

We evaluated the distributions of SBP response during the first year for the intensive and standard treatment groups of SPRINT and ACCORD using growth mixture models. We assessed whether significant differences existed between trials in the distributions of SBP achieved at 1 year and the treatment‐independent relationships of achieved SBP with risks of primary outcomes defined in each trial, heart failure, stroke, and all‐cause death. We examined whether visit‐to‐visit variability was associated with heterogeneous treatment effects. Among the included 9027 SPRINT and 4575 ACCORD participants, the difference in mean SBP achieved between treatment groups was 15.7 mm Hg in SPRINT and 14.2 mm Hg in ACCORD, but SPRINT had significantly less between‐group overlap in the achieved SBP (standard deviations of intensive and standard groups, respectively: 6.7 and 5.9 mm Hg in SPRINT versus 8.8 and 8.2 mm Hg in ACCORD; P<0.001). The relationship between achieved SBP and outcomes was consistent across trials except for stroke and all‐cause death. Higher visit‐to‐visit variability was more common in SPRINT but without treatment‐effect heterogeneity.

Conclusions

SPRINT and ACCORD had different degrees of separation in achieved SBP between treatment groups, even as they had similar mean differences. The greater between‐group overlap of achieved SBP may have contributed to the discordant trial results.

Detection of gene-environment interactions in the presence of linkage disequilibrium and noise by using genetic risk scores with internal weights from elastic net regression

BACKGROUND

For the analysis of gene-environment (GxE) interactions commonly single nucleotide polymorphisms (SNPs) are used to characterize genetic susceptibility, an approach that mostly lacks power and has poor reproducibility. One promising approach to overcome this problem might be the use of weighted genetic risk scores (GRS), which are defined as weighted sums of risk alleles of gene variants. The gold-standard is to use external weights from published meta-analyses.

METHODS

In this study, we used internal weights from the marginal genetic effects of the SNPs estimated by a multivariate elastic net regression and thereby provided a method that can be used if there are no external weights available. We conducted a simulation study for the detection of GxE interactions and compared power and type I error of single SNPs analyses with Bonferroni correction and corresponding analysis with unweighted and our weighted GRS approach in scenarios with six risk SNPs and an increasing number of highly correlated (up to 210) and noise SNPs (up to 840).

RESULTS

Applying weighted GRS increased the power enormously in comparison to the common single SNPs approach (e.g. 94.2% vs. 35.4%, respectively, to detect a weak interaction with an OR ≈ 1.04 for six uncorrelated risk SNPs and n = 700 with a well-controlled type I error). Furthermore, weighted GRS outperformed the unweighted GRS, in particular in the presence of SNPs without any effect on the phenotype (e.g. 90.1% vs. 43.9%, respectively, when 20 noise SNPs were added to the six risk SNPs). This outperforming of the weighted GRS was confirmed in a real data application on lung inflammation in the SALIA cohort (n = 402). However, in scenarios with a high number of noise SNPs (>200 vs. 6 risk SNPs), larger sample sizes are needed to avoid an increased type I error, whereas a high number of correlated SNPs can be handled even in small samples (e.g. n = 400).

CONCLUSION

In conclusion, weighted GRS with weights from the marginal genetic effects of the SNPs estimated by a multivariate elastic net regression were shown to be a powerful tool to detect gene-environment interactions in scenarios of high Linkage disequilibrium and noise.

Association between NF-κB Pathway Gene Variants and sICAM1 Levels in Taiwanese

Intercellular adhesion molecule–1 (ICAM1) is crucial to the development and progression of atherosclerosis. Recent genome-wide association studies (GWAS) have revealed that single nucleotide polymorphisms (SNPs) in two of the nuclear factor-κB (NF-κB) pathway genes, NFKBIK and RELA, are associated with soluble ICAM1 (sICAM1) levels. However, neither of these two gene variants is found in the Asian populations. This study aimed to elucidate whether other candidate gene variants involved in the NF-κB pathway may be associated with sICAM1 levels in Taiwanese. After excluding carriers of the ICAM1 rs5491-T allele, three SNPs in the ICAM1 gene and eight SNPs in six of the NF-κB pathway genes (NFKB1, PDCD11, TNFAIP3, NKAPL, IKBKE, and PRKCB) were analyzed for their association with sICAM1 levels in 480 individuals. Our data showed that two SNPs, rs5498 of ICAM1 and rs1635 of NKAPL, were significantly associated with sICAM1 levels (P = 0.002 and 0.004, respectively) in the Taiwanese population. Using a multivariate analysis, rs5498 and rs1635 as well as the previously reported ABO genotypes and rs12051272 of the CDH13 gene were independently associated with sICAM1 levels (P = 0.001, 0.001, 0.006 and 0.031, respectively). An analysis with combined risk alleles of four candidate SNPs in the ICAM1, NKAPL, ABO, and CDH13 genes showed an increase in sICAM1 levels with added numbers of risk alleles and weighted genetic risk score. Our findings thus expanded the repertoire of gene variants responsible for the regulation of sICAM1 levels in the Asian populations.

Increasing insulin resistance accentuates the effect of triglyceride-associated loci on serum triglycerides during 5 years

Blood concentrations of triglycerides are influenced by genetic factors as well as a number of environmental factors, including adiposity and glucose homeostasis. The aim was to investigate the association between a serum triglyceride weighted genetic risk score (wGRS) and changes in fasting serum triglyceride level over 5 years and to test whether the effect of the wGRS was modified by 5 year changes of adiposity, insulin resistance, and lifestyle factors. A total of 3,474 nondiabetic individuals from the Danish Inter99 cohort participated in both the baseline and 5 year follow-up physical examinations and had information on the wGRS comprising 39 genetic variants. In a linear regression model adjusted for age, sex, and baseline serum triglyceride, the wGRS was associated with increased serum triglyceride levels over 5 years [per allele effect = 1.3% (1.0-1.6%); P = 1.0 × 10-17]. This triglyceride-increasing effect of the wGRS interacted with changes in insulin resistance (Pinteraction = 1.5 × 10-6). This interaction indicated that the effect of the wGRS was stronger in individuals who became more insulin resistant over 5 years. In conclusion, our findings suggest that increased genetic risk load is associated with a larger increase in fasting serum triglyceride levels in nondiabetic individuals during 5 years of follow-up. This effect of the wGRS is accentuated by increasing insulin resistance.

Interaction of Insulin Resistance and Related Genetic Variants With Triglyceride-Associated Genetic Variants

BACKGROUND

Several studies suggest that some triglyceride-associated single-nucleotide polymorphisms (SNPs) have pleiotropic and opposite effects on glycemic traits. This potentially implicates them in pathways such as de novo lipogenesis, which is presumably upregulated in the context of insulin resistance. We therefore tested whether the association of triglyceride-associated SNPs with triglyceride levels differs according to one's level of insulin resistance.

METHODS AND RESULTS

In 3 cohort studies (combined n=12 487), we tested the interaction of established triglyceride-associated SNPs (individually and collectively) with several traits related to insulin resistance, on triglyceride levels. We also tested the interaction of triglyceride SNPs with fasting insulin-associated SNPs, individually and collectively, on triglyceride levels. We find significant interactions of a weighted genetic risk score for triglycerides with insulin resistance on triglyceride levels (Pinteraction=2.73×10(-11) and Pinteraction=2.48×10(-11) for fasting insulin and homeostasis model assessment of insulin resistance, respectively). The association of the triglyceride genetic risk score with triglyceride levels is >60% stronger among those in the highest tertile of homeostasis model assessment of insulin resistance compared with those in the lowest tertile. Individual SNPs contributing to this trend include those in/near GCKR, CILP2, and IRS1, whereas PIGV-NROB2 and LRPAP1 display an opposite trend of interaction. In the pooled data set, we also identify a SNP-by-SNP interaction involving a triglyceride-associated SNP, rs4722551 near MIR148A, with a fasting insulin-associated SNP, rs4865796 in ARL15 (Pinteraction=4.1×10(-5)).

CONCLUSIONS

Our findings may thus provide genetic evidence for the upregulation of triglyceride levels in insulin-resistant individuals, in addition to identifying specific genetic loci and a SNP-by-SNP interaction implicated in this process.