Stochastic responses are not left to pure "chance".

Many coregulated genes assemble in multigene complexes via stochastic inter- and intrachromosomal interactions. In this issue, Fanucchi et al. report that chromatin loop formation governs hierarchical cotranscription within a multigene complex.

Identification of a dynamic gene regulatory network required for pluripotency factor-induced reprogramming of mouse fibroblasts and hepatocytes.

The generation of induced pluripotent stem cells (iPSCs) from somatic cells provides an excellent model to study mechanisms of transcription factor-induced global alterations of the epigenome and genome function. Here, we have investigated the early transcriptional events of cellular reprogramming triggered by the co-expression of Oct4, Sox2, Klf4, and c-Myc (OSKM) in mouse embryonic fibroblasts (MEFs) and mouse hepatocytes (mHeps). In this analysis, we identified a gene regulatory network composed of nine transcriptional regulators (9TR; Cbfa2t3, Gli2, Irf6, Nanog, Ovol1, Rcan1, Taf1c, Tead4, and Tfap4), which are directly targeted by OSKM, in vivo. Functional studies using single and double shRNA knockdowns of any of these factors caused disruption of the network and dramatic reductions in reprogramming efficiency, indicating that this network is essential for the induction and establishment of pluripotency. We demonstrate that the stochastic co-expression of 9TR network components occurs in a remarkably small number of cells, approximating the percentage of terminally reprogrammed cells as a result of dynamic molecular events. Thus, the early DNA-binding patterns of OSKM and the subsequent probabilistic co-expression of essential 9TR components in subpopulations of cells undergoing reprogramming steer the reconstruction of a gene regulatory network marking the transition to pluripotency.

Decoding a gene expression program that accompanies the phenotype of sporadic and basal cell nevus syndrome-associated odontogenic keratocyst.

BACKGROUND: Odontogenic keratocyst is characterized by local aggressive behavior and a high recurrence rate, as well as its potential to develop in association with the basal cell nevus syndrome. The aim of this study was to decode the gene expression program accompanying odontogenic keratocyst phenotype. METHODS: 150-bp paired-end RNA-sequencing was applied on six sporadic and six basal cell nevus syndrome-associated whole-tissue odontogenic keratocyst samples in comparison to six dental follicles, coupled with bioinformatics and complemented by immunohistochemistry. RESULTS: 2654 and 2427 differentially expressed genes were captured to characterize the transcriptome of sporadic and basal cell nevus syndrome-associated odontogenic keratocysts, respectively. Gene ontologies related to "epidermis/skin development" and "keratinocyte/epidermal cell differentiation" were enriched among the upregulated genes (KRT10, NCCRP1, TP63, GRHL3, SOX21), while "extracellular matrix organization" (ITGA5, LOXL2) and "odontogenesis" (MSX1, LHX8) gene ontologies were overrepresented among the downregulated genes in odontogenic keratocyst. Interestingly, upregulation of various embryonic stem cells markers (EPHA1, SCNN1A) and genes committed in cellular reprogramming (SOX2, KLF4, OVOL1, IRF6, TACSTD2, CDH1) was found in odontogenic keratocyst. These findings were highly shared between sporadic and basal cell nevus syndrome-associated odontogenic keratocysts. Immunohistochemistry verified SOX2, KLF4, OVOL1, IRF6, TACSTD2/TROP2, CDH1/E-cadherin, and p63 expression predominantly in the odontogenic keratocyst suprabasal epithelial layers. CONCLUSION: The odontogenic keratocyst transcriptomic profile is characterized by a prominent epidermal and dental epithelial fate, a repressed dental mesenchyme fate combined with deregulated extracellular matrix organization, and enhanced stemness gene signatures. Thus, we propose a developed epidermis-like phenotype in the odontogenic keratocyst suprabasal epithelial cells, established in parallel to a significant upregulation of marker genes related to embryonic stem cells and cellular reprogramming.

The immunohistochemical profile of basal cell nevus syndrome-associated and sporadic odontogenic keratocysts: a systematic review and meta-analysis.

OBJECTIVES: To provide a systematic review of the literature on studies comparing the immunoprofile of nevoid basal cell carcinoma syndrome (BCNS)-associated and sporadic odontogenic keratocysts (OKCs), in order to identify markers that could accurately distinguish the two OKC subtypes. MATERIALS AND METHODS: We searched MEDLINE/Pubmed, Web of Science, EMBASE via OVID, and grey literature for publications until December 28th, 2019, that compared the immunohistochemical expression of the two OKC subtypes. The studies were qualitatively assessed using the Critical Appraisal Tool for Case Series (Joana Briggs Institute). Sensitivity and specificity, positive and negative likelihood ratio, diagnostic odds ratio and area under the curve, and pooled estimates were calculated, using a random-effects model. RESULTS: Seventy-one studies were qualitatively analyzed; 61 markers were evaluated in one study and 32 in ≥ 2 studies. Twenty-five studies reported differential expression of 29 markers in the form of higher number of positive cells or greater staining intensity usually in BCNS-associated OKCs. Meta-analysis for bcl-2, Cyclin D1, CD56, CK18, p53, and PCNA showed that none of those markers is distinguishable between BCNS-associated and sporadic OKCs, in a 95% confidence interval. The risk of bias was high in 34 studies, moderate in 22, and low in 15. CONCLUSIONS: The present systematic review and meta-analysis uncovered that, although several immunohistochemical markers might characterize the OKC phenotype, they cannot discriminate between the BCNS-associated and sporadic OKCs. CLINICAL RELEVANCE: This study highlighted the requirement for additional screening for markers by immunohistochemistry, preferentially coupled to alternative diagnostic applications such as genomics technologies.

Transcription Factor Ets-2 Acts as a Preinduction Repressor of Interleukin-2 (IL-2) Transcription in Naive T Helper Lymphocytes.

IL-2 is the first cytokine produced when naive T helper (Th) cells are activated and differentiate into dividing pre-Th0 proliferating precursors. IL-2 expression is blocked in naive, but not activated or memory, Th cells by the transcription factor Ets-2 that binds to the antigen receptor response element (ARRE)-2 of the proximal IL-2 promoter. Ets-2 acts as an independent preinduction repressor in naive Th cells and does not interact physically with the transcription factor NFAT (nuclear factor of activated T-cells) that binds to the ARRE-2 in activated Th cells. In naive Th cells, Ets-2 mRNA expression, Ets-2 protein levels, and Ets-2 binding to ARRE-2 decrease upon cell activation followed by the concomitant expression of IL-2. Cyclosporine A stabilizes Ets-2 mRNA and protein when the cells are activated. Ets-2 silences directly constitutive or induced IL-2 expression through the ARRE-2. Conversely, Ets-2 silencing allows for constitutive IL-2 expression in unstimulated cells. Ets-2 binding to ARRE-2 in chromatin is stronger in naive compared with activated or memory Th cells; in the latter, Ets-2 participates in a change of the IL-2 promoter architecture, possibly to facilitate a quick response when the cells re-encounter antigen. We propose that Ets-2 expression and protein binding to the ARRE-2 of the IL-2 promoter are part of a strictly regulated process that results in a physiological transition of naive Th cells to Th0 cells upon antigenic stimulation. Malfunction of such a repression mechanism at the molecular level could lead to a disturbance of later events in Th cell plasticity, leading to autoimmune diseases or other pathological conditions.

Neutrophil-fibroblast crosstalk drives immunofibrosis in Crohn's disease through IFNα pathway.

Introduction: Crohn's disease (CD) is characterized by chronic inflammation and intestinal fibrosis leading to lifelong complications. However, the disease pathogenesis remains elusive, and the therapeutic options are limited. Here, we investigated the interaction between neutrophils and intestinal fibroblasts in the development of CD immunofibrosis, a disease mechanism predisposing to inflammatory and fibrotic complications. Methods: Peripheral neutrophils, enriched neutrophil extracellular traps (eNETs), serum, primary intestinal fibroblasts (PIFs) and intestinal biopsies from CD, ulcerative colitis (UC) patients, and healthy individuals (HI), were studied. Transcriptome analysis of neutrophils, multi-cytokine profiling and cell-based functional assays at mRNA/protein level were performed. Results: Compared to UC, PIFs from CD patients, independently to the presence of strictures, displayed a distinct pro-fibrotic phenotype characterized by negative Krüppellike Factor-2 (KLF2) and increased cellular communication network factor-2 (CCN2) expression leading to collagen production. In both UC and CD, PIFs-derived IL-8 acted as a culprit chemoattractant for neutrophils in the intestine, where CD neutrophils were accumulated close to fibrotic lesions. Functionally, only CD neutrophils via eNETs induced a CD-like phenotype in HI PIFs, suggesting their fibrotic plasticity. High IFNa in serum and IFΝ-responsive signature in peripheral neutrophils were observed in CD, distinguishing it from UC. Moreover, CD serum stimulated the release of fibrogenic eNETs from neutrophils in an IFNa-dependent manner, suggesting the priming role of IFNa in circulating neutrophils. Inhibition of eNETs or JAK signaling in neutrophils or PIFs prevented the neutrophil-mediated fibrotic effect on PIFs. Furthermore, both serum IFNa levels and mRNA levels of key IFN signaling components in neutrophils were wellcorrelated with CD severity. Conclusions: This study reveals the important role of the IFNa/neutrophil/fibroblast axis in CD immunofibrosis, suggesting candidate biomarkers and putative therapeutic targets.

MYCN Amplifications and Metabolic Rewiring in Neuroblastoma.

Cancer is a disease caused by (epi)genomic and gene expression abnormalities and characterized by metabolic phenotypes that are substantially different from the normal phenotypes of the tissues of origin. Metabolic reprogramming is one of the key features of tumors, including those established in the human nervous system. In this work, we emphasize a well-known cancerous genomic alteration: the amplification of MYCN and its downstream effects in neuroblastoma phenotype evolution. Herein, we extend our previous computational biology investigations by conducting an integrative workflow applied to published genomics datasets and comprehensively assess the impact of MYCN amplification in the upregulation of metabolism-related transcription factor (TF)-encoding genes in neuroblastoma cells. The results obtained first emphasized overexpressed TFs, and subsequently those committed in metabolic cellular processes, as validated by gene ontology analyses (GOs) and literature curation. Several genes encoding for those TFs were investigated at the mechanistic and regulatory levels by conducting further omics-based computational biology assessments applied on published ChIP-seq datasets retrieved from MYCN-amplified- and MYCN-enforced-overexpression within in vivo systems of study. Hence, we approached the mechanistic interrelationship between amplified MYCN and overexpression of metabolism-related TFs in neuroblastoma and showed that many are direct targets of MYCN in an amplification-inducible fashion. These results illuminate how MYCN executes its regulatory underpinnings on metabolic processes in neuroblastoma.

Antiphospholipid antibodies induce proinflammatory and procoagulant pathways in endothelial cells.

Antiphospholipid syndrome (APS) is an autoimmune thrombophilia characterized by recurrent thrombotic events and/or pregnancy morbidity in the presence of antiphospholipid antibodies detected either as anti-cardiolipin, anti-ß2 Glycoprotein I (anti-ß2GPI) or Lupus anticoagulant (LA). Endothelial deregulation characterizes the syndrome. To address gene expression changes accompanying the development of autoimmune phenotype in endothelial cells in the context of APS, we performed transcriptomics analysis in Human Umbilical Vein Endothelial Cells (HUVECs) stimulated with IgG from APS patients and ß2GPI, followed by intersection of RNA-seq data with published microarray and ChIP-seq results (Chromatin Immunoprecipitation). Our strategy revealed that during HUVEC activation diverse signaling pathways such as TNF-α, TGF-ß, MAPK38, and Hippo are triggered as indicated by Gene Ontology (GO) classification and pathway analysis. Finally, cell biology approaches performed side-by-side in naïve and stimulated cultured HUVECs, as well as, in placenta specimens derived from Healthy donors (HDs) and APS-patients verified the evolution of an APS-characteristic gene expression program in endothelial cells during the initial stages of the disease's development.

HGCA2.0: An RNA-Seq Based Webtool for Gene Coexpression Analysis in Homo sapiens.

Genes with similar expression patterns in a set of diverse samples may be considered coexpressed. Human Gene Coexpression Analysis 2.0 (HGCA2.0) is a webtool which studies the global coexpression landscape of human genes. The website is based on the hierarchical clustering of 55,431 Homo sapiens genes based on a large-scale coexpression analysis of 3500 GTEx bulk RNA-Seq samples of healthy individuals, which were selected as the best representative samples of each tissue type. HGCA2.0 presents subclades of coexpressed genes to a gene of interest, and performs various built-in gene term enrichment analyses on the coexpressed genes, including gene ontologies, biological pathways, protein families, and diseases, while also being unique in revealing enriched transcription factors driving coexpression. HGCA2.0 has been successful in identifying not only genes with ubiquitous expression patterns, but also tissue-specific genes. Benchmarking showed that HGCA2.0 belongs to the top performing coexpression webtools, as shown by STRING analysis. HGCA2.0 creates working hypotheses for the discovery of gene partners or common biological processes that can be experimentally validated. It offers a simple and intuitive website design and user interface, as well as an API endpoint.

Typical Enhancers, Super-Enhancers, and Cancers.

Non-coding segments of the human genome are enriched in cis-regulatory modules that constitute functional elements, such as transcriptional enhancers and Super-enhancers. A hallmark of cancer pathogenesis is the dramatic dysregulation of the "archetype" gene expression profiles of normal human cells. Genomic variations can promote such deficiencies when occurring across enhancers and Super-enhancers, since they affect their mechanistic principles, their functional capacity and specificity, and the epigenomic features of the chromatin microenvironment across which these regulatory elements reside. Here, we comprehensively describe: fundamental mechanisms of gene expression dysregulation in cancers that involve genomic abnormalities within enhancers' and Super-enhancers' (SEs) sequences, which alter the expression of oncogenic transcription factors (TFs); cutting-edge technologies applied for the analysis of variation-enriched hotspots of the cancer genome; and pharmacological approaches for the treatment of Super-enhancers' aberrant function. Finally, we provide an intratumor meta-analysis, which highlights that genomic variations in transcription-factor-driven tumors are accompanied overexpression of genes, a portion of which encodes for additional cancer-related transcription factors.

The TΑp63/BCL2 axis represents a novel mechanism of clinical aggressiveness in chronic lymphocytic leukemia.

The TA-isoform of the p63 transcription factor (TAp63) has been reported to contribute to clinical aggressiveness in chronic lymphocytic leukemia (CLL) in a hitherto elusive way. Here, we sought to further understand and define the role of TAp63 in the pathophysiology of CLL. First, we found that elevated TAp63 expression levels are linked with adverse clinical outcomes, including disease relapse and shorter time-to-first treatment and overall survival. Next, prompted by the fact that TAp63 participates in an NF-κB/TAp63/BCL2 antiapoptotic axis in activated mature, normal B cells, we explored molecular links between TAp63 and BCL2 also in CLL. We documented a strong correlation at both the protein and the messenger RNA (mRNA) levels, alluding to the potential prosurvival role of TAp63. This claim was supported by inducible downregulation of TAp63 expression in the MEC1 CLL cell line using clustered regularly interspaced short palindromic repeats (CRISPR) system, which resulted in downregulation of BCL2 expression. Next, using chromatin immunoprecipitation (ChIP) sequencing, we examined whether BCL2 might constitute a transcriptional target of TAp63 and identified a significant binding profile of TAp63 in the BCL2 gene locus, across a genomic region previously characterized as a super enhancer in CLL. Moreover, we identified high-confidence TAp63 binding regions in genes mainly implicated in immune response and DNA-damage procedures. Finally, we found that upregulated TAp63 expression levels render CLL cells less responsive to apoptosis induction with the BCL2 inhibitor venetoclax. On these grounds, TAp63 appears to act as a positive modulator of BCL2, hence contributing to the antiapoptotic phenotype that underlies clinical aggressiveness and treatment resistance in CLL.

The Causes and Consequences of Spatial Organization of the Genome in Regulation of Gene Expression.

Regulation of gene expression in time, space and quantity is orchestrated by the functional interplay of cis-acting elements and trans-acting factors. Our current view postulates that transcription factors recognize enhancer DNA and read the transcriptional regulatory code by cooperative DNA binding to specific DNA motifs, thus instructing the recruitment of transcriptional regulatory complexes forming a plethora of higher-ordered multi-protein-DNA and protein-protein complexes. Here, we reviewed the formation of multi-dimensional chromatin assemblies implicated in gene expression with emphasis on the regulatory role of enhancer hubs as coordinators of stochastic gene expression. Enhancer hubs contain many interacting regulatory elements and represent a remarkably dynamic and heterogeneous network of multivalent interactions. A functional consequence of such complex interaction networks could be that individual enhancers function synergistically to ensure coordination, tight control and robustness in regulation of expression of spatially connected genes. In this review, we discuss fundamental paradigms of such inter- and intra- chromosomal associations both in the context of immune-related genes and beyond.

The Role of Coronavirus RNA-Processing Enzymes in Innate Immune Evasion.

Viral RNA sensing triggers innate antiviral responses in humans by stimulating signaling pathways that include crucial antiviral genes such as interferon. RNA viruses have evolved strategies to inhibit or escape these mechanisms. Coronaviruses use multiple enzymes to synthesize, modify, and process their genomic RNA and sub-genomic RNAs. These include Nsp15 and Nsp16, whose respective roles in RNA capping and dsRNA degradation play a crucial role in coronavirus escape from immune surveillance. Evolutionary studies on coronaviruses demonstrate that genome expansion in Nidoviruses was promoted by the emergence of Nsp14-ExoN activity and led to the acquisition of Nsp15- and Nsp16-RNA-processing activities. In this review, we discuss the main RNA-sensing mechanisms in humans as well as recent structural, functional, and evolutionary insights into coronavirus Nsp15 and Nsp16 with a view to potential antiviral strategies.

Composite macroH2A/NRF-1 Nucleosomes Suppress Noise and Generate Robustness in Gene Expression.

The histone variant macroH2A (mH2A) has been implicated in transcriptional repression, but the molecular mechanisms that contribute to global mH2A-dependent genome regulation remain elusive. Using chromatin immunoprecipitation sequencing (ChIP-seq) coupled with transcriptional profiling in mH2A knockdown cells, we demonstrate that singular mH2A nucleosomes occupy transcription start sites of subsets of both expressed and repressed genes, with opposing regulatory consequences. Specifically, mH2A nucleosomes mask repressor binding sites in expressed genes but activator binding sites in repressed genes, thus generating distinct chromatin landscapes that limit genetic or extracellular inductive signals. We show that composite nucleosomes containing mH2A and NRF-1 are stably positioned on gene regulatory regions and can buffer transcriptional noise associated with antiviral responses. In contrast, mH2A nucleosomes without NRF-1 bind promoters weakly and mark genes with noisier gene expression patterns. Thus, the strategic position and stabilization of mH2A nucleosomes in human promoters defines robust gene expression patterns.

cgChIP: a cell type- and gene-specific method for chromatin analysis.

Hox and other homeobox-containing genes encode critical transcriptional regulators of animal development. Although these genes are well known for their roles in the body axis and appendage development, little is known regarding the mechanisms by which these factors influence chromatin landscapes. Chromatin structure can have a profound influence on gene expression during animal body formation. However, when applied to developing embryos, conventional chromatin analysis of genes and cis-regulatory modules (CRMs) typically lacks the required cell type-specific resolution due to the heterogeneous nature of animal bodies. Here we present a strategy to analyze both the composition and conformation of in vivo-tagged CRM sequences in a cell type-specific manner, using as a system Drosophila embryos. We term this method cgChIP (cell- and gene-specific Chromatin Immunoprecipitation) by which we access and analyze regulatory chromatin in specific cell types. cgChIP is an in vivo method designed to analyze genetic elements derived from limited cell populations. cgChIP can be used for both the analysis of chromatin structure (e.g., long-distance interactions between DNA elements) and the composition of histones and histone modifications and the occupancy of transcription factors and chromatin modifiers. This method was applied to the Hox target gene Distalless (Dll), which encodes for a homeodomain-containing transcription factor critical for the formation of appendages in Drosophila. However, cgChIP can be applied in diverse animal models to better dissect CRM-dependent gene regulation and body formation in developing animals.

Developmental regulation of chromatin conformation by Hox proteins in Drosophila.

We present a strategy to examine the chromatin conformation of individual loci in specific cell types during Drosophila embryogenesis. Regulatory DNA is tagged with binding sites (lacO) for LacI, which is used to immunoprecipitate the tagged chromatin from specific cell types. We applied this approach to Distalless (Dll), a gene required for limb development in Drosophila. We show that the local chromatin conformation at Dll depends on the cell type: in cells that express Dll, the 5' regulatory region is in close proximity to the Dll promoter. In Dll-nonexpressing cells this DNA is in a more extended configuration. In addition, transcriptional activators and repressors are bound to Dll regulatory DNA in a cell type-specific manner. The pattern of binding by GAGA factor and the variant histone H2Av suggest that they play a role in the regulation of Dll chromatin conformation in expressing and nonexpressing cell types, respectively.

Epigenetic determination of a cell-specific gene expression program by ATF-2 and the histone variant macroH2A.

Transcriptional activation of the interleukin-8 (IL-8) gene is restricted to distinct cell types, although the transcriptional regulatory proteins controlling IL-8 gene expression are ubiquitous. We show that cell-specific transcription of IL-8 is due to the distinct chromatin architecture on the enhancer/promoter before the arrival of the inducing signal. In expressing epithelial cells the enhancer/promoter is nucleosome-free, whereas in non-expressing B cells a nucleosome masks the entire regulatory region. The B-cell-specific nucleosome contains the histone variant macroH2A, which is responsible for preventing transcription factor binding. Recruitment of the repressive macroH2A nucleosome requires direct interactions between ATF-2 bound to the nearby AP1 site and macroH2A and it is regulated by DNA-induced protein allostery. siRNA against ATF-2 or macroH2A rescues IL-8 transcription in B cells. Thus, a transcription factor can work as a transcriptional repressor by orchestrating and maintaining the assembly of specialized local chromatin architectures.