ENIGMA CHEK2gether Project: A Comprehensive Study Identifies Functionally Impaired CHEK2 Germline Missense Variants Associated with Increased Breast Cancer Risk

PURPOSE

Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT).

EXPERIMENTAL DESIGN

We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1-CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case-control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls.

RESULTS

A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)-like (N = 226). We then examined their association with breast cancer risk in the case-control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35-3.41), 1.57 (95% CI, 1.41-1.75), and 1.19 (95% CI, 1.08-1.31), respectively. The meta-analysis of population-specific datasets showed similar results.

CONCLUSIONS

We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

Performance Evaluation of Three DNA Sample Tracking Tools in a Whole Exome Sequencing Workflow

INTRODUCTION

Next-generation sequencing applications are becoming indispensable for clinical diagnostics. These experiments require numerous wet- and dry-laboratory steps, each one increasing the probability of a sample swap or contamination. Therefore, identity confirmation at the end of the process is recommended to ensure the right data are used for each patient.

METHODS

We tested three commercially available, single nucleotide polymorphism (SNP)-based sample tracking kits in a diagnostic workflow to evaluate their ease of use and performance. The coverage uniformity, on-target specificity, sample identification, and genotyping performance were determined to assess the reliability and cost effectiveness of each kit.

RESULTS AND DISCUSSION

Hands-on time and manual steps are almost identical for the kits from pxlence and Nimagen. The Swift kit has an extra purification step, making it the longest and most demanding protocol. Furthermore, the Swift kit failed to correctly genotype 26 of the 46 samples. The Nimagen kit identified all but one sample and the pxlence kit unambiguously identified all samples, making it the most reliable and robust kit of this evaluation. The Nimagen kit showed poor on-target mapping rates, resulting in deeper sequencing needs and higher sequencing costs compared with the other two kits.

CONCLUSION

Our conclusion is that the Human Sample ID kit from pxlence is the most cost effective of the three tested tools for DNA sample tracking and identification.

Breast and Prostate Cancer Risks for Male BRCA1 and BRCA2 Pathogenic Variant Carriers Using Polygenic Risk Scores

BACKGROUND

Recent population-based female breast cancer and prostate cancer polygenic risk scores (PRS) have been developed. We assessed the associations of these PRS with breast and prostate cancer risks for male BRCA1 and BRCA2 pathogenic variant carriers.

METHODS

483 BRCA1 and 1318 BRCA2 European ancestry male carriers were available from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). A 147-single nucleotide polymorphism (SNP) prostate cancer PRS (PRSPC) and a 313-SNP breast cancer PRS were evaluated. There were 3 versions of the breast cancer PRS, optimized to predict overall (PRSBC), estrogen receptor (ER)-negative (PRSER-), or ER-positive (PRSER+) breast cancer risk.

RESULTS

PRSER+ yielded the strongest association with breast cancer risk. The odds ratios (ORs) per PRSER+ standard deviation estimates were 1.40 (95% confidence interval [CI] =1.07 to 1.83) for BRCA1 and 1.33 (95% CI = 1.16 to 1.52) for BRCA2 carriers. PRSPC was associated with prostate cancer risk for BRCA1 (OR = 1.73, 95% CI = 1.28 to 2.33) and BRCA2 (OR = 1.60, 95% CI = 1.34 to 1.91) carriers. The estimated breast cancer odds ratios were larger after adjusting for female relative breast cancer family history. By age 85 years, for BRCA2 carriers, the breast cancer risk varied from 7.7% to 18.4% and prostate cancer risk from 34.1% to 87.6% between the 5th and 95th percentiles of the PRS distributions.

CONCLUSIONS

Population-based prostate and female breast cancer PRS are associated with a wide range of absolute breast and prostate cancer risks for male BRCA1 and BRCA2 carriers. These findings warrant further investigation aimed at providing personalized cancer risks for male carriers and informing clinical management.

Dealing with Pseudogenes in Molecular Diagnostics in the Next Generation Sequencing Era

Presence of pseudogenes is a dreadful issue in next generation sequencing (NGS), because their contamination can interfere with the detection of variants in the genuine gene and generate false positive and false negative variants.In this chapter we focus on issues related to the application of NGS strategies for analysis of genes with pseudogenes in a clinical setting. The degree to which a pseudogene impacts the ability to accurately detect and map variants in its parent gene depends on the degree of similarity (homology) with the parent gene itself. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences. For target enrichment, we describe advantages and disadvantages of PCR- and capture-based strategies. For mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some examples of genes associated with Mendelian disorders, for which interesting NGS approaches are described to avoid interference with pseudogene sequences.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

ABCA4-associated disease as a model for missing heritability in autosomal recessive disorders: novel noncoding splice, cis-regulatory, structural, and recurrent hypomorphic variants

Purpose

ABCA4-associated disease, a recessive retinal dystrophy, is hallmarked by a large proportion of patients with only one pathogenic ABCA4 variant, suggestive for missing heritability.

Methods

By locus-specific analysis of ABCA4, combined with extensive functional studies, we aimed to unravel the missing alleles in a cohort of 67 patients (p), with one (p = 64) or no (p = 3) identified coding pathogenic variants of ABCA4.

Results

We identified eight pathogenic (deep-)intronic ABCA4 splice variants, of which five are novel and six structural variants, four of which are novel, including two duplications. Together, these variants account for the missing alleles in 40.3% of patients. Furthermore, two novel variants with a putative cis-regulatory effect were identified. The common hypomorphic variant c.5603A>T p.(Asn1868Ile) was found as a candidate second allele in 43.3% of patients. Overall, we have elucidated the missing heritability in 83.6% of our cohort. In addition, we successfully rescued three deep-intronic variants using antisense oligonucleotide (AON)-mediated treatment in HEK 293-T cells and in patient-derived fibroblast cells.

Conclusion

Noncoding pathogenic variants, novel structural variants, and a common hypomorphic allele of the ABCA4 gene explain the majority of unsolved cases with ABCA4-associated disease, rendering this retinopathy a model for missing heritability in autosomal recessive disorders.

BRCA1 and BRCA2 5' noncoding region variants identified in breast cancer patients alter promoter activity and protein binding

The widespread use of next generation sequencing for clinical testing is detecting an escalating number of variants in noncoding regions of the genome. The clinical significance of the majority of these variants is currently unknown, which presents a significant clinical challenge. We have screened over 6,000 early-onset and/or familial breast cancer (BC) cases collected by the ENIGMA consortium for sequence variants in the 5' noncoding regions of BC susceptibility genes BRCA1 and BRCA2, and identified 141 rare variants with global minor allele frequency < 0.01, 76 of which have not been reported previously. Bioinformatic analysis identified a set of 21 variants most likely to impact transcriptional regulation, and luciferase reporter assays detected altered promoter activity for four of these variants. Electrophoretic mobility shift assays demonstrated that three of these altered the binding of proteins to the respective BRCA1 or BRCA2 promoter regions, including NFYA binding to BRCA1:c.-287C>T and PAX5 binding to BRCA2:c.-296C>T. Clinical classification of variants affecting promoter activity, using existing prediction models, found no evidence to suggest that these variants confer a high risk of disease. Further studies are required to determine if such variation may be associated with a moderate or low risk of BC.

Mutational spectrum in a worldwide study of 29,700 families with BRCA1 or BRCA2 mutations

The prevalence and spectrum of germline mutations in BRCA1 and BRCA2 have been reported in single populations, with the majority of reports focused on White in Europe and North America. The Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) has assembled data on 18,435 families with BRCA1 mutations and 11,351 families with BRCA2 mutations ascertained from 69 centers in 49 countries on six continents. This study comprehensively describes the characteristics of the 1,650 unique BRCA1 and 1,731 unique BRCA2 deleterious (disease-associated) mutations identified in the CIMBA database. We observed substantial variation in mutation type and frequency by geographical region and race/ethnicity. In addition to known founder mutations, mutations of relatively high frequency were identified in specific racial/ethnic or geographic groups that may reflect founder mutations and which could be used in targeted (panel) first pass genotyping for specific populations. Knowledge of the population-specific mutational spectrum in BRCA1 and BRCA2 could inform efficient strategies for genetic testing and may justify a more broad-based oncogenetic testing in some populations.

Accurate detection and quantification of epigenetic and genetic second hits in BRCA1 and BRCA2-associated hereditary breast and ovarian cancer reveals multiple co-acting second hits

BACKGROUND

This study characterizes the second hit spectrum in BRCA1 and BRCA2-associated breast and ovarian cancers at both gene loci to investigate if second hit mechanisms are mutually exclusive or able to coincide within the same tumor.

METHODS

Loss of heterozygosity, somatic point mutations and copy number alterations along with promoter methylation were studied in 56 breast and 15 ovarian cancers from BRCA1 and BRCA2 germline mutation carriers. A mathematical methodology was introduced to quantify the tumor cell population carrying a second hit.

RESULTS

Copy neutral LOH was the most prevalent LOH mechanism in this cohort (BC 69%, OC 67%). However, only 36% of BC and 47% of OC showed LOH in all cancerous cells. Somatic intragenic deletions and methylated subclones were also found in combination with (partial) loss of heterozygosity. Unequivocal deleterious somatic point mutations were not identified in this cohort.

CONCLUSION

Different mechanisms inactivating the wild type allele are present within the same tumor sample at various extents. Results indicate that BRCA1/2-linked breast and ovarian cancer cells are predominantly characterized by LOH, but harbor a complex combination of second hits at various frequencies.

Thorough in silico and in vitro cDNA analysis of 21 putative BRCA1 and BRCA2 splice variants and a complex tandem duplication in BRCA2 allowing the identification of activated cryptic splice donor sites in BRCA2 exon 11

For 21 putative BRCA1 and BRCA2 splice site variants, the concordance between mRNA analysis and predictions by in silico programs was evaluated. Aberrant splicing was confirmed for 12 alterations. In silico prediction tools were helpful to determine for which variants cDNA analysis is warranted, however, predictions for variants in the Cartegni consensus region but outside the canonical sites, were less reliable. Learning algorithms like Adaboost and Random Forest outperformed the classical tools. Further validations are warranted prior to implementation of these novel tools in clinical settings. Additionally, we report here for the first time activated cryptic donor sites in the large exon 11 of BRCA2 by evaluating the effect at the cDNA level of a novel tandem duplication (5' breakpoint in intron 4; 3' breakpoint in exon 11) and of a variant disrupting the splice donor site of exon 11 (c.6841+1G > C). Additional sites were predicted, but not activated. These sites warrant further research to increase our knowledge on cis and trans acting factors involved in the conservation of correct transcription of this large exon. This may contribute to adequate design of ASOs (antisense oligonucleotides), an emerging therapy to render cancer cells sensitive to PARP inhibitor and platinum therapies.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Prediction of Breast and Prostate Cancer Risks in Male BRCA1 and BRCA2 Mutation Carriers Using Polygenic Risk Scores

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

Analysis of chromosomal radiosensitivity of healthy BRCA2 mutation carriers and non-carriers in BRCA families with the G2 micronucleus assay

Breast cancer risk drastically increases in individuals with a heterozygous germline BRCA1 or BRCA2 mutation, while it is estimated to equal the population risk for relatives without the familial mutation (non-carriers). The aim of the present study was to use a G2 phase-specific micronucleus assay to investigate whether lymphocytes of healthy BRCA2 mutation carriers are characterized by increased radiosensitivity compared to controls without a family history of breast/ovarian cancer and how this relates to healthy non-carrier relatives. BRCA2 is active in homologous recombination, a DNA damage repair pathway, specifically active in the late S/G2 phase of the cell cycle. We found a significantly increased radiosensitivity in a cohort of healthy BRCA2 mutation carriers compared to individuals without a familial history of breast cancer (P=0.046; Mann-Whitney U test). At the individual level, 50% of healthy BRCA2 mutation carriers showed a radiosensitive phenotype (radiosensitivity score of 1 or 2), whereas 83% of the controls showed no radiosensitivity (P=0.038; one-tailed Fishers exact test). An odds ratio of 5 (95% CI, 1.07–23.47) indicated an association between the BRCA2 mutation and radiosensitivity in healthy mutation carriers. These results indicate the need for the gentle use of ionizing radiation for either diagnostic or therapeutic use in BRCA2 mutation carriers. We detected no increased radiosensitivity in the non-carrier relatives.

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression: identification of a modifier of breast cancer risk at locus 11q22.3

PURPOSE

Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways.

METHODS

Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of ~320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2.

RESULTS

We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 × 10-6). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance.

CONCLUSION

We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND

Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk.

METHOD

We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation.

RESULTS

Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05.

CONCLUSION

This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

Increased chromosomal radiosensitivity in asymptomatic carriers of a heterozygous BRCA1 mutation

Background

Breast cancer risk increases drastically in individuals carrying a germline BRCA1 mutation. The exposure to ionizing radiation for diagnostic or therapeutic purposes of BRCA1 mutation carriers is counterintuitive, since BRCA1 is active in the DNA damage response pathway. The aim of this study was to investigate whether healthy BRCA1 mutations carriers demonstrate an increased radiosensitivity compared with healthy individuals.

Methods

We defined a novel radiosensitivity indicator (RIND) based on two endpoints measured by the G₂ micronucleus assay, reflecting defects in DNA repair and G₂ arrest capacity after exposure to doses of 2 or 4 Gy. We investigated if a correlation between the RIND score and nonsense-mediated decay (NMD) could be established.

Results

We found significantly increased radiosensitivity in the cohort of healthy BRCA1 mutation carriers compared with healthy controls. In addition, our analysis showed a significantly different distribution over the RIND scores (p = 0.034, Fisher’s exact test) for healthy BRCA1 mutation carriers compared with non-carriers: 72 % of mutation carriers showed a radiosensitive phenotype (RIND score 1–4), whereas 72 % of the healthy volunteers showed no radiosensitivity (RIND score 0). Furthermore, 28 % of BRCA1 mutation carriers had a RIND score of 3 or 4 (not observed in control subjects). The radiosensitive phenotype was similar for relatives within several families, but not for unrelated individuals carrying the same mutation. The median RIND score was higher in patients with a mutation leading to a premature termination codon (PTC) located in the central part of the gene than in patients with a germline mutation in the 5′ end of the gene.

Conclusions

We show that BRCA1 mutations are associated with a radiosensitive phenotype related to a compromised DNA repair and G₂ arrest capacity after exposure to either 2 or 4 Gy. Our study confirms that haploinsufficiency is the mechanism involved in radiosensitivity in patients with a PTC allele, but it suggests that further research is needed to evaluate alternative mechanisms for mutations not subjected to NMD.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-016-0709-1) contains supplementary material, which is available to authorized users.

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.

Targeted resequencing and variant validation using pxlence PCR assays

The advent of next-generation sequencing technologies had a profound impact on molecular diagnostics. PCR is a popular method for target enrichment of disease gene panels. Using our proprietary primer-design pipeline, primerXL, we have created almost one million assays covering over 98% of the human exome. Here we describe the assay specification and both in silico and wet-lab validation of a selected set of 2294 assays using both next-generation sequencing and Sanger sequencing. Using a universal PCR protocol without optimization, these assays result in high coverage uniformity and limited non-specific coverage. In addition, data indicates a positive correlation between the predictive in silico specificity score and the amount of assay non-specific coverage.

Flexible, scalable, and efficient targeted resequencing on a benchtop sequencer for variant detection in clinical practice

The release of benchtop next-generation sequencing (NGS) instruments has paved the way to implement the technology in clinical setting. The need for flexible, qualitative, and cost-efficient workflows is high. We used singleplex-PCR for highly efficient target enrichment, allowing us to reach the quality standards set in Sanger sequencing-based diagnostics. For the library preparation, a modified NexteraXT protocol was used, followed by sequencing on a MiSeq instrument. With an innovative pooling strategy, high flexibility, scalability, and cost-efficiency were obtained, independent of the availability of commercial kits. The approach was validated for ∼250 genes associated with monogenic disorders. An overall sensitivity (>99%) similar to Sanger sequencing was observed in combination with a positive predictive value of >98%. The distribution of coverage was highly uniform, guaranteeing a minimal number of gaps to be filled with alternative methods. ISO15189-accreditation was obtained for the workflow. A major asset of the singleplex PCR-based enrichment is that new targets can be easily implemented. Diagnostic laboratories have validated assays available ensuring that the proposed workflow can easily be adopted. Although our platform was optimized for constitutional variant detection of monogenic disease genes, it is now also used as a model for somatic mutation detection in acquired diseases.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

INTRODUCTION

Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

METHODS

We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

RESULTS

We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

CONCLUSIONS

This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Association of type and location of BRCA1 and BRCA2 mutations with risk of breast and ovarian cancer

IMPORTANCE

Limited information about the relationship between specific mutations in BRCA1 or BRCA2 (BRCA1/2) and cancer risk exists.

OBJECTIVE

To identify mutation-specific cancer risks for carriers of BRCA1/2.

DESIGN, SETTING, AND PARTICIPANTS

Observational study of women who were ascertained between 1937 and 2011 (median, 1999) and found to carry disease-associated BRCA1 or BRCA2 mutations. The international sample comprised 19,581 carriers of BRCA1 mutations and 11,900 carriers of BRCA2 mutations from 55 centers in 33 countries on 6 continents. We estimated hazard ratios for breast and ovarian cancer based on mutation type, function, and nucleotide position. We also estimated RHR, the ratio of breast vs ovarian cancer hazard ratios. A value of RHR greater than 1 indicated elevated breast cancer risk; a value of RHR less than 1 indicated elevated ovarian cancer risk.

EXPOSURES

Mutations of BRCA1 or BRCA2.

MAIN OUTCOMES AND MEASURES

Breast and ovarian cancer risks.

RESULTS

Among BRCA1 mutation carriers, 9052 women (46%) were diagnosed with breast cancer, 2317 (12%) with ovarian cancer, 1041 (5%) with breast and ovarian cancer, and 7171 (37%) without cancer. Among BRCA2 mutation carriers, 6180 women (52%) were diagnosed with breast cancer, 682 (6%) with ovarian cancer, 272 (2%) with breast and ovarian cancer, and 4766 (40%) without cancer. In BRCA1, we identified 3 breast cancer cluster regions (BCCRs) located at c.179 to c.505 (BCCR1; RHR = 1.46; 95% CI, 1.22-1.74; P = 2 × 10(-6)), c.4328 to c.4945 (BCCR2; RHR = 1.34; 95% CI, 1.01-1.78; P = .04), and c. 5261 to c.5563 (BCCR2', RHR = 1.38; 95% CI, 1.22-1.55; P = 6 × 10(-9)). We also identified an ovarian cancer cluster region (OCCR) from c.1380 to c.4062 (approximately exon 11) with RHR = 0.62 (95% CI, 0.56-0.70; P = 9 × 10(-17)). In BRCA2, we observed multiple BCCRs spanning c.1 to c.596 (BCCR1; RHR = 1.71; 95% CI, 1.06-2.78; P = .03), c.772 to c.1806 (BCCR1'; RHR = 1.63; 95% CI, 1.10-2.40; P = .01), and c.7394 to c.8904 (BCCR2; RHR = 2.31; 95% CI, 1.69-3.16; P = .00002). We also identified 3 OCCRs: the first (OCCR1) spanned c.3249 to c.5681 that was adjacent to c.5946delT (6174delT; RHR = 0.51; 95% CI, 0.44-0.60; P = 6 × 10(-17)). The second OCCR spanned c.6645 to c.7471 (OCCR2; RHR = 0.57; 95% CI, 0.41-0.80; P = .001). Mutations conferring nonsense-mediated decay were associated with differential breast or ovarian cancer risks and an earlier age of breast cancer diagnosis for both BRCA1 and BRCA2 mutation carriers.

CONCLUSIONS AND RELEVANCE

Breast and ovarian cancer risks varied by type and location of BRCA1/2 mutations. With appropriate validation, these data may have implications for risk assessment and cancer prevention decision making for carriers of BRCA1 and BRCA2 mutations.

Novel FRMD7 Mutations and Genomic Rearrangement Expand the Molecular Pathogenesis of X-Linked Idiopathic Infantile Nystagmus

PURPOSE

Idiopathic infantile nystagmus (IIN; OMIM 31700) with X-linked inheritance is one of the most common forms of infantile nystagmus. Up to date, three X-linked loci have been identified, Xp11.4-p11.3 (calcium/calmodulin-dependent serine protein kinase [CASK]), Xp22 (GPR143), and Xq26-q27 (FRMD7), respectively. Here, we investigated the role of mutations and copy number variations (CNV) of FRMD7 and GPR143 in the molecular pathogenesis of IIN in 49 unrelated Belgian probands.

METHODS

We set up a comprehensive molecular genetic workflow based on Sanger sequencing, targeted next generation sequencing (NGS) and CNV analysis using multiplex ligation-dependent probe amplification (MLPA) for FRMD7 (NM_194277.2) and GPR143 (NM_000273.2).

RESULTS

In 11/49 probands, nine unique FRMD7 changes were found, five of which are novel: frameshift mutation c.2036del, missense mutations c.801C>A and c.875T>C, splice-site mutation c.497+5G>A, and one genomic rearrangement (1.29 Mb deletion) in a syndromic case. Additionally, four known mutations were found: c.70G>A, c.886G>C, c.910C>T, and c.660del. The latter was found in three independent families. In silico predictions and segregation testing of the novel mutations support their pathogenic effect. No GPR143 mutations or CNVs were found in the remainder of the probands (38/49).

CONCLUSIONS

Overall, genetic defects of FRMD7 were found in 11/49 (22.4%) probands, including the first reported genomic rearrangement of FRMD7 in IIN, expanding its mutational spectrum. Finally, we generate a discovery cohort of IIN patients potentially harboring either hidden a variation of FRMD7 or mutations in genes at known or novel loci sustaining the genetic heterogeneity of IIN.

Identification of six new susceptibility loci for invasive epithelial ovarian cancer

Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.

Novel pathogenic COL11A1/COL11A2 variants in Stickler syndrome detected by targeted NGS and exome sequencing

INTRODUCTION

Stickler syndrome is caused by mutations in genes encoding type II and type XI collagens. About 85% of the pathogenic variants is found in COL2A1 (Stickler type 1), whereas a minority of mutations has been reported in COL11A1 (Stickler type 2) and COL11A2 (Stickler type 3). Beside the typical skeletal and orofacial manifestations, ocular anomalies are predominantly present in type 1 and type 2, while hearing loss is more pronounced in type 2 and type 3.

METHODS

We performed COL11A1 mutation analysis for 40 type 2 Stickler patients and COL11A2 mutation analysis for five type 3 Stickler patients, previously all COL2A1 mutation-negative, using targeted next-generation sequencing (NGS) whereas whole-exome sequencing (WES) was performed in parallel for two patients. Three patients were analyzed for both genes due to unclear ocular findings.

RESULTS

In total 14 COL11A1 and two COL11A2 mutations could be identified, seven of which are novel. Splice site alterations are the most frequent mutation type, followed by glycine substitutions. In addition, six variants of unknown significance (VUS) have been found. Identical mutations and variants were identified with both NGS techniques.

CONCLUSION

We expand the mutation spectrum of COL11A1 and COL11A2 in Stickler syndrome patients and show that targeted NGS is an efficient and cost-effective molecular tool in the genetic diagnosis of Stickler syndrome, whereas the more standardized WES might be an alternative approach.

Breast-cancer risk in families with mutations in PALB2

BACKGROUND

Germline loss-of-function mutations in PALB2 are known to confer a predisposition to breast cancer. However, the lifetime risk of breast cancer that is conferred by such mutations remains unknown.

METHODS

We analyzed the risk of breast cancer among 362 members of 154 families who had deleterious truncating, splice, or deletion mutations in PALB2. The age-specific breast-cancer risk for mutation carriers was estimated with the use of a modified segregation-analysis approach that allowed for the effects of PALB2 genotype and residual familial aggregation.

RESULTS

The risk of breast cancer for female PALB2 mutation carriers, as compared with the general population, was eight to nine times as high among those younger than 40 years of age, six to eight times as high among those 40 to 60 years of age, and five times as high among those older than 60 years of age. The estimated cumulative risk of breast cancer among female mutation carriers was 14% (95% confidence interval [CI], 9 to 20) by 50 years of age and 35% (95% CI, 26 to 46) by 70 years of age. Breast-cancer risk was also significantly influenced by birth cohort (P<0.001) and by other familial factors (P=0.04). The absolute breast-cancer risk for PALB2 female mutation carriers by 70 years of age ranged from 33% (95% CI, 25 to 44) for those with no family history of breast cancer to 58% (95% CI, 50 to 66) for those with two or more first-degree relatives with breast cancer at 50 years of age.

CONCLUSIONS

Loss-of-function mutations in PALB2 are an important cause of hereditary breast cancer, with respect both to the frequency of cancer-predisposing mutations and to the risk associated with them. Our data suggest the breast-cancer risk for PALB2 mutation carriers may overlap with that for BRCA2 mutation carriers. (Funded by the European Research Council and others.).

Dealing with pseudogenes in molecular diagnostics in the next-generation sequencing era

In this chapter, we focus on issues related to the application of next-generation sequencing (NGS) strategies for the analysis of genes with pseudogenes in a clinical setting. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences, which are characterized by higher mutation rates than their functional parental genes. For the target enrichment strategies, we describe advantages and disadvantages of PCR- and capture-based enrichment methodologies. For the mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some concrete examples of genes with known pseudogenes and associated with Mendelian disorders that were analyzed by NGS on various platforms and starting from different library preparations.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing

BACKGROUND

Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting.

METHODS

We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design.

RESULTS

PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp).

CONCLUSIONS

We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

Variant ataxia telangiectasia: clinical and molecular findings and evaluation of radiosensitive phenotypes in a patient and relatives

Variant ataxia telangiectasia (A-T) may be an underdiagnosed entity. We correlate data from radiosensitivity and kinase assays with clinical and molecular data from a patient with variant A-T and relatives. The coding region of ATM was sequenced. To evaluate the functional effect of the mutations, we performed kinase assays and developed a novel S-G2 micronucleus test. Our patient presented with mild dystonia, moderately dysarthric speech, increased serum α-fetoprotein but no ataxia nor telangiectasias, no nystagmus or oculomotor dyspraxia. She has a severe IgA deficiency, but does not have recurrent infections. She is compound heterozygote for ATM c.8122G>A (p.Asp2708Asn) and c.8851-1G>T, leading to in frame loss of 63 nucleotides at the cDNA level. A trace amount of ATM protein is translated from both alleles. Residual kinase activity is derived only from the p.Asp2708Asn allele. The conventional G0 micronucleus test, based on irradiation of resting lymphocytes, revealed a radiosensitive phenotype for the patient, but not for the heterozygous relatives. As ATM is involved in homologous recombination and G2/M cell cycle checkpoint, we optimized an S-G2 micronucleus assay, allowing to evaluate micronuclei in lymphocytes irradiated in the S and G2 phases. This test showed increased radiosensitivity for both the patient and the heterozygous carriers. Intriguingly, heterozygous carriers of c.8851-1G>T (mutation associated with absence of kinase activity) showed a stronger radiosensitive phenotype with this assay than heterozygous carriers of p.Asp2708Asn (mutation associated with residual kinase activity). The modified S-G2 micronucleus assay provided phenotypic insight into complement the diagnosis of this atypical A-T patient.

Analysis of the novel fanconi anemia gene SLX4/FANCP in familial breast cancer cases

SLX4/FANCP is a recently discovered novel disease gene for Fanconi anemia (FA), a rare recessive disorder characterized by chromosomal instability and increased cancer susceptibility. Three of the 15 FA genes are breast cancer susceptibility genes in heterozygous mutation carriers--BRCA2, PALB2, and BRIP1. To investigate if defects in SLX4 also predispose to breast cancer, the gene was sequenced in a cohort of 729 BRCA1/BRCA2-negative familial breast cancer cases. We identified a single splice site mutation (c.2013+2T>A), which causes a frameshift by skipping of exon 8. We also identified 39 missense variants, four of which were selected for functional testing in a Mitomycin C-induced growth inhibition assay, and appeared indistinguishable from wild type. Although this is the first study that describes a truncating SLX4 mutation in breast cancer patients, our data indicate that germline mutations in SLX4 are very rare and are unlikely to make a significant contribution to familial breast cancer.

Practical tools to implement massive parallel pyrosequencing of PCR products in next generation molecular diagnostics

Despite improvements in terms of sequence quality and price per basepair, Sanger sequencing remains restricted to screening of individual disease genes. The development of massively parallel sequencing (MPS) technologies heralded an era in which molecular diagnostics for multigenic disorders becomes reality. Here, we outline different PCR amplification based strategies for the screening of a multitude of genes in a patient cohort. We performed a thorough evaluation in terms of set-up, coverage and sequencing variants on the data of 10 GS-FLX experiments (over 200 patients). Crucially, we determined the actual coverage that is required for reliable diagnostic results using MPS, and provide a tool to calculate the number of patients that can be screened in a single run. Finally, we provide an overview of factors contributing to false negative or false positive mutation calls and suggest ways to maximize sensitivity and specificity, both important in a routine setting. By describing practical strategies for screening of multigenic disorders in a multitude of samples and providing answers to questions about minimum required coverage, the number of patients that can be screened in a single run and the factors that may affect sensitivity and specificity we hope to facilitate the implementation of MPS technology in molecular diagnostics.

Applying massive parallel sequencing to molecular diagnosis of Marfan and Loeys-Dietz syndromes

The Marfan (MFS) and Loeys-Dietz (LDS) syndromes are caused by mutations in the fibrillin-1 (FBN1) and Transforming Growth Factor Beta Receptor 1 and 2 (TGFBR1 and TGFBR2) genes, respectively. With the current conventional mutation screening technologies, analysis of this set of genes is time consuming and expensive. We have tailored a cost-effective and reliable mutation discovery strategy using multiplex PCR followed by Next Generation Sequencing (NGS). In a first stage, genomic DNA from five MFS or LDS patient samples with previously identified mutations and/or polymorphisms in FBN1 and TGFBR1 and 2 were analyzed and revealed all expected variants. In a second stage, we validated the technique on 87 samples from MFS patients fulfilling the Ghent criteria. This resulted in the identification of 75 FBN1 mutations, of which 67 were unique. Subsequent Multiplex Ligation-dependent Probe Amplification (MLPA) analysis of the remaining negative samples identified four large deletions/insertions. Finally, Sanger sequencing identified a missense mutation in FBN1 exon 1 that was not included in the NGS workflow. In total, there was an overall mutation identification rate of 92%, which is in agreement with data published previously. We conclude that multiplex PCR of all coding exons of FBN1 and TGFBR1/2 followed by NGS analysis and MLPA is a robust strategy for time- and cost-effective identification of mutations.

Massive parallel amplicon sequencing of the breast cancer genes BRCA1 and BRCA2: opportunities, challenges, and limitations

This study describes how the new massive parallel sequencing technology can be implemented in a diagnostic setting for the breast cancer susceptibility genes (BRCA1 and BRCA2). The throughput was maximized by increasing uniformity in coverage, obtained by a multiplex approach, which outperformed pooling of singleplex PCRs. We evaluated the sensitivity by analysis of 133 distinct sequence variants; three (2%) deletions or duplications in homopolymers of greater than or equal to seven nucleotides remained undetected, illustrating a limitation of pyrosequencing. Furthermore, other limitations like nonrandom sequencing errors, pseudogene amplification, and failure to detect multiexon deletions are thoroughly described. Our workflow illustrates the potential of massive parallel sequencing of large genes in a diagnostic setting, which is of great importance to meet the increasing expectations of genetic testing. Implementation of this approach will hopefully lead to a strong reduction in turnaround times. As a consequence a wider spectrum of at risk women will be able to benefit from therapeutic interventions and prophylactic interventions.

Analysing 454 amplicon resequencing experiments using the modular and database oriented Variant Identification Pipeline

BACKGROUND

Next-generation amplicon sequencing enables high-throughput genetic diagnostics, sequencing multiple genes in several patients together in one sequencing run. Currently, no open-source out-of-the-box software solution exists that reliably reports detected genetic variations and that can be used to improve future sequencing effectiveness by analyzing the PCR reactions.

RESULTS

We developed an integrated database oriented software pipeline for analysis of 454/Roche GS-FLX amplicon resequencing experiments using Perl and a relational database. The pipeline enables variation detection, variation detection validation, and advanced data analysis, which provides information that can be used to optimize PCR efficiency using traditional means. The modular approach enables customization of the pipeline where needed and allows researchers to adopt their analysis pipeline to their experiments. Clear documentation and training data is available to test and validate the pipeline prior to using it on real sequencing data.

CONCLUSIONS

We designed an open-source database oriented pipeline that enables advanced analysis of 454/Roche GS-FLX amplicon resequencing experiments using SQL-statements. This modular database approach allows easy coupling with other pipeline modules such as variant interpretation or a LIMS system. There is also a set of standard reporting scripts available.

Genotyping of frequent BRCA1/2 SNPs with unlabeled probes: a supplement to HRMCA mutation scanning, allowing the strong reduction of sequencing burden

We previously validated mutation scanning for BRCA1 and 2 using high-resolution melting curve analysis (HRMCA). Due to recurrent single nucleotide polymorphisms (SNPs), a considerable amount of sequencing work remains after HRMCA, as melting curves for SNPs and deleterious mutations may be similar. Here, we present a simple approach for the optimization of SNP genotyping with HRMCA using unlabeled probes. Protocols were optimized for 14 frequent SNPs in BRCA1 and 2. Two probes contained an additional mismatch to detect a rare polymorphism a few nucleotides upstream. PCR was performed in the presence of LCgreenPlus and analyzed on a Lightscanner. Genotyping assays were optimized with five wild-type, heterozygous, and homozygous mutant samples. Sensitivity and specificity of the assays were evaluated with a blind screening of 95 samples. All unlabeled probes correctly genotyped the SNPs. A 1:5 asymmetric primer ratio produced sufficient probe-strand duplexes to accurately genotype the SNP of interest. The most important parameter to optimize was the number of PCR cycles. By complementing our BRCA1/2 HRMCA with 14 unlabeled probe assays, we reduced the sequencing burden by three-fold. Our simple approach for optimization can be used as a blueprint to design genotyping assays for other genes. This is one of the largest studies reported to date and the first that presents an approach combining genotyping and mutation scanning of two large polymorphic genes.

Rapid and Sensitive Detection of BRCA1/2 Mutations in a Diagnostic Setting: Comparison of Two High-Resolution Melting Platforms

Background: High-resolution melting is an emerging technique for detection of nucleic acid sequence variations. Developments in instrumentation and saturating intercalating dyes have made accurate high-resolution melting analysis possible and created opportunities to use this technology in diagnostic settings. We evaluated 2 high-resolution melting instruments for screening BRCA1 and BRCA2 mutations.

Methods: To cover the complete coding region and splice sites, we designed 112 PCR amplicons (136–435 bp), amplifiable with a single PCR program. LCGreen® Plus was used as the intercalating dye. High-resolution melting analysis was performed on the 96-well Lightscanner™ (Idaho Technology Inc.) and the 96-well LightCycler® 480 (Roche) instruments. We evaluated sensitivity by analyzing 212 positive controls scattered over almost all amplicons and specificity by blind screening of 22 patients for BRCA1 and BRCA2. In total, we scanned 3521 fragments.

Results: All 212 known heterozygous sequence variants were detected on the Lightscanner by analysis on normal sensitivity setting. On the LightCycler 480, the standard instrument sensitivity setting of 0.3 had to be increased to 0.7 to detect all variants, decreasing the specificity to 95.9% (vs 98.7% for the Lightscanner).

Conclusions: Previously, we screened BRCA1/2 by direct sequencing of the large exon 11 and denaturing gel gradient electrophoresis (DGGE) for all other coding exons. Since the introduction of high-resolution melting, our turnaround time has been one third of that with direct sequencing and DGGE, as post-PCR handling is no longer required and the software allows fast analyses. High-resolution melting is a rapid, cost-efficient, sensitive method simple enough to be readily implemented in a diagnostic laboratory.