Search

Background: Despite a large body of literature evaluating the association between recreational physical activity and epithelial ovarian cancer (EOC) risk, the extant evidence is inconclusive, and little is known about the independent association between recreational physical inactivity and EOC risk. We conducted a pooled analysis of nine studies from the Ovarian Cancer Association Consortium to investigate the association between chronic recreational physical inactivity and EOC risk., Methods: In accordance with the 2008 Physical Activity Guidelines for Americans, women reporting no regular, weekly recreational physical activity were classified as inactive. Multivariable logistic regression was utilized to estimate the ORs and 95% confidence intervals (CI) for the association between inactivity and EOC risk overall and by subgroups based upon histotype, menopausal status, race, and body mass index., Results: The current analysis included data from 8,309 EOC patients and 12,612 controls. We observed a significant positive association between inactivity and EOC risk (OR = 1.34; 95% CI, 1.14-1.57), and similar associations were observed for each histotype., Conclusions: In this large pooled analysis examining the association between recreational physical inactivity and EOC risk, we observed consistent evidence of an association between chronic inactivity and all EOC histotypes., Impact: These data add to the growing body of evidence suggesting that inactivity is an independent risk factor for cancer. If the apparent association between inactivity and EOC risk is substantiated, additional work via targeted interventions should be pursued to characterize the dose of activity required to mitigate the risk of this highly fatal disease. Cancer Epidemiol Biomarkers Prev; 25(7); 1114-24. ©2016 AACR., (©2016 American Association for Cancer Research.)

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.

Background: Little is known about modifiable behaviours that may be associated with epithelial ovarian cancer (EOC) survival. We conducted a pooled analysis of 12 studies from the Ovarian Cancer Association Consortium to investigate the association between pre-diagnostic physical inactivity and mortality., Methods: Participants included 6806 women with a primary diagnosis of invasive EOC. In accordance with the Physical Activity Guidelines for Americans, women reporting no regular, weekly recreational physical activity were classified as inactive. We utilised Cox proportional hazard models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) representing the associations of inactivity with mortality censored at 5 years., Results: In multivariate analysis, inactive women had significantly higher mortality risks, with (HR=1.34, 95% CI: 1.18-1.52) and without (HR=1.22, 95% CI: 1.12-1.33) further adjustment for residual disease, respectively., Conclusion: In this large pooled analysis, lack of recreational physical activity was associated with increased mortality among women with invasive EOC.

Data presentation for scientific publications in small sample size studies has not changed substantially in decades. It relies on static figures and tables that may not provide sufficient information for critical evaluation, particularly of the results from small sample size studies. Interactive graphics have the potential to transform scientific publications from static reports of experiments into interactive datasets. We designed an interactive line graph that demonstrates how dynamic alternatives to static graphics for small sample size studies allow for additional exploration of empirical datasets. This simple, free, web-based tool (http://statistika.mfub.bg.ac.rs/interactive-graph/) demonstrates the overall concept and may promote widespread use of interactive graphics.

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer., Competing Interests: The authors declare no competing financial interests.

Carboplatin/taxane combination is first-line therapy for ovarian cancer. However, patients can encounter treatment delays, impaired quality of life, even death because of chemotherapy-induced gastrointestinal (GI) toxicity. A candidate gene study was conducted to assess potential association of genetic variants with GI toxicity in 808 patients who received carboplatin/taxane in the Scottish Randomized Trial in Ovarian Cancer 1 (SCOTROC1). Patients were randomized into discovery and validation cohorts consisting of 404 patients each. Clinical covariates and genetic variants associated with grade III/IV GI toxicity in discovery cohort were evaluated in replication cohort. Chemotherapy-induced GI toxicity was significantly associated with seven single-nucleotide polymorphisms in the ATP7B, GSR, VEGFA and SCN10A genes. Patients with risk genotypes were at 1.53 to 18.01 higher odds to develop carboplatin/taxane-induced GI toxicity (P<0.01). Variants in the VEGF gene were marginally associated with survival time. Our data provide potential targets for modulation/inhibition of GI toxicity in ovarian cancer patients.

Background: Observational studies have reported a positive association between body mass index (BMI) and ovarian cancer risk. However, questions remain as to whether this represents a causal effect, or holds for all histological subtypes. The lack of association observed for serous cancers may, for instance, be due to disease-associated weight loss. Mendelian randomization (MR) uses genetic markers as proxies for risk factors to overcome limitations of observational studies. We used MR to elucidate the relationship between BMI and ovarian cancer, hypothesizing that genetically predicted BMI would be associated with increased risk of non-high grade serous ovarian cancers (non-HGSC) but not HGSC., Methods: We pooled data from 39 studies (14 047 cases, 23 003 controls) in the Ovarian Cancer Association Consortium. We constructed a weighted genetic risk score (GRS, partial F-statistic = 172), summing alleles at 87 single nucleotide polymorphisms previously associated with BMI, weighting by their published strength of association with BMI. Applying two-stage predictor-substitution MR, we used logistic regression to estimate study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted BMI and risk, and pooled these using random-effects meta-analysis., Results: Higher genetically predicted BMI was associated with increased risk of non-HGSC (pooled OR = 1.29, 95% CI 1.03-1.61 per 5 units BMI) but not HGSC (pooled OR = 1.06, 95% CI 0.88-1.27). Secondary analyses stratified by behaviour/subtype suggested that, consistent with observational data, the association was strongest for low-grade/borderline serous cancers (OR = 1.93, 95% CI 1.33-2.81)., Conclusions: Our data suggest that higher BMI increases risk of non-HGSC, but not the more common and aggressive HGSC subtype, confirming the observational evidence., (© The Author 2016; all rights reserved. Published by Oxford University Press on behalf of the International Epidemiological Association.)

Numerous studies demonstrating that statistical errors are common in basic science publications have led to calls to improve statistical training for basic scientists. In this article, we sought to evaluate statistical requirements for PhD training and to identify opportunities for improving biostatistics education in the basic sciences. We provide recommendations for improving statistics training for basic biomedical scientists, including: 1. Encouraging departments to require statistics training, 2. Tailoring coursework to the students' fields of research, and 3. Developing tools and strategies to promote education and dissemination of statistical knowledge. We also provide a list of statistical considerations that should be addressed in statistics education for basic scientists.

Purpose: Socioeconomic status (SES) is a known predictor of survival for several cancers and it has been suggested that SES differences affecting tumour stage at diagnosis may be the most important explanatory factor for this. However, only a limited number of studies have investigated SES differences in tumour stage at diagnosis of ovarian cancer. In a pooled analysis, we investigated whether SES as represented by level of education is predictive for advanced tumour stage at diagnosis of ovarian cancer, overall and by histotype. The effect of cigarette smoking and body mass index (BMI) on the association was also evaluated., Methods: From 18 case-control studies, we obtained information on 10,601 women diagnosed with epithelial ovarian cancer. Study specific odds ratios (ORs) with corresponding 95% confidence intervals (CI) were obtained from logistic regression models and combined into a pooled odds ratio (pOR) using a random effects model., Results: Overall, women who completed ≤high school had an increased risk of advanced tumour stage at diagnosis compared with women who completed >high school (pOR 1.15; 95% CI 1.03-1.28). The risk estimates for the different histotypes of ovarian cancer resembled that observed for ovarian cancers combined but did not reach statistical significance. Our results were unchanged when we included BMI and cigarette smoking., Conclusion: Lower level of education was associated with an increased risk of advanced tumour stage at diagnosis of ovarian cancer. The observed socioeconomic difference in stage at diagnosis of ovarian cancer calls for further studies on how to reduce this diagnostic delay., Competing Interests: None, (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC)., Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped)., Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E-6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E-6; Pcorrected = 0.01)., Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed., Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446-54. ©2016 AACR., (©2016 American Association for Cancer Research.)

Objectives: Bipolar disorder (BD) is a complex disease associated with various hereditary traits, including a higher body mass index (BMI). In a prior genome-wide association study, we found that BMI modified the association of rs12772424 - a common variant in the gene encoding transcription factor 7-like 2 (TCF7L2) - with risk for BD. TCF7L2 is a transcription factor in the canonical Wnt pathway, involved in multiple disorders, including diabetes, cancer and psychiatric conditions. Here, using an independent sample, we evaluated 26 TCF7L2 single nucleotide polymorphisms (SNPs) to explore further the association of BD with the TCF7L2-BMI interaction., Methods: Using a sample of 662 BD cases and 616 controls, we conducted SNP-level and gene-level tests to assess the evidence for an association between BD and the interaction of BMI and genetic variation in TCF7L2. We also explored the potential mechanism behind the detected associations using human brain expression quantitative trait loci (eQTL) analysis., Results: The analysis provided independent evidence of an rs12772424-BMI interaction (p = 0.011). Furthermore, while overall there was no evidence for SNP marginal effects on BD, the TCF7L2-BMI interaction was significant at the gene level (p = 0.042), with seven of the 26 SNPs showing SNP-BMI interaction effects with p < 0.05. The strongest evidence of interaction was observed for rs7895307 (p = 0.006). TCF7L2 expression showed a significant enrichment of association with the expression of other genes in the Wnt canonical pathway., Conclusions: The current study provides further evidence suggesting that TCF7L2 involvement in BD risk may be regulated by BMI. Detailed, prospective assessment of BMI, comorbidity, and other possible contributing factors is necessary to explain fully the mechanisms underlying this association., (© 2016 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Mucocele-like lesions (MLLs) of the breast are characterized by cystic architecture with stromal mucin and frequent atypia, but it is unknown whether they convey long-term breast cancer risk. We evaluated 102 MLLs that were derived from a single-institution benign breast disease cohort of 13412 women who underwent biopsy from 1967 to 2001. MLLs were histologically characterized by type of lining epithelium, architecture of the lesion, associated atypical hyperplasia (AH), and incidence of breast cancer (14.8-year median follow-up). A relatively large proportion of MLLs (42%) were diagnosed in women older than 55 years. AH was significantly more frequent in MLL patient compared to the cohort overall (27% versus 5%; P < .001). Breast cancer has developed in 13 patients with MLL. This frequency is only slightly higher than population expected rates overall (standardized incidence ratio, 2.28; 95% confidence interval, 1.21-3.91) and not significantly different from women in the cohort with (nonatypical) proliferative breast lesions. Younger women (<45) with MLL had a nonsignificant increase in risk of cancer compared to the general population (standardized incidence ratio, 5.16; 95% confidence interval, 1.41-13.23). We conclude that MLL is an uncommon breast lesion that is often associated with coexisting AH. However, in women older than 45 years, MLLs do not convey additional risk of breast cancer beyond that associated with the presence of proliferative disease., Competing Interests: The authors declare that they have no financial or nonfinancial relationships to disclose., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci., Results: At nine of eleven previously published EOC susceptibility regions (2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13), novel variants were identified that were more strongly associated with risk than previously reported variants. Beyond known susceptibility regions, no variants were found to be associated with EOC risk at genome-wide statistical significance (p <5x10(-8)), nor were any significant after Bonferroni correction for 17,000 variants (p< 3x10-6)., Methods: A customized genotyping array was used to assess over 17,000 variants in coding, non-coding, regulatory, and known susceptibility regions in 4,973 EOC cases and 5,640 controls from 13 independent studies. Susceptibility for EOC overall and for select histotypes was evaluated using logistic regression adjusted for age, study site, and population substructure., Conclusion: Given the novel variants identified within the 2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13 regions, larger follow-up genotyping studies, using imputation where necessary, are needed for fine-mapping and confirmation of low frequency variants that fall below statistical significance.

Background: Women with benign breast disease (BBD) have an increased risk of developing breast cancer (BC). Nearly 30% of all BCs develop in women with prior BBD. Information regarding features of the expected number of BCs after BBD would enhance individualized surveillance and prevention strategies for these women. In the current study, the authors sought to characterize BCs developing in a large cohort of women with BBD., Methods: The current study cohort included 13,485 women who underwent breast biopsy for mammographic or palpable concerns between 1967 and 2001. Biopsy slides were reviewed and classified as nonproliferative disease, proliferative disease without atypia, or atypical hyperplasia. BCs were identified by follow-up questionnaires, medical records, and Tumor Registry data. BC tissues were obtained and reviewed., Results: With median follow-up of 15.8 years, 1273 women developed BC. The majority of BCs were invasive (81%), of which 61% were ductal, 13% were mixed ductal/lobular, and 14% were lobular. Approximately two-thirds of the BC cases were intermediate or high grade, and 29% were lymph node positive. Cancer characteristics were similar across the 3 histologic categories of BBD, with a similar frequency of ductal carcinoma in situ, invasive disease, tumor size, time to invasive BC, histologic type of BC, lymph node positivity, and human epidermal growth factor receptor 2 positivity. Women with atypical hyperplasia were found to have a higher frequency of estrogen receptor-positive BC (91%) compared with women with proliferative disease without atypia (80%) or nonproliferative disease (85%) (P = .02)., Conclusions: A substantial percentage of all BCs develop in women with prior BBD. The majority of BCs after BBD are invasive tumors of ductal type, with a substantial number demonstrating lymph node positivity. Of all the BCs in the current study, 84% were estrogen receptor positive. Prevention therapy should be strongly encouraged in higher-risk women with BBD., (© 2015 American Cancer Society.)

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8×10(-11)). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4×10(-8)). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11-1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03-1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction., (© 2016 The authors.)

Machine learning methods, including Random Forests (RF), are increasingly used for genetic data analysis. However, the standard RF algorithm does not correctly model the effects of X chromosome single nucleotide polymorphisms (SNPs), leading to biased estimates of variable importance. We propose extensions of RF to correctly model X SNPs, including a stratified approach and an approach based on the process of X chromosome inactivation. We applied the new and standard RF approaches to case-control alcohol dependence data from the Study of Addiction: Genes and Environment (SAGE), and compared the performance of the alternative approaches via a simulation study. Standard RF applied to a case-control study of alcohol dependence yielded inflated variable importance estimates for X SNPs, even when sex was included as a variable, but the results of the new RF methods were consistent with univariate regression-based approaches that correctly model X chromosome data. Simulations showed that the new RF methods eliminate the bias in standard RF variable importance for X SNPs when sex is associated with the trait, and are able to detect causal autosomal and X SNPs. Even in the absence of sex effects, the new extensions perform similarly to standard RF. Thus, we provide a powerful multimarker approach for genetic analysis that accommodates X chromosome data in an unbiased way. This method is implemented in the freely available R package "snpRF" (http://www.cran.r-project.org/web/packages/snpRF/)., (© 2015 WILEY PERIODICALS, INC.)

Objective: To determine prevalence rates and clinical correlates of current DSM-5 eating disorders in patients with bipolar disorder (BP)., Methods: Prevalence rates of current DSM-5- and DSM-IV-defined binge eating disorder (BED), bulimia nervosa (BN), and anorexia nervosa (AN) were assessed with the Eating Disorder Diagnostic Scale (EDDS) in 1092 patients with BP. Psychiatric illness burden was evaluated with five proxy measures of BP illness severity. Medical illness burden was evaluated with the Cumulative Index Rating Scale (CIRS)., Results: Twenty-seven percent of patients had a current DSM-5 eating disorder: 12% had BED, 15% had BN, and 0.2% had AN. Rates of DSM-5-defined BED and BN were higher than clinical diagnosis rates and rates of DSM-IV-defined BED and BN. Compared with BP patients without an eating disorder, BP patients with a DSM-5 eating disorder were younger and more likely to be women; had an earlier age of onset of BP; had higher EDDS composite scores and higher degrees of suicidality, mood instability, and anxiety disorder comorbidity; and had a higher mean BMI, higher rate of obesity, and higher CIRS total scores. In a logistic regression model controlling for previously identified correlates of an eating disorder, younger age, female gender, and higher BMI remained significantly associated with an eating disorder., Limitations: The EDDS has not been validated in BP patients., Conclusion: DSM-5-defined BED and BN are common in BP patients, possibly more common than DSM-IV-defined BED and BN, and associated with greater psychiatric and general medical illness burden. Further studies assessing DSM-5 eating disorders in people with BP are greatly needed., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Purpose: Chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that germline polymorphisms might be associated with clinical outcome., Experimental Design: We analyzed approximately 2.8 million genotyped and imputed SNPs from the iCOGS experiment for progression-free survival (PFS) and overall survival (OS) in 2,901 European epithelial ovarian cancer (EOC) patients who underwent first-line treatment of cytoreductive surgery and chemotherapy regardless of regimen, and in a subset of 1,098 patients treated with ≥ 4 cycles of paclitaxel and carboplatin at standard doses. We evaluated the top SNPs in 4,434 EOC patients, including patients from The Cancer Genome Atlas. In addition, we conducted pathway analysis of all intragenic SNPs and tested their association with PFS and OS using gene set enrichment analysis., Results: Five SNPs were significantly associated (P ≤ 1.0 × 10(-5)) with poorer outcomes in at least one of the four analyses, three of which, rs4910232 (11p15.3), rs2549714 (16q23), and rs6674079 (1q22), were located in long noncoding RNAs (lncRNAs) RP11-179A10.1, RP11-314O13.1, and RP11-284F21.8, respectively (P ≤ 7.1 × 10(-6)). ENCODE ChIP-seq data at 1q22 for normal ovary show evidence of histone modification around RP11-284F21.8, and rs6674079 is perfectly correlated with another SNP within the super-enhancer MEF2D, expression levels of which were reportedly associated with prognosis in another solid tumor. YAP1- and WWTR1 (TAZ)-stimulated gene expression and high-density lipoprotein (HDL)-mediated lipid transport pathways were associated with PFS and OS, respectively, in the cohort who had standard chemotherapy (pGSEA ≤ 6 × 10(-3))., Conclusions: We have identified SNPs in three lncRNAs that might be important targets for novel EOC therapies., (©2015 American Association for Cancer Research.)

High-risk mutations in several genes predispose to both colorectal cancer (CRC) and endometrial cancer (EC). We therefore hypothesised that some lower-risk genetic variants might also predispose to both CRC and EC. Using CRC and EC genome-wide association series, totalling 13,265 cancer cases and 40,245 controls, we found that the protective allele [G] at one previously-identified CRC polymorphism, rs2736100 near TERT, was associated with EC risk (odds ratio (OR) = 1.08, P = 0.000167); this polymorphism influences the risk of several other cancers. A further CRC polymorphism near TERC also showed evidence of association with EC (OR = 0.92; P = 0.03). Overall, however, there was no good evidence that the set of CRC polymorphisms was associated with EC risk, and neither of two previously-reported EC polymorphisms was associated with CRC risk. A combined analysis revealed one genome-wide significant polymorphism, rs3184504, on chromosome 12q24 (OR = 1.10, P = 7.23 × 10(-9)) with shared effects on CRC and EC risk. This polymorphism, a missense variant in the gene SH2B3, is also associated with haematological and autoimmune disorders, suggesting that it influences cancer risk through the immune response. Another polymorphism, rs12970291 near gene TSHZ1, was associated with both CRC and EC (OR = 1.26, P = 4.82 × 10(-8)), with the alleles showing opposite effects on the risks of the two cancers.

We previously demonstrated that prodynorphin (PDYN) haplotypes and single nucleotide polymorphism (SNP) rs2281285 are associated with alcohol dependence and the propensity to drink in negative emotional states, and recent studies suggest that PDYN gene effects on substance dependence risk may be sex-related. We examined sex-dependent associations of PDYN variation with alcohol dependence and related phenotypes, including negative craving, time until relapse after treatment and the length of sobriety episodes before seeking treatment, in discovery and validation cohorts of European ancestry. We found a significant haplotype-by-sex interaction (p  =  0.03), suggesting association with alcohol dependence in males (p = 1E-4) but not females. The rs2281285 G allele increased risk for alcohol dependence in males in the discovery cohort (OR = 1.49, p = 0.002), with a similar trend in the validation cohort (OR = 1.35, p = 0.086). However, rs2281285 showed a trend towards association with increased negative craving in females in both the discovery (beta = 10.16, p = 0.045) and validation samples (OR = 7.11, p = 0.066). In the discovery cohort, rs2281285 was associated with time until relapse after treatment in females (HR = 1.72, p = 0.037); in the validation cohort, it was associated with increased length of sobriety episodes before treatment in males (beta = 13.49, p = 0.001). Our findings suggest that sex-dependent effects of PDYN variants in alcohol dependence are phenotype-specific.

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Excessive exposure to estrogen is a well-established risk factor for endometrial cancer (EC), particularly for cancers of endometrioid histology. The physiological function of estrogen is primarily mediated by estrogen receptor alpha, encoded by ESR1. Consequently, several studies have investigated whether variation at the ESR1 locus is associated with risk of EC, with conflicting results. We performed comprehensive fine-mapping analyses of 3633 genotyped and imputed single nucleotide polymorphisms (SNPs) in 6607 EC cases and 37 925 controls. There was evidence of an EC risk signal located at a potential alternative promoter of the ESR1 gene (lead SNP rs79575945, P=1.86×10(-5)), which was stronger for cancers of endometrioid subtype (P=3.76×10(-6)). Bioinformatic analysis suggests that this risk signal is in a functionally important region targeting ESR1, and eQTL analysis found that rs79575945 was associated with expression of SYNE1, a neighbouring gene. In summary, we have identified a single EC risk signal located at ESR1, at study-wide significance. Given SNPs located at this locus have been associated with risk for breast cancer, also a hormonally driven cancer, this study adds weight to the rationale for performing informed candidate fine-scale genetic studies across cancer types., (© 2015 Society for Endocrinology.)

Objectives: To examine the independent effects of sex on the risk of rapid cycling and other indicators of adverse illness course in patients with bipolar I disorder (BP-I) or bipolar II disorder (BP-II)., Methods: We analyzed data from the first 1,225 patients enrolled in the Mayo Clinic Individualized Medicine Biobank for Bipolar Disorder. Demographic and clinical variables were ascertained using standardized questionnaires; height and weight were assessed to determine body mass index (BMI). Rates of rapid cycling, cycle acceleration, and increased severity of mood episodes over time were compared between women and men overall and within subgroups defined by bipolar disorder subtype (BP-I or BP-II). Multiple logistic regression analysis was used to assess the independent effect of sex on the risk of these indicators of adverse illness course., Results: Women had significantly higher rates of rapid cycling than men. Overall rates of rapid cycling were higher in patients with BP-II than BP-I; and sex differences in the rate of rapid cycling were more pronounced in patients with BP-II than BP-I, although the power to detect statistically significant differences was reduced due to the lower sample size of subjects with BP-II. Female sex was a significant predictor of rapid cycling, cycle acceleration, and increased severity of mood episodes over time after adjusting for age, bipolar disorder subtype, BMI, having any comorbid psychiatric disorder, and current antidepressant use., Conclusions: Female sex was associated with significantly higher risk of rapid cycling, cycle acceleration, and increased severity of mood episodes over time in a sample of 1,225 patients with bipolar disorders., (© 2015 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Benign breast disease (BBD) is diagnosed in 1-2 million women/year in the US, and while these patients are known to be at substantially increased risk for subsequent development of breast cancer, existing models for risk assessment perform poorly at the individual level. Here, we describe a DNA-microarray-based transcriptional model for breast cancer risk prediction for patients with sclerosing adenosis (SA), which represent ¼ of all BBD patients. A training set was developed from 86 patients diagnosed with SA, of which 27 subsequently developed cancer within 10 years (cases) and 59 remained cancer-free at 10 years (controls). An diagonal linear discriminate analysis-prediction model for prediction of cancer within 10 years (SA TTC10) was generated from transcriptional profiles of FFPE biopsy-derived RNA. This model was tested on a separate validation case-control set composed of 65 SA patients. The SA TTC10 gene signature model, composed of 35 gene features, achieved a clear and significant separation between case and control with receiver operating characteristic area under the curve of 0.913 in the training set and 0.836 in the validation set. Our results provide the first demonstration that benign breast tissue contains transcriptional alterations that indicate risk of breast cancer development, demonstrating that essential precursor biomarkers of malignancy are present many years prior to cancer development. Furthermore, the SA TTC10 gene signature model, which can be assessed on FFPE biopsies, constitutes a novel prognostic biomarker for patients with SA.

Sex differences in incidence and prevalence of and morbidity and mortality from cardiovascular disease are well documented. However, many studies examining the genetic basis for cardiovascular disease fail to consider sex as a variable in the study design, in part, because there is an inherent difficulty in studying the contribution of the sex chromosomes in women due to X chromosome inactivation. This paper will provide general background on the X and Y chromosomes (including gene content, the pseudoautosomal regions, and X chromosome inactivation), discuss how sex chromosomes have been ignored in Genome-wide Association Studies (GWAS) of cardiovascular diseases, and discuss genetics influencing development of cardiovascular risk factors and atherosclerosis with particular attention to carotid intima-medial thickness, and coronary arterial calcification based on sex-specific studies. In addition, a brief discussion of how ethnicity and hormonal status act as confounding variables in sex-based analysis will be considered along with methods for statistical analysis to account for sex in cardiovascular disease., (Copyright © 2015 The Authors. Published by Elsevier Ireland Ltd.. All rights reserved.)

Figures in scientific publications are critically important because they often show the data supporting key findings. Our systematic review of research articles published in top physiology journals (n = 703) suggests that, as scientists, we urgently need to change our practices for presenting continuous data in small sample size studies. Papers rarely included scatterplots, box plots, and histograms that allow readers to critically evaluate continuous data. Most papers presented continuous data in bar and line graphs. This is problematic, as many different data distributions can lead to the same bar or line graph. The full data may suggest different conclusions from the summary statistics. We recommend training investigators in data presentation, encouraging a more complete presentation of data, and changing journal editorial policies. Investigators can quickly make univariate scatterplots for small sample size studies using our Excel templates.

Bipolar disorder (BD) is associated with higher body mass index (BMI) and increased metabolic comorbidity. Considering the associated phenotypic traits in genetic studies of complex diseases, either by adjusting for covariates or by investigating interactions between genetic variants and covariates, may help to uncover the missing heritability. However, obesity-related traits have not been incorporated in prior genome-wide analyses of BD as covariates or potential interacting factors. To investigate the genetic factors underlying BD while considering BMI, we conducted genome-wide analyses using data from the Genetic Association Information Network BD study. We analyzed 729,454 genotyped single-nucleotide polymorphism (SNP) markers on 388 European-American BD cases and 1020 healthy controls with available data for maximum BMI. We performed genome-wide association analyses of the genetic effects while accounting for the effect of maximum BMI, and also evaluated SNP-BMI interactions. A joint test of main and interaction effects demonstrated significant evidence of association at the genome-wide level with rs12772424 in an intron of TCF7L2 (P=2.85E-8). This SNP exhibited interaction effects, indicating that the bipolar susceptibility risk of this SNP is dependent on BMI. TCF7L2 codes for the transcription factor TCF/LF, part of the Wnt canonical pathway, and is one of the strongest genetic risk variants for type 2 diabetes (T2D). This is consistent with BD pathophysiology, as the Wnt pathway has crucial implications in neurodevelopment, neurogenesis and neuroplasticity, and is involved in the mechanisms of action of BD and depression treatments. We hypothesize that genetic risk for BD is BMI dependent, possibly related to common genetic risk with T2D.

Background: Bipolar disorder (BD) is a highly heritable disease. While genome-wide association (GWA) studies have identified several genetic risk factors for BD, few of these studies have investigated the genetic etiology of specific disease subtypes. In particular, BD is positively associated with eating dysregulation traits such as binge eating behavior (BE), yet the genetic risk factors underlying BD with comorbid BE have not been investigated., Methods: Utilizing data from the Genetic Association Information Network study of BD, which included 729,454 single nucleotide polymorphisms (SNPs) genotyped in 1001 European American bipolar cases and 1034 controls, we performed GWA analyses of bipolar subtypes defined by the presence or absence of BE history, and performed a case-only analysis comparing BD subjects with and without BE history. Association signals were refined using imputation, and network analysis was performed with Ingenuity Pathway Analysis software. Based on these results, candidate SNPs were selected for replication in an independent sample of 855 cases and 857 controls., Results: Top ranking SNPs in the discovery set included rs6006893 in PRR5, rs17045162 in ANK2, rs13233490 near PER4, rs4665788 and rs10198175 downstream of APOB, rs2367911 in CACNA2D1, and rs7249968 near ZNF536. Rs10198175 in APOB also demonstrated evidence of association in the replication sample and a meta-analysis of the two samples., Limitations: Without information of BE history in controls, it is not possible to determine whether the observed association with APOB reflects a risk factor for BE behavior in general or a risk factor for a subtype of BD with BE. Further longitudinal and functional studies are needed to determine the causal pathways underlying the observed associations., Conclusions: This study identified new potential BD-susceptibility genes, highlighting the advantages of phenotypic sub-classification in genetic research and clinical practice., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Due to its potential as a biomarker for early cancer detection, blood-based DNA methylation (DNAm) is of interest in cancer research. Specifically, highly predictive mechanisms for early detection of epithelial ovarian cancer (EOC) are desired, so previous studies have compared DNAm between EOC cases and controls. However, case-control studies are confounded by the distribution of white blood cell types through an immune response induced by the cancer. Rather than determining the distribution of the cell types manually or investigating isolated cell types, an alternative approach involves the use of complete blood count (CBC), which is routinely collected. In the analysis of an EOC case-control study of DNAm, we incorporate CBC measures to adjust for this confounding and compare DNAm between 242 EOC cases and 181 age-matched controls (assayed on the Illumina Infinium HumanMethylation27 or HumanMethylation450 Beadchips), at both the individual CpG and CpG island levels. We found that adjustment for leukocyte distribution using CBC measurements dramatically reduced confounding, with 62 single CpG sites found to be associated with EOC status after adjustment (P < 5E-8). Additionally, regional DNAm was assessed by applying principal components analysis to CpG islands. The top associated CpG island (P = 7E-6) was located in the promoter/transcription start site of the human basonuclin 2 gene (BNC2), a known susceptibility gene for EOC risk identified through GWAS. Follow-up studies are necessary to establish the role of BNC2 in blood-based DNA and EOC, including prospective studies to validate this region as a potential biomarker and predictor of EOC susceptibility., (© 2014 WILEY PERIODICALS, INC.)

Background: Genome-wide interrogation of DNA methylation (DNAm) in blood-derived leukocytes has become feasible with the advent of CpG genotyping arrays. In epithelial ovarian cancer (EOC), one report found substantial DNAm differences between cases and controls; however, many of these disease-associated CpGs were attributed to differences in white blood cell type distributions., Methods: We examined blood-based DNAm in 336 EOC cases and 398 controls; we included only high-quality CpG loci that did not show evidence of association with white blood cell type distributions to evaluate association with case status and overall survival., Results: Of 13,816 CpGs, no significant associations were observed with survival, although eight CpGs associated with survival at p < 10-3, including methylation within a CpG island located in the promoter region of GABRE (p = 5.38 x 10-5, HR = 0.95). In contrast, 53 CpG methylation sites were significantly associated with EOC risk (p <5 x10-6). The top association was observed for the methylation probe cg04834572 located approximately 315 kb upstream of DUSP13 (p = 1.6 x10-14). Other disease-associated CpGs included those near or within HHIP (cg14580567; p =5.6x10-11), HDAC3 (cg10414058; p = 6.3x10-12), and SCR (cg05498681; p = 4.8x10-7)., Conclusions: We have identified several CpGs in leukocytes that are differentially methylated by case-control status. Since a retrospective study design was used, we cannot differentiate whether DNAm was etiologic or resulting from EOC; thus, prospective studies of EOC-associated loci are the critical next step.

To further explore reports of association of alcohol dependence and response to acamprosate treatment with the GATA4 rs13273672 single nucleotide polymorphism (SNP), we genotyped this and 10 other GATA4 SNPs in 816 alcohol-dependent cases and 1248 controls. We tested for association of alcohol dependence with the 11 SNPs individually and performed a global test for association using a principle components analysis. Our analyses demonstrate significant association between GATA4 and alcohol dependence at the gene level (P = 0.009) but no association with rs13273672. Further studies are needed to identify potential causal GATA4 variation(s) and the functional mechanism(s) contributing to this association., (© 2012 The Authors, Addiction Biology © 2012 Society for the Study of Addiction.)

Identifying genetic variants associated with complex disease in high-dimensional data is a challenging problem, and complicated etiologies such as gene-gene interactions are often ignored in analyses. The data-mining method Random Forests (RF) can handle high-dimensions; however, in high-dimensional data, RF is not an effective filter for identifying risk factors associated with the disease trait via complex genetic models such as gene-gene interactions without strong marginal components. Here we propose an extension called Weighted Random Forests (wRF), which incorporates tree-level weights to emphasize more accurate trees in prediction and calculation of variable importance. We demonstrate through simulation and application to data from a genetic study of addiction that wRF can outperform RF in high-dimensional data, although the improvements are modest and limited to situations with effect sizes that are larger than what is realistic in genetics of complex disease. Thus, the current implementation of wRF is unlikely to improve detection of relevant predictors in high-dimensional genetic data, but may be applicable in other situations where larger effect sizes are anticipated.

Introduction: One of the proposed psychobiological pathways of craving attributes the desire for drinking in the context of tension, discomfort or unpleasant emotions, to "negative" (or "relief") craving. The aim of this study was to replicate a previously reported association of the PDYN rs2281285 variant with negative craving using a different phenotyping approach., Methods: The TaqMan® Genotyping Assay was used to genotype the rs2281285 variant in 417 German alcohol-dependent subjects. The presence of negative/relief craving was assessed by asking if participants ever ingested alcohol to avoid unwanted emotional or somatic discomfort., Results: The minor allele of rs2281285 was associated with an increased risk of drinking to avoid/escape unwanted emotional or somatic events (OR=2.29, 95% CI=1.08-4.85, p=0.0298)., Discussion: Despite the use of a different phenotyping approach to the measurement of negative craving, our results confirm the association between negative craving and PDYN rs2281285. Genetic markers of negative craving may help to identify subgroups of alcohol-dependent individuals vulnerable to relapse in the context of negative emotions or somatic discomfort, leading to the development of specifically tailored treatment strategies.

Purpose: The combination of a platinum and taxane are standard of care for many cancers, but the utility is often limited due to debilitating neurotoxicity. We examined whether single-nucleotide polymorphisms (SNP) from annotated candidate genes will identify genetic risk for chemotherapy-induced neurotoxicity., Patients and Methods: A candidate-gene association study was conducted to validate the relevance of 1,261 SNPs within 60 candidate genes in 404 ovarian cancer patients receiving platinum/taxane chemotherapy on the SCOTROC1 trial. Statistically significant variants were then assessed for replication in a separate 404 patient replication cohort from SCOTROC1., Results: Significant associations with chemotherapy-induced neurotoxicity were identified and replicated for four SNPs in SOX10, BCL2, OPRM1, and TRPV1. The population attributable risk for each of the four SNPs ranged from 5% to 35%, with a cumulative risk of 62%. According to the multiplicative model, the odds of developing neurotoxicity increase by a factor of 1.64 for every risk genotype. Patients possessing three risk variants have an estimated OR of 4.49 (2.36-8.54) compared to individuals with 0 risk variants. Neither the four SNPs nor the risk score were associated with progression-free survival or overall survival., Conclusions: This study shows that SNPs in four genes have a significant cumulative association with increased risk for the development of chemotherapy-induced neurotoxicity, independent of patient survival., (©2013 AACR.)

Background: Complex psychiatric traits have long been thought to be the result of a combination of genetic and environmental factors, and gene-environment interactions are thought to play a crucial role in behavioral phenotypes and the susceptibility and progression of psychiatric disorders. Candidate gene studies to investigate hypothesized gene-environment interactions are now fairly common in human genetic research, and with the shift toward genome-wide association studies, genome-wide gene-environment interaction studies are beginning to emerge., Methods: We summarize the basic ideas behind gene-environment interaction, and provide an overview of possible study designs and traditional analysis methods in the context of genome-wide analysis. We then discuss novel approaches beyond the traditional strategy of analyzing the interaction between the environmental factor and each polymorphism individually., Results: Two-step filtering approaches that reduce the number of polymorphisms tested for interactions can substantially increase the power of genome-wide gene-environment studies. New analytical methods including data-mining approaches, and gene-level and pathway-level analyses, also have the capacity to improve our understanding of how complex genetic and environmental factors interact to influence psychologic and psychiatric traits. Such methods, however, have not yet been utilized much in behavioral and mental health research., Conclusions: Although methods to investigate gene-environment interactions are available, there is a need for further development and extension of these methods to identify gene-environment interactions in the context of genome-wide association studies. These novel approaches need to be applied in studies of psychology and psychiatry., (© 2013 The Authors. Journal of Child Psychology and Psychiatry © 2013 Association for Child and Adolescent Mental Health.)

Ovarian cancer remains the leading cause of death in women with gynecologic malignancies, despite surgical advances and the development of more effective chemotherapeutics. As increasing evidence indicates that clear-cell ovarian cancer may have unique pathogenesis, further understanding of molecular features may enable us to begin to understand the underlying biology and histology-specific information for improved outcomes. To study epigenetics in clear-cell ovarian cancer, fresh frozen tumor DNA (n = 485) was assayed on Illumina Infinium HumanMethylation450 BeadChips. We identified a clear-cell ovarian cancer tumor methylation profile (n = 163) which we validated in two independent replication sets (set 1, n = 163; set 2, n = 159), highlighting 22 CpG loci associated with nine genes (VWA1, FOXP1, FGFRL1, LINC00340, KCNH2, ANK1, ATXN2, NDRG21 and SLC16A11). Nearly all of the differentially methylated CpGs showed a propensity toward hypermethylation among clear-cell cases. Several loci methylation inversely correlated with tumor gene expression, most notably KCNH2 (HERG, a potassium channel) (P = 9.5 × 10(-7)), indicating epigenetic silencing. In addition, a predicted methylation class mainly represented by the clear-cell cases (20 clear cell out of 23 cases) had improved survival time. Although these analyses included only 30 clear-cell carcinomas, results suggest that loss of expression of KCNH2 (HERG) by methylation could be a good prognostic marker, given that overexpression of the potassium (K(+)) channel Eag family members promotes increased proliferation and results in poor prognosis. Validation in a bigger cohort of clear-cell tumors of the ovary is warranted.

Synthetic κ-opioid receptor (KOR) agonists induce dysphoric and pro-depressive effects and variations in the KOR (OPRK1) and prodynorphin (PDYN) genes have been shown to be associated with alcohol dependence. We genotyped 23 single nucleotide polymorphisms (SNPs) in the PDYN and OPRK1 genes in 816 alcohol-dependent subjects and investigated their association with: (1) negative craving measured by a subscale of the Inventory of Drug Taking Situations; (2) a self-reported history of depression; (3) the intensity of depressive symptoms measured by the Beck Depression Inventory-II. In addition, 13 of the 23 PDYN and OPRK1 SNPs, which were previously genotyped in a set of 1248 controls, were used to evaluate association with alcohol dependence. SNP and haplotype tests of association were performed. Analysis of a haplotype spanning the PDYN gene (rs6045784, rs910080, rs2235751, rs2281285) revealed significant association with alcohol dependence (p = 0.00079) and with negative craving (p = 0.0499). A candidate haplotype containing the PDYN rs2281285-rs1997794 SNPs that was previously associated with alcohol dependence was also associated with negative craving (p = 0.024) and alcohol dependence (p = 0.0008) in this study. A trend for association between depression severity and PDYN variation was detected. No associations of OPRK1 gene variation with alcohol dependence or other studied phenotypes were found. These findings support the hypothesis that sequence variation in the PDYN gene contributes to both alcohol dependence and the induction of negative craving in alcohol-dependent subjects.

Background: Identifying variants associated with complex human traits in high-dimensional data is a central goal of genome-wide association studies. However, complicated etiologies such as gene-gene interactions are ignored by the univariate analysis usually applied in these studies. Random Forests (RF) are a popular data-mining technique that can accommodate a large number of predictor variables and allow for complex models with interactions. RF analysis produces measures of variable importance that can be used to rank the predictor variables. Thus, single nucleotide polymorphism (SNP) analysis using RFs is gaining popularity as a potential filter approach that considers interactions in high-dimensional data. However, the impact of data dimensionality on the power of RF to identify interactions has not been thoroughly explored. We investigate the ability of rankings from variable importance measures to detect gene-gene interaction effects and their potential effectiveness as filters compared to p-values from univariate logistic regression, particularly as the data becomes increasingly high-dimensional., Results: RF effectively identifies interactions in low dimensional data. As the total number of predictor variables increases, probability of detection declines more rapidly for interacting SNPs than for non-interacting SNPs, indicating that in high-dimensional data the RF variable importance measures are capturing marginal effects rather than capturing the effects of interactions., Conclusions: While RF remains a promising data-mining technique that extends univariate methods to condition on multiple variables simultaneously, RF variable importance measures fail to detect interaction effects in high-dimensional data in the absence of a strong marginal component, and therefore may not be useful as a filter technique that allows for interaction effects in genome-wide data.

Paclitaxel is one of the most frequently used chemotherapeutic agents for the treatment of breast cancer patients. Using a candidate gene approach, we hypothesized that polymorphisms in genes relevant to the metabolism and transport of paclitaxel are associated with treatment efficacy and toxicity. Patient and tumor characteristics and treatment outcomes were collected prospectively for breast cancer patients treated with paclitaxel-containing regimens in the neoadjuvant setting. Treatment response was measured before and after each phase of treatment by clinical tumor measurement and categorized according to RECIST criteria, while toxicity data were collected from physician notes. The primary endpoint was achievement of clinical complete response (cCR) and secondary endpoints included clinical response rate (complete response+partial response) and grade 3+ peripheral neuropathy. The genotypes and haplotypes assessed were CYP1B1*3, CYP2C8*3, CYP3A4*1B/CYP3A5*3C, and ABCB1*2. A total of 111 patients were included in this study. Overall, cCR was 30.1% to the paclitaxel component. CYP2C8*3 carriers (23/111, 20.7%) had higher rates of cCR (55% vs. 23%; OR=3.92 [95% CI: 1.46-10.48], corrected p=0.046). In the secondary toxicity analysis, we observed a trend toward greater risk of severe neuropathy (22% vs. 8%; OR=3.13 [95% CI: 0.89-11.01], uncorrected p=0.075) in subjects carrying the CYP2C8*3 variant. Other polymorphisms interrogated were not significantly associated with response or toxicity. Patients carrying CYP2C8*3 are more likely to achieve clinical complete response from neoadjuvant paclitaxel treatment, but may also be at increased risk of experiencing severe peripheral neurotoxicity.

Background: A breadth of high-dimensional data is now available with unprecedented numbers of genetic markers and data-mining approaches to variable selection are increasingly being utilized to uncover associations, including potential gene-gene and gene-environment interactions. One of the most commonly used data-mining methods for case-control data is Multifactor Dimensionality Reduction (MDR), which has displayed success in both simulations and real data applications. Additional software applications in alternative programming languages can improve the availability and usefulness of the method for a broader range of users., Results: We introduce a package for the R statistical language to implement the Multifactor Dimensionality Reduction (MDR) method for nonparametric variable selection of interactions. This package is designed to provide an alternative implementation for R users, with great flexibility and utility for both data analysis and research. The 'MDR' package is freely available online at http://www.r-project.org/. We also provide data examples to illustrate the use and functionality of the package., Conclusions: MDR is a frequently-used data-mining method to identify potential gene-gene interactions, and alternative implementations will further increase this usage. We introduce a flexible software package for R users.

The standard in genetic association studies of complex diseases is replication and validation of positive results, with an emphasis on assessing the predictive value of associations. In response to this need, a number of analytical approaches have been developed to identify predictive models that account for complex genetic etiologies. Multifactor Dimensionality Reduction (MDR) is a commonly used, highly successful method designed to evaluate potential gene-gene interactions. MDR relies on classification error in a cross-validation framework to rank and evaluate potentially predictive models. Previous work has demonstrated the high power of MDR, but has not considered the accuracy and variance of the MDR prediction error estimate. Currently, we evaluate the bias and variance of the MDR error estimate as both a retrospective and prospective estimator and show that MDR can both underestimate and overestimate error. We argue that a prospective error estimate is necessary if MDR models are used for prediction, and propose a bootstrap resampling estimate, integrating population prevalence, to accurately estimate prospective error. We demonstrate that this bootstrap estimate is preferable for prediction to the error estimate currently produced by MDR. While demonstrated with MDR, the proposed estimation is applicable to all data-mining methods that use similar estimates., (© 2010 The Authors Annals of Human Genetics © 2010 Blackwell Publishing Ltd/University College London.)

Background: A fundamental goal of human genetics is the discovery of polymorphisms that predict common, complex diseases. It is hypothesized that complex diseases are due to a myriad of factors including environmental exposures and complex genetic risk models, including gene-gene interactions. Such epistatic models present an important analytical challenge, requiring that methods perform not only statistical modeling, but also variable selection to generate testable genetic model hypotheses. This challenge is amplified by recent advances in genotyping technology, as the number of potential predictor variables is rapidly increasing., Methods: Decision trees are a highly successful, easily interpretable data-mining method that are typically optimized with a hierarchical model building approach, which limits their potential to identify interacting effects. To overcome this limitation, we utilize evolutionary computation, specifically grammatical evolution, to build decision trees to detect and model gene-gene interactions. In the current study, we introduce the Grammatical Evolution Decision Trees (GEDT) method and software and evaluate this approach on simulated data representing gene-gene interaction models of a range of effect sizes. We compare the performance of the method to a traditional decision tree algorithm and a random search approach and demonstrate the improved performance of the method to detect purely epistatic interactions., Results: The results of our simulations demonstrate that GEDT has high power to detect even very moderate genetic risk models. GEDT has high power to detect interactions with and without main effects., Conclusions: GEDT, while still in its initial stages of development, is a promising new approach for identifying gene-gene interactions in genetic association studies.

Background: It is hypothesized that common, complex diseases may be due to complex interactions between genetic and environmental factors, which are difficult to detect in high-dimensional data using traditional statistical approaches. Multifactor Dimensionality Reduction (MDR) is the most commonly used data-mining method to detect epistatic interactions. In all data-mining methods, it is important to consider internal validation procedures to obtain prediction estimates to prevent model over-fitting and reduce potential false positive findings. Currently, MDR utilizes cross-validation for internal validation. In this study, we incorporate the use of a three-way split (3WS) of the data in combination with a post-hoc pruning procedure as an alternative to cross-validation for internal model validation to reduce computation time without impairing performance. We compare the power to detect true disease causing loci using MDR with both 5- and 10-fold cross-validation to MDR with 3WS for a range of single-locus and epistatic disease models. Additionally, we analyze a dataset in HIV immunogenetics to demonstrate the results of the two strategies on real data., Results: MDR with 3WS is computationally approximately five times faster than 5-fold cross-validation. The power to find the exact true disease loci without detecting false positive loci is higher with 5-fold cross-validation than with 3WS before pruning. However, the power to find the true disease causing loci in addition to false positive loci is equivalent to the 3WS. With the incorporation of a pruning procedure after the 3WS, the power of the 3WS approach to detect only the exact disease loci is equivalent to that of MDR with cross-validation. In the real data application, the cross-validation and 3WS analyses indicate the same two-locus model., Conclusions: Our results reveal that the performance of the two internal validation methods is equivalent with the use of pruning procedures. The specific pruning procedure should be chosen understanding the trade-off between identifying all relevant genetic effects but including false positives and missing important genetic factors. This implies 3WS may be a powerful and computationally efficient approach to screen for epistatic effects, and could be used to identify candidate interactions in large-scale genetic studies.