Integration of gene expression and GWAS results supports involvement of calcium signaling in Schizophrenia.

The number of Genome Wide Association Studies (GWAS) of schizophrenia is rapidly growing. However, the small effect of individual genes limits the number of reliably implicated genes, which are too few and too diverse to perform reliable pathway analysis; hence the biological roles of the genes implicated in schizophrenia are unclear. To overcome these limitations we combine GWAS with genome-wide expression data from human post-mortem brain samples of schizophrenia patients and controls, taking these steps: 1) Identify 36 GWAS-based genes which are expressed in our dataset. 2) Find a cluster of 19 genes with highly correlated expression. We show that this correlation pattern is robust and statistically significant. 3) GO-enrichment analysis of these 19 genes reveals significant enrichment of ion channels and calcium-related processes. This finding (based on analyzing a small number of coherently expressed genes) is validated and enhanced in two ways: First, the emergence of calcium channels and calcium signaling is corroborated by identifying proteins that interact with those encoded by the cluster of 19. Second, extend the 19 cluster genes into 1028 genes, whose expression is highly correlated with the cluster's average profile. When GO-enrichment analysis is performed on this extended set, many schizophrenia related pathways appear, with calcium-related processes enriched with high statistical significance. Our results give further, expression-based validation to GWAS results, support a central role of calcium-signaling in the pathogenesis of schizophrenia, and point to additional pathways potentially related to the disease.

Reactivation of epigenetically silenced miR-512 and miR-373 sensitizes lung cancer cells to cisplatin and restricts tumor growth.

MicroRNAs (miRs) regulate a variety of cellular processes, and their impaired expression is involved in cancer. Silencing of tumor-suppressive miRs in cancer can occur through epigenetic modifications, including DNA methylation and histone deacetylation. We performed comparative miR profiling on cultured lung cancer cells before and after treatment with 5'aza-deoxycytidine plus Trichostatin A to reverse DNA methylation and histone deacetylation, respectively. Several tens of miRs were strongly induced by such 'epigenetic therapy'. Two representatives, miR-512-5p (miR-512) and miR-373, were selected for further analysis. Both miRs were secreted in exosomes. Re-expression of both miRs augmented cisplatin-induced apoptosis and inhibited cell migration; miR-512 also reduced cell proliferation. TEAD4 mRNA was confirmed as a direct target of miR-512; likewise, miR-373 was found to target RelA and PIK3CA mRNA directly. Our results imply that miR-512 and miR-373 exert cell-autonomous and non-autonomous tumor-suppressive effects in lung cancer cells, where their re-expression may benefit epigenetic cancer therapy.

Expression of TP53 mutation-associated microRNAs predicts clinical outcome in head and neck squamous cell carcinoma patients.

BACKGROUND: TP53 mutation is associated with decreased survival rate in head and neck squamous cell carcinoma (HNSCC) patients. We set out to identify microRNAs (miRNAs) whose expression associates with TP53 mutation and survival in HNSCC. PATIENTS AND METHODS: We analyzed TP53 status by direct sequencing of exons 2 through 11 of a prospective series of 121 HNSCC samples and assessed its association with outcome in 109 followed-up patients. We carried out miRNA expression profiling on 121 HNSCC samples and 66 normal counterparts. miRNA associations with TP53 mutations and outcome were evaluated. RESULTS: A TP53 mutation was present in 58% of the tumors and TP53 mutations were significantly associated with a shorter recurrence-free survival. This association was stronger in the clinical subgroup of patients subjected to adjuvant therapy after surgery. The expression of 49 miRNAs was significantly associated with TP53 status. Among these 49, we identified a group of 12 miRNAs whose expression correlates with recurrence-free survival and a group of 4 miRNAs that correlates with cancer-specific survival. The two groups share three miRNAs. Importantly, miRNAs that correlate with survival are independent prognostic factors either when considered individually or as signatures. CONCLUSIONS: miRNAs expression associates with TP53 status and with reduced survival after surgical treatment of squamous cell carcinoma of the head and neck.

Modeling ductal carcinoma in situ: a HER2-Notch3 collaboration enables luminal filling.

A large fraction of ductal carcinoma in situ (DCIS), a non-invasive precursor lesion of invasive breast cancer, overexpresses the HER2/neu oncogene. The ducts of DCIS are abnormally filled with cells that evade apoptosis, but the underlying mechanisms remain incompletely understood. We overexpressed HER2 in mammary epithelial cells and observed growth factor-independent proliferation. When grown in extracellular matrix as three-dimensional spheroids, control cells developed a hollow lumen, but HER2-overexpressing cells populated the lumen by evading apoptosis. We demonstrate that HER2 overexpression in this cellular model of DCIS drives transcriptional upregulation of multiple components of the Notch survival pathway. Importantly, luminal filling required upregulation of a signaling pathway comprising Notch3, its cleaved intracellular domain and the transcriptional regulator HES1, resulting in elevated levels of c-MYC and cyclin D1. In line with HER2-Notch3 collaboration, drugs intercepting either arm reverted the DCIS-like phenotype. In addition, we report upregulation of Notch3 in hyperplastic lesions of HER2 transgenic animals, as well as an association between HER2 levels and expression levels of components of the Notch pathway in tumor specimens of breast cancer patients. Therefore, it is conceivable that the integration of the Notch and HER2 signaling pathways contributes to the pathophysiology of DCIS.

Modeling invasive breast cancer: growth factors propel progression of HER2-positive premalignant lesions.

The HER2/neu oncogene encodes a receptor-like tyrosine kinase whose overexpression in breast cancer predicts poor prognosis and resistance to conventional therapies. However, the mechanisms underlying aggressiveness of HER2 (human epidermal growth factor receptor 2)-overexpressing tumors remain incompletely understood. Because it assists epidermal growth factor (EGF) and neuregulin receptors, we overexpressed HER2 in MCF10A mammary cells and applied growth factors. HER2-overexpressing cells grown in extracellular matrix formed filled spheroids, which protruded outgrowths upon growth factor stimulation. Our transcriptome analyses imply a two-hit model for invasive growth: HER2-induced proliferation and evasion from anoikis generate filled structures, which are morphologically and transcriptionally analogous to preinvasive patients' lesions. In the second hit, EGF escalates signaling and transcriptional responses leading to invasive growth. Consistent with clinical relevance, a gene expression signature based on the HER2/EGF-activated transcriptional program can predict poorer prognosis of a subgroup of HER2-overexpressing patients. In conclusion, the integration of a three-dimensional cellular model and clinical data attributes progression of HER2-overexpressing lesions to EGF-like growth factors acting in the context of the tumor's microenvironment.

Gene expression signature is shared by patients with Alzheimer's disease and schizophrenia at the superior temporal gyrus.

BACKGROUND: Alzheimer's disease and Schizophrenia are two common diseases of the brain with significant differences in neuropathology, etiology and symptoms. This dissimilarity in the two diseases makes a comparison of the two ideal for detecting molecular substrates that are common to brain disorders in general. METHODS: In this study, we compared gene expression profiles across multiple brain areas, taken postmortem from patients with well-characterized Alzheimer's disease and Schizophrenia, and from cognitively normal control group with no neuro- or psychopathology. RESULTS: Although the totality of gene expression changes in the two diseases is dissimilar, a subset of genes appears to play a role in both diseases in specific brain regions. We find at Brodmann area 22, the superior temporal gyrus, a statistically significant number of genes with apparently disregulated expression in both diseases. Furthermore, we found genes that differentiate the two diseases from the control across multiple brain regions, and note that these genes were usually down-regulated. Brodmann area 8, part of the superior frontal cortex, is relatively abundant with them. CONCLUSION: We show overwhelming statistical evidence for Alzheimer's and Schizophrenia sharing a specific molecular background at the superior temporal gyrus. We suggest that impairment of the regulation of autophagy pathway is shared, in BA 22, by the two diseases.

Mutant p53(R175H) upregulates Twist1 expression and promotes epithelial-mesenchymal transition in immortalized prostate cells.

A mutation within one allele of the p53 tumor suppressor gene can inactivate the remaining wild-type allele in a dominant-negative manner and in some cases can exert an additional oncogenic activity, known as mutant p53 'gain of function' (GOF). To study the role of p53 mutations in prostate cancer and to discriminate between the dominant-negative effect and the GOF activity of mutant p53, we measured, using microarrays, the expression profiles of three immortalized prostate epithelial cultures expressing wild-type, inactivated p53 or mutated p53. Analysis of these gene expression profiles showed that both inactivated p53 and p53(R175H) mutant expression resulted in the upregulation of cell cycle progression genes. A second group, which was upregulated exclusively by mutant p53(R175H), was predominantly enriched in developmental genes. This group of genes included the Twist1, a regulator of metastasis and epithelial-mesenchymal transition (EMT). Twist1 levels were also elevated in metastatic prostate cancer-derived cell line DU145, in immortalized lung fibroblasts and in a subset of lung cancer samples, all in a mutant p53-dependent manner. p53(R175H) mutant bearing immortalized epithelial cells showed typical features of EMT, such as higher expression of mesenchymal markers, lower expression of epithelial markers and enhanced invasive properties in vitro. The mechanism by which p53(R175H) mutant induces Twist1 expression involves alleviation of the epigenetic repression. Our data suggest that Twist1 expression might be upregulated following p53 mutation in cancer cells.

Age-dependent spatial memory loss can be partially restored by immune activation.

Aging is often associated with a decline in hippocampus-dependent spatial memory. Here, we show that functional cell-mediated immunity is required for the maintenance of hippocampus-dependent spatial memory. Sudden imposition of immune compromise in young mice caused spatial memory impairment, whereas immune reconstitution reversed memory deficit in immune-deficient mice. Analysis of hippocampal gene expression suggested that immune-dependent spatial memory performance was associated with the expression of insulin-like growth factor (Igf1) and of genes encoding proteins related to presynaptic activity (Syt10, Cplx2). We further showed that memory loss in aged mice could be attributed to age-related attenuation of the immune response and could be reversed by immune system activation. Homeostatic-driven proliferation of lymphocytes, which expands the existing T cell repertoire, restored spatial memory deficits in aged mice. Thus, our results identify a novel function of the immune system in the maintenance of spatial memory and suggest an original approach for arresting or reversing age-associated memory loss.

TGFbeta-dependent gene expression profile during maturation of dendritic cells.

Primary immune response to pathogens involves the maturation of antigen-presenting dendritic cells (DC). Bacterial lipopolysacharride (LPS) is a potent inducer of DC maturation, whereas the transforming growth factor beta (TGFbeta) attenuates much of this process. Here, we analyzed the global gene expression pattern in LPS-treated bone marrow derived DC during inhibition of their maturation process by TGFbeta. Exposure of DC to LPS induces a pronounced cell response, manifested in altered expression of a large number of genes. Interestingly, TGFbeta did not affect most of the LPS responding genes. Nevertheless, analysis identified a subset of genes that did respond to TGFbeta, among them the two inflammatory cytokines interleukin (IL)-12 and IL-18. Expression of IL-12, the major proinflammatory cytokine secreted by mature DC, was downregulated by TGFbeta, whereas the expression level of the proinflammatory cytokine IL-18, known to potentiate the IL-12 effect, was upregulated. Expression of the peroxisome proliferator-activated receptor gamma (PPARgamma) increased in response to TGFbeta, concomitantly with reduced expression of chemokine receptor 7 (CCR7). This finding supports the possibility that TGFbeta-dependent inhibition of CCR7 expression in DC is mediated by PPARgamma.

Gene expression profiles of AML derived stem cells; similarity to hematopoietic stem cells.

Tumors contain a fraction of cancer stem cells that maintain the propagation of the disease. The CD34(+)CD38(-) cells, isolated from acute myeloid leukemia (AML), were shown to be enriched leukemic stem cells (LSC). We isolated the CD34(+)CD38(-) cell fraction from AML and compared their gene expression profiles to the CD34(+)CD38(+) cell fraction, using microarrays. We found 409 genes that were at least twofold over- or underexpressed between the two cell populations. These include underexpression of DNA repair, signal transduction and cell cycle genes, consistent with the relative quiescence of stem cells, and chromosomal aberrations and mutations of leukemic cells. Comparison of the LSC expression data to that of normal hematopoietic stem cells (HSC) revealed that 34% of the modulated genes are shared by both LSC and HSC, supporting the suggestion that the LSC originated within the HSC progenitors. We focused on the Notch pathway since Jagged-2, a Notch ligand was found to be overexpressed in the LSC samples. We show that DAPT, an inhibitor of gamma-secretase, a protease that is involved in Jagged and Notch signaling, inhibits LSC growth in colony formation assays. Identification of additional genes that regulate LSC self-renewal may provide new targets for therapy.

Antigen-chip technology for accessing global information about the state of the body.

Traditionally, immunologic diagnosis has been based on an attempt to correlate each disease with a specific immune reactivity, such as an antibody or a T-cell response to a single antigen specific for the disease entity. The state of the body, however, appears to be encoded by the immune system in collectives of reactivities and not by single reactivities. Here we describe our use of microarray technology and informatics to develop an antigen chip capable of detecting global patterns of antibodies binding to hundreds of antigens simultaneously. The patterns fashion diagnostic signatures.

Sorting points into neighborhoods (SPIN): data analysis and visualization by ordering distance matrices.

SUMMARY: We introduce a novel unsupervised approach for the organization and visualization of multidimensional data. At the heart of the method is a presentation of the full pairwise distance matrix of the data points, viewed in pseudocolor. The ordering of points is iteratively permuted in search of a linear ordering, which can be used to study embedded shapes. Several examples indicate how the shapes of certain structures in the data (elongated, circular and compact) manifest themselves visually in our permuted distance matrix. It is important to identify the elongated objects since they are often associated with a set of hidden variables, underlying continuous variation in the data. The problem of determining an optimal linear ordering is shown to be NP-Complete, and therefore an iterative search algorithm with O(n3) step-complexity is suggested. By using sorting points into neighborhoods, i.e. SPIN to analyze colon cancer expression data we were able to address the serious problem of sample heterogeneity, which hinders identification of metastasis related genes in our data. Our methodology brings to light the continuous variation of heterogeneity--starting with homogeneous tumor samples and gradually increasing the amount of another tissue. Ordering the samples according to their degree of contamination by unrelated tissue allows the separation of genes associated with irrelevant contamination from those related to cancer progression. AVAILABILITY: Software package will be available for academic users upon request.

Expression profiles of acute lymphoblastic and myeloblastic leukemias with ALL-1 rearrangements.

The ALL-1 gene is directly involved in 5-10% of acute lymphoblastic leukemias (ALLs) and acute myeloid leukemias (AMLs) by fusion to other genes or through internal rearrangements. DNA microarrays were used to determine expression profiles of ALLs and AMLs with ALL-1 rearrangements. These profiles distinguish those tumors from other ALLs and AMLs. The expression patterns of ALL-1-associated tumors, in particular ALLs, involve oncogenes, tumor suppressors, antiapoptotic genes, drug-resistance genes, etc., and correlate with the aggressive nature of the tumors. The genes whose expression differentiates between ALLs with and without ALL-1 rearrangement were further divided into several groups, enabling separation of ALL-1-associated ALLs into two subclasses. One of the groups included 43 genes that exhibited expression profiles closely linked to ALLs with ALL-1 rearrangements. Further, there were evident differences between the expression profiles of AMLs in which ALL-1 had undergone fusion to other genes and AMLs with partial duplication of ALL-1. The extensive analysis described here pinpointed genes that might have a direct role in pathogenesis.

Computational capacity of an odorant discriminator: the linear separability of curves.

We introduce and study an artificial neural network inspired by the probabilistic receptor affinity distribution model of olfaction. Our system consists of N sensory neurons whose outputs converge on a single processing linear threshold element. The system's aim is to model discrimination of a single target odorant from a large number p of background odorants within a range of odorant concentrations. We show that this is possible provided p does not exceed a critical value p(c) and calculate the critical capacity alpha(c) = p(c)/N. The critical capacity depends on the range of concentrations in which the discrimination is to be accomplished. If the olfactory bulb may be thought of as a collection of such processing elements, each responsible for the discrimination of a single odorant, our study provides a quantitative analysis of the potential computational properties of the olfactory bulb. The mathematical formulation of the problem we consider is one of determining the capacity for linear separability of continuous curves, embedded in a large-dimensional space. This is accomplished here by a numerical study, using a method that signals whether the discrimination task is realizable, together with a finite-size scaling analysis.

Statistical properties of contact vectors.

We study the statistical properties of contact vectors, a construct to characterize a protein's structure. The contact vector of an N-residue protein is a list of N integers n(i), representing the number of residues in contact with residue i. We study analytically (at mean-field level) and numerically the amount of structural information contained in a contact vector. Analytical calculations reveal that a large variance in the contact numbers reduces the degeneracy of the mapping between contact vectors and structures. Exact enumeration for lengths up to N=16 on the three-dimensional cubic lattice indicates that the growth rate of number of contact vectors as a function of N is only 3% less than that for contact maps. In particular, for compact structures we present numerical evidence that, practically, each contact vector corresponds to only a handful of structures. We discuss how this information can be used for better structure prediction.

Resampling method for unsupervised estimation of cluster validity.

We introduce a method for validation of results obtained by clustering analysis of data. The method is based on resampling the available data. A figure of merit that measures the stability of clustering solutions against resampling is introduced. Clusters that are stable against resampling give rise to local maxima of this figure of merit. This is presented first for a one-dimensional data set, for which an analytic approximation for the figure of merit is derived and compared with numerical measurements. Next, the applicability of the method is demonstrated for higher-dimensional data, including gene microarray expression data.

DNA microarrays identification of primary and secondary target genes regulated by p53.

The transcriptional program regulated by the tumor suppressor p53 was analysed using oligonucleotide microarrays. A human lung cancer cell line that expresses the temperature sensitive murine p53 was utilized to quantitate mRNA levels of various genes at different time points after shifting the temperature to 32 degrees C. Inhibition of protein synthesis by cycloheximide (CHX) was used to distinguish between primary and secondary target genes regulated by p53. In the absence of CHX, 259 and 125 genes were up or down-regulated respectively; only 38 and 24 of these genes were up and down-regulated by p53 also in the presence of CHX and are considered primary targets in this cell line. Cluster analysis of these data using the super paramagnetic clustering (SPC) algorithm demonstrate that the primary genes can be distinguished as a single cluster among a large pool of p53 regulated genes. This procedure identified additional genes that co-cluster with the primary targets and can also be classified as such genes. In addition to cell cycle (e.g. p21, TGF-beta, Cyclin E) and apoptosis (e.g. Fas, Bak, IAP) related genes, the primary targets of p53 include genes involved in many aspects of cell function, including cell adhesion (e.g. Thymosin, Smoothelin), signaling (e.g. H-Ras, Diacylglycerol kinase), transcription (e.g. ATF3, LISCH7), neuronal growth (e.g. Ninjurin, NSCL2) and DNA repair (e.g. BTG2, DDB2). The results suggest that p53 activates concerted opposing signals and exerts its effect through a diverse network of transcriptional changes that collectively alter the cell phenotype in response to stress.

Spin domains generate hierarchical ground state structure in J = +/-1 spin glasses.

Unbiased samples of ground states were generated for the short-range Ising spin glass with J(ij) = +/-1, in three dimensions. Clustering the ground states revealed their hierarchical structure, which is explained by correlated spin domains, serving as cores for macroscopic zero energy "excitations."

Coupled two-way clustering analysis of gene microarray data.

We present a coupled two-way clustering approach to gene microarray data analysis. The main idea is to identify subsets of the genes and samples, such that when one of these is used to cluster the other, stable and significant partitions emerge. The search for such subsets is a computationally complex task. We present an algorithm, based on iterative clustering, that performs such a search. This analysis is especially suitable for gene microarray data, where the contributions of a variety of biological mechanisms to the gene expression levels are entangled in a large body of experimental data. The method was applied to two gene microarray data sets, on colon cancer and leukemia. By identifying relevant subsets of the data and focusing on them we were able to discover partitions and correlations that were masked and hidden when the full dataset was used in the analysis. Some of these partitions have clear biological interpretation; others can serve to identify possible directions for future research.