1. Introduction
Women who carry a pathogenic mutation in BRCA1 or BRCA2 are at increased risk of developing breast and ovarian cancer. Kuchenbaecker et al. estimated the cumulative breast cancer risk to age 80 years to be 72% (95% confidence interval (CI) = 65–79%) and 69% (95% CI = 61–77%) for BRCA1 and BRCA2 pathogenic variant carriers, respectively (Kuchenbaecker et al., Reference Kuchenbaecker, Hopper and Barnes2017). For ovarian cancer, the cumulative cancer risk to age 80 years is estimated to be 44% (95% CI = 36–53%) and 17% (95% CI = 11–25%) for BRCA1 and BRCA2 pathogenic variant carriers, respectively (Kuchenbaecker et al., Reference Kuchenbaecker, Hopper and Barnes2017).
In many countries, genetic testing for BRCA1 and BRCA2 has shown clear clinical utility and validity. Evidence-based best practice guidelines are available to inform the clinical management of women who carry a BRCA1 or BRCA2 pathogenic variant. These guidelines support personalized risk assessment, targeted treatment regimens and informed decision-making about the use of risk-reducing medications, bilateral salpingo-oophorectomy, mammography, risk-reducing mastectomy, magnetic resonance imaging and other screening modalities.
It is well established that, as a consequence of a founder effect, different ethnic and geographical regions can have different BRCA1 and BRCA2 mutation spectra and prevalence rates.
Cybulski et al. recently reported on the mutation spectrum in BRCA1, BRCA2 and other genes associated, or putatively associated, with increased risk of breast cancer in 1018 probands from multiple-case breast cancer families from Poland. In their study, three founder mutations were identified with high prevalence: BRCA1:c.5266dup (20%, 204/1018), BRCA1:c.181T>G (8.3%, 84/1018) and BRCA1:c.4035del (1.5%, 15/1018). Other mutations reported at lower prevalence (≤1.0%) in the 1108 familial breast cancer cases included BRCA1:c.3700_3704del (1.0% 10/1018), BRCA1:c.68_69del (0.9%, 9/1018), BRCA1:c.5251C>T (0.6%, 6/1018) and BRCA1:c.5346G>A (0.5%, 5/1108) (Cybulski et al., Reference Cybulski, Kluzniak and Huzarski2019). Recurrent mutations were also reported in BRCA2 by Cybulski et al., all at a prevalence of 0.5% or below in the Polish women with familial breast cancer (Cybulski et al., Reference Cybulski, Kluzniak and Huzarski2019).
The BRCA1 mutations c.5266dupC and c.4035delA have been confirmed by haplotype analysis as founder mutations in Eastern Europe (Hamel et al., Reference Hamel, Feng and Foretova2011; Janavicius et al., Reference Janavicius, Rudaitis and Feng2013). BRCA1:c.181T>G mutation carriers of Polish and Jewish ancestry have been shown to present the same haplotype (Kaufman et al., Reference Kaufman, Laitman, Gronwald, Lubinski and Friedman2009). Haplotype analysis of BRCA1:c.68_69delAG showed a common haplotype among Ashkenazi Jews (Laitman et al., Reference Laitman, Feng and Zamir2013). Currently, genetic testing in the Polish population mainly relies on testing of the founder mutations. Comprehensive BRCA1 and BRCA2 genetic testing could, however, identify more women with pathogenic variants, thus leading to improved cancer prevention for more women at high risk of breast and ovarian cancer.
Participants in this study were unrelated women affected with breast or ovarian cancer from south-west Poland and west Ukraine. The Polish participants had previously been genotyped for BRCA1:c.5266dup, BRCA1:c.181T>G, BRCA1:c.4035del, BRCA1:c.68_69del and BRCA2:c.5946del. Thirty-four women identified as carrying one of these mutations were excluded from further testing in this study. Participants recruited in Ukraine had not undergone prior genetic testing and have thus all been included in the genetic testing reported in this study.
2. Materials and methods
2.1. Study participants
The women participating in this study were unrelated women diagnosed with breast and/or ovarian cancer recruited after or during oncological treatment from Wroclaw Medical University, Lower Silesia, Poland, between 2004 and 2008, or Lviv State Oncology Regional Treatment and Diagnostic Center, Lviv, Ukraine, between 2004 and 2010, as described previously (Myszka et al., Reference Myszka, Nguyen-Dumont and Karpinski2018). The Polish cohort consisted of 238 women affected with breast cancer, 95 women affected with ovarian cancer and 4 women affected with breast and ovarian cancer. Of the 242 women with breast cancer, 95 had hereditary breast cancer, 18 had familial breast cancer and 125 were sporadic cases, according to the criteria described by Berliner et al. (Reference Berliner and Fay2007). Of the 95 Polish women with ovarian cancer, 28 had hereditary ovarian cancer, 10 had familial ovarian cancer and 57 were sporadic ovarian cancer cases. All four women with breast and ovarian cancers met the criteria for hereditary disease. The Polish cohort thus consisted of 337 women, all of whom had previously been genotyped for four mutations in BRCA1 (c.5266dup, c.181T>G, c.4035del, c.68_69del) (Table 1) and one mutation in BRCA2 (c.5946delT). Thirty-four women were identified as carriers of one of these BRCA1 founder mutations and were thus not included in the targeted sequence screening described in this study. No BRCA2 c.5946delT carrier was observed.
a Variant nomenclature according to the Human Genome Variation Society (HGVS). Transcript sequences is BRCA1: NM_007294.3.
b Variant nomenclature according to the Breast Cancer Information Core (BIC) (https://research.nhgri.nih.gov/bic).
c Minor allele frequency (MAF) in the ExAC database, in the non-Finnish European population minus The Cancer Genome Atlas (Lek et al., Reference Lek, Karczewski and Minikel2016).
The Ukrainian cohort consisted of 112 women with breast cancer, 10 women with ovarian cancer and 1 woman with breast and ovarian cancer. Seventy-three women affected with breast cancer met the hereditary cancer criteria and 38 women met the familial cancer criteria. For one Ukrainian participant with breast cancer, insufficient information was available to classify her cancer. Of the women with ovarian cancer, six had hereditary ovarian cancer and four had familial ovarian cancer. The woman with breast and ovarian cancer met the hereditary cancer criteria. There had been no previous testing for mutations in BRCA1 and BRCA2 conducted in the Ukrainian participants.
All participants provided informed consent for participation in this research programme, which was approved by the Commission of Bioethics of the Institute of Hereditary Pathology of the National Academy of Medical Sciences of Ukraine, the Ethics Committee of Wroclaw Medical University (Poland), the Ethics Committee of the University of Rzeszow (Poland) and the University of Melbourne Human Research Ethics Committee (Australia).
2.2. Mutation screening
Amplicon-based massively parallel sequencing of the protein-coding regions and proximal intron–exon junctions of BRCA1 (NM_007294.3) and BRCA2 (NM_000059.3) was performed using lymphocyte-derived germline DNA and the Hi-Plex protocol (Nguyen-Dumont et al., Reference Nguyen-Dumont, Hammet and Mahmoodi2015). All oligonucleotides were synthesized by Integrated DNA Technologies (Coralville, IA, USA). Gene-specific primers and adapter primers were purified to standard desalting and high-performance liquid chromatography grade, respectively. All oligonucleotide sequences are available upon request. Massively parallel sequencing (150 bp paired-end) was performed on the MiSeq system (Illumina, San Diego, CA, USA). Mapping to the human reference build GRCh37 was performed using bwa-mem 0.7.17 (Li & Durbin, Reference Li and Durbin2009). Variant calling was performed using the Java version of VarDict in single-sample, amplicon mode (Lai et al., Reference Lai, Markovets and Ahdesmaki2016).
2.3. Annotation and classification of variants
Variants were annotated using VEP (v.90) and loaded into GEMINI, according to the authors’ recommendations (Paila et al., Reference Paila, Chapman, Kirchner and Quinlan2013; McLaren et al., Reference McLaren, Gil and Hunt2016). Classification of genetic variants in BRCA1 and BRCA2 was then performed in accordance with the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium's recommendations (Spurdle et al., Reference Spurdle, Healey and Devereau2012). Pathogenicity calls were retrieved from the BRCA Exchange portal (http://brcaexchange.org; accessed 17 June 2019). When no information was available (i.e., variants reported as ‘no data’, ‘not available’ or ‘not yet reviewed’), variants were classified in accordance with the ENIGMA criteria. Specifically, truncating variants (i.e., nonsense, frameshift insertions or deletions and variants occurring in consensus splice junctions) were classified as pathogenic/likely pathogenic (P/LP). Missense substitutions that had not been reviewed by ENIGMA were classified as variants of unknown significance (VUS).
3. Results
Targeted sequencing identified a total of 31 and 18 women as carriers of a P/LP variant in BRCA1 and BRCA2, respectively (Tables 2 & 3). The clinical characteristics of these women are available in Supplementary Table S1.
a Genetic variants in BRCA1 and BRCA2 that are classified as pathogenic by the expert panel Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) (Spurdle et al., Reference Spurdle, Healey and Devereau2012) as reported on the BRCA Exchange portal (http://brcaexchange.org), unreported truncating variants and variants occurring in consensus splice sites.
b Variant nomenclature based on +1 as A of ATG start codon, according to the Human Genome Variation Society (HGVS), HGVS_c for coding DNA and HGVS_p for protein variants.
c Classification according to the ENIGMA expert panel (Spurdle et al., Reference Spurdle, Healey and Devereau2012), available from the BRCA Exchange portal (http://brcaexchange.org).
d Pol = Polish; Ukr = Ukrainian.
e BC = breast cancer; OC = ovarian cancer. Some women were diagnosed with both BC and OC.
f Transcript sequences are BRCA1: NM_007294.3 and BRCA2: NM_00059.3.
g Founder mutation – identified via targeted sequencing (this study) or via Sanger sequencing (prior testing) (Table 1).
h No DNA was available for validation by Sanger sequencing.
a Genetic variants in BRCA1 and BRCA2 that are classified as pathogenic/likely pathogenic (P/LP) by the expert panel Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) (Spurdle et al., Reference Spurdle, Healey and Devereau2012) as reported on the BRCA Exchange portal (http://brcaexchange.org), unreported truncating variants and variants occurring in consensus splice sites.
b Classification as per Berliner et al. (Reference Berliner and Fay2007).
c Founder mutations in Table 1.
The 49 P/LP variants with a variant allele fraction >0.2 and total depth ≥10× were verified by Sanger sequencing, except for the case that was identified as carrying the BRCA1 variant NM_007294.3:c.4357+1G>C, as there was insufficient DNA available. That position was covered by 1302 reads, and the variant allele fraction was 0.53 and thus highly unlikely to be a sequencing artefact.
The Polish women participating in this study were previously genotyped for mutations in BRCA1 and BRCA2. The prevalence of BRCA1 founder mutations (c.5266dup, c.181T>G, c.4035del, c.68_69del) identified prior to this study was 10.1% (34/337) (Tables 2 & 3). Targeted sequencing did not identify any additional carrier of the mutations previously tested. We detected five carriers of a non-founder P/LP BRCA1 variant and eight carriers of a P/LP BRCA2 variant (3.9% combined prevalence of non-founder mutations) (Table 3).
Among the Ukrainian women, the prevalence of founder mutations was 13% (16/123) (Table 3). We identified 11 carriers of BRCA1:c.5266dup, 3 carriers of BRCA1:c.181T>G and 2 carriers of BRCA1:c.68_69del. We did not observe BRCA1:c.4035del in the Ukrainian cohort. There were 10/123 carriers of a P/LP variant in BRCA1 and 10/123 carriers of a P/LP variant in BRCA2 (16.3% combined prevalence of non-founder P/LP variants) (Table 3). BRCA1:c.5030_5033del was observed in three unrelated women with breast cancer from Ukraine, all of whom having a family history of cancer (one had a sister with breast cancer and two had mothers with ovarian cancer) (Supplementary Table S1).
Targeted sequencing also identified 11 and 27 rare missense substitutions in BRCA1 and BRCA2, respectively, which are currently classified as VUS (Table 4).
a Missense substitutions in BRCA1 and BRCA2 that are present at less than 1% in ExAC, that have not been reviewed yet or are classified as variants of unknown significance on BRCA Exchange by the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) expert panel (Spurdle et al., Reference Spurdle, Healey and Devereau2012).
b Variant nomenclature based on +1 as A of ATG start codon, according to the Human Genome Variation Society (HGVS), HGVS_c for coding DNA and HGVS_p for protein variants.
c ExAC non-Finnish European population minus The Cancer Genome Atlas (Lek et al., Reference Lek, Karczewski and Minikel2016).
d CADD (Kircher et al., Reference Kircher, Witten, Jain, O'Roak, Cooper and Shendure2014); REVEL (Ioannidis et al., Reference Ioannidis, Rothstein and Pejaver2016).
e Pol = Polish; Ukr = Ukrainian.
f BC = breast cancer; OC = ovarian cancer.
g Transcript sequences are BRCA1: NM_007294.3 and BRCA2: NM_00059.3.
4. Discussion
The overall prevalence of P/LP BRCA1 and BRCA2 variants in the Polish women in this study was 13.9% (47/337) (Table 3). Of these, over a quarter were non-founder mutations (27.7%, 13/47). Our findings are consistent with a recent report by Kowalik et al., who screened for BRCA1 and BRCA2 in Polish women who qualified for genetic testing and identified 161 P/LP variants, 64% (103/161) of which were founder mutations and 24.8% (40/161) of which were non-founder P/LP variants (Kowalik et al., Reference Kowalik, Siolek and Kopczynski2018).
To the best of our knowledge, our study is the first to report the mutation screening of the complete coding regions of BRCA1 and BRCA2 in Ukrainian women affected with breast and/or ovarian cancer. In the Ukrainian women, the overall prevalence of P/LP BRCA1 and BRCA2 variants was 29.3% (36/123) (Table 3). There was no difference in the proportion of founder and non-founder mutations (44.4%, 16/36, and 55.5%, 20/36, respectively).
The higher overall prevalence of P/LP variants observed in the Ukrainian participants (29.3% versus 13.9% in the Polish participants) is likely due to differences in selection criteria. Among the Ukrainian women, 99% (122/123) of participants met the criteria for hereditary (65%, 80/123) or familial cancer (34%, 42/123), whereas the Polish cohort included a majority of sporadic cancers (54%, 182/337) (Berliner et al., Reference Berliner and Fay2007). Hereditary and familial cancers in the Polish cohort accounted for only 38% (127/337) and 8% (28/337) of all participants, respectively (Table 3).
Overall, non-founder variants represented 43.9% (18/41), 62.5% (5/8) and 66.7% (2/3) of all P/LP variants observed in women affected with hereditary, familial and sporadic breast cancer, respectively (Table 3). In women affected with ovarian cancer, non-founder variants represented 14.3% (3/21) and 100% (3/3) of all P/LP variants observed in hereditary and familial cancer, respectively. Expanding genetic testing beyond genotyping for founder variants has thus enabled us to identify 33 women carrying a clinically actionable variant who will be able to receive personalized clinical advice for themselves and their family. These results support the utility of comprehensive gene testing of BRCA1 and BRCA2 in Polish and Ukrainian patients, especially in women with hereditary and familial cancers.
In addition to P/LP variants, our study identified 38/427 carriers (7.7%) of rare missense variants of unknown clinical significance. Missense substitutions may result in variant proteins with functions ranging from normal to severely altered. Therefore, this group of variants is highly likely to be made up of variants with differing levels of associated risks (including none). There are substantial ongoing efforts by ENIGMA to classify VUS in BRCA1 and BRCA2 (Vallee et al., Reference Vallee, Francy and Judkins2012). However, methods such as calibrated functional assays are essential for enabling variant classification, and we currently lack the evidence base from which to interpret and report most missense substitutions.
Our study applied an amplicon-based targeted sequencing methodology that is not designed to detect copy number variations (CNVs). The inherent nature of targeted sequencing poses substantial challenges for the detection of these variants. A number of software tools for CNV detection from targeted sequencing data have recently emerged to try to address this gap (Li et al., Reference Li, Lupat and Amarasinghe2012; Ellingford et al., Reference Ellingford, Campbell and Barton2017; Kerkhof et al., Reference Kerkhof, Schenkel and Reilly2017). However, they are developed for probe-based enrichment rather than amplicon-based methodologies and, to date, multiplex ligation-dependent probe amplification remains the gold standard method for the clinical identification of such events.
5. Conclusions
Our results show that performing comprehensive genetic testing of BRCA1 and BRCA2 instead of testing for founder mutations only will be highly valuable in Poland and Ukraine. Massively parallel sequencing is an effective way of performing comprehensive genetic testing of BRCA1 and BRCA2 that will increase the detection rate of clinically actionable variants and thus enhance risk assessment and management for these women and their families.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0016672320000075
Acknowledgements
We would like to thank Magdalena Puniach for her assistance with translating the Polish data.
Author contributions
TN-D contributed to the study design, performed the sequencing data analyses and drafted the manuscript. PK, MMS, RS HA, AgS, KP and NK collected the study sample and provided clinical information. JAS and DT contributed to the processing of the sequencing data. FH contributed to the genetic testing. DJP and BJP designed the Hi-Plex technology for the panel testing. PK, AlS and HT prepared the DNA samples. AM performed the genetic testing and validations. MCS and AM were responsible for the overall study design and contributed to the data analysis and drafting of the manuscript.
Financial support
This work was supported by the Australian National Health and Medical Research Council (APP1029974 and APP1074383) and by a Victorian Life Sciences Computation Initiative grant (number VR0182) on its Peak Computing Facility, an initiative of the Victorian Government. TN-D is a Career Development Fellow of the National Breast Cancer Foundation (Australia, ECF-17-001). BJP is the recipient of a Victorian Health and Medical Research Fellowship (Australia). MCS is a National Health and Medical Research Council (Australia) Senior Research Fellow. AM was supported by a Research Fellowship within ‘UR – modernity and future of region’ from the European Social Fund, Human Capital, National Cohesion Strategy (contract no. UDA-POKL.04.01.01-00-068/10-00).
Conflicts of interest
The authors declare that they have no conflicts of interest.
Availability of data and material
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.