Research Papers:

Determination of genetic predisposition to early breast cancer in women of Kazakh ethnicity

Metrics: PDF 191 views  |   Full Text 1539 views  |   ?  

Gulnur Zhunussova _, Nazgul Omarbayeva _, Dilyara Kaidarova, Saltanat Abdikerim, Natalya Mit, Ilya Kisselev, Kanagat Yergali, Aigul Zhunussova, Tatyana Goncharova, Aliya Abdrakhmanova and Leyla Djansugurova


Gulnur Zhunussova1,3, Nazgul Omarbayeva2,4, Dilyara Kaidarova2,4, Saltanat Abdikerim1,3, Natalya Mit1, Ilya Kisselev1, Kanagat Yergali1, Aigul Zhunussova1, Tatyana Goncharova2, Aliya Abdrakhmanova2,4 and Leyla Djansugurova1,3

1 Laboratory of Molecular Genetics, Institute of Genetics and Physiology, Almaty 050060, Kazakhstan

2 Kazakh Institute of Oncology and Radiology, Almaty 050060, Kazakhstan

3 Al-Farabi Kazakh National University, Almaty 050060, Kazakhstan

4 Asfendiyarov Kazakh National Medical University, Almaty 050060, Kazakhstan

Correspondence to:

Gulnur Zhunussova, email: [email protected]
Nazgul Omarbayeva, email: [email protected]

Keywords: early-onset breast cancer; triple negative breast cancer; next-generation sequencing; pathogenic variant; Kazakh population

Received: May 06, 2023     Accepted: September 21, 2023     Published: October 04, 2023

Copyright: © 2023 Zhunussova et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Breast cancer (BC) is the most common type of cancer among women in Kazakhstan. To date, little data are available on the spectrum of genetic variation in Kazakh women with BC. We aimed to identify population-specific genetic markers associated with the risk of developing early-onset BC and test their association with clinical and prognostic factors. The study included 224 Kazakh women diagnosed with BC (≤40 age). Entire coding regions (>1700 exons) and the flanking noncoding regions of 94 cancer-associated genes were sequenced from blood DNA using MiSeq platform. We identified 38 unique pathogenic variants (PVs) in 13 different cancer-predisposing genes among 57 patients (25.4%), of which 6 variants were novel. In total, 12 of the 38 distinct PVs were detected recurrently, including BRCA1 c.5266dup, c.5278-2del, and c.2T>C, and BRCA2 c.9409dup and c.9253del that may be founder in this population. BRCA1 carriers were significantly more likely to develop triple-negative BC (OR = 6.61, 95% CI 2.44–17.91, p = 0.0002) and have family history of BC (OR = 3.17, 95% CI 1.14–8.76, p = 0.03) compared to non-carriers. This study allowed the identification of PVs specific to early-onset BC, which may be used as a foundation to develop regional expertise and diagnostic tools for early detection of BC in young Kazakh women.


In Kazakhstan, breast cancer (BC) is the most common malignancy in women. According to the latest national study on the epidemiology of BC in Kazakhstan, the number of reported cases increased by 14% between 2017 and 2021 [1]. By nationality, most BC patients with a confirmed primary diagnosis are Kazakh (48.1%), followed by Russian (33.1%) [1]. A detailed study within ethnic groups shows that BC accounts for 26.3% of all cancer cases and 8.7% of cancer-related deaths in Kazakh women [2]. At the same time, although the prevalence of BC in Kazakhstan shows a steady upward trend in the age range from 42.5 years to 62.5 years, the prevalence of BC in young women has also increased from 309 cases in 2017 to 330 cases in 2021 [1]. Thus, combined with a poorer prognosis, more aggressive histologic features, and more frequent recurrences [3, 4], BC in young Kazakh women in Kazakhstan is a growing threat. The problem is exacerbated by the fact that unlike an effective golden standard of mammography screening for women aged 40–70 years, screening programs for BC in young women are still underdeveloped in Kazakhstan. Targeted development of molecular genetic biomarkers for genetic predisposition to early BC in Kazakh women could provide a useful choice of efficient and noninvasive diagnostic methods.

Genetic factors significantly contribute to the risk of early-onset BC. In the United States, approximately 5 to 10% of BC cases are believed to have a hereditary risk [5]. However, these percentages are higher in a subset of women with early-onset BC, familial history, or triple-negative BC (TNBC), i.e., estrogen and progesterone receptors (ER and PR, respectively) and human epidermal growth factor receptor 2 (HER2) negative, based on immunohistochemistry (IHC) [6, 7]. Moreover, the cumulative risk of developing BC by age 50 in sisters and mothers of patients under 40 years old with BC was 6 and 3 times higher than the average population, respectively [7].

BC-associated variants have variable penetrance. Rare, highly penetrant genes that confer a >5-fold relative risk [8] tend to cluster in families and are associated with hereditary BC subtypes. Moderate-penetrance variants that confer two- to five-fold increased risks [8] also contribute to a minority of BC cases. According to epidemiological studies, 15–20% of familial BC cases and less than 5% of all BCs are caused by germline variants in the BRCA1 and BRCA2 genes [9]. Carriers of pathogenic or likely pathogenic variants (PVs) in these genes have a 60–85% risk of developing BC by the age of 70 [10]. Results from population studies have shown that PVs in BRCA1 and BRCA2 are responsible for the development of BC in only 5–10% of cases [11]. Thus, other high-to-moderate-penetrance variants might also contribute to the onset of BC in the remaining cases. These include the rare highly penetrant variants in the CDH1, PTEN, STK11, TP53 genes, and the moderately penetrant variants in the ATM, CHEK2, PALB2 genes, among others. Genetic variants in these genes are associated with a 2–5-fold increase in the risk of developing BC [1214]. Currently, according to the guidelines of the National Comprehensive Cancer Network (NCCN), preventive diagnostic and therapeutic measures are recommended for carriers of pathogenic variants in the ATM, BARD1, BRCA1/BRCA2, CDH1, CHEK2, NF1, PALB2, PTEN, RAD51C/RAD51D, STK11, and TP53 genes.

The existing differences in genetic structure and history between different populations require identifying and analyzing the spectrum and prevalence of PVs in genes involved in carcinogenesis in specific ethnic groups. To date, little data are available on the hereditary cancer risk among Kazakh women with BC. To characterize the inherited BC risk among Kazakh women with early onset BC, we used next-generation sequencing (NGS) to identify PVs involved in early-onset BC, and assessed their relationship with clinical and prognostic factors. Our study may reveal previously uncharacterized population-specific variants that may increase the risk of BC in the Kazakh population.


Clinical characteristics of the study population

Our study enrolled 224 unrelated Kazakh women diagnosed with early onset BC. The median age at diagnosis was 34.6 with a range from 19 to 40 years. Of these, 35 patients (15.6%) were diagnosed with BC under the age of 30 and 31 patients (13.8%) had family history of BC. Pathological examination showed predominantly invasive ductal carcinoma (90.6%) and moderately (64.7%) and poorly differentiated grade (30.4%) of the tumor. At the time of diagnosis, most patients (71.4%) had stage II of the disease and 14.6% had metastases in the axillary lymph node. The majority of cases (88.4%) had high (≥14%) Ki-67 expression, indicating rapid proliferation of cancer cells [15]. Immunohistochemistry (IHC) revealed that most patients had the aggressive subtypes of BC: TNBC (32.6%) and Luminal B (Her2 negative) (27.3%). The detailed clinicopathological characteristics of patients are summarized in Table 1.

Table 1: Characteristics of the study cohort and comparison of carriers with non-carriers of a pathogenic variants

Breast cancer (BC)All patients
No (%)
carriers No (%)
carriers No (%)
No (%)
No of patients22438 (17.0)19 (8.5)167 (74.6)
Age at diagnosis, years
 Median Age (Min-Max)34.6 (19–40)34.9 (27–40)34.7 (23–40)34.7 (19–4)0.110
 Age 19–3043 (19.20)9 (23.7)7 (36.8)27 (16.2)0.079
 Age 31–40181 (80.80)29 (76.3)12 (63.1)140 (83.8)
Family history of cancer
 No family history of cancer178 (79.5)26 (68.4)19 (100.0)133 (79.6)0.041
 Breast cancer31 (13.8)10 (26.3)0 (0.0)21 (12.6)
 Other cancers15 (6.7)2 (5.3)0 (0.0)13 (7.8)
 Invasive ductal carcinoma203 (90.6)35 (92.1)14 (73.7)154 (92.2)0.034
 Invasive lobular carcinoma12 (5.4)3 (7.9)3 (15.8)6 (3.6)
 Others9 (4.0)0 (0.0)2 (9.1)7 (4.2)
 Well differentiated8 (3.6)1 (2.6)0 (0.0)7 (4.2)0.821
 Moderately differentiated145 (64.7)24 (63.2)11 (57.9)110 (65.9)
 Poorly differentiated68 (30.4)12 (31.6)8 (42.1)48 (28.7)
 Not applicable3 (1.4)1 (2.7)0 (0.0)2 (1.2)
In situ2 (0.9)0 (0.0)0 (0.0)2 (1.2)0.09
 I15 (6.7)2 (5.3)1 (5.3)12 (7.2)
 II160 (71.4)29 (76.3)9 (47.4)122 (73.1)
 III38 (17.0)6 (15.8)6 (31.6)26 (15.6)
 IV9 (4.0)1 (2.6)3 (15.8)5 (3.0)
Axillary lymph node metastases (stages II–III)
 Yes29 (14.6)5 (14.3)6 (40.0)18 (12.2)0.023
 No169 (85.4)30 (85.7)9 (60.0)130 (87.8)
 High (≥14)198 (88.4)34 (89.5)18 (94.7)146 (87.4)0.872
 Low (<14)24 (10.7)4 (10.5)1 (5.3)19 (11.4)
 Not applicable2 (0.9)0 (0.0)0 (0.0)2 (1.2)
Molecular subtype
 Luminal A (ER+ and/or PR+ HER2− KI67−)27 (12.1)3 (7.9)2 (10.5)22 (13.2)0.005
 Luminal B− (ER+ and/or PR+ HER2− KI67+)61 (27.3)9 (23.7)7 (36.8)45 (27.0)
 Luminal B+ (ER+ and/or PR+ HER2+ KI67+)38 (17.1)6 (15.8)1 (5.3)31 (18.6)
 HER2 overexpression (ER− and/or PR− HER2+)23 (10.3)1 (2.6)3 (15.8)19 (11.4)
 TNBC (ER− and/or PR− HER2−)73 (32.6)19 (50.0)6 (31.6)48 (28.7)
 Not applicable2 (0.9)0 (0.0)0 (0.0)2 (1.2)

Spectrum of germline PVs

A total of 23,013 unique variants in 87 genes were detected across 224 unrelated patients with early-onset BC. The identified variants were 7,720 missense, 12,727 synonymous, 1,230 3′-UTR, 668 5′-UTR, 31 frameshift, 16 nonsense, 10 in-frame, 9 splice site, 3 start-lost, and 599 intron or non-coding. Of these, 1,093 variants were found at frequency ≤1%. The spectrum of rare variants in the 224 patients is represented in Figure 1. To filter out rare polymorphisms (1093 variants, Figure 1), we eliminated variants that occur at a frequency of more than 1% in population databases 1000 Genomes database (1000G) [16], the Exome Sequencing Project database (ESP6500, https://esp.gs.washington.edu), and the Exome Aggregation Consortium database (ExAC).

Figure 1: Flowchart illustrating the filtering and spectrum of rare variants detected by NGS.

In total, we identified 59 PVs (38 unique positions) located in the 13 genes, which occurred in 57 patients (25.4%) (Supplementary Table 1). The prevalence of PVs out of the total PVs identified in the 13 genes was as follows: BRCA1 (37.3%), BRCA2 (27.1%), TP53 (8.5%), CHEK2 (6.8%), PALB2 (5.1%), SDHB (3.4%), ATM (1.7%), BLM (1.7%), FANCM (1.7%), NBN (1.7%), PMS1 (1.7%), PMS2 (1.7%), and XPA (1.7%). Thus, the most common pathogenic germline variants in 57 young patients with PVs were the BRCA1 (39%, 22/57) gene, followed by the BRCA2 (28%, 16/57), TP53 (9%, 5/57), CHEK2 (7%, 4/57), PALB2 (5%, 3/57), and SDHB (4%, 2/57) genes. Other genes ATM, BLM, FANCM, NBN, PMS1, PMS2, and XPA were present in 14% of patients with PVs (2%, 11/57 for each gene).

Patients with multiple PVs were also found in our study. Specifically, two patients were carriers of two PVs in different cancer predisposition genes. One patient diagnosed with hormone-positive BC at the age of 31 years and a first-degree relative with BC (mother diagnosed with BC at the age of 45) had BRCA2 c.9409dup in combination with ATM c.2465T>G. The other patient diagnosed with TNBC at the age of 39 was carrier of BRCA1 c.2T>C and BLM c.320dup.

The 59 PVs identified in the 57 patients included 29 frameshift, 11 nonsense, 8 splice site, 7 missense, 3 start-loss, and 1 synonymous. As shown in Figure 2, frameshift variants were the predominant type in both the BRCA1/2 and non-BRCA groups. Spice site and start-loss variants were detected only in the BRCA1/2 group. In the non-BRCA group, nonsense and missense variants were evenly distributed. All the detected PVs existed in a heterozygous state.

Figure 2: The spectrum of pathogenic variants in the early-onset breast cancer cohort.

In addition to the 59 PVs, we also identified additional variants (5 in-frame, 2 nonsense, 1 frameshift and 1 splice site) in 16 patients that were not annotated neither as (likely) pathogenic nor as (likely) benign in LOVD, ClinVar, and ARUP’ databases. Thus, these were labeled as variants of unknown significance (VUS, Supplementary Table 2).

Identification of recurrent PVs

In total, 12 (31.6%) of the 38 unique PVs were seen recurrently. We identified five PVs that each occurred in ≥3 patients in our cohort (Figure 3). Of these, three variants were observed in the BRCA1 gene, and two were detected in the BRCA2 gene. The most common PVs were the frameshift insertion c.5266dup (also known as c.5382ins) and the splice acceptor site variant c.5278-2del (c.5341-2del) in BRCA1 gene, each accounting for 2.7% (n = 6) of the cohort (Figure 3). Two frameshift variants (c.9409dup and c.9253del) in BRCA2 gene and one nonsense variant in BRCA1 (c.2T>C) were each observed in three patients. Five additional PVs (BRCA2 c.2808_2811del, BRCA2 c.7567_7568del, BRCA1 c.2498del (novel), CHEK2 c.470T>C, and SDHB c.725G>A) each occurred in two patients, whereas 25 PVs (65.8%) occurred in single patients.

Figure 3: Lollipop charts of the genes with the highest number of variants in this study.

The variant type is shown by colored balls. Protein domains are indicated as colored squares. BRCT1/2, BRCA1/2 C-terminus. Abbreviations: TAD: transactivation domain; FHA: forkhead-associated domain.

Identification of novel PVs

Overall, of the 59 PVs, we identified 6 novel PVs (BRCA1 c.2498del; PALB2 c.1034T>G, c.18_22del; TP53 c.154C>T, XPA c.20del, PMS1 c.1258del, Supplementary Table 1) that were not previously reported in any of the databases (LOVD, ClinVar, and ARUP). Among these, the novel frameshift deletion BRCA1 c.2498del was found in two patients diagnosed with unilateral TNBC at young age (31 and 28 years) and with first-degree relatives with cancer.

Rare missense variants and variants of unknown significance

A total of 559 rare missense variants (minor allele frequency (MAF) <1%) in the population frequency databases) were detected in this study. Except for 1.3% variants reported as pathogenic in both ClinVar and LOVD, most of the variants were classified as VUS. All of them were in a heterozygous state. In ClinVar, 46.3% of the rare missense variants were registered as VUS, 23.6% were not reported, and 28.8% were considered to be benign. When we evaluated the impact of rare missense variants on protein structure and function using in silico algorithms (SIFT and PolyPhen-2) [17, 18], 291 variants in 161 patients were identified as potentially deleterious and probably damaging according to at least one prediction tool (Supplementary Table 3), while 53 variants in 46 patients were reported to be deleterious by both prediction tools (Supplementary Table 4). Six patients were carriers of at least two VUS each, all found in non-BRCA1/2 genes.

Association between germline PVs and clinical factors

To analyze the impact of the identified PVs on the clinical and pathological characteristics of early onset BC patients, we divided the cohort into BRCA1/2 carriers, non-BRCA1/2 carriers, and non-carriers (of PVs) (Table 1). The number of patients with at least one first-degree relative with BC was significantly higher among BRCA1/2 carriers (26.3%) compared to non-carriers (12.6%), while no family history of BC was reported among BRCA1/2 non-carriers (p = 0.041). While then number of patients with invasive ductal carcinoma was high in all three groups, a significantly higher proportion was observed in BRCA1/2 carriers (92.1%) and non-carriers (92.2%) compared to BRCA1/2 non-carriers (73.7%, p = 0.034). A significantly higher proportion of non-BRCA1/2 carriers (40.0%) had axillary lymph node metastases (stages II–III) compared with both BRCA1/2 carriers (14.3%) and non-carriers (12.2%, p = 0.023). Molecular subtype assessment indicated that 50% of BRCA1/2 carriers had TNBC, which was significantly higher than the frequency observed in non-BRCA1/2 carriers (31.6%) and non-carriers (28.7%, p = 0.005).

Biomarkers and prognostic factors

To investigate further the contribution of the identified PVs to clinical and pathological characteristics of patients with early onset BC, we divided the cohort into two groups (carriers and non-carriers), and tested their association with hormonal status (Table 2) and prognostic factors (Table 3). Since most of the PVs affect the BRCA1 and BRCA2 genes, we focused on associations with these genes.

Table 2: Association between BRCA1/2 pathogenic variants and hormonal status

Hormonal statusBRCA1BRCA2
Carriers (%)OR (95% CI)P-valueCarriers (%)OR (95% CI)P-value
ER status
 Negative17 (77.3)1.000.0043 (20)1.000.08
 Positive5 (22.7)0.22 (0.08–0.62)12 (80)2.96 (0.80–10.87)
PR status
 Negative17 (77.3)1.000.0065 (33.33)1.000.40
 Positive5 (22.7)0.23 (0.08–0.66)10 (66.67)1.59 (0.52–4.85)
Her2 status
 Negative21 (95.5)1.000.03610 (66.67)1.000.77
 Positive1 (4.5)0.11 (0.02–0.87)5 (33.33)1.20 (0.39–3.70)
 No6 (27.3)1.000.000212 (80)1.000.56
 Yes16 (72.7)6.61 (2.44–17.91)3 (20)0.62 (0.16–2.29)

Table 3: Association of BRCA1/2 pathogenic variant carriers for prognostic factors compared to noncarriers

Clinical statusBRCA1
Carriers (%)
OR (95% CI)P-valueBRCA2
Carriers (%)
OR (95% CI)P-value
Bone metastasis
 Negative21 (95.5)153 (91.6)1.000.5414 (87.5)153 (91.6)1.000.58
 Positive1 (4.5)14 (8.4)0.52 (0.06–4.16)2 (12.5)14 (8.4)1.56 (0.32–7.57)
Lymph node metastasis
 No19 (86.4)147 (88.0)1.000.8213 (81.3)147 (88.0)1.000.44
 Yes3 (13.6)20 (12.0)1.16 (0.31–4.27)3 (18.7)20 (12.0)1.70 (0.44–6.47)
Lung metastasis
 No20 (90.9)152 (91.0)1.000.9914 (87.5)152 (91.0)1.000.64
 Yes2 (9.1)15 (9.0)1.01 (0.21–4.76)2 (12.5)15 (9.0)1.45 (0.30–6.98)
Liver metastasis
 No21 (95.5)159 (95.2)1.000.9614 (87.5)159 (95.2)1.000.21
 Yes1 (4.5)8 (4.8)0.95 (0.11–7.94)2 (12.5)8 (4.8)2.84 (0.54–14.68)
CNS metastasisv
 No21 (95.5)157 (94.0)1.000.7915 (93.7)157 (94.0)1.000.97
 Yes1 (4.5)10 (6.0)0.75 (0.09–6.14)1 (6.3)10 (6.0)1.05 (0.12–8.74)
Family history of cancer
 Negative14 (63.6)133 (79.6)1.0011 (73.3)133 (79.6)1.00
 Positive for BC7 (31.8)21 (12.6)3.17 (1.14–8.76)0.033 (20.0)21 (12.6)1.73 (0.44–6.71)0.43
 Positive for other cancers1 (4.6)13 (7.8)0.73 (0.09-6.01)0.771 (6.7)13 (7.8)0.93 (0.11–7.79)0.95
Histological grade
 I–II17 (77.3)134 (81.3)1.000.6614 (82.4)134 (81.3)1.000.91
 III–IV5 (22.7)31 (18.7)1.27 (0.43–3.71)3 (17.6)31 (18.7)0.93 (0.25–3.42)

As shown in Table 2, the association with TNBC in BRCA1 carriers was significant (OR = 6.61, 95% CI 2.44–17.91, p = 0.0002). We found a significant association between BRCA1 carriers and a family history of BC (OR = 3.17, 95% CI 1.14–8.76, p = 0.03) whereas BRCA2 carriers had worse prognostic outcomes of developing a metastatic phenotype but not statistically significant as shown in Table 3.


Pathogenic variants prevalence

In clinical practice, genetic testing for BC risks has been based chiefly on BRCA1/2 gene analysis, despite new evidence suggesting the clinical significance of a broader number of cancer-related genes [19]. In this study, we used a targeted panel including a vast number of genes implicated in hereditary cancer syndromes and overall BC cancer predisposition.

To our knowledge, this is the first study using NGS technology to study the genetic predisposition to early-onset BC women from Kazakhstan and assess their impact on the patients’ clinical outcomes. Overall, 25.4% of Kazakh women were carriers of one or more PVs in the 94 cancer-associated genes analyzed. PVs in the 12 established BC predisposition genes (ATM, BARD1, BRCA1, BRCA2, CDH1, CHEK2, NF1, PALB2, PTEN, RAD51C, RAD51D, and TP53) were detected in 22.3% of our patient population, which is higher than that reported in the CARRIERS consortium that was a population-based study of hereditary BC genes in the United States (5.03%) [20]. When compared to other studies of early-onset BC, the proportion of PVs found in our study was lower than reported in the Greek population (31.5%), but similar to other studies in diverse populations [21, 22].

The genetic landscape of the identified PVs was dominated by the BRCA1 gene, followed by the BRCA2, TP53, CHEK2, PALB2, and SDHB genes, and was consistent with the profile of early-onset BC and/or ovarian cancer in the Indian population and in BC patients in the Brazilian miscegenated population, where PVs in the BRCA1 gene were also dominant, followed by the BRCA2 and TP53 genes [23, 24]. In contrast, the most common pathogenic germline variants in the Greek patients with early-onset BC were PVs in CHEK2, followed by BRCA1/2 and TP53 genes [21]. However, compared to the Kazakh women, no PVs were detected in the PALB2, SDHB, ATM, BLM, FANCM, NBN, PMS1, PMS2, and XPA genes among the 94 cancer-associated genes analyzed in the Greek and Indian populations, which may be due to the smaller samples. Interestingly, when the same panel of 94 cancer genes was examined in the group of TNBC BRCA-negative patients from Cyprus, germline PVs of the PALB2 gene were found in older women but not in young [25]. At the same time, in a large-scale study of BC women aged 35–59 years, who were mainly of Northern and Western European descent, PVs in CHEK2, ATM, PALB2, and PMS2 genes were the most frequent among the 25 cancer genes tested, after PVs in BRCA1/2 [26].

BRCA1/2 pathogenic variants

Our finding suggests that germline PVs in the BRCA1 and BRCA2 genes are major contributors to early-onset BC in the young Kazakh female population. Among the 59 PVs detected in 57 patients (25.4% of patients), the majority of the PVs were found in BRCA genes (64.4% of total PVs or 67% of total PV patients), which reflects the high proportion of BRCA1/2 variants in early-onset BC women reported in other ethnic groups [21, 26]. In our study, the five PVs in BRCA1 identified in 15 patients may lead to defects in BRCT domains, while the five PVs in BRCA2 occurring in 7 patients were found in functional domains that may affect interactions with RAD51, NPM1, and SEM1 (Figure 3). It is known that PVs in these genes are most prevalent in patients with hereditary breast and/or ovarian cancer [21, 2729]. However, some studies using a targeted sequencing approach have discovered a higher proportion (62–73%) of PVs in other cancer susceptibility genes other than BRCA1/2 [23, 30], with nearly equal contribution to hereditary BC [26, 31].

Non-BRCA1/2 pathogenic variants

In our study, 35.6% of all PVs found in non-BRCA genes might also contribute to early onset BC. Overall, 21 PVs in 11 non-BRCA genes (ATM, BLM, CHEK2, FANCM, NBN, PALB2, PMS1, PMS2, SDHB, TP53, and XPA) were found in 21 patients. Two patients had both BRCA and non-BRCA mutated genes (BRCA2 c.9409dup plus ATM c.2465T>G, and BRCA1 c.2T>C plus BLM c.320dup), with both having family history of BC and poor clinical outcomes.

The most frequently mutated genes after BRCA1/2 genes were TP53 (5/224; 2.2%), followed by CHEK2 (4/224; 1.8%) and PALB2 (3/224; 1.3%) (Figure 3 and Supplementary Table 1). Among the studied cohort, we did not find any constitutional mosaicism in TP53 variants. Sixteen out of 21 PVs detected in this study were protein-truncating (the majority being frameshift), potentially resulting in loss of protein function.

Even though the frequency and spectrum of PVs in Kazakh women resembled those reported in Guindalini R.S.C et al. [24], most BC-associated genes and variants differ between races and populations [32, 33]. Work is evolving to better understand the implications of multiple pathogenic variants in an individual in regard to how it impacts cancer phenotype, age of cancer onset and tumor biology [34]. It remains to be determined whether these variants modify the BC risk. There is still limited data on tumor pathogenesis among patients with concurrent PVs in cancer genes, and future studies may provide valuable insight in counseling and management of individuals with multiple germline PVs.

Recurrent and founder pathogenic variants

A founder variant is a PV that is observed at high frequency in a given population due to the presence of the variant in a single ancestor or a small number of ancestors [35]. Previous studies have shown that different populations have founder genetic variants that affect BC risk [3639].

We detected 12 (31.6%) PVs that occurred in multiple patients and affected 13.8% (31/224) of the patient cohort. The most prevalent PVs were BRCA1 c.5278-2del and c.5266dup, each found in 2.7% of total case (Supplementary Table 1). The c.5278-2del found within the BRCT (BRCA1 C-terminal) domain affects an acceptor splice site in intron 19 of the BRCA1 gene, which disrupts RNA splicing by generating two aberrant transcripts (one lacking exon 21 and the other skipping 8 bp at the 5′-end of exon 21) that lead to the formation of a truncated protein [40]. The frameshift insertion c.5266dup in the last exon of the BRCA1 mRNA transcript is predicted to escape nonsense-mediated decay and be expressed as a truncated protein (p.Gln1756Profs*74) that lacks the C-terminal BRCT domain.

There is no evidence to date about the founder effect of BRCA1 c.5278-2del within different populations, whereas BRCA1 c.5266dup was reported as one of the most frequent PVs in Eastern and Central Europe [41] and North Africa [42]. It is also discussed as a founder mutation for the development of BC in Uzbek, Ashkenazi Jewish, Tatar and Russian populations [4346]. Previous studies of the BRCA1 gene showed the absence of the c.5266dup PV in a cohort of 121 Kazakh patients with BC, mean age of 50.3 ± 11.5 years [47]. A previous study using Sanger sequencing to identify BRCA variants in sporadic BC women from Kazakhstan with an average age of 50–51 years did not reveal any recurrent PVs, including BRCA1 c.5266dup and c.5278-2del [48]. It is worth noting the advantage of using NGS technology, which significantly expands the possibilities for detecting relevant variants, rather than the conventional genotyping methods. Since BRCA1 c.5266dup and c.5278-2del were found at high frequency among patients under 40 years of age (mean age 34 years) in our cohort, we can assume that they are characteristic of the early development of BC, although future studies investigating the frequency distribution of these variants in patients subgroups by age are needed for confirmation. BRCA1 c.5278-2del has been detected in one Italian patient with a personal and family history of BC and reported as a rare variant [49]. This variant is currently not reported in the 1000G, ESP6500 and ExAC databases. However, it was also previously identified in a colorectal cancer patient of Kazakh ethnicity [50]. Phelan et al. in their prospective cohort study showed four-fold increased risk of early-onset colorectal cancer in female BRCA1 PVs carriers [51]. Our findings suggest that BRCA1 c.5266dup and c.5278-2del may be founder PVs in the Kazakh population. This could have important implications for preventive screening and targeted treatments for BC in the Kazakhstan population.

Other common PVs identified in our study are BRCA2 c.9253del (also known as c.9481del) and c.9409dup, and BRCA1 c.2T>C, each accounting for 1.3% of the entire cohort. Both frameshift variants in exons 24 and 25 of the BRCA2 gene trigger the nonsense-mediated decay by introducing a premature termination codon (p.Thr3085GlnfsTer19 and p.Thr3137AsnfsTer13) in the DNA-binding OB-fold C-terminal domain (OB3). BRCA2 c.9409dup was reported in LOVD and ARUP, but not in ClinVar. This PV has also been identified at a frequency of 0.2% in a Chinese BC cohort, which is much lower than the frequency reported in our study [52]. BRCA2 c.9253del was also detected at a low frequency in Korean BC patients [5255]. A sequencing study on a cohort of 32 early-onset BC patients from multiple ethnic regions of China also found this variant in one woman of Kazakh ethnicity [56]. Thus, these results suggest that BRCA2 c.9409dup and c.9253del could be specific to Asian populations. The BRCA1 c.2T>C variant modifies the methionine residue at the initiation codon. This variant is not reported in any population frequency database, and no functional characterizations are available yet. However, in our study, c.2T>C occurred in three patients with TNBC, and a neighbor PV (c.1A>G) resulting in the loss of the initial 17 amino acids of BRCA1 was also frequently identified in breast and ovarian cancer families [57, 58].

Overall, the profile of founder mutations for the Kazakh population appears to be represented by BRCA1 c.5266dup, which is common in European, African, and Asian ethnic groups; BRCA2 c.9409dup and c.9253del, which are specific to Asian populations; and BRCA1 c.5278-2del and c.2T>C. The last two PVs could be ethnicity-specific founder mutations for Kazakhs compared to other Central Asian ethnicities, as they have not been discussed as founder PVs in these ethnicities. For example, the founder mutation profile in Tatar women with BC includes BRCA1 c.5161C>T and c.300T>G and BRCA2 c.468dup in addition to BRCA1 c.5266dup [45]. The most used test for Russian BC and ovarian cancer patients includes BRCA1 c.4153delA, c.185delAG and BRCA2 c.6174delT in addition to BRCA1 c.5266dup. At the same time, the inclusion of four additional recurrent BRCA1 mutations c.C61G, c.2080delA, c.3819del5, and c.3875del4 in the assay led to a significant increase in the number of identified Russian mutation carriers [46, 59]. The spectrum of mutations identified in Uighur women with early-onset and sporadic BC was completely different from that of Kazakhs and includes BRCA1 c.3180insA, c.3538insT, c.3694insAA, c.1963insT, c.3948G>C, c.3182A>C, and c.3538G>T [60].

Further investigation of the revealed founder variants in other cancers in Kazakhs could lead to overlapping results, as these PVs cause truncated or absent BRCA1 and BRCA2 proteins and have already been found in other malignancies where loss of function is a known disease mechanism. BRCA1 c.5266dup has previously been found in prostate cancer [61] and pancreatic cancer [62], in addition to its known role in ovarian cancer [6365]. c.5278-2del has been detected in colorectal cancer [50] and intrahepatic cholangiocarcinoma [66] while BRCA2 c.9253del has also been reported for prostate cancer [67].

Novel pathogenic variants

We also identified six novel PVs: BRCA1 c.2498del; PALB2 c.1034T>G, c.18_22del; TP53 c.154C>T, XPA c.20del, and PMS1 c.1258del. The regions affected by these variants appear to be critical for the function of the corresponding proteins. The BRCA1 c.2498del variant results in a frameshift that alters the amino acid sequence of the protein starting at position 2498 and leads to a premature termination codon 13 amino acids downstream. The other frameshift variants XPA c.20del and PMS1 c.1258del generate a premature translational stop signal (p.Ala7ValfsTer8 and p.His420IlefsTer22, respectively) and result in a premature termination codon 8 and 22 amino acids downstream, respectively. The nonsense mutation PALB2 c.1034T>G results from a substitution of T to G at nucleotide position 1034, changing the amino acid from a leucine to a stop codon in the protein, and the TP53 c.154C>T nonsense variant produces a premature translation stop signal (p.Gln52Ter). All of these variants are expected to result in an absent or disrupted protein product. Because these variants have not previously been described in population frequency databases, we hypothesise that they may be specific to the Kazakh population. Our further studies will involve the clinical scrutiny of carriers with these PVs. It will help clarify the role of these pathogenic variants in the development and progression of breast cancer.

Association between genetic profiles and clinicopathological features

Population-based studies indicate that BC patients with specific clinical features are more likely to have PVs in BC susceptibility genes [68]. In this study, BRCA1 PVs carriers had a greater risk of having family history (statistically significant) whereas BRCA2 carriers had worse prognostic outcomes of developing a metastatic phenotype (not statistically significant, Table 3). In our study, the proportion of BRCA1/2 carriers with TNBC was 8.4%, which is within the range of those previously reported in other populations [69]. In line with the existing evidence [68], tumors from BRCA1/2 carriers were more frequently ductal carcinomas, whereas non-BRCA1/2 carriers variants were more likely to have lobular BC. Interestingly, family history of BC was not observed in non-BRCA1/2 carriers, whereas 26.3% of BRCA1/2 carriers and 12.6% of non-carriers had positive BC family history. The lack of family history of BC in non-BRCA1/2 carriers could be explained by either the limited sample size of our study or the moderate/low penetrance of non-BRCA1/2 genes.

Logistic regression analysis showed that the association with TNBC was almost 7-fold greater in BRCA1 carriers, while BRCA2 status was not significantly associated with TNBC. Patients with TNBC were highly enriched in BRCA1 c.5278-2del, c.5382ins, c.2T>C, and c.2498del. It should be noted that BRCA1 carriers were less likely to develop ER, PR, and Her2-positive BC than carriers of BRCA2 variants (Table 2). These results are consistent with the findings from a previous study in Chinese women with BC, which also reported a strong association between BRCA1 variants, but not BRCA2 variants, and TNBC [33]. Unfortunately, we did not have sufficient clinical data or statistical power to perform a subpopulation analysis. The limited sample size would also prevent from obtaining reliable estimates from subpopulation analyses.

Although our efforts focused on BC family history, we were interested in how the identified strong PVs might be related to cancer development on the male side. In our study, this association was traced for cancer of the digestive system. The father of the patient with BRCA2 c.6468_6469delTC/BIC: 6696delTC had laryngeal cancer, and the brother and father of the carrier of PV BRCA1 c.2498delT had gastric and esophageal cancer. Numerous studies confirm that pathogenic BRCA germline variants play a predisposing role for gastrointestinal malignancies [7073]. Interestingly, the analysis of personal data did not reveal a single case of cancer of the male reproductive organs among the patients’ relatives. However, cases of cancer of the male reproductive organs, generally prostate cancer, have been reported in the literature in carriers of germline mutations in the BRCA1, TP53, NBN, PALB2, CHEK2, and ATM genes [61, 7478].

In conclusion, our study comprehensively analyzed the spectrum and prevalence of PVs with early-onset BC in unrelated young women of the Kazakh population. Using NGS-based multigene panel testing, we were able to identify recurrent, possible founder and novel PVs that were undetected in earlier studies due to a less comprehensive and sensitive methodology. Moreover, the identified PVs were differentially associated with biomarkers and prognostic factors. Since the average incidence in young women of all ethnicities in Kazakhstan is 329.6 cases per year [1], we have every reason to discuss the applicability of the results to the general population and the sufficient representativeness of our sample, consisting of 224 young Kazakh women.

We demonstrated the remarkable efficacy of an NGS-based panel to identify rare germline variants in early onset BC patients. These findings could contribute to the development of population-specific multigene panels for more rapid and cost-effective testing. Expanding routine genetic testing for hereditary BC from traditional BRCA testing to multigene panels could improve diagnostic yield, increase cancer prevention options for both identified carriers and their relatives, and reduce the likelihood that hereditary families with early onset BC will be overlooked [79]. Nevertheless, further studies are needed to validate the clinical utility of these panels. The significant gaps in knowledge regarding genotype–phenotype correlations, expressivity, and penetrance should be addressed [80]. To maximize clinical utility and minimize potentially harmful effects, the efficacy of each testing strategy should be evaluated very specifically, from test panels containing only founder PVs to test panels combining genes with high and moderate penetrance or even VUS [81]. A result with a high number of VUS is likely to be more expensive than useful and increases patient anxiety because the likelihood of detecting VUS increases when many genes are tested [82]. At last, the problem of disclosing genetic test results to young patients and family members and the problem of standardizing counseling and clinical management of carriers of gene variants should be reviewed [81]. In this way, the multi-gene panel test for hereditary BC will be adopted by health care providers as a screening tool once clear guidelines are available. Otherwise, inadequate implementation of genetic testing may result in high health care expenditures, wasted time, and other resources without a positive impact on health outcomes [82].

With this in mind, we will focus in the future on segregation analyzes of family members and functional analyzes to evaluate the inheritance pattern and pathogenicity of the identified recurrent and novel BC variants. Retrospective analyses of their possible association with progression-free, metastasis-free, and overall survival are also an exciting direction for future research. No less interesting would be the study of these variants regarding the chemosensitivity and efficacy of specific targeted therapies for their carriers.

Materials and Methods

Study population and data collection

The study included 224 unrelated female patients from the Kazakh ethnic group diagnosed with BC at or under 40 years of age. All patients were treated and followed up at the Kazakh Institute of Oncology and Radiology (KazIOR, Almaty, Kazakhstan) from August 2017 to October 2019. The study protocol was approved by the Ethics Committee of the KazIOR and was conducted in accordance with the Declaration of Helsinki. All the participants signed the informed consent and donated 5 ml of blood for molecular-genetic analysis. The patients also gave consent for publishing the results of their molecular-genetic and clinical data anonymously under a specific ID code assigned for the study. Collected blood in EDTA tubes was transported to the Institute of Genetics and Physiology in a portable refrigerated container within several hours after the collection and frozen at −20°C for further molecular-genetic studies. Information about sociodemographic status, clinicopathological, and family history were obtained from all patients by the ordering physician. ER, PR, HER2, and antigen Ki-67 status were evaluated by IHC using standard procedures and manufacturer protocols. Histological grade was evaluated using the Nottingham grading system.

DNA preparation

Rigid quality control procedures were carried out at several stages of the sequencing process throughout our work. The QIAamp DNA Blood Mini Kit (QIAGEN) was used to isolate DNA from peripheral blood samples while according to manufacturer recommendations. The quantification and quality evaluation of the isolated genomic DNA (gDNA) samples were then performed using spectrophotometry (Biophotometer plus, Eppendorf, Germany) and fluorometry (Qubit fluorometer, Invitrogen, USA). At least 50 ng of double-stranded nuclear DNA per sample were obtained. In order to prepare libraries, gDNA samples had to be broken up, and then tags had to be attached to the resulting pieces. The removal of undesirable components was ensured by subsequent purifying stages, which produced libraries of the highest caliber. Qubit was used to quantify these libraries, and an Agilent Bioanalyzer was used to determine their size distribution.

NGS – library preparation and sequencing

According to the manufacturer’s recommended protocol, we performed NGS using the TruSight Rapid Capture Kit (Illumina) in combination with the TruSight Cancer Sequencing Panel (Illumina). TruSight Cancer Sequencing Panel is a specific enrichment system targeted to the whole coding regions (>1700 exons) and the flanking noncoding regions (on average 50 bp upstream and downstream each exon) of 94 high-risk genes associated with various types of cancer predisposition according to genome-wide association studies (GWAS). The panel covers a total of 255 kb of the human genome. The selection of this panel enables a comprehensive evaluation of putative cancer-related genetic variations.

For library preparation, 50 ng of each gDNA sample was fragmented by Nextera transposome and adapter sequences (tags) were attached to the ends of the fragments. The “tagmented” DNA fragments were cleaned up from the Nextera transposome taking into account the individual bar-code for each patient, and standard adapters required for cluster amplification were added. The “tagmented” DNA was amplified by PCR followed by purification of the amplified fragments. Purification of DNA library was carried out for removing unwanted amplification products using magnetic particles (Sample purification beads, SPB). The library was then qualitatively and quantitatively evaluated using a fluorometric method, such as Qubit Fluorometer (Invitrogen, USA) and by using an Agilent High Sensitivity DNA Kit (Agilent Technologies) on 2100 Bioanalyzer (Agilent Technologies, USA). Next, 500 ng of individual DNA libraries were combined into a single pool in batches of 24 samples followed by double hybridization to capture labeled probes specific to the targeted genomic regions of interest. Nonhybridized material was removed by washing. Before PCR amplification, the captured library was purified using SPB for removing nonhybridized material, and then the enriched library was amplified and re-purified with SPB. The full library was quantified using the Qubit Fluorometer (Invitrogen, USA) and library quality was assessed using Agilent High Sensitivity DNA Kit (Agilent Technologies, USA) on 2100 Bioanalyzer (Agilent Technologies, USA) as previously described. Finally, DNA libraries with 4 mM molarities were subjected to cluster generation using a standard flow cell.

Sequencing was carried out using the MiSeq platform and MiSeq Reagent Kit v3 (Illumina) in paired-end runs of 24 samples of each separate library as recommended by the manufacturer. To analyze 224 samples, 10 MiSeq Reagent Kits v3 were used, with 24 pooled samples each. The MiSeq Reagent Kit v3 provide the highest output among all MiSeq kits. It utilize the same pre-filled, ready-to-use reagent cartridges as the v1 and v2 kits, but incorporate enhanced chemistry to boost cluster density, extend read lengths up to 2 × 300 bp using the 600-cycle kit, and improve quality scores. We used PhiX library as a quality and calibration control for cluster generation, sequencing runs, alignment, and cross-talk matrix generation.

Mean target coverage was 170,7X, and a 95.4% of the targeted regions was covered by at least 20 reads. These values are sufficient to consider the sequencing results reliable and valid.

NGS – bioinformatics data analysis

The bioinformatics data analysis was performed using two methodological approaches accompanying different algorithms. In the first approach, NGS data analysis employed MiSeq Reporter v.3.0 software (Illumina, USA) using the Burrows-Wheeler Aligner (BWA) algorithm. We mapped and aligned the raw demultiplexing reads and compared to the reference human genome (GRCH37/hg19). By utilizing the Genome Analysis Toolkit from (GATK, Broad Institute, Cambridge, USA), we searched and detected variants for specific regions of the genome.

In the second approach, we used various bioinformatics methods and software packages to optimize the workflow. By utilizing the Bowtie2 algorithm with a very sensitive local parameter, we mapped and aligned the sequence reads against the human reference genome version GRCH37/hg19. The quality of reading sequences was evaluated using FastQC and MultiQC. For both the tools default parameters were used. The SAM files were converted into BAM files by Picard Tools and SAMTools/BCF tools. Afterward, we used Java Runtime Environment and R Bioconductor software scripts to sort and index the mapped reads, remove intermediate files, merge the BAM files, and identify duplicates. For re-alignment, the mapped reads around the areas with insertions/deletions GATK were used, and its “Haplotype caller” strategy was used to filter and detect genomic variants. Ultimately, all variant call format (VCF) files obtained from GATK contained only alternative variants imported into the Variant Studio 3.0 software (Illumina, USA). In our study, we have only considered the variants identified by both approaches to increase the robustness of the analysis.

Annotation, interpretation, and classification of variants

The genetic variants were annotated and interpreted by using the Variant Studio 3.0 software (Illumina, USA). The acceptance threshold value was selected in terms of a Q-score of 30, corresponding to an error rate of 0.1%. The poor-quality nucleotide variants were excluded by utilizing the following filtering parameters: (i) a read depth of greater than 50×, (ii) alternative read depth of greater than 20×, and (iii) quality value of greater than 100. In this study, only variants that matched the filtering parameters were selected.

The annotation of genetic variants was performed in accordance with the nomenclature of the Human Genome Variation Society (HGVS) [83], and the interpretation was accomplished using the Single Nucleotide Polymorphism Database (dbSNP, http://www.ncbi.nlm.nih.gov/projects/SNP/), ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar/), LOVD, (https://www.lovd.nl/), BRCA Exchange (https://brcaexchange.org/), and ARUP’ (https://arup.utah.edu/database/). The Integrative Genomics Viewer (IGV) [8485] was employed to visualize both the quality and variance of the BAM and VCF files. Variants were considered rare if the minor allele frequency (MAF) was ≤1% according to 1000G, ESP6500, and ExAC.

Following the guidelines of the American College of Medical Genetics and Genomics (ACMG) [86], the identified variants were categorized as pathogenic, likely pathogenic, VUS, benign, and likely benign. According to the results of functional studies, pathogenic variants included those resulting in a premature stop codon (frameshift and nonsense), variants with uncorrected splicing, and variants affecting protein function. Novel frameshift and nonsense variants leading to protein truncation and variants without well-established functional effects were classified as likely pathogenic.

The ClinVar database was also used to classify variants based on previously established pathogenic or benign effects. Synonymous and intronic variants not affecting splicing were regarded as benign/probably benign. The remaining variants that did not exhibit functional effect and did not satisfy the classification criteria as pathogenic/likely pathogenic or benign/likely benign, or with inconclusive evidence of their benign and pathogenic nature, were defined as VUS.

In silico, bioinformatics tools such as SIFT [17] and PolyPhen-2 [18] were used to predict potential pathogenic effects of missense variants on protein structure and function. We considered variants with a SIFT score of less than 0.05 and a PolyPhen-2 score of more than 0.95 strongly suspected of being deleterious [87].

Statistical analysis

We analyzed the distribution of patients’ characteristics, such as demographics, clinical-pathological characteristics, and personal and family history, between patient groups. The descriptive statistics were summarized as frequency distributions for categorical variables and as means for continuous variables. Kruskal-Wallis test and Fisher’s exact test were used to compare continuous and categorical variables between groups, respectively. The odds ratio (OR) and 95% confidence interval (CI) were calculated through logistical regression. Results were considered statistically significant at a p-value < 0.05. All statistical tests are two-sided at the 5% level of significance. The statistical data were processed and analyzed using the R v4.1.2 software.


BC: breast cancer; TNBC: triple-negative breast cancer; NGS: next-generation sequencing; PV: pathogenic variant; VUS: variants of unknown significance; OR: odds ratio; CI: confidence interval; ER: estrogen receptor; PR: progesterone receptor; HER2: human epidermal growth factor receptor 2; NCCN: National Comprehensive Cancer Network; ACMG: American College of Medical Genetics and Genomics; HGVS: Human Genome Variation Society; LOVD: Leiden Open Variation Database; ARUP’: Associated Regional and University Pathologists; MAF: minor allele frequency; BWA: Burrows-Wheeler Aligner; ExAC: Exome Aggregation Consortium database; ESP6500: Exome Sequencing Project database; 1000G: 1000 Genomes database; GATK: Genome Analysis Toolkit; BAM: Binary Alignment Map; SAM: Sequence Alignment Map; gDNA: genomic DNA; GWAS: genome-wide association studies; IHC: immunohistochemistry; CNS: central nervous system; EDTA: ethylenediaminetetraacetic acid; PCR: polymerase chain reaction.

Author contributions

Gulnur Zhunussova, Dilyara Kaidarova, Leyla Djansugurova: Conceptualization; Gulnur Zhunussova, Nazgul Omarbayeva: Data curation; Gulnur Zhunussova, Nazgul Omarbayeva, Saltanat Abdikerim: Formal analysis; Gulnur Zhunussova, Dilyara Kaidarova: Funding acquisition; Gulnur Zhunussova, Nazgul Omarbayeva, Saltanat Abdikerim, Ilya Kisselev, Kanagat Yergali, Aigul Zhunussova: Investigation; Gulnur Zhunussova, Dilyara Kaidarova, Nazgul Omarbayeva, Saltanat Abdikerim: Methodology; Gulnur Zhunussova, Dilyara Kaidarova: Project administration; Gulnur Zhunussova, Dilyara Kaidarova, Nazgul Omarbayeva, Tatyana Goncharova, Aliya Abdrakhmanova: Resources; Gulnur Zhunussova, Dilyara Kaidarova, Leyla Djansugurova: Supervision; Gulnur Zhunussova, Dilyara Kaidarova, Nazgul Omarbayeva, Saltanat Abdikerim, Natalya Mit, Leyla Djansugurova: Validation; Gulnur Zhunussova, Dilyara Kaidarova, Nazgul Omarbayeva, Saltanat Abdikerim, Tatyana Goncharova, Natalya Mit, Leyla Djansugurova: Visualization; Gulnur Zhunussova: Writing – original draft; Nazgul Omarbayeva, Dilyara Kaidarova, Saltanat Abdikerim, Ilya Kisselev, Aigul Zhunussova, Natalya Mit, Leyla Djansugurova: Writing – review and editing. The work reported in the article has been performed by the authors, unless clearly specified in the text.


The authors are grateful to all the patients enrolled in the study.


Authors have no conflicts of interest to declare.

Ethical statement and consent

The study protocol was approved by the Ethics Committee of the KazIOR and was conducted in accordance with the Declaration of Helsinki. All the participants signed the informed consent and donated 5 ml of blood for molecular-genetic analysis.


This research was funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan (Grant No. AP09259034) and the Ministry of Healthcare of the Republic of Kazakhstan (programs No. OR12165486, No. 49019/ПЦФ-M3CP-OT-18). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Shertaeva A, Ospanova D, Grjibovsky A, Shamsutdinova A, Rakhmetov N, Dushimova Z, Salimgereeva B, Yermentayeva Z, Kaketaeva I, Kuandykov Y, Tanabayeva S, Fakhradiyev I, Zharmenov S. Study on Breast Cancer in Kazakhstan Using the Functional Time Series. Asian Pac J Cancer Prev. 2023; 24:1037–46. https://doi.org/10.31557/APJCP.2023.24.3.1037. [PubMed].

2. Kaidarova D, Shatkovskaya O, Ongarbayev B, Seisenbayeva G, Azhmagambetova A, Zhylkaidarova A, Lavrentyeva I, Sagi M. The indicators of the Oncology Service of the Republic of Kazakhstan, 2021 (Statistical and Analytical Materials)/eds. KIOR JSC. 2022; 384. https://doi.org/10.52532/1-11-2021-1-384.

3. Fernandes U, Guidi G, Martins D, Vieira B, Leal C, Marques C, Freitas F, Dupont M, Ribeiro J, Gomes C, Marques R, Avelar P, Esteves AS, Pinto-de-Sousa J. Breast cancer in young women: a rising threat: A 5-year follow-up comparative study. Porto Biomed J. 2023; 8:e213. https://doi.org/10.1097/j.pbj.0000000000000213. [PubMed].

4. Erić I, Petek Erić A, Kristek J, Koprivčić I, Babić M. Breast cancer in young women: pathologic and immunohistochemical features. Acta Clin Croat. 2018; 57:497–502. https://doi.org/10.20471/acc.2018.57.03.13. [PubMed].

5. Shen L, Zhang S, Wang K, Wang X. Familial Breast Cancer: Disease Related Gene Mutations and Screening Strategies for Chinese Population. Front Oncol. 2021; 11:740227. https://doi.org/10.3389/fonc.2021.740227. [PubMed].

6. Afghahi A, Telli ML, Kurian AW. Genetics of triple-negative breast cancer: Implications for patient care. Curr Probl Cancer. 2016; 40:130–40. https://doi.org/10.1016/j.currproblcancer.2016.09.007. [PubMed].

7. Saleem M, Ghazali MB, Wahab MAMA, Yusoff NM, Mahsin H, Seng CE, Khalid IA, Rahman MNG, Yahaya BH. The BRCA1 and BRCA2 Genes in Early-Onset Breast Cancer Patients. Adv Exp Med Biol. 2020; 1292:1–12. https://doi.org/10.1007/5584_2018_147. [PubMed].

8. Robson M. Management of Women With Breast Cancer and Pathogenic Variants in Genes Other Than BRCA1 or BRCA2. J Clin Oncol. 2021; 39:2528–34. https://doi.org/10.1200/JCO.21.00999. [PubMed].

9. Kim H, Choi DH, Park W. Germline BRCA Mutation and Clinical Outcomes in Breast Cancer Patients Focusing on Survival and Failure Patterns: A Long-Term Follow-Up Study of Koreans. Medicina (Kaunas). 2020; 56:514. https://doi.org/10.3390/medicina56100514. [PubMed].

10. Ripperger T, Gadzicki D, Meindl A, Schlegelberger B. Breast cancer susceptibility: current knowledge and implications for genetic counselling. Eur J Hum Genet. 2009; 17:722–31. https://doi.org/10.1038/ejhg.2008.212. [PubMed].

11. Casaubon JT, Kashyap S, Regan JP. BRCA1 and BRCA2 Mutations. In: StatPearls. Treasure Island (FL): StatPearls Publishing; 2023. [PubMed].

12. Yadav S, Boddicker NJ, Na J, Polley EC, Hu C, Hart SN, Gnanaolivu RD, Larson N, Holtegaard S, Huang H, Dunn CA, Teras LR, Patel AV, et al. Contralateral Breast Cancer Risk Among Carriers of Germline Pathogenic Variants in ATM, BRCA1, BRCA2, CHEK2, and PALB2. J Clin Oncol. 2023; 41:1703–13. https://doi.org/10.1200/JCO.22.01239. [PubMed].

13. Kleiblova P, Stolarova L, Krizova K, Lhota F, Hojny J, Zemankova P, Havranek O, Vocka M, Cerna M, Lhotova K, Borecka M, Janatova M, Soukupova J, et al. Identification of deleterious germline CHEK2 mutations and their association with breast and ovarian cancer. Int J Cancer. 2019; 145:1782–97. https://doi.org/10.1002/ijc.32385. [PubMed].

14. Stucci LS, Internò V, Tucci M, Perrone M, Mannavola F, Palmirotta R, Porta C. The ATM Gene in Breast Cancer: Its Relevance in Clinical Practice. Genes (Basel). 2021; 12:727. https://doi.org/10.3390/genes12050727. [PubMed].

15. Kreipe H, Harbeck N, Christgen M. Clinical validity and clinical utility of Ki67 in early breast cancer. Ther Adv Med Oncol. 2022; 14:17588359221122725. https://doi.org/10.1177/17588359221122725. [PubMed].

16. Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, Hurles ME, McVean GA, and 1000 Genomes Project Consortium. A map of human genome variation from population-scale sequencing. Nature. 2010; 467:1061–73. https://doi.org/10.1038/nature09534. [PubMed].

17. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009; 4:1073–81. https://doi.org/10.1038/nprot.2009.86. [PubMed].

18. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet. 2013; Chapter 7:Unit7.20. https://doi.org/10.1002/0471142905.hg0720s76. [PubMed].

19. Rostami P, Zendehdel K, Shirkoohi R, Ebrahimi E, Ataei M, Imanian H, Najmabadi H, Akbari MR, Sanati MH. Gene Panel Testing in Hereditary Breast Cancer. Arch Iran Med. 2020; 23:155–62. [PubMed].

20. Hu C, Hart SN, Gnanaolivu R, Huang H, Lee KY, Na J, Gao C, Lilyquist J, Yadav S, Boddicker NJ, Samara R, Klebba J, Ambrosone CB, et al. A Population-Based Study of Genes Previously Implicated in Breast Cancer. N Engl J Med. 2021; 384:440–51. https://doi.org/10.1056/NEJMoa2005936. [PubMed].

21. Fostira F, Kostantopoulou I, Apostolou P, Papamentzelopoulou MS, Papadimitriou C, Faliakou E, Christodoulou C, Boukovinas I, Razis E, Tryfonopoulos D, Barbounis V, Vagena A, Vlachos IS, et al. One in three highly selected Greek patients with breast cancer carries a loss-of-function variant in a cancer susceptibility gene. J Med Genet. 2020; 57:53–61. https://doi.org/10.1136/jmedgenet-2019-106189. [PubMed].

22. Bychkovsky B, Rana HQ, Ademuyiwa F, Plichta J, Anderson K, Nogueira-Rodrigues A, Santa-Maria CA, Coffman LG, Marquez C, Das A, Taghian A, Koeller DR, Sandoval RL, et al. Call for action: expanding global access to hereditary cancer genetic testing. Lancet Oncol. 2022; 23:1124–26. https://doi.org/10.1016/S1470-2045(22)00378-3. [PubMed].

23. Mannan AU, Singh J, Lakshmikeshava R, Thota N, Singh S, Sowmya TS, Mishra A, Sinha A, Deshwal S, Soni MR, Chandrasekar A, Ramesh B, Ramamurthy B, et al. Detection of high frequency of mutations in a breast and/or ovarian cancer cohort: implications of embracing a multi-gene panel in molecular diagnosis in India. J Hum Genet. 2016; 61:515–22. https://doi.org/10.1038/jhg.2016.4. [PubMed].

24. Guindalini RSC, Viana DV, Kitajima JPFW, Rocha VM, López RVM, Zheng Y, Freitas É, Monteiro FPM, Valim A, Schlesinger D, Kok F, Olopade OI, Folgueira MAAK. Detection of germline variants in Brazilian breast cancer patients using multigene panel testing. Sci Rep. 2022; 12:4190. https://doi.org/10.1038/s41598-022-07383-1. [PubMed].

25. Zanti M, Loizidou MA, Michailidou K, Pirpa P, Machattou C, Marcou Y, Kyriakou F, Kakouri E, Tanteles GA, Spanou E, Spyrou GM, Kyriacou K, Hadjisavvas A. NGS Panel Testing of Triple-Negative Breast Cancer Patients in Cyprus: A Study of BRCA-Negative Cases. Cancers (Basel). 2020; 12:3140. https://doi.org/10.3390/cancers12113140. [PubMed].

26. Buys SS, Sandbach JF, Gammon A, Patel G, Kidd J, Brown KL, Sharma L, Saam J, Lancaster J, Daly MB. A study of over 35,000 women with breast cancer tested with a 25-gene panel of hereditary cancer genes. Cancer. 2017; 123:1721–30. https://doi.org/10.1002/cncr.30498. [PubMed].

27. Singh J, Thota N, Singh S, Padhi S, Mohan P, Deshwal S, Sur S, Ghosh M, Agarwal A, Sarin R, Ahmed R, Almel S, Chakraborti B, et al. Screening of over 1000 Indian patients with breast and/or ovarian cancer with a multi-gene panel: prevalence of BRCA1/2 and non-BRCA mutations. Breast Cancer Res Treat. 2018; 170:189–96. https://doi.org/10.1007/s10549-018-4726-x. [PubMed].

28. Lincoln SE, Kobayashi Y, Anderson MJ, Yang S, Desmond AJ, Mills MA, Nilsen GB, Jacobs KB, Monzon FA, Kurian AW, Ford JM, Ellisen LW. A Systematic Comparison of Traditional and Multigene Panel Testing for Hereditary Breast and Ovarian Cancer Genes in More Than 1000 Patients. J Mol Diagn. 2015; 17:533–44. https://doi.org/10.1016/j.jmoldx.2015.04.009. [PubMed].

29. Schroeder C, Faust U, Sturm M, Hackmann K, Grundmann K, Harmuth F, Bosse K, Kehrer M, Benkert T, Klink B, Mackenroth L, Betcheva-Krajcir E, Wimberger P, et al. HBOC multi-gene panel testing: comparison of two sequencing centers. Breast Cancer Res Treat. 2015; 152:129–36. https://doi.org/10.1007/s10549-015-3429-9. [PubMed].

30. Susswein LR, Marshall ML, Nusbaum R, Vogel Postula KJ, Weissman SM, Yackowski L, Vaccari EM, Bissonnette J, Booker JK, Cremona ML, Gibellini F, Murphy PD, Pineda-Alvarez DE, et al. Pathogenic and likely pathogenic variant prevalence among the first 10,000 patients referred for next-generation cancer panel testing. Genet Med. 2016; 18:823–32. https://doi.org/10.1038/gim.2015.166. [PubMed].

31. Tung N, Lin NU, Kidd J, Allen BA, Singh N, Wenstrup RJ, Hartman AR, Winer EP, Garber JE. Frequency of Germline Mutations in 25 Cancer Susceptibility Genes in a Sequential Series of Patients With Breast Cancer. J Clin Oncol. 2016; 34:1460–68. https://doi.org/10.1200/JCO.2015.65.0747. [PubMed].

32. Caswell-Jin JL, Gupta T, Hall E, Petrovchich IM, Mills MA, Kingham KE, Koff R, Chun NM, Levonian P, Lebensohn AP, Ford JM, Kurian AW. Racial/ethnic differences in multiple-gene sequencing results for hereditary cancer risk. Genet Med. 2018; 20:234–39. https://doi.org/10.1038/gim.2017.96. [PubMed].

33. Hata C, Nakaoka H, Xiang Y, Wang D, Yang A, Liu D, Liu F, Zou Q, Wei L, Zheng K, Inoue I, You H. Germline mutations of multiple breast cancer-related genes are differentially associated with triple-negative breast cancers and prognostic factors. J Hum Genet. 2020; 65:577–87. https://doi.org/10.1038/s10038-020-0729-7. [PubMed].

34. Horton C, Lo M, Yussuf A, Young C, Carraway C, Agaoglu NB, Bychkovsky BL, Rana HQ. Double jeopardy? A closer look at cancer histories of individuals with multiple germline pathogenic variants. J Clin Oncol. 2022; 40:10514. https://doi.org/10.1200/JCO.2022.40.16_suppl.10514.

35. Wallace SE, Bean LJH. Resources for Genetics Professionals -- Genetic Disorders Associated with Founder Variants Common in the Cree and/or Ojibway Population. In: Adam MP, Mirzaa GM, Pagon RA, et al, editors. GeneReviews®. Seattle (WA): University of Washington, Seattle; 1993-2023; NBK541231. https://www.ncbi.nlm.nih.gov/books/.

36. Papamentzelopoulou M, Apostolou P, Fostira F, Dimitrakakis C, Loutradis D, Fountzilas G, Yannoukakos D, Konstantopoulou I. Prevalence and founder effect of the BRCA1 p.(Val1833Met) variant in the Greek population, with further evidence for pathogenicity and risk modification. Cancer Genet. 2019; 237:90–96. https://doi.org/10.1016/j.cancergen.2019.06.006. [PubMed].

37. Meng H, Yao L, Yuan H, Xu Y, Ouyang T, Li J, Wang T, Fan Z, Fan T, Lin B, Xie Y. BRCA1 c.5470_5477del, a founder mutation in Chinese Han breast cancer patients. Int J Cancer. 2020; 146:3044–52. https://doi.org/10.1002/ijc.32877. [PubMed].

38. Cybulski C, Kluźniak W, Huzarski T, Wokołorczyk D, Kashyap A, Rusak B, Stempa K, Gronwald J, Szymiczek A, Bagherzadeh M, Jakubowska A, Dębniak T, Lener M, et al. The spectrum of mutations predisposing to familial breast cancer in Poland. Int J Cancer. 2019; 145:3311–20. https://doi.org/10.1002/ijc.32492. [PubMed].

39. Lieberman S, Chen-Shtoyerman R, Levi Z, Shkedi-Rafid S, Zuckerman S, Bernstein-Molho R, Levi GR, Shachar SS, Flugelman A, Libman V, Kedar I, Naftaly-Nathan S, Lagovsky I, et al. Common founder BRCA2 pathogenic variants and breast cancer characteristics in Ethiopian Jews. Breast Cancer Res Treat. 2022; 193:217–24. https://doi.org/10.1007/s10549-022-06557-9. [PubMed].

40. Colombo M, De Vecchi G, Caleca L, Foglia C, Ripamonti CB, Ficarazzi F, Barile M, Varesco L, Peissel B, Manoukian S, Radice P. Comparative in vitro and in silico analyses of variants in splicing regions of BRCA1 and BRCA2 genes and characterization of novel pathogenic mutations. PLoS One. 2013; 8:e57173. https://doi.org/10.1371/journal.pone.0057173. [PubMed].

41. Laitman Y, Friebel TM, Yannoukakos D, Fostira F, Konstantopoulou I, Figlioli G, Bonanni B, Manoukian S, Zuradelli M, Tondini C, Pasini B, Peterlongo P, Plaseska-Karanfilska D, et al. The spectrum of BRCA1 and BRCA2 pathogenic sequence variants in Middle Eastern, North African, and South European countries. Hum Mutat. 2019; 40:e1–23. https://doi.org/10.1002/humu.23842. [PubMed].

42. ElBiad O, Laraqui A, El Boukhrissi F, Mounjid C, Lamsisi M, Bajjou T, Elannaz H, Lahlou AI, Kouach J, Benchekroune K, Oukabli M, Chahdi H, Ennaji MM, et al. Prevalence of specific and recurrent/founder pathogenic variants in BRCA genes in breast and ovarian cancer in North Africa. BMC Cancer. 2022; 22:208. https://doi.org/10.1186/s12885-022-09181-4. [PubMed].

43. Abdikhakimov A, Tukhtaboeva M, Adilov B, Turdikulova S. The Potential Contribution of BRCA Mutations to Early Onset and Familial Breast Cancer in Uzbekistan. Cent Asian J Glob Health. 2016; 5:228. https://doi.org/10.5195/cajgh.2016.228. [PubMed].

44. Qin Z, Li J, Tam B, Sinha S, Zhao B, Bhaskaran SP, Huang T, Wu X, Chian JS, Guo M, Kou SH, Lei H, Zhang L, et al. Ethnic-specificity, evolution origin and deleteriousness of Asian BRCA variation revealed by over 7500 BRCA variants derived from Asian population. Int J Cancer. 2023; 152:1159–73. https://doi.org/10.1002/ijc.34359. [PubMed].

45. Shagimardanova E, Shigapova L, Gusev O, Nikitin A, Druzhkov M, Enikeev R, Gordiev M. Germline BRCA screening in breast cancer patients in Tatar women population. Ann Oncol. 2016 (Suppl 6); 27:vi43–67. https://doi.org/10.1093/annonc/mdw364.26.

46. Anisimenko MS, Mitrofanov DV, Chasovnikova OB, Voevoda MI, Kovalenko SP. BRCA1 gene mutations frequency estimation by allele-specific real-time PCR of pooled genomic DNA samples. Breast. 2013; 22:532–36. https://doi.org/10.1016/j.breast.2012.12.007. [PubMed].

47. Balmukhanov T, Hanseitova A, Nigmatova V, Ashirbekov E, Popova I, Prnazarova A, Cherusheva A, Hodaeva A, Talaeva SH, Aithojina N. The prevalence of some polymorphisms in the BRCA1 and BRCA2 genes in breast cancer among the population of the Republic of Kazakhstan. Reports of the National Academy of Sciences of the Republic of Kazakhstan. 2012; 4:58–63.

48. Akilzhanova AR, Nyshanbekkyzy B, Nurkina ZM, Shtephanov II, Makishev AK, Adylkhanov TA, Rakhypbekov TK, Ramanculov EM, Momynaliev KT. BRCA1 and BRCA2 Gene Mutations Screening In Sporadic Breast Cancer Patients In Kazakhstan. Cent Asian J Glob Health. 2013; 2:29. https://doi.org/10.5195/cajgh.2013.29. [PubMed].

49. Rebbeck TR, Friebel TM, Friedman E, Hamann U, Huo D, Kwong A, Olah E, Olopade OI, Solano AR, Teo SH, Thomassen M, Weitzel JN, Chan TL, et al. Mutational spectrum in a worldwide study of 29,700 families with BRCA1 or BRCA2 mutations. Hum Mutat. 2018; 39:593–620. https://doi.org/10.1002/humu.23406. [PubMed].

50. Zhunussova G, Afonin G, Abdikerim S, Jumanov A, Perfilyeva A, Kaidarova D, Djansugurova L. Mutation Spectrum of Cancer-Associated Genes in Patients With Early Onset of Colorectal Cancer. Front Oncol. 2019; 9:673. https://doi.org/10.3389/fonc.2019.00673. [PubMed].

51. Phelan CM, Iqbal J, Lynch HT, Lubinski J, Gronwald J, Moller P, Ghadirian P, Foulkes WD, Armel S, Eisen A, Neuhausen SL, Senter L, Singer CF, et al. Incidence of colorectal cancer in BRCA1 and BRCA2 mutation carriers: results from a follow-up study. Br J Cancer. 2014; 110:530–34. https://doi.org/10.1038/bjc.2013.741. [PubMed].

52. Kwong A, Ho JCW, Shin VY, Kurian AW, Tai E, Esserman LJ, Weitzel JN, Lin PH, Field M, Domchek SM, Lo J, Ngan HYS, Ma ESK, et al. Rapid detection of BRCA1/2 recurrent mutations in Chinese breast and ovarian cancer patients with multiplex SNaPshot genotyping panels. Oncotarget. 2018; 9:7832–43. https://doi.org/10.18632/oncotarget.23471. [PubMed].

53. Kang E, Seong MW, Park SK, Lee JW, Lee J, Kim LS, Lee JE, Kim SY, Jeong J, Han SA, Kim SW, and Korean Hereditary Breast Cancer Study Group. The prevalence and spectrum of BRCA1 and BRCA2 mutations in Korean population: recent update of the Korean Hereditary Breast Cancer (KOHBRA) study. Breast Cancer Res Treat. 2015; 151:157–68. https://doi.org/10.1007/s10549-015-3377-4. [PubMed].

54. Son BH, Ahn SH, Kim SW, Kang E, Park SK, Lee MH, Noh WC, Kim LS, Jung Y, Kim KS, Noh DY, Moon BI, Suh YJ, et al. Prevalence of BRCA1 and BRCA2 mutations in non-familial breast cancer patients with high risks in Korea: the Korean Hereditary Breast Cancer (KOHBRA) Study. Breast Cancer Res Treat. 2012; 133:1143–52. https://doi.org/10.1007/s10549-012-2001-0. [PubMed].

55. Kim H, Cho DY, Choi DH, Choi SY, Shin I, Park W, Huh SJ, Han SH, Lee MH, Ahn SH, Son BH, Kim SW, Haffty BG, and Korean Breast Cancer Study Group. Characteristics and spectrum of BRCA1 and BRCA2 mutations in 3,922 Korean patients with breast and ovarian cancer. Breast Cancer Res Treat. 2012; 134:1315–26. https://doi.org/10.1007/s10549-012-2159-5. [PubMed].

56. Ou J, Wu T, Sijmons R, Ni D, Xu W, Upur H. Prevalence of BRCA1 and BRCA2 Germline Mutations in Breast Cancer Women of Multiple Ethnic Region in Northwest China. J Breast Cancer. 2013; 16:50–54. https://doi.org/10.4048/jbc.2013.16.1.50. [PubMed].

57. Meindl A, and German Consortium for Hereditary Breast and Ovarian Cancer. Comprehensive analysis of 989 patients with breast or ovarian cancer provides BRCA1 and BRCA2 mutation profiles and frequencies for the German population. Int J Cancer. 2002; 97:472–80. https://doi.org/10.1002/ijc.1626. [PubMed].

58. Cunningham JM, Cicek MS, Larson NB, Davila J, Wang C, Larson MC, Song H, Dicks EM, Harrington P, Wick M, Winterhoff BJ, Hamidi H, Konecny GE, et al. Clinical characteristics of ovarian cancer classified by BRCA1, BRCA2, and RAD51C status. Sci Rep. 2014; 4:4026. https://doi.org/10.1038/srep04026. [PubMed].

59. Sokolenko AP, Sokolova TN, Ni VI, Preobrazhenskaya EV, Iyevleva AG, Aleksakhina SN, Romanko AA, Bessonov AA, Gorodnova TV, Anisimova EI, Savonevich EL, Bizin IV, Stepanov IA, et al. Frequency and spectrum of founder and non-founder BRCA1 and BRCA2 mutations in a large series of Russian breast cancer and ovarian cancer patients. Breast Cancer Res Treat. 2020; 184:229–35. https://doi.org/10.1007/s10549-020-05827-8. [PubMed].

60. Cao YW, Fu XG, Wan GX, Yu SY, Cui XB, Li L, Jiang JF, Zheng YQ, Zhang WJ, Li F. BRCA1 gene exon 11 mutations in Uighur and Han women with early-onset sporadic breast cancer in the northwest region of China. Asian Pac J Cancer Prev. 2014; 15:4513–18. https://doi.org/10.7314/apjcp.2014.15.11.4513. [PubMed].

61. Boilève A, Lavaud P, Caron O. Germline BRCA1 Mutation and Prostate Cancer: Be Careful on Causality. Eur Urol Oncol. 2021; 4:674–75. https://doi.org/10.1016/j.euo.2021.06.004. [PubMed].

62. Vietri MT, D’Elia G, Caliendo G, Albanese L, Signoriello G, Napoli C, Molinari AM. Pancreatic Cancer with Mutation in BRCA1/2, MLH1, and APC Genes: Phenotype Correlation and Detection of a Novel Germline BRCA2 Mutation. Genes (Basel). 2022; 13:321. https://doi.org/10.3390/genes13020321. [PubMed].

63. Janavičius R. Founder BRCA1/2 mutations in the Europe: implications for hereditary breast-ovarian cancer prevention and control. EPMA J. 2010; 1:397–412. https://doi.org/10.1007/s13167-010-0037-y. [PubMed].

64. Gorodetska I, Serga S, Levkovich N, Lahuta T, Ostapchenko L, Demydov S, Anikusko N, Cheshuk V, Smolanka I, Sklyar S, Polenkov S, Boichenko O, Kozeretska I. The frequency of BRCA1 founder mutation c.5266dupC (5382insC) in breast cancer patients from Ukraine. Hered Cancer Clin Pract. 2015; 13:19. https://doi.org/10.1186/s13053-015-0040-3. [PubMed].

65. Bahsi T, Erdem HB. Spectrum of BRCA1/BRCA2 variants in 1419 Turkish breast and ovarian cancer patients: a single center study. Turk J Biochem. 2019; 45:83–90. https://doi.org/10.1515/tjb-2019-0424.

66. De Bonis M, Mannavola F, Salvatore L, De Paolis E, Nero C, Giacò L, Tortora G, Giuliante F, Urbani A, Scambia G, Normanno N, Minucci A. Unexpected finding of a rare pathogenic germline BRCA1 variant in an intrahepatic cholangiocarcinoma using the Oncomine Focus DNA assay: clinical and diagnostic implications. Mol Biol Rep. 2023; 50:7113–20. https://doi.org/10.1007/s11033-023-08588-x. [PubMed].

67. Brandão A, Paulo P, Teixeira MR. Hereditary Predisposition to Prostate Cancer: From Genetics to Clinical Implications. Int J Mol Sci. 2020; 21:5036. https://doi.org/10.3390/ijms21145036. [PubMed].

68. Larsen MJ, Thomassen M, Gerdes AM, Kruse TA. Hereditary breast cancer: clinical, pathological and molecular characteristics. Breast Cancer (Auckl). 2014; 8:145–55. https://doi.org/10.4137/BCBCR.S18715. [PubMed].

69. Shimelis H, LaDuca H, Hu C, Hart SN, Na J, Thomas A, Akinhanmi M, Moore RM, Brauch H, Cox A, Eccles DM, Ewart-Toland A, Fasching PA, et al. Triple-Negative Breast Cancer Risk Genes Identified by Multigene Hereditary Cancer Panel Testing. J Natl Cancer Inst. 2018; 110:855–62. https://doi.org/10.1093/jnci/djy106. [PubMed].

70. Correa T, Laux DE, Hoffman HT. A pathogenic germline BRCA2 variant in a patient with hypopharyngeal squamous cell carcinoma. Clin Case Rep. 2021; 9:429–32. https://doi.org/10.1002/ccr3.3548. [PubMed].

71. Premnath N, O’Reilly EM. BReast CAncer (BRCA) gene mutations as an emerging biomarker for the treatment of gastrointestinal malignancies. Chin Clin Oncol. 2020; 9:64. https://doi.org/10.21037/cco-2019-ddp-05. [PubMed].

72. Maccaroni E, Giampieri R, Lenci E, Scortichini L, Bianchi F, Belvederesi L, Brugiati C, Pagliaretta S, Ambrosini E, Berardi R. BRCA mutations and gastrointestinal cancers: When to expect the unexpected? World J Clin Oncol. 2021; 12:565–80. https://doi.org/10.5306/wjco.v12.i7.565. [PubMed].

73. Molinaro E, Andrikou K, Casadei-Gardini A, Rovesti G. BRCA in Gastrointestinal Cancers: Current Treatments and Future Perspectives. Cancers (Basel). 2020; 12:3346. https://doi.org/10.3390/cancers12113346. [PubMed].

74. Rusak B, Kluźniak W, Wokołorczykv D, Stempa K, Kashyap A, Gronwald J, Huzarski T, Dębniak T, Jakubowska A, Masojć B, Akbari MR, Narodv SA, Lubiński J, Cybulski C. Inherited NBN Mutations and Prostate Cancer Risk and Survival. Cancer Res Treat. 2019; 51:1180–87. https://doi.org/10.4143/crt.2018.532. [PubMed].

75. Maxwell KN, Cheng HH, Powers J, Gulati R, Ledet EM, Morrison C, Le A, Hausler R, Stopfer J, Hyman S, Kohlmann W, Naumer A, Vagher J, et al. Inherited TP53 Variants and Risk of Prostate Cancer. Eur Urol. 2022; 81:243–50. https://doi.org/10.1016/j.eururo.2021.10.036. [PubMed].

76. Wokołorczyk D, Kluźniak W, Stempa K, Rusak B, Huzarski T, Gronwald J, Gliniewicz K, Kashyap A, Morawska S, Dębniak T, Jakubowska A, Szwiec M, Domagała P, et al. PALB2 mutations and prostate cancer risk and survival. Br J Cancer. 2021; 125:569–75. https://doi.org/10.1038/s41416-021-01410-0. [PubMed].

77. Karlsson Q, Brook MN, Dadaev T, Wakerell S, Saunders EJ, Muir K, Neal DE, Giles GG, MacInnis RJ, Thibodeau SN, McDonnell SK, Cannon-Albright L, Teixeira MR, et al. Rare Germline Variants in ATM Predispose to Prostate Cancer: A PRACTICAL Consortium Study. Eur Urol Oncol. 2021; 4:570–79. https://doi.org/10.1016/j.euo.2020.12.001. [PubMed].

78. Wu Y, Yu H, Zheng SL, Na R, Mamawala M, Landis T, Wiley K, Petkewicz J, Shah S, Shi Z, Novakovic K, McGuire M, Brendler CB, et al. A comprehensive evaluation of CHEK2 germline mutations in men with prostate cancer. Prostate. 2018; 78:607–15. https://doi.org/10.1002/pros.23505. [PubMed].

79. Asphaug L, Melberg HO. The Cost-Effectiveness of Multigene Panel Testing for Hereditary Breast and Ovarian Cancer in Norway. MDM Policy Pract. 2019; 4:2381468318821103. https://doi.org/10.1177/2381468318821103. [PubMed].

80. Slavin TP, Niell-Swiller M, Solomon I, Nehoray B, Rybak C, Blazer KR, Weitzel JN. Clinical Application of Multigene Panels: Challenges of Next-Generation Counseling and Cancer Risk Management. Front Oncol. 2015; 5:208. https://doi.org/10.3389/fonc.2015.00208. [PubMed].

81. Reid S, Pal T. Update on multi-gene panel testing and communication of genetic test results. Breast J. 2020; 26:1513–19. https://doi.org/10.1111/tbj.13971. [PubMed].

82. Catana A, Apostu AP, Antemie RG. Multi gene panel testing for hereditary breast cancer - is it ready to be used? Med Pharm Rep. 2019; 92:220–25. https://doi.org/10.15386/mpr-1083. [PubMed].

83. Nomenclature for the description of sequence variants. Human Genome Variation Society. 2013. https://www.hgvs.org/mutnomen.

84. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011; 29:24–26. https://doi.org/10.1038/nbt.1754. [PubMed].

85. Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013; 14:178–92. https://doi.org/10.1093/bib/bbs017. [PubMed].

86. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, Voelkerding K, Rehm HL, and ACMG Laboratory Quality Assurance Committee. 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:405–24. https://doi.org/10.1038/gim.2015.30. [PubMed].

87. Lin PH, Kuo WH, Huang AC, Lu YS, Lin CH, Kuo SH, Wang MY, Liu CY, Cheng FT, Yeh MH, Li HY, Yang YH, Hsu YH, et al. Multiple gene sequencing for risk assessment in patients with early-onset or familial breast cancer. Oncotarget. 2016; 7:8310–20. https://doi.org/10.18632/oncotarget.7027. [PubMed].

Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 28518