Skip to main content

Advertisement

Comprehensive mutation detection of BRCA1/2 genes reveals large genomic rearrangements contribute to hereditary breast and ovarian cancer in Chinese women

Article metrics

Abstract

Background

Mutated BRCA1/2 genes are associated with hereditary breast and ovarian cancer (HBOC). So far most of the identified BRCA1/2 pathogenic variants are single nucleotide variants (SNVs) or insertions/deletions (Indels). However, large genomic rearrangements (LGRs) such as copy number variants (CNVs) are also playing an important role in HBOC predisposition. Their frequency and spectrum have been well studied in western populations but remain largely unknown for Chinese population.

Methods

Peripheral blood samples were collected from 218 unrelated familial breast and/or ovarian cancer (FBOC) patients living in Eastern China. PCR-based Sanger sequencing and panel-based next-generation sequencing (NGS) were performed to detect pathogenic SNVs and Indels in BRCA1/2 genes. For the patients lacking small pathogenic variants, multiplex ligation dependent probe amplification (MLPA) assay was conducted to screen for LGRs.

Results

In total, we identified 44 samples (20.1%) carrying small pathogenic variants (26 in BRCA1 and 18 in BRCA2, respectively). Among the rest of 174 samples, five were found carrying novel deleterious LGRs in BRCA1 which are exon5-7dup (1 patient), exon13-14dup (2 patients), and exon1-22del (2 patients). No LGR was found in BRCA2. Overall, LGRs accounted for 16.1% (5/31) of BRCA1 pathogenic variants, and were detected in 2.3% (5/218) of all FBOC patients.,

Conclusions

LGR variants in BRCA1 gene play a significant role in Chinese HBOC patients. MLPA or other similar LGR-detecting methods should be recommended along with nucleotide sequencing as the initial screening approach for Chinese HBOC women.

Background

According to National Central Cancer Registry of China, breast cancer ranks No.1 in cancer incidence and sixth in cancer-associated death for Chinese women, with over 250,000 newly diagnosed cases and 70,000 breast cancer-associated death in 2015 [1]. The average onset age of breast cancer is 45–55 years old for Chinese women, which is also younger than observed for Caucasian women [2]. While majority of breast cancer cases are sporadic, patients with familial history or other risk factors such as early onset age have been frequently observed in clinic, suggesting an important role of genetic factors in the disease development. Indeed, germline pathogenic variants in the two major breast cancer susceptibility genes BRCA1/2 have been detected within Chinese patients [3,4,5,6,7,8,9].

Studying BRCA1/2 pathogenic variants requires accurate and comprehensive testing methods. Short-read DNA sequencing methods, including both Sanger and next-generation sequencing (NGS), are only capable of reliably detecting small variants such as single nucleotide variants (SNVs) or insertion/ deletion (Indels), but not suitable for detecting large genomic rearrangements (LGRs), which involve deletions or duplications of multiple exons [i.e. copy number variants (CNVs)]. Therefore, sequencing alone may lead to underestimated frequency of pathogenic variants. Southern blotting could be used to detect LGRs [10], but is labor intensive and generally low-throughput. SNP or CGH arrays can detect copy number variants but their unit cost is high and resolution is usually over hundreds of Kb. Several multiplex PCR-based techniques have been recently developed to achieve higher processing and cost efficiency. For instance, multiplex ligation dependent probe amplification (MLPA) assay and multiplex amplicon quantification (MAQ) have been developed as fast and reproducible methods for CNV detection [11]. At the present, MLPA remains to be the most commonly used method for LGRs, and has detected 82.7 and 53% LGRs in BRCA1 and BRCA2, respectively [12].

As of today BRCA1/2 LGR studies have been mostly conducted in western countries, showing different prevalence with ethnicity and geography. For example, there was no BRCA1/2 LGR variants detected in Ashkenazi Jewish familial breast cancer patients [13, 14], but in non-Ashkenazi Jewish, the frequency of LGRs was 6% [14]. Very limited research has been conducted for Chinese, and only 12 BRCA1/2 LGRs have been so far reported. Those studies were conducted in Hong Kong [15], Singapore [16] and Malaysian [17]. The frequency and spectrum of BRCA1/2 LGRs in familial breast cancer patients from China mainland remain largely unknown.

Methods

Patient subjects

A total of 218 unrelated familial breast cancer patients were enrolled into this study between 2008 and 2017. All patients were diagnosed in Zhejiang Cancer Hospital in Eastern China and had a family history of at least one first- or second-degree relatives affected with breast cancer and/or ovarian cancer, regardless of age. Peripheral blood samples from the patients were collected in EDTA tubes and stored at − 80 °C. SNVs and Indels variants of BRCA1/2 were firstly determined for all patients using sequencing methods (PCR-based Sanger sequencing and panel-based NGS). The patients with negative sequencing finding were further screened for LGRs by MLPA. The written informed consents were obtained from all participating patients prior to clinical data and peripheral blood collection. This study was approved by the Research and Ethical Committee of Zhejiang Cancer Hospital, China. All experiments were performed in accordance with the approved guidelines.

DNA extraction

Genomic DNA was extracted from peripheral blood samples using QIAamp DNA Blood Mini kit (Qiagen, Hilden, Germany) by following the manufacture’s manual. DNA purity and concentration were measured by NanoDrop 2000 Spectrophotometer and Qubit 3.0 (Thermo Fisher Scientific, Waltham, USA), and DNA integrity was determined by agrose gel electrophoresis.

Nucleotide sequencing

The present study used both PCR-based Sanger sequencing and panel-based NGS to interrogate small nucleotide variants including SNVs and Indels. In the first phase of the project, Sanger sequencing was performed on 133 unrelated FBOC cases using a total of 72 pairs of oligos to cover all coding exons and intron-exon boundaries of BRCA1/2. The primer oligo sequences were listed in Additional file 1: Table S1 and Additional file 2: Table S2. In the second phase of the project, in order to achieve high-throughput and cost-effective sequencing, we designed a NGS panel by adopting the NEBNext Direct sequencing technology developed by New England Biolabs (Ipswich, MA). The panel contains BRCA1/2 genes as well as other 96 known cancer risk-associated genes. We performed panel NGS on all of the 133 Sanger cases along with 85 new cases newly collected. Individually prepared libraries were pooled for Hiseq X sequencing (Illumina, CA, USA) to achieve a minimum 500x mean coverage for the included panel genes. Raw FASTQ data run through in house bioinformatic pipeline with variant calling generated for BRCA1/2 genes. Variant filtering and final interpretation were conducted by following the ACMG Standards and Guideline for the Interpretation of Sequence Variants [18] and based on a set of criteria such as allele frequency as well as information from clinical genome databases including ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/), Online Mendelian Inheritance in Man (OMIM) (http://www.omim.org/) and Human Gene and Mutation Database (HGMD) (http://www.hgmd.cf.ac.uk/ac/index.php).

LGRs analysis

BRCA1/2 LGRs was screened by (Multiplex ligation dependent probe assay) MLPA assay using the SALSA P002 kit and P045 kit for BRCA1 and BRCA2 genes, respectively (MRC-Holland, Amsterdam, the Netherlands). The MLPA reactions were performed according to the manufacturer’s instruction. Five normal control samples were included as reference within each MLPA run. Fragment analysis of the PCR products were performed on an ABI 3130xl Genetic Analyzer (Applied Biosystems, Foster City, CA). The data was analyzed by using the Coffalyser software v.9 (Applied Biosystems, Foster City, CA). All of the peak heights were normalized, and the ratio value between 0.7–1.3 was considered as normal. A ratio value ≤0.7 or ≥ 1.3 was threshold suggestive of a deletion or duplication, respectively. All patients with a value ≤0.7 or ≥ 1.3 were confirmed by independent experiments.

Two primer oligos were designed to validate BRCA1 Exon 5–7 duplication. The forward primer sequence was CCGTGCCAAAAGACTTCTACA (Exon 7) and the reverse primer sequence was TTGCTTCCAACCTAGCATCA (Exon 5). Long range PCR amplification was performed with Takara LA Taq DNA polymerase (Takara Bio, USA) by following the manufacturer’s manual. The amplified product was run on 0.8% Agrose gel electrophoresis with EB (i.e. ethidium bromide) and visualized under UV light. The purified amplicons were subjected to Sanger sequencing to confirm amplification fidelity.

Results

Small pathogenic variants in BRCA1/2 genes

Overall, we identified a total of 31 BRCA1 or BRCA2 pathogenic SNVs and Indels in 44 unrelated patients by combining Sanger sequencing and the 98-gene panel NGS assay. Table 1 lists all these small variants. In summary, nearly 59% (26 of 44) of patients had BRCA1 pathogenic variants, and 41% (18 of 44) had BRCA2 pathogenic variants. Two recurrent pathogenic variants (c.5154G > A and c.5468-1_5474del GCAATTGG) in BRCA1 were reported as putative founder mutations [19]. In total, frequency of BRCA1/2 small pathogenic variants was 20.2% (44/218) in the studied cohort.

Table 1 Small pathogenic variants of BRCA1 and BRCA2 in 218 familial breast and/or ovarian cancer patients

Novel LGRs identified in BRCA1

Among the 174 patients lacking BRCA1/2 small pathogenic variants, three unique BRCA1 LGRs were detected in 5 (2.9%) cases by MLPA assay (Fig. 1a, b, c, d). These include one case with exon5-7dup, two cases with exon13-14dup, and two cases with exon1-22del (Table 2). To our knowledge, these three LGRs have not been reported in Chinese HBOC patients. To confirm MLPA results, we validated exon 5-7dup by designing oligo primers surrounding the putative junction. We obtained a clear and strong 6-8Kb PCR amplicon (Fig. 1e), whose sequence identity was confirmed by Sanger sequencing (data not shown) supporting a tandem duplication event. Overall, BRCA1 LGRs accounted for 16.1% (5/31) of all patients with BRCA1 pathogenic variants. No LGR was identified in BRCA2. Combining small nucleotide variants and LGRs, we obtained a frequency of 22.5% (49/218) for this cohort, with BRCA1 LGRs accounting for 10.2% (5/49) of all BRCA1/2 pathogenic variants.

Fig. 1
figure1

Three BRCA1 LGRs detected in the HBOC cohort. MLPA result for BRCA1 LGRs in bar chart format generated by Coffalyser software v.9. BRCA1 exons and intra normalized ratio are given on the X axis and Y axis, respectively. Exons having reduced or increased peak ratio are denoted by red or blue dot, respectively. a MLPA result for BRCA1 wildtype. b exon5-7dup. c exon13-14dup. d exon1-22del. e exon5-7dup was confirmed by LR-PCR

Table 2 BRCA1 LGRs in 174 familial breast and/or ovarian cancer patients

Disease pathology associated with BRCA1 LGRs

The characteristics and familial cancer history of patients with BRCA1 LGRs were listed in Table 2 and their family pedigrees were shown in Fig. 2. In the five breast cancer patients with BRCA1 LGRs, the most common tumor type was invasive ductal carcinoma. Moreover, Micropapillary carcinoma, mucinous carcinoma as well as medullary carcinoma was found in these patients with BRCA1 LGRs. Although triple negative (ER−/PR−/HER2-) subtype was the most common subtype, luminal subtype (ER+) and HER2-overexpression subtype (ER−/PR−/HER2+) were also exited in patients with BRCA1 LGRs.

Fig. 2
figure2

Pedigree of families with proband carrying BRCA1 LGRs. a 147th family with proband carrying exon5-7dup. b 10th family with proband carrying exon13-14dup. c 213th family with proband carrying exon13-14dup. d 113th family with proband carrying exon1-22del. e 203th family with proband carrying exon1-22del. (BC, breast cancer; OC, ovarian cancer; Bi-BC, bilateral breast cancer)

Discussion

In this study, we performed nucleotide sequencing for 218 unrelated FBOC patients living in Eastern China and observed a 20.2% overall pathogenic variant frequency for BRCA1/2 genes. MLPA assay on the patients lacking small pathogenic variants (174 of 218) further identified 3 unique LCRs in 5 patients, increasing the total BRCA1/2 pathogenic variant frequency to 22.5%. All of the three BRCA1 LGRs were not previously reported by Chinese population studies [15,16,17]. Interestingly, no LGR was identified in BRCA2 gene within our cohort, consistent to the knowledge that LGRs are more frequently observed in BRCA1 than BRCA2. It was revealed that Alu-mediated unequal homologous recombination could be the most common mechanism of LGRs found in BRCA1/2, as 72.84% (59/81) and 52.94% (9/17) LGRs in BRCA1 and BRCA2 respectively were mediated by this manner [20]. The reason behind higher LGRs frequency in BRCA1 than in BRCA2 might be due to the higher Alu density (41.5%) in the BRCA1 gene than in the BRCA2 gene [21].

It has been reported that frequency of LGRs ranges from approximately 6–27% of all detected BRCA1 pathogenic variants, and BRCA2 LGRs play a less role in hereditary breast cancer patients [20]. Thirty five out of three hundreds (12%) of BRCA1/2-sequencing negative familial breast cancer patients from non-Ashkenazi Jewish in US were found carrying BRCA LCRs, with 10% (31/300) LCRs detected in BRCA1 and 1% (4/300) detected in BRCA2, respectively [22]. A nationwide study conducted in South Korea showed that LGRs were detected in 3.7% (3/81) of patients bearing BRCA1/2 pathogenic variants and 7.5% (3/40) of patients bearing only BRCA1 pathogenic variants [23]. A large sample screening of high-risk breast cancer patients from Hong Kong showed that LGRs accounted for 6.67% (8/120) of all BRCA1/2 pathogenic variants, involving 8.77% (5/57) of BRCA1 and 4.76% (3/63) of BRCA2, respectively [15]. In our present study, BRCA1 LGRs account for 16.1% (5/31) of all BRCA1 pathogenic variants and 10.2% (5/49) of all BRCA1/2 pathogenic variants. For 174 cases with negative sequencing results, five (2.9%) were identified with BRCA1 LCRs.

Conclusions

To conclude, our study has provided evidence that BRCA1/2 LGRs contribute significantly to the development of HBOC in Chinese mainland population and LGRs screening should be taken into consideration in hereditary breast cancer consulting. It is imperative that the frequency and spectrum of BRCA1/2 should be investigated in the context of both small nucleotide variants and LGRs.

Availability of data and materials

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

Abbreviations

BC:

Breast cancer

IHC:

Immunohistochemistry

LGRs:

Large genomic rearrangements

MAQ:

Multiplex amplicon quantification

MLPA:

Multiplex ligation dependent probe amplification

NGS:

Next generation sequencing

OC:

Ovarian cancer

PCR:

Polymerase chain reaction

References

  1. 1.

    Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115–32.

  2. 2.

    Fan L, Strasser-Weippl K, Li JJ, St Louis J, Finkelstein DM, Yu KD, et al. Breast cancer in China. Lancet Oncol. 2014;15:e279–89.

  3. 3.

    Tang NL, Pang CP, Yeo W, Choy KW, Lam PK, Suen M, et al. Prevalence of mutations in the BRCA1 gene among Chinese patients with breast cancer. J Natl Cancer Inst. 1999;91:882–5.

  4. 4.

    Cao W, Wang X, Li JC. Hereditary breast cancer in the Han Chinese population. J Epidemiol. 2013;23:75–84.

  5. 5.

    Pang Z, Yao L, Zhang J, Ouyang T, Li J, Wang T, et al. RAD51C germline mutations in Chinese women with familial breast cancer. Breast Cancer Res Treat. 2011;129:1019–20.

  6. 6.

    Sun J, Wang Y, Xia Y, Xu Y, Ouyang T, Li J, et al. Mutations in RECQL gene are associated with predisposition to breast Cancer. PLoS Genet. 2015;11:e1005228.

  7. 7.

    Sun J, Meng H, Yao L, Lv M, Bai J, Zhang J, et al. Germline mutations in Cancer susceptibility genes in a large series of unselected breast Cancer patients. Clin Cancer Res. 2017;23:6113–9.

  8. 8.

    Lang GT, Shi JX, Hu X, Zhang CH, Shan L, Song CG, et al. The spectrum of BRCA mutations and characteristics of BRCA-associated breast cancers in China: screening of 2,991 patients and 1,043 controls by next-generation sequencing. Int J Cancer. 2017;141:129–42.

  9. 9.

    Li J, Jing R, Wei H, Wang M, Qi X, Liu H, et al. Germline mutations in 40 cancer susceptibility genes among Chinese patients with high hereditary risk breast cancer. Int J Cancer. 2018. https://doi.org/10.1002/ijc.31601.

  10. 10.

    Puget N, Torchard D, Serova-Sinilnikova OM, Lynch HT, Feunteun J, Lenoir GM, et al. A 1-kb Alu-mediated germ-line deletion removing BRCA1 exon 17. Cancer Res. 1997;57:828–31.

  11. 11.

    Concolino P, Mello E, Minucci A, Santonocito C, Scambia G, Giardina B, et al. Advanced tools for BRCA1/2 mutational screening: comparison between two methods for large genomic rearrangements (LGRs) detection. Clin Chem Lab Med. 2014;52:1119–27.

  12. 12.

    Ruiz De Garibay G, Gutiérrez-Enríquez S, Garre P, Bonache S, Romero A, Palomo L, et al. Characterization of four novel BRCA2 large genomic rearrangements in Spanish breast/ovarian cancer families: review of the literature, and reevaluation of the genetic mechanisms involved in their origin. Breast Cancer Res Treat. 2012;133:273–83.

  13. 13.

    Stadler ZK, Saloustros E, Hansen NA, Schluger AE, Kauff ND, Offit K, et al. Absence of genomic BRCA1 and BRCA2 rearrangements in Ashkenazi breast and ovarian cancer families. Breast Cancer Res Treat. 2010;123:581–5.

  14. 14.

    Palma MD, Domchek SM, Stopfer J, Erlichman J, Siegfried JD, Tigges-Cardwell J, et al. The relative contribution of point mutations and genomic rearrangements in BRCA1 and BRCA2 in high-risk breast cancer families. Cancer Res. 2008;68:7006–14.

  15. 15.

    Kwong A, Chen J, Shin VY, Ho JCW, Law FBF, Au CH, et al. The importance of analysis of long-range rearrangement of BRCA1 and BRCA2 in genetic diagnosis of familial breast cancer. Cancer Genet. 2015;208:448–54.

  16. 16.

    Lim YK, Iau PTC, Ali AB, Lee SC, Wong JE, Putti TC, et al. Identification of novel BRCA large genomic rearrangements in Singapore Asian breast and ovarian patients with cancer. Clin Genet. 2007;71:331–42.

  17. 17.

    Kang P, Mariapun S, Phuah SY, Lim LS, Liu J, Yoon SY, et al. Large BRCA1 and BRCA2 genomic rearrangements in Malaysian high risk breast-ovarian cancer families. Breast Cancer Res Treat. 2010;124:579–84.

  18. 18.

    Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.

  19. 19.

    Cao WM, Gao Y, Yang HJ, Xie SN, Ding XW, Pan ZW, et al. Novel germline mutations and unclassified variants of BRCA1 and BRCA2 genes in Chinese women with familial breast/ovarian cancer. BMC Cancer. 2016;16:64.

  20. 20.

    Sluiter MD, Van Rensburg EJ. Large genomic rearrangements of the BRCA1 and BRCA2 genes: review of the literature and report of a novel BRCA1 mutation. Breast Cancer Res Treat. 2011;125:325–49.

  21. 21.

    Ewald IP, Ribeiro PL, Palmero EI, Cossio SL, Giugliani R, Ashton-Prolla P. Genomic rearrangements in BRCA1 and BRCA2: a literature review. Genet Mol Biol. 2009;32:437–46.

  22. 22.

    Walsh T, Casadei S, Coats KH, Swisher E, Higgins J, Roach KC, et al. Spectrum of mutations in BRCA1, BRCA2, CHEK2, and TP53 in families at high risk of breast cancer. JAMA. 2006;295:1379–88.

  23. 23.

    Seong MW, Cho SI, Kim KH, Chung IY, Kang E, Lee JW, et al. A multi-institutional study of the prevalence of BRCA1 and BRCA2 large genomic rearrangements in familial breast cancer patients. BMC Cancer. 2014;14:645.

Download references

Acknowledgements

Not applicable.

Funding

This research was supported by the grants from National Natural Science Foundation of China (grant number: 81672597), Natural Science Foundation of Zhejiang Province, China (grant number: LY17H160038), Key Research-Development Program of Zhejiang Province (grant number: 2017C03013, 2019C04001), Science and Technology Program offered by Health Bureau of Zhejiang Province, China (2012RCB006 and 2014KYA006). The funding bodies had no role in study design, collection, analysis, or interpretation of data, or in writing the manuscript.

Author information

WMC: designed the study, analyzed the MLPA analysis data and drafted the manuscript. YBZ: performed the analysis of MLPA. YG and ZWP: extracted DNA and performed PCR. XWD, YS, YH and CJL: contributed samples and patient information. GP: reviewed the data and drafted the manuscript. XJW: conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Correspondence to Xiao-Jia Wang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Research and Ethical Committee of Zhejiang Cancer Hospital, China. The written informed consents were obtained from all of the participating patients prior to clinical data and peripheral blood collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Primers for entire coding exons and intron-exon boundaries of BRCA1. (DOCX 30 kb)

Additional file 2:

Table S2. Primers for entire coding exons and intron-exon boundaries of BRCA2. (DOCX 31 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Chinese
  • Familial breast cancer
  • Familial ovarian cancer
  • BRCA1
  • BRCA2
  • Rearrangement