NF-Y Binding Site Architecture Defines a C-Fos Targeted Promoter Class

On the NF-Y regulome as in ENCODE (2019)

NF-Y is a trimeric Transcription Factor -TF- which binds with high selectivity to the conserved CCAAT element. Individual ChIP-seq analysis as well as ENCODE have progressively identified locations shared by other TFs. Here, we have analyzed data introduced by ENCODE over the last five years in K562, HeLa-S3 and GM12878, including several chromatin features, as well RNA-seq profiling of HeLa cells after NF-Y inactivation. We double the number of sequence-specific TFs and co-factors reported. We catalogue them in 4 classes based on co-association criteria, infer target genes categorizations, identify positional bias of binding sites and gene expression changes. Larger and novel co-associations emerge, specifically concerning subunits of repressive complexes as well as RNA-binding proteins. On the one hand, these data better define NF-Y association with single members of major classes of TFs, on the other, they suggest that it might have a wider role in the control of mRNA production.

The ongoing ENCODE consortium represents a useful compendium of locations of TFs, chromatin marks, gene expression data. In previous reports, we identified modules of CCAAT-binding NF-Y with individual TFs. Here, we analyzed all 363 factors currently present: 68 with enrichment of CCAAT in their locations, 38 with overlap of peaks. New sequence-specific TFs, co-activators and co-repressors are reported. Co-association patterns correspond to specific targeted genes categorizations and gene expression changes, as assessed by RNA-seq after NF-Y inactivation. These data widen and better define a coherent model of synergy of NF-Y with selected groups of TFs and co-factors.

MicroRNA‑21 serves an important role during PAOO‑facilitated orthodontic tooth movement

Periodontal accelerate osteogenesis orthodontics (PAOO) is an extension of described techniques that surgically alter the alveolar bone; however, the specific mechanism underlying the technique is not completely understood. The aim of the present study was to evaluate the roles of microRNA (miR)-21 during PAOO. Sprague-Dawley rats were divided into the following four groups: i) Group tooth movement (TM), underwent TM and were administered normal saline (NS); ii) Group PAOO, underwent PAOO + TM and were administered NS; iii) Group agomiR-21, underwent PAOO + TM and were administered agomiR-21; and iv) Group antagomiR-21, underwent PAOO + TM and were administered antagomiR-21. To validate the rat model of PAOO, morphological analyses were performed and measurements were collected. Reverse transcription-quantitative PCR, western blotting and immunohistochemical staining were performed to examine the expression levels of programmed cell death 4 (PDCD4), activin A receptor type 2B (ACVR2b), receptor activator of NF-κΒ ligand (RANKL) and C-Fos. Dual-luciferase reporter assays were performed to validate PDCD4 as a target of miR-21 in vitro. Following 7 days of treatment, the TM distance of group PAOO was longer compared with groups TM and antagomiR-21 (P<0.05), but shorter compared with group agomiR-21 (P<0.05). Tartrate-resistant acid phosphatase staining indicated that following treatment with agomiR-21, osteoclast activity was notably increased, whereas the mRNA and protein expression levels of PDCD4 were notably decreased compared with group PAOO. The mRNA and protein expression levels of RANKL and C-Fos in group agomiR-21 were notably increased compared with group PAOO, whereas group antagomiR-21 displayed the opposite pattern (P<0.05). With regard to ACVR2b, no significant differences were observed among the group agomiR-21 and antagomiR-21 compared with group PAOO. Bioinformatics analysis predicted that PDCD4 was a potential target gene of miR-21, and dual-luciferase reporter assays demonstrated that miR-21 directly targeted PDCD4. In conclusion, the present study demonstrated that miR-21 serves an important role during PAOO-mediated orthodontic TM.

Integrating Peak Colocalization and Motif Enrichment Analysis for the Discovery of Genome-Wide Regulatory Modules and Transcription Factor Recruitment Rules

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-Seq) has opened new avenues of research in the genome-wide characterization of regulatory DNA-protein interactions at the genetic and epigenetic level. As a consequence, it has become the de facto standard for studies on the regulation of transcription, and literally thousands of data sets for transcription factors and cofactors in different conditions and species are now available to the scientific community. However, while pipelines and best practices have been established for the analysis of a single experiment, there is still no consensus on the best way to perform an integrated analysis of multiple datasets in the same condition, in order to identify the most relevant and widespread regulatory modules composed by different transcription factors and cofactors. We present here a computational pipeline for this task, that integrates peak summit colocalization, a novel statistical framework for the evaluation of its significance, and motif enrichment analysis. We show examples of its application to ENCODE data, that led to the identification of relevant regulatory modules composed of different factors, as well as the organization on DNA of the binding motifs responsible for their recruitment.

NF-YA overexpression protects from glutamine deprivation

The heterotrimeric transcription factor NF-Y binds to CCAAT boxes of genes of glutamine metabolism. We set out to study the role of the regulatory NF-YA subunit in this pathway. We produced U2OS and A549 clones stably overexpressing -OE- the two splicing isoforms of NF-YA. NF-YA OE cells show normal growth and colony formation rates, but they become resistant to cell death upon glutamine deprivation. Increased mRNA and protein expression of the key biosynthetic enzyme GLUL in U2OS entails increased production of endogenous glutamine upon deprivation. The use of GLUL inhibitors dampens the NF-YA-mediated effect. NF-YA OE prevents activation of the pro-apoptotic transcription factor CHOP/DDIT3. Elevated basal levels of SERCA1/2, coding for the molecular target of Thapsigargin, correlate with resistance of NF-YA OE cells to the drug. The work represents a proof-of-principle that elevated levels of NF-YA, as found in some tumor types, helps altering cancer metabolic pathways.

The Plant NF-Y DNA Matrix In Vitro and In Vivo

Nuclear Factor Y (NF-Y) is an evolutionarily conserved trimer formed by a Histone-Fold Domain (HFD) heterodimeric module shared by core histones, and the sequence-specific NF-YA subunit. In plants, the genes encoding each of the three subunits have expanded in number, giving rise to hundreds of potential trimers. While in mammals NF-Y binds a well-characterized motif, with a defined matrix centered on the CCAAT box, the specificity of the plant trimers has yet to be determined. Here we report that Arabidopsis thaliana NF-Y trimeric complexes, containing two different NF-YA subunits, bind DNA in vitro with similar affinities. We assayed precisely sequence-specificity by saturation mutagenesis, and analyzed genomic DNA sites bound in vivo by selected HFDs. The plant NF-Y CCAAT matrix is different in nucleotides flanking CCAAT with respect to the mammalian matrix, in vitro and in vivo. Our data point to flexible DNA-binding rules by plant NF-Ys, serving the scope of adapting to a diverse audience of genomic motifs.

Disentangling transcription factor binding site complexity

The binding motifs of many transcription factors (TFs) comprise a higher degree of complexity than a single position weight matrix model permits. Additional complexity is typically taken into account either as intra-motif dependencies via more sophisticated probabilistic models or as heterogeneities via multiple weight matrices. However, both orthogonal approaches have limitations when learning from in vivo data where binding sites of other factors in close proximity can interfere with motif discovery for the protein of interest. In this work, we demonstrate how intra-motif complexity can, purely by analyzing the statistical properties of a given set of TF-binding sites, be distinguished from complexity arising from an intermix with motifs of co-binding TFs or other artifacts. In addition, we study the related question whether intra-motif complexity is represented more effectively by dependencies, heterogeneities or variants in between. Benchmarks demonstrate the effectiveness of both methods for their respective tasks and applications on motif discovery output from recent tools detect and correct many undesirable artifacts. These results further suggest that the prevalence of intra-motif dependencies may have been overestimated in previous studies on in vivo data and should thus be reassessed.

Epigenetic alterations of CYLD promoter modulate its expression in gastric adenocarcinoma: A footprint of infections

Gastric cancer (GC) is one of the most common causes of cancer‐related death in the world, with multiple genetic and epigenetic alterations involved in disease development. CYLD tumor suppressor gene encodes a multifunctional deubiquitinase which negatively regulates various signaling pathways. Deregulation of this gene has been found in different types of cancer. This study aimed to evaluate for the first time the CpG island methylation pattern of CYLD gene promoter, and its expression level in gastric adenocarcinoma. CYLD messenger RNA expression and promoter methylation in 53 tumoral and their non‐neoplastic counterpart tissues were assessed using quantitative polymerase chain reaction and bisulfite sequencing. Also, we investigated the impacts of the infectious agents including Helicobacter pylori (H. pylori), EBV, and CMV on CYLD expression and promoter methylation in GC. Results showed that the expression level of CYLD was downregulated in GC, and was significantly associated with gender (female), patient’s age (<60), high grade, and no lymph‐node metastasis (p = 0.001, 0.002, 0.03, and 0.003, respectively). Among the 31 analyzed CpG sites located in about 600 bp region within the promoter, two CpG sites were hypermethylated in GC tissues. We also found a significant inverse association between DNA promoter methylation and CYLD expression (p = 0.02). Furthermore, a direct association between H. pylori, EBV, and CMV infections with hypermethylation and reduced CYLD expression was observed (p = 0.04, 0.03, and 0.03, respectively). Our findings indicate that CYLD is downregulated in GC. Infectious agents may influence CYLD expression.

Genome wide features, distribution and correlations of NF-Y binding sites

NF-Y is a trimeric transcription factor that binds on DNA the CCAAT-box motif. In this article we reviewed and complemented with additional bioinformatic analysis existing data on genome-wide NF-Y binding characterization in human, reaching the following main conclusions: (1) about half of NF-Y binding sites are located at promoters, about 60-80 base pairs from transcription start sites; NF-Y binding to distal genomic regions takes place at inactive chromatin loci and/or DNA repetitive elements more often than active enhancers; (2) on almost half of its binding sites, regardless of their genomic localization (promoters or distal regions), NF-Y finds on DNA more than one CCAAT-box, and most of those multiple CCAAT binding loci present precise spacing and organization of the elements composing them; (3) there exists a well defined class of transcription factors that show genome-wide co-localization with NF-Y. Some of them lack their canonical binding site in binding regions overlapping with NF-Y, hence hinting at NF-Y mediated recruitment, while others show a precise positioning on DNA of their binding sites with respect to the CCAAT box bound by NF-Y. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.