Exome Sequencing Identifies 2 Rare Variants for Low High-Density Lipoprotein Cholesterol in an Extended FamilyClinical Perspective
Background—Exome sequencing is a recently implemented method to discover rare mutations for Mendelian disorders. Less is known about its feasibility to identify genes for complex traits. We used exome sequencing to search for rare variants responsible for a complex trait, low levels of serum high-density lipoprotein cholesterol (HDL-C).
Methods and Results—We conducted exome sequencing in a large French-Canadian family with 75 subjects available for study, of which 27 had HDL-C values less than the fifth age-sex-specific population percentile. We captured ≈50 Mb of exonic and transcribed sequences of 3 closely related family members with HDL-C levels less than the fifth age-sex percentile and sequenced the captured DNA. Approximately 82 000 variants were detected in each individual, of which 41 rare nonsynonymous variants were shared by the sequenced affected individuals after filtering steps. Two rare nonsynonymous variants in the ATP-binding cassette, subfamily A (ABC1), member 1 (ABCA1), and lipoprotein lipase genes predicted to be damaging were investigated for cosegregation with the low HDL-C trait in the entire extended family. The carriers of either variant had low HDL-C levels, and the individuals carrying both variants had the lowest HDL-C values. Interestingly, the ABCA1 variant exhibited a sex effect which was first functionally identified, and, subsequently, statistically demonstrated using additional French-Canadian families with ABCA1 mutations.
Conclusions—This complex combination of 2 rare variants causing low HDL-C in the extended family would not have been identified using traditional linkage analysis, emphasizing the need for exome sequencing of complex lipid traits in unexplained familial cases.
Low high-density lipoprotein cholesterol (HDL-C) is the most common lipoprotein abnormality and established risk factor of coronary heart disease (CHD). Low HDL-C is caused by multiple genetic factors, common and rare, interacting with one another and with the environment and behavior. In the last 2 decades, significant effort has been devoted to the identification of low HDL-C susceptibility genes. This was initially done using the genome-wide linkage analysis.1,2 However, progress in identification of the actual disease genes was very slow despite the discovery of many linked intervals. More recently, genome-wide association studies (GWAS) have successfully identified multiple common variants associated with decreased levels of HDL-C.3
Clinical Perspective on p 546
However, the sum of common variants identified so far through GWAS explains only a small fraction (10%–15%) of the variance in the HDL-C levels.3 Hence, it has become evident that other types of DNA variants must contribute substantially to HDL-C levels as well. To identify new rare and low-frequency variants underlying low HDL-C, massive parallel sequencing technologies can be used. The whole-genome sequencing is the most complete approach, but it remains significantly more expensive than exome sequencing that only analyzes coding and transcribed regions, which constitutes <5% of the whole genome sequence.4 It is estimated that the protein-coding regions of the human genome constitute about 85% of the disease-causing mutations.4
We used whole-exome sequencing to search for rare variants conferring susceptibility to low HDL-C. We sequenced the exomes of closely related family members with low HDL-C from a large multigenerational French-Canadian family with 75 subjects available for study and followed up the candidate variants by examining the cooccurrence patterns in the entire extended family.
The study sample consists of a large multigenerational French-Canadian family collected in the Cardiovascular Genetics Laboratory, McGill University Health Centre, Royal Victoria Hospital, Montreal, Canada, as described previously.5 There are 75 family members (35 males and 40 females) with both DNA and extensive demographic and clinical phenotype information available for study in this family. We selected 3 closely related family members with HDL-C levels less than the fifth age-sex percentile from the uppermost generations (Figure 1) for exome sequencing to focus on most severe cases and avoid genetic heterogeneity typical for complex lipid traits.
For a gene–sex interaction analysis, 10 additional French-Canadian families with previously identified mutations in ABCA16–8 comprising 125 individuals were also included in the study. The affection status in all families was determined using the fifth age-sex specific population percentile of HDL-C.5 Family members were sampled (blood collection for lipoprotein analyses, DNA isolation for genetic studies, and skin biopsy for culture of skin fibroblasts used in cellular cholesterol efflux assays) after a 12-hour fast and discontinuation of lipid modifying medications for >4 weeks. Lipids and lipoproteins were measured using standardized techniques as described previously.6,9 The research protocol was approved by the Research Ethics Board of the McGill University Health Center, and all subjects gave informed consent.
Library Construction and Sequencing
Library construction was performed using 3 μg of genomic DNA and Agilent SureSelect All Exon Kit (50-Mb design) according to the manufacturer’s instructions. Further details of library construction and sequencing are given in the online-only Data Supplement Methods.
We converted the qseq files into a Sanger-formatted FASTQ files that were aligned to a reference sequence (hg19) using the default options of the Burroughs-Wheeler Aligner.10 Duplicates were removed and a pileup file was generated using SAM tools.11 The pileup file was used to run the quality control metrics including: a minimum read depth of 4, a maximum read depth of 600, a maximum of 2 single nucleotide polymorphisms (SNPs) per a window size of 10 bases, and a minimum indel score of 25 for filtering nearby SNPs and Phred quality >40. The BED file supplied by Agilent was used to filter only those reads corresponding to the 50 Mb targets.
Annovar was used for functional annotation, dividing the variants into coding and noncoding variants.12 The coding variants were further divided into synonymous, nonsynonymous (missense), and stop gain or stop loss variants. The synonymous variants were subsequently discarded because they are less likely to be causal. The variants were filtered against the variants present in the HapMap,13 The 1000 Genomes Project,14 and dbSNP13214 databases. Along with novel variants, we selected known rare variants with a minor allele frequency <5%. These variants were classified into damaging and benign based on their predicted protein effect using PolyPhen15 and SIFT.16
Parametric Linkage and Association Analysis
Two-point parametric linkage analysis was performed in the extended family using the Location-Score option of the Mendel software17 as described in detail in the online-only Data Supplement Methods. Association analysis was performed using a measured genotype approach utilizing the Polygenic-QTL option of Mendel,18 using continuous HDL-C levels with age and sex as covariates and allele counts of either the ABCA1 variant, lipoprotein lipase (LPL) variant, both variants or none (ie, null model). The heritability and variance explained were calculated as the percent change in total and genetic variance between the null model and the models, including the genotypes as covariates. The LPL variant was further tested for association with log transformed triglyceride (TG) values in a similar fashion.
Genotype by Sex Interaction
We included the extended family together with 10 additional families with previously identified mutations in ABCA16–8 in a gene–sex interaction analysis, comprising 200 individuals and 9 different mutations in ABCA1 (DelED1893, G616V, K776N, N1800H, Q2210H, R1851X, R2084X, R909X, and S1731C). Genotype by sex interaction was tested by the SOLAR program19 using variance-component analysis for discrete traits. We compared models with and without the gene–sex interaction term while keeping the ABCA1 genotypes in both the null and interaction model. We assumed a dominant genetic inheritance, classifying carriers of a mutation as 1 and 0 otherwise, and a multiplicative interaction term, multiplying the genotype score by sex (men=1 and women=0). We also coded a sex-interaction term in which men and postmenopausal women (≥50 years of age) were coded as 1 and premenopausal women (<50 years of age) were coded as 0. Subjects with HDL-C levels less than the age-sex specific 10th percentiles were classified as affected and subjects with HDL-C levels more than the age-sex specific 20th percentiles as unaffected. P values were generated by comparing the 2 models using a likelihood ratio statistic with 1 degree of freedom. Because the affection status is adjusted for sex, the inclusion of the main effect of sex in the model was no longer necessary. The binary HDL-C affection was tested because the variance of HDL-C levels in these ascertained families is reduced and thus limited for effective quantitative analysis.1
Human skin fibroblasts were obtained from 3.0-mm punch biopsies of the forearm of a healthy control subject and the affected proband homozygous for the ABCA1 S1731C variant. The fibroblasts were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 0.1% nonessential amino acids, penicillin (100 U/mL), streptomycin (100 μg/mL), and 10% fetal bovine serum.
To identify rare genetic variants underlying low HDL-C, we sequenced the entire exomes (≈50 Mb) of 3 family members with HDL-C less than the fifth age-sex specific population percentile from a large multigenerational French-Canadian family. The 3 sequenced family members were closely related: 1 affected, his sibling, and his child were sequenced. Exome capture and sequencing were performed using the Agilent Sure Select in-solution method and Illumina Hiseq2000 platform as described in the Methods. We obtained an average of 90 million reads per person and successfully mapped ≈90% of these reads to the reference sequence (Table 1). After quality control the mean coverage was 50X.
Filtering of Identified Variants
On average, 82 000 single-nucleotide variants were detected in each individual. We focused on variants shared by all 3 exome sequenced subjects and filtered the variants based on their type, frequency, and functional predictions. Filtering for missense and stop gain or stop loss variants that were shared by all 3 affected individuals resulted in 3428 nonsynonymous variants and 31 stop gain or loss variants (Table 2). The transition/transversion (Ti/Tv) ratio of the coding variants was 3.4, whereas the Ti/Tv ratio of the noncoding was 2.5, in good agreement with the expected ratios.22
The identified variants were further filtered against variants present in the HapMap,13 1000 Genomes Project,14 and dbSNP13214 databases, resulting in 332 novel variants and known variants with minor allele frequency <5%. These variants were further filtered by selecting variants predicted to affect protein function using PolyPhen15 and SIFT16 and expressed in a relevant tissue including liver, adipose, and heart, resulting in 41 shared potentially functional variants that were either novel or known but relatively rare (Table 2). Among the shared variants there were 2 rare functional variants in the ABCA1 and LPL genes that are excellent susceptibility candidates as their key role in HDL-C metabolism is well established.23 We confirmed their presence by both Sanger sequencing and genotyping.
Rare Missense Variant in ABCA1
The ABCA1 (S1731C) variant is not present in dbSNP13214 or The 1000 Genomes Project14 data, and is located in exon 38. This previously reported rare variant is changing a conserved amino acid from serine to cysteine, and is known to result in decreased cholesterol efflux.7,24–26
To further determine the effect of the S1731C variant on cholesterol efflux, we used human fibroblasts from the affected proband homozygous for the S1731C variant, and compared these cells to a normal control. Assays were performed in 22OH/9CRA stimulated fibroblasts (to induce ABCA1 expression), and unstimulated cells, in the presence or absence of lipid free ApoA-I (Figure 2A). We observed a significant decrease (≈40%) in apoA-I-mediated cellular cholesterol efflux in the proband, compared with the control without the ABCA1 variant (P=1.23×10–4 using Student t-test and P=0.0495 using a nonparametric 2-sample Wilcoxon rank sum test). These results are in agreement with previously documented findings.7,25,26 Low efflux levels were also observed in unstimulated cells, presumably due to basal levels of ABCA1 expression and the presence of other apoA-I binding sites at the cell surface. Also, as expected, background basal conditions of passive diffusion of cellular cholesterol were not affected by mutations at the ABCA1 gene locus.
Because the lipid levels of the ABCA1 S1731C variant carriers suggested a possible gene-sex effect (Tables 3 and 4), we further investigated whether exposure to 17β-estradiol steroid hormone endogenously expressed in females, possibly corrects the cholesterol efflux defect in fibroblasts from the S1731C male ABCA1 carrier during the 22OH/9CRA ABCA1 stimulation phase of 17 hours (Figure 2B). Interestingly, after adjusting for basal cholesterol diffusion, we observed that upon treatment with elevated doses of estradiol (>20 nM), efflux in the S1731C proband significantly increased (P=7.2×10–6, r=0.78 using a nonparametric Spearman trend test), whereas that in the wild-type control remained constant (P=0.2, r=0.25) (Figure 2B). Taken together, these results support a genotype–sex interaction effect, as hormonal regulation with 17β-estradiol partially restored the low efflux observed in the S1731C male proband but had no significant effect on the efflux of a wild-type control.
Rare Missense Variant in LPL
The identified LPL variant rs118204060 is present in the dbSNP13214 and The 1000 Genomes Project14 data with an unknown frequency. The rs118204060 located in exon 5 changes a conserved amino acid from proline to leucine (P234L). This variant was initially identified in familial chylomicronemia and was reported as P207L27–29 because of the differences in genome builds. Upon sequence comparisons, we confirmed that they are indeed the same variant.
Investigation of the ABCA1 and LPL Variants in the Entire Family
We examined the pedigree members for cooccurrence of nonsynonymous ABCA1 and LPL variants with low HDL-C. By stratifying individuals by their HDL percentiles, we can see that all the affected family members with HDL-C less than the fifth percentile carry a risk allele for either 1 or both of the variants (3 P234L, 11 S1731C, and 8 P234L/S1731C), except in 1 separate branch of the extended family in which the low HDL-C traits appears to be inherited from the affected spouse’s side (Figure 1). No family member with an HDL-C value greater than the fifth percentile had both LPL and ABCA1 variants. Furthermore, no family member with HDL-C values greater than the 15th percentile had the LPL variant. Seven subjects (1 male and 6 females) had the ABCA1 variant with the HDL-C percentile of 22% for the male and with an average HDL-percentile of 35% for the 6 females (Figure 1). Two of the 3 exome sequenced subjects were heterozygous for both variants and 1 was homozygous for the ABCA1 variant and heterozygous for the LPL variant. There were 4 homozygous subjects for the ABCA1 variant, 2 of which were also heterozygous for the LPL variant, whereas the LPL variant was heterozygous in all 14 family members it was observed (Figure 1). Thus, a heterozygous, milder form of LPL deficiency exists in this family. Accordingly, the LPL variant P234L is also associated with elevated levels of TGs (1.65±0.27 mmol/L, P=6.14×10–3). In addition, we observed that the subjects with both variants have a lower HDL-C than the subjects with only 1 variant, and that the subjects heterozygous or homozygous for the ABCA1 variant do not differ in the HDL-C levels (Table 4).
Explained Variance and Heritability
We estimated that the effect of the ABCA1 and LPL variants on continuous HDL-C measurements in the extended family is −0.17±0.08 mmol/L (P=0.025) and −0.27±0.09 mmol/L (P=0.006), respectively. Together, these 2 variants explain 60% of the genetic variance in this family and 26% of the total (genetic + environment) variance in this family, which amounts to 46% of the heritability explained as assessed in a measured genotype analysis.18 We also repeated the analysis while excluding the 3 affected subjects that were exome sequenced to reduce the potential for ascertainment bias. In this analysis, the effect sizes of ABCA1 and LPL remained the same (−0.18±0.08 and −0.27±0.10 mmol/L, respectively), and the additive and total variance explained were 50% and 24%, respectively, with 34% of the heritably explained. Importantly, if the subfamily with the bilineal introduction of the low HDL-C trait through the affected spouse is excluded from these analyses, virtually all of the additive variance of HDL-C and virtually all of the heritability of HDL-C is explained by the ABCA1 and LPL variants, suggesting that the ABCA1 and LPL variants can explain the low HDL-C in the non-bilineal part of the extended family.
Genome-Wide Linkage Analyses
To further investigate that we did not miss a major susceptibility variant, we performed a whole-genome 2-point linkage analysis for low HDL-C using a dominant mode of inheritance. We first estimated using the SLINK simulation program30 that under the assumption of homogeneity, the maximum lod score this family can provide is 4.34. However, none of the actual 553 microsatellite markers reached this lod score, most probably due to the existence of multiple low HDL-C variants in the family (ie, heterogeneity). In more detail, no lod scores >3 were observed anywhere in the genome. The only lod score >2.0 was observed on chromosome 21 for marker D21S1255. However, we noticed that this signal on chr 21 seems to arise from the bilineal branch (Figure 1), as the signal diminishes to lod score of 0.6 when we excluded this subfamily from the analysis and increases to 2.5 when we analyzed this bilineal branch of the family alone. Hence the genome-wide linkage data suggest that there might be another susceptibility variant on chr 21 that accounts for the low HDL-C in the bilineal subfamily branch of the extended family. However, because we did not sequence any family members from this branch, none of the 3459 filtered-out variants would be good candidates. Importantly, we observed lod scores >1 near the LPL and ABCA1 genes (lod scores of 1.62 and 1.28 10.3 Mb and 5.8 Mb from LPL and ABCA1, respectively). Without the bilineal branch these lod scores increased to 2.14 and 1.45, respectively.
Genotype by Sex Interaction
The effect of the ABCA1 S1731C variant on low HDL-C levels appears more profound in the males than in females in the extended family (Table 4). Furthermore, our efflux study also suggested a gene × sex interaction (Figure 2B). Although the frequency of the S1731C variant may be individually too rare for testing genetic interactions (as large sample sizes are necessary), rare variants with large phenotypic effects are collectively common in low HDL-C families.26 We hypothesized that the apparent sex effect may not be restricted to the S1731C allele, but rather it may generally extend to ABCA1 alleles with major phenotypic effects. Thus, to further investigate this intriguing relationship between ABCA1 and sex, we examined the collective effect of multiple rare variants in ABCA1 by sex on HDL-C affection. All in all, 10 additional low HDL-C French-Canadian families with known mutations in ABCA16–8 were included in the sex-interaction analysis using the SOLAR program,19 comprising to a total of 93 males and 107 females. The percentage of mutation carriers was 42% and 53% in males and females, respectively (Figure 3). The S1731C variant was present in 3 of these additional low HDL-C pedigrees,7 and together with the exome sequenced family, the association signal for the main effect of S1731C on low HDL-C status resulted in a P value of 0.008. In all 11 families we observed, as expected, a highly significant main effect for ABCA1 genotypes (P=1×10–09), as well as a significant ABCA1 genotype × sex interaction on the qualitative HDL-C affection (P=0.03) (Figure 3). Furthermore, the interaction effect appeared to be more pronounced when comparing premenopausal women (aged <50 years) to men and postmenopausal women (P=0.003).
By using exome sequencing we identified 2 functional rare variants in the ABCA1 and LPL genes, cosegregating with low HDL-C and explaining a major proportion of the HDL-C variance and heritability in an extended family. We also observed a sex effect for ABCA1 variants, male carriers exhibiting significantly lower HDL-C levels than females. Furthermore, none of the unaffected family members had the LPL variant or both variants. Our study exemplifies how use of exome sequencing was critical to reveal the complex combination of 2 variants of which 1 is less severe in females. Traditional linkage analysis was unable to elucidate this type of complex pattern of variants in this extended family,31 suggesting that many such combinations have been missed in previous linkage analyses of complex traits.
ABCA1 and LPL are major players of lipid metabolism. The ABCA1 is a key protein involved in reverse cholesterol transport that transports cellular cholesterol to lipid-poor acceptor apolipoproteins, such as apolipoproteinA-I.20,32 As a result, the apolipoprotein is released with the extracted phospholipid and cholesterol, forming nascent HDL particles. Mutations disrupting the normal function of ABCA1 result in little or no circulating HDL.33 Previous studies have shown that cell lines with the identified ABCA1 variant S1731C exhibit low levels of protein expression,25 and that cells transfected with the S1731C allele express abundant ABCA1 mRNA but fail to generate significant amounts of ABCA1 protein.25 Furthermore, the cholesterol efflux of S1731C has been shown to be reduced to 12.3% to 68.0% of the wild-type.7,25,26 Here, we observed a ≈40% cholesterol efflux reduction in the proband homozygous for the S1731C variant as compared with a normal control, in line with the earlier findings.7,25,26 In our previous article,7 we showed that 3 heterozygous subjects with the S1731C variant have cholesterol efflux values of 63%, 66%, and 68% of the wild type. Thus, about the same 40% decrease is observed in the heterozygous subjects as in the homozygous subject with the S173C variant, in line with their similar HDL-C values (Table 4). Although the phenotype data suggest that the S1731C variant does not have a gene dose effect, this conclusion warrants additional functional studies in future, because there are only 4 homozygotes in the family 2 of which are also heterozygous for the LPL variant, and furthermore the variant has a large range (12.3%–68.0% of the wild-type) in its effect on the cholesterol efflux.7,25,26
The main function of LPL is to hydrolyze TGs to deliver fatty acids to the tissue. LPL also hydrolyzes very-low-density lipoproteins. Sequence variation in LPL has been reported to be associated with the risk of coronary heart disease, TGs, and HDL-C.3 An efficient LPL function is associated with lower TG and low-density lipoproteins and higher HDL. Regarding the identified P207L variant, individuals with this mutation have reduced HDL particles compared with the control subjects.34 Previous studies have also shown that missense mutations in exon 5 of the LPL gene where the P207L variant resides are the most common cause of LPL deficiency.35,36 Importantly, Ma et al reported that upon site-directed in vitro mutagenesis this variant produces a catalytically inactive LPL protein, which is the cause of the LPL deficiency in the patients.29 Taken together, these previous data, along with our PolyPhen15 and SIFT16 predictions, show that it is highly likely that both identified variants S1731C and P207L affect protein function.
The 2 identified variants, S1731C and P207L, have been reported previously in French-Canadian dyslipidemic individuals but not in normal controls.7,24–29 The S1731C ABCA1 variant was present in 3 French-Canadian dyslipidemic families with low HDL-C levels,7 but not in 528 chromosomes from French-Canadian subjects with normal HDL-C levels.24 It was also absent in 108 French-Canadian subjects with high HDL-C.26 The P207L LPL variant was previously observed in 37 unrelated French-Canadian patients with LPL deficiency with 54 mutant alleles present in that study sample.29 In the same study, the variant was also genotyped in 34 unrelated patients with LPL deficiency from ancestries other than French-Canadian. Only 1 German patient was found to be heterozygous for the risk allele. Furthermore, 11 out of 180 French-Canadian hyperlipidemic cases were heterozygous for the P207L variant, whereas none of the 170 normolipidemic controls had the P207L variant.29
It is important to study the effect of sex on lipid traits to better understand the sex-specific differences in incidence of dyslipidemia and cardiovascular disease. The results of an earlier study demonstrated that ABCA1 has a sex-specific effect, because elevated levels of ABCA1 were observed in females,37 which is in line with the higher HDL-C levels and the lower risk of females for coronary artery disease. In this study, we observed that functional mutations in ABCA1 affecting the cholesterol efflux7,25,26 have a larger effect on HDL-C levels in male than female carriers of these variants. It is possible that the observed genotype–sex interaction results from the previously observed sex differences in ABCA1 expression levels,37 because if males have lower baseline levels of ABCA1, the effect of the mutations could be even more profound in males. These interesting sex-specific mechanisms of ABCA1 may involve hormonal regulation of ABCA1, a hypothesis supported by our efflux experiment (Figure 2B), demonstrating that exposure of fibroblasts of a male proband with the ABCA1 S1731C variant to increasing concentrations of 17β-estradiol led to a significantly increased efflux in the male proband with the ABCA1 variant. These intriguing findings warrant further investigation in future studies.
Our results demonstrate that 2 relatively rare functional ABCA1 and LPL variants contribute to the risk of low HDL-C in a unique combination involving a sex-effect in an extended family. We first identified a set of variants by filtering the variants shared by the exome-sequenced affected family members for variant type, frequency, and functional predictions. Because filtering has limitations caused by heterogeneity of complex traits,4 we then used the extended family structure for statistical analysis exploring how much of the trait variance and heritability the 2 key ABCA1 and LPL variants explain. Thus, our study highlights the fact that the filtering strategy used in exome studies of Mendelian disorders4 is not directly applicable for complex disorders, and that new methodologies that incorporate multiple susceptibility variants within a family are warranted. Because the 2 variants explain a major part of the variance in HDL-C and are shown to be functional,7,25,26,29 they represent the key underlying HDL-C variants in this family, though other rare and common variants are likely to explain the remaining portion of the variance. Because both ABCA1 and LPL are known to affect HDL-C, our study did not reveal a novel HDL gene. However, our study does highlight the importance of exome sequencing of dyslipidemic families, because traditional linkage or haplotype analysis cannot detect complex segregation of several functional rare variants due to the inherent parametric restrictions of linkage analysis. This type of underlying biological complexity must have contributed to the low lod scores and weak success of linkage analysis in gene identification of complex lipid traits. In this study, we demonstrate for how family-based exome sequencing can successfully identify multiple rare variants to be followed up utilizing the effective cosegregation information available in extended dyslipidemic families. To the best of our knowledge, our study is the first described example of 2 functional rare variants conferring the susceptibility to low HDL-C in an extended family.
We thank the family members who participated in this study. We also thank Cindy Montes and UCLA core facilities for laboratory technical assistance.
Sources of Funding
This research was supported by the grants HL095056 and HL-28481 from the National Institutes of Health (Drs Pajukanta and Sinsheimer); grants MOP 97752 from the Canadian Institutes of Health Research (CIHR) and from the Heart and Stroke Foundation of Canada (Dr Genest); and by the American Heart Association grant 11POST7380028 (Dr Linga Reddy).
The online-only Data Supplement is available at http://circgenetics.ahajournals.org/lookup/suppl/doi:10.1161/CIRCGENETICS.112.963264/-/DC1.
- Received January 4, 2012.
- Accepted July 23, 2012.
- © 2012 American Heart Association, Inc.
- Dastani Z,
- Quiogue L,
- Plaisier C,
- Engert JC,
- Marcil M,
- Genest J,
- et al
- Marcil M,
- Yu L,
- Krimbou L,
- Boucher B,
- Oram JF,
- Cohn JS Jr.,
- et al
- Li H,
- Handsaker B,
- Wysoker A,
- Fennell T,
- Ruan J,
- Homer N,
- et al
- Wang K,
- Li M,
- Hakonarson H
- Ramensky V,
- Bork P,
- Sunyaev S
- Ng PC,
- Henikoff S
- Marcil M,
- Bissonnette R,
- Vincent J,
- Krimbou L,
- Genest J
- Denis M,
- Haidar B,
- Marcil M,
- Bouvier M,
- Krimbou L,
- Genest J Jr.
- Weissglas-Volkov D,
- Pajukanta P
- Clee SM,
- Zwinderman AH,
- Engert JC,
- Zwarts KY,
- Molhuizen HO,
- Roomp K,
- et al
- Cohen JC,
- Kiss RS,
- Pertsemlidis A,
- Marcel YL,
- McPherson R,
- Hobbs HH
- Brisson D,
- Ledoux K,
- Bossé Y,
- St-Pierre J,
- Julien P,
- Perron P,
- et al
- Weeks DE
- Oram JF,
- Lawn RM
- Peterson J,
- Ayyobi AF,
- Ma Y,
- Henderson H,
- Reina M,
- Deeb SS,
- et al
It is important to elucidate the genetic background of low levels of high-density lipoprotein cholesterol (HDL-C), because low serum HDL-C predisposes to coronary artery disease, the major cause of mortality. Recent genome-wide association studies have successfully identified common DNA sequence variants for HDL-C. However, these common variants explain only a small proportion of variation in the HDL-C levels, indicating that currently still undiscovered rare variants contribute to HDL-C levels as well. Although exome sequencing has been successfully applied to discover rare mutations for Mendelian disorders, not much is known about its feasibility to identify genes for complex traits. We hypothesized that exome sequencing can be used to identify rare variants in unexplained familial cases with low HDL-C because the family-based approach allows the examination of the co-occurrence patterns in the entire family. We sequenced the exomes of subjects with low HDL-C from a large multigenerational family and identified 2 rare nonsynonymous variants in ABCA1 and lipoprotein lipase. The family members carrying either variant had low HDL-C levels, and the individuals carrying both variants had the lowest HDL-C values. Our data also demonstrated a sex effect for the ABCA1 variant. This type of complex combination of 2 rare variants causing low HDL-C in an extended family would not have been identified using traditional linkage analysis, emphasizing the conclusion that exome sequencing can effectively discover rare variants in familial unexplained forms of complex lipid traits.