Spectrum and prevalence of BRCA1/2 germline mutations in Pakistani breast cancer patients: results from a large comprehensive study

Background

Pathogenic germline mutations in BRCA1 and BRCA2 (BRCA1/2) account for the majority of hereditary breast and/or ovarian cancers worldwide. To refine the spectrum of BRCA1/2 mutations and to accurately estimate the prevalence of mutation in the Pakistani population, we studied 539 breast cancer patients selected for family history and age of diagnosis.

Methods

Comprehensive screening for BRCA1/2 germline mutations was performed using state-of-the-art technologies.

Results

A total of 133 deleterious mutations were identified in 539 families (24.7%), comprising 110 in BRCA1 and 23 in BRCA2. The prevalence of BRCA1/2 small-range mutations and large genomic rearrangements was 55.4% (36/65) for families with breast and ovarian cancer, 27.4% (67/244) for families with two or more cases of breast cancer, 18.5% (5/27) for families with male breast cancer, and 12.3% (25/203) for families with a single case of early-onset breast cancer. Nine mutations were specific to the Pakistani population. Eighteen mutations in BRCA1 and three in BRCA2 were recurrent and accounted for 68.2% (75/110) and 34.8% (8/23) of all identified mutations in BRCA1 and BRCA2, respectively. Most of these mutations were exclusive to a specific ethnic group and may result from founder effects.

Conclusions

Our findings show that BRCA1/2 mutations account for one in four cases of hereditary breast/ovarian cancer, one in five cases of male breast cancer, and one in eight cases of early-onset breast cancer in Pakistan. Our study suggests genetic testing of an extended panel of 21 recurrent BRCA1/2 mutations for appropriately selected patients and their families in Pakistan.

Individuals harboring BRCA1/2 germline mutations have high lifetime risks of breast and ovarian cancer. The identification of individuals harboring BRCA1/2 mutations is crucial to assess their cancer risk, consider preventive measures and tailor cancer management strategies.

Several studies have investigated the prevalence of BRCA1/2 small-range mutations and/or large genomic rearrangements (LGRs) with frequencies varying from 17.6% to 29.8% in white populations from Europe and Australia [1,2,3,4,5] and 9.4% to 21.7% in non-whites from Asia [6,7,8]. The prevalence and distribution of BRCA1/2 mutations vary across populations, mainly due to population-specific recurrent or founder mutations. Accurate identification of the population-specific mutation spectrum is therefore the first step towards incorporating appropriate genetic BRCA1/2 testing into clinical practice in a particular population. This information is not fully elucidated in Pakistan, a country with one of the highest rates of breast cancer in Asia.

To date, no large comprehensive studies evaluating the BRCA1/2 mutations have been reported in the Pakistani population and mutations in males have not been identified so far. Small-range mutations were previously reported in 341 unselected breast and 120 ovarian cancer patients, in which the analysis was restricted to a few exons only [9]. We conducted two studies in early-onset and familial breast/ovarian cancer patients from Pakistan. In the initial study the complete coding regions and exon-intron boundaries of BRCA1/2 were screened for small-range mutations in 176 patients [10]. In the other study 120 BRCA1/2 small-range mutations negative patients were screened for LGRs [11]. Other Asian studies also had small sample sizes [12, 13], reported small-range mutations only [14, 15], and/or restricted LGR analyses to a small number of study participants [6, 16, 17].

Here, we refined the spectrum of BRCA1/2 mutations and more precisely estimated the mutation frequencies including small-range mutations and LGRs in 539 early-onset and familial breast cancer patients from Pakistan.

Enrollment of families

Five hundred and ninety-three breast cancer only or breast and ovarian cancer families were enrolled through index breast and/or ovarian cancer patients who presented at the Shaukat Khanum Memorial Cancer Hospital and Research Centre (SKMCH&RC) in Lahore, Pakistan, from September 2004 to August 2015. The recruited families were classified into five risk groups based on family history of breast/ovarian cancer or age at diagnosis (Table 1) as described previously [19]. After enrollment, 54 families were excluded (Fig. 1), leaving 539 families in the study.

Flow diagram for the study participant’s enrollment, screening methods used and BRCA1/2 mutations detected

Full size image

Clinical and histopathological data and comprehensive information on personal and family history of cancer(s), and ethnicity were obtained from all study participants. The Institutional Review Board of the SKMCH&RC approved the study. All study participants signed an informed written consent before providing a blood sample.

BRCA1/2 mutation screening

Genomic DNA was extracted from 9 to 18 ml of whole blood samples, as described previously [20]. The entire coding regions of the BRCA1 (Genbank accession number U14680) and BRCA2 (Genbank accession number U43746) genes including exon-intron boundaries were screened in 139 patient DNA samples for small-range mutations using denaturing high-performance liquid chromatography (DHPLC) analysis as previously described [21, 22]. Each sample revealing variants was sequenced as described elsewhere [18]. Four-hundred female patients selected based on triple-negative vs. non triple-negative breast cancer phenotype had been screened for small-range BRCA1/2 mutations and 33 mutations were described [18]. All small-range mutation-negative patients had been screened for LGRs using multiplex ligation-dependent probe amplification and three LGRs were described [11]. For the current study, families were selected on the basis of family history of breast/ovarian cancer, male breast cancer or age at diagnosis.

Mutation classification

All BRCA1/2 alterations identified in the current study were classified into pathogenic mutations, variants of unknown significance, or polymorphisms. Pathogenic mutations were defined as (i) small-range mutations which affect one or a few nucleotides including frameshift, nonsense, or splice-site mutations and generate a premature termination codon, except BRCA2 exon 27 variants generating a premature termination codon after codon 3010 [23] and (ii) LGRs that span one or more exons. Mutations were designated using the Human Genome Variation Society (HGVS) and the Breast Cancer Information Core (BIC) committee nomenclature.

All identified mutations were searched in various mutation databases including BIC (https://research.nhgri.nih.gov/bic/), ClinVar (http://www.ncbi.nim.nih.gov/clinvar/), LOVD (http://databases.lovd.nl/shared/genes/BRCA2), ARUP (http://arup.utah.edu/database/BRCA/), and BRCA Exchange (http://brcaexchange.org/). Mutations not reported in these databases were considered as novel and specific to Pakistani population.

Statistical analyses

Distribution of clinical and histopathological characteristics between BRCA1/2 carriers and non-carriers were estimated using Fisher’s exact test for categorical variables and the Wilcoxon rank-sum test for quantitative variables. All statistical tests were two-sided. Results were considered significant at a p value of <0.05. All statistical analyses were done using StatXact 4 for Windows (Cytel, Cambridge, USA) and R, version 2.1.

A total of 539 index patients from unrelated families were enrolled and stratified into five risk groups (Table 1). The mean age of disease onset was 35.4 years (range 18-78) for female breast cancer (n=502), 45.4 years (range 23-66) for ovarian cancer (n=30) and 54.5 years (range 27-76) for male breast cancer (n=27) patients.

Spectrum of BRCA1/2 mutations

Evaluation of pooled data from 539 patients yielded 71 distinct pathogenic mutations in 133 families (24.7%) (Table 1). Fifty-three BRCA1 mutations were detected in 110 families (20.4%) and 18 BRCA2 mutations in 23 families (4.3%). Five mutations in BRCA1 (9.4%) and four mutations in BRCA2 (22.2%) were novel (Table 2). The phenotypes of all families carrying BRCA1/2 mutations are presented in Table 3.

Twenty-one (21/71; 29.6%) mutations including 18 in BRCA1 and three in BRCA2 occurred more than once (Fig. 2a, b). These mutations were identified in 83 unrelated families and accounted for 62.4% (83/133) of all families with mutations. The most common BRCA1 mutation was c.3770_3771del (ten Punjabi families), followed by c.5503C>T (nine Punjabi families), exon 1-2 deletion (seven Punjabi families), c.685del (five Punjabi and two Balochi families), c.4485-1G>A (six Punjabi families), c.2269del (one Punjabi and four Pathan families), c.2405_2406del (five Punjabi families), c.4065_4068del (three Punjabi and one Pathan families), c.1793T>G and c.4508C>A (three Punjabi families each), exon 21-24 deletion, c.2603C>G, c.3339_3341del, c.3598C>T, c.5035del, c.5074+1G>A, and c.5361_5362del (two Punjabi families each) and c.2340_2343del (one Pathan and one Mohajir families). The most common BRCA2 mutation was c.5222_5225del (one Punjabi and three Mohajir families), followed by c.92G>A (two Punjabi families) and c.5682C>A (two Pathan families).

Distribution of deleterious germline mutations identified in Pakistani breast/ovarian cancer families across the BRCA1 and BRCA2 genes. Recurrent mutations are marked with (RM). The distribution of breast cancer (BC) and ovarian cancer (OC) in families according to the position of the mutations in BRCA1 (a) and BRCA2 (b) is shown. Regions inferred to be breast cancer cluster regions (BCCRs) and ovarian cancer cluster regions (OCCRs) according to Rebbeck and colleagues [24] are shown at the bottom

Full size image

In addition to the deleterious mutations, 153 (28.4%) distinct BRCA1/2 sequence variants were detected: 79 missense variants, 48 non-coding variants, 24 synonymous variants, one in-frame deletion, and one polymorphic nonsense variant in exon 27 of BRCA2 (data not shown).

BRCA1/2 mutation frequencies

The frequencies of BRCA1/2 mutations by risk group are provided in Table 1. For BRCA1, the highest mutation frequency was noted in families with breast and ovarian cancer (53.8%), followed by families with at least three breast cancer cases (24.8%), families with two breast cancer cases (18.3%), or families with one early-onset breast cancer case (<30 years) (10.8%). For BRCA2, the highest frequency was observed in families with male breast cancer (14.8%).

Patient and tumors characteristics by BRCA1/2 status

BRCA1 carriers (n=110) were more often identified among female patients (99.1% vs. 94.6%, p=0.039) and belonged to the Punjabi ethnic group (81.8% vs. 68.7%, p=0.030) compared to non-carriers (n=406). In contrast, BRCA2 carriers (n=23) were more common among male patients (17.4% vs. 5.4%, p=0.043) and more often belonged to Pathan ethnic group (34.8% vs. 15.5%, p=0.009).

Female breast cancer patients with mutations in BRCA1 (n=106) or BRCA2 (n=19) had a similar mean age of diagnosis (34.0 years (range 21–72) and 37.7 years (range 23-56), respectively, p=0.073, Wilcoxon rank-sum test), which did not differ to that of non-carriers (n=377) (35.7 years (range 18-78). In contrast, male breast cancer patients harboring BRCA2 mutations (n=4) had an older mean age of diagnosis than non-carriers (n=22) (66.5 years (range 54-76) and 52.5 years (range 27-69) years, respectively, p=0.039, Wilcoxon rank-sum test).

BRCA1-associated breast tumors more often were invasive ductal carcinomas (99.0% vs. 91.4%, p=0.004), triple-negative (60.8% vs. 22.6%, p=<0.0001), and of higher tumor grade (grade 3: 94.9% vs. 63.2%, p=<0.0001) compared to tumors of non-carriers. BRCA2-associated breast tumors more often were PR positive compared to tumors of non-carriers (81.8% vs. 57.2%, p=0.025) (data not shown).

To our knowledge, this is the largest Pakistani study that comprehensively investigated the spectrum of BRCA1/2 small-range mutations and LGRs and prevalence of mutations in 539 high-risk families. Mutations were identified in 24.7% (133/539) of families. Eighteen BRCA1 and three BRCA2 mutations were recurrent and accounted for 68.2% and 34.8% of all mutations in BRCA1 and BRCA2, respectively. Nine mutations were specific to the Pakistani population, whereas other mutations had been reported elsewhere.

The most common type of identified mutations were frameshift mutations (60.6%) followed by nonsense mutations (25.4%). These data are consistent with our previous report [10] and a recent worldwide study [25]. In Pakistani patients, BRCA1 mutations were about 5-times more frequent than BRCA2 mutations. A similar distribution was observed in two Asian studies from South India [26] and Saudi Arabia [27] and most studies among white populations [3,4,5, 28]. This is likely due to the predominance of recurrent BRCA1 mutations in these populations. Contradictory results were reported in other Asian studies from China, Hong Kong, Korea, and Indonesia, where BRCA2 mutations were observed at an equal or a higher frequency than BRCA1 mutations [6, 12, 15,16,17].

Among the 133 mutations identified in our study, 18 BRCA1 and three BRCA2 mutations were recurrent, accounting for 68.2% and 34.8% of all mutations in BRCA1 and BRCA2, respectively. The proportion of recurrent BRCA1 mutations to the total number of identified BRCA1 mutations is higher than our previous report [10], which is likely due to the larger size of the present study. Of the identified recurrent mutations, the majority was also reported as recurrent mutations in other populations [1, 4, 25], while few were exclusively identified in a specific ethnic group of Pakistan. Fourteen BRCA1 mutations (c.3770_3771del, c.5503C>T, c.4485-1G>A, c.2405_2406del, c.1793T>G, c.4508C>A, c.2603C>G, c.3339_3341del, c.3598C>T, c.5035del, c.5074+1G>A, c.5361_5362del, exon 1-2 deletion, and exon 21-24 deletion) and one BRCA2 mutation (c.92G>A) were identified only in the Punjabi ethnic group. One BRCA2 mutation (c.5682C>A) was found only in the Pathan ethnic group. Five other recurrent mutations were identified in more than one ethnic group. Our findings imply that a panel of ethnic specific recurrent mutations may be useful for initial screening of high-risk patients from these ethnic groups. Founder effects were previously suggested for six of the 18 recurrent BRCA1 mutations (c.3770_3771del, c.4065_4068del, c.4485-1G>A, c.4508C>A, c.5503C>T, exon 1-2 deletion) [9,10,11], while haplotype analyses of the remaining recurrent mutations have not been performed so far. The high percentage of recurrent BRCA1 mutations facilitates the development of a local, economical, and efficient ethnic-specific genetic testing strategy in Pakistan.

BRCA1/2 mutations were identified in 24.7% of Pakistani breast cancer families. This frequency is higher than that from our initial report (17%) [10], probably due to the larger study size and comprehensive mutation analyses of both genes. This frequency is also higher than those from other Asian reports from Hong Kong, Malaysia, and Korea, ranging from 9.4% to 21.7% [6,7,8, 16, 17]. These findings further support the notion that the BRCA1/2 mutation frequencies vary among different populations. Our data are similar to those reported in white populations, ranging from 17.6% to 29.8% [1, 2, 4]. We found the highest mutation frequency in breast and ovarian cancer families (55.4%), in agreement with previous studies from Pakistan [10], Korea [16], and studies in white populations [4, 28]. We observed a 2.52 fold (53.8% vs. 21.3%) increased occurrence of BRCA1 mutations in breast and ovarian cancer families compared to breast cancer only families, in line with previous reports [1, 4, 6, 28]. Our findings support the notion that the presence of ovarian cancer in Pakistani breast cancer families increases the likelihood for the occurrence of BRCA1 mutation.

In the present study on 27 families with male breast cancer, a BRCA1/2 mutation frequency of approximately 19% was observed, with BRCA2 mutations being about 4-times more common than BRCA1 mutations. Our data are in line with previous studies [4, 14]. This observed frequency is higher than that reported in our initial much smaller study, in which no mutations were identified [10]. In agreement with the National Comprehensive Cancer Network (NCCN) guidelines, our data also warrant BRCA1/2 testing in families with male breast cancer (NCCN Guidelines Version 2.2019).

The main strength of this study is its large size of 539 high-risk families, the comprehensive screening of both genes for small-range mutations and LGRs using highly sensitive methods (allowing the identification of recurrent BRCA1/2 mutations in the Pakistani population and the more accurate estimation of BRCA1/2 mutation frequencies among high-risk families), and the confirmations of mutations in an independent patient’s sample. However, our study also has some limitations. Participants were recruited at one tertiary care cancer center in Lahore, which may have introduced selection bias. Families belonging to Punjabi and Pathan ethnic groups are over-represented and, therefore, mutations in these groups may be over-represented. Nevertheless, Punjabi (44.7%) and Pathan (15.4%) are the most common ethnic groups reported in Pakistan (The World Factbook). Further, our data are based on self-reported ethnicity of study participants, which may lead to a misclassification of the ethnic origin of some of them.

In summary, our study showed that BRCA1/2 mutations account for 24.7% of high-risk breast cancer patients in Pakistan. Our results have important clinical implications, such as personalized treatment with platinum-based or PARP-inhibitor therapy for breast/ovarian cancer patients carrying a pathogenic BRCA1/2 mutation and early detection, surgical prevention, and chemoprevention strategies for their unaffected BRCA1/2 mutation positive relatives. Overall, BRCA1/2 mutations account for one in four patients with a family history of breast cancer/breast and ovarian cancer, one in five patients with male breast cancer, and one in eight patients with early-onset breast cancer. Eighteen mutations in BRCA1 and three in BRCA2 were recurrent and accounted for 68.2% and 34.8% of all identified mutations in BRCA1 and BRCA2, respectively. Our data suggest that BRCA1 testing may be justified for families with multiple female breast cancers, breast and ovarian cancer or early-onset breast cancer and BRCA2 testing for families with male breast cancer from Pakistan. Our findings will help in tailoring cost-effective genetic testing approach for the high-risk Pakistani population or for individuals of Pakistani origin residing in other countries.

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

BIC:

Breast cancer information core

HGVS:

Human genome variation society

LGRs:

Large genomic rearrangements

LOVD:

The Leiden open variation database

SKMCH&RC:

Shaukat Khanum Memorial Cancer Hospital and Research Centre

We are grateful to all study subjects for their participation in this study. We thank the clinicians (Neelam Siddiqui, Mazhar Ali Shah, Narjis Muzaffar, Usman Ahmad, Umm E Kalsoom, Amir Ali Syed, Huma Majeed, Zulqarnain Chaudhry, Muhammad Asad Parvaiz, and Amina Khan) for their help in recruitment of study participants. We thank Jörg Hoheisel for critical reading of the manuscript.

Funding

The study was supported by the Shaukat Khanum Memorial Cancer Hospital and Research Centre (grant number ONC-BRCA-002) and the German Cancer Research Center.

Affiliations

1. Basic Sciences Research, Shaukat Khanum Memorial Cancer Hospital and Research Centre (SKMCH&RC), Lahore, Pakistan
   • Muhammad Usman Rashid
   • , Noor Muhammad
   • , Humaira Naeemi
   • , Faiz Ali Khan
   •  & Sidra Gull
2. Molecular Genetics of Breast Cancer, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 580, 69120, Heidelberg, Germany
   • Muhammad Usman Rashid
   •  & Ute Hamann
3. Clinical Research Office, SKMCH&RC, Lahore, Pakistan
   • Mariam Hassan
   •  & Saima Faisal
4. Leibniz Institute on Aging - Fritz Lipmann Institute, Jena, Germany
   • Sidra Gull
5. Levine Cancer Institute, Charlotte, USA
   • Asim Amin
6. Pathology Department, SKMCH&RC, Lahore, Pakistan
   • Asif Loya

Contributions

MUR contributed to conception and design of the study, patient recruitment and data acquisition. In addition, he was involved in data analysis, interpretation and in drafting and revising the manuscript. NM performed the molecular analyses and contributed to data analysis and interpretation. He was also involved in writing the first draft of the manuscript. HN, FAK and SG performed the molecular analysis. MH and SF were involved in patient recruitment and clinical data acquisition. AA was involved in the recruitment of study subjects, clinical data collection and revising the manuscript. AL was involved in the pathological data acquisition and interpretation. UH contributed to conception and design of the study, data analysis and interpretation and led the writing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ute Hamann.

Ethics approval and consent to participate

This study was conducted in accordance with the Declaration of Helsinki and Good Clinical Practice guidelines. It was approved by the ethics committee of Shaukat Khanum Memorial Cancer Hospital and Research Centre (SKMCH&RC), Lahore Pakistan. The ethics committee name is the “Institutional Review Board”. The approval number is ONC-BRCA-002. Written informed consent was obtained from all study participants.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing of interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions