The Parkinson's disease variant rs356182 regulates neuronal differentiation independently from alpha-synuclein.

One of the most significant risk variants for Parkinson's disease (PD), rs356182, is located at the PD-associated locus near the alpha-synuclein (α-syn) encoding gene, SNCA. SNCA-proximal variants, including rs356182, are thought to function in PD risk through enhancers via allele-specific regulatory effects on SNCA expression. However, this interpretation discounts the complex activity of genetic enhancers and possible non-conical functions of α-syn. Here we investigated a novel risk mechanism for rs356182. We use CRISPR-Cas9 in LUHMES cells, a model for dopaminergic midbrain neurons, to generate precise hemizygous lesions at rs356182. The PD-protective (A/-), PD-risk (G/-) and wild-type (A/G) clones were neuronally differentiated and then compared transcriptionally and morphologically. Among the affected genes was SNCA, whose expression was promoted by the PD-protective allele (A) and repressed in its absence. In addition to SNCA, hundreds of genes were differentially expressed and associated with neurogenesis and axonogenesis-an effect not typically ascribed to α-syn. We also found that the transcription factor FOXO3 specifically binds to the rs356182 A-allele in differentiated LUHMES cells. Finally, we compared the results from the rs356182-edited cells to our previously published knockouts of SNCA and found only minimal overlap between the sets of significant differentially expressed genes. Together, the data implicate a risk mechanism for rs356182 in which the risk-allele (G) is associated with abnormal neuron development, independent of SNCA expression. We speculate that these pathological effects manifest as a diminished population of dopaminergic neurons during development leading to the predisposition for PD later in life.

Alpha-synuclein negatively controls cell proliferation in dopaminergic neurons.

As researchers grapple with the mechanisms and implications of alpha-synuclein (α-syn) in neuropathology, it is often forgotten that the function(s) of α-syn in healthy cells remain largely elusive. Previous work has relied on observing α-syn localization in the cell or using knockout mouse models. Here, we address the specific role of α-syn in human dopaminergic neurons by disrupting its gene (SNCA) in the human dopaminergic neuron cell line, LUHMES. SNCA-null cells were able to differentiate grossly normally and showed modest effects on gene expression. The effects on gene expression were monodirectional, resulting primarily in the significant decrease of expression for 401 genes, implicating them as direct, or indirect positive targets of α-syn. Gene ontological analysis of these genes showed enrichment in terms associated with proliferation, differentiation, and synapse activity. These results add to the tapestry of α-syn biological functions. SIGNIFICANCE STATEMENT: The normal functions of α-syn have remained controversial, despite its clear importance in Parkinson's Disease pathology, where it accumulates in Lewy bodies and contributes to neurodegeneration. Its name implies synaptic and nuclear functions, but how it participates at these locations has not been resolved. Via knock-out experiments in dopaminergic neurons, we implicate α-syn as a functional participant in synapse activity and in proliferation/differentiation, the latter being novel and provide insight into α-syn's role in neuronal development.

Understanding Non-Mendelian Genetic Risk.

This opinion paper highlights strategies for a better understanding of non-Mendelian genetic risk that was revealed by genome-wide association studies (GWAS) of complex diseases. The genetic risk resides predominantly in non-coding regulatory DNA, such as in enhancers. The identification of mechanisms, the causal variants (mainly SNPs), and their target genes are, however, not always apparent but are likely involved in a network of risk determinants; the identification presents a bottle-neck in the full understanding of the genetics of complex phenotypes. Here, we propose strategies to identify functional SNPs and link risk enhancers with their target genes. The strategies are 1) identifying fine-mapped SNPs that break/form response elements within chromatin bio-features in relevant cell types 2) considering the nearest gene on linear DNA, 3) analyzing eQTLs, 4) mapping differential DNA methylation regions and relating them to gene expression, 5) employing genomic editing with CRISPR/cas9 and 6) identifying topological associated chromatin domains using chromatin conformation capture.

Parkinson's disease genetic risk in a midbrain neuronal cell line.

In genome-wide association studies of complex diseases, many risk polymorphisms are found to lie in non-coding DNA and likely confer risk through allele-dependent differences in gene regulatory elements. However, because distal regulatory elements can alter gene expression at various distances on linear DNA, the identity of relevant genes is unknown for most risk loci. In Parkinson's disease, at least some genetic risk is likely intrinsic to a neuronal subpopulation of cells in the brain regions affected. In order to compare neuron-relevant methods of pairing risk polymorphisms to target genes as well as to further characterize a single-cell model of a neurodegenerative disease, we used the portionally-dopaminergic, neuronal, mesencephalic-derived cell line LUHMES to dissect differentiation-specific mechanisms of gene expression. We compared genome-wide gene expression in undifferentiated and differentiated cells with genome-wide histone H3K27ac and CTCF-bound regions. Whereas promoters and CTCF binding were largely consistent between differentiated and undifferentiated cells, enhancers were mostly unique. We matched the differentiation-specific appearance or disappearance of enhancers with changes in gene expression and identified 22,057 enhancers paired with 6388 differentially expressed genes by proximity. These enhancers are enriched with at least 13 transcription factor response elements, driving a cluster of genes involved in neurogenesis. We show that differentiated LUHMES cells, but not undifferentiated cells, show enrichment for PD-risk SNPs. Candidate genes for these loci are largely unrelated, though a subset is linked to synaptic vesicle cycling and transport, implying that PD-related disruption of these pathways is intrinsic to dopaminergic neurons.

Cell-type-specific enrichment of risk-associated regulatory elements at ovarian cancer susceptibility loci.

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.

Comprehensive functional annotation of 77 prostate cancer risk loci.

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.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1.

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

motifbreakR: an R/Bioconductor package for predicting variant effects at transcription factor binding sites.

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.

Biological implications and regulatory mechanisms of long-range chromosomal interactions.

Development of high-throughput sequencing-based methods has enabled us to examine nuclear architecture at unprecedented resolution, allowing further examination of the function of long-range chromosomal interactions. Here, we review methods used to investigate novel long-range chromosomal interactions and genome-wide organization of chromatin. We further discuss transcriptional activation and silencing in relation to organization and positioning of gene loci and regulation of chromatin organization through protein complexes and noncoding RNAs.

Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells.

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.

Identification of Genetic Susceptibility Loci for Colorectal Tumors in a Genome-Wide Meta-analysis.

BACKGROUND & AIMS: Heritable factors contribute to the development of colorectal cancer. Identifying the genetic loci associated with colorectal tumor formation could elucidate the mechanisms of pathogenesis. METHODS: We conducted a genome-wide association study that included 14 studies, 12,696 cases of colorectal tumors (11,870 cancer, 826 adenoma), and 15,113 controls of European descent. The 10 most statistically significant, previously unreported findings were followed up in 6 studies; these included 3056 colorectal tumor cases (2098 cancer, 958 adenoma) and 6658 controls of European and Asian descent. RESULTS: Based on the combined analysis, we identified a locus that reached the conventional genome-wide significance level at less than 5.0 × 10(-8): an intergenic region on chromosome 2q32.3, close to nucleic acid binding protein 1 (most significant single nucleotide polymorphism: rs11903757; odds ratio [OR], 1.15 per risk allele; P = 3.7 × 10(-8)). We also found evidence for 3 additional loci with P values less than 5.0 × 10(-7): a locus within the laminin gamma 1 gene on chromosome 1q25.3 (rs10911251; OR, 1.10 per risk allele; P = 9.5 × 10(-8)), a locus within the cyclin D2 gene on chromosome 12p13.32 (rs3217810 per risk allele; OR, 0.84; P = 5.9 × 10(-8)), and a locus in the T-box 3 gene on chromosome 12q24.21 (rs59336; OR, 0.91 per risk allele; P = 3.7 × 10(-7)). CONCLUSIONS: In a large genome-wide association study, we associated polymorphisms close to nucleic acid binding protein 1 (which encodes a DNA-binding protein involved in DNA repair) with colorectal tumor risk. We also provided evidence for an association between colorectal tumor risk and polymorphisms in laminin gamma 1 (this is the second gene in the laminin family to be associated with colorectal cancers), cyclin D2 (which encodes for cyclin D2), and T-box 3 (which encodes a T-box transcription factor and is a target of Wnt signaling to ß-catenin). The roles of these genes and their products in cancer pathogenesis warrant further investigation.

Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells.

Expression profiling of T24 cells revealed that 17 out of 313 human miRNAs were upregulated more than 3-fold by simultaneous treatment with the chromatin-modifying drugs 5-aza-2'-deoxycytidine and 4-phenylbutyric acid. One of these, miR-127, is embedded in a CpG island and is highly induced from its own promoter after treatment. miR-127 is usually expressed as part of a miRNA cluster in normal cells but not in cancer cells, suggesting that it is subject to epigenetic silencing. In addition, the proto-oncogene BCL6, a potential target of miR-127, was translationally downregulated after treatment. These results suggest that DNA demethylation and histone deacetylase inhibition can activate expression of miRNAs that may act as tumor suppressors.

Genome-wide Runx2 occupancy in prostate cancer cells suggests a role in regulating secretion.

Runx2 is a metastatic transcription factor (TF) increasingly expressed during prostate cancer (PCa) progression. Using PCa cells conditionally expressing Runx2, we previously identified Runx2-regulated genes with known roles in epithelial-mesenchymal transition, invasiveness, angiogenesis, extracellular matrix proteolysis and osteolysis. To map Runx2-occupied regions (R2ORs) in PCa cells, we first analyzed regions predicted to bind Runx2 based on the expression data, and found that recruitment to sites upstream of the KLK2 and CSF2 genes was cyclical over time. Genome-wide ChIP-seq analysis at a time of maximum occupancy at these sites revealed 1603 high-confidence R2ORs, enriched with cognate motifs for RUNX, GATA and ETS TFs. The R2ORs were distributed with little regard to annotated transcription start sites (TSSs), mainly in introns and intergenic regions. Runx2-upregulated genes, however, displayed enrichment for R2ORs within 40 kb of their TSSs. The main annotated functions enriched in 98 Runx2-upregulated genes with nearby R2ORs were related to invasiveness and membrane trafficking/secretion. Indeed, using SDS-PAGE, mass spectrometry and western analyses, we show that Runx2 enhances secretion of several proteins, including fatty acid synthase and metastasis-associated laminins. Thus, combined analysis of Runx2's transcriptome and genomic occupancy in PCa cells lead to defining its novel role in regulating protein secretion.

FunciSNP: an R/bioconductor tool integrating functional non-coding data sets with genetic association studies to identify candidate regulatory SNPs.

UNLABELLED: Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance. AVAILABILITY: FunciSNP is available from Bioconductor (bioconductor.org).

Parkinson's disease risk enhancers in microglia.

Genome-wide association studies have identified thousands of single nucleotide polymorphisms that associate with increased risk for Parkinson's disease (PD), but the functions of most of them are unknown. Using assay for transposase-accessible chromatin (ATAC) and H3K27ac chromatin immunoprecipitation (ChIP) sequencing data, we identified 73 regulatory elements in microglia that overlap PD risk SNPs. To determine the target genes of a "risk enhancer" within intron two of SNCA, we used CRISPR-Cas9 to delete the open chromatin region where two PD risk SNPs reside. The loss of the enhancer led to reduced expression of multiple genes including SNCA and the adjacent gene MMRN1. It also led to expression changes of genes involved in glucose metabolism, a process that is known to be altered in PD patients. Our work expands the role of SNCA in PD and provides a connection between PD-associated genetic variants and underlying biology that points to a risk mechanism in microglia.

Gene expression asymmetry in Parkinson's Disease; variation of CCT and BEX gene expression levels are correlated with hemisphere specific severity.

Parkinson's Disease (PD) develops unilaterally, which may be related to brain hemispheric differences in gene expression. Here we measured bulk RNA-seq levels in neuronal nuclei obtained from prefrontal cortex postmortem brain samples from males and females with PD and from healthy controls. Left and right hemispheres from each brain were related the side of symptom onset and compared. We employed two a priori approaches; first we identified genes differentially expressed between PD and controls and between left vs right PD brain hemispheres. Second, we examined the presence of, and correlates to, variable asymmetry seen in candidate PD differentially expressed genes. We found large variation among individuals with PD, and PD stratification by gene expression similarity was required for patterns of genetic asymmetry to emerge. For a subset of PD brains, hemispherical variation of CCT and BEX gene levels correlated with the side of PD symptom onset.

Knockdown of the cochaperone SGTA results in the suppression of androgen and PI3K/Akt signaling and inhibition of prostate cancer cell proliferation.

Solid tumors have an increased reliance on Hsp70/Hsp90 molecular chaperones for proliferation, survival and maintenance of intracellular signaling systems. An underinvestigated component of the chaperone system is the tetratricopeptide repeat (TPR)-containing cochaperone, which coordinates Hsp70/Hsp90 involvement on client proteins as well as having diverse individual actions. A potentially important cochaperone in prostate cancer (PCa) is small glutamine-rich TPR-containing protein alpha (SGTA), which interacts with the androgen receptor (AR) and other critical cancer-related client proteins. In this study, the authors used small interfering RNA coupled with genome-wide expression profiling to investigate the biological significance of SGTA in PCa and its influence on AR signaling. Knockdown of SGTA for 72 hr in PCa C4-2B cells significantly altered expression of >1,900 genes (58% decreased) and reduced cell proliferation (p < 0.05). The regulation of 35% of 5α-dihydrotestosterone (DHT) target genes was affected by SGTA knockdown, with gene-specific effects on basal or DHT-induced expression or both. Pathway analysis revealed a role for SGTA in p53, generic PCa and phosphoinositol kinase (PI3K) signaling pathways; the latter evident by a reduction in PI3K subunit p100ß levels and decreased phosphorylated Akt. Immunohistochemical analysis of 64 primary advanced PCa samples showed a significant increase in the AR:SGTA ratio in cancerous lesions compared to patient-matched benign prostatic hyperplasia tissue (p < 0.02). This study not only provides insight into the biological actions of SGTA and its effect on genome-wide AR transcriptional activity and other therapeutically targeted intracellular signaling pathways but also provides evidence for PCa-specific alterations in SGTA expression.

8q24 prostate, breast, and colon cancer risk loci show tissue-specific long-range interaction with MYC.

The 8q24 gene desert contains risk loci for multiple epithelial cancers, including colon, breast, and prostate. Recent evidence suggests these risk loci contain enhancers. In this study, data are presented showing that each risk locus bears epigenetic marks consistent with enhancer elements and forms a long-range chromatin loop with the MYC proto-oncogene located several hundred kilobases telomeric and that these interactions are tissue-specific. We therefore propose that the 8q24 risk loci operate through a common mechanism-as tissue-specific enhancers of MYC.

Functional enhancers at the gene-poor 8q24 cancer-linked locus.

Multiple discrete regions at 8q24 were recently shown to contain alleles that predispose to many cancers including prostate, breast, and colon. These regions are far from any annotated gene and their biological activities have been unknown. Here we profiled a 5-megabase chromatin segment encompassing all the risk regions for RNA expression, histone modifications, and locations occupied by RNA polymerase II and androgen receptor (AR). This led to the identification of several transcriptional enhancers, which were verified using reporter assays. Two enhancers in one risk region were occupied by AR and responded to androgen treatment; one contained a single nucleotide polymorphism (rs11986220) that resides within a FoxA1 binding site, with the prostate cancer risk allele facilitating both stronger FoxA1 binding and stronger androgen responsiveness. The study reported here exemplifies an approach that may be applied to any risk-associated allele in non-protein coding regions as it emerges from genome-wide association studies to better understand the genetic predisposition of complex diseases.