Up until recently, determining penetrance of pathogenic variants required large observational cohort studies. The availability of large population cohort data in the Exome Aggregate Consortium (ExAC) now makes it possible to use a Bayesian approach to calculate penetrance, in that population frequencies of pathogenic germline variants should be inversely proportional to their penetrance for disease. We tested this hypothesis using data from our genetic testing for medullary thyroid cancer (MTC) and pheochromocytoma/paraganglioma (PC/PGL) susceptibility genes. Frequency of each pathogenic variant found in RET, SDH subunit genes A-D and VHL in our cases was compared with its frequency in ExAC.
Genetic testing was performed in 437 Australian subjects presenting with MTC and 575 subjects presenting with PC/PGL. We identified pathogenic RET variants in 58 subjects (13.3%) presenting with MTC, and pathogenic RET, SDHx, or VHL variants in 162 subjects (28.3%) presenting with pheochromocytoma or paraganglioma. Many of these variants are reported in ExAC (for RET 38%, for SDHA 29%, for SDHB 20%, for SDHC 20% and for VHL 7%). None of our pathogenic SDHD variants are reported in ExAC. Cumulative frequency of these individually rare pathogenic variants in ExAC was 0.00015 (RET) and 0.0002 (SDHA, SDHB, SDHC), four-fold higher than expected from observed population prevalence of MTC or PC/PGL. For RET, the three American Thyroid Association (ATA) risk categories correlated well with groups of variants predicted to have high, moderate, or low/negligible risk of disease penetrance. For SDHB variants present in ExAC, our estimates predict lifetime penetrance of PC/PGL of 23% (95% confidence interval [CI] 12-39%) and lower for SDHA (1.5%, CI 0.5-4.2%) or SDHC (8%, CI 1-33%). Our findings have important implications for genetic counseling and surveillance of subjects carrying pathogenic variants in these genes.