Investigation of monogenic causes of familial breast cancer: data from the BEACCON case-control study

Association between missense variants of uncertain significance in the CHEK2 gene and hereditary breast cancer: a cosegregation and bioinformatics analysis

Inherited mutations in the CHEK2 gene have been associated with an increased lifetime risk of developing breast cancer (BC). We aim to identify in the study population the prevalence of mutations in the CHEK2 gene in diagnosed BC patients, evaluate the phenotypic characteristics of the tumor and family history, and predict the deleteriousness of the variants of uncertain significance (VUS). A genetic study was performed, from May 2016 to April 2020, in 396 patients diagnosed with BC at the University Hospital Lozano Blesa of Zaragoza, Spain. Patients with a genetic variant in the CHEK2 gene were selected for the study. We performed a descriptive analysis of the clinical variables, a bibliographic review of the variants, and a cosegregation study when possible. Moreover, an in-depth bioinformatics analysis of CHEK2 VUS was carried out. We identified nine genetic variants in the CHEK2 gene in 10 patients (two pathogenic variants and seven VUS). This supposes a prevalence of 0.75% and 1.77%, respectively. In all cases, there was a family history of BC in first- and/or second-degree relatives. We carried out a cosegregation study in two families, being positive in one of them. The bioinformatics analyses predicted the pathogenicity of six of the VUS. In conclusion, CHEK2 mutations have been associated with an increased risk for BC. This risk is well-established for foundation variants. However, the risk assessment for other variants is unclear. The incorporation of bioinformatics analysis provided supporting evidence of the pathogenicity of VUS.

Integrating a Polygenic Risk Score into a clinical setting would impact risk predictions in familial breast cancer

BACKGROUND

Low-impact genetic variants identified in population-based genetic studies are not routinely measured as part of clinical genetic testing in familial breast cancer (BC). We studied the consequences of integrating an established Polygenic Risk Score (PRS) (BCAC 313, PRS313) into clinical sequencing of women with familial BC in Sweden.

METHODS

We developed an add-on sequencing panel to capture 313 risk variants in addition to the clinical screening of hereditary BC genes. Index patients with no pathogenic variant from 87 families, and 1000 population controls, were included in comparative PRS calculations. Including detailed family history, sequencing results and tumour pathology information, we used BOADICEA (Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm) V.6 to estimate contralateral and lifetime risks without and with PRS313.

RESULTS

Women with BC but no pathogenic variants in hereditary BC genes have a higher PRS313 compared with population controls (mean+0.78 SD, p<3e-9). Implementing PRS313 in the clinical risk estimation before their BC diagnosis would have changed the recommended follow-up in 24%-45% of women.

CONCLUSIONS

Our results show the potential impact of incorporating PRS313 directly in the clinical genomic investigation of women with familial BC.

Detection Rate and Spectrum of Pathogenic Variations in a Cohort of 83 Patients with Suspected Hereditary Risk of Kidney Cancer

Hereditary predisposition to cancer affects about 3-5% of renal cancers. Testing criteria have been proposed in France for genetic testing of non-syndromic renal cancer. Our study explores the detection rates associated with our testing criteria. Using a comprehensive gene panel including 8 genes related to renal cancer and 50 genes related to hereditary predisposition to other cancers, we evaluated the detection rate of pathogenic variants in a cohort of 83 patients with suspected renal cancer predisposition. The detection rate was 7.2% for the renal cancer genes, which was 2.41-fold higher than the estimated 3% proportion of unselected kidney cases with inherited risk. Pathogenic variants in renal cancer genes were observed in 44.5% of syndromic cases, and in 2.7% of non-syndromic cases. Incidental findings were observed in CHEK2, MSH2, MUTYH and WRN. CHEK2 was associated with renal cancer (OR at 7.14; 95% CI 1.74-29.6; p < 0.003) in our study in comparison to the gnomAD control population. The detection rate in renal cancer genes was low in non-syndromic cases. Additional causal mechanisms are probably involved, and further research is required to find them. A study of the management of renal cancer risk for CHEK2 pathogenic variant carriers is needed.

The CHK2 kinase is recurrently mutated and functionally impaired in the germline of pediatric cancer patients

Predisposing CHEK2 germline variants are associated with various adult-type malignancies, whereas their impact on cancer susceptibility in childhood cancer is unclear. To understand the frequency of germline variants in the CHEK2 gene and their impact on pediatric malignancies, we used whole-exome sequencing to search for CHEK2 variants in the germlines of 418 children diagnosed with cancer in our clinics. Moreover, we performed functional analysis of the pathogenic CHEK2 variants to analyze the effect of the alterations on CHK2 protein function. We detected a CHEK2 germline variant in 32/418 (7.7%) pediatric cancer patients and 46.8% of them had leukemia. Functional analysis of the pathogenic variants revealed that 5 pathogenic variants impaired CHK2 protein function. 6/32 patients carried one of these clearly damaging CHEK2 variants and two of them harbored a matching family history of cancer. In conclusion, we detected germline CHEK2 variants in 7.7% of all pediatric cancer patients, of which a minority but still relevant fraction of approximately 20% had a profound impact on protein expression or its phosphorylation after irradiation-induced DNA damage. Accordingly, we conclude that CHEK2 variants increase the risk for not only adult-onset but also pediatric cancer.

Contribution of large genomic rearrangements in PALB2 to familial breast cancer: implications for genetic testing

BACKGROUND

PALB2 is the most important contributor to familial breast cancer after BRCA1 and BRCA2. Large genomic rearrangements (LGRs) in BRCA1 and BRCA2 are routinely assessed in clinical testing and are a significant contributor to the yield of actionable findings. In contrast, the contribution of LGRs in PALB2 has not been systematically studied.

METHODS

We performed targeted sequencing and real-time qPCR validation to identify LGRs in PALB2 in 5770 unrelated patients with familial breast cancer and 5741 cancer-free control women from the same Australian population.

RESULTS

Seven large deletions ranging in size from 0.96 kbp to 18.07 kbp involving PALB2 were identified in seven cases, while no LGRs were identified in any of the controls. Six LGRs were considered pathogenic as they included one or more exons of PALB2 and disrupted the WD40 domain at the C terminal end of the PALB2 protein while one LGR only involved a partial region of intron 10 and was considered a variant of unknown significance. Altogether, pathogenic LGRs identified in this study accounted for 10.3% (6 of 58) of the pathogenic PALB2 variants detected among the 5770 families with familial breast cancer.

CONCLUSIONS

Our data show that a clinically important proportion of PALB2 pathogenic mutations in Australian patients with familial breast cancer are LGRs. Such observations have provided strong support for inclusion of PALB2 LGRs in routine clinical genetic testing.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Genetic analyses of DNA repair pathway associated genes implicate new candidate cancer predisposing genes in ancestrally defined ovarian cancer cases

Not all familial ovarian cancer (OC) cases are explained by pathogenic germline variants in known risk genes. A candidate gene approach involving DNA repair pathway genes was applied to identify rare recurring pathogenic variants in familial OC cases not associated with known OC risk genes from a population exhibiting genetic drift. Whole exome sequencing (WES) data of 15 OC cases from 13 families tested negative for pathogenic variants in known OC risk genes were investigated for candidate variants in 468 DNA repair pathway genes. Filtering and prioritization criteria were applied to WES data to select top candidates for further analyses. Candidates were genotyped in ancestry defined study groups of 214 familial and 998 sporadic OC or breast cancer (BC) cases and 1025 population-matched controls and screened for additional carriers in 605 population-matched OC cases. The candidate genes were also analyzed in WES data from 937 familial or sporadic OC cases of diverse ancestries. Top candidate variants in ERCC5, EXO1, FANCC, NEIL1 and NTHL1 were identified in 5/13 (39%) OC families. Collectively, candidate variants were identified in 7/435 (1.6%) sporadic OC cases and 1/566 (0.2%) sporadic BC cases versus 1/1025 (0.1%) controls. Additional carriers were identified in 6/605 (0.9%) OC cases. Tumour DNA from ERCC5, NEIL1 and NTHL1 variant carriers exhibited loss of the wild-type allele. Carriers of various candidate variants in these genes were identified in 31/937 (3.3%) OC cases of diverse ancestries versus 0-0.004% in cancer-free controls. The strategy of applying a candidate gene approach in a population exhibiting genetic drift identified new candidate OC predisposition variants in DNA repair pathway genes.

Somatic inactivation of breast cancer predisposition genes in tumors associated with pathogenic germline variants

BACKGROUND

Breast cancers (BCs) that arise in individuals heterozygous for a germline pathogenic variant in a susceptibility gene, such as BRCA1 and BRCA2, PALB2, and RAD51C, have been shown to exhibit biallelic loss in the respective genes and be associated with triple-negative breast cancer (TNBC) and distinctive somatic mutational signatures. Tumor sequencing thus presents an orthogonal approach to assess the role of candidate genes in BC development.

METHODS

Exome sequencing was performed on paired normal-breast tumor DNA from 124 carriers of germline loss-of-function (LoF) or missense variant carriers in 15 known and candidate BC predisposition genes identified in the BEACCON case-control study. Biallelic inactivation and association with tumor genome features including mutational signatures and homologous recombination deficiency (HRD) score were investigated.

RESULTS

BARD1-carrying TNBC (4 of 5) displayed biallelic loss and associated high HRD scores and mutational signature 3, as did a RAD51D-carrying TNBC and ovarian cancer. Biallelic loss was less frequent in BRIP1 BCs (4 of 13) and had low HRD scores. In contrast to other established BC genes, BCs from carriers of CHEK2 LoF (6 of 17) or missense (2 of 20) variant had low rates of biallelic loss. Exploratory analysis of BC from carriers of LoF variants in candidate genes such as BLM, FANCM, PARP2, and RAD50 found little evidence of biallelic inactivation.

CONCLUSIONS

BARD1 and RAD51D behave as classic BRCA-like predisposition genes with biallelic inactivation, but this was not observed for any of the candidate genes. However, as demonstrated for CHEK2, the absence of biallelic inactivation does not provide definitive evidence against the gene's involvement in BC predisposition.

Estimating the proportion of pathogenic variants from breast cancer case-control data: Application to calibration of ACMG/AMP variant classification criteria

For genes with reliable estimates of disease risk associated with loss-of-function variants, case-control data can be used to estimate the proportion of variants of typical risk effect for defined groups of variants, of relevance for variant classification. A calculation was derived for a maximum likelihood estimate of the proportion of pathogenic variants of typical effect from case-control data and applied to rare variant counts for ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, and RAD51D from published breast cancer studies: BEACCON (5770 familial cases and 5741 controls) and breast cancer risk after diagnostic sequencing (60,466 familial and population-based cases and 53,461 controls). There was significant evidence of pathogenic variants among rare noncoding variants, in particular deeper intronic variants, for BRCA1 (13%, p = 8.3 × 10^-7 ), BRCA2 (6%, p = 0.016) and PALB2 (13%, p = 0.001). The estimated proportion of pathogenic missense variants varied markedly between genes, generally with enrichment in familial cases, for example, 9% for BRCA2 versus 60%-90% for CHEK2. Stratifying missense variants by position indicated that, for most genes, location within a functional domain significantly predicted pathogenicity, whereas location outside domains provided robust evidence against pathogenicity. Our approach provides novel insights into the spectrum of pathogenic variants of specific breast cancer genes and has wider application to inform gene-focused specifications of American College of Medical Genetics and Genomics (ACMG)/Association of Molecular Pathology (AMP) codes for variant curation.

Integration of tumour sequencing and case-control data to assess pathogenicity of RAD51C missense variants in familial breast cancer

While protein-truncating variants in RAD51C have been shown to predispose to triple-negative (TN) breast cancer (BC) and ovarian cancer, little is known about the pathogenicity of missense (MS) variants. The frequency of rare RAD51C MS variants was assessed in the BEACCON study of 5734 familial BC cases and 14,382 population controls, and findings were integrated with tumour sequencing data from 21 cases carrying a candidate variant. Collectively, a significant enrichment of rare MS variants was detected in cases (MAF < 0.001, OR 1.57, 95% CI 1.00-2.44, p = 0.05), particularly for variants with a REVEL score >0.5 (OR 3.95, 95% CI 1.40-12.01, p = 0.006). Sequencing of 21 tumours from 20 heterozygous and 1 homozygous carriers of nine candidate MS variants identified four cases with biallelic inactivation through loss of the wild-type allele, while six lost the variant allele and ten that remained heterozygous. Biallelic loss of the wild-type alleles corresponded strongly with ER- and TN breast tumours, high homologous recombination deficiency scores and mutational signature 3. Using this approach, the p.Gly264Ser variant, which was previously suspected to be pathogenic based on small case-control analyses and loss of activity in in vitro functional assays, was shown to be benign with similar prevalence in cases and controls and seven out of eight tumours showing no biallelic inactivation or characteristic mutational signature. Conversely, evaluation of case-control findings and tumour sequencing data identified p.Ile144Thr, p.Arg212His, p.Gln143Arg and p.Gly114Arg as variants warranting further investigation.

Association of AXIN2 s2240308 C>T, rs1133683 C>T, rs7224837 A>G Polymorphisms with Susceptibility to Breast Cancer

BACKGROUND

The association of genetic polymorphisms at Axis inhibition protein 2 (AXIN2) gene and susceptibility to different cancer has attracted much interest. The present study aimed to evaluate the association between AXIN2 rs2240308 C>T, rs1133683 C>T, rs7224837 A>G polymorphisms with susceptibility to breast cancer.

METHODS

A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay was designed to genotype the AXIN2 rs2240308 C>T, rs1133683 C>T, rs7224837 A>G polymorphisms among 150 breast cancer patients and 150 healthy subjects.

RESULTS

The frequencies of these genetic variants were in agreement with Hardy-Weinberg equilibrium in healthy controls (p>0.05). The frequencies of AXIN2 rs2240308 C>T, rs1133683 C>T, rs7224837 A>G genotypes were similar in breast cancer patients and controls. There was no a significant association between the AXIN2 SNP and risk of breast cancer.

CONCLUSION

The impact of AXIN2 polymorphisms in the breast cancer development remains unclear. Our results indicated that AXIN2 rs2240308, rs7224837 and rs1133683 polymorphisms did not contribute to increased risk of breast cancer. More studies with larger sample sizes and diverse ethnicities are warranted to verify our finding.

Investigation of monogenic causes of familial breast cancer: data from the BEACCON case-control study

Breast cancer (BC) has a significant heritable component but the genetic contribution remains unresolved in the majority of high-risk BC families. This study aims to investigate the monogenic causes underlying the familial aggregation of BC beyond BRCA1 and BRCA2, including the identification of new predisposing genes. A total of 11,511 non-BRCA familial BC cases and population-matched cancer-free female controls in the BEACCON study were investigated in two sequencing phases: 1303 candidate genes in up to 3892 cases and controls, followed by validation of 145 shortlisted genes in an additional 7619 subjects. The coding regions and exon–intron boundaries of all candidate genes and 14 previously proposed BC genes were sequenced using custom designed sequencing panels. Pedigree and pathology data were analysed to identify genotype-specific associations. The contribution of ATM, PALB2 and CHEK2 to BC predisposition was confirmed, but not RAD50 and NBN. An overall excess of loss-of-function (LoF) (OR 1.27, p = 9.05 × 10⁻⁹) and missense (OR 1.27, p = 3.96 × 10⁻⁷³) variants was observed in the cases for the 145 candidate genes. Leading candidates harbored LoF variants with observed ORs of 2–4 and individually accounted for no more than 0.79% of the cases. New genes proposed by this study include NTHL1, WRN, PARP2, CTH and CDK9. The new candidate BC predisposition genes identified in BEACCON indicate that much of the remaining genetic causes of high-risk BC families are due to genes in which pathogenic variants are both very rare and convey only low to moderate risk.