Carm1-arginine methylation of the transcription factor C/EBPα regulates transdifferentiation velocity.

Here, we describe how the speed of C/EBPα-induced B cell to macrophage transdifferentiation (BMT) can be regulated, using both mouse and human models. The identification of a mutant of C/EBPα (C/EBPαR35A) that greatly accelerates BMT helped to illuminate the mechanism. Thus, incoming C/EBPα binds to PU.1, an obligate partner expressed in B cells, leading to the release of PU.1 from B cell enhancers, chromatin closing and silencing of the B cell program. Released PU.1 redistributes to macrophage enhancers newly occupied by C/EBPα, causing chromatin opening and activation of macrophage genes. All these steps are accelerated by C/EBPαR35A, initiated by its increased affinity for PU.1. Wild-type C/EBPα is methylated by Carm1 at arginine 35 and the enzyme's perturbations modulate BMT velocity as predicted from the observations with the mutant. Increasing the proportion of unmethylated C/EBPα in granulocyte/macrophage progenitors by inhibiting Carm1 biases the cell's differentiation toward macrophages, suggesting that cell fate decision velocity and lineage directionality are closely linked processes.

The trophectoderm acts as a niche for the inner cell mass through C/EBPα-regulated IL-6 signaling.

IL-6 has been shown to be required for somatic cell reprogramming into induced pluripotent stem cells (iPSCs). However, how Il6 expression is regulated and whether it plays a role during embryo development remains unknown. Here, we describe that IL-6 is necessary for C/EBPα-enhanced reprogramming of B cells into iPSCs but not for B cell to macrophage transdifferentiation. C/EBPα overexpression activates both Il6 and Il6ra genes in B cells and in PSCs. In embryo development, Cebpa is enriched in the trophectoderm of blastocysts together with Il6, while Il6ra is mostly expressed in the inner cell mass (ICM). In addition, Il6 expression in blastocysts requires Cebpa. Blastocysts secrete IL-6 and neutralization of the cytokine delays the morula to blastocyst transition. The observed requirement of C/EBPα-regulated IL-6 signaling for pluripotency during somatic cell reprogramming thus recapitulates a physiologic mechanism in which the trophectoderm acts as niche for the ICM through the secretion of IL-6.

Real age prediction from the transcriptome with RAPToR.

Transcriptomic data is often affected by uncontrolled variation among samples that can obscure and confound the effects of interest. This variation is frequently due to unintended differences in developmental stages between samples. The transcriptome itself can be used to estimate developmental progression, but existing methods require many samples and do not estimate a specimen's real age. Here we present real-age prediction from transcriptome staging on reference (RAPToR), a computational method that precisely estimates the real age of a sample from its transcriptome, exploiting existing time-series data as reference. RAPToR works with whole animal, dissected tissue and single-cell data for the most common animal models, humans and even for non-model organisms lacking reference data. We show that RAPToR can be used to remove age as a confounding factor and allow recovery of a signal of interest in differential expression analysis. RAPToR will be especially useful in large-scale single-organism profiling because it eliminates the need for accurate staging or synchronisation before profiling.

Neuronal perception of the social environment generates an inherited memory that controls the development and generation time of C. elegans.

An old and controversial question in biology is whether information perceived by the nervous system of an animal can "cross the Weismann barrier" to alter the phenotypes and fitness of their progeny. Here, we show that such intergenerational transmission of sensory information occurs in the model organism, C. elegans, with a major effect on fitness. Specifically, that perception of social pheromones by chemosensory neurons controls the post-embryonic timing of the development of one tissue, the germline, relative to others in the progeny of an animal. Neuronal perception of the social environment thus intergenerationally controls the generation time of this animal.

Single cell RNA-seq identifies the origins of heterogeneity in efficient cell transdifferentiation and reprogramming.

Forced transcription factor expression can transdifferentiate somatic cells into other specialised cell types or reprogram them into induced pluripotent stem cells (iPSCs) with variable efficiency. To better understand the heterogeneity of these processes, we used single-cell RNA sequencing to follow the transdifferentation of murine pre-B cells into macrophages as well as their reprogramming into iPSCs. Even in these highly efficient systems, there was substantial variation in the speed and path of fate conversion. We predicted and validated that these differences are inversely coupled and arise in the starting cell population, with Mychigh large pre-BII cells transdifferentiating slowly but reprogramming efficiently and Myclow small pre-BII cells transdifferentiating rapidly but failing to reprogram. Strikingly, differences in Myc activity predict the efficiency of reprogramming across a wide range of somatic cell types. These results illustrate how single cell expression and computational analyses can identify the origins of heterogeneity in cell fate conversion processes.

Maternal age generates phenotypic variation in Caenorhabditis elegans.

Genetically identical individuals that grow in the same environment often show substantial phenotypic variation within populations of organisms as diverse as bacteria, nematodes, rodents and humans. With some exceptions, the causes are poorly understood. Here we show that isogenic Caenorhabditis elegans nematodes vary in their size at hatching, speed of development, growth rate, starvation resistance, fecundity, and also in the rate of development of their germline relative to that of somatic tissues. We show that the primary cause of this variation is the age of an individual's mother, with the progeny of young mothers exhibiting several phenotypic impairments. We identify age-dependent changes in the maternal provisioning of the lipoprotein complex vitellogenin to embryos as the molecular mechanism that underlies the variation in multiple traits throughout the life of an animal. The production of sub-optimal progeny by young mothers may reflect a trade-off between the competing fitness traits of a short generation time and the survival and fecundity of the progeny.

C/EBPα creates elite cells for iPSC reprogramming by upregulating Klf4 and increasing the levels of Lsd1 and Brd4.

Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) is typically inefficient and has been explained by elite-cell and stochastic models. We recently reported that B cells exposed to a pulse of C/EBPα (Bα' cells) behave as elite cells, in that they can be rapidly and efficiently reprogrammed into iPSCs by the Yamanaka factors OSKM. Here we show that C/EBPα post-transcriptionally increases the abundance of several hundred proteins, including Lsd1, Hdac1, Brd4, Med1 and Cdk9, components of chromatin-modifying complexes present at super-enhancers. Lsd1 was found to be required for B cell gene silencing and Brd4 for the activation of the pluripotency program. C/EBPα also promotes chromatin accessibility in pluripotent cells and upregulates Klf4 by binding to two haematopoietic enhancers. Bα' cells share many properties with granulocyte/macrophage progenitors, naturally occurring elite cells that are obligate targets for leukaemic transformation, whose formation strictly requires C/EBPα.

Adenomatous polyposis coli (APC) regulates miR17-92 cluster through ß-catenin pathway in colorectal cancer.

Adenomatous polyposis coli (APC) mutation is the most common genetic change in sporadic colorectal cancer (CRC). Although deregulations of miRNAs have been frequently reported in this malignancy, APC-regulated miRNAs have not been extensively documented. Here, by using an APC-inducible cell line and array analysis, we identified a total of 26 deregulated miRNAs. Among them, members of miR-17-92 cluster were dramatically inhibited by APC and induced by enforced expression of ß-catenin. Furthermore, we demonstrate that activated ß-catenin resulted from APC loss binds to and activates the miR-17-92 promoter. Notably, enforced expression of miR-19a overrides APC tumor suppressor activity, and knockdown of miR-19a in cancer cells with compromised APC function reduced their aggressive features in vitro. Finally, we observed that expression of miR-19a significantly correlates with ß-catenin levels in colorectal cancer specimens, and it is associated to the aggressive stage of tumor progression. Thus, our study reveals that miR-17-92 cluster is directly regulated by APC/ß-catenin pathway and could be a potential therapeutic target in colon cancers with aberrant APC/ß-catenin signaling.

The propagation of perturbations in rewired bacterial gene networks.

What happens to gene expression when you add new links to a gene regulatory network? To answer this question, we profile 85 network rewirings in E. coli. Here we report that concerted patterns of differential expression propagate from reconnected hub genes. The rewirings link promoter regions to different transcription factor and σ-factor genes, resulting in perturbations that span four orders of magnitude, changing up to ∼ 70% of the transcriptome. Importantly, factor connectivity and promoter activity both associate with perturbation size. Perturbations from related rewirings have more similar transcription profiles and a statistical analysis reveals ∼ 20 underlying states of the system, associating particular gene groups with rewiring constructs. We examine two large clusters (ribosomal and flagellar genes) in detail. These represent alternative global outcomes from different rewirings because of antagonism between these major cell states. This data set of systematically related perturbations enables reverse engineering and discovery of underlying network interactions.

Reconstructing and analysing cellular states, space and time from gene expression profiles of many cells and single cells.

Genome-wide gene expression profiling is a fast, cheap and standardised analysis that provides a high dimensional measurement of the state of a biological sample. In this review we describe computational methods that can be applied to identify and interpret sources of variance in gene expression in whole organisms, organs, tissues or single cells. This allows the identification of constituent cell types and states in complex mixtures, the reconstruction of temporal trajectories of development, differentiation and progression, and the reconstruction of spatial patterning. When applied to genetically variable samples, these methods allow the efficient investigation of how genetic variation influences gene expression and biological processes in space and time.

The DREAM complex promotes gene body H2A.Z for target repression.

The DREAM (DP, Retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell cycle genes, but its mechanism of action is poorly understood. Here we show that Caenorhabditis elegans DREAM targets have an unusual pattern of high gene body HTZ-1/H2A.Z. In mutants of lin-35, the sole p130/Rb-like gene in C. elegans, DREAM targets have reduced gene body HTZ-1/H2A.Z and increased expression. Consistent with a repressive role for gene body H2A.Z, many DREAM targets are up-regulated in htz-1/H2A.Z mutants. Our results indicate that the DREAM complex facilitates high gene body HTZ-1/H2A.Z, which plays a role in target gene repression.

The effects of genetic variation on gene expression dynamics during development.

The development of a multicellular organism and physiological responses require massive coordinated changes in gene expression across several cell and tissue types. Polymorphic regions of the genome that influence gene expression levels have been identified by expression quantitative trait locus (eQTL) mapping in many species, including loci that have cell-dependent, tissue-dependent and age-dependent effects. However, there has been no comprehensive characterization of how polymorphisms influence the complex dynamic patterns of gene expression that occur during development and in physiology. Here we describe an efficient experimental design to infer gene expression dynamics from single expression profiles in different genotypes, and apply it to characterize the effect of local (cis) and distant (trans) genetic variation on gene expression at high temporal resolution throughout a 12-hour period of the development of Caenorhabditis elegans. Taking dynamic variation into account identifies >50% more cis-eQTLs, including more than 900 that alter the dynamics of expression during this period. Local sequence polymorphisms extensively affect the timing, rate, magnitude and shape of expression changes. Indeed, many local sequence variants both increase and decrease gene expression, depending on the time-point profiled. Expression dynamics during this 12-hour period are also influenced extensively in trans by distal loci. In particular, several trans loci influence genes with quite diverse dynamic expression patterns, but they do so primarily during a common time interval. Trans loci can therefore act as modifiers of expression during a particular period of development. This study provides the first characterization, to our knowledge, of the effect of local and distant genetic variation on the dynamics of gene expression throughout an extensive time period. Moreover, the approach developed here should facilitate the genetic dissection of other dynamic processes, including potentially development, physiology and disease progression in humans.

Integrated genome-scale prediction of detrimental mutations in transcription networks.

A central challenge in genetics is to understand when and why mutations alter the phenotype of an organism. The consequences of gene inhibition have been systematically studied and can be predicted reasonably well across a genome. However, many sequence variants important for disease and evolution may alter gene regulation rather than gene function. The consequences of altering a regulatory interaction (or "edge") rather than a gene (or "node") in a network have not been as extensively studied. Here we use an integrative analysis and evolutionary conservation to identify features that predict when the loss of a regulatory interaction is detrimental in the extensively mapped transcription network of budding yeast. Properties such as the strength of an interaction, location and context in a promoter, regulator and target gene importance, and the potential for compensation (redundancy) associate to some extent with interaction importance. Combined, however, these features predict quite well whether the loss of a regulatory interaction is detrimental across many promoters and for many different transcription factors. Thus, despite the potential for regulatory diversity, common principles can be used to understand and predict when changes in regulation are most harmful to an organism.

Complex patterns of gene expression in human T cells during in vivo aging.

Human aging is associated with complex alterations that contribute to remodelling of physiological processes and ultimately manifests in loss of tissue/organ function. Peripheral blood T cells do not escape this phenomenon and undergo profound remodelling with aging. Thus, investigating the effects of aging on T cells transcriptomics and identifying the underlying regulatory mechanisms can be of extreme importance to understand the aging process in the Immune System (IS). To this aim, we performed an analysis of gene expression data of T cells collected from peripheral blood of 25 healthy human donors of different age from 25 to more than 95 years, in order to characterize changes that occur throughout the entire adult lifespan. By means of microarray analysis, we observed large groups of genes exhibiting non-monotonic expression patterns over time: such behaviour, that could not be observed in typical "two-group" experiments (e.g. young vs. old people) highlights similarities in gene expression profiles of young and "successfully aged" individuals. Genes whose expression profiles change during lifespan were grouped into three main patterns (eigenmodes) to which different biological functions were significantly associated. The analysis of KEGG pathways to which these genes belong indicated that the biological processes altered in T cell aging are not only those typically associated with immune cells (Jak-STAT signalling, T cell receptor signalling, cytokine-cytokine receptor interactions, etc.) but also some not specific of immune cells, such as long-term depression, PPAR and mTOR signalling, glucose and glutathione metabolism, suggesting that T cell aging may be representative of a more generalised aging phenomenon. Thus, the T cell may represent a useful cellular model to study organismal aging. We further searched for over-represented transcription factor binding sites (TFBSs) in the promoter regions of genes clustered by similarity of their age-related patterns to evidence possible co-regulation. A comparison between over-representation of TFBSs and the time course of the corresponding transcription factor (TF) expression levels revealed that a restricted group of TFs may play a central role in driving aging-specific changes in gene expression of T cells.

Regulation of gene expression in hepatic cells by the mammalian Target of Rapamycin (mTOR).

BACKGROUND: We investigated mTOR regulation of gene expression by studying rapamycin effect in two hepatic cell lines, the non-tumorigenic WB-F344 cells and the tumorigenic WB311 cells. The latter are resistant to the growth inhibitory effects of rapamycin, thus providing us with an opportunity to study the gene expression effects of rapamycin without confounding effects on cell proliferation. METHODOLOGY/PRINCIPAL FINDINGS: The hepatic cells were exposed to rapamycin for 24 hr. Microarray analysis on total RNA preparations identified genes that were affected by rapamycin in both cell lines and, therefore, modulated independent of growth arrest. Further studies showed that the promoter regions of these genes included E-box-containing transcription factor binding sites at higher than expected rates. Based on this, we tested the hypothesis that c-Myc is involved in regulation of gene expression by mTOR by comparing genes altered by rapamycin in the hepatic cells and by c-Myc induction in fibroblasts engineered to express c-myc in an inducible manner. Results showed enrichment for c-Myc targets among rapamycin sensitive genes in both hepatic cell lines. However, microarray analyses on wild type and c-myc null fibroblasts showed similar rapamycin effect, with the set of rapamycin-sensitive genes being enriched for c-Myc targets in both cases. CONCLUSIONS/SIGNIFICANCE: There is considerable overlap in the regulation of gene expression by mTOR and c-Myc. However, regulation of gene expression through mTOR is c-Myc-independent and cannot be attributed to the involvement of specific transcription factors regulated by the rapamycin-sensitive mTOR Complex 1.

Rapamycin response in tumorigenic and non-tumorigenic hepatic cell lines.

BACKGROUND: The mTOR inhibitor rapamycin has anti-tumor activity across a variety of human cancers, including hepatocellular carcinoma. However, resistance to its growth inhibitory effects is common. We hypothesized that hepatic cell lines with varying rapamycin responsiveness would show common characteristics accounting for resistance to the drug. METHODOLOGY/PRINCIPAL FINDINGS: We profiled a total of 13 cell lines for rapamycin-induced growth inhibition. The non-tumorigenic rat liver epithelial cell line WB-F344 was highly sensitive while the tumorigenic WB311 cell line, originally derived from the WB-F344 line, was highly resistant. The other 11 cell lines showed a wide range of sensitivities. Rapamycin induced inhibition of cyclin E-dependent kinase activity in some cell lines, but the ability to do so did not correlate with sensitivity. Inhibition of cyclin E-dependent kinase activity was related to incorporation of p27(Kip1) into cyclin E-containing complexes in some but not all cell lines. Similarly, sensitivity of global protein synthesis to rapamycin did not correlate with its anti-proliferative effect. However, rapamycin potently inhibited phosphorylation of two key substrates, ribosomal protein S6 and 4E-BP1, in all cases, indicating that the locus of rapamycin resistance was downstream from inhibition of mTOR Complex 1. Microarray analysis did not disclose a unifying mechanism for rapamycin resistance, although the glycolytic pathway was downregulated in all four cell lines studied. CONCLUSIONS/SIGNIFICANCE: We conclude that the mechanisms of rapamycin resistance in hepatic cells involve alterations of signaling downstream from mTOR and that the mechanisms are highly heterogeneous, thus predicting that maintaining or promoting sensitivity will be highly challenging.

Overcoming resistance to conventional drugs in Ewing sarcoma and identification of molecular predictors of outcome.

PURPOSE: The improvement of Ewing sarcoma (EWS) therapy is currently linked to the discovery of strategies to select patients with poor and good prognosis and of modified treatment regimens. In this study, we analyzed the molecular factors governing EWS response to chemotherapy to identify genetic signatures to be used for risk-adapted therapy. PATIENTS AND METHODS: Microarray technology was used for profiling 30 primary tumors and seven metastases of patients who were classified according to event-free survival. For selected genes, real-time polymerase chain reaction was applied in 42 EWS primary tumors as validation assay. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide test was used to evaluate in vitro drug sensitivity. RESULTS: We identified molecular signatures that reflect tumor resistance to chemotherapy. Annotation analysis was applied to reveal the biologic functions that critically influenced clinical outcome. The prognostic relevance of glutathione metabolism pathway was validated. The expression of MGST1, the microsomal glutathione S-transferase (GST), was found to clearly predict EWS prognosis. MGST1 expression was associated with doxorubicin chemosensitivity. This prompted us to assess the in vitro effectiveness of 6-(7-nitro-2,1,3-benzoxadiazol-4-ylthio)hexanol (NBDHEX), a new anticancer agent that efficiently inhibits GST enzymes. Six cell lines were found to be sensitive to this new drug. CONCLUSION: Classification of EWS patients into high- and low-risk groups is feasible with restricted molecular signatures that may have practical value at diagnosis for selecting patients with EWS who are unresponsive to current treatments. Glutathione metabolism pathway emerged as one of the most significantly altered prognosis-associated pathway. NBDHEX is proposed as a new potential therapeutic possibility.

Displayed correlation between gene expression profiles and submicroscopic alterations in response to cetuximab, gefitinib and EGF in human colon cancer cell lines.

BACKGROUND: EGFR is frequently overexpressed in colon cancer. We characterized HT-29 and Caco-2, human colon cancer cell lines, untreated and treated with cetuximab or gefitinib alone and in combination with EGF. METHODS: Cell growth was determined using a variation on the MTT assay. Cell-cycle analysis was conducted by flow cytometry. Immunohistochemistry was performed to evaluate EGFR expression and scanning electron microscopy (SEM) evidenced the ultrastructural morphology. Gene expression profiling was performed using hybridization of the microarray Ocimum Pan Human 40 K array A. RESULTS: Caco-2 and HT-29 were respectively 66.25 and 59.24 % in G0/G1. They maintained this level of cell cycle distribution after treatment, suggesting a predominantly differentiated state. Treatment of Caco-2 with EGF or the two EGFR inhibitors produced a significant reduction in their viability. SEM clearly showed morphological cellular transformations in the direction of cellular death in both cell lines treated with EGFR inhibitors. HT-29 and Caco-2 displayed an important reduction of the microvilli (which also lose their erect position in Caco-2), possibly invalidating microvilli absorption function. HT-29 treated with cetuximab lost their boundary contacts and showed filipodi; when treated with gefitinib, they showed some vesicles: generally membrane reshaping is evident. Both cell lines showed a similar behavior in terms of on/off switched genes upon treatment with cetuximab. The gefitinib global gene expression pattern was different for the 2 cell lines; gefitinib treatment induced more changes, but directly correlated with EGF treatment. In cetuximab or gefitinib plus EGF treatments there was possible summation of the morphological effects: cells seemed more weakly affected by the transformation towards apoptosis. The genes appeared to be less stimulated than for single drug cases. CONCLUSION: This is the first study to have systematically investigated the effect of cetuximab or gefitinib, alone and in combination with EGF, on human colon cancer cell lines. The EGFR inhibitors have a weaker effect in the presence of EGF that binds EGFR. Cetuximab treatment showed an expression pattern that inversely correlates with EGF treatment. We found interesting cyto-morphological features closely relating to gene expression profile. Both drugs have an effect on differentiation towards cellular death.

Reconstructing networks of pathways via significance analysis of their intersections.

BACKGROUND: Significance analysis at single gene level may suffer from the limited number of samples and experimental noise that can severely limit the power of the chosen statistical test. This problem is typically approached by applying post hoc corrections to control the false discovery rate, without taking into account prior biological knowledge. Pathway or gene ontology analysis can provide an alternative way to relax the significance threshold applied to single genes and may lead to a better biological interpretation. RESULTS: Here we propose a new analysis method based on the study of networks of pathways. These networks are reconstructed considering both the significance of single pathways (network nodes) and the intersection between them (links). We apply this method for the reconstruction of networks of pathways to two gene expression datasets: the first one obtained from a c-Myc rat fibroblast cell line expressing a conditional Myc-estrogen receptor oncoprotein; the second one obtained from the comparison of Acute Myeloid Leukemia and Acute Lymphoblastic Leukemia derived from bone marrow samples. CONCLUSION: Our method extends statistical models that have been recently adopted for the significance analysis of functional groups of genes to infer links between these groups. We show that groups of genes at the interface between different pathways can be considered as relevant even if the pathways they belong to are not significant by themselves.