- Research article
- Open Access
Six years’ experience with LipidSeq: clinical and research learnings from a hybrid, targeted sequencing panel for dyslipidemias
BMC Medical Genomics volume 13, Article number: 23 (2020)
In 2013, our laboratory designed a targeted sequencing panel, “LipidSeq”, to study the genetic determinants of dyslipidemia and metabolic disorders. Over the last 6 years, we have analyzed 3262 patient samples obtained from our own Lipid Genetics Clinic and international colleagues. Here, we highlight our findings and discuss research benefits and clinical implications of our panel.
LipidSeq targets 69 genes and 185 single-nucleotide polymorphisms (SNPs) either causally related or associated with dyslipidemia and metabolic disorders. This design allows us to simultaneously evaluate monogenic—caused by rare single-nucleotide variants (SNVs) or copy-number variants (CNVs)—and polygenic forms of dyslipidemia. Polygenic determinants were assessed using three polygenic scores, one each for low-density lipoprotein cholesterol, triglyceride, and high-density lipoprotein cholesterol.
Among 3262 patient samples evaluated, the majority had hypertriglyceridemia (40.1%) and familial hypercholesterolemia (28.3%). Across all samples, we identified 24,931 unique SNVs, including 2205 rare variants predicted disruptive to protein function, and 77 unique CNVs. Considering our own 1466 clinic patients, LipidSeq results have helped in diagnosis and improving treatment options.
Our LipidSeq design based on ontology of lipid disorders has enabled robust detection of variants underlying monogenic and polygenic dyslipidemias. In more than 50 publications related to LipidSeq, we have described novel variants, the polygenic nature of many dyslipidemias—some previously thought to be primarily monogenic—and have uncovered novel mechanisms of disease. We further demonstrate several tangible clinical benefits of its use.
Dyslipidemias, defined as extreme deviations of plasma lipids or lipoproteins, are commonly encountered clinically . They are often associated with increased risk of cardiovascular disease and other complications such as acute pancreatitis [2, 3]. There are 25 monogenic dyslipidemias caused by variants in 25 genes [1, 4, 5], most of which were identified > 10 years ago using classical biochemical and genetic mapping methods . With the exception of heterozygous familial hypercholesterolemia (FH), monogenic dyslipidemias are rare disorders  and can sometimes display multisystem syndromic features [1, 4]. Most show recessive inheritance and typically result from pathogenic rare variants—either single-nucleotide variants (SNVs) or copy-number variants (CNVs)—in well-established causal genes. Further, some dyslipidemias are polygenic, resulting from contributions of several types of genetic determinants including incompletely penetrant rare variants and small-effect common variants [2, 7]. Accumulated variants within an individual’s genome can predispose to more severe phenotypic expression . In addition to genetic determinants, several secondary factors—diet, obesity, activity level, other medical conditions such as diabetes or hypothyroidism, and certain medications—can exacerbate the clinical presentation of both monogenic and polygenic dyslipidemias [8, 9].
For over 25 years, our laboratory has studied both monogenic and polygenic dyslipidemias. Patient care and genetic analysis have coexisted through fortuitous geographic convergence of our lipid clinic, genetics research laboratory and genomic core facility, and through uninterrupted funding for a genetics research program. At the time the clinic and research program were established, our ethics review panel stipulated that genetic results were to be shared with patients, and this became our practice. Patient samples come from both local clinical practices and international colleagues; virtually all patients seen in the clinic have consented to provide DNA samples for research. Between 1998 and 2012, DNA analysis was performed by automated Sanger sequencing. In 2013, we transitioned to next-generation sequencing using the custom-designed “LipidSeq” panel; results from the latter are reported here.
Because our clinical practice spans all dyslipidemias, we have focused on their ontology [1, 10, 11] and on documenting dyslipidemia-associated gene variants . Our molecular classification system ultimately informed the design of the LipidSeq panel for genes underlying monogenic dyslipidemias [12, 13]. We also designed the panel to target genes causing monogenic disorders for which dyslipidemia is a secondary manifestation, such as inherited forms of diabetes. A benefit of the high depth of coverage afforded by our panel is the ability to concurrently identify CNVs along with SNVs. Furthermore, our longstanding interest in the polygenic basis of plasma lipids [14,15,16] motivated us to simultaneously genotype common single-nucleotide polymorphisms (SNPs) . We easily accommodated 185 “micro-sequencing” reactions to genotype lipid trait-associated SNPs from the Global Lipid Genetics Consortium genome-wide association studies (GWAS) of plasma lipids [17,18,19].
Thus, LipidSeq is a hybrid panel that detects: 1) functionally relevant rare SNVs and CNVs in genes underlying monogenic dyslipidemias, and 2) common variants, particularly SNPs, that we use to build polygenic scores . Since 2014, LipidSeq results have helped clarify the genetic basis for hundreds of dyslipidemic patients and have been reported in > 50 peer-reviewed publications. In this report, we briefly describe our aggregated research findings and discuss the clinical benefit of our LipidSeq panel.
The LipidSeq panel
LipidSeq was designed for clinical resequencing of genomic loci associated with dyslipidemia and related metabolic traits [4, 12]. It targets exons plus 250 bp into each flanking intron and the 5′- and 3′-untranslated regions of 69 genes, including: 1) 25 causative genes for monogenic dyslipidemias; 2) 16 causative genes for inherited lipodystrophies; 3) 13 for subtypes of maturity-onset diabetes of the young (MODY) and inherited diabetes; and 4) 15 candidate genes in lipoprotein metabolism for which no pathogenic mutations have yet been found (Additional file 1: Table S1). LipidSeq also “micro-sequences” the area surrounding 185 GWAS SNPs, from which we use certain SNP subsets to build small-scale polygenic scores for low-density lipoprotein (LDL) cholesterol, triglyceride, and high-density lipoprotein (HDL) cholesterol . An expanded rationale for the LipidSeq design are reported elsewhere , as well as quality assessment and validation of the panel .
Figure 1 shows the breakdown of samples studied using LipidSeq. The Lipid Genetics Clinic at the London Health Sciences Centre, University Hospital (London ON, Canada) was established in 1998 and operates a half-day each week, serving a region with a population of ~ 2 million people and providing care to outpatients referred from family practitioners and specialists. Because the main clinic physician (R.A.H.) also holds peer-reviewed research funding to study genetics of dyslipidemia and related disorders, patients are invited to provide DNA samples for research. The journey of a patient and their consented DNA sample are outlined in Fig. 2. On initial consultation (Visit 1), after taking a medical history and performing a physical examination, informed consent is obtained and the patient provides a fasting blood sample for: 1) determination of routine lipid profile (total, LDL and HDL cholesterol, and triglyceride); 2) advanced lipid profiling (including apolipoprotein [apo] B and A1, and lipoprotein [a]); 3) screening for secondary causes of dyslipidemia (including hypothyroidism, and hepatic and renal conditions); 4) screening non-traditional cardiovascular risk markers (including high sensitivity C-reactive protein and carotid intima-medial thickness); and 5) DNA extraction. After sample processing and reporting, results are discussed with the patient at Visit 2. The project was first approved in 1998 by the Western University ethics review board (#07290E) and has been updated and renewed continuously.
We also accept patient samples referred from colleagues provided consent is obtained following appropriate institutional standards. We also receive contracted samples from pharmaceutical companies to genotype de novo or validate previous diagnostic results for clinical trials. External samples follow the same processing flow as internal samples (Fig. 2).
DNA extraction and isolation
Genomic DNA is extracted from blood (95% of samples) using the Puregene® DNA Blood Kit (Gentra Systems, Qiagen Inc., Mississauga ON, Canada) (Cat No. 158389) or saliva (5% of samples) using the Oragene DNA kit (DNA Genotek Inc., Ottawa ON, Canada; Cat No. OG-500).
Sample preparation and sequencing
DNA samples prepared in batches of 24 are indexed and enriched using the Nextera® Rapid Capture Custom Enrichment Kit (Cat No. FC-140-1009) “LipidSeq” design . Sequencing is performed for each batch at the London Regional Genomics Centre (www.lrgc.on.ca; London ON, Canada) on an Illumina MiSeq personal sequencer (Illumina, San Diego CA, USA).
Bioinformatic processing and quality assurance
Paired FASTQ files are generated for each sample after sequencing and are imported into CLC Bio Genomics Workbench (CLC Bio, Aarhus, Denmark) for bioinformatic processing. Sequencing reads are aligned to the human reference genome (hg19/GRCh37) and undergo local realignment to improve alignment quality. From the assembled reads, variants are called if there are discrepancies between the reference genome and the sample’s sequence. Following this, VCF and BAM files are created for each sample; these files contain information on the genomic position and zygosity of identified variants, as well as the depth of coverage for each sequencing read. A detailed explanation of the bioinformatic and quality assurance processes have already been reported .
The initial tool used for the annotation and analysis of variants was the open-sourced tool, ANNOVAR . Recently, we have switched to the commercially available software, VarSeq® (Golden Helix, Inc., Bozeman MT, USA) for variant annotation and analysis. Our software upgrade allowed us to assess for CNVs, which was not previously accessible using ANNOVAR alone.
Single-nucleotide variants and indels
Rare variants with potential for protein-altering effects are of primary interest. We consider variants with a minor allele frequency of ≤1% or absent from publicly available genotype databases. Our reference database has changed over the years as more comprehensive databases became available; starting with the 1000 Genomes Project (http://browser.1000genomes.org/index.html) , we progressed to the Exome Aggregation Consortium (ExAC; http://exac.broadinstitute.org/) , and finally to the Genome Aggregation Database (gnomAD; https://gnomad.broadinstitute.org/) . Only rare variants that impact the amino acid sequence or canonical splice sites are considered, including missense, nonsense, insertions or deletions (indels), splice-donor, and splice-acceptor variants; the rationale for this is that a change to the encoded protein will likely have a phenotypic impact. In an attempt to avoid benign variants, multiple in silico prediction tools are used to identify rare variants with possible damaging or deleterious effects, including Combined Annotation Dependent Depletion (CADD; http://cadd.gs.washington.edu/score) [26, 27], Polymorphism Phenotyping version 2 (PolyPhen2; http://genetics.bwh.harvard.edu/pph2/) , Sorting Intolerant From Tolerant (SIFT; http://sift.jcvi.org/) , and MutationTaster (http://www.mutationtaster.org/) . Since its introduction in 2015, we also consider ACMG classifications for each variant of interest, and are in the process of reannotating our entire variant database using these criteria . We have been utilizing Franklin by Genoox (https://franklin.genoox.com/home), a web tool for variant interpretation for this process.
CNVs are detected using the VarSeq-CNV® caller algorithm. Using BAM files, this algorithm detects differences in read depth of a sample compared to a group of “reference” samples without CNVs. More details of this method and our standard parameters were reported previously .
We calculate small polygenic scores using lipid-altering alleles from a subset of SNP loci captured by our panel; 10, 16, and 9 SNPs comprise the LDL cholesterol, triglyceride and HDL cholesterol scores, respectively. All targeted SNPs were reported by the Global Lipids Genetics Consortium as having a statistically significant association with at least one of the three lipid traits [17,18,19]. Our weighted polygenic score calculation considers the total number of trait-raising alleles at a single locus (0, 1, or 2) multiplied by that allele’s beta coefficient determined from GWAS [17,18,19]. Each product is summed to produce the overall weighted polygenic score for the trait. A more detailed explanation behind polygenic scores and their calculations is available . Each individual sequenced by LipidSeq has polygenic scores calculated for each of LDL cholesterol, triglyceride and HDL cholesterol, regardless of their referral phenotype. An extreme accumulation of common SNPs was defined as an extreme polygenic score, classified as a score greater than or equal to the 90th percentile previously determined using a normolipidemic reference group .
Reporting clinically relevant genetic determinants of interest
The preliminary list of computationally prioritized rare variants from each patient sample—either SNVs, indels, or CNVs—with potentially damaging or deleterious effects derived from the pipeline is first checked and reviewed manually by two laboratory personnel (A.D.M. and J.W.). Polygenic scores are also shown on the draft patient report. Prior to Visit 2, the patient’s draft report is reviewed by the laboratory scientist (A.D.M.) and the physician (R.A.H.) before finalizing the report of both rare variant results and polygenic scores, with interpretations derived by consensus from the three reviewers (A.D.M., J.W. and R.A.H.). We put particular emphasis on rare variants disrupting genes with a direct relation to the phenotype of interest. With respect to polygenic scores, for brevity and simplicity, we only report to the patient the score associated with their referral phenotype (i.e. a patient with hypertriglyceridemia will only receive results from the triglyceride risk score). Based on a patient’s report, they can receive either: a) a genetically based diagnosis; b) a genetic confirmation of a previously received diagnosis; c) a “relevant” genetic result (i.e. a variant that has not been functionally confirmed to cause the phenotype, but is still predicted to be damaging and occurs within a phenotypically associated gene); or d) a negative result, indicating that we were unable to identify any sort of genetic determinant related to the phenotype. Each report is proofread and signed by the laboratory director (R.A.H.). Hard copies of reports for Lipid Genetics Clinic patients are added to patient paper charts and the findings are reported to the patient at Visit 2. Hard copies of reports for externally referred samples are mailed to the referring physician.
Characterization of sequenced samples
To date, we have sequenced 3262 samples from both internal and external sources (Fig. 1), of which 1466 (45.0%) were from the Lipid Genetics Clinic and 1796 (55.0%) were received from external sources for diagnostic (32.0%), research (16.2%), and industry-contracted (6.8%) purposes. Demographic and clinical information from our cohort of internal patient samples is presented in Table 1. The phenotypic breakdown of our total sample cohort is illustrated in Fig. 3.
The most prevalent phenotype is hypertriglyceridemia (40.1%), followed by FH (28.3%). Briefly, patients with hypertriglyceridemia have elevated triglyceride levels (≥ 1.8 mmol/L) and can present with different clinical features depending on whether the patient has a mild-to-moderate (> 1.8 and < 10 mmol/L) or severe (≥ 10 mmol/L) deviation [33, 34]. These patients are referred to clinic to identify a possible genetic basis for their condition, and for recommendation of treatment options. In contrast, patients referred with “FH” have high prior clinical suspicion of this condition, mainly due to extremely elevated LDL cholesterol levels (> 5.0 mmol/L). With our LipidSeq panel, we are often able to provide a genetic confirmation of the FH diagnosis, which in turn may support the use of more intensive therapeutic strategies to lower LDL cholesterol levels and decrease risk for cardiovascular disease.
The remaining 31.6% of samples include: low or high HDL cholesterol levels (i.e. hypo- and hyperalphalipoproteinemia, respectively), low LDL cholesterol levels (hypobetalipoproteinemia/abetalipoproteinemia), familial combined hyperlipidemia, diabetes, lipodystrophy, and miscellaneous conditions including elevated levels of lipoprotein(a), atypical hyperlipidemia, and severe obesity (Additional file 1: Table S2). Patients referred with low HDL cholesterol levels may be at an increased risk for cardiovascular disease [35,36,37]. Conversely, patients with high HDL cholesterol levels were previously thought to be at a decreased risk for cardiovascular disease; however, in some instances the causative molecular mechanism increases risk due to impaired clearance of HDL particles . Patients with familial combined hyperlipidemia have elevations of both LDL cholesterol and triglyceride levels, which increases risk for cardiovascular disease. Meanwhile, patients with diabetes, insulin resistance, or uncontrolled glucose are usually referred to clinic for assistance in management of the dyslipidemic component of their phenotype. Patients referred to us with a clinical suspicion of partial lipodystrophy are often able to receive a genetic confirmation of this diagnosis. Meanwhile, patients with lipoprotein(a) levels in the top 5th percentile of the population (i.e. ≥ 36 mg/dL) are referred to our clinic for assistance in managing modifiable cardiovascular risk factors, since no treatment is presently available. “Atypical hyperlipidemia” is used to describe patients with multiple perturbations of lipid variables that do not fit in the “familial combined hyperlipidemia” category. Patients referred due to severe obesity often have dyslipidemia and diabetes-related complications requiring management.
Rare variant analysis
A total of 24,931 unique variants were identified across 3262 samples sequenced with LipidSeq (Additional file 1: Table S3). After applying our rare variant filtering criteria (Fig. 4a), 2205 variants were of potential interest based on their disruptive sequence ontology and predictions of producing a deleterious or damaging protein product. Of these, 289 are predicted to be loss-of-function variants based on ontologies of either frameshift, splice acceptor, splice donor, nonsense (stop gain), or stop loss (Fig. 4b). After filtering these variants further based on a stricter CADD PHRED score of ≥20, 258 variants would likely be of clinical interest given their predicted level of having a damaging effect (Additional file 1: Table S4), especially if the variant disrupts a gene with a known relationship to the patient’s phenotype (ex. an LDLR variant in a patient with FH). Lastly, from our rare variant findings, we identified 191 total CNVs, of which 77 were unique (Table 2). The distributions of CADD PHRED-scaled scores for all variants predicted to be possibly deleterious or damaging are presented in Fig. 4c.
Genetic characterization of familial hypercholesterolemia and hypertriglyceridemia
We show the relevant genetic determinants present in the patient subsets for the two most prevalent phenotypes encountered in the clinic (Table 3). Rare SNVs and indels were considered if they occurred in phenotypically relevant genes (i.e. LDLR, PCSK9, or APOB for patients with FH; LPL, APOA5, LMF1, GPIHBP1, or APOC2 for patients with hypertriglyceridemia) had a CADD PHRED-scaled score ≥ 10 plus a predicted deleterious or damaging outcomes by SIFT, PolyPhen2, or MutationTaster, and resulted in a change to the encoded protein’s amino acid sequence. The CNVs described in Table 2 were also used in this characterization.
When considering rare variants—both SNVs and CNVs—and extreme common SNP accumulation, FH patients were more likely to carry a rare variant compared to hypertriglyceridemia patients (46.3% vs. 23.9%), while hypertriglyceridemia patients were more likely to have an extreme accumulation of common SNPs, represented as an extremely high polygenic score compared to FH patients (32.7% vs. 12.4%). Overall, each cohort had ~ 50% of patients with an identifiable, relevant genetic determinant, although the most prominent determinant for FH patients was the presence of a rare variant, while an extreme polygenic score was the most prominent feature in hypertriglyceridemia patients. We are involved in updating ACMG pathogenicity criteria for FH-causing variants and will adjust our diagnostic process when these criteria are published.
We report our clinical and research experience with LipidSeq, a targeted hybrid panel designed for clinical resequencing of genomic loci known to be associated with dyslipidemia and related metabolic traits and disorders. Since 2014, the results from this panel have contributed to 39 publications reporting original scientific findings, including seven on FH [32, 43, 44, 46,47,48,49], seven on hypertriglyceridemia [42, 45, 50,51,52,53,54], four on extremes of HDL cholesterol [39, 55,56,57], and 21 case reports [40, 41, 58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76]. We have published an additional 15 reviews and methods articles related to this work [4, 5, 7, 11,12,13, 20, 34, 77,78,79,80,81,82,83]. Some highlights of outcomes from the use of LipidSeq are summarized in Tables 4 and 5. Several insights emerged, particularly from the 1466 samples acquired from patients of the Lipid Genetics Clinic who were referred for medical care and treatment advice. Sometimes, the research results could be applied directly to inform and guide patient management, especially when confirming a clinical diagnosis of FH and for other less common dyslipidemias (Table 4).
Perhaps the largest impact of DNA-based diagnosis has been upon patients with suspected FH; our laboratory is among the largest contributors of FH variants to the ClinVar database . In contrast to the low yield of FH-causing variants in population-based samples with hypercholesterolemia , we find that ~ 50% of referred patients suspected to have FH with LDL cholesterol > 5 mmol/L (> 190 mg/dL) had likely or definite pathogenic variants, which rose to > 90% for patients with LDL cholesterol > 8 mmol/L (> 310 mg/dL) . Furthermore, by simultaneously assessing for CNVs, we increased the diagnostic yield of likely pathogenic LDLR variants by ~ 10% [32, 77, 79, 80, 83]. When rare variants were absent, we found at least 20% of patients with suspected heterozygous FH had a high polygenic SNP score [7, 20, 44], indicating accumulated trait-raising alleles at SNP loci associated with LDL cholesterol.
In contrast to FH, most defined cases of severe hypertriglyceridemia (> 30%) were not monogenic, while only 1–2% of cases were diagnosed as familial chylomicronemia syndrome due to biallelic pathogenic variants affecting lipolysis . Among individuals with this monogenic, autosomal recessive condition, there are minimal phenotypic differences when stratifying by causative gene or type of genetic determinant . Among patients with monogenic chylomicronemia, ~ 5% of causative variants were CNVs in the GPIHBP1 gene . While individuals with monogenic hypertriglyceridemia had higher relative risk of acute pancreatitis than those with multifactorial or polygenic hypertriglyceridemia , the absolute number of cases was larger in the latter group, since it is much more prevalent . We showed how the clinical phenotype in some patients with multifactorial hypertriglyceridemia can be as severe as in those with monogenic hypertriglyceridemia [63, 66, 71].
Among patients with severely lowered HDL cholesterol, 2–3% have monogenic disorders (i.e. recessive Tangier disease, LCAT deficiency or apo A-I deficiency) . As with severe hypertriglyceridemia, polygenic factors like heterozygous rare variants with incomplete penetrance and extreme polygenic SNP scores, were much more common among individuals with very low HDL cholesterol . Also, we detected heterozygous large-scale deletions of ABCA1 in four patients with severely lowered HDL cholesterol, the first report of ABCA1 CNVs in the context of this phenotype .
Beyond characterizing the genetic determinants underlying our patients’ phenotypes, we have also uncovered new mechanisms of disease. In two families with severe FH, we discovered a heterozygous whole-gene duplication of PCSK9 with extremely high circulating PCSK9 levels . As well, a gain-of-function mutation in APOC3 was revealed as a new cause for hypertriglyceridemia .
Our findings have also been individually impactful for our patients. DNA-based confirmation of the diagnosis of heterozygous FH has helped > 50 patients to obtain private coverage for PCSK9 inhibitor drugs. A pilot pharmacogenetic analysis showed that these agents are equally effective in monogenic and polygenic severe hypercholesterolemia . As well, we diagnosed several patients whose severe hypertriglyceridemia was due to subclinical undiagnosed partial lipodystrophy , which altered monitoring and management. Other examples of positive clinical outcomes from use of LipidSeq include: 1) ending the protracted diagnostic odyssey endured by some patients [72, 74, 75]; 2) increasing the diagnostic yield in MODY diabetes by ~ 6% through simultaneous screening for CNVs [40, 76]; 3) switching some patients with GCK CNVs (diagnosed with MODY2) from insulin to oral hypoglycemic agents ; 4) diagnosing sitosterolemia in patients who were initially diagnosed with homozygous FH, resulting in a dramatic change in management ; and 5) ruling out genetic contributions in several patients with severe dyslipidemias due to secondary causes [52, 60, 65].
Sharing research results with patients during follow-up visits has become routine in our practice and has allowed some general observations. We noted early that research findings were frequently illuminating in a clinical context, and as a result we routinely inform patients of their DNA findings on follow-up visits, and cautiously interpret these. The advice many years ago from our ethics review board seems to have anticipated the current importance of disclosure of results to research patients . Since LipidSeq’s content is restricted to lipid disorders, there is no possibility of identifying secondary or incidental findings, except those related to other dyslipidemias.
In the course of reporting findings, we have observed in FH patients that: 1) knowing that there is a causative or relevant genetic finding seems to improve adherence to advice, particularly medication use; 2) when reporting polygenic effects, patients seem to understand the concept better when a simple tally of alleles is presented, compared to a weighted normalized statistic; 3) compliance seems unaffected by whether a monogenic or polygenic basis is communicated (we are undertaking a formal research project to address this issue); and 4) biochemical cascade screening is recommended regardless as to whether the hypercholesterolemia has a monogenic or polygenic basis, since multiple small-effect variants also tend to cluster in closely related family members.
Following from the positive experiences with LipidSeq, we have used it as the template to develop a similar panel for neurodegenerative conditions . Given our experiences, we hope that more individuals, both physicians and researchers, will be able to use similar sequencing strategies for their clinical work and research, respectively. We note that the overall design and procedure used for several commercial dyslipidemia panels offered since about 2016–17 resemble LipidSeq quite closely.
We have extensively applied our targeted sequencing panel for lipid disorders in a clinical context over several years. LipidSeq has enabled robust detection not just of rare variants underlying monogenic dyslipidemias, but also of CNVs due to high depth of coverage, and polygenic effects because of flexibility to detect common SNPs used in polygenic scores. This additional ability for assessing multiple genetic determinants across multiple genes simultaneously reduces genetic non-diagnoses that might otherwise result from overreliance on methods such as exome sequencing, which are optimized to uncover rare coding SNVs. In contrast, genome sequencing can potentially detect this wide range of variant types, but optimization of bioinformatic and ethical issues are needed first. Our accumulated observations, anecdotes and small case series suggest the value of genetic diagnosis for certain patients and clinical situations. But whether genetic diagnosis alters treatment decisions for the majority of dyslipidemic patients, above and beyond clinical and biochemical criteria alone, requires further study.
Availability of data and materials
Patient-level data is unavailable and cannot be shared due to patient privacy and our ethics form. Instead, annotated variant information for each identified rare variant from our total cohort (N = 3,262) is available in Additional file 1: Table S3 and Table S4.
Combined Annotation Dependent Depletion
Exome Aggregation Consortium
Genome Aggregation Database
Genome-wide association studies
Insertions or deletions
Maturity-onset diabetes of the young
Polymorphism Phenotyping version 2
Sorting Intolerant from Tolerant
Hegele RA. Plasma lipoproteins: genetic influences and clinical implications. Nat Rev Genet. 2009;10(2):109–21.
Hegele RA, Ginsberg HN, Chapman MJ, Nordestgaard BG, Kuivenhoven JA, Averna M, Boren J, Bruckert E, Catapano AL, Descamps OS, et al. The polygenic nature of hypertriglyceridaemia: implications for definition, diagnosis, and management. Lancet Diabetes Endocrinol. 2014;2(8):655–66.
Defesche JC, Gidding SS, Harada-Shiba M, Hegele RA, Santos RD, Wierzbicki AS. Familial hypercholesterolaemia. Nat Rev Dis Primers. 2017;3:17093.
Hegele RA, Ban MR, Cao H, McIntyre AD, Robinson JF, Wang J. Targeted next-generation sequencing in monogenic dyslipidemias. Curr Opin Lipidol. 2015;26(2):103–13.
Dron JS, Hegele RA. Genetics of lipid and lipoprotein disorders and traits. Curr Genet Med Rep. 2016;4(3):130–41.
Farhan SM, Hegele RA. Exome sequencing: new insights into lipoprotein disorders. Curr Cardiol Rep. 2014;16(7):507.
Dron JS, Hegele RA. Polygenic influences on dyslipidemias. Curr Opin Lipidol. 2018;29(2):133–43.
Johansen CT, Hegele RA. Genetic bases of hypertriglyceridemic phenotypes. Curr Opin Lipidol. 2011;22(4):247–53.
Brahm AJ, Hegele RA. Combined hyperlipidemia: familial but not (usually) monogenic. Curr Opin Lipidol. 2016;27(2):131–40.
Rahalkar AR, Hegele RA. Monogenic pediatric dyslipidemias: classification, genetics and clinical spectrum. Mol Genet Metab. 2008;93(3):282–94.
Fu J, Kwok S, Sinai L, Abdel-Razek O, Babula J, Chen D, Farago E, Fernandopulle N, Leith S, Loyzer M, et al. Western database of lipid variants (WDLV): a catalogue of genetic variants in monogenic dyslipidemias. Can J Cardiol. 2013;29(8):934–9.
Johansen CT, Dube JB, Loyzer MN, MacDonald A, Carter DE, McIntyre AD, Cao H, Wang J, Robinson JF, Hegele RA. LipidSeq: a next-generation clinical resequencing panel for monogenic dyslipidemias. J Lipid Res. 2014;55(4):765–72.
Hegele RA. Editorial: designing targeted sequencing panels for dyslipidemia. Curr Opin Lipidol. 2019;30(2):53–5.
Hegele RA, Brunt JH, Connelly PW. Multiple genetic determinants of variation of plasma lipoproteins in Alberta Hutterites. Arterioscler Thromb Vasc Biol. 1995;15(7):861–71.
Wang J, Ban MR, Zou GY, Cao H, Lin T, Kennedy BA, Anand S, Yusuf S, Huff MW, Pollex RL, et al. Polygenic determinants of severe hypertriglyceridemia. Hum Mol Genet. 2008;17(18):2894–9.
Johansen CT, Wang J, Lanktree MB, McIntyre AD, Ban MR, Martins RA, Kennedy BA, Hassell RG, Visser ME, Schwartz SM, 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.
Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466(7307):707–13.
Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45(11):1274–83.
Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, Schadt EE, Kaplan L, Bennett D, Li Y, Tanaka T, et al. Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet. 2009;41(1):56–65.
Dron JS, Hegele RA. The evolution of genetic-based risk scores for lipids and cardiovascular disease. Curr Opin Lipidol. 2019;30(2):71–81.
Dilliott AA, Farhan SMK, Ghani M, Sato C, Liang E, Zhang M, McIntyre AD, Cao H, Racacho L, Robinson JF, et al. Targeted next-generation sequencing and bioinformatics pipeline to evaluate genetic determinants of constitutional disease. J Vis Exp. 2018;134:57266.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.
Genomes Project C, Auton a, brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285–91.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. BioRxiv. 2019:531210.
Kircher M, Witten DM, Jain P, O'Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014;46(3):310–5.
Rentzsch P, Witten D, Cooper GM, Shendure J, Kircher M. CADD: predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019;47(D1):D886–94.
Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013, Chapter 7:Unit7:20.
Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4.
Schwarz JM, Rodelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010;7(8):575–6.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.
Iacocca MA, Wang J, Dron JS, Robinson JF, McIntyre AD, Cao H, Hegele RA. Use of next-generation sequencing to detect LDLR gene copy number variation in familial hypercholesterolemia. J Lipid Res. 2017;58(11):2202–9.
Brahm AJ, Hegele RA. Chylomicronaemia--current diagnosis and future therapies. Nat Rev Endocrinol. 2015;11(6):352–62.
Dron JS, Ho R, Hegele RA. Recent advances in the genetics of Atherothrombotic disease and its determinants. Arterioscler Thromb Vasc Biol. 2017;37(10):e158–66.
Di Angelantonio E, Gao P, Pennells L, Kaptoge S, Caslake M, Thompson A, Butterworth AS, Sarwar N, Wormser D, Saleheen D, et al. Lipid-related markers and cardiovascular disease prediction. JAMA. 2012;307(23):2499–506.
Perk J, De Backer G, Gohlke H, Graham I, Reiner Z, Verschuren M, Albus C, Benlian P, Boysen G, Cifkova R, et al. European Guidelines on cardiovascular disease prevention in clinical practice (version 2012). The fifth joint task force of the European Society of Cardiology and Other Societies on cardiovascular disease prevention in clinical practice (constituted by representatives of nine societies and by invited experts). Eur Heart J. 2012;33(13):1635–701.
Parish S, Offer A, Clarke R, Hopewell JC, Hill MR, Otvos JD, Armitage J, Collins R. Lipids and lipoproteins and risk of different vascular events in the MRC/BHF heart protection study. Circulation. 2012;125(20):2469–78.
Zanoni P, Khetarpal SA, Larach DB, Hancock-Cerutti WF, Millar JS, Cuchel M, DerOhannessian S, Kontush A, Surendran P, Saleheen D, et al. Rare variant in scavenger receptor BI raises HDL cholesterol and increases risk of coronary heart disease. Science. 2016;351(6278):1166–71.
Dron JS, Wang J, Berberich AJ, Iacocca MA, Cao H, Yang P, Knoll J, Tremblay K, Brisson D, Netzer C, et al. Large-scale deletions of the ABCA1 gene in patients with hypoalphalipoproteinemia. J Lipid Res. 2018;59(8):1529–35.
Berberich AJ, Huot C, Cao H, AD MI, Robinson JF, Wang J, Hegele RA. Copy Number Variation in GCK in patients with Maturity-Onset Diabetes of the Young. J Clin Endocrinol Metab. 2019;104:3428–36.
Berberich AJ, Mokashi A, McIntyre AD, Robinson JF, Cao H, Wang J, Hegele RA. Bioinformatic detection of copy number variation in HNF4A causing maturity onset diabetes of the young. Clin Genet. 2019;96(4):376–7.
Dron JS, Wang J, McIntyre AD, Cao H, Robinson JF, Duell B, Manjoo P, Feng J, Movsesyan I, Malloy MJ, et al. Partial LPL deletions: rare copy-number variants contributing towards the polygenic form of severe hypertriglyceridemia. J Lipid Res. 2019;60:1953–8.
Iacocca MA, Wang J, Sarkar S, Dron JS, Lagace T, McIntyre AD, Lau P, Robinson JF, Yang P, Knoll JH, et al. Whole-gene duplication of PCSK9 as a novel genetic mechanism for severe familial hypercholesterolemia. Can J Cardiol. 2018;34(10):1316–24.
Wang J, Dron JS, Ban MR, Robinson JF, McIntyre AD, Alazzam M, Zhao PJ, Dilliott AA, Cao H, Huff MW, et al. Polygenic versus monogenic causes of hypercholesterolemia ascertained clinically. Arterioscler Thromb Vasc Biol. 2016;36(12):2439–45.
Dron JS, Wang J, Cao H, McIntyre AD, Iacocca MA, Menard JR, Movsesyan I, Malloy MJ, Pullinger CR, Kane JP, et al. Severe hypertriglyceridemia is primarily polygenic. J Clin Lipidol. 2019;13(1):80–8.
Corral P, Geller AS, Polisecki EY, Lopez GI, Banares VG, Cacciagiu L, Berg G, Hegele RA, Schaefer EJ, Schreier LE. Unusual genetic variants associated with hypercholesterolemia in Argentina. Atherosclerosis. 2018;277:256–61.
Iacocca MA, Chora JR, Carrie A, Freiberger T, Leigh SE, Defesche JC, Kurtz CL, DiStefano MT, Santos RD, Humphries SE, et al. ClinVar database of global familial hypercholesterolemia-associated DNA variants. Hum Mutat. 2018;39(11):1631–40.
Brinton EA, Hopkins PN, Hegele RA, Geller AS, Polisecki EY, Diffenderfer MR, Schaefer EJ. The association between hypercholesterolemia and sitosterolemia, and report of a sitosterolemia kindred. J Clin Lipidol. 2018;12(1):152–61.
Lee T, Iacocca MA, Ban MR, Hegele RA. Efficacy of evolocumab in monogenic versus polygenic hypercholesterolemia. CJC Open. 2019:115–8.
Hegele RA, Berberich AJ, Ban MR, Wang J, Digenio A, Alexander VJ, D'Erasmo L, Arca M, Jones A, Bruckert E, et al. Clinical and biochemical features of different molecular etiologies of familial chylomicronemia. J Clin Lipidol. 2018;12(4):920–7 e924.
Paquette M, Bernard S, Hegele RA, Baass A. Chylomicronemia: differences between familial chylomicronemia syndrome and multifactorial chylomicronemia. Atherosclerosis. 2019;283:137–42.
Eguchi J, Miyashita K, Fukamachi I, Nakajima K, Murakami M, Kawahara Y, Yamashita T, Ohta Y, Abe K, Nakatsuka A, et al. GPIHBP1 autoantibody syndrome during interferon beta1a treatment. J Clin Lipidol. 2019;13(1):62–9.
Chyzhyk V, Kozmic S, Brown AS, Hudgins LC, Starc TJ, Davila AD, Blevins TC, Diffenderfer MR, He L, Geller AS, et al. Extreme hypertriglyceridemia: genetic diversity, pancreatitis, pregnancy, and prevalence. J Clin Lipidol. 2019;13(1):89–99.
Berberich AJ, Ziada A, Zou GY, Hegele RA. Conservative management in hypertriglyceridemia-associated pancreatitis. J Intern Med. 2019;286:644–50.
Dron JS, Wang J, Low-Kam C, Khetarpal SA, Robinson JF, McIntyre AD, Ban MR, Cao H, Rhainds D, Dube MP, et al. Polygenic determinants in extremes of high-density lipoprotein cholesterol. J Lipid Res. 2017;58(11):2162–70.
Low-Kam C, Rhainds D, Lo KS, Barhdadi A, Boule M, Alem S, Pedneault-Gagnon V, Rheaume E, Dube MP, Busseuil D, et al. Variants at the APOE /C1/C2/C4 locus modulate cholesterol efflux capacity independently of high-density lipoprotein cholesterol. J Am Heart Assoc. 2018;7(16):e009545.
Geller AS, Polisecki EY, Diffenderfer MR, Asztalos BF, Karathanasis SK, Hegele RA, Schaefer EJ. Genetic and secondary causes of severe HDL deficiency and cardiovascular disease. J Lipid Res. 2018;59(12):2421–35.
Paquette M, Dufour R, Hegele RA, Baass A. A tale of 2 cousins: an atypical and a typical case of abetalipoproteinemia. J Clin Lipidol. 2016;10(4):1030–4.
Chan D, McIntyre AD, Hegele RA, Don-Wauchope AC. Familial partial lipodystrophy presenting as metabolic syndrome. J Clin Lipidol. 2016;10(6):1488–91.
Xiao C, Hegele RA, Lewis GF. Treatment of severe hypercholesterolemia in a woman with advanced primary Sclerosing cholangitis. JAMA Cardiol. 2017;2(5):575–6.
Jackson C, Geraghty MT, Hegele RA, El Demellawy D, Jimenez-Rivera C. Fatty liver in a non-obese patient. Paediatr Child Health. 2017;22(2):59–60.
Varma S, McIntyre AD, Hegele RA. Atypical presentation and treatment response in a child with familial hypercholesterolemia having a novel LDLR mutation. JIMD Rep. 2017;35:67–70.
Soto AG, McIntyre A, Agrawal S, Bialo SR, Hegele RA, Boney CM. Severe Hypertriglyceridemia due to a novel p.Q240H mutation in the Lipoprotein Lipase gene. Lipids Health Dis. 2015;14:102.
Rosenson RS, Najera SD, Hegele RA. Heterozygous familial hypercholesterolemia presenting as chylomicronemia syndrome. J Clin Lipidol. 2017;11(1):294–6.
Aljenedil S, Hegele RA, Genest J, Awan Z. Estrogen-associated severe hypertriglyceridemia with pancreatitis. J Clin Lipidol. 2017;11(1):297–300.
Stahel P, Xiao C, Hegele RA, Lewis GF. Polygenic risk for hypertriglyceridemia can mimic a major monogenic mutation. Ann Intern Med. 2017;167(5):360–1.
Ueda M, Dunbar RL, Wolska A, Sikora TU, Escobar MDR, Seliktar N, deGoma E, DerOhannessian S, Morrell L, AD MI, et al. A novel APOC2 missense mutation causing Apolipoprotein C-II deficiency with severe Triglyceridemia and pancreatitis. J Clin Endocrinol Metab. 2017;102(5):1454–7.
Sundaram M, Curtis KR, Amir Alipour M, LeBlond ND, Margison KD, Yaworski RA, Parks RJ, McIntyre AD, Hegele RA, Fullerton MD, et al. The apolipoprotein C-III (Gln38Lys) variant associated with human hypertriglyceridemia is a gain-of-function mutation. J Lipid Res. 2017;58(11):2188–96.
Paquette M, Hegele RA, Pare G, Baass A. A novel mutation in GPIHBP1 causes familial chylomicronemia syndrome. J Clin Lipidol. 2018;12(2):506–10.
Morise AP, Hegele RA. Atypical familial dysbetalipoproteinemia associated with high polygenic cholesterol and triglyceride scores treated with ezetimibe and evolocumab. J Clin Lipidol. 2019.
Wojcik C, Fazio S, McIntyre AD, Hegele RA. Co-occurrence of heterozygous CREB3L3 and APOA5 nonsense variants and polygenic risk in a patient with severe hypertriglyceridemia exacerbated by estrogen administration. J Clin Lipidol. 2018;12(5):1146–50.
Wang LR, McIntyre AD, Hegele RA. Complex genetic architecture in severe hypobetalipoproteinemia. Lipids Health Dis. 2018;17(1):48.
Hussain I, Patni N, Ueda M, Sorkina E, Valerio CM, Cochran E, Brown RJ, Peeden J, Tikhonovich Y, Tiulpakov A, et al. A novel generalized Lipodystrophy-associated Progeroid syndrome due to recurrent heterozygous LMNA p.T10I mutation. J Clin Endocrinol Metab. 2018;103(3):1005–14.
Perera U, Kennedy BA, Hegele RA. Multiple Symmetric Lipomatosis (Madelung Disease) in a Large Canadian Family With the Mitochondrial MTTK c.8344A>G Variant. J Investig Med High Impact Case Rep. 2018;6:2324709618802867.
Wang LR, Radonjic A, Dilliott AA, McIntyre AD, Hegele RA. A De novo POLD1 mutation associated with mandibular hypoplasia, deafness, Progeroid features, and Lipodystrophy syndrome in a family with Werner syndrome. J Investig Med High Impact Case Rep. 2018;6:2324709618786770.
Brahm AJ, Wang G, Wang J, McIntyre AD, Cao H, Ban MR, Hegele RA. Genetic confirmation rate in clinically suspected maturity-onset diabetes of the young. Can J Diabetes. 2016;40(6):555–60.
Iacocca MA, Hegele RA. Recent advances in genetic testing for familial hypercholesterolemia. Expert Rev Mol Diagn. 2017;17(7):641–51.
Dron JS, Hegele RA. Genetics of triglycerides and the risk of atherosclerosis. Curr Atheroscler Rep. 2017;19(7):31.
Iacocca MA, Hegele RA. Role of DNA copy number variation in dyslipidemias. Curr Opin Lipidol. 2018;29(2):125–32.
Iacocca MA, Dron JS, Hegele RA. Progress in finding pathogenic DNA copy number variations in dyslipidemia. Curr Opin Lipidol. 2019;30(2):63–70.
Berberich AJ, Hegele RA. The role of genetic testing in dyslipidaemia. Pathology. 2019;51(2):184–92.
Farhan SMK, Dilliott AA, Ghani M, Sato C, Liang E, Zhang M, McIntyre AD, Cao H, Racacho L, Robinson JF, et al. The ONDRISeq panel: custom-designed next-generation sequencing of genes related to neurodegeneration. NPJ Genom Med. 2016;1:16032.
Kerkhof J, Schenkel LC, Reilly J, McRobbie S, Aref-Eshghi E, Stuart A, Rupar CA, Adams P, Hegele RA, Lin H, et al. Clinical validation of copy number variant detection from targeted next-generation sequencing panels. J Mol Diagn. 2017;19(6):905–20.
Chora JR, Alves AC, Medeiros AM, Mariano C, Lobarinhas G, Guerra A, Mansilha H, Cortez-Pinto H, Bourbon M. Lysosomal acid lipase deficiency: a hidden disease among cohorts of familial hypercholesterolemia? J Clin Lipidol. 2017;11(2):477–84 e472.
Khera AV, Won HH, Peloso GM, Lawson KS, Bartz TM, Deng X, van Leeuwen EM, Natarajan P, Emdin CA, Bick AG, et al. Diagnostic yield and clinical utility of sequencing familial hypercholesterolemia genes in patients with severe hypercholesterolemia. J Am Coll Cardiol. 2016;67(22):2578–89.
Bombard Y, Brothers KB, Fitzgerald-Butt S, Garrison NA, Jamal L, James CA, Jarvik GP, McCormick JB, Nelson TN, Ormond KE, et al. The responsibility to Recontact research participants after reinterpretation of genetic and genomic research results. Am J Hum Genet. 2019;104(4):578–95.
We would like to acknowledge and thank the patients involved in this study. We thank Ericka Simon, Brooke Kennedy, Lauren Thomson and Lydia Vermeer for valuable contributions to clinical databases. We also thank Shehan Perera, Laney Beaulieu, and Ziyad Khatab for their efforts in classifying variants using the ACMG guidelines.
J.S.D. is supported by the Canadian Institutes of Health Research (Doctoral Research Award). R.A.H. is supported by the Jacob J. Wolfe Distinguished Medical Research Chair, the Edith Schulich Vinet Research Chair in Human Genetics, and the Martha G. Blackburn Chair in Cardiovascular Research. R.A.H. has also received operating grants from the Canadian Institutes of Health Research (Foundation award), and the Heart and Stroke Foundation of Ontario (G-15-0009214). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Patients provided informed, signed consent with ethics approval from the Western University ethics review board (#07290E).
Consent for publication
R.A.H. reports consulting fees from Acasti, Aegerion, Akcea/Ionis, Amgen and Sanofi. The other authors have no disclosures.
This additional file contains four supplemental tables, each in its own labelled tab. Table S1. Contains a description of each gene on the LipidSeq panel and what disorder or trait it is related to, as well its chromosmal coordinates. Table S2. Contains information on our phenotypic categorization of samples that have been sequenced using the LipidSeq panel design. Table S3. Contains annotated information on every unique variant we have identified from our entire cohort. Table S4. Contains a subset of annotated variants with CADD PHRED scores greater than or equal to 20.
About this article
Cite this article
Dron, J.S., Wang, J., McIntyre, A.D. et al. Six years’ experience with LipidSeq: clinical and research learnings from a hybrid, targeted sequencing panel for dyslipidemias. BMC Med Genomics 13, 23 (2020). https://doi.org/10.1186/s12920-020-0669-2
- Targeted next-generation sequencing panel
- Familial hypercholesterolemia
- Metabolic disorder