Proteomics and in silico approaches to extend understanding of the glutathione transferase superfamily of the tropical liver fluke Fasciola gigantica

Fasciolosis is an important foodborne, zoonotic disease of livestock and humans, with global annual health and economic losses estimated at several billion US$. Fasciola hepatica is the major species in temperate regions, while F. gigantica dominates in the tropics. In the absence of commercially available vaccines to control fasciolosis, increasing reports of resistance to current chemotherapeutic strategies and the spread of fasciolosis into new areas, new functional genomics approaches are being used to identify potential new drug targets and vaccine candidates. The glutathione transferase (GST) superfamily is both a candidate drug and vaccine target. This study reports the identification of a putatively novel Sigma class GST, present in a water-soluble cytosol extract from the tropical liver fluke F. gigantica. The GST was cloned and expressed as an enzymically active recombinant protein. This GST shares a greater identity with the human schistosomiasis GST vaccine currently at Phase II clinical trials than previously discovered F. gigantica GSTs, stimulating interest in its immuno-protective properties. In addition, in silico analysis of the GST superfamily of both F. gigantica and F. hepatica has revealed an additional Mu class GST, Omega class GSTs, and for the first time, a Zeta class member.

Expression of the NF-kappaB targets BCL2 and BIRC5/Survivin characterizes small B-cell and aggressive B-cell lymphomas, respectively

Nuclear factor kappa B (NF-kappaB) activation has been proposed as a cardinal feature of tumourigenesis, although the precise mechanism, frequency, relevance, and extent of NF-kappaB activation in lymphomas remain to be fully elucidated. In this study, expression profiling and tissue microarray studies of 209 and 323 non-Hodgkin's lymphomas (NHLs) respectively, including the most frequent sub-types of NHL, were employed to generate a hypothesis concerning the most common NF-kappaB targets in NHL. These analyses showed that NF-kappaB activation is a common phenomenon in NHL, resulting in the expression of distinct sets of NF-kappaB target genes, depending on the cell context. BCL2 and BIRC5/Survivin were identified as key NF-kappaB targets and their expression distinguished small and aggressive B-cell lymphomas, respectively. Interestingly, in the vast majority of B-cell lymphomas, the expression of these markers was mutually exclusive. A set of genes was identified whose expression correlates either with BIRC5/Survivin or with BCL2. BIRC5/Survivin expression, in contrast to BCL2, was associated with a signature of cell proliferation (overexpression of cell cycle control, DNA repair, and polymerase genes), which may contribute to the aggressive phenotype and poor prognosis of these lymphomas. Strikingly, mantle cell lymphoma and chronic lymphocytic leukaemia expressed highly elevated levels of BCL2 protein and mRNA, higher than that observed in reactive mantle zone cells or even in follicular lymphomas, where BCL2 expression is deregulated through the t(14;18) translocation. In parallel with this observation, BIRC5/Survivin expression was higher in Burkitt's lymphoma and diffuse large B-cell lymphoma than in non-tumoural germinal centre cells. In vitro studies confirmed that NF-kappaB activation contributes to the expression of both markers. In cell lines representing aggressive lymphomas, NF-kappaB inhibition resulted in a decrease in BIRC5/Survivin expression. Meanwhile, in chronic lymphocytic leukaemia (CLL)-derived lymphocytes, NF-kappaB inhibition resulted in a marked decrease in BCL2 expression.

Transcriptional response of T cells to IFN-alpha:changes induced in IFN-alpha-sensitive and resistant cutaneous T cell lymphoma

Interferon-alpha (IFN-alpha) therapy is commonly used in the treatment of neoplastic and autoimmune diseases, including cutaneous T cell lymphoma (CTCL). However, the IFN-alpha response is unpredictable, and the IFN-alpha cell targets and pathways are only partially understood. To delineate the molecular mechanisms of IFN-alpha activity, gene expression profiling was performed in a time-course experiment of both IFN-alpha sensitive and IFN-alpha-resistant variants of a CTCL cell line. These experiments revealed that IFN-alpha is responsible for the regulation of hundreds of genes in both variants and predominantly involves genes implicated in signal transduction, cell cycle control, apoptosis, and transcription regulation. Specifically, the IFN-alpha response of tumoral T cells is due to a combination of induction of apoptosis in which TNFSF10 and HSXIAPAF1 may play an important role and cell cycle arrest achieved by downregulation of CDK4 and CCNG2 and upregulation of CDKN2C and tumor suppressor genes (TSGs). Resistance to IFN-alpha appears to be associated with failure to induce IRF1 and IRF7 and deregulation of the apoptotic signals of HSXIAPAF1, TRADD, BAD, and BNIP3. Additionally, cell cycle progression is heralded by upregulation of CDC25A and CDC42. A critical role of NF-kappaB in promoting cell survival in IFN-alpha-resistant cells is indicated by the upregulation of RELB and LTB.

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.