Efficient computation of the joint probability of multiple inherited risk alleles from pedigree data

Variant-specific Mendelian Risk Prediction Model

Many pathogenic sequence variants (PSVs) have been associated with increased risk of cancers. Mendelian risk prediction models use Mendelian laws of inheritance to predict the probability of having a PSV based on family history, as well as specified PSV frequency and penetrance (agespecific probability of developing cancer given genotype). Most existing models assume penetrance is the same for any PSVs in a certain gene. However, for some genes (for example, BRCA1/2), cancer risk does vary by PSV. We propose an extension of Mendelian risk prediction models to relax the assumption that risk is the same for any PSVs in a certain gene by incorporating variant-specific penetrances and illustrating these extensions on two existing Mendelian risk prediction models, BRCAPRO and PanelPRO. Our proposed BRCAPRO-variant and PanelPRO-variant models incorporate variant-specific BRCA1/2 PSVs through the region classifications. Due to the sparsity of the variant information we classify BRCA1/2 PSVs into three regions; the breast cancer clustering region (BCCR), the ovarian cancer clustering region (OCCR), and an other region. Simulations were conducted to evaluate the performance of the proposed BRCAPRO-variant model compared to the existing BRCAPRO model which assumes the penetrance is the same for any PSVs in BRCA1 (and respectively BRCA2). Simulation results showed that the BRCAPRO-variant model was well calibrated to predict region-specific BRCA1/2 carrier status with high discrimination and accuracy on the region-specific level. In addition, we showed that the BRCAPRO-variant model achieved performance gains over the existing risk prediction models in terms of calibration without loss in discrimination and accuracy. We also evaluated the performance of the two proposed models, BRCAPRO-variant and PanelPRO-variant, on a cohort of 1,961 families from the Cancer Genetics Network (CGN). We showed that our proposed models provide region-specific PSV carrier probabilities with high accuracy, while the calibration, discrimination and accuracy of gene-specific PSV carrier probabilities were comparable to the existing gene-specific models. As more variant-specific PSV penetrances become available, we have shown that Mendelian risk prediction models can be extended to integrate the additional information, providing precise variant or region-specific PSV carrier probabilities and improving future cancer risk predictions.

Statistical methods for Mendelian models with multiple genes and cancers

Risk evaluation to identify individuals who are at greater risk of cancer as a result of heritable pathogenic variants is a valuable component of individualized clinical management. Using principles of Mendelian genetics, Bayesian probability theory, and variant-specific knowledge, Mendelian models derive the probability of carrying a pathogenic variant and developing cancer in the future, based on family history. Existing Mendelian models are widely employed, but are generally limited to specific genes and syndromes. However, the upsurge of multigene panel germline testing has spurred the discovery of many new gene-cancer associations that are not presently accounted for in these models. We have developed PanelPRO, a flexible, efficient Mendelian risk prediction framework that can incorporate an arbitrary number of genes and cancers, overcoming the computational challenges that arise because of the increased model complexity. We implement an 11-gene, 11-cancer model, the largest Mendelian model created thus far, based on this framework. Using simulations and a clinical cohort with germline panel testing data, we evaluate model performance, validate the reverse-compatibility of our approach with existing Mendelian models, and illustrate its usage. Our implementation is freely available for research use in the PanelPRO R package.

Multi-syndrome, multi-gene risk modeling for individuals with a family history of cancer with the novel R package PanelPRO

Identifying individuals who are at high risk of cancer due to inherited germline mutations is critical for effective implementation of personalized prevention strategies. Most existing models focus on a few specific syndromes; however, recent evidence from multi-gene panel testing shows that many syndromes are overlapping, motivating the development of models that incorporate family history on several cancers and predict mutations for a comprehensive panel of genes.

We present PanelPRO, a new, open-source R package providing a fast, flexible back-end for multi-gene, multi-cancer risk modeling with pedigree data. It includes a customizable database with default parameter values estimated from published studies and allows users to select any combinations of genes and cancers for their models, including well-established single syndrome BayesMendel models (BRCAPRO and MMRPRO). This leads to more accurate risk predictions and ultimately has a high impact on prevention strategies for cancer and clinical decision making. The package is available for download for research purposes at https://projects.iq.harvard.edu/bayesmendel/panelpro.