HES6 drives a critical AR transcriptional programme to induce castration-resistant prostate cancer through activation of an E2F1-mediated cell cycle network

Castrate-resistant prostate cancer (CRPC) is poorly characterized and heterogeneous and while the androgen receptor (AR) is of singular importance, other factors such as c-Myc and the E2F family also play a role in later stage disease. HES6 is a transcription co-factor associated with stem cell characteristics in neural tissue. Here we show that HES6 is up-regulated in aggressive human prostate cancer and drives castration-resistant tumour growth in the absence of ligand binding by enhancing the transcriptional activity of the AR, which is preferentially directed to a regulatory network enriched for transcription factors such as E2F1. In the clinical setting, we have uncovered a HES6-associated signature that predicts poor outcome in prostate cancer, which can be pharmacologically targeted by inhibition of PLK1 with restoration of sensitivity to castration. We have therefore shown for the first time the critical role of HES6 in the development of CRPC and identified its potential in patient-specific therapeutic strategies.

The mitochondrial and autosomal mutation landscapes of prostate cancer

BACKGROUND: Prostate cancer (PCa) is the most common cancer in men. PCa is strongly age associated; low death rates in surveillance cohorts call into question the widespread use of surgery, which leads to overtreatment and a reduction in quality of life. There is a great need to increase the understanding of tumor characteristics in the context of disease progression.

OBJECTIVE: To perform the first multigenome investigation of PCa through analysis of both autosomal and mitochondrial DNA, and to integrate exome sequencing data, and RNA sequencing and copy-number alteration (CNA) data to investigate how various different tumor characteristics, commonly analyzed separately, are interconnected.

DESIGN, SETTING, AND PARTICIPANTS: Exome sequencing was applied to 64 tumor samples from 55 PCa patients with varying stage and grade. Integrated analysis was performed on a core set of 50 tumors from which exome sequencing, CNA, and RNA sequencing data were available.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Genes, mutated at a significantly higher rate relative to a genomic background, were identified. In addition, mitochondrial and autosomal mutation rates were correlated to CNAs and proliferation, assessed as a cell cycle gene expression signature.

RESULTS AND LIMITATIONS: Genes not previously reported to be significantly mutated in PCa, such as cell division cycle 27 homolog (Saccharomyces cerevisiae) (CDC27), myeloid/lymphoid or mixed-lineage leukemia 3 (MLL3), lysine (K)-specific demethylase 6A (KDM6A), and kinesin family member 5A (KIF5A) were identified. The mutation rate in the mitochondrial genome was 55 times higher than that of the autosomes. Multilevel analysis demonstrated a tight correlation between high reactive-oxygen exposure, chromosomal damage, high proliferation, and in parallel, a transition from multiclonal indolent primary PCa to monoclonal aggressive disease. As we only performed targeted sequence analysis; copy-number neutral rearrangements recently described for PCa were not accounted for.

CONCLUSIONS: The mitochondrial genome displays an elevated mutation rate compared to the autosomal chromosomes. By integrated analysis, we demonstrated that different tumor characteristics are interconnected, providing an increased understanding of PCa etiology.

ER stress-mediated autophagy promotes Myc-dependent transformation and tumor growth

The proto-oncogene c-Myc paradoxically activates both proliferation and apoptosis. In the pathogenic state, c-Myc-induced apoptosis is bypassed via a critical, yet poorly understood escape mechanism that promotes cellular transformation and tumorigenesis. The accumulation of unfolded proteins in the ER initiates a cellular stress program termed the unfolded protein response (UPR) to support cell survival. Analysis of spontaneous mouse and human lymphomas demonstrated significantly higher levels of UPR activation compared with normal tissues. Using multiple genetic models, we demonstrated that c-Myc and N-Myc activated the PERK/eIF2α/ATF4 arm of the UPR, leading to increased cell survival via the induction of cytoprotective autophagy. Inhibition of PERK significantly reduced Myc-induced autophagy, colony formation, and tumor formation. Moreover, pharmacologic or genetic inhibition of autophagy resulted in increased Myc-dependent apoptosis. Mechanistically, we demonstrated an important link between Myc-dependent increases in protein synthesis and UPR activation. Specifically, by employing a mouse minute (L24+/-) mutant, which resulted in wild-type levels of protein synthesis and attenuation of Myc-induced lymphomagenesis, we showed that Myc-induced UPR activation was reversed. Our findings establish a role for UPR as an enhancer of c-Myc-induced transformation and suggest that UPR inhibition may be particularly effective against malignancies characterized by c-Myc overexpression.

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.

Global identification of androgen response elements

Chromatin immunoprecipitation (ChIP) is an invaluable tool in the study of transcriptional regulation. ChIP methods require both a priori knowledge of the transcriptional regulators which are important for a given biological system and high-quality specific antibodies for these targets. The androgen receptor (AR) is known to play essential roles in male sexual development, in prostate cancer and in the function of many other AR-expressing cell types (e.g. neurons and myocytes). As a ligand-activated transcription factor the AR also represents an endogenous, inducible system to study transcriptional biology. Therefore, ChIP studies of the AR can make use of treatment contrast experiments to define its transcriptional targets. To date several studies have mapped AR binding sites using ChIP in combination with genome tiling microarrays (ChIP-chip) or direct sequencing (ChIP-seq), mainly in prostate cancer cell lines and with varying degrees of genomic coverage. These studies have provided new insights into the DNA sequences to which the AR can bind, identified AR cooperating transcription factors, mapped thousands of potential AR regulated genes and provided insights into the biological processes regulated by the AR. However, further ChIP studies will be required to fully characterise the dynamics of the AR-regulated transcriptional programme, to map the occupancy of different AR transcriptional complexes which result in different transcriptional output and to delineate the transcriptional networks downstream of the AR.

Terminal and progenitor lineage-survival oncogenes as cancer markers

Tumour classification has traditionally focused on differentiation and cellular morphology, and latterly on the application of genomic approaches. By combining chromatin immunoprecipitation with expression array, it has been possible to identify direct gene targets for transcription factors for nuclear hormone receptors. At the same time, there have been great strides in deriving stem and progenitor cells from tissues. It is therefore timely to propose that pairing the isolation of these cell subpopulations from tissues and tumours with these genomics approaches will reveal conserved gene targets for transcription factors. By focusing on transcription factors (lineage-survival oncogenes) with roles in both organogenesis and tumourigenesis at multiple organ sites, we suggest that this comparative genomics approach will enable developmental biology to be used more fully in relation to understanding tumour progression and will reveal new cancer markers. We focus here on neurogenesis and neuroendocrine differentiation in tumours.

New androgen receptor genomic targets show an interaction with the ETS1 transcription factor

The androgen receptor (AR) initiates important developmental and oncogenic transcriptional pathways. The AR is known to bind as a homodimer to 15-base pair bipartite palindromic androgen-response elements; however, few direct AR gene targets are known. To identify AR promoter targets, we used chromatin immunoprecipitation with on-chip detection of genomic fragments. We identified 1,532 potential AR-binding sites, including previously known AR gene targets. Many of the new AR target genes show altered expression in prostate cancer. Analysis of sequences underlying AR-binding sites showed that more than 50% of AR-binding sites did not contain the established 15 bp AR-binding element. Unbiased sequence analysis showed 6-bp motifs, which were significantly enriched and were bound directly by the AR in vitro. Binding sequences for the avian erythroblastosis virus E26 homologue (ETS) transcription factor family were also highly enriched, and we uncovered an interaction between the AR and ETS1 at a subset of AR promoter targets.

A role for neurotensin in bicalutamide resistant prostate cancer cells

BACKGROUND: Anti-androgens are administered as a principal treatment for prostate cancer. Aggressive hormone refractory disease is characterized in some cases by the development of a neuroendocrine phenotype. However little attention has been paid to resistance pathways selected for by long-term treatment with non-steroidal anti-androgens.

METHODS: Using a resistant sub-line, LNCaP-Bic, we performed a comparative gene expression profiling using cDNA microarrays and target validation by qRT-PCR. Targets were then explored using cell proliferation, cell cycle analysis and in vitro invasion assays using siRNA technology.

RESULTS: Neurotensin/Neuromedin N (NTS) was upregulated in the LNCaP-Bic line at both the transcript and protein level. The resistant line was found to have an increased proliferation rate, more rapid cell cycle progression and increased invasiveness through Matrigel. Each phenotypic difference could be reduced using siRNA knockdown of NT.

CONCLUSION: Increased expression of NT in bicalutamide resistant prostate cancer cells induces cell proliferation and invasion suggesting that this peptide may contribute to the development of bicalutamide resistant prostate cancer.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

Identification of novel Factor XI mutations in a patient requiring dental extractions

Factor XI is a serine protease that participates in the intrinsic pathway of blood coagulation. Patients deficient in factor XI exhibit varying degrees of post operative bleeding following invasive surgical procedures such as dental extractions. Objectives: The aim of the study was to identify the specific mutations in a patient from a family with known factor XI deficiency. Methods: Samples were obtained from the patient, his mother and his father and subjected to DNA sequencing. Each protein coding exon 2-15 of the factor XI gene was amplified by polymerase chain reaction (PCR) followed by bidirectional sequencing utilizing di-deoxy chain termination chemistry. Results: The patient had a factor XI level of 20% of normal. Initial sequencing of factor XI from the patient identified a point mutation (646G>A) and a putative splice site mutation (1567+4A>T) in intron 13. These are novel previously unreported mutations. DNA sequence analysis of the mother revealed the 1567+4A>T mutation and the father exhibited the 646G>A mutation. As a consequence the treatment proceeded without serious bleeding complication and required administration only of transexamic acid though factor XI was available as haemostatic cover. Conclusion: The two mutations identified in this family are novel; further laboratory investigation of the functional consequences of those mutations is currently underway. Although factor XI deficiency is rare in the Northern Irish population this study highlights the techniques available to sequence and analyse this and similar haematological disorders.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Environmental selection on transcriptome-derived SNPs in a high gene flow marine fish, the Atlantic herring (Clupea harengus)

High gene flow is considered the norm for most marine organisms and is expected to limit their ability to adapt to local environments. Few studies have directly compared the patterns of differentiation at neutral and selected gene loci in marine organisms. We analysed a transcriptome-derived panel of 281 SNPs in Atlantic herring (Clupea harengus), a highly migratory small pelagic fish, for elucidating neutral and selected genetic variation among populations and to identify candidate genes for environmental adaptation. We analysed 607 individuals from 18 spawning locations in the northeast Atlantic, including two temperature clines (5-12 °C) and two salinity clines (5-35‰). By combining genome scan and landscape genetic analyses, four genetically distinct groups of herring were identified: Baltic Sea, Baltic-North Sea transition area, North Sea/British Isles and North Atlantic; notably, samples exhibited divergent clustering patterns for neutral and selected loci. We found statistically strong evidence for divergent selection at 16 outlier loci on a global scale, and significant correlations with temperature and salinity at nine loci. On regional scales, we identified two outlier loci with parallel patterns across temperature clines and five loci associated with temperature in the North Sea/North Atlantic. Likewise, we found seven replicated outliers, of which five were significantly associated with low salinity across both salinity clines. Our results reveal a complex pattern of varying spatial genetic variation among outlier loci, likely reflecting adaptations to local environments. In addition to disclosing the fine scale of local adaptation in a highly vagile species, our data emphasize the need to preserve functionally important biodiversity.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Production of extended-spectrum β-lactamases and the potential indirect pathogenic role of Prevotella isolates from the cystic fibrosis respiratory microbiota

Extended-spectrum β-lactamase (ESBL) production and the prevalence of the β-lactamase-encoding gene blaTEM were determined in Prevotella isolates (n=50) cultured from the respiratory tract of adults and young people with cystic fibrosis (CF). Time-kill studies were used to investigate the concept of passive antibiotic resistance and to ascertain whether a β-lactamase-positive Prevotella isolate can protect a recognised CF pathogen from the action of ceftazidime in vitro. The results indicated that approximately three-quarters (38/50; 76%) of Prevotella isolates produced ESBLs. Isolates positive for ESBL production had higher minimum inhibitory concentrations (MICs) of β-lactam antibiotics compared with isolates negative for production of ESBLs (P<0.001). The blaTEM gene was detected more frequently in CF Prevotella isolates from paediatric patients compared with isolates from adults (P=0.002), with sequence analysis demonstrating that 21/22 (95%) partial blaTEM genes detected were identical to blaTEM-116. Furthermore, a β-lactamase-positive Prevotella isolate protected Pseudomonas aeruginosa from the antimicrobial effects of ceftazidime (P=0.03). Prevotella isolated from the CF respiratory microbiota produce ESBLs and may influence the pathogenesis of chronic lung infection via indirect methods, including shielding recognised pathogens from the action of ceftazidime.