==== Front PLoS One PLoS One plos plosone PLoS ONE 1932-6203 Public Library of Science San Francisco, CA USA 10.1371/journal.pone.0243919 PONE-D-20-29487 Research Article Biology and Life Sciences Computational Biology Genome Complexity Introns Biology and Life Sciences Genetics Genomics Genome Complexity Introns Biology and Life Sciences Biochemistry Lipids Biology and Life Sciences Biochemistry Enzymology Enzymes Hydrolases Lipases Biology and Life Sciences Biochemistry Proteins Enzymes Hydrolases Lipases Research and Analysis Methods Database and Informatics Methods Biological Databases Protein Structure Databases Biology and Life Sciences Molecular Biology Macromolecular Structure Analysis Protein Structure Protein Structure Databases Biology and Life Sciences Biochemistry Proteins Protein Structure Protein Structure Databases Biology and Life Sciences Genetics Genetic Loci Alleles Biology and Life Sciences Biochemistry Lipids Cholesterol Biology and Life Sciences Molecular Biology Macromolecular Structure Analysis Protein Structure Biology and Life Sciences Biochemistry Proteins Protein Structure Biology and life sciences Biochemistry Nucleic acids RNA Non-coding RNA Long non-coding RNA Hepatic lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels LIPC sequencing in extreme lipid phenotypeshttps://orcid.org/0000-0002-0522-9432Pirim Dilek ConceptualizationData curationFormal analysisInvestigationMethodologyValidationWriting – original draftWriting – review & editing12 Bunker Clareann H. ResourcesWriting – review & editing3 Hokanson John E. ResourcesWriting – review & editing4 Hamman Richard F. ResourcesWriting – review & editing4 Demirci F. Yesim InvestigationMethodologyProject administrationSupervisionWriting – review & editing1 Kamboh M. Ilyas ConceptualizationFormal analysisFunding acquisitionInvestigationMethodologyProject administrationResourcesSupervisionValidationVisualizationWriting – original draftWriting – review & editing1* 1 Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America 2 Department of Molecular Biology and Genetics, Faculty of Arts & Science, Bursa Uludag University, Gorukle, Bursa, Turkey 3 Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America 4 Department of Epidemiology, Colorado School of Public Health, University of Colorado Denver, Aurora, Colorado, United States of America Baroni Marco Giorgio Editor Universita degli Studi di Roma La Sapienza, ITALY Competing Interests: The authors have declared that no competing interests exist. * E-mail: kamboh@pitt.edu 16 12 2020 2020 15 12 e024391918 9 2020 1 12 2020 This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.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. http://dx.doi.org/10.13039/100000050National Heart, Lung, and Blood InstituteHL084613Kamboh M. Ilyas This study was supported by the National Heart, Lung and Blood Institute (NHLBI) grant, HL084613 (M. Ilyas Kamboh). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study. Data AvailabilityAll relevant data are within the paper and its Supporting information files.Data Availability All relevant data are within the paper and its Supporting information files. ==== Body 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. 10.1371/journal.pone.0243919.t001Table 1 Biometric and quantitative data (mean± SD) of the 95 NHW and 95 African black samples with extreme* HDL-C levels. NHWs (n = 95) Entire Sample (n = 623) African blacks (n = 95) Entire Sample (n = 788) High HDL-C (n = 47) Low HDL–C (n = 48) P-value High HDL-C (n = 48) Low HDL-C (n = 47) P-value (HDL-C range: 58–106 mg/dL) (HDL-C range: 20–40 mg/dL) (HDL-C range: 68.30–99 mg/dL) (HDL-C range: 10.30–35 mg/dL) Sex (M/F) 24/23 24/24 1 295/328 24/24 23/24 1 493/293 Age (years) 55.45 ± 9.80 53.03 ± 10.54 0.25 52.83 ± 11.41 41.29 ± 8.72 40.87 ± 7.12 0.80 40.95 ± 8.39 BMI (kg/m2) 23.17 ± 3.17 27.35 ± 3.90 1.2E-07 25.51 ± 4.06 22.06 ± 4.70 23.91 ± 5.51 0.08 22.87 ± 4.04 TC (mg/dl) 227.34 ± 51.76 208.81± 44.65 0.07 216.99 ± 43.55 201 ± 39.68 141.68 ± 31.03 2.4E-12 172.01 ± 38.47 LDL-C(mg/dl) 126.84 ± 46.95 125.54 ± 54.97 0.90 136.99 ± 40.80 112.55 ± 39.75 95.04 ± 28.28 0.02 109.25 ± 34.40 HDL-C (mg/dl) 77.68 ± 13.32 31.81 ± 4.37 2.2E-16 50.76 ± 14.35 76.05 ± 7.53 25.51 ± 5.66 2.2E-16 47.88 ± 12.87 TG (mg/dl) 114.09 ± 60.88 240.21 ± 153.22 1.7E-06 142.72 ± 93.49 61.98 ± 19.85 95.79 ± 73.21 0.004 72.96 ± 39.32 ApoB(mg/dl) 87.88 ± 25.49 89.61± 25.18 0.80 149.62 ± 33.33 66.00 ± 20.22 69.64 ± 21.46 0.40 66.98 ± 22.19 ApoA1(mg/dl) 174.08 ± 3.57 130.20 ± 2.78 1.4E-06 87.72 ± 24.27 166.04 ± 28.19 103.84 ± 27.23 2.2E-16 137.03 ± 28.46 P-values were calculated based on the original values by using t-test. No covariates were included. *Adjusted for sex and age: High/Low HDL-C groups correspond to ≥90th % tile and ≤10th % tile of the HDL-C distribution. TC: Total cholesterol; LDL-C: Low-density lipoprotein cholesterol; HDL-C; high-density lipoprotein cholesterol; TG: Triglycerides. 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). 10.1371/journal.pone.0243919.g001Fig 1 Locations of the LIPC variants identified in NHWs (n = 95) and ABs (n = 95). 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. 10.1371/journal.pone.0243919.g002Fig 2 Distributions of the minor allele frequencies of LIPC variants in NHWs (n = 95) and ABs (n = 95). 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. 10.1371/journal.pone.0243919.g003Fig 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). 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. 10.1371/journal.pone.0243919.t002Table 2 LIPC coding variants identified in 95 NHWs and 95 African blacks with extreme HDL-C levels. African blacks NHWs gnomAD (population)* Ref SNP ID Location Amino acid change Alleles High HDL-C MAF Low HDL-C MAF MAF Alleles High HDL-C MAF Low HDL-C MAF MAF SIFT PP2 MT MAF rs755990193 Exon 2 p.T44T G>A 0 0.011 0.005 G>A - - - P 1.09E-04 (Asian) rs113174258 Exon 2 p.T71T G>A 0 0.011 0.005 G>A - - - DC 7.18E-04 rs7175412 Exon 2 p.H88H C>T 0.042 0.087 0.064 C>T - - - DC 0.004 rs6078 Exon 3 p.V95M G>A 0.076 0.06 0.067 G>A 0.011 0 0.005 T B P 0.07 rs776118661 Exon 3 p.H127H C>T 0.011 0.00 0.006 C>T - - - P 2.01E-05 rs690 Exon 4 p.V155V T>G 0.5 0.44 0.473 T>G 0.5 0.362 0.428 P 0.5 rs6082 Exon 5 p.G197G A>G 0.115 0.02 0.07 A>G 0.076 0.064 0.07 P 0.459 rs6083 Exon 5 p.N215S G>A 0.219 0.244 0.231 A>G 0.38 0.34 0.36 T B P 0.474 rs6084 Exon 5 p.T224T C>G 0.188 0.200 0.194 G>C 0.424 0.5 0.462 P 0.424 rs146299102 Exon 6 p.H279H C>T 0.021 0.022 0.022 C>T - - - DC 0.001 rs121912502 Exon 6 p.S289F C>T - - - C>T 0.011 0 0.005 D PD DC 0.001 rs145811475 Exon 6 p.V342I G>A 0 0.011 0.005 G>A - - - T B P 1.06E-04 rs3829462 Exon 7 p.F356L A>C 0.032 0.03 0.031 A>C - - - T B P 0.031 rs3829461 Exon 7 p.T366T G>A 0.032 0.03 0.031 G>A - - - P 0.031 rs113298164 Exon 8 p.T405M C>T - - - C>T 0.011 0 0.005 D PD DC 0.003 rs75983069 Exon 8 p.P438P A>G 0.032 0.02 0.027 A>G - - - T DC 7.46E-04 rs6074 Exon 9 p.T479T C>A 0.062 0.085 0.074 C>A 0.174 0.104 0.138 T P 0.188 T: Tolerated, B: Benign, PD: Probably damaging, DC: Disease causing, P: Polymorphism, Bold indicates non-synonymous variants. *Study-wide MAF is given if variant is found more than one population in the gnomAD (Genome Aggregation Database). 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). 10.1371/journal.pone.0243919.t003Table 3 Distributions of the regulatory variants (rSNPs) identified in 190 individuals (95 NHWs and 95 ABs) with extreme HDL-C levels. RefSNP ID Alleles Location Call rate HW-P High HDL-C MAF Low HDL-C MAF Total MAF RegulomeDB Score Related regulatory elements non-Hispanic Whites (n = 95) rs143186931 A>G Intron 7 95.8 1 0.023 0 0.011 5 circRNA region rs11071390 A>G Intron 7 95.8 1 0 0.011 0.005 5 circRNA region rs144831345 G>A Intron 7 95.8 1 0 0.011 0.005 5 circRNA region rs117852639 C>A Intron 7 95.8 1 0 0.021 0.011 5 circRNA region rs113298164 C>T Exon 8 96.8 1 0.011 0 0.005 5 circRNA region rs11631342 A>G Intron 1 97.9 1 0.022 0 0.011 4 lncRNA region rs4774305 G>C Intron 7 97.9 1 0 0.01 0.005 4 circRNA region rs139878091 A>G Intron 7 98.9 1 0.011 0 0.005 5 circRNA region rs117911817 G>A Intron 7 98.9 1 0.033 0 0.016 5 circRNA region rs1869129 T>C Intron 7 100 1 0 0.01 0.005 3a circRNA region rs1869130 C>T Intron 7 100 1 0 0.01 0.005 4 circRNA region rs12438032 G>A Intron 7 100 1 0 0.01 0.005 5 circRNA region rs34964641 T>G Intron 7 100 1 0 0.01 0.005 5 circRNA region rs35925692 Ins1 Intron 7 100 1 0 0.01 0.005 5 circRNA region rs1869131 A>T Intron 7 100 1 0 0.01 0.005 5 circRNA region rs4775079 C>T Intron 7 100 1 0 0.01 0.005 5 circRNA region rs8026372 A>G Intron 7 100 1 0 0.01 0.005 5 circRNA region rs1839928 A>G Intron 7 100 1 0 0.01 0.005 5 circRNA region rs1839927 A>G Intron 7 100 1 0 0.01 0.005 5 circRNA region rs8030903 T>C Intron 7 100 1 0 0.01 0.005 4 circRNA region rs10851636 C>T Intron 7 100 1 0 0.01 0.005 4 circRNA region rs7170227 G>A Intron 7 100 1 0 0.01 0.005 3a circRNA region rs35412158 G>A Intron 7 100 1 0 0.01 0.005 5 circRNA region African Blacks (n = 95) rs35631005 C:T Intron 6 97.9 1 0.011 0 0.005 4 rs533300601 G:A Intron 7 97.9 0.0324 0.033 0 0.016 3a circRNA region rs16940468 T:G Intron 2 96.8 1 0.00 0.011 0.005 5 circRNA region rs12909325 G:A Intron 2 96.8 1 0.00 0.021 0.011 5 circRNA region rs115408618 G:C Intron 7 96.8 1 0.023 0 0.011 3a circRNA region rs190375050 C:T Intron 7 98.9 1 0.00 0.021 0.011 5 circRNA region rs181084356 C:G Intron 7 98.9 1 0.00 0.021 0.011 5 circRNA region rs568646677 G:T Intron 7 98.9 1 0.011 0 0.005 4 circRNA region rs143889538 G:A Intron 7 98.9 1 0.021 0 0.011 2b circRNA region 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 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. (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. 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). (TIF) Click here for additional data file. S1 Table Primers used in DNA sequencing and PCR. (DOCX) Click here for additional data file. S2 Table Sequencing results for the LIPC gene in NHWs (n = 95). (DOCX) Click here for additional data file. S3 Table Sequencing results for the LIPC gene in African blacks (n = 95). (DOCX) Click here for additional data file. S4 Table Tagger results for 161 LIPC variants (MAF≥0.05, r2≥0.8) in ABs. (DOCX) Click here for additional data file. S5 Table Tagger results for the 78 LIPC variants (MAF≥0.05, r2≥0.8) in NHWs. (DOCX) Click here for additional data file. 10.1371/journal.pone.0243919.r001 Decision Letter 0 Baroni Marco Giorgio Academic Editor © 2020 Marco Giorgio Baroni2020Marco Giorgio BaroniThis 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.Submission Version0 19 Oct 2020 PONE-D-20-29487 Hepatic Lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels PLOS ONE Dear Dr. Pirim, Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process. ============================== ACADEMIC EDITOR: The paper is interesting, with good methodologies. The reviewers have raised some points of concern, that the authors must address in their rebuttal. Please submit your revised manuscript by Nov 28 2020 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file. Please include the following items when submitting your revised manuscript: A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'. A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'. An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'. If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter. If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols We look forward to receiving your revised manuscript. Kind regards, Marco Giorgio Baroni, MD, PhD Academic Editor PLOS ONE Journal Requirements: When submitting your revision, we need you to address these additional requirements. 1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf 2. Thank you for stating in your Funding Statement: "This study was supported by the National Heart, Lung and Blood Institute (NHLBI) grant, HL084613 (M. Ilyas Kamboh). ". i) Please provide an amended statement that declares *all* the funding or sources of support (whether external or internal to your organization) received during this study, as detailed online in our guide for authors at http://journals.plos.org/plosone/s/submit-now.  Please also include the statement “There was no additional external funding received for this study.” in your updated Funding Statement. ii) Please include your amended Funding Statement within your cover letter. We will change the online submission form on your behalf. [Note: HTML markup is below. Please do not edit.] Reviewers' comments: Reviewer's Responses to Questions Comments to the Author 1. Is the manuscript technically sound, and do the data support the conclusions? The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. Reviewer #1: Yes Reviewer #2: Yes ********** 2. Has the statistical analysis been performed appropriately and rigorously? Reviewer #1: Yes Reviewer #2: Yes ********** 3. Have the authors made all data underlying the findings in their manuscript fully available? The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified. Reviewer #1: Yes Reviewer #2: Yes ********** 4. Is the manuscript presented in an intelligible fashion and written in standard English? PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here. Reviewer #1: Yes Reviewer #2: Yes ********** 5. Review Comments to the Author Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters) Reviewer #1: In this study, Pirim D and colleagues performed a resequencing of the LIPC gene in 95 non-Hispanic Whites (NHWs) and 95 African blacks (ABs) selected from the upper and lower 10th percentile of HDL cholesterol distribution. The aim of the study was to identify and understand the role of common, uncommon and rare variants that may influence plasma HDL-C levels. They identified a total of 464 variants, including 43 novel. To assess their functional relevance, in silico functional analyses were executed. Two nonsynonymous variants (p.S289F, p.T405M), found in NHWs with higher HDL-C levels were predicted to have damaging effect on LIPC protein. Furthermore, the authors found several non-coding variants that possibly reside in the circRNA and lncRNA binding sites and may have regulatory function. The authors conclude that this study highlights the importance of LIPC polymorphisms, common and rare, in influencing plasma HDL-C levels and that functional studies are needed to further evaluate the roles of non-coding regulatory variants in the lipid metabolism. The results presented in the paper are well written and easy to follow and provide additional evidences on the importance of the LIPC gene polymorphisms in the modulation of plasma HDL-C levels. Nonetheless, I would have some comments to the authors. 1) It’s not clear to me why the authors state that they identified a total of 464 variants if summing the 122 shared in both populations, the 156 unique to NHWs and the 234 unique to Abs, the variants are 512 in total. 2) I found a discrepancy between S2 Table where are reported 180 variants in NHWs of which 88 are common and the test in page 8 “Of the 178 variant identified in NHWs, 86 were common (MAF≥0.05)….”. In paragraph 3.2 (line 4) the authors state they found 88 common variants. And the same for S3 Table: 358 variants of which 174 common in AB cohort and in the test on page 8 “Among the 356 variants found in ABs, 172 were common (MAF≥0.05)……” Could the authors explain the differences? 3) in paragraph 3.2, it’s not clear why 161 variants in AB cohort were analysed (S4 Table) if the common variants were 172. I wonder if the 3 variants deviating from H-W equilibrium and the 10 with low call rate have been eliminated. However, in this case the variants analysed in S4 Table should be 159. This should be clarified. 4) I’ve a doubt how the authors calculated the RegulomeDB Score in Table 3. Indeed, I tried to calculate the Score selecting randomly few snps. However, the results are different from those reported in table 3. May the authors explain these different results? I wonder if the RegulomeDB scores calculated in S2 and S3 Tables present the same discrepancy. NHWs RefSNP ID RegulomeDB rank rs117911817 5 rs1869129 3a rs1869130 4 rs12438032 5 rs34964641 5 rs35925692 5 rs1869131 5 rs4775079 5 ABs RefSNP ID RegulomeDB rank rs35631005 4 rs533300601 3a rs16940468 5 rs12909325 5 rs115408618 3a rs190375050 5 rs181084356 5 rs568646677 4 rs143889538 2b Reviewer #2: Pirim and collaborator sequenced LIPC gene in 190 subjects selected from the upper and lower 10th percentile of HDL levels in two cohorts, African Blacks and non-Hispanic Whites. 464 variants were recorded, 43 of which were novel. SNPs were evaluated using some predictive online software for coding and non-coding (regulatory) variants. Among all, 17 coding variants were identified: 6 cause aminoacid change while 11 were synonymous. Authors analysis were focused on difference in MAF distribution within the extremes of HDL levels. Overall, 3 variants showed a significant different distribution. In particularly, MAF were lower in low-HDL than high-HDL group. Moreover, authors found 3 common variants associated with extremes of HDL levels. Other findings are in line with previous published studies. Language and exposition are very clear. Also the amount of data, figures, tables and supplementary materials provided make easy to follow authors argument. Interestingly, the analysis considers cohorts from different ethnic group, non-Hispanic White and African Black. Minor revision: Lacks of power calculation (even if cited at page 15). As it would be useful for readers and for fully understand study design, could authors add a power calculation in statistical analysis paragraph, 2.4? Conclusion paragraphs seems too short. Could authors add a brief summary of findings, or most promising variants, or some hypothesis for observed effects in coding and regulatory variants, or some future prospective of the study? at page 8. On paragraph 3.1 the common variants found in NWHs are n=86, while in the same page, but in paragraph 3.2, are n=88. Please correct or explain. ********** 6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files. If you choose “no”, your identity will remain anonymous but your review may still be made public. Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy. Reviewer #1: No Reviewer #2: No [NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.] While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step. 10.1371/journal.pone.0243919.r002 Author response to Decision Letter 0 Submission Version1 19 Nov 2020 Dear Editors: Thank you for the opportunity to revise our manuscript. We are pleased with the overall positive review and found the comments very helpful in improving our manuscript. We have revised the manuscript accordingly (highlighted in red in the manuscript) and below we provide a point-by-point response to each comment in red. The Funding Statement has been also updated as "This study was supported by the National Heart, Lung and Blood Institute (NHLBI) grant, HL084613 (M. Ilyas Kamboh). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study ". Reviewer #1: In this study, Pirim D and colleagues performed a resequencing of the LIPC gene in 95 non-Hispanic Whites (NHWs) and 95 African blacks (ABs) selected from the upper and lower 10th percentile of HDL cholesterol distribution. The aim of the study was to identify and understand the role of common, uncommon and rare variants that may influence plasma HDL-C levels. They identified a total of 464 variants, including 43 novel. To assess their functional relevance, in silico functional analyses were executed. Two nonsynonymous variants (p.S289F, p.T405M), found in NHWs with higher HDL-C levels were predicted to have damaging effect on LIPC protein. Furthermore, the authors found several non-coding variants that possibly reside in the circRNA and lncRNA binding sites and may have regulatory function. The authors conclude that this study highlights the importance of LIPC polymorphisms, common and rare, in influencing plasma HDL-C levels and that functional studies are needed to further evaluate the roles of non-coding regulatory variants in the lipid metabolism. The results presented in the paper are well written and easy to follow and provide additional evidences on the importance of the LIPC gene polymorphisms in the modulation of plasma HDL-C levels. Nonetheless, I would have some comments to the authors. Response: We greatly appreciate the reviewer for noting that our manuscript is well written, easy to follow and the importance of novel findings. 1) It’s not clear to me why the authors state that they identified a total of 464 variants if summing the 122 shared in both populations, the 156 unique to NHWs and the 234 unique to Abs, the variants are 512 in total. Response: We apologize for this typo as the number of variants unique to NHWs should have been 56 not "156", so the total number of identified variants is 412. Abstract, results and discussion sections have now been edited accordingly. 2) I found a discrepancy between S2 Table where are reported 180 variants in NHWs of which 88 are common and the test in page 8 “Of the 178 variant identified in NHWs, 86 were common (MAF≥0.05)….”. In paragraph 3.2 (line 4) the authors state they found 88 common variants. And the same for S3 Table: 358 variants of which 174 common in AB cohort and in the test on page 8 “Among the 356 variants found in ABs, 172 were common (MAF≥0.05)……” Could the authors explain the differences? Response: The discrepancy was due to a triallelic variant (rs7171818) that was found in both populations. Including the triallelic variant, the total number of common variants identified in NHWs and ABs were 86 and 172, respectively. However, since the three alleles for the rs7171818 were presented in 3 separate rows, the number of total variants falsely got added to 88 and 174 in NHWs and ABs, respectively. In order to avoid this apparent mis-representation, now the triallelic variant in Table S2 and S3 is presented in a single row. Of note, Haploview analyses was conducted for all variants, including three alleles for the triallelic variant. For further clarification, we have edited the text on Page 8, Section "TagSNP identification and Linkage Disequilibrium (LD) analyses", line 2-6 as follows: “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).” 3) in paragraph 3.2, it’s not clear why 161 variants in AB cohort were analysed (S4 Table) if the common variants were 172. I wonder if the 3 variants deviating from H-W equilibrium and the 10 with low call rate have been eliminated. However, in this case the variants analysed in S4 Table should be 159. This should be clarified. Response: In ABs, 174 common variants (including the triallelic variant rs7171818) were analyzed. "172" was a typo, which is now corrected on Page 8, Section 3.2, line 2-4, as indicated above. 4) I’ve a doubt how the authors calculated the RegulomeDB Score in Table 3. Indeed, I tried to calculate the Score selecting randomly few snps. However, the results are different from those reported in table 3. May the authors explain these different results? I wonder if the RegulomeDB scores calculated in S2 and S3 Tables present the same discrepancy. NHWs RefSNP ID RegulomeDB rank rs117911817 5 rs1869129 3a rs1869130 4 rs12438032 5 rs34964641 5 rs35925692 5 rs1869131 5 rs4775079 5 ABs RefSNP ID RegulomeDB rank rs35631005 4 rs533300601 3a rs16940468 5 rs12909325 5 rs115408618 3a rs190375050 5 rs181084356 5 rs568646677 4 rs143889538 2b Response: Thank you for pointing this out. We used the RegulomeDB version 1 in our analyses and since then the database has been updated, resulting in difference scores for some variants than the previous version. We have now reanalyzed all variants by using RegulomeDB v2 and updated all scores in the tables and related text accordingly. Reviewer #2: Pirim and collaborator sequenced LIPC gene in 190 subjects selected from the upper and lower 10th percentile of HDL levels in two cohorts, African Blacks and non-Hispanic Whites. 464 variants were recorded, 43 of which were novel. SNPs were evaluated using some predictive online software for coding and non-coding (regulatory) variants. Among all, 17 coding variants were identified: 6 cause aminoacid change while 11 were synonymous. Authors analysis were focused on difference in MAF distribution within the extremes of HDL levels. Overall, 3 variants showed a significant different distribution. In particularly, MAF were lower in low-HDL than high-HDL group. Moreover, authors found 3 common variants associated with extremes of HDL levels. Other findings are in line with previous published studies. Language and exposition are very clear. Also the amount of data, figures, tables and supplementary materials provided make easy to follow authors argument. Interestingly, the analysis considers cohorts from different ethnic group, non-Hispanic White and African Black. Response: We greatly appreciate the reviewer for noting multiple strengths in our study. Minor revision: Lacks of calculation (even if cited at page 15). As it would be useful for readers and for fully understand study design, could authors add a power calculation in statistical analysis paragraph, 2.4? Response: We have now calculated the power by using G power based on our sample size and the text has been revised as below on Page 6, Section Statistical and Bioinformatic Analysis, line 4-6: “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).” Conclusion paragraphs seems too short. Could authors add a brief summary of findings, or most promising variants, or some hypothesis for observed effects in coding and regulatory variants, or some future prospective of the study? Response: Since the summary of the findings is already highlighted in the Abstract and the promising variants discussed in the Discussion section, we avoided their repetitions in the Conclusions. However, in deference to the reviewer's suggestion, we have added the the below text in the Conclusions section on Page 16, Section 5, line 3-6: “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....” at page 8. On paragraph 3.1 the common variants found in NWHs are n=86, while in the same page, but in paragraph 3.2, are n=88. Please correct or explain. Response: As responded to under Reviewer 1, in NHWs, we identified 86 common variants, including one triallelic SNP (rs7171818) and for clarification, we have edited the text on Page 8, Section TagSNP identification and Linkage Disequilibrium (LD) analyses, line 2-6. We hope that we have satisfactorily addressed reviewers’ concerns and look forward to learning a favorable response. Sincerely, M. Ilyas Kamboh, PhD, FAHA On behalf of all authors Attachment Submitted filename: Response to Reviewers.docx Click here for additional data file. 10.1371/journal.pone.0243919.r003 Decision Letter 1 Baroni Marco Giorgio Academic Editor © 2020 Marco Giorgio Baroni2020Marco Giorgio BaroniThis 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.Submission Version1 1 Dec 2020 Hepatic Lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels PONE-D-20-29487R1 Dear Dr. Pirim, We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements. Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication. An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org. If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org. Kind regards, Marco Giorgio Baroni, MD, PhD Academic Editor PLOS ONE Additional Editor Comments (optional): Reviewers' comments: Reviewer's Responses to Questions Comments to the Author 1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation. Reviewer #1: All comments have been addressed Reviewer #2: All comments have been addressed ********** 2. Is the manuscript technically sound, and do the data support the conclusions? The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented. Reviewer #1: Yes Reviewer #2: Yes ********** 3. Has the statistical analysis been performed appropriately and rigorously? Reviewer #1: Yes Reviewer #2: Yes ********** 4. Have the authors made all data underlying the findings in their manuscript fully available? The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified. Reviewer #1: Yes Reviewer #2: Yes ********** 5. Is the manuscript presented in an intelligible fashion and written in standard English? PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here. Reviewer #1: Yes Reviewer #2: Yes ********** 6. Review Comments to the Author Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters) Reviewer #1: Dear authors, I read the paper carefully and the current form is suitable to be published. I really appreciated your collaboration in improving your manuscript. Reviewer #2: (No Response) ********** 7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files. If you choose “no”, your identity will remain anonymous but your review may still be made public. Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy. Reviewer #1: No Reviewer #2: No 10.1371/journal.pone.0243919.r004 Acceptance letter Baroni Marco Giorgio Academic Editor © 2020 Marco Giorgio Baroni2020Marco Giorgio BaroniThis 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. 4 Dec 2020 PONE-D-20-29487R1 Hepatic Lipase (LIPC) sequencing in individuals with extremely high and low high-density lipoprotein cholesterol levels Dear Dr. Pirim: I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department. If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org. If we can help with anything else, please email us at plosone@plos.org. Thank you for submitting your work to PLOS ONE and supporting open access. Kind regards, PLOS ONE Editorial Office Staff on behalf of Prof Marco Giorgio Baroni Academic Editor PLOS ONE ==== Refs References 1 Connelly PW . The role of hepatic lipase in lipoprotein metabolism . Clinica Chimica Acta . 1999 ;286 (1–2 ):243 –55 . 10.1016/s0009-8981(99)00105-9 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 . 10.1194/jlr.r200008-jlr200 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 . 10.1097/00041433-199806000-00005 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 . 2324091 5 Cai SJ , Wong DM , Chen SH , Chan L . Structure of the human hepatic triglyceride lipase gene . Biochemistry . 1989 ;28 (23 ):8966 –71 . 10.1021/bi00449a002 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 . 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 . 10.1172/JCI117603 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 . 10.1097/00041433-200304000-00010 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 . 10.1161/hc5001.100795 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 10.1038/ng.75 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 . 10.1093/hmg/ddx062 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 10.1038/s41588-018-0222-9 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 10.1038/s41588-018-0047-6 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 10.18632/aging.101407 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 10.1038/srep19429 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 . 10.1161/CIRCGENETICS.115.001096 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 . 10.1136/jmedgenet-2012-101461 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 10.1371/journal.pgen.1001300 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 . 10.2337/db10-1011 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 10.1038/ng.3300 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 . 10.1126/science.1099870 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 . 10.1016/j.atherosclerosis.2010.11.026 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 . 10.1161/CIRCGENETICS.111.960864 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 . 10.1161/ATVBAHA.107.150680 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 . 10.1016/j.amjcard.2012.02.052 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 . 10.1016/j.atherosclerosis.2016.04.011 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 .30626917 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 . 10.1093/oxfordjournals.aje.a115134 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 . 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 . 10.1161/01.hyp.26.4.616 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 . 10.1016/0021-9150(96)05810-8 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 . 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 . 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 . 10.1016/s0021-9150(97)00230-x 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 . 10.1093/bioinformatics/bth457 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 10.1038/nprot.2009.86 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 10.1038/nmeth0810-575 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 10.1038/nmeth0410-248 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 . 10.1101/gr.137323.112 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 .15292318 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 . 10.1194/jlr.P001578 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 . 10.1111/j.1399-0004.2006.00659.x 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 10.1371/journal.pone.0067805 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 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 . 10.1038/nrcardio.2017.115 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 10.7860/JCDR/2016/20496.8378 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 10.1186/s12944-018-0884-4 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 10.1038/s41598-018-37832-9 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 10.1038/ng.2926 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 10.1371/journal.pone.0214060 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 . 10.1016/j.metabol.2015.09.020 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 . 10.1194/jlr.M043265 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 10.1186/s12881-015-0250-6 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 . 10.1161/CIRCGENETICS.114.000559 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 . 10.1038/ejhg.2014.268 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 10.1371/journal.pone.0114618 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 10.1002/mgg3.840 31257730