PAVIS: a tool for Peak Annotation and Visualization

RAG1/2 induces double-stranded DNA breaks at non-Ig loci in the proximity of single sequence repeats in developing B cells

In developing B cells, V(D)J gene recombination is initiated by the RAG1/2 endonuclease complex, introducing double-stranded DNA breaks (DSBs) in V, D, and J genes and resulting in the formation of the hypervariable parts of immunoglobulins (Ig). Persistent or aberrant RAG1/2 targeting is a potential threat to genome integrity. While RAG1 and RAG2 have been shown to bind various regions genome-wide, the in vivo off-target DNA damage instigated by RAG1/2 endonuclease remains less well understood. In the current study, we identified regions containing RAG1/2-induced DNA breaks in mouse pre-B cells on a genome-wide scale using a global DNA DSB detection strategy. We detected 1489 putative RAG1/2-dependent DSBs, most of which were located outside the Ig loci. DNA sequence motif analysis showed a specific enrichment of RAG1/2-induced DNA DSBs at GA- and CA-repeats and GC-rich motifs. These findings provide further insights into RAG1/2 off-target activity. The ability of RAG1/2 to introduce DSBs on the non-Ig loci during the endogenous V(D)J recombination emphasizes its genotoxic potential in developing lymphocytes.

ERK signalling eliminates Nanog and maintains Oct4 to drive the formative pluripotency transition

Naïve epiblast cells in the embryo and pluripotent stem cells in vitro undergo developmental progression to a formative state competent for lineage specification. During this transition, transcription factors and chromatin are rewired to encode new functional features. Here, we examine the role of mitogen-activated protein kinase (ERK1/2) signalling in pluripotent state transition. We show that a primary consequence of ERK activation in mouse embryonic stem cells is elimination of Nanog, which precipitates breakdown of the naïve state gene regulatory network. Variability in pERK dynamics results in heterogeneous loss of Nanog and metachronous state transition. Knockdown of Nanog allows exit without ERK activation. However, transition to formative pluripotency does not proceed and cells collapse to an indeterminate identity. This outcome is due to failure to maintain expression of the central pluripotency factor Oct4. Thus, during formative transition ERK signalling both dismantles the naïve state and preserves pluripotency. These results illustrate how a single signalling pathway can both initiate and secure transition between cell states.

DiPRO1 distinctly reprograms muscle and mesenchymal cancer cells

We have recently identified the uncharacterized ZNF555 protein as a component of a productive complex involved in the morbid function of the 4qA locus in facioscapulohumeral dystrophy. Subsequently named DiPRO1 (Death, Differentiation, and PROliferation related PROtein 1), our study provides substantial evidence of its role in the differentiation and proliferation of human myoblasts. DiPRO1 operates through the regulatory binding regions of SIX1, a master regulator of myogenesis. Its relevance extends to mesenchymal tumors, such as rhabdomyosarcoma (RMS) and Ewing sarcoma, where DiPRO1 acts as a repressor via the epigenetic regulators TIF1B and UHRF1, maintaining methylation of cis-regulatory elements and gene promoters. Loss of DiPRO1 mimics the host defense response to virus, awakening retrotransposable repeats and the ZNF/KZFP gene family. This enables the eradication of cancer cells, reprogramming the cellular decision balance towards inflammation and/or apoptosis by controlling TNF-α via NF-kappaB signaling. Finally, our results highlight the vulnerability of mesenchymal cancer tumors to si/shDiPRO1-based nanomedicines, positioning DiPRO1 as a potential therapeutic target.

Control of endothelial cell function and arteriogenesis by MEG3:EZH2 epigenetic regulation of integrin expression

Epigenetic processes involving long non-coding RNAs regulate endothelial gene expression. However, the underlying regulatory mechanisms causing endothelial dysfunction remain to be elucidated. Enhancer of zeste homolog 2 (EZH2) is an important rheostat of histone H3K27 trimethylation (H3K27me3) that represses endothelial targets, but EZH2 RNA binding capacity and EZH2:RNA functional interactions have not been explored in post-ischemic angiogenesis. We used formaldehyde/UV-assisted crosslinking ligation and sequencing of hybrids and identified a new role for maternally expressed gene 3 (MEG3). MEG3 formed the predominant RNA:RNA hybrid structures in endothelial cells. Moreover, MEG3:EZH2 assists recruitment onto chromatin. By EZH2-chromatin immunoprecipitation, following MEG3 depletion, we demonstrated that MEG3 controls recruitment of EZH2/H3K27me3 onto integrin subunit alpha4 (ITGA4) promoter. Both MEG3 knockdown or EZH2 inhibition (A-395) promoted ITGA4 expression and improved endothelial cell migration and adhesion to fibronectin in vitro. The A-395 inhibitor re-directed MEG3-assisted chromatin remodeling, offering a direct therapeutic benefit by increasing endothelial function and resilience. This approach subsequently increased the expression of ITGA4 in arterioles following ischemic injury in mice, thus promoting arteriogenesis. Our findings show a context-specific role for MEG3 in guiding EZH2 to repress ITGA4. Novel therapeutic strategies could antagonize MEG3:EZH2 interaction for pre-clinical studies.

DNA double-stranded breaks, a hallmark of aging, defined at the nucleotide resolution, are increased and associated with transcription in the cardiac myocytes in LMNA-cardiomyopathy

AIMS

An intrinsic feature of gene transcription is the formation of DNA superhelices near the transcription bubble, which are resolved upon induction of transient double-stranded breaks (DSBs) by topoisomerases. Unrepaired DSBs are pathogenic as they lead to cell cycle arrest, senescence, inflammation, and organ dysfunction. We posit that DSBs would be more prevalent at the genomic sites that are associated with gene expression. The objectives were to identify and characterize genome-wide DSBs at the nucleotide resolution and determine the association of DSBs with transcription in cardiac myocytes.

METHODS AND RESULTS

We identified the genome-wide DSBs in ∼1 million cardiac myocytes per heart in three wild-type and three myocyte-specific LMNA-deficient (Myh6-Cre:LmnaF/F) mice by END-Sequencing. The prevalence of DSBs was 0.8% and 2.2% in the wild-type and Myh6-Cre:LmnaF/F myocytes, respectively. The END-Seq signals were enriched for 8 and 6764 DSBs in the wild-type and Myh6-Cre:LmnaF/F myocytes, respectively (q < 0.05). The DSBs were preferentially localized to the gene regions, transcription initiation sites, cardiac transcription factor motifs, and the G quadruplex forming structures. Because LMNA regulates transcription through the lamin-associated domains (LADs), we defined the LADs in cardiac myocytes by a Cleavage Under Targets & Release Using Nuclease (CUT&RUN) assay (N = 5). On average there were 818 LADs per myocyte. Constitutive LADs (cLADs), defined as LADs that were shared by at least three genomes (N = 2572), comprised about a third of the mouse cardiac myocyte genomes. Transcript levels of the protein-coding genes located at the cLADs (N = 3975) were ∼16-fold lower than those at the non-LAD regions (N = ∼17 778). The prevalence of DSBs was higher in the non-LAD as compared to the cLAD regions. Likewise, DSBs were more common in the loss-of-LAD regions, defined as the genomic regions in the Myh6-Cre:LmnaF/F that were juxtaposed to the LAD regions in the wild-type myocytes.

CONCLUSION

To our knowledge, this is the first identification of the DSBs, at the nucleotide resolution in the cardiovascular system. The prevalence of DSBs was higher in the genomic regions associated with transcription. Because transcription is pervasive, DSBs are expected to be common and pathogenic in various states and aging.

Machine learning unveils RNA polymerase II binding as a predictor for SMAD2-dependent transcription dynamics in response to Actvin signalling

The transforming growth factor-β (TGF-β) superfamily, including Nodal and Activin, plays a critical role in various cellular processes. Understanding the intricate regulation and gene expression dynamics of TGF-β signalling is of interest due to its diverse biological roles. A machine learning approach is used to predict gene expression patterns induced by Activin using features, such as histone modifications, RNA polymerase II binding, SMAD2-binding, and mRNA half-life. RNA sequencing and ChIP sequencing datasets were analysed and differentially expressed SMAD2-binding genes were identified. These genes were classified into activated and repressed categories based on their expression patterns. The predictive power of different features and combinations was evaluated using logistic regression models and their performances were assessed. Results showed that RNA polymerase II binding was the most informative feature for predicting the expression patterns of SMAD2-binding genes. The authors provide insights into the interplay between transcriptional regulation and Activin signalling and offers a computational framework for predicting gene expression patterns in response to cell signalling.

Risk SNP in a transcript of RP11-638I2.4 increases lncRNA-YY1 interaction and pancreatic cancer susceptibility

Tens of thousands of long non-coding RNAs (lncRNAs) have been identified through RNA-seq analysis, but the biological and pathological significance remains unclear. By integrating the genome-wide lncRNA data with a cross-ancestry meta-analysis of PDAC GWASs, we depicted a comprehensive atlas of pancreatic ductal adenocarcinoma (PDAC)-associated lncRNAs, containing 1,204 lncRNA (445 novel lncRNAs and 759 GENCODE annotated lncRNAs) and 4,368 variants. Furthermore, we found that PDAC-associated lncRNAs could function by altering chromatin activity, transcription factors, and RNA-binding proteins binding affinity. Importantly, genetic variants linked to PDAC are preferentially found at PDAC-associated lncRNA regions, supporting the biological and clinical relevance of PDAC-associated lncRNAs. Finally, we prioritized a novel transcript (MICT00000110172.1) of RP11-638I2.4 as a potential tumor promoter. MICT00000110172.1 is able to reinforce the interaction with YY1, which could reverse the effect of YY1 on pancreatic cancer cell cycle arrest to promote the pancreatic cancer growth. G > A change at rs2757535 in the second exon of MICT00000110172.1 induces a spatial structural change and creates a target region for YY1 binding, which enforces the effect of MICT00000110172.1 in an allele-specific manner, and thus confers susceptibility to tumorigenesis. In summary, our results extend the repertoire of PDAC-associated lncRNAs that could act as a starting point for future functional explorations, and the identification of lncRNA-based target therapy.

Chromatin architectural factor CTCF is essential for progesterone-dependent uterine maturation

Receptors for estrogen and progesterone frequently interact, via Cohesin/CTCF loop extrusion, at enhancers distal from regulated genes. Loss-of-function CTCF mutation in >20% of human endometrial tumors indicates its importance in uterine homeostasis. To better understand how CTCF-mediated enhancer-gene interactions impact endometrial development and function, the Ctcf gene was selectively deleted in female reproductive tissues of mice. Prepubertal Ctcfd/d uterine tissue exhibited a marked reduction in the number of uterine glands compared to those without Ctcf deletion (Ctcff/f mice). Post-pubertal Ctcfd/d uteri were hypoplastic with significant reduction in both the amount of the endometrial stroma and number of glands. Transcriptional profiling revealed increased expression of stem cell molecules Lif, EOMES, and Lgr5, and enhanced inflammation pathways following Ctcf deletion. Analysis of the response of the uterus to steroid hormone stimulation showed that CTCF deletion affects a subset of progesterone-responsive genes. This finding indicates (1) Progesterone-mediated signaling remains functional following Ctcf deletion and (2) certain progesterone-regulated genes are sensitive to Ctcf deletion, suggesting they depend on gene-enhancer interactions that require CTCF. The progesterone-responsive genes altered by CTCF ablation included Ihh, Fst, and Errfi1. CTCF-dependent progesterone-responsive uterine genes enhance critical processes including anti-tumorigenesis, which is relevant to the known effectiveness of progesterone in inhibiting progression of early-stage endometrial tumors. Overall, our findings reveal that uterine Ctcf plays a key role in progesterone-dependent expression of uterine genes underlying optimal post-pubertal uterine development.

Ezh2 emerges as an epigenetic checkpoint regulator during monocyte differentiation limiting cardiac dysfunction post-MI

Epigenetic regulation of histone H3K27 methylation has recently emerged as a key step during alternative immunoregulatory M2-like macrophage polarization; known to impact cardiac repair after Myocardial Infarction (MI). We hypothesized that EZH2, responsible for H3K27 methylation, could act as an epigenetic checkpoint regulator during this process. We demonstrate for the first time an ectopic EZH2, and putative, cytoplasmic inactive localization of the epigenetic enzyme, during monocyte differentiation into M2 macrophages in vitro as well as in immunomodulatory cardiac macrophages in vivo in the post-MI acute inflammatory phase. Moreover, we show that pharmacological EZH2 inhibition, with GSK-343, resolves H3K27 methylation of bivalent gene promoters, thus enhancing their expression to promote human monocyte repair functions. In line with this protective effect, GSK-343 treatment accelerated cardiac inflammatory resolution preventing infarct expansion and subsequent cardiac dysfunction in female mice post-MI in vivo. In conclusion, our study reveals that pharmacological epigenetic modulation of cardiac-infiltrating immune cells may hold promise to limit adverse cardiac remodeling after MI.

p53 governs an AT1 differentiation programme in lung cancer suppression

Lung cancer is the leading cause of cancer deaths worldwide1. Mutations in the tumour suppressor gene TP53 occur in 50% of lung adenocarcinomas (LUADs) and are linked to poor prognosis1-4, but how p53 suppresses LUAD development remains enigmatic. We show here that p53 suppresses LUAD by governing cell state, specifically by promoting alveolar type 1 (AT1) differentiation. Using mice that express oncogenic Kras and null, wild-type or hypermorphic Trp53 alleles in alveolar type 2 (AT2) cells, we observed graded effects of p53 on LUAD initiation and progression. RNA sequencing and ATAC sequencing of LUAD cells uncovered a p53-induced AT1 differentiation programme during tumour suppression in vivo through direct DNA binding, chromatin remodelling and induction of genes characteristic of AT1 cells. Single-cell transcriptomics analyses revealed that during LUAD evolution, p53 promotes AT1 differentiation through action in a transitional cell state analogous to a transient intermediary seen during AT2-to-AT1 cell differentiation in alveolar injury repair. Notably, p53 inactivation results in the inappropriate persistence of these transitional cancer cells accompanied by upregulated growth signalling and divergence from lung lineage identity, characteristics associated with LUAD progression. Analysis of Trp53 wild-type and Trp53-null mice showed that p53 also directs alveolar regeneration after injury by regulating AT2 cell self-renewal and promoting transitional cell differentiation into AT1 cells. Collectively, these findings illuminate mechanisms of p53-mediated LUAD suppression, in which p53 governs alveolar differentiation, and suggest that tumour suppression reflects a fundamental role of p53 in orchestrating tissue repair after injury.

Light and temperature regulate m6A-RNA modification to regulate growth in plants

N6-methyladenosine is a highly dynamic, abundant mRNA modification which is an excellent potential mechanism for fine tuning gene expression. Plants adapt to their surrounding light and temperature environment using complex gene regulatory networks. The role of m6A in controlling gene expression in response to variable environmental conditions has so far been unexplored. Here, we map the transcriptome-wide m6A landscape under various light and temperature environments. Identified m6A-modifications show a highly specific spatial distribution along transcripts with enrichment occurring in 5'UTR regions and around transcriptional end sites. We show that the position of m6A modifications on transcripts might influence cellular transcript localization and the presence of m6A-modifications is associated with alternative polyadenylation, a process which results in multiple RNA isoforms with varying 3'UTR lengths. RNA with m6A-modifications exhibit a higher preference for shorter 3'UTRs. These shorter 3'UTR regions might directly influence transcript abundance and localization by including or excluding cis-regulatory elements. We propose that environmental stimuli might change the m6A landscape of plants as one possible way of fine tuning gene regulation through alternative polyadenylation and transcript localization.

A Systems Biology Approach on the Regulatory Footprint of Human Endogenous Retroviruses (HERVs)

Human endogenous retroviruses (HERVs) are a family of endogenous retroviruses that comprise the ~8.93% of the human genome sequence, with a high proportion being human specific. The recent expansion of repeated HERV sequences has offered a framework for genetic and epigenetic innovation. In the current report, a systematic approach is implemented to catalogue regulatory elements within HERVs, as a roadmap to potential functions of HERV sequences in gene networks. ENCODE Project has offered a wealth of epigenetic data based on omics technologies. I analyzed the presence of HERV sequences on consensus cis-regulatory elements (cCREs) from ENCODE data. On the one side, HERVs are in 1 out of 9 cCREs (>100.000 cCREs in total), dispersed within the genome and present in cis-regulatory regions of ~81% of human genes, as calculated following gene enrichment analysis. On the other side, promoter-associated HERV cCREs are present adjacent to (in a 200 bp window) the transcription start sites of 256 human genes. Regulatory network production, followed by centrality analysis led to the discovery of 90 core genes containing HERV-associated promoters. Pathway analysis on the core network genes and their immediate neighbors revealed a regulatory footprint that, among others, is associated with inflammation, chemokine signaling and response to viral infection. Collectively, these results support the concept that the expansion of regulatory sequences derived from HERVs is critical for epigenetic innovation that may have wired together genes into novel transcriptional networks with critical roles in cellular physiology and pathology.

Exploring Epigenomic Datasets by ChIPseeker

In many aspects of life, epigenetics, or the altering of phenotype without changes in sequences, play an essential role in biological function. A vast number of epigenomic datasets are emerging as a result of the advent of next-generation sequencing. Annotation, comparison, visualization, and interpretation of epigenomic datasets remain key aspects of computational biology. ChIPseeker is a Bioconductor package for performing these analyses among variable epigenomic datasets. The fundamental functions of ChIPseeker, including data preparation, annotation, comparison, and visualization, are explained in this article. ChIPseeker is a freely available open-source package that may be found at https://www.bioconductor.org/packages/ChIPseeker. © 2022 Wiley Periodicals LLC. Basic Protocol 1: ChIPseeker and epigenomic dataset preparation Basic Protocol 2: Annotation of epigenomic datasets Basic Protocol 3: Comparison of epigenomic datasets Basic Protocol 4: Visualization of annotated results Basic Protocol 5: Functional analysis of epigenomic datasets Basic Protocol 6: Genome-wide and locus-specific distribution of epigenomic datasets Basic Protocol 7: Heatmaps and metaplots of epigenomic datasets.

The Estrogen Receptor α Cistrome in Human Endometrium and Epithelial Organoids

Endometrial health is affected by molecular processes that underlie estrogen responses. We assessed estrogen regulation of endometrial function by integrating the estrogen receptor α (ESR1) cistromes and transcriptomes of endometrial biopsies taken from the proliferative and mid-secretory phases of the menstrual cycle together with hormonally stimulated endometrial epithelial organoids. The cycle stage-specific ESR1 binding sites were determined by chromatin immunoprecipitation and next-generation sequencing and then integrated with changes in gene expression from RNA sequencing data to infer candidate ESR1 targets in normal endometrium. Genes with ESR1 binding in whole endometrium were enriched for chromatin modification and regulation of cell proliferation. The distribution of ESR1 binding sites in organoids was more distal from gene promoters when compared to primary endometrium and was more similar to the proliferative than the mid-secretory phase ESR1 cistrome. Inferred organoid estrogen/ESR1 candidate target genes affected formation of cellular protrusions and chromatin modification. Comparison of signaling effected by candidate ESR1 target genes in endometrium vs organoids reveals enrichment of both overlapping and distinct responses. Our analysis of the ESR1 cistromes and transcriptomes from endometrium and organoids provides important resources for understanding how estrogen affects endometrial health and function.

The BASP1 transcriptional corepressor modifies chromatin through lipid-dependent and lipid-independent mechanisms

The transcriptional corepressor BASP1 requires N-terminal myristoylation for its activity and functions through interactions with nuclear lipids. Here we determine the role of BASP1 lipidation in histone modification and the modulation of chromatin accessibility. We find that the removal of the active histone modifications H3K9ac and H3K4me3 by BASP1 requires the N-terminal myristoylation of BASP1. In contrast, the placement of the repressive histone modification, H3K27me3, by BASP1 does not require BASP1 lipidation. RNA-seq and ATAC-seq analysis finds that BASP1 regulates the activity of multiple transcription factors and induces extensive changes in chromatin accessibility. We find that ∼50% of BASP1 target genes show lipidation-dependent chromatin compaction and transcriptional repression. Our results suggest that BASP1 elicits both lipid-dependent and lipid-independent functions in histone modification and transcriptional repression. In accordance with this, we find that the tumor suppressor activity of BASP1 is also partially dependent on its myristoylation.

Progesterone Signaling in Endometrial Epithelial Organoids

For pregnancy to be established, uterine cells respond to the ovarian hormones, estrogen, and progesterone, via their nuclear receptors, the estrogen receptor (ESR1) and progesterone receptor (PGR). ESR1 and PGR regulate genes by binding chromatin at genes and at distal enhancer regions, which interact via dynamic 3-dimensional chromatin structures. Endometrial epithelial cells are the initial site of embryo attachment and invasion, and thus understanding the processes that yield their receptive state is important. Here, we cultured and treated organoids derived from human epithelial cells, isolated from endometrial biopsies, with estrogen and progesterone and evaluated their transcriptional profiles, their PGR cistrome, and their chromatin conformation. Progesterone attenuated estrogen-dependent gene responses but otherwise minimally impacted the organoid transcriptome. PGR ChIPseq peaks were co-localized with previously described organoid ESR1 peaks, and most PGR and ESR1 peaks were in B (inactive) compartment regions of chromatin. Significantly more ESR1 peaks were assigned to estrogen-regulated genes by considering chromatin loops identified using HiC than were identified using ESR1 peak location relative to closest genes. Overall, the organoids model allowed a definition of the chromatin regulatory components governing hormone responsiveness.

ASC proneural factors are necessary for chromatin remodeling during neuroectodermal to neuroblast fate transition to ensure the timely initiation of the neural stem cell program

Background

In both Drosophila and mammals, the achaete-scute (ASC/ASCL) proneural bHLH transcription factors are expressed in the developing central and peripheral nervous systems, where they function during specification and maintenance of the neural stem cells in opposition to Notch signaling. In addition to their role in nervous system development, ASC transcription factors are oncogenic and exhibit chromatin reprogramming activity; however, the impact of ASC on chromatin dynamics during neural stem cell generation remains elusive. Here, we investigate the chromatin changes accompanying neural commitment using an integrative genetics and genomics methodology.

Results

We found that ASC factors bind equally strongly to two distinct classes of cis-regulatory elements: open regions remodeled earlier during maternal to zygotic transition by Zelda and less accessible, Zelda-independent regions. Both classes of cis-elements exhibit enhanced chromatin accessibility during neural specification and correlate with transcriptional regulation of genes involved in a variety of biological processes necessary for neuroblast function/homeostasis. We identified an ASC-Notch regulated TF network that includes likely prime regulators of neuroblast function. Using a cohort of ASC target genes, we report that ASC null neuroblasts are defectively specified, remaining initially stalled, unable to divide, and lacking expression of many proneural targets. When mutant neuroblasts eventually start proliferating, they produce compromised progeny. Reporter lines driven by proneural-bound enhancers display ASC dependency, suggesting that the partial neuroblast identity seen in the absence of ASC genes is likely driven by other, proneural-independent, cis-elements. Neuroblast impairment and the late differentiation defects of ASC mutants are corrected by ectodermal induction of individual ASC genes but not by individual members of the TF network downstream of ASC. However, in wild-type embryos, the induction of individual members of this network induces CNS hyperplasia, suggesting that they synergize with the activating function of ASC to consolidate the chromatin dynamics that promote neural specification.

Conclusions

We demonstrate that ASC proneural transcription factors are indispensable for the timely initiation of the neural stem cell program at the chromatin level by regulating a large number of enhancers in the vicinity of neural genes. This early chromatin remodeling is crucial for both neuroblast homeostasis as well as future progeny fidelity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01300-8.

R-loop mapping and characterization during Drosophila embryogenesis reveals developmental plasticity in R-loop signatures

R-loops are involved in transcriptional regulation, DNA and histone post-translational modifications, genome replication and genome stability. To what extent R-loop abundance and genome-wide localization is actively regulated during metazoan embryogenesis is unknown. Drosophila embryogenesis provides a powerful system to address these questions due to its well-characterized developmental program, the sudden onset of zygotic transcription and available genome-wide data sets. Here, we measure the overall abundance and genome localization of R-loops in early and late-stage embryos relative to Drosophila cultured cells. We demonstrate that absolute R-loop levels change during embryogenesis and that RNaseH1 catalytic activity is critical for embryonic development. R-loop mapping by strand-specific DRIP-seq reveals that R-loop localization is plastic across development, both in the genes which form R-loops and where they localize relative to gene bodies. Importantly, these changes are not driven by changes in the transcriptional program. Negative GC skew and absolute changes in AT skew are associated with R-loop formation in Drosophila. Furthermore, we demonstrate that while some chromatin binding proteins and histone modifications such as H3K27me3 are associated with R-loops throughout development, other chromatin factors associated with R-loops in a developmental specific manner. Our findings highlight the importance and developmental plasticity of R-loops during Drosophila embryogenesis.

Jasmonates and Histone deacetylase 6 activate Arabidopsis genome-wide histone acetylation and methylation during the early acute stress response

BACKGROUND

Jasmonates (JAs) mediate trade-off between responses to both biotic and abiotic stress and growth in plants. The Arabidopsis thaliana HISTONE DEACETYLASE 6 is part of the CORONATINE INSENSITIVE1 receptor complex, co-repressing the HDA6/COI1-dependent acetic acid-JA pathway that confers plant drought tolerance. The decrease in HDA6 binding to target DNA mirrors histone H4 acetylation (H4Ac) changes during JA-mediated drought response, and mutations in HDA6 also cause depletion in the constitutive repressive marker H3 lysine 27 trimethylation (H3K27me3). However, the genome-wide effect of HDA6 on H4Ac and much of the impact of JAs on histone modifications and chromatin remodelling remain elusive.

RESULTS

We performed high-throughput ChIP-Seq on the HDA6 mutant, axe1-5, and wild-type plants with or without methyl jasmonate (MeJA) treatment to assess changes in active H4ac and repressive H3K27me3 histone markers. Transcriptional regulation was investigated in parallel by microarray analysis in the same conditions. MeJA- and HDA6-dependent histone modifications on genes for specialized metabolism; linolenic acid and phenylpropanoid pathways; and abiotic and biotic stress responses were identified. H4ac and H3K27me3 enrichment also differentially affects JAs and HDA6-mediated genome integrity and gene regulatory networks, substantiating the role of HDA6 interacting with specific families of transposable elements in planta and highlighting further specificity of action as well as novel targets of HDA6 in the context of JA signalling for abiotic and biotic stress responses.

CONCLUSIONS

The findings demonstrate functional overlap for MeJA and HDA6 in tuning plant developmental plasticity and response to stress at the histone modification level. MeJA and HDA6, nonetheless, maintain distinct activities on histone modifications to modulate genetic variability and to allow adaptation to environmental challenges.

Gibberellins regulate ovule number through a DELLA-CUC2 complex in Arabidopsis

Ovule development is a key process for plant reproduction, helping to ensure correct seed production. Several molecular factors and plant hormones such as gibberellins are involved in ovule initiation and development. Gibberellins control ovule development by the destabilization of DELLA proteins, whereas DELLA activity has been shown to act as a positive factor for ovule primordia emergence. But the molecular mechanism by which DELLA acts in ovule primordia initiation remained unknown. In this study we report that DELLA proteins participate in ovule initiation by the formation of a protein complex with the CUC2 transcription factor. The DELLA protein GAI requires CUC2 to promote ovule primordia formation, through the direct GAI-CUC2 interaction in placental cells that would determine the boundary regions between ovules during pistil development. Analysis of GAI-CUC2 interaction and co-localization in the placenta supports this hypothesis. Moreover, molecular analysis identified a subset of the loci for which the GAI protein may act as a transcriptional co-regulator in a CUC2-dependent manner. The DELLA-CUC2 complex is a component of the gene regulatory network controlling ovule primordia initiation in Arabidopsis.

Structural and genome-wide analyses suggest that transposon-derived protein SETMAR alters transcription and splicing

Extensive portions of the human genome have unknown function, including those derived from transposable elements. One such element, the DNA transposon Hsmar1, entered the primate lineage approximately 50 million years ago leaving behind terminal inverted repeat (TIR) sequences and a single intact copy of the Hsmar1 transposase, which retains its ancestral TIR-DNA-binding activity, and is fused with a lysine methyltransferase SET domain to constitute the chimeric SETMAR gene. Here, we provide a structural basis for recognition of TIRs by SETMAR and investigate the function of SETMAR through genome-wide approaches. As elucidated in our 2.37 Å crystal structure, SETMAR forms a dimeric complex with each DNA-binding domain bound specifically to TIR-DNA through the formation of 32 hydrogen bonds. We found that SETMAR recognizes primarily TIR sequences (∼5000 sites) within the human genome as assessed by chromatin immunoprecipitation sequencing analysis. In two SETMAR KO cell lines, we identified 163 shared differentially expressed genes and 233 shared alternative splicing events. Among these genes are several pre–mRNA-splicing factors, transcription factors, and genes associated with neuronal function, and one alternatively spliced primate-specific gene, TMEM14B, which has been identified as a marker for neocortex expansion associated with brain evolution. Taken together, our results suggest a model in which SETMAR impacts differential expression and alternative splicing of genes associated with transcription and neuronal function, potentially through both its TIR-specific DNA-binding and lysine methyltransferase activities, consistent with a role for SETMAR in simian primate development.

BRN2 as a key gene drives the early primate telencephalon development

Evolutionary mutations in primate-specific genes drove primate cortex expansion. However, whether conserved genes with previously unidentified functions also play a key role in primate brain expansion remains unknown. Here, we focus on BRN2 (POU3F2), a gene encoding a neural transcription factor commonly expressed in both primates and mice. Compared to the limited effects on mouse brain development, BRN2 biallelic knockout in cynomolgus monkeys (Macaca fascicularis) is lethal before midgestation. Histology analysis and single-cell transcriptome show that BRN2 deficiency decreases RGC expansion, induces precocious differentiation, and alters the trajectory of neurogenesis in the telencephalon. BRN2, serving as an upstream factor, controls specification and differentiation of ganglionic eminences. In addition, we identified the conserved function of BRN2 in cynomolgus monkeys to human RGCs. BRN2 may function by directly regulating SOX2 and STAT3 and maintaining HOPX. Our findings reveal a previously unknown mechanism that BRN2, a conserved gene, drives early primate telencephalon development by gaining novel mechanistic functions.

DNA-Protein Cross-Linking Sequencing for Genome-Wide Mapping of Thymidine Glycol

Thymidine glycol (Tg) is the most prevalent form of oxidatively induced pyrimidine lesions in DNA. Tg can arise from direct oxidation of thymidine in DNA. In addition, 5-methyl-2'-deoxycytidine (5-mdC) can be oxidized to 5-mdC glycol, and its subsequent deamination also yields Tg. However, Tg's distribution in the human genome remains unknown. Here, we presented a DNA-protein cross-linking sequencing (DPC-Seq) method for genome-wide mapping of Tg in human cells. Our approach capitalizes on the specificity of a bifunctional DNA glycosylase, i.e., NTHL1, for the covalent labeling, as well as DPC pulldown, SDS-PAGE fractionation, and membrane transfer for highly efficient and selective enrichment of Tg-bearing DNA. By employing DPC-Seq, we detected thousands of Tg sites in the human genome, where dual ablation of NTHL1 and NEIL1, the major DNA glycosylases responsible for Tg repair, led to pronounced increases in the number of Tg peaks. In addition, Tg is depleted in genomic regions associated with active transcription but enriched at nucleosome-binding sites, especially at heterochromatin sites marked with H3K9me2. Collectively, we developed a DPC-Seq method for highly efficient enrichment of Tg-containing DNA and for genome-wide mapping of Tg in human cells. Our work offers a robust tool for future functional studies of Tg in DNA, and we envision that the method can also be adapted for mapping other modified nucleosides in genomic DNA in the future.

Characterization of Chromatin Accessibility and Gene Expression upon Cold Stress Reveals that the RAV1 Transcription Factor Functions in Cold Response in Vitis Amurensis

Cold tolerance is regulated by a variety of transcription factors (TFs) and their target genes. Except for the well-characterized C-repeat binding factors (CBFs)-dependent transcriptional cascade, the mechanisms of cold tolerance mediated by other transcriptional regulatory networks are still largely unknown. Here, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq to identify cold responsive TFs in Vitis amurensis, a grape species with high cold hardiness. Nine TFs, including CBF4, RAV1 and ERF104, were identified after cold treatment. Weighted gene co-expression network analysis (WGCNA) and gene ontology (GO) analysis revealed that these TFs may regulate cold response through different pathways. As a prime candidate TF, overexpression of VaRAV1 in grape cells improved its cold tolerance. The transgenic cells exhibited low electrolyte leakage and malondialdehyde content and high peroxidase activity. Moreover, the TF gene TCP8 and a gene involving in homogalacturonan biosynthesis were found to be regulated by VaRAV1, suggesting that the contribution of VaRAV1 to cold tolerance may be achieved by enhancing the stability of cell membrane and regulating the expression of target genes involved in plant cell wall composition. Our work provides novel insights into plant response to cold stress and demonstrates the utility of ATAC-seq and RNA-seq for the rapid identification of TFs in response to cold stress in grapevine. VaRAV1 may play an important role in adaption to cold stress.

Discrete cis-acting element regulates developmentally timed gene-lamina relocation and neural progenitor competence in vivo

The nuclear lamina is typically associated with transcriptional silencing, and peripheral relocation of genes highly correlates with repression. However, the DNA sequences and proteins regulating gene-lamina interactions are largely unknown. Exploiting the developmentally timed hunchback gene movement to the lamina in Drosophila neuroblasts, we identified a 250 bp intronic element (IE) both necessary and sufficient for relocation. The IE can target a reporter transgene to the lamina and silence it. Endogenously, however, hunchback is already repressed prior to relocation. Instead, IE-mediated relocation confers a heritably silenced gene state refractory to activation in descendent neurons, which terminates neuroblast competence to specify early-born identity. Surprisingly, we found that the Polycomb group chromatin factors bind the IE and are required for lamina relocation, revealing a nuclear architectural role distinct from their well-known function in transcriptional repression. Together, our results uncover in vivo mechanisms underlying neuroblast competence and lamina association in heritable gene silencing.

In Drosophila neuroblasts, relocation of the hunchback gene locus to the nuclear lamina confers heritable silencing in daughter neurons. Lucas et al. identify a genomic element necessary and sufficient for hunchback gene movement in vivo. Polycomb proteins target this element for lamina relocation, thereby regulating competence, but not hunchback expression.

Neonatal Exposure to BPA, BDE-99, and PCB Produces Persistent Changes in Hepatic Transcriptome Associated With Gut Dysbiosis in Adult Mouse Livers

Recent evidence suggests that complex diseases can result from early life exposure to environmental toxicants. Polybrominated diphenyl ethers (PBDEs), and polychlorinated biphenyls (PCBs) are persistent organic pollutants (POPs) and remain a continuing risk to human health despite being banned from production. Developmental BPA exposure mediated-adult onset of liver cancer via epigenetic reprogramming mechanisms has been identified. Here, we investigated whether the gut microbiome and liver can be persistently reprogrammed following neonatal exposure to POPs, and the associations between microbial biomarkers and disease-prone changes in the hepatic transcriptome in adulthood, compared with BPA. C57BL/6 male and female mouse pups were orally administered vehicle, BPA, BDE-99 (a breast milk-enriched PBDE congener), or the Fox River PCB mixture (PCBs), once daily for three consecutive days (postnatal days [PND] 2-4). Tissues were collected at PND5 and PND60. Among the three chemicals investigated, early life exposure to BDE-99 produced the most prominent developmental reprogramming of the gut-liver axis, including hepatic inflammatory and cancer-prone signatures. In adulthood, neonatal BDE-99 exposure resulted in a persistent increase in Akkermansia muciniphila throughout the intestine, accompanied by increased hepatic levels of acetate and succinate, the known products of A. muciniphila. In males, this was positively associated with permissive epigenetic marks H3K4me1 and H3K27, which were enriched in loci near liver cancer-related genes that were dysregulated following neonatal exposure to BDE-99. Our findings provide novel insights that early life exposure to POPs can have a life-long impact on disease risk, which may partly be regulated by the gut microbiome.

AnnoMiner is a new web-tool to integrate epigenetics, transcription factor occupancy and transcriptomics data to predict transcriptional regulators

Gene expression regulation requires precise transcriptional programs, led by transcription factors in combination with epigenetic events. Recent advances in epigenomic and transcriptomic techniques provided insight into different gene regulation mechanisms. However, to date it remains challenging to understand how combinations of transcription factors together with epigenetic events control cell-type specific gene expression. We have developed the AnnoMiner web-server, an innovative and flexible tool to annotate and integrate epigenetic, and transcription factor occupancy data. First, AnnoMiner annotates user-provided peaks with gene features. Second, AnnoMiner can integrate genome binding data from two different transcriptional regulators together with gene features. Third, AnnoMiner offers to explore the transcriptional deregulation of genes nearby, or within a specified genomic region surrounding a user-provided peak. AnnoMiner’s fourth function performs transcription factor or histone modification enrichment analysis for user-provided gene lists by utilizing hundreds of public, high-quality datasets from ENCODE for the model organisms human, mouse, Drosophila and C. elegans. Thus, AnnoMiner can predict transcriptional regulators for a studied process without the strict need for chromatin data from the same process. We compared AnnoMiner to existing tools and experimentally validated several transcriptional regulators predicted by AnnoMiner to indeed contribute to muscle morphogenesis in Drosophila. AnnoMiner is freely available at http://chimborazo.ibdm.univ-mrs.fr/AnnoMiner/.

Reactivation of the Hedgehog pathway in esophageal progenitors turns on an embryonic-like program to initiate columnar metaplasia

Columnar metaplasia of the esophagus is the main risk factor for esophageal adenocarcinoma. There is a lack of evidence to demonstrate that esophageal progenitors can be the source of columnar metaplasia. In this study, using transgenic mouse models, lineage tracing, single-cell RNA sequencing, and transcriptomic and epigenetic profiling, we found that the activation of the Hedgehog pathway in esophageal cells modifies their differentiation status in vivo. This process involves an initial step of dedifferentiation into embryonic-like esophageal progenitors. Moreover, a subset of these cells undergoes full squamous-to-columnar conversion and expresses selected intestinal markers. These modifications of cell fate are associated with remodeling of the chromatin and the appearance of Sox9. Using a conditional knockout mouse, we show that Sox9 is required for columnar conversion but not for the step of dedifferentiation. These results provide insight into the mechanisms by which esophageal cells might initiate columnar metaplasia.

Cooperation of LIM domain-binding 2 (LDB2) with EGR in the pathogenesis of schizophrenia

Genomic defects with large effect size can help elucidate unknown pathologic architecture of mental disorders. We previously reported on a patient with schizophrenia and a balanced translocation between chromosomes 4 and 13 and found that the breakpoint within chromosome 4 is located near the LDB2 gene. We show here that Ldb2 knockout (KO) mice displayed multiple deficits relevant to mental disorders. In particular, Ldb2 KO mice exhibited deficits in the fear-conditioning paradigm. Analysis of the amygdala suggested that dysregulation of synaptic activities controlled by the immediate early gene Arc is involved in the phenotypes. We show that LDB2 forms protein complexes with known transcription factors. Consistently, ChIP-seq analyses indicated that LDB2 binds to > 10,000 genomic sites in human neurospheres. We found that many of those sites, including the promoter region of ARC, are occupied by EGR transcription factors. Our previous study showed an association of the EGR family genes with schizophrenia. Collectively, the findings suggest that dysregulation in the gene expression controlled by the LDB2-EGR axis underlies a pathogenesis of subset of mental disorders.

Molecular mechanism of cytokinin-activated cell division in Arabidopsis

Mitogens trigger cell division in animals. In plants, cytokinins, a group of phytohormones derived from adenine, stimulate cell proliferation. Cytokinin signaling is initiated by membrane-associated histidine kinase receptors and transduced through a phosphorelay system. We show that in the Arabidopsis shoot apical meristem (SAM), cytokinin regulates cell division by promoting nuclear shuttling of Myb-domain protein 3R4 (MYB3R4), a transcription factor that activates mitotic gene expression. Newly synthesized MYB3R4 protein resides predominantly in the cytoplasm. At the G2-to-M transition, rapid nuclear accumulation of MYB3R4-consistent with an associated transient peak in cytokinin concentration-feeds a positive feedback loop involving importins and initiates a transcriptional cascade that drives mitosis and cytokinesis. An engineered nuclear-restricted MYB3R4 mimics the cytokinin effects of enhanced cell proliferation and meristem growth.

OsFD4 promotes the rice floral transition via florigen activation complex formation in the shoot apical meristem

In rice, the florigens Heading Date 3a (Hd3a) and Rice Flowering Locus T 1 (RFT1), OsFD-like basic leucine zipper (bZIP) transcription factors, and Gf14 proteins assemble into florigen activation/repressor complexes (FACs/FRCs), which regulate transition to flowering in leaves and apical meristem. Only OsFD1 has been described as part of complexes promoting flowering at the meristem, and little is known about the role of other bZIP transcription factors, the combinatorial complexity of FAC formation, and their DNA-binding properties. Here, we used mutant analysis, protein-protein interaction assays and DNA affinity purification (DAP) sequencing coupled to in silico prediction of binding syntaxes to study several bZIP proteins that assemble into FACs or FRCs. We identified OsFD4 as a component of a FAC promoting flowering at the shoot apical meristem, downstream of OsFD1. The osfd4 mutants are late flowering and delay expression of genes promoting inflorescence development. Protein-protein interactions indicate an extensive network of contacts between several bZIPs and Gf14 proteins. Finally, we identified genomic regions bound by bZIPs with promotive and repressive effects on flowering. We conclude that distinct bZIPs orchestrate floral induction at the meristem and that FAC formation is largely combinatorial. While binding to the same consensus motif, their DNA-binding syntax is different, suggesting discriminatory functions.

Dissecting Hes-centred transcriptional networks in neural stem cell maintenance and tumorigenesis in Drosophila

Neural stem cells divide during embryogenesis and juvenile life to generate the entire complement of neurons and glia in the nervous system of vertebrates and invertebrates. Studies of the mechanisms controlling the fine balance between neural stem cells and more differentiated progenitors have shown that, in every asymmetric cell division, progenitors send a Delta-Notch signal to their sibling stem cells. Here, we show that excessive activation of Notch or overexpression of its direct targets of the Hes family causes stem-cell hyperplasias in the Drosophila larval central nervous system, which can progress to malignant tumours after allografting to adult hosts. We combined transcriptomic data from these hyperplasias with chromatin occupancy data for Dpn, a Hes transcription factor, to identify genes regulated by Hes factors in this process. We show that the Notch/Hes axis represses a cohort of transcription factor genes. These are excluded from the stem cells and promote early differentiation steps, most likely by preventing the reversion of immature progenitors to a stem-cell fate. We describe the impact of two of these 'anti-stemness' factors, Zfh1 and Gcm, on Notch/Hes-triggered tumorigenesis.

Loss of histone H3.3 results in DNA replication defects and altered origin dynamics in C. elegans

Histone H3.3 is a replication-independent variant of histone H3 with important roles in development, differentiation, and fertility. Here, we show that loss of H3.3 results in replication defects in Caenorhabditis elegans embryos at elevated temperatures. To characterize these defects, we adapt methods to determine replication timing, map replication origins, and examine replication fork progression. Our analysis of the spatiotemporal regulation of DNA replication shows that despite the very rapid embryonic cell cycle, the genome is replicated from early and late firing origins and is partitioned into domains of early and late replication. We find that under temperature stress conditions, additional replication origins become activated. Moreover, loss of H3.3 results in altered replication fork progression around origins, which is particularly evident at stress-activated origins. These replication defects are accompanied by replication checkpoint activation, a delayed cell cycle, and increased lethality in checkpoint-compromised embryos. Our comprehensive analysis of DNA replication in C. elegans reveals the genomic location of replication origins and the dynamics of their firing, and uncovers a role of H3.3 in the regulation of replication origins under stress conditions.

Disturbed mitochondrial dynamics in CD8+ TILs reinforce T cell exhaustion

The metabolic challenges present in tumors attenuate the metabolic fitness and antitumor activity of tumor-infiltrating T lymphocytes (TILs). However, it remains unclear whether persistent metabolic insufficiency can imprint permanent T cell dysfunction. We found that TILs accumulated depolarized mitochondria as a result of decreased mitophagy activity and displayed functional, transcriptomic and epigenetic characteristics of terminally exhausted T cells. Mechanistically, reduced mitochondrial fitness in TILs was induced by the coordination of T cell receptor stimulation, microenvironmental stressors and PD-1 signaling. Enforced accumulation of depolarized mitochondria with pharmacological inhibitors induced epigenetic reprogramming toward terminal exhaustion, indicating that mitochondrial deregulation caused T cell exhaustion. Furthermore, supplementation with nicotinamide riboside enhanced T cell mitochondrial fitness and improved responsiveness to anti-PD-1 treatment. Together, our results reveal insights into how mitochondrial dynamics and quality orchestrate T cell antitumor responses and commitment to the exhaustion program.

Diet Modifies Pioglitazone's Influence on Hepatic PPARγ-Regulated Mitochondrial Gene Expression

Pioglitazone (Pio) is a thiazolidinedione (TZD) insulin-sensitizing drug whose effects result predominantly from its modulation of the transcriptional activity of peroxisome proliferator-activated-receptor-gamma (PPARγ). Pio is used to treat human insulin-resistant diabetes and also frequently considered for treatment of nonalcoholic steatohepatitis (NASH). In both settings, Pio's beneficial effects are believed to result primarily from its actions on adipose PPARγ activity, which improves insulin sensitivity and reduces the delivery of fatty acids to the liver. Nevertheless, a recent clinical trial showed variable efficacy of Pio in human NASH. Hepatocytes also express PPARγ, and such expression increases with insulin resistance and in nonalcoholic fatty liver disease (NAFLD). Furthermore, mice that overexpress hepatocellular PPARγ and Pio-treated mice with extrahepatic PPARγ gene disruption develop features of NAFLD. Thus, Pio's direct impact on hepatocellular gene expression might also be a determinant of this drug's ultimate influence on insulin resistance and NAFLD. Previous studies have characterized Pio's PPARγ-dependent effects on hepatic expression of specific adipogenic, lipogenic, and other metabolic genes. However, such transcriptional regulation has not been comprehensively assessed. The studies reported here address that consideration by genome-wide comparisons of Pio's hepatic transcriptional effects in wildtype (WT) and liver-specific PPARγ-knockout (KO) mice given either control or high-fat (HFD) diets. The results identify a large set of hepatic genes for which Pio's liver PPARγ-dependent transcriptional effects are concordant with its effects on RXR-DNA binding in WT mice. These data also show that HFD modifies Pio's influence on a subset of such transcriptional regulation. Finally, our findings reveal a broader influence of Pio on PPARγ-dependent hepatic expression of nuclear genes encoding mitochondrial proteins than previously recognized. Taken together, these studies provide new insights about the tissue-specific mechanisms by which Pio affects hepatic gene expression and the broad scope of this drug's influence on such regulation.

Landscape of G-quadruplex DNA structural regions in breast cancer

Response and resistance to anticancer therapies vary due to intertumor and intratumor heterogeneity1. Here, we map differentially enriched G-quadruplex (G4) DNA structure-forming regions (∆G4Rs) in 22 breast cancer patient-derived tumor xenograft (PDTX) models. ∆G4Rs are associated with the promoters of highly amplified genes showing high expression, and with somatic single-nucleotide variants. Differences in ΔG4R landscapes reveal seven transcription factor programs across PDTXs. ∆G4R abundance and locations stratify PDTXs into at least three G4-based subtypes. ∆G4Rs in most PDTXs (14 of 22) were found to associate with more than one breast cancer subtype, which we also call an integrative cluster (IC)2. This suggests the frequent coexistence of multiple breast cancer states within a PDTX model, the majority of which display aggressive triple-negative IC10 gene activity. Short-term cultures of PDTX models with increased ∆G4R levels are more sensitive to small molecules targeting G4 DNA. Thus, G4 landscapes reveal additional IC-related intratumor heterogeneity in PDTX biopsies, improving breast cancer stratification and potentially identifying new treatment strategies.

Acetylation-dependent regulation of PD-L1 nuclear translocation dictates the efficacy of anti-PD-1 immunotherapy

Immunotherapies that target programmed cell death protein 1 (PD-1) and its ligand PD-L1 as well as cytotoxic T-lymphocyte-associated protein 4 (CTLA4) have shown impressive clinical outcomes for multiple tumours. However, only a subset of patients achieves durable responses, suggesting that the mechanisms of the immune checkpoint pathways are not completely understood. Here, we report that PD-L1 translocates from the plasma membrane into the nucleus through interactions with components of the endocytosis and nucleocytoplasmic transport pathways, regulated by p300-mediated acetylation and HDAC2-dependent deacetylation of PD-L1. Moreover, PD-L1 deficiency leads to compromised expression of multiple immune-response-related genes. Genetically or pharmacologically modulating PD-L1 acetylation blocks its nuclear translocation, reprograms the expression of immune-response-related genes and, as a consequence, enhances the anti-tumour response to PD-1 blockade. Thus, our results reveal an acetylation-dependent regulation of PD-L1 nuclear localization that governs immune-response gene expression, and thereby advocate targeting PD-L1 translocation to enhance the efficacy of PD-1/PD-L1 blockade.

Induction of ICER is superseded by smICER, challenging the impact of ICER under chronic beta-adrenergic stimulation

ICER (inducible cAMP early repressor) isoforms are transcriptional repressors encoded by the Crem (cAMP responsive element modulator) gene. They were linked to the regulation of a multitude of cellular processes and pathophysiological mechanisms. Here, we show for the first time that two independent induction patterns for CREM repressor isoforms exist in the heart, namely for ICER and smICER (small ICER), which are induced in response to β-adrenergic stimulation in a transient- and saturation-like manner, respectively. This time-shifted induction pattern, driven by two internal promoters in the Crem gene, leads to the predominant transcription of smIcer after prolonged β-adrenergic stimulation. Using an ICER knockout mouse model with preserved smICER induction, we show that the transient-like induction of Icer itself has minor effects on gene regulation, cardiac hypertrophy or contractile function in the heart. We conclude that the functions previously linked to ICER may be rather attributed to smICER, also beyond the cardiac background.

RNA-DNA strand exchange by the Drosophila Polycomb complex PRC2

Polycomb Group (PcG) proteins form memory of transient transcriptional repression that is necessary for development. In Drosophila, DNA elements termed Polycomb Response Elements (PREs) recruit PcG proteins. How PcG activities are targeted to PREs to maintain repressed states only in appropriate developmental contexts has been difficult to elucidate. PcG complexes modify chromatin, but also interact with both RNA and DNA, and RNA is implicated in PcG targeting and function. Here we show that R-loops form at many PREs in Drosophila embryos, and correlate with repressive states. In vitro, both PRC1 and PRC2 can recognize R-loops and open DNA bubbles. Unexpectedly, we find that PRC2 drives formation of RNA-DNA hybrids, the key component of R-loops, from RNA and dsDNA. Our results identify R-loop formation as a feature of Drosophila PREs that can be recognized by PcG complexes, and RNA-DNA strand exchange as a PRC2 activity that could contribute to R-loop formation.

Computational Oncology in the Multi-Omics Era: State of the Art

Cancer is the quintessential complex disease. As technologies evolve faster each day, we are able to quantify the different layers of biological elements that contribute to the emergence and development of malignancies. In this multi-omics context, the use of integrative approaches is mandatory in order to gain further insights on oncological phenomena, and to move forward toward the precision medicine paradigm. In this review, we will focus on computational oncology as an integrative discipline that incorporates knowledge from the mathematical, physical, and computational fields to further the biomedical understanding of cancer. We will discuss the current roles of computation in oncology in the context of multi-omic technologies, which include: data acquisition and processing; data management in the clinical and research settings; classification, diagnosis, and prognosis; and the development of models in the research setting, including their use for therapeutic target identification. We will discuss the machine learning and network approaches as two of the most promising emerging paradigms, in computational oncology. These approaches provide a foundation on how to integrate different layers of biological description into coherent frameworks that allow advances both in the basic and clinical settings.

Human Endometrial Transcriptome and Progesterone Receptor Cistrome Reveal Important Pathways and Epithelial Regulators

CONTEXT

Poor uterine receptivity is one major factor leading to pregnancy loss and infertility. Understanding the molecular events governing successful implantation is hence critical in combating infertility.

OBJECTIVE

To define Progesterone Receptor (PGR)-regulated molecular mechanisms and epithelial roles in receptivity.

DESIGN

RNA-sequencing and PGR-ChIP-seq were conducted in parallel to identify PGR-regulated pathways during the Window of implantation (WOI) in endometrium of fertile women.

SETTING

Endometrial biopsies from the proliferative and mid-secretory phases were analyzed.

PATIENTS OR OTHER PARTICIPANTS

Participants were fertile, reproductive aged (18-37 years) women with normal cycle length, and without any history of dysmenorrhea, infertility, or irregular cycles. In total, 42 endometrial biopsies obtained from 42 women were analyzed in this study.

INTERVENTIONS

There were no interventions during this study.

MAIN OUTCOME MEASURES

Here we measured the alterations in gene expression and PGR occupancy in the genome during the WOI, based on the hypothesis that PGR binds uterine chromatin cycle dependently to regulate genes involved in uterine cell differentiation and function.

RESULTS

653 genes were identified with regulated PGR binding and differential expression during the WOI. These were involved in regulating inflammatory response, xenobiotic metabolism, epithelial mesenchymal transition, cell death, interleukin/Signal Transducer And Activator Of Transcription (STAT) signaling, estrogen response, and Mammalian target of rapamycin complex 1 (MTORC1) response. Transcriptome of the epithelium identified 3052 differentially expressed genes, of which 658 were uniquely regulated. Transcription factors Interferon Regulatory Factor 8 (IRF8) and Myocyte Enhancer Factor 2C (MEF2C) were found to be regulated in the epithelium during the WOI at the protein level, suggesting potentially important functions that are previously unrecognized.

CONCLUSION

PGR binds the genomic regions of genes regulating critical processes in uterine receptivity and function.

BedSect: An Integrated Web Server Application to Perform Intersection, Visualization, and Functional Annotation of Genomic Regions From Multiple Datasets

A large number of genomic regions, such as transcription factor binding sites (TFBSs) captured from next generation sequencing (NGS) data analyses or those available from the public resource database ENCODE, are generally overlapped to answer a variety of biological questions. Though several command-line tools are available to perform such an analysis, there is a notable lack of an integrated webserver application with which to identify genomic region intersections, generate publication-ready plots depicting subsets of the overlapped regions, and perform functional annotation. Thus, there is an ardent need for a comprehensive and user-friendly webserver application that allows the users to either upload multiple datasets or select from the integrated Gene Transcription Regulation Database (GTRD). We thus introduce BedSect (http://imgsb.org/bedsect/.), which not only fulfils the above criteria but also performs intersection analysis along with visualization of the intersection regions as an UpSet and correlation plot using the integrated Shiny application. Moreover, analyses, including functional annotation, gene ontology, and biological pathways enrichment for the identified unique and intersected genomic regions, can also be performed using the integrated GREAT tool. To view the genomic regions in the genome browser, the inbuilt hyperlink for UCSC can redirect the user to visualize the results as custom tracks.

5-hydroxymethylcytosine and gene activity in mouse intestinal differentiation

Cytosine hydroxymethylation (5hmC) in mammalian DNA is the product of oxidation of methylated cytosines (5mC) by Ten-Eleven-Translocation (TET) enzymes. While it has been shown that the TETs influence 5mC metabolism, pluripotency and differentiation during early embryonic development, the functional relationship between gene expression and 5hmC in adult (somatic) stem cell differentiation is still unknown. Here we report that 5hmC levels undergo highly dynamic changes during adult stem cell differentiation from intestinal progenitors to differentiated intestinal epithelium. We profiled 5hmC and gene activity in purified mouse intestinal progenitors and differentiated progeny to identify 43425 differentially hydroxymethylated regions and 5325 differentially expressed genes. These differentially marked regions showed both losses and gains of 5hmC after differentiation, despite lower global levels of 5hmC in progenitor cells. In progenitors, 5hmC did not correlate with gene transcript levels, however, upon differentiation the global increase in 5hmC content showed an overall positive correlation with gene expression level as well as prominent associations with histone modifications that typify active genes and enhancer elements. Our data support a gene regulatory role for 5hmC that is predominant over its role in controlling DNA methylation states.

Cell lineage inference from SNP and scRNA-Seq data

Several recent studies focus on the inference of developmental and response trajectories from single cell RNA-Seq (scRNA-Seq) data. A number of computational methods, often referred to as pseudo-time ordering, have been developed for this task. Recently, CRISPR has also been used to reconstruct lineage trees by inserting random mutations. However, both approaches suffer from drawbacks that limit their use. Here, we develop a method to detect significant, cell type specific, sequence mutations from scRNA-Seq data. We show that only a few mutations are enough for reconstructing good branching models. Integrating these mutations with expression data further improves the accuracy of the reconstructed models. As we show, the majority of mutations we identify are likely RNA editing events indicating that such information can be used to distinguish cell types.

WUSCHEL acts as an auxin response rheostat to maintain apical stem cells in Arabidopsis

To maintain the balance between long-term stem cell self-renewal and differentiation, dynamic signals need to be translated into spatially precise and temporally stable gene expression states. In the apical plant stem cell system, local accumulation of the small, highly mobile phytohormone auxin triggers differentiation while at the same time, pluripotent stem cells are maintained throughout the entire life-cycle. We find that stem cells are resistant to auxin mediated differentiation, but require low levels of signaling for their maintenance. We demonstrate that the WUSCHEL transcription factor confers this behavior by rheostatically controlling the auxin signaling and response pathway. Finally, we show that WUSCHEL acts via regulation of histone acetylation at target loci, including those with functions in the auxin pathway. Our results reveal an important mechanism that allows cells to differentially translate a potent and highly dynamic developmental signal into stable cell behavior with high spatial precision and temporal robustness.

Spatial control of auxin signaling maintains a balance between stem-cell self-renewal and differentiation at the plant shoot apex. Here Ma et al. show that rheostatic control of auxin response by the WUSCHEL transcription factor maintains stem cells by conferring resistance to auxin mediated differentiation.

Gene regulatory network inference resources: A practical overview

Transcriptional regulation is a fundamental molecular mechanism involved in almost every aspect of life, from homeostasis to development, from metabolism to behavior, from reaction to stimuli to disease progression. In recent years, the concept of Gene Regulatory Networks (GRNs) has grown popular as an effective applied biology approach for describing the complex and highly dynamic set of transcriptional interactions, due to its easy-to-interpret features. Since cataloguing, predicting and understanding every GRN connection in all species and cellular contexts remains a great challenge for biology, researchers have developed numerous tools and methods to infer regulatory processes. In this review, we catalogue these methods in six major areas, based on the dominant underlying information leveraged to infer GRNs: Coexpression, Sequence Motifs, Chromatin Immunoprecipitation (ChIP), Orthology, Literature and Protein-Protein Interaction (PPI) specifically focused on transcriptional complexes. The methods described here cover a wide range of user-friendliness: from web tools that require no prior computational expertise to command line programs and algorithms for large scale GRN inferences. Each method for GRN inference described herein effectively illustrates a type of transcriptional relationship, with many methods being complementary to others. While a truly holistic approach for inferring and displaying GRNs remains one of the greatest challenges in the field of systems biology, we believe that the integration of multiple methods described herein provides an effective means with which experimental and computational biologists alike may obtain the most complete pictures of transcriptional relationships. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Recapitulation and Reversal of Schizophrenia-Related Phenotypes in Setd1a-Deficient Mice

SETD1A, a lysine-methyltransferase, is a key schizophrenia susceptibility gene. Mice carrying a heterozygous loss-of-function mutation of the orthologous gene exhibit alterations in axonal branching and cortical synaptic dynamics accompanied by working memory deficits. We show that Setd1a binds both promoters and enhancers with a striking overlap between Setd1a and Mef2 on enhancers. Setd1a targets are highly expressed in pyramidal neurons and display a complex pattern of transcriptional up- and downregulations shaped by presumed opposing functions of Setd1a on promoters and Mef2-bound enhancers. Notably, evolutionarily conserved Setd1a targets are associated with neuropsychiatric genetic risk burden. Reinstating Setd1a expression in adulthood rescues cognitive deficits. Finally, we identify LSD1 as a major counteracting demethylase for Setd1a and show that its pharmacological antagonism results in a full rescue of the behavioral and morphological deficits in Setd1a-deficient mice. Our findings advance understanding of how SETD1A mutations predispose to schizophrenia (SCZ) and point to novel therapeutic interventions.

Genome-Wide Transcription Factor Binding in Leaves from C3 and C4 Grasses

The majority of plants use C3 photosynthesis, but over 60 independent lineages of angiosperms have evolved the C4 pathway. In most C4 species, photosynthesis gene expression is compartmented between mesophyll and bundle-sheath cells. We performed DNaseI sequencing to identify genome-wide profiles of transcription factor binding in leaves of the C4 grasses Zea mays, Sorghum bicolor, and Setaria italica as well as C3 Brachypodium distachyon In C4 species, while bundle-sheath strands and whole leaves shared similarity in the broad regions of DNA accessible to transcription factors, the short sequences bound varied. Transcription factor binding was prevalent in gene bodies as well as promoters, and many of these sites could represent duons that influence gene regulation in addition to amino acid sequence. Although globally there was little correlation between any individual DNaseI footprint and cell-specific gene expression, within individual species transcription factor binding to the same motifs in multiple genes provided evidence for shared mechanisms governing C4 photosynthesis gene expression. Furthermore, interspecific comparisons identified a small number of highly conserved transcription factor binding sites associated with leaves from species that diverged around 60 million years ago. These data therefore provide insight into the architecture associated with C4 photosynthesis gene expression in particular and characteristics of transcription factor binding in cereal crops in general.

Cooperation of cancer drivers with regulatory germline variants shapes clinical outcomes

Pediatric malignancies including Ewing sarcoma (EwS) feature a paucity of somatic alterations except for pathognomonic driver-mutations that cannot explain overt variations in clinical outcome. Here, we demonstrate in EwS how cooperation of dominant oncogenes and regulatory germline variants determine tumor growth, patient survival and drug response. Binding of the oncogenic EWSR1-FLI1 fusion transcription factor to a polymorphic enhancer-like DNA element controls expression of the transcription factor MYBL2 mediating these phenotypes. Whole-genome and RNA sequencing reveals that variability at this locus is inherited via the germline and is associated with variable inter-tumoral MYBL2 expression. High MYBL2 levels sensitize EwS cells for inhibition of its upstream activating kinase CDK2 in vitro and in vivo, suggesting MYBL2 as a putative biomarker for anti-CDK2-therapy. Collectively, we establish cooperation of somatic mutations and regulatory germline variants as a major determinant of tumor progression and highlight the importance of integrating the regulatory genome in precision medicine.