Hepatic lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels

Abstract

Common variants in the hepatic lipase (LIPC) gene have been shown to be associated with plasma lipid levels; however, the distribution and functional features of rare and regulatory LIPC variants contributing to the extreme lipid phenotypes are not well known. This study was aimed to catalogue LIPC variants by resequencing the entire LIPC gene in 95 non-Hispanic Whites (NHWs) and 95 African blacks (ABs) with extreme HDL-C levels followed by in silico functional analyses. A total of 412 variants, including 43 novel variants were identified; 56 were unique to NHWs and 234 were unique to ABs. Seventy-eight variants in NHWs and 89 variants in ABs were present either in high HDL-C group or low HDL-C group. Two non-synonymous variants (p.S289F, p.T405M), found in NHWs with high HDL-C group were predicted to have damaging effect on LIPC protein by SIFT, MT2 and PP2. We also found several non-coding variants that possibly reside in the circRNA and lncRNA binding sites and may have regulatory potential, as identified in rSNPbase and RegulomeDB databases. Our results shed light on the regulatory nature of rare and non-coding LIPC variants as well as suggest their important contributions in affecting the extreme HDL-C phenotypes.

About the Authors:

Dilek Pirim

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

Affiliations Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Department of Molecular Biology and Genetics, Faculty of Arts & Science, Bursa Uludag University, Gorukle, Bursa, Turkey

ORCID logo https://orcid.org/0000-0002-0522-9432

Clareann H. Bunker

Roles Resources, Writing – review & editing

Affiliation: Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

John E. Hokanson

Roles Resources, Writing – review & editing

Affiliation: Department of Epidemiology, Colorado School of Public Health, University of Colorado Denver, Aurora, Colorado, United States of America

Richard F. Hamman

Roles Resources, Writing – review & editing

Affiliation: Department of Epidemiology, Colorado School of Public Health, University of Colorado Denver, Aurora, Colorado, United States of America

F. Yesim Demirci

Roles Investigation, Methodology, Project administration, Supervision, Writing – review & editing

Affiliation: Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

M. Ilyas Kamboh

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

Introduction

The hepatic lipase (HL) gene, also known as LIPC, is one of the member of the lipase gene family that functions in the uptake of high density lipoprotein cholesterol (HDL-C) and plays an essential multifunctional role in lipoprotein and lipid metabolism [1–3]. The LIPC gene is located on chromosome 15q21-q23 spanning approximately 138kb and consists of a total 9 exons coding 499 amino acid (~55kDa) glycoprotein HL or LIPC [4–6]. Due to its critical role in lipid metabolism, the LIPC activity dramatically causes fluctuating in plasma lipoprotein-lipid levels [7–9]. Previous studies have shown genetic variants in the LIPC affect the activity of LIPC and contribute to the risk of several diseases, including coronary artery disease (CAD), type 2 diabetes, metabolic syndrome and HL deficiency. Genome-wide association studies (GWASs) have revealed several common [minor allele frequency (MAF>5%)] LIPC variants that were associated with plasma levels of HDL-C, triglycerides (TG), and total cholesterol (TC) [10–19]. In addition to common variants, rare and low frequency variants (MAF≤5%) also impact phenotypic variation in plasma lipid levels [20]. Thus, it is imperative to identify and elucidate the roles of rare variants in lipid genes to explain the missing heritability of lipid phenotypes. Indeed, resequencing individuals with extreme lipid profiles has shown to confer advantage for detecting rare variants contributing to plasma lipid variation [21–27], which warrant further investigation in diverse population.

In this study, we resequenced the LIPC gene in 190 subjects with extreme HDL-C levels in order to identify the distribution of LIPC sequence variants in extreme phenotypes as well as to assess their functional relevance.

Materials and methods

Subjects

The study sample comprised of 190 samples, including 95 African Blacks (ABs) and 95 non-Hispanic Whites (NHWs) with extreme HDL-C levels. These small subset of samples were selected from the upper and lower 10th percentile of HDL-C distribution from two epidemiological well-characterized samples including 623 NHWs and 788 ABs for variant discovery purposes [28–31]. The biometric and quantitative data of 95 NHW and 95 AB samples with extreme HDL-C groups are presented in Table 1. The research was conducted in accordance with the relevant ethical guidelines/regulations and approved by the Institutional Ethical Review Boards of the University of Pittsburgh and University of Colorado Denver. Written informed consent was obtained from all participants.

[Figure omitted. See PDF.]

Table 1. Biometric and quantitative data (mean± SD) of the 95 NHW and 95 African black samples with extreme* HDL-C levels.

https://doi.org/10.1371/journal.pone.0243919.t001

Lipid measurements

Friedewald equation was used to calculate the plasma low-density lipoprotein cholesterol (LDL-C) levels and esterase-oxidase method was employed for measuring fasting total cholesterol levels [32, 33]. Enzymatic procedures described in Harris et al. [34] were used to determine serum HDL-C and triglyceride concentrations.

DNA sequencing

The entire LIPC gene (NG_011465, NM_000236) (excluding the large intron 1) (see S1 Fig), plus 1kb in 3' and 5' flanking sequence comprising ~ 33kb was sequenced in both direction using a total of 46 resequencing amplicons. DNA isolations were conducted using buffy coat and blood clots from of NHW and AB samples, respectively, following the standard DNA isolation protocols. PCR primers were prepared by using the Primer 3 software (http://frodo.wi.mit.edu/primer3/) (see S1 Table). PCR conditions are available upon request. Sanger sequencing was performed on the ABI 3730x1DNA analyzers in a commercial lab (Beckman Coulter Genomics, Danvers, MA). Variant Reporter (Applied Biosystems, Foster City, CA) and Sequencher (Gene Codes, Ann Arbor, MI) softwares were used for analyzing the sequence chromatograms and variant detection.

Statistical and bioinformatic analysis

Haploview software (www.broadinstitute.org/haploview) was used to calculate allele frequencies. Departure from Hardy-Weinberg equilibrium and difference of allele frequencies between high and low HDL groups were determined using a chi-squared test [35]. A p-value of <0.05 was considered as "suggestive evidence". G power software (version 3.1) was used for power analysis, which showed that our sample size has sufficient power (97% at 5% alpha level) to detect low frequency variants (MAF<0.05) with moderate effect sizes (Cohen’s d = 0.05). TagSNPs [(MAF)≥0.05, r2≥0.8)] and linkage disequilibrium (LD) analyses were also performed for common variants [minor allele frequency (MAF) ≥5%] in Haploview.

Functional effects of the coding variants were predicted using common in silico tools [Sorting Intolerant From Tolerant (SIFT) [36], Mutation Taster 2 (MT2) [37] and Polymorphism Phenotyping v.2 (PP2) [38]. Functional annotations of variants located in non-coding regions were determined using RegulomeDB database (http://regulome.stanford.edu/) and rSNPBase 3.1 database (http://rsnp3.psych.ac.cn) [39, 40]. rSNPBase 3.1 helps to identify SNP-related regulatory elements [transcription factor binding regions (TFBRs), chromatin interactive regions (CIRs), miRNA target sites, lncRNA regions, TADs and circRNAs (circular RNAs)] and their target regulatory genes incorporating mostly experimental data from the Encyclopaedia of DNA Elements (ENCODE) project and other sources such as LNCipedia, CircNet, miRBase, miRNAda and TargetScan. It also provides information related eQTL (expression Quantitative Trait Loci) and disease associations for each query variant. Data annotations for epigenetic marks including active chromatin state regions, histone binding regions, and methylation sites are not available in rSNPBase 3.1, thus we also used RegulomeDB database which uses a scoring scheme to annotate variants based on their regulatory impact related epigenetic marks as well as other regulatory information.

Results

DNA sequencing

We identified a total of 412 variants by resequencing the LIPC gene in 190 individuals from two distinct population; of which 122 [4 indel, 1 triallelic single nucleotide variant (SNV) (rs7171818) and 65 diallelic SNVs] were shared in both populations, 56 (4 indels, 152 SNVs) were unique to NHWs and 234 (9 indels, 225 SNVs) were unique to ABs (S2 and S3 Tables). A total of 43 novel variants (not previously reported in any public databases) were observed; of which 26 were present only in ABs and 14 were found only in NHWs. Novel variants were submitted to dbSNP database using handle ID: KAMBOH and distinct dbSNP IDs were assigned for each novel variant (http://www.ncbi.nlm.nih.gov/SNP/snp_viewTable.cgi?handle=KAMBOH). While majority of the identified variants were located in introns, 17 were present in the coding regions of the LIPC gene (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Locations of the LIPC variants identified in NHWs (n = 95) and ABs (n = 95).

https://doi.org/10.1371/journal.pone.0243919.g001

Of the 178 variant identified in NHWs, 86 were common (MAF≥0.05), 21 were uncommon (0.01≤MAF<0.05) and 71 were rare (MAF<0.01). Among the 356 variants found in ABs, 172 were common (MAF≥0.05), 116 were uncommon (0.01≤MAF<0.05) and 68 were rare (MAF<0.01) (Fig 2). MAFs of all 178 variants in the total NHW samples and between the two extreme HDL-C groups are shown in S2 Table. Likewise, the same information for the 356 variants in ABs, is shown in S3 Table.

[Figure omitted. See PDF.]

Fig 2. Distributions of the minor allele frequencies of LIPC variants in NHWs (n = 95) and ABs (n = 95).

https://doi.org/10.1371/journal.pone.0243919.g002

We identified 17 coding SNVs in both populations; of which 11 were synonymous and 6 were non-synonymous. The 6 non-synonymous variants were: valine to methionine (V95M) in exon 3, asparagine to serine (N215S) in exon 5, serine to phenylalanine (S289F) in exon 6, valine to isoleucine (V342I) in exon 6, phenylalanine to leucine (F356L) in exon 7, threonine to methionine (T405M) in exon 8.

TagSNP identification and linkage disequilibrium (LD) analyses

By using Haploview software, we identified tagSNPs and determined LD between common variants (MAF≥0.05, r2≥0.80) in NHWs and ABs. Of the 174 common variants in ABs (including 3 alleles of the triallelic variant rs7171818, See S3), 3 deviated from Hardy-Weinberg equilibrium (HW-P-value<0.0001) and 10 had low call rate (<80%). Of the 88 common variants in NHWs (including 3 alleles of the triallelic variant rs7171818, See S2), 9 had low call rate (<80%) and 1 deviated from Hardy-Weinberg equilibrium (HW-P-value<0.0003). Accordingly, LD structures for 78 variants in NHWs and 161 variants in ABs were analyzed and tagSNPs were identified by Tagger (see S4 and S5 Tables). Among 78 common variants in NHWs and 161 common variants in ABs, 58 biallelic variants were shared between the two populations. Fig 3 displays the pairwise LD structure of common LIPC variants which were highly different in NHWs and ABs.

[Figure omitted. See PDF.]

Fig 3. LD structure of the LIPC variants identified in both populations.

(a) NHWs, (b) African blacks. Shade intensity indicates the degree of LD (r2 between 0 and 1). Black indicates complete LD (r2 = 1), white indicates no LD (r2 = 0).

https://doi.org/10.1371/journal.pone.0243919.g003

Distribution of LIPC variants in NHWs and ABs with extreme HDL levels

We analyzed the MAF distributions of quality control (QC) passed 78 common variants in NHWs and 161 variants in ABs between the two extreme HDL-C groups (see S2 and S3 Tables). Two variants [rs6082 (P = 0.0136), rs12592139 (P = 0.046)] in ABs (see S3 Table) and one variant [rs143731122 (P = 0.028)] in NHWs (see S2 Table) show suggestive evidence for different MAF distributions between extreme HDL-C groups. MAFs of these three variants were lower in the low HDL-C group compared to the high HDL-C group. There were also some variants in both populations that were found either in individuals with the low HDL-C group or in the high HDL-C group. In NHWs, 78 variants (0.005≤MAF≤0.022) were present in one or the other group, including 53 only in low HDL-C group and 25 in only high HDL-C group (S2 Table). In ABs, a total of 89 variants (0.005≤MAF≤0.039) were observed in one or the other group, including 47 in low HDL-C group and 42 in high HDL-C group (S3 Table).

Functional annotations of LIPC variants

Of the 17 LIPC coding variants identified, 6 were associated with amino acid changes and 11 were synonymous. SIFT and PP2 tools only predict the effects of non-synonymous variants, however, MT2 evaluates the possible effects of both synonymous and non-synonymous variants. Two rare (MAF = 0.005) LIPC variants [rs121912502 (p.Ser289Phe), rs113298164 (p.Thr405Met)], which were unique to NHW population and seen in only individuals with high HDL-C levels were found to be damaging, probably damaging and disease causing in SIFT, PP2 and MT2, respectively (Table 2). The locations of the missense variants on protein structure of the LIPC were shown in S2–S7 Figs and MAFs of all identified coding LIPC variants reported in gnomAD were listed in Table 2.

[Figure omitted. See PDF.]

Table 2. LIPC coding variants identified in 95 NHWs and 95 African blacks with extreme HDL-C levels.

https://doi.org/10.1371/journal.pone.0243919.t002

Moreover, 4 synonymous LIPC variants that were observed in only AB individuals were predicted as disease causing by MT2 computational predictions. Among them, one [rs113174258 (p.Thr71Thr) was a rare variant (MAF = 0.005) and observed in only one individual with low HDL-C level (29.3mg/dl).

Regulatory relevance of the intronic variants using rSNPbase and RegulomeDB databases revealed several regulatory variants (rSNP) that have potential to affect the binding of regulatory elements. A total of 37 and 46 rSNPs were identified in NHWs and ABs, respectively, that have putative binding regions for circRNA or lncRNA. Among the 37 rSNP identified in NHWss, 23 were unique to one extreme HDL-C group and they were found to be located in circRNA and lncRNA binding regions. Of the 46 rSNPs identified in ABs, 9 were seen in only one extreme HDL-C group with MAF≤0.016 and all were determined to be located in circRNA regions and had RegulomeDB score≥3a indicating their regulatory impact (Table 3).

[Figure omitted. See PDF.]

Table 3. Distributions of the regulatory variants (rSNPs) identified in 190 individuals (95 NHWs and 95 ABs) with extreme HDL-C levels.

https://doi.org/10.1371/journal.pone.0243919.t003

We also assessed the regulatory features of all LIPC variants by submitting all refSNP IDs to the RegulomeDB database (v 2.0) and variants with available data were assigned a score ranging from 1 to 6. S2 and S3 Tables list the RegulomeDB score for all variants identified in NHWs and ABs, respectively. Lower score 1 indicates higher evidence that the variant may affect the binding region of regulatory elements and have impact on the expression of the target of the gene, whereas score 2 indicates strong evidence that the variant may reside in the binding site of the regulatory elements and proteins. RegulomeDB analyses revealed 36 variants in ABs and 20 variants in NHWs with strong evidence (RegulomeDB score<3) indicating their regulatory roles, and some were also identified as rSNP based on rSNPbase results (S2 and S3 Tables).

Discussion

Accumulating evidence supports the influence of genetic variants on plasma concentration of HDL-C levels. Majority of the reported genetic variants associated with plasma HDL-C levels reside in genes coding for key enzymes in lipid metabolism. LIPC is one of the key lipid genes where its common sequence variants have been reported to be associated with variation in lipoprotein-lipid levels and CAD risk [41–45]. Moreover, individuals with rare loss of functional variants in the LIPC gene were observed to have deficient hepatic lipase, resulting in lipid disturbance and affecting CHD risk [46].

In order to further understand the role of common and rare LIPC variants in affecting plasma HDL-C levels, we resequenced the LIPC gene in individuals with extreme HDL-C levels in two ethnic groups and identified 412 variants, including 290 population-specific and several rare variants with MAF<0.01 (Fig 2), of which some were only present in one of the two extreme HDL-C groups (see S2 and S3 Tables). Our LD analyses showed population-specific LD structure of common LIPC variants where there was high LD between variants in NHWs compared to ABs (Fig 3).

To the best of our knowledge, this is the first study that provides a catalog of common and rare LIPC variants by resequencing the entire gene in two well-characterized population based samples with extreme HDL-C levels. In a recent study, only exons and exon-intron junctions of LIPC were resequenced in Koreans with extreme HDL-C levels (n = 42) where no rare variants were identified [27]. Three non-synonymous variants (rs3829462, rs6078, rs6083) were identified in Korean subjects with extreme HDL-C levels of which rs3829462 and rs6078 were seen in all studied subjects (n = 42) and rs6083 were found in only 16 individuals. In our study, the rs3829462 SNP was observed only in the AB sample and its distribution along with rs6078 and rs6083 were similar between the two HDL-C groups in ABs. In NHWs, the distribution of rs6083 was also similar between the two extreme groups, and rs6078 was observed in only one individual with high HDL-C (see S3 Table).

Moreover, we found three common variants [rs6082 (P = 0.0136), rs12592139 (P = 0.0463), rs143731122 (P = 0.028)] that showed evidence of association with extreme HDL-C levels (see S3 and S4 Tables). The rs6082 (p.G197G) is a synonymous coding variant that was previously identified in individuals with HL deficiency and is listed as a benign variant in the ClinVar database. The most extensively studied common LIPC variants are located in the promoter, -250G/A (rs2070895) and -514C/T (rs1800588), that have been shown to be associated with HL activity, HDL-C and with the risk of metabolic diseases [19, 41, 45, 47–49]. The location of -514C/T (rs1800588) was out of our sequenced region range and thus we did not have data for this variant On the other hand, -250G/A (rs2070895) was observed in both populations and the frequency of the A allele was higher in individuals with high HDL-C group than in low HDL-C group in both NHWs (20.7% vs 12.8%) and ABs (65.6% vs 54.3%). Our results are in good agreement with a recent GWAS finding where the rs2070895-A was associated with increased levels of HDL-C levels (p = 4 x 10−24) [50].

We also observed 78 uncommon or rare variants (MAF<0.05) in NHWs and 89 in ABs that were present in only one or the other extreme HDL-C group. Of the variants that were identified in only one extreme group, seven were coding including four [rs6078 (p.V95M), rs121912502 (p.S289F), rs145811475 (p.V342I), rs113298164 (p.T405M) causing amino acid change (Table 2). The rs121912502 (p.Ser289Phe) variant was observed in only one NHW individual in the high HDL-C group and this was reported to be of uncertain significance associated with HL deficiency in ClinVar database. Another coding variant, rs113298164 (p.Thr405Met), which was reported as pathogenic and associated with HL deficiency in ClinVar, was also found in only one NHW individual with high HDL-C levels. Our finding is in accordance with a recent study which suggested the association of the minor T allele of rs113298164 with high HDL-C levels [51]. Our results indicate that the sequencing subset of samples having extreme HDL-C levels selected from the larger population-based samples have sufficient power to detect variants with strong effect size. Notably, our sequencing study design on extreme lipid phenotypes using the same sample size has enabled us to detect multiple rare and novel variants in addition to common variants in other candidate lipid genes [52–58].

Although the effects of coding LIPC variants on lipid traits have been widely investigated, the contribution of LIPC non-coding regulatory variants in interindividual variation in plasma lipids is still unclear. Thus, we also assessed the regulatory nature of non-coding LIPC variants by using two databases (RegulomeDB and rSNPbase) to prioritize the non-coding variants as regulatory variants. None of the identified 43 novel non-coding variants were found to be located in regions that have regulatory significance. However, 42 LIPC variants, seen in one or the other extreme HDL-C group, were implicated to reside in binding regions of regulatory elements, including one exonic variant [rs113298164 (p.Thr405Met)]. Our analyses suggest that non-coding regulatory LIPC variants have the potential to disrupt the regulatory functions of circRNA and lncRNA as well as highlight the possible regulatory role of exonic rs113298164 variant in splicing dysregulation. The impact of exonic variants in splicing mechanisms and their contribution to phenotypes highlight the imperative role of investigating the roles of exonic variants in splicing dysregulation [59]. Interestingly, majority of the regulatory LIPC variants were seen in individuals with the low HDL-C in both populations (Table 3).

Our study has some limitations. Our sequencing did not cover the large intron 1 due to technical reasons and thus we missed the potential rare and common LIPC sequence variants in this region. Our results depend on the data of 190 chromosomes with extreme phenotypes in each ethnic group and our sample size is small to claim associations between variants and HDL-C levels. Thus variants with suggestive evidence of association with HDL-C levels should be tested in a larger sample. Nevertheless, our effort produced a large catalog of LIPC sequence variants in samples from two ethnic groups and yielded several novel variants as well as rare population-specific variants located in the LIPC regulatory regions that could be followed up in larger future studies. Also, our findings suggest that individuals with extreme HDL-C levels carrying variants that might affect the activity of LIPC should be followed up for their risk cardiovascular disease.

Conclusions

In conclusion, our study reaffirms the considerable contribution of LIPC variants to plasma concentrations of HDL-C levels in the general population and also highlights the high prevalence of the rare variants that may play key role in the regulation of the HDL-C levels. Our results also emphasize the possible regulatory impact of non-coding LIPC variants in the determination of HDL-C levels. However, associations between suggested regulatory LIPC variants and lipid traits should be tested in large samples and functional studies need to be conducted to further evaluate the roles of non-coding regulatory variants in the lipid metabolism.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Exon and intron structure of the LIPC gene.

The image was retrieved from http://www.ensembl.org/. Boxes and lines between boxes indicate exons and introns, respectively. Unfilled box indicates untranslated region.

https://doi.org/10.1371/journal.pone.0243919.s001

(TIF)

S2 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs6078 (p.V95M) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s002

(TIF)

S3 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs6083 (p.N215S) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s003

(TIF)

S4 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs121912502 (p.S289F) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s004

(TIF)

S5 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs145811475 (p. V342I) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s005

(TIF)

S6 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs145811475 (p.F356L) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s006

(TIF)

S7 Fig. Diagram depicts a part of the hepatic lipase protein structure where rs145811475 (p. T405M) is located.

The image was retrieved from the VarSite database (https://www.ebi.ac.uk/thornton-srv/databases/VarSite).

https://doi.org/10.1371/journal.pone.0243919.s007

(TIF)

S1 Table. Primers used in DNA sequencing and PCR.

https://doi.org/10.1371/journal.pone.0243919.s008

(DOCX)

S2 Table. Sequencing results for the LIPC gene in NHWs (n = 95).

https://doi.org/10.1371/journal.pone.0243919.s009

(DOCX)

S3 Table. Sequencing results for the LIPC gene in African blacks (n = 95).

https://doi.org/10.1371/journal.pone.0243919.s010

(DOCX)

S4 Table. Tagger results for 161 LIPC variants (MAF≥0.05, r2≥0.8) in ABs.

https://doi.org/10.1371/journal.pone.0243919.s011

(DOCX)

S5 Table. Tagger results for the 78 LIPC variants (MAF≥0.05, r2≥0.8) in NHWs.

https://doi.org/10.1371/journal.pone.0243919.s012

(DOCX)

Citation: Pirim D, Bunker CH, Hokanson JE, Hamman RF, Demirci FY, Kamboh MI (2020) Hepatic lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels. PLoS ONE 15(12): e0243919. https://doi.org/10.1371/journal.pone.0243919

References

1. Connelly PW. The role of hepatic lipase in lipoprotein metabolism. Clinica Chimica Acta. 1999;286(1–2):243–55. pmid:10511296

2. Jansen H, Verhoeven AJ, Sijbrands EJ. Hepatic lipase a pro-or anti-atherogenic protein? Journal of lipid research. 2002;43(9):1352–62. pmid:12235167

3. Santamarina-Fojo S, Haudenschild C, Amar M. The role of hepatic lipase in lipoprotein metabolism and atherosclerosis. Current opinion in lipidology. 1998;9(3):211–9. pmid:9645503

4. Ameis D, Stahnke G, Kobayashi J, McLean J, Lee G, Büscher M, et al. Isolation and characterization of the human hepatic lipase gene. Journal of Biological Chemistry. 1990;265(12):6552–5. pmid:2324091

5. Cai SJ, Wong DM, Chen SH, Chan L. Structure of the human hepatic triglyceride lipase gene. Biochemistry. 1989;28(23):8966–71. pmid:2605236

6. Perret B, Mabile L, Martinez L, Tercé F, Barbaras R, Collet X. Hepatic lipase structure/function relationship, synthesis, and regulation. Journal of lipid research. 2002;43(8):1163–9. pmid:12177160

7. Cohen JC, Wang Z, Grundy SM, Stoesz MR, Guerra R. Variation at the hepatic lipase and apolipoprotein AI/CIII/AIV loci is a major cause of genetically determined variation in plasma HDL cholesterol levels. The Journal of clinical investigation. 1994;94(6):2377–84. pmid:7989594

8. Zambon A, Deeb SS, Pauletto P, Crepaldi G, Brunzell JD. Hepatic lipase: a marker for cardiovascular disease risk and response to therapy. Current opinion in lipidology. 2003;14(2):179–89. pmid:12642787

9. Dugi KA, Brandauer K, Schmidt N, Nau B, Schneider JG, Mentz S, et al. Low hepatic lipase activity is a novel risk factor for coronary artery disease. Circulation. 2001;104(25):3057–62. pmid:11748100

10. Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, et al. Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nature genetics. 2008;40(2):189. pmid:18193044

11. Spracklen CN, Chen P, Kim YJ, Wang X, Cai H, Li S, et al. Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels. Human molecular genetics. 2017;26(9):1770–84. pmid:28334899

12. Klarin D, Damrauer SM, Cho K, Sun YV, Teslovich TM, Honerlaw J, et al. Genetics of blood lipids among~ 300,000 multi-ethnic participants of the Million Veteran Program. Nature genetics. 2018;50(11):1514. pmid:30275531

13. Kanai M, Akiyama M, Takahashi A, Matoba N, Momozawa Y, Ikeda M, et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nature genetics. 2018;50(3):390. pmid:29403010

14. Kulminski AM, Huang J, Loika Y, Arbeev KG, Bagley O, Yashkin A, et al. Strong impact of natural-selection–free heterogeneity in genetics of age-related phenotypes. Aging (Albany NY). 2018;10(3):492. pmid:29615537

15. Below JE, Parra EJ, Gamazon ER, Torres J, Krithika S, Candille S, et al. Meta-analysis of lipid-traits in Hispanics identifies novel loci, population-specific effects, and tissue-specific enrichment of eQTLs. Scientific reports. 2016;6:19429. pmid:26780889

16. Lu X, Huang J, Mo Z, He J, Wang L, Yang X, et al. Genetic susceptibility to lipid levels and lipid change over time and risk of incident hyperlipidemia in Chinese populations. Circulation: Cardiovascular Genetics. 2016;9(1):37–44. pmid:26582766

17. Weissglas-Volkov D, Aguilar-Salinas CA, Nikkola E, Deere KA, Cruz-Bautista I, Arellano-Campos O, et al. Genomic study in Mexicans identifies a new locus for triglycerides and refines European lipid loci. Journal of medical genetics. 2013;50(5):298–308. pmid:23505323

18. Lettre G, Palmer CD, Young T, Ejebe KG, Allayee H, Benjamin EJ, et al. Genome-wide association study of coronary heart disease and its risk factors in 8,090 African Americans: the NHLBI CARe Project. PLoS genetics. 2011;7(2):e1001300. pmid:21347282

19. Kraja AT, Vaidya D, Pankow JS, Goodarzi MO, Assimes TL, Kullo IJ, et al. A bivariate genome-wide approach to metabolic syndrome: STAMPEED consortium. Diabetes. 2011;60(4):1329–39. pmid:21386085

20. Surakka I, Horikoshi M, Mägi R, Sarin A-P, Mahajan A, Lagou V, et al. The impact of low-frequency and rare variants on lipid levels. Nature genetics. 2015;47(6):589. pmid:25961943

21. Cohen JC, Kiss RS, Pertsemlidis A, Marcel YL, McPherson R, Hobbs HH. Multiple rare alleles contribute to low plasma levels of HDL cholesterol. Science. 2004;305(5685):869–72. pmid:15297675

22. Evans D, Arzer J, Aberle J, Beil F. Rare variants in the lipoprotein lipase (LPL) gene are common in hypertriglyceridemia but rare in Type III hyperlipidemia. Atherosclerosis. 2011;214(2):386–90. pmid:21159338

23. Johansen CT, Wang J, McIntyre AD, Martins RA, Ban MR, Lanktree MB, et al. Excess of rare variants in non–genome-wide association study candidate genes in patients with hypertriglyceridemia. Circulation: Cardiovascular Genetics. 2012;5(1):66–72. pmid:22135386

24. Wang J, Cao H, Ban MR, Kennedy BA, Zhu S, Anand S, et al. Resequencing genomic DNA of patients with severe hypertriglyceridemia (MIM 144650). Arteriosclerosis, thrombosis, and vascular biology. 2007;27(11):2450–5. pmid:17717288

25. Khovidhunkit W, Chartyingcharoen P, Siriwong S, Limumpornpetch P, Plengpanich W. Resequencing CETP, LIPC and LIPG genes in Thai subjects with hyperalphalipoproteinemia. The American journal of cardiology. 2012;110(1):62–6. pmid:22464213

26. Patel AP, Peloso GM, Pirruccello JP, Johansen CT, Dubé JB, Larach DB, et al. Targeted exonic sequencing of GWAS loci in the high extremes of the plasma lipids distribution. Atherosclerosis. 2016;250:63–8. pmid:27182959

27. Lee CJ, Park MS, Kim M, Ann S-j, Lee J, Park S, et al. CETP, LIPC, and SCARB1 variants in individuals with extremely high high-density lipoprotein-cholesterol levels. Scientific reports. 2019;9(1):1–7.

28. Hamman RF, Marshall JA, Baxter J, Kahn LB, Mayer EJ, Orleans M, et al. Methods and prevalence of non-insulin-dependent diabetes mellitus in a biethnic Colorado population: the San Luis Valley Diabetes Study. American Journal of Epidemiology. 1989;129(2):295–311. pmid:2912042

29. Rewers M, Shetterly S, Hoag S, Baxter J, Marshall J, Hamman R. Is the risk of coronary heart disease lower in Hispanics than in non-Hispanic whites? The San Luis Valley Diabetes Study. Ethnicity & disease. 1993;3(1):44–54. pmid:8508104

30. Bunker CH, Ukoli FA, Matthews KA, Kriska AM, Huston SL, Kuller LH. Weight threshold and blood pressure in a lean black population. Hypertension. 1995;26(4):616–23. pmid:7558221

31. Bunker CH, Ukoli FA, Okoro FI, Olomu AB, Kriska AM, Huston SL, et al. Correlates of serum lipids in a lean black population. Atherosclerosis. 1996;123(1–2):215–25. pmid:8782852

32. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clinical chemistry. 1972;18(6):499–502. pmid:4337382

33. Richmond W. Preparation and properties of a cholesterol oxidase from Nocardia sp. and its application to the enzymatic assay of total cholesterol in serum. Clinical chemistry. 1973;19(12):1350–6. pmid:4757363

34. Harris MR, Bunker CH, Hamman RF, Sanghera DK, Aston CE, Kamboh MI. Racial differences in the distribution of a low density lipoprotein receptor-related protein (LRP) polymorphism and its association with serum lipoprotein, lipid and apolipoprotein levels. Atherosclerosis. 1998;137(1):187–95. pmid:9568751

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

36. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nature protocols. 2009;4(7):1073. pmid:19561590

37. Schwarz JM, Rödelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nature methods. 2010;7(8):575. pmid:20676075

38. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nature methods. 2010;7(4):248. pmid:20354512

39. Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome research. 2012;22(9):1790–7. pmid:22955989

40. Guo L, Wang J. rSNPBase 3.0: an updated database of SNP-related regulatory elements, element-gene pairs and SNP-based gene regulatory networks. Nucleic acids research. 2017;46(D1):D1111–D6.

41. Isaacs A, Sayed-Tabatabaei FA, Njajou OT, Witteman JC, van Duijn CM. The− 514 C→ T hepatic lipase promoter region polymorphism and plasma lipids: a meta-analysis. The Journal of Clinical Endocrinology & Metabolism. 2004;89(8):3858–63.

42. Hodoğlugil U, Williamson DW, Mahley RW. Polymorphisms in the hepatic lipase gene affect plasma HDL-cholesterol levels in a Turkish population. Journal of lipid research. 2010;51(2):422–30. pmid:19734193

43. McCaskie P, Cadby G, Hung J, McQuillan B, Chapman C, Carter K, et al. The C-480T hepatic lipase polymorphism is associated with HDL-C but not with risk of coronary heart disease. Clinical genetics. 2006;70(2):114–21. pmid:16879193

44. Verdier C, Ruidavets J-B, Bongard V, Taraszkiewicz D, Martinez LO, Elbaz M, et al. Association of hepatic lipase-514T allele with coronary artery disease and ankle-brachial index, dependence on the lipoprotein phenotype: the GENES study. PloS one. 2013;8(7):e67805. pmid:23874450

45. Goodarzynejad H, Boroumand M, Behmanesh M, Ziaee S, Jalali A, Pourgholi L. Association between the hepatic lipase promoter region polymorphism (-514 C/T) and the presence and severity of premature coronary artery disease. The Journal of Tehran University Heart Center. 2017;12(3):119. pmid:29062379

46. Rosenson R.S., Brewer H.B., Barter P.J. Jr, Björkegren J.L.M., Chapman M.J., Gaudet D., et al. HDL and atherosclerotic cardiovascular disease: genetic insights into complex biology. Nat. Rev. Cardiol. 15 (2018) 9–19. pmid:28795686

47. Verma P, Verma DK, Sethi R, Singh S, Krishna A. The rs2070895 (-250G/A) single nucleotide polymorphism in hepatic lipase (HL) gene and the risk of coronary artery disease in North Indian population: a case-control study. Journal of clinical and diagnostic research: JCDR. 2016;10(8):GC01. pmid:27656463

48. Mohammadzadeh G, Ghaffari M-A, Bazyar M, Kheirollah A. Association between two common polymorphisms (single nucleotide polymorphism-250G/A and-514C/T) of the hepatic lipase gene and coronary artery disease in type 2 diabetic patients. Advanced biomedical research. 2016;5.

49. Zhao X, Ren Y, Li H, Wu Y. Association of LIPC-250G/A and-514C/T polymorphisms and hypertension: a systematic review and meta-analysis. Lipids in health and disease. 2018;17(1):238. pmid:30322388

50. Moon S, Kim YJ, Han S, Hwang MY, Shin DM, Park MY, et al. The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits. Scientific reports. 2019;9(1):1382. pmid:30718733

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

52. Pirim D, Radwan ZH, Wang X, et al. Apolipoprotein E-C1-C4-C2 gene cluster region and inter-individual variation in plasma lipoprotein levels: a comprehensive genetic association study in two ethnic groups. PLoS One. 2019;14(3):e0214060. pmid:30913229

53. Pirim D, Wang X, Niemsiri V, et al. Resequencing of the CETP gene in American whites and African blacks: Association of rare and common variants with HDL-cholesterol levels. Metabolism. 2016;65(1):36–47. pmid:26683795

54. Pirim D, Wang X, Radwan ZH, et al. Lipoprotein lipase gene sequencing and plasma lipid profile. J Lipid Res. 2014;55(1):85–93. pmid:24212298

55. Niemsiri V, Wang X, Pirim D, et al. Genetic contribution of SCARB1 variants to lipid traits in African Blacks: a candidate gene association study. BMC Med Genet. 2015;16:106. pmid:26563154

56. Niemsiri V, Wang X, Pirim D, et al. Impact of genetic variants in human scavenger receptor class B type I (SCARB1) on plasma lipid traits [published correction appears in Circ Cardiovasc Genet. 2015 Feb;8(1):245] [published correction appears in Circ Cardiovasc Genet. 2015 Feb;8(1):245]. Circ Cardiovasc Genet. 2014;7(6):838–847. pmid:25245032

57. Pirim D, Wang X, Radwan ZH, et al. Resequencing of LPL in African Blacks and associations with lipoprotein-lipid levels. Eur J Hum Genet. 2015;23(9):1244–1253. pmid:25626708

58. Radwan ZH, Wang X, Waqar F, et al. Comprehensive evaluation of the association of APOE genetic variation with plasma lipoprotein traits in U.S. whites and African blacks. PLoS One. 2014;9(12):e114618. pmid:25502880

59. Katneni UK, Liss A, Holcomb D, Katagiri NH, Hunt R, Bar H, et al. Splicing dysregulation contributes to the pathogenicity of several F9 exonic point variants. Mol Genet Genomic Med. 2019; 7(8):e840. pmid:31257730