Functional analysis and clinical classification of 462 germline BRCA2 missense variants affecting the DNA binding domain

Variants of uncertain significance (VUSs) in BRCA2 are a common result of hereditary cancer genetic testing. While more than 4,000 unique VUSs, comprised of missense or intronic variants, have been identified in BRCA2, the few missense variants now classified clinically as pathogenic or likely pathogenic are predominantly located in the region encoding the C-terminal DNA binding domain (DBD). We report on functional evaluation of the influence of 462 BRCA2 missense variants affecting the DBD on DNA repair activity of BRCA2 using a homology-directed DNA double-strand break repair assay. Of these, 137 were functionally abnormal, 313 were functionally normal, and 12 demonstrated intermediate function. Comparisons with other functional studies of BRCA2 missense variants yielded strong correlations. Sequence-based in silico prediction models had high sensitivity, but limited specificity, relative to the homology-directed repair assay. Combining the functional results with clinical and genetic data in an American College of Medical Genetics (ACMG)/Association for Molecular Pathology (AMP)-like variant classification framework from a clinical testing laboratory, after excluding known splicing variants and functionally intermediate variants, classified 431 of 442 (97.5%) missense variants (129 as pathogenic/likely pathogenic and 302 as benign/likely benign). Functionally abnormal variants classified as pathogenic by ACMG/AMP rules were associated with a slightly lower risk of breast cancer (odds ratio [OR] 5.15, 95% confidence interval [CI] 3.43-7.83) than BRCA2 DBD protein truncating variants (OR 8.56, 95% CI 6.03-12.36). Overall, functional studies of BRCA2 variants using validated assays substantially improved the variant classification yield from ACMG/AMP models and are expected to improve clinical management of many individuals found to harbor germline BRCA2 missense VUS.

Saturation genome editing-based functional evaluation and clinical classification of BRCA2 single nucleotide variants

Germline BRCA2 loss-of function (LOF) variants identified by clinical genetic testing predispose to breast, ovarian, prostate and pancreatic cancer. However, variants of uncertain significance (VUS) (n>4000) limit the clinical use of testing results. Thus, there is an urgent need for functional characterization and clinical classification of all BRCA2 variants. Here we report on comprehensive saturation genome editing-based functional characterization of 97% of all possible single nucleotide variants (SNVs) in the BRCA2 DNA Binding Domain hotspot for pathogenic missense variants that is encoded by exons 15 to 26. The assay was based on deep sequence analysis of surviving endogenously targeted haploid cells. A total of 7013 SNVs were characterized as functionally abnormal (n=955), intermediate/uncertain, or functionally normal (n=5224) based on 95% agreement with ClinVar known pathogenic and benign standards. Results were validated relative to batches of nonsense and synonymous variants and variants evaluated using a homology directed repair (HDR) functional assay. Breast cancer case-control association studies showed that pooled SNVs encoding functionally abnormal missense variants were associated with increased risk of breast cancer (odds ratio (OR) 3.89, 95%CI: 2.77-5.51). In addition, 86% of tumors associated with abnormal missense SNVs displayed loss of heterozygosity (LOH), whereas 26% of tumors with normal variants had LOH. The functional data were added to other sources of information in a ClinGen/ACMG/AMP-like model and 700 functionally abnormal SNVs, including 220 missense SNVs, were classified as pathogenic or likely pathogenic, while 4862 functionally normal SNVs, including 3084 missense SNVs, were classified as benign or likely benign. These classified variants can now be used for risk assessment and clinical care of variant carriers and the remaining functional scores can be used directly for clinical classification and interpretation of many additional variants.

Summary

Germline BRCA2 loss-of function (LOF) variants identified by clinical genetic testing predispose to several types of cancer. However, variants of uncertain significance (VUS) limit the clinical use of testing results. Thus, there is an urgent need for functional characterization and clinical classification of all BRCA2 variants to facilitate current and future clinical management of individuals with these variants. Here we show the results from a saturation genome editing (SGE) and functional analysis of all possible single nucleotide variants (SNVs) from exons 15 to 26 that encode the BRCA2 DNA Binding Domain hotspot for pathogenic missense variants. The assay was based on deep sequence analysis of surviving endogenously targeted human haploid HAP1 cells. The assay was calibrated relative to ClinVar known pathogenic and benign missense standards and 95% prevalence thresholds for functionally abnormal and normal variants were identified. Thresholds were validated based on nonsense and synonymous variants. SNVs encoding functionally abnormal missense variants were associated with increased risks of breast and ovarian cancer. The functional assay results were integrated into a ClinGen/ACMG/AMP-like model for clinical classification of the majority of BRCA2 SNVs as pathogenic/likely pathogenic or benign/likely benign. The classified variants can be used for improved clinical management of variant carriers.

BRCA1 frameshift variants leading to extended incorrect protein C termini

Carriers of BRCA1 germline pathogenic variants are at substantially higher risk of developing breast and ovarian cancer than the general population. Accurate identification of at-risk individuals is crucial for risk stratification and the implementation of targeted preventive and therapeutic interventions. Despite significant progress in variant classification efforts, a sizable portion of reported BRCA1 variants remain as variants of uncertain clinical significance (VUSs). Variants leading to premature protein termination and loss of essential functional domains are typically classified as pathogenic. However, the impact of frameshift variants that result in an extended incorrect terminus is not clear. Using validated functional assays, we conducted a systematic functional assessment of 17 previously reported BRCA1 extended incorrect terminus variants (EITs) and concluded that 16 constitute loss-of-function variants. This suggests that most EITs are likely to be pathogenic. However, one variant, c.5578dup, displayed a protein expression level, affinity to known binding partners, and activity in transcription and homologous recombination assays comparable to the wild-type BRCA1 protein. Twenty-three additional carriers of c.5578dup were identified at a US clinical diagnostic lab and assessed using a family history likelihood model providing, in combination with the functional data, a likely benign interpretation. These results, consistent with family history data in the current study and available data from ClinVar, indicate that most, but not all, BRCA1 variants leading to an extended incorrect terminus constitute loss-of-function variants and underscore the need for comprehensive assessment of individual variants.

Evolution of the triplet BRCT domain

Organisms have evolved a complex system, called the DNA damage response (DDR), which maintains genome integrity. The DDR is responsible for identifying and repairing a variety of lesions and alterations in DNA. DDR proteins coordinate DNA damage detection, cell cycle arrest, and repair, with many of these events regulated by protein phosphorylation. In the human proteome, 23 proteins contain the BRCT (BRCA1 C-Terminus domain) domain, a modular signaling domain that can bind phosphopeptides and mediate protein-protein interactions. BRCTs can be found as functional single units, tandem (tBRCT), triplet (tpBRCT), and quartet. Here we examine the evolution of the tpBRCT architecture present in TOPBP1 (DNA topoisomerase II binding protein 1) and ECT2 (epithelial cell transforming 2), and their respective interaction partners RAD9 (Cell cycle checkpoint control protein RAD9) and CYK-4 (Rac GTPase-activating protein 1), with a focus on the conservation of the phosphopeptide-binding residues. The pair TOPBP1-RAD9 arose with the Eukaryotes and ECT2-CYK-4 with the Eumetazoans. Triplet structural and functional characteristics were conserved in almost all organisms. The first unit of the triplet (BRCT0) is different from the other two BRCTs but conserved between orthologs for both TOPBP1 and ECT2. BRCT domain evolution simulations suggest a trend to retain the singlet or towards two or three BRCT copies per protein consistent with functional tBRCT and tpBRCT architectures. Our results shed light on the emergence of the function and architecture of multiple BRCT domain organizations and provide information about the evolution of the BRCT triplet. Knowledge of BRCT domain evolution can improve the understanding of DNA damage response mechanisms and signal transduction in DDR.

Assessment of small in-frame indels and C-terminal nonsense variants of BRCA1 using a validated functional assay

BRCA1 (Breast Cancer 1, early onset) is linked to breast and ovarian cancer predisposition. Still, the risks conferred by a significant portion of BRCA1 variants identified in the population remains unknown. Most of these variants of uncertain significance are missense alterations. However, the functional implications of small in-frame deletions and/or insertions (indels) are also difficult to predict. Our group has previously evaluated the functional impact of 347 missense variants using an extensively validated transcriptional activity assay. Here we show a systematic assessment of 30 naturally occurring in-frame indels located at the C-terminal region of BRCA1. We identified positions sensitive and tolerant to alterations, expanding the knowledge of structural determinants of BRCA1 function. We further designed and assessed the impact of four single codon deletions in the tBRCT linker region and six nonsense variants at the C-terminus end of BRCA1. Amino acid substitutions, deletions or insertions in the disordered region do not significantly impact activity and are not likely to constitute pathogenic alleles. On the other hand, a sizeable fraction of in-frame indels at the BRCT domain significantly impact function. We then use a Bayesian integrative statistical model to derive the probability of pathogenicity for each variant. Our data highlights the importance of assessing the impact of small in-frame indels in BRCA1 to improve risk assessment and clinical decisions for carriers.

Identification of Two Genetic Loci Associated with Leukopenia after Chemotherapy in Patients with Breast Cancer

PURPOSE

To identify molecular predictors of grade 3/4 neutropenic or leukopenic events (NLE) after chemotherapy using a genome-wide association study (GWAS).

EXPERIMENTAL DESIGN

A GWAS was performed on patients in the phase III chemotherapy study SUCCESS-A (n = 3,322). Genotyping was done using the Illumina HumanOmniExpress-12v1 array. Findings were functionally validated with cell culture models and the genotypes and gene expression of possible causative genes were correlated with clinical treatment response and prognostic outcomes.

RESULTS

One locus on chromosome 16 (rs4784750; NLRC5; P = 1.56E-8) and another locus on chromosome 13 (rs16972207; TNFSF13B; P = 3.42E-8) were identified at a genome-wide significance level. Functional validation revealed that expression of these two genes is altered by genotype-dependent and chemotherapy-dependent activity of two transcription factors. Genotypes also showed an association with disease-free survival in patients with an NLE.

CONCLUSIONS

Two loci in NLRC5 and TNFSF13B are associated with NLEs. The involvement of the MHC I regulator NLRC5 implies the possible involvement of immuno-oncological pathways.

Functional Restoration of BRCA1 Nonsense Mutations by Aminoglycoside-Induced Readthrough

BRCA1 is a major tumor suppressor that functions in the accurate repair of DNA double-strand breaks via homologous recombination (HR). Nonsense mutations in BRCA1 lead to inactive truncated protein products and are associated with high risk of breast and ovarian cancer. These mutations generate premature termination codons (PTCs). Different studies have shown that aminoglycosides can induce PTC suppression by promoting stop codon readthrough and restoring full-length (FL) protein expression. The use of these compounds has been studied in clinical trials for genetic diseases such as cystic fibrosis and Duchenne muscular dystrophy, with encouraging results. Here we show proof-of-concept data demonstrating that the aminoglycoside G418 can induce BRCA1 PTC readthrough and restore FL protein synthesis and function. We first demonstrate that G418 treatment restores BRCA1 FL protein synthesis in HCC1395, a human breast tumor cell line carrying the R1751X mutation. HCC1395 cells treated with G418 also recover HR DNA repair and restore cell cycle checkpoint activation. A set of naturally occurring BRCA1 nonsense variants encoding different PTCs was evaluated in a GFP C-terminal BRCA1 construct model and BRCA1 PTC readthrough levels vary depending on the stop codon context. Because PTC readthrough could generate FL protein carrying pathogenic missense mutations, variants representing the most probable acquired amino acid substitutions in consequence of readthrough were functionally assessed by a validated transcription activation assay. Overall, this is the first study that evaluates the readthrough of PTC variants with clinical relevance in the breast and ovarian cancer-predisposing gene BRCA1.

An integrative model for the comprehensive classification of BRCA1 and BRCA2 variants of uncertain clinical significance

Loss-of-function variants in the BRCA1 and BRCA2 susceptibility genes predispose carriers to breast and/or ovarian cancer. The use of germline testing panels containing these genes has grown dramatically, but the interpretation of the results has been complicated by the identification of many sequence variants of undefined cancer relevance, termed "Variants of Uncertain Significance (VUS)." We have developed functional assays and a statistical model called VarCall for classifying BRCA1 and BRCA2 VUS. Here we describe a multifactorial extension of VarCall, called VarCall XT, that allows for co-analysis of multiple forms of genetic evidence. We evaluated the accuracy of models defined by the combinations of functional, in silico protein predictors, and family data for VUS classification. VarCall XT classified variants of known pathogenicity status with high sensitivity and specificity, with the functional assays contributing the greatest predictive power. This approach could be used to identify more patients that would benefit from personalized cancer risk assessment and management.

Two distinct mechanisms underlie estrogen-receptor-negative breast cancer susceptibility at the 2p23.2 locus

Genome wide-association studies (GWAS) have established over 400 breast cancer risk loci defined by common single nucleotide polymorphisms (SNPs), including several associated with estrogen-receptor (ER)-negative disease. Most of these loci have not been studied systematically and the mechanistic underpinnings of risk are largely unknown. Here we explored the landscape of genomic features at an ER-negative breast cancer susceptibility locus at chromosome 2p23.2 and assessed the functionality of 81 SNPs with strong evidence of association from previous fine mapping. Five candidate regulatory regions containing risk-associated SNPs were identified. Regulatory Region 1 in the first intron of WDR43 contains SNP rs4407214, which showed allele-specific interaction with the transcription factor USF1 in in vitro assays. CRISPR-mediated disruption of Regulatory Region 1 led to expression changes in the neighboring PLB1 gene, suggesting that the region acts as a distal enhancer. Regulatory Regions 2, 4, and 5 did not provide sufficient evidence for functionality in in silico and experimental analyses. Two SNPs (rs11680458 and rs1131880) in Regulatory Region 3, mapping to the seed region for miRNA-recognition sites in the 3' untranslated region of WDR43, showed allele-specific effects of ectopic expression of miR-376 on WDR43 expression levels. Taken together, our data suggest that risk of ER-negative breast cancer associated with the 2p23.2 locus is likely driven by a combinatorial effect on the regulation of WDR43 and PLB1.

PALB2 Variants: Protein Domains and Cancer Susceptibility

Since its discovery, partner and localizer of breast cancer 2 (BRCA2) (PALB2) has emerged as a major tumor suppressor gene linked to breast cancer (BC), pancreatic cancer (PC), and ovarian cancer (OC) susceptibility. Its protein product plays a pivotal role in the maintenance of genome integrity. Here we discuss the first functional evaluation of a large set of PALB2 missense variants of uncertain significance (VUSs). Assessment of 136 VUSs interrogating a range of PALB2 biological functions resulted in the identification of 15 variants with consistent loss of function across different assays. All loss-of-function variants are located at the PALB2 coiled coil (CC) or at the WD40 domain, highlighting the importance of modular domains mechanistically involved in the DNA damage response (DDR) and pinpointing their roles in tumor suppression.

Integration of functional assay data results provides strong evidence for classification of hundreds of BRCA1 variants of uncertain significance

Purpose

BRCA1 pathogenic variant heterozygotes are at a substantially increased risk for breast and ovarian cancer. The widespread uptake of testing has led to a significant increase in the detection of missense variants in BRCA1, the vast majority of which are variants of uncertain clinical significance (VUS), posing a challenge to genetic counseling. Here, we harness a wealth of functional data for thousands of variants to aid in variant classification.

Methods

We have collected, curated, and harmonized functional data for 2701 missense variants representing 24.5% of possible missense variants in BRCA1. Results were harmonized across studies by converting data into binary categorical variables (functional impact versus no functional impact). Using a panel of reference variants we identified a subset of assays with high sensitivity and specificity (≥80%) and apply the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) variant interpretation guidelines to assign evidence criteria for classification.

Results

Integration of data from validated assays provided ACMG/AMP evidence criteria in favor of pathogenicity for 297 variants or against pathogenicity for 2058 representing 96.2% of current VUS functionally assessed. We also explore discordant results and identify limitations in the approach.

Conclusion

High quality functional data are available for BRCA1 missense variants and provide evidence for classification of 2355 VUS according to their pathogenicity.

DNA damage response and repair in perspective: Aedes aegypti, Drosophila melanogaster and Homo sapiens

Background

The maintenance of genomic integrity is the responsibility of a complex network, denominated the DNA damage response (DDR), which controls the lesion detection and DNA repair. The main repair pathways are base excision repair (BER), nucleotide excision repair (NER), mismatch repair (MMR), homologous recombination repair (HR) and non-homologous end joining repair (NHEJ). They correct double-strand breaks (DSB), single-strand breaks, mismatches and others, or when the damage is quite extensive and repair insufficient, apoptosis is activated.

Methods

In this study we used the BLAST reciprocal best-hit methodology to search for DDR orthologs proteins in Aedes aegypti. We also provided a comparison between Ae. aegypti, D. melanogaster and human DDR network.

Results

Our analysis revealed the presence of ATR and ATM signaling, including the H2AX ortholog, in Ae. aegypti. Key DDR proteins (orthologs to RAD51, Ku and MRN complexes, XP-components, MutS and MutL) were also identified in this insect. Other proteins were not identified in both Ae. aegypti and D. melanogaster, including BRCA1 and its partners from BRCA1-A complex, TP53BP1, PALB2, POLk, CSA, CSB and POLβ. In humans, their absence affects DSB signaling, HR and sub-pathways of NER and BER. Seven orthologs not known in D. melanogaster were found in Ae. aegypti (RNF168, RIF1, WRN, RAD54B, RMI1, DNAPKcs, ARTEMIS).

Conclusions

The presence of key DDR proteins in Ae. aegypti suggests that the main DDR pathways are functional in this insect, and the identification of proteins not known in D. melanogaster can help fill gaps in the DDR network. The mapping of the DDR network in Ae. aegypti can support mosquito biology studies and inform genetic manipulation approaches applied to this vector.

CTDP1 regulates breast cancer survival and DNA repair through BRCT-specific interactions with FANCI

BRCA1 C-terminal domains are found in a specialized group of 23 proteins that function in the DNA damage response to protect genomic integrity. C-terminal domain phosphatase 1 (CTDP1) is the only phosphatase with a BRCA1 C-terminal domain in the human proteome, yet direct participation in the DNA damage response has not been reported. Examination of the CTDP1 BRCA1 C-terminal domain-specific protein interaction network revealed 103 high confidence interactions enriched in DNA damage response proteins, including FANCA and FANCI that are central to the Fanconi anemia DNA repair pathway necessary for the resolution of DNA interstrand crosslink damage. CTDP1 expression promotes DNA damage-induced FANCA and FANCD2 foci formation and enhances homologous recombination repair efficiency. CTDP1 was found to regulate multiple aspects of FANCI activity, including chromatin localization, interaction with γ-H2AX, and SQ motif phosphorylations. Knockdown of CTDP1 increases MCF-10A sensitivity to DNA interstrand crosslinks and double-strand breaks, but not ultraviolet radiation. In addition, CTDP1 knockdown impairs in vitro and in vivo growth of breast cancer cell lines. These results elucidate the molecular functions of CTDP1 in Fanconi anemia interstrand crosslink repair and identify this protein as a potential target for breast cancer therapy.

Impact of amino acid substitutions at secondary structures in the BRCT domains of the tumor suppressor BRCA1: Implications for clinical annotation

Genetic testing for BRCA1, a DNA repair protein, can identify carriers of pathogenic variants associated with a substantially increased risk for breast and ovarian cancers. However, an association with increased risk is unclear for a large fraction of BRCA1 variants present in the human population. Most of these variants of uncertain clinical significance lead to amino acid changes in the BRCA1 protein. Functional assays are valuable tools to assess the potential pathogenicity of these variants. Here, we systematically probed the effects of substitutions in the C terminus of BRCA1: the N- and C-terminal borders of its tandem BRCT domain, the BRCT-[N-C] linker region, and the α1 and α'1 helices in BRCT-[N] and -[C]. Using a validated transcriptional assay based on a fusion of the GAL4 DNA-binding domain to the BRCA1 C terminus (amino acids 1396-1863), we assessed the functional impact of 99 missense variants of BRCA1. We include the data obtained for these 99 missense variants in a joint analysis to generate the likelihood of pathogenicity for 347 missense variants in BRCA1 using VarCall, a Bayesian integrative statistical model. The results from this analysis increase our understanding of BRCA1 regions less tolerant to changes, identify functional borders of structural domains, and predict the likelihood of pathogenicity for 98% of all BRCA1 missense variants in this region recorded in the population. This knowledge will be critical for improving risk assessment and clinical treatment of carriers of BRCA1 variants.

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

Variants in genes encoding small GTPases and association with epithelial ovarian cancer susceptibility

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

Assessment of the Clinical Relevance of BRCA2 Missense Variants by Functional and Computational Approaches

Many variants of uncertain significance (VUS) have been identified in BRCA2 through clinical genetic testing. VUS pose a significant clinical challenge because the contribution of these variants to cancer risk has not been determined. We conducted a comprehensive assessment of VUS in the BRCA2 C-terminal DNA binding domain (DBD) by using a validated functional assay of BRCA2 homologous recombination (HR) DNA-repair activity and defined a classifier of variant pathogenicity. Among 139 variants evaluated, 54 had ?99% probability of pathogenicity, and 73 had ?95% probability of neutrality. Functional assay results were compared with predictions of variant pathogenicity from the Align-GVGD protein-sequence-based prediction algorithm, which has been used for variant classification. Relative to the HR assay, Align-GVGD significantly (p < 0.05) over-predicted pathogenic variants. We subsequently combined functional and Align-GVGD prediction results in a Bayesian hierarchical model (VarCall) to estimate the overall probability of pathogenicity for each VUS. In addition, to predict the effects of all other BRCA2 DBD variants and to prioritize variants for functional studies, we used the endoPhenotype-Optimized Sequence Ensemble (ePOSE) algorithm to train classifiers for BRCA2 variants by using data from the HR functional assay. Together, the results show that systematic functional assays in combination with in silico predictors of pathogenicity provide robust tools for clinical annotation of BRCA2 VUS.

Mutational heterogeneity in non-serous ovarian cancers

Epithelial ovarian cancer is a leading cause of death in gynecological cancers. While several systematic studies have revealed the mutation landscape of serous epithelial ovarian cancer, other non-serous subtypes of the disease have not been explored as extensively. Here we conduct exome sequencing of nine non-serous epithelial ovarian tumors (six endometrioid and three mucinous) and their corresponding normal DNA as well as a tumor-only granulosa cell sample. We integrated the exome data with targeted gene sequencing for 1,321 genes selected for their involvement in cancer from additional 28 non-serous ovarian tumors and compared our results to TCGA ovarian serous cystadenocarcinoma and uterine corpus endometrial carcinomas. Prevalence of TP53 mutations in non-serous was much lower than in serous epithelial OC, whereas the prevalence of PIK3CA, PIK3R1, PTEN, CTNNB1, ARID1A, and KRAS was higher. We confirmed the high prevalence of FOXL2 and KRAS mutations in granulosa cell tumors and in mucinous tumors, respectively. We also identified POLE proofreading domain mutations in three endometrioid ovarian tumors. These results highlight mutational differences between serous and non-serous ovarian cancers, and further distinguish different non-serous subtypes.

Mutational heterogeneity in non-serous ovarian cancers

Epithelial ovarian cancer is a leading cause of death in gynecological cancers. While several systematic studies have revealed the mutation landscape of serous epithelial ovarian cancer, other non-serous subtypes of the disease have not been explored as extensively. Here we conduct exome sequencing of nine non-serous epithelial ovarian tumors (six endometrioid and three mucinous) and their corresponding normal DNA as well as a tumor-only granulosa cell sample. We integrated the exome data with targeted gene sequencing for 1,321 genes selected for their involvement in cancer from additional 28 non-serous ovarian tumors and compared our results to TCGA ovarian serous cystadenocarcinoma and uterine corpus endometrial carcinomas. Prevalence of TP53 mutations in non-serous was much lower than in serous epithelial OC, whereas the prevalence of PIK3CA, PIK3R1, PTEN, CTNNB1, ARID1A, and KRAS was higher. We confirmed the high prevalence of FOXL2 and KRAS mutations in granulosa cell tumors and in mucinous tumors, respectively. We also identified POLE proofreading domain mutations in three endometrioid ovarian tumors. These results highlight mutational differences between serous and non-serous ovarian cancers, and further distinguish different non-serous subtypes.

BRCA2 Hypomorphic Missense Variants Confer Moderate Risks of Breast Cancer

Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case-control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A>C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G>A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G>A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G>T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789-99. ©2017 AACR.

A targeted genetic association study of epithelial ovarian cancer susceptibility

BACKGROUND

Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci.

RESULTS

At nine of eleven previously published EOC susceptibility regions (2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13), novel variants were identified that were more strongly associated with risk than previously reported variants. Beyond known susceptibility regions, no variants were found to be associated with EOC risk at genome-wide statistical significance (p <5x10(-8)), nor were any significant after Bonferroni correction for 17,000 variants (p< 3x10-6).

METHODS

A customized genotyping array was used to assess over 17,000 variants in coding, non-coding, regulatory, and known susceptibility regions in 4,973 EOC cases and 5,640 controls from 13 independent studies. Susceptibility for EOC overall and for select histotypes was evaluated using logistic regression adjusted for age, study site, and population substructure.

CONCLUSION

Given the novel variants identified within the 2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13 regions, larger follow-up genotyping studies, using imputation where necessary, are needed for fine-mapping and confirmation of low frequency variants that fall below statistical significance.

Germline missense pathogenic variants in the BRCA1 BRCT domain, p.Gly1706Glu and p.Ala1708Glu, increase cellular sensitivity to PARP inhibitor olaparib by a dominant negative effect

BRCA1-deficient cells show defects in DNA repair and rely on other members of the DNA repair machinery, which makes them sensitive to PARP inhibitors (PARPi). Although carrying a germline pathogenic variant in BRCA1/2 is the best determinant of response to PARPi, a significant percentage of the patients do not show sensitivity and/or display increased toxicity to the agent. Considering previously suggested mutation-specific BRCA1 haploinsufficiency, we aimed to investigate whether there are any differences in cellular response to PARPi olaparib depending on the BRCA1 mutation type. Lymphoblastoid cell lines derived from carriers of missense pathogenic variants in the BRCA1 BRCT domain (c.5117G > A, p.Gly1706Glu and c.5123C > A, p.Ala1708Glu) showed higher sensitivity to olaparib than cells with truncating variants or wild types (WT). Response to olaparib depended on a basal PARP enzymatic activity, but did not correlate with PARP1 expression. Interestingly, cellular sensitivity to the agent was associated with the level of BRCA1 recruitment into γH2AX foci, being the lowest in cells with missense variants. Since these variants lead to partially stable protein mutants, we propose a model in which the mutant protein acts in a dominant negative manner on the WT BRCA1, impairing the recruitment of BRCA1 into DNA damage sites and, consequently, increasing cellular sensitivity to PARPi. Taken together, our results indicate that carriers of different BRCA1 mutations could benefit from olaparib in a distinct way and show different toxicities to the agent, which could be especially relevant for a potential future use of PARPi as prophylactic agents in BRCA1 mutation carriers.

A multigene mutation classification of 468 colorectal cancers reveals a prognostic role for APC

Colorectal cancer (CRC) is a highly heterogeneous disease, for which prognosis has been relegated to clinicopathologic staging for decades. There is a need to stratify subpopulations of CRC on a molecular basis to better predict outcome and assign therapies. Here we report targeted exome-sequencing of 1,321 cancer-related genes on 468 tumour specimens, which identified a subset of 17 genes that best classify CRC, with APC playing a central role in predicting overall survival. APC may assume 0, 1 or 2 truncating mutations, each with a striking differential impact on survival. Tumours lacking any APC mutation carry a worse prognosis than single APC mutation tumours; however, two APC mutation tumours with mutant KRAS and TP53 confer the poorest survival among all the subgroups examined. Our study demonstrates a prognostic role for APC and suggests that sequencing of APC may have clinical utility in the routine staging and potential therapeutic assignment for CRC.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers

BACKGROUND

The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers.

METHODS

Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided.

RESULTS

The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed.

CONCLUSIONS

Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

Functional analysis of the 11q23.3 glioma susceptibility locus implicates PHLDB1 and DDX6 in glioma susceptibility

Glioma is the most common malignant primary brain tumor and is associated with poor prognosis. Genetic factors contributing to glioma risk have recently been investigated through genome-wide association studies (GWAS), implicating seven independent glioma risk loci in six chromosomal regions. Here, we performed an in-depth functional analysis of the risk locus proximal to the PHLDB1 gene on 11q23.3. We retrieved all SNPs in linkage disequilibrium (r² ≥ 0.2) with the glioma-associated SNP (rs498872) and performed a comprehensive bioinformatics and experimental functional analysis for the region. After testing candidate SNPs for allele-specific activity in a luciferase-based enhancer scanning assay, we established a subset of 10 functional SNPs in the promoters of PHLDB1 and DDX6, and in a putative enhancer element. Chromatin conformation capture (3C) identified a physical interaction between the enhancer element containing a functional SNP (rs73001406) and the promoter of the DDX6 gene. Knockdown experiments in cell culture and 3D assays to evaluate the role of PHLDB1 and DDX6 suggest that both genes may contribute to the phenotype. These studies reveal the functional landscape of the 11q23.3 glioma susceptibility locus and identify a network of functional SNPs in regulatory elements and two target genes as a possible mechanism driving glioma risk association.

Common variants at the CHEK2 gene locus and risk of epithelial ovarian cancer

Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

METHODS

We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS

Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION

We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

Genetic determinants of telomere length and risk of common cancers: a Mendelian randomization study

Epidemiological studies have reported inconsistent associations between telomere length (TL) and risk for various cancers. These inconsistencies are likely attributable, in part, to biases that arise due to post-diagnostic and post-treatment TL measurement. To avoid such biases, we used a Mendelian randomization approach and estimated associations between nine TL-associated SNPs and risk for five common cancer types (breast, lung, colorectal, ovarian and prostate cancer, including subtypes) using data on 51 725 cases and 62 035 controls. We then used an inverse-variance weighted average of the SNP-specific associations to estimate the association between a genetic score representing long TL and cancer risk. The long TL genetic score was significantly associated with increased risk of lung adenocarcinoma (P = 6.3 × 10(-15)), even after exclusion of a SNP residing in a known lung cancer susceptibility region (TERT-CLPTM1L) P = 6.6 × 10(-6)). Under Mendelian randomization assumptions, the association estimate [odds ratio (OR) = 2.78] is interpreted as the OR for lung adenocarcinoma corresponding to a 1000 bp increase in TL. The weighted TL SNP score was not associated with other cancer types or subtypes. Our finding that genetic determinants of long TL increase lung adenocarcinoma risk avoids issues with reverse causality and residual confounding that arise in observational studies of TL and disease risk. Under Mendelian randomization assumptions, our finding suggests that longer TL increases lung adenocarcinoma risk. However, caution regarding this causal interpretation is warranted in light of the potential issue of pleiotropy, and a more general interpretation is that SNPs influencing telomere biology are also implicated in lung adenocarcinoma risk.

Common Genetic Variation in Circadian Rhythm Genes and Risk of Epithelial Ovarian Cancer (EOC)

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10^-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

BRCA1 and BRCA2 genetic testing-pitfalls and recommendations for managing variants of uncertain clinical significance

BACKGROUND

Increasing use of BRCA1/2 testing for tailoring cancer treatment and extension of testing to tumour tissue for somatic mutation is moving BRCA1/2 mutation screening from a primarily prevention arena delivered by specialist genetic services into mainstream oncology practice. A considerable number of gene tests will identify rare variants where clinical significance cannot be inferred from sequence information alone. The proportion of variants of uncertain clinical significance (VUS) is likely to grow with lower thresholds for testing and laboratory providers with less experience of BRCA. Most VUS will not be associated with a high risk of cancer but a misinterpreted VUS has the potential to lead to mismanagement of both the patient and their relatives.

DESIGN

Members of the Clinical Working Group of ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) global consortium (www.enigmaconsortium.org) observed wide variation in practices in reporting, disclosure and clinical management of patients with a VUS. Examples from current clinical practice are presented and discussed to illustrate potential pitfalls, explore factors contributing to misinterpretation, and propose approaches to improving clarity.

RESULTS AND CONCLUSION

Clinicians, patients and their relatives would all benefit from an improved level of genetic literacy. Genetic laboratories working with clinical geneticists need to agree on a clinically clear and uniform format for reporting BRCA test results to non-geneticists. An international consortium of experts, collecting and integrating all available lines of evidence and classifying variants according to an internationally recognized system, will facilitate reclassification of variants for clinical use.

Cell-type-specific enrichment of risk-associated regulatory elements at ovarian cancer susceptibility loci

Understanding the regulatory landscape of the human genome is a central question in complex trait genetics. Most single-nucleotide polymorphisms (SNPs) associated with cancer risk lie in non-protein-coding regions, implicating regulatory DNA elements as functional targets of susceptibility variants. Here, we describe genome-wide annotation of regions of open chromatin and histone modification in fallopian tube and ovarian surface epithelial cells (FTSECs, OSECs), the debated cellular origins of high-grade serous ovarian cancers (HGSOCs) and in endometriosis epithelial cells (EECs), the likely precursor of clear cell ovarian carcinomas (CCOCs). The regulatory architecture of these cell types was compared with normal human mammary epithelial cells and LNCaP prostate cancer cells. We observed similar positional patterns of global enhancer signatures across the three different ovarian cancer precursor cell types, and evidence of tissue-specific regulatory signatures compared to non-gynecological cell types. We found significant enrichment for risk-associated SNPs intersecting regulatory biofeatures at 17 known HGSOC susceptibility loci in FTSECs (P = 3.8 × 10(-30)), OSECs (P = 2.4 × 10(-23)) and HMECs (P = 6.7 × 10(-15)) but not for EECs (P = 0.45) or LNCaP cells (P = 0.88). Hierarchical clustering of risk SNPs conditioned on the six different cell types indicates FTSECs and OSECs are highly related (96% of samples using multi-scale bootstrapping) suggesting both cell types may be precursors of HGSOC. These data represent the first description of regulatory catalogues of normal precursor cells for different ovarian cancer subtypes, and provide unique insights into the tissue specific regulatory variation with respect to the likely functional targets of germline genetic susceptibility variants for ovarian cancer.

Common Genetic Variation In Cellular Transport Genes and Epithelial Ovarian Cancer (EOC) Risk

BACKGROUND

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

METHODS

In total, DNA samples were obtained from 14,525 case subjects with invasive EOC and from 23,447 controls from 43 sites in the Ovarian Cancer Association Consortium (OCAC). Two hundred seventy nine SNPs, representing 131 genes, were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). SNP analyses were conducted using unconditional logistic regression under a log-additive model, and the FDR q<0.2 was applied to adjust for multiple comparisons.

RESULTS

The most significant evidence of an association for all invasive cancers combined and for the serous subtype was observed for SNP rs17216603 in the iron transporter gene HEPH (invasive: OR = 0.85, P = 0.00026; serous: OR = 0.81, P = 0.00020); this SNP was also associated with the borderline/low malignant potential (LMP) tumors (P = 0.021). Other genes significantly associated with EOC histological subtypes (p<0.05) included the UGT1A (endometrioid), SLC25A45 (mucinous), SLC39A11 (low malignant potential), and SERPINA7 (clear cell carcinoma). In addition, 1785 SNPs in six genes (HEPH, MGST1, SERPINA, SLC25A45, SLC39A11 and UGT1A) were imputed from the 1000 Genomes Project and examined for association with INV EOC in white-European subjects. The most significant imputed SNP was rs117729793 in SLC39A11 (per allele, OR = 2.55, 95% CI = 1.5-4.35, p = 5.66x10-4).

CONCLUSION

These results, generated on a large cohort of women, revealed associations between inherited cellular transport gene variants and risk of EOC histologic subtypes.

Incorporating computational resources in a cancer research program

Recent technological advances have transformed cancer genetics research. These advances have served as the basis for the generation of a number of richly annotated datasets relevant to the cancer geneticist. In addition, many of these technologies are now within reach of smaller laboratories to answer specific biological questions. Thus, one of the most pressing issues facing an experimental cancer biology research program in genetics is incorporating data from multiple sources to annotate, visualize, and analyze the system under study. Fortunately, there are several computational resources to aid in this process. However, a significant effort is required to adapt a molecular biology-based research program to take advantage of these datasets. Here, we discuss the lessons learned in our laboratory and share several recommendations to make this transition effective. This article is not meant to be a comprehensive evaluation of all the available resources, but rather highlight those that we have incorporated into our laboratory and how to choose the most appropriate ones for your research program.

Understanding Germ-Line Mutations inBRCA1

Germ-line mutations in BRCA1 account for the majority of familial breast and ovarian cancer cases and development of cancer in individuals who carry such mutations requires somatic inactivation of the normal allele. BRCA1 is highly polymorphic with more than 1,200 distinct documented variants. Approximately 70% of reported variants lead to absence of full-length BRCA1 protein, through loss of expression or protein truncation, and are suspected to predispose to cancer. These include regulatory mutations, splice site alterations, large rearrangements, large and small deletions or insertions, and nonsense mutations. However, characterizing the remaining missense alterations as either deleterious (cancer-associated mutations) or neutral variants is more complex, as the functional significance of the respective amino acid substitution is not straightforward to evaluate. In addition, many missense variants have been identified only once in defined ethnic groups and represent alleles with very low frequency. Most often, little information is available about segregation of the variant with disease in families, and assessment of disease risk for low frequency alleles through association studies is problematic, requiring a large number of samples stratified and matched by ethnicity. The fact that a significant proportion of BRCA1 variants remain unclassified represents a gap in risk assessment, such that individuals undergoing genetic testing will receive noninformative test results. An approach for assessing the potential clinical significance of missense variants is to combine available genetic data with functional and structural studies. Here we review the available information on BRCA1 variants and explore ways in which we can analyze unclassified variants.

DNA damage response: Determining the fate of phosphorylated histone H2AX

Phosphorylation of histone H2AX at Serine 139 is one of the earliest events after DNA damage and is required for the retention of factors involved in repair at the site of the break. Intriguingly, H2AX phosphorylation spreads from the vicinity of the break to both directions spanning large chromosomal regions. Phosphorylated H2AX (also known as gamma-H2AX) then progressively disappears with kinetics that correlates with the completion of DNA repair. Despite intense investigation on the kinases and stimuli involved in gamma-H2AX formation, the mechanism of gamma-H2AX disappearance has remained obscure. Three recent papers shed light on this process and suggest that H2AX may serve as a signaling platform that integrates repair and cell cycle checkpoints.

The accidental cancer geneticist: Hilário de Gouvêa and hereditary retinoblastoma

In 1872 a Brazilian ophthalmologist performed an enucleation in a young boy with retinoblastoma. This boy survived and married a woman without any family history of cancer. The couple had two daughters with bilateral retinoblastoma also seen by the same ophthalmologist, Hilário de Gouvêa. This case became the first documented report of a family with retinoblastoma in more than one generation. Here we examine the life of de Gouvêa and his contribution which raised the possibility that cancer had a genetic basis. We discuss how de Gouvêa's mind had been prepared to realize the importance of this observation. We attempt to define the conditions that allowed not only his discovery, but also the report of the findings and a dogged pursuit for credit over many years in a country which had virtually no research tradition and was still grappling with its recent colonial history.

Functional annotation signatures of disease susceptibility loci improve SNP association analysis

BACKGROUND

Genetic association studies are conducted to discover genetic loci that contribute to an inherited trait, identify the variants behind these associations and ascertain their functional role in determining the phenotype. To date, functional annotations of the genetic variants have rarely played more than an indirect role in assessing evidence for association. Here, we demonstrate how these data can be systematically integrated into an association study's analysis plan.

RESULTS

We developed a Bayesian statistical model for the prior probability of phenotype-genotype association that incorporates data from past association studies and publicly available functional annotation data regarding the susceptibility variants under study. The model takes the form of a binary regression of association status on a set of annotation variables whose coefficients were estimated through an analysis of associated SNPs in the GWAS Catalog (GC). The functional predictors examined included measures that have been demonstrated to correlate with the association status of SNPs in the GC and some whose utility in this regard is speculative: summaries of the UCSC Human Genome Browser ENCODE super-track data, dbSNP function class, sequence conservation summaries, proximity to genomic variants in the Database of Genomic Variants and known regulatory elements in the Open Regulatory Annotation database, PolyPhen-2 probabilities and RegulomeDB categories. Because we expected that only a fraction of the annotations would contribute to predicting association, we employed a penalized likelihood method to reduce the impact of non-informative predictors and evaluated the model's ability to predict GC SNPs not used to construct the model. We show that the functional data alone are predictive of a SNP's presence in the GC. Further, using data from a genome-wide study of ovarian cancer, we demonstrate that their use as prior data when testing for association is practical at the genome-wide scale and improves power to detect associations.

CONCLUSIONS

We show how diverse functional annotations can be efficiently combined to create 'functional signatures' that predict the a priori odds of a variant's association to a trait and how these signatures can be integrated into a standard genome-wide-scale association analysis, resulting in improved power to detect truly associated variants.

Probing structure-function relationships in missense variants in the carboxy-terminal region of BRCA1

Germline inactivating variants in BRCA1 lead to a significantly increased risk of breast and ovarian cancers in carriers. While the functional effect of many variants can be inferred from the DNA sequence, determining the effect of missense variants present a significant challenge. A series of biochemical and cell biological assays have been successfully used to explore the impact of these variants on the function of BRCA1, which contribute to assessing their likelihood of pathogenicity. It has been determined that variants that co-localize with structural or functional motifs are more likely to disrupt the stability and function of BRCA1. Here we assess the functional impact of 37 variants chosen to probe the functional impact of variants in phosphorylation sites and in the BRCT domains. In addition, we perform a meta-analysis of 170 unique variants tested by the transcription activation assays in the carboxy-terminal domain of BRCA1 using a recently developed computation model to provide assessment for functional impact and their likelihood of pathogenicity.

Functional assays for analysis of variants of uncertain significance in BRCA2

Missense variants in the BRCA2 gene are routinely detected during clinical screening for pathogenic mutations in patients with a family history of breast and ovarian cancer. These subtle changes frequently remain of unknown clinical significance because of the lack of genetic information that may help establish a direct correlation with cancer predisposition. Therefore, alternative ways of predicting the pathogenicity of these variants are urgently needed. Since BRCA2 is a protein involved in important cellular mechanisms such as DNA repair, replication, and cell cycle control, functional assays have been developed that exploit these cellular activities to explore the impact of the variants on protein function. In this review, we summarize assays developed and currently utilized for studying missense variants in BRCA2. We specifically depict details of each assay, including variants of uncertain significance analyzed, and describe a validation set of (genetically) proven pathogenic and neutral missense variants to serve as a golden standard for the validation of each assay. Guidelines are proposed to enable implementation of laboratory-based methods to assess the impact of the variant on cancer risk.

Lessons from postgenome-wide association studies: functional analysis of cancer predisposition loci

In the last few years, genome-wide association studies (GWASs) have identified hundreds of predisposition loci for several types of human cancers. Recent progress has been made in determining the underlying mechanisms through which different single-nucleotide polymorphisms (SNPs) affect predisposition to cancer. Although there has been much debate about the clinical utility of GWASs, less attention has been paid to how GWASs and post-GWASs functional analysis have contributed to understanding the aetiology of cancer. Most common variants associated with cancer risk are localized in nonprotein-coding regions highlighting transcriptional regulation as a common theme in the mechanism of cancer predisposition. Here, we outline strategies to functionally dissect predisposition loci and discuss their limitations as well as challenges for future studies.

GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer

Genome-wide association studies (GWAS) have identified four susceptibility loci for epithelial ovarian cancer (EOC), with another two suggestive loci reaching near genome-wide significance. We pooled data from a GWAS conducted in North America with another GWAS from the UK. We selected the top 24,551 SNPs for inclusion on the iCOGS custom genotyping array. We performed follow-up genotyping in 18,174 individuals with EOC (cases) and 26,134 controls from 43 studies from the Ovarian Cancer Association Consortium. We validated the two loci at 3q25 and 17q21 that were previously found to have associations close to genome-wide significance and identified three loci newly associated with risk: two loci associated with all EOC subtypes at 8q21 (rs11782652, P = 5.5 × 10(-9)) and 10p12 (rs1243180, P = 1.8 × 10(-8)) and another locus specific to the serous subtype at 17q12 (rs757210, P = 8.1 × 10(-10)). An integrated molecular analysis of genes and regulatory regions at these loci provided evidence for functional mechanisms underlying susceptibility and implicated CHMP4C in the pathogenesis of ovarian cancer.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Modes of delivery of genetic testing services and the uptake of cancer risk management strategies in BRCA1 and BRCA2 carriers

BRCA testing services are now offered by various healthcare providers, thus it is important to evaluate whether the implementation of cancer risk management (CRM) strategies varies by service provider. Using a registry-based sample of 795 female BRCA mutation carriers, we explored the association between uptake of CRM strategies with duration of genetic counseling (GC) sessions, provider type, and other demographic and clinical variables. All participants completed a baseline questionnaire. Information about uptake of CRM strategies was collected for a subset of 438 participants who completed additional questions. Summary statistics and Pearson chi-squared analysis were used to examine the associations between demographic and clinical variables with service delivery factors and with the uptake of various CRM strategies. Overall uptake of CRM strategies was high across all provider types. However, GC sessions were longer when provided by a genetics professional than by another provider (p < 0.001). Furthermore, higher frequencies of uptake of most CRM strategies were associated with longer GC sessions and when testing was performed by a genetics professional. Identification of factors to optimize delivery of these specialized GC services is important to maximize implementation of CRM strategies in BRCA carriers.

Biallelic deleterious BRCA1 mutations in a woman with early-onset ovarian cancer

BRCA1 and BRCA2 are the most important breast and ovarian cancer susceptibility genes. Biallelic mutations in BRCA2 can lead to Fanconi anemia and predisposition to cancers, whereas biallelic BRCA1 mutations have not been confirmed, presumably because one wild-type BRCA1 allele is required during embryogenesis. This study describes an individual who was diagnosed with ovarian carcinoma at age 28 and found to have one allele with a deleterious mutation in BRCA1, c.2457delC (p.Asp821Ilefs*25), and a second allele with a variant of unknown significance in BRCA1, c.5207T>C (p.Val1736Ala). Medical records revealed short stature, microcephaly, developmental delay, and significant toxicity from chemotherapy. BRCA1 p.Val1736Ala cosegregated with cancer in multiple families, associated tumors showed loss of wild-type BRCA1, and BRCA1 p.Val1736Ala showed reduced DNA damage localization. These findings represent the first validated example of biallelic deleterious human BRCA1 mutations and have implications for the interpretation of genetic test results.

SIGNIFICANCE

Accurate assessment of genetic testing data for BRCA1 mutations is essential for clinical monitoring and treatment strategies. Here, we report the fi rst validated example of an individual with biallelic BRCA1 mutations, early-onset ovarian cancer, and clinically significant hypersensitivity to chemotherapy.

Charting the landscape of tandem BRCT domain-mediated protein interactions

Eukaryotic cells have evolved an intricate system to resolve DNA damage to prevent its transmission to daughter cells. This system, collectively known as the DNA damage response (DDR) network, includes many proteins that detect DNA damage, promote repair, and coordinate progression through the cell cycle. Because defects in this network can lead to cancer, this network constitutes a barrier against tumorigenesis. The modular BRCA1 carboxyl-terminal (BRCT) domain is frequently present in proteins involved in the DDR, can exist either as an individual domain or as tandem domains (tBRCT), and can bind phosphorylated peptides. We performed a systematic analysis of protein-protein interactions involving tBRCT in the DDR by combining literature curation, yeast two-hybrid screens, and tandem affinity purification coupled to mass spectrometry. We identified 23 proteins containing conserved BRCT domains and generated a human protein-protein interaction network for seven proteins with tBRCT. This study also revealed previously unknown components in DNA damage signaling, such as COMMD1 and the target of rapamycin complex mTORC2. Additionally, integration of tBRCT domain interactions with DDR phosphoprotein studies and analysis of kinase-substrate interactions revealed signaling subnetworks that may aid in understanding the involvement of tBRCT in disease and DNA repair.