Abstract
OBJECTIVES: In breast cancer (BC) patients, the frequency of germline BRCA mutations (gBRCA) may vary according to the ethnic background, age, and family history of cancer. Phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) is the second most common somatic mutated gene in BC; however, the association of mutations in both genes with cancer has not been thoroughly investigated. Thus, our aims were to investigate gBRCA mutation frequency in a cohort of postmenopausal Brazilian BC patients and the association of gBRCA1/BRCA2 and PIK3CA somatic mutations.
METHODS: Forty-nine postmenopausal (>55 years) and forty-one young (≤35 years) BC patients were included in this study. The postmenopausal group included patients who reported a positive family history of cancer. For these patients, gBRCA1/BRCA2 were sequenced using next-generation sequencing (NGS) or Sanger sequencing. Data for gBRCA in young patients were already available from a previous study. DNA from formalin-fixed, paraffin-embedded (FFPE) tumors was obtained from 27 postmenopausal and 41 young patients for analyzing exons 9 and 20 of PIK3CA. The association between gBRCA1/BRCA2 and somatic mutations in PIK3CA was investigated.
RESULTS: The overall frequency of gBRCA1/BRCA2 among the 49 postmenopausal patients was 10.2%. The frequencies of somatic mutations in PIK3CA in the postmenopausal and young patients were 37% and 17%, respectively (ns). The most common PIK3CA mutation was found to be E454A. Nonsense and frameshift mutations, which may counteract the oncogenic potential of PIK3CA were also detected. Regardless of age, 25% of BRCA1/BRCA2 mutation carriers and non-carriers , each, had PIK3CA somatic mutations.
CONCLUSIONS: Data obtained indicate that BRCA1/BRCA2 gene testing may be considered for postmenopausal patients with BC who have a family history of cancer. Although some of them are not considered pathogenic, somatic variants of PIK3CA are frequently observed in BC patients, especially in postmenopausal patients.
Breast Cancer; Germline Mutation; Somatic Mutation; BRCA1; BRCA2; PIK3CA
BACKGROUND
Breast cancer affects women of all ages; however, the incidence of breast cancer increases with age, and the peak incidence occurs between 45-64 years (1). In addition, breast cancer is the most prevalent cancer in women aged 30-39 years (2). The main risk factors for breast cancer are a) age, b) positive family history of breast and ovarian cancer, and c) hormone exposure (3).
A positive family history is observed in approximately 10-20% of the breast cancer patients, but mutations in predisposing genes have been identified in <30% of these cases (4). BRCA1/BRCA2—both related to the homologous repair of DNA double-strand breaks—are the major breast/ovarian cancer susceptibility genes. Generally, women who harbor BRCA1/BRCA2 mutations are more frequently diagnosed with breast cancer at an early age (≤40 years) or with ovarian cancer at any age. In addition, women who develop breast cancer at an older age and report a strong family history of breast/ovarian cancer mainly in close relatives—first, second, or third degree—may also be BRCA1/BRCA2 mutation carriers (5). However, the majority of breast cancer cases are sporadic, i.e., not related to genetic syndromes. In this case, somatic mutations accumulate over an individual’s lifetime, similar to an ‘evolutionary’ process, a phenomenon that makes age itself a risk factor for cancer (6). In this process, some cells acquire mutations that are advantageous from a tumoral perspective, which allows aberrant proliferation, invasion, and metastasis.
In breast cancer, somatic mutations in the PIK3CA gene are the most frequent, just after TP53 (7). The PIK3CA gene encodes the p110 catalytic subunit of a heterodimeric lipid kinase called PI3K that is activated in response to various extracellular signals that are transduced through receptor tyrosine kinases. After activation, PI3K phosphorylates phosphatidylinositol-4,5-bisphosphate (PI-4,5-P2), generating phosphatidylinositol-3,4,5-trisphosphate (PIP3), which functions as a second messenger and recruits proteins that harbor pleckstrin homology (PH) domains (e.g., AKT) (8). Mutations in the helical or kinase domain of PIK3CA resulted in the activation of the p110a kinase, with the subsequent downstream activation of mediators that culminates in cell proliferation, angiogenesis, and promotion of metastasis (9,10).
In breast cancer, an association between somatic mutations in PIK3CA and the positive expression of the estrogen receptor (ER) has been reported (11-14). However, the association between the frequency of somatic mutations in PIK3CA and age is unclear (15,16). Moreover, it seems likely that the frequency of somatic mutations in PIK3CA increases in ER-positive tumors in aging patients (7).
Thus, BRCA1 and BRCA2 are the most common germline mutated genes, while PIK3CA is the second most common somatic mutated gene in breast cancer patients; however, subtle frequency differences may be related to the age of onset of the disease. Carcinogenic mechanisms elicited by BRCA1/BRCA2 loss of function and PIK3CA gain of function may be targeted for therapy. There is evidence that combination therapies targeting tumors harboring BRCA mutations—such as PARP inhibitors—with PI3K pathway inhibition therapies may exhibit synergy in vivo for the treatment of endogenous BRCA1-related breast cancer mouse model (17). However, it has been previously reported that the frequency of PIK3CA mutations may be different in breast cancer patients based on the presence of germline mutations in BRCA1/BRCA2 (in both women and men) (18,19). Thus, our aim was to investigate the frequency of BRCA mutations in a cohort of postmenopausal Brazilian breast cancer patients, for whom scarce information is available. The secondary exploratory aim of this study was to evaluate the association of germline BRCA1/BRCA2 mutations with somatic PIK3CA mutations in a cohort of young and postmenopausal patients with breast cancer.
METHODS
Patients
Patients were recruited at the Instituto do Câncer do Estado de São Paulo (ICESP), the cancer treatment branch of Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, the largest public hospital complex in Latin America, São Paulo, Brazil. This study was approved by the Institutional Ethics Committee (Comitê de Ética da Faculdade de Medicina da Universidade de São Paulo; protocol 397/11). All patients signed informed consent forms.
The inclusion criteria were 1) histopathological diagnosis of invasive breast carcinoma in patients aged <36 years or >54 years; 2) patients aged 55 years or older with at least one relative having first, second, or third degrees and diagnosed with breast, ovarian pancreatic, or prostate cancer; 3) triple-negative tumor and age ≤60 years. The expression of hormone receptor was classified as positive if at least 1% of the malignant cells were stained with antibodies against estrogen or progesterone receptor; HER2 positivity was defined as immunohistochemistry scores of 3(+) or 2(+), the latter, associated with fluorescence in situ hybridization (FISH)-amplification. HER2 immunohistochemistry and FISH were scored according to the ASCO/CAP guidelines (20). The Ki67 expression cut-off was set at >14% for a high proliferation index. The molecular subtypes were classified using previously established criteria (21).
Personal and familial cancer histories were collected through a structured questionnaire. Patients were also asked about their ancestry to obtain information about the country or continent where their parents and grandparents (at least) were born. A pedigree that reached up to third-degree relatives was designed. Clinical and pathological data were retrieved from hospital files.
In a previous study, 79 very young breast cancer patients (≤35 years) were evaluated for the presence of germline mutations in BRCA1 and BRCA2, among whom, four harbored BRCA1 mutations (c.66_67insA; c.211A>G; c.3331_3334delCAAG; c.5263_5264insC) and nine harbored BRCA2 mutations (c.483T>A; c.1138_1138delA; c.2808_2811delACAA (n=2); c.3956_3959delATGA; c.6656C>G; c.6990_6994delTACCT; c.9154C>T; c.9382C>T) (22). For detecting PIK3CA mutations, tumor samples were available for 41 patients (among the 79 patients) and were included in the present analysis. Clinical data and tumor subtypes based on ER, PR, HER2, and Ki67 expression levels (as described above) are summarized in Table 4 (22). Six of these forty-one patients harbored BRCA1 or BRCA2 mutations.
DNA Extraction from the Blood and Tumor Tissue
Genomic DNA from peripheral blood samples was extracted using the Illustra Blood Genomic Prep Mini Spin Kit (GE Healthcare Bio-Sciences, Pittsburgh, PA, USA), and from cancer cell-enriched areas from the formalin-fixed, paraffin-embedded (FFPE) tumor samples using the QIAamp® DNA FFPE Tissue (Qiagen, Valencia, CA, USA), as per the manufacturer's protocol.
DNA concentration and purity were determined using a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific, Massachusetts, USA), and the absorbance260/280 ratio varied from 1.42 to 2.2. DNA concentration from samples analyzed using next-generation sequencing (NGS) was also evaluated using a Qubit® dsDNA BR Assay kit on a Qubit® 3.0 Fluorometer (Invitrogen, Carlsbad, California, USA).
Analysis of Germline Mutations in BRCA1/BRCA2
The entire coding regions of BRCA1 and BRCA2, including exon-intron boundaries, were sequenced by NGS using the Ion Torrent Personal Genome Machine (PGM) platform (n=38) or by Sanger sequencing (n=11), to determine the presence of germline mutations.
Next-Generation Sequencing
BRCA1 and BRCA2 were sequenced using the Ion AmpliSeq™ BRCA1 and BRCA2 Panel (Life Technologies, Carlsbad, CA, USA) consisting of three primer pools, covering the target regions in 167 amplicons that target the entire coding region, including 10-20 bp of non-coding sequences, flanking the 5’ and 3’ ends of each exon, for both genes. Libraries containing the PCR product were sequenced on a 314 v2 Ion Chip, which allows the simultaneous analysis of 12 samples per chip on a PGM sequencer (Ion Torrent™), and the Ion PGM Sequencing 200 Kit version 2 (Life Technologies, Carlsbad, CA, USA). Data analysis was performed using the Ion Reporter™ Server System (Thermo Fisher Scientific, Massachusetts, USA). Sequence data were also visually evaluated using the Integrative Genomics Viewer (IGV). Amplicons with coverage less than 30x, pathogenic variants, and new variants were confirmed by PCR followed by conventional bidirectional Sanger sequencing. Full details of the methods are provided in the Appendix.
PCR and Sanger Sequencing
All coding regions, including the intron-exon boundaries of BRCA1 (NM_7294.3) and BRCA2 (NM_000059.3) were amplified by PCR. Primers and conditions are described in the Appendix. The amplicons were purified (Illustra™ ExoStar™ 1-Step-GE Healthcare Bio-Sciences, Pittsburgh, PA, USA) and were sequenced using the BigDyeTM Terminator v3.1 Cycle Sequencing kit (Applied BiosystemsTM, Foster City, California, USA), as described previously (22). Following purification, samples were analyzed on a 3500 Genetic Analyzer or ABI 3730 DNA Analyzer (Applied Biosystems™, Foster City, California, USA) in both forward and reverse directions (Appendix). The results were analyzed using Mutation Surveyor DNA Variant Analysis Software (v3.30, SoftGenetics LLC). All pathogenic mutations were confirmed using Sanger sequencing.
Analysis of Copy Number Variation in BRCA1 and BRCA2
For the analysis of large deletions and duplications—that would have provided comprehensive information regarding germline mutations—patient DNA was subjected to BRCA1 and BRCA2 multiplex ligation-dependent probe amplification (MLPA) analysis (BRCA1: SALSA® MLPA® P002 and P087 Probemix; BRCA2: SALSA® MLPA® P045 BRCA2/CHEK2 Probemix; MRC-Holland, Amsterdam, The Netherlands), as per the manufacturer's protocols (Appendix), as described previously (22,23).
Mutation Nomenclature and Classification
BRCA1 and BRCA2 variants were named according to the Human Genome Variation Society (HGVS) nomenclature (24) and were searched in publicly accessible databases, i.e., BRCA Share™, BRCA Exchange, BRCA Mutation Database, and ClinVar. The search was performed in 2020 (between January and June). In silico analyzes were performed using the following prediction tools: Polymorphism Phenotyping (PolyPhen; v2.2.2), Sorting Intolerant From Tolerant (SIFT; v1.0.3), Align-GVGD, Protein Variation Effect Analyzer (Provean; v1.1), and Human Splicing Finder to analyze variants of unknown clinical significance. Minor allele frequency (MAF) was checked on the 1000 Genomes Project database, Exome Aggregation Consortium (ExAC), Global MAF dbSNP, Exome Variant Server, NHLBI GO Exome Sequencing Project (ESP), Genome Aggregation Database (gnomAD), Trans-Omics for Precision Medicine (TOPMed), and Brazilian genomic variants (ABraOM). More details are provided in the Appendix.
The variants were then classified as pathogenic, likely pathogenic, benign, likely benign, and variant of uncertain significance (VUS) based on the recommendations of the American College of Medical Genetics and Genomics (25). VUS for BRCA was also checked for co-occurrence with known pathogenic mutations in the same patient. For some variants, we considered that consensus information in ≥2 databases was strong enough to classify them as benign or VUS.
Analysis of Somatic Mutations in PIK3CA
Among the 49 postmenopausal patients, 27 FFPE tumor samples were available for analysis. Tumor samples from another 22 patients were not available because they had been operated on at another service. Tumor samples from all 41 young patients were used for further analysis (Figure 1).
PIK3CA (NM_006218.2) exons 9 (helical domain) and 20 (kinase domain), which are the regions with the highest mutation frequency (26), were amplified by PCR and were analyzed by Sanger sequencing. Primer sets were designed using software Primer3 (http://bioinfo.ut.ee/primer3/). To avoid non-specific product formation, BLAST (http://www.ncbi.nlm.nih.gov/blast) and BLAT (https://genome.ucsc.edu/cgi-bin/hgBlat) were performed. Primers and conditions are described in the Appendix.
Statistical Analysis and Sample Size Calculation
To detect the frequency of germline BRCA mutations in postmenopausal breast cancer patients (varying from 2% to 17%), a sample size of 50 was estimated (27,28). For analyzing the frequency of PIK3CA mutations in young and postmenopausal patients; this was a convenient sample size, because only 55% of tumor samples were available for the latter. Assuming that the frequency of PIK3CA mutations in young and postmenopausal patients was 7% and 35%, respectively (7) and the correlation of two postmenopausal patients for every three young patients, the estimated sample size to detect a difference with 0.05 one-sided significance level and 80% power would be 31 young and 21 postmenopausal patients.
Pearson’s chi-square test was used to evaluate the association between variables, and a two-sided significance level of 0.05 was considered.
RESULTS
Patients
Forty-nine elderly women aged ≥55 years who were diagnosed with invasive ductal breast carcinoma were included between May 2014 and May 2015 and evaluated for the presence of germline mutations in BRCA. FFPE tumor samples of 27 patients were analyzed for the presence of somatic mutations in PIK3CA. The median ages at the time of diagnosis and enrollment in the study were 61 years (55-80 years) and 64 years (56-87 years), respectively. The majority of the patients had Nottingham histological grade II tumors (63.3%) and clinical stage I/II tumors (67.4%). With respect to the tumor subtype, most tumors were luminal B (44.9%) or luminal A (22.4%), followed by HER2+ and triple-negative tumors (10.2% each) (Table 1; Additional Table 1). Most patients (95.9%)—except for two patients (one with a triple-negative tumor and age ≤60 years)—reported a positive family history of breast, ovarian, pancreatic, or prostate cancers. A large proportion of the patients (69.4%) reported at least one affected first-degree family member with breast and/or ovarian cancer. Most women were born in the Southeast (67.3%)—followed by the Northeast (18.4%)—regions of Brazil. With respect to ancestry, 28.6% of the patients reported Brazilian and European ancestries, 26.5% reported only Brazilian ancestry, and 18.4% and 8.4% reported European-only or Asian ancestry, respectively (Table 1).
Clinical and pathological features of breast cancer patients according to deleterious BRCA1 and BRCA2 mutations.
Another 41 young patients, aged ≤35 years, had their tumor samples evaluated for the presence of somatic mutations in PIK3CA. This is a subgroup of patients whose clinical data, as well as germline BRCA1 and BRCA2 sequencing results had already been reported in a previous study (22). The cohort of patients now reported comprehends those young patients who had FFPE tumor samples available for PIK3CA analysis. The median age at the time of diagnosis was 32 years (range, 23-35 years). Most patients presented tumors with histological grade II (43.9%) or III (48.8%), and disease clinical stage I/II (65.7%). Luminal B (46.3%) was the most frequent tumor subtype, followed by triple-negative (24.4%) and HER2 (+) (12.2%) tumors (Table 4). Among these patients, 14.6% and 12.2% reported first-or second-degree relatives diagnosed with breast and/or ovarian cancer, respectively, while 39% reported a negative family history of breast and/or ovarian cancer, and 24.4% were not able to describe their family history. Six out of the forty-one patients harbored pathogenic mutations (14.6%) in BRCA1 or BRCA2, as previously reported (22).
Germline Mutations in BRCA1 and BRCA2 in Postmenopausal Patients
Among 49 postmenopausal unrelated women, 5 (10.2%) were identified to harbor mutations of clinical significance, 3 in BRCA1 and 2 in BRCA2 (Table 2; Additional Tables 2-3). All five BRCA mutations were identified among 47 patients who reported a positive family history of breast, ovarian, prostate, and pancreatic cancers in close relatives (10.6%), including four mutations detected among 34 patients reporting first-degree relatives affected by these types of cancer (11.76%) (Table 1).
Mutations in BRCA1 comprised one splice-site variant (c.5074+2T>C, in exon 17), one missense mutation (c.5123C>A), and one BRCA1 rearrangement generating a large deletion encompassing exons 1-19. The two pathogenic mutations in BRCA2 included one missense variant (c.2T>G) and one nonsense variant (c.5645C>A) (Table 2). The presence of CHEK2 c.1100delC mutation was investigated in 47 out of the 49 patients; however, no mutations were detected.
Eight VUS were detected, five in BRCA1 and three in BRCA2. Among the VUS, four distinct missense variants were identified, two in each gene (BRCA1: c.3305A>G and c.3752G>A; BRCA2: c.3371A>G and c.8942A>G), among which BRCA2 c.3371A>G was predicted to be deleterious by at least three out of four mutation function prediction models (SIFT, Polyphen-2, Align-GVGD, or Provean) (Table 3). The remaining VUS were located in the intronic regions, at least 36 nucleotides away from the intron-exon boundary.
In silico analysis of VUS identified in BRCA1 and BRCA2 using mutation function prediction models.
Presence of Somatic Mutations in PIK3CA in Postmenopausal and Young Patients
Tumor sequencing was performed on samples from 27 elderly patients to identify PIK3CA mutations. Fourteen tumors (51.8%) were found to harbor mutations in exons 9 or 20; however, only ten (37%) harbored meaningful deleterious or possibly deleterious variants (pathogenic in at least one out of four function prediction tests). Recurrent mutations were E545A (observed in four samples) and H1047L (in the other two samples). Among these 27 elderly patients, two were BRCA1 mutation carriers, both of whom harbored somatic pathogenic (E545A) or possibly pathogenic PIK3CA mutations (Additional Table 4).
Another three tumors (11.1%), all luminal A subtypes, harbored synonymous variants (in one case, associated with an intronic variant) (sample 39). In addition, tumors from another six patients harbored multiple PIK3CA variants; however, two tumors harbored (samples 26 and 39) a combination of non-pathogenic variants represented by missense non-pathogenic and nonsense variants (sample 26) or a combination of a deep intronic and two synonymous variants (sample 39). In the third tumor, PIK3CA double mutation (sample 47) (S541P and E1037V) was considered pathogenic in at least three function prediction tests, even though none of them were located in a hotspot. In the fourth and fifth tumors (samples 36 and 46), the contribution of the mutations were difficult to define because the PIK3CA pathogenic missense variant (E545A) was accompanied by a frameshift (FS) mutation (S553FS). If it occurs in cis, FS S553FS might counteract the oncogenic potential of E545A. The sixth tumor (sample 8) harbored a pathogenic hotspot (H1047L) and a synonymous variant (Additional Table 4).
In a cohort of young patients, PIK3CA variants were observed in 12 tumors, including synonymous variants—detected in two tumors (one luminal B, sample 484, and one HER2+ sample 503)—and missense non-pathogenic variants detected in another two samples (samples 455 and 478). In addition, a nonsense variant, W552* was detected in a luminal A tumor (sample 468). Hence, pathogenic or possibly pathogenic PIK3CA mutations were detected in seven out of forty-one young patients (17.1%) (Additional Table 5).
Among the young patients, E545A was the most frequent mutation (detected in three different samples, one luminal B and two triple-negative tumors). In one of these triple-negative tumors, E545A occurred concomitantly with N1068T, another pathogenic variant. The variant P539S, considered pathogenic in the prediction models, was detected in two luminal B samples, in one of these cases, in combination with R555K, which is also a pathogenic variant.
We then compared frequency of pathogenic PIK3CA mutation in tumors from postmenopausal and young patients (37% vs. 17%); however, we could not find a significant difference (Table 4). Using our data with a sample size of 27 postmenopausal and 41 young women and the reported frequency of PIK3CA mutation, the power to detect a difference with a one-sided significance level of 0.05% was 58.51%.
The frequency of PIK3CA is enriched in ER-positive tumors, and in a previous study we detected a trend toward a higher frequency of PIK3CA mutations in ER-positive tumor from elderly women compared to that observed in younger women (7). Upon considering the characteristics of the patients in the present series, we observed differences between the two groups, reflecting a higher proportion of luminal tumors in postmenopausal women. We then analyzed the frequency of PIK3CA mutations in luminal tumors and observed that eight out of the twenty-four samples (33.3%) from postmenopausal patients and five out of the twenty-five samples (20%) from young patients harbored pathogenic PIK3CA mutations (Table 4; p=0.291). A future meta-analysis including more recent data may help to clarify this aspect.
We next considered a total of 68 patients, postmenopausal as well as young, who were tested for the presence of germline mutations in BRCA1/BRCA2 and somatic mutations in PIK3CA. Upon simultaneously considering patients from both age groups, two out of eight germline BRCA1/BRCA2 mutant carriers (25%) were also found to harbor somatic mutations in PIK3CA. Among the 60 patients who were BRCA1 and BRCA2 wild type, 15 manifested tumors harboring PIK3CA mutations (25%).
DISCUSSION
In this cohort of postmenopausal breast cancer patients, 10.2% harbored pathogenic germline BRCA1/BRCA2 variants; 11.7% of these patients had at least one family member who was affected with breast, ovarian, prostate, or pancreatic cancer.
Age at the onset of breast cancer and a family history of breast and ovarian cancer are important factors associated with the frequency of germline BRCA mutations (29). For elderly patients who were not selected for a family history of cancer, the frequency of BRCA mutations tended to be relatively low. Accordingly, a recent nested case-control study conducted in the USA revealed that only 1.18% of the unselected postmenopausal breast cancer patients were BRCA1/BRCA2 mutation carriers (27). In a large cohort comprising 1554 Brazilian breast cancer patients referred for genetic testing at a single clinical diagnostic laboratory in Brazil, 9.84% were found to be BRCA1 or BRCA2 mutation carriers independent of age (30). Higher BRCA mutation frequencies (varying from 15% to 22%) have been reported among young Brazilian breast cancer patients with ages up to 35 years (22,31,30). However, specifically for postmenopausal Brazilian patients with breast cancer, little data are available. Our study indicates that 10.6% of the breast cancer patients with at least one close relative affected by the disease (until third degree) harbor germline BRCA1/BRCA2 mutations. A previous study evaluated 39 breast cancer patients aged more than 50 years, among whom 17.9% were BRCA mutation carriers (32). These latter patients reported a strong family history based on the early age of cancer onset or multiple relatives with breast cancer and/or ovarian cancer at any age, which may explain the higher BRCA mutation frequency.
An important issue to take into consideration is the cost-effectiveness of the diagnostic program for germline mutations in BRCA1/BRCA2 genes and preventative strategies for relatives of patients diagnosed with the mutation. In the scenario of Brazilian ovarian cancer patients, for whom BRCA1/BRCA2 mutation frequency is 20%, performing genetic testing and adopting prophylactic measures for family members was considered a cost-effective measure (33). In a more inclusive model, BRCA testing may be offered to women of the general population to avoid missing mutation carriers, owing to test indications based on clinical criteria and family history. In this context, population-based BRCA testing was estimated to be cost-effective for the Brazilian population and to prevent a large number of breast and ovarian cancer cases (34). Although direct studies for postmenopausal Brazilian breast cancer patients are necessary, the previous two studies might suggest that genetic testing may be valuable for these women in the context of a positive family history.
The variants detected in the present study were not among the most frequent mutations in BRCA1 and BRCA2 in Brazilian patients with breast cancer. Variants BRCA1 c.5074+2T>C, BRCA1 c.5123C>A, and BRCA2 c.2T>G respectively represent 2.2%, 0.5%, and 1.2% of the BRCA1/BRCA2 mutations previously reported (28).
The other two BRCA mutations, BRCA1 large rearrangement (del exons 1-19) and BRCA2 c.5645C>A, have not been previously reported in the Brazilian population. The variant, BRCA2 c.5645C>A has been reported in breast cancer patients from Japan, China, and the Czech Republic (35,36,37), and in prostate cancer patients (38). Interestingly, our patient who harbored this variant was also born in Japan.
Somatic mutations in PIK3CA gene are the second most common mutations in breast cancer, just after TP53 (7). The PIK3CA mutation hotspots were clustered in exon 9 in nucleotides corresponding to codons E542K and E545K (helical domain) and in exon 20 in nucleotides corresponding to codon H1047R (kinase domain) (39,40).
In the present series, the most frequent mutation in PIK3CA in tumors from both postmenopausal and young patients was E545A, a variant with intermediate oncogenic potency, located in the helical domain (39). In agreement with our data, studies on breast cancer patients from Singapore and Peru have also found E545A to be the most frequent PIK3CA variant in tumor samples (41,42). Nevertheless, a method was developed to specifically enhance the detection of E454A (43). In contrast, data from another cohort of Brazilian patients with sporadic breast cancer have reported that the most frequent PIK3CA hotspot mutations were E542K, E545K, and H1047R (13).
The second most commonly found mutations in elderly patients were H1047L and S553FS. H1047L is located in the kinase domain and is associated with high oncogenic potential (39). Further, the frameshift mutation S553FS might counteract the proto-oncogene potential of PIK3CA. In addition, nonsense mutations were detected in tumors from both elderly and young patients, which might also neutralize the proto-oncogenic activity of PIK3CA. However, another study has specified that nonsense mutations in PIK3CA are not frequently encountered (44).
Six tumors were found to harbor double or triple PIK3CA variants (four from elderly patients and two from young patients). It has been previously shown that approximately 13% of all the PIK3CA mutations correspond to multiple variants occurring in the same tumor. It has also been reported that most double mutations occur in cis and induce the activation of the downstream PI3K pathway (compared to single-hotspot mutants) (40). However, in the present study, among the four tumors in elderly patients harboring double or triple variants, only one might be deleterious, involving a combination of S541P and E1037V. In the other three tumors, concomitant variants included nonsense, frameshift, synonymous, and intronic variants, in addition to missense variants with pathogenic or non-pathogenic potential. The determination of whether these variants were in cis might have helped to determine the oncogenic potential of the combinations because if a driver mutation occurred in trans, the effect of the driver mutation might have prevailed.
In the present cohort of patients, somatic mutations in PIK3CA were detected in 25% of the patients harboring germline BRCA1/BRCA2 mutations (two of the eight postmenopausal patients were analyzed for the presence of both gene mutations). This finding may be attributed to the small sample size. In other studies, the frequency of the combination of both mutations appeared to be less than that of individual mutations. In Chinese breast cancer patients, PIK3CA somatic mutations were detected in 14% and 43% of the patients harboring germline BRCA1/BRCA2 mutations (vs. wild type carriers), respectively (18). PIK3CA somatic mutations were not detected in male patients with breast cancer who harbored BRCA2 mutations (19).
Although we were not able to identify any associations between the germline BRCA and somatic PIK3CA mutations because of the small number of patients involved in this study, this is an intriguing situation involving two genes that are treatment targets; therefore, this information may be aggregated in future studies.
The limitations of our study are the small sample size and the sequencing of hotspots (but not all exons of PIK3CA), which may have resulted in the underestimation of the mutation frequency. The strengths of this study are the combined analysis of germline BRCA1/BRCA2 and somatic PIK3CA mutations in a group of postmenopausal and young patients with breast cancer.
In conclusion, the present data indicate that BRCA1/BRCA2 sequencing may be considered for postmenopausal breast cancer patients having a family history of cancer. In addition, although the frequency of PIK3CA variants in exons 9 and 20 is high in both elderly and young patients, some of these variants may not be pathogenic in the context of breast cancer.
ACKNOWLEDGMENTS
We acknowledge the helpful assistance of Dr. Rossana Veronica Mendoza Lopez for the statistical analysis. The São Paulo Research Foundation (FAPESP, grant #2012/12306-4) supported this work. Simone Maistro received a postdoctoral scholarship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, PVE #029/2012), and received a PhD scholarship grant from the São Paulo Research Foundation (FAPESP, #2011/09572-1). Tauana Rodrigues Nagy and Gláucia Fernanda de Lima Pereira received a scholarship grant from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior. Maria Aparecida Azevedo Koike Folgueira received a research grant from the Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazil (CNPq-308876/2017-2).
REFERENCES
-
1 DeSantis CE, Ma J, Goding Sauer A, Newman LA, Jemal A. Breast cancer statistics, 2017, racial disparity in mortality by state. CA Cancer J Clin. 2017;67(6):439-48. https://doi.org/10.3322/caac.21412
» https://doi.org/10.3322/caac.21412 -
2 Ward EM, Sherman RL, Henley SJ, Jemal A, Siegel DA, Feuer EJ, et al. Annual Report to the Nation on the Status of Cancer, Featuring Cancer in Men and Women Age 20-49 Years. J Natl Cancer Inst. 2019;111(12):1279-97. https://doi.org/10.1093/jnci/djz106
» https://doi.org/10.1093/jnci/djz106 -
3 Sun YS, Zhao Z, Yang ZN, Xu F, Lu HJ, Zhu ZY, et al. Risk Factors and Preventions of Breast Cancer. Int J Biol Sci. 2017;13(11):1387-97. https://doi.org/10.7150/ijbs.21635
» https://doi.org/10.7150/ijbs.21635 -
4 Wooster R, Weber BL. Breast and ovarian cancer. N Engl J Med. 2003;348(23):2339-47. https://doi.org/10.1056/NEJMra012284
» https://doi.org/10.1056/NEJMra012284 -
5 Malone KE, Daling JR, Doody DR, Hsu L, Bernstein L, Coates RJ, et al. Prevalence and predictors of BRCA1 and BRCA2 mutations in a population-based study of breast cancer in white and black American women ages 35 to 64 years. Cancer Res. 2006;66(16):8297-308. https://doi.org/10.1158/0008-5472.CAN-06-0503
» https://doi.org/10.1158/0008-5472.CAN-06-0503 -
6 Stratton MR, Campbell PJ, Futreal PA. The cancer genome. Nature. 2009;458(7239):719-24. https://doi.org/10.1038/nature07943
» https://doi.org/10.1038/nature07943 - 7 Encinas G, Maistro S, Pasini FS, Katayama ML, Brentani MM, Bock GH, et al. Somatic mutations in breast and serous ovarian cancer young patients: a systematic review and meta-analysis. Rev Assoc Med Bras (1992). 2015;61(5):474-83.
-
8 Lien EC, Dibble CC, Toker A. PI3K signaling in cancer: beyond AKT. Curr Opin Cell Biol. 2017;45:62-71. https://doi.org/10.1016/j.ceb.2017.02.007
» https://doi.org/10.1016/j.ceb.2017.02.007 -
9 Ikenoue T, Kanai F, Hikiba Y, Obata T, Tanaka Y, Imamura J, et al. Functional analysis of PIK3CA gene mutations in human colorectal cancer. Cancer Res. 2005;65(11):4562-7. https://doi.org/10.1158/0008-5472.CAN-04-4114
» https://doi.org/10.1158/0008-5472.CAN-04-4114 -
10 Isakoff SJ, Engelman JA, Irie HY, Luo J, Brachmann SM, Pearline RV, et al. Breast cancer-associated PIK3CA mutations are oncogenic in mammary epithelial cells. Cancer Res. 2005;65(23):10992-1000. https://doi.org/10.1158/0008-5472.CAN-05-2612
» https://doi.org/10.1158/0008-5472.CAN-05-2612 -
11 Maruyama N, Miyoshi Y, Taguchi T, Tamaki Y, Monden M, Noguchi S. Clinicopathologic analysis of breast cancers with PIK3CA mutations in Japanese women. Clin Cancer Res. 2007;13(2 Pt 1):408-14. https://doi.org/10.1158/1078-0432.CCR-06-0267
» https://doi.org/10.1158/1078-0432.CCR-06-0267 -
12 Cizkova M, Susini A, Vacher S, Cizeron-Clairac G, Andrieu C, Driouch K, et al. PIK3CA mutation impact on survival in breast cancer patients and in ERα, PR and ERBB2-based subgroups. Breast Cancer Res. 2012;14(1):R28. https://doi.org/10.1186/bcr3113
» https://doi.org/10.1186/bcr3113 - 13 Mangone FR, Bobrovnitchaia IG, Salaorni S, Manuli E, Nagai MA. PIK3CA exon 20 mutations are associated with poor prognosis in breast cancer patients. Clinics (Sao Paulo). 2012;67(11):1285-90.
-
14 Liedtke C, Cardone L, Tordai A, Yan K, Gomez HL, Figureoa LJ, et al. PIK3CA-activating mutations and chemotherapy sensitivity in stage II-III breast cancer. Breast Cancer Res. 2008;10(2):R27. https://doi.org/10.1186/bcr1984
» https://doi.org/10.1186/bcr1984 -
15 Kalinsky K, Jacks LM, Heguy A, Patil S, Drobnjak M, Bhanot UK, et al. PIK3CA mutation associates with improved outcome in breast cancer. Clin Cancer Res. 2009;15(16):5049-59. https://doi.org/10.1158/1078-0432.CCR-09-0632
» https://doi.org/10.1158/1078-0432.CCR-09-0632 -
16 Loi S, Michiels S, Lambrechts D, Fumagalli D, Claes B, Kellokumpu-Lehtinen PL, et al. Somatic mutation profiling and associations with prognosis and trastuzumab benefit in early breast cancer. J Natl Cancer Inst. 2013;105(13):960-7. https://doi.org/10.1093/jnci/djt121
» https://doi.org/10.1093/jnci/djt121 -
17 Juvekar A, Burga LN, Hu H, Lunsford EP, Ibrahim YH, Balmaãè J, et al. Combining a PI3K inhibitor with a PARP inhibitor provides an effective therapy for BRCA1-related breast cancer. Cancer Discov. 2012;2(11):1048-63. https://doi.org/10.1158/2159-8290.CD-11-0336
» https://doi.org/10.1158/2159-8290.CD-11-0336 - 18 Chen B, Zhang G, Li X, Ren C, Wang Y, Li K, et al. Comparison of BRCA versus non-BRCA germline mutations and associated somatic mutation profiles in patients with unselected breast cancer. Aging (Albany NY). 2020;12(4):3140-55.
-
19 Deb S, Do H, Byrne D, Jene N; kConFab Investigators, Dobrovic A, et al. PIK3CA mutations are frequently observed in BRCAX but not BRCA2-associated male breast cancer. Breast Cancer Res. 2013;15(4):R69. https://doi.org/10.1186/bcr3463
» https://doi.org/10.1186/bcr3463 -
20 Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J Clin Oncol. 2013;31(31):3997-4013. https://doi.org/10.1200/JCO.2013.50.9984
» https://doi.org/10.1200/JCO.2013.50.9984 -
21 Azim HA Jr, Michiels S, Bedard PL, Singhal SK, Criscitiello C, Ignatiadis M, et al. Elucidating prognosis and biology of breast cancer arising in young women using gene expression profiling. Clin Cancer Res. 2012;18(5):1341-51. https://doi.org/10.1158/1078-0432.CCR-11-2599
» https://doi.org/10.1158/1078-0432.CCR-11-2599 -
22 Encinas G, Sabelnykova VY, de Lyra EC, Hirata Katayama ML, Maistro S, de Vasconcellos Valle PWM, et al. Somatic mutations in early onset luminal breast cancer. Oncotarget. 2018;9(32):22460-79. https://doi.org/10.18632/oncotarget.25123
» https://doi.org/10.18632/oncotarget.25123 -
23 Maistro S, Teixeira N, Encinas G, Katayama ML, Niewiadonski VD, Cabral LG, et al. Germline mutations in BRCA1 and BRCA2 in epithelial ovarian cancer patients in Brazil. BMC Cancer. 2016;16(1):934. https://doi.org/10.1186/s12885-016-2966-x
» https://doi.org/10.1186/s12885-016-2966-x -
24 den Dunnen JT, Dalgleish R, Maglott DR, Hart RK, Greenblatt MS, McGowan-Jordan J, et al. HGVS Recommendations for the Description of Sequence Variants: 2016 Update. Hum Mutat. 2016;37(6):564-9. https://doi.org/10.1002/humu.22981
» https://doi.org/10.1002/humu.22981 -
25 Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405-24. https://doi.org/10.1038/gim.2015.30
» https://doi.org/10.1038/gim.2015.30 -
26 Dirican E, Akkiprik M, Özer A. Mutation distributions and clinical correlations of PIK3CA gene mutations in breast cancer. Tumour Biol. 2016;37(6):7033-45. https://doi.org/10.1007/s13277-016-4924-2
» https://doi.org/10.1007/s13277-016-4924-2 -
27 Kurian AW, Bernhisel R, Larson K, Caswell-Jin JL, Shadyab AH, Ochs-Balcom H, et al. Prevalence of Pathogenic Variants in Cancer Susceptibility Genes Among Women With Postmenopausal Breast Cancer. JAMA. 2020;323(10):995-7. https://doi.org/10.1001/jama.2020.0229
» https://doi.org/10.1001/jama.2020.0229 -
28 Palmero EI, Carraro DM, Alemar B, Moreira MAM, Ribeiro-Dos-Santos Â, Abe-Sandes K, et al. The germline mutational landscape of BRCA1 and BRCA2 in Brazil. Sci Rep. 2018;8(1):9188. https://doi.org/10.1038/s41598-018-27315-2
» https://doi.org/10.1038/s41598-018-27315-2 -
29 Frank TS, Deffenbaugh AM, Reid JE, Hulick M, Ward BE, Lingenfelter B, et al. Clinical characteristics of individuals with germline mutations in BRCA1 and BRCA2: analysis of 10,000 individuals. J Clin Oncol. 2002;20(6):1480-90. https://doi.org/10.1200/JCO.2002.20.6.1480
» https://doi.org/10.1200/JCO.2002.20.6.1480 -
30 Guindalini RSC, Viana D, Kitajima JP, Valim A, Schlesinger D, Kok F, et al. Detection of inherited mutations in Brazilian breast cancer patients using multi-gene panel testing. In: 2018 American Society of Clinical Oncology Annual Meeting, 2018, Chicago. J Clin Oncol. 2018;36(15_suppl):e13610. https://doi.org/10.1200/JCO.2018.36.15_suppl.e13610
» https://doi.org/10.1200/JCO.2018.36.15_suppl.e13610 -
31 Carraro DM, Koike Folgueira MA, Garcia Lisboa BC, Ribeiro Olivieri EH, Vitorino Krepischi AC, de Carvalho AF, et al. Comprehensive analysis of BRCA1, BRCA2 and TP53 germline mutation and tumor characterization: a portrait of early-onset breast cancer in Brazil. PLoS One. 2013;8(3):e57581. https://doi.org/10.1371/journal.pone.0057581
» https://doi.org/10.1371/journal.pone.0057581 -
32 Fernandes GC, Michelli RA, Galvão HC, Paula AE, Pereira R, Andrade CE, et al. Prevalence of BRCA1/BRCA2 mutations in a Brazilian population sample at-risk for hereditary breast cancer and characterization of its genetic ancestry. Oncotarget. 2016;7(49):80465-81. https://doi.org/10.18632/oncotarget.12610
» https://doi.org/10.18632/oncotarget.12610 -
33 Ramos MCA, Folgueira MAAK, Maistro S, Campolina AG, Soárez PC, Bock GH, et al. Cost effectiveness of the cancer prevention program for carriers of the BRCA1/2 mutation. Rev Saude Publica. 2018;52:94. https://doi.org/10.11606/S1518-8787.2018052000643
» https://doi.org/10.11606/S1518-8787.2018052000643 - 34 Manchanda R, Sun L, Patel S, Evans O, Wilschut J, De Freitas Lopes AC, et al. Economic Evaluation of Population-Based BRCA1/BRCA2 Mutation Testing across Multiple Countries and Health Systems. Cancers (Basel). 2020;12(7):1929.
-
35 Momozawa Y, Iwasaki Y, Parsons MT, Kamatani Y, Takahashi A, Tamura C, et al. Germline pathogenic variants of 11 breast cancer genes in 7,051 Japanese patients and 11,241 controls. Nat Commun. 2018;9(1):4083. https://doi.org/10.1038/s41467-018-06581-8
» https://doi.org/10.1038/s41467-018-06581-8 -
36 Sun J, Meng H, Yao L, Lv M, Bai J, Zhang J, et al. Germline Mutations in Cancer Susceptibility Genes in a Large Series of Unselected Breast Cancer Patients. Clin Cancer Res. 2017;23(20):6113-9. https://doi.org/10.1158/1078-0432.CCR-16-3227
» https://doi.org/10.1158/1078-0432.CCR-16-3227 -
37 Rebbeck TR, Friebel TM, Friedman E, Hamann U, Huo D, Kwong A, et al. Mutational spectrum in a worldwide study of 29,700 families with BRCA1 or BRCA2 mutations. Hum Mutat. 2018;39(5):593-620. https://doi.org/10.1002/humu.23406
» https://doi.org/10.1002/humu.23406 -
38 Pritchard CC, Mateo J, Walsh MF, De Sarkar N, Abida W, Beltran H, et al. Inherited DNA-Repair Gene Mutations in Men with Metastatic Prostate Cancer. N Engl J Med. 2016;375(5):443-53. https://doi.org/10.1056/NEJMoa1603144
» https://doi.org/10.1056/NEJMoa1603144 -
39 Gymnopoulos M, Elsliger MA, Vogt PK. Rare cancer-specific mutations in PIK3CA show gain of function. Proc Natl Acad Sci U S A. 2007;104(13):5569-74. https://doi.org/10.1073/pnas.0701005104
» https://doi.org/10.1073/pnas.0701005104 -
40 Vasan N, Razavi P, Johnson JL, Shao H, Shah H, Antoine A, et al. Double PIK3CA mutations in cis increase oncogenicity and sensitivity to PI3Kα inhibitors. Science. 2019;366(6466):714-23. https://doi.org/10.1126/science.aaw9032
» https://doi.org/10.1126/science.aaw9032 -
41 Liang X, Lau QC, Salto-Tellez M, Putti TC, Loh M, Sukumar S. Mutational hotspot in exon 20 of PIK3CA in breast cancer among Singapore Chinese. Cancer Biol Ther. 2006;5(5):544-8. https://doi.org/10.4161/cbt.5.5.2656
» https://doi.org/10.4161/cbt.5.5.2656 -
42 Castaneda CA, Lopez-Ilasaca M, Pinto JA, Chirinos-Arias M, Doimi F, Neciosup SP, et al. PIK3CA mutations in Peruvian patients with HER2-amplified and triple negative non-metastatic breast cancers. Hematol Oncol Stem Cell Ther. 2014;7(4):142-8. https://doi.org/10.1016/j.hemonc.2014.09.007
» https://doi.org/10.1016/j.hemonc.2014.09.007 -
43 Desriani, Al-Ahwani F. The sensitivity and efficacy method of PIK3CA exon 9 E545A as a high diagnostic accuracy in breast cancer. J Genet Eng Biotechnol. 2018;16(1):71-6. https://doi.org/10.1016/j.jgeb.2017.10.002
» https://doi.org/10.1016/j.jgeb.2017.10.002 -
44 Samuels Y, Velculescu VE. Oncogenic mutations of PIK3CA in human cancers. Cell Cycle. 2004;3(10):1221-4. https://doi.org/10.4161/cc.3.10.1164
» https://doi.org/10.4161/cc.3.10.1164
APPENDIX
ADDITIONAL METHODSNGS
BRCA1 and BRCA2 were analyzed for mutations using the Ion AmpliSeq™ BRCA1 and BRCA2 panels (Thermo Fisher Scientific). This panel consists of three primer pools (167 amplicons) covering the entire coding region, including 10-20 bp of non-coding sequences flanking the 5’ and 3’ ends of each exon. Library preparation was performed using the Ion AmpliSeq™ Library Kit 2.0 and Ion Xpress™ Barcode Adapter 1-96 kit. DNA amplification was performed using 30 ng of DNA with three primer pools and 5x Ion AmpliSeq™ HiFi Master Mix. The PCR cycle included the following: 2 min at 99°C, followed by 19 cycles of 99°C for 15s and 60°C for 4 min, ending with a hold step at 10°C on a Veriti Thermal Cycler (Thermo Fisher Scientific). Next, the three PCR amplicons were mixed (30 µL), and 20 µL was treated with 2 µL FuPa Reagent to partially digest the primer sequences and phosphorylate the amplicons at 50°C for 10 min, followed by 55°C for 10 min, then 60°C for 20 min, and then held at 10°C. Next, sequencing adaptors (A: conjugated to biotin and P1) and barcodes (consisting of short stretches of index sequences that enable sample multiplexing) were ligated to the amplicons using the Ion Xpress™ Barcode Adapters kit (Thermo Fisher Scientific) for 30 min at 22°C, 5 min at 68°C, and 5 min at 72°C, ending with a hold at 10°C. The adaptor-ligated amplicon (libraries) were purified with 45 μL of the Agencourt® AMPure® XP Reagents (Beckman Coulter) and incubated for 5 min at room temperature(22-25°C). The tube was placed in a magnetic rack such that it was incubated for 2 min or until the solution became clear. After the supernatant was removed carefully and discarded without disturbing the pellet, freshly prepared 70% ethanol (150 μL) was added, and the tube was moved side‐to‐side of the magnet to wash the beads, and then the supernatant was discarded (two rounds of purification were repeated). The tube was placed on the magnet, and the beads were air-dried at room temperature for 5 min. The library was subjected to a second round of amplification using 50 μL of Platinum® PCR SuperMix HiFi and 2 μL of Equalizer™ Primers (added to each bead-pellet); the PCR cycles included 98°C for 2 min, followed by 9 cycles of 98°C for 15s and 64°C for 1 min, ending with a hold at 10°C. Then, 10 μL of Equalizer™ Capture was added to each amplified library, mixed by pipetting, and incubated at room temperature for 5 min. Next, 6 μL of washed Equalizer™ beads was added to each tube containing the captured library, mixed, and incubated at room temperature for 5 min. The tube was then placed in the magnet and incubated for 2 min or until the solution became clear. After the supernatant was removed carefully without disturbing the pellet, the Equalizer™ Wash Buffer (150 μL) was added to each reaction, to wash the beads, the tube was moved side‐to‐side of the magnet, and then the supernatant was removed and discarded (two rounds of purification were repeated). Next, the tube was removed from the magnet and 100 μL of Equalizer™ Elution Buffer was added to each pellet, mixed, and incubated at 32°C in a thermal cycler for 5 min. The tube was placed in the magnet and incubated at room temperature for 5 min or until the solution became clear. The supernatant contained the equalized library at ∼100 pM, and the same amount of the 12 libraries was pooled to perform the emulsion PCR. Next, emulsion PCR was performed using the Ion OneTouch™ System and Ion OneTouch™ 200 Template Kit v2 (Thermo Fisher Scientific). Template-positive Ion Sphere™ Particles (ISPs) were enriched using Dynabeads MyOne™ Streptavidin C1 beads (Invitrogen) and were washed with Ion OneTouch Wash Solution. This process was performed on an Ion OneTouch™ ES System (Thermo Fisher Scientific). The quality of the ISPs was evaluated using a Qubit 2.0 Fluorometer (Invitrogen). The enriched ISPs were sequenced on a 314 v2 Ion Chip (12 samples per chip) using an Ion Torrent Personal Genome Machine (PGM) sequencer system (Thermo Fisher Scientific) using the Ion PGM Sequencing 200 Kit version 2 (Thermo Fisher Scientific). Sequencing was performed using 500 flow runs, which generated approximately 200 bp. The PGM sequencing run outputs were directly loaded to the Torrent Server and stored as ‘.dat’ files. Data analysis comprising annotation of single-nucleotide variants, insertions, deletions, and splice-site alterations was performed using the Ion Reporter™ Server System (Life Technologies). Sequence data were also visually examined and verified using the Integrative Genomics Viewer (IGV). Sequencing generated an average of 302,346 reads per patient, and 96.14% of these regions were mapped to the BRCA1 and BRCA2 loci. Amplicons with coverage less than 30x on the Ion Torrent™ platform, as well as pathogenic variants, and new variants were reanalyzed by Sanger sequencing.
PCR amplification and Sanger sequencing
The complete coding regions of BRCA1 (NM_7294.3) and BRCA2 (NM_000059.3), including 50-100 base pairs (bp) of non-coding sequences flanking the 5’ and 3’ ends of each exon, were amplified by PCR using 33 pairs of primers (for BRCA1) (1-2) (Table 1), and 48 pairs of primers (for BRCA2) (3) (Table 2). Exons 9 (helical domain) and 20 (kinase domain) of PIK3CA (NM_006218.2) were amplified by PCR (Table 3). The PCR products were analyzed by Sanger sequencing in both the forward and reverse directions.
The reaction mixture (total volume, 20 μL) contained AmpliTaq Gold enzyme 250 U (final concentration, 0.04 U/μL) (Applied Biosystems, Foster City, CA, USA), 1× AmpliTaq Gold buffer, 1.5-3.0 mM AmpliTaq Gold magnesium chloride, 0.16 mM deoxynucleotides (Invitrogen, Carlsbad, CA, USA-AM8200), primers (0.4 μM each pair), and 50 ng DNA. PCR was performed on a Veriti® 96-well Thermal Cycler (Applied Biosystems™). The PCR cycle consisted of 1 cycle at 95°C for 10 min, 40 cycles at 94°C for 50 s, 54-66°C for 50s, 72°C for 50s, and 1 cycle at 72°C for 7 min. The BRCA2 exon 11 fragments were amplified by touchdown PCR, with annealing temperatures decreasing from 63°C to 56°C for fragments corresponding to the beginning of this exon to nucleotide 4526, and annealing temperatures from 68°C to 61°C for fragments corresponding to the end of the exon. PCR products were loaded onto a 1.5% agarose gel, stained with GelRed Nucleic Acid Stain (Biotium, Hayward, CA, USA), and evaluated. PCR products were treated with Illustra™ ExoStar™ 1-Step (GE Healthcare Bio-Sciences, Pittsburgh, PA, USA) and incubated at 37°C for 15 min, followed by incubation at 80°C for 15 min. All PCR products were sequenced in both forward and reverse directions using BigDye® Terminator v3.1 (Applied Biosystems, Foster City, CA, USA-4337456), according to the manufacturer’s instructions. The final product was sequenced on a 3500 Genetic Analyzer (Applied Biosystems™) or ABI 3730 DNA Analyzer (Applied Biosystems™). Sequences obtained were analyzed using Mutation Surveyor DNA Variant Analysis Software (v3.30, SoftGenetics LLC, State College, PA, USA). All pathogenic mutations were confirmed by Sanger sequencing.
MLPA
All patients were investigated for large rearrangements, specifically deletions and duplications, using the MLPA commercial kits, SALSA® MLPA® P002 BRCA1 probemix (P002-100R) and SALSA® MLPA® P045 BRCA2/CHEK2 probemix (P045-100R) (MRC-Holland, Amsterdam, The Netherlands). At first, 80 ng of genomic DNA resuspended in 2.5 μL ultrapure water was denatured for 10 min at 98°C after which 1.5 μL of the probemix mixture was added (0.75 μL of MLPA probe and 0.75 μL of MLPA buffer). The sample DNA and probemix mixture were heated at 95°C for 1 min and then incubated overnight at 60°C (17h). Afterward, ligation was performed using 1.5 μL of ligase buffer A, 1.5 μL of ligase buffer B, 0.5 μL of Ligase-65, and 12.5 μL of water and the reaction was incubated at 54°C for 15 min. The ligase was then inactivated by incubating the reaction at 98°C for 5 min. Amplification was performed by adding 5 μL of Polymerase mix (1 μL of SALSA PCR primers, 0.25 μL of SALSA polymerase, and 3.75 μL of water) and heated at 95°C for 1 min. PCR was carried out for 35 cycles (30s at 95°C, 30s at 60°C, and 60s at 72°C), followed by 20 min at 72°C on a GeneAmp 9700 Thermal Cycler (Applied Biosystems). Then, 1 μL of the PCR product was diluted 1:10 in water, mixed with 0.075 μL GeneScan™ 600 LIZ® dye Size Standard v2.0 (Applied Biosystems-4408399) and 9 μL Hi-Di Formamide (Applied Biosystems-4440753), and incubated at 80°C for 2 min on a GeneAmp 9700 Thermal Cycler (Applied Biosystems). The fragments were analyzed on an Applied Biosystems 3500 Genetic Analyzer (Applied Biosystems™) or ABI 3730 DNA Analyzer (Applied Biosystems™), and analysis was performed using Coffalyser.Net MLPA Analysis Software (MRC-Holland, Amsterdam, Netherlands). To normalize the data, at least three genomic DNA samples obtained from the peripheral blood cells of healthy donors were used as controls in each analysis. Normal values were considered when the ratio was between 0.8 and 1.2.
Nomenclature and classification of mutations
Variants were named according to the Human Genome Variation Society (HGVS) nomenclature (4). BRCA1 and BRCA2 variants were searched in publicly accessible databases, i.e., BRCA Share™ (5,6), BRCA Exchange (7), BRCA Mutation Database (8), and ClinVar (9); this search was performed between January and June 2020. Gene variants were evaluated using the following in silico prediction models: Polymorphism Phenotyping (PolyPhen; v2.2.2) (10), Sorting Intolerant From Tolerant (SIFT; v1.0.3) (11), Align-GVGD (12,13), Protein Variation Effect Analyzer (Provean; v1.1) (14), and Human Splicing Finder (15) to identify variants of unknown clinical significance. Minor allele frequency (MAF) was checked using the 1000 Genomes Project database (16), the Exome Aggregation Consortium (ExAC) (17,18), Global MAF dbSNP (19), Exome Variant Server, NHLBI GO Exome Sequencing Project (ESP) (20), Genome Aggregation Database (gnomAD) (21), Trans-Omics for Precision Medicine (TOPMed) (22), and Brazilian genomic variants (ABraOM) (23).
The variants were then classified according to the recommendations of the American College of Medical Genetics and Genomics in pathogenic, likely pathogenic, benign, likely benign, and variant of uncertain significance (VUS) (24). VUS for BRCA was also checked for co-occurrence with known pathogenic mutations in the same patient. For some variants, we considered that consensus information in ≥2 databases was strong enough to classify them as benign or VUS.
Clinical and pathological characteristics of breast cancer patients, BRCA sequencing, and the multiplex ligation-dependent probe amplification (MLPA) results.
In silico analysis of the alterations in exons 9 and 20 of PIK3CA in postmenopausal patients with breast cancer.
REFERENCES FOR ADDITIONAL MATERIAL
-
1 Arnold N, Gross E, Schwarz-Boeger U, Pfisterer J, Jonat W, Kiechle M. A highly sensitive, fast, and economical technique for mutation analysis in hereditary breast and ovarian cancers. Hum Mutat. 1999;14(4):333-9. https://doi.org/10.1002/(SICI)1098-1004(199910)14:4<333::AID-HUMU9>3.0.CO;2-C
» https://doi.org/10.1002/(SICI)1098-1004(199910)14:4<333::AID-HUMU9>3.0.CO;2-C -
2 Friedman LS, Ostermeyer EA, Szabo CI, Dowd P, Lynch ED, Rowell SE, et al. Confirmation of BRCA1 by analysis of germline mutations linked to breast and ovarian cancer in ten families. Nat Genet. 1994;8(4):399-404. https://doi.org/10.1038/ng1294-399
» https://doi.org/10.1038/ng1294-399 -
3 Wagner TM, Hirtenlehner K, Shen P, Moeslinger R, Muhr D, Fleischmann E, et al. Global sequence diversity of BRCA2: analysis of 71 breast cancer families and 95 control individuals of worldwide populations. Hum Mol Genet. 1999;8(3):413-23. https://doi.org/10.1093/hmg/8.3.413
» https://doi.org/10.1093/hmg/8.3.413 -
4 den Dunnen JT, Dalgleish R, Maglott DR, Hart RK, Greenblatt MS, McGowan-Jordan J, et al. HGVS Recommendations for the Description of Sequence Variants: 2016 Update. Hum Mutat. 2016;37(6):564-9. https://doi.org/10.1002/humu.22981
» https://doi.org/10.1002/humu.22981 -
5 Béroud C, Collod-Béroud G, Boileau C, Soussi T, Junien C. UMD (Universal mutation database): a generic software to build and analyze locus-specific databases. Hum Mutat. 2000;15(1):86-94. https://doi.org/10.1002/(SICI)1098-1004(200001)15:1<86::AID-HUMU16>3.0.CO;2-4
» https://doi.org/10.1002/(SICI)1098-1004(200001)15:1<86::AID-HUMU16>3.0.CO;2-4 -
6 Caputo S, Benboudjema L, Sinilnikova O, Rouleau E, Béroud C, Lidereau R, et al. Description and analysis of genetic variants in French hereditary breast and ovarian cancer families recorded in the UMD-BRCA1/BRCA2 databases. Nucleic Acids Res. 2012;40(Database issue):D992-1002. https://doi.org/10.1093/nar/gkr1160
» https://doi.org/10.1093/nar/gkr1160 -
7 BRCA Exchange. Available from: https://brcaexchange.org/ [Accessed June 28th, 2020]
» https://brcaexchange.org/ -
8 Arup Laboratories. Available from: https://arup.utah.edu/database/BRCA/Home/BRCA1_landing.php [Accessed June 28th, 2020]
» https://arup.utah.edu/database/BRCA/Home/BRCA1_landing.php -
9 Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, et al. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42(Database issue):D980-5. https://doi.org/10.1093/nar/gkt1113
» https://doi.org/10.1093/nar/gkt1113 -
10 Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248-9. https://doi.org/10.1038/nmeth0410-248
» https://doi.org/10.1038/nmeth0410-248 -
11 Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11(5):863-74. https://doi.org/10.1101/gr.176601
» https://doi.org/10.1101/gr.176601 -
12 Tavtigian SV, Deffenbaugh AM, Yin L, Judkins T, Scholl T, Samollow PB, et al. Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral. J Med Genet. 2006;43(4):295-305. https://doi.org/10.1136/jmg.2005.033878
» https://doi.org/10.1136/jmg.2005.033878 -
13 Mathe E, Olivier M, Kato S, Ishioka C, Hainaut P, Tavtigian SV. Computational approaches for predicting the biological effect of p53 missense mutations: a comparison of three sequence analysis based methods. Nucleic Acids Res. 2006;34(5):1317-25. https://doi.org/10.1093/nar/gkj518
» https://doi.org/10.1093/nar/gkj518 -
14 Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7(10):e46688. https://doi.org/10.1371/journal.pone.0046688
» https://doi.org/10.1371/journal.pone.0046688 -
15 Desmet FO, Hamroun D, Lalande M, Collod-Béroud G, Claustres M, Béroud C. Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009;37(9):e67. https://doi.org/10.1093/nar/gkp215
» https://doi.org/10.1093/nar/gkp215 -
16 1000 Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56-65. https://doi.org/10.1038/nature11632
» https://doi.org/10.1038/nature11632 -
17 http://exac.broadinstitute.org/variant/9-139413097-T-G [Accessed June 28th, 2020]
» http://exac.broadinstitute.org/variant/9-139413097-T-G -
18 Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285-91. https://doi.org/10.1038/nature19057
» https://doi.org/10.1038/nature19057 -
19 https://pubchem.ncbi.nlm.nih.gov/
» https://pubchem.ncbi.nlm.nih.gov/ -
20 http://evs.gs.washington.edu/EVS/
» http://evs.gs.washington.edu/EVS/ -
21 Genome Aggregation Database (gnomAD). Availabel from: https://gnomad.broadinstitute.org/
» https://gnomad.broadinstitute.org/ - 22 Trans-Omics for Precision Medicine(TOPMed)
-
23 Naslavsky MS, Yamamoto GL, de Almeida TF, Ezquina SAM, Sunaga DY, Pho N, et al. Exomic variants of an elderly cohort of Brazilians in the ABraOM database. Hum Mutat. 2017;38(7):751-63. https://doi.org/10.1002/humu.23220
» https://doi.org/10.1002/humu.23220 -
24 Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405-24. https://doi.org/10.1038/gim.2015.30
» https://doi.org/10.1038/gim.2015.30
Publication Dates
-
Publication in this collection
16 July 2021 -
Date of issue
2021
History
-
Received
3 Feb 2021 -
Accepted
21 May 2021