From microarray to biology: an integrated experimental, statistical and in silico analysis of how the extracellular matrix modulates the phenotype of cancer cells

A statistically robust and biologically-based approach for analysis of microarray data is described that integrates independent biological knowledge and data with a global F-test for finding genes of interest that minimizes the need for replicates when used for hypothesis generation. First, each microarray is normalized to its noise level around zero. The microarray dataset is then globally adjusted by robust linear regression. Second, genes of interest that capture significant responses to experimental conditions are selected by finding those that express significantly higher variance than those expressing only technical variability. Clustering expression data and identifying expression-independent properties of genes of interest including upstream transcriptional regulatory elements (TREs), ontologies and networks or pathways organizes the data into a biologically meaningful system. We demonstrate that when the number of genes of interest is inconveniently large, identifying a subset of "beacon genes" representing the largest changes will identify pathways or networks altered by biological manipulation. The entire dataset is then used to complete the picture outlined by the "beacon genes." This allow construction of a structured model of a system that can generate biologically testable hypotheses. We illustrate this approach by comparing cells cultured on plastic or an extracellular matrix which organizes a dataset of over 2,000 genes of interest from a genome wide scan of transcription. The resulting model was confirmed by comparing the predicted pattern of TREs with experimental determination of active transcription factors

C. elegans SWAN-1 Binds to EGL-9 and Regulates HIF-1-Mediated Resistance to the Bacterial Pathogen Pseudomonas aeruginosa PAO1

Pseudomonas aeruginosa is a nearly ubiquitous human pathogen, and infections can be lethal to patients with impaired respiratory and immune systems. Prior studies have established that strong loss-of-function mutations in the egl-9 gene protect the nematode C. elegans from P. aeruginosa PAO1 fast killing. EGL-9 inhibits the HIF-1 transcription factor via two pathways. First, EGL-9 is the enzyme that targets HIF-1 for oxygen-dependent degradation via the VHL-1 E3 ligase. Second, EGL-9 inhibits HIF-1-mediated gene expression through a VHL-1-independent mechanism. Here, we show that a loss-of-function mutation in hif-1 suppresses P. aeruginosa PAO1 resistance in egl-9 mutants. Importantly, we find stabilization of HIF-1 protein is not sufficient to protect C. elegans from P. aeruginosa PAO1 fast killing. However, mutations that inhibit both EGL-9 pathways result in higher levels of HIF-1 activity and confer resistance to the pathogen. Using forward genetic screens, we identify additional mutations that confer resistance to P. aeruginosa. In genetic backgrounds that stabilize C. elegans HIF-1 protein, loss-of-function mutations in swan-1 increase the expression of hypoxia response genes and protect C. elegans from P. aeruginosa fast killing. SWAN-1 is an evolutionarily conserved WD-repeat protein belonging to the AN11 family. Yeast two-hybrid and co-immunoprecipitation assays show that EGL-9 forms a complex with SWAN-1. Additionally, we present genetic evidence that the DYRK kinase MBK-1 acts downstream of SWAN-1 to promote HIF-1-mediated transcription and to increase resistance to P. aeruginosa. These data support a model in which SWAN-1, MBK-1 and EGL-9 regulate HIF-1 transcriptional activity and modulate resistance to P. aeruginosa PAO1 fast killing

Robust Metabolic Responses to Varied Carbon Sources in Natural and Laboratory Strains of Saccharomyces cerevisiae

Understanding factors that regulate the metabolism and growth of an organism is of fundamental biologic interest. This study compared the influence of two different carbon substrates, dextrose and galactose, on the metabolic and growth rates of the yeast Saccharomyces cerevisiae. Yeast metabolic and growth rates varied widely depending on the metabolic substrate supplied. The metabolic and growth rates of a yeast strain maintained under long-term laboratory conditions was compared to strain isolated from natural condition when grown on different substrates. Previous studies had determined that there are numerous genetic differences between these two strains. However, the overall metabolic and growth rates of a wild isolate of yeast was very similar to that of a strain that had been maintained under laboratory conditions for many decades. This indicates that, at in least this case, metabolism and growth appear to be well buffered against genetic differences. Metabolic rate and cell number did not co-vary in a simple linear manner. When grown in either dextrose or galactose, both strains showed a growth pattern in which the number of cells continued to increase well after the metabolic rate began a sharp decline. Previous studied have reported that O2 consumption in S. cerevisiae grown in reduced dextrose levels were elevated compared to higher levels. Low dextrose levels have been proposed to induce caloric restriction and increase life span in yeast. However, there was no evidence that reduced levels of dextrose increased metabolic rates, measured by either O2 consumption or CO2 production, in the strains used in this study

The Unfolded Protein Response Is Not Necessary for the G1/S Transition, but It Is Required for Chromosome Maintenance in Saccharomyces cerevisiae

BACKGROUND: The unfolded protein response (UPR) is a eukaryotic signaling pathway, from the endoplasmic reticulum (ER) to the nucleus. Protein misfolding in the ER triggers the UPR. Accumulating evidence links the UPR in diverse aspects of cellular homeostasis. The UPR responds to the overall protein synthesis capacity and metabolic fluxes of the cell. Because the coupling of metabolism with cell division governs when cells start dividing, here we examined the role of UPR signaling in the timing of initiation of cell division and cell cycle progression, in the yeast Saccharomyces cerevisiae. METHODOLOGY/PRINCIPAL FINDINGS: We report that cells lacking the ER-resident stress sensor Ire1p, which cannot trigger the UPR, nonetheless completed the G1/S transition on time. Furthermore, loss of UPR signaling neither affected the nutrient and growth rate dependence of the G1/S transition, nor the metabolic oscillations that yeast cells display in defined steady-state conditions. Remarkably, however, loss of UPR signaling led to hypersensitivity to genotoxic stress and a ten-fold increase in chromosome loss. CONCLUSIONS/SIGNIFICANCE: Taken together, our results strongly suggest that UPR signaling is not necessary for the normal coupling of metabolism with cell division, but it has a role in genome maintenance. These results add to previous work that linked the UPR with cytokinesis in yeast. UPR signaling is conserved in all eukaryotes, and it malfunctions in a variety of diseases, including cancer. Therefore, our findings may be relevant to other systems, including humans

Discovery and Annotation of Functional Chromatin Signatures in the Human Genome

Transcriptional regulation in human cells is a complex process involving a multitude of regulatory elements encoded by the genome. Recent studies have shown that distinct chromatin signatures mark a variety of functional genomic elements and that subtle variations of these signatures mark elements with different functions. To identify novel chromatin signatures in the human genome, we apply a de novo pattern-finding algorithm to genome-wide maps of histone modifications. We recover previously known chromatin signatures associated with promoters and enhancers. We also observe several chromatin signatures with strong enrichment of H3K36me3 marking exons. Closer examination reveals that H3K36me3 is found on well-positioned nucleosomes at exon 5′ ends, and that this modification is a global mark of exon expression that also correlates with alternative splicing. Additionally, we observe strong enrichment of H2BK5me1 and H4K20me1 at highly expressed exons near the 5′ end, in contrast to the opposite distribution of H3K36me3-marked exons. Finally, we also recover frequently occurring chromatin signatures displaying enrichment of repressive histone modifications. These signatures mark distinct repeat sequences and are associated with distinct modes of gene repression. Together, these results highlight the rich information embedded in the human epigenome and underscore its value in studying gene regulation

Mapping the Hsp90 Genetic Interaction Network in Candida albicans Reveals Environmental Contingency and Rewired Circuitry

The molecular chaperone Hsp90 regulates the folding of diverse signal transducers in all eukaryotes, profoundly affecting cellular circuitry. In fungi, Hsp90 influences development, drug resistance, and evolution. Hsp90 interacts with ∼10% of the proteome in the model yeast Saccharomyces cerevisiae, while only two interactions have been identified in Candida albicans, the leading fungal pathogen of humans. Utilizing a chemical genomic approach, we mapped the C. albicans Hsp90 interaction network under diverse stress conditions. The chaperone network is environmentally contingent, and most of the 226 genetic interactors are important for growth only under specific conditions, suggesting that they operate downstream of Hsp90, as with the MAPK Hog1. Few interactors are important for growth in many environments, and these are poised to operate upstream of Hsp90, as with the protein kinase CK2 and the transcription factor Ahr1. We establish environmental contingency in the first chaperone network of a fungal pathogen, novel effectors upstream and downstream of Hsp90, and network rewiring over evolutionary time

CHD7 Targets Active Gene Enhancer Elements to Modulate ES Cell-Specific Gene Expression

CHD7 is one of nine members of the chromodomain helicase DNA–binding domain family of ATP–dependent chromatin remodeling enzymes found in mammalian cells. De novo mutation of CHD7 is a major cause of CHARGE syndrome, a genetic condition characterized by multiple congenital anomalies. To gain insights to the function of CHD7, we used the technique of chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP–Seq) to map CHD7 sites in mouse ES cells. We identified 10,483 sites on chromatin bound by CHD7 at high confidence. Most of the CHD7 sites show features of gene enhancer elements. Specifically, CHD7 sites are predominantly located distal to transcription start sites, contain high levels of H3K4 mono-methylation, found within open chromatin that is hypersensitive to DNase I digestion, and correlate with ES cell-specific gene expression. Moreover, CHD7 co-localizes with P300, a known enhancer-binding protein and strong predictor of enhancer activity. Correlations with 18 other factors mapped by ChIP–seq in mouse ES cells indicate that CHD7 also co-localizes with ES cell master regulators OCT4, SOX2, and NANOG. Correlations between CHD7 sites and global gene expression profiles obtained from Chd7+/+, Chd7+/−, and Chd7−/− ES cells indicate that CHD7 functions at enhancers as a transcriptional rheostat to modulate, or fine-tune the expression levels of ES–specific genes. CHD7 can modulate genes in either the positive or negative direction, although negative regulation appears to be the more direct effect of CHD7 binding. These data indicate that enhancer-binding proteins can limit gene expression and are not necessarily co-activators. Although ES cells are not likely to be affected in CHARGE syndrome, we propose that enhancer-mediated gene dysregulation contributes to disease pathogenesis and that the critical CHD7 target genes may be subject to positive or negative regulation

SHANK proteins limit integrin activation by directly interacting with Rap1 and R-Ras

SHANK3, a synaptic scaffold protein and actin regulator, is widely expressed outside of the central nervous system with predominantly unknown function. Solving the structure of the SHANK3 N-terminal region revealed that the SPN domain is an unexpected Ras-association domain with high affinity for GTP-bound Ras and Rap G-proteins. The role of Rap1 in integrin activation is well established but the mechanisms to antagonize it remain largely unknown. Here, we show that SHANK1 and SHANK3 act as integrin activation inhibitors by sequestering active Rap1 and R-Ras via the SPN domain and thus limiting their bioavailability at the plasma membrane. Consistently, SHANK3 silencing triggers increased plasma membrane Rap1 activity, cell spreading, migration and invasion. Autism-related mutations within the SHANK3 SPN domain (R12C and L68P) disrupt G-protein interaction and fail to counteract integrin activation along the Rap1-RIAM-talin axis in cancer cells and neurons. Altogether, we establish SHANKs as critical regulators of G-protein signalling and integrin-dependent processes

The Genomic Distribution and Function of Histone Variant HTZ-1 during C. elegans Embryogenesis

In all eukaryotes, histone variants are incorporated into a subset of nucleosomes to create functionally specialized regions of chromatin. One such variant, H2A.Z, replaces histone H2A and is required for development and viability in all animals tested to date. However, the function of H2A.Z in development remains unclear. Here, we use ChIP-chip, genetic mutation, RNAi, and immunofluorescence microscopy to interrogate the function of H2A.Z (HTZ-1) during embryogenesis in Caenorhabditis elegans, a key model of metazoan development. We find that HTZ-1 is expressed in every cell of the developing embryo and is essential for normal development. The sites of HTZ-1 incorporation during embryogenesis reveal a genome wrought by developmental processes. HTZ-1 is incorporated upstream of 23% of C. elegans genes. While these genes tend to be required for development and occupied by RNA polymerase II, HTZ-1 incorporation does not specify a stereotypic transcription program. The data also provide evidence for unexpectedly widespread independent regulation of genes within operons during development; in 37% of operons, HTZ-1 is incorporated upstream of internally encoded genes. Fewer sites of HTZ-1 incorporation occur on the X chromosome relative to autosomes, which our data suggest is due to a paucity of developmentally important genes on X, rather than a direct function for HTZ-1 in dosage compensation. Our experiments indicate that HTZ-1 functions in establishing or maintaining an essential chromatin state at promoters regulated dynamically during C. elegans embryogenesis

Epigenetic Analysis of KSHV Latent and Lytic Genomes

Epigenetic modifications of the herpesviral genome play a key role in the transcriptional control of latent and lytic genes during a productive viral lifecycle. In this study, we describe for the first time a comprehensive genome-wide ChIP-on-Chip analysis of the chromatin associated with the Kaposi's sarcoma-associated herpesvirus (KSHV) genome during latency and lytic reactivation. Depending on the gene expression class, different combinations of activating [acetylated H3 (AcH3) and H3K4me3] and repressive [H3K9me3 and H3K27me3] histone modifications are associated with the viral latent genome, which changes upon reactivation in a manner that is correlated with their expression. Specifically, both the activating marks co-localize on the KSHV latent genome, as do the repressive marks. However, the activating and repressive histone modifications are mutually exclusive of each other on the bulk of the latent KSHV genome. The genomic region encoding the IE genes ORF50 and ORF48 possesses the features of a bivalent chromatin structure characterized by the concomitant presence of the activating H3K4me3 and the repressive H3K27me3 marks during latency, which rapidly changes upon reactivation with increasing AcH3 and H3K4me3 marks and decreasing H3K27me3. Furthermore, EZH2, the H3K27me3 histone methyltransferase of the Polycomb group proteins (PcG), colocalizes with the H3K27me3 mark on the entire KSHV genome during latency, whereas RTA-mediated reactivation induces EZH2 dissociation from the genomic regions encoding IE and E genes concurrent with decreasing H3K27me3 level and increasing IE/E lytic gene expression. Moreover, either the inhibition of EZH2 expression by a small molecule inhibitor DZNep and RNAi knockdown, or the expression of H3K27me3-specific histone demethylases apparently induced the KSHV lytic gene expression cascade. These data indicate that histone modifications associated with the KSHV latent genome are involved in the regulation of latency and ultimately in the control of the temporal and sequential expression of the lytic gene cascade. In addition, the PcG proteins play a critical role in the control of KSHV latency by maintaining a reversible heterochromatin on the KSHV lytic genes. Thus, the regulation of the spatial and temporal association of the PcG proteins with the KSHV genome may be crucial for propagating the KSHV lifecycle