The transcriptional network that controls growth arrest and differentiation in a human myeloid leukemia cell line

Identification of TNF-α-responsive promoters and enhancers in the intestinal epithelial cell model Caco-2

The Caco-2 cell line is one of the most important in vitro models for enterocytes, and is used to study drug absorption and disease, including inflammatory bowel disease and cancer. In order to use the model optimally, it is necessary to map its functional entities. In this study, we have generated genome-wide maps of active transcription start sites (TSSs), and active enhancers in Caco-2 cells with or without tumour necrosis factor (TNF)-α stimulation to mimic an inflammatory state. We found 520 promoters that significantly changed their usage level upon TNF-α stimulation; of these, 52% are not annotated. A subset of these has the potential to confer change in protein function due to protein domain exclusion. Moreover, we locate 890 transcribed enhancer candidates, where ∼50% are changing in usage after TNF-α stimulation. These enhancers share motif enrichments with similarly responding gene promoters. As a case example, we characterize an enhancer regulating the laminin-5 γ2-chain (LAMC2) gene by nuclear factor (NF)-κB binding. This report is the first to present comprehensive TSS and enhancer maps over Caco-2 cells, and highlights many novel inflammation-specific promoters and enhancers.

Personalized medicine in hematology - A landmark from bench to bed

Personalized medicine is the cornerstone of medical practice. It tailors treatments for specific conditions of an affected individual. The borders of personalized medicine are defined by limitations in technology and our understanding of biology, physiology and pathology of various conditions. Current advances in technology have provided physicians with the tools to investigate the molecular makeup of the disease. Translating these molecular make-ups to actionable targets has led to the development of small molecular inhibitors. Also, detailed understanding of genetic makeup has allowed us to develop prognostic markers, better known as companion diagnostics. Current attempts in the development of drug delivery systems offer the opportunity of delivering specific inhibitors to affected cells in an attempt to reduce the unwanted side effects of drugs.

Early B cell factor 1 regulates adipocyte morphology and lipolysis in white adipose tissue

White adipose tissue (WAT) morphology characterized by hypertrophy (i.e., fewer but larger adipocytes) associates with increased adipose inflammation, lipolysis, insulin resistance, and risk of diabetes. However, the causal relationships and the mechanisms controlling WAT morphology are unclear. Herein, we identified EBF1 as an adipocyte-expressed transcription factor with decreased expression/activity in WAT hypertrophy. In human adipocytes, the regulatory targets of EBF1 were enriched for genes controlling lipolysis and adipocyte morphology/differentiation, and in both humans and murine models, reduced EBF1 levels associated with increased lipolysis and adipose hypertrophy. Although EBF1 did not affect adipose inflammation, TNFα reduced EBF1 gene expression. High-fat diet intervention in Ebf1(+/-) mice resulted in more pronounced WAT hypertrophy and attenuated insulin sensitivity compared with wild-type littermate controls. We conclude that EBF1 is an important regulator of adipose morphology and fat cell lipolysis and may constitute a link between WAT inflammation, altered lipid metabolism, adipose hypertrophy, and insulin resistance.

RNA sequencing: from tag-based profiling to resolving complete transcript structure

Technological advances in the sequencing field support in-depth characterization of the transcriptome. Here, we review genome-wide RNA sequencing methods used to investigate specific aspects of gene expression and its regulation, from transcription to RNA processing and translation. We discuss tag-based methods for studying transcription, alternative initiation and polyadenylation events, shotgun methods for detection of alternative splicing, full-length RNA sequencing for the determination of complete transcript structures, and targeted methods for studying the process of transcription and translation. With the ensemble of technologies available, it is now possible to obtain a comprehensive view on transcriptome complexity and the regulation of transcript diversity.

The genome-wide molecular signature of transcription factors in leukemia

Transcription factors control expression of genes essential for the normal functioning of the hematopoietic system and regulate development of distinct blood cell types. During leukemogenesis, aberrant regulation of transcription factors such as RUNX1, CBFβ, MLL, C/EBPα, SPI1, GATA, and TAL1 is central to the disease. Here, we will discuss the mechanisms of transcription factor deregulation in leukemia and how in recent years next-generation sequencing approaches have helped to elucidate the molecular role of many of these aberrantly expressed transcription factors. We will focus on the complexes in which these factors reside, the role of posttranslational modification of these factors, their involvement in setting up higher order chromatin structures, and their influence on the local epigenetic environment. We suggest that only comprehensive knowledge on all these aspects will increase our understanding of aberrant gene expression in leukemia as well as open new entry points for therapeutic intervention.

Transcriptional profiling of the human fibrillin/LTBP gene family, key regulators of mesenchymal cell functions

The fibrillins and latent transforming growth factor binding proteins (LTBPs) form a superfamily of extracellular matrix (ECM) proteins characterized by the presence of a unique domain, the 8-cysteine transforming growth factor beta (TGFβ) binding domain. These proteins are involved in the structure of the extracellular matrix and controlling the bioavailability of TGFβ family members. Genes encoding these proteins show differential expression in mesenchymal cell types which synthesize the extracellular matrix. We have investigated the promoter regions of the seven gene family members using the FANTOM5 CAGE database for human. While the protein and nucleotide sequences show considerable sequence similarity, the promoter regions were quite diverse. Most genes had a single predominant transcription start site region but LTBP1 and LTBP4 had two regions initiating different transcripts. Most of the family members were expressed in a range of mesenchymal and other cell types, often associated with use of alternative promoters or transcription start sites within a promoter in different cell types. FBN3 was the lowest expressed gene, and was found only in embryonic and fetal tissues. The different promoters for one gene were more similar to each other in expression than to promoters of the other family members. Notably expression of all 22 LTBP2 promoters was tightly correlated and quite distinct from all other family members. We located candidate enhancer regions likely to be involved in expression of the genes. Each gene was associated with a unique subset of transcription factors across multiple promoters although several motifs including MAZ, SP1, GTF2I and KLF4 showed overrepresentation across the gene family. FBN1 and FBN2, which had similar expression patterns, were regulated by different transcription factors. This study highlights the role of alternative transcription start sites in regulating the tissue specificity of closely related genes and suggests that this important class of extracellular matrix proteins is subject to subtle regulatory variations that explain the differential roles of members of this gene family.

Patterns of genome-wide VDR locations

The genome-wide analysis of the binding sites of the transcription factor vitamin D receptor (VDR) is essential for a global appreciation the physiological impact of the nuclear hormone 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). Genome-wide analysis of lipopolysaccharide (LPS)-polarized THP-1 human monocytic leukemia cells via chromatin immunoprecipitation sequencing (ChIP-seq) resulted in 1,318 high-confidence VDR binding sites, of which 789 and 364 occurred uniquely with and without 1,25(OH)₂D₃ stimulation, while only 165 were common. We re-analyzed five public VDR ChIP-seq datasets with identical peak calling settings (MACS, version 2) and found, using a novel consensus summit identification strategy, in total 23,409 non-overlapping VDR binding sites, 75% of which are unique within the six analyzed cellular models. LPS-differentiated THP-1 cells have 22% more genomic VDR locations than undifferentiated cells and both cell types display more overlap in their VDR locations than the other investigated cell types. In general, the intersection of VDR binding profiles of ligand-stimulated cells is higher than those of unstimulated cells. De novo binding site searches and HOMER screening for binding motifs formed by direct repeats spaced by three nucleotides (DR3) suggest for all six VDR ChIP-seq datasets that these sequences are found preferentially at highly ligand responsive VDR loci. Importantly, all VDR ChIP-seq datasets display the same relationship between the VDR occupancy and the percentage of DR3-type sequences below the peak summits. The comparative analysis of six VDR ChIP-seq datasets demonstrated that the mechanistic basis for the action of the VDR is independent of the cell type. Only the minority of genome-wide VDR binding sites contains a DR3-type sequence. Moreover, the total number of identified VDR binding sites in each ligand-stimulated cell line inversely correlates with the percentage of peak summits with DR3 sites.

RECLU: a pipeline to discover reproducible transcriptional start sites and their alternative regulation using capped analysis of gene expression (CAGE)

BACKGROUND

Next generation sequencing based technologies are being extensively used to study transcriptomes. Among these, cap analysis of gene expression (CAGE) is specialized in detecting the most 5' ends of RNA molecules. After mapping the sequenced reads back to a reference genome CAGE data highlights the transcriptional start sites (TSSs) and their usage at a single nucleotide resolution.

RESULTS

We propose a pipeline to group the single nucleotide TSS into larger reproducible peaks and compare their usage across biological states. Importantly, our pipeline discovers broad peaks as well as the fine structure of individual transcriptional start sites embedded within them. We assess the performance of our approach on a large CAGE datasets including 156 primary cell types and two cell lines with biological replicas. We demonstrate that genes have complicated structures of transcription initiation events. In particular, we discover that narrow peaks embedded in broader regions of transcriptional activity can be differentially used even if the larger region is not.

CONCLUSIONS

By examining the reproducible fine scaled organization of TSS we can detect many differentially regulated peaks undetected by previous approaches.

Transcription and enhancer profiling in human monocyte subsets

In-depth regulome analysis of human monocyte subsets, including transcription and enhancer profiling. Description of metabolomic differences in human monocyte subsets.

Characterization of genomic vitamin D receptor binding sites through chromatin looping and opening

The vitamin D receptor (VDR) is a transcription factor that mediates the genomic effects of 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). Genome-wide there are several thousand binding sites and hundreds of primary 1,25(OH)₂D₃ target genes, but their functional relation is largely elusive. In this study, we used ChIA-PET data of the transcription factor CTCF in combination with VDR ChIP-seq data, in order to map chromatin domains containing VDR binding sites. In total, we found 1,599 such VDR containing chromatin domains and studied in THP-1 human monocytic leukemia cells four representatives of them. Our combined ChIP-seq and FAIRE-seq time course data showed that each of these four domains contained a master VDR binding site, where an increase of VDR binding pairs with 1,25(OH)₂D₃-promoted chromatin opening and the presence of a highly significant DR3-type sequence below the peak summit. These sites differed in their relative VDR binding but not in their kinetics, while other loci either had a weaker and delayed VDR association or could not be confirmed at all. All studied chromatin domains contained at least one primary 1,25(OH)₂D₃ target gene demonstrating a characteristic slope of mRNA increase, while neighboring genes responded delayed, if at all. In conclusion, the observation of ligand-inducible VDR binding and chromatin opening combined with a DR3-type sequence highlighted genome-wide 160 VDR loci that have within their chromatin domain a more than 4-fold increased likelihood to identify a primary 1,25(OH)₂D₃ target gene than in the vicinity of other genomic VDR binding sites.

Additive effects of microRNAs and transcription factors on CCL2 production in human white adipose tissue

Adipose tissue inflammation is present in insulin-resistant conditions. We recently proposed a network of microRNAs (miRNAs) and transcription factors (TFs) regulating the production of the proinflammatory chemokine (C-C motif) ligand-2 (CCL2) in adipose tissue. We presently extended and further validated this network and investigated if the circuits controlling CCL2 can interact in human adipocytes and macrophages. The updated subnetwork predicted that miR-126/-193b/-92a control CCL2 production by several TFs, including v-ets erythroblastosis virus E26 oncogene homolog 1 (avian) (ETS1), MYC-associated factor X (MAX), and specificity protein 12 (SP1). This was confirmed in human adipocytes by the observation that gene silencing of ETS1, MAX, or SP1 attenuated CCL2 production. Combined gene silencing of ETS1 and MAX resulted in an additive reduction in CCL2 production. Moreover, overexpression of miR-126/-193b/-92a in different pairwise combinations reduced CCL2 secretion more efficiently than either miRNA alone. However, although effects on CCL2 secretion by co-overexpression of miR-92a/-193b and miR-92a/-126 were additive in adipocytes, the combination of miR-126/-193b was primarily additive in macrophages. Signals for miR-92a and -193b converged on the nuclear factor-κB pathway. In conclusion, TF and miRNA-mediated regulation of CCL2 production is additive and partly relayed by cell-specific networks in human adipose tissue that may be important for the development of insulin resistance/type 2 diabetes.

Comparison of CAGE and RNA-seq transcriptome profiling using clonally amplified and single-molecule next-generation sequencing

CAGE (cap analysis gene expression) and RNA-seq are two major technologies used to identify transcript abundances as well as structures. They measure expression by sequencing from either the 5' end of capped molecules (CAGE) or tags randomly distributed along the length of a transcript (RNA-seq). Library protocols for clonally amplified (Illumina, SOLiD, 454 Life Sciences [Roche], Ion Torrent), second-generation sequencing platforms typically employ PCR preamplification prior to clonal amplification, while third-generation, single-molecule sequencers can sequence unamplified libraries. Although these transcriptome profiling platforms have been demonstrated to be individually reproducible, no systematic comparison has been carried out between them. Here we compare CAGE, using both second- and third-generation sequencers, and RNA-seq, using a second-generation sequencer based on a panel of RNA mixtures from two human cell lines to examine power in the discrimination of biological states, detection of differentially expressed genes, linearity of measurements, and quantification reproducibility. We found that the quantified levels of gene expression are largely comparable across platforms and conclude that CAGE and RNA-seq are complementary technologies that can be used to improve incomplete gene models. We also found systematic bias in the second- and third-generation platforms, which is likely due to steps such as linker ligation, cleavage by restriction enzymes, and PCR amplification. This study provides a perspective on the performance of these platforms, which will be a baseline in the design of further experiments to tackle complex transcriptomes uncovered in a wide range of cell types.

A promoter-level mammalian expression atlas

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

Redefinition of the human mast cell transcriptome by deep-CAGE sequencing

Mast cells (MCs) mature exclusively in peripheral tissues, hampering research into their developmental and functional programs. Here, we employed deep cap analysis of gene expression on skin-derived MCs to generate the most comprehensive view of the human MC transcriptome ever reported. An advantage is that MCs were embedded in the FANTOM5 project, giving the opportunity to contrast their molecular signature against a multitude of human samples. We demonstrate that MCs possess a unique and surprising transcriptional landscape, combining hematopoietic genes with those exclusively active in MCs and genes not previously reported as expressed by MCs (several of them markers of unrelated tissues). We also found functional bone morphogenetic protein receptors transducing activatory signals in MCs. Conversely, several immune-related genes frequently studied in MCs were not expressed or were weakly expressed. Comparing MCs ex vivo with cultured counterparts revealed profound changes in the MC transcriptome in in vitro surroundings. We also determined the promoter usage of MC-expressed genes and identified associated motifs active in the lineage. Befitting their uniqueness, MCs had no close relative in the hematopoietic network (also only distantly related with basophils). This rich data set reveals that our knowledge of human MCs is still limited, but with this resource, novel functional programs of MCs may soon be discovered.

High-throughput transcription profiling identifies putative epigenetic regulators of hematopoiesis

Hematopoietic differentiation is governed by a complex regulatory program controlling the generation of different lineages of blood cells from multipotent hematopoietic stem cells. The transcriptional program that dictates hematopoietic cell fate and differentiation requires an epigenetic memory function provided by a network of epigenetic factors regulating DNA methylation, posttranslational histone modifications, and chromatin structure. Aberrant interactions between epigenetic factors and transcription factors cause perturbations in the blood cell differentiation program that result in various types of hematopoietic disorders. To elucidate the contributions of different epigenetic factors in human hematopoiesis, high-throughput cap analysis of gene expression was used to build transcription profiles of 199 epigenetic factors in a wide range of blood cells. Our epigenetic transcriptome analysis revealed cell type- (eg, HELLS and ACTL6A), lineage- (eg, MLL), and/or leukemia- (eg, CHD2, CBX8, and EPC1) specific expression of several epigenetic factors. In addition, we show that several epigenetic factors use alternative transcription start sites in different cell types. This analysis could serve as a resource for the scientific community for further characterization of the role of these epigenetic factors in blood development.

Analysis of the DNA methylome and transcriptome in granulopoiesis reveals timed changes and dynamic enhancer methylation

In development, epigenetic mechanisms such as DNA methylation have been suggested to provide a cellular memory to maintain multipotency but also stabilize cell fate decisions and direct lineage restriction. In this study, we set out to characterize changes in DNA methylation and gene expression during granulopoiesis using 4 distinct cell populations ranging from the oligopotent common myeloid progenitor stage to terminally differentiated neutrophils. We observed that differentially methylated sites (DMSs) generally show decreased methylation during granulopoiesis. Methylation appears to change at specific differentiation stages and overlap with changes in transcription and activity of key hematopoietic transcription factors. DMSs were preferentially located in areas distal to CpG islands and shores. Also, DMSs were overrepresented in enhancer elements and enriched in enhancers that become active during differentiation. Overall, this study depicts in detail the epigenetic and transcriptional changes that occur during granulopoiesis and supports the role of DNA methylation as a regulatory mechanism in blood cell differentiation.

Basonuclin-1 modulates epithelial plasticity and TGF-β1-induced loss of epithelial cell integrity

Transforming growth factor-β1 (TGF-β1) is a multifunctional cytokine and critically involved in the progression of a variety of cancers. TGF-β1 signaling can impair tumor development by its anti-proliferative and pro-apoptotic features. In contrast, it may actively promote tumor progression and cancer cell dissemination by inducing a gradual switch from epithelial towards mesenchymal-like cell features (EMT-like), including decreased intercellular adhesion. Here, we show that expression of the transcription factor Basonuclin-1 (Bnc1) modulates TGF-β1-induced epithelial dedifferentiation of mammary epithelial cells. RNAi-mediated repression of Bnc1 resulted in enhanced intercellular adhesion and strongly impaired TGF-β1-dependent sheet disintegration and cell scattering. In contrast, forced expression of Bnc1 modifies plasma membrane/cytoskeletal dynamics and seemingly interferes with the initiation of sustainable cell-cell contacts. Follow-up analyses revealed that Bnc1 affects the expression of numerous TGF-β1-responsive genes including distinct EMT-related transcription factors, some of which modulate the expression of Bnc1 themselves. These results suggest that Bnc1 is part of a transcription factor network related to epithelial plasticity with reciprocal feedback-loop connections on which Smad-factors integrate TGF-β1 signaling. Our study demonstrates that Bnc1 regulates epithelial plasticity of mammary epithelial cells and influences outcome of TGF-β1 signaling.

Data integration in the era of omics: current and future challenges

To integrate heterogeneous and large omics data constitutes not only a conceptual challenge but a practical hurdle in the daily analysis of omics data. With the rise of novel omics technologies and through large-scale consortia projects, biological systems are being further investigated at an unprecedented scale generating heterogeneous and often large data sets. These data-sets encourage researchers to develop novel data integration methodologies. In this introduction we review the definition and characterize current efforts on data integration in the life sciences. We have used a web-survey to assess current research projects on data-integration to tap into the views, needs and challenges as currently perceived by parts of the research community.

The functional consequences of variation in transcription factor binding

One goal of human genetics is to understand how the information for precise and dynamic gene expression programs is encoded in the genome. The interactions of transcription factors (TFs) with DNA regulatory elements clearly play an important role in determining gene expression outputs, yet the regulatory logic underlying functional transcription factor binding is poorly understood. Many studies have focused on characterizing the genomic locations of TF binding, yet it is unclear to what extent TF binding at any specific locus has functional consequences with respect to gene expression output. To evaluate the context of functional TF binding we knocked down 59 TFs and chromatin modifiers in one HapMap lymphoblastoid cell line. We then identified genes whose expression was affected by the knockdowns. We intersected the gene expression data with transcription factor binding data (based on ChIP-seq and DNase-seq) within 10 kb of the transcription start sites of expressed genes. This combination of data allowed us to infer functional TF binding. Using this approach, we found that only a small subset of genes bound by a factor were differentially expressed following the knockdown of that factor, suggesting that most interactions between TF and chromatin do not result in measurable changes in gene expression levels of putative target genes. We found that functional TF binding is enriched in regulatory elements that harbor a large number of TF binding sites, at sites with predicted higher binding affinity, and at sites that are enriched in genomic regions annotated as “active enhancers.”

An important question in genomics is to understand how a class of proteins called “transcription factors” controls the expression level of other genes in the genome in a cell-type-specific manner – a process that is essential to human development. One major approach to this problem is to study where these transcription factors bind in the genome, but this does not tell us about the effect of that binding on gene expression levels and it is generally accepted that much of the binding does not strongly influence gene expression. To address this issue, we artificially reduced the concentration of 59 different transcription factors in the cell and then examined which genes were impacted by the reduced transcription factor level. Our results implicate some attributes that might influence what binding is functional, but they also suggest that a simple model of functional vs. non-functional binding may not suffice.

Protein synthesis rate is the predominant regulator of protein expression during differentiation

The contribution of transcription, protein synthesis and degradation rates to the control of protein expression during differentiation was analyzed using quantitative proteomics and transcriptomics data. Protein synthesis rate was identified as the main determinant of protein expression.

The lack of correlation usually observed between transcript and protein levels can be fully explained when correcting for the synthesis and degradation rates of the individual proteins.

Synthesis rates for individual proteins are extensively regulated, in contrast to degradation rates that mostly remain constant in response to differentiation.

The modularity of macromolecular complexes is maintained during synthesis and degradation of the complexes.

External perturbations, by forcing cells to adapt to a new environment, often elicit large-scale changes in gene expression resulting in an altered proteome that improves the cell's fitness in the new conditions. Steady-state levels of a proteome depend on transcription, the levels of transcripts, translation and protein degradation but system-level contribution that each of these processes make to the final protein expression change has yet to be explored. We therefore applied a systems biology approach to characterize the regulation of protein expression during cellular differentiation using quantitative proteomics. As a general rule, it seems that protein expression during cellular differentiation is largely controlled by changes in the relative synthesis rate, whereas the relative degradation rate of the majority of proteins stays constant. In these data, we also observe that the proteins in defined sub-structures of larger protein complexes tend to have highly correlated synthesis and degradation rates but that this does not necessarily extend to the holo-complex. Finally, we provide strong evidence that the generally poor correlation observed between transcript and protein levels can fully be explained once the protein synthesis and degradation rates are taken into account.

Pegylated IFN-α regulates hepatic gene expression through transient Jak/STAT activation

The use of pegylated interferon-α (pegIFN-α) has replaced unmodified recombinant IFN-α for the treatment of chronic viral hepatitis. While the superior antiviral efficacy of pegIFN-α is generally attributed to improved pharmacokinetic properties, the pharmacodynamic effects of pegIFN-α in the liver have not been studied. Here, we analyzed pegIFN-α-induced signaling and gene regulation in paired liver biopsies obtained prior to treatment and during the first week following pegIFN-α injection in 18 patients with chronic hepatitis C. Despite sustained high concentrations of pegIFN-α in serum, the Jak/STAT pathway was activated in hepatocytes only on the first day after pegIFN-α administration. Evaluation of liver biopsies revealed that pegIFN-α induces hundreds of genes that can be classified into four clusters based on different temporal expression profiles. In all clusters, gene transcription was mainly driven by IFN-stimulated gene factor 3 (ISGF3). Compared with conventional IFN-α therapy, pegIFN-α induced a broader spectrum of gene expression, including many genes involved in cellular immunity. IFN-induced secondary transcription factors did not result in additional waves of gene expression. Our data indicate that the superior antiviral efficacy of pegIFN-α is not the result of prolonged Jak/STAT pathway activation in hepatocytes, but rather is due to induction of additional genes that are involved in cellular immune responses.

Transcription factor binding sites prediction based on modified nucleosomes

In computational methods, position weight matrices (PWMs) are commonly applied for transcription factor binding site (TFBS) prediction. Although these matrices are more accurate than simple consensus sequences to predict actual binding sites, they usually produce a large number of false positive (FP) predictions and so are impoverished sources of information. Several studies have employed additional sources of information such as sequence conservation or the vicinity to transcription start sites to distinguish true binding regions from random ones. Recently, the spatial distribution of modified nucleosomes has been shown to be associated with different promoter architectures. These aligned patterns can facilitate DNA accessibility for transcription factors. We hypothesize that using data from these aligned and periodic patterns can improve the performance of binding region prediction. In this study, we propose two effective features, “modified nucleosomes neighboring” and “modified nucleosomes occupancy”, to decrease FP in binding site discovery. Based on these features, we designed a logistic regression classifier which estimates the probability of a region as a TFBS. Our model learned each feature based on Sp1 binding sites on Chromosome 1 and was tested on the other chromosomes in human CD4+T cells. In this work, we investigated 21 histone modifications and found that only 8 out of 21 marks are strongly correlated with transcription factor binding regions. To prove that these features are not specific to Sp1, we combined the logistic regression classifier with the PWM, and created a new model to search TFBSs on the genome. We tested the model using transcription factors MAZ, PU.1 and ELF1 and compared the results to those using only the PWM. The results show that our model can predict Transcription factor binding regions more successfully. The relative simplicity of the model and capability of integrating other features make it a superior method for TFBS prediction.

Monocyte subpopulations in angiogenesis

Growing understanding of the role of the tumor microenvironment in angiogenesis has brought monocyte-derived cells into focus. Monocyte subpopulations are an increasingly attractive therapeutic target in many pathologic states, including cancer. Before monocyte-directed therapies can be fully harnessed for clinical use, understanding of monocyte-driven angiogenesis in tissue development and homeostasis, as well as malignancy, is required. Here, we provide an overview of the mechanisms by which monocytic subpopulations contribute to angiogenesis in tissue and tumor development, highlight gaps in our existing knowledge, and discuss opportunities to exploit these cells for clinical benefit.

Tead2 expression levels control the subcellular distribution of Yap and Taz, zyxin expression and epithelial-mesenchymal transition

The cellular changes during an epithelial-mesenchymal transition (EMT) largely rely on global changes in gene expression orchestrated by transcription factors. Tead transcription factors and their transcriptional co-activators Yap and Taz have been previously implicated in promoting an EMT; however, their direct transcriptional target genes and their functional role during EMT have remained elusive. We have uncovered a previously unanticipated role of the transcription factor Tead2 during EMT. During EMT in mammary gland epithelial cells and breast cancer cells, levels of Tead2 increase in the nucleus of cells, thereby directing a predominant nuclear localization of its co-factors Yap and Taz via the formation of Tead2-Yap-Taz complexes. Genome-wide chromatin immunoprecipitation and next generation sequencing in combination with gene expression profiling revealed the transcriptional targets of Tead2 during EMT. Among these, zyxin contributes to the migratory and invasive phenotype evoked by Tead2. The results demonstrate that Tead transcription factors are crucial regulators of the cellular distribution of Yap and Taz, and together they control the expression of genes critical for EMT and metastasis.

ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

Heterozygous loss-of-function mutations in YAP1 cause both isolated and syndromic optic fissure closure defects

Exome sequence analysis of affected individuals from two families with autosomal-dominant inheritance of coloboma identified two different cosegregating heterozygous nonsense mutations (c.370C>T [p.Arg124*] and c. 1066G>T [p.Glu356*]) in YAP1. The phenotypes of the affected families differed in that one included no extraocular features and the other manifested with highly variable multisystem involvement, including hearing loss, intellectual disability, hematuria, and orofacial clefting. A combined LOD score of 4.2 was obtained for the association between YAP1 loss-of-function mutations and the phenotype in these families. YAP1 encodes an effector of the HIPPO-pathway-induced growth response, and whole-mount in situ hybridization in mouse embryos has shown that Yap1 is strongly expressed in the eye, brain, and fusing facial processes. RT-PCR showed that an alternative transcription start site (TSS) in intron 1 of YAP1 and Yap1 is widely used in human and mouse development, respectively. Transcripts from the alternative TSS are predicted to initiate at codon Met179 relative to the canonical transcript (RefSeq NM_001130145). In these alternative transcripts, the c.370C>T mutation in family 1305 is within the 5' UTR and cannot result in nonsense-mediated decay (NMD). The c. 1066G>T mutation in family 132 should result in NMD in transcripts from either TSS. Amelioration of the phenotype by the alternative transcripts provides a plausible explanation for the phenotypic differences between the families.

Hox in motion: tracking HoxA cluster conformation during differentiation

Three-dimensional genome organization is an important higher order transcription regulation mechanism that can be studied with the chromosome conformation capture techniques. Here, we combined chromatin organization analysis by chromosome conformation capture-carbon copy, computational modeling and epigenomics to achieve the first integrated view, through time, of a connection between chromatin state and its architecture. We used this approach to examine the chromatin dynamics of the HoxA cluster in a human myeloid leukemia cell line at various stages of differentiation. We found that cellular differentiation involves a transient activation of the 5'-end HoxA genes coinciding with a loss of contacts throughout the cluster, and by specific silencing at the 3'-end with H3K27 methylation. The 3D modeling of the data revealed an extensive reorganization of the cluster between the two previously reported topologically associated domains in differentiated cells. Our results support a model whereby silencing by polycomb group proteins and reconfiguration of CTCF interactions at a topologically associated domain boundary participate in changing the HoxA cluster topology, which compartmentalizes the genes following differentiation.

Computational modeling identifies key gene regulatory interactions underlying phenobarbital-mediated tumor promotion

Gene regulatory interactions underlying the early stages of non-genotoxic carcinogenesis are poorly understood. Here, we have identified key candidate regulators of phenobarbital (PB)-mediated mouse liver tumorigenesis, a well-characterized model of non-genotoxic carcinogenesis, by applying a new computational modeling approach to a comprehensive collection of in vivo gene expression studies. We have combined our previously developed motif activity response analysis (MARA), which models gene expression patterns in terms of computationally predicted transcription factor binding sites with singular value decomposition (SVD) of the inferred motif activities, to disentangle the roles that different transcriptional regulators play in specific biological pathways of tumor promotion. Furthermore, transgenic mouse models enabled us to identify which of these regulatory activities was downstream of constitutive androstane receptor and β-catenin signaling, both crucial components of PB-mediated liver tumorigenesis. We propose novel roles for E2F and ZFP161 in PB-mediated hepatocyte proliferation and suggest that PB-mediated suppression of ESR1 activity contributes to the development of a tumor-prone environment. Our study shows that combining MARA with SVD allows for automated identification of independent transcription regulatory programs within a complex in vivo tissue environment and provides novel mechanistic insights into PB-mediated hepatocarcinogenesis.

S5a binds death receptor-6 to induce THP-1 monocytes differentiation via NF-κB pathway

Analyses of apoptotic cell supernatants have helped identify many signals that modulate the states of activation and differentiation in the congeneric or other cells. However, the current knowledge about these soluble factors that are released during apoptosis is rather limited. Previous studies have shown that S5a/Angiocidin induced human acute monocytic leukemia cells (THP-1 cells) to differentiation into macrophages, but the cell surface receptor of S5a has not been identified. In this study we show that apoptotic THP-1 cells released endogenous S5a, and S5a bound with death receptor-6, which was identified as an orphan receptor, to induce THP-1 cells differentiation. Furthermore, we found NF-κB pathway was activated and the transcription factor WT1 and c-myb mediated THP-1 differentiation induced by S5a. And we also show that the differentiation was blocked after anti-DR6 antibody, DR6 siRNA, DR6-Fc, NF-κB inhibitor, or WT1 siRNA treatment. Our finding indicated that the interaction between cells can determine their destination. And we provided evidence for a functional interaction between S5a and DR6, which provides a novel target that can induce the differentiation of cancer cells especially for biotherapy of leukemia.

Detecting expressed genes using CAGE

Cap analysis of gene expression (CAGE) provides accurate high-throughput measurement of RNA expression. By the large-scale analysis of 5' end of transcripts using CAGE method, it enables not only determination of the transcription start site but also prediction of promoter region. Here we provide a protocol for the construction of no-amplification non-tagging CAGE libraries for Illumina next-generation sequencers (nAnT-iCAGE). We have excluded the commonly used PCR amplification and cleavage of restriction enzyme to eliminate any potential biases. As a result, we achieved less biased simple preparation process.

Chromatin signatures at transcriptional start sites separate two equally populated yet distinct classes of intergenic long noncoding RNAs

BACKGROUND

Mammalian transcriptomes contain thousands of long noncoding RNAs (lncRNAs). Some lncRNAs originate from intragenic enhancers which, when active, behave as alternative promoters producing transcripts that are processed using the canonical signals of their host gene. We have followed up this observation by analyzing intergenic lncRNAs to determine the extent to which they might also originate from intergenic enhancers.

RESULTS

We integrated high-resolution maps of transcriptional initiation and transcription to annotate a conservative set of intergenic lncRNAs expressed in mouse erythroblasts. We subclassified intergenic lncRNAs according to chromatin status at transcriptional initiation regions, defined by relative levels of histone H3K4 mono- and trimethylation. These transcripts are almost evenly divided between those arising from enhancer-associated (elncRNA) or promoter-associated (plncRNA) elements. These two classes of 5' capped and polyadenylated RNA transcripts are indistinguishable with regard to their length, number of exons or transcriptional orientation relative to their closest neighboring gene. Nevertheless, elncRNAs are more tissue-restricted, less highly expressed and less well conserved during evolution. Of considerable interest, we found that expression of elncRNAs, but not plncRNAs, is associated with enhanced expression of neighboring protein-coding genes during erythropoiesis.

CONCLUSIONS

We have determined globally the sites of initiation of intergenic lncRNAs in erythroid cells, allowing us to distinguish two similarly abundant classes of transcripts. Different correlations between the levels of elncRNAs, plncRNAs and expression of neighboring genes suggest that functional lncRNAs from the two classes may play contrasting roles in regulating the transcript abundance of local or distal loci.

Understanding the transcriptional regulation of cervix cancer using microarray gene expression data and promoter sequence analysis of a curated gene set

Cervical cancer, the malignant neoplasm of the cervix uteri is the second most common cancer among women worldwide and the top-most cancer in India. Several factors are responsible for causing cervical cancer, which alter the expression of oncogenic genes resulting in up or down-regulation of gene expression and inactivation of tumor-suppressor genes/gene products. Gene expression is regulated by interactions between transcription factors (TFs) and specific regulatory elements in the promoter regions of target genes. Thus, it is important to decipher and analyze TFs that bind to regulatory regions of diseased genes and regulate their expression. In the present study, computational methods involving the combination of gene expression data from microarray experiments and promoter sequence analysis of a curated gene set involved in the cervical cancer causation have been utilized for identifying potential regulatory elements. Consensus predictions of two approaches led to the identification of twelve TFs that might be crucial to the regulation of cervical cancer progression. Subsequently, TF enrichment and oncomine expression analysis suggested that the transcription factor family E2F played an important role for the regulation of genes involve in cervical carcinogenesis. Our results suggest that E2F possesses diagnostic/prognostic value and can act as a potential drug target in cervical cancer.

Changes in the expression of miR-381 and miR-495 are inversely associated with the expression of the MDR1 gene and development of multi-drug resistance

Multidrug resistance (MDR) frequently develops in cancer patients exposed to chemotherapeutic agents and is usually brought about by over-expression of P-glycoprotein (P-gp) which acts as a drug efflux pump to reduce the intracellular concentration of the drug(s). Thus, inhibiting P-gp expression might assist in overcoming MDR in cancer chemotherapy. MiRNAome profiling using next-generation sequencing identified differentially expressed microRNAs (miRs) between parental K562 cells and MDR K562 cells (K562/ADM) induced by adriamycin treatment. Two miRs, miR-381 and miR-495, that were strongly down-regulated in K562/ADM cells, are validated to target the 3'-UTR of the MDR1 gene. These miRs are located within a miR cluster located at chromosome region 14q32.31, and all miRs in this cluster appear to be down-regulated in K562/ADM cells. Functional analysis indicated that restoring expression of miR-381 or miR-495 in K562/ADM cells was correlated with reduced expression of the MDR1 gene and its protein product, P-gp, and increased drug uptake by the cells. Thus, we have demonstrated that changing the levels of certain miR species modulates the MDR phenotype in leukemia cells, and propose further exploration of the use of miR-based therapies to overcome MDR.

Means to an end: mechanisms of alternative polyadenylation of messenger RNA precursors

Expression of mature messenger RNAs (mRNAs) requires appropriate transcription initiation and termination, as well as pre-mRNA processing by capping, splicing, cleavage, and polyadenylation. A core 3'-end processing complex carries out the cleavage and polyadenylation reactions, but many proteins have been implicated in the selection of polyadenylation sites among the multiple alternatives that eukaryotic genes typically have. In recent years, high-throughput approaches to map both the 3'-end processing sites as well as the binding sites of proteins that are involved in the selection of cleavage sites and in the processing reactions have been developed. Here, we review these approaches as well as the insights into the mechanisms of polyadenylation that emerged from genome-wide studies of polyadenylation across a range of cell types and states.

Dynamic modelling of microRNA regulation during mesenchymal stem cell differentiation

Background

Network inference from gene expression data is a typical approach to reconstruct gene regulatory networks. During chondrogenic differentiation of human mesenchymal stem cells (hMSCs), a complex transcriptional network is active and regulates the temporal differentiation progress. As modulators of transcriptional regulation, microRNAs (miRNAs) play a critical role in stem cell differentiation. Integrated network inference aimes at determining interrelations between miRNAs and mRNAs on the basis of expression data as well as miRNA target predictions. We applied the NetGenerator tool in order to infer an integrated gene regulatory network.

Results

Time series experiments were performed to measure mRNA and miRNA abundances of TGF-beta1+BMP2 stimulated hMSCs. Network nodes were identified by analysing temporal expression changes, miRNA target gene predictions, time series correlation and literature knowledge. Network inference was performed using NetGenerator to reconstruct a dynamical regulatory model based on the measured data and prior knowledge. The resulting model is robust against noise and shows an optimal trade-off between fitting precision and inclusion of prior knowledge. It predicts the influence of miRNAs on the expression of chondrogenic marker genes and therefore proposes novel regulatory relations in differentiation control. By analysing the inferred network, we identified a previously unknown regulatory effect of miR-524-5p on the expression of the transcription factor SOX9 and the chondrogenic marker genes COL2A1, ACAN and COL10A1.

Conclusions

Genome-wide exploration of miRNA-mRNA regulatory relationships is a reasonable approach to identify miRNAs which have so far not been associated with the investigated differentiation process. The NetGenerator tool is able to identify valid gene regulatory networks on the basis of miRNA and mRNA time series data.

Proteome changes induced by c-myb silencing in human chronic myeloid leukemia cells suggest molecular mechanisms and putative biomarkers of hematopoietic malignancies

To shed light on the molecular mechanisms associated with aberrant accumulation of c-Myb in chronic myeloid leukemia, comparative proteomic analysis was performed on c-myb RNAi-specifically silenced K562 cells, sampled on a time-course basis. 2D-DIGE technology highlighted 37 differentially-represented proteins that were further characterized by nLC-ESI-LIT-MS/MS and validated by western blotting and qRT-PCR analysis. Most of the deregulated proteins were related to protein folding, energy/primary metabolism, transcription/translation regulation and oxidative stress response. Protein network analysis suggested that glycolysis, gluconeogenesis and protein ubiquitination biosynthesis pathways were highly represented, confirming also the pivotal role of c-Myc. A specific reduced representation was observed for glyceraldehyde-3-phosphate-dehydrogenase and α-enolase, suggesting a possible role of c-Myb in the activation of aerobic glycolysis. A reduced amount was also observed for stress responsive heat shock 70kDa protein and 78kDa glucose-regulated protein, previously identified as direct targets of c-Myb. Among over-represented proteins, worth mentioning is the chromatin modifier chromobox protein homolog 3 that contributes to silencing of E2F- and Myc-responsive genes in quiescent G0 cells. Data here presented, while providing novel insights onto the molecular mechanisms underlying c-Myb activity, indicate potential protein biomarkers for monitoring the progression of chronic myeloid leukemia.

BIOLOGICAL SIGNIFICANCE

Myeloid leukemia is a malignant disease of the hematopoietic system in which cells of myeloid lineages accumulate to an undifferentiated state. In particular, it was shown that an aberrant accumulation of the c-Myb transcriptional factor is associated with the suppression of normal differentiation processes promoting the development of the hematopoietic malignancies. Many efforts have been recently made to identify novel genes directly targeted by c-Myb at a transcriptome level. In this work, we originally describe a differential proteomic approach to facilitate the comprehension of the regulation of the protein networks exerted by c-Myb. Our study reveals a complex network of proteins regulated by c-Myb. The functional heterogeneity of these proteins emphasizes the pleiotropic role of c-Myb as a regulator of genes that are crucial for energy production and stress response in leukemia. In fact, variations in glyceraldehyde-3-phosphate-dehydrogenase and α-enolase suggest a possible role of c-Myb in the activation of aerobic glycolysis. Moreover, significant differences were found for heat shock 70kDa protein and 78kDa glucose-regulated protein known as direct c-Myb targets. This work highlights potential protein biomarkers to look into disease progression and to develop translational medicine approaches in myeloid leukemia.

Identification of direction in gene networks from expression and methylation

BACKGROUND

Reverse-engineering gene regulatory networks from expression data is difficult, especially without temporal measurements or interventional experiments. In particular, the causal direction of an edge is generally not statistically identifiable, i.e., cannot be inferred as a statistical parameter, even from an unlimited amount of non-time series observational mRNA expression data. Some additional evidence is required and high-throughput methylation data can viewed as a natural multifactorial gene perturbation experiment.

RESULTS

We introduce IDEM (Identifying Direction from Expression and Methylation), a method for identifying the causal direction of edges by combining DNA methylation and mRNA transcription data. We describe the circumstances under which edge directions become identifiable and experiments with both real and synthetic data demonstrate that the accuracy of IDEM for inferring both edge placement and edge direction in gene regulatory networks is significantly improved relative to other methods.

CONCLUSION

Reverse-engineering directed gene regulatory networks from static observational data becomes feasible by exploiting the context provided by high-throughput DNA methylation data.An implementation of the algorithm described is available at http://code.google.com/p/idem/.

Functional transcriptomics in the post-ENCODE era

The last decade has seen tremendous effort committed to the annotation of the human genome sequence, most notably perhaps in the form of the ENCODE project. One of the major findings of ENCODE, and other genome analysis projects, is that the human transcriptome is far larger and more complex than previously thought. This complexity manifests, for example, as alternative splicing within protein-coding genes, as well as in the discovery of thousands of long noncoding RNAs. It is also possible that significant numbers of human transcripts have not yet been described by annotation projects, while existing transcript models are frequently incomplete. The question as to what proportion of this complexity is truly functional remains open, however, and this ambiguity presents a serious challenge to genome scientists. In this article, we will discuss the current state of human transcriptome annotation, drawing on our experience gained in generating the GENCODE gene annotation set. We highlight the gaps in our knowledge of transcript functionality that remain, and consider the potential computational and experimental strategies that can be used to help close them. We propose that an understanding of the true overlap between transcriptional complexity and functionality will not be gained in the short term. However, significant steps toward obtaining this knowledge can now be taken by using an integrated strategy, combining all of the experimental resources at our disposal.

The MYB proto-oncogene suppresses monocytic differentiation of acute myeloid leukemia cells via transcriptional activation of its target gene GFI1

The MYB gene is a master regulator of hematopoiesis and contributes to leukemogenesis in several species including humans. Although it is clear that MYB can promote proliferation, suppress apoptosis and block differentiation, the identities of the MYB target genes that mediate these effects have only been partially elucidated. Several studies, including our own, have collectively identified substantial numbers of MYB target genes, including candidates for each of these activities; however, functional validation, particularly in the case of differentiation suppression, has lagged well behind. Here we show that GFI1, which encodes an important regulator of hematopoietic stem cell (HSC) function and granulocytic differentiation, is a direct target of MYB in myeloid leukemia cells. Chromatin immunoprecipitation and reporter studies identified a functional MYB-binding site in the promoter region of GFI, whereas ectopic expression and small hairpin RNA-mediated knockdown of MYB resulted in concomitant increases and decreases, respectively, in GFI1 expression. We also demonstrate that GFI1, like MYB, can block the induced monocytic differentiation of a human acute myeloid leukemia cell line, and most importantly, that GFI1 is essential for MYB's ability to block monocytic differentiation. Thus, we have identified a target of MYB that is a likely mediator of its myeloid differentiation-blocking activity, and which may also be involved in MYB's activities in regulating normal HSC function and myeloid differentiation.

Dynamic analysis of gene expression and genome-wide transcription factor binding during lineage specification of multipotent progenitors

We used the paradigmatic GATA-PU.1 axis to explore, at the systems level, dynamic relationships between transcription factor (TF) binding and global gene expression programs as multipotent cells differentiate. We combined global ChIP-seq of GATA1, GATA2, and PU.1 with expression profiling during differentiation to erythroid and neutrophil lineages. Our analysis reveals (1) differential complexity of sequence motifs bound by GATA1, GATA2, and PU.1; (2) the scope and interplay of GATA1 and GATA2 programs within, and during transitions between, different cell compartments, and the extent of their hard-wiring by DNA motifs; (3) the potential to predict gene expression trajectories based on global associations between TF-binding data and target gene expression; and (4) how dynamic modeling of DNA-binding and gene expression data can be used to infer regulatory logic of TF circuitry. This rubric exemplifies the utility of this cross-platform resource for deconvoluting the complexity of transcriptional programs controlling stem/progenitor cell fate in hematopoiesis.

•

Cross-platform resource for TF-network regulation of multipotent blood cell fate

•

DNA motif dependence and changing specificity of GATA factors in lineage choice

•

Modeling-based inference identifies GATA2 repression of PU.1 in multipotent cells

•

Priming, recruitment, and switching modes of GATA interplay during differentiation

Comparison of RNA- or LNA-hybrid oligonucleotides in template-switching reactions for high-speed sequencing library preparation

Background

Analyzing the RNA pool or transcription start sites requires effective means to convert RNA into cDNA libraries for digital expression counting. With current high-speed sequencers, it is necessary to flank the cDNAs with specific adapters. Adding template-switching oligonucleotides to reverse transcription reactions is the most commonly used approach when working with very small quantities of RNA even from single cells.

Results

Here we compared the performance of DNA-RNA, DNA-LNA and DNA oligonucleotides in template-switching during nanoCAGE library preparation. Test libraries from rat muscle and HeLa cell RNA were prepared in technical triplicates and sequenced for comparison of the gene coverage and distribution of the reads within transcripts. The DNA-RNA oligonucleotide showed the highest specificity for capped 5′ ends of mRNA, whereas the DNA-LNA provided similar gene coverage with more reads falling within exons.

Conclusions

While confirming the cap-specific preference of DNA-RNA oligonucleotides in template-switching reactions, our data indicate that DNA-LNA hybrid oligonucleotides could potentially find other applications in random RNA sequencing.

Polo-like kinase 4 transcription is activated via CRE and NRF1 elements, repressed by DREAM through CDE/CHR sites and deregulated by HPV E7 protein

Infection by oncogenic viruses is a frequent cause for tumor formation as observed in cervical cancer. Viral oncoproteins cause inactivation of p53 function and false transcriptional regulation of central cell cycle genes. Here we analyze the regulation of Plk4, serving as an example of many cell cycle- and p53-regulated genes. Cell cycle genes are often repressed via CDE and CHR elements in their promoters and activated by NF-Y binding to CCAAT-boxes. In contrast, general activation of Plk4 depends on NRF1 and CRE sites. Bioinformatic analyses imply that NRF1 and CRE are central elements of the transcriptional network controlling cell cycle genes. We identify CDE and CHR sites in the Plk4 promoter, which are necessary for binding of the DREAM (DP, RB-like, E2F4 and MuvB) complex and for mediating repression in G₀/G₁. When cells progress to G₂ and mitosis, DREAM is replaced by the MMB (Myb-MuvB) complex that only requires the CHR element for binding. Plk4 expression is downregulated by the p53-p21^WAF1/CIP1-DREAM signaling pathway through the CDE and CHR sites. Cell cycle- and p53-dependent repression is abrogated by HPV E7 oncoprotein. Together with genome-wide analyses our results imply that many cell cycle genes upregulated in tumors by viral infection are bound by DREAM through CDE/CHR sites.

Discerning the subfibrillar structure of mineralized collagen fibrils: a model for the ultrastructure of bone

Biomineralization templated by organic molecules to produce inorganic-organic nanocomposites is a fascinating example of nature using bottom-up strategies at nanoscale to accomplish highly ordered multifunctional materials. One such nanocomposite is bone, composed primarily of hydroxyapatite (HA) nanocrystals that are embedded within collagen fibrils with their c-axes arranged roughly parallel to the long axis of the fibrils. Here we discern the ultra-structure of biomimetic mineralized collagen fibrils (MCFs) as consisting of bundles of subfibrils with approximately 10 nm diameter; each one with an organic-inorganic core-shell structure. Through an amorphous calcium phosphate precursor phase the HA nanocrystals were specifically grown along the longitudinal direction of the collagen microfibrils and encapsulated them within the crystal lattice. They intercalated throughout the collagen fibrils such that the mineral phase surrounded the surface of collagen microfibrils forming an interdigitated network. It appears that this arrangement of collagen microfibrils in collagen fibrils is responsible for the observed ultrastructure. Such a subfibrillar nanostructure in MCFs was identified in both synthetic and natural bone, suggesting this is the basic building block of collagen-based hard tissues. Insights into the ultrastructure of mineralized collagen fibrils have the potential to advance our understanding on the biomineralization principles and the relationship between bone’s structure and mechanical properties, including fracture toughness mechanisms. We anticipate that these principles from biological systems can be applied to the rational design of new nanocomposites with improved performance.

Sox4 is a master regulator of epithelial-mesenchymal transition by controlling Ezh2 expression and epigenetic reprogramming

Gene expression profiling has uncovered the transcription factor Sox4 with upregulated activity during TGF-β-induced epithelial-mesenchymal transition (EMT) in normal and cancerous breast epithelial cells. Sox4 is indispensable for EMT and cell survival in vitro and for primary tumor growth and metastasis in vivo. Among several EMT-relevant genes, Sox4 directly regulates the expression of Ezh2, encoding the Polycomb group histone methyltransferase that trimethylates histone 3 lysine 27 (H3K27me3) for gene repression. Ablation of Ezh2 expression prevents EMT, whereas forced expression of Ezh2 restores EMT in Sox4-deficient cells. Ezh2-mediated H3K27me3 marks associate with key EMT genes, representing an epigenetic EMT signature that predicts patient survival. Our results identify Sox4 as a master regulator of EMT by governing the expression of the epigenetic modifier Ezh2.

Analysis and prediction of pathways in HeLa cells by integrating biological levels of organization with systems-biology approaches

It has recently begun to be considered that cancer is a systemic disease and that it must be studied at every level of complexity using many of the currently available approaches, including high-throughput technologies and bioinformatics. To achieve such understanding in cervical cancer, we collected information on gene, protein and phosphoprotein expression of the HeLa cell line and performed a comprehensive analysis of the different signaling pathways, transcription networks and metabolic events in which they participate. A total expression analysis by RNA-Seq of the HeLa cell line showed that 19,974 genes were transcribed. Of these, 3,360 were over-expressed, and 2,129 under-expressed when compared to the NHEK cell line. A protein-protein interaction network was derived from the over-expressed genes and used to identify central elements and, together with the analysis of over-represented transcription factor motifs, to predict active signaling and regulatory pathways. This was further validated by Metal-Oxide Affinity Chromatography (MOAC) and Tandem Mass Spectrometry (MS/MS) assays which retrieved phosphorylated proteins. The 14-3-3 family members emerge as important regulators in carcinogenesis and as possible clinical targets. We observed that the different over- and under-regulated pathways in cervical cancer could be interrelated through elements that participate in crosstalks, therefore belong to what we term "meta-pathways". Additionally, we highlighted the relations of each one of the differentially represented pathways to one or more of the ten hallmarks of cancer. These features could be maintained in many other types of cancer, regardless of mutations or genomic rearrangements, and favor their robustness, adaptations and the evasion of tissue control. Probably, this could explain why cancer cells are not eliminated by selective pressure and why therapy trials directed against molecular targets are not as effective as expected.