Large-scale genetic perturbations reveal regulatory networks and an abundance of gene-specific repressors.

To understand regulatory systems, it would be useful to uniformly determine how different components contribute to the expression of all other genes. We therefore monitored mRNA expression genome-wide, for individual deletions of one-quarter of yeast genes, focusing on (putative) regulators. The resulting genetic perturbation signatures reflect many different properties. These include the architecture of protein complexes and pathways, identification of expression changes compatible with viability, and the varying responsiveness to genetic perturbation. The data are assembled into a genetic perturbation network that shows different connectivities for different classes of regulators. Four feed-forward loop (FFL) types are overrepresented, including incoherent type 2 FFLs that likely represent feedback. Systematic transcription factor classification shows a surprisingly high abundance of gene-specific repressors, suggesting that yeast chromatin is not as generally restrictive to transcription as is often assumed. The data set is useful for studying individual genes and for discovering properties of an entire regulatory system.

Functional overlap and regulatory links shape genetic interactions between signaling pathways.

To understand relationships between phosphorylation-based signaling pathways, we analyzed 150 deletion mutants of protein kinases and phosphatases in S. cerevisiae using DNA microarrays. Downstream changes in gene expression were treated as a phenotypic readout. Double mutants with synthetic genetic interactions were included to investigate genetic buffering relationships such as redundancy. Three types of genetic buffering relationships are identified: mixed epistasis, complete redundancy, and quantitative redundancy. In mixed epistasis, the most common buffering relationship, different gene sets respond in different epistatic ways. Mixed epistasis arises from pairs of regulators that have only partial overlap in function and that are coupled by additional regulatory links such as repression of one by the other. Such regulatory modules confer the ability to control different combinations of processes depending on condition or context. These properties likely contribute to the evolutionary maintenance of paralogs and indicate a way in which signaling pathways connect for multiprocess control.

Molecular mechanisms that distinguish TFIID housekeeping from regulatable SAGA promoters.

An important distinction is frequently made between constitutively expressed housekeeping genes versus regulated genes. Although generally characterized by different DNA elements, chromatin architecture and cofactors, it is not known to what degree promoter classes strictly follow regulatability rules and which molecular mechanisms dictate such differences. We show that SAGA-dominated/TATA-box promoters are more responsive to changes in the amount of activator, even compared to TFIID/TATA-like promoters that depend on the same activator Hsf1. Regulatability is therefore an inherent property of promoter class. Further analyses show that SAGA/TATA-box promoters are more dynamic because TATA-binding protein recruitment through SAGA is susceptible to removal by Mot1. In addition, the nucleosome configuration upon activator depletion shifts on SAGA/TATA-box promoters and seems less amenable to preinitiation complex formation. The results explain the fundamental difference between housekeeping and regulatable genes, revealing an additional facet of combinatorial control: an activator can elicit a different response dependent on core promoter class.

EBV MicroRNA BART16 Suppresses Type I IFN Signaling.

Type I IFNs play critical roles in orchestrating the antiviral defense by inducing direct antiviral activities and shaping the adaptive immune response. Viruses have evolved numerous strategies to specifically interfere with IFN production or its downstream mediators, thereby allowing successful infection of the host to occur. The prototypic human gammaherpesvirus EBV, which is associated with infectious mononucleosis and malignant tumors, harbors many immune-evasion proteins that manipulate the adaptive and innate immune systems. In addition to proteins, the virus encodes >40 mature microRNAs for which the functions remain largely unknown. In this article, we identify EBV-encoded miR-BART16 as a novel viral immune-evasion factor that interferes with the type I IFN signaling pathway. miR-BART16 directly targets CREB-binding protein, a key transcriptional coactivator in IFN signaling, thereby inducing CREB-binding protein downregulation in EBV-transformed B cells and gastric carcinoma cells. miR-BART16 abrogates the production of IFN-stimulated genes in response to IFN-α stimulation and it inhibits the antiproliferative effect of IFN-α on latently infected BL cells. By obstructing the type I IFN-induced antiviral response, miR-BART16 provides a means to facilitate the establishment of latent EBV infection and enhance viral replication.

The specificity and topology of chromatin interaction pathways in yeast.

Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data are assembled into a network of chromatin interaction pathways. The network is function based, has a branched, interconnected topology, and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts additional interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.

A consensus of core protein complex compositions for Saccharomyces cerevisiae.

Analyses of biological processes would benefit from accurate definitions of protein complexes. High-throughput mass spectrometry data offer the possibility of systematically defining protein complexes; however, the predicted compositions vary substantially depending on the algorithm applied. We determine consensus compositions for 409 core protein complexes from Saccharomyces cerevisiae by merging previous predictions with a new approach. Various analyses indicate that the consensus is comprehensive and of high quality. For 85 out of 259 complexes not recorded in GO, literature search revealed strong support in the form of coprecipitation. New complexes were verified by an independent interaction assay and by gene expression profiling of strains with deleted subunits, often revealing which cellular processes are affected. The consensus complexes are available in various formats, including a merge with GO, resulting in 518 protein complex compositions. The utility is further demonstrated by comparison with binary interaction data to reveal interactions between core complexes.

FOXO target gene CTDSP2 regulates cell cycle progression through Ras and p21(Cip1/Waf1).

Activity of FOXO (forkhead box O) transcription factors is inhibited by growth factor-PI3K (phosphoinositide 3-kinase)-PKB (protein kinase B)/Akt signalling to control a variety of cellular processes including cell cycle progression. Through comparative analysis of a number of microarray datasets we identified a set of genes commonly regulated by FOXO proteins and PI3K-PKB/Akt, which includes CTDSP2 (C-terminal domain small phosphatase 2). We validated CTDSP2 as a genuine FOXO target gene and show that ectopic CTDSP2 can induce cell cycle arrest. We analysed transcriptional regulation after CTDSP2 expression and identified extensive regulation of genes involved in cell cycle progression, which depends on the phosphatase activity of CTDSP2. The most notably regulated gene is the CDK (cyclin-dependent kinase) inhibitor p21(Cip1/Waf1) and in the present study we show that p21(Cip1/Waf1) is partially responsible for the cell cycle arrest through decreasing cyclin-CDK activity. Our data suggest that CTDSP2 induces p21(Cip1/Waf1) through increasing the activity of Ras. As has been described previously, Ras induces p21(Cip1/Waf1) through p53-dependent and p53-independent pathways and indeed both p53 and MEK inhibition can mitigate the CTDSP2-induced p21(Cip1/Waf1) mRNA up-regulation. In support of Ras activation by CTDSP2, depletion of endogenous CTDSP2 results in reduced Ras activity and thus CTDSP2 seems to be part of a larger set of genes regulated by FOXO proteins, which increase growth factor signalling upon FOXO activation.

A high-resolution gene expression atlas of epistasis between gene-specific transcription factors exposes potential mechanisms for genetic interactions.

BACKGROUND: Genetic interactions, or non-additive effects between genes, play a crucial role in many cellular processes and disease. Which mechanisms underlie these genetic interactions has hardly been characterized. Understanding the molecular basis of genetic interactions is crucial in deciphering pathway organization and understanding the relationship between genotype, phenotype and disease. RESULTS: To investigate the nature of genetic interactions between gene-specific transcription factors (GSTFs) in Saccharomyces cerevisiae, we systematically analyzed 72 GSTF pairs by gene expression profiling double and single deletion mutants. These pairs were selected through previously published growth-based genetic interactions as well as through similarity in DNA binding properties. The result is a high-resolution atlas of gene expression-based genetic interactions that provides systems-level insight into GSTF epistasis. The atlas confirms known genetic interactions and exposes new ones. Importantly, the data can be used to investigate mechanisms that underlie individual genetic interactions. Two molecular mechanisms are proposed, "buffering by induced dependency" and "alleviation by derepression". CONCLUSIONS: These mechanisms indicate how negative genetic interactions can occur between seemingly unrelated parallel pathways and how positive genetic interactions can indirectly expose parallel rather than same-pathway relationships. The focus on GSTFs is important for understanding the transcription regulatory network of yeast as it uncovers details behind many redundancy relationships, some of which are completely new. In addition, the study provides general insight into the complex nature of epistasis and proposes mechanistic models for genetic interactions, the majority of which do not fall into easily recognizable within- or between-pathway relationships.

Cell cycle population effects in perturbation studies.

Growth condition perturbation or gene function disruption are commonly used strategies to study cellular systems. Although it is widely appreciated that such experiments may involve indirect effects, these frequently remain uncharacterized. Here, analysis of functionally unrelated Saccharyomyces cerevisiae deletion strains reveals a common gene expression signature. One property shared by these strains is slower growth, with increased presence of the signature in more slowly growing strains. The slow growth signature is highly similar to the environmental stress response (ESR), an expression response common to diverse environmental perturbations. Both environmental and genetic perturbations result in growth rate changes. These are accompanied by a change in the distribution of cells over different cell cycle phases. Rather than representing a direct expression response in single cells, both the slow growth signature and ESR mainly reflect a redistribution of cells over different cell cycle phases, primarily characterized by an increase in the G1 population. The findings have implications for any study of perturbation that is accompanied by growth rate changes. Strategies to counter these effects are presented and discussed.

Two distinct repressive mechanisms for histone 3 lysine 4 methylation through promoting 3'-end antisense transcription.

Histone H3 di- and trimethylation on lysine 4 are major chromatin marks that correlate with active transcription. The influence of these modifications on transcription itself is, however, poorly understood. We have investigated the roles of H3K4 methylation in Saccharomyces cerevisiae by determining genome-wide expression-profiles of mutants in the Set1 complex, COMPASS, that lays down these marks. Loss of H3K4 trimethylation has virtually no effect on steady-state or dynamically-changing mRNA levels. Combined loss of H3K4 tri- and dimethylation results in steady-state mRNA upregulation and delays in the repression kinetics of specific groups of genes. COMPASS-repressed genes have distinct H3K4 methylation patterns, with enrichment of H3K4me3 at the 3'-end, indicating that repression is coupled to 3'-end antisense transcription. Further analyses reveal that repression is mediated by H3K4me3-dependent 3'-end antisense transcription in two ways. For a small group of genes including PHO84, repression is mediated by a previously reported trans-effect that requires the antisense transcript itself. For the majority of COMPASS-repressed genes, however, it is the process of 3'-end antisense transcription itself that is the important factor for repression. Strand-specific qPCR analyses of various mutants indicate that this more prevalent mechanism of COMPASS-mediated repression requires H3K4me3-dependent 3'-end antisense transcription to lay down H3K4me2, which seems to serve as the actual repressive mark. Removal of the 3'-end antisense promoter also results in derepression of sense transcription and renders sense transcription insensitive to the additional loss of SET1. The derepression observed in COMPASS mutants is mimicked by reduction of global histone H3 and H4 levels, suggesting that the H3K4me2 repressive effect is linked to establishment of a repressive chromatin structure. These results indicate that in S. cerevisiae, the non-redundant role of H3K4 methylation by Set1 is repression, achieved through promotion of 3'-end antisense transcription to achieve specific rather than global effects through two distinct mechanisms.

Retinoic acid-dependent and -independent gene-regulatory pathways of Pitx3 in meso-diencephalic dopaminergic neurons.

Development of meso-diencephalic dopamine (mdDA) neurons requires the combined actions of the orphan nuclear receptor Nurr1 and the paired-like homeobox transcription factor Pitx3. Whereas all mdDA neurons require Nurr1 for expression of Th and survival, dependence on Pitx3 is displayed only by the mdDA subpopulation that will form the substantia nigra (SNc). Previously, we have demonstrated that Pitx3(-/-) embryos lack the expression of the retinoic acid (RA)-generating enzyme Ahd2, which is normally selectively expressed in the Pitx3-dependent DA neurons of the SNc. Restoring RA signaling in Pitx3(-/-) embryos revealed a selective dependence of SNc neurons on the presence of RA for differentiation into Th-positive neurons and maintenance throughout embryonic development. Whereas these data are suggestive of an important developmental role for RA in neurons of the SNc, it remained unclear whether other Nurr1 and Pitx3 target genes depend on RA signaling in a manner similar to Th. In the search for genes that were affected in Pitx3-deficient mdDA neurons and restored upon embryonic RA treatment, we provide evidence that Delta-like 1, D2R (Drd2) and Th are regulated by Pitx3 and RA signaling, which influences the mdDA terminal differentiated phenotype. Furthermore, we show that regulation of Ahd2-mediated RA signaling represents only one aspect of the Pitx3 downstream cascade, as Vmat2, Dat, Ahd2 (Aldh1a1), En1, En2 and Cck were unaffected by RA treatment and are (subset) specifically modulated by Pitx3. In conclusion, our data reveal several RA-dependent and -independent aspects of the Pitx3-regulated gene cascade, suggesting that Pitx3 acts on multiple levels in the molecular subset-specification of mdDA neurons.

Yeast glucose pathways converge on the transcriptional regulation of trehalose biosynthesis.

BACKGROUND: Cellular glucose availability is crucial for the functioning of most biological processes. Our understanding of the glucose regulatory system has been greatly advanced by studying the model organism Saccharomyces cerevisiae, but many aspects of this system remain elusive. To understand the organisation of the glucose regulatory system, we analysed 91 deletion mutants of the different glucose signalling and metabolic pathways in Saccharomyces cerevisiae using DNA microarrays. RESULTS: In general, the mutations do not induce pathway-specific transcriptional responses. Instead, one main transcriptional response is discerned, which varies in direction to mimic either a high or a low glucose response. Detailed analysis uncovers established and new relationships within and between individual pathways and their members. In contrast to signalling components, metabolic components of the glucose regulatory system are transcriptionally more frequently affected. A new network approach is applied that exposes the hierarchical organisation of the glucose regulatory system. CONCLUSIONS: The tight interconnection between the different pathways of the glucose regulatory system is reflected by the main transcriptional response observed. Tps2 and Tsl1, two enzymes involved in the biosynthesis of the storage carbohydrate trehalose, are predicted to be the most downstream transcriptional components. Epistasis analysis of tps2Δ double mutants supports this prediction. Although based on transcriptional changes only, these results suggest that all changes in perceived glucose levels ultimately lead to a shift in trehalose biosynthesis.

Redox-sensitive cysteines bridge p300/CBP-mediated acetylation and FoxO4 activity.

Cellular damage invoked by reactive oxygen species plays a key role in the pathobiology of cancer and aging. Forkhead box class O (FoxO) transcription factors are involved in various cellular processes including cell cycle regulation, apoptosis and resistance to reactive oxygen species, and studies in animal models have shown that these transcription factors are of vital importance in tumor suppression, stem cell maintenance and lifespan extension. Here we report that the activity of FoxO in human cells is directly regulated by the cellular redox state through a unique mechanism in signal transduction. We show that reactive oxygen species induce the formation of cysteine-thiol disulfide-dependent complexes of FoxO and the p300/CBP acetyltransferase, and that modulation of FoxO biological activity by p300/CBP-mediated acetylation is fully dependent on the formation of this redox-dependent complex. These findings directly link cellular redox status to the activity of the longevity protein FoxO.

Adaptable gene-specific dye bias correction for two-channel DNA microarrays.

DNA microarray technology is a powerful tool for monitoring gene expression or for finding the location of DNA-bound proteins. DNA microarrays can suffer from gene-specific dye bias (GSDB), causing some probes to be affected more by the dye than by the sample. This results in large measurement errors, which vary considerably for different probes and also across different hybridizations. GSDB is not corrected by conventional normalization and has been difficult to address systematically because of its variance. We show that GSDB is influenced by label incorporation efficiency, explaining the variation of GSDB across different hybridizations. A correction method (Gene- And Slide-Specific Correction, GASSCO) is presented, whereby sequence-specific corrections are modulated by the overall bias of individual hybridizations. GASSCO outperforms earlier methods and works well on a variety of publically available datasets covering a range of platforms, organisms and applications, including ChIP on chip. A sequence-based model is also presented, which predicts which probes will suffer most from GSDB, useful for microarray probe design and correction of individual hybridizations. Software implementing the method is publicly available.

Improved genome-wide localization by ChIP-chip using double-round T7 RNA polymerase-based amplification.

Chromatin immunoprecipitation combined with DNA microarrays (ChIP-chip) is a powerful technique to detect in vivo protein-DNA interactions. Due to low yields, ChIP assays of transcription factors generally require amplification of immunoprecipitated genomic DNA. Here, we present an adapted linear amplification method that involves two rounds of T7 RNA polymerase amplification (double-T7). Using this we could successfully amplify as little as 0.4 ng of ChIP DNA to sufficient amounts for microarray analysis. In addition, we compared the double-T7 method to the ligation-mediated polymerase chain reaction (LM-PCR) method in a ChIP-chip of the yeast transcription factor Gsm1p. The double-T7 protocol showed lower noise levels and stronger binding signals compared to LM-PCR. Both LM-PCR and double-T7 identified strongly bound genomic regions, but the double-T7 method increased sensitivity and specificity to allow detection of weaker binding sites.

Growth plate expression profiling: Large and small breed dogs provide new insights in endochondral bone formation.

The difference in the adult height of mammals, and hence in endochondral bone formation, is not yet fully understood and may serve to identify targets for bone and cartilage regeneration. In line with this hypothesis, the intra-species disparity between the adult height of Great Danes and Miniature Poodles was investigated at a transcriptional level. Microarray analysis of the growth plate of five Great Danes and five Miniature Poodles revealed 2,981 unique genes that were differentially expressed, including many genes with an unknown role in skeletal development. A signaling pathway impact analysis indicated activation of the cell cycle, extracellular matrix receptor interaction and the tight junction pathway, and inhibition of pathways associated with inflammation and the complement cascade. In additional validation steps, the gene expression profile of the separate growth plate zones for both dog breeds were determined. Given that the BMP signaling is known for its crucial role in skeletal development and fracture healing, and BMP-2 is used in orthopaedic and spine procedures for bone augmentation, further investigations concentrated on the BMP pathway.The canonical BMP-2 and BMP-6 signaling pathway was activated in the Great Danes compared to Miniature Poodles. In conclusion, investigating the differential expression of genes involved in endochondral bone formation in small and large breed dogs, could be a game changing strategy to provide new insights in growth plate development and identify new targets for bone and cartilage regeneration. © 2017 The Authors. Journal of Orthopaedic Research® published by Wiley Periodicals, Inc. on behalf of the Orthopaedic Research Society. J Orthop Res 36:138-148, 2018.

FOXO Transcription Factors Both Suppress and Support Breast Cancer Progression.

FOXO transcription factors are regulators of cellular homeostasis and putative tumor suppressors, yet the role of FOXO in cancer progression remains to be determined. The data on FOXO function, particularly for epithelial cancers, are fragmentary and come from studies that focused on isolated aspects of cancer. To clarify the role of FOXO in epithelial cancer progression, we characterized the effects of inducible FOXO activation and loss in a mouse model of metastatic invasive lobular carcinoma. Strikingly, either activation or loss of FOXO function suppressed tumor growth and metastasis. We show that the multitude of cellular processes critically affected by FOXO function include proliferation, survival, redox homeostasis, and PI3K signaling, all of which must be carefully balanced for tumor cells to thrive.Significance: FOXO proteins are not solely tumor suppressors, but also support tumor growth and metastasis by regulating a multitude of cellular processes essential for tumorigenesis. Cancer Res; 78(9); 2356-69. ©2018 AACR.

An endometrial gene expression signature accurately predicts recurrent implantation failure after IVF.

The primary limiting factor for effective IVF treatment is successful embryo implantation. Recurrent implantation failure (RIF) is a condition whereby couples fail to achieve pregnancy despite consecutive embryo transfers. Here we describe the collection of gene expression profiles from mid-luteal phase endometrial biopsies (n = 115) from women experiencing RIF and healthy controls. Using a signature discovery set (n = 81) we identify a signature containing 303 genes predictive of RIF. Independent validation in 34 samples shows that the gene signature predicts RIF with 100% positive predictive value (PPV). The strength of the RIF associated expression signature also stratifies RIF patients into distinct groups with different subsequent implantation success rates. Exploration of the expression changes suggests that RIF is primarily associated with reduced cellular proliferation. The gene signature will be of value in counselling and guiding further treatment of women who fail to conceive upon IVF and suggests new avenues for developing intervention.

Twenty-two novel mutations in the lysosomal alpha-glucosidase gene (GAA) underscore the genotype-phenotype correlation in glycogen storage disease type II.

Patients with glycogen storage disease type II (GSDII, Pompe disease) suffer from progressive muscle weakness due to acid alpha-glucosidase deficiency. The disease is inherited as an autosomal recessive trait with a spectrum of clinical phenotypes. We have investigated 29 cases of GSDII and thereby identified 55 pathogenic mutations of the acid alpha-glucosidase gene (GAA) encoding acid maltase. There were 34 different mutations identified, 22 of which were novel. All of the missense mutations and two other mutations with an unpredictable effect on acid alpha-glucosidase synthesis and function were transiently expressed in COS cells. The effect of a novel splice-site mutation was investigated by real-time PCR analysis. The outcome of our analysis underscores the notion that the clinical phenotype of GSDII is largely dictated by the nature of the mutations in the GAA alleles. This genotype-phenotype correlation makes DNA analysis a valuable tool to help predict the clinical course of the disease.

A Low-Cost Computerized System to Monitor Running Performance and Circadian Rhythms of Twenty Mice Simultaneously.

This paper describes the design and functioning of a low-cost computerized system for monitoring the voluntary activity of mice in running wheels. The required software is written in Turbo Pascal(r) and provided via the Internet (http://www.eur.nl/fgg/ch1/rodent.html). The system accommodates the simultaneous monitoring of 20 animals over a virtually unlimited period. Two applications of the system are presented; one monitors the circadian rhythm of mice, and the other tests muscle strength and endurance.