Search

© 2014 Nature America, Inc. All rights reserved.Genome-wide association studies (GWAS) have identified 76 variants associated with prostate cancer risk predominantly in populations of European ancestry. To identify additional susceptibility loci for this c

Motivation: motifbreakR scans genetic variants against position weight matrices of transcription factors (TFs) to determine the potential for the disruption of binding at the site of the variant. It leverages the Bioconductor suite of software packages and annotations to query a diverse array of genomes and motif databases. Initially developed to interrogate the effect of single-nucleotide variants on TF binding sites, in motifbreakR v2, we have updated the functionality., Results: New features include the ability to query other types of complex genetic variants, such as short insertions and deletions. This capability allows modeling a more extensive array of variants that may have significant effects on TF binding. Additionally, predictions based on sequence preference alone can indicate many more potential binding events than observed. Adding information from DNA-binding sequencing datasets lends confidence to motif disruption prediction by demonstrating TF binding in cell lines and tissue types. Therefore, motifbreakR can directly query the ReMap2022 database for evidence that a TF matching the disrupted motif binds over the disrupting variant. Finally, in motifbreakR , in addition to the existing interface, we implemented an R/Shiny graphical user interface to simplify and enhance access to researchers with different skill sets., Availability and Implementation: motifbreakR is implemented in R. Source code, documentation, and tutorials are available on Bioconductor at https://bioconductor.org/packages/release/bioc/html/motifbreakR.html and GitHub at https://github.com/Simon-Coetzee/motifBreakR., Competing Interests: None declared., (© The Author(s) 2024. Published by Oxford University Press.)

Background: With recent advances in single cell technology, high-throughput methods provide unique insight into disease mechanisms and more importantly, cell type origin. Here, we used multi-omics data to understand how genetic variants from genome-wide association studies influence development of disease. We show in principle how to use genetic algorithms with normal, matching pairs of single-nucleus RNA- and ATAC-seq, genome annotations, and protein-protein interaction data to describe the genes and cell types collectively and their contribution to increased risk., Results: We used genetic algorithms to measure fitness of gene-cell set proposals against a series of objective functions that capture data and annotations. The highest information objective function captured protein-protein interactions. We observed significantly greater fitness scores and subgraph sizes in foreground vs. matching sets of control variants. Furthermore, our model reliably identified known targets and ligand-receptor pairs, consistent with prior studies., Conclusions: Our findings suggested that application of genetic algorithms to association studies can generate a coherent cellular model of risk from a set of susceptibility variants. Further, we showed, using breast cancer as an example, that such variants have a greater number of physical interactions than expected due to chance., Competing Interests: Competing interests There are no financial and non-financial competing interests declared by authors. Additional Declarations: No competing interests reported.

MotifbreakR is a software tool that scans genetic variants against position weight matrices of transcription factors (TF) to determine the potential for the disruption of TF binding at the site of the variant. It leverages the Bioconductor suite of software packages and annotations to operate across a diverse array of genomes and motif databases. Initially developed to interrogate the effect of single nucleotide variants (common and rare SNVs) on potential TF binding sites, in motifbreakR v2, we have updated the functionality. New features include the ability to query other types of more complex genetic variants, such as short insertions and deletions (indels). This function allows modeling a more extensive array of variants that may have more significant effects on TF binding. Additionally, while TF binding is based partly on sequence preference, predictions of TF binding based on sequence preference alone can indicate many more potential binding events than observed. Adding information from DNA-binding sequencing datasets lends confidence to motif disruption prediction by demonstrating TF binding in cell lines and tissue types. Therefore, motifbreakR implements querying the ReMap2022 database for evidence that a TF matching the disrupted motif binds over the disrupting variant. Finally, in motifbreakR , in addition to the existing interface, we have implemented an R/Shiny graphical user interface to simplify and enhance access to researchers with different skill sets., Competing Interests: Declaration of competing interests None to Declare.

Cancer is an organism-level disease, impacting processes from cellular metabolism and the microenvironment to systemic immune response. Nevertheless, efforts to distinguish overarching mutational processes from interactions with the cell of origin for a tumor have seen limited success, presenting a barrier to individualized medicine. Here we present a pathway-centric approach, extracting somatic mutational profiles within and between tissues, largely orthogonal to cell of origin, mutational burden, or stage. Known predisposition variants are equally distributed among clusters, and largely independent of molecular subtype. Prognosis and risk of death vary jointly by cancer type and cluster. Analysis of metastatic tumors reveals that differences are largely cluster-specific and complementary, implicating convergent mechanisms that combine familiar driver genes with diverse low-frequency lesions in tumor-promoting pathways, ultimately producing distinct molecular phenotypes. The results shed new light on the interplay between organism-level dysfunction and tissue-specific lesions., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Quantifying the functional effects of complex disease risk variants can provide insights into mechanisms underlying disease biology. Genome-wide association studies have identified 39 regions associated with risk of epithelial ovarian cancer (EOC). The vast majority of these variants lie in the non-coding genome, where they likely function through interaction with gene regulatory elements. In this study we first estimated the heritability explained by known common low penetrance risk alleles for EOC. The narrow sense heritability (h g 2 ) of EOC overall and high-grade serous ovarian cancer (HGSOCs) were estimated to be 5%-6%. Partitioned SNP heritability across broad functional categories indicated a significant contribution of regulatory elements to EOC heritability. We collated epigenomic profiling data for 77 cell and tissue types from Roadmap Epigenomics and ENCODE, and from H3K27Ac ChIP-seq data generated in 26 ovarian cancer and precursor-related cell and tissue types. We identified significant enrichment of risk single-nucleotide polymorphisms (SNPs) in active regulatory elements marked by H3K27Ac in HGSOCs. To further investigate how risk SNPs in active regulatory elements influence predisposition to ovarian cancer, we used motifbreakR to predict the disruption of transcription factor binding sites. We identified 469 candidate causal risk variants in H3K27Ac peaks that are predicted to significantly break transcription factor (TF) motifs. The most frequently broken motif was REST (p value = 0.0028), which has been reported as both a tumor suppressor and an oncogene. Overall, these systematic functional annotations with epigenomic data improve interpretation of EOC risk variants and shed light on likely cells of origin., (Copyright © 2020. Published by Elsevier Inc.)

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (p adj  < 8 × 10 -4 ). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (p adj  < 2 × 10 -4 ), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10 -16 ). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: The development of next generation sequencing (NGS) methods led to a rapid rise in the generation of large genomic datasets, but the development of user-friendly tools to analyze and visualize these datasets has not developed at the same pace. This presents a two-fold challenge to biologists; the expertise to select an appropriate data analysis pipeline, and the need for bioinformatics or programming skills to apply this pipeline. The development of graphical user interface (GUI) applications hosted on web-based servers such as Shiny can make complex workflows accessible across operating systems and internet browsers to those without programming knowledge., Results: We have developed GENAVi (Gene Expression Normalization Analysis and Visualization) to provide a user-friendly interface for normalization and differential expression analysis (DEA) of human or mouse feature count level RNA-Seq data. GENAVi is a GUI based tool that combines Bioconductor packages in a format for scientists without bioinformatics expertise. We provide a panel of 20 cell lines commonly used for the study of breast and ovarian cancer within GENAVi as a foundation for users to bring their own data to the application. Users can visualize expression across samples, cluster samples based on gene expression or correlation, calculate and plot the results of principal components analysis, perform DEA and gene set enrichment and produce plots for each of these analyses. To allow scalability for large datasets we have provided local install via three methods. We improve on available tools by offering a range of normalization methods and a simple to use interface that provides clear and complete session reporting and for reproducible analysis., Conclusion: The development of tools using a GUI makes them practical and accessible to scientists without bioinformatics expertise, or access to a data analyst with relevant skills. While several GUI based tools are currently available for RNA-Seq analysis we improve on these existing tools. This user-friendly application provides a convenient platform for the normalization, analysis and visualization of gene expression data for scientists without bioinformatics expertise.

Motivation: DNA methylation has been used to identify functional changes at transcriptional enhancers and other cis-regulatory modules (CRMs) in tumors and other disease tissues. Our R/Bioconductor package ELMER (Enhancer Linking by Methylation/Expression Relationships) provides a systematic approach that reconstructs altered gene regulatory networks (GRNs) by combining enhancer methylation and gene expression data derived from the same sample set., Results: We present a completely revised version 2 of ELMER that provides numerous new features including an optional web-based interface and a new Supervised Analysis mode to use pre-defined sample groupings. We show that Supervised mode significantly increases statistical power and identifies additional GRNs and associated Master Regulators, such as SOX11 and KLF5 in Basal-like breast cancer., Availability and Implementation: ELMER v.2 is available as an R/Bioconductor package at http://bioconductor.org/packages/ELMER/., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press.)

Objective: Genome-wide association studies (GWASs) for epithelial ovarian cancer (EOC) have focused largely on populations of European ancestry. We aimed to identify common germline variants associated with EOC risk in Asian women., Methods: Genotyping was performed as part of the OncoArray project. Samples with >60% Asian ancestry were included in the analysis. Genotyping was performed on 533,631 SNPs in 3238 Asian subjects diagnosed with invasive or borderline EOC and 4083 unaffected controls. After imputation, genotypes were available for 11,595,112 SNPs to identify associations., Results: At chromosome 6p25.2, SNP rs7748275 was associated with risk of serous EOC (odds ratio [OR] = 1.34, P = 8.7 × 10 -9 ) and high-grade serous EOC (HGSOC) (OR = 1.34, P = 4.3 × 10 -9 ). SNP rs6902488 at 6p25.2 (r 2  = 0.97 with rs7748275) lies in an active enhancer and is predicted to impact binding of STAT3, P300 and ELF1. We identified additional risk loci with low Bayesian false discovery probability (BFDP) scores, indicating they are likely to be true risk associations (BFDP <10%). At chromosome 20q11.22, rs74272064 was associated with HGSOC risk (OR = 1.27, P = 9.0 × 10 -8 ). Overall EOC risk was associated with rs10260419 at chromosome 7p21.3 (OR = 1.33, P = 1.2 × 10 -7 ) and rs74917072 at chromosome 2q37.3 (OR = 1.25, P = 4.7 × 10 -7 ). At 2q37.3, expression quantitative trait locus analysis in 404 HGSOC tissues identified ESPNL as a putative candidate susceptibility gene (P = 1.2 × 10 -7 )., Conclusion: While some risk loci were shared between East Asian and European populations, others were population-specific, indicating that the landscape of EOC risk in Asian women has both shared and unique features compared to women of European ancestry., (Copyright © 2019. Published by Elsevier Inc.)

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro , verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439 ., (©2018 American Association for Cancer Research.)

Treatment of prostate cancer (PC) by androgen suppression promotes the emergence of aggressive variants that are androgen receptor (AR) independent. Here we identify the transcription factor ONECUT2 (OC2) as a master regulator of AR networks in metastatic castration-resistant prostate cancer (mCRPC). OC2 acts as a survival factor in mCRPC models, suppresses the AR transcriptional program by direct regulation of AR target genes and the AR licensing factor FOXA1, and activates genes associated with neural differentiation and progression to lethal disease. OC2 appears active in a substantial subset of human prostate adenocarcinoma and neuroendocrine tumors. Inhibition of OC2 by a newly identified small molecule suppresses metastasis in mice. These findings suggest that OC2 displaces AR-dependent growth and survival mechanisms in many cases where AR remains expressed, but where its activity is bypassed. OC2 is also a potential drug target in the metastatic phase of aggressive PC.

Background: Recent genome-wide association studies (GWAS) have identified more than 100 loci associated with increased risk of prostate cancer, most of which are in non-coding regions of the genome. Understanding the function of these non-coding risk loci is critical to elucidate the genetic susceptibility to prostate cancer., Results: We generate genome-wide regulatory element maps and performed genome-wide chromosome confirmation capture assays (in situ Hi-C) in normal and tumorigenic prostate cells. Using this information, we annotate the regulatory potential of 2,181 fine-mapped prostate cancer risk-associated SNPs and predict a set of target genes that are regulated by prostate cancer risk-related H3K27Ac-mediated loops. We next identify prostate cancer risk-associated CTCF sites involved in long-range chromatin loops. We use CRISPR-mediated deletion to remove prostate cancer risk-associated CTCF anchor regions and the CTCF anchor regions looped to the prostate cancer risk-associated CTCF sites, and we observe up to 100-fold increases in expression of genes within the loops when the prostate cancer risk-associated CTCF anchor regions are deleted., Conclusions: We identify GWAS risk loci involved in long-range loops that function to repress gene expression within chromatin loops. Our studies provide new insights into the genetic susceptibility to prostate cancer.

PAX8 is a lineage-restricted transcription factor that is expressed in epithelial ovarian cancer (EOC) precursor tissues, and in the major EOC histotypes. Frequent overexpression of PAX8 in primary EOCs suggests this factor functions as an oncogene during tumorigenesis, however, the biological role of PAX8 in EOC development is poorly understood. We found that stable knockdown of PAX8 in EOC models significantly reduced cell proliferation and anchorage dependent growth in vitro, and attenuated tumorigenicity in vivo . Chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq) and transcriptional profiling were used to create genome-wide maps of PAX8 binding and putative target genes. PAX8 binding sites were significantly enriched in promoter regions ( p < 0.05) and superenhancers ( p < 0.05). MEME-ChIP analysis revealed that PAX8 binding sites overlapping superenhancers or enhancers, but not promoters, were enriched for JUND/B and ARNT/AHR motifs. Integrating PAX8 ChIP-seq and gene expression data identified PAX8 target genes through their associations within shared topological association domains. Across two EOC models we identified 62 direct regulatory targets based on PAX8 binding in promoters and 1,330 putative enhancer regulatory targets . SEPW1, which is involved in oxidation-reduction, was identified as a PAX8 target gene in both cell line models. While the PAX8 cistrome exhibits a high degree of cell-type specificity, analyses of PAX8 target genes and putative cofactors identified common molecular targets and partners as candidate therapeutic targets for EOC., Competing Interests: CONFLICT OF INTEREST None to declare.

Granulocyte-monocyte progenitors (GMPs) and monocyte-dendritic cell progenitors (MDPs) produce monocytes during homeostasis and in response to increased demand during infection. Both progenitor populations are thought to derive from common myeloid progenitors (CMPs), and a hierarchical relationship (CMP-GMP-MDP-monocyte) is presumed to underlie monocyte differentiation. Here, however, we demonstrate that mouse MDPs arose from CMPs independently of GMPs, and that GMPs and MDPs produced monocytes via similar but distinct monocyte-committed progenitors. GMPs and MDPs yielded classical (Ly6C hi ) monocytes with gene expression signatures that were defined by their origins and impacted their function. GMPs produced a subset of "neutrophil-like" monocytes, whereas MDPs gave rise to a subset of monocytes that yielded monocyte-derived dendritic cells. GMPs and MDPs were also independently mobilized to produce specific combinations of myeloid cell types following the injection of microbial components. Thus, the balance of GMP and MDP differentiation shapes the myeloid cell repertoire during homeostasis and following infection., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Prostate cancer incidence is 1.6-fold higher in African Americans than in other populations. The risk factors that drive this disparity are unknown and potentially consist of social, environmental, and genetic influences. To investigate the genetic basis of prostate cancer in men of African ancestry, we performed a genome-wide association meta-analysis using two-sided statistical tests in 10 202 case subjects and 10 810 control subjects. We identified novel signals on chromosomes 13q34 and 22q12, with the risk-associated alleles found only in men of African ancestry (13q34: rs75823044, risk allele frequency = 2.2%, odds ratio [OR] = 1.55, 95% confidence interval [CI] = 1.37 to 1.76, P = 6.10 × 10-12; 22q12.1: rs78554043, risk allele frequency = 1.5%, OR = 1.62, 95% CI = 1.39 to 1.89, P = 7.50 × 10-10). At 13q34, the signal is located 5' of the gene IRS2 and 3' of a long noncoding RNA, while at 22q12 the candidate functional allele is a missense variant in the CHEK2 gene. These findings provide further support for the role of ancestry-specific germline variation in contributing to population differences in prostate cancer risk., (© The Author 2017. Published by Oxford University Press.)

Background: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits., Methods: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background., Results: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis., Conclusions: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures., Impact: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR., Competing Interests: There are no conflicts of interest, (©2016 American Association for Cancer Research.)

Background: Genome-wide association studies (GWAS) in European populations have identified genetic risk variants associated with multiple myeloma., Methods: We performed association testing of common variation in eight regions in 1,318 patients with multiple myeloma and 1,480 controls of European ancestry and 1,305 patients with multiple myeloma and 7,078 controls of African ancestry and conducted a meta-analysis to localize the signals, with epigenetic annotation used to predict functionality., Results: We found that variants in 7p15.3, 17p11.2, 22q13.1 were statistically significantly (P < 0.05) associated with multiple myeloma risk in persons of African ancestry and persons of European ancestry, and the variant in 3p22.1 was associated in European ancestry only. In a combined African ancestry-European ancestry meta-analysis, variation in five regions (2p23.3, 3p22.1, 7p15.3, 17p11.2, 22q13.1) was statistically significantly associated with multiple myeloma risk. In 3p22.1, the correlated variants clustered within the gene body of ULK4 Correlated variants in 7p15.3 clustered around an enhancer at the 3' end of the CDCA7L transcription termination site. A missense variant at 17p11.2 (rs34562254, Pro251Leu, OR, 1.32; P = 2.93 × 10 -7 ) in TNFRSF13B encodes a lymphocyte-specific protein in the TNF receptor family that interacts with the NF-κB pathway. SNPs correlated with the index signal in 22q13.1 cluster around the promoter and enhancer regions of CBX7 CONCLUSIONS: We found that reported multiple myeloma susceptibility regions contain risk variants important across populations, supporting the use of multiple racial/ethnic groups with different underlying genetic architecture to enhance the localization and identification of putatively functional alleles., Impact: A subset of reported risk loci for multiple myeloma has consistent effects across populations and is likely to be functional. Cancer Epidemiol Biomarkers Prev; 25(12); 1609-18. ©2016 AACR., (©2016 American Association for Cancer Research.)

Recent genome-wide association studies (GWAS) of Parkinson's disease (PD) revealed at least 26 risk loci, with associated single nucleotide polymorphisms (SNPs) located in non-coding DNA having unknown functions in risk. In order to explore in which cell types these SNPs (and their correlated surrogates at r(2) ≥ 0.8) could alter cellular function, we assessed their location overlap with histone modification regions that indicate transcription regulation in 77 diverse cell types. We found statistically significant enrichment of risk SNPs at 12 loci in active enhancers or promoters. We investigated 4 risk loci in depth that were most significantly enriched (-logeP > 14) and contained 8 putative enhancers in the different cell types. These enriched loci, along with eQTL associations, were unexpectedly present in non-neuronal cell types. These included lymphocytes, mesendoderm, liver- and fat-cells, indicating that cell types outside the brain are involved in the genetic predisposition to PD. Annotating regulatory risk regions within specific cell types may unravel new putative risk mechanisms and molecular pathways that contribute to PD development.

The 8q24 region harbors multiple risk variants for distinct cancers, including >8 for prostate cancer. In this study, we conducted fine mapping of the 8q24 risk region (127.8-128.8Mb) in search of novel associations with common and rare variation in 4853 prostate cancer case patients and 4678 control subjects of African ancestry. All statistical tests were two-sided. We identified three independent associations at P values of less than 5.00×10(-8), all of which were replicated in studies from Ghana and Uganda (combined sample = 5869 case patients, 5615 control subjects; rs114798100: risk allele frequency [RAF] = 0.04, per-allele odds ratio [OR] = 2.31, 95% confidence interval [CI] = 2.04 to 2.61, P = 2.38×10(-40); rs72725879: RAF = 0.33, OR = 1.37, 95% CI = 1.30 to 1.45, P = 3.04×10(-27); and rs111906932: RAF = 0.03, OR = 1.79, 95% CI = 1.53 to 2.08, P = 1.39×10(-13)). Risk variants rs114798100 and rs111906923 are only found in men of African ancestry, with rs111906923 representing a novel association signal. The three variants are located within or near a number of prostate cancer-associated long noncoding RNAs (lncRNAs), including PRNCR1, PCAT1, and PCAT2. These findings highlight ancestry-specific risk variation and implicate prostate-specific lncRNAs at the 8q24 prostate cancer susceptibility region., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Genome-wide association studies (GWAS) have revealed numerous genomic 'hits' associated with complex phenotypes. In most cases these hits, along with surrogate genetic variation as measure by numerous single nucleotide polymorphisms (SNPs) that are in linkage disequilibrium, are not in coding genes making assignment of functionality or causality intractable. Here we propose that fine-mapping along with the matching of risk SNPs at chromatin biofeatures lessen this complexity by reducing the number of candidate functional/causal SNPs. For example, we show here that only on average 2 SNPs per prostate cancer risk locus are likely candidates for functionality/causality; we further propose that this manageable number should be taken forward in mechanistic studies. The candidate SNPs can be looked up for each prostate cancer risk region in 2 recent publications in 2015 (1,2) from our groups.

Unlabelled: Functional annotation represents a key step toward the understanding and interpretation of germline and somatic variation as revealed by genome-wide association studies (GWAS) and The Cancer Genome Atlas (TCGA), respectively. GWAS have revealed numerous genetic risk variants residing in non-coding DNA associated with complex diseases. For sequences that lie within enhancers or promoters of transcription, it is not straightforward to assess the effects of variants on likely transcription factor binding sites. Consequently we introduce motifbreakR, which allows the biologist to judge whether the sequence surrounding a polymorphism or mutation is a good match, and how much information is gained or lost in one allele of the polymorphism or mutation relative to the other. MotifbreakR is flexible, giving a choice of algorithms for interrogation of genomes with motifs from many public sources that users can choose from. MotifbreakR can predict effects for novel or previously described variants in public databases, making it suitable for tasks beyond the scope of its original design. Lastly, it can be used to interrogate any genome curated within bioconductor., Availability and Implementation: https://github.com/Simon-Coetzee/MotifBreakR, www.bioconductor.org., Contact: dennis.hazelett@cshs.org., (© The Author 2015. Published by Oxford University Press.)

Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region., (© The Author 2015. Published by Oxford University Press.)

Understanding the regulatory landscape of the human genome is a central question in complex trait genetics. Most single-nucleotide polymorphisms (SNPs) associated with cancer risk lie in non-protein-coding regions, implicating regulatory DNA elements as functional targets of susceptibility variants. Here, we describe genome-wide annotation of regions of open chromatin and histone modification in fallopian tube and ovarian surface epithelial cells (FTSECs, OSECs), the debated cellular origins of high-grade serous ovarian cancers (HGSOCs) and in endometriosis epithelial cells (EECs), the likely precursor of clear cell ovarian carcinomas (CCOCs). The regulatory architecture of these cell types was compared with normal human mammary epithelial cells and LNCaP prostate cancer cells. We observed similar positional patterns of global enhancer signatures across the three different ovarian cancer precursor cell types, and evidence of tissue-specific regulatory signatures compared to non-gynecological cell types. We found significant enrichment for risk-associated SNPs intersecting regulatory biofeatures at 17 known HGSOC susceptibility loci in FTSECs (P = 3.8 × 10(-30)), OSECs (P = 2.4 × 10(-23)) and HMECs (P = 6.7 × 10(-15)) but not for EECs (P = 0.45) or LNCaP cells (P = 0.88). Hierarchical clustering of risk SNPs conditioned on the six different cell types indicates FTSECs and OSECs are highly related (96% of samples using multi-scale bootstrapping) suggesting both cell types may be precursors of HGSOC. These data represent the first description of regulatory catalogues of normal precursor cells for different ovarian cancer subtypes, and provide unique insights into the tissue specific regulatory variation with respect to the likely functional targets of germline genetic susceptibility variants for ovarian cancer., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Purpose: The pathophysiology of ocular hypertension (OH) leading to primary open-angle glaucoma shares many features with a secondary form of OH caused by treatment with glucocorticoids, but also exhibits distinct differences. In this study, a pharmacogenomics approach was taken to discover candidate genes for this disorder., Methods: A genome-wide association study was performed, followed by an independent candidate gene study, using a cohort enrolled from patients treated with off-label intravitreal triamcinolone, and handling change in IOP as a quantitative trait., Results: An intergenic quantitative trait locus (QTL) was identified at chromosome 6p21.33 near the 5' end of HCG22 that attained the accepted statistical threshold for genome-level significance. The HCG22 transcript, encoding a novel mucin protein, was expressed in trabecular meshwork cells, and expression was stimulated by IL-1, and inhibited by triamcinolone acetate and TGF-β. Bioinformatic analysis defined the QTL as an approximately 4 kilobase (kb) linkage disequilibrium block containing 10 common single nucleotide polymorphisms (SNPs). Four of these SNPs were identified in the National Center for Biotechnology Information (NCBI) GTEx eQTL browser as modifiers of HCG22 expression. Most are predicted to disrupt or improve motifs for transcription factor binding, the most relevant being disruption of the glucocorticoid receptor binding motif. A second QTL was identified within the predicted signal peptide of the HCG22 encoded protein that could affect its secretion. Translation, O-glycosylation, and secretion of the predicted HCG22 protein was verified in cultured trabecular meshwork cells., Conclusions: Identification of two independent QTLs that could affect expression of the HCG22 mucin gene product via two different mechanisms (transcription or secretion) is highly suggestive of a role in steroid-induced OH.

Dysfunction of the RNA-binding protein, TDP-43, is strongly implicated as a causative event in many neurodegenerative diseases including amyotrophic lateral sclerosis (ALS). TDP-43 is normally found in the nucleus and pathological hallmarks of ALS include the presence of cytoplasmic protein aggregates containing TDP-43 and an associated loss of TDP-43 from the nucleus. Loss of nuclear TDP-43 likely contributes to neurodegeneration. Using Drosophila melanogaster to model TDP-43 loss of function, we show that reduced levels of the voltage-gated calcium channel, cacophony, mediate some of the physiological effects of TDP-43 loss. Null mutations in the Drosophila orthologue of TDP-43, named TBPH, resulted in defective larval locomotion and reduced levels of cacophony protein in whole animals and at the neuromuscular junction. Restoring the levels of cacophony in all neurons or selectively in motor neurons rescued these locomotion defects. Using TBPH immunoprecipitation, we showed that TBPH associates with cacophony transcript, indicating that it is likely to be a direct target for TBPH. Loss of TBPH leads to reduced levels of cacophony transcript, possibly due to increased degradation. In addition, TBPH also appears to regulate the inclusion of some alternatively spliced exons of cacophony. If similar effects of cacophony or related calcium channels are found in human ALS patients, these could be targets for the development of pharmacological therapies for ALS., (© 2013 Published by Elsevier B.V.)

Genome-wide association studies (GWAS) have identified 76 variants associated with prostate cancer risk predominantly in populations of European ancestry. To identify additional susceptibility loci for this common cancer, we conducted a meta-analysis of > 10 million SNPs in 43,303 prostate cancer cases and 43,737 controls from studies in populations of European, African, Japanese and Latino ancestry. Twenty-three new susceptibility loci were identified at association P < 5 × 10(-8); 15 variants were identified among men of European ancestry, 7 were identified in multi-ancestry analyses and 1 was associated with early-onset prostate cancer. These 23 variants, in combination with known prostate cancer risk variants, explain 33% of the familial risk for this disease in European-ancestry populations. These findings provide new regions for investigation into the pathogenesis of prostate cancer and demonstrate the usefulness of combining ancestrally diverse populations to discover risk loci for disease.

Background: The precise nature of how cell type specific chromatin structures at enhancer sites affect gene expression is largely unknown. Here we identified cell type specific enhancers coupled with gene expression in two different types of breast epithelial cells, HMEC (normal breast epithelial cells) and MDAMB231 (triple negative breast cancer cell line)., Results: Enhancers were defined by modified neighboring histones [using chromatin immunoprecipitation followed by sequencing (ChIP-seq)] and nucleosome depletion [using formaldehyde-assisted isolation of regulatory elements followed by sequencing (FAIRE-seq)]. Histone modifications at enhancers were related to the expression levels of nearby genes up to 750 kb away. These expression levels were correlated with enhancer status (poised or active), defined by surrounding histone marks. Furthermore, about fifty percent of poised and active enhancers contained nucleosome-depleted regions. We also identified response element motifs enriched at these enhancer sites that revealed key transcription factors (e.g. TP63) likely involved in regulating breast epithelial enhancer-mediated gene expression. By utilizing expression data, potential target genes of more than 600 active enhancers were identified. These genes were involved in proteolysis, epidermis development, cell adhesion, mitosis, cell cycle, and DNA replication., Conclusions: These findings facilitate the understanding of epigenetic regulation specifically, such as the relationships between regulatory elements and gene expression and generally, how breast epithelial cellular phenotypes are determined by cell type specific enhancers.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process.

The human Tar-DNA binding protein, TDP-43, is associated with amyotrophic lateral sclerosis (ALS) and other neurodegenerative disorders. TDP-43 contains two conserved RNA-binding motifs and has documented roles in RNA metabolism, including pre-mRNA splicing and repression of transcription. Here, using Drosophila melanogaster as a model, we generated loss-of-function and overexpression genotypes of Tar-DNA binding protein homolog (TBPH) to study their effect on the transcriptome of the central nervous system (CNS). By using massively parallel sequencing methods (RNA-seq) to profile the CNS, we find that loss of TBPH results in widespread gene activation and altered splicing, much of which are reversed by rescue of TBPH expression. Conversely, TBPH overexpression results in decreased gene expression. Although previous studies implicated both absence and mis-expression of TDP-43 in ALS, our data exhibit little overlap in the gene expression between them, suggesting that the bulk of genes affected by TBPH loss-of-function and overexpression are different. In combination with computational approaches to identify likely TBPH targets and orthologs of previously identified vertebrate TDP-43 targets, we provide a comprehensive analysis of enriched gene ontologies. Our data suggest that TDP-43 plays a role in synaptic transmission, synaptic release, and endocytosis. We also uncovered a potential novel regulation of the Wnt and BMP pathways, many of whose targets appear to be conserved.

Pde1c is a calcium/calmodulin-regulated, dual-specificity cyclic nucleotide phosphodiesterase. We have used a transposon insertion line to investigate the physiological function of Pde1c in Drosophila melanogaster and to show that the insertion leads to male sterility and male mating behavior defects that include reduced copulation rates. Sterility appears to be primarily due to elimination of sperm from the female reproductive system. The male mating behavior defects were fully rescued by expression of exogenous Pde1c under the control of either a Pde1c or a pan-neuronal promoter, whereas the sterility could be only partially rescued by expression of exogenous Pde1c under the control of these promoters. We also show that Pde1c has a male-specific expression pattern in the CNS with an increased number of Pde1c-expressing neurons in the abdominal ganglion in males.

Methods: A new method was developed for identifying novel transcription factor regulatory targets based on calculating Local Affinity Density. Techniques from the signal-processing field were used, in particular the Hann digital filter, to calculate the relative binding affinity of different regions based on previously published in vitro binding data. To illustrate this approach, the complete genomes of Drosophila melanogaster and D.pseudoobscura were analyzed for binding sites of the homeodomain proteinc Tinman, an essential heart development gene in both Drosophila and Mouse. The significant binding regions were identified relative to genomic background and assigned to putative target genes. Valid candidates common to both species of Drosophila were selected as a test of conservation., Results: The new method was more sensitive than cluster searches for conserved binding motifs with respect to positive identification of known Tinman targets. Our Local Affinity Density method also identified a significantly greater proportion of Tinman-coexpressed genes than equivalent, optimized cluster searching. In addition, this new method predicted a significantly greater than expected number of genes with previously published RNAi phenotypes in the heart., Availability: Algorithms were implemented in Python, LISP, R and maxima, using MySQL to access locally mirrored sequence data from Ensembl (D.melanogaster release 4.3) and flybase (D.pseudoobscura). All code is licensed under GPL and freely available at http://www.ohsu.edu/cellbio/dev&#95;biol&#95;prog/affinitydensity/.

During metamorphosis of the hawkmoth, Manduca sexta, some larval muscles degenerate while others are respecified for new functions. In larvae, accessory planta retractor muscles (APRMs) are present in abdominal segments 1 to 6 (A1 to A6). APRMs serve as proleg retractors in A3 to A6 and body wall muscles in A1 and A2. At pupation, all APRMs degenerate except those in A2 and A3, which are respecified to circulate hemolymph in pupae. The motoneurons that innervate APRMs, the APRs, likewise undergo segment-specific programmed cell death (PCD), as a direct, cell-autonomous response to the prepupal peak of ecdysteroids. The segment-specific patterns of APR and APRM death differ. The present study tested the hypothesis that APRM death is a direct, cell-autonomous response to the prepupal peak of ecdysteroids. Prevention of the prepupal peak prevented APRM degeneration, and replacement of the peak by infusion of 20-hydroxyecdysone restored the correct segment-specific pattern of APRM degeneration. Surgical denervation of APRMs did not perturb their segment-specific degeneration at pupation, indicating that signals from APRs are not required for the muscles' segment-specific responses to ecdysteroids. The possibility that instructive signals originate from APRMs' epidermal attachment points was tested by treating the epidermis with a juvenile hormone analog to prevent pupal development. This manipulation likewise did not alter APRM fate. We conclude that both the muscles and motoneurons in this motor system respond directly and cell-autonomously to prepupal ecdysteroids to produce a segment-specific pattern of PCD that is matched to the functional requirements of the pupal body., (2004 Wiley Periodicals, Inc)

The wave of differentiation that traverses the Drosophila eye disc requires rapid transitions in gene expression that are controlled by a number of signaling molecules also required in other developmental processes. We have used a mosaic genetic screen to systematically identify autosomal genes required for the normal pattern of photoreceptor differentiation, independent of their requirements for viability. In addition to genes known to be important for eye development and to known and novel components of the Hedgehog, Decapentaplegic, Wingless, Epidermal growth factor receptor, and Notch signaling pathways, we identified several members of the Polycomb and trithorax classes of genes encoding general transcriptional regulators. Mutations in these genes disrupt the transitions between zones along the anterior-posterior axis of the eye disc that express different combinations of transcription factors. Different trithorax group genes have very different mutant phenotypes, indicating that target genes differ in their requirements for chromatin remodeling, histone modification, and coactivation factors.

Cells in the morphogenetic furrow of the Drosophila eye disc undergo a striking shape change immediately prior to their neuronal differentiation. We have isolated mutations in a novel gene, act up (acu), that is required for this shape change. acu encodes a homolog of yeast cyclase-associated protein, which sequesters monomeric actin; we show that acu is required to prevent actin filament polymerization in the eye disc. In contrast, profilin promotes actin filament polymerization, acting epistatically to acu. However, both acu and profilin are required to prevent premature Hedgehog-induced photoreceptor differentiation ahead of the morphogenetic furrow. These findings suggest that dynamic changes in actin filaments alter cell shape to control the movement of signals that coordinate a wave of differentiation.

Signaling by the secreted hedgehog, decapentaplegic and wingless proteins organizes the pattern of photoreceptor differentiation within the Drosophila eye imaginal disc; hedgehog and decapentaplegic are required for differentiation to initiate at the posterior margin and progress across the disc, while wingless prevents it from initiating at the lateral margins. Our analysis of these interactions has shown that initiation requires both the presence of decapentaplegic and the absence of wingless, which inhibits photoreceptor differentiation downstream of the reception of the decapentaplegic signal. However, wingless is unable to inhibit differentiation driven by activation of the epidermal growth factor receptor pathway. The effect of wingless is subject to regional variations in control, as the anterior margin of the disc is insensitive to wingless inhibition. The eyes absent and eyegone genes encode members of a group of nuclear proteins required to specify the fate of the eye imaginal disc. We show that both eyes absent and eyegone are required for normal activation of decapentaplegic expression at the posterior and lateral margins of the disc, and repression of wingless expression in presumptive retinal tissue. The requirement for eyegone can be alleviated by inhibition of the wingless signaling pathway, suggesting that eyegone promotes eye development primarily by repressing wingless. These results provide a link between the early specification and later differentiation of the eye disc.