A total of 772 patients with breast and/or ovarian cancers were included in this study and additional 149 patients were screened only for the 185delAG hotspot mutation. The characteristics of the study subjects are presented in Table 1. We identified 181 (23.4%) patients with clinically significant variants and based on ACMG 2015 guidelines, these variants were classified: 111 pathogenic or likely pathogenic and 70 VUS. The overall diagnostic yield (mutation positivity), including all cases, was 23%. There was a higher yield seen in patients with family history (34.9%) than those without family history (20.4%), (Table 2). The BRCA1 mutation load in this study was comparatively higher with 118 variants, while 68 variants were identified in the BRCA2 gene. We identified five novel variants (two pathogenic, two likely pathogenic, and one VUS). A list of the identified pathogenic and likely pathogenic variants is given in Table 3.
Multiplex PCR based target enrichments followed by sequencing yielded a mean coverage of >500X with more than 98% of the bases covered at a minimum coverage of 20X (yielding coverage of coding regions and splice site regions without any gaps). Sanger sequencing validation results were consistent with the NGS readouts. However, the challenge with the Ion Torrent sequencing technology is the presence of homopolymer errors (false calls at the repeats of the same nucleotide) (31). Therefore, the variants identified in the homopolymer regions were filtered by cross-comparing similar allele frequency false calls across the samples. Also, the unique variant calls in the homopolymer region specific to a sample, even with < 40% allele frequency calls, were further confirmed by Sanger sequencing. With these strategies, we were able to identify and validate variants even in the homopolymers regions in BRCA1 and BRCA2 genes.
In this study, out of 772 patients, 181 patients were identified with clinically significant variants, which consisted of 63 frameshift (34.8%), 26 nonsense (14.4%) and 64 missense (35.4%) mutations. Splice site variants were identified in 23 (12.7%) patients (Figure 3).
Description of Novel Variants
In this study, we have identified five novel BRCA variants.
Subject BC294 was found to be positive for a heterozygous BRCA1 deletion c.5334_5335del (p.Asp1778GlufsTer72) that results in a frameshift at codon 1778, which creates a premature translational stop signal at codon 1849 instead of 1885. The patient had bilateral metachronous breast cancer with a strong family history of first- and second-degree relatives with breast and ovarian cancers.
Subject BC219, with serous cystadenocarcinoma of the ovary, was positive for a novel heterozygous BRCA1 nonsense variant c.1369G>T (p. Glu457Ter), resulting in premature termination. Also, the patient had recurrent metastatic disease and has been counselled to use PARP inhibitors. She had no first or second-degree family members with HBOC (Hereditary Breast and Ovarian Cancer). According to ACMG 2015 guidelines, these novel BRCA1 loss-of-function variants are likely pathogenic. The families of both these patients have been referred for genetic counselling and are on regular monitoring.
The subject BC704 diagnosed with ovarian cancer at 57 years of age, was found to be positive for a novel frameshift BRCA2(NM_000059.4):c.7219del (p.Val2407PhefsTer60) variant, resulting in premature truncation of BRCA2 protein at 2466 instead of 3419. There was no family history of cancers.
Another novel likely pathogenic frameshift BRCA2(NM_000059.4):c.521del (p.Arg174LeufsTer11) variant was identified in the subject BC724 who was diagnosed with estrogen-receptor positive breast cancer and the age of onset was 29. The patient, paternal grandfather’s, brother had prostate cancer.
Subject BC152, is a 13 years old girl from Bangladesh who underwent BRCA testing because of a strong family history with her mother, who was diagnosed with carcinoma breast at the age of 33 years. The proband was found to be positive for a novel BRCA2 variant c.6661A>G (p.Asn2221Asp). The identified BRCA2 variant is not reported in GnomAD exomes or 1000genome databases. However, the mutation status of other family members is unknown and based on ACMG 2015 guidelines this variant has been classified as a variant with uncertain significance.
These results were further validated with sanger sequencing and NGS (Supplementary material Figure 2a, 2b, 3 and 4).
Recurrent Mutations
We identified six recurrent variants in the BRCA1 gene-BRCA1(NM_007300.4):c.68_69del, c.3995G>T, c.5098del, c.5429del, c.55C>T and c.4738G>C and three recurrent variants in the BRCA2 gene- BRCA2(NM_000059.4):c.2593G>C, c.4570_4573del and c.5604_5605del. The recurrent variants with ACMG classification have been listed in Table 4.
The BRCA1(NM_007300.4):c.68_69delAG on the BRCA1 gene, also known as the 185delAG, was identified 29 times out of 115 deleterious BRCA1/2 variants (25.2%). Of these 29 patients, 11 were diagnosed with carcinoma breast and 13 with ovary. Interestingly, out of 67 patients screened due to strong family history, five were positive for 185delAG mutation. This mutation accounted for 25.2% of all pathogenic variants in our cohort.
There were four patients (BC086, BC094, BC116 and BC757 ) in our cohort positive for BRCA1(NM_007300.4):c.3995G>T (p.Gly1332Val) variant in exon 10 of BRCA1 gene; BC086 diagnosed with breast cancer at the age of 26 whose mother had breast cancer and BC094 and BC757 had ovarian cancer at the age of 68 and 70 respectively. BC094with subsequent recurrence at 70 years of age and BC757 had a strong family history of cancers. This missense variant results in the substitution of amino acid Glycine to Valine at position 1332. The majority of the in-silico predictions suggest a damaging outcome and also the minor allele frequency (0.0000653) of 1 in 15306 (in South Asians) and only two heterozygote alleles have been reported in gnomAD (13,14). This variant has not been previously described in affected individuals(15,16). However, in this study, four positive alleles among a cohort of 772 patients suggest that the variant may have a likely pathogenic role in Breast and Ovarian cancer and require further investigations.
There were two patients (BC012 and BC175) positive for BRCA1(NM_007300.4):c.5098del (p.Leu1700Ter) nonsense mutation, which results in the deletion of a single nucleotide from exon 16 causing a translational premature stop signal at codon 1700. Both these patients were diagnosed with early onset (age 36 and 37 respectively) triple negative breast cancer with a strong family history of cancer.
Two patients (BC169 and BC001) were positive for a pathogenic BRCA1(NM_007300.4):c.5429del (p.Ala1810ValfsTer4) mutation in exon 22. BC169 was diagnosed with epithelial ovarian cancer at 44 while BC001 had triple negative breast cancer at 37 . No family history of breast, ovarian, pancreatic or prostate cancer is known in these families.
The pathogenic variant BRCA1(NM_007300.4):c.55C>T (p.Gln19Ter) in exon 2 was identified in two patients (BC023 and BC328) diagnosed with triple negative breast cancer at ages 52 and 37 respectively. The patient BC023 had recurrent breast cancer after ten years and her daughter expired due to bilateral breast cancer. The patient BC328, as per the records has no family history.
Two patients (BC120 & BC172) tested positive for BRCA2(NM_000059.4):c.2593G>C (Glu865Gln) variant in exon 11 with uncertain significance. BC120 (age 38 years) was screened due to a family history of breast cancer, whereas BC172 was diagnosed with estrogen positive breast cancer at 37 years, however, the family history (as per the record) is not known.
The pathogenic variant BRCA1:c.4738G>C was identified in two patients (Subject ID: BC569 and BC587). The patient BC587 was diagnosed with ovarian cancer at the age of 62 years and had a strong family history of cancer. BC569 also had ovarian cancer at 49 and had a family history of cancers.
Another recurrent variant in our cohort was BRCA2:c.4570_4573delTTTC (p.Phe1524IlefsTer18), identified in patients with subject ID BC502 and BC579. BC 502 was diagnosed with Ca ovary at 33 and has no family history known to date while the BC579 has a family history of cancer and is diagnosed with Ca ovary at 60.
Two patients, BC484 and BC571 were positive for a pathogenic BRCA2:c.5604_5605del (p.Asp1868GlufsTer4) frameshift mutation leading to protein truncation. The patient BC484 was diagnosed with Ca ovary at 54 with a history of cancer in the family and BC571 had estrogen-positive breast cancer at 31 with a family history of cancer.
Factors predicting deleterious BRCA1/2 Variants and 185delAG mutation:
We performed a univariate analysis of the factors predictive of the presence of BRCA1 or BRCA2 mutations. Triple negative breast cancer, hormone positive breast cancer, and a family history of BOPP malignancy were significantly associated with higher risks of detecting BRCA1 or BRCA2 mutations. After multivariate binary logistic regression, only a positive family history of BOPP malignancy remained significantly associated with being BRCA mutation positive (Out of 181 positive patients, 58 had a family history of BOPP cancers and 123 did not have any family history (OR 2.15 [95% C.I 1.46 to 3.16] p= <0.0001).
In a logistic regression model for factors predictive of a 185delAG variant, being of South Indian origin was the only significant factor identified (OR 4.93 [95% C.I 95%=1.86 to 13.04, p=0.0013). Out of 29 delAG positive patients, 24 are from south India and 5 are from northern India.
Characteristics of patients positive for variants
There were 111 patients identified with pathogenic/likely pathogenic mutations and the mean age of these patients was (41.2 ±14) years. Out of these 111 patients, 51 were diagnosed with breast cancer while 54 with ovarian cancer, and six were considered for genetic screening due to a strong family history of BOPP cancer.
Interestingly, 89/111 (80%) pathogenic/likely pathogenic mutations were found in the BRCA1 gene and remaining 22/111 (20%) in BRCA2 gene (Table 4). Twenty-seven patients were positive for mutations in exon 2 (3 missense, 22 185delAG mutation and other dels) and 37 patients with mutations in exon 10 of BRCA1 gene in the study cohort. Therefore, BRCA1 exons 2 and 10 can be considered as hotspot exons and by screening these exons alone can yield, identification of 57% of pathogenic/likely pathogenic mutations.
There were 70 patients identified to have a VUS, requiring further investigation to understand the role of these variants in HBOC. The clinical features of patients positive for VUS variants varied with age and family history (Supplementary Material Table 2).