Abstract
Background/Aim: The aim of this study was to analyse the genetic profiles of metastatic castration-resistant prostate cancer (mCRPC) by using next generation sequencing to identify variants with pathogenic potential in nine DNA repair genes – BRCA1, BRCA2, RAD50, RAD51, RAD51C/D, ATM and ATR. Materials and Methods: Isolated genomic DNA from peripheral blood of 50 patients with mCRPC was used for the sequencing of 111 genes associated with hereditary cancer on an Illumina platform. Identified variants were tested in Integrative Genomic Viewer, their clinical significance confirmed in databases and their potential impact on protein function predicted by in silico tools. Results: From nine analysed DNA repair genes, we identified 14 relevant variants; three pathogenic variants – BRCA2 (rs80359306), RAD50 (rs786201531) and ATM (rs1555099760), and eleven other variants with pathogenic potential. Conclusion: The pathogenic variants identified in this study are located in evolutionarily conserved regions of proteins and are highly likely to affect DNA repair efficiency.
Prostate cancer (PC) is the second most common oncological disease in Slovakia, one of the major medical problems of the male population occurring in ages over 50 years. Despite extensive studies of PC, the exact cause is not yet known.
Based on DNA change, certain cells in prostate start to grow and proliferate without control which mostly results in tumour formation. Upon DNA damage, DNA repair mechanisms are activated to repair damaged DNA and thus prevent cancer. Some DNA changes can either be inherited or acquired during a lifetime. According to American Cancer Society, more than 10% of PC cases represent hereditary cancer, in which inherited genes play an important role. It has long been known that heredity plays a significant role in PC. Men with a family history of PC have a greater risk of developing PC (1). Germline mutations in DNA repair genes are linked to poor prognosis, early onset, increased aggressiveness, drug treatment response, higher possibility of early development of CRPC and the course of the disease itself (2).
In our study, we focused on changes in DNA repair genes, that may greatly influence DNA repair mechanisms as well as genomic integrity. Deregulation of repair mechanisms can cause disruption of gene function, accelerated aging and higher risk for cancer (3). There are five major DNA repair pathways whose role is to repair and eliminate single-stranded (SSB) and double-stranded (DSB) DNA breaks. We focused on genes involved in homologous recombination (HR). Approximately 20-25% of mCRPC are mostly linked with defects in HR genes (4), implying the importance of HR as a major DNA damage repairing pathway, monitoring genomic stability and preventing cells from being transformed into cancerous cells (5).
The first activated complex involved in the recognition of DSB is the MRE11-RAD50-NBS1 (MRN) complex, which due to its DNA-binding ability adheres to DNA damage sites. DSB also activates the ATM and ATR kinases whose role lies in coordinating cell cycle checkpoint and DNA repair pathways. For the final homology-directed repair (HDR) of DNA breaks is crucial the formation of RAD51 nucleofilaments mediated by checkpoint regulation factors PALB2 and BRCA2, which cooperate with BRCA1 (6). The RAD51 nucleofilaments, along with paralogues RAD51C/D, mediate the homology search, prepare and use sister chromatid as a template for a new DNA strand synthesis (4).
Materials and Methods
The study group included 50 mCRPC patients with unknown records of the hereditary origin of prostate cancer [mean age=71.64, median age=71.5; mean prostate-specific antigen (PSA)=6.06 ng/ml, median PSA=19.29 ng/ml]. This study was approved by the Ethical Committee at Jessenius Faculty in Medicine in Martin, Comenius University in Bratislava.
Genomic DNA from buffy coat (approx. 200 μl) was purified using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). Quality and quantity of the DNA were measured on NanoDrop™2000 Spectrophotometer (Thermo Fisher Scientific, MA, USA). The Qubit® 3.0 Fluorometer (Thermo Fisher Scientific) was used for subsequently dilution of DNA samples to concentration 20 ng/μl required for the preparation of NGS library. NGS was performed by the external supplier Sistemas Genómicos (Valencia, Spain). The analysed panel GeneSGKit® contained 111 genes related to hereditary cancer of heterogeneous origin. Their protein-coding exons, flanking regions and introns were sequenced using HiSeq 2500 System (Illumina, CA, USA). Raw sequencing data were analysed through GeneSystems© (Sistemas Genómicos, Valencia, Spain) and a VCF file was created for every patient. Identified variants with potential clinical significance were tested in Integrative Genomic Viewer (IGV) and their clinical significance was confirmed in VarSome search engine (7). To predict the potential impact of amino acid substitution on protein function we used in silico predictor tools SIFT (8) and PolyPhen-2 (9). The source of protein sequence expressed from individual genes used for in silico predictions was UniProt (10). To evaluate indel variants for their disease-causing potential we used MutationTaster (11).
Results
We examined gene profiles of 50 patients with mCRPC for the presence and frequency of germline mutations in nine DNA-repair genes – BRCA1, BRCA2, PALB2, RAD50, RAD51 and its paralogues RAD51C and RAD51D, ATM and ATR. With respect to the reference genome (GRCh38.p7 assembly), individual genes contained multiple variants with mean value 17.44 variants per gene: missense 48 (30.6%), frameshift 5 (3.2%), in frame 1 (0.6%), synonymous 26 (16.6%) and intron 77 (49.0%). Average sequencing depth was 365 reads per gene (minimum 171 reads for ATM, maximum 579 reads for BRCA1). From all 157 variants, 14 variants with pathogenic potential according to American College of Medical Genetics and Genomics (ACMG) classification implemented into the automated criteria and rules in VarSome, were filtered. The Table I lists Human Genome Variation Society (HGVS) names and other characteristics of all filter passed variants in our patient sample, and in the Table II we present the characteristics of the pathogenic potential and clinical significance of these variants.
BRCA1. In BRCA1 gene, we detected 10 missense variants; eight in exon 10 (rs1799950, rs4986850, rs28897683, rs799917, rs16941, rs4986852, rs16942 and rs28897689) and two in exon 15 (rs1799966, rs1799967). The first nine variants mentioned do not lie in any of the important BRCA1 domains, but in mutational hot spots of 31 base-pairs length, which have been connected to early onset and high histological grade of cancer (12). The pathogenicity of variants in these areas varies between 54.9-81.0% and include up to 28 classified variants mostly with pathogenic status (7). Only the benign missense variant p.Met1652Ile (rs1799967) is located in BRCT1 domain, one of the integral signalling modules in DNA damage response (13). According to VarSome and ClinVar all BRCA1 variants identified in this study are benign or likely benign, deleterious effect (SIFT and PolyPhen-2) was predicted for only three variants, p.Thr826Lys, p.Ser1040Asn and p.Ser1613Cys, but MutationTaster evaluated them as frequent polymorphisms.
BRCA2. In BRCA2 gene, we analysed one pathogenic frameshift variant (rs80359306) in exon 10, localised in the PALB2-binding domain, which is essential for many functions, such as G2/M checkpoint, protection of replication fork and DNA breaks repair by homologous recombination (14). In exon 11, we detected three missense variants (rs41293503, rs1799944, rs4987117), which are located in the RAD51-binding domain at the N-terminus of the protein. RAD51 and BRCA2 provide exchange of genetic information between damaged and sister chromatids (15). In silico tools predicted that only the variant p.Glu2236Lys (rs41293503) may be associated with disease cause in the Human Gene Mutation Database (HGMD), while other databases cannot determine the exact clinical significance of this variant. In exon 14, we found another missense variant p.Gly2353Arg (rs80358935), located in DNA-binding domain, which is an integral part of homologous recombination facilitating the binding of the recombination protein RAD51 at DNA breaks (15). Regarding the frameshift indel variant (rs276174803) in exon 27 at position c.10095 with deletion of one cytosine and insertion of 11 bases (C/GAATTATATCT), the formation of a stop codon and subsequent termination of protein synthesis is predicted.
PALB2. In each of exons 8 and 9 of the PALB2 gene (partner and localizer of BRCA2), we found two benign missense variants: p.Val932Met (rs45624036) and p.Gly998Glu (rs45551636).
RAD50. In the coiled-coil domain of RAD50, we detected the pathogenic frameshift variant p.Ser181PhefsTer7 (rs786201531), which terminates the protein at position 188 and the missense variant of uncertain significance (VUS) p.Glu245Lys (rs797044753).
RAD51 and its paralogues. We did not find any potentially pathogenic variants in the RAD51 gene and its paralogue RAD51D. However, in exon 6 of the paralogue RAD51C, we found the likely benign missense variant p.Thr287Ala (rs28363317) that is located at the C-terminal Walker B motif ATPase domain.
ATM. In the HEAT repeat segment of ATM, we detected the pathogenic frameshift variant p.Met1484ArgfsTer15 (rs1555099760) that terminates ATM protein at position 1499 and thus significantly truncates the protein by 1557 amino acids and two benign missense variants, p.Phe858Leu (rs1800056) and p.Pro1054Arg (rs1800057). In the important regulatory FAT domain, we found one missense VUS p.Gly2023Arg (rs11212587). In the PI3K domain, we found the likely benign missense variant p.Leu2492Arg (rs56399857) with deleterious prediction in SIFT, PolyPhen-2 and also in MutationTaster.
ATR. We found two variants in the N-terminal HEAT repeat segment of ATR; the pathogenic frameshift variant p.Val327HisfsTer3 (rs759554602), which truncates the protein from the original 2644 amino acids to 330 and the missense VUS p.His1531Arg (rs200070057) with deleterious predictions by SIFT, PolyPhen-2 and also in MutationTaster.
Discussion
Despite the rapid development of new technologies that allow us to observe the profile of the human genome, the genes predisposing to prostate cancer are largely unknown. In this study, we examine profiles of DNA damage repair (DDR) genes in which defects are more prevalent than previously predicted. Recent insights showed that approx. 23% of mCRPC patients are carriers of somatic/germline mutations in DDR genes, such as BRCA1/2, PALB2, CHEK2, ATM (16) and FANCA and RAD51 paralogues (17). It has been suggested that DDR defects may be involved in the development of carcinogenesis, disease progression and higher aggressiveness (16, 18). In terms of treatment, for mCRPC patients with germline mutations in DDR genes, drugs like docetaxel, abiraterone, enzalutamide and others have been developed and rapidly improved survival, but there remains a need to develop therapies to improve clinical outcomes (16). Several studies have demonstrated that the hereditary component of PC is one of the strongest risk factors. Numerous germline mutations in BRCA1 and BRCA2 genes may increase the risk of prostate cancer (1). Approximately 5.3% of prostate cancer patients carrying germline mutations in BRCA2 are linked to poor prognosis and early onset (age <55 years), higher Gleason score (>8), elevated PSA, more aggressive phenotype and poor clinical outcomes (1, 19, 20).
In our study, we identified 14 variants having a pathogenic status in VarSome, but only three of them have pathogenic status in ClinVar. Each of the pathogenic frameshift variants is located in a different gene; p.Ile605AsnfsTer11 (BRCA2), p.Ser181PhefsTer7 (RAD50) and p.Met1484ArgfsTer15 (ATM). The frameshift variant p.Ile605AsnfsTer11 (rs80359306) was multiple times associated with familial prostate cancer (19, 21, 22), as well as with breast and ovarian cancer. Through the creation of a stop codon at position 616 and protein truncation, there has been predicted to cause loss of protein (23). Nykamp et al. have found that the truncating variant p.Ser181PhefsTer7 in RAD50 that we detected in heterozygous status, is lethal or severe (24). These variants may also have the potential to disrupt the communication between the enzymatic domains and functions of the highly conserved MRN protein complex. However, they have not been associated with RAD50-related disease. In exon 30 of ATM, we identified the pathogenic frameshift variant p.Met1484ArgfsTer15 (rs1555099760) that causes a premature truncation of ATM protein, which binds to the MRN complex (25).
In addition to these pathogenic variants, we identified other variants with the potential to affect protein functions in DNA repair mechanisms. The variants identified in genes BRCA1, RAD51 and RAD51D were not included in any results table because their pathogenicity was not demonstrated by the in silico tools or listed in databases. In BRCA2, we identified one frameshift indel variant in heterozygous status. It starts at position 10095 by deleting one cytosine and inserting 11 bases GAATTATATCT resulting in premature termination of translation in the very C-terminus of BRCA2 protein. Wojcik et al. have classified this variant as non-pathogenic, however, Balabas et al. and Ratajska et al. have suggested an association of the variant p.Ser3366AsnfsTer5 with familial aggregation of breast and/or ovarian cancer (26-28).
We identified two missense variants in PALB2; p.Gly998Glu and p.Val932Met. Both are located at the C-terminus of the WD40 domain, whose role is direct binding of RAD51C and BRCA2 in the process of HR. The protein PALB2 is indispensable for the HR pathway, and it is also considered as a linker between BRCA1 and BRCA2 (29).
There have not been many studies dealing with the association between germline mutations and prostate cancer in RAD51, but RAD51 and XRCC3 polymorphisms may be important for the development of PC which by affecting the DNA repair process and increasing the susceptibility to neoplastic transformation (30). We identified one missense variant p.Thr287Ala (rs28363317) in RAD51C associated with PC (31) with possible loss of protein function prediction, generally considered as non-pathogenic and likely benign (31-33). This variant is located at the C-terminal Walker B motif ATPase domain. Alterations in this domain can decrease the ability to repair DNA interstrand crosslinks by RAD51 paralogues and gene XRCC2 (34).
We identified two missense variants in ATM; the p.Gly2023Arg which involves the alteration of a conserved nucleotide with unknown clinical significance (35, 36), and the p.Leu2492Arg with a predicted deleterious effect (SIFT, PolyPhen-2) on the ATM protein henceforth referred to as a VUS (37). This variant is located in the PI3K domain responsible for the phosphorylation of BRCA1, p53 and other proteins (38).
In the ATR, which is considered to be a DNA damage sensor, we detected the p.Val327HisfsTer3 variant, which is likely a pathogenic frameshift that causes premature truncation of the protein and loss of the main protein function; and the missense variant p.His1531Arg with deleterious outcome and prediction of possible disruption of protein function. There are no records of these two variants in the literature.
To conclude, we identified multiple germline mutations in DNA repair genes involved in homologous recombination pathway. These variants may have the potential to reduce or cause protein dysfunction, and have been associated with more aggressive phenotype of PC and earlier onset of mCRPC. Detection of relevant pathogenic variants in these genes would lead to earlier diagnosis of PC and would allow clinicians to set earlier, more targeted treatment for patients with a family predisposition or patients in the first stages of PC.
Acknowledgements
This study was supported by the Agency of Ministry of Education, Science, Research and Sport of Slovakia by the Slovak Research and Development Agency under the contract no. APVV-15-0181 and by the Grant of the Comenius University no. UK/26/2019.
Footnotes
Authors' Contributions
KH and KB conceived the study and were in charge of overall direction and planning. KH analysed the data, performed literature reviews and drafted the manuscript in consultation with MH, LM, JK and MS. JK provided samples from Urology Clinic. All Authors discussed the results and contributed to the final manuscript.
This article is freely accessible online.
Conflicts of Interest
The Authors confirm that there are no conflicts of interest in regard to this study.
- Received March 19, 2020.
- Revision received March 30, 2020.
- Accepted March 31, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved