If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Center for Molecular Medicine and Genetics, Wayne State University School of Medicine, Detroit, MI, United StatesDepartment of Oncology, Wayne State University School of Medicine, Detroit, MI, United StatesMolecular Therapeutics Program, Karmanos Cancer Institute at Wayne State University School of Medicine, Detroit, MI, United States
Extensive in vitro analysis of several variants of unknown significance in ATM reveals several likely pathogenic variants.
•
Results of in vitro assays differ from predictions made via in silico algorithms.
•
Provides a method for identification and assessment of rare SNPs in hereditary cancer syndromes.
•
Wet lab functional assays provide useful data in determining individual SNP pathogenicity.
Abstract
Variants of unknown significance (VUS) remain a constant challenge in the diagnosis of hereditary cancer and the counseling of patients with pedigrees suggestive of such a syndrome. In order to assess some of this limitation, several variants in the DNA repair gene ATM were selected from a cohort of high risk individuals with negative genetic diagnoses. ATM has proven a challenge in the counseling of patients due to its nature as a moderate penetrance gene. In this study, six ATM missense mutations with a high likelihood for pathogenicity were assessed through a battery of experiments to yield high fidelity information on their biochemical effect on ATM activity. We report that several of these variants show signs of reduced ATM function indicative of likely pathogenicity. With further study, this data may be used in clinic, improving diagnosis, surveillance, and outcome for patients carrying these mutations.
In 2021, an estimated 21,410 individuals will be diagnosed with ovarian cancer in the USA. In addition, 13,700 patients with ovarian cancer will die from the disease [
]. The five year survival rate for ovarian cancer is currently 47%, a number reflecting high survival rates for early stage disease (up to 92% for stage 1A and 1B). However, unlike breast cancer, few ovarian cancers are diagnosed at early stage. Only 15% of patients present early, and the survival rates decline rapidly as the disease becomes more advanced [
The pattern of disease presentation can be addressed by early surveillance and prophylactic surgery for at risk patients but is not recommended for the general population, or even many patients at moderate risk [
]. Cancer genetic testing can play a key role in identifying patients for whom early surveillance and enhanced testing can be beneficial. This, in turn can drastically increase the survival of the most vulnerable patients by enhancing cancer detection at early stages when it is most treatable [
Multigene Panel Testing Detects Equal Rates of Pathogenic BRCA1/2 Mutations and has a Higher Diagnostic Yield Compared to Limited BRCA1/2 Analysis Alone in Patients at Risk for Hereditary Breast Cancer.
]. Panel testing has been able to identify pathogenic mutations in 17% of BRCA1/2 negative high risk families and 24% of unselected ovarian cancer patients [
]. The results of gene panel testing, when unambiguous, identify at risk individuals and change their clinical management, enabling increased surveillance [
]. Despite this success, at least 49% of the heritable component of hereditary breast and ovarian cancer remains unknown, with many seemingly familial cases testing negative with the current panels [
]. Failure to offer panel testing has led to the underdiagnoses of HBOC in Medicare and Medicaid patients, for whom insurance does not cover most testing panels and methods [
]. Overall, there is an increasing trend toward the discovery of new cancer predisposing mutations, and the inclusion of these on testing panels would further increase the panel's utility in the clinic [
The protein kinase ATM plays an essential role in the repair of double strand breaks, a mechanism often associated with HBOC. Along with the related protein, ATR, ATM is a key activator of the DNA damage response, and the cellular response to other types of stress, including oxidative stress and reactive oxygen species, replication fork stress, telomere dysfunction, hypoxia, and hyperthermia, phosphorylating over 150 target proteins leading to widespread pathway activation, chromatin relaxation, and nucleosome remodeling [
]. Upon DNA damage, auto-phosphorylation of ATM on S1981 in the FAT domain causes the inactive multimer form to disassociate, and ATM is stabilized in its active monomeric form [
]. Additional phosphorylation occurs on S367 and S1893 along with acetylation of K3016 to further stabilize active ATM, which then interacts with its numerous targets [
ATM mutations have long been implicated in the development of breast cancer, conferring a 60% risk by age 80 to individuals carrying germline mutations [
]. Damaging ATM mutations are present in the heterozygous state in up to 2% of the general population, conferring a 2–4 fold risk increase for general cancers, and a 5–9 fold risk increase for cancers in women [
]. ATM is hypothesized to be a significant component of genetic risk specifically for hereditary breast and ovarian cancer syndrome (HBOC), and mutations have been found in many families with HBOC even as it presents a challenge to genetic counselors as a moderate penetrance gene [
]. The increasing trend of ATM mutation discovery in HBOC is strongly suggestive of a role in the development of cancer, and individual ATM mutations may confer moderate risk for these cancers. These women, therefore, undergo enhanced screening, and they may be more responsive to PARP inhibitors for cancer treatment [
], we examined, using whole exome sequencing, a cohort of 48 women with a personal history of ovarian cancer and a significant family history of cancer but who were negative for pathogenic mutations in BRCA1 and BRCA2. Hypothesizing that variants of unknown significance in genes related to DNA repair and cell cycle control may account for their missing heritability, we found several single nucleotide polymorphisms (SNPs) classified as variants of unknown significance with a moderate to high predicted impact on function in our cohort. In addition, five families were re-contacted for clinical genetic testing and counseling based on new relevant mutations detected by our study. The majority of the SNPs in our cohort were found to be low frequency mutations of uncertain significance and poorly annotated in both literature and clinical databases [
]. In accordance with our hypothesis and data, others have found that retesting patients who were previously negative when assessed for BRCA1/2 specific mutations yields important information about risk, pointing to the further relevance of panel testing in the clinic [
Because of the likelihood of hereditary cancer in our patient cohort, we hypothesized that these SNPs could represent pathogenic variants. In our cohort, several missense mutations in ATM were found to be significantly overrepresented (Table 1) in addition to two rare truncations in ATM [
]. All mutations identified by NGS whole exome sequencing were validated by Sanger sequencing. Here, we report the in-depth functional analysis of these variants. The analysis of these variants demonstrates that by combining functional studies with small curated cohorts, our process can facilitate the discovery of variants likely responsible for missing heritability.
Table 1Variants of unknown significance in ATM from a cohort of high risk ovarian cancer. Variants are listed in amino acid order from N to C-Terminus with rsID. All variants have a low minor allele frequency (MAF) and are listed in ClinVar as conflicting interpretations of pathogenicity. REVEL scores are shown, which combine several variant classification algorithms into one higher fidelity score [
]. REVEL scores are closer to 0 than 1, indicating that most algorithms predict that these variants are unlikely to be pathogenic. Scores from the EVE algorithm
Cohort and IRB: The sample consists of DNA from 48 BRCA1/2 mutation negative Caucasian OVCA patients from 47 families (one mother-daughter pair). Informed consent was signed and permission was obtained for the collection of blood samples and for access to medical records for all subjects. The protocol (HIC#024199MP2F(5R)) was approved following Full Board Review by the Human Investigation Committee at Wayne State University, Detroit, Michigan.
ATM variants: pcDNA3.1(+)Flag-His-ATM WT was a gift from Michael Kastan (Addgene plasmid # 31,985; http://n2t.net/addgene:31985 ; RRID:Addgene_31,985) [
]. ATM variants were constructed using the NEB Q5 site directed mutagenesis kit. Following completion of the PCR, ATM variants were transformed into STBL2 Max Efficiency competent cells (Invitrogen), and grown for 2 days on LB agar + 100 µg/mL Carbenicillin at 30°C. Small colonies were selected and transferred to liquid culture in LB broth + 100 µg/mL Carbenicillin, and grown for 2 days at 30°C. Presence of the desired variants were confirmed by Sanger sequencing (Genewiz), and clones carrying the correct variant sequence were then assessed by full plasmid Sanger sequencing via several primers resulting in overlapping 500 bp high quality regions. Clones with no other alterations except the VUS were used for subsequent experiments. ATM-S1981A was expressed from hATMS1981A, which was a gift from Michael Kastan (Addgene plasmid # 32,300 ; http://n2t.net/addgene:32300 ; RRID:Addgene_32,300) [
]. The same cell line was used for all experiments. Quality control for mycoplasma was completed once every 6 months throughout the experimental workflow. Cell line authentication was performed by the Biobanking and Correlative Sciences Core Lab of the Karmanos Cancer Institute.
Lentiviral shRNA knockdown: ATM knockdown was accomplished using shRNA expressed from Lentiviral vectors targeting the 3′UTR (Sigma-Aldrich, TRCN0000010299, ready to use lentiviral particles). Hela-DRGFP cells were transduced with lentivirus at a MOI of 2 on the day of plating for experiments. Freshly thawed lentiviral supernatants were combined with counted cells in suspension and then plated as one unit. Knockdown efficacy of ∼.6 was confirmed by Western Blot. Verification of ATM knockdown and ATM VUS models is shown in Supplemental Figure S-1 and Fig. 1d-E.
Fig. 1ATM variants in HDR analysis. A: Diagram of HDR analysis in Hela-DRGFP. Treatment with iSce1 creates a double strand break, which, if repaired by HDR, results in the expression of GFP from a now functional allele. B: Comparison of GFP expression after pCMV3 (empty vector) and pCBASce1 transfection in Hela DRGFP. C: Six ATM variants were assessed for fluorescence after treatment with iSce1. S1891A is a known pathogenic variant lacking kinase activity. pATM-WT expressed wild-type ATM and was used as a positive control. Sce1 only shows baseline Hela-DRGFP fluorescence and is used for calibration of WT-ATM rescue. The group pCMV3 shows baseline Hela fluorescence without iSce1 treatment. All conditions receiving ATM shRNA (shATM) also receive iSce1. Assessment is completed using FACS (BD-LSR II). Statistics use alpha = 0.05, 1 tailed t-test, assuming unequal variances, comparisons made to pATM WT. D: Relative expression of ATM as a ratio to vinculin assessed via western blot and normalized to mock transfection with pCMV3 (loading control). E. Western Blots showing ATM expression and Vinculin.
Transfections: Plasmid transfections were carried out using JetPrime (Polyplus transfection) according to manufacturer's specification. For a 12-well plate, 70,000 Hela-DRGFP cells were seeded 1 day prior to transfection and checked for appropriate patterning before transfection. Transfections were carried out using 1 µg of plasmid (pCMV3 for control, or pATM-Flag-His (variant or WT)), 100 µL JetPrime buffer, and 4 µL of JetPrime per well. For HDR assays, an additional 500 ng or iSce1 or pCMV3 was included. This formula was scaled for kinase assay, colony survival assay, and protein harvesting for western blotting.
Western Blotting: For ATM, protein lysates were collected on ice using RIPA buffer with 1.5% PITC and 1% SDS. Samples were incubated in RIPA 15 min on ice before shearing with a 21 ga. needle. Following shearing, samples were incubated on ice for 1 hour, then centrifuged at 15,000 rpm at 4°C for 15 min. Protein concentration was assessed via Bradford assay. Two hundred µg of each sample were loaded for gel electrophoresis, which was carried out using 6% polyacrylamide gels at 120 V x 20 min + 220 V x 45 min. Transfer was accomplished in Bjerrum Schafer-Neilsen transfer buffer to nitrocellulose at 180 mA for 18 hrs in cold room. Membranes were air dried for 15 min at room temperature to fix protein, then blocked in 5% milk-TBST overnight at 4°C. After 3 × 10 min wash in TBST, ATM primary antibody (1:500 ATM-AM9, Millipore-Sigma, 05–513) in 5% BSA-TBST was used, and incubated overnight at 4°C. Following another 3 × 10 min TBST wash, anti-mouse secondary (AF 680, Abcam, ab175774) was incubated for 1 hour at room temperature in 5% BSA-TBST. After a final 3 × 10 min TBST wash, imaging was carried out using a Li-cor Odyssey imager using Odyssey software (Li-Cor) for visualization and quantification. Vinculin was used as a loading control and was probed after ATM imaging to avoid oversaturation. The membrane was incubated with vinculin primary antibody (Invitrogen, 700,062) for 2 h at room temperature in 5% BSA-TBST, washed as before, then incubated for 1 hour at room temperature with anti-rabbit IgG secondary antibody (AF-790 DxR, Invitrogen, A11374) in BSA-TBST, washed again, and imaged on the Licor. p-CHK2-T68, p-p53-S15, and p-H2AX-S139 were probed concurrently. For these proteins, 100 ug of lysate was loaded onto a 12% polyacrylamide gel, then transferred to nitrocellulose membranes. After 15 min air drying, the membranes were blocked in 5% BSA-TBST for 1 hour, then probed for 2 h in 5%BSA-TBST with three antibodies (P-CHK2 Thr68 [Cell Signaling, 2661S] at 1:1000, P-P53 Ser15 [Cell Signaling, 9284S] at 1:1000 and P-H2AX Ser139 [Cell Signaling, 2577S] at 1:800) for 2 h at room temperature. After wash, anti-rabbit IgG secondary was used to probe for 1 hour at 1:10,000 in 5%BSA-TBST. The membranes were then washed and visualized on a Li-Cor Odyssey imager. The membranes were then stripped and re-blocked overnight at 4°C in 5% BSA-TBST and probed for loading control (beta actin, (anti-ACTB, Sigma, HPA041264) as before. pBRCA1 was blotted from a 6% gel run and transferred as ATM above. The membrane was blocked overnight at 4°C in 5% BSA-TBST, then probed in 5% BSA-TBST using pBRCA1-S1387 (Bethyl, A300–007A) at 1:1000 and vinculin at 1:10,000 at 4°C overnight. Following wash, the membrane was incubated with the anti-rabbit IgG secondary antibody for 1 hour at room temperature, washed again, and imaged on the Li-Cor Odyssey as above. p-p53-S20 and p21 were probed using a 12% gel as above for p-CHK2, p-p53-S15, and p-H2AX. The membrane was incubated with anti-p-p53-S20 (Cell Signaling, 9287S) in 5% BSA-TBST overnight at 4°C after blocking, then washed, probed with anti-rabbit secondary antibody and imaged. The membrane was then stripped, re-blocked for 1 hour using BSA-TBST, and probed for actin and p21 (Cell signaling, 2947S) for 2 h at room temperature, after which it was washed, probed with anti-rabbit IgG secondary antibody, washed again, and imaged as before.
CDC25A-S178 and p-CDC25C-S216 were probed on separate 12% gels as above. Each gel was incubated separately with either CDC25A-S178 antibody (Bioss, bs-3095R) or CDC25C-S216 antibody (Cell Signaling, 2577S) for 2 h at room temperature after blocking in BSA-TBST. Anti-rabbit IgG secondary antibody was used for both, and actin was used as a loading control.
Kinase assay: Hela cells were transfected with all ATM variants and controls as described above. Three days post transfection, 1 plate of cells for each group was treated with 2.5 µM Bleomycin (Selleck Chemical), with a duplicate kept as a control. After 4 h incubation, protein was purified from both plates and used for western blotting as described above.
HDR assay: To obtain relative GFP expression in experimental groups, Hela-DRGFP cells were plated on day 1 at 70k/well in a 12 well plate, and transduced with lentiviral shRNA as described above. On day 2, transfection was performed using JetPrime as described above to introduce pCMV3, iSce1, and pATM variants. On days 3 and 4, the medium in the wells was changed to remove dead cells. Finally, on day 5 the cells were harvested in PBS for use in Flow Cytometry. Flow cytometry was performed using the Wayne State University Flow Cytometry Core as described before [
], except using a BD-LSR instrument. Gates were used to select for singlets, GFP+ (FITC), and DAPI+ (live) cells.
Colony Survival Assay: Hela Cells were transfected as described above. Two days post transfection, cells were re-suspended and counted, and then 300 cells were seeded onto each well of 6 well plates. A total of 5 concentrations with three replicates each were plated for each condition. One day post plating, cells were treated with bleomycin (0–10 µM) for 4 h. Media was then replaced with standard culture media and the colonies were allowed to grow for 7 days. After 7 days, colonies were fixed and stained using crystal violet. Counting was performed with a dissection miscrscope. Colonies of at least 20 cells were scored, each as 1 colony. IC50s and 95% confidence intervals were calculated using GraphPad Prism. Survival was calculated after controlling for plating efficiency.
Results
Assessment of ATM variants in homology directed repair
ATM variants were initially assessed via homology directed repair (HDR) assay to provide a global view of decreased function in HDR. To model the ATM VUS, Hela DRGFP cells were transduced with lentiviral shRNA targeting the 3′UTR of ATM to perform knockdown. ATM cDNA from plasmids expressing either wild-type ATM or one of the six VUS was then used to restore ATM expression to previous levels (Fig. 1). An empty vector control was included to visualize HDR activity in unaltered Hela cells and to verify restoration of ATM signaling to wild-type levels. Transfection with iSce1 was performed to cut the DRGFP signaling plasmid was done, and the relative HDR activity of the models was then assessed after three days by flow cytometry using GFP expression as a marker for successful HDR. A negative control with no iSce1 was included and used to correct for any background HDR activity at the DRGFP allele. Restoration of complete HDR function was observed by the transfected pATM-WT plasmid in comparison to the unaltered control (pCMV3). The ATM knockdown control (shATM+ pCMV3), modeling ATM dysfunction, showed marked loss of HDR activity (p = 0.00067 vs. pATM-WT). Several of the variants we tested exhibited defects in HDR function compared to the wild type plasmid, with two variants exhibiting HDR function equivalent to the known pathogenic control allele, S1981A (p = 0.00078 vs. pATM-WT). S49C showed a significant decrease in HDR function (p = 0.03), but this variant likely retains some HDR function (see discussion below). S333F and D1853V were similar in exhibiting a significant decrease in HDR activity (p = 0.007 and p= 0.0009, respectively), but retain some HDR function. Some mechanism conferred by these mutations likely reduces the efficiency of ATM's kinase activity. F1463C was interesting in its appearance in our previous study in three patients from unrelated families. Thus, it was a high priority for analysis. This variant possessed lower HDR activity than the known pathogenic control, suggesting that it causes a near complete loss of function to the affected allele (p = 0.00003 vs. pATM-WT). Similarly, V2540I exhibited a near complete loss of HDR activity in this assay (p = 0.000014 vs. pATM-WT)). Because they showed such a large decrease, these two VUS are of prime interest as potential pathogenic variants Finally, L2307F showed complete retention of function, with levels of HDR activity similar to the wild-type control allele (p = 0.18). These data strongly suggest that there is severe loss of ATM function for several of the variants present in our high risk patient cohort.
Effect of atm variants on colony formation
Successful repair of DNA damage increases the cell's ability to survive mutagenic stimuli, including DNA damaging drugs commonly employed in chemotherapy. We used the radiomimetic bleomycin [
Ionizing radiation induces ataxia telangiectasia mutated-dependent checkpoint signaling and G(2) but not G(1) cell cycle arrest in pluripotent human embryonic stem cells.
] (Fig. 2). To ensure comparable data, the same cell model system was used as in the HDR assay, but without transfection of the iSce1 plasmid to ensure that double strand breaks were strictly proportional to the amount of bleomycin used. As before [
], the resulting variant models were first assessed for plating efficiency, a baseline measure of their ability to form a colony from a single cell. Plating efficiency can be reduced in cells with dysfunction of critical genes, indicating relative fragility. The assessment of plating efficiency is also used to control for the effect of the drug treatments on each variant model. Here, the ATM knockdown cells and the cells expressing the kinase deficient control variant S1981A showed significant reduction in their viability after plating but without any treatment with bleomycin (p = 0.0003 and p = 0.0004, respectively). It is noteworthy that several of the variants showed significant reductions in viability after plating when compared to the ATM wild type model, although all of these were similar or above the plating efficiency of the empty vector control transfected cells (F1463C p = 0.01, D1853V p = 0.025, L2307F p = 0.003, V2540I p = 0.004). Bleomycin IC50 was assessed from the colony survival assay data, and none of the VUS were found to differ significantly from the wild type model in this metric, although the knockdown model did exhibit significantly reduced IC50 (p= <0.05, IC50 = 1.11 [95% CI 0.82 – 1.45] vs. IC50 3.51 [95% Ci 2.34–5.5] for ATM-WT). Because there was widespread cell death at higher concentrations, we considered whether the surviving percentage at the mid-point concentrations of bleomycin (1.5 and 2.5 µM, before high lethality is uniformly present) might be more illuminating for specific defects, even if the overall IC50 did not vary significantly between variants. In this analysis, it was found that at both concentrations of bleomycin, the ATM knockdown model showed significantly reduced survival compared to the wild type model, which in turn did not significantly vary from the empty vector control, indicating that the overall experiment was sensitive enough to detect a defect in ATM signaling, and that the wild type plasmid did fully restore ATM function, rendering analysis of these concentrations useful. At 1.5 µM bleomycin, the variants V2540I and F1463C showed significantly reduced survival (p = 0.046, and p = 0.033, 47% and 46% survival relative to plating, respectively). Similar early reduced viability was present in the ATM knockdown cells and those expressing S1981A (p = 0.02 and p = 0.008, 37% and 32% survival relative to plating, respectively). The reduced viability of the F1463C was even more apparent at the higher concentration of 2.5 µM bleomycin (p = 0.002 vs. ATM-WT, 13% survival relative to plating) again comparable to the ATM knockdown cells (p = 0.003 vs. ATM-WT, 11% survival relative to plating). These results are notable in light of the severe reduction in HDR activity seen for these two variants. The fact that they are also the least likely to recover from more mild doses of bleomycin compared to the ATM WT control indicates that these ATM variants may be unable to repair minimal DNA damage as efficiently, thus cells carrying them may be more susceptible to accumulating unrepaired mutations as a result of fewer double strand breaks. Such inability to repair and recover from smaller amounts of DNA damage could cause the elevated genetic risk seen in HBOC families.
Fig. 2Colony Survival Analysis of ATM Variants. A: Plating efficiency as a percent of the control group, pCMV3, for each variant model. B. IC50 for each variant model calculated using GraphPad Prism Software using the [inhibitor] vs. normalized response (variable slope) model. Significance calculated by GraphPad for 95% confidence. C-D: Percentage survival for each model group at the specified concentration of Bleomycin, as expressed as a percentage of plating efficiency. For all panels, N = 3, statistics use alpha = 0.05, 1 tailed t-test assuming unequal variance, comparisons made to pATM-WT. .
The major biochemical function of ATM in the DNA damage response is as a protein kinase. Upon activation, ATM phosphorylates and modulates the activity of over 150 unique targets, many with redundant activation from different pathways [
]. In order to investigate a potential mechanism by which our variants of interest may disrupt the DNA damage signaling cascade, we analyzed eight such targets. These targets were BRCA1 at S1387 [
Fig. 3Analysis of the Kinase Activity of ATM Variants. A: Representative western blots of all phosphorylation targets analyzed and loading controls. Blots are grouped by membrane, with loading controls present for each separate membrane-target. B: Graph of relative induction for each target. Significance uses alpha =0.05, 1 tailed t-test assuming unequal variants. .
Fig. 3Analysis of the Kinase Activity of ATM Variants. A: Representative western blots of all phosphorylation targets analyzed and loading controls. Blots are grouped by membrane, with loading controls present for each separate membrane-target. B: Graph of relative induction for each target. Significance uses alpha =0.05, 1 tailed t-test assuming unequal variants. .
Anova analysis of our measurements of phosphorylation of BRCA1-S1387 showed significant differences between variants (Single Factor Anova, p = 0.01, F = 3.22, Fcrit=2.39). Compared to the ATM-WT control, both ATM knockdown and ATM-S1981A exhibited markedly reduced phosphorylation in response to DNA damage (p = 0.046 and p = 0.015, respectively). Replacement with the VUS V2540I showed a reduction in induced phosphorylation similar to the negative controls (p = 0.037). Surprisingly, the L2307F variant showed a notable and significant increase in phosphorylation of BRCA1-S1387 (p = 0.021). Visible but non-significant increases were also observed for cells expressing the variants S49C and F1463C (p = 0.26 and p = 0.15, respectively). This target on the BRCA1 protein is shared with two other kinases, and the resultant increase may reflect compensatory activation of those kinases (ATR and DNA-PKcs) in response to changes in signaling by this variant. Unfortunately, we were unable to acquire antibodies for the more ATM specific residues in BRCA1, S1423 and S1524 [
Induction of p-H2AX-s139 was significantly reduced in the ATM knockdown cells, and those expressing the S1981A and S333F variants (p = 0.04, p = 0.035, and p = 0.027, respectively) when compared directly to ATM-WT. However, these comparisons risk error in that the overall dataset didn't register as significant by Single Factor Anova (p = 0.17).
Analysis of p-p53-S15, likewise, did not register as significant overall by Anova (p = 0.056). However, with the much narrower gap to significance, we think that the significant pairwise comparisons hold more weight. For this target, ATM knockdown, S1981A, and V2540I were significant (p = 0.038, p = 0.032, and p = 0.032, respectively). Additionally, phosphorylation at this site was notably reduced for S49C (p = 0.051), and D1853V (p = 0.07).
A major phosphorylation target for ATM is CHK2-T68. We tested p-CHK2-T68 for our variants, which showed significant differences among variants in activation (Single Factor ANOVA, p = 0.000799, F = 543, Fcrit =2.39). The ATM knockdown cells, and those expressing S1981A and S333F variants showed significant decreases in induction (p = 0.022, 0.026, and 0.041 respectively). Additionally, compared to ATM-WT, phosphorylation by the variant ATM S49C was decreased but not significantly (p = 0.057) as was phosphorylation by the variant V2540I (p = 0.056).
No differences were detected among the experimental conditions in p-p53-S20 induction (Single Factor ANOVA, p = 0.72, F = 0.67, Fcrit=2.39), including the variants affected for CHK2 phosphorylation above, indicating that the reduction in p-CHK2 protein kinase activity for S333F was not sufficient to register any change in p-p53-S20 induction.
Given its central role in the DNA damage response of the activation of TP53, we investigated whether any of the ATM variants affected the expression of p21. (Despite expressing the HPV E6 protein, HeLa cells do express detectable, functional p53 [
Repression of endogenous p53 transactivation function in HeLa cervical carcinoma cells by human papillomavirus type 16 E6, human mdm-2, and mutant p53.
]) The induction of p21 was assessed and found not to differ significantly among any of the variants (Single factor Anova, p = 0.22, F = 1.48, Fcrit=2.39). Visually, more p21 accumulation was present, though not by a significant factor in the ATM knockdown cells, and those expressing the S1981A, S49C, and V2540I variants (p = 0.17, 0.13, 0.14, and 0.23, respectively). It is possible that this is the result of more unrepaired double strand breaks in these models, causing more cell cycle arrest as reflected in the accumulation of p21. However, accumulation of p21 is not currently known to be evidence of pathogenicity for ATM variants.
Considering that some of the variants may affect CHK2 signaling, we also investigated CDC25A and CDC25C phosphorylation by the variants S49C, V2540I, and F1463C (Supplemental figure S-2, S-3). While there was no difference in CDC25C phosphorylation, CDC25A phosphorylation was significantly reduced for cells expressing both variants S49C and V2540I, but not for the ATM knockdown cells with the empty vector, a surprising result.
Discussion
In this study, we assessed the functional deficiencies conferred by the ATM variants S49C, S333F, F1483C, D1853V, L2307F, and V2540I (Table 2), all of which currently have conflicting interpretations of pathogenicity. Compared to our negative controls (ATM knockdown + empty vector and the kinase deficient variant ATM-S1981A), which showed broad defects across assays, we found that the variants S333F, F1463C, D1853V, and V2540I exhibited markedly reduced activity. Cells expressing S333F intriguingly were defective in HDR and also had reduced phosphorylation of two targets in response to DNA damage, suggesting a potential mechanism for this variant. By contrast, the variant L2307F showed almost entirely wild type activity, with reduction only in plating efficiency during colony survival assay, but a surprising increase in BRCA1 phosphorylation. S49C showed an intermediate phenotype, with S49C expressing cells exhibiting a significant reduction in HDR activity only, while many phosphorylation targets had defects approaching significance. Interestingly, despite low to intermediate REVEL scores for all variants (Table 1), the in silico assessments both alone [
] and combined failed to identify likely damaging variants. Rather than assessing the effects of many variants against a panel of control variants in a high throughput study, we studied the individual effects of several variants in detail using specific in vitro biochemical assessments of the variant activity. Multiple replications of each assay combined with the ability to then compare the variant effects across assays and enable the assessment of variant phenotype with high confidence.
Table 2Comparison of ATM Variant activity in functional assays showing defects by significance. Values approaching significance are listed as numerical p values. S1981A is a known pathogenic variant, and shATM+pCMV3 represents a knockdown of ATM only. Both negative controls show defects across assays. Comparisons are made to ATM-WT, thus ATM-WT control is not shown. Magnitude effect is shown for measurements at or near significance and represents a fraction of the measurement for ATM-WT. ns = not significant. .
Buchholz, T.A., M.M. Weil, C.L. Ashorn, E.A. Strom, A. Sigurdson, M. Bondy, et al., A Ser49Cys variant in the ataxia telangiectasia, mutated, gene that is more common in patients with breast carcinoma compared with population controls. Cancer, 2004. 100(7):1345–51 DOI: https://doi.org/10.1002/cncr.20133.
]. The intermediate nature of this variant is further evident in the difference in OR for breast cancer between heterozygotic carriers and homozygotic carriers [
Mangone, F.R., E.C. Miracca, H.E. Feilotter, L.M. Mulligan, and M.A. Nagai, ATM gene mutations in sporadic breast cancer patients from Brazil. Springerplus, 2015. 4:23 DOI: 10.1186/s40064-015-0787-z.
], but was benign in our analyses. Previously, the variant D1853V was studied with other variants in an in vitro assessment, and the authors concluded that the variant was fully functional [
]. However, this study only assessed p21 induction as a functional endpoint, finding no difference in this metric. Likewise, we did not see any difference in p21 expression, but our study was much more extensive, and we showed defects across multiple assays for this variant. Furthermore, in contrast to our results [
]. The elevated frequency was surprising considering the rarity of this mutation. This same study found one patient carrying the F1463C mutation as well, but with the P604S variant on the same allele, hampering assessment of the individual variants. F1463C had also previously been noted in a study of patient derived acute lymphoblastic leukemia cell lines, which patients with ataxia telangiectasia are at elevated risk for [
]. However, these studies are statistical studies of small populations, and the predictive ability of any variant may differ widely from population to population and between studies. Additionally, rare variants are harder to study in populations at the statistical level, thus in vitro assessment of the variant is much more likely to yield accurate, biochemically relevant data about any variant. Assessment of variants in ATM in particular is hampered by ATM's nature as a moderate penetrance gene. Compared to genes like BRCA1 or BRCA2, however, ATM mutations are low penetrance in terms of cancer risk [
Germline Pathogenic Variants in the Ataxia Telangiectasia Mutated (<em>ATM</em>) Gene are Associated with High and Moderate Risks for Multiple Cancers.
The frequencies of the alleles that showed functional deficits in our assays are also of note. In gnomAD, a database aggregating the allelic frequencies of healthy individuals, all five of the candidates for pathogenic variants are present at very low frequencies, slightly lower than other databases (Table 1). It is very intriguing that for several of the variants, homozygotes are present in the gnomAD database and none of these individuals are known to have ataxia telangiectasia (AT). This leads to the conclusion that none of the variants for which non-AT homozygotes are documented are complete loss of function alleles. Indeed in our data we observed each of these variants presenting with a significant, but not complete loss of function. The variant S49C, for which our experiments suggest at worst an intermediate phenotype, is present in 12 individuals in the homozygous state. S333F is present in one individual as a homozygote, F1463C in 7 homozygotic individuals, and D1853V is present in 6 homozygotic individuals. Interestingly, the variant for which our study suggested the worst phenotype, V2540I, is present in only 14 individuals in gnomAD, with no homozygotes present.
There are several possible reasons for this discrepancy that are able to incorporate all of the data. The first is that ATM is a moderate penetrance gene, and deleterious mutations in ATM, especially missense mutations that hamper ATM's functions rather than a null allele, will likely retain enough activity to prevent AT. The disease AT itself has a highly variable presentation, with variant AT coined for patients with mild symptoms [
]. Age of symptom onset and presentation are known to be influenced by the genotype of AT patients and whether activity of the ATM protein is retained [
]. Numerous case reports exist for such patients who are either diagnosed late or otherwise do not follow the natural history of AT, including patients with mild or missing symptoms as a result of a mutation with partial function [
ATM Mutations and Phenotypes in Ataxia-Telangiectasia Families in the British Isles: expression of Mutant ATM and the Risk of Leukemia, Lymphoma, and Breast Cancer.
Residual ataxia telangiectasia mutated protein function in cells from ataxia telangiectasia patients, with 5762ins137 and 7271T–>G mutations, showing a less severe phenotype.
]. Frequently, cases of variant AT occur when a patient is the carrier of a hypomorphic ATM allele, which retains some function but otherwise is pathogenic [
]. However, in comparison to AT patients, who display mostly null mutations, missense mutations in ATM are vastly overrepresented in The Cancer Genome Atlas [
It stands to reason, then, that a missense ATM allele may have enough activity to prevent the phenotype of AT while at the same time conferring an elevated lifetime risk of cancer through some impairment of ATM activity. Such an intermediate phenotype is exactly what is suggested by our experiments for the variants S333F, F1463C, D1853V, and possibly S49C, though this latter variant is the most likely to have a benign clinical phenotype, based both on our data, allele frequency, and the data of others reviewed in this discussion. Of our variants, V2540I showed the broadest defects across assays, and possesses the lowest allele frequency with no known homozygotic carriers in gnomAD. Thus, this variant is a clear candidate for pathogenicity.
That these variants have different levels of impairment as seen in the assays conducted likely indicates that the cancer predisposition for patients carrying these mutations will be different. However, how different is difficult to quantify at this stage, and will require further study to explore.
Mechanistically, we explored the effect of each of the six ATM VUS on the phosphorylation of several targets. Our study of kinase targets revealed significant impairment of several phosphorylation events for both the ATM knockdown and kinase dead (S1981A) controls, and similar defects for one of our most likely pathogenic variants from other assays, V2540I. While we did not see any defects in the phosphorylation of downstream targets for F1463C, the specific event causing impairment for this variant may be present in an unstudied phosphorylation target, in the activation of the ATM protein itself, or some other mechanism. Defective kinase activity was also revealed in several assays for S333F, possibly explaining the mechanism for the reduction in HDR activity for this variant.
Other variants showed non-significant changes in phosphorylation activity that was nonetheless visible, such as the reduction in CHK2-T68, p53-S15, and H2AX-139 phosphorylation seen for S49C. It is entirely possible that western blotting is not sensitive enough to discriminate between the small differences in kinase activity seen across assays for S49C. Alternatively, these and possibly additional small but not statistically significant reductions in signaling may be responsible for the reduction in DNA repair ability observed for S49C in the HDR assay. The significant reduction was seen in the phosphorylation of CDC25A, an end stage event in the DNA damage response pathway, may be reflective of this occurrence.
Of the kinase targets studied, only p-p53-S20 was completely without utility in visualizing differences among the VUS. This is likely due to its downstream nature, being phosphorylated by CHK2 after CHK2 is activated by ATM.
Overall, these data lead us to hypothesize that the mechanism by which most, if not all, of the variants which showed reduced function in the HDR and colony survival assays is impaired phosphorylation of one or more downstream targets of ATM in the DNA damage response pathway. Importantly, the specific deficiency in kinase activity appears to vary among the germline ATM variants observed in our cohort.
Exploring the exact mechanism of function loss for each variant would be an extensive process, as ATM has many targets for its kinase activity, and these in turn have even more downstream targets that may be the exact mechanism of the defect. Additionally, some of the loss of function could stem from a small, but target wide defect in ATM's kinase activity, affecting each target only slightly but with an additive affect over the many targets. Finally, the defect may not be related to the kinase function of ATM at all, but could affect folding, multimerization, disassociation from the inactive multimer, or lead to early degradation of the misfolded protein. Such extensive testing for defects is beyond the scope of this work to evaluate the clinical significance of germline missense mutations. Our data here provide strong evidence according to the ACMG guidelines, PS3 (“Well-established in vitro or in vivo functional studies supportive of a damaging effect on the gene or gene product”), and PM2 (“Absent from controls [or at extremely low frequency if recessive]…in Exome Sequencing Project, 1000 Genomes Project, or Exome Aggregation Consortium”) for pathogenicity for V2540I (Table 1) [
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.
]. Replication of this data in clinical laboratories could provide additional justification for the inclusion of these variants as risk alleles.
Even with more knowledge of damaging mutations, moderate penetrance genes like ATM will continue to pose a challenge to physicians and genetic counselors because of the variability in phenotypic presentation. Prophylactic surgery is currently not recommended for carriers of pathogenic mutations in ATM, further limiting prevention options [
]. The best strategy to address this limitation will continue to be the use of the patient's pedigree, individual preference, and the ability of the counselor to interpret the variant data for the benefit of the patient such that the individual can make an informed decision on the direction of her treatment. This further underscores the importance of tests like those we have presented here. Functional data on the VUS enables clinicians to reach conclusions about their individual patient's risk.
Most importantly, this work, along with our previous studies [
], indicates that a cohort of patients with negative genetic diagnoses but pedigrees suggestive of hereditary cancer syndromes can be effectively utilized to curate high yield variant pools for analysis of potentially pathogenic mutations. Such analyses provide an approach that is capable of finding and interrogating the “missing” hereditary mutations that continue to elude clinicians and yield false negative results on cancer gene testing panels, thereby addressing the problem directly. By using arrays of relevant in vitro functional studies, these studies avoid the high false positive and false negative error rates associated with high throughput variant studies and reduce the risk of random error in the results through replication and exclusion of unnecessary assessments by curating the test pool to likely variants found in the patient population. Further research is needed to increase the scale and efficiency of such work. Based on minor allele frequencies and estimated incidence, conclusive data on the six variants investigated in this study could help inform the decisions of 534 ovarian cancer patients and their offspring and a further 8254 breast cancer patients and their offspring. Therefore, based on the cost-benefit of identifying a single pathogenic mutation, small scale but thorough in vitro analysis of variants of unknown significance is useful in promoting human health and increasing the accuracy and utility of cancer genetic testing.
CRediT authorship contribution statement
Scott L. Baughan: Methodology, Project administration, Resources, Investigation, Validation, Formal analysis, Writing – original draft, Writing – review & editing, Visualization. Fatima Darwiche: Investigation, Validation. Michael A. Tainsky: Conceptualization, Funding acquisition, Supervision, Data curation, Methodology, Project administration, Resources, Writing – review & editing, Visualization.
Acknowledgements
This project was supported by The Barbara and Fred Erb Endowed Chair in Cancer Genetics to MAT and from the Graduate School of Wayne State University, SLB.
Multigene Panel Testing Detects Equal Rates of Pathogenic BRCA1/2 Mutations and has a Higher Diagnostic Yield Compared to Limited BRCA1/2 Analysis Alone in Patients at Risk for Hereditary Breast Cancer.
Ionizing radiation induces ataxia telangiectasia mutated-dependent checkpoint signaling and G(2) but not G(1) cell cycle arrest in pluripotent human embryonic stem cells.
Repression of endogenous p53 transactivation function in HeLa cervical carcinoma cells by human papillomavirus type 16 E6, human mdm-2, and mutant p53.
Buchholz, T.A., M.M. Weil, C.L. Ashorn, E.A. Strom, A. Sigurdson, M. Bondy, et al., A Ser49Cys variant in the ataxia telangiectasia, mutated, gene that is more common in patients with breast carcinoma compared with population controls. Cancer, 2004. 100(7):1345–51 DOI: https://doi.org/10.1002/cncr.20133.
Mangone, F.R., E.C. Miracca, H.E. Feilotter, L.M. Mulligan, and M.A. Nagai, ATM gene mutations in sporadic breast cancer patients from Brazil. Springerplus, 2015. 4:23 DOI: 10.1186/s40064-015-0787-z.
Germline Pathogenic Variants in the Ataxia Telangiectasia Mutated (<em>ATM</em>) Gene are Associated with High and Moderate Risks for Multiple Cancers.
ATM Mutations and Phenotypes in Ataxia-Telangiectasia Families in the British Isles: expression of Mutant ATM and the Risk of Leukemia, Lymphoma, and Breast Cancer.
Residual ataxia telangiectasia mutated protein function in cells from ataxia telangiectasia patients, with 5762ins137 and 7271T–>G mutations, showing a less severe phenotype.
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.