Assessing the Variations in Breast/Ovarian Cancer Risk for Chinese <i>BRCA1/2</i> Carriers

Li, Ang; Hao, Shuai; Luo, Jiaqi; Zi, Yi; Lan, Zhaoji; Zhou, Tianliangwen; Zhi, Qihuan; Zhan, Jiamin; Sun, Gang; Shi, Yujian; Luo, Donglin

doi:https://doi.org/10.1155/2022/9390539

Journal of Oncology

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Additional Points Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 9390539 | https://doi.org/10.1155/2022/9390539

Assessing the Variations in Breast/Ovarian Cancer Risk for Chinese BRCA1/2 Carriers

Ang Li,¹Shuai Hao,^2,3Jiaqi Luo,¹Yi Zi,¹Zhaoji Lan,¹Tianliangwen Zhou,¹Qihuan Zhi,¹Jiamin Zhan,¹Gang Sun,¹Yujian Shi,¹and Donglin Luo²

Academic Editor: Francesca De Felice

Received23 Jan 2022

Accepted28 Feb 2022

Published26 Mar 2022

Abstract

Background. Cancer risks vary in different BRCA1/2 mutations. We are interested in identifying regions associated with elevated/reduced risks of breast/ovarian cancers in the Chinese population and comparing with previously reported Caucasian-based breast/ovarian cancer cluster regions (OCCR/BCCR). We also aim to characterize the distribution and estimate the cancer risks of different Chinese recurrent mutations. Methods. A total of 3,641 cancer-free women and 4,278 female cancer patients were included in the study. Germline BRCA1/2 status was detected with amplicon-based next-generation sequencing. We calculated the odds ratio (OR) of breast cancer and OR of ovarian cancer, and their ratio of the two ORs (ROR) for each region. ROR >1 indicated elevated odds of breast cancer and/or decreasing odds of ovarian cancer, and vice versa. The frequency, distribution, and penetrance of six known Chinese founder mutations were characterized, respectively. Haplotype analysis and age estimation were performed on the most prevalent founder mutation BRCA1: c.5470_5477del. Results. A total of 729 subjects were detected with germline BRCA1/2 deleterious mutations. The putative Chinese OCCR/BCCR partially overlapped with Caucasian-based OCCR/BCCR and shared structural-functional characteristics. The six known Chinese founder mutations greatly vary in both distribution and penetrance. The two widely spread mutations are estimated to convey low penetrance, while the area-restricted founder mutations seemed to confer higher/complete penetrance. BRCA1: c.5470_5477del is estimated to have emerged ∼2,090 years ago (70 B.C.) during the Han dynasty. Conclusions. BRCA1/2 carriers with different genotypes have significantly different cancer risks. An optimal risk assessment should be mutation specific, rather than concerning a single figure.

1. Introduction

Since the establishment of predisposing effects of BRCA1 (MIM:113705) and BRCA2 (MIM: 114480) to breast/ovarian cancers, we have accumulated understandings of the biological roles of BRCA1/2 and their prevalence of mutations in different ethnic populations through a lot of functional experiments and sample-based studies. It is now generally believed that deleterious mutations in BRCA1/2 or other related genes result in the production of malformed proteins that are unable to function during the error-free DNA repair process mediated by homologous recombination (HR) upon DNA double-strand break, which will lead to genomic instabilities and eventually the development of malignancies. It is also known that race/ethnic differences are presented in the mutation spectrum, the prevalence of mutations, and in the recurrently mutated positions, which are likely to reflect the so-called “founder effects.” For example, Ashkenazi Jews generally have a higher risk of being BRCA1/2 carriers because of the highly prevalent founder mutations BRCA1 185delAG, BRCA1 5382insC, and BRCA2 6174delT [1]. Several Chinese founder mutations have also been previously reported [2–6].

Moreover, it is recognized that great differences in cancer risks present within BRCA1/2 mutation carriers, depending on the location and type of mutation they bear. According to previous observations, mutations within certain regions (defined as ovarian cancer cluster regions, OCCRs) are associated with higher ovarian cancer and/or lower breast cancer risks than other regions and certain regions (defined as breast cancer cluster regions, BCCRs) with higher breast cancer and/or lower ovarian cancer risks. The biological impact caused by the mutation is the determining factor of cancer risk, which apparently affects the position of the OCCRs and BCCRs; the race/ethnicity of the studied population is another variable as it affects the mutation spectrum (e.g., the position of mutation hotspots). Several studies have calculated putative OCCR/BCCRs of BRCA1/2 using samples mostly of Caucasian origin, and the estimated regions are considerably overlapping [7–9]. Interestingly, for both BRCA1 and BRCA2, the estimated OCCRs seemed to locate in the center of the coding sequence (CDS) and significantly overlap with the largest exon (Exon10 of BRCA1 and Exon11 of BRCA2), while the BCCRs seemed to occupy the 5’ and the 3’ ends of the CDS. As we expect race/ethnic difference to cause some degree of variability in OCCR/BCCRs when studying different populations, and since Asian populations only represented 1% of the total samples in previous studies, we are interested in defining OCCR/BCCRs in Chinese BRCA1/2 carriers and finding out to what extent the position of OCCR/BCCRs could be varied. Moreover, it has been widely acknowledged that the identification and screening of founder mutations are highly cost-effective measures for cancer risk management [10]. Hence, we also aim to characterize the distribution of Chinese recurrent mutations (or founder mutations) and estimate the cancer risks they confer, in order to provide a reference to the precision management of genetic risks for Chinese BRCA1/2 carriers.

2. Methods

2.1. Samples

A total of 3,641 cancer-free women and 4,278 female cancer patients (breast cancer, ovarian cancer, colon cancers, and pancreatic adenocarcinoma) in the nationwide of China volunteered to enroll in this study. The participants selected for further study should meet the following criteria: (1) ≥18 years old with unambiguous clinical status (cancer-free/cancer); (2) have taken the BRCA1/2 gene test; and (3) have unambiguous pathology and a sufficient amount of materials for the experiment if they were diagnosed with cancer. Clinical characteristics of samples in relation to BRCA1/2 are presented in Table 1. Among the 3,641 cancer-free individuals, 2,615 were negative of family history and further selected as a healthy control for penetrance estimation and haplotype analysis. All the subjects involved had given their written informed consent in accordance with the Chinese ethical standards and the 2008 Helsinki Declaration.

2.2. BRCA1/2 Genetic Testing and Mutation Classification

Blood samples were collected from all enrolled individuals and subjected to NGS-based BRCA1/2 whole-exon sequencing (all coding regions and exon-intron boundaries20 bp) with average sequencing depth >1000x on the NextSeq CN500 platform (Berry Genomics, China). The analysis pipeline and mutation classification followed protocols as described in our previous research [11]. The human genome hg19/GRC37 was used as reference; the NCBI reference sequences NM_007294.3 and NM_000059.3 were used for annotations of BRCA1 and BRCA2 variants, respectively. Nomenclature of mutations followed the latest version of Human Genome Variation Society Sequence Variant Nomenclature (HGVS, https://varnomen.hgvs.org/) and Mutalyzer Name Checker (https://mutalyzer.nl). Deleterious mutations in this study include likely pathogenic and pathogenic mutations. Further validations of all deleterious mutations and variants of uncertain significance were performed by Sanger sequencing. Visualization of the variants was presented with circos plots [12].

2.3. Mutation Grouping and Statistical Analysis

To estimate OCCR/BCCR, segments of regions containing all deleterious mutations (regardless of mutation type or function) need to be created for statistical calculations. We divided the coding sequence (CDS) of BRCA1 and BRCA2 into bins by base pair location so that each bin contains a roughly equal number of carriers (Table 2). Large rearrangements were excluded to avoid spanning multiple bins.

We then calculate the odds ratio (OR) of breast cancer and the OR of ovarian cancer, respectively, for each bin. We computed a statistical measure ROR, defined as the ratio of breast versus ovarian cancer OR. The value of ROR is associated with elevated/reduced breast or ovarian cancer odds (ROR >1, increasing odds of breast cancer and/or decreasing odds of ovarian cancer; ROR <1, increasing odds of ovarian cancer and/or decreasing breast cancer odds). RORs of all bins were compared to identify significant outliers as putative OCCR/BCCRs. values were adjusted with the Benjamini-Hochberg procedure to reduce the possibility of false positive. Adjusted value <0.05 was considered significant.

2.4. Characterization of BRCA1: c.5470_5477delATTGGGCA

Among the recurrent mutations identified in this study, BRCA1: c.5470_5477delATTGGGCA occurred in 9.5% (50/527) of BRCA1 carriers and had enough materials for further experimental validation. We then performed haplotype analysis on this variant. A total of 81 subjects (31 patients with the same mutations from independent families and 50 unrelated controls without the mutation) were included in the haplotype analysis. According to marker selection methods previously described [4, 13, 14], we selected nine polymorphic markers flanking the BRCA1 gene-D17S800, D17S1320, D17S1321, D17S855, D17S1323, D17S1327, D17S1326, D17S1325, and D17S791 (short tandem repeats, STR, spanning approximately 5.8 Mbp, Supplementary Figure S1). We detected STR lengths using fluorescently end-labeled PCR primers and ABI 3730xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). We reconstructed all possible haplotypes using PHASE v.2.1.1 [15] (Supplementary Table S1). A chi-square test was used to evaluate the differences in allele frequencies between patients and healthy subjects for each STR involved.

We also performed an age estimation of BRCA1: c.5470_5477delATTGGGCA. We used DMLE + v2.3 to estimate the original time when the mutation emerged in the BRCA1 gene [16]. Based on a Markov Chain Monte Carlo algorithm, this program enabled Bayesian inference to estimate the age of the specific mutation with the knowledge of the patients’ and controls’ haplotypes or genotypes observed, physical distances between markers (cM), and the estimated population growth rate. The population statistics of ancient time points were taken from historical data [17]; the population figure of the year 1949 was collected from the National Bureau of Statistics of China (https://data.stats.gov.cn/index.htm). The equation (1) shown below was performed to calculate a list of average growth rates (‰ per year) from official population records [17] of different historical time points ranging from 684 B.C. to 1949 A.D (Table 3). The calculated growth rates slightly differ from each other, with a maximum of 7.81‰ (1741 A.D-present) and a minimum of 1.54‰ (2 A.D-present). For reference, the estimated growth rate of the year 2019 is 3.34‰ (https://data.stats.gov.cn/index.htm).X = average growth rate, n = time till now (years), = population as recorded n years ago, and = Chinese population in 2020 = 1.4 billion.

2.5. Penetrance Estimations on Recurrent Mutations for Breast Cancer

We estimated the penetrance of the recurrent mutations of breast cancer using the allelic model by Bayes’ theorem [18]. This formula required the knowledge of the lifetime risk of breast cancer in the Chinese population. Restricted by the access to the official data, we used the estimated mean value of lifetime risk (0.053; Figure 1, gray vertical line) for breast cancer derived from the Gail model among Chinese women [19] as the population baseline risk. The approximate lower bound (0.35; Figure 1, black dashed line) of breast/ovarian cancer risk in BRCA1/2 carriers (according to previous studies [20, 21]) was used as a minimum expected penetrance of all deleterious mutations throughout the whole gene.

3. Results

Out of the total 7,919 samples containing 4,278 cancer patients and 3,641 cancer-free individuals, 729 were detected with BRCA1/2 deleterious germline mutations, of which 539 individuals carried BRCA1 mutations and 213 carried BRCA2 mutations. The median age of breast cancer diagnosis was 43, 44, and 46 in BRCA1 carriers, BRCA2 carriers, and noncarriers; the median age of ovarian cancer diagnosis was 51, 54, and 52 in BRCA1, BRCA2 carriers, and noncarriers, respectively.

A total of 236 deleterious germline mutations in BRCA1 (Figure 2) and 122 in BRCA2 (Figure 3) were detected and verified in this study. Several recurrent mutations are identified (Figures 2 and 3 highlighted with red, bigger font). BRCA1: c.5470_5477delATTGGGCA (Ile1824AspfsTer3) was the top hit, accounting for 9.5% (50/527) BRCA1 carriers; the second hit for BRCA1 was c.981_982delAT, with 4.0% (21/527) prevalence in BRCA1 carriers; next were c.3770_3771delAG (12/527, 2.3% of BRCA1 carriers), c.5521delA (11/527, 2.1%), and c.4801A > T (8/527, 1.5%). The BRCA2: c.5722_5723delCT was the top hit for BRCA2 (10/202, 5.0%), followed by BRCA2: c.3109C > T (9/202, 4.5%).

Figure 2

Circos plot representation of all BRCA1 deleterious mutations identified in this study. The outmost ring displays the 23 exons of BRCA1 and each variant at corresponding positions. Recurrent variants are highlighted in a red and enlarged font. The second circle is the mutation density graph, each dot corresponds to a variant, and its distance to the outmost ring represents the frequency (the closer the higher; the most prevalent mutation is colored in red). The ROR is represented by a heatmap in the intermediate ring. The next circle displays the BRCA1 functional domains. The innermost arcs represent previously reported Caucasian-based OCCR/BCCRs. The areas enclosed by the arches indicate the estimated Chinese OCCR/BCCR, which are statistically significant (FDR value <0.05).

Figure 3

Circos plot representation of all BRCA2 deleterious mutations identified in this study. The outmost ring displays the 27 exons of BRCA2 and each variant at corresponding positions. Recurrent variants are highlighted in a red and enlarged font. The second circle is the mutation density graph, each dot corresponds to a variant, and its distance to the outmost ring represents the frequency (the closer the higher; the most prevalent mutation is colored in red). The ROR is represented by the heatmap in the intermediate ring. The next circle displays the BRCA2 functional domains. The innermost arcs represent previously reported Caucasian-based OCCR/BCCRs. The areas enclosed by the arches indicate the estimated Chinese OCCR/BCCR, which are statistically significant (FDR value <0.05).

3.1. Estimated Chinese OCCR/BCCR in BRCA1

We predicted an OCCR at c.1154–c.2111 (Table 2; Figure 2, blue arch) with a relative decrease in breast cancer risk and a relative increase in ovarian cancer risk (ROR = 0.29; 95% CI = 0.15–0.57; FDR-corrected value = 1.2 × 10⁻³). The putative OCCR lies within the largest exon (Exon 10) of BRCA1 and is partially overlapped with previously reported OCCR (c.1380–c.4062) [8] (Figure 2, blue arc). The putative OCCR entirely or partially spans several functional domains, including the binding sites for Rb, Rad50, c-Myc, and the nuclear localization sequence (NLS) [22]. A putative BCCR was found at c.5470–c.5524 near the 3’ end of the CDS (Table 2; Figure 2, orange arch), explained by a relative decrease in ovarian cancer risk (ROR = 3.12, 95% CI = 1.56–5.88, FDR-corrected = 1.2 × 10⁻³). The BCCR lies within the second BRCT domain (c.5268–c.5526) and is also partially overlapped with previously reported BCCR (c.5261–c.5563) [8] (Figure 2, orange arc).

3.2. Estimated Chinese OCCR/BCCR in BRCA2

We predicted an OCCR at c.5745–c.7805 + 1 (Table 2; Figure 3, blue arch) with a relative decrease in breast cancer risk and a relative increase in ovarian cancer risk (ROR = 0.38; 95% CI = 0.16–0.86; FDR-corrected value = 0.044). The putative OCCR is partially overlapped with previously reported OCCR (c.6645–c.7471) [8] (Figure 3, blue arc), spanning from the last two BRC repeats within exon 11 (the largest exon of BRCA2) to approximately the boundary of exons 17 and 18, which contains the most part of the helical DNA-binding domain (c.7437–c.8001). A putative BCCR was found at c.10–c.2176 (Table 2; Figure 3, orange arch), explained by a relative increase in breast cancer risk and a relative decrease in ovarian cancer risk (ROR = 3.47, 95% CI = 1.51–7.96, FDR-corrected ); it is also partially overlapped with previously reported BCCR (c.1–c.596, c.772–c.1806) [8] (Figure 3, orange arcs). The putative BCCR spans from the 5’ end of the CDS to approximately the boundary of exons 10 and 11.

3.3. Estimated Penetrance of Chinese Founder Mutations Greatly Differ

Of the six founder mutations examined, the two most prevalent mutations (BRCA1: c.5470_5477del and BRCA1: c.981_982del) were estimated to have penetrance lower than 0.35; one (BRCA2: c.3109C > T) with either above or lower than 0.35 penetrance subjected to family history (high penetrance only occurs in carriers with positive family history), and three with small sample sizes (BRCA1: c.3342_3345del, BRCA1: c.5154G > A, and BRCA1: c.4801A > T) showed complete penetrance (100%) due to no carriers found in control. Our results demonstrated that the penetrance of different BRCA1/2 deleterious mutations greatly varies and shows a large deviation from the expected value.

3.4. The Founder BRCA1: c.5470_5477del Is Estimated to Have Emerged More than 2000 Years Ago

As mentioned above, the BRCA1: c.5470_5477del was the most recurrent deleterious mutation, accounting for 9.5% of BRCA1 mutation carriers. Haplotype analysis was carried out on 31 unrelated patients and 50 unrelated controls without the mutation. The haplotype analysis was performed independent of the work of Meng et al. [6]; similar to their findings, our haplotype analysis suggested a strong founder effect (Supplementary Table S1) of the mutation. Moreover, carriers of this variant are distributed throughout the country (Figure 4), except provinces with sampling size <100 (regions colored in gray; total sample size: 7,919). Compared with other known Chinese founder mutations, BRCA1: c.5470_5477del (orange dots) showed the highest allele frequency and the most thorough spread in terms of geographical location, indicating a relatively early emergence of the mutation. We used DMLE + 2.3 to estimate the distribution of possible mutation age (years) under each growth rate (Supplementary Figure S2), and then, we compared the estimated mutation age (peak) and the actual time from which the growth rate is drawn. Among the seven growth rates used, the estimated mutation age () and the actual time (n) showed the best correlation with 1.54‰ (2 A.D-present), for which = 2,090 and n = 2,018. The estimated time (∼2,090 years ago, i.e., ∼70 B.C.) of the emergence of the mutation lies within the period of the Han dynasty (206 B.C–220 A.D.), currently known as the second imperial dynasty of China.

4. Discussion

We have presented a large-scale BRCA1/2 screening in both cancer patients and cancer-free individuals throughout China. We preliminarily defined the Chinese OCCR/BCCRs in BRCA1/2 and compared them with the previously reported OCCR/BCCRs based on mainly Caucasian population. We estimated penetrance for each of the six known Chinese founder mutations and demonstrated great variations between them. Unlike previous studies that were based on patients and/or their family members, this study included a large number of samples from the normal population, allowing us to assess the risk of a mutation carrier without a family history, which can significantly differ from the risk of those with family history. We also performed haplotype analysis on the most recurrent founder mutation BRCA1: c.5470_5477del and estimated its time of emergence to be ∼2,090 years ago within the Han dynasty.

Our OCCR/BCCRs partially overlap with previously reported OCCR/BCCRs drawn from the White/Jewish (Ashkenazi) population (Figures 2 and 3) [8]. We observed both in the Chinese and the White/Jewish studies the OCCRs’ position at the middle of the gene where the largest exon with the highest mutation rate is located and with binding sites for key proteins involved in DNA repair processes such as RAD51 and PALB2 [22, 23], whereas BCCRs tend to position at the ends (5’ and/or 3’) of the gene where smaller exons with secondary peaks of mutation rate are located and include transcription activation domain (TAD) and DNA-binding domain (DBD). OCCRs tend to be longer than BCCRs. The location of OCCR/BCCRs is thought to be mainly determined by two factors: the biological impact caused by mutations within a certain region and the frequency of these mutations detected within the studied population. While the former is likely not to be affected by ethnicity and may explain the concordant parts of the OCCR/BCCRs, the latter is known to be considerably affected by ethnicity and may explain the parts with discordance (38% of Chinese BRCA1/2 variants have not been reported in other populations [24]). Indeed, the Caucasian OCCR/BCCRs and Chinese OCCR/BCCRs show concordance with hotspot exons (i.e., frequently mutated exons) drawn from BIC and dbBRCA-Chinese databases [24], respectively. The difference in mutation frequency of each exon between the two databases may explain some of the differences shown between our Chinese OCCR/BCCRs and the previously reported Caucasian-based OCCR/BCCRs. For example, the reported mutation rate of BRCA1 Exon5 in BIC is threefold of that in the dbBRCA-Chinese database, which may explain the presence of a BCCR in Caucasian data and its absence in our Chinese data (note that due to historical reasons, the first exon was named exon 2 so that exon 11 is the largest exon, as seen in many literature and databases and that exons 2, 5, and 11 are equal to exons 1, 4, and 10 in our study). However, the exact mechanism of how the biological impact is affected by the position of the mutation is unclear, since the vast majority of BRCA2 mutations are truncating and able to trigger nonsense-mediated mRNA decay (NMD) [25]; truncated BRCA2 proteins are cytoplasmic [26] and unable to enter the nucleus due to lack of an intact nuclear localization signal. One possible explanation might be that NMD is not removing all truncating transcripts, which translate into truncated proteins and might compete with intact proteins for binding partners at the cytoplasm, or even become able to enter the nucleus through the carriage by a binding partner.

With the awareness of the great heterogeneity in cancer risks for mutations in different regions as evidenced by the OCCR/BCCRs, we focused on estimating independent penetrance for each recurrent mutation instead of a single penetrance for the whole gene. The two most prevalent mutations, which have been shown to spread through the country, are estimated to have relatively low penetrance, while those less prevalent ones residing in a local scope tend to have higher or complete penetrance. Our results demonstrated the “one risk does not fit all” as previously suggested by De Bock et al. [27]. Not every BRCA1/2 carrier has the same risk of developing cancers; some may never develop cancer throughout their lives. It is, therefore, important to separately consider when assessing cancer risks for BRCA1/2 carriers with different genotypes. For carriers of low-risk genotypes, regular follow-up surveillance is required. Follow-up screening for breast tumors is based on clinical palpation combined with imaging examination (MRI examination is preferred if available or ultrasound combined with mammography), while screening for ovarian cancer includes CA125 examination and ultrasound examination [28]. A more preventive strategy such as risk-reducing resection should be considered for carriers with high-risk genotype [29]. In BRCA mutation carriers, risk-reducing salpingo-oophorectomy (RRSO) alone or combined with bilateral risk-reducing mastectomy (BRRM) significantly reduces all-cause mortality, risk of ovarian cancer, and risk of breast cancer [30, 31]. When is the optimal time to perform RRSO? It is recommended to be performed at the age of 35–40 years after completion of childbearing for those with the BRCA1 mutation and at the age of 40–45 years for those with BRCA2 mutation [32]. In addition, complications associated with prophylactic surgical resection should not be ignored. Among them, the occurrence of early menopause is the most feared by both patient and physician after RRSO. Therefore, before RRSO is performed, patients should be informed of the common sequelae of medically induced menopause and informed of the benefits and risks of the appropriate remedies. Recent studies have found that HRT seems to be a safe therapeutic option in BRCA 1 and BRCA 2 mutation carriers undergoing RRSO and that the risk of BC treated with HRT after receiving RRSO was not significantly elevated [33]. Furthermore, assisted reproductive technology [34] should be considered for carriers who are of childbearing age and with high-risk genotype to prevent the passage of the high-risk allele. Women at high risk have a relatively short window of “reproductive age” and if they have not yet completed childbearing and have to consider RRSO, assisted reproduction with oocyte and/or embryo cryopreservation is recommended [35]. Therefore, the management of high-risk women with BRCA germline mutations requires the multidisciplinary involvement of geneticists, gynecologic oncologists, breast surgeons, and psychologists. For those putative high-risk mutations that seemed to be limited within a relatively small local area, it is necessary to carry out concentrated screening in order to further reveal the true penetrance in a larger sample size with standard incidence ratio (SIR), and if the high penetrance is confirmed, to identify more carriers of the mutation and take early preventative measures.

We estimated the time of emergence of the currently most recurrent Chinese founder mutation (also the most widely spread and probably the oldest), BRCA1: c.5470_5477del. Our estimation indicated that the emergence of this mutation may have happened ∼2,090 years ago (∼70 B.C.) during the Han dynasty, or more specifically, during the Western Han (Xi Han, 206 B.C.-8 A.D.). The thorough spread of the mutant allele throughout the country is most likely to be explained by multiple large-scale population migration events caused by frequent wars. There have been three major waves of population migration in Chinese history: the first started from the end of the Han dynasty (220 A.D.), i.e., the start of the Three Kingdoms period when the country was divided into three and suffered from long-lasting wars, until the end of the Southern and Northern dynasties (also known as the start of the Sui dynasty, 589 A.D.); the second occurred during the An Shi Rebellion in Tang dynasty (755 A.D.); the third happened during the Jingkang Incident (1127 A.D.), which has led to the end of the Northern Song dynasty. The three major migrations all happened after the estimated time of emergence of BRCA1: c.5470_5477del. There are 290 years (∼14.5 generation) between the estimated emergence (70 B.C.) and the first wave of migration (220 A.D.), which would be enough for the initial accumulation of mutant alleles so that the founder mutation can survive wars and natural disasters and inherit for more than two thousand years.

There are several limitations of this study. It is retrospective, and the number of subjects may be insufficient for a comprehensive estimation of the 1.4 billion Chinese population; follow-up time for BRCA1/2 carriers is short; most of the cancer-free carriers are under 70 years old, which can cause underestimation of the cancer risk; the portion of BRCA1/2 carriers in cancer-free individuals is higher than expected, which is likely due to the fact that those with family histories are more willing to participate in testing; we did not include large rearrangement events of BRCA1/2 in our research, which may account for more than 5% of all BRCA1/2 mutations [36]; due to the lack of genome-wide single-nucleotide polymorphism of all subjects included, we have no access to the population stratification considering the geographical difference and 56 ethnic origins of the Chinese population.

5. Conclusions

We preliminarily defined the OCCR/BCCRs based on a large number of Chinese samples. The Chinese OCCR/BCCRs partially overlap with the previously defined OCCR/BCCRs based on Caucasian samples. We estimated the penetrance of the six major Chinese founder mutations, respectively, and demonstrated great variations between them, which strongly suggests that cancer risks should be calculated and separately considered depending on the genotype rather than looking at a fixed risk figure. From a rigorous perspective, it is still necessary to calculate the true penetrance in a larger study cohort with population-based standardized incidences of breast and ovarian cancers. Finally, we investigated the most prevalent and nationally spread Chinese founder mutation BRCA1: c.5470_5477del and estimated that it has more than two thousand years of history.

Data Availability

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Additional Points

This manuscript has previously been released as a preprint at https://www.medrxiv.org/content/10.1101/2020.07.20.20135202v1 [37].

Ethical Approval

All patients provided informed consent, and the study was approved by the Institutional Review Board (IRB) of Daping Hospital, Army Medical University (IRB No. 2020(100)). The study protocol adhered to the principles of the Declaration of Helsinki.

Conflicts of Interest

Authors Ang Li, Jiaqi Luo, Yi Zi, Zhaoji Lan, Tianliangwen Zhou, Qihuan Zhi, Jiamin Zhan, Gang Sun, and Yujian Shi were employed by the company Top Gene Tech (Guangzhou) Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest. In addition, Ang Li, Yi Zi, Jiamin Zhan, and Yujian Shi were employed by Top Gene Tech (Guangzhou) Co., Ltd., when the work was conducted but no longer by the time of manuscript submission.

Authors’ Contributions

All authors contributed to the study conception and design, material preparation, data collection, and analysis. The draft of the manuscript was written by Ang Li, Shuai Hao, Jiaqi Luo, and Yi Zi, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. Ang Li, Shuai Hao, Jiaqi Luo, and Yi Zi equally contributed to this work.

Acknowledgments

The authors would like to thank all participants of this study for their contributions to scientific research. The authors also thank Dr. Xiaobin You, Yangming Wu, Huijun Su, Mei Zhu, Siwen Xu, Yun Gao, Zaixuan Zhong, Ling Xie, Yuanqin Wang, Qiuping Lin, Xiaoting Li, Hui Weng, Dan Li, and Shulan Xu’s suggestions about the manuscript. This study was funded by the Chongqing Science and Technology Bureau (no. 2021MSXM338).

Supplementary Materials

Table S1: most likely haplotype analysis of BRCA1 c.5470 5477delATTGGGCA mutation patient carriers. Figure S1: genetic map of the nine STR markers used in the haplotype analysis of BRCA1: c.5470_5477del (p.Ile1824AspfsTer3). Figure S2: age estimation of BRCA1: c.5470_5477del. Seven different growth rates (gr, %) were used together with the haplotype profiles as input to DMLE + 2.3, and a distribution of mutation age vs probability (frequency) was estimated for each simulation. The highest probability was achieved at gr = 0.154%, for which the estimated mutation age was 2090 years. (Supplementary Materials)

References

E. Levy-Lahad, R. Catane, S. Eisenberg et al., “Founder BRCA1 and BRCA2 mutations in Ashkenazi Jews in Israel: frequency and differential penetrance in ovarian cancer and in breast-ovarian cancer families,” The American Journal of Human Genetics, vol. 60, pp. 1059–1067, 1997.
View at: Google Scholar
U.-S. Khoo, K. Y. K. Chan, A. N. Y. Cheung et al., “Recurrent BRCA1 and BRCA2 germline mutations in ovarian cancer: a founder mutation of BRCA1 identified in the Chinese population,” Human Mutation, vol. 19, no. 3, pp. 307-308, 2002.
View at: Publisher Site | Google Scholar
A. Kwong, L. P. Wong, H. N. Wong et al., “A BRCA2 founder mutation and seven novel deleterious BRCA mutations in southern Chinese women with breast and ovarian cancer,” Breast Cancer Research and Treatment, vol. 117, no. 3, pp. 683–686, 2009.
View at: Publisher Site | Google Scholar
A. Kwong, E. K. Ng, C. L. Wong et al., “Identification of BRCA1/2 founder mutations in Southern Chinese breast cancer patients using gene sequencing and high resolution DNA melting analysis,” PLoS One, vol. 7, Article ID e43994, 2012.
View at: Publisher Site | Google Scholar
W.-M. Cao, Y. Gao, H.-J. Yang et al., “Novel germline mutations and unclassified variants of BRCA1 and BRCA2 genes in Chinese women with familial breast/ovarian cancer,” BMC Cancer, vol. 16, no. 1, p. 64, 2016.
View at: Publisher Site | Google Scholar
H. Meng, L. Yao, H. Yuan et al., “BRCA1 c.5470_5477del, a founder mutation in Chinese Han breast cancer patients,” International Journal of Cancer, vol. 146, no. 11, pp. 3044–3052, 2020.
View at: Publisher Site | Google Scholar
D. Gallardo-Rincón, R. M. Álvarez-Gómez, E. Montes-Servín et al., “Clinical evaluation of BRCA1/2 mutation in Mexican ovarian cancer patients,” Translational Oncology, vol. 13, no. 2, pp. 212–220, 2020.
View at: Publisher Site | Google Scholar
T. R. Rebbeck, N. Mitra, F. Wan et al., “Association of type and location of BRCA1 and BRCA2 mutations with risk of breast and ovarian cancer,” JAMA, vol. 313, pp. 1347–1361, 2015.
View at: Publisher Site | Google Scholar
K. B. Kuchenbaecker, J. L. Hopper, D. R. Barnes et al., “Risks of breast, ovarian, and contralateral breast cancer for BRCA1 and BRCA2 mutation carriers,” JAMA, vol. 317, no. 23, pp. 2402–2416, 2017.
View at: Publisher Site | Google Scholar
R. Ferla, V. Calo, S. Cascio et al., “Founder mutations in BRCA1 and BRCA2 genes,” Annals of Oncology, vol. 18, 2007.
View at: Publisher Site | Google Scholar
A. Li, R. Xie, Q. Zhi et al., “BRCA germline mutations in an unselected nationwide cohort of Chinese patients with ovarian cancer and healthy controls,” Gynecologic Oncology, vol. 151, no. 1, pp. 145–152, 2018.
View at: Publisher Site | Google Scholar
M. Krzywinski, J. Schein, İ. Birol et al., “Circos: an information aesthetic for comparative genomics,” Genome Research, vol. 19, no. 9, pp. 1639–1645, 2009.
View at: Publisher Site | Google Scholar
G. Cini, M. Mezzavilla, L. Della Puppa et al., “Tracking of the origin of recurrent mutations of the BRCA1 and BRCA2 genes in the north-east of Italy and improved mutation analysis strategy,” BMC Medical Genetics, vol. 17, no. 1, 2016.
View at: Publisher Site | Google Scholar
S. L. Neuhausen, S. Mazoyer, L. Friedman et al., “Haplotype and phenotype analysis of six recurrent BRCA1 mutations in 61 families: results of an international study,” The American Journal of Human Genetics, vol. 58, pp. 271–280, 1996.
View at: Google Scholar
M. Stephens, N. J. Smith, and P. Donnelly, “A new statistical method for haplotype reconstruction from population data,” The American Journal of Human Genetics, vol. 68, no. 4, pp. 978–989, 2001.
View at: Publisher Site | Google Scholar
J. P. Reeve and B. Rannala, “DMLE+: Bayesian linkage disequilibrium gene mapping,” Bioinformatics, vol. 18, no. 6, pp. 894-895, 2002.
View at: Publisher Site | Google Scholar
Z. G. Wu, J. Y. Zhou, X. D. Feng, Y. R. Fang, Y. S. Li, and D. G. Wu, Dictionary of Population Science, Southwest University of Finance and Economics Press, Chengdu, China, 1997.
E. V. Minikel, S. M. Vallabh, M. Lek et al., “Quantifying prion disease penetrance using large population control cohorts,” Science Translational Medicine, vol. 8, Article ID 322ra9ra9, 2016.
View at: Publisher Site | Google Scholar
T. Y. Wu, C. Wu, S. L. Chen, J. G. Shi, and K. H. Lin, “Breast cancer risk assessment and gail model among Chinese women,” Journal of Epidemiology and Preventive Medicine, vol. 1, no. 1, p. 107, 2015.
View at: Publisher Site | Google Scholar
M. S. Brose, T. R. Rebbeck, K. A. Calzone, J. E. Stopfer, K. L. Nathanson, and B. L. Weber, “Cancer risk estimates for BRCA1 mutation carriers identified in a risk evaluation program,” Journal of National Cancer Institute, vol. 94, no. 18, pp. 1365–1372, 2002.
View at: Publisher Site | Google Scholar
L. Yao, J. Sun, J. Zhang et al., “Breast cancer risk in Chinese women with BRCA1 or BRCA2 mutations,” Breast Cancer Research and Treatment, vol. 156, no. 3, pp. 441–445, 2016.
View at: Publisher Site | Google Scholar
S. L. Clark, A. M. Rodriguez, R. R. Snyder, G. D. Hankins, and D. Boehning, “Structure-function of the tumor suppressor BRCA1,” Computational and Structural Biotechnology Journal, vol. 1, Article ID e201204005, 2012.
View at: Publisher Site | Google Scholar
J. S. Martinez, C. Baldeyron, and A. Carreira, “Molding BRCA2 function through its interacting partners,” Cell Cycle, vol. 14, no. 21, pp. 3389–3395, 2015.
View at: Publisher Site | Google Scholar
S. P. Bhaskaran, K. Chandratre, H. Gupta et al., “Germline variation in BRCA1/2 is highly ethnic‐specific: evidence from over 30,000 Chinese hereditary breast and ovarian cancer patients,” International Journal of Cancer, vol. 145, no. 4, pp. 962–973, 2019.
View at: Publisher Site | Google Scholar
M. D. Ware, D. DeSilva, O. M. Sinilnikova, D. Stoppa-Lyonnet, S. V. Tavtigian, and S. Mazoyer, “Does nonsense-mediated mRNA decay explain the ovarian cancer cluster region of the BRCA2 gene?” Oncogene, vol. 25, no. 2, pp. 323–328, 2006.
View at: Publisher Site | Google Scholar
B. H. Spain, C. J. Larson, L. S. Shihabuddin, F. H. Gage, and I. M. Verma, “Truncated BRCA2 is cytoplasmic: implications for cancer-linked mutations,” Proceedings of the National Academy of Sciences, vol. 96, no. 24, pp. 13920–13925, 1999.
View at: Publisher Site | Google Scholar
G. H. De Bock, M. J. E. Mourits, and J. C. Oosterwijk, “One risk fits all?” Journal of Clinical Oncology, vol. 25, no. 22, pp. 3383-3384, 2007.
View at: Publisher Site | Google Scholar
E. Warner, K. Hill, P. Causer et al., “Prospective study of breast cancer incidence in women with a BRCA1 or BRCA2 mutation under surveillance with and without magnetic resonance imaging,” Journal of Clinical Oncology, vol. 29, no. 13, pp. 1664–1669, 2011.
View at: Publisher Site | Google Scholar
H. W. Tuffaha, A. Mitchell, R. L. Ward et al., “Cost-effectiveness analysis of germ-line BRCA testing in women with breast cancer and cascade testing in family members of mutation carriers,” Genetics in Medicine, vol. 20, no. 9, pp. 985–994, 2018.
View at: Publisher Site | Google Scholar
S. M. Domchek, T. M. Friebel, C. F. Singer et al., “Association of risk-reducing surgery in BRCA1 or BRCA2 mutation carriers with cancer risk and mortality,” JAMA, vol. 304, no. 9, pp. 967–975, 2010.
View at: Publisher Site | Google Scholar
C. Marchetti, F. De Felice, I. Palaia et al., “Risk-reducing salpingo-oophorectomy: a meta-analysis on impact on ovarian cancer risk and all cause mortality in BRCA 1 and BRCA 2 mutation carriers,” BMC Women’s Health, vol. 14, no. 1, p. 150, 2014.
View at: Publisher Site | Google Scholar
F. D. Felica, C. Marchetti, S. M. Boccia et al., “Risk-reducing salpingo-oophorectomy in BRCA1 and BRCA2 mutated patients: an evidence-based approach on what women should know,” Cancer Treatment Reviews, vol. 61, pp. 1–5, 2017.
View at: Google Scholar
C. Marchetti, F. De Felice, S. Boccia et al., “Hormone replacement therapy after prophylactic risk-reducing salpingo-oophorectomy and breast cancer risk in BRCA1 and BRCA2 mutation carriers: a meta-analysis,” Critical Reviews in Oncology, vol. 132, pp. 111–115, 2018.
View at: Publisher Site | Google Scholar
Q. Wang, J. F. Chow, W. S. Yeung et al., “Preimplantation genetic diagnosis using combined strategies on a breast cancer patient with a novel genomic deletion in BRCA2,” Journal of Assisted Reproduction and Genetics, vol. 31, no. 12, pp. 1719–1726, 2014.
View at: Publisher Site | Google Scholar
G. Corrado, C. Marchetti, R. Trozzi, G. Scambia, and A. Fagotti, “Fertility preservation in patients with BRCA mutations or lynch syndrome,” International Journal of Gynecological Cancer, vol. 31, no. 3, pp. 332–338, 2021.
View at: Publisher Site | Google Scholar
A. Kwong, J. Chen, V. Y. Shin et al., “The importance of analysis of long-range rearrangement of BRCA1 and BRCA2 in genetic diagnosis of familial breast cancer,” Cancer Genetics, vol. 208, no. 9, pp. 448–454, 2015.
View at: Publisher Site | Google Scholar
A. Li, Y. Zi, J. Luo et al., “Assessing the variations in breast/ovarian cancer risk for Chinese BRCA1/2 carriers.medRxiv,” 2020, https://www.medrxiv.org/content/10.1101/2020.07.20.20135202v1.
View at: Google Scholar

Copyright

Copyright © 2022 Ang Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

493

Downloads

526

Citations