Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Oxidised, glycated LDL selectively influences tissue inhibitor of metalloproteinase-3 gene expression and protein production in human retinal capillary pericytes

Aims/hypothesis: Matrix metalloproteinases (MMPs) and their natural inhibitors, tissue inhibitor of metalloproteinases (TIMPs), regulate important biological processes including the homeostasis of the extracellular matrix, proteolysis of cell surface proteins, proteinase zymogen activation, angiogenesis and inflammation. Studies have shown that their balance is altered in retinal microvascular tissues in diabetes. Since LDLs modified by oxidation/glycation are implicated in the pathogenesis of diabetic vascular complications, we examined the effects of modified LDL on the gene expression and protein production of MMPs and TIMPs in retinal pericytes. Methods: Quiescent human retinal pericytes were exposed to native LDL (N-LDL), glycated LDL (G-LDL) and heavily oxidised and glycated LDL (HOG-LDL) for 24 h. We studied the expression of the genes encoding MMPs and TIMPs mRNAs by analysis of microarray data and quantitative PCR, and protein levels by immunoblotting and ELISA. Results: Microarray analysis showed that MMP1, MMP2, MMP11, MMP14 and MMP25 and TIMP1, TIMP2, TIMP3 and TIMP4 were expressed in pericytes. Of these, only TIMP3 mRNA showed altered regulation, being expressed at significantly lower levels in response to HOG- vs N-LDL. Quantitative PCR and immunoblotting of cell/matrix proteins confirmed the reduction in TIMP3 mRNA and protein in response to HOG-LDL. In contrast to cellular TIMP3 protein, analysis of secreted TIMP1, TIMP2, MMP1 and collagenase activity indicated no changes in their production in response to modified LDL. Combined treatment with N- and HOG-LDL restored TIMP3 mRNA expression to a level comparable with that after N-LDL alone. Conclusions/interpretation: Among the genes encoding for MMPs and TIMPs expressed in retinal pericytes, TIMP3 is uniquely regulated by HOG-LDL. Reduced TIMP3 expression might contribute to microvascular abnormalities in diabetic retinopathy. © 2007 Springer-Verlag.

Fosfomycin and tobramycin in combination downregulate nitrate reductase genes narG and narH, resulting in increased activity against Pseudomonas aeruginosa under anaerobic conditions

The activity of aminoglycosides, used to treat Pseudomonas aeruginosa respiratory infection in cystic fibrosis (CF) patients, is reduced under the anaerobic conditions that reflect the CF lung in vivo. In contrast, a 4:1 (w/w) combination of fosfomycin and tobramycin (F:T), under investigation for use in the treatment of CF lung infection, has increased activity against P. aeruginosa under anaerobic conditions. The aim of this study was to elucidate the mechanisms underlying the increased activity of F:T under anaerobic conditions. Microarray analysis was used to identify the transcriptional basis of increased F:T activity under anaerobic conditions, and key findings were confirmed by microbiological tests including nitrate utilization assays, growth curves and susceptibility testing. Notably, growth in sub-inhibitory concentrations of F:T, but not tobramycin or fosfomycin alone, significantly downregulated (p <0.05) nitrate reductase genes narG and narH, essential for normal anaerobic growth of P. aeruginosa. Under anaerobic conditions, F:T significantly decreased (p

Regulator of G-protein signaling 4 (RGS4) gene is associated with schizophrenia in Irish high density families

The regulator of the G-protein signaling 4 (RGS4) gene was shown to have a different expression pattern in schizophrenia patients in a microarray study. A family-based study subsequently implicated the association of this gene with schizophrenia. We replicated the study with our sample from the Irish Study of High Density Schizophrenia Families (ISHDSF). Single marker transmission disequilibrium tests (TDT) for the four core SNPs showed modest association for SNP 18 (using a narrow diagnostic approach with FBAT P = 0.044; with PDT P = 0.0073) and a trend for SNP 4 (with FBAT P = 0.1098; with PDT P = 0.0249). For SNP 1 and 7, alleles overtransmitted to affected subjects were the same as previously reported. Haplotype analyses suggested that haplotype G-G-G for SNP1-4-18, which is the most abundant haplotype (42.3%) in the Irish families, was associated with the disease (narrow diagnosis, FBAT P = 0.0061, PDT P = 0.0498). This was the same haplotype implicated in the original study. While P values were not corrected for multiple testing because of the clear prior hypothesis, these results could be interpreted as supporting evidence for the association between RGS4 and schizophrenia.

Novel hydrolysis-probe based qPCR assay to detect saxitoxin transcripts of dinoflagellates in environmental samples

Paralytic Shellfish Poisoning (PSP) is a serious human illness caused by ingestion of seafood enriched with paralytic shellfish toxins (PSTs). PSTs are neurotoxic compounds produced by marine dinoflagellates, specifically by Alexandrium spp., Gymnodinium catenatum and Pyrodinium bahamense. Every year, massive monitoring of PSTs and their producers is undertaken worldwide to avoid PSP incidences. Here we developed a sensitive, hydrolysis probe-based quantitative PCR (qPCR) assay to detect a gene essential for PST synthesis across different dinoflagellate species and genera and tested it on cDNA generated from environmental samples spiked with Alexandrium minutum or Alexandrium fundyense cells. The assay was then applied to two environmental sample series from Norway and Spain and the results were complemented with cell counts, LSU-based microarray data and toxin measurements (enzyme-linked immunosorbent assay (ELISA) and surface plasmon resonance (SPR) biosensor method). The overall agreement between the results of the qPCR assay and the complementary data was good. The assay reliably detected sxtA transcripts from Alexandrium spp. and G. catenatum, even though Alexandrium spp. cell concentrations were mostly so low that they could not be quantified microscopically. Agreement between the novel assay and toxin measurements or cell counts was generally good; the few inconsistencies observed were most likely due to disparate residence times of sxtA transcripts and PSTs in seawater, or, in the case of cell counts, to dissimilar sxtA4 transcript numbers per cell in different dinoflagellate strains or species. © 2013 Elsevier B.V.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Identification and Validation of an Anthracycline/Cyclophosphamide–Based Chemotherapy Response Assay in Breast Cancer

Background: There is no method routinely used to predict response to anthracycline and cyclophosphamide–based chemotherapy in the clinic; therefore patients often receive treatment for breast cancer with no benefit. Loss of the Fanconi anemia/BRCA (FA/BRCA) DNA damage response (DDR) pathway occurs in approximately 25% of breast cancer patients through several mechanisms and results in sensitization to DNA-damaging agents. The aim of this study was to develop an assay to detect DDR-deficient tumors associated with loss of the FA/BRCA pathway, for the purpose of treatment selection.

Methods: DNA microarray data from 21 FA patients and 11 control subjects were analyzed to identify genetic processes associated with a deficiency in DDR. Unsupervised hierarchical clustering was then performed using 60 BRCA1/2 mutant and 47 sporadic tumor samples, and a molecular subgroup was identified that was defined by the molecular processes represented within FA patients. A 44-gene microarray-based assay (the DDR deficiency assay) was developed to prospectively identify this subgroup from formalin-fixed, paraffin-embedded samples. All statistical tests were two-sided.

Results: In a publicly available independent cohort of 203 patients, the assay predicted complete pathologic response vs residual disease after neoadjuvant DNA-damaging chemotherapy (5-fluorouracil, anthracycline, and cyclophosphamide) with an odds ratio of 3.96 (95% confidence interval [Cl] =1.67 to 9.41; P = .002). In a new independent cohort of 191 breast cancer patients treated with adjuvant 5-fluorouracil, epirubicin, and cyclophosphamide, a positive assay result predicted 5-year relapse-free survival with a hazard ratio of 0.37 (95% Cl = 0.15 to 0.88; P = .03) compared with the assay negative population.

Conclusions: A formalin-fixed, paraffin-embedded tissue-based assay has been developed and independently validated as a predictor of response and prognosis after anthracycline/cyclophosphamide–based chemotherapy in the neoadjuvant and adjuvant settings. These findings warrant further validation in a prospective clinical study.

Macrophage migration inhibitory factor (MIF) expression in human malignant gliomas contributes to immune escape and tumour progression

Macrophage migration inhibitory factor (MIF), which inhibits apoptosis and promotes angiogenesis, is expressed in cancers suppressing immune surveillance. Its biological role in human glioblastoma is, however, only poorly understood. We examined in-vivo expression of MIF in 166 gliomas and 23 normal control brains by immunohistochemistry. MIF immunoreactivity was enhanced in neoplastic astrocytes in WHO grade II glioma and increased significantly in higher tumour grades (III-IV). MIF expression was further assessed in 12 glioma cell lines in vitro. Quantitative RT-PCR showed that MIF mRNA expression was elevated up to 800-fold in malignant glioma cells compared with normal brain. This translated into high protein levels as assessed by immunoblotting of total cell lysates and by ELISA-based measurement of secreted MIF. Wild-type p53-retaining glioma cell lines expressed higher levels of MIF, which may be connected with the previously described role of MIF as a negative regulator of wild-type p53 signalling in tumour cells. Stable knockdown of MIF by shRNA in glioma cells significantly increased tumour cell susceptibility towards NK cell-mediated cytotoxicity. Furthermore, supernatant from mock-transfected cells, but not from MIF knockdown cells, induced downregulation of the activating immune receptor NKG2D on NK and CD8+ T cells. We thus propose that human glioma cell-derived MIF contributes to the immune escape of malignant gliomas by counteracting NK and cytotoxic T-cell-mediated tumour immune surveillance. Considering its further cell-intrinsic and extrinsic tumour-promoting effects and the availability of small molecule inhibitors, MIF seems to be a promising candidate for future glioma therapy.

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Gene expression in normal Urothelium depends on location within the bladder:A possible link to bladder carcinogenesis

Objectives: Clinical studies have shown that more than 70% of primary bladder tumours arise in the area around the ureteric orifice. In this study a genomic approach was taken to explore the molecular mechanisms that may influence this phenomenon.

Methods: RNA was isolated from each individual normal ureteric orifice and the dome biopsy from 33 male patients. Equal amounts of the pooled ureteric orifice and dome mRNAs were labelled with Cy3 and Cy5, respectively before hybridising to the gene chip (UniGEM 2.0, Incyte Genomics Inc., Wilmington, Delaware, USA). Results: Significant changes (more than a twofold difference) in gene expression were observed in 3.1% (312) of the 10,176 gene array: 211 genes upregulated and 101 downregulated. Analysis of Cdc25B, TK1, PKM, and PDGFra with RT-PCR supported the reliability of the microarray result. Seladin-1 was the most upregulated gene in the ureteric orifice: 8.3-fold on the microarray and 11.4-fold by real time PCR.

Conclusions: Overall, this study suggests significant altered gene expression between these two anatomically distinct areas of the normal human bladder. Of particular note is Seladin-1, whose significance in cancer is yet to be clarified. Further studies of the genes discovered by this work will help clarify which of these differences influence primary bladder carcinogenesis. (c) 2006 European Association of Urology. Published by Elsevier B.V. All rights reserved.

Digital pathology and image analysis in tissue biomarker research

Digital pathology and the adoption of image analysis have grown rapidly in the last few years. This is largely due to the implementation of whole slide scanning, advances in software and computer processing capacity and the increasing importance of tissue-based research for biomarker discovery and stratified medicine. This review sets out the key application areas for digital pathology and image analysis, with a particular focus on research and biomarker discovery. A variety of image analysis applications are reviewed including nuclear morphometry and tissue architecture analysis, but with emphasis on immunohistochemistry and fluorescence analysis of tissue biomarkers. Digital pathology and image analysis have important roles across the drug/companion diagnostic development pipeline including biobanking, molecular pathology, tissue microarray analysis, molecular profiling of tissue and these important developments are reviewed. Underpinning all of these important developments is the need for high quality tissue samples and the impact of pre-analytical variables on tissue research is discussed. This requirement is combined with practical advice on setting up and running a digital pathology laboratory. Finally, we discuss the need to integrate digital image analysis data with epidemiological, clinical and genomic data in order to fully understand the relationship between genotype and phenotype and to drive discovery and the delivery of personalized medicine.

Spindle assembly checkpoint protein expression correlates with cellular proliferation and shorter time to recurrence in ovarian cancer.

Ovarian carcinoma (OC) is the most lethal of the gynecological malignancies, often presenting at an advanced stage. Treatment is hampered by high levels of drug resistance. The taxanes are microtubule stabilizing agents, used as first-line agents in the treatment of OC that exert their apoptotic effects through the spindle assembly checkpoint. BUB1-related protein kinase (BUBR1) and mitotic arrest deficient 2 (MAD2), essential spindle assembly checkpoint components, play a key role in response to taxanes. BUBR1, MAD2, and Ki-67 were assessed on an OC tissue microarray platform representing 72 OC tumors of varying histologic subtypes. Sixty-one of these patients received paclitaxel and platinum agents combined; 11 received platinum alone. Overall survival was available for all 72 patients, whereas recurrence-free survival (RFS) was available for 66 patients. Increased BUBR1 expression was seen in serous carcinomas, compared with other histologies (P = .03). Increased BUBR1 was significantly associated with tumors of advanced stage (P = .05). Increased MAD2 and BUBR1 expression also correlated with increased cellular proliferation (P < .0002 and P = .02, respectively). Reduced MAD2 nuclear intensity was associated with a shorter RFS (P = .03), in ovarian tumors of differing histologic subtype (n = 66). In this subgroup, for those women who received paclitaxel and platinum agents combined (n = 57), reduced MAD2 intensity also identified women with a shorter RFS (P < .007). For the entire cohort of patients, irrespective of histologic subtype or treatment, MAD2 nuclear intensity retained independent significance in a multivariate model, with tumors showing reduced nuclear MAD2 intensity identifying patients with a poorer RFS (P = .05).

Discovering subgroups of patients from DNA copy number data using NMF on compacted matrices

In the study of complex genetic diseases, the identification of subgroups of patients sharing similar genetic characteristics represents a challenging task, for example, to improve treatment decision. One type of genetic lesion, frequently investigated in such disorders, is the change of the DNA copy number (CN) at specific genomic traits. Non-negative Matrix Factorization (NMF) is a standard technique to reduce the dimensionality of a data set and to cluster data samples, while keeping its most relevant information in meaningful components. Thus, it can be used to discover subgroups of patients from CN profiles. It is however computationally impractical for very high dimensional data, such as CN microarray data. Deciding the most suitable number of subgroups is also a challenging problem. The aim of this work is to derive a procedure to compact high dimensional data, in order to improve NMF applicability without compromising the quality of the clustering. This is particularly important for analyzing high-resolution microarray data. Many commonly used quality measures, as well as our own measures, are employed to decide the number of subgroups and to assess the quality of the results. Our measures are based on the idea of identifying robust subgroups, inspired by biologically/clinically relevance instead of simply aiming at well-separated clusters. We evaluate our procedure using four real independent data sets. In these data sets, our method was able to find accurate subgroups with individual molecular and clinical features and outperformed the standard NMF in terms of accuracy in the factorization fitness function. Hence, it can be useful for the discovery of subgroups of patients with similar CN profiles in the study of heterogeneous diseases.

Genomic lesions associated with a different clinical outcome in diffuse large B-Cell lymphoma treated with R-CHOP-21

Despite recent therapeutic improvements, the clinical course of diffuse large B-cell lymphoma (DLBCL) still differs considerably among patients. We conducted this retrospective multi-centre study to evaluate the impact of genomic aberrations detected using a high-density genome wide-single nucleotide polymorphism-based array on clinical outcome in a population of DLBCL patients treated with R-CHOP-21 (rituximab, cyclophosphamide, doxorubicine, vincristine and prednisone repeated every 21_d). 166 DNA samples were analysed using the GeneChip Human Mapping 250K NspI. Genomic anomalies were analysed regarding their impact on the clinical course of 124 patients treated with R-CHOP-21. Unsupervised clustering was performed to identify genetically related subgroups of patients with different clinical outcomes. Twenty recurrent genetic lesions showed an impact on the clinical course. Loss of genomic material at 8p23.1 showed the strongest statistical significance and was associated with additional aberrations, such as 17p- and 15q-. Unsupervised clustering identified five DLBCL clusters with distinct genetic profiles, clinical characteristics and outcomes. Genetic features and clusters, associated with a different outcome in patients treated with R-CHOP, have been identified by arrayCGH.