Discovering Genetic Interactions in Large-Scale Association Studies by Stage-wise Likelihood Ratio Tests

Despite the success of genome-wide association studies in medical genetics, the underlying genetics of many complex diseases remains enigmatic. One plausible reason for this could be the failure to account for the presence of genetic interactions in current analyses. Exhaustive investigations of interactions are typically infeasible because the vast number of possible interactions impose hard statistical and computational challenges. There is, therefore, a need for computationally efficient methods that build on models appropriately capturing interaction. We introduce a new methodology where we augment the interaction hypothesis with a set of simpler hypotheses that are tested, in order of their complexity, against a saturated alternative hypothesis representing interaction. This sequential testing provides an efficient way to reduce the number of non-interacting variant pairs before the final interaction test. We devise two different methods, one that relies on a priori estimated numbers of marginally associated variants to correct for multiple tests, and a second that does this adaptively. We show that our methodology in general has an improved statistical power in comparison to seven other methods, and, using the idea of closed testing, that it controls the family-wise error rate. We apply our methodology to genetic data from the PRO-CARDIS coronary artery disease case/control cohort and discover three distinct interactions. While analyses on simulated data suggest that the statistical power may suffice for an exhaustive search of all variant pairs in ideal cases, we explore strategies for a priori selecting subsets of variant pairs to test. Our new methodology facilitates identification of new disease-relevant interactions from existing and future genome-wide association data, which may involve genes with previously unknown association to the disease. Moreover, it enables construction of interaction networks that provide a systems biology view of complex diseases, serving as a basis for more comprehensive understanding of disease pathophysiology and its clinical consequences.</p

A simple method for co-segregation analysis to evaluate the pathogenicity of unclassified variants; BRCA1 and BRCA2 as an example

BACKGROUND: Assessment of the clinical significance of unclassified variants (UVs) identified in BRCA1 and BRCA2 is very important for genetic counselling. The analysis of co-segregation of the variant with the disease in families is a powerful tool for the classification of these variants. Statistical methods have been described in literature but these methods are not always easy to apply in a diagnostic setting. METHODS: We have developed an easy to use method which calculates the likelihood ratio (LR) of an UV being deleterious, with penetrance as a function of age of onset, thereby avoiding the use of liability classes. The application of this algorithm is publicly available http://www.msbi.nl/cosegregation. It can easily be used in a diagnostic setting since it requires only information on gender, genotype, present age and/or age of onset for breast and/or ovarian cancer. RESULTS: We have used the algorithm to calculate the likelihood ratio in favour of causality for 3 UVs in BRCA1 (p.M18T, p.S1655F and p.R1699Q) and 5 in BRCA2 (p.E462G p.Y2660D, p.R2784Q, p.R3052W and p.R3052Q). Likelihood ratios varied from 0.097 (BRCA2, p.E462G) to 230.69 (BRCA2, p.Y2660D). Typing distantly related individuals with extreme phenotypes (i.e. very early onset cancer or old healthy individuals) are most informative and give the strongest likelihood ratios for or against causality. CONCLUSION: Although co-segregation analysis on itself is in most cases insufficient to prove pathogenicity of an UV, this method simplifies the use of co-segregation as one of the key features in a multifactorial approach considerably

Somatic genomic alterations in retinoblastoma beyond RB1 are rare and limited to copy number changes

Retinoblastoma is a rare childhood cancer initiated by RB1 mutation or MYCN amplification, while additional alterations may be required for tumor development. However, the view on single nucleotide variants is very limited. To better understand oncogenesis, we determined the genomic landscape of retinoblastoma. We performed exome sequencing of 71 retinoblastomas and matched blood DNA. Next, we determined the presence of single nucleotide variants, copy number alterations and viruses. Aside from RB1, recurrent gene mutations were very rare. Only a limited fraction of tumors showed BCOR (7/71, 10%) or CREBBP alterations (3/71, 4%). No evidence was found for the presence of viruses. Instead, specific somatic copy number alterations were more common, particularly in patients diagnosed at later age. Recurrent alterations of chromosomal arms often involved less than one copy, also in highly pure tumor samples, suggesting within-tumor heterogeneity. Our results show that retinoblastoma is among the least mutated cancers and signify the extreme sensitivity of the childhood retina for RB1 loss. We hypothesize that retinoblastomas arising later in retinal development benefit more from subclonal secondary alterations and therefore, these alterations are more selected for in these tumors. Targeted therapy based on these subclonal events might be insufficient for complete tumor control

A method to assess the clinical significance of unclassified variants in the BRCA1 and BRCA2 genes based on cancer family history

Introduction Unclassified variants (UVs) in the BRCA1/BRCA2 genes are a frequent problem in counseling breast cancer and/or ovarian cancer families. Information about cancer family history is usually available, but has rarely been used to evaluate UVs. The aim of the present study was to identify which is the best combination of clinical parameters that can predict whether a UV is deleterious, to be used for the classification of UVs. Methods We developed logistic regression models with the best combination of clinical features that distinguished a positive control of BRCA pathogenic variants (115 families) from a negative control population of BRCA variants initially classified as UVs and later considered neutral (38 families). Results The models included a combination of BRCAPRO scores, Myriad scores, number of ovarian cancers in the family, the age at diagnosis, and the number of persons with ovarian tumors and/ or breast tumors. The areas under the receiver operating characteristic curves were respectively 0.935 and 0.836 for the BRCA1 and BRCA2 models. For each model, the minimum receiver operating characteristic distance (respectively 90% and 78% specificity for BRCA1 and BRCA2) was chosen as the cutoff value to predict which UVs are deleterious from a study population of 12 UVs, present in 59 Dutch families. The p. S1655F, p. R1699W, and p. R1699Q variants in BRCA1 and the p. Y2660D, p. R2784Q, and p. R3052W variants in BRCA2 are classified as deleterious according to our models. The predictions of the p. L246V variant in BRCA1 and of the p. Y42C, p. E462G, p. R2888C, and p. R3052Q variants in BRCA2 are in agreement with published information of them being neutral. The p. R2784W variant in BRCA2 remains uncertain. Conclusions The present study shows that these developed models are useful to classify UVs in clinical genetic practic

Identification of a BRCA2-Specific modifier locus at 6p24 related to breast cancer risk

Common genetic variants contribute to the observed variation in breast cancer risk for BRCA2 mutation carriers; those known to date have all been found through population-based genome-wide association studies (GWAS). To comprehensively identify breast cancer risk modifying loci for BRCA2 mutation carriers, we conducted a deep replication of an ongoing GWAS discovery study. Using the ranked P-values of the breast cancer associations with the imputed genotype of 1.4 M SNPs, 19,029 SNPs were selected and designed for inclusion on a custom Illumina array that included a total of 211,155 SNPs as part of a multi-consortial project. DNA samples from 3,881 breast cancer affected and 4,330 unaffected BRCA2 mutation carriers from 47 studies belonging to the Consortium of Investigators of Modifiers of BRCA1/2 were genotyped and available for analysis. We replicated previously reported breast cancer susceptibility alleles in these BRCA2 mutation carriers and for several regions (including FGFR2, MAP3K1, CDKN2A/B, and PTHLH) identified SNPs that have stronger evidence of association than those previously published. We also identified a novel susceptibility allele at 6p24 that was inversely associated with risk in BRCA2 mutation carriers (rs9348512; per allele HR = 0.85, 95% CI 0.80-0.90, P = 3.9×10−8). This SNP was not associated with breast cancer risk either in the general population or in BRCA1 mutation carriers. The locus lies within a region containing TFAP2A, which encodes a transcriptional activation protein that interacts with several tumor suppressor genes. This report identifies the first breast cancer risk locus specific to a BRCA2 mutation background. This comprehensive update of novel and previously reported breast cancer susceptibility loci contributes to the establishment of a panel of SNPs that modify breast cancer risk in BRCA2 mutation carriers. This panel may have clinical utility for women with BRCA2 mutations weighing options for medical prevention of breast cancer

Assessing associations between the AURKAHMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood appr

Identification of a BRCA2-Specific Modifier Locus at 6p24 Related to Breast Cancer Risk

Peer reviewe