1078 0432 - CCR 17 1397v2
1078 0432 - CCR 17 1397v2
Cancer
Research
Biochemical, Molecular, and Clinical
Characterization of Succinate Dehydrogenase
Subunit A Variants of Unknown Significance
Amber E. Bannon1, Jason Kent1, Isaac Forquer2, Ajia Town3, Lillian R. Klug4,
Kelly McCann5, Carol Beadling6, Oliver Harismendy7, Jason K. Sicklick8,
Christopher Corless9, Ujwal Shinde10, and Michael C. Heinrich11
www.aacrjournals.org OF1
Bannon et al.
process known as flavination. We also measured the protein marked decrease in Sdh1 protein abundance (R179W, G251R,
abundance of Sdh1 and Sdh2 (Fig. 2). Sdh1 flavination is critical H287Y, R303C, Y399C, G410R, C431F, G432E, R444H, A447E,
for the catalytic activity of Sdh1; without covalent attachment of R458W, and R582G). All of the loss-of-function variants resulted
flavin, the Sdh1 protein is unable to oxidize succinate (14). As a in a decrease in Sdh2 protein abundance, consistent with loss of a
positive control for a variant that inhibited insertion of FAD into functional SDH complex as described previously (25). In a subset
Sdh1, we used the previously described Sdh1 H90S variant (14). of no effect and loss-of-function variants, we confirmed that Sdh1
Variants that complemented sdh1D yeast in our functional assays mRNA expression was not significantly different from WT.
had similar levels of flavinated Sdh1, total Sdh1, and total Sdh2 as Based on all of the above results, we characterized 16 variants
yeast complemented with WT Sdh1 (Fig. 2D). However, variants (73%) as causing loss of function and six variants (27%) as having
that were unable to perform oxidative phosphorylation had no effect on SDH function (Table 2). All the loss-of-function
differential effects on the levels of flavinated Sdh1 and total Sdh1 variants were unable to grow on glycerol, had decreased oxygen
protein (Fig. 2A–C). Some loss-of-function variants inhibited consumption, and showed decreased Sdh2 protein abundance. In
flavination without affecting the abundance of Sdh1 (H90S, contrast, the no effect group had similar results to cells comple-
G97R, T134M, R444C, and H601Y), whereas others caused a mented with WT Sdh1 protein in these assays.
A GLUCOSE GLYCEROL
B GLUCOSE GLYCEROL
WT WT
-SDH1 -SDH1
+SDH1 +SDH1
R19X N109S
H90S R162H
G97R R186W
R179W W605F
Figure 1.
ySdh1 variants causing loss of function G251R V633I
are unable to use glycerol as a sole H287Y
carbon source or consume oxygen.
R303C
Yeasts were grown to a confluency
of 2 107/mL and then serially diluted Y399C
and plated onto either glucose or
G410R
glycerol media plates (left to right:
highest to lowest). All variants, C431F
regardless of SDH complex function,
G432E
can grow on glucose, but only those
variants that complement sdh1D and R444C
restore SDH complex activity are able
R444H
to grow on glycerol and consume
oxygen. A, Loss-of-function variants A447E
that are unable to grow with glycerol R458W
as their sole carbon source. B, Variants
with no effect are able to grow with R582W
glycerol as their sole carbon source. H601Y
C, Sdh1 variants that consume
significantly less oxygen than cells
complemented with WT Sdh1 (þSdh1
C D
1.5 1.5
control). D, Variants that consume
oxygen at a similar rate to the þSdh1
0.5 0.5
0.0 0.0
T D 1 X S R R Y C C R F E C H E G Y T D 1 S H I F I
W h1 dh 19 90 97 34M79W51 87 03 99 10 31 32 44 44 47 58W 82 01 W h1 dh 09 62 86W 501 05 633
4 1 6
Sd +S R H G T1 R1 G2 H2 R3 Y3 G4 C4 G4 R4 R4 A R4 R5 H6 Sd +S N R1 R1 T W V
Computational modeling function for most of these amino acids (31). In addition to the
We next performed computational modeling to predict the amino acids directly interacting with the flavin-binding pocket, it
potential structural implications of each variant. Using the two has been previously shown that the C-terminal domain of Sdh1 is
groups (loss-of-function and no effect) identified in our yeast crucial for the flavination of Sdh1 (25). Two variants located in the
model, we further characterized each variant by visualizing the C-terminal domain of Sdh1, yR582G and yH601Y, both inhibited
location in the Sdh1 protein. Both groups, loss-of-function and no flavination of Sdh1. To visualize the substrate (succinate) and
effect, were located throughout the four domains of Sdh1, sug- cofactor (FAD), which were not visualized using the yeast struc-
gesting that there are not "hotspot" areas or domains for path- ture, these Sdh1 variants were modeled using the E. Coli crystal
ogenic variants, unlike the situation for some proteins. Within the structure of Sdh1.
loss-of-function group, 12 variants [yH90S (loss-of-function con- The other six loss-of-function variants [yR19X (control—not
trol), yG97R, yG251R, yH287Y, yR303C, yY399C, yG432E, pictured in computational models because it is an early truncation
yR444C, yR444H, yA447E, yR582G, yH601Y] were identified as of the protein), yT134M, yR179W, yG410R, yG431F, and
being involved in cofactor (FAD) or substrate (succinate) binding yR458W] are not located in the active site of Sdh1. Changes in
in the Sdh1-active site (Fig. 3, Table 2, and Supplementary Video protein structure based on each of these variants are visualized
S1). Structural studies from the homologous protein quinol: in Fig. 4A and Supplementary Video S1. Further, FoldX analyses of
fumarate reductase in E. coli suggests an important catalytic the variants compared with WT provided insight into how the
R179W
G410R
R303C
G432E
G251R
Sdh1D
H287Y
Y399C
+Sdh1
T134M
C431F
Sdh1D
+Sdh1
B
G97R
R19X
H90S
A
WT
WT
FAD-Sdh1 (72 kDa) FAD-Sdh1 (72 kDa)
C D
R458W
R582G
R444C
R444H
R186W
Sdh1D
A447E
H601Y
+Sdh1
W605F
R162H
Sdh1D
N109S
+Sdh1
V633I
T501I
WT
WT
E 1.4
F 1.2
1.2
Normalized Sdh1/Sdh2
1
1
Normalized Sdh2
0.8
0.8
0.6
0.6
0.4 0.4
0.2 0.2
0 0
-Sdh1D
WT
Sdh1
R19X
H90S
G97R
T134M
R179W
G251R
H287Y
R303C
Y399C
G410R
C431F
G432E
R444C
R444H
A447E
R458W
R582G
H601Y
N109S
R162H
R186W
T501I
W605F
V633I
-Sdh1D
WT
Sdh1
N109S
R162H
R186W
T501I
W605F
V633I
Figure 2.
Abundance of flavinated Sdh1, total Sdh1, and total Sdh2 in mitochondrial lysates from yeast expressing each variant. Western blotting of mitochondrial
lysates from each variant is shown. Positive and negative control lanes are included on each Western blot (WT, SDH1delta, þSDH1). WT and þSdh1 show
bands at expected size for each protein. Sdh1D does not show bands for Sdh1-FAD, Sdh1, or Sdh2, but does have normal expression of the mitochondrial
marker, porin. A–C, Loss-of-function variants are shown in numerical order by altered amino acid residue. D, No effect variants are shown in numerical
order by altered amino acid residue and have with no reduction in flavin, Sdh1, or Sdh2 expression compared with WT controls. E, Sdh2 protein abundance
was calculated relative to Sdh1 and then normalized to porin (control for mitochondrial protein loading). The no-effect variants are labeled in black, whereas
the loss-of-function variants are in gray. No-effect variants have Sdh2 protein abundance more similar to complemented Sdh1, and loss-of-function
variants have decreased Sdh2 similar to sdh1D. F, Ratio of Sdh1/Sdh2 protein abundance in no effect variants.
variants affect stability (Supplementary Table S2; ref. 28). problem is exacerbated in GIST as loss of SDHA protein is the
yT134M, yR179W, yG410R, and yG431F all have positive DDG most common cause of SDH deficiency in GIST, but SDHA
(change in free energy), indicating these changes are highly variants have not been extensively characterized at the func-
destabilizing. Some of these also were associated with decreased tional level. Finally, many SDHA VUS are also found in the
Sdh1 protein abundance (Fig. 2). germline, meaning that they can be inherited. Knowing a
Interestingly, all the variants with no functional effects were patient has a loss-of-function SDHA variant in a tumor dra-
localized to the surface of the protein (Fig. 4B; Supplementary matically changes clinical decision making concerning screen-
Video S1) where it would be predicted that they would be less ing recommendations for tumor syndromes, which affects both
likely to affect protein structure and function. Notably, these the patient and other family members. Because there are no
mutations did not map into any of the known critical domains effective medical treatments for advanced SDH-deficient
of Sdh1 or predicted Sdh2 interaction sites. tumors, early detection of tumors allowing curative surgical
resection is crucial.
A major problem in the field has been the lack of a validated
Discussion model system to assess the biochemical function of SDHA VUS.
Molecular testing plays an important role in the clinical Currently, there are no human SDHA-deficient cell lines, and
management of GIST, including decision making about appro- assessing SDH complex activity in tumor samples is difficult.
priate medical and surgical therapy. SDH deficiency is a poten- Previously, it has been shown that a yeast model can be used to
tial cause of GIST in tumors lacking known oncogenic drivers evaluate the pathogenic significance of SDHB mutations
such as gain-of-function KIT mutations. To attempt to identify (17, 18). Extending these studies, we report the successful use
the cause of SDH deficiency in such GIST, it is necessary to use of a yeast model to characterize SDHA variants found in GIST
multi-gene sequencing panels to search for SDHx pathogenic tumors as either causing loss of function or having no bio-
mutations; however, as the routine use of multi-gene sequenc- chemical consequence. Together with computational model-
ing panels has expanded, there has also been an increase in ing, structural homology, and patient data, we have drawn
reported VUS. Due to the multiple potential causes of SDH- conclusions on how and why the SDHA variants we studied
deficiency, it is difficult to determine if a newly discovered VUS affect SDH function, providing clinicians with information to
is responsible for the defect, or if a different genetic mechanism guide genetic counseling of patients and family members who
is responsible for the loss of SDH complex activity. This harbor one of these germline VUS.
Table 2. Consistent findings across multiple functional assays for each variant in our yeast model
Growth Oxygen consump- Sdh2 Protein
hSDHA ySdh1 phenotype on tion (% of (% of Structural implications of
Mutation Mutation glycerol complemented) complemented) variants Group name
R31X R19X Unable to grow Decreased (20%) Decreased (0%) Truncates mitochondrial Loss-of-function
(control) (control) targeting sequence (47) (control)
H99S H90S Unable to grow Decreased (2%) Decreased (15%) Inhibits covalent bond to Loss-of-function
(control) (control) FAD (14) (control)
G106R G97R Unable to grow Decreased (7%) Decreased (14%) Distorts the active site (31) Loss-of-function
N118S N109S Growth Similar (108%) Similar (83%) Surface of protein No effect
T143M T134M Unable to grow Decreased (3%) Decreased (35%) No obvious disturbances Loss-of-function
Using our yeast model, we characterized variants as either the loss-of-function variants were associated with a dramatic
causing loss of protein function or having no effect. All variants decrease in Sdh2 protein abundance, which is in concordance
in the loss-of-function group were unable to grow on a nonfer- with current clinical testing to determine SDH-deficiency in
mentable carbon source (glycerol). In addition, these variants also tumors—assessment of SDHB (human equivalent to ySdh2)
failed to consume oxygen, indicating defective oxidative phos- expression using IHC (2, 32–36). Our data provide further
phorylation. Interestingly, these variants had differential effects evidence supporting the role of SDHB IHC as the most reliable
on Sdh1 flavination and/or protein abundance. As a group, all of clinical test to identify SDH-deficiency in a tumor. Notably, we
identified four novel loss-of-function variants that did not affect The conclusions drawn from our yeast model are largely
Sdh1 protein abundance (yG97R, yT134M, yR444C, yH601Y). supported by our computational modeling results. For many
These variants prevent flavination of Sdh1, making Sdh1 dysfunc- variants, we could find confirmatory structural modeling evi-
tional, but still expressed at normal levels. Our findings are dence that the amino acids affected by these mutations played a
consistent with previous reports where SDHA expression was role in the catalytic site of Sdh1. The computational model and
retained in some SDHA-mutant tumors that lack SDHB expres- energy minimization also allowed us to explain why some
sion (37), confirming that not all dysfunctional SDHA (or Sdh1) variants could reduce Sdh1 protein abundance. Furthermore,
variants lead to decreased protein expression of this subunit. all the variants that did not affect protein function were on the
There are three clinical samples (R188W, G260R, A454E) that surface of the protein as visualized by our computational
we characterized as loss-of-function variants but are associated model, suggesting that they would not interact with the cata-
with normal sdh1 in our yeast model. Based on these results and lytic mechanism of Sdh1.
coupled with the fact the SDHA IHC is not universally available, Based on our biochemical data, there are still several tumors
we advocate for the use of SDHB IHC to identify tumors that are that have no known mechanism for loss of SDH activity
SDH-deficient. (Table 1). One of these tumors had the hR195W (yR186W)
variant with an allele frequency of 66.7% and negative SDHB sequencing panels (15, 16). Given the limitations of using
IHC. This tumor is likely driven by a different mechanism of deidentified results from clinical testing, we did not have access
SDH-deficiency. The rest of the unexplained tumors had no to residual tumor samples to perform additional IHC or geno-
IHC data available, so it is possible that these GISTs are not mic testing.
driven by SDH-deficiency and instead belong to a different We used a yeast model to characterize 22 SDHA VUS. These
subtype of GIST (11). Alternatively, these GISTs may be SDH- data revealed 16 (73%) of SDHA VUS as loss of function
deficient by a different genetic mechanism than the SDHA (and therefore pathogenic), highlighting the importance of
variant identified using targeted exome sequencing. There are understanding such variants to provide better clinical recommen-
several weaknesses of various targeted exome sequencing dations for genetic counselors concerning family screening
panels including potential lack of coverage for regions of and early detection protocols. However, our approach using a
certain SDHx genes, failure to detect large genomic deletions functional yeast model paired with computational modeling
involving SDHx subunits, as well as the inability to measure can distinguish between SDHA VUS that cause loss of
hypermethylation of the SDHC promoter. In addition, muta- function and those that have no biochemical effect, allowing
tions that inactivate newly identified SDH assembly factors us to discriminate between pathogenic and nonpathogenic
(SDHAF3, 4) are typically not included in clinical cancer gene variants.
References
1. Bezawork-Geleta A, Rohlena J, Dong L, Pacak K, Neuzil J. Mitochondrial 14. Hao HX, Khalimonchuk O, Schraders M, Dephoure N, Bayley JP, Kunst H,
complex II: at the crossroads. Trends Biochem Sci 2017;42:312–25. et al. SDH5, a gene required for flavination of succinate dehydrogenase, is
2. Janeway KA, Kim SY, Lodish M, Nose V, Rustin P, Gaal J, et al. Defects mutated in paraganglioma. Science 2009;325:1139–42.
in succinate dehydrogenase in gastrointestinal stromal tumors lacking 15. Van Vranken JG, Bricker DK, Dephoure N, Gygi SP, Cox JE, Thummel CS,
KIT and PDGFRA mutations. Proc Natl Acad Sci U S A 2011;108: et al. SDHAF4 promotes mitochondrial succinate dehydrogenase activity
314–8. and prevents neurodegeneration. Cell Metab 2014;20:241–52.
3. Burnichon N, Abermil N, Buffet A, Favier J, Gimenez-Roqueplo AP. The 16. Na U, Yu W, Cox J, Bricker DK, Brockmann K, Rutter J, et al. The LYR factors
genetics of paragangliomas. Eur Ann Otorhinolaryngol Head Neck Dis SDHAF1 and SDHAF3 mediate maturation of the iron-sulfur subunit of
2012;129:315–8. succinate dehydrogenase. Cell Metab 2014;20:253–66.
4. Baysal BE. A recurrent stop-codon mutation in succinate dehydrogenase 17. Panizza E, Ercolino T, Mori L, Rapizzi E, Castellano M, Opocher G, et al.
subunit B gene in normal peripheral blood and childhood T-cell acute Yeast model for evaluating the pathogenic significance of SDHB, SDHC
leukemia. PLoS One 2007;2:e436. and SDHD mutations in PHEO-PGL syndrome. Hum Mol Genet 2013;
5. Beamer LC. Cowden syndrome: what oncology nurses need to know about 22:804–15.
increased risk of developing certain cancers. Oncol Nurs Forum 18. Goffrini P, Ercolino T, Panizza E, Giache V, Cavone L, Chiarugi A, et al.
2014;41:555–7. Functional study in a yeast model of a novel succinate dehydrogenase
6. Niemeijer ND, Papathomas TG, Korpershoek E, de Krijger RR, Oudijk L, subunit B gene germline missense mutation (C191Y) diagnosed in a
Morreau H, et al. Succinate dehydrogenase (SDH)-deficient pancreatic patient affected by a glomus tumor. Hum Mol Genet 2009;18:1860–8.
neuroendocrine tumor expands the SDH-related tumor spectrum. J Clin 19. Szeto SS, Reinke SN, Sykes BD, Lemire BD. Ubiquinone-binding site
Endocrinol Metab 2015;100:E1386–93. mutations in the Saccharomyces cerevisiae succinate dehydrogenase gen-
7. van Nederveen FH, Gaal J, Favier J, Korpershoek E, Oldenburg RA, de Bruyn erate superoxide and lead to the accumulation of succinate. J Biol Chem
EM, et al. An immunohistochemical procedure to detect patients with 2007;282:27518–26.
paraganglioma and phaeochromocytoma with germline SDHB, SDHC, or 20. Alston CL, Ceccatelli Berti C, Blakely EL, Olahova M, He L, McMahon CJ,
SDHD gene mutations: a retrospective and prospective analysis. Lancet et al. A recessive homozygous p.Asp92Gly SDHD mutation causes prenatal
Oncol 2009;10:764–71. cardiomyopathy and a severe mitochondrial complex II deficiency. Hum
8. Killian JK, Miettinen M, Walker RL, Wang Y, Zhu YJ, Waterfall JJ, et al. Genet 2015;134:869–79.
Recurrent epimutation of SDHC in gastrointestinal stromal tumors. Sci 21. Fontanesi F, Diaz F, Barrientos A. Evaluation of the mitochondrial respi-
Transl Med 2014;6:268ra177. ratory chain and oxidative phosphorylation system using yeast models of
9. Ricketts CJ, Forman JR, Rattenberry E, Bradshaw N, Lalloo F, Izatt L, et al. OXPHOS deficiencies. Curr Protoc Hum Genet 2009;Chapter 19:Unit19.5.
Tumor risks and genotype-phenotype-proteotype analysis in 358 patients 22. Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B,
with germline mutations in SDHB and SDHD. Hum Mutat 2010;31:41–51. et al. Functional characterization of the S. cerevisiae genome by gene
10. Neumann HP, Pawlu C, Peczkowska M, Bausch B, McWhinney SR, Mur- deletion and parallel analysis. Science 1999;285:901–6.
esan M, et al. Distinct clinical features of paraganglioma syndromes 23. Sikorski RS, Hieter P. A system of shuttle vectors and yeast host strains
associated with SDHB and SDHD gene mutations. JAMA 2004;292: designed for efficient manipulation of DNA in Saccharomyces cerevisiae.
943–51. Genetics 1989;122:19–27.
11. Bannon AE, Klug LR, Corless CL, Heinrich MC. Using molecular diagnostic 24. Muller F, Crofts AR, Kramer DM. Multiple Q-cycle bypass reactions at the
testing to personalize the treatment of patients with gastrointestinal Qo site of the cytochrome bc1 complex. Biochemistry 2002;41:7866–74.
stromal tumors. Expert Rev Mol Diagn 2017;17:445–57. 25. Kim HJ, Jeong MY, Na U, Winge DR. Flavinylation and assembly of
12. Evenepoel L, Papathomas TG, Krol N, Korpershoek E, de Krijger RR, Persu succinate dehydrogenase are dependent on the C-terminal tail of the
A, et al. Toward an improved definition of the genetic and tumor spectrum flavoprotein subunit. J Biol Chem 2012;287:40670–9.
associated with SDH germ-line mutations. Genet Med 2015;17:610–20. 26. Bafunno V, Giancaspero TA, Brizio C, Bufano D, Passarella S, Boles E, et al.
13. Shi E, Chmielecki J, Tang CM, Wang K, Heinrich MC, Kang G, et al. FGFR1 Riboflavin uptake and FAD synthesis in Saccharomyces cerevisiae mito-
and NTRK3 actionable alterations in "Wild-Type" gastrointestinal stromal chondria: involvement of the Flx1p carrier in FAD export. J Biol Chem
tumors. J Transl Med 2016;14:339. 2004;279:95–102.
27. Kitanovic A, Bonowski F, Heigwer F, Ruoff P, Kitanovic I, Ungewiss C, et al. 38. Italiano A, Chen CL, Sung YS, Singer S, DeMatteo RP, LaQuaglia MP, et al.
Acetic acid treatment in S. cerevisiae creates significant energy deficiency SDHA loss of function mutations in a subset of young adult wild-type
and nutrient starvation that is dependent on the activity of the mitochon- gastrointestinal stromal tumors. BMC Cancer 2012;12:408.
drial transcriptional complex Hap2–3-4–5. Front Oncol 2012;2:118. 39. Pantaleo MA, Nannini M, Astolfi A, Biasco G. A distinct pediatric-type
28. Schymkowitz J, Borg J, Stricher F, Nys R, Rousseau F, Serrano L. The FoldX gastrointestinal stromal tumor in adults: potential role of succinate
web server: an online force field. Nucleic Acids Res 2005;33(Web Server dehydrogenase subunit A mutations. Am J Surg Pathol 2011;35:
issue):W382–8. 1750–2.
29. Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. 40. Korpershoek E, Favier J, Gaal J, Burnichon N, van Gessel B, Oudijk L, et al.
Analysis of protein-coding genetic variation in 60,706 humans. Nature SDHA immunohistochemistry detects germline SDHA gene mutations in
2016;536:285–91. apparently sporadic paragangliomas and pheochromocytomas. J Clin
30. Landrum MJ, Lee JM, Benson M, Brown G, Chao C, Chitipiralla S, et al. Endocrinol Metab 2011;96:E1472–6.