Signatures of defective DNA repair and replication in early-onset renal cancer patients referred for germline genetic testing

eoRCC patients at the Fox Chase Cancer Center (FCCC) and family history of cancer

We analyzed the personal and family history of the probands in a cohort of 22 eoRCC patients. Multiple probands (6/22, 27%) had a second primary cancer, with breast cancer diagnosed in 3 probands (3/22, 14%) prior to diagnosis of RCC (Table 1). Here, 73% (n=16/22) of probands had a family history of RCC, with 50% (n=11/22) of probands having a first-degree relative with RCC. Intriguingly, 64% (n=14/22) of probands had a family history of cancers of the prostate, bladder, and thyroid, and melanoma, which have been associated with an RCC diagnosis (9).

Analysis of whole-exome sequencing data reveals enrichment of germline variation in DNA repair and replication genes in eoRCC patients

We performed WES on lymphocyte DNA from the 22 eoRCC probands, detecting candidate PGVs based on a candidate gene list that was developed in our prior studies ((10); see (Supplementary Table 1) and Supplementary Materials and Methods). Briefly, the candidate gene list was developed by a comprehensive hypothesis-driven framework with the following assumptions: 1) genes involved in genome stability (using Gene Ontology terms such as DNA repair, DNA replication, DNA damage checkpoints, cell cycle, mitotic machinery, replication stress, DNA damage response, chromatin remodeling) would be important for hereditary cancer risk (5, 11-13); and 2) an expanded network of genes relevant to renal biology (such as cellular metabolism) and genes somatically mutated in RCC that might be relevant for eoRCC-predisposition (5, 11, 13).

Novel candidate PGVs were stringently defined as those that are predicted to disrupt protein function by the consensus of at least 4 protein predictor algorithms; are rare (gnomAD allele frequency<0.01); and are nonsynonymous variants (frameshifts, stop gains, and splicing) (see Supplementary Methods for variant prioritization). We identified candidate PGVs in 18/22 eoRCC patients in the study, yielding a total of 44 variants in 39 genes (Table 1, and Supplementary Table 2). Gene Ontology analysis confirmed that the candidate PGVs were enriched in DNA repair and replication pathway genes (Figure 1A and Supplementary Figure 1, WebGestalt) (14). Here, 10 patients (46%; 10/22) had 17 candidate PGVs in 14 genes currently associated with hereditary cancers across major organ systems (ATM, BRCA2, POLD1, POLE, FH, MITF, MSH3, MUTYH, PDGFRA, RET, SDHB, SMARCA4, SMARCE1, TSC2). Only 4 patients had candidate PGVs in RCC familial syndrome genes (18%, 4/22 – FH, MITF, SDHB, TSC2). Finally, a total of 14 patients (64%; 14/22) had candidate PGVs from our expanded candidate gene list, from genes not currently defined as RCC-predisposing.

• Download figure
• Open in new tab

Figure 1. Enrichment of predicted pathogenic variants in DNA repair genes in the cohort.

A. Summary of variants in genes and pathways, identified in the cohort. In color - number of variants identified for each gene. For detailed information, see Supplementary tables 1 and 2.

Among the DNA repair-associated genes, candidate PGVs were found in BRCA2 (2/22 patients; 9%, Table 1, Supplementary Table 2) and in ATM (2/22 patients; 9%, Table 1, Supplementary Table 2). In addition, 5 candidate variants in DNA replication-repair genes (4/22 patients; 18%, Table 1, Supplementary Table 2, POLD1 and POLH (Pt #1), POLE (Pt #2), POLK (Pt #16), and RRM2B (Pt #4)). Pt #1 had candidate missense variants in PolD1, a catalytic subunit of the replicative DNA polymerase, Pol δ, and in the translesion synthesis DNA polymerase, Pol η. POLD1 G2275A p.V759I is in a highly conserved region of PolD1 subunit of the Pol δ (coded by POLD1 gene) (15) and occurs at a high allele frequency in the Ashkenazi Jewish population (0.0213) reported in gnomAD (versus 0.0018 in the complete gnomAD dataset (16)). POLH G626T (p.G209V) is in the Pol η catalytic core (17, 18). Pt #2 has a candidate stop-gain G4872A (p.W1624X) variant in the POLE gene, coding replicative DNA polymerase Pol ε, in the conserved C-terminal domain (19). Pt #4 had a splice site (intronic) variant in RRM2B, coding a subunit of p53-inducible ribonucleotide reductase, which performs de novo conversion of ribonucleotide diphosphates into the corresponding deoxyribonucleotide diphosphates for DNA synthesis (20), in genome position #103237248. Pt#16 has a missense variant in the highly conserved N-terminus domain of another translesion synthesis polymerase, POLK G85A (p.E29K); this variant has been described previously as compromising enzyme activity (21). ClinVar classifications of the variants are presented in Supplementary Table 2, with most variants currently classified as variants of uncertain significance (VUS).

Primary lymphocytes from eoRCC patients have reduced capacity to suppress DNA double strand breaks (DSBs)

To begin to assess the functional effect of candidate PGVs in genes linked to DNA replication and repair, we assessed the numbers of γ(phospho)-H2AX foci (a marker of DSBs, (22)) in patient peripheral blood monocytes (PBMCs) at baseline and after treatment with the DNA polymerase inhibitor aphidicolin (Figure 2A). In PBMCs from both matched cancer-free controls (by age and gender) and eoRCC patients, aphidicolin significantly elevated the number of γH2AX foci; however, aphidicolin-treated cells from eoRCC patients had markedly higher numbers of γH2AX foci than those from similarly treated controls on treatment, indicating reduced DSB repair mechanism in eoRCC patient cells (Figure 2A, P<0.001). In complementary work, we tested whether the genes bearing candidate PGVs were specifically needed to suppress DNA DSBs in RCC cells. For this, we used siRNA to deplete the POLD1, POLE, POLH, POLK, RRM2B, and ATM genes in the Caki RCC cell line. For each gene, knockdown significantly increased γH2AX foci relative to control (Supplementary Figure 2) further supporting a role for these proteins in DSB repair in renal cells.

• Download figure
• Open in new tab

Figure 2. Cell-based functional analysis revealed defects in DNA repair and DNA replication in lymphocytes from eoRCC patients.

A. γH2AX foci immune fluorescence staining in primary PBMCs from eoRCC patients versus matched controls, at baseline or post treatment with aphidicolin (2h). PBMCs from patients showed statistically significant elevation of γH2AX foci post treatment with aphidicolin. Data were normalized and are presented as percent of positive γH2AX foci. B. Representative Western blots showing expression of PolD1, Pol η, Pol ε and Pol κ in EBV-transformed cell lines carrying variants versus matched controls (without the variants). Data quantification was performed based on 3 independent biological repeats, technical repeats are presented on gels. C-E. Relative viability of EBV-transformed cell lines was assessed by CTB assay at baseline or after treatment with aphidicolin or UV. Data were normalized to CTB values for controls and are presented as percent cellular viability for POLD1/POLH (C), POLE (D), and POLK (E) cell lines. Data from 3 independent biological repeats are presented. F-G. Difference in DNA replication fork elongation/restoration in EBV-transformed cell lines (F - POLD1, POLH; G - POLE lines) at the baseline and post replications stress was assessed using DNA fiber assay. At baseline the EBV-transformed cells were labeled with IdU for 20 min, for fork restoration cells then were treated with 100 μM aphidicolin or 1 uM MNNG for 2 h, and then labeled with CldU for 40 min. For all conditions, post labeling, cells were lysed, and DNA fibers stretched onto glass-slides, fixed, denatured, blocked, and stained with corresponding antibodies. Fiber images were captured using the Nikon TS2R Inverted Microscope and analyzed in ImageJ software. Data for 3 independent repeats are presented as IdU tract length or CldU/IdU tract length ratio. For all graphs: *** for p<0.001, ** for p<0.01, * for p<0.05 and NS for p>0.05, unpaired, non-parametric t-test, Mann-Whitney criteria.

Patient-derived cell lines with candidate PGVs in DNA polymerases have DNA replication defects

We prepared EBV-transformed cell lines from the primary lymphocytes of 3 patients bearing candidate PGVs (henceforth referred to as the POLD1/POLH cell line, POLE cell line, and POLK cell line) and from several age- and gender-matched cancer-free controls. The POLD1/POLH cell line had significantly reduced levels of the Pol η; for the other candidate PGVs, polymerase level was not affected (Figure 2B). Cell Titer Blue (CTB) cellular assays showed significantly better metabolic capacity or cellular viability than control-derived cell lines when treated either with aphidicolin, or with ultraviolet light (which causes bulky adducts in DNA), suggesting that cell lines from patients had better ability to tolerate DNA damage (Figures 2C-E). Analysis of cell cycle did not show any significant differences in patients and matched control cell lines (Supplementary Figure 3A).

• Download figure
• Open in new tab

Figure 3. Structural and biochemical assays revealed altered enzymatic activities of the PolD1 and Pol η variants.

A. Pol δ complex and primer extension assay. On the left - representative gel image of purified wt and variant Pol δ protein complexes, containing 4 subunits: PolD1 (125 kDa), PolD2 (50 kDa), PolD3 (66 kDa) and PolD4 (12.5 kDa). Center and right – Pol δ complex primer extension assay with quantification. Representative gel image showing reactions performed with 20 nM Cy-3 labeled DNA-duplex template (SA#1), 20 nM of indicated proteins and 500 uM dNTPs. PolD1 V759I complex extended DNA-template less efficiently comparing to wt protein complex. Data for 3 independent repeats are presented. B. Pol η and primer extension assay. On the left - representative gel image of purified wt and variant Pol η catalytic cores (432 amino acids), molecular weight ∼56 kDa. Center and right – Pol η primer extension assay with quantification. Representative gel image showing reactions performed with 20 nM Cy-3 labeled DNA-duplex template (SA#1), 20 nM of indicated proteins and 500 uM dNTPs. Data for 3 independent repeats are presented. C. Pol η lesion (8-oxoG) bypass assay with quantification. Representative gel image showing reactions performed with 20 nM Cy-3 labeled DNA-duplex template with 8-oxoG in the position, opposite to 3‵-OH group (SA#4), and template of the same sequence without lesion (SA#3), 20 nM of indicated proteins and 500 uM dNTPs. Data for 3 independent repeats are presented. For A-C: *** for p<0.001, ** for p<0.01, * for p<0.05 and NS for p>0.05, unpaired, non-parametric t-test, Mann-Whitney criteria. All template sequences may be found in Supplementary Table 4. D-E. Homology modeling structures using yeast protein templates for human PolD1 V759I (D) and for human Pol η G209V (E).

To further expand on these cell-based findings, we used a DNA fiber assay (see Supplementary Methods) and directly compared DNA replication in patient-derived versus control cells, either untreated or following treatment with aphidicolin or with a DNA-alkylating agent, 1-Methyl-3-nitro-1-nitrosoguanidine (MNNG) (Figure 2F-G Supplementary Figure 4A-D). The POLD1/POLH cell line and the POLE cell line exhibited a significantly lower rate of DNA replication in untreated cells (∼1.4-fold decrease, p<0.001 for the POLD1/POLH cell line, and ∼1.9-fold decrease, p<0.05 for the POLE cell line versus controls). We also observed a significantly lower replication fork recovery after 2h treatment with aphidicolin (∼1.44-fold decrease, p<0.001 for the POLD1/POLH cell line, and ∼1.88-fold decrease, p<0.01 for the POLE cell line versus controls) (Figure 2G-H). Intriguingly, the POLD1/POLH cell line showed defective replication fork restoration (∼1.2-fold decrease, p<0.001 versus control line) 2h post-treatment with MNNG (Figure 2G). The POLK cell line did not show any defects in DNA replication and replication recovery under the conditions tested (Supplementary Figure 4A). A complete summary of results for the DNA polymerase variants is provided in Supplementary Table 3.

• Download figure
• Open in new tab

Figure 4. Renal tumors carrying polymerase variants showed high TMB, MSS, and no LOH.

A. Percent alteration frequency in 897 tumors from TCGA in different histological types of RCC: chromophobe (n=66), ccRCC - clear cell renal cell carcinoma (n=538), ccRCC (hyper) - hypermutated samples (n=12), papillary (n=293). B. TMB and MSS data are presented for Pt #1 (POLD1 V759I, POLH G209V) and Pt #2 (POLE W1624X). C. Tumor and normal Sanger sequencing for variants in Pt #1 (POLD1 V759I, POLH G209V) and Pt #2 (POLE W1624X) showing no LOH. Arrows show variants of interest on sequencing tracks.

Altered enzymatic activity of Pol δ and Pol η variant proteins

Among the candidate PGVs detected in polymerases, the Pol κ variant E29K has previously been biochemically shown to possess not only a significantly reduced catalytic efficiency but also reduced replication fidelity (21). E29K is in a conserved region of the Pol κ N-terminus, the N-clasp subdomain (1-32 aa), which is essential to maintaining the stability of the open conformation of the Pol κ active site (23). Intriguingly, a previous study showed that deletion of the first 67 amino acids reduces Pol κ activity during translesion synthesis (TLS, i.e., replication by efficient bypass of bulky lesions in DNA) (24).

To directly test effects of the other candidate PGVs on polymerase activity, we first purified the polymerase delta (Pol δ) protein complex, with the PolD1 (POLD1), PolD2 (POLD2), PolD3 (POLD3), and PolD4 (POLD4) subunits from recombinant protein co-expressed in E. coli, and with preparations containing either wild type (wt) PolD1 or PolD1 V759I variant (Figure 3). Both the wt and the variant-containing Pol δ complexes extended a Cy3-labeled DNA primer-template; however, the V579I variant complex had significantly less robust polymerase activity than the wt complex (Figure 3A, p<0.001). Furthermore, when Pol δ complexes containing PolD1 wt or PolD1 V759I proteins were mixed in a ratio of 1:1, the appearance of the extended primer-template was significantly more robust than the variant alone but significantly less robust than the wt alone. This result suggests that the variant is not only impaired for function but has a partial dominance over the wt in this assay (Supplementary Figure 5, p<0.001).

• Download figure
• Open in new tab

Figure 5. Structure mapping of the novel PolD1 and Pol ε variants from hypermutated ccRCCs in TCGA.

A-E. DNA-bound PolD1 3D-model was refined from PDB:3IAY. The colored functional domains are exonuclease (light blue, residues 131-477) and polymerase (green, residues 550-978). A. Red spheres represent known cancer drivers. B. Blue spheres represent variants of uncertain significance in ccRCC. C-E. Fragments of PolD1 model showing variants: A810S and R823G (C), D893N, A810S, and R978C (D), P151S and E245K (E). F. DNA-bound N-terminal domain of Pol ε was refined from PDB:4M8O. The colored functional domains are N-terminal subdomain (dark grey, residues 31-281), exonuclease (wheat, residues 282-527), polymerase (light pink, palm: 528-950; cyan, fingers: 769-833; lime, thumb: 951-1186). Red spheres represent known cancer drivers (structure above). Blue spheres represent variants of uncertain significance in ccRCC (structure below). G. A fragment of Pol ε model showing variant P696L. H. 3D-model of whole-length Pol ε (without DNA). Structure was refined as described in Methods based on (26). The colored functional domains are N-terminal subdomain (dark grey, residues 31-281), exonuclease (wheat, residues 282-527), polymerase (light pink, palm: 528-950; cyan, fingers: 769-833; lime, thumb: 951-1186), C-terminal domain (light grey, residues 1308-2222). Blue spheres represent variants of uncertain significance in ccRCC. I. Fragments of Pol ε model showing variants S803L and F753.

Pol η is a low fidelity polymerase, which contributes to its ability to perform TLS (25). Hence, G209V variant and wt Pol η prepared in E. coli were assessed for their ability to extend labeled DNA primer-template duplexes (Figures 3B-C). In the absence of DNA damage (e.g., in a normally base-paired template), the wt and variant proteins both extended the template (Figure 3B, p>0.05), but the observed bands suggest higher processivity for the variant on template without lesions compared to wt (Figure 3B). To evaluate repair of DNA damage, TLS activity was also tested using a template containing an 8-oxoGuanine (8-oxoG) DNA lesion (25). Pol η wt bypassed the 8-oxoG lesion robustly compared to the Pol η variant (p<0.05, Figure 3C), suggesting better processivity for the wt protein on template with DNA lesion (25).

Finally, biochemical analysis of a purified Pol ε variant, W1624X, was not performed as it is a stopgain variant in the C-terminal domain or CTD, truncating 662 amino acids of the protein. The CTD region is not well-studied, but is thought to be essential for stability of the Pol ε holoenzyme (26).

Structural modeling of DNA polymerase variants in eoRCC suggests impact on polymerase function

The PolD1 V759I variant is targets a residue two amino acids away from residue D757 (Figure 3D, in dark green). In the PolD1 active site, D757 coordinates the Mg²⁺ ions (neon green spheres) required for DNA synthesis and plays a direct role in the catalytic mechanism and binding of DNA (15). Structural modeling indicates that a substitution of the valine (V) 759 to isoleucine (I) could plausibly alter the position of D757 and disrupt the efficiency of DNA polymerization. To further understand the structural changes that each polymerase variant might induce, we calculated the change in stability of each amino acid substitution as described in Supplemental Methods and Supplementary Table 4. Interestingly, the PolD1 I759 yielded a mild stabilization (ΔΔG of -1.36 kcal/mole) relative to the wildtype V759. The I759 residue makes twice the number of hydrophobic contacts as the wildtype V759 with the long helix below. The variant could lock this strand in an overly rigid position that compromises DNA polymerization steps that involve flexibility, consistent with the biochemical results observed in Figure 3A.

The Pol η G209 residue is in the catalytic core of the polymerase (Figure 3E, residues 1-432, colored green) at a position often called the C-cap, i.e. the residue in this position is proximal to the C-terminal end of the α-helix (in orange) (27). Typically, valine, threonine, and isoleucine residues are not preferred in the C-cap, due to poor solvation at the C-terminus of the helix when the side chains are bulky (27). Structural modeling of the G209V substitution showed that the valine with a bulkier side chain could not only alter the stability of the α-helix, but also the nearby β-strands of the catalytic active site (in pink). Rosetta modeling shows the variant G209V is capable of making 5 hydrophobic contacts across the cleft with R24, which is just downstream from the key residue D13, which coordinates active site metals that are central to the catalytic mechanism and the binding of the incoming NTP. The predicted change in stability of the Pol η G209V variant revealed a significant destabilization (ΔΔG of 5.67 kcal/mole) relative to the wildtype G209, consistent with lower levels observed in protein purification (Figure 2B, Figure 3B (left gel), and Supplementary Table 4). Interestingly, the activity of the G209V compares well with wt for a normal primer template and might even be more processive (Figure 3B, middle gel). However, the variant appears defective for TLS when the template contains an 8oxoG (Figure 3C). Consideration of these data together show that changes in stability and conformation may be subtle and even result in an alteration of substrate preference. Thus, a careful combination of both computational and biochemical methods is required to gain a clear understanding of the role a given polymerase variant might have in DNA replication, repair, and cancer initiation and/or progression.

EoRCCs carrying candidate PGVs in DNA polymerases are hypermutant and microsatellite stable (MSS)

To extend these functional tests, we next explored tumor mutation burden (TMB) in tumors from RCC patients in TCGA and from the FCCC eoRCC patients. Previous studies have shown that colorectal and endometrial tumors carrying mutations in POLE exonuclease domain (ExoD) and in POLD1 exhibit a high burden of mutations, are typically MSS but few cases with microsatellite instability (MSI) have been reported, and do not exhibit loss of heterozygosity (LOH) (15, 28-36). RCCs are typically non-hypermutated, with an average TMB of ∼1 mut/Mb (37); however, rare hypermutated (≥10 mut/Mb) and rarer ultra-hypermutated (>100 mut/Mb) RCCs carrying polymerase mutations, with or without MSI, have been reported (37). Analysis of renal tumors in TCGA found that several candidate PGV genes from this study (Figure 4A), including specifically DNA polymerase genes were mutated in these hypermutant clear cell RCCs (ccRCCs) s (Figure 4A). We analyzed the TMB and MSI/MSS status of tumors from the FCCC eoRCC patients and found hypermutation without MSI (POLD1/POLH tumor with 12.85 mut/Mb; POLE tumor 14.44 mut/Mb) (Figure 4B). The tumors from the FCCC eoRCC patients did not exhibit LOH of the polymerase genes, as has been observed in other polymerase-mutated tumors (37, 38) (Figure 4C).

To expand the analysis of DNA polymerases in RCC, we next modeled the structural consequences of the PolD1 and Pol ε variants in hypermutated ccRCCs from TCGA. Supplementary Table 5 shows the predicted changes in stability for 20 different PolD1 variants and 23 Pol ε variants from hypermutated ccRCCs in TCGA. A broad range of both stabilizing and destabilizing variants was found in all domains of each of these polymerases. Figure 5 shows these variants in relation to known pathogenic variants in POLD1 (Figure 5A-E) and POLE (Figure 5F-I). PolD1 R823G/L is in a β-sheet of the polymerase domain close to the DNA (Figure 5C). The substitution of a positively charged arginine (R) to a hydrophobic glycine (G) or leucine (L) could destabilize the β-sheet and impact DNA binding. PolD1 D893N is positioned close to the DNA and a variant in this region may destabilize DNA binding or position for DNA-protein interactions (Figure 5D). PolD1 P151S is in a β-sheet in the ExoD and the change from proline (P) to serine (S) could destabilize the β-sheet geometry (Figure 5E). P151 is close to a known cancer driver mutation (37, 39), E245K, in the unfolded region of the ExoD (Figure 5A).

Pol ε P696R is in the palm region of the polymerase domain, which is highly conserved among replicative polymerases (Figure 5G). Arginine (R) has a very large positively charged side chain when compared to smaller proline (P), suggesting this variant may disrupt the polymerase structure and impact DNA synthesis. Pol ε S803L and F753L are in the flexible region of the polymerase domain or the fingers (in cyan, Figure 5I). This finger region shifts (27° tilt) on DNA binding (40), and thus plays an essential role in polymerase function. S803 is close to the positively charged lysine (K) 431 in the ExoD, and serine (S) is polar and a smaller residue than the hydrophobic leucine (L). Pol ε S803L is near a site of known cancer driving mutations, C810 (37), suggesting that specific alterations in this α-helix could impact polymerase functioning. F753L is on the border of the fingers and the polymerase domain, close to the ExoD and could be important in the coordination of these regions with or without DNA binding. Finally, several variants were found in the C-terminal domain (Figure 5H, in light grey), which is currently not well-studied, but is known for stabilizing the Pol2 (human Pol ε) complex in yeast (26).

eoRCC patient population, and peripheral blood DNA analysis

Case-only eoRCC probands that underwent clinical germline genetic testing between 2010-2016 were included in this study (n=22). Patients were followed by the Genitourinary Program at the Fox Chase Cancer Center and had undergone evaluation for inherited cancer risk at the FCCC Family Risk Assessment Program (RAP). Each participant had a strong family cancer history as shown in Table 1, with either multiple first-degree or second-degree relatives with RC, RC-associated cancers, or other cancers. The mean age at eoRCC diagnosis was 48 years (range 36–59 years). No pathogenic mutations were identified from sequencing the following RC-specific genes: VHL, MET, FLCN, TSC1, TSC2, FH, SDHx, PTEN and BAP1. The patients had consented to the FCCC RAP Registry, which allowed further research, genomic sequencing, and banking of their biospecimens in the BioSample Repository Facility. The patients reported here were self-reported white, non-Hispanic. Family histories were obtained by trained licensed genetic counselors and verified by attending physicians. The analyses performed and publication of deidentified information was under the approval of the FCCC Institutional Review Board Committee protocol number 14-831. See Supplementary Methods for more methods.