The Contribution of GWAS Loci in Familial Dyslipidemias

Abstract

Familial combined hyperlipidemia (FCH) is a complex and common familial dyslipidemia characterized by elevated total cholesterol and/or triglyceride levels with over five-fold risk of coronary heart disease. The genetic architecture and contribution of rare Mendelian and common variants to FCH susceptibility is unknown. In 53 Finnish FCH families, we genotyped and imputed nine million variants in 715 family members with DNA available. We studied the enrichment of variants previously implicated with monogenic dyslipidemias and/or lipid levels in the general population by comparing allele frequencies between the FCH families and population samples. We also constructed weighted polygenic scores using 212 lipid-associated SNPs and estimated the relative contributions of Mendelian variants and polygenic scores to the risk of FCH in the families. We identified, across the whole allele frequency spectrum, an enrichment of variants known to elevate, and a deficiency of variants known to lower LDL-C and/or TG levels among both probands and affected FCH individuals. The score based on TG associated SNPs was particularly high among affected individuals compared to non-affected family members. Out of 234 affected FCH individuals across the families, seven (3%) carried Mendelian variants and 83 (35%) showed high accumulation of either known LDL-C or TG elevating variants by having either polygenic score over the 90th percentile in the population. The positive predictive value of high score was much higher for affected FCH individuals than for similar sporadic cases in the population. FCH is highly polygenic, supporting the hypothesis that variants across the whole allele frequency spectrum contribute to this complex familial trait. Polygenic SNP panels improve identification of individuals affected with FCH, but their clinical utility remains to be defined.

Familial combined hyperlipidemia (FCH) is a complex and common familial dyslipidemia characterized by elevated total cholesterol and/or triglyceride levels with over five-fold risk of coronary heart disease. The genetic architecture and contribution of rare Mendelian and common variants to FCH susceptibility is unknown. In 53 Finnish FCH families, we genotyped and imputed nine million variants in 715 family members with DNA available. We studied the enrichment of variants previously implicated with monogenic dyslipidemias and/or lipid levels in the general population by comparing allele frequencies between the FCH families and population samples. We also constructed weighted polygenic scores using 212 lipid-associated SNPs and estimated the relative contributions of Mendelian variants and polygenic scores to the risk of FCH in the families. We identified, across the whole allele frequency spectrum, an enrichment of variants known to elevate, and a deficiency of variants known to lower LDL-C and/or TG levels among both probands and affected FCH individuals. The score based on TG associated SNPs was particularly high among affected individuals compared to non-affected family members. Out of 234 affected FCH individuals across the families, seven (3%) carried Mendelian variants and 83 (35%) showed high accumulation of either known LDL-C or TG elevating variants by having either polygenic score over the 90th percentile in the population. The positive predictive value of high score was much higher for affected FCH individuals than for similar sporadic cases in the population. FCH is highly polygenic, supporting the hypothesis that variants across the whole allele frequency spectrum contribute to this complex familial trait. Polygenic SNP panels improve identification of individuals affected with FCH, but their clinical utility remains to be defined.

Author Summary

Familial combined hyperlipidemia (FCH) is a familial dyslipidemia and the most common familial risk factor for premature coronary heart disease. Its genetic architecture is poorly understood. Rare high-impact variants have been identified in some patients, but have not explained a substantial portion of the trait. FCH has previously been speculated to be a polygenic disorder, but genetic data supporting this hypothesis have so far been incomplete. We provide experimental evidence for the polygenicity and heterogeneity of FCH in a large set of affected families using comprehensive genome-wide variant data. Approximately a third of the affected FCH individuals in our sample had high polygenic burden, and only a minority carried high-impact variants identifiable by genotyping. We show that the polygenic burden of affected FCH family members is comparable to that observed in individuals with similar lipid phenotypes in the general population. Genetic variants identified in large-scale population studies can also underlie the typical phenotypes observed in complex familial diseases such as FCH. Advances in genetic diagnosis based on population samples may thus also benefit FCH families. Families without high polygenic burden are good candidates for sequencing studies to identify rare variants not observable with genotyping.

Citation: Ripatti P, Rämö JT, Söderlund S, Surakka I, Matikainen N, Pirinen M, et al. (2016) The Contribution of GWAS Loci in Familial Dyslipidemias. PLoS Genet 12(5): e1006078. doi:10.1371/journal.pgen.1006078

Editor: Greg Gibson, Georgia Institute of Technology, UNITED STATES

Received: January 22, 2016; Accepted: May 3, 2016; Published: May 26, 2016

Copyright: © 2016 Ripatti et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Due to the consent given by study participants, data cannot be made publicly available. Data are available through the Institute for Molecular Medicine Finland FIMM Data Access Committee (DAC) for authorized researchers who have IRB/ethics approval and an institutionally approved study plan. For more details, please contact the FIMM DAC ([email protected]).

Funding: Genotyping was funded by National Institutes of Health (NIH) (http://www.nih.gov) (grant HL113315). SR, MRT, and AP were supported by NIH. PP was supported by the NIH grants HL095056 and HL28481. SR, VS, MRT and AP were supported by Finnish Foundation for Cardiovascular Research (http://www.sydantutkimussaatio.fi). SR was additionally supported by Sigrid Juselius Foundation (http://www.sigridjuselius.fi) and Biocentrum Helsinki (http://www.helsinki.fi/biocentrum/). MRT was supported by Sigrid Juselius Foundation, EU-project RESOLVE (EU 7th Framework Program, grant 305707) (http://cordis.europa.eu/fp7/home_en.html) and Helsinki University Central Hospital Research Funds (http://www.med.helsinki.fi/english/research/). PR was funded by travel grants from the University of Helsinki Jubilee Fund (http://www.helsinki.fi/insight/rahastot/juhlarahasto.html) and Doctoral Programme in Population Health (http://www.helsinki.fi/docpop/index.html). JTR was supported by the MD/PhD Program of the Faculty of Medicine, University of Helsinki, (http://www.med.helsinki.fi/peruskoulutus/perustutkinnot/tutkijalaakari_ohjeet.html), the Doctoral Programme in Population Health, University of Helsinki, and by grants from the Finnish Medical Foundation (http://www.laaketieteensaatio.fi), Emil Aaltonen Foundation (http://www.emilaaltonen.fi) and Biomedicum Helsinki Foundation (http://www.biomedicum.com/index.php?page=112). We acknowledge support from the Academy of Finland (grant 251217) (http://www.aka.fi). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Familial combined hyperlipidemia (FCH), classically defined by elevations in serum total cholesterol (TC), triglycerides (TG), or both, in two or more first degree relatives, displays a prevalence of greater than 1% in Western populations [1, 2]. It is the most common familial risk factor for premature coronary heart disease (CHD), occurring in 11-14% of individuals with this condition, and raising by up to five-fold the CHD risk in first- and second-degree relatives of affected individuals [3, 4]. In clinical practice, FCH is characterized by elevations of low-density lipoprotein cholesterol (LDL-C), TG, or both [5]. The phenotype within a family shows high inter- and intraindividual variability of lipid values (TG, LDL-C, high-density lipoprotein cholesterol [HDL-C], and apolipoprotein B) and therefore the diagnosis is commonly missed. Despite attempts to identify rare high-impact variants underlying FCH, no such variant has explained a substantial proportion of the trait [6]. Also, so far the genetics of FCH has not been addressed using high-density genotyping panels [7].

Here, we present a comprehensive evaluation of the genetic background of FCH using a dense-marker genotyping panel. We hypothesized that FCH risk could derive in part from a combination of common variants associated, in population cohorts, with LDL-C and/or TG as well as uncommon variants in genes that have been implicated in Mendelian lipid syndromes. To test our hypothesis we evaluated Finnish FCH families to determine if they demonstrated an enrichment of known lipid-associated variants, compared to Finnish general population samples, and assessed whether such variants could account for differences in lipid levels and degree of aggregation of disease observed between the FCH families. Additionally, following the logic of recent studies that have demonstrated the combined impact of multiple small effect variants (polygenic scores) on risk for familial hypercholesterolemia (FH) and other lipid disorders [8-11], we measured, in these families, the relative contributions of LDL-C and TG associated polygenic scores to the lipid traits and the FCH phenotype.

Results

We evaluated the contribution of lipid-level associated genetic variation to FCH in 715 genotyped individuals (234 of whom were considered affected by FCH; Table 1, S1 Table, and S1 Fig) from 53 FCH families (average family size 13 individuals). Of the 234 FCH-affected individuals, 78 (33%) were identified because of elevated (≥ 90th age- and sex-specific population percentile) TC without elevated TG, 76 (32%) had elevated TG without elevated TC, and 80 (34%) had both lipids elevated. FCH probands (n = 48, because five probands had DNA unavailable) and other FCH affected individuals showed higher mean levels of both LDL-C (4.45 and 4.29 mmol/l, respectively) and TG (4.06 and 2.51 mmol/l, respectively) compared to all FCH family members (LDL-C and TG 3.64 and 1.60 mmol/l, respectively) and to the general Finnish population (LDL-C and TG 3.43 and 1.43 mmol/l, respectively) (Fig 1). The median LDL-C and TG values of the genotyped family members ranged from 1.08 to 4.69 mmol/l for LDL-C and from 0.70 to 4.44 mmol/l for TG demonstrating considerable variation among families (S2 Fig).

Fig 1. Distributions of lipid levels in subsets of the Finnish general population and the FCH samples.

Distributions of (a) LDL-C and (b) TG are shown for the Finnish FINRISK population cohort (FINRISK all, blue), hyperlipidemic Finnish population samples (FINRISK hyperlipidemic, green), all FCH family members (FCH all, brown), affected family members (FCH affected, red), and proband individuals (FCH probands, purple). Hyperlipidemia in the population samples is defined as TC or TG ≥ 90th age- and sex-specific population percentile, analogously with the FCH diagnostic criteria. In (b) the x-axis is cut at 9 mmol/l. FCH, familial combined hyperlipidemia; FINRISK, The National FINRISK Study.

http://dx.doi.org/10.1371/journal.pgen.1006078.g001

Table 1. Clinical and metabolic characteristics of genotyped individuals.

http://dx.doi.org/10.1371/journal.pgen.1006078.t001

Compared to 18,715 random Finnish population samples, alleles that elevated LDL-C and TG levels in the population were overall enriched and alleles that lowered LDL-C and TG levels were depleted in 234 affected FCH individuals (sign test p = 0.0025 for LDL-C and p = 0.0016 for TG). We included in this analysis all SNPs with at least one affected carrier (194 of the 212 SNPs previously associated with either LDL-C or TG in population genome-wide association (GWA) studies or with monogenic dyslipidemias). In total, 60 out of 95 LDL-C elevating SNPs were more common in EUFAM than in population samples, and 57 out of 99 of LDL-C lowering SNPs were less common. Similarly, 61 out of 96 TG elevating SNPs were enriched and 57 out of 98 TG lowering SNPs were depleted (Fig 2, S3 Table). S3 Fig shows the enrichment ratios for probands and all family members compared to the population samples.

Fig 2. Enrichment of LDL-C or TG associated SNPs in FCH affected individuals by their frequency and effect on LDL-C and TG levels.

Enrichment ratio is the ratio of the allele frequencies in the affected individuals (n = 234) to the allele frequencies in the Finnish FINRISK population cohort (n = 18,715). Only individuals without diabetes or other relevant confounders were included (S1 Text). Under the null hypothesis of no enrichment, a 95% credible interval (shaded area) was estimated by calculating the enrichment statistic (enrichment ratio) for all variants with MAF > 0.001% across the genome excluding the loci of the 212 SNPs. Variants are designated as either lipid level elevating (red) or lowering (blue) for (a) LDL-C and (b) TG based on [Beta] estimates from linear regression in the FINRISK samples. Point size and color intensity reflect the magnitude of the effect. Only SNPs with at least one heterozygous carrier are shown (n = 194). *The enrichment ratio for LPL rs1801177 fell within the 95% credible interval.

http://dx.doi.org/10.1371/journal.pgen.1006078.g002

We then examined if high-impact variants implicated in monogenic dyslipidemias would be enriched and contribute to FCH. APOA5 rs3135506, which predisposes to hypertriglyceridemia in homozygous form [12], was 1.8-fold more frequent in the affected individuals of the FCH families than in the general population (minor allele frequency [MAF] in affected FCH individuals = 0.11, MAF in population = 0.062). In all FCH family members, there were 105 heterozygotes (43 of them affected) and six homozygotes (five of them affected). Homozygotes had a mean TG level (mmol/l) of 2.61, heterozygotes 1.97 and wild type carriers 1.55 (S4 Fig). APOE variations have a well-established role in dyslipidemias yet incompletely defined contribution to FCH. By genotyping rs7412 and imputing rs429358 we were able to determine the phenotyped apolipoprotein E isoform with a minimum accuracy of 95% (Table 2). The APOE [varepsilon]2[varepsilon]2 haplotype that predisposes to type III hyperlipoproteinemia [13], was observed in three FCH family members (two of them affected), all with elevated cholesterol and TG concentrations in very low-density lipoprotein and intermediate-density lipoprotein fractions (Table 2). Its frequency was 0.0023 in the Finnish population, 0.0045 in all FCH families and 0.0085 among FCH affected individuals. LIPC rs28933094 predisposes to hepatic lipase deficiency, the cardiovascular effects of which are unclear but may be dependent on the underlying lipid phenotype [14, 15]. It is 4.8-fold more frequent in Finns compared to other Europeans and additionally 2.6-fold more frequent in affected FCH individuals compared to the Finnish population (MAF in affected FCH individuals = 0.041, MAF in Finnish population = 0.016, MAF in Non-Finnish Europeans = 0.0033). Finally, it is of special interest whether the FCH families carry any of the classical FH variants. In the Finnish population there are five major extensively documented, FH-associated LDLR variants [16]. As the genotyping array does not capture these, we successfully imputed one out of the five variants (FH-Pogosta), but observed no carriers. In summary, variants that are strongly linked to known forms of dyslipidemia (in APOA5 and APOE), could at most explain only 7 (3.0%) of all 234 affected FCH individuals, thus being of minor importance in our FCH family sample.

Table 2. Apolipoprotein E phenotypes, imputed haplotypes, and the carriers' lipid profiles.

http://dx.doi.org/10.1371/journal.pgen.1006078.t002

We then estimated what proportion of the affected FCH individuals had a high cumulative burden of enriched LDL-C and TG elevating variants or carried high-impact variants. Out of the 234 affected FCH individuals, 83 individuals (35%) showed high accumulation of known LDL-C or TG elevating variants, having either one or both polygenic lipid scores over the 90th percentile in the population (Table 3). Three additional FCH affected individuals carried Mendelian variants and did not have elevated polygenic scores. In 14 out of the 53 (25%) families, over half of the affected individuals had high polygenic scores or carried known Mendelian variants (Fig 3). In six out of the 53 (11%) families, all affected members had high polygenic scores or carried known Mendelian variants. This did not appear to be driven by differences in genetic correlation between affected individuals in the families (S5 Fig).

Fig 3. Number of affected and unaffected individuals with high polygenic lipid scores or carriers of high-impact Mendelian variants.

The black shading presents the number of (a) affected, and (b) unaffected individuals with high polygenic lipid scores (LDL-C, TG, or both polygenic scores over the 90th percentile in the population) or carriers of a high-impact Mendelian variant (APOE [varepsilon]2[varepsilon]2, n = 2; or homozygosity for APOA5 rs3135506). The families are sorted by the number of affected individuals with high polygenic lipid scores or a high-impact Mendelian variant in (a). The grey shading presents the number of other (a) affected, and (b) unaffected individuals in the family.

http://dx.doi.org/10.1371/journal.pgen.1006078.g003

Table 3. Number of affected individuals with high polygenic lipid score or Mendelian SNPs.

http://dx.doi.org/10.1371/journal.pgen.1006078.t003

Finally, we evaluated how well the high polygenic scores predicted FCH affectedness in the families or among hyperlipidemic individuals in the population (Table 4). Having a high polygenic score (≥ the 90th population percentile) for either LDL-C or TG had a positive predictive value (PPV) of 0.45 in the FCH families and 0.23 in the general population. Negative predictive values (NPV) were 0.71 and 0.89, respectively.

Table 4. Positive and negative predictive values of high polygenic lipid scores.

http://dx.doi.org/10.1371/journal.pgen.1006078.t004

In summary, we observed enrichment of both common and uncommon lipid-elevating variants in the FCH families. Only a minority (7 out of 234, 3%) of affected individuals were carriers of high-impact variants of Mendelian dyslipidemias. The polygenic lipid scores contribute to FCH in over one third of 234 FCH subjects in these families, with considerable heterogeneity between families. In more than half of the families we did not observe either an increased load of common variants or carriers of high-impact variants.

Discussion

Our results demonstrate an enrichment of many known lipid-level elevating variants in FCH families. This enrichment was observed for both uncommon and common variants known to affect either LDL-C or TG levels in populations. When the known variants were combined into polygenic lipid scores for each individual, 17% and 25% of the affected individuals had polygenic scores that were higher than the 90th percentile of the population for LDL-C and TG, respectively. In 13 families three or more affected individuals had high polygenic score and in 22 families none of the affected individuals had polygenic score above the 90th percentile of the population.

Our results allow us to draw several conclusions about the genetic background of FCH. First, variants across the whole frequency spectrum were enriched and contribute to the high LDL-C and/or TG levels in the affected FCH individuals. This emphasizes the polygenic nature of FCH and is in line with previous results with familial hypercholesterolemia demonstrating that this less prevalent syndrome has a polygenic component [9]. Genes in which we observed enriched variants include APOE, LIPC and APOA5, whose role in dyslipidemias is already established, but also include several genes whose function in lipid metabolism is not yet clear (e.g. UBR1, MTHFD2L and the PIGV-NR0B2 region).

Second, the majority of the enriched variants were originally identified in random population samples. This finding confirms that variants identified in population screens play a considerable role also in the complex familial disease FCH, and highlights the potential of using variants identified in population-based GWA studies to characterize familial dyslipidemia cases genetically.

Third, the observed enrichment was stronger for TG loci than for LDL-C loci. Many of the enriched TG SNPs were located in genes known to contribute to hypertriglyceridemia, such as the APOA1-C3-A4-A5 cluster, showing the central role of genetically driven TG in FCH [7, 17].

Fourth, over a third of the affected FCH individuals had a high load (polygenic lipid scores over the 90th percentile of the population) of either known LDL-C or TG associated variants. This proportion with high polygenic score was comparable to the population samples with similar levels of high LDL-C or TG in the population. This observation suggests that the role of the known variants is similar in familial hyperlipidemias and randomly ascertained high LDL-C and/or TG individuals. However, the accumulation of LDL-C or TG-elevating alleles in families was highlighted by the prediction analysis results. High load of known LDL-C or TG associated variants predicted affectedness in FCH families almost two times better than comparable hyperlipidemia in the population.

There were large differences among families in the proportion of affected individuals with high polygenic scores. While in nine families more than two thirds of the affected individuals had high polygenic score, in over a third of the families, none of the cases had high score and the underlying genetic architecture remains unexplained in the majority of families. As genotyping and imputation allow for the evaluation of only a portion of catalogued high-impact variants, future sequencing studies might pinpoint rare causal variants in these families.

Our study supports the hypothesis that FCH is a genetically heterogeneous disease, reflected by its heterogeneous lipid phenotype [18-20]. large number of LDL-C associated variants are mainly responsible for the elevation of LDL-C in the FCH lipid profile. In contrast, a handful of low- to moderate-impact TG variants drive the elevation of triglycerides. Our study provides support for the polygenic rather than monogenic nature of FCH, and highlights the central role of genetically driven TG in FCH.

Materials and Methods

Ethics statement

Written informed consent was obtained from all study participants. All samples were collected in accordance with the Helsinki declaration and study protocols were approved by the ethics committees of the participating centers (The Hospital District of Helsinki and Uusimaa Coordinating Ethics Committee, approval number 184/13/03/00/12).

Subjects and measurements

As part of the European Multicenter Study on Familial Dyslipidemias in Patients with Premature Coronary Heart Disease (EUFAM), the Finnish FCH families were identified from patients admitted to university hospitals with a diagnosis of premature CHD who demonstrated levels of TC, TG, or both that were ≥ 90th Finnish age- and sex-specific population percentile (S4 Table) [21, 22]. Families that contained at least one other first-degree relative affected with hyperlipidemia (according to the same lipid-level criteria as used for the probands) were included in the study. Additionally, in the included families, at least one of the affected individuals had high TG. All family members with hyperlipidemia were then considered affected by FCH. Probands with a diagnosis of FH (screened with a functional LDL receptor test), and any subjects with diabetes or other chronic diseases were considered unaffected and did not contribute to establishing the family's FCH status (S1 Text).

For analyses of continuous lipid traits, individuals using lipid-lowering or estrogen medication at the time of sampling were excluded. Samples from the Finnish National FINRISK study were used as a Finnish population-specific comparison group, and individuals with known diabetes or cancer were excluded from the analyses (S1 Text). For the FCH families, venous blood samples were obtained after an overnight fast. FINRISK participants were advised to fast for four hours before the examination and avoid heavy meals earlier during the day. For both the EUFAM and FINRISK samples, circulating biochemical markers were measured from the venous blood samples using standard methods, as described in the S1 Text.

Genotyping and imputation

All FCH individuals with DNA available (715 out of 1161 from all 53 FCH families) and the FINRISK samples (n = 20,626) were genotyped with the HumanCoreExome BeadChip (Illumina Inc., San Diego, CA, USA) using standard methods (S1 Text). Additionally, over nine million variants were imputed using a combined reference panel of 1000 Genomes Phase I integrated haplotypes and 1943 Finnish genomes (S1 Text).

Polygenic lipid score calculation

To construct a polygenic lipid score, we catalogued and used 212 SNPs, representing either the lead SNPs for genome-wide significant associations to LDL-C or TG from GWA studies of these traits (171 SNPs), or SNPs catalogued in the Online Mendelian Inheritance in Man (OMIM) database as located in genes implicated in primary and secondary monogenic dyslipidemic syndromes (44 SNPs, four of which overlap with the GWA lead-SNPs, S2 Table) [23-26]. The SNPs were assigned to LDL-C and TG scores based on their previously reported associations (S1 Text). The scores were calculated as the sum of the risk alleles weighted by their effect estimates drawn from a multiple linear model estimating all SNP effects on the trait (LDL-C or TG) at the same time in the FINRISK samples (S1 Text). The FINRISK samples used to estimate the weights were independent from the FCH samples.

Statistical analysis

We estimated allele frequencies for all 212 SNPs in the 234 affected FCH individuals and in 18,715 FINRISK samples after excluding individuals with known diabetes or cancer similarly to the EUFAM exclusion criteria. The minor allele in FINRISK was designated as the effect allele. We then calculated enrichment ratios by dividing the effect allele frequencies in affected FCH individuals with the effect allele frequencies in the FINRISK samples. Under the null hypothesis of no enrichment, we estimated the null distribution for enrichment testing by calculating the enrichment statistic (enrichment ratio) for all variants with MAF > 0.001% (12,234,754 SNPs) across the genome excluding regions within 50 000 base pairs from the lipid-associated SNPs. We estimated the null distribution in MAF bins across the genome while keeping the family structure fixed (FINRISK minor allele frequency intervals [0.1%,0.5%); [0.5%,1%); [1%,1.5%); [1.5%,2%); [2%,2.5%), [2.5%,5%); [5%,10%); [10%,15%); [15%,20%); [20%,25%); [25%,30%); [30%,35%); [35%,40%); [40%,45%); and [45%,50%)). The observed lipid SNP allele frequency enrichment ratios were compared to the null distribution to provide one-sided p-values for each SNP. We used sign test for all SNPs together to estimate whether the direction of enrichment (more or less common in the affected FCH individuals than in the population) was associated with the SNPs' effects on the lipid in question (elevating or lowering).

We applied linear mixed models to test for differences in metabolic and clinical characteristics between FCH affected and unaffected individuals (S1 Text). An empirical genetic correlation matrix between individuals was included as the covariance structure of a random effect. Linear mixed models were applied with MMM (version 1.01) [27], and the other statistical analyses were performed using R (version 3.2.1) [28].

Supporting Information

S1 Text. Supplementary Materials and Methods.

The recruitment, assessment, and genotyping of study subjects, the calculation of polygenic lipid scores, and the statistical analyses performed are described in more detail.

(DOCX)

S1 Fig. Number of affected and unaffected individuals in FCH families.

The number of genotyped affected (black) and unaffected individuals (white) is presented for each FCH family. The families are sorted by the number of affected individuals in each family, but the family numbers on the x-axis correspond to those presented in Fig 3.

(PDF)

S2 Fig. Lipid level median values and ranges in FCH families in genotyped individuals.

Median values of (a) LDL-C and (b) TG for genotyped individuals in all families (n = 53). Vertical lines represent the range of values within the families. The families are sorted by the medians of LDL-C and TG levels, respectively, but the family numbers on the x-axis correspond to those presented in Fig 3. In (b) the y-axis is cut at 7 mmol/l.

(PDF)

S3 Fig. Enrichment of LDL-C or TG associated SNPs in FCH probands and all FCH family members by their frequency and effect on LDL-C and TG levels.

Enrichment ratio is the ratio of effect allele frequency in FCH probands (n = 48, (a), (b)), or in all FCH family members (n = 715, (c), (d)), to allele frequencies in the Finnish FINRISK population cohort (n = 20,626 in c) and d); n = 18,715 in a) and b) after excluding individuals with diabetes and cancer). Under the null hypothesis of no enrichment, a 95% credible interval (shaded area) was estimated by calculating the enrichment statistic (enrichment ratio) for all variants with MAF > 0.001% across the genome excluding the loci of the 212 SNPs. Variants are designated as either lipid level elevating (red) or lowering (blue) for LDL-C ((a), (c)) and TG ((b), (d)) based on [Beta] estimates from linear regression in the FINRISK samples. Point size and color intensity reflect the magnitude of the effect. Only SNPs with at least one heterozygous carrier are shown (n = 191 for probands and n = 196 for all family members). FINRISK, The National FINRISK Study. *The enrichment ratio for LPL rs1801177 fell within the 95% credible interval.

(PDF)

S4 Fig. Triglyceride levels in carriers and wild type individuals of APOA5 rs3135506 and APOE rs7412.

Homozygosity for (a) APOA5 rs3135506 predisposes to hypertriglyceridemia and homozygosity for (b) APOE rs7412 predisposes to type III hyperlipoproteinemia, which typically presents in elevated levels of VLDL and TG species. Only individuals without diabetes and other relevant confounders were included in this comparison. In (b) the y-axis is log-scaled.

(PDF)

S5 Fig. Mean genetic correlation between affected family members in FCH families.

For each family with at least two genotyped affected members (n = 51), we calculated mean empiric genetic correlation of all combinations of pairs of affected subjects. Families are ranked on the x-axis as in Fig 3, with higher ranking representing a higher number of affected subjects with high polygenic lipid scores or high-impact Mendelian variants. A significance estimate for the relationship was derived from simple linear regression.

(PDF)

S1 Table. Clinical and metabolic characteristics of all individuals.

P-values were calculated using Wald test by a linear mixed model correcting for sample relatedness. FCH, familial combined hyperlipidemia.

(PDF)

S2 Table. LDL-C or TG associated variants included in the study.

NFE, non-Finnish Europeans; FINRISK, The National FINRISK Study; FCH, familial combined hyperlipidemia. *Effect allele. [dagger]Global minor allele. [dagger]MAF in non-Finnish Europeans calculated from the public ExAC data set (version 0.3) for coding SNPs (those present in the ExAC data set), and from the public 1000 Genomes Phase 3 data set for non-coding SNPs (S1 Text). §MAF in the FINRISK cohort after excluding those with known diabetes or cancer. ||Common: MAF > 5%, low-frequency: 0.5% < MAF ≤ 5%, rare: MAF ≤ 0.5%. ¶SNPs were selected based on OMIM entries for genes implicated in monogenic dyslipidemia, and recent genome-wide association studies (S1 Text). [diamonds]Effect allele weight not estimated due to only one copy of minor allele present in FINRISK.

(PDF)

S3 Table. Enrichment ratios and observed p-values for all score SNPs with at least one affected carrier.

Under the null hypothesis of no enrichment, we established the null distribution for enrichment testing by calculating the enrichment statistic (enrichment ratio) for all variants with effect allele (minor allele in the FINRISK cohort) frequency > 0.001% across the genome excluding the loci of the 212 SNPs. Enrichment ratios were calculated for each of the 212 SNPs if at least one carrier was identified (194 SNPs for affected individuals, 191 SNPs for probands), and one-tailed p-values were estimated from the null distribution. Enrichment ratios are presented for the affected individuals. The carrier counts were calculated from all genotyped FCH family members who passed exclusion criteria and did not have diabetes or other relevant confounders (n = 661), and did not account for non-independence within family. The + and--symbols denote elevating and lowering effects on the lipid in question, respectively. FCH, familial combined hyperlipidemia; FINRISK, The National FINRISK Study. J.A.K, J. A. Kuivenhoven; C.J.W., C. J. Willer; T.M.T., T. M. Teslovich; I.S., I. Surakka (S1 Text). [dagger]Effect allele frequencies in non-Finnish Europeans were calculated from the public ExAC data set (version 0.3) for coding SNPs (those present in the ExAC data set), and from the public 1000 Genomes Phase 3 data set for non-coding SNPs (S1 Text). [dagger]Effect allele frequencies in the FINRISK cohort after excluding those with known diabetes or cancer. §Enrichment ratio is the ratio of the effect allele frequencies in affected individuals (n = 234) to the Finnish population (n = 18,715).

(PDF)

S4 Table. The 90th age- and sex-specific Finnish population percentiles for total cholesterol and triglycerides.

In the EUFAM study, dyslipidemia was established based on levels of total cholesterol, triglycerides, or both that were ≥ 90th Finnish age- and sex-specific population percentile. The population percentiles were derived from FINMONICA, a large population survey performed in 1992. The percentile estimates and a description of the polynomial regression analyses used to establish them have been reported in detail previously [21].

(PDF)

Acknowledgments

We would like to thank Sari Kivikko and Huei-Yi Shen for management assistance.

Author Contributions

Conceived and designed the experiments: SR MRT NBF AP RKW VS CE PP. Performed the experiments: SR AP MRT NBF PPL APS SKS RKW VS. Analyzed the data: PR JTR IS. Contributed reagents/materials/analysis tools: MRT SS NM VS MP SR SKS AP NBF. Wrote the paper: PR JTR SR MRT AP NBF SS SKS. Provided support in writing of the manuscript: IS MP APS.

References

Goldstein JL, Schrott HG, Hazzard WR, Bierman EL, Motulsky AG. Hyperlipidemia in coronary heart disease. II. Genetic analysis of lipid levels in 176 families and delineation of a new inherited disorder, combined hyperlipidemia. J Clin Invest. 1973;52(7):1544-68. doi: 10.1172/JCI107332. pmid:4718953

Nikkila EA, Aro A. Family study of serum lipids and lipoproteins in coronary heart-disease. Lancet. 1973;1(7810):954-9. pmid:4121585 doi: 10.1016/s0140-6736(73)91598-5

Voors-Pette C, de Bruin TW. Excess coronary heart disease in Familial Combined Hyperlipidemia, in relation to genetic factors and central obesity. Atherosclerosis. 2001;157(2):481-9. pmid:11472750 doi: 10.1016/s0021-9150(00)00752-8

Genest JJ Jr., Martin-Munley SS, McNamara JR, Ordovas JM, Jenner J, Myers RH, et al. Familial lipoprotein disorders in patients with premature coronary artery disease. Circulation. 1992;85(6):2025-33. pmid:1534286 doi: 10.1161/01.cir.85.6.2025

European Association for Cardiovascular Prevention & Rehabilitation, Reiner Z, Catapano AL, De Backer G, Graham I, Taskinen MR, et al. ESC/EAS Guidelines for the management of dyslipidaemias: the Task Force for the management of dyslipidaemias of the European Society of Cardiology (ESC) and the European Atherosclerosis Society (EAS). Eur Heart J. 2011;32(14):1769-818. doi: 10.1093/eurheartj/ehr158. pmid:21712404

Brahm AJ, Hegele RA. Combined hyperlipidemia: familial but not (usually) monogenic. Curr Opin Lipidol. 2015. doi: 10.1097/MOL.0000000000000270

Brouwers MC, van Greevenbroek MM, Stehouwer CD, de Graaf J, Stalenhoef AF. The genetics of familial combined hyperlipidaemia. Nat Rev Endocrinol. 2012;8(6):352-62. doi: 10.1038/nrendo.2012.15. pmid:22330738

Johansen CT, Wang J, Lanktree MB, McIntyre AD, Ban MR, Martins RA, et al. An increased burden of common and rare lipid-associated risk alleles contributes to the phenotypic spectrum of hypertriglyceridemia. Arterioscler Thromb Vasc Biol. 2011;31(8):1916-26. doi: 10.1161/ATVBAHA.111.226365. pmid:21597005

Talmud PJ, Shah S, Whittall R, Futema M, Howard P, Cooper JA, et al. Use of low-density lipoprotein cholesterol gene score to distinguish patients with polygenic and monogenic familial hypercholesterolaemia: a case-control study. Lancet. 2013;381(9874):1293-301. doi: 10.1016/S0140-6736(12)62127-8. pmid:23433573

Hegele RA, Ginsberg HN, Chapman MJ, Nordestgaard BG, Kuivenhoven JA, Averna M, et al. The polygenic nature of hypertriglyceridaemia: implications for definition, diagnosis, and management. Lancet Diabetes Endocrinol. 2014;2(8):655-66. doi: 10.1016/S2213-8587(13)70191-8. pmid:24731657

Stitziel NO, Peloso GM, Abifadel M, Cefalu AB, Fouchier S, Motazacker MM, et al. Exome sequencing in suspected monogenic dyslipidemias. Circ Cardiovasc Genet. 2015;8(2):343-50. doi: 10.1161/CIRCGENETICS.114.000776. pmid:25632026

Talmud PJ, Hawe E, Martin S, Olivier M, Miller GJ, Rubin EM, et al. Relative contribution of variation within the APOC3/A4/A5 gene cluster in determining plasma triglycerides. Hum Mol Genet. 2002;11(24):3039-46. pmid:12417525 doi: 10.1093/hmg/11.24.3039

Zhao SP, Smelt AH, Van den Maagdenberg AM, Van Tol A, Vroom TF, Gevers Leuven JA, et al. Plasma lipoprotein profiles of normocholesterolemic and hypercholesterolemic homozygotes for apolipoprotein E2(Arg158-->Cys) compared. Clin Chem. 1994;40(8):1559-66. pmid:8044997

Brunzell JD, Zambon A, Deeb SS. The effect of hepatic lipase on coronary artery disease in humans is influenced by the underlying lipoprotein phenotype. Biochim Biophys Acta. 2012;1821(3):365-72. doi: 10.1016/j.bbalip.2011.09.008. pmid:21986251

Hegele RA, Vezina C, Moorjani S, Lupien PJ, Gagne C, Brun LD, et al. A hepatic lipase gene mutation associated with heritable lipolytic deficiency. J Clin Endocrinol Metab. 1991;72(3):730-2. doi: 10.1210/jcem-72-3-730. pmid:1671786

Vuorio AF, Aalto-Setälä K, Koivisto UM, Turtola H, Nissen H, Kovanen PT, et al. Familial hypercholesterolaemia in Finland: common, rare and mild mutations of the LDL receptor and their clinical consequences Finnish FH-group. Ann Med. 2001;33:410-21. pmid:11585102 doi: 10.3109/07853890108995954

Calandra S, Priore Oliva C, Tarugi P, Bertolini S. APOA5 and triglyceride metabolism, lesson from human APOA5 deficiency. Curr Opin Lipidol. 2006;17(2):122-7. doi: 10.1097/01.mol.0000217892.00618.54. pmid:16531747

Jarvik GP, Brunzell JD, Austin MA, Krauss RM, Motulsky AG, Wijsman E. Genetic predictors of FCHL in four large pedigrees. Influence of ApoB level major locus predicted genotype and LDL subclass phenotype. Arterioscler Thromb. 1994;14(11):1687-94. pmid:7947591 doi: 10.1161/01.atv.14.11.1687

van Greevenbroek MM, Stalenhoef AF, de Graaf J, Brouwers MC. Familial combined hyperlipidemia: from molecular insights to tailored therapy. Curr Opin Lipidol. 2014;25(3):176-82. doi: 10.1097/MOL.0000000000000068. pmid:24811296

Lewis GF, Xiao C, Hegele RA. Hypertriglyceridemia in the genomic era: a new paradigm. Endocr Rev. 2015;36(1):131-47. doi: 10.1210/er.2014-1062. pmid:25554923

Porkka KV, Nuotio I, Pajukanta P, Ehnholm C, Suurinkeroinen L, Syvanne M, et al. Phenotype expression in familial combined hyperlipidemia. Atherosclerosis. 1997;133(2):245-53. pmid:9298685 doi: 10.1016/s0021-9150(97)00134-2

Borodulin K, Vartiainen E, Peltonen M, Jousilahti P, Juolevi A, Laatikainen T, et al. Forty-year trends in cardiovascular risk factors in Finland. Eur J Public Health. 2015;25(3):539-46. doi: 10.1093/eurpub/cku174. pmid:25422363

Kuivenhoven JA, Hegele RA. Mining the genome for lipid genes. Biochim Biophys Acta. 2014;1842(10):1993-2009. doi: 10.1016/j.bbadis.2014.04.028. pmid:24798233

Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466(7307):707-13. Epub 2010/08/06. doi: 10.1038/nature09270. pmid:20686565

Surakka I, Horikoshi M, Magi R, Sarin AP, Mahajan A, Lagou V, et al. The impact of low-frequency and rare variants on lipid levels. Nat Genet. 2015;47(6):589-97. doi: 10.1038/ng.3300. pmid:25961943

Global Lipids Genetics Consortium, Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45:1274-83. doi: 10.1038/ng.2797. pmid:24097068

Pirinen M, Donnelly P, Spencer CCA. Efficient computation with a linear mixed model on large-scale data sets with applications to genetic studies. Ann Appl Stat. 2013;7(1):369-90. doi: 10.1214/12-AOAS586

R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; 2014