NF-Y coassociates with FOS at promoters, enhancers, repetitive elements, and inactive chromatin regions, and is stereo-positioned with growth-controlling transcription factors

PhpCNF-Y transcription factor infiltrates heterochromatin to generate cryptic intron-containing transcripts crucial for small RNA production

The assembly of repressive heterochromatin in eukaryotic genomes is crucial for silencing lineage-inappropriate genes and repetitive DNA elements. Paradoxically, transcription of repetitive elements within constitutive heterochromatin domains is required for RNA-based mechanisms, such as the RNAi pathway, to target heterochromatin assembly proteins. However, the mechanism by which heterochromatic repeats are transcribed has been unclear. Using fission yeast, we show that the conserved trimeric transcription factor (TF) PhpCNF-Y complex can infiltrate constitutive heterochromatin via its histone-fold domains to transcribe repeat elements. PhpCNF-Y collaborates with a Zn-finger containing TF to bind repeat promoter regions with CCAAT boxes. Mutating either the TFs or the CCAAT binding site disrupts the transcription of heterochromatic repeats. Although repeat elements are transcribed from both strands, PhpCNF-Y-dependent transcripts originate from only one strand. These TF-driven transcripts contain multiple cryptic introns which are required for the generation of small interfering RNAs (siRNAs) via a mechanism involving the spliceosome and RNAi machinery. Our analyses show that siRNA production by this TF-mediated transcription pathway is critical for heterochromatin nucleation at target repeat loci. This study reveals a mechanism by which heterochromatic repeats are transcribed, initiating their own silencing by triggering a primary cascade that produces siRNAs necessary for heterochromatin nucleation.

The role(s) of NF-Y in development and differentiation

NF-Y is a conserved sequence-specific trimeric Transcription Factor -TF- binding to the CCAAT element. We review here the role(s) in development, from pre-implantation embryo to terminally differentiated tissues, by rationalizing and commenting on genetic, genomic, epigenetic and biochemical studies. This effort brings to light the impact of NF-YA isoforms on stemness and differentiation, as well as binding to distal vs promoter proximal sites and connections with selected TFs.

Novel function of U7 snRNA in the repression of HERV1/LTR12s and lincRNAs in human cells

U7 snRNA is part of the U7 snRNP complex, required for the 3' end processing of replication-dependent histone pre-mRNAs in S phase of the cell cycle. Here, we show that U7 snRNA plays another function in inhibiting the expression of a subset of long terminal repeats of human endogenous retroviruses (HERV1/LTR12s) and LTR12-containing long intergenic noncoding RNAs (lincRNAs), both bearing sequence motifs that perfectly match the 5' end of U7 snRNA. We demonstrate that U7 snRNA inhibits LTR12 and lincRNA transcription and propose a mechanism in which U7 snRNA hampers the binding/activity of the NF-Y transcription factor to CCAAT motifs within LTR12 elements. Thereby, U7 snRNA plays a protective role in maintaining the silencing of deleterious genetic elements in selected types of cells.

Genomic binding of NF-Y in mouse and human cells

NF-Y is a Transcription Factor that regulates transcription through binding to the CCAAT-box. To understand its strategy, we analyzed 16 ChIP-seq datasets from human and mouse cells. Shared loci, mostly located in promoters of expressed genes of cell cycle, metabolism and gene expression pathways, are associated with histone marks of active chromatin and specific modules of TFs. Other peaks are in enhancers and Transposable Elements -TE- of retroviral origin in human and mouse. We evaluated the relationship with USF1, a common synergistic partner in promoters and MLT1 TEs, upon NF-YB inactivation: USF1 binding decreases in promoters, modestly in MLT1, suggesting a pioneering role of NF-Y in formers, not in the latters. These data define a common set of NF-Y functional targets across different mammalian cell types, suggesting a pioneering role in promoters with respect to TEs.

BATF2 is a regulator of interferon-γ signaling in astrocytes during neuroinflammation

Astrocytic interferon (IFN)γ signaling is associated with a reduction in neuroinflammation. We have previously shown that the benefits of astrocytic IFNγ arise from a variety of mechanisms; however, downstream effectors responsible for regulating this protection are unknown. We address this by identifying a specific transcription factor that may play a key role in modulating the consequences of IFNγ signaling. RNA-sequencing of primary human astrocytes treated with IFNγ revealed basic leucine zipper ATF-like transcription factor ( BATF )2 as a highly expressed interferon-specific gene. Primarily studied in the periphery, BATF2 has been shown to exert both inflammatory and protective functions; however, its function in the central nervous system (CNS) is unknown. Here, we demonstrate that human spinal cord astrocytes upregulate BATF2 transcript and protein in an IFNγ-specific manner. Additionally, we found that BATF2 prevents overexpression of interferon regulatory factor (IRF)1 and IRF1 targets such as Caspase-1, which are known downstream pro-inflammatory mediators. We also show that Batf2 -/- mice exhibit exacerbated clinical disease severity in a murine model of CNS autoimmunity, characterized by an increase in both CNS immune cell infiltration and demyelination. Batf2 -/- mice also exhibit increased astrocyte-specific expression of IRF1 and Caspase-1, suggesting an amplified interferon response in vivo . Further, we demonstrate that BATF2 is expressed primarily in astrocytes in MS lesions and that this expression is co-localized with IRF1. Collectively, our results further support a protective role for IFNγ and implicate BATF2 as a key suppressor of overactive immune signaling in astrocytes during neuroinflammation.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

Transcription factor NFYa controls cardiomyocyte metabolism and proliferation during mouse fetal heart development

Cardiomyocytes are highly metabolic cells responsible for generating the contractile force in the heart. During fetal development and regeneration, these cells actively divide but lose their proliferative activity in adulthood. The mechanisms that coordinate their metabolism and proliferation are not fully understood. Here, we study the role of the transcription factor NFYa in developing mouse hearts. Loss of NFYa alters cardiomyocyte composition, causing a decrease in immature regenerative cells and an increase in trabecular and mature cardiomyocytes, as identified by spatial and single-cell transcriptome analyses. NFYa-deleted cardiomyocytes exhibited reduced proliferation and impaired mitochondrial metabolism, leading to cardiac growth defects and embryonic death. NFYa, interacting with cofactor SP2, activates genes linking metabolism and proliferation at the transcription level. Our study identifies a nodal role of NFYa in regulating prenatal cardiac growth and a previously unrecognized transcriptional control mechanism of heart metabolism, highlighting the importance of mitochondrial metabolism during heart development and regeneration.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

RNA exosome ribonuclease DIS3 degrades Pou6f1 to promote mouse pre-implantation cell differentiation

Mammalian development is precisely controlled by cell differentiation. Identifying new regulators and investigating their interactions provide insight into genetic networks defining pre-implantation development. We established a knockout mouse model of Dis3, an exosome associated ribonuclease. Homozygous Dis3 null embryos arrest at the morula stage of development. Using single-embryo RNA sequencing (RNA-seq), we observed persistence of Pou6f1 mRNA in homozygous null Dis3 embryos and that the cognate protein represses transcription of Nanog and Cdx2. The resultant defects in cell differentiation disrupt the morula-to-blastocyst transition and are embryonic lethal. Microinjection of Dis3 mRNA into zygotes rescues the phenotype. Point mutations of Dis3 ribonuclease in individual blastomeres prevents their incorporation into embryos. To overcome the paucity of embryos, we derived homozygous Dis3 null mouse embryonic stem cells to identify additional gene targets of POU6F1. Our findings delineate a regulatory pathway of DIS3-POU6F1 in pre-implantation mammalian embryogenesis.

LEAFY COTYLEDONs: Connecting different stages of plant development

The life of higher plants progresses successively through embryonic, juvenile, adult, and reproductive stages. LEAFY COTYLEDON (LEC) transcription factors, first discovered in Arabidopsis thaliana several decades ago, play a key role in regulating plant embryonic development, seed maturation, and subsequent growth. Existing studies have demonstrated that LECs together with other transcription factors form a huge and complex regulatory network to regulate many aspects of plant growth and development and respond to environmental stresses. Here, we focus on the role that has received little attention about the LECs linking different developmental stages and generational cycles in plants. We summarize the current fragmented research progress on the LECs role and molecular mechanism in connecting embryonic and vegetative growth periods and the reproductive stage. Furthermore, the possibility of LECs controlling the maintenance and transition of plant growth stages through epigenetic modifications is discussed.

Drosophila transcription factor NF-Y suppresses transcription of the lipase 4 gene, a key gene for lipid storage

The CCAAT motif-binding factor NF-Y consists of three different subunits, NF-YA, NF-YB, and NF-YC. Although it is suggested that NF-Y activity is essential for normal tissue homeostasis, survival, and metabolic function, its precise role in lipid metabolism is not clarified yet. In Drosophila, eye disc specific knockdown of Drosophila NF-YA (dNF-YA) induced aberrant morphology of the compound eye, the rough eye phenotype in adults and mutation of the lipase 4 (lip4) gene suppressed the rough eye phenotype. RNA-seq analyses with dNF-YA knockdown third instar larvae identified the lip4 gene as one of the genes that are up-regulated by the dNF-YA knockdown. We identified three dNF-Y-binding consensuses in the 5'flanking region of the lip4 gene, and a chromatin immunoprecipitation assay with the specific anti-dNF-YA IgG demonstrated dNF-Y binding to this genomic region. The luciferase transient expression assay with cultured Drosophila S2 cells and the lip4 promoter-luciferase fusion genes with and without mutations in the dNF-Y-binding consensuses showed that each of the three dNF-Y consensus sequences negatively regulated lip4 gene promoter activity. Consistent with these results, qRT-PCR analysis with the dNF-YA knockdown third instar larvae revealed that endogenous lip4 mRNA levels were increased by the knockdown of dNF-YA in vivo. The specific knockdown of dNF-YA in the fat body with the collagen-GAL4 driver resulted in smaller oil droplets in the fat body cells. Collectively, these results suggest that dNF-Y is involved in lipid storage through its negative regulation of lip4 gene transcription.

The transcription factor NFYC positively regulates expression of MHCIa in the red-spotted grouper (Epinephelus akaara)

Mammalian studies have shown that the nuclear transcription factor Y (NFYC) regulates the expression of major histocompatibility complex (MHC) by binding to CCAAT-box on promoters. However, few studies have focused on the regulatory mechanisms of NFYC in MHC pathway in fish. To explore the transcriptional regulatory mechanism of MHCIa in fish, we characterized NFYC and MHCIa of red-spotted grouper (Epinephelus akaara) (named EaNFYC and EaMHCIa, respectively). The EaNFYC genome sequence is 13,796 bp and contains 1,065 bp open reading frame. It is composed of ten exons and nine introns and encode a 354 amino acid sequence. The putative EaNFYC protein sequence shared 67.2-99.4% identity to vertebrate NFYC and possesses a typically conserved domain (histone- or haem-associated protein 5 domain (HAP5)) at the N-terminus. Transcripts of both EaNFYC and EaMHCIa were ubiquitously expressed in all detect tissues, and higher mRNA levels were detected in immune-relevant tissues (middle-kidney). EaNFYC expression increased after treatment with polyinosinic: polycytidylic acid, lipopolysaccharide, nervous necrosis virus, zymosan A, and Singapore grouper iridovirus. Analysis of subcellular localization indicated that EaNFYC was localized at the cell nucleus only. Furthermore, overexpression of EaNFYC significantly stimulated the expression of EaMHCIa, interferon signalling molecules and inflammatory cytokine. The region -878 bp to +82 bp of EaMHCIa promoter was identified to be the core promoter which EaNFYC take effect on. Additionally, point mutations and electrophoretic mobility shift assays verified that NFYC activate MHCIa expression by binding at the M1 and M2 binding sites that do not contain CCAAT-box. These results contribute to elucidating the function of fish NFYC on MHC transcriptional mechanisms, and provide the first evidence of positive regulation of MHCIa expression by NFYC in fish.

The USR domain of USF1 mediates NF-Y interactions and cooperative DNA binding

The trimeric CCAAT-binding NF-Y is a "pioneer" Transcription Factor -TF- known to cooperate with neighboring TFs to regulate gene expression. Genome-wide analyses detected a precise stereo-alignment -10/12 bp- of CCAAT with E-box elements and corresponding colocalization of NF-Y with basic-Helix-Loop-Helix (bHLH) TFs. We dissected here NF-Y interactions with USF1 and MAX. USF1, but not MAX, cooperates in DNA binding with NF-Y. NF-Y and USF1 synergize to activate target promoters. Reconstruction of complexes by structural means shows independent DNA binding of MAX, whereas USF1 has extended contacts with NF-Y, involving the USR, a USF-specific amino acid sequence stretch required for trans-activation. The USR is an intrinsically disordered domain and adopts different conformations based on E-box-CCAAT distances. Deletion of the USR abolishes cooperative DNA binding with NF-Y. Our data indicate that the functionality of certain unstructured domains involves adapting to small variation in stereo-alignments of the multimeric TFs sites.

PRMT6 physically associates with nuclear factor Y to regulate photoperiodic flowering in Arabidopsis

The timing of floral transition is critical for reproductive success in flowering plants. In long-day (LD) plant Arabidopsis, the floral regulator gene FLOWERING LOCUS T (FT) is a major component of the mobile florigen. FT expression is rhythmically activated by CONSTANS (CO), and specifically accumulated at dusk of LDs. However, the underlying mechanism of adequate regulation of FT transcription in response to day-length cues to warrant flowering time still remains to be investigated. Here, we identify a homolog of human protein arginine methyltransferases 6 (HsPRMT6) in Arabidopsis, and confirm AtPRMT6 physically interacts with three positive regulators of flowering Nuclear Factors YC3 (NF-YC3), NF-YC9, and NF-YB3. Further investigations find that AtPRMT6 and its encoding protein accumulate at dusk of LDs. PRMT6-mediated H3R2me2a modification enhances the promotion of NF-YCs on FT transcription in response to inductive LD signals. Moreover, AtPRMT6 and its homologues proteins AtPRMT4a and AtPRMT4b coordinately inhibit the expression of FLOWERING LOCUS C, a suppressor of FT. Taken together, our study reveals the role of arginine methylation in photoperiodic pathway and how the PRMT6-mediating H3R2me2a system interacts with NF-CO module to dynamically control FT expression and facilitate flowering time.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s42994-021-00065-y.

TCP7 interacts with Nuclear Factor-Ys to promote flowering by directly regulating SOC1 in Arabidopsis

The success of plant reproduction depends on the timely transition from the vegetative phase to reproductive growth, a process often referred to as flowering. Although several plant-specific transcription factors belonging to the Teosinte Branched 1/Cycloidea/Proliferating Cell Factor (TCP) family are reportedly involved in the regulation of flowering in Arabidopsis, the molecular mechanisms, especially for Class I TCP members, are poorly understood. Here, we genetically identified Class I TCP7 as a positive regulator of flowering time. Protein interaction analysis indicated that TCP7 interacted with several Nuclear Factor-Ys (NF-Ys), known as the 'pioneer' transcription factors; CONSTANS (CO), a main photoperiod regulator of flowering. SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was differentially expressed in the dominant-negative mutant of TCP7 (lcu) and the loss-of-function mutant of Class I TCP members (septuple). Additionally, we obtained genetic and molecular evidence showing that TCP7 directly activates the flowering integrator gene, SOC1. Moreover, TCP7 synergistically activated SOC1 expression upon interacting with CO and NF-Ys in vivo. Collectively, our results provide compelling evidence that TCP7 synergistically interacts with NF-Ys to activate the transcriptional expression of the flowering integrator gene SOC1.

Proteasome inhibition alters mitotic progression through the upregulation of centromeric α-Satellite RNAs

Cell cycle progression requires control of the abundance of several proteins and RNAs over space and time to properly transit from one phase to the next and to ensure faithful genomic inheritance in daughter cells. The proteasome, the main protein degradation system of the cell, facilitates the establishment of a proteome specific to each phase of the cell cycle. Its activity also strongly influences transcription. Here, we detected the upregulation of repetitive RNAs upon proteasome inhibition in human cancer cells using RNA‐seq. The effect of proteasome inhibition on centromeres was remarkable, especially on α‐Satellite RNAs. We showed that α‐Satellite RNAs fluctuate along the cell cycle and interact with members of the cohesin ring, suggesting that these transcripts may take part in the regulation of mitotic progression. Next, we forced exogenous overexpression and used gapmer oligonucleotide targeting to demonstrate that α‐Sat RNAs have regulatory roles in mitosis. Finally, we explored the transcriptional regulation of α‐Satellite DNA. Through in silico analyses, we detected the presence of CCAAT transcription factor‐binding motifs within α‐Satellite centromeric arrays. Using high‐resolution three‐dimensional immuno‐FISH and ChIP‐qPCR, we showed an association between the α‐Satellite upregulation and the recruitment of the transcription factor NFY‐A to the centromere upon MG132‐induced proteasome inhibition. Together, our results show that the proteasome controls α‐Satellite RNAs associated with the regulation of mitosis.

YAP and TAZ are transcriptional co-activators of AP-1 proteins and STAT3 during breast cellular transformation

The YAP and TAZ paralogs are transcriptional co-activators recruited to target sites by TEAD proteins. Here, we show that YAP and TAZ are also recruited by JUNB (a member of the AP-1 family) and STAT3, key transcription factors that mediate an epigenetic switch linking inflammation to cellular transformation. YAP and TAZ directly interact with JUNB and STAT3 via a WW domain important for transformation, and they stimulate transcriptional activation by AP-1 proteins. JUNB, STAT3, and TEAD co-localize at virtually all YAP/TAZ target sites, yet many target sites only contain individual AP-1, TEAD, or STAT3 motifs. This observation and differences in relative crosslinking efficiencies of JUNB, TEAD, and STAT3 at YAP/TAZ target sites suggest that YAP/TAZ is recruited by different forms of an AP-1/STAT3/TEAD complex depending on the recruiting motif. The different classes of YAP/TAZ target sites are associated with largely non-overlapping genes with distinct functions. A small minority of target sites are YAP- or TAZ-specific, and they are associated with different sequence motifs and gene classes from shared YAP/TAZ target sites. Genes containing either the AP-1 or TEAD class of YAP/TAZ sites are associated with poor survival of breast cancer patients with the triple-negative form of the disease.

NF-Y Subunits Overexpression in HNSCC

Simple Summary

Cancer cells have altered gene expression profiles. This is ultimately elicited by altered structure, expression or binding of transcription factors to regulatory regions of genomes. The CCAAT-binding trimer is a pioneer transcription factor involved in the activation of “cancer” genes. We and others have shown that the regulatory NF-YA subunit is overexpressed in epithelial cancers. Here, we examined large datasets of bulk gene expression profiles, as well as single-cell data, in head and neck squamous cell carcinomas by bioinformatic methods. We partitioned tumors according to molecular subtypes, mutations and positivity for HPV. We came to the conclusion that high levels of the histone-like subunits and the “short” NF-YAs isoform are protective in HPV-positive tumors. On the other hand, high levels of the “long” NF-YAl were found in the recently identified aggressive and metastasis-prone cell population undergoing partial epithelial to mesenchymal transition, p-EMT.

Abstract

NF-Y is the CCAAT-binding trimer formed by the histone fold domain (HFD), NF-YB/NF-YC and NF-YA. The CCAAT box is generally prevalent in promoters of “cancer” genes. We reported the overexpression of NF-YA in BRCA, LUAD and LUSC, and of all subunits in HCC. Altered splicing of NF-YA was found in breast and lung cancer. We analyzed RNA-seq datasets of TCGA and cell lines of head and neck squamous cell carcinomas (HNSCC). We partitioned all TCGA data into four subtypes, deconvoluted single-cell RNA-seq of tumors and derived survival curves. The CCAAT box was enriched in the promoters of overexpressed genes. The “short” NF-YAs was overexpressed in all subtypes and the “long” NF-YAl in Mesenchymal. The HFD subunits are overexpressed, except Basal (NF-YB) and Atypical (NF-YC); NF-YAl is increased in p53 mutated tumors. In HPV-positive tumors, high levels of NF-YAs, p16 and ΔNp63 correlate with better prognosis. Deconvolution of single cell RNA-seq (scRNA-seq) found a correlation of NF-YAl with Cancer Associated Fibroblasts (CAFs) and p-EMT cells, a population endowed with metastatic potential. We conclude that overexpression of HFD subunits and NF-YAs is protective in HPV-positive tumors; expression of NF-YAl is largely confined to mutp53 tumors and malignant p-EMT cells.

mTOR Driven Gene Transcription Is Required for Cholesterol Production in Neurons of the Developing Cerebral Cortex

Dysregulated mammalian target of rapamycin (mTOR) activity is associated with various neurodevelopmental disorders ranging from idiopathic autism spectrum disorders (ASD) to syndromes caused by single gene defects. This suggests that maintaining mTOR activity levels in a physiological range is essential for brain development and functioning. Upon activation, mTOR regulates a variety of cellular processes such as cell growth, autophagy, and metabolism. On a molecular level, however, the consequences of mTOR activation in the brain are not well understood. Low levels of cholesterol are associated with a wide variety of neurodevelopmental disorders. We here describe numerous genes of the sterol/cholesterol biosynthesis pathway to be transcriptionally regulated by mTOR complex 1 (mTORC1) signaling in vitro in primary neurons and in vivo in the developing cerebral cortex of the mouse. We find that these genes are shared targets of the transcription factors SREBP, SP1, and NF-Y. Prenatal as well as postnatal mTORC1 inhibition downregulated expression of these genes which directly translated into reduced cholesterol levels, pointing towards a substantial metabolic function of the mTORC1 signaling cascade. Altogether, our results indicate that mTORC1 is an essential transcriptional regulator of the expression of sterol/cholesterol biosynthesis genes in the developing brain. Altered expression of these genes may be an important factor contributing to the pathogenesis of neurodevelopmental disorders associated with dysregulated mTOR signaling.

Live cell dynamics of the NF-Y transcription factor

Transcription factors (TFs) are core players in the control of gene expression, evolutionarily selected to recognise a subset of specific DNA sequences and nucleate the recruitment of the transcriptional machinery. How TFs assemble and move in the nucleus to locate and bind their DNA targets and cause a transcriptional response, remains mostly unclear. NF-Y is a highly conserved, heterotrimeric TF with important roles in both housekeeping and lineage-specific gene expression, functioning as a promoter organiser. Despite a large number of biochemical, structural and genomic studies of NF-Y, there is a lack of experiments in single living cells; therefore, basic assumptions of NF-Y biology remain unproven in vivo. Here we employ a series of dynamic fluorescence microscopy methods (FLIM-FRET, NB, RICS and FRAP) to study NF-Y dynamics and complex formation in live cells. Specifically, we provide quantitative measurement of NF-Y subunit association and diffusion kinetics in the nucleus that collectively suggest NF-Y to move and bind chromatin as a trimeric complex in vivo.

Transcription factor competition at the γ-globin promoters controls hemoglobin switching

BCL11A, the major regulator of fetal hemoglobin (HbF, α2γ2) level, represses γ-globin expression through direct promoter binding in adult erythroid cells in a switch to adult hemoglobin (HbA, α2β2). To uncover how BCL11A initiates repression, we used CRISPR-Cas9, dCas9, dCas9-KRAB and dCas9-VP64 screens to dissect the γ-globin promoters and identified an activator element near the BCL11A-binding site. Using CUT&RUN and base editing, we demonstrate that a proximal CCAAT box is occupied by the activator NF-Y. BCL11A competes with NF-Y binding through steric hindrance to initiate repression. Occupancy of NF-Y is rapidly established following BCL11A depletion, and precedes γ-globin derepression and locus control region (LCR)-globin loop formation. Our findings reveal that the switch from fetal to adult globin gene expression within the >50-kb β-globin gene cluster is initiated by competition between a stage-selective repressor and a ubiquitous activating factor within a remarkably discrete region of the γ-globin promoters.

TopicNet: a framework for measuring transcriptional regulatory network change

Motivation

Recently, many chromatin immunoprecipitation sequencing experiments have been carried out for a diverse group of transcription factors (TFs) in many different types of human cells. These experiments manifest large-scale and dynamic changes in regulatory network connectivity (i.e. network ‘rewiring’), highlighting the different regulatory programs operating in disparate cellular states. However, due to the dense and noisy nature of current regulatory networks, directly comparing the gains and losses of targets of key TFs across cell states is often not informative. Thus, here, we seek an abstracted, low-dimensional representation to understand the main features of network change.

Results

We propose a method called TopicNet that applies latent Dirichlet allocation to extract functional topics for a collection of genes regulated by a given TF. We then define a rewiring score to quantify regulatory-network changes in terms of the topic changes for this TF. Using this framework, we can pinpoint particular TFs that change greatly in network connectivity between different cellular states (such as observed in oncogenesis). Also, incorporating gene expression data, we define a topic activity score that measures the degree to which a given topic is active in a particular cellular state. And we show how activity differences can indicate differential survival in various cancers.

Availability and Implementation

The TopicNet framework and related analysis were implemented using R and all codes are available at https://github.com/gersteinlab/topicnet.

Supplementary information

Supplementary data are available at Bioinformatics online.

NHP2 downregulation counteracts hTR-mediated activation of the DNA damage response at ALT telomeres

About 10% of cancer cells employ the "alternative lengthening of telomeres" (ALT) pathway instead of re-activating the hTERT subunit of human telomerase. The hTR RNA subunit is also abnormally silenced in some ALT⁺ cells not expressing hTERT, suggesting a possible negative non-canonical impact of hTR on ALT. Indeed, we show that ectopically expressed hTR reduces phosphorylation of ssDNA-binding protein RPA (p-RPA^S33 ) at ALT telomeres by promoting the hnRNPA1- and DNA-PK-dependent depletion of RPA. The resulting defective ATR checkpoint signaling at telomeres impairs recruitment of the homologous recombination protein, RAD51. This induces ALT telomere fragility, increases POLD3-dependent C-circle production, and promotes the recruitment of the DNA damage marker 53BP1. In ALT⁺ cells that naturally retain hTR expression, NHP2 H/ACA ribonucleoprotein levels are downregulated, likely in order to restrain DNA damage response (DDR) activation at telomeres through reduced 53BP1 recruitment. This unexpected role of NHP2 is independent from hTR's non-canonical function in modulating telomeric p-RPA^S33 . Collectively, our study shines new light on the interference between telomerase- and ALT-dependent pathways and unravels a crucial role for hTR and NHP2 in DDR regulation at ALT telomeres.

Insights into epithelial cell senescence from transcriptome and secretome analysis of human oral keratinocytes

Senescent cells produce chronic inflammation that contributes to the diseases and debilities of aging. How this process is orchestrated in epithelial cells, the origin of human carcinomas, is poorly understood. We used human normal oral keratinocytes (NOKs) to elucidate senescence programs in a prototype primary mucosal epithelial cell that senesces spontaneously. While NOKs exhibit several typical facets of senescence, they also display distinct characteristics. These include expression of p21WAF1/CIP1 at early passages, making this common marker of senescence unreliable in NOKs. Transcriptome analysis by RNA-seq revealed specific commonalities with and differences from cancer cells, explicating the tumor avoidance role of senescence. Repression of DNA repair genes that correlated with downregulation of E2F1 mRNA and protein was observed for two donors; a divergent result was seen for the third. Using proteomic profiling of soluble (non-vesicular) and extracellular vesicle (EV) associated secretions, we propose additions to the senescence associated secretory phenotype, including HSP60, which localizes to the surface of EVs. Finally, EVs from senescent NOKs activate interferon pathway signaling in THP-1 monocytes in a STING-dependent manner and associate with mitochondrial and nuclear DNA. Our results highlight senescence changes in epithelial cells and how they might contribute to chronic inflammation and age-related diseases.

LEAFY is a pioneer transcription factor and licenses cell reprogramming to floral fate

Master transcription factors reprogram cell fate in multicellular eukaryotes. Pioneer transcription factors have prominent roles in this process because of their ability to contact their cognate binding motifs in closed chromatin. Reprogramming is pervasive in plants, whose development is plastic and tuned by the environment, yet little is known about pioneer transcription factors in this kingdom. Here, we show that the master transcription factor LEAFY (LFY), which promotes floral fate through upregulation of the floral commitment factor APETALA1 (AP1), is a pioneer transcription factor. In vitro, LFY binds to the endogenous AP1 target locus DNA assembled into a nucleosome. In vivo, LFY associates with nucleosome occupied binding sites at the majority of its target loci, including AP1. Upon binding, LFY 'unlocks' chromatin locally by displacing the H1 linker histone and by recruiting SWI/SNF chromatin remodelers, but broad changes in chromatin accessibility occur later. Our study provides a mechanistic framework for patterning of inflorescence architecture and uncovers striking similarities between LFY and animal pioneer transcription factor.

Resetting FLOWERING LOCUS C Expression After Vernalization Is Just Activation in the Early Embryo by a Different Name

The reproductive success of many plants depends on their capacity to respond appropriately to their environment. One environmental cue that triggers flowering is the extended cold of winter, which promotes the transition from vegetative to reproductive growth in a response known as vernalization. In annual plants of the Brassicaceae, the floral repressor, FLOWERING LOCUS C (FLC), is downregulated by exposure to low temperatures. Repression is initiated during winter cold and then maintained as the temperature rises, allowing plants to complete their life cycle during spring and summer. The two stages of FLC repression, initiation and maintenance, are distinguished by different chromatin states at the FLC locus. Initiation involves the removal of active chromatin marks and the deposition of the repressive mark H3K27me3 over a few nucleosomes in the initiation zone, also known as the nucleation region. H3K27me3 then spreads to cover the entire locus, in a replication dependent manner, to maintain FLC repression. FLC is released from repression in the next generation, allowing progeny of a vernalized plant to respond to winter. Activation of FLC in this generation has been termed resetting to denote the restoration of the pre-vernalized state in the progeny of a vernalized plant. It has been assumed that resetting must differ from the activation of FLC expression in progeny of plants that have not experienced winter cold. Considering that there is now strong evidence indicating that chromatin undergoes major modifications during both male and female gametogenesis, it is time to challenge this assumption.

MITO-Luc/GFP zebrafish model to assess spatial and temporal evolution of cell proliferation in vivo

We developed a novel reporter transgenic zebrafish model called MITO-Luc/GFP zebrafish in which GFP and luciferase expression are under the control of the master regulator of proliferation NF-Y. In MITO-Luc/GFP zebrafish it is possible to visualize cell proliferation in vivo by fluorescence and bioluminescence. In this animal model, GFP and luciferase expression occur in early living embryos, becoming tissue specific in juvenile and adult zebrafish. By in vitro and ex vivo experiments we demonstrate that luciferase activity in adult animals occurs in intestine, kidney and gonads, where detectable proliferating cells are located. Further, by time lapse experiments in live embryos, we observed a wave of GFP positive cells following fin clip. In adult zebrafish, in addition to a bright bioluminescence signal on the regenerating tail, an early unexpected signal coming from the kidney occurs indicating not only a fin cell proliferation, but also a systemic response to tissue damage. Finally, we observed that luciferase activity was inhibited by anti-proliferative interventions, i.e. 5FU, cell cycle inhibitors and X-Rays. In conclusion, MITO-Luc/GFP zebrafish is a novel animal model that may be crucial to assess the spatial and temporal evolution of cell proliferation in vivo.

Nonlinear relationship between chromatin accessibility and estradiol-regulated gene expression

Chromatin accessibility is central to basal and inducible gene expression. Through ATAC-seq experiments in estrogen receptor-positive (ER+) breast cancer cell line MCF-7 and integration with multi-omics data, we found estradiol (E2) induced chromatin accessibility changes in a small number of breast cancer-relevant E2-regulated genes. As expected, open chromatin regions associated with E2-inducible gene expression showed enrichment of estrogen response element (ERE) and those associated with E2-repressible gene expression were enriched for ERE, PBX1, and PBX3. While a significant number of open chromatin regions showed pioneer factor FOXA1 occupancy in the absence of E2, E2-treatment further enhanced FOXA1 occupancy suggesting that ER-E2 enhances chromatin occupancy of FOXA1 to a subset of E2-regulated genes. Surprisingly, promoters of 80% and enhancers of 60% of E2-inducible genes displayed closed chromatin configuration both in the absence and presence of E2. Integration of ATAC-seq data with ERα ChIP-seq data revealed that ~40% ERα binding sites in the genome are found in chromatin regions that are not accessible as per ATAC-seq. Such ERα binding regions were enriched for binding sites of multiple nuclear receptors including ER, ESRRB, ERRγ, COUP-TFII (NR2F2), RARα, EAR2 as well as traditional pioneer factors FOXA1 and GATA3. Similar data were also obtained when ERα ChIP-seq data were integrated with MNase-seq and DNase-seq data sets. In summation, our results reveal complex mechanisms of ER-E2 interaction with nucleosomes. Notably, "closed chromatin" configuration as defined by ATAC-seq or by other techniques is not necessarily associated with lack of gene expression and technical limitations may preclude ATAC-seq to demonstrate accessibility of chromatin regions that are bound by ERα.

Structural determinants for NF-Y subunit organization and NF-Y/DNA association in plants

NF-Y transcription factor comprises three subunits: NF-YA, NF-YB and NF-YC. NF-YB and NF-YC dimerize through their histone fold domain (HFD), which can bind DNA in a non-sequence-specific fashion while serving as a scaffold for NF-YA trimerization. Upon trimerization, NF-YA specifically recognizes the CCAAT box sequence on promoters and enhancers. In plants, each NF-Y subunit is encoded by several genes giving rise to hundreds of potential heterotrimeric combinations. In addition, plant NF-YBs and NF-YCs interact with other protein partners to recognize a plethora of genomic motifs, as the CCT protein family that binds CORE sites. The NF-Y subunit organization and its DNA-binding properties, together with the NF-Y HFD capacity to adapt different protein modules, represent plant-specific features that play a key role in development, growth and reproduction. Despite their relevance, these features are still poorly understood at the molecular level. Here, we present the structures of Arabidopsis and rice NF-YB/NF-YC dimers, and of an Arabidopsis NF-Y trimer in complex with the FT CCAAT box, together with biochemical data on NF-Y mutants. The dimeric structures identify the key residues for NF-Y HFD stabilization. The NF-Y/DNA structure and the mutation experiments shed light on HFD trimerization interface properties and the NF-YA sequence appetite for the bases flanking the CCAAT motif. These data explain the logic of plant NF-Y gene expansion: the trimerization adaptability and the flexible DNA-binding rules serve the scopes of accommodating the large number of NF-YAs, CCTs and possibly other NF-Y HFD binding partners and a diverse audience of genomic motifs.

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.

Gene expression profiling in neuronal cells identifies a different type of transcriptome modulated by NF-Y

A heterotrimeric transcription factor NF-Y is crucial for cell-cycle progression in various types of cells. In contrast, studies using NF-YA knockout mice have unveiled its essential role in endoplasmic reticulum (ER) homeostasis in neuronal cells. However, whether NF-Y modulates a different transcriptome to mediate distinct cellular functions remains obscure. Here, we knocked down NF-Y in two types of neuronal cells, neuro2a neuroblastoma cells and mouse brain striatal cells, and performed gene expression profiling. We found that down-regulated genes preferentially contained NF-Y-binding motifs in their proximal promoters, and notably enriched genes related to ER functions rather than those for cell cycle. This contrasts with the profiling data of HeLa and embryonic stem cells in which distinct down-regulation of cell cycle-related genes was observed. Clustering analysis further identified several functional clusters where populations of the down-regulated genes were highly distinct. Further analyses using chromatin immunoprecipitation and RNA-seq data revealed that the transcriptomic difference was not correlated with DNA binding of NF-Y but with splicing of NF-YA. These data suggest that neuronal cells have a different type of transcriptome in which ER-related genes are dominantly modulated by NF-Y, and imply that NF-YA splicing alteration could be involved in this cell type-specific gene modulation.

NF-Y Overexpression in Liver Hepatocellular Carcinoma (HCC)

NF-Y is a pioneer trimeric transcription factor formed by the Histone Fold Domain (HFD) NF-YB/NF-YC subunits and NF-YA. Three subunits are required for DNA binding. CCAAT-specificity resides in NF-YA and transactivation resides in Q-rich domains of NF-YA and NF-YC. They are involved in alternative splicing (AS). We recently showed that NF-YA is overexpressed in breast and lung carcinomas. We report here on the overexpression of all subunits in the liver hepatocellular carcinoma (HCC) TCGA database, specifically the short NF-YAs and NF-YC2 (37 kDa) isoforms. This is observed at all tumor stages, in viral-infected samples and independently from the inflammatory status. Up-regulation of NF-YAs and NF-YC, but not NF-YB, is associated to tumors with mutant p53. We used a deep-learning-based method (DeepCC) to extend the partitioning of the three molecular clusters to all HCC TCGA tumors. In iCluster3, CCAAT is a primary matrix found in promoters of up-regulated genes, and cell-cycle pathways are enriched. Finally, clinical data indicate that, globally, only NF-YAs, but not HFD subunits, correlate with the worst prognosis; in iCluster1 patients, however, all subunits correlate. The data show a difference with other epithelial cancers, in that global overexpression of the three subunits is reported and clinically relevant in a subset of patients; yet, they further reinstate the regulatory role of the sequence-specific subunit.

Transcription factor NFYA promotes G1/S cell cycle transition and cell proliferation by transactivating cyclin D1 and CDK4 in clear cell renal cell carcinoma

NFYA (nuclear transcription factor Y, subunit A) is a CCAAT-binding transcription factor. Accumulating evidence suggests that NFYA plays an important role in breast, ovarian, lung and gastric cancer. However, the role of NFYA in clear cell renal cell carcinoma (ccRCC) remains unclear. In this study, it was discovered that the expression of NFYA is elevated in tissues of ccRCC patient and high NFYA expression is linked to poor overall survival in ccRCC patient. Inhibition of G1/S cell cycle transition and decreased cell proliferation were observed upon NFYA knockdown in ccRCC cells. Moreover, further investigation revealed that NFYA binds directly to the promoter region of both CDK4 and cyclin D1 (CCND1) thus transactivating their expression, resulting in RB phosphorylation and the activation of subsequent E2F pathway activation. Taken together, these findings imply the oncogenic role of NFYA in ccRCC progression and its potential as a target for ccRCC therapy.

Influenza A virus directly modulates mouse eosinophil responses

Allergic asthma and influenza are common respiratory diseases with a high probability of co-occurrence. During the 2009 influenza pandemic, hospitalized patients with influenza experienced lower morbidity if asthma was an underlying condition. We have previously demonstrated that acute allergic asthma protects mice from severe influenza and have implicated eosinophils in the airways of mice with allergic asthma as participants in the antiviral response. However, very little is known about how eosinophils respond to direct exposure to influenza A virus (IAV) or the microenvironment in which the viral burden is high. We hypothesized that eosinophils would dynamically respond to the presence of IAV through phenotypic, transcriptomic, and physiologic changes. Using our mouse model of acute fungal asthma and influenza, we showed that eosinophils in lymphoid tissues were responsive to IAV infection in the lungs and altered surface expression of various markers necessary for cell activation in a niche-specific manner. Siglec-F expression was altered in a subset of eosinophils after virus exposure, and those expressing high Siglec-F were more active (IL-5Rα^hi CD62L^lo ). While eosinophils exposed to IAV decreased their overall transcriptional activity and mitochondrial oxygen consumption, transcription of genes encoding viral recognition proteins, Ddx58 (RIG-I), Tlr3, and Ifih1 (MDA5), were up-regulated. CD8⁺ T cells from IAV-infected mice expanded in response to IAV PB1 peptide-pulsed eosinophils, and CpG methylation in the Tbx21 promoter was reduced in these T cells. These data offer insight into how eosinophils respond to IAV and help elucidate alternative mechanisms by which they regulate antiviral immune responses during IAV infection.

The nuclear oncoprotein Fra-1: a transcription factor knocking on therapeutic applications' door

Among the FOS-related members of the AP-1 dimeric complex, the transcription factor Fra-1, encoded by FOSL1, is crucially involved in human tumor progression and metastasis, thus representing a promising therapeutic target. Here we review the state of the art and discuss the emerging topics and perspectives on FOSL1 and its gene product. First, we summarize the present knowledge on the FOSL1 transcriptional and epigenetic controls, driving Fra-1 accumulation in a variety of human solid tumors. We also present a model on the regulatory interactions between Fra-1, p53, and miRNAs. Then, we outline the multiple roles of Fra-1 posttranslational modifications and transactivation mechanisms of select Fra-1 target genes. In addition to summarizing the Fra-1-dependent gene networks controlling proliferation, survival, and epithelial-mesenchymal transitions (EMT) in multiple cancer cell types, we highlight the roles played by Fra-1 in nonneoplastic cell populations recruited to the tumor microenvironment, and in mouse models of tumorigenesis. Next, we review the prognostic power of the Fra-1-associated gene signatures, and envisage potential strategies aimed at Fra-1 therapeutic inhibition. Finally, we discuss several recent reports showing the emerging roles of Fra-1 in the mechanisms of both resistance and addiction to targeted therapies.

Single-nucleus RNA-seq identifies divergent populations of FSHD2 myotube nuclei

FSHD is characterized by the misexpression of DUX4 in skeletal muscle. Although DUX4 upregulation is thought to be the pathogenic cause of FSHD, DUX4 is lowly expressed in patient samples, and analysis of the consequences of DUX4 expression has largely relied on artificial overexpression. To better understand the native expression profile of DUX4 and its targets, we performed bulk RNA-seq on a 6-day differentiation time-course in primary FSHD2 patient myoblasts. We identify a set of 54 genes upregulated in FSHD2 cells, termed FSHD-induced genes. Using single-cell and single-nucleus RNA-seq on myoblasts and differentiated myotubes, respectively, we captured, for the first time, DUX4 expressed at the single-nucleus level in a native state. We identified two populations of FSHD myotube nuclei based on low or high enrichment of DUX4 and FSHD-induced genes ("FSHD-Lo" and "FSHD Hi", respectively). FSHD-Hi myotube nuclei coexpress multiple DUX4 target genes including DUXA, LEUTX and ZSCAN4, and also upregulate cell cycle-related genes with significant enrichment of E2F target genes and p53 signaling activation. We found more FSHD-Hi nuclei than DUX4-positive nuclei, and confirmed with in situ RNA/protein detection that DUX4 transcribed in only one or two nuclei is sufficient for DUX4 protein to activate target genes across multiple nuclei within the same myotube. DUXA (the DUX4 paralog) is more widely expressed than DUX4, and depletion of DUXA suppressed the expression of LEUTX and ZSCAN4 in late, but not early, differentiation. The results suggest that the DUXA can take over the role of DUX4 to maintain target gene expression. These results provide a possible explanation as to why it is easier to detect DUX4 target genes than DUX4 itself in patient cells and raise the possibility of a self-sustaining network of gene dysregulation triggered by the limited DUX4 expression.

Epigenetic Regulation of Hepatic Lipogenesis: Role in Hepatosteatosis and Diabetes

Hepatosteatosis, which is frequently associated with development of metabolic syndrome and insulin resistance, manifests when triglyceride (TG) input in the liver is greater than TG output, resulting in the excess accumulation of TG. Dysregulation of lipogenesis therefore has the potential to increase lipid accumulation in the liver, leading to insulin resistance and type 2 diabetes. Recently, efforts have been made to examine the epigenetic regulation of metabolism by histone-modifying enzymes that alter chromatin accessibility for activation or repression of transcription. For regulation of lipogenic gene transcription, various known lipogenic transcription factors, such as USF1, ChREBP, and LXR, interact with and recruit specific histone modifiers, directing specificity toward lipogenesis. Alteration or impairment of the functions of these histone modifiers can lead to dysregulation of lipogenesis and thus hepatosteatosis leading to insulin resistance and type 2 diabetes.

Dynamic Transcriptional Responses to Injury of Regenerative and Non-regenerative Cardiomyocytes Revealed by Single-Nucleus RNA Sequencing

The adult mammalian heart is incapable of regeneration following injury. In contrast, the neonatal mouse heart can efficiently regenerate during the first week of life. The molecular mechanisms that mediate the regenerative response and its blockade in later life are not understood. Here, by single-nucleus RNA sequencing, we map the dynamic transcriptional landscape of five distinct cardiomyocyte populations in healthy, injured, and regenerating mouse hearts. We identify immature cardiomyocytes that enter the cell cycle following injury and disappear as the heart loses the ability to regenerate. These proliferative neonatal cardiomyocytes display a unique transcriptional program dependent on nuclear transcription factor Y subunit alpha (NFYa) and nuclear factor erythroid 2-like 1 (NFE2L1) transcription factors, which exert proliferative and protective functions, respectively. Cardiac overexpression of these two factors conferred protection against ischemic injury in mature mouse hearts that were otherwise non-regenerative. These findings advance our understanding of the cellular basis of neonatal heart regeneration and reveal a transcriptional landscape for heart repair following injury.

The Switch from NF-YAl to NF-YAs Isoform Impairs Myotubes Formation

NF-YA, the regulatory subunit of the trimeric transcription factor (TF) NF-Y, is regulated by alternative splicing (AS) generating two major isoforms, “long” (NF-YAl) and “short” (NF-YAs). Muscle cells express NF-YAl. We ablated exon 3 in mouse C2C12 cells by a four-guide CRISPR/Cas9n strategy, obtaining clones expressing exclusively NF-YAs (C2-YAl-KO). C2-YAl-KO cells grow normally, but are unable to differentiate. Myogenin and—to a lesser extent, MyoD— levels are substantially lower in C2-YAl-KO, before and after differentiation. Expression of the fusogenic Myomaker and Myomixer genes, crucial for the early phases of the process, is not induced. Myomaker and Myomixer promoters are bound by MyoD and Myogenin, and Myogenin overexpression induces their expression in C2-YAl-KO. NF-Y inactivation reduces MyoD and Myogenin, but not directly: the Myogenin promoter is CCAAT-less, and the canonical CCAAT of the MyoD promoter is not bound by NF-Y in vivo. We propose that NF-YAl, but not NF-YAs, maintains muscle commitment by indirectly regulating Myogenin and MyoD expression in C2C12 cells. These experiments are the first genetic evidence that the two NF-YA isoforms have functionally distinct roles.

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 in Lung Cancer: LUAD

The trimeric transcription factor (TF) NF-Y regulates the CCAAT box, a DNA element enriched in promoters of genes overexpressed in many types of cancer. The regulatory NF-YA is present in two major isoforms, NF-YAl ("long") and NF-YAs ("short"). There is growing indication that NF-YA levels are increased in tumors. Here, we report interrogation of RNA-Seq TCGA (The Cancer Genome Atlas)-all 576 samples-and GEO (Gene Expression Ominibus) datasets of lung adenocarcinoma (LUAD). NF-YAs is overexpressed in the three subtypes, proliferative, inflammatory, and TRU (terminal respiratory unit). CCAAT is enriched in promoters of tumor differently expressed genes (DEG) and in the proliferative/inflammatory intersection, matching with KEGG (Kyoto Encyclopedia of Genes and Genomes) terms cell-cycle and signaling. Increasing levels of NF-YAs are observed from low to high CpG island methylator phenotypes (CIMP). We identified 166 genes overexpressed in LUAD cell lines with low NF-YAs/NF-YAl ratios: applying this centroid to TCGA samples faithfully predicted tumors' isoform ratio. This signature lacks CCAAT in promoters. Finally, progression-free intervals and hazard ratios concurred with the worst prognosis of patients with either a low or high NF-YAs/NF-YAl ratio. In conclusion, global overexpression of NF-YAs is documented in LUAD and is associated with aggressive tumor behavior; however, a similar prognosis is recorded in tumors with high levels of NF-YAl and overexpressed CCAAT-less genes.

NF-YA Overexpression in Lung Cancer: LUSC

The CCAAT box is recognized by the trimeric transcription factor NF-Y, whose NF-YA subunit is present in two major splicing isoforms, NF-YAl (“long”) and NF-YAs (“short”). Little is known about the expression levels of NF-Y subunits in tumors, and nothing in lung cancer. By interrogating RNA-seq TCGA and GEO datasets, we found that, unlike NF-YB/NF-YC, NF-YAs is overexpressed in lung squamous cell carcinomas (LUSC). The ratio of the two isoforms changes from normal to cancer cells, with NF-YAs becoming predominant in the latter. NF-YA increased expression correlates with common proliferation markers. We partitioned all 501 TCGA LUSC tumors in the four molecular cohorts and verified that NF-YAs is similarly overexpressed. We analyzed global and subtype-specific RNA-seq data and found that CCAAT is the most abundant DNA matrix in promoters of genes overexpressed in all subtypes. Enriched Gene Ontology terms are cell-cycle and signaling. Survival curves indicate a worse clinical outcome for patients with increasing global amounts of NF-YA; same with hazard ratios with very high and, surprisingly, very low NF-YAs/NF-YAl ratios. We then analyzed gene expression in this latter cohort and identified a different, pro-migration signature devoid of CCAAT. We conclude that overexpression of the NF-Y regulatory subunit in LUSC has the scope of increasing CCAAT-dependent, proliferative (NF-YAs^high) or CCAAT-less, pro-migration (NF-YAl^high) genes. The data further reinstate the importance of analysis of single isoforms of TFs involved in tumor development.

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.

Overexpression and alternative splicing of NF-YA in breast cancer

NF-Y is a CCAAT-binding trimeric transcription factor, whose regulome, interactome and oncogenic potential point to direct involvement in cellular transformation. Yet little is known about the levels of NF-Y subunits in tumors. We focused on breast carcinomas, and analyzed RNA-Seq datasets of TCGA and 54 BRCA cell lines at gene and isoforms level. We partitioned all tumors in the four major subclasses. NF-YA, but not histone-fold subunits NF-YB/NF-YC, is globally overexpressed, correlating with the proliferative Ki67 marker and a common set of 840 genes, with cell-cycle, metabolism GO terms. Their promoters are enriched in NF-Y, GC-rich and E2F sites. Surprisingly, there is an isoform switch, with the “short” isoform -NF-YAs- becoming predominant in tumors. E2F genes are also overexpressed in BRCA, but no switch in isoforms is observed. In Basal-like Claudin^low cell lines and tumors, expression of NF-YAl -long- isoform is high, together with 11 typical EMT markers and low levels of basal Keratins. Analysis of Progression-Free-Intervals indicates that tumors with unbalance of NF-YA isoforms ratios have worst clinical outcomes. The data suggest that NF-YA overexpression increases CCAAT-dependent, pro-growth genes in BRCA. NF-YAs is associated with a proliferative signature, but high levels of NF-YAl signal loss of epithelial features, EMT and acquisition of a more aggressive behavior in a subset of Claudin^low Basal-like tumors.

NF-Y controls fidelity of transcription initiation at gene promoters through maintenance of the nucleosome-depleted region

Faithful transcription initiation is critical for accurate gene expression, yet the mechanisms underlying specific transcription start site (TSS) selection in mammals remain unclear. Here, we show that the histone-fold domain protein NF-Y, a ubiquitously expressed transcription factor, controls the fidelity of transcription initiation at gene promoters in mouse embryonic stem cells. We report that NF-Y maintains the region upstream of TSSs in a nucleosome-depleted state while simultaneously protecting this accessible region against aberrant and/or ectopic transcription initiation. We find that loss of NF-Y binding in mammalian cells disrupts the promoter chromatin landscape, leading to nucleosomal encroachment over the canonical TSS. Importantly, this chromatin rearrangement is accompanied by upstream relocation of the transcription pre-initiation complex and ectopic transcription initiation. Further, this phenomenon generates aberrant extended transcripts that undergo translation, disrupting gene expression profiles. These results suggest NF-Y is a central player in TSS selection in metazoans and highlight the deleterious consequences of inaccurate transcription initiation.

The mechanisms underlying specific TSS selection in mammals remain unclear. Here the authors show that the ubiquitously expressed transcription factor NF-Y regulate fidelity of transcription initiation at gene promoters, maintaining the region upstream of TSSs in a nucleosome-depleted state, while protecting this region from ectopic transcription initiation.

Inflammatory regulatory network mediated by the joint action of NF-kB, STAT3, and AP-1 factors is involved in many human cancers

Using an inducible, inflammatory model of breast cellular transformation, we describe the transcriptional regulatory network mediated by STAT3, NF-κB, and AP-1 factors on a genomic scale. These proinflammatory regulators form transcriptional complexes that directly regulate the expression of hundreds of genes in oncogenic pathways via a positive feedback loop. This transcriptional feedback loop and associated network functions to various extents in many types of cancer cells and patient tumors, and it is the basis for a cancer inflammation index that defines cancer types by functional criteria. We identify a network of noninflammatory genes whose expression is well correlated with the cancer inflammatory index. Conversely, the cancer inflammation index is negatively correlated with the expression of genes involved in DNA metabolism, and transformation is associated with genome instability. We identify drugs whose efficacy in cell lines is correlated with the cancer inflammation index, suggesting the possibility of using this index for personalized cancer therapy. Inflammatory tumors are preferentially associated with infiltrating immune cells that might be recruited to the site of the tumor via inflammatory molecules produced by the cancer cells.

The phosphorylatable Ser320 of NF-YA is involved in DNA binding of the NF-Y trimer

Nuclear factor Y (NF-Y) is a transcription factor trimer binding to the functionally important CCAAT box, present in promoters of growth-promoting and cell cycle-regulated genes. The regulatory nuclear factor YA (NF-YA) subunit confers sequence-specificity to the histone-like nuclear factor YB/YC dimer. NF-YA harbors 2 serines-Ser320 and Ser326-shown to be phosphorylated by cyclin-dependent kinase 2. High-throughput proteomics data indicate that they are phosphorylated in vivo. Specifically, Ser320 makes structural contacts with the DNA phosphate backbone; Ser320-P is the major NF-YA phosphorylation isoform following overexpression in HeLa cells, increasing upon mitotic arrest. EMSA with recombinant Ala and Glu mutants confirm a role of Ser320, but not Ser326, in stabilization of DNA binding. Transactivation assays of the CCAAT-dependent MDR1 and RHOB promoters show loss in transcription function for Ser320Glu and Ser320Ala NF-YA mutants. Phylogenetic analysis of NF-YA proteins indicates that Ser320 is indeed evolutionarily conserved. We conclude that phosphorylation of this residue belongs to the core mechanisms of DNA-binding control, possibly driven by the necessity to unfasten binding of or to evict NF-Y from CCAAT sites under specific conditions of growth regulation.-Bernardini, A., Lorenzo, M., Nardini, M., Mantovani, R., Gnesutta, N. The phosphorylatable Ser320 of NF-YA is involved in DNA binding of the NF-Y trimer.