BRCA1 is an important mediator of the DNA damage response, which promotes homologous recombination (HR) and antagonizes 53BP1-dependent non-homologous end joining in S/G2 phase. But how this is achieved remains unclear. Here, we report that the E3 ubiquitin ligase UHRF1 (Ubiquitin-like, with PHD and RING finger domains 1) directly participates in the interplay between BRCA1 and 53BP1. Mechanistically, UHRF1 is recruited to DNA double-strand breaks (DSBs) by BRCA1 in S phase, which requires the BRCT domain of BRCA1 and phosphorylated Ser674 of UHRF1. Subsequently, UHRF1 mediates K63-linked polyubiquitination of RIF1, and results in its dissociation from 53BP1 and DSBs thereby facilitating HR initiation. Thus, UHRF1 is a key regulator of DSB repair choice, which is separate from its role in heterochromatin formation and epigenetic regulator.

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.

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

Ovarian cancer represents the most lethal tumor type among malignancies of the female reproductive system. Overall survival rates remain low. In this study, we identify the serine protease inhibitor Kazal type 1 (SPINK1) as a potential therapeutic target for a subset of ovarian cancers. We show that SPINK1 drives ovarian cancer cell proliferation through activation of epidermal growth factor receptor (EGFR) signaling, and that SPINK1 promotes resistance to anoikis through a distinct mechanism involving protease inhibition. In analyses of ovarian tumor specimens from a Mayo Clinic cohort of 490 patients, we further find that SPINK1 immunostaining represents an independent prognostic factor for poor survival, with the strongest association in patients with nonserous histological tumor subtypes (endometrioid, clear cell, and mucinous). This study provides novel insight into the fundamental processes underlying ovarian cancer progression, and also suggests new avenues for development of molecularly targeted therapies.

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

COMMD1 deficiency results in defective copper homeostasis, but the mechanism for this has remained elusive. Here we report that COMMD1 is directly linked to early endosomes through its interaction with a protein complex containing CCDC22, CCDC93, and C16orf62. This COMMD/CCDC22/CCDC93 (CCC) complex interacts with the multisubunit WASH complex, an evolutionarily conserved system, which is required for endosomal deposition of F-actin and cargo trafficking in conjunction with the retromer. Interactions between the WASH complex subunit FAM21, and the carboxyl-terminal ends of CCDC22 and CCDC93 are responsible for CCC complex recruitment to endosomes. We show that depletion of CCC complex components leads to lack of copper-dependent movement of the copper transporter ATP7A from endosomes, resulting in intracellular copper accumulation and modest alterations in copper homeostasis in humans with CCDC22 mutations. This work provides a mechanistic explanation for the role of COMMD1 in copper homeostasis and uncovers additional genes involved in the regulation of copper transporter recycling.

Accurately identifying patients with high-grade serous ovarian carcinoma (HGSOC) who respond to poly(ADP-ribose) polymerase inhibitor (PARPi) therapy is of great clinical importance. Here we show that quantitative BRCA1 methylation analysis provides new insight into PARPi response in preclinical models and ovarian cancer patients. The response of 12 HGSOC patient-derived xenografts (PDX) to the PARPi rucaparib was assessed, with variable dose-dependent responses observed in chemo-naive BRCA1/2-mutated PDX, and no responses in PDX lacking DNA repair pathway defects. Among BRCA1-methylated PDX, silencing of all BRCA1 copies predicts rucaparib response, whilst heterozygous methylation is associated with resistance. Analysis of 21 BRCA1-methylated platinum-sensitive recurrent HGSOC (ARIEL2 Part 1 trial) confirmed that homozygous or hemizygous BRCA1 methylation predicts rucaparib clinical response, and that methylation loss can occur after exposure to chemotherapy. Accordingly, quantitative BRCA1 methylation analysis in a pre-treatment biopsy could allow identification of patients most likely to benefit, and facilitate tailoring of PARPi therapy.

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

Previous studies have shown that human immunodeficiency virus (HIV) protease cleaves procaspase 8 to a fragment, termed Casp8p41, that lacks caspase activity but nonetheless contributes to T cell apoptosis. Herein, we show that Casp8p41 contains a domain that interacts with the BH3-binding groove of pro-apoptotic Bak to cause Bak oligomerization, Bak-mediated membrane permeabilization, and cell death. Levels of active Bak are higher in HIV-infected T cells that express Casp8p41. Conversely, targeted mutations in the Bak-interacting domain diminish Bak binding and Casp8p41-mediated cell death. Similar mutations in procaspase 8 impair the ability of HIV to kill infected T cells. These observations support a novel paradigm in which HIV converts a normal cellular constituent into a direct activator that functions like a BH3-only protein.

Mammalian target of rapamycin (mTOR) represents a key downstream intermediate for a myriad of oncogenic receptor tyrosine kinases. In the case of the insulin-like growth factor (IGF) pathway, the mTOR complex (mTORC1) mediates IGF-1 receptor (IGF-1R)-induced estrogen receptor alpha (ERα) phosphorylation/activation and leads to increased proliferation and growth in breast cancer cells. As a result, the prevalence of mTOR inhibitors combined with hormonal therapy has increased in recent years. Conversely, activated mTORC1 provides negative feedback regulation of IGF signaling via insulin receptor substrate (IRS)-1/2 serine phosphorylation and subsequent proteasomal degradation. Thus, the IGF pathway may provide escape (e.g. de novo or acquired resistance) from mTORC1 inhibitors. It is therefore plausible that combined inhibition of mTORC1 and IGF-1R for select subsets of ER-positive breast cancer patients presents as a viable therapeutic option.

Interest in preclinical drug development for ovarian cancer has stimulated development of patient-derived xenograft (PDX) or tumorgraft models. However, the unintended formation of human lymphoma in severe combined immunodeficiency (SCID) mice from Epstein-Barr virus (EBV)-infected human lymphocytes can be problematic. In this study, we have characterized ovarian cancer PDXs which developed human lymphomas and explore methods to suppress lymphoproliferative growth. Fresh human ovarian tumors from 568 patients were transplanted intraperitoneally in SCID mice. A subset of PDX models demonstrated atypical patterns of dissemination with mediastinal masses, hepatosplenomegaly, and CD45-positive lymphoblastic atypia without ovarian tumor engraftment. Expression of human CD20 but not CD3 supported a B-cell lineage, and EBV genomes were detected in all lymphoproliferative tumors. Immunophenotyping confirmed monoclonal gene rearrangements consistent with B-cell lymphoma, and global gene expression patterns correlated well with other human lymphomas. The ability of rituximab, an anti-CD20 antibody, to suppress human lymphoproliferation from a patient's ovarian tumor in SCID mice and prevent growth of an established lymphoma led to a practice change with a goal to reduce the incidence of lymphomas. A single dose of rituximab during the primary tumor heterotransplantation process reduced the incidence of CD45-positive cells in subsequent PDX lines from 86.3% (n = 117 without rituximab) to 5.6% (n = 160 with rituximab), and the lymphoma rate declined from 11.1% to 1.88%. Taken together, investigators utilizing PDX models for research should routinely monitor for lymphoproliferative tumors and consider implementing methods to suppress their growth.

The mechanism by which the proapoptotic Bcl-2 family members Bax and Bak release cytochrome c from mitochondria is incompletely understood. In this paper, we show that activator BH3-only proteins bind tightly but transiently to the Bak hydrophobic BH3-binding groove to induce Bak oligomerization, liposome permeabilization, mitochondrial cytochrome c release, and cell death. Analysis by surface plasmon resonance indicated that the initial binding of BH3-only proteins to Bak occurred with similar kinetics with or without detergent or mitochondrial lipids, but these reagents increase the strength of the Bak-BH3-only protein interaction. Point mutations in Bak and reciprocal mutations in the BH3-only proteins not only confirmed the identity of the interacting residues at the Bak-BH3-only protein interface but also demonstrated specificity of complex formation in vitro and in a cellular context. These observations indicate that transient protein-protein interactions involving the Bak BH3-binding groove initiate Bak oligomerization and activation.

Gene set methods aim to assess the overall evidence of association of a set of genes with a phenotype, such as disease or a quantitative trait. Multiple approaches for gene set analysis of expression data have been proposed. They can be divided into two types: competitive and self-contained. Benefits of self-contained methods include that they can be used for genome-wide, candidate gene, or pathway studies, and have been reported to be more powerful than competitive methods. We therefore investigated ten self-contained methods that can be used for continuous, discrete and time-to-event phenotypes. To assess the power and type I error rate for the various previously proposed and novel approaches, an extensive simulation study was completed in which the scenarios varied according to: number of genes in a gene set, number of genes associated with the phenotype, effect sizes, correlation between expression of genes within a gene set, and the sample size. In addition to the simulated data, the various methods were applied to a pharmacogenomic study of the drug gemcitabine. Simulation results demonstrated that overall Fisher's method and the global model with random effects have the highest power for a wide range of scenarios, while the analysis based on the first principal component and Kolmogorov-Smirnov test tended to have lowest power. The methods investigated here are likely to play an important role in identifying pathways that contribute to complex traits.

Several methods have been proposed to impute genotypes at untyped markers using observed genotypes and genetic data from a reference panel. We used the Genetic Analysis Workshop 16 rheumatoid arthritis case-control dataset to compare the performance of four of these imputation methods: IMPUTE, MACH, PLINK, and fastPHASE. We compared the methods' imputation error rates and performance of association tests using the imputed data, in the context of imputing completely untyped markers as well as imputing missing genotypes to combine two datasets genotyped at different sets of markers. As expected, all methods performed better for single-nucleotide polymorphisms (SNPs) in high linkage disequilibrium with genotyped SNPs. However, MACH and IMPUTE generated lower imputation error rates than fastPHASE and PLINK. Association tests based on allele "dosage" from MACH and tests based on the posterior probabilities from IMPUTE provided results closest to those based on complete data. However, in both situations, none of the imputation-based tests provide the same level of evidence of association as the complete data at SNPs strongly associated with disease.

Genome-wide association studies (GWASs) assess correlation between traits and DNA sequence variation using large numbers of genetic variants such as single nucleotide polymorphisms (SNPs) distributed across the genome. A GWAS produces many trait-SNP associations with low p-values, but few are replicated in subsequent studies. We sought to determine if characteristics of the genomic loci associated with a trait could be used to identify initial associations with a higher chance of replication in a second cohort. Data from the age-related eye disease study (AREDS) of 100,000 SNPs on 395 subjects with and 198 without age-related macular degeneration (AMD) were employed. Loci highly associated with AMD were characterized based on the distribution of genotypes, level of significance, and clustering of adjacent SNPs also associated with AMD suggesting linkage disequilibrium or multiple effects. Forty nine loci were highly associated with AMD, including 3 loci (CFH, C2/BF, LOC387715/HTRA1) already known to contain important genetic risks for AMD. One additional locus (C3) reported during the course of this study was identified and replicated in an additional study group. Tag-SNPs and haplotypes for each locus were evaluated for association with AMD in additional cohorts to account for population differences between discovery and replication subjects, but no additional clearly significant associations were identified. Relying on a significant genotype tests using a log-additive model would have excluded 57% of the non-replicated and none of the replicated loci, while use of other SNP features and clustering might have missed true associations.

Akt is a central regulator of cell growth. Its activity can be negatively regulated by the phosphatase PHLPP that specifically dephosphorylates the hydrophobic motif of Akt (Ser473 in Akt1). However, how PHLPP is targeted to Akt is not clear. Here we show that FKBP51 (FK506-binding protein 51) acts as a scaffolding protein for Akt and PHLPP and promotes dephosphorylation of Akt. Furthermore, FKBP51 is downregulated in pancreatic cancer tissue samples and several cancer cell lines. Decreased FKBP51 expression in cancer cells results in hyperphosphorylation of Akt and decreased cell death following genotoxic stress. Overall, our findings identify FKBP51 as a negative regulator of the Akt pathway, with potentially important implications for cancer etiology and response to chemotherapy.

The nuclear factor-kappaB (NF-kappaB) family of transcription factors regulates the expression of a variety of genes involved in apoptosis and immune response. We examined relationships between genotypes at five NF-kappaB subunits (NFKB1, NFKB2, REL, RELA, and RELB) and variable expression levels of 15 NF-kappaB regulated proteins with heritability greater than 0.40: BCL2A1, BIRC2, CD40, CD44, CD80, CFLAR, CR2, FAS, ICAM1, IL15, IRF1, JUNB, MYC, SLC2A5, and VCAM1. SNP genotypes and expression phenotypes from pedigrees of Utah residents with ancestry from northern and western Europe were provided by Genetic Analysis Workshop 15 and supplemented with additional genotype data from the International HapMap Consortium. We conducted association, linkage, and family-based association analyses between each candidate gene and the 15 heritable expression phenotypes. We observed consistent results in association and linkage analyses of the NFKB1 region (encoding p50) and levels of FAS and IRF1 expression. FAS is a cell surface protein that also belongs to the TNF-receptor family; signals through FAS are able to induce apoptosis. IRF1 is a member of the interferon regulatory transcription factor family, which has been shown to regulate apoptosis and tumor-suppression. Analyses in the REL region (encoding c-Rel) revealed linkage and association with CD40 phenotype. CD40 proteins belong to the tumor necrosis factor (TNF)-receptor family, which mediates a broad variety of immune and inflammatory responses. We conclude that variation in the genes encoding p50 and c-Rel may play a role in NF-kappaB-related transcription of FAS, IRF1, and CD40.

Several methods to identify tagging single-nucleotide polymorphisms (SNPs) are in common use for genetic epidemiologic studies; however, there may be loss of information when using only a subset of SNPs. We sought to compare the ability of commonly used pairwise, multimarker, and haplotype-based tagging SNP selection methods to detect known associations with quantitative expression phenotypes. Using data from HapMap release 21 on unrelated Utah residents with ancestors from northern and western Europe (CEPH-Utah, CEU), we selected tagging SNPs in five chromosomal regions using ldSelect, Tagger, and TagSNPs. We found that SNP subsets did not substantially overlap, and that the use of trio data did not greatly impact SNP selection. We then tested associations between HapMap genotypes and expression phenotypes on 28 CEU individuals as part of Genetic Analysis Workshop 15. Relative to the use of all SNPs (n = 210 SNPs across all regions), most subset methods were able to detect single-SNP and haplotype associations. Generally, pairwise selection approaches worked extremely well, relative to use of all SNPs, with marked reductions in the number of SNPs required. Haplotype-based approaches, which had identified smaller SNP subsets, missed associations in some regions. We conclude that the optimal tagging SNP method depends on the true model of the genetic association (i.e., whether a SNP or haplotype is responsible); unfortunately, this is often unknown at the time of SNP selection. Additional evaluations using empirical and simulated data are needed.