The use of next-generation sequencing in clinical diagnosis of familial hypercholesterolemia

Purpose:Familial hypercholesterolemia is a common Mendelian disorder associated with early-onset coronary heart disease that can be treated by cholesterol-lowering drugs. The majority of cases in the United Kingdom are currently without a molecular diagnosis, which is partly due to the cost and time associated with standard screening techniques. The main purpose of this study was to test the sensitivity and specificity of two next-generation sequencing protocols for genetic diagnosis of familial hypercholesterolemia.Methods:Libraries were prepared for next-generation sequencing by two target enrichment protocols; one using the SureSelect Target Enrichment System and the other using the PCR-based Access Array platform.Results:In the validation cohort, both protocols showed 100% specificity, whereas the sensitivity for short variant detection was 100% for the SureSelect Target Enrichment and 98% for the Access Array protocol. Large deletions/duplications were only detected using the SureSelect Target Enrichment protocol. In the prospective cohort, the mutation detection rate using the Access Array was highest in patients with clinically definite familial hypercholesterolemia (67%), followed by patients with possible familial hypercholesterolemia (26%).Conclusion:We have shown the potential of target enrichment methods combined with next-generation sequencing for molecular diagnosis of familial hypercholesterolemia. Adopting these assays for patients with suspected familial hypercholesterolemia could improve cost-effectiveness and increase the overall number of patients with a molecular diagnosis.Genet Med 15 12, 948–957.Genetics in Medicine (2013); 15 12, 948–957. doi:10.1038/gim.2013.55

Original research article specific mutations and populations and therefore do not detect less common mutations or any novel mutations. 10,11 During the past few years, high-throughput next-generation sequencing (NGS)-based methods have become available for DNA analysis. They have not only proven successful in new disease gene identification, but following the availability of economic "benchtop" sequencers, they have become more easily applicable to targeted diagnostic sequencing. The combination of high-throughput and relatively small DNA target selection allows for many genes and samples to be processed simultaneously, making it an attractive solution for the processing of large sample numbers in a diagnostic laboratory. Moreover, the flexible design of these assays enables molecular screening to be extended to other relevant genes and polymorphisms that are not covered by routine tests.
In this study, two target enrichment protocols, the hybridization-based SureSelect Target Enrichment System and the PCR-based Access Array System platform, followed by NGS sequencing, have been validated for the detection of FH-causing mutations. We show excellent performance for both approaches and discuss their potential for clinical screening programs and discovery of novel FH-causing genes.

Patients
We studied two groups of patients, a validation group including patients with known molecular diagnosis and a prospective cohort of previously unscreened FH patients. In the validation group, DNA from 104 patients from the Hammersmith Hospital Lipid Clinic were analyzed. All 104 samples had previously undergone complete or partial screening of LDLR-coding regions and exon/intron boundaries, p.Asp374Tyr in PCSK9, and p.Arg3527Gln in APOB exon 26 using Sanger sequencing. In addition, MLPA was used to detect large deletions and duplications in a subset (>80%) of the patients in the validation cohort, as described in the study by Tosi et al. 9 Forty-five samples were known to carry heterozygous FH-causing point mutations or small insertions/deletions (<20 bp) in LDLR, PCSK9, or APOB, and six had large heterozygous insertions/ deletions in LDLR. One sample was homozygous for LDLR p.Gln384Pro, and one sample was a compound heterozygote with one missense mutation and one large deletion in LDLR. The remaining 51 samples in the validation cohort had negative molecular diagnoses. Of the 104 samples, 29 were processed using the SureSelect Target Enrichment System alone, 42 using the Access Array System microfluidic platform alone, and the remaining 33 using both platforms. In the prospective cohort, 84 consecutive patients referred for molecular testing by the Hammersmith Hospital Lipid Clinic over a period of 1 year were studied by the PCR-based Access Array System microfluidic platform. Fifteen of the 84 samples were also analyzed using the SureSelect Target Enrichment System. One individual was referred with suspected homozygous FH. All the remaining 83 were suspected to be heterozygous. Six patients had a diagnosis of definite FH based on Simon Broome criteria. The median highest total cholesterol in this group was 9.7 mmol/l (minimum: 7.9 mmol/l; maximum: 15.9 mmol/l). Sixty-five individuals had a diagnosis of possible FH, with a median highest total cholesterol level of 8.7 mmol/l (minimum: 5.2 mmol/l; maximum: 13 mmol/l), and 13 patients did not fulfill Simon Broome criteria. The median highest total cholesterol in this group was 8.3 mmol/l (minimum: 6.3 mmol/l; maximum: 9.6 mmol/l). The mean age at measurement was 43 years (minimum: 2 years; maximum: 68 years).
The study was approved under ethics committee references REC2002/6451 and REC11/LO/0883. All patients provided informed consent. dna extraction DNA was extracted from blood using a standard phenol-chloroform protocol or the Maxwell 16 system (Promega, Madison, WI) or from saliva samples using the Oragene (Genotek, Ottawa, Ontario, Canada) protocol. Both protocols followed the manufacturer's recommendations.

Rna extraction and reverse-transcription PCR
Blood for RNA extraction was collected in Tempus blood RNA tubes (Applied Biosystems, Foster City, CA) and extracted using the Paxgene blood RNA kit (Qiagen, Hilden, Germany) following the manufacturer's protocol. Reverse transcription was performed using the iScript cDNA synthesis kit (BioRad, Hercules, CA). Primers (forward: TCGAGTTCCACTGCCTAAGTG and reverse: GTTGTTGTCCAAGCATTCGTT) were designed to amplify exons 4-7 of LDLR.
The genomic coordinates of the 18 targeted genes were determined using the March 2006 build (NCBI36/hg18) of the human genome in the Ensembl genome browser. 12 The density of bait tiling was fivefold, and the baits were allowed to overlap into repeat regions by 30 bp. The total targeted DNA length was 399 kb. All libraries were generated from sheared DNA (Covaris, Woburn, MA) with an average insert size of 200 bp following the SureSelect Target Enrichment System XT protocol for Illumina multiplexed sequencing version 1.2 (Illumina, San Diego, CA). After dilution to 2 nmol/l, up to 30 libraries were pooled and sequenced on one lane on the HiSeq2000 platform (Illumina) to generate 2 × 100 bp paired-end reads.

PCR-based access array system
The design included 43 amplicons covering all exons of LDLR, with the majority of the coding sequence covered by more than one overlapping fragment. Amplicons were also designed to cover exons 2, 4, 7, and 9 of PCSK9, containing the most common gain-of-function pathogenic mutations; APOB (one amplicon covering the most common familial defective apolipoprotein B-100 mutation, p.Arg3527Gln); APOE (one amplicon covering the APOE E2 variant site, rs7412); and SLCO1B1 (one amplicon covering rs4149056, the myopathy-associated variant). The average amplicon length was 184 bp, with 57% GC content. Primer sequences are shown in Supplementary Table S1 online. Samples were processed using the Access Array System (Fluidigm, South San Francisco, CA) according to the manufacturer's 4-Primer amplicon tagging protocol generating paired-end libraries for Illumina sequencing. Purified pooled products of 47 barcoded samples and one negative control were sequenced in one run using the MiSeq sequencer (Illumina). One amplicon, a GC-rich (75%) amplicon in APOE targeting rs7412, failed to amplify and was excluded from subsequent data analysis.

data analysis
Sequences were mapped to the GRCh37/hg19 human reference sequence using Burrows Wheeler Aligner v0.6.1. 13 PCR duplicate reads were removed from SureSelect Target Enrichment System data (Picard tools v1.35; http://picard.sourceforge.net). Sequence reads from both data sets were further processed, and variants were called using GATK v1.0.6 14 with hard filtering options. Variant annotation was carried out with Ensembl's Variant Effect Predictor tool 12 and was based on the transcripts ENST00000558518 (LDLR), ENST00000233242 (APOB), and ENST00000302118 (PCSK9). The annotation included Sorting Intolerant From Tolerant, Condel, and PolyPhen. Conservation scores (Genomic Evolutionary Rate Profiling scores) were obtained from Ensembl 12 version 68. To identify potentially pathogenic mutations, all nonsynonymous, splice site, frameshift, and truncating mutations were examined and compared with an FH locus-specific database 6,15 as a guide to interpretation of variant pathogenicity. Synonymous and intronic variants located outside exon/intron boundaries as well as single-nucleotide polymorphisms with minor allele frequency >1% in the International HapMap Project 16 or the 1000 Genomes Project 17 were excluded from further analysis. All single-nucleotide variants and short insertions/deletions that were potentially disease causing were verified by conventional Sanger sequencing. Novel variants were followed up by segregation analysis wherever possible.

Copy-number variant analysis
Copy-number variant (CNV) analysis from NGS data was performed for the samples sequenced using targeted capture. A read depth-based method, 18 as implemented in R package ExomeDepth, was used to identify deletions and duplications spanning at least one exon. Each sequencing batch of samples was processed separately to increase the quality of a reference set for each sample and therefore to maximize the power to detect CNVs. Read depth was assessed for each exon in the target region, and the ratio of expected and observed read count was obtained, as well as a Bayes factor for the CNV calls, as implemented in the ExomeDepth method.
As an independent method of CNV analysis, MLPA was performed using the kit LDLR-P062 (MRC-Holland, Amsterdam, The Netherlands) following the manufacturer's protocol. The novel exon 16 deletion was confirmed by PCR. The previously described large deletions and duplications are detailed in the study by Tosi et al. 9

statistical analysis
The sensitivity of an assay was defined as the percentage of pathogenic mutations correctly identified with respect to previous or new Sanger sequencing and MLPA. The specificity is defined as the percentage of mutation-negative samples correctly identified as negative with respect to previous or new Sanger sequencing and MLPA.

Myopathy-associated variant in SLCO1B1
Genotypes of the SLCO1B1 myopathy-associated variant rs4149056 were scored from the Access Array System microfluidic platform and SureSelect Targeted Enrichment System data in all patients, and any history of adverse effects was obtained by review of medical records. Side effects were defined biochemically (transaminase or creatine kinase levels more than three times the upper limit of the normal range) or symptomatically for myalgia and other side effects that coincided temporally with statin treatment. Tests for deviation from Hardy-Weinberg equilibrium and association tests were performed using the DeFinetti software (http://ihg.gsf.de/cgi-bin/hw/ hwa1.pl).

Validation study
Custom SureSelect Target Enrichment System. To validate the assay, DNA was analyzed from 62 previously screened individuals who either had a confirmed molecular diagnosis of FH (n = 28) or were mutation-negative (n = 34). An average of 330 Gb (191-381 Gb) of sequence was obtained per sample, with an average coverage of ×827. Overall, 64% of total mapped reads aligned within the target region, and 99.8 and 98.8% of nucleotides were covered at ×4 and ×25, respectively. The insufficiently covered regions are consistent among runs and were found mostly outside coding regions. Hybridization-based capture is known to target GC-rich regions poorly; however, all regions containing known FH-causing mutations were covered sufficiently (more than ×25) for confident variant calling. The initial analysis of the sequencing results was carried out by a researcher blinded to the gene and mutation details for each sample. All 20 heterozygous and 1 homozygous short pathogenic mutations, including point mutations and insertions/ deletions of <15 bp, were detected ( Table 1). In addition, one compound heterozygote and six large insertions/deletions in the LDLR gene were detected, resulting in 100% sensitivity for this assay. No false positives were detected, and specificity for this assay was also 100%. A further three mutations were identified in LDLR in patients without a previous molecular diagnosis. One was a known single-nucleotide substitution, p.Asp227Glu, and two were large deletions (deletion of exon 16 and deletion of promoter/exon 1). In addition, a novel variant was identified in the LDLR promoter (c. -227G>T; Genomic Evolutionary Rate Profiling score: 3.29; located in the highly conserved footprint 1 site), and two rare variants were identified in APOB in patients without a previous molecular diagnosis.

Prospective cohort
To test the feasibility of an NGS-based mutation screening in a clinical setting, a consecutive cohort of 84 unrelated patients referred for genetic screening over a single year was analyzed. Sixty-nine of these patients were screened using the Access Array System assay alone, and 15 were processed using both protocols. Sequencing and coverage metrics are provided in Supplementary Table S2 online. In total, 23 variants that have previously been reported as pathogenic in FH patients were identified in 22 individuals ( Table 2), including 13 single-nucleotide changes, 4 variants predicted to affect splicing, and 6 short insertions/deletions. In addition, two novel variants (p.Ala521Thr and p.Asn316Lysfs*54) were identified that had not been previously reported ( Table 2). There was a good correlation between the results obtained from the Access Array System assay and SureSelect Target Enrichment System, with five of six mutations concordant between the data sets. The discordant variant, p.Lys730Hisfs*48, is the same one that was not detected using the Access Array System protocol in the validation cohort. Despite the presence of the same variant, these patients in the validation and prospective cohort were not known to be related. No additional mutations were identified using MLPA.
Of the two previously unreported variants ( Table 2)  Two (c.817+1G>A and c.941-4G>A) of the four splice site variants identified ( Table 2) were previously detected in FH patients, 6 but their pathogenicity had not been experimentally validated before. An RNA sample was available for one of them (patient ID: 35) and reverse-transcriptase PCR in this patient showed that the c.817+1G>A variant disrupted correct splicing (data not shown). Alternative splice products were generated by partial intron retention that led to truncation of the protein, by exon skipping that caused the loss of exon 5 sequence, and by a predicted frameshift of the remaining protein. Of the remaining variants previously identified in FH patients, one missense variant, p.Val827Ile, was conservative, without a published proof of pathogenicity and therefore was considered to be of unknown significance. Two children of this index case were available for segregation analysis, and the results showed that an unaffected daughter (total cholesterol: 4.6 mmol/l; age: 27 years) of the proband had inherited the p.Val827Ile variant.

Mutation detection rates
The prospective cohort included six patients with definite FH as defined by the Simon Broome criteria, 65 with possible FH, and Original research article 13 hypercholesterolemic patients not fulfilling Simon Broome criteria for FH. The highest detection rate of clearly pathogenic mutations was in the group of patients with a definite FH diagnosis (4 of 6, 67%) followed by the group with a diagnosis of possible FH (17 of 65, 26%). One mutation (1 of 13, 8%) was identified among the 13 hypercholesterolemic patients who did not fulfill FH diagnostic criteria.

Variants in genes involved in cholesterol metabolism
In addition to the three known FH-causing genes and SLCO1B1, the SureSelect Target Enrichment System assay included 13 genes that reside on cholesterol metabolism pathways. Rare variants identified in these genes, particularly in hypercholesterolemic individuals who had been screened negative for mutations in LDLR, PCSK9, and APOB, are potentially responsible for patients' raised cholesterol levels. Rare variants (minor allele frequency < 0.01) identified in such individuals are listed in Supplementary Table S3 online. In patients with no previously known molecular diagnosis, seven rare nonsynonymous variants were found, of which six were not predicted to be functionally significant by Sorting Intolerant From Tolerant and PolyPhen. The single variant that was most likely to be of functional significance and therefore potentially pathogenic, was p.Val809Met in SREBF1. However, this was excluded from further analysis because it did not segregate with the phenotype in the family (data not shown).

SLCO1B1 myopathy-associated variant
Data on exposure to statins were available for 149 patients. Twenty-seven individuals were heterozygous for SLCO1B1 rs4149056, and two were homozygous. One patient (a heterozygote) had suffered severe biopsy-proven hepatitis while taking

disCUssion
We investigated the sensitivity and specificity of two target enrichment protocols, combined with NGS, for the detection of disease-causing mutations in patients with proven or suspected FH. Using the SureSelect Target Enrichment System and the HiSeq 2000 system, 98.8% of targeted regions were covered at more than ×25 as compared with the PCR-based Access Array System followed by MiSeq sequencing, for which on average 10% of nucleotides failed to reach this coverage. Regions with low coverage were found outside the coding sequence of LDLR and mutation hotspots in APOB and PCSK9 and therefore did not affect the overall success of mutation detection. Both techniques showed high sensitivity and specificity for the detection of single base substitutions and short insertions/deletions. In the validation part of the study, 100% of previously detected mutations were correctly identified using the SureSelect/HiSeq protocol. In comparison, the Access Array System/MiSeq approach led to correct identification of 98% of all variants detected previously by Sanger sequencing, although large insertions/deletions could not be detected (see below). The single variant that was not detected by the PCR-based Access Array System protocol was an 11 bp deletion (p.Lys730Hisfs*48) that overlapped a primer site. However, the recent availability of longer sequencing reads, currently up to 250 bp when using MiSeq, will be anticipated in most cases to eliminate the need to have primer sites within coding regions or allow design modification to include longer overlaps between amplicons. Three mutations in the LDLR coding sequence (one singlenucleotide variant, p.Asp227Glu, and two large deletions) and one LDLR promoter variant (c.-227G>T) were identified in patients who were previously classified as negative for mutations. Inspection of the previous laboratory database indicated that the p.Asp227Glu variant and the two large deletions had not been detected because of incomplete capillary sequencing in the case of the p.Asp227Glu and promoter variants and absence of MLPA data in the case of the two deletions.
A small proportion of FH cases are caused by large deletions or duplications of LDLR. 19 The current standard screening is based on MLPA, a technique that is highly reliable but that is a costly and time-consuming addition to sequencing. The detection of large variants from NGS data has been shown previously, but its use in FH diagnostics has not yet been investigated. 20,21 Here, the combination of SureSelect Target Enrichment System/ HiSeq and data analysis using ExomeDepth software 18 led to correct identification of all eight large deletions and one large duplication. This shows the potential of using hybrid capture for the detection of both short and large sequence variants in FH. The read depth from the PCR-based Access Array System assay was not used for this analysis method because read depth in this study did not correlate with exon deletions, although it has recently been shown that amplicon multiplex PCR can be optimized for CNV detection. 22 When testing a new molecular workflow for routine diagnostics, in addition to sensitivity and specificity, it is important to consider the cost and time efficiency of the protocol. The SureSelect Target Enrichment System protocol allows for comprehensive coverage of targeted regions (current custom design up to 24 Mb), whereas the Access Array System protocol is limited to 480 amplicons with a maximum length of 400 bp when sequenced using the latest MiSeq system. Sequence capture also allows the detection of all types of variants, including large deletions and duplications, which was not possible in this study using the PCR-based Access Array System. On the other hand, the Access Array System is considerably cheaper to run, 20,23 with reagent costs approximately 10-fold lower than those for the SureSelect Target Enrichment System protocol, 23 and, in addition, the library preparation turnaround time is shorter. In our hands, 96 samples can be processed within a day using the Access Array System protocol as compared with at least 3 days needed for the SureSelect Target Enrichment System insolution capture. Initial DNA requirements of the Access Array System/MiSeq protocol are also low (50 ng vs. 1 ug), and most of the process is automated, reducing the risks of contamination and human error to a minimum. Recently, Hollants et al. 24 Original research article published a protocol for FH mutation detection that was also based on the Access Array System, but the amplicons were sequenced using the pyrosequencing-based GS-FLX (Roche, Branford, CT) platform. They successfully identified all short variants but, as with our study, could not detect large rearrangements. Our Access Array System protocol is here combined with the MiSeq sequencing platform (a benchtop personal NGS system specifically developed to suit the needs of the routine laboratory setting) and offers similar quality to GS-FLX but at a lower sequencing cost. 25 The use of the MiSeq system is also reported to limit errors of variant calls within homopolymer regions that are known to occur in pyrosequencing. 26

novel variants identified in the validation study
The LDLR promoter variant c.-227G>T has not previously been reported in FH patients, or in other populations, but was investigated previously as part of a study delineating the conserved footprint 1 site. A luciferase assay showed that the c. -227G>T variant had ~75% transcription levels as compared with the wild-type site. 27 Although a 25% reduction is not a definitive decrease in promoter activity, it cannot be excluded that this change is sufficient to cause raised cholesterol levels in this patient. Family members were unavailable for segregation analysis, and we therefore classify this variant as being of uncertain significance pending further functional and segregation data. The remaining two novel variants in APOB, rs12714097 and rs72653087, were located outside LDLR-binding sites, regions currently not associated with hypercholesterolemia, and therefore their pathogenicity remains to be elucidated.

Genotype-phenotype correlation in the prospective cohort
In the United Kingdom, FH diagnosis is made based on the Simon Broome criteria that identify patients with definite or possible FH. Following the National Institute for Health and Clinical Excellence guidelines, 6 genetic testing is recommended for all suspected index cases and should be followed by family cascade screening if a mutation is found. Our study focused on a group of consecutive patients referred for genetic screening by a local lipid clinic over the course of 1 year. To assess clinical significance, all rare variants (minor allele frequency < 0.01) that were identified in LDLR, APOB, or PCSK9 were compared against locus-specific FH databases. 6,15 All nonsense variants and insertion/deletion variants that have been found previously among FH patients were also considered pathogenic. Missense variants were considered to be pathogenic if they were predicted by Sorting Intolerant From Tolerant and PolyPhen to be deleterious and probably damaging, respectively, and were previously identified in FH patients. On the basis of these criteria, a missense change, p.Val827Ile, was included among the list of pathogenic variants, although the conservative valine to isoleucine substitution may not affect protein function. This residue is the third amino acid within the internalization signal for LDLR, NPVY, a position generally not necessary for correct internalization. 28 In addition, our segregation data suggest that this variant may not be disease causing, and we therefore classified this variant as being of unknown significance. Of the four splice site variants, three (c.313+1G>A, c.190+4A>T, and c.817+1G>A) were shown experimentally to disrupt splicing, two in published data 29,30 and one within this study, and are therefore listed as pathogenic. The remaining splice variant, c.941-4G>A, was described previously in FH patients, but its pathogenicity remains to be experimentally confirmed. An RNA sample from this individual was not available for confirmation within this study.
The mutation detection rates among patients with definite (66%) and possible (26%) FH are comparable to those of previously published studies. 7 One mutation was identified among 13 patients who did not fulfill Simon Broome criteria. This suggests that a number of patients with high cholesterol in the general population who do not fulfill formal clinical criteria for FH diagnosis may have FH-causing mutations. Such patients may therefore be diagnosed by the assays developed here, and, because these mutations are dominantly inherited, family cascade screening would probably identify the same mutations in 50% of their first-degree relatives.
The types of mutation identified here reflect the distribution of variants published in the LDLR locus-specific database, 15 with exonic substitutions (56%) being the most common, followed by short insertions/deletions (28%). Most of the mutations identified here were unique, with only 20% common among UK FH patients that would be identified using the Elucigene FH20 commercial ARMS kit. In our prospective study, we did not identify any mutations in PCSK9 or large rearrangements of LDLR. These variants are generally rare among patients with FH, 7 and we therefore consider it unlikely that our results are biased in any way.
Recently, the importance of screening additional regions of APOB has been highlighted by Motazacker et al., 31 who identified two patients with novel mutations located outside the commonly sequenced region of APOB exon 26. Screening of entire coding regions instead of focusing on mutation hot spots is therefore likely to enhance the discovery of FH-causing mutations. The format of next-generation assays is flexible and can be readily extended to include full coding regions of APOB and PCSK9, as well as coding regions of other medically relevant genes, such as APOE and SLCO1B1, which currently need to be genotyped separately.
Our current Access Array System covers the same regions that would be screened using conventional Sanger sequencing, and therefore the number of variants identified would be no larger than that identified after Sanger sequencing. The SureSelect Target Enrichment System design also included genes in cholesterol metabolism pathways and regions of APOB and PCSK9 outside those known to be causative of FH. To establish the pathogenicity of such variants, further in silico and experimental analyses will be required. For detection of variants outside known FH-causing genes or gene regions, we focused on variants that are rare (minor allele frequency < 0.01), are probably pathogenic (i.e., located in coding regions, promoters, and exon/intron boundaries), and are present only in individuals without an existing molecular diagnosis of FH. Using these filters, we were left with a manageable list of variants (<10) that we decided to follow up. When screening larger regions in more patients, the identification of large numbers of variants of unknown significance will necessitate further extensive family segregation and functional assays. However, the increase of data available in public databases and improvements of bioinformatics tools will allow more efficient filtering of variants than is available at present.

statin-induced myopathy
Previous published data on larger cohorts show an odds ratio of 4.5 for association of rs4149056 with myopathy, 32 and on that basis, guidelines have recommended avoiding higher doses of simvastatin in heterozygotes or homozygotes for rs4149056. 33 Association testing restricted to more severe side effects (transaminases or creatine kinase greater than three times the upper limit of normal) was not informative in our cohort due to the small sample size. However, less-severe patient-reported adverse effects to statin treatment such as myalgia have an effect on patient adherence to treatment in this high-risk population. Therefore, the association identified here between SLCO1B1 genotype and a wider range of milder adverse effects in this lipid clinic cohort gives an insight into the potential benefits of prospective genotyping before the initiation of statin treatment.
In conclusion, we have shown the potential utility of sequence target enrichment methods in combination with NGS in molecular diagnostics of FH. Owing to the comprehensive coverage, SureSelect Target Enrichment System protocols (either whole-exome or region-specific) may provide the most benefit for studies that aim to identify novel disease-causing genes and for diseases for which the number of genes needing to be screened is very large. In FH diagnostics for which <5 genes need to be analyzed, PCR-based enrichment techniques offer more streamlined protocols that provide high sensitivity for mutation detection and may offer an affordable option for clinical screening of large numbers of patients with suspected FH. If adopted, greater numbers of FH patients and their relatives may potentially benefit from early diagnosis and treatment.

SUPPLEMENTARY MATERIAL
Supplementary material is linked to the online version of the paper at http://www.nature.com/gim