Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

Genome-wide gene pathway analysis of psychotic illness symptom dimensions based on a new schizophrenia-specific model of the OPCRIT

Empirically derived phenotypic measurements have the potential to enhance gene-finding efforts in schizophrenia. Previous research based on factor analyses of symptoms has typically included schizoaffective cases. Deriving factor loadings from analysis of only narrowly defined schizophrenia cases could yield more sensitive factor scores for gene pathway and gene ontology analyses. Using an Irish family sample, this study 1) factor analyzed clinician-rated Operational Criteria Checklist items in cases with schizophrenia only, 2) scored the full sample based on these factor loadings, and 3) implemented genome-wide association, gene-based, and gene-pathway analysis of these SCZ-based symptom factors (final N= 507). Three factors emerged from the analysis of the schizophrenia cases: a manic, a depressive, and a positive symptom factor. In gene-based analyses of these factors, multiple genes had q<. 0.01. Of particular interest are findings for PTPRG and WBP1L, both of which were previously implicated by the Psychiatric Genomics Consortium study of SCZ; results from this study suggest that variants in these genes might also act as modifiers of SCZ symptoms. Gene pathway analyses of the first factor indicated over-representation of glutamatergic transmission, GABA-A receptor, and cyclic GMP pathways. Results suggest that these pathways may have differential influence on affective symptom presentation in schizophrenia.

Convergent genetic and expression data implicate immunity in Alzheimer's disease

Background: Late-onset Alzheimer's disease (AD) is heritable with 20 genes showing genome-wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease, we extended these genetic data in a pathway analysis.

Methods: The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain.

Results: ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response (P = 3.27 X 10(-12) after multiple testing correction for pathways), regulation of endocytosis (P = 1.31 X 10(-11)), cholesterol transport (P = 2.96 X 10(-9)), and proteasome-ubiquitin activity (P = 1.34 X 10(-6)). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected P = .002-.05).

Conclusions: The immime response, regulation of endocytosis, cholesterol transport, and protein ubiquitination represent prime targets for AD therapeutics. (C) 2015 Published by Elsevier Inc. on behalf of The Alzheimer's Association.

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.

Autosomal dominant retinitis pigmentosa:linkage to rhodopsin and evidence for genetic heterogeneity

Retinitis pigmentosa (RP) is the most prevalent human retinopathy of genetic origin. Chromosomal locations for X-linked RP and autosomal dominant RP genes have recently been established. Multipoint analyses with ADRP and seven markers on the long arm of chromosome 3 demonstrate that the gene for rhodopsin, the pigment of the rod photoreceptors, cosegregates with the disease locus with a maximum lod score of approximately 19, implicating rhodopsin as a causative gene. Recent studies have indicated the presence of a point mutation at codon 23 in exon 1 of rhodopsin which results in the substitution of histidine for the highly conserved amino acid proline, suggesting that this mutation is a cause of rhodopsin-linked ADRP. This mutation is not present in the Irish pedigree in which ADRP has been mapped close to rhodopsin. Another mutation in the rhodopsin gene or in a gene closely linked to rhodopsin may be involved. Moreover, the gene in a second ADRP pedigree, with Type II late onset ADRP, does not segregate with chromosome 3q markers, indicating that nonallelic as well as perhaps allelic genetic heterogeneity exists in the autosomal dominant form of this disease.

Epidermolysis bullosa:evidence for linkage to genetic markers on chromosome 1 in a family with the autosomal dominant simplex form

DNA from members of a three-generation pedigree of Irish origin, displaying an autosomal dominant simplex form of epidermolysis bullosa of the epidermolytic, simplex, or Koebner variety (EBS2), was analyzed for linkage with a set of markers derived from the long arm of chromosome 1. Two-point analysis revealed positive lod scores for five of these markers, AT3 (Z = 2.107, theta = 0), APOA2 (Z = 1.939, theta = 0.15), D1S66 (Z = 1.204, theta = 0), D1S13 (Z = 1.026, theta = 0.15), and D1S65 (Z = 0.329, theta = 0.15). Multilocus analysis, incorporating the markers D1S19, D1S16, D1S13, APOA2, D1S66, AT3, and D1S65, resulted in a lod score of 3 maximizing at AT3. These data strongly support previous tentative indications of linkage between EBS2 and genetic markers on the long arm of chromosome 1.

Autosomal dominant retinitis pigmentosa (ADRP):localization of an ADRP gene to the long arm of chromosome 3

Members of a large pedigree of Irish origin presenting with early onset Type I autosomal dominant retinitis pigmentosa (ADRP) have been typed for D3S47 (C17), a polymorphic marker from the long arm of chromosome 3. Significant, tight linkage of ADRP to D3S47, with a lod score of 14.7 maximizing at 0.00 recombination, has been obtained, hence localizing the ADRP gene (RP1) segregating in this pedigree to 3q.

Report from the 2nd Summer School in Computational Biology organized by the Queen's University of Belfast

In this paper, we present a meeting report for the 2nd Summer School in Computational Biology organized by the Queen's University of Belfast. We describe the organization of the summer school, its underlying concept and student feedback we received after the completion of the summer school.

Living long and ageing well: is epigenomics the missing link between nature and nurture?

Human longevity is a complex trait and increasingly we understand that both genes and lifestyle interact in the longevity phenotype. Non-genetic factors, including diet, physical activity, health habits, and psychosocial factors contribute approximately 50 % of the variability in human lifespan with another 25 % explained by genetic differences. Family clusters of nonagenarian and centenarian siblings, who show both exceptional age-span and health-span, are likely to have inherited facilitatory gene groups, but also have nine decades of life experiences and behaviours which have interacted with their genetic profiles. Identification of their shared genes is just one small step in the link from genes to their physical and psychological profiles. Behavioural genomics is beginning to demonstrate links to biological mechanisms through regulation of gene expression, which directs the proteome and influences the personal phenotype. Epigenetics has been considered the missing link between nature and nurture. Although there is much that remains to be discovered, this article will discuss some of genetic and environmental factors which appear important in good quality longevity and link known epigenetic mechanisms to themes identified by nonagenarians themselves related to their longevity. Here we suggest that exceptional 90-year old siblings have adopted a range of behaviours and life-styles which have contributed to their ageing-well-phenotype and which link with important public health messages.

Meta-analysis of prostate cancer gene expression data identifies a novel discriminatory signature enriched for glycosylating enzymes

BACKGROUND: Tumorigenesis is characterised by changes in transcriptional control. Extensive transcript expression data have been acquired over the last decade and used to classify prostate cancers. Prostate cancer is, however, a heterogeneous multifocal cancer and this poses challenges in identifying robust transcript biomarkers.

METHODS: In this study, we have undertaken a meta-analysis of publicly available transcriptomic data spanning datasets and technologies from the last decade and encompassing laser capture microdissected and macrodissected sample sets.

RESULTS: We identified a 33 gene signature that can discriminate between benign tissue controls and localised prostate cancers irrespective of detection platform or dissection status. These genes were significantly overexpressed in localised prostate cancer versus benign tissue in at least three datasets within the Oncomine Compendium of Expression Array Data. In addition, they were also overexpressed in a recent exon-array dataset as well a prostate cancer RNA-seq dataset generated as part of the The Cancer Genomics Atlas (TCGA) initiative. Biologically, glycosylation was the single enriched process associated with this 33 gene signature, encompassing four glycosylating enzymes. We went on to evaluate the performance of this signature against three individual markers of prostate cancer, v-ets avian erythroblastosis virus E26 oncogene homolog (ERG) expression, prostate specific antigen (PSA) expression and androgen receptor (AR) expression in an additional independent dataset. Our signature had greater discriminatory power than these markers both for localised cancer and metastatic disease relative to benign tissue, or in the case of metastasis, also localised prostate cancer.

CONCLUSION: In conclusion, robust transcript biomarkers are present within datasets assembled over many years and cohorts and our study provides both examples and a strategy for refining and comparing datasets to obtain additional markers as more data are generated.

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.

Taking risks with translational research

The recent bankruptcy filing by deCODE, a company with an exceptional pedigree in associating genetic variance with disease onset, highlights the commercial risks of translational research. Indeed, deCODE's approach was similar to that adapted by academic researchers who seek to connect genetics and disease. We argue here that neither a purely corporate nor purely academic model is entirely appropriate for such research. Instead, we suggest that the private sector undertake the high-throughput elements of translational research, while the public sector and governments assume the role of providing long-term funding to develop gifted scientists with the confidence to attempt to use genetic data as a stepping stone to a truly mechanistic understanding of complex disease.

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.

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.

Pro-neural transcription factors as cancer markers

BACKGROUND: The aberrant transcription in cancer of genes normally associated with embryonic tissue differentiation at various organ sites may be a hallmark of tumour progression. For example, neuroendocrine differentiation is found more commonly in cancers destined to progress, including prostate and lung. We sought to identify proteins which are involved in neuroendocrine differentiation and differentially expressed in aggressive/metastatic tumours.

RESULTS: Expression arrays were used to identify up-regulated transcripts in a neuroendocrine (NE) transgenic mouse model of prostate cancer. Amongst these were several genes normally expressed in neural tissues, including the pro-neural transcription factors Ascl1 and Hes6. Using quantitative RT-PCR and immuno-histochemistry we showed that these same genes were highly expressed in castrate resistant, metastatic LNCaP cell-lines. Finally we performed a meta-analysis on expression array datasets from human clinical material. The expression of these pro-neural transcripts effectively segregates metastatic from localised prostate cancer and benign tissue as well as sub-clustering a variety of other human cancers.

CONCLUSION: By focussing on transcription factors known to drive normal tissue development and comparing expression signatures for normal and malignant mouse tissues we have identified two transcription factors, Ascl1 and Hes6, which appear effective markers for an aggressive phenotype in all prostate models and tissues examined. We suggest that the aberrant initiation of differentiation programs may confer a selective advantage on cells in all contexts and this approach to identify biomarkers therefore has the potential to uncover proteins equally applicable to pre-clinical and clinical cancer biology.