The relevance of many BRCA2 variants of uncertain significance (VUS) to breast cancer has not been determined due to limited genetic information from families carrying these alterations. Here, we classified six new variants as pathogenic or nonpathogenic by analysis of genetic information from families carrying 64 individual BRCA2 DNA binding domain (DBD) missense mutations using a multifactorial likelihood model of cancer causality. Next, we evaluated the use of a homology-directed DNA break repair (HDR) functional assay as a method for inferring the clinical relevance of VUS in the DBD of BRCA2 using 18 established nonpathogenic missense variants and all 13 established pathogenic missense mutations from the BRCA2 DBD. Compared with the known status of these variants based on the multifactorial likelihood model, the sensitivity of the HDR assay for pathogenic mutations was estimated at 100% [95% confidence interval (CI): 75.3%–100%] and specificity was estimated at 100% (95% CI: 81.5%–100%). A statistical classifier for predicting the probability of pathogenicity of BRCA2 DBD variants was developed using these functional results. When applied to 33 additional VUS, the classifier identified eight with 99% or more probability of nonpathogenicity and 18 with 99% or more probability of pathogenicity. Thus, in the absence of genetic evidence, a cell-based HDR assay can provide a probability of pathogenicity for all VUS in the BRCA2 DBD, suggesting that the assay can be used in combination with other information to determine the cancer relevance of BRCA2 VUS. Cancer Res; 73(1); 265–75. ©2012 AACR.

The tumor suppressor gene BRCA2 encodes a large protein (3,418 amino acids) involved in the repair of DNA double strand breaks (DSB; refs. 1, 2). Inactivating germline mutations in the BRCA2 gene predispose to breast, ovarian, pancreas, and several other types of cancers (3) and confer an average cumulative lifetime risk to age 80 of 51% for breast cancer and 11%–12% for ovarian cancer (4). Because mutations that truncate or inactivate BRCA2 are associated with an elevated risk of cancer, identification of germline mutations through clinical genetic testing is widely used to identify individuals who can benefit from risk management strategies.

The majority of known pathogenic mutations in BRCA2 result in protein truncation. However, hundreds of missense BRCA2 variants (5) influence on cancer risk has also been identified in the BRCA2 gene. These variants of uncertain significance (VUS) represent a significant challenge for cancer risk assessment. Characterization of the cancer relevance of VUS has relied on a combination of genetic and in silico approaches. Specifically, a multifactorial likelihood model has been developed that estimates an overall likelihood ratio (LR) of pathogenicity for each VUS based on cosegregation of variants with cancer in families, co-occurrence of variants with known pathogenic mutations, empirical evaluation of personal and family history of cancer, and tumor histopathology (6–9). The likelihood of pathogenicity can be combined with the prior probability of pathogenicity, based on cross-species sequence conservation and physicochemical properties of mutated residues (Align GVGD; ref. 10), to generate a posterior probability of pathogenicity (11, 12). The posterior probability of pathogenicity is used for categorizing VUS according to a 5-tier classification system introduced by an International Agency for Research on Cancer (IARC) Working Group (12). Briefly, Class 1 (posterior probability < 0.001) and Class 2 (0.001 < posterior probability < 0.049) VUS are nonpathogenic and likely nonpathogenic VUS. Class 3 (0.05 < posterior probability < 0.949) VUS remain unclassified because of lack of sufficient family information for classification. Class 4 (0.95 < posterior probability < 0.99) and Class 5 (posterior probability > 0.99) VUS are likely pathogenic and pathogenic, respectively. Many VUS in the BRCA1 and BRCA2 genes lack sufficient family information for classification. As a result, alternative methods are needed to interpret VUS pathogenicity. One approach is to use functional assays that assess the impact of genetic variants on the activity of the protein (13, 14). In this study, we evaluate the ability of a cell-based homology-directed repair (HDR) assay to assess the pathogenicity of missense variants in the DNA binding domain (DBD) of BRCA2. While this assay has been used to identify missense mutations that inactivate BRCA2 (15, 16), the sensitivity and specificity of the assay for classified (12) pathogenic BRCA2 mutations has not been established. Here, we identified a series of pathogenic and nonpathogenic BRCA2 DBD mutations and used these standards to define the sensitivity and specificity of the HDR assay. On the basis of the results, we also developed a statistical classifier for prediction of the pathogenicity of all BRCA2 DBD variants and applied this model to 33 additional VUS from the BRCA2 DBD.

VUS selection

A total of 64 missense variants localized to the DBD of BRCA2 were selected for functional analysis by the HDR assay. All VUS were identified in patients who underwent clinical genetic testing for mutations in the BRCA1 and BRCA2 genes. Of these, 63 are reported in the Breast Cancer Information Core (BIC) database (5). Details of selection criteria are provided in the Results section.

Posterior probability approach

LR of pathogenicity for each variant were estimated on the basis of cosegregation of the variant with breast cancer in families of mutation carriers; co-occurrence (in trans) with known pathogenic mutations, and the family and personal history of cancer of individuals found to carry the variant (6, 9). Family information was obtained by Myriad Genetics Laboratories Inc. (Salt Lake City, Utah), cancer risk assessment clinics at the Mayo Clinic, the University of Pennsylvania (Philadelphia, PA), The Research Institute of the McGill University Health Center (Quebec, Canada), the German Hereditary Breast Cancer Consortium (GC-HBOC), the Austrian Hereditary Breast and Ovarian Cancer Group, and the Instituto Nazionale dei Tumori in Milan (Milan, Italy). Prior probabilities of pathogenicity based on sequence analysis were derived from Align GVGD (17). A prior probability of 0.01 [95% confidence interval (CI): 0.00–0.06] was assigned to substitutions at evolutionary variable positions (18), located within the BRCA2 DBD functional domain (11).

HDR assay

The selected missense variants were cloned into a 3× FLAG-tagged full-length BRCA2 cDNA expression plasmid using the QuickChange Site-Directed Mutagenesis Kit XL (Stratagene; ref. 13). Mutations were verified by DNA sequencing. Wild-type and mutant constructs were cotransfected with an I-Sce1–expressing pcBASce plasmid (13) into V-C8 cells containing the DR-GFP reporter plasmid. The BRCA2-deficient V-C8 hamster lung fibroblast cell line (XRCC11) displays chromosomal instability and abnormal centrosomes, reduced nuclear localization of the RAD51 protein that is central in the HDR processes, and sensitivity to cross-linking agents such as methyl methanesulfonate (19). Cells expressing GFP following HDR-dependent repair of an I-Sce1–induced DNA double strand break in the DR-GFP reporter plasmid were quantified by fluorescence-activated cell sorting (FACS) after 72 hours. Transfection efficiency was evaluated by immunofluorescence-based counting of GFP-positive cells. Expression of wild-type and mutant BRCA2 was evaluated by immunoblotting with rabbit Anti-Flag antibody (1:1,000 F7425; Sigma) of mouse Anti-Flag M2 (F1804; Sigma) immunoprecipitates from protein lysates.

Statistical analysis

Data normalization.

Mutants were evaluated by the HDR assay in batches. In each batch, wild-type BRCA2 and the D2723H known pathogenic mutant were included as positive and negative controls. The basal HDR activity due to the transfection of the vector alone was also measured. All experiments were carried out in duplicate. The mean fold change in GFP-positive cells for each mutant relative to vector alone was calculated. These values were normalized relative to the ratio of D2723H to wild-type BRCA2 in each experiment and rescaled to reference fold-change values of 5.0 for wild-type and 1.0 for D2723H, which reflected the prenormalization HDR fold changes while providing simple integer values as anchors for known nonpathogenic or pathogenic variants, respectively. Using these normalized values, the average normalized log (HDR) value and SE for each variant were computed while pooling data across all replicated experiments. Because the HDR assay results on the fold-difference scale are positively skewed, the SE increases as a function of the mean HDR fold difference. This relationship is not evident when computing the summary statistics from log-transformed data. Therefore, all statistical analyses were conducted after transforming the fold-difference measurements to the log scale. A linear mixed model was also used to analyze all data available for each variant and to extract model-based estimates of the mean and SE of the HDR activity for each VUS.

Sensitivity and specificity.

A receiver operating characteristics (ROC) analysis (20) of the normalized mean HDR assay results of the known pathogenic and nonpathogenic variants was conducted to identify an HDR assay cutoff that simultaneously maximized both the sensitivity and specificity of the assay for pathogenic variants. The midpoint between the extremes of the ranges of HDR assay values for pathogenic and nonpathogenic variants was selected. The sensitivity and specificity of the assay was calculated together with exact 95% CIs based on this cut-point. To determine the degree to which variability in the assay might influence the results, the sensitivity and specificity of the HDR assay was also estimated using the most extreme assay reading for each variant: the highest observed values for the pathogenic variants and the lowest observed measurements for the nonpathogenic variants.

VUS classification via the HDR assay.

A statistical classifier for determining the clinical relevance of BRCA2 variants was developed using the HDR data from variants already classified as pathogenic or nonpathogenic by means of the multifactorial likelihood method. A discriminant analysis approach was used to develop the classifier (21, 22). After verifying that the distribution of the normalized and log transformed HDR assay results were well approximated by Gaussian distributions, the data observed for the known pathogenic and nonpathogenic variants were used to estimate the parameters of the 2 Gaussian distributions that describe the HDR activity within each of these 2 classes of variant. Linear mixed models that incorporated random per-variant and per-experiment intercepts to allow for the presence of repeated measurements were used to estimate means and variance components that quantified the behavior of the HDR assay when applied to known pathogenic and nonpathogenic variants. Separate models were fit to the normalized and loge-transformed data from the known pathogenic and nonpathogenic variants and the model parameters were extracted. The estimated Gaussian distribution of the nonpathogenic variants, fN(x), had a mean of 1.57 and a variance equal to the sum of the variance components from the linear mixed model or 0.06. The equivalent distribution of the pathogenic variants, fP(x), had a mean of 0.05 and a variance of 0.07.

Given the distinct Gaussian distributions for the nonpathogenic and pathogenic variants, fN(x) and fP(x), respectively, the posterior probability that a new variant, with HDR assay result equal to exp(x), is pathogenic is:

where pP is the prior probability that the VUS is pathogenic. It is possible to use external estimates of prior probability in this equation, such as prior probabilities based on the Align GVGD model (17). However, for this study, we chose to set equal prior probabilities that a variant was either pathogenic or nonpathogenic because the results obtained with the HDR-based classifications are meant to be independent of other data. Consequently, it is possible to compare the HDR-based posterior probabilities with the posterior probabilities obtained using other methods, such as the multifactorial likelihood model or its components. Assuming equal prior probabilities, we computed posterior probabilities that each variant was pathogenic based on the HDR assay result. We elected to consider the assay result as inconclusive unless the posterior probability of pathogenicity was either ≥0.99 or ≤0.01. To assess the influence of the assumed prior probabilities on this model, we also estimated posterior probabilities of pathogenicity using prior probabilities from the Align GVGD model.

Identification of pathogenic and nonpathogenic variants using a multifactorial likelihood model

In this study, we focused on variants in the C-terminal DBD of BRCA2 because the region has been implicated in the homologous recombination (HDR) activity of BRCA2 and contains all of the known pathogenic BRCA2 missense mutations (18, 23). In addition, many of the highly conserved amino acids in the region, that are likely important for function, have been found to contain missense mutations in individuals with a family history of breast and/or ovarian cancer. Likelihood ratios of pathogenicity for many variants in the BRCA2 DBD have been estimated on the basis of a large collection of family data provided by Myriad Genetics Laboratories (6). Posterior probabilities of pathogenicity have also been estimated by combining the likelihood ratios with prior probabilities of pathogenicity based on Align GVGD modeling of protein sequences (18, 23).

Collection of family information for BRCA2 DBD Class 3 VUS with posterior probabilities of pathogenicity approaching Class 2 (<0.05) or Class 4 (>0.95) thresholds yielded 18 pedigrees that were informative for segregation analysis. Reestimation of the posterior probability of pathogenicity for the associated variants resulted in classification of 3 variants (S3020C, S2414L, and T3013I) as Class 2 (likely nonpathogenic) and 3 (G2609D, L2688P, and N3124I) as Class 4 (likely pathogenic). When these Class 4 missense substitutions were combined with all known Class 4 and Class 5 missense mutations (18, 23), a total of 13 known pathogenic missense mutations from the BRCA2 DBD were available for further study (Table 1). Of these, 8 were Class 5 and 5 were Class 4. In addition, we selected the 3 new Class 2 variants, 13 previously classified Class 1 and Class 2 variants (23), and the polymorphisms A2466V (24) and A2951T (25, 26) as nonpathogenic standards for our HDR studies (Table 1).

Table 1.

Classification of BRCA2 DBD variants as pathogenic and nonpathogenic by a multifactorial probability-based model

LRs derived from each component
VUSAlign GVGDaCosegregation oddsCo-occurence oddsFamily history oddsCombined odds in favor of causalityPosteriorPbIARC classc
K2411T 0.01 1.00 0.06 0.15 0.01 8.49 × 10−5 
S2414L 0.01 1.00 1.02 2.79 2.86 0.03 
A2466V 0.01 — — — — — — 
P2589H 0.01 1.00 1.02 1.44 1.48 0.01 
G2609D 0.81 1.00 1.02 4.02 4.12 0.95 
W2626C 0.81 1.00 2.13 22.65 48.14 1.00 
I2627F 0.29 2.01 1.28 407.62 1.05 × 103 1.00 
L2647P 0.81 1.00 1.19 3.86 4.58 0.95 
L2653P 0.81 0.58 1.16 35.96 24.06 0.99 
D2665G 0.81 0.01 0.13 0.24 3.34 × 10−4 1.42 × 10−3 
M2676T 0.01 1.00 1.16 0.83 0.96 0.01 
L2688P 0.81 1.00 1.19 3.44 4.08 0.95 
T2722R 0.81 10.80 1.24 6.95 93.47 1.00 
D2723G 0.81 1.00 1.59 75.96 120.57 1.00 
D2723H 0.81 1.38 × 104 4.17 9.77 × 106 5.62 × 1011 1.00 
K2729N 0.01 0.04 6.35 0.01 1.63 × 10−3 1.65 × 10−5 
G2748D 0.81 6.73 1.44 257.49 2.49 × 103 1.00 
R2787H 0.01 1.00 1.05 1.25 1.31 0.01 
E2856A 0.01 NA NA 3.89 × 10−11 3.89 × 10−11 3.93 × 10−13 
V2908G 0.66 1.00 0.06 0.07 3.91 × 10−3 0.01 
A2951T 0.01 — — — — — — 
D2965H 0.01 1.00 1.08 1.07 1.15 0.01 
R2973C 0.01 2.36 × 10−3 1.84 0.04 1.76 × 10−4 1.78 × 10−6 
T3013I 0.01 1.82 NA NA 1.82 0.02 
S3020C 0.01 1.00 1.02 1.94 1.98 0.02 
R3052W 0.81 3.98 × 103 1.48 1.61 9.55 × 103 0.99 
P3063S 0.01 1.00 1.10 0.67 0.73 0.01 
Y3092C 0.81 1.02 1.28 4.62 × 10−3 0.01 0.03 
D3095E 0.66 0.43 1.67 31.82 22.59 0.98 
Y3098H 0.01 5.03 0.19 3.54 × 10−4 3.46 × 10−4 3.49 × 10−6 
N3124I 0.81 0.04 2.07 81.64 7.23 0.97 
LRs derived from each component
VUSAlign GVGDaCosegregation oddsCo-occurence oddsFamily history oddsCombined odds in favor of causalityPosteriorPbIARC classc
K2411T 0.01 1.00 0.06 0.15 0.01 8.49 × 10−5 
S2414L 0.01 1.00 1.02 2.79 2.86 0.03 
A2466V 0.01 — — — — — — 
P2589H 0.01 1.00 1.02 1.44 1.48 0.01 
G2609D 0.81 1.00 1.02 4.02 4.12 0.95 
W2626C 0.81 1.00 2.13 22.65 48.14 1.00 
I2627F 0.29 2.01 1.28 407.62 1.05 × 103 1.00 
L2647P 0.81 1.00 1.19 3.86 4.58 0.95 
L2653P 0.81 0.58 1.16 35.96 24.06 0.99 
D2665G 0.81 0.01 0.13 0.24 3.34 × 10−4 1.42 × 10−3 
M2676T 0.01 1.00 1.16 0.83 0.96 0.01 
L2688P 0.81 1.00 1.19 3.44 4.08 0.95 
T2722R 0.81 10.80 1.24 6.95 93.47 1.00 
D2723G 0.81 1.00 1.59 75.96 120.57 1.00 
D2723H 0.81 1.38 × 104 4.17 9.77 × 106 5.62 × 1011 1.00 
K2729N 0.01 0.04 6.35 0.01 1.63 × 10−3 1.65 × 10−5 
G2748D 0.81 6.73 1.44 257.49 2.49 × 103 1.00 
R2787H 0.01 1.00 1.05 1.25 1.31 0.01 
E2856A 0.01 NA NA 3.89 × 10−11 3.89 × 10−11 3.93 × 10−13 
V2908G 0.66 1.00 0.06 0.07 3.91 × 10−3 0.01 
A2951T 0.01 — — — — — — 
D2965H 0.01 1.00 1.08 1.07 1.15 0.01 
R2973C 0.01 2.36 × 10−3 1.84 0.04 1.76 × 10−4 1.78 × 10−6 
T3013I 0.01 1.82 NA NA 1.82 0.02 
S3020C 0.01 1.00 1.02 1.94 1.98 0.02 
R3052W 0.81 3.98 × 103 1.48 1.61 9.55 × 103 0.99 
P3063S 0.01 1.00 1.10 0.67 0.73 0.01 
Y3092C 0.81 1.02 1.28 4.62 × 10−3 0.01 0.03 
D3095E 0.66 0.43 1.67 31.82 22.59 0.98 
Y3098H 0.01 5.03 0.19 3.54 × 10−4 3.46 × 10−4 3.49 × 10−6 
N3124I 0.81 0.04 2.07 81.64 7.23 0.97 

NOTE: GenBank Reference BRCA2 NM_000059.3. —, not applicable because the variant is a well-defined polymorphism.

Abbreviation: NA, not available.

aPrior probability based on Align GVGD (17).

bPosterior probability = Posterior odds/(posterior odds + 1), where the posterior odds = odds for causality × [prior probability/(1-prior probability)].

cClass 1/2: nonpathogenic/likely nonpathogenic; Class 4/5: pathogenic/likely pathogenic (12).

A further 27 VUS from the BRCA2 DBD were selected for evaluation based on high cross-species sequence conservation of the parental residues (Align GVGD-based prior probability of 0.81). Six other VUS with prior probabilities ranging from 0.01 to 0.61 were included as an initial investigation of the correlation between probabilities derived from Align GVGD and the HDR assay. Family data were available for 22 of these 33 variants, but the information was insufficient for classification of the VUS as pathogenic or nonpathogenic using the multifactorial likelihood model. Thus, the 22 variants with family data and the 11 without data remained classified as Class 3 VUS (Table 2).

Table 2.

Genetic and functional characterization of 64 BRCA2 DBD variants

Multifactorial likelihood modelFunctional model
VUSAlign GVGDaPosteriorPbIARC classcVUSHDR normalizedHDR SEPr(deleterious)dPr(neutral)d
T3013I 0.01 0.02 T3013I 5.96 0.48 1.37 ×10−7 
K2411T 0.01 8.49 × 10−5 K2411T 5.74 0.43 2.25 × 10−7 
M2676T 0.01 0.01 M2676T 5.68 0.43 2.67 × 10−7 
S2414L 0.01 0.03 S2414L 5.58 0.52 5.55 × 10−7 
E2856A 0.01 3.93 × 10−13 E2856A 5.54 0.42 3.99 × 10−7 
Y3098H 0.01 3.49 × 10−6 Y3098H 5.53 0.16 1.84 × 10−7 
D2665G 0.81 1.42 × 10−3 D2665G 5.30 0.65 2.99 × 10−6 
D2965H 0.01 0.01 D2965H 5.23 0.39 1.01 × 10−6 
A2466V 0.01 — — A2466V 4.68 0.43 8.64 × 10−6 
P3063S 0.01 0.01 P3063S 4.43 0.21 9.96 × 10−6 
P2589H 0.01 0.01 P2589H 4.35 0.39 2.67 × 10−6 
Y3092C 0.81 0.03 Y3092C 4.29 0.36 2.94 × 10−5 
A2951T 0.01 — — A2951T 4.14 0.38 6.02 × 10−5 
V2908G 0.66 0.01 V2908G 3.94 0.50 2.60 × 10−4 
S3020C 0.01 0.02 S3020C 3.69 0.24 2.79 × 10−4 
R2787H 0.01 0.01 R2787H 3.61 0.27 4.55 × 10−4 
K2729N 0.01 1.65 × 10−5 K2729N 3.30 0.42 3.7 × 10−3 
R2973C 0.01 1.78 × 10−6 R2973C 3.15 0.40 7.57 × 10−3 0.99 
F2406L 0.01 NA F2406L 6.58 0.60 4.05 × 10−8 
L2865V 0.29 0.25 L2865V 4.96 0.64 9.27 × 10−6 
G2813E 0.81 0.77 G2813E 4.66 0.42 9.02 × 10−6 
H2440R 0.01 NA H2440R 4.30 0.34 2.72 × 10−5 
D2679G 0.81 0.73 D2679G 4.25 0.38 3.84 × 10−5 
D2679Y 0.81 0.83 D2679Y 4.09 0.25 4.64 × 10−5 
S2807L 0.81 NA S2807L 3.84 0.35 2.01 × 10−4 
G2901D 0.81 0.73 G2901D 3.48 0.23 7.43 × 10−4 
T3033I 0.81 0.9 T3033I 2.80 0.21 0.03 0.97 
R2842L 0.81 NA R2842L 2.78 0.21 0.04 0.96 
E3002D 0.66 0.75 E3002D 2.69 0.25 0.06 0.94 
Y3092S 0.81 0.64 Y3092S 2.51 0.12 0.17 0.83 
G2812E 0.81 0.75 G2812E 2.51 0.18 0.17 0.83 
P2800S 0.81 NA P2800S 2.43 0.18 0.27 0.73 
P2800R 0.81 NA P2800R 2.43 0.15 0.27 0.73 
R2842C 0.81 0.84 R2842C 1.76 0.16 0.99 0.01 
Y2726C 0.81 0.76 Y2726C 1.70 0.16 0.99 6.33 × 10−3 
D2611G 0.81 NA D2611G 1.67 0.15 4.39 × 10−3 
D3073G 0.81 0.86 D3073G 1.52 0.14 8.93 × 10−4 
L2654P 0.81 0.79 L2654P 1.45 0.13 4.06 × 10−4 
G2793E 0.81 0.75 G2793E 1.32 0.12 7.78 × 10−5 
Y2660D 0.81 NA Y2660D 1.30 0.10 4.54 × 10−5 
G2585R 0.81 NA G2585R 1.26 0.12 3.55 × 10−5 
G2793R 0.81 0.61 G2793R 1.22 0.11 1.94 × 10−5 
L2792P 0.81 0.92 L2792P 1.20 0.11 1.35 × 10−5 
R2659G 0.81 NA R2659G 1.12 0.07 2.51 × 10−6 
R2784Q 0.66 0.12 R2784Q 1.11 0.10 3.97 × 10−6 
G3076V 0.81 0.94 G3076V 1.09 0.10 2.80 × 10−6 
E3002K 0.66 0.66 E3002K 1.05 0.10 1.61 × 10−6 
G3076E 0.81 0.92 G3076E 1.04 0.09 1.14 × 10−6 
L2510P 0.81 NA L2510P 1.04 0.10 1.16 × 10−6 
D2723A 0.81 0.85 D2723A 0.88 0.08 6.40 × 10−6 
R2784W 0.81 0.83 R2784W 0.83 0.08 2.42 × 10−8 
W2626C 0.81 W2626C 1.30 0.12 6.13 × 10−5 
G2609D 0.81 0.95 G2609D 1.28 0.08 2.91 × 10−5 
T2722R 0.81 T2722R 1.26 0.12 3.55 × 10−5 
I2627F 0.29 I2627F 1.09 0.10 2.79 × 10−6 
D2723G 0.81 D2723G 1.03 0.10 1.00 × 10−6 
D3095E 0.66 0.98 D3095E 1.02 0.09 8.42 × 10−7 
L2688P 0.81 0.95 L2688P 1.01 0.09 6.90 × 10−7 
D2723H 0.81 D2723H 1.00 0.02 1.39 × 10−7 
L2647P 0.81 0.95 L2647P 0.99 0.09 5.19 × 10−7 
N3124I 0.81 0.97 N3124I 0.96 0.10 3.91 × 10−7 
L2653P 0.81 0.99 L2653P 0.96 0.09 2.77 × 10−7 
G2748D 0.81 G2748D 0.83 0.08 2.32 × 10−8 
R3052W 0.81 0.99 R3052W 0.81 0.07 1.38 × 10−8 
Multifactorial likelihood modelFunctional model
VUSAlign GVGDaPosteriorPbIARC classcVUSHDR normalizedHDR SEPr(deleterious)dPr(neutral)d
T3013I 0.01 0.02 T3013I 5.96 0.48 1.37 ×10−7 
K2411T 0.01 8.49 × 10−5 K2411T 5.74 0.43 2.25 × 10−7 
M2676T 0.01 0.01 M2676T 5.68 0.43 2.67 × 10−7 
S2414L 0.01 0.03 S2414L 5.58 0.52 5.55 × 10−7 
E2856A 0.01 3.93 × 10−13 E2856A 5.54 0.42 3.99 × 10−7 
Y3098H 0.01 3.49 × 10−6 Y3098H 5.53 0.16 1.84 × 10−7 
D2665G 0.81 1.42 × 10−3 D2665G 5.30 0.65 2.99 × 10−6 
D2965H 0.01 0.01 D2965H 5.23 0.39 1.01 × 10−6 
A2466V 0.01 — — A2466V 4.68 0.43 8.64 × 10−6 
P3063S 0.01 0.01 P3063S 4.43 0.21 9.96 × 10−6 
P2589H 0.01 0.01 P2589H 4.35 0.39 2.67 × 10−6 
Y3092C 0.81 0.03 Y3092C 4.29 0.36 2.94 × 10−5 
A2951T 0.01 — — A2951T 4.14 0.38 6.02 × 10−5 
V2908G 0.66 0.01 V2908G 3.94 0.50 2.60 × 10−4 
S3020C 0.01 0.02 S3020C 3.69 0.24 2.79 × 10−4 
R2787H 0.01 0.01 R2787H 3.61 0.27 4.55 × 10−4 
K2729N 0.01 1.65 × 10−5 K2729N 3.30 0.42 3.7 × 10−3 
R2973C 0.01 1.78 × 10−6 R2973C 3.15 0.40 7.57 × 10−3 0.99 
F2406L 0.01 NA F2406L 6.58 0.60 4.05 × 10−8 
L2865V 0.29 0.25 L2865V 4.96 0.64 9.27 × 10−6 
G2813E 0.81 0.77 G2813E 4.66 0.42 9.02 × 10−6 
H2440R 0.01 NA H2440R 4.30 0.34 2.72 × 10−5 
D2679G 0.81 0.73 D2679G 4.25 0.38 3.84 × 10−5 
D2679Y 0.81 0.83 D2679Y 4.09 0.25 4.64 × 10−5 
S2807L 0.81 NA S2807L 3.84 0.35 2.01 × 10−4 
G2901D 0.81 0.73 G2901D 3.48 0.23 7.43 × 10−4 
T3033I 0.81 0.9 T3033I 2.80 0.21 0.03 0.97 
R2842L 0.81 NA R2842L 2.78 0.21 0.04 0.96 
E3002D 0.66 0.75 E3002D 2.69 0.25 0.06 0.94 
Y3092S 0.81 0.64 Y3092S 2.51 0.12 0.17 0.83 
G2812E 0.81 0.75 G2812E 2.51 0.18 0.17 0.83 
P2800S 0.81 NA P2800S 2.43 0.18 0.27 0.73 
P2800R 0.81 NA P2800R 2.43 0.15 0.27 0.73 
R2842C 0.81 0.84 R2842C 1.76 0.16 0.99 0.01 
Y2726C 0.81 0.76 Y2726C 1.70 0.16 0.99 6.33 × 10−3 
D2611G 0.81 NA D2611G 1.67 0.15 4.39 × 10−3 
D3073G 0.81 0.86 D3073G 1.52 0.14 8.93 × 10−4 
L2654P 0.81 0.79 L2654P 1.45 0.13 4.06 × 10−4 
G2793E 0.81 0.75 G2793E 1.32 0.12 7.78 × 10−5 
Y2660D 0.81 NA Y2660D 1.30 0.10 4.54 × 10−5 
G2585R 0.81 NA G2585R 1.26 0.12 3.55 × 10−5 
G2793R 0.81 0.61 G2793R 1.22 0.11 1.94 × 10−5 
L2792P 0.81 0.92 L2792P 1.20 0.11 1.35 × 10−5 
R2659G 0.81 NA R2659G 1.12 0.07 2.51 × 10−6 
R2784Q 0.66 0.12 R2784Q 1.11 0.10 3.97 × 10−6 
G3076V 0.81 0.94 G3076V 1.09 0.10 2.80 × 10−6 
E3002K 0.66 0.66 E3002K 1.05 0.10 1.61 × 10−6 
G3076E 0.81 0.92 G3076E 1.04 0.09 1.14 × 10−6 
L2510P 0.81 NA L2510P 1.04 0.10 1.16 × 10−6 
D2723A 0.81 0.85 D2723A 0.88 0.08 6.40 × 10−6 
R2784W 0.81 0.83 R2784W 0.83 0.08 2.42 × 10−8 
W2626C 0.81 W2626C 1.30 0.12 6.13 × 10−5 
G2609D 0.81 0.95 G2609D 1.28 0.08 2.91 × 10−5 
T2722R 0.81 T2722R 1.26 0.12 3.55 × 10−5 
I2627F 0.29 I2627F 1.09 0.10 2.79 × 10−6 
D2723G 0.81 D2723G 1.03 0.10 1.00 × 10−6 
D3095E 0.66 0.98 D3095E 1.02 0.09 8.42 × 10−7 
L2688P 0.81 0.95 L2688P 1.01 0.09 6.90 × 10−7 
D2723H 0.81 D2723H 1.00 0.02 1.39 × 10−7 
L2647P 0.81 0.95 L2647P 0.99 0.09 5.19 × 10−7 
N3124I 0.81 0.97 N3124I 0.96 0.10 3.91 × 10−7 
L2653P 0.81 0.99 L2653P 0.96 0.09 2.77 × 10−7 
G2748D 0.81 G2748D 0.83 0.08 2.32 × 10−8 
R3052W 0.81 0.99 R3052W 0.81 0.07 1.38 × 10−8 

NOTE: GenBank Reference BRCA2 NM_000059.3.

Abbreviation: NA, not available.

—, not applicable because the variant is a well-defined polymorphism.

aPrior probability based on Align GVGD (17).

bPosterior probability = Posterior odds/(posterior odds + 1), where the posterior odds = odds for causality × [prior probability/(1-prior probability)].

cClass 1/2: nonpathogenic/likely nonpathogenic; Class 3: unclassified; and Class 4/5: pathogenic/likely pathogenic (12).

dModel-based estimates of conditional probabilities in favor or against pathogenicity for VUS based on HDR activity.

Sensitivity and specificity of the HDR assay

To assess the sensitivity and specificity of the HDR assay for pathogenic mutations, the influence of the 13 known pathogenic and 18 known nonpathogenic mutations on the HDR activity of BRCA2 was evaluated. Wild-type and mutant forms of BRCA2 were coexpressed with I-Sce1 in BRCA2-deficient V-C8 cells containing the DR-GFP reporter. The BRCA2 constructs exhibited similar transfection efficiencies and produced similar levels of full-length BRCA2 protein (Supplementary Fig. S1; ref. 13). Cells expressing GFP after BRCA2-dependent repair of I-Sce1 DNA DSBs were quantified by flow cytometry. Normalized fold changes in GFP-positive cells, which equate to HDR activity, are shown in Fig. 1 for the known pathogenic and nonpathogenic variants. Importantly, the HDR activity of the Class 1 and 2 variants did not overlap with the HDR activity of Class 4 or 5 mutations (Fig. 1). For the purpose of assessing sensitivity and specificity of the HDR assay, the midpoint between the mean of the normalized fold change for Class 1/2 variants and the mean for the Class 4/5 variants was selected as a cut-point for discriminating between pathogenic and nonpathogenic variants. Using this threshold of 2.25, the 13 pathogenic mutations and 18 nonpathogenic variants were correctly classified. Sensitivity was estimated at 100% (95% CI: 75.3%–100%), and specificity was estimated at 100% (95% CI: 81.5%–100%). Importantly, these estimates were not altered when considering the lowest observed HDR assay values for the nonpathogenic variants and the highest observed values for the pathogenic variants instead of mean values. Furthermore, when restricting to the 7 Class 1 and 8 Class 5 variants, the sensitivity (100%; 95% CI: 63.1%–100%) and the specificity (100%; 95% CI: 71.5%–100%) of the assay were again estimated at 100%. Likewise the 5 Class 4 variants and 11 Class 2 variants yielded a sensitivity of 100% (95% CI: 47.8%–100%) and specificity of 100% (95% CI: 71.5%–100%). Given these results, we combined the Class 1 and Class 2 variants as nonpathogenic and the Class 4 and Class 5 variants as pathogenic for the purpose of developing a HDR-based probability model of pathogenicity.

Figure 1.

HDR activity of 31 classified nonpathogenic variants and pathogenic mutations. The normalized HDR fold increase with SE is displayed on a logarithmic scale for the combined Class 1/Class 2 variants and combined Class 4/Class 5 mutations.

Figure 1.

HDR activity of 31 classified nonpathogenic variants and pathogenic mutations. The normalized HDR fold increase with SE is displayed on a logarithmic scale for the combined Class 1/Class 2 variants and combined Class 4/Class 5 mutations.

Close modal

Evaluation of VUS with the HDR assay

Given the high sensitivity and specificity of the HDR assay for pathogenic BRCA2 DBD mutations, we evaluated the influence of 33 Class 3 VUS from the BRCA2 DBD on BRCA2 HDR activity in an effort to determine the functional and disease relevance of these VUS. The transfection efficiencies and expression levels of the mutants were similar to wild-type BRCA2 (Supplementary Fig. S1; ref. 13). The normalized HDR results for each VUS are shown in Supplementary Table S1. Thirteen VUS displayed low HDR activity within the range of the known pathogenic mutations (Table 2, Fig. 2). In addition, 8 VUS displayed high HDR activity within the range of activity of the known nonpathogenic mutations (Table 2, Fig. 2). A further 5 VUS displayed activity that was marginally above the activity of known pathogenic mutations. However, 7 VUS seemed to have intermediate HDR activity outside the ranges of the known pathogenic or nonpathogenic mutations (Fig. 2). Recloning and revaluation of these variants did not alter these results.

Figure 2.

HDR activity of 64 BRCA2 DBD mutants. The discriminant model-based normalized HDR fold change with SE is displayed on a logarithmic scale for the analyzed nonpathogenic, unclassified, and pathogenic variants. Solid lines represent 99.9% and 0.1% probability of pathogenicity and dotted lines represent 99% and 1% probability of pathogenicity. The SE is included only as a measure of the reproducibility of the HDR assay for each variant.

Figure 2.

HDR activity of 64 BRCA2 DBD mutants. The discriminant model-based normalized HDR fold change with SE is displayed on a logarithmic scale for the analyzed nonpathogenic, unclassified, and pathogenic variants. Solid lines represent 99.9% and 0.1% probability of pathogenicity and dotted lines represent 99% and 1% probability of pathogenicity. The SE is included only as a measure of the reproducibility of the HDR assay for each variant.

Close modal

We also evaluated possible influences of the VUS on mRNA splicing. Here, we calculated consensus values of potential splice sites using the MaxEnt (27) algorithm and Human Splicing Finder (HSF; refs. 9, 28). The c.8035G>T mutation (p.D2679Y) had an HSF score of 87.91 (wild-type 61.07) and a MaxEnt score of 9.11 (wild-type 1.46) and a moderate probability in favor of a de novo donor site. Aberrant in vivo splicing was not assessed because of lack of specimens. In addition, the D2723G variant was previously shown to cause aberrant splicing (deletion of 163 bp from exon 18) of a small proportion of the mutant allele (9). In contrast, we found low probabilities of de novo splice acceptor or donor formation for all other VUS in this study. This suggested that the variants evaluated in this study likely do not influence breast cancer risk through aberrant splicing.

A statistical classifier for BRCA2 DBD VUS based on the HDR assay

Although the HDR assay results for the 33 VUS suggested that a number of the VUS were functionally relevant, a quantitative model that accounted for the variance in the assay was needed to make predictions of pathogenicity or nonpathogenicity. Therefore, we developed a statistical classifier to estimate the probability of pathogenicity for any BRCA2 DBD evaluated by the HDR assay. Specifically, we estimated the mean and variances of the distributions of the HDR assay results for the known pathogenic and nonpathogenic mutation groups and used these results to compute the probability of each tested VUS being pathogenic, conditional on the observed HDR assay result (Table 2). Variants with 99% or more probability of pathogenicity were classified as pathogenic and variants with 1% or less probability of pathogenicity were classified as nonpathogenic. Thus, 18 unclassified variants were considered pathogenic, 8 were considered nonpathogenic, and 7 VUS remained unclassified.

Because each variant displays a unique degree of variability in the HDR assay, the final posterior probability that any specific variant is pathogenic or nonpathogenic will depend on the results of individual assays. However, we also estimated cut-points in HDR activity that could be used for visual inspection of the probability and classification of any DBD variant. Specifically, the minimum values of the HDR assay that grouped all variants with probabilities of at least 99% (≥3.02) and 99.9% (≥3.48) of nonpathogenicity were identified (Fig. 2). Similarly, the maximum values of the HDR assay that grouped all variants with probabilities of 99% (≤1.75) and 99.9% (≤1.53) of pathogenicity were identified (Fig. 2). Importantly, the SE associated with each variant is only shown to provide information on the reproducibility of the HDR results for each variant because the variability of the mean HDR result for each variant is incorporated into the computation of the selected cut-points. Therefore, a variant is considered classified even when the SE of the normalized fold change in HDR activity crosses these visual cut-points.

Structural interpretation of predicted pathogenic variants

Thirteen BRCA2 DBD variants were classified as pathogenic using the multifactorial likelihood model and a further 18 variants were classified as pathogenic by the HDR assay alone. In an attempt to understand the mechanism by which these variants altered BRCA2 function the location of the mutated residues in the predicted 3-dimensional structure of the mouse BRCA2-DSS1-ssDNA domain (29) was determined and the impact of the mutations on the local protein structure was assessed (Fig. 3). The missense variants mapped to 4 regions of the BRCA2 DBD, region 1: Helical domain (Fig. 3A); region 2: Helical-OB1-DSS1 interface (Fig. 3B); region 3: OB1-OB2-DSS1 interface (Fig. 3C), and region 4: OB2-OB3 interface (Fig. 3D). In region 1, several mutated residues were associated with the terminal helix of the helical domain. Residues Y2581 (human Y2660), L2574 (human L2653), and L2575 (human L2654) are located on one face of the helix, confined to the helical domain, whereas R2580 (human R2659) is located on the opposite face of the helix that interacts with OB1 and DSS1. Thus Y2581, L2575, and L2576 are likely more important for the stability of the helix and the helical domain, whereas R2580 may influence the interaction of the helical domain with the OB1 domain and DSS1. This is supported by the observation that R2580 is potentially involved in hydrogen bond contacts with H2544 (human H2623) within the helical domain and G2645 (human G2724) of the OB1 domain. In region 2, a cluster of residues containing likely pathogenic variants is connected through a network of hydrogen bonds. Disruption of this network by mutations is potentially detrimental to the integrity of the helical-OB1 interface. Specifically, W2547C (human W2626C) in the helical domain may disrupt a direct hydrogen bond to the backbone of W2646 (human W2725) in the OB1 domain. Similarly, Y2647C (human Y2726C) in OB1 may disrupt the hydrogen bonding to H2544 (human H2623) and E2584 (human E2663) in the helical domain. In addition, the T2643R mutation (human T2722R) in OB1 potentially leads to a steric clash with W2547 in the helical domain. Furthermore, OB1 residue D2644 (human D2723) associated with the human D2723H, D2723G, and D2723A mutations is hydrogen bonded to the OB1 residues S2591 (human S2670), the L2593 (human I2672) and the W2646 (human W2725), suggesting that these mutations may disrupt the helical-OB1 interface. Separately, disruption of the loop residue G2669 (human G2748) by the G2748D mutation may result in steric clashes with the DSS1 residues Q55 and A58, and destabilization of the DSS1-OB1 interaction. In region 3, the R2705 (human R2784) in OB1 forms hydrogen bonds with the side chains of N2702 (human N2781) in OB1, S2728 (human S2807) in OB2, and the D41 side chain of DSS1. This critical residue serves to clamp three flexible loops through polar interactions. The R2784W and R2784Q mutation of this residue may affect tertiary interactions between OB1 and OB2 and the interaction between BRCA2 and DSS1. In region 4, R2971 (human R3052) in OB3 is positioned at the junction of the OB2 and OB3 folds (Fig. 2D). Remarkably, this residue seems to interact with residues E2921 (human E3002) and V2746 (human A2825) in OB2 and with D2992 (human D3073) and N3042 (human N3124) in OB3. Thus, the R3052W, E3002K, D3073G, and N3124I pathogenic mutations likely disrupt this interaction network and the structure and stability of the OB2-OB3 interface.

Figure 3.

Diagram of the 3-dimensional structure of the BRCA2 DBD and the 4 regions containing HDR-inactivating variants. Variants classified as pathogenic by both the multifactorial likelihood model and the HDR assay are shown in red. Variants classified as pathogenic by the HDR assay but unclassified by the multifactorial likelihood model because of lack of genetic evidence are shown in purple. Helical domain of BRCA2 (A); Helical-OB1-DSS1 interface (B); OB1-OB2-DSS1 interface (C); and OB2-OB3 interface (D). The mouse BRCA2-DSS1-ssDNA complex structure (pdb ID: 1MJE) served as the model for analysis of the variants. All figures were generated using Pymol Molecular Graphics version 1.3.

Figure 3.

Diagram of the 3-dimensional structure of the BRCA2 DBD and the 4 regions containing HDR-inactivating variants. Variants classified as pathogenic by both the multifactorial likelihood model and the HDR assay are shown in red. Variants classified as pathogenic by the HDR assay but unclassified by the multifactorial likelihood model because of lack of genetic evidence are shown in purple. Helical domain of BRCA2 (A); Helical-OB1-DSS1 interface (B); OB1-OB2-DSS1 interface (C); and OB2-OB3 interface (D). The mouse BRCA2-DSS1-ssDNA complex structure (pdb ID: 1MJE) served as the model for analysis of the variants. All figures were generated using Pymol Molecular Graphics version 1.3.

Close modal

Because of the paucity of genetic data that can be used for the classification of BRCA1 and BRCA2 VUS (23), alternative methods that can predict the clinical relevance of these mutations are needed. Because disruption or depletion of BRCA2 causes defective HDR-mediated DNA DSB repair, we evaluated the use of a cell-based HDR assay for classification of VUS in the DBD of BRCA2. Once established, we asked whether the HDR-based functional assay could independently validate 31 variants previously classified as pathogenic or nonpathogenic. The HDR assay clearly segregated these mutations into 2 groups, with all of the known pathogenic mutations defined by genetic studies displaying impaired HDR activity, suggesting that the assay could be used for characterization of germline variants in the BRCA2 DBD.

In an effort to gain insight into the clinical relevance of variants that could not be classified by genetic studies due to the limited availability of family data, we developed a model for estimating the probability of pathogenicity of BRCA2 DBD VUS. This classifier was based on the HDR assay results for known pathologic and nonpathologic variants. Although the number of variants in each of the known classes was limited, the degree of separation observed in the HDR assay between the 2 classes was very large relative to the variability of the assay. Therefore, we were able to estimate with sufficient accuracy the parameters used in the model to discriminate between pathogenic and nonpathogenic variants. To test the robustness of the model, we also conducted a sensitivity analysis that perturbed the parameters relative to the SE of the measurements and verified that the overall conclusions were unchanged. The resulting estimates of the conditional probabilities that a variant is pathogenic accurately classified the 13 known pathogenic and 18 known nonpathogenic mutations. Furthermore, using probability of pathogenicity thresholds of 99% and 1%, the model predicted 18 other VUS as pathogenic and 8 as nonpathogenic. In total, the multifactorial likelihood model and the HDR assay-based probability model used in this study reclassified 21 VUS from the BRCA2 DBD as pathogenic or likely pathogenic. Given that only 10 BRCA2 DBD missense variants were previously identified as pathogenic, this study represents a substantial improvement in VUS classification. The large proportion of predicted pathogenic variants in the DBD is consistent with previous findings implicating the domain in the tumor suppression function of BRCA2. Indeed, Easton and colleagues estimated that 35% of missense mutations in the BRCA2 DBD were pathogenic (6). Overall, our report suggests that in the absence of detailed genetic background needed for the multifactorial likelihood model, the HDR assay presents a viable alternative for assessing the clinical relevance of VUS in the BRCA2 DBD.

While the HDR-based classifier effectively identified 18 VUS as pathogenic and suggested that 8 other VUS were nonpathogenic, 7 VUS displayed intermediate HDR activity. Whether additional VUS in the BRCA2 DBD also exhibit intermediate HDR activity remains to be determined through additional studies. Because of limited family data, it is not known if these VUS with intermediate activity influence breast cancer risk. To investigate whether intermediate function equates to intermediate risk, it may be necessary to identify all intermediate function VUS in the DBD, collect family information for all of these variants, and combine the information for the entire group to have sufficient data for estimation of risk and/or penetrance. The recently established ENIGMA consortium (Evidence-based Network for the Interpretation of Germline Mutant Alleles; ref. 30) for large-scale collaborative studies of BRCA1 and BRCA2 sequence variants has begun this process.

We also evaluated the influence of Align GVGD–based prior probabilities of pathogenicity rather than equal prior probabilities for each variant on the posterior probability of pathogenicity. The results were similar, with a correlation of 0.987 (95% CI: 0.979–0.992) between the posterior probabilities from the 2 models. However, VUS with high (0.81) and low (0.01) Align GVGD-based prior probabilities had slightly higher and lower posterior probabilities of pathogenicity, respectively, than those obtained from the equal prior probability model (Supplementary Table S2). Importantly, the 7 VUS with intermediate HDR activity remained unclassified with either prior probability model (1% and 99% thresholds).

To ascertain if a structure–function correlation exists for the pathogenic missense mutations, we mapped the variants on the predicted 3-dimensional structure of a murine BRCA2-DSS1 complex. We found that many of the parental residues of the pathogenic missense mutants were located at the interdomain interfaces of the BRCA2 DBD. These data suggest that the preservation of these interdomain interactions are important for retention of BRCA2 HDR activity. However, 5 nonpathogenic variants were also found to be located at interdomain interfaces (Supplementary Fig. S2). Specifically, R2787H (mouse R2708) is at the OB1-DSS1 interface, S3020C (mouse F2939) is at the OB2-OB3 and OB2-DNA interface, and T3013I (mouse V2932) is at the OB1-OB2 interface. Furthermore, the S2807L (mouse S2728) and G2813E (mouse G2734) variants, classified as nonpathogenic by the HDR assay but not by the multifactorial likelihood model because of lack of genetic evidence, are both at the OB1-OB2 interface. In addition, 4 of 7 intermediate activity variants were located at interfaces (Supplementary Fig. S2). Specifically, P2800R/S (mouse P2721) and G2812E (mouse G2733) are located at the OB1-OB2 interface and E3002D (mouse E2921) is at the OB2-OB3 interface. However, P2800R/S and G2812E are solvent exposed and can be accommodated in the structure. In addition, while the pathogenic E3002K substitutes a negatively charged amino acid with a positively charged amino acid and is predicted to induce a dramatic change in the stability of the structure, the E3002D variant at the same residue maintains the negative charge and only affects hydrogen bonds with surrounding residues. The remaining intermediate activity VUS R2842C/L (mouse R2763), T3033I (mouse T2952), and Y3092S (mouse Y3011S) are not located at interdomain interfaces but may influence DNA binding. Overall, these results show that not all residues at interdomain interfaces are required for BRCA2 HDR activity. Mutation of certain interface residues may not affect the function because the specific change can be tolerated within the DBD structure or because the parental residues are not important for the maintenance of interdomain interactions. These findings suggest that the location of mutant residues alone is insufficient for predicting the effects of VUS on protein function or on cancer risk.

We chose a cell-based HDR assay rather than an in vitro biochemical assay to ensure that the in vivo effects of VUS on HDR activity were measured. In addition, we based the HDR assay on the use of the full-length BRCA2 protein, to understand how BRCA2 VUS affect the HDR function of BRCA2 in the context of all other functional domains. Because several functional properties have been attributed to BRCA2, it is important to account for the effects of these characteristics, some of which have unknown relevance to cancer risk on BRCA2 HDR activity. The use of full-length protein also limits the potential for misinterpretation of effects observed in partial proteins. Importantly, even though the assay is conducted using a full-length protein, only variants in domains directed implicated in HDR repair should be evaluated. Here, we focused on the BRCA2 DBD, but it is possible that variants in other domains, such as the N-terminal PALB2-interacting domain, can also be effectively evaluated. It is also possible that there are other functions of BRCA2 associated with the DBD of BRCA2 that are not detected by the HDR assay. If these alternative functions influence breast cancer risk, then VUS with no effect or limited influence on HDR activity may actually influence breast cancer risk. For this reason, VUS with wild-type HDR activity should be treated with caution.

F.J. Couch has ownership interest (including patents) in Myriad Genetics Laboratories Inc. No potential conflicts of interest were disclosed by the other authors.

Conception and design: L. Guidugli, F.J. Couch

Development of methodology: L. Guidugli, V.S. Pankratz, F.J. Couch

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): L. Guidugli, C.A. Erding, C. Engel, R.K. Schmutzler, S. Domchek, K. Nathanson, P. Radice, C.F. Singer, P.N. Tonin, N.M. Lindor, F.J. Couch

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): L. Guidugli, V.S. Pankratz, N. Singh, J. Thompson, C.F. Singer, D.E. Goldgar

Writing, review, and/or revision of the manuscript: L. Guidugli, V.S. Pankratz, N. Singh, J. Thompson, C. Engel, R.K. Schmutzler, S. Domchek, K. Nathanson, P. Radice, P.N. Tonin, D.E. Goldgar, F.J. Couch

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): C.A. Erding, C. Engel, R.K. Schmutzler, S. Domchek, F.J. Couch

Study supervision: F.J. Couch

This research was supported by an NIH specialized program of research excellence (SPORE) in breast cancer award to the Mayo Clinic (CA116201) NIH grants R01 CA116167 and R01 CA128978. L. Guidugli is supported by a fellowship from the Komen Foundation for the Cure. S. Domchek is supported by funding from the Komen Foundation for the Cure. F.J. Couch and K. Nathanson are supported by funding from the Breast Cancer Research Foundation.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Moynahan
ME
,
Pierce
AJ
,
Jasin
M
. 
BRCA2 is required for homology-directed repair of chromosomal breaks
.
Mol Cell
2001
;
7
:
263
72
.
2.
Jensen
RB
,
Carreira
A
,
Kowalczykowski
SC
. 
Purified human BRCA2 stimulates RAD51-mediated recombination
.
Nature
2010
;
467
:
678
83
.
3.
Easton
DF
,
Steele
L
,
Fields
P
,
Ormiston
W
,
Averill
D
,
Daly
PA
, et al
Cancer risks in two large breast cancer families linked to BRCA2 on chromosome 13q12-13
.
Am J Hum Genet
1997
;
61
:
120
8
.
4.
Antoniou
AC
,
Cunningham
AP
,
Peto
J
,
Evans
DG
,
Lalloo
F
,
Narod
SA
, et al
The BOADICEA model of genetic susceptibility to breast and ovarian cancers: updates and extensions
.
Br J Cancer
2008
;
98
:
1457
66
.
5.
The Breast Cancer Information Core
. 
An open access on-line breast cancer mutation data base
. 
1996
;
[cited 1996]. Available from
: http://research.nhgri.nih.gov/bic/
6.
Easton
DF
,
Deffenbaugh
AM
,
Pruss
D
,
Frye
C
,
Wenstrup
RJ
,
Allen-Brady
K
, et al
A systematic genetic assessment of 1,433 sequence variants of unknown clinical significance in the BRCA1 and BRCA2 breast cancer-predisposition genes
.
Am J Hum Genet
2007
;
81
:
873
83
.
7.
Goldgar
DE
,
Easton
DF
,
Deffenbaugh
AM
,
Monteiro
AN
,
Tavtigian
SV
,
Couch
FJ
. 
Integrated evaluation of DNA sequence variants of unknown clinical significance: application to BRCA1 and BRCA2
.
Am J Hum Genet
2004
;
75
:
535
44
.
8.
Chenevix-Trench
G
,
Healey
S
,
Lakhani
S
,
Waring
P
,
Cummings
M
,
Brinkworth
R
, et al
Genetic and histopathologic evaluation of BRCA1 and BRCA2 DNA sequence variants of unknown clinical significance
.
Cancer Res
2006
;
66
:
2019
27
.
9.
Walker
LC
,
Whiley
PJ
,
Couch
FJ
,
Farrugia
DJ
,
Healey
S
,
Eccles
DM
, et al
Detection of splicing aberrations caused by BRCA1 and BRCA2 sequence variants encoding missense substitutions: implications for prediction of pathogenicity
.
Hum Mutat
2010
;
31
:
E1484
505
.
10.
Tavtigian
SV
,
Deffenbaugh
AM
,
Yin
L
,
Judkins
T
,
Scholl
T
,
Samollow
PB
, et al
Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral
.
J Med Genet
2006
;
43
:
295
305
.
11.
Tavtigian
SV
,
Greenblatt
MS
,
Lesueur
F
,
Byrnes
GB
. 
In silico analysis of missense substitutions using sequence-alignment based methods
.
Hum Mutat
2008
;
29
:
1327
36
.
12.
Plon
SE
,
Eccles
DM
,
Easton
D
,
Foulkes
WD
,
Genuardi
M
,
Greenblatt
MS
, et al
Sequence variant classification and reporting: recommendations for improving the interpretation of cancer susceptibility genetic test results
.
Hum Mutat
2008
;
29
:
1282
91
.
13.
Farrugia
DJ
,
Agarwal
MK
,
Pankratz
VS
,
Deffenbaugh
AM
,
Pruss
D
,
Frye
C
, et al
Functional assays for classification of BRCA2 variants of uncertain significance
.
Cancer Res
2008
;
68
:
3523
31
.
14.
Lee
MS
,
Green
R
,
Marsillac
SM
,
Coquelle
N
,
Williams
RS
,
Yeung
T
, et al
Comprehensive analysis of missense variations in the BRCT domain of BRCA1 by structural and functional assays
.
Cancer Res
2010
;
70
:
4880
90
.
15.
Wu
K
,
Hinson
SR
,
Ohashi
A
,
Farrugia
D
,
Wendt
P
,
Tavtigian
SV
, et al
Functional evaluation and cancer risk assessment of BRCA2 unclassified variants
.
Cancer Res
2005
;
65
:
417
26
.
16.
Couch
FJ
,
Rasmussen
LJ
,
Hofstra
R
,
Monteiro
AN
,
Greenblatt
MS
,
de Wind
N
. 
Assessment of functional effects of unclassified genetic variants
.
Hum Mutat
2008
;
29
:
1314
26
.
17.
Tavtigian
SV
,
Byrnes
GB
,
Goldgar
DE
,
Thomas
A
. 
Classification of rare missense substitutions, using risk surfaces, with genetic- and molecular-epidemiology applications
.
Hum Mutat
2008
;
29
:
1342
54
.
18.
Vallee
MP
,
Francy
TC
,
Judkins
MK
,
Babikyan
D
,
Lesueur
F
,
Gammon
A
, et al
Classification of missense substitutions in the BRCA genes: a database dedicated to Ex-UVs
.
Hum Mutat
2012
;
33
:
22
8
.
19.
Kraakman-van der Zwet
M
,
Overkamp
WJ
,
van Lange
RE
,
Essers
J
,
van Duijn-Goedhart
A
,
Wiggers
I
, et al
Brca2 (XRCC11) deficiency results in radioresistant DNA synthesis and a higher frequency of spontaneous deletions
.
Mol Cell Biol
2002
;
22
:
669
79
.
20.
Fawcett
T
. 
An introduction to ROC analysis
.
Pattern Recognit Lett
2006
;
27
:
861
74
.
21.
McLachlan
GJ
. 
Discriminant analysis and statistical pattern recognition
.
New York, NY
:
Wiley
; 
1992
.
22.
Fraley
C
,
Raftery
AE
. 
Model-based clustering, discriminant analysis, and density estimation
.
J Am Stat Assoc
2002
;
97
:
611
31
.
23.
Lindor
NM
,
Guidugli
L
,
Wang
X
,
Vallee
MP
,
Monteiro
AN
,
Tavtigian
S
, et al
A review of a multifactorial probability-based model for classification of BRCA1 and BRCA2 variants of uncertain significance (VUS)
.
Hum Mutat
2012
;
33
:
8
21
.
24.
Freedman
ML
,
Penney
KL
,
Stram
DO
,
Le Marchand
L
,
Hirschhorn
JN
,
Kolonel
LN
, et al
Common variation in BRCA2 and breast cancer risk: a haplotype-based analysis in the Multiethnic Cohort
.
Hum Mol Genet
2004
;
13
:
2431
41
.
25.
Deffenbaugh
AM
,
Frank
TS
,
Hoffman
M
,
Cannon-Albright
L
,
Neuhausen
SL
. 
Characterization of common BRCA1 and BRCA2 variants
.
Genet Test
2002
;
6
:
119
21
.
26.
Wagner
TM
,
Hirtenlehner
K
,
Shen
P
,
Moeslinger
R
,
Muhr
D
,
Fleischmann
E
, et al
Global sequence diversity of BRCA2: analysis of 71 breast cancer families and 95 control individuals of worldwide populations
.
Hum Mol Genet
1999
;
8
:
413
23
.
27.
Yeo
G
,
Burge
CB
. 
Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals
.
J Comput Biol
2004
;
11
:
377
94
.
28.
Desmet
FO
,
Hamroun
D
,
Lalande
M
,
Collod-Beroud
G
,
Claustres
M
,
Beroud
C
. 
Human Splicing Finder: an online bioinformatics tool to predict splicing signals
.
Nucleic Acids Res
2009
;
37
:
e67
.
29.
Yang
H
,
Jeffrey
PD
,
Miller
J
,
Kinnucan
E
,
Sun
Y
,
Thoma
NH
, et al
BRCA2 function in DNA binding and recombination from a BRCA2-DSS1-ssDNA structure
.
Science
2002
;
297
:
1837
48
.
30.
Spurdle
AB
,
Healey
S
,
Devereau
A
,
Hogervorst
FB
,
Monteiro
AN
,
Nathanson
KL
, et al
ENIGMA–evidence-based network for the interpretation of germline mutant alleles: an international initiative to evaluate risk and clinical significance associated with sequence variation in BRCA1 and BRCA2 genes
.
Hum Mutat
2012
;
33
:
2
7
.