Search

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., (© The Author(s) 2022. 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.)

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.

Genome-wide association studies (GWAS) of breast cancer defined by hormone receptor status have revealed loci contributing to susceptibility of estrogen receptor (ER)-negative subtypes. To identify additional genetic variants for ER-negative breast cancer, we conducted the largest meta-analysis of ER-negative disease to date, comprising 4754 ER-negative cases and 31 663 controls from three GWAS: NCI Breast and Prostate Cancer Cohort Consortium (BPC3) (2188 ER-negative cases; 25 519 controls of European ancestry), Triple Negative Breast Cancer Consortium (TNBCC) (1562 triple negative cases; 3399 controls of European ancestry) and African American Breast Cancer Consortium (AABC) (1004 ER-negative cases; 2745 controls). We performed in silico replication of 86 SNPs at P ≤ 1 × 10(-5) in an additional 11 209 breast cancer cases (946 with ER-negative disease) and 16 057 controls of Japanese, Latino and European ancestry. We identified two novel loci for breast cancer at 20q11 and 6q14. SNP rs2284378 at 20q11 was associated with ER-negative breast cancer (combined two-stage OR = 1.16; P = 1.1 × 10(-8)) but showed a weaker association with overall breast cancer (OR = 1.08, P = 1.3 × 10(-6)) based on 17 869 cases and 43 745 controls and no association with ER-positive disease (OR = 1.01, P = 0.67) based on 9965 cases and 22 902 controls. Similarly, rs17530068 at 6q14 was associated with breast cancer (OR = 1.12; P = 1.1 × 10(-9)), and with both ER-positive (OR = 1.09; P = 1.5 × 10(-5)) and ER-negative (OR = 1.16, P = 2.5 × 10(-7)) disease. We also confirmed three known loci associated with ER-negative (19p13) and both ER-negative and ER-positive breast cancer (6q25 and 12p11). Our results highlight the value of large-scale collaborative studies to identify novel breast cancer risk loci.

The androgen receptor (AR) gene contains a polymorphic trinucleotide repeat region, (CAG)(n), in its N-terminal transactivation domain (NTD) that encodes a polyglutamine (polyQ) tract in the receptor protein. Whereas the length of the CAG repeat ranges from 6 to 39 in healthy individuals, the variations in repeat length both within and outside the normal range are associated with disease, including impaired spermatogenesis and Kennedy's disease, and with the risk of developing breast and prostate cancer. Whereas it has been proposed that the inverse relationship between polyQ tract length within the normal range and AR transactivation potential may be responsible for altered risk of disease, the molecular mechanisms underlying polyQ length modulation of AR function have not been elucidated. In this study, we provide detailed characterization of a somatic AR gene mutation detected in a human prostate tumor that results in interruption of the polyQ tract by two non-consecutive leucine residues (AR-polyQ2L). Compared with wtAR, AR-polyQ2L exhibits disrupted inter-domain communication (N/C interaction) and a lower protein level, but paradoxically has markedly increased transactivation activity. Molecular modeling and the response to cofactors indicate that the increased activity of AR-polyQ2L results from the presentation of a more stable platform for the recruitment of accessory proteins than wild-type AR. Analysis of the relationship between polyQ tract length and AR function revealed a critical size (Q16-Q29) for maintenance of N/C interaction. That between 91 and 99% of AR alleles in different racial-ethnic groups encode a polyQ tract in the range of Q16-Q29 suggests that N/C interaction has been preserved as an essential component of androgen-induced AR signaling.

Normal polymorphic size variation of the exon 1 CAG microsatellite of the androgen receptor (AR) is associated with prostate cancer, benign prostatic hyperplasia and male infertility. Furthermore, abnormal expansion of the satellite leads to Kennedy's disease. We have shown recently that the AR N-terminal domain (NTD), which contains the polyglutamine (polyQ) stretch (encoded by the CAG repeat), functionally interacts with the C-termini of p160 coactivators. In the present study we explored possible AR CAG size effects on the p160 coactivator-mediated transactivation activity of the receptor. First, we mapped the p160 coactivator interaction on the AR NTD and found an interaction surface between amino acids 351 and 537. Although this region is 'downstream' from the polyQ stretch, it is still within the AR NTD, is implicated in constitutive transactivation activity of the receptor, and thus might be subject to polyQ size modulation. Indeed, cotrans- fection experiments in cultured prostate epithelial cells, using AR constructs of varying CAG sizes and p160 coactivator expression vectors, revealed that increased polyQ length, up to a size of 42 repeats, inhibited both basal and coactivator-mediated AR transactivation activity. AR expression in these cells, on the other hand, was unaffected by the same increased CAG repeat size range. We conclude that the AR NTD contributes to AR transactivation activity via functional interactions with p160 coactivators and that increasing polyQ length negatively affects p160-mediated coactivation of the AR. This molecular mechanism thus might explain, at least in part, the observed phenotypic effects of the AR CAG size polymorphism.

A subgroup of trinucleotide repeat diseases result from abnormal expansions of CAG repeats which are translated into polyglutamine stretches. As yet there is little understanding of how the polyglutamines function either normally, or when expanded. We have investigated these sequences in the Machado-Joseph disease, androgen receptor and spinocerebellar ataxia type 1 genes in humans and other primates. None of the 748 normal chromosomes that were examined had more than 34 uninterrupted glutamine codons in the Machado-Joseph disease gene. Similarly, no normal alleles with more than 39 uninterrupted glutamine codons have been reported for the other disease genes associated with polyglutamine expansions. Sequence analyses of the repeats in primates revealed shorter polyglutamine stretches in some of the non-human primates at all three loci and marked diversions from the expected polyglutamines in the orang-utan Machado-Joseph gene and in the marmoset spinocerebellar ataxia type 1 gene. These data suggest that conservation of these polyglutamine stretches may not always be necessary for normal gene function.