GT198 (PSMC3IP) germline variants in early-onset breast cancer patients from hereditary breast and ovarian cancer families

ABSTRACT

GT198, located 470 kb downstream of BRCA1, encodes for the nuclear PSMC3-interacting protein, which functions as co-activator of steroid hormone-mediated gene expression, and is involved in RAD51 and DMC1-mediated homologous recombination during DNA repair of double-strand breaks. Recently, germline variants in GT198 have been identified in hereditary breast and ovarian cancer (HBOC) patients, mainly in cases with early-onset. We screened a cohort of 166 BRCA1/2 mutation-negative HBOC patients, of which 56 developed early-onset breast cancer before the age of 36 years, for GT198 variants. We identified 7 novel or rare GT198 variants in 8 out of 166 index patients: c.-115G>A (rs191843707); c.-70T>A (rs752276800); c.-37A>T (rs199620968); c.-24C>G (rs200359709); c.519G>A p.(Trp173*); c.537+51G>C (rs375509656); c.*24G>A. Three out of 7 identified variants (c.-115G>A, c.519G>A and c.*24G>A) with putative pathogenic impact were found in HBOC patients with breast cancer onset at ≤ 36 years. The nonsense mutation c.519G>A p.(Trp173*) was located within the DNA binding domain of GT198 and is predicted to induce nonsense-mediated mRNA decay. Functional analyses of c.-115G>A, and c.*24A>G indicated an influence of these variants on gene expression. This is the second study that gives evidence for an association between pathogenic GT198 germline variants and early-onset breast cancer in HBOC.

INTRODUCTION

Roughly 5-10% of all breast and ovarian cancers occur in the context of genetic predisposition [1]. Pathogenic mutations in BRCA1 and BRCA2 account for approximately 25% of all cases of hereditary breast and ovarian cancer (HBOC) [1]. HBOC is characterized by an autosomal dominant inheritance pattern with incomplete, age-dependent penetrance, variable expressivity, an early age of breast cancer onset, and/or a positive family history with first and second degree relatives affected with breast and/or ovarian cancer [2].

Since the identification of BRCA1 [3] and BRCA2 [4] in 1994 as the two main HBOC-susceptibility genes, further genes involved in DNA-repair mechanism and/or predisposing for other rare cancer predisposing syndromes have been identified as additional low or moderate risk factors for HBOC (e.g. ATM, BARD1, BRIP1, CDH1, CHEK2, NBN1, PALB2. RAD51C, RAD51D, RECQL and TP53) [1, 5-15]. However, genetic testing of these HBOC risk genes, including BRCA1 and BRCA2, detects causative variants in no more than 40% of the families with HBOC [1, 16-17]. Thus, in the majority of cases, the HBOC predisposing genetic factors remain unknown.

GT198 has been described as a novel potential candidate gene for early-onset breast and ovarian cancer by Peng et al. [18]. GT198, also known as PSMC3IP, TBPIP (Tat binding protein interacting protein), and HOP2 (ortholog of S. cerevisiae Hop2), has been mapped 470 kb proximal of BRCA1 on chromosome 17q21 [19, 20]. It encodes for the PSMC3 (proteasome 26S subunit, ATPase, 3)-interacting protein, which is strongly expressed in adult testis and, at much lower levels, in other tissues, such as ovary and mammary gland. It acts as a transcriptional coactivator by interacting with the DNA-binding domains of nuclear receptors, such as estrogen receptor alpha and beta, thyroid hormone receptor beta 1, androgen receptor, glucocorticoid receptor, and progesterone receptor [21]. Furthermore, GT198 has been shown to stimulate RAD51 or meiotic DMC1-mediated DNA strand exchange during repair of DNA double-strand breaks [22-24]. GT198 also has an anti-apoptotic role by repressing caspase 8 activity in estrogen receptor-positive and triple-negative breast cancer cells [25].

In 2011, GT198 has been described as a novel candidate gene for primary ovarian insufficiency, when a homozygous 3 bp in-frame deletion in exon 8 (NM_016556.3, c.600_602del; p.Glu201del) was found in five affected females of a consanguineous Palestinian family with XX-female gonadal dysgenesis [26]. However, no association of mutated GT198 with primary ovarian insufficiency has been found in a cohort of 50 patients with Swedish ethnicity [27]. Subsequently, potential pathogenic germline variants in GT198 were identified at a low frequency in patients with HBOC, mostly with early cancer onset and in one patient with apparently sporadic early-onset breast cancer [18]. Deleterious somatic variants, which often cluster in the 5´-UTR and at the exon 4/intron 4 border of GT198, are abundantly detectable in breast and ovarian cancers and in fallopian tube tumors [18, 24, 28, 29]. In order to evaluate the role of GT198 in HBOC, we screened 166 BRCA1/2 mutation-negative patients, who fulfilled the diagnostic criteria of the German Consortium of Familial Breast and Ovarian Cancer (criteria details see Supplementary Table 1). Fifty-six of them developed breast cancer before the age of 36 years and, thus, were regarded as early-onset breast cancer patients (≤ 35 years). GT198 variants were investigated regarding their functional impairment.

RESULTS

A germ line nonsense mutation in GT198 has been identified in a family with hereditary breast and ovarian cancer and early-onset breast cancer and in another unrelated case with early-onset breast cancer [18]. This report prompted us to screen 166 HBOC-affected index patients, 56 of them showing early-onset breast cancer, for GT198 variants (Supplementary Table 1). We found rare or novel GT198 variants with possible pathogenic significance in 8 unrelated index cases with a family history of breast and/or ovarian cancer. (Table 1, Figure 1, Figure 2). Seven patients with GT198 variants were affected with breast cancer with a median age of cancer onset of 36 years, and one heterozygous index case was diagnosed with ovarian cancer at the age of 35 years. GT198 variants were identified in 2 out of 56 early-onset breast cancer cases (3.6%) and in 6 out of 110 breast and ovarian cancer patients (5.5%) with suspected HBOC diagnosis without early-onset (Table 1).

We identified 1 common (rs2292752, c.338-15C>G) and 5 rare nucleotide substitutions (rs191843707 (c.-115G>A), rs752276800 (c.-70T>A), rs199620968 (c.-37A>T), rs200359709 (c.-24C>G) and rs375509656 (c.537+51G>C)). These 5 variants were listed in the European population in public databases with allele frequencies of <1% (Exome Aggregation Consortium and the NCBI data base, including the 1000 Genome Project). We observed a Hardy-Weinberg equilibrium for all detected variants in cases and controls, with the exception of the common variant rs2292752 (c.338-15C>G), which was in disequilibrium in controls. We found a significant difference for the allele frequencies of rs752276800 (c.-70T>A) and rs375509656 (c.537+51G>C) between cases and controls (Table 2).

We also detected a heterozygous nonsense mutation (c.519G>A; p.(Trp173*)) in exon 6 (Table 1), which was classified as disease causing by MutationTaster, presumably by inducing nonsense-mediated mRNA decay (NMD), and its localization within the DNA-binding domain of GT198 [18]. Since there are many GT198 isoforms the prediction was made also for the protein coding transcript variants ENST00000253789 (c.483G>A; p.(Trp161*), ENST00000587209 (c.330G>A; p.(Trp110*) and ENST00000590760 (c.144G>A; p.(Trp48*). The affected amino acid tryptophan is highly conserved among vertebrates (PhyloP:6,302, PhasCons:1). This truncating mutation is also listed in the COSMIC database (mutation ID 4431647) and has been detected as a somatic variant by exome sequencing in one patient with esophagus squamous cell carcinoma [30]. In our cohort, the nonsense mutation p.(Trp173*) was found in two sisters (F16 and F17), which were both diagnosed with unilateral breast cancer (invasive ductal carcinomas) at 33 years of age (Table 1, Figure 1, Figure 2). One of the sisters was also heterozygous for the c.-37A>T variant (rs199620968). The c.-37A>T substitution was also detected in two other unrelated cancer patients, each in the heterozygous state: in patient C4, which was diagnosed with ovarian cancer at the age of 35 years followed by unilateral invasive ductal breast cancer at the age of 61 years, and in index case D13, which was affected by unilateral breast cancer at the age of 68 (Table 1). The heterozygous substitution c.-37A>T was also present in the 43-year-old meningioma affected niece (D17) of index patient D13.

As copy number gains of a mutated GT198 allele with a nonsense mutation have recently been reported in a breast cancer affected patient [18], we additionally screened the index case F16 and her sister F17 for copy number changes by a custom-made 60k eArray. No copy number gains or losses of GT198 and no further structural rearrangements were detected by array-CGH. From all identified substitutions in the 5’-UTR (i.e. rs191843707 (c.-115G>A), rs752276800 (c.-70T>A), rs199620968 (c.-37A>T) and rs200359709 (c.-24C>G)), the variant c.-115G>A was found once in our own study in a female that developed unilateral breast cancer at the age of 33 years (Table 1, Figure 1, Figure 2) and has previously been described in hereditary breast and ovarian cancer [18]. It was predicted that c.-115G>A deleteriously alters the binding site for the ETS domain-containing factor ELK1 at positions c.-111_-120 (Supplementary Figure 1). To evaluate whether the predicted effects on transcription factor-binding might influence GT198 expression, luciferase assays were performed for all identified 5´-UTR variants in HEK293T cells (Figure 3). Transfection of the c.-115A construct showed a significant decrease of relative luciferase activity by 22% compared to the wild type allele (Figure 3). In contrast, for c.-70T>A, c.-37A>T, and c.-24C>G, no significant differences of luciferase activities were detected. In silico analysis of the intron 6 variant c.537+51G>C (Table 1) provided no evidence for altered splicing.

We further detected a novel nucleotide substitution within the 3’-UTR (c.*24G>A) of GT198 in patient H16 diagnosed with unilateral breast cancer at the age of 36 years (Table 1, Figure 1, Figure 2). The variant is also present in her 32-year-old unaffected sister and in her cousin H17 diseased from breast cancer at the age of 25 years (Table 1, Figure 2). In silico analyses indicated that only less conserved microRNA-binding sites (hsa-miR-1224-3p, hsa-miR-1280, hsa-miR-2114, hsa-miR-2355-5p, and hsa-miR-4286) were affected by c.*24G>A. To investigate whether this position is important for microRNA-binding, luciferase assays according to Buurman et al. [31] were performed (Figure 3). A significant decrease of luciferase activity of approximately 42% was observed when the wildtype 3’-UTR was introduced downstream to the reporter gene (Figure 3). Luciferase activity in HEK293T cells transfected with the c.*24A construct was decreased by only 26% compared to the control, suggesting a negative effect of c.*24G>A on microRNA-binding and, thus, to an increased gene expression.

We further confirmed the variants c.-70T>A, c.-37A>T, c.-24C>G, c.537+51 G>C and the nonsense mutation p.(Trp173*) in FFPE-breast cancer samples of the respective index cases. For the variants c.-115G>A, and c.*24G>A, no tumor material was available (Table 1). As somatic variants are frequently observed in breast and ovarian cancers [28, 29], we sequenced GT198 in the available tumor samples (Table 1, Supplementary Table 2). No additional second hits were detected.

We further screened all 8 index cases carrying GT198 variants for pathogenic variant in additional low or moderate risk genes for HBOC (e.g ATM, CDH1, CHEK2, NBN1, PALB2, RAD51C, RAD51D, and TP53) using the TruSight Cancer panel (Illumina, San Diego CA). No pathogenic variants were detected in any of these additional risk genes. In order to screen for HBOC predisposing copy number changes in the 8 carriers of GT198 variants, a customized high-resolution 8x60k array (Design:069100, HBOC-2, Agilent technologies) for comparative genomic hybridization CGH covering the 94 genes of the TruSight cancer panel was used [32]. We identified no aberrant copy number changes in the HBOC cancer risk genes [5, 8-13, 15], including BRCA1 and BRCA2, by high resolution array CGH.

DISCUSSION

GT198’s location in a genomic region on 17q21, previously linked to hereditary breast and ovarian cancer, makes an association of GT198 disease-causing changes with HBOC and sporadic early-onset breast cancer likely [18-20, 24, 33]. Germline variants with possible pathogenic impact have been found in HBOC cases with mostly early onsets (median age 35 years) and in an apparently sporadic case of breast cancer with an onset age of 30 years [18].

Eight out of 166 unrelated index cases (4.8%) in our study were heterozygous for rare or novel GT198 variants with yet unknown impact on GT198 function, which is similar to the detection frequency of the first report, in which 8 out of 212 index patients (3.8%) have been heterozygous for putative pathogenic GT198 germline variants [18]. Three out of 8 heterozygous carriers of GT198 variants in the present study showed early-onset of breast or ovarian cancer (≤35 years), which is, albeit in lower frequency, congruent with the former findings of Peng et al. [18], who reported that 6 out of 8 index cases carrying GT198 variants were affected by early-onset breast or ovarian cancer.

We identified two germ line variants in GT198 (c.519G>A p.(Trp173*) and c.*24G>A), that were neither listed in the NCBI database nor in SNP data bases of the ExAc Browser and of EVS. The nonsense mutation c.519G>A; p.(Trp173*) is reported once in the COSMIC database and was detected in a human esophagus squamous cell carcinoma [30]. Our own in silico analysis suggest a negative impact of c.519G>A on GT198 expression and function. The induced premature translation termination codon of c.519G>A is located in the DNA binding domain of GT198, and is predicted to induce nonsense-mediated mRNA decay of the aberrant transcript [34]. This stop codon is located downstream of an alternative translation initiation codon within exon 5, that leads to the expression of a truncated protein isoform, harboring the DNA binding domain and the C-terminus [24]. The DNA binding domain is able to bind both single- and double-stranded DNA and is important for GT198 DNA repair activity after DNA double-strand breaks [18]. It has been shown by in vitro assays that especially the amino acid residues 171-178 of murine Hop2 at the C-terminus, which are 100% identical to the human orthologue, have a high affinity for single-stranded DNA [22], but it still remains unknown, whether amino acid residue 173 is indispensable for RAD51 single-stranded DNA presynaptic filament stabilization or homologous DNA pairing, which are both important for RAD51-mediated homologous recombination of damaged chromosomes [23]. Its presence in two breast cancer affected sisters that were both diagnosed having cancer at 33 years of age (Table 1, Figure 2), makes a positive association of c.519G>A with early-onset breast cancer likely. A pathogenic nonsense mutation (c.310C>T; p.(Q104*)) has already been identified in a former study by Peng et al. [18] in two unrelated female breast cancer patients that were diagnosed with breast cancer at the age of 30 and 33 years. The induced premature stop codon affects the leucine zipper dimerization domain of GT198, which is also required for protein-protein interaction and transcriptional regulation and has been shown in in vitro cell culture experiments to abolish RAD51-mediated DNA repair activity after γ-irradiation [18]. It is assumed that mutated and alternate transcripts are counteracting in a dominant negative manner with wildtype GT198 [18, 24].

Some microRNA-binding sites are predicted in the 3’-UTR of GT198. The substitution c.24*G>A is located within a weakly-conserved microRNA-binding site for hsa-miR-1224-3p, hsa-miR-1280, hsa-miR-2114, hsa-miR-2355-5p, and hsa-miR-4286. However, it is currently unknown whether GT198 expression is regulated by one of these in vivo. Our own in vitro data points to an impaired effect of the variant c.24*G>A on microRNA-binding. Whether this would also affect GT198 expression in vivo requires further elucidation. The c.24*G>A substitution was found to segregate with early-onset breast cancer in an affected female maternal cousin, but was also present in a 32-year-old healthy sister of the index case (Table 1, Figure 2).

Interestingly, rare variants within the GT198 5’-UTR and 3’-UTR were frequently found in breast, ovarian and fallopian tube cancers [18, 24, 28, 29]. Six out of 166 cancer patients of our cohort were heterozygous for rare or novel nucleotide substitutions located in the 5´-UTR or 3’-UTR of GT198. Two of these variants, c.-115G>A and c.-37A>T, have already been identified as hereditary variants in familial cases of breast and ovarian cancer [18]. Interestingly, one variant in the 5’-UTR (c.-37A>T) has been found in the germline of two unrelated familial cases in our own and in a former study [18]. One of the carriers of c.-37A>T was affected by ovarian cancer at early-onset (35 years). This variant has previously been presented as a somatic variant in serous ovarian carcinoma, fallopian tube cancer, and endometrial carcinoma [24]. All identified GT198 5’-UTR variants are also detectable in the general population, albeit at low allele frequency, and a possible disease association still remains unknown. Especially the substitution c.-115G>A is listed in the European population of the Exome Aggregation Consortium with an allele frequency of 0.64%. Our own in silico and in vitro data suggest no influence of c.-70T>A, c.-37A>T and c.-24C>G on GT198 expression. In contrast, the variant c.-115G>A induced a slightly, albeit significant, decrease of reporter gene expression in transiently transfected HEK293T cells. Our own in silico predictions led us to speculate if this effect could be mediated by destroying a binding site for the transcription factor ELK1 (ETS domain-containing protein Elk-1), a member of the ETS family of transcription factors, which regulates the expression of genes involved in cell proliferation, chromatin modelling and apoptosis [35, 36]. ELK1 overexpression is frequently observed in many carcinomas, including breast cancer [36]. Recently, the MZF1/Elk-1 complex has been identified as mediator of protein kinase C alpha (PKCα) expression in triple-negative breast cancer, which induces cell migration and invasion of triple negative breast cancer cells and poor outcome [37]. However, the pathogenic effect of altered GT198 expression in general and the influence of c.-115G>A on GT198 expression in breast or ovary in vivo is still unknown and requires further elucidation.

The minor alleles c.-70T>A and c.537+51G>C were observed significantly more frequent in cases than in European controls of the ExAC database (Table 2). Our own analyses suggest a benign effect of both variants on GT198 function, and we, therefore, ascribe the discrepancies in allele frequencies to the small sample size of the analyzed case cohort.

We here present a screen for pathogenic changes in GT198 in patients with BRCA1/2-negative HBOC. We identified seven different rare or novel GT198 variants with yet unknown impact on GT198 function, of which six were absent or extremely rare in the ExAC database in Europeans. Three variants (c.-115G>A, c.519G>A and c*24G>A) found in familial breast cancer patients with early-onset at ≤ 36 years seem to have an impact on GT198 function and may contribute to breast cancer predisposition. GT198 participation in steroid hormone receptor-mediated gene expression, its function in DNA recombination, and its ability to stimulate RAD51 mediated DNA strand exchange [18, 23, 29], makes its implication in oncogenesis conceivable. Further, comprehensive mutation screenings in multi-national case and control collectives are required to evaluate the role of GT198 in breast and ovarian cancer predisposition.

MATERIAL AND METHODS

Study cohort

The analyzed cohort was composed of 166 unrelated, female breast and/or ovarian cancer patients of mixed Caucasian, mostly German origin, which were referred to our outpatient clinic between 2004 and 2014. All selected individuals, including familial cases (n=158) and early-onset breast cancer cases without a positive family history (n=8), fulfilled the inclusion criteria for BRCA1 and BRCA2 testing of the German Consortium for HBOC (Supplementary Table 1) [38]. All women have given their informed consent for participating in the study, which was approved by the hospital´s ethics committee (Hannover Medical School, ethic votum 4121). Comprehensive data about the family history, including data about breast and ovarian cancer development over at least three generations, tumor pathology, and about BRCA1/2 mutational status were available for each case. All samples were previously shown to be negative for deleterious variants within BRCA1 and BRCA2 using routine diagnostic methods, including sequencing and multiplex ligation probe-dependent amplification (MRC Holland, Amsterdam, the Netherlands) for BRCA1. For 162 out of 166 samples structural aberrations affecting BRCA2 were also excluded by MLPA and/or array-CGH analyses (Agilent Technologies, Santa Clara, United States). However, 5 samples were heterozygous for missense variants of unknown significance in BRCA2 (ENST00000380152, c.831T>G, p.(Asn277Lys); c.995T>A, p.(Ile332Asn); c.4782G>T, p.(Met1594Ile); c.6101G>A, p.(Arg2034His); c.7562T>C p.(Ile2521Thr)).

The cohort encompasses 155 individuals with breast cancer (135 unilaterally and 20 bilaterally affected women), 9 individuals with ovarian cancer and 2 individuals with breast and ovarian cancer. Age at diagnosis of breast cancer ranged from 17 to 68 years (median 39 years), while the age of onset for ovarian cancer vary from 21 to 67 years (median 33.5 years). In 56 out of 155 breast cancer patients, diagnosis was made before 36 years. The majority of index patients (n=90) originated from families with at least two first and/or second degree female relatives affected by breast cancer, of whom one individual was diagnosed before age of 51 years. For 32 patients, there was a family history of at least one female breast and one ovarian cancer or one woman with breast and ovarian cancer. Thirty one samples were derived from families with at least one woman diagnosed with bilateral breast cancer before age of 51 years. Eight patients developed early-onset (≤ 35 years) breast cancer and had no familial history for HBOC. Three patients had a family history with at least 3 first or second degree relatives with breast cancer and two recruited index patients had a family history with at least two ovarian cancers (Supplementary Table 1).

DNA extraction and sequencing

Genomic DNA from EDTA blood samples and buccal mucosa smears were extracted using QIAamp DNA Blood Midi and Mini Kit (Qiagen, Hilden, Germany), respectively. From selected breast cancer cases, genomic DNA from formalin-fixed and paraffin-embedded (FFPE) tumor tissues was extracted from 5 µm serial sections using GeneRead DNA FFPE Kit (Qiagen) according to the manufacturer’s instructions. All tissue samples were previously histologically examined by pathologists. A hematoxylin and eosin-stained section of each tumor paraffin block was histologically examined to define the area with 15-80 % tumour cells to be macro-dissected for DNA extraction.

All 8 GT198 exons encompassing the entire coding exons, adjacent intronic regions and parts of the flanking 5´-UTR and 3´-UTR were PCR amplified and subsequently sequenced using an ABI genetic analyzer 3130xl (Applied Biosystems, Darmstadt, Germany). For FFPE samples a set of additional primers was used. Primers were designed using the software Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/primer3/) (Supplementary Table 2). Illustra ExoProStar 1-Step (GE Healthcare, Munich, Germany) was used for PCR product purification before sequencing PCRs. Finally, sequencing products were cleaned up by Sephadex G-50 purification (Sigma-Aldrich Chemie GmbH, Steinheim, Germany). Variant analyses were performed using Sequence Pilot version 4.3.0 (JSI Medical Systems, Ettenheim, Germany) and the NCBI sequence NM_016556.3 as reference.

The TruSight Cancer panel (Illumina, San Diego CA) was used for target enrichment of ATM, CHEK2, CDH1, NBN1, PALB2, RAD51C, RAD51D and TP53 of selected index cases. Sequence analyses for the additional HBOC risk genes were performed with the module NextSeq of Sequence Pilot version 4.3.0 (JSI Medical Systems, Ettenheim, Germany) and the NCBI transcripts NM_000051 (ATM), NM_007194.3 (CHEK2), NM_002485 (NBN1), and the Ensembl transcripts ENST00000261769 (CDH1), ENST00000261584 (PALB2), ENST00000337432 (RAD51C), ENST00000345365 (RAD51D) and ENST00000269305 (TP53) as references.

Variant and statistical analyses

All identified rare variants were validated by PCR amplification and sequencing of an independent DNA sample.

Detailed in silico predictions for identified GT198 variants were made with tools implemented in alamut visual (Align GVGD, SIFT, MutationTaster, Polyphen2) (interactive biosoftware, Version 2.7 rev.1, Rouen, France) [39-42]. Data from the European population retrieved from the ExAC Browser (http://exac.broadinstitute.org/) or the NCBI database (https://www.ncbi.nlm.nih.gov/nuccore) and the European-African population from the exome variant server (http://evs.gs.washington.edu/EVS/) served as controls. Variants within the 5´-UTR were analyzed for putative effects on transcription factor binding using “JASPAR” (version 5.0_ALPHA, http://jaspar.binf.ku.dk/) and PROMO (http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3). For splice analyses the programs Berkeley Drosophila Genome Project (http://www.fruitfly.org/seq_tools/splice.html), and ASSP (http://wangcomputing.com/assp/index.html) were used. Additionally, a splice analysis using alamut visual (interactive biosoftware) with the included tools SpliceSiteFinder-like, MaxEntScan, NNSPLICE, GeneSplicer and Human Splicing Finder (HSF2.4.1) was also performed. MicroRNA-binding site predictions were made with the program Alamut visual based on miRanda predictions and microRNA.org targets. The effect of the nonsense mutation was evaluated by the program MutationTaster (www.mutationtaster.org) [42].

To evaluate whether the genotype distribution of GT198 variants of cases and controls were in Hardy-Weinberg equilibrium (HWE), a chi-squared goodness-of-fit test with one degree of freedom was performed. For HWE analysis the online software tool http://www.koonec.com/k-blog/2010/06/20/ and hardy-weinberg-equilibrium-calculator/from Strom and Wienker were used (http://ihg2.helmholtz-muenchen.de/cgi-bin/hw/hwa1.pl). Allele frequencies between cases and controls were compared with Fisher´s exact test for low number of individuals carrying the rare allele (≤5) or with Pearson’s chi squared goodness-of-fit test (R-software,https://www.r-project.org/). P-values of these comparisons were assessed descriptively, and defined to be statistically significant if p<0.05. As control, datasets (genotype and allele frequencies) for the European non-Finnish population retrieved from the ExAC Browser (http://exac.broadinstitute.org/) were used.

Cloning and Luciferase assays

To investigate the effects of c.-115G>A, c.-70T>A, c.-37A>T, c.-24C>G and c.*24G>A on gene expression, luciferase reporter assays were performed in HEK293T cells. To subclone genomic regions immediately 5’ of the start and 3’ of the stop codon of GT198, specific primers tagged with sequences containing restriction enzyme sites were used for amplification (Supplementary Table 2). Primers GT198_5’UTR_XhoI_f and GT198_5’UTR_HindIII_r were used to amplify a 564 bp genomic fragment on the corresponding patients’ DNA (A13, B19, D13, E10, Table1, Figure 2). Following an XhoI/HindIII double digestion, the fragment was integrated upstream of the firefly luciferase gene of the pGL3.basic vector (Promega, Mannheim, Germany). Primers GT198_3’UTR_XbaI_f and GT198_3’UTR_SpeI_r were used for amplification of a 1,070 bp amplicon on the corresponding patients’ DNA (H16, Table1, Figure 2). Amplicon and pGL3.promoter vector (Promega, Mannheim, Germany) were cut with XbaI and the amplicon was subcloned downstream of the firefly luciferase gene. Inserts of the obtained vectors containing mutant or wildtype alleles were evaluated by sequencing using an ABI genetic analyzer 3130xl (Applied Biosystems). Thus, seven different vectors containing either the wildtype (p.GL3.basic.5’UTR.WT and p.GL3.promoter.3’UTR.WT) or the rare GT198 alleles (p.GL3.basic.5’UTR.c.-115G>A, p.GL3.basic.5’UTR.c.-70T>A, p.GL3.basic.5’UTR.c.-37A>T, p.GL3.basic.5’UTR.c.-24C>G, and p.GL3.promoter.3’UTR.c.*24G>A) were generated.

HEK293T cells were cultured in Dulbecco´s modified Eagle´s medium supplemented with 1 mM sodium pyruvate, 10% heat-inactivated fetal bovine serum, 100 units/ml penicillin and 100 µg/ml streptomycin in a humidified atmosphere with 5% CO₂ at 37 ° C. For each of the seven firefly luciferase reporter constructs triplicates of 8,000 HEK293T cells were seeded in a 96-well plate in 100 µl medium. After 24 hours, cells were cotransfected with 25 ng luciferase reporter plasmid and 2.5 ng pGL-4.70 (Promega) using Lipofectamine 2000 (Invitrogen, Paisley, UK). Twenty-four hours after transfection, cells were lysed and firefly and renilla luciferase activity were measured using the Dual-Glo^® Luciferase Assay System (Promega) and a Synergy 2 Multi-Mode Microplate Reader (BioTek, Winooski, VT) in accordance to the manufacturer´s instructions.

eArray-CGH analysis

EDTA blood-derived genomic DNA (average DLRS value 0,1) of F16 and F17, carrying the GT198 truncating mutation p.(Trp173*) (Table 1), was further screened for copy number changes and structural rearrangements affecting GT198 using a custom-made 60k eArray (Design:081270, brca1-2region, Agilent Technologies) with high resolution in the genomic regions on chromosome 17 (chr17:40,600,000-41,756,000, 1.2 Mb, GRCh37/hg19) and chromosome 13 (chr13:32,414,344-33,490,000, 1.1 Mb, GRCh37/hg19), with an average probe spacing of 300 bp. The GT198 gene was covered by 85 probes.

Selected index cases (n=8) carrying putative pathogenic GT198 variants were also screened for HBOC predisposing copy number variations by a customized high resolution 60k eArray (Design:069100, HBOC-2, Agilent technologies) [32]. Array CGH analysis was performed as recommended by the manufacturer. The female human DNA EA-100F was used as control (Kreatech Biotechnology, Amsterdam, The Netherlands). Fluorescence signals were scanned using a Dual Laser Scanner G2565CA (Agilent Technologies). Raw data analysis was performed using Feature extraction version 11.0.1.1 (Agilent Technologies). For further data analysis, Genomic Workbench 7.0.4.0 (Agilent Technologies) was used: ADM-2 algorithm, threshold 6, and no aberration filter for the brca1-2region, while a 2log0.2 filter was used for the HBOC-2 design.

AUTHORS CONTRIBUTIONS

Stephanie Schubert: Conception and design of the study, analysis and interpretation of the data, drafting of the manuscript

Tim Ripperger: Cloning and Luciferase assays in HEK293T cells, scientific adviser and critically revised the manuscript.

Melanie Rood: acquisition of data (PCR and Sanger sequencing and in silico analyses)

Anthoy Petkidis: acquisition of data (PCR and Sanger sequencing and in silico analyses)

Winfried Hofmann: statistical analyses and scientific advise

Hildegard Frye-Boukhriss: acquisition of data (assisted in sample preparation, PCR and Sanger-sequencing)

Marcel Tauscher: acquisition of data (eArray-CGH analyses)

Bernd Auber: acquisition of data (TruSight Cancer panel)

Ursula Hille-Betz: acquisition of data (examined participants, assisted interpretation of data, critically revised the manuscript

Thomas Illig: study design and critical revision of the manuscript

Brigitte Schlegelberger: supervised the work, critical reviewed and edited the manuscript

Doris Steinemann: study concept and design of the study, scientific adviser, interpretation of results, critical reviewed and edited the manuscript

ACKNOWLEDGMENTS

The authors thank Claudia Davenport for critical reading and revising a first version of the manuscript.

CONFICTS OF INTEREST

The authors declared that none of the contributors has commercial, professional or financial affiliations that may be perceived as having biased the presentation. Furthermore all authors have agreed that no conflict of interest with respect to the research, authorship and publication exists. The data presented here have not been submitted, published, and are not being submitted elsewhere for publication.

FUNDING

This work is supported in part by the Claudia-von-Schilling foundation (to D.S. and B.S.) and the German research Foundation (DFG, Cluster of Excellence, REBIRTH).

• 1. Nielsen FC, van Overeem Hansen T, Sørensen CS. Hereditary breast and ovarian cancer: new genes in confined pathways. Nat Rev Cancer. 2016; 16:599-612. [PubMed]
• 2. Kobayashi H, Ohno S, Sasaki Y, Matsuura M. Hereditary breast and ovarian cancer susceptibility genes. Oncol Rep. 2013; 30:1019-1029. [PubMed]
• 3. Miki Y, Swensen J, Shattuck-Eidens D, Futreal PA, Harshman K, Tavtigian S, Liu Q, Cochran C, Bennett LM, Ding W, Bell R, Rosenthal J, Hussei C et al. A strong candidate for the breast and ovarian cancer susceptibility gene BRCA1. Science. 1994; 266:66-71. [PubMed]
• 4. Wooster R, Neuhausen SL, Mangion J, Quirk Y, Ford D, Collins N, Nguyen K, Seal S, Tran T, Averill D, Fields P, Marshall G, Narod S et al. Localization of a breast cancer susceptibility gene, BRCA2, to chromosome 13q12-12. Science. 1994; 265:2088-2090. [PubMed]
• 5. Athma P, Rappaport R, Swift M. Molecular genotyping shows that ataxia-telangiectasia heterozygotes are predisposed to breast cancer. Cancer Genet Cytogenet. 1996; 92:130-134. [PubMed]
• 6. De Brakeleer S, De Grève J, Loris R, Janin N, Lissens W, Sermijn E, Teugels E. Cancer predisposing missense and protein truncating BARD1 mutations in non-BRCA1 or BRCA2 breast cancer families. Hum Mutat. 2010; 31:E1175-1185. [PubMed]
• 7. Seal S, Thompson D, Renwick A, Elliott A, Kelly P, Barfoot R, Chagtai T, Jayatilake H, Ahmed M, Spanova K, North B, McGuffog L, Evans DG et al. Truncating mutations in the Fanconi anemia J gene BRIP1 are low-penetrance breast cancer susceptibility alleles. Nat Genet. 2006; 38:12391241. [PubMed]
• 8. Guilford PJ, Hopkins JB, Grady WM, Markowitz SD, Willis J, Lynch H, Rajput A, Wiesner GL, Lindor NM, Burgart LJ, Toro TT, Lee D, Limacher JM, et al. E-cadherin germline mutations define an inherited cancer syndrome dominated by diffuse gastric cancer. Hum Mutat. 1999; 14:249-255. [PubMed]
• 9. Meijers-Heijboer H, van den Ouweland A, Klijn J, Wasielewski M, de Snoo A, Oldenburg R, Hollestelle A, Houben M, Crepin E, van Veghel-Plandsoen M, Elstrodt F, van Duijn C, Bartels C et al. Low-penetrance susceptibility to breast cancer due to CHEK2(*)1100delC in noncarriers of BRCA1 or BRCA2 mutations. Nat Genet. 2002; 31:5559. [PubMed]
• 10. Górski B, Debniak T, Masojć B, Mierzejewski M, Medrek K, Cybulski C, Jakubowska A, Kurzawski G, Chosia M, Scott R, Lubiński J. Germline 657del5 mutation in the NBS1 gene in breast cancer patients. Int J Cancer. 2003; 106:379-381. [PubMed]
• 11. Rahman N, Seal S, Thompson D, Kelly P, Renwick A, Elliott A, Reid S, Spanova K, Barfoot R, Chagtai T, Jayatilake H, McGuffog L, Hanks S et al. PALB2, which encodes a BRCA2-interacting protein, is a breast cancer susceptibility gene. Nat Genet. 2007; 39:165-167. [PubMed] https://doi.org/10.1038/ng1959.
• 12. Meindl A, Hellebrand H, Wiek C, Erven V, Wappenschmidt B, Niederacher D, Freund M, Lichtner P, Hartmann L, Schaal H. Ramser J, Honisch E, Kubisch C et al. Germline mutations in breast and ovarian cancer pedigrees establish RAD51C as a human cancer susceptibility gene. Nat Genet. 2010; 42:410-414. [PubMed]
• 13. Loveday C, Turnbull C, Ramsay E, Hughes D, Ruark E, Frankum JR, Bowden G, Kalmyrzaev B, Warren-Perry M, Snape K, Adlard JW, Barwell J, Berg J et al. Germline mutations in RAD51D confer susceptibility to ovarian cancer. Nat Genet. 2011; 43:879-882. [PubMed] https://doi.org/10.1038/ng.893.
• 14. Cybulski C, Carrot-Zhang J, Kluźniak W, Rivera B, Kashyap A, Wokołorczyk D, Giroux S, Nadaf J, Hamel N, Zhang S, Huzarski T, Gronwald J, Byrski T et al. Germline RECQL mutations are associated with breast cancer susceptibility. Nat Genet. 2015; 47:643-646. [PubMed]
• 15. Malkin D, Li FP, Strong LC, Fraumeni JF Jr, Nelson CE, Kim DH, Kassel J, Gryka MA, Bischoff FZ, Tainsky MA, Friend SH. Germ line p53 mutations in a familial syndrome of breast cancer, sarcomas, and other neoplasms. Science. 1990; 250:1233-1238. [PubMed]
• 16. Meindl A, Ditsch N, Kast, K, Rhiem K, Schmutzler RK. Hereditary breast and ovarian cancer: new genes, new treatments, new concepts. Dtsch Arztebl Int. 2011; 108:323330. [PubMed] https://doi.org/10.3238/arztebl.2011.0323.
• 17. Kurian AW, Hare EE, Mills MA, Kingham KE, McPherson L, Whittemore AS, McGuire V, Ladabaum U, Kobayashi Y, Lincoln SE, Cargill M, Ford JM. Clinical evaluation of a multiple-gene sequencing panel for hereditary cancer risk assessment. J Clin Oncol. 2014; 32:2001-2009. [PubMed] https://doi.org/10.1200/JCO.2013.53.6607.
• 18. Peng M, Bakker JL, Dicioccio RA, Gille JJ, Zhao H, Odunsi K, Sucheston L, Jaafar L, Mivechi NF, Waisfisz Q, Ko L. Inactivating Mutations in GT198 in Familial and EarlyOnset Breast and Ovarian Cancers. Genes Cancer. 2013; 4:15-25. [PubMed] https://doi.org/10.1177/1947601913486344.
• 19. Ijichi H, Tanaka T, Nakamura T, Yagi H, Hakuba A, Sato M. Molecular cloning and characterization of a human homologue of TBPIP, a BRCA1 locus-related gene. Gene. 2000; 248:99-107. [PubMed]
• 20. Rommens JM, Durocher F, McArthur J, Tonin P, LeBlanc JF, Allen T, Samson C, Ferri L, Narod S, Morgan K, Simard J. Generation of a transcription map at the HSD17B locus centromeric to BRCA1 at 17q21. Genomics. 1995; 28:530542. [PubMed]
• 21. Ko L, Cardona GR, Henrion-Caude A, Chin WW. Identification and characterization of a tissue-specific coactivator, GT198, that interacts with the DNA-binding domains of nuclear receptors.Mol Cell Biol. 2002; 22:357369.. 1995; 28:530542. [PubMed] https://doi.org/10.1128/MCB.22.1.357-369.2002.
• 22. Zhao W, Saro D , Hammel M, Kwon Y, Xu Y, Rambo RP, Williams GJ, Chi P, Lu L, Pezza RJ, Camerini-Otero RD, Tainer JA, Wang HW et al. Mechanistic insights into the role of Hop2-Mnd1 in meiotic homologous DNA pairing. Nucleic Acids Res. 2014; 42:906-917. [PubMed] https://doi.org/10.1093/nar/gkt924.
• 23. Zhao W, Sung P. Significance of ligand interactions involving Hop2-Mnd1 and the RAD51 and DMC1 recombinases in homologous DNA repair and XX ovarian dysgenesis. Nucleic Acids Res. 2015; 43:4055-4066. [PubMed] https://doi.org/10.1093/nar/gkv259.
• 24. Peng M, Yang Z, Zhang H, Jaafar L, Wang G, Liu M, Flores-Rozas H, Xu J, Mivechi NF, Ko L. GT198 Splice Variants Display Dominant-Negative Activities and Are Induced by Inactivating Mutations. Genes Cancer. 2013; 4:26-38. [PubMed] https://doi.org/10.1177/1947601913486345.
• 25. Capdevila-Busquets E, Badiola N, Arroyo R, Alcalde V, Soler-López M, Aloy P. Breast cancer genes PSMC3IP and EPSTI1 play a role in apoptosis regulation. PLoS One. 2015; 10:e0115352. [PubMed] https://doi.org/10.1371/journal.pone.0115352.
• 26. Zangen D, Kaufman Y, Zeligson S, Perlberg S, Fridman H, Kanaan M, Abdulhadi-Atwan M, Abu Libdeh A, Gussow A, Kisslov I, Carmel L, Renbaum P, Levy-Lahad E. XX ovarian dysgenesis is caused by a PSMC3IP/HOP2 mutation that abolishes coactivation of estrogen-driven transcription. Am J Hum Genet. 2011; 89:572-579. [PubMed] https://doi.org/10.1016/j.ajhg.2011.09.006.
• 27. Norling A, Hirschberg AL, Karlsson L, Rodriguez-Wallberg KA, Iwarsson E, Wedell A, Barbaro M. No mutations in the PSMC3IP gene identified in a Swedish cohort of women with primary ovarian insufficiency. Sex Dev. 2014; 8:146150. [PubMed]
• 28. Peng M, Zhang H, Jaafar L, Risinger JI, Huang S, Mivechi NF, Ko L. Human ovarian cancer stroma contains luteinized theca cells harboring tumor suppressor gene GT198 mutations. J Biol Chem. 2013; 288:33387-33397. [PubMed] https://doi.org/10.1074/jbc.M113.485581.
• 29. Yang Z, Peng M, Cheng L, Jones K, Maihle NJ, Mivechi NF, Ko L. GT198 Expression Defines Mutant Tumor Stroma in Human Breast Cancer. Am J Pathol. 2016; 186:1340-1350. [PubMed] https://doi.org/10.1016/j.ajpath.2016.01.006.
• 30. Gao YB, Chen ZL, Li JG, Hu XD, Shi XJ, Sun ZM, Zhang F, Zhao ZR, Li ZT, Liu ZY, Zhao YD, Sun J, Zhou CC et al. Genetic landscape of esophageal squamous cell carcinoma. Nat Genet. 2014; 46:1097-1102. [PubMed]
• 31. Buurman R, Gürlevik E, Schäffer V, Eilers M, Sandbothe M, Kreipe H, Wilkens L, Schlegelberger B, Kühnel F, Skawran B. Histone deacetylases activate hepatocyte growth factor signaling by repressing microRNA-449 in hepatocellular carcinoma cells. Gastroenterology. 2012; 143:811-820. [PubMed]
• 32. Hackmann K, Kuhlee F, Betcheva-Krajcir E, Kahlert AK, Mackenroth L, Klink B, Di Donato N, Tzschach A, Kast K, Wimberger P, Schrock E, Rump A. Ready to clone: CNV detection and breakpoint fine-mapping in breast and ovarian cancer susceptibility genes by high-resolution array CGH. Breast Cancer Res Treat. 2016; 159:585-590. [PubMed]
• 33. KingMC. “The race” to clone BRCA1. Science. 2014; 343:1462-1465. [PubMed]
• 34. Frischmeyer PA, Dietz HC. Nonsense-mediated mRNA decay in health and disease. Hum Mol Genet. 1999; 8:18931900. [PubMed]
• 35. Niida A, Smith AD, Imoto S, Tsutsumi S, Aburatani H, Zhang MQ, Akiyama T. Integrative bioinformatics analysis of transcriptional regulatory programs in breast cancer cells. BMC Bioinformatics. 2008; 9:404. [PubMed] https://doi.org/10.1186/1471-2105-9-404.
• 36. Laliotis A, Vrekoussis T, Kafousi M, Sanidas E, Askoxilakis J, Melissas J, Mavroudis D, Castanas E, Stathopoulos EN. Immunohistochemical study of pElk-1 expression in human breast cancer: association with breast cancer biologic profile and clinicopathologic features. Breast. 2013; 22:89-95. [PubMed]
• 37. Lee CJ, Hsu LS, Yue CH, Lin H, Chiu YW, Lin YY, Huang CY, Hung MC, Liu JY. MZF-1/Elk-1 interaction domain as therapeutic target for protein kinase Cα-based triplenegative breast cancer cells. Oncotarget. 2016; 7:5984559859. [PubMed] https://doi.org/10.18632/oncotarget.11337.
• 38. Kast K, Rhiem K, Wappenschmidt B, Hahnen E, Hauke J, Bluemcke B, Zarghooni V, Herold N, Ditsch N, Kiechle M, Braun M, Fischer C, Dikow N et al. German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC). Prevalence of BRCA1/2 germline mutations in 21401 families with breast and ovarian cancer. J Med Genet. 2016; 53:465-471. [PubMed]
• 39. Tavtigian SV, Deffenbaugh AM, Yin L, Judkins T, Scholl T, Samollow PB, de Silva D, Zharkikh A, Thomas A. Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral. J Med Genet. 2006; 43:295-305. [PubMed] https://doi.org/10.1136/jmg.2005.033878.
• 40. Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001; 11:863-874. [PubMed] https://doi.org/10.1101/gr.176601.
• 41. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods 2010; 7:248-249.. 2001; 11:863-874. [PubMed] https://doi.org/10.1038/nmeth0410-248.
• 42. Schwarz JM, Rödelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010; 7:575-576. [PubMed]