The functional impact of BRCA1 BRCT domain variants using multiplexed DNA double-strand break repair assays

Pathogenic variants in BRCA1 are associated with a greatly increased risk of hereditary breast and ovarian cancer (HBOC). With the increased availability and affordability of genetic testing, many individuals have been identified with BRCA1 variants of uncertain significance (VUSs), which are individually detected in the population too infrequently to ascertain a clinical risk. Functional assays can be used to experimentally assess the effects of these variants. In this study, we used multiplexed DNA repair assays of variants in the BRCA1 carboxyl terminus to functionally characterize 2,271 variants for homology-directed repair function (HDR) and 1,427 variants for cisplatin resistance (CR). We found a high level of consistent results (Pearson’s r = 0.74) in the two multiplexed functional assays with non-functional variants located within regions of the BRCA1 protein necessary for its tumor suppression activity. In addition, functional categorizations of variants tested in the multiplex HDR and CR assays correlated with known clinical significance and with other functional assays for BRCA1 (Pearson’s r = 0.53 to 0.71). The results of the multiplex HDR and CR assays are useful resources for characterizing large numbers of BRCA1 VUSs.     

Molding BRCA2 function through its interacting partners

The role of the tumor suppressor BRCA2 has been shaped over 2 decades thanks to the discovery of its protein and nucleic acid partners, biochemical and structural studies of the protein, and the functional evaluation of germline variants identified in breast cancer patients. Yet, the pathogenic and functional effect of many germline mutations in BRCA2 remains undetermined, and the heterogeneity of BRCA2-associated tumors challenges the identification of causative variants that drive tumorigenesis.In this review, we propose an overview of the established and emerging interacting partners and functional pathways attributed to BRCA2, and we speculate on how variants altering these functions may contribute to cancer susceptibility.     

A new bioinformatics tool to help assess the significance of <Emphasis Type="Italic">BRCA1</Emphasis> variants

Background

Germline pathogenic variants in the breast cancer type 1 susceptibility gene BRCA1 are associated with a 60% lifetime risk for breast and ovarian cancer. This overall risk estimate is for all BRCA1 variants; obviously, not all variants confer the same risk of developing a disease. In cancer patients, loss of BRCA1 function in tumor tissue has been associated with an increased sensitivity to platinum agents and to poly-(ADP-ribose) polymerase (PARP) inhibitors. For clinical management of both at-risk individuals and cancer patients, it would be important that each identified genetic variant be associated with clinical significance. Unfortunately for the vast majority of variants, the clinical impact is unknown. The availability of results from studies assessing the impact of variants on protein function may provide insight of crucial importance.

Results and conclusion

We have collected, curated, and structured the molecular and cellular phenotypic impact of 3654 distinct BRCA1 variants. The data was modeled in triple format, using the variant as a subject, the studied function as the object, and a predicate describing the relation between the two. Each annotation is supported by a fully traceable evidence. The data was captured using standard ontologies to ensure consistency, and enhance searchability and interoperability. We have assessed the extent to which functional defects at the molecular and cellular levels correlate with the clinical interpretation of variants by ClinVar submitters. Approximately 30% of the ClinVar BRCA1 missense variants have some molecular or cellular assay available in the literature. Pathogenic variants (as assigned by ClinVar) have at least some significant functional defect in 94% of testable cases. For benign variants, 77% of ClinVar benign variants, for which neXtProt Cancer variant portal has data, shows either no or mild experimental functional defects. While this does not provide evidence for clinical interpretation of variants, it may provide some guidance for variants of unknown significance, in the absence of more reliable data.The neXtProt Cancer variant portal (https://www.nextprot.org/portals/breast-cancer) contains over 6300 observations at the molecular and/or cellular level for BRCA1 variants.

     

DNA methylation profiling to assess pathogenicity of BRCA1 unclassified variants in breast cancer

Germline pathogenic mutations in BRCA1 increase risk of developing breast cancer. Screening for mutations in BRCA1 frequently identifies sequence variants of unknown pathogenicity and recent work has aimed to develop methods for determining pathogenicity. We previously observed that tumor DNA methylation can differentiate BRCA1-mutated from BRCA1-wild type tumors. We hypothesized that we could predict pathogenicity of variants based on DNA methylation profiles of tumors that had arisen in carriers of unclassified variants. We selected 150 FFPE breast tumor DNA samples [47 BRCA1 pathogenic mutation carriers, 65 BRCAx (BRCA1-wild type), 38 BRCA1 test variants] and analyzed a subset (n=54) using the Illumina 450K methylation platform, using the remaining samples for bisulphite pyrosequencing validation. Three validated markers (BACH2, C8orf31, and LOC654342) were combined with sequence bioinformatics in a model to predict pathogenicity of 27 variants (independent test set).  Predictions were compared with standard multifactorial likelihood analysis. Prediction was consistent for c.5194-12G>A (IVS 19-12 G>A) (P>0.99); 13 variants were considered not pathogenic or likely not pathogenic using both approaches. We conclude that tumor DNA methylation data alone has potential to be used in prediction of BRCA1 variant pathogenicity but is not independent of estrogen receptor status and grade, which are used in current multifactorial models to predict pathogenicity.     

BRCA1 and BRCA2 missense variants of high and low clinical significance influence lymphoblastoid cell line post-irradiation gene expression

The functional consequences of missense variants in disease genes are difficult to predict. We assessed if gene expression profiles could distinguish between BRCA1 or BRCA2 pathogenic truncating and missense mutation carriers and familial breast cancer cases whose disease was not attributable to BRCA1 or BRCA2 mutations (BRCAX cases). 72 cell lines from affected women in high-risk breast ovarian families were assayed after exposure to ionising irradiation, including 23 BRCA1 carriers, 22 BRCA2 carriers, and 27 BRCAX individuals. A subset of 10 BRCAX individuals carried rare BRCA1/2 sequence variants considered to be of low clinical significance (LCS). BRCA1 and BRCA2 mutation carriers had similar expression profiles, with some subclustering of missense mutation carriers. The majority of BRCAX individuals formed a distinct cluster, but BRCAX individuals with LCS variants had expression profiles similar to BRCA1/2 mutation carriers. Gaussian Process Classifier predicted BRCA1, BRCA2 and BRCAX status, with a maximum of 62% accuracy, and prediction accuracy decreased with inclusion of BRCAX samples carrying an LCS variant, and inclusion of pathogenic missense carriers. Similarly, prediction of mutation status with gene lists derived using Support Vector Machines was good for BRCAX samples without an LCS variant (82–94%), poor for BRCAX with an LCS (40–50%), and improved for pathogenic BRCA1/2 mutation carriers when the gene list used for prediction was appropriate to mutation effect being tested (71–100%). This study indicates that mutation effect, and presence of rare variants possibly associated with a low risk of cancer, must be considered in the development of array-based assays of variant pathogenicity.     

A systematic genetic assessment of 1,433 sequence variants of unknown clinical significance in the BRCA1 and BRCA2 breast cancer-predisposition genes

Mutation screening of the breast and ovarian cancer-predisposition genes BRCA1 and BRCA2 is becoming an increasingly important part of clinical practice. Classification of rare nontruncating sequence variants in these genes is problematic, because it is not known whether these subtle changes alter function sufficiently to predispose cells to cancer development. Using data from the Myriad Genetic Laboratories database of nearly 70,000 full-sequence tests, we assessed the clinical significance of 1,433 sequence variants of unknown significance (VUSs) in the BRCA genes. Three independent measures were employed in the assessment: co-occurrence in trans of a VUS with known deleterious mutations; detailed analysis, by logistic regression, of personal and family history of cancer in VUS-carrying probands; and, in a subset of probands, an analysis of cosegregation with disease in pedigrees. For each of these factors, a likelihood ratio was computed under the hypothesis that the VUSs were equivalent to an "average" deleterious mutation, compared with neutral, with respect to risk. The likelihood ratios derived from each component were combined to provide an overall assessment for each VUS. A total of 133 VUSs had odds of at least 100 : 1 in favor of neutrality with respect to risk, whereas 43 had odds of at least 20 : 1 in favor of being deleterious. VUSs with evidence in favor of causality were those that were predicted to affect splicing, fell at positions that are highly conserved among BRCA orthologs, and were more likely to be located in specific domains of the proteins. In addition to their utility for improved genetics counseling of patients and their families, the global assessment reported here will be invaluable for validation of functional assays, structural models, and in silico analyses.     

Detection of New Point Mutations of BRCA1 and BRCA2 in Breast Cancer Patients

This study included 20 selected female patients with breast cancer, 30 of their female relatives (sisters and daughters), and 10 healthy females as a control group. Genomic DNA was extracted from peripheral blood lymphocytes of all the subjects, and the polymerase chain reaction was carried out using specific primers for BRCA1 (exons 2 and 8) and BRCA2 (exons 9, 11, and 21). The mutations were detected using a single-strand conformation polymorphism assay and heteroduplex analysis. Finally, the sample variants and their controls were sequenced. Mutations were detected in 44% of the study population, with 18% found in the BRCA1 gene and 26% attributed to BRCA2. Five sequence variants were identified, including two frameshift mutations, one nonsense mutation, and two missense mutations. Therefore, we conclude that germline mutations in two major genes, BRCA1 and BRCA2, may have an important influence on the predisposition and development of familial breast cancer.     

Functional profiling of nucleotide Excision repair in breast cancer

Homologous recombination deficiency conferred by alterations in BRCA1 or BRCA2 are common in breast tumors and can drive sensitivity to platinum chemotherapy and PARP inhibitors. Alterations in nucleotide excision repair (NER) activity can also impact sensitivity to DNA damaging agents, but NER activity in breast cancer has been poorly characterized. Here, we apply a novel immunofluorescence-based cellular NER assay to screen a large panel of breast epithelial and cancer cell lines. Although the majority of breast cancer models are NER proficient, we identify an example of a breast cancer cell line with profound NER deficiency. We show that NER deficiency in this model is driven by epigenetic silencing of the ERCC4 gene, leading to lack of expression of the NER nuclease XPF, and that ERCC4 methylation is also strongly correlated with ERCC4 mRNA and XPF protein expression in primary breast tumors. Re-expression of XPF in the ERCC4-deficient breast cancer rescues NER deficiency and cisplatin sensitivity, but does not impact PARP inhibitor sensitivity. These findings demonstrate the potential to use functional assays to identify novel mechanisms of DNA repair deficiency and nominate NER deficiency as a platinum sensitivity biomarker in breast cancer.     

A new assay for functional screening of BRCA2 linker region mutations identifies variants that alter chemoresistance to cisplatin

Variants of unknown significance (VUS) complicate the assignment of risk to new DNA sequence variants found in at-risk populations. This study focused on the poorly studied linker region of the cancer-associated BRCA2 protein encoded by exons twelve through fourteen of BRCA2. To develop a new method to characterize VUS in this region of BRCA2, we first chose to study 4 reported VUS occurring on evolutionarily conserved residues within the linker region. To determine if these VUS represent neutral changes or if they impact the function of the BRCA2 protein, we stably transfected expression plasmids encoding wild-type or each mutant peptide into T47D breast cancer cells, which are wild-type for BRCA2. Four mutant peptide expressing cell lines and a wild-type linker region expressing cell line next were studied by challenging transfected cell lines with the DNA crosslinking compound cisplatin (10 μM) for 5 days. Expression of the wild-type linker region and certain mutant linker peptides (N2452D and I2285V) decreased apoptosis (as demonstrated by cell death detection assay) in transfected cell lines, indicating that the linker region peptide directly or indirectly affects the DNA damage repair pathway. By determining the cell survival and assaying the apoptotic index of treated cell lines, one could potentially use this screen to determine that a particular VUS has a functional impact on BRCA2 function, and hence is of functional significance. We conclude that this method is useful for screening the effect of linker region VUS on BRCA2 function, and to identify mutations for further testing. We also conclude that mutations in the linker region may have heretofore unappreciated roles in BRCA2 function.     

Rare BRCA1 haplotypes including 3′UTR SNPs associated with breast cancer risk

Genetic markers identifying women at an increased risk of developing breast cancer exist, yet the majority of inherited risk remains elusive. While numerous BRCA1 coding sequence mutations are associated with breast cancer risk, BRCA1 mutations account for less then 5% of breast cancer risk. Since 3′ untranslated region (3′UTR) polymorphisms disrupting microRNA (miRNA) binding can be functional and can act as genetic markers of cancer risk, we tested the hypothesis that such polymorphisms in the 3′UTR of BRCA1 and haplotypes containing these functional polymorphisms may be associated with breast cancer risk. We sequenced the BRCA1 3′UTR from breast cancer patients to identify miRNA disrupting polymorphisms. We further evaluated haplotypes of this region including the identified 3′UTR variants in a large population of controls and breast cancer patients (n = 221) with known breast cancer subtypes and ethnicities. We identified three 3′UTR variants in BRCA1 that are polymorphic in breast cancer populations, and haplotype analysis including these variants revealed that breast cancer patients harbor five rare haplotypes not generally found among controls (9.50% for breast cancer chromosomes, 0.11% for control chromosomes, p = 0.0001). Three of these rare haplotypes contain the rs8176318 BRCA1 3′UTR functional variant. These haplotypes are not biomarkers for BRCA1 coding mutations, as they are found rarely in BRCA1 mutant breast cancer patients (1/129 patients = 0.78%). These rare BRCA1 haplotypes and 3′UTR SNPs may represent new genetic markers of breast cancer risk.Key words: BRCA1, haplotype, microRNA, SNP, 3′UTR, breast cancer, triple negative breast cancer     

Whole Exome Sequencing Suggests Much of Non-BRCA1/BRCA2 Familial Breast Cancer Is Due to Moderate and Low Penetrance Susceptibility Alleles

The identification of the two most prevalent susceptibility genes in breast cancer, BRCA1 and BRCA2, was the beginning of a sustained effort to uncover new genes explaining the missing heritability in this disease. Today, additional high, moderate and low penetrance genes have been identified in breast cancer, such as P53, PTEN, STK11, PALB2 or ATM, globally accounting for around 35 percent of the familial cases. In the present study we used massively parallel sequencing to analyze 7 BRCA1/BRCA2 negative families, each having at least 6 affected women with breast cancer (between 6 and 10) diagnosed under the age of 60 across generations. After extensive filtering, Sanger sequencing validation and co-segregation studies, variants were prioritized through either control-population studies, including up to 750 healthy individuals, or case-control assays comprising approximately 5300 samples. As a result, a known moderate susceptibility indel variant (CHEK2 1100delC) and a catalogue of 11 rare variants presenting signs of association with breast cancer were identified. All the affected genes are involved in important cellular mechanisms like DNA repair, cell proliferation and survival or cell cycle regulation. This study highlights the need to investigate the role of rare variants in familial cancer development by means of novel high throughput analysis strategies optimized for genetically heterogeneous scenarios. Even considering the intrinsic limitations of exome resequencing studies, our findings support the hypothesis that the majority of non-BRCA1/BRCA2 breast cancer families might be explained by the action of moderate and/or low penetrance susceptibility alleles.     

Targeted Sequencing of the Mitochondrial Genome of Women at High Risk of Breast Cancer without Detectable Mutations in BRCA1/2

Breast Cancer is a complex multifactorial disease for which high-penetrance mutations have been identified. Approaches used to date have identified genomic features explaining about 50% of breast cancer heritability. A number of low- to medium penetrance alleles (per-allele odds ratio < 1.5 and 4.0, respectively) have been identified, suggesting that the remaining heritability is likely to be explained by the cumulative effect of such alleles and/or by rare high-penetrance alleles. Relatively few studies have specifically explored the mitochondrial genome for variants potentially implicated in breast cancer risk. For these reasons, we propose an exploration of the variability of the mitochondrial genome in individuals diagnosed with breast cancer, having a positive breast cancer family history but testing negative for BRCA1/2 pathogenic mutations. We sequenced the mitochondrial genome of 436 index breast cancer cases from the GENESIS study. As expected, no pathogenic genomic pattern common to the 436 women included in our study was observed. The mitochondrial genes MT-ATP6 and MT-CYB were observed to carry the highest number of variants in the study. The proteins encoded by these genes are involved in the structure of the mitochondrial respiration chain, and variants in these genes may impact reactive oxygen species production contributing to carcinogenesis. More functional and epidemiological studies are needed to further investigate to what extent variants identified may influence familial breast cancer risk.     

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient’s RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZ_EI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.     

The BRCA1 Variant p.Ser36Tyr Abrogates BRCA1 Protein Function and Potentially Confers a Moderate Risk of Breast Cancer

The identification of variants of unknown clinical significance (VUS) in the BRCA1 gene complicates genetic counselling and causes additional anxiety to carriers. In silico approaches currently used for VUS pathogenicity assessment are predictive and often produce conflicting data. Furthermore, functional assays are either domain or function specific, thus they do not examine the entire spectrum of BRCA1 functions and interpretation of individual assay results can be misleading. PolyPhen algorithm predicted that the BRCA1 p.Ser36Tyr VUS identified in the Cypriot population was damaging, whereas Align-GVGD predicted that it was possibly of no significance. In addition the BRCA1 p.Ser36Tyr variant was found to be associated with increased risk (OR = 3.47, 95% CI 1.13-10.67, P = 0.02) in a single case-control series of 1174 cases and 1109 controls. We describe a cellular system for examining the function of exogenous full-length BRCA1 and for classifying VUS. We achieved strong protein expression of full-length BRCA1 in transiently transfected HEK293T cells. The p.Ser36Tyr VUS exhibited low protein expression similar to the known pathogenic variant p.Cys61Gly. Co-precipitation analysis further demonstrated that it has a reduced ability to interact with BARD1. Further, co-precipitation analysis of nuclear and cytosolic extracts as well as immunofluorescence studies showed that a high proportion of the p.Ser36Tyr variant is withheld in the cytoplasm contrary to wild type protein. In addition the ability of p.Ser36Tyr to co-localize with conjugated ubiquitin foci in the nuclei of S-phase synchronized cells following genotoxic stress with hydroxyurea is impaired at more pronounced levels than that of the p.Cys61Gly pathogenic variant. The p.Ser36Tyr variant demonstrates abrogated function, and based on epidemiological, genetic, and clinical data we conclude that the p.Ser36Tyr variant is probably associated with a moderate breast cancer risk.     

Germline Variation in Cancer-Susceptibility Genes in a Healthy,Ancestrally Diverse Cohort: Implications for Individual Genome Sequencing

Technological advances coupled with decreasing costs are bringing whole genome and whole exome sequencing closer to routine clinical use. One of the hurdles to clinical implementation is the high number of variants of unknown significance. For cancer-susceptibility genes, the difficulty in interpreting the clinical relevance of the genomic variants is compounded by the fact that most of what is known about these variants comes from the study of highly selected populations, such as cancer patients or individuals with a family history of cancer. The genetic variation in known cancer-susceptibility genes in the general population has not been well characterized to date. To address this gap, we profiled the nonsynonymous genomic variation in 158 genes causally implicated in carcinogenesis using high-quality whole genome sequences from an ancestrally diverse cohort of 681 healthy individuals. We found that all individuals carry multiple variants that may impact cancer susceptibility, with an average of 68 variants per individual. Of the 2,688 allelic variants identified within the cohort, most are very rare, with 75% found in only 1 or 2 individuals in our population. Allele frequencies vary between ancestral groups, and there are 21 variants for which the minor allele in one population is the major allele in another. Detailed analysis of a selected subset of 5 clinically important cancer genes, BRCA1, BRCA2, KRAS, TP53, and PTEN, highlights differences between germline variants and reported somatic mutations. The dataset can serve a resource of genetic variation in cancer-susceptibility genes in 6 ancestry groups, an important foundation for the interpretation of cancer risk from personal genome sequences.     

Accurate prediction of the functional significance of single nucleotide polymorphisms and mutations in the ABCA1 gene

The human genome contains an estimated 100,000 to 300,000 DNA variants that alter an amino acid in an encoded protein. However, our ability to predict which of these variants are functionally significant is limited. We used a bioinformatics approach to define the functional significance of genetic variation in the ABCA1 gene, a cholesterol transporter crucial for the metabolism of high density lipoprotein cholesterol. To predict the functional consequence of each coding single nucleotide polymorphism and mutation in this gene, we calculated a substitution position-specific evolutionary conservation score for each variant, which considers site-specific variation among evolutionarily related proteins. To test the bioinformatics predictions experimentally, we evaluated the biochemical consequence of these sequence variants by examining the ability of cell lines stably transfected with the ABCA1 alleles to elicit cholesterol efflux. Our bioinformatics approach correctly predicted the functional impact of greater than 94% of the naturally occurring variants we assessed. The bioinformatics predictions were significantly correlated with the degree of functional impairment of ABCA1 mutations (r² = 0.62, p = 0.0008). These results have allowed us to define the impact of genetic variation on ABCA1 function and to suggest that the in silico evolutionary approach we used may be a useful tool in general for predicting the effects of DNA variation on gene function. In addition, our data suggest that considering patterns of positive selection, along with patterns of negative selection such as evolutionary conservation, may improve our ability to predict the functional effects of amino acid variation.