Shared common variants in prostate cancer and blood lipids

BACKGROUND: Epidemiological and clinical studies suggest comorbidity between prostate cancer (PCA) and cardiovascular disease (CVD) risk factors. However, the relationship between these two phenotypes is still not well understood. Here we sought to identify shared genetic loci between PCA and CVD risk factors.

METHODS: We applied a genetic epidemiology method based on conjunction false discovery rate (FDR) that combines summary statistics from different genome-wide association studies (GWAS), and allows identification of genetic overlap between two phenotypes. We evaluated summary statistics from large, multi-centre GWA studies of PCA (n=50 000) and CVD risk factors (n=200 000) [triglycerides (TG), low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, body mass index, waist-hip ratio and type 2 diabetes (T2D)]. Enrichment of single nucleotide polymorphisms (SNPs) associated with PCA and CVD risk factors was assessed with conditional quantile-quantile plots and the Anderson-Darling test. Moreover, we pinpointed shared loci using conjunction FDR.

RESULTS: We found the strongest enrichment of P-values in PCA was conditional on LDL and conditional on TG. In contrast, we found only weak enrichment conditional on HDL or conditional on the other traits investigated. Conjunction FDR identified altogether 17 loci; 10 loci were associated with PCA and LDL, 3 loci were associated with PCA and TG and additionally 4 loci were associated with PCA, LDL and TG jointly (conjunction FDR <0.01). For T2D, we detected one locus adjacent to HNF1B.

CONCLUSIONS: We found polygenic overlap between PCA predisposition and blood lipids, in particular LDL and TG, and identified 17 pleiotropic gene loci between PCA and LDL, and PCA and TG, respectively. These findings provide novel pathobiological insights and may have implications for trials using targeting lipid-lowering agents in a prevention or cancer setting.

Exome sequencing of prostate cancer supports the hypothesis of independent tumour origins

BACKGROUND: Prostate cancer (PCa) is a clinically and pathologically heterogeneous disease. The rapid development of sequencing technology has the potential to deliver new biomarkers with emphasis on aggressive disease and to revolutionise personalised cancer treatment. However, a prostate harbouring cancer commonly contains multiple separate tumour foci, with the potential to aggravate tumour sampling. The level of intraprostatic tumour heterogeneity remains to be determined.

OBJECTIVE: To determine the level of intraprostatic tumour heterogeneity through genome-wide, high-resolution profiling of multiple tumour samples from the same individual.

DESIGN, SETTINGS, AND PARTICIPANTS: Multiple tumour samples were obtained from four individuals following radical prostatectomy. One individual (SWE-1) contained >70% cancer cells in all tumour samples, whereas the other three (SWE-2 to SWE-4) required the use of laser capture microdissection for tumour cell enrichment. Subsequently, DNA was extracted from all tissue samples, and exome sequencing was performed. All tumour foci of SWE-1 were also profiled using a high-resolution array for the identification of copy number alterations (CNA).

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Shared somatic high-frequency single nucleotide variants (SNV) and CNAs were used to infer the level of intraprostatic tumour heterogeneity.

RESULTS AND LIMITATIONS: No high-frequency mutations, common for the three tumour samples of SWE-1, were identified. Ten randomly chosen positions were validated with Sanger sequencing in all foci, which verified the exome data. The high level of intraprostatic heterogeneity was consistent in all individuals. In total, three out of four individuals harboured tumours without an apparent common somatic denominator. Although we cannot exclude the presence of common structural rearrangements, a high-density array was used for the detection of deletions and amplifications in SWE-1, which agreed with the exome data.

CONCLUSIONS: We present evidence for the presence of somatically independent tumours within the same prostate. This finding will have implications for personalised cancer treatment and biomarker discovery.

Chromatin binding by the androgen receptor in prostate cancer

Alterations in transcriptional programs are fundamental to the development of cancers. The androgen receptor is central to the normal development of the prostate gland and to the development of prostate cancer. To a large extent this is believed to be due to the control of gene expression through the interaction of the androgen receptor with chromatin and subsequently with coregulators and the transcriptional machinery. Unbiased genome-wide studies have recently uncovered the recruitment sites that are gene-distal and intragenic rather than associated with proximal promoter regions. Whilst expression profiles from AR-positive primary prostate tumours and cell lines can directly relate to the AR cistrome in prostate cancer cells, this distribution raises significant challenges in making direct mechanistic connections. Furthermore, extrapolating from datasets assembled in one model to other model systems or clinical samples poses challenges if we are to use the AR-directed transcriptome to guide the development of novel biomarkers or treatment decisions. This review will provide an overview of the androgen receptor before addressing the challenges and opportunities created by whole-genome studies of the interplay between the androgen receptor and chromatin.

Elevated NCOR1 disrupts a network of dietary-sensing nuclear receptors in bladder cancer cells

Increasingly invasive bladder cancer cells lines displayed insensitivity toward a panel of dietary-derived ligands for members of the nuclear receptor superfamily. Insensitivity was defined through altered gene regulatory actions and cell proliferation and reflected both reduced receptor expression and elevated nuclear receptor corepressor 1 (NCOR1) expression. Stable overexpression of NCOR1 in sensitive cells (RT4) resulted in a panel of clones that recapitulated the resistant phenotype in terms of gene regulatory actions and proliferative responses toward ligand. Similarly, silencing RNA approaches to NCOR1 in resistant cells (EJ28) enhanced ligand gene regulatory and proliferation responses, including those mediated by peroxisome proliferator-activated receptor (PPAR) gamma and vitamin D receptor (VDR) receptors. Elevated NCOR1 levels generate an epigenetic lesion to target in resistant cells using the histone deacetylase inhibitor vorinostat, in combination with nuclear receptor ligands. Such treatments revealed strong-additive interactions toward the PPARgamma, VDR and Farnesoid X-activated receptors. Genome-wide microarray and microfluidic quantitative real-time, reverse transcription-polymerase chain reaction approaches, following the targeting of NCOR1 activity and expression, revealed the selective capacity of this corepressor to govern common transcriptional events of underlying networks. Combined these findings suggest that NCOR1 is a selective regulator of nuclear receptors, notably PPARgamma and VDR, and contributes to their loss of sensitivity. Combinations of epigenetic therapies that target NCOR1 may prove effective, even when receptor expression is reduced.

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

ChIPping away at gene regulation

The coordinated regulation of gene expression in higher eukaryotes is complex and poorly understood. Recent technological advances have allowed the first insights into these networks on a genome-wide scale. These investigations have identified transcription factor target sites in the genome and successfully predicted cooperative interactions with other factors. However, a detailed understanding of the processes that coordinate gene expression remains elusive. Here, we highlight the advances that have been made using current methods, and the need for new technologies to address the gaps in our knowledge and to map these complex pathways further.

A new look towards BAC-based array CGH through a comprehensive comparison with oligo-based array CGH

BACKGROUND: Currently, two main technologies are used for screening of DNA copy number; the BAC (Bacterial Artificial Chromosome) and the recently developed oligonucleotide-based CGH (Chromosomal Comparative Genomic Hybridization) arrays which are capable of detecting small genomic regions with amplification or deletion. The correlation as well as the discriminative power of these platforms has never been compared statistically on a significant set of human patient samples.

RESULTS: In this paper, we present an exhaustive comparison between the two CGH platforms, undertaken at two independent sites using the same batch of DNA from 19 advanced prostate cancers. The comparison was performed directly on the raw data and a significant correlation was found between the two platforms. The correlation was greatly improved when the data were averaged over large chromosomic regions using a segmentation algorithm. In addition, this analysis has enabled the development of a statistical model to discriminate BAC outliers that might indicate microevents. These microevents were validated by the oligo platform results.

CONCLUSION: This article presents a genome-wide statistical validation of the oligo array platform on a large set of patient samples and demonstrates statistically its superiority over the BAC platform for the Identification of chromosomic events. Taking advantage of a large set of human samples treated by the two technologies, a statistical model has been developed to show that the BAC platform could also detect microevents.

Meta-analysis of Positive and Negative Symptoms Reveals Schizophrenia Modifier Genes

BACKGROUND: Evidence suggests that genetic factors may influence both schizophrenia (Scz) and its clinical presentation. In recent years, genome-wide association studies (GWAS) have demonstrated considerable success in identifying risk loci. Detection of "modifier loci" has the potential to further elucidate underlying disease processes.

METHODS: We performed GWAS of empirically derived positive and negative symptom scales in Irish cases from multiply affected pedigrees and a larger, independent case-control sample, subsequently combining these into a large Irish meta-analysis. In addition to single-SNP associations, we considered gene-based and pathway analyses to better capture convergent genetic effects, and to facilitate biological interpretation of these findings. Replication and testing of aggregate genetic effects was conducted using an independent European-American sample.

RESULTS: Though no single marker met the genome-wide significance threshold, genes and ontologies/pathways were significantly associated with negative and positive symptoms; notably, NKAIN2 and NRG1, respectively. We observed limited overlap in ontologies/pathways associated with different symptom profiles, with immune-related categories over-represented for negative symptoms, and addiction-related categories for positive symptoms. Replication analyses suggested that genes associated with clinical presentation are generalizable to non-Irish samples.

CONCLUSIONS: These findings strongly support the hypothesis that modifier loci contribute to the etiology of distinct Scz symptom profiles. The finding that previously implicated "risk loci" actually influence particular symptom dimensions has the potential to better delineate the roles of these genes in Scz etiology. Furthermore, the over-representation of distinct gene ontologies/pathways across symptom profiles suggests that the clinical heterogeneity of Scz is due in part to complex and diverse genetic factors.

Choline Kinase Alpha as an Androgen Receptor Chaperone and Prostate Cancer Therapeutic Target

BACKGROUND: The androgen receptor (AR) is a major drug target in prostate cancer (PCa). We profiled the AR-regulated kinome to identify clinically relevant and druggable effectors of AR signaling.

METHODS: Using genome-wide approaches, we interrogated all AR regulated kinases. Among these, choline kinase alpha (CHKA) expression was evaluated in benign (n = 195), prostatic intraepithelial neoplasia (PIN) (n = 153) and prostate cancer (PCa) lesions (n = 359). We interrogated how CHKA regulates AR signaling using biochemical assays and investigated androgen regulation of CHKA expression in men with PCa, both untreated (n = 20) and treated with an androgen biosynthesis inhibitor degarelix (n = 27). We studied the effect of CHKA inhibition on the PCa transcriptome using RNA sequencing and tested the effect of CHKA inhibition on cell growth, clonogenic survival and invasion. Tumor xenografts (n = 6 per group) were generated in mice using genetically engineered prostate cancer cells with inducible CHKA knockdown. Data were analyzed with χ(2) tests, Cox regression analysis, and Kaplan-Meier methods. All statistical tests were two-sided.

RESULTS: CHKA expression was shown to be androgen regulated in cell lines, xenografts, and human tissue (log fold change from 6.75 to 6.59, P = .002) and was positively associated with tumor stage. CHKA binds directly to the ligand-binding domain (LBD) of AR, enhancing its stability. As such, CHKA is the first kinase identified as an AR chaperone. Inhibition of CHKA repressed the AR transcriptional program including pathways enriched for regulation of protein folding, decreased AR protein levels, and inhibited the growth of PCa cell lines, human PCa explants, and tumor xenografts.

CONCLUSIONS: CHKA can act as an AR chaperone, providing, to our knowledge, the first evidence for kinases as molecular chaperones, making CHKA both a marker of tumor progression and a potential therapeutic target for PCa.

Genetic Sharing with Cardiovascular Disease Risk Factors and Diabetes Reveals Novel Bone Mineral Density Loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

Pleiotropic analysis of cancer risk loci on esophageal adenocarcinoma risk

Background: Several cancer-associated loci identified from genome-wide association studies (GWAS) have been associated with risks of multiple cancer sites, suggesting pleiotropic effects. We investigated whether GWAS-identified risk variants for other common cancers are associated with risk of esophageal adenocarcinoma (EA) or its precursor, Barrett's esophagus. 

Methods: We examined the associations between risks of EA and Barrett's esophagus and 387 SNPs that have been associated with risks of other cancers, by using genotype imputation data on 2,163 control participants and 3,885 (1,501 EA and 2,384 Barrett's esophagus) case patients from the Barrett's and Esophageal Adenocarcinoma Genetic Susceptibility Study, and investigated effect modification by smoking history, body mass index (BMI), and reflux/heartburn. 

Results: After correcting for multiple testing, none of the tested 387 SNPs were statistically significantly associated with risk of EA or Barrett's esophagus. No evidence of effect modification by smoking, BMI, or reflux/heartburn was observed. 

Conclusions: Genetic risk variants for common cancers identified from GWAS appear not to be associated with risks of EA or Barrett's esophagus. 

Impact: To our knowledge, this is the first investigation of pleiotropic genetic associations with risks of EA and Barrett's esophagus.

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F(CT) = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F(CT) range 0.275-0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.

A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Detection of quantitative trait loci for androstenone, skatole and boar taint in a cross between Large White and Meishan pigs

'Boar taint' is a strong perspiration-like, urine-like unpleasant odour given off upon heating or cooking of meat from some intact (uncastrated) male pigs. Data from the F(2) generation of a Large White (LW) x Meishan (MS) crossbred population were analysed to detect quantitative trait loci (QTL) for traits associated with boar taint. Fat samples from 178 intact male pigs slaughtered at 85 +/- 5 kg were analysed for the major contributors to boar taint (androstenone, indole and skatole). Fat and lean samples from cooked meat were scored for boar, abnormal and pork flavour and odour by a trained sensory panel (SP). A scan with 117 markers covering the whole genome was performed in the F(2) individuals, together with their F(1) parents and purebred grandparents. At the 5% chromosomal significance threshold (approximately equal to the genome-wide suggestive significance threshold), QTL were detected for the laboratory estimate of androstenone on chromosomes 2, 4, 6, 7 and 9. However, only on chromosome 6 were there QTL for boar flavour (BF) traits in the same or adjacent marker intervals as a QTL for the laboratory estimate of androstenone. On chromosome 14, QTL were detected for the laboratory estimates of indole and skatole, the SP score for skatole and the scores for BF in lean and BF in fat. In all five cases, the MS allele generally increased the estimate or score, compared with the LW allele, but it appeared that desirable and undesirable alleles were present in both breeds. This locus on chromosome 14 has considerable potential for use to reduce the incidence of boar taint, especially if further research can identify the causative polymorphism or strongly associated markers.