The Landscape of L1 Retrotransposons in the Human Genome Is Shaped by Pre-insertion Sequence Biases and Post-insertion Selection.

L1 retrotransposons are transposable elements and major contributors of genetic variation in humans. Where L1 integrates into the genome can directly impact human evolution and disease. Here, we experimentally induced L1 retrotransposition in cells and mapped integration sites at nucleotide resolution. At local scales, L1 integration is mostly restricted by genome sequence biases and the specificity of the L1 machinery. At regional scales, L1 shows a broad capacity for integration into all chromatin states, in contrast to other known mobile genetic elements. However, integration is influenced by the replication timing of target regions, suggesting a link to host DNA replication. The distribution of new L1 integrations differs from those of preexisting L1 copies, which are significantly reshaped by natural selection. Our findings reveal that the L1 machinery has evolved to efficiently target all genomic regions and underline a predominant role for post-integrative processes on the distribution of endogenous L1 elements.

Tyrosine-1 of RNA Polymerase II CTD Controls Global Termination of Gene Transcription in Mammals.

The carboxy-terminal domain (CTD) of RNA polymerase (Pol) II is composed of a repetition of YSPTSPS heptads and functions as a loading platform for protein complexes that regulate transcription, splicing, and maturation of RNAs. Here, we studied mammalian CTD mutants to analyze the function of tyrosine1 residues in the transcription cycle. Mutation of 3/4 of the tyrosine residues (YFFF mutant) resulted in a massive read-through transcription phenotype in the antisense direction of promoters as well as in the 3' direction several hundred kilobases downstream of genes. The YFFF mutant shows reduced Pol II at promoter-proximal pause sites, a loss of interaction with the Mediator and Integrator complexes, and impaired recruitment of these complexes to chromatin. Consistent with these observations, Pol II loading at enhancers and maturation of snRNAs are altered in the YFFF context genome-wide. We conclude that tyrosine1 residues of the CTD control termination of transcription by Pol II.

RNA polymerase II CTD is dispensable for transcription and required for termination in human cells.

The largest subunit of RNA polymerase (Pol) II harbors an evolutionarily conserved C-terminal domain (CTD), composed of heptapeptide repeats, central to the transcriptional process. Here, we analyze the transcriptional phenotypes of a CTD-Δ5 mutant that carries a large CTD truncation in human cells. Our data show that this mutant can transcribe genes in living cells but displays a pervasive phenotype with impaired termination, similar to but more severe than previously characterized mutations of CTD tyrosine residues. The CTD-Δ5 mutant does not interact with the Mediator and Integrator complexes involved in the activation of transcription and processing of RNAs. Examination of long-distance interactions and CTCF-binding patterns in CTD-Δ5 mutant cells reveals no changes in TAD domains or borders. Our data demonstrate that the CTD is largely dispensable for the act of transcription in living cells. We propose a model in which CTD-depleted Pol II has a lower entry rate onto DNA but becomes pervasive once engaged in transcription, resulting in a defect in termination.

Notch-dependent and -independent functions of transcription factor RBPJ.

Signal transduction pathways often involve transcription factors that promote activation of defined target gene sets. The transcription factor RBPJ is the central player in Notch signaling and either forms an activator complex with the Notch intracellular domain (NICD) or a repressor complex with corepressors like KYOT2/FHL1. The balance between these two antagonizing RBPJ-complexes depends on the activation state of the Notch receptor regulated by cell-to-cell interaction, ligand binding and proteolytic cleavage events. Here, we depleted RBPJ in mature T-cells lacking active Notch signaling and performed RNA-Seq, ChIP-Seq and ATAC-seq analyses. RBPJ depletion leads to upregulation of many Notch target genes. Ectopic expression of NICD1 activates several Notch target genes and enhances RBPJ occupancy. Based on gene expression changes and RBPJ occupancy we define four different clusters, either RBPJ- and/or Notch-regulated genes. Importantly, we identify early (Hes1 and Hey1) and late Notch-responsive genes (IL2ra). Similarly, to RBPJ depletion, interfering with transcriptional repression by squelching with cofactor KYOT2/FHL1, leads to upregulation of Notch target genes. Taken together, RBPJ is not only an essential part of the Notch co-activator complex but also functions as a repressor in a Notch-independent manner.

Evaluation of DNA Methylation Episignatures for Diagnosis and Phenotype Correlations in 42 Mendelian Neurodevelopmental Disorders.

Genetic syndromes frequently present with overlapping clinical features and inconclusive or ambiguous genetic findings which can confound accurate diagnosis and clinical management. An expanding number of genetic syndromes have been shown to have unique genomic DNA methylation patterns (called "episignatures"). Peripheral blood episignatures can be used for diagnostic testing as well as for the interpretation of ambiguous genetic test results. We present here an approach to episignature mapping in 42 genetic syndromes, which has allowed the identification of 34 robust disease-specific episignatures. We examine emerging patterns of overlap, as well as similarities and hierarchical relationships across these episignatures, to highlight their key features as they are related to genetic heterogeneity, dosage effect, unaffected carrier status, and incomplete penetrance. We demonstrate the necessity of multiclass modeling for accurate genetic variant classification and show how disease classification using a single episignature at a time can sometimes lead to classification errors in closely related episignatures. We demonstrate the utility of this tool in resolving ambiguous clinical cases and identification of previously undiagnosed cases through mass screening of a large cohort of subjects with developmental delays and congenital anomalies. This study more than doubles the number of published syndromes with DNA methylation episignatures and, most significantly, opens new avenues for accurate diagnosis and clinical assessment in individuals affected by these disorders.

Fra-1 regulates its target genes via binding to remote enhancers without exerting major control on chromatin architecture in triple negative breast cancers.

The ubiquitous family of dimeric transcription factors AP-1 is made up of Fos and Jun family proteins. It has long been thought to operate principally at gene promoters and how it controls transcription is still ill-understood. The Fos family protein Fra-1 is overexpressed in triple negative breast cancers (TNBCs) where it contributes to tumor aggressiveness. To address its transcriptional actions in TNBCs, we combined transcriptomics, ChIP-seqs, machine learning and NG Capture-C. Additionally, we studied its Fos family kin Fra-2 also expressed in TNBCs, albeit much less. Consistently with their pleiotropic effects, Fra-1 and Fra-2 up- and downregulate individually, together or redundantly many genes associated with a wide range of biological processes. Target gene regulation is principally due to binding of Fra-1 and Fra-2 at regulatory elements located distantly from cognate promoters where Fra-1 modulates the recruitment of the transcriptional co-regulator p300/CBP and where differences in AP-1 variant motif recognition can underlie preferential Fra-1- or Fra-2 bindings. Our work also shows no major role for Fra-1 in chromatin architecture control at target gene loci, but suggests collaboration between Fra-1-bound and -unbound enhancers within chromatin hubs sometimes including promoters for other Fra-1-regulated genes. Our work impacts our view of AP-1.

Critical role for TRIM28 and HP1ß/γ in the epigenetic control of T cell metabolic reprograming and effector differentiation.

Naive CD4+ T lymphocytes differentiate into different effector types, including helper and regulatory cells (Th and Treg, respectively). Heritable gene expression programs that define these effector types are established during differentiation, but little is known about the epigenetic mechanisms that install and maintain these programs. Here, we use mice defective for different components of heterochromatin-dependent gene silencing to investigate the epigenetic control of CD4+ T cell plasticity. We show that, upon T cell receptor (TCR) engagement, naive and regulatory T cells defective for TRIM28 (an epigenetic adaptor for histone binding modules) or for heterochromatin protein 1 ß and γ isoforms (HP1ß/γ, 2 histone-binding factors involved in gene silencing) fail to effectively signal through the PI3K-AKT-mTOR axis and switch to glycolysis. While differentiation of naive TRIM28-/- T cells into cytokine-producing effector T cells is impaired, resulting in reduced induction of autoimmune colitis, TRIM28-/- regulatory T cells also fail to expand in vivo and to suppress autoimmunity effectively. Using a combination of transcriptome and chromatin immunoprecipitation-sequencing (ChIP-seq) analyses for H3K9me3, H3K9Ac, and RNA polymerase II, we show that reduced effector differentiation correlates with impaired transcriptional silencing at distal regulatory regions of a defined set of Treg-associated genes, including, for example, NRP1 or Snai3. We conclude that TRIM28 and HP1ß/γ control metabolic reprograming through epigenetic silencing of a defined set of Treg-characteristic genes, thus allowing effective T cell expansion and differentiation into helper and regulatory phenotypes.

MIR sequences recruit zinc finger protein ZNF768 to expressed genes.

Mammalian-wide interspersed repeats (MIRs) are retrotransposed elements of mammalian genomes. Here, we report the specific binding of zinc finger protein ZNF768 to the sequence motif GCTGTGTG (N20) CCTCTCTG in the core region of MIRs. ZNF768 binding is preferentially associated with euchromatin and promoter regions of genes. Binding was observed for genes expressed in a cell type-specific manner in human B cell line Raji and osteosarcoma U2OS cells. Mass spectrometric analysis revealed binding of ZNF768 to Elongator components Elp1, Elp2 and Elp3 and other nuclear factors. The N-terminus of ZNF768 contains a heptad repeat array structurally related to the C-terminal domain (CTD) of RNA polymerase II. This array evolved in placental animals but not marsupials and monotreme species, displays species-specific length variations, and possibly fulfills CTD related functions in gene regulation. We propose that the evolution of MIRs and ZNF768 has extended the repertoire of gene regulatory mechanisms in mammals and that ZNF768 binding is associated with cell type-specific gene expression.

Noncoding transcription at enhancers: general principles and functional models.

Mammalian genomes are extensively transcribed outside the borders of protein-coding genes. Genome-wide studies recently demonstrated that cis-regulatory genomic elements implicated in transcriptional control, such as enhancers and locus-control regions, represent major sites of extragenic noncoding transcription. Enhancer-templated transcripts provide a quantitatively small contribution to the total amount of cellular nonribosomal RNA; nevertheless, the possibility that enhancer transcription and the resulting enhancer RNAs may, in some cases, have functional roles, rather than represent mere transcriptional noise at accessible genomic regions, is supported by an increasing amount of experimental data. In this article we review the current knowledge on enhancer transcription and its functional implications.

Immunopathological manifestations in Kabuki syndrome: a registry study of 177 individuals.

PURPOSE: Kabuki syndrome (KS) (OMIM 147920 and 300867) is a rare genetic disorder characterized by specific facial features, intellectual disability, and various malformations. Immunopathological manifestations seem prevalent and increase the morbimortality. To assess the frequency and severity of the manifestations, we measured the prevalence of immunopathological manifestations as well as genotype-phenotype correlations in KS individuals from a registry. METHODS: Data were for 177 KS individuals with KDM6A or KMT2D pathogenic variants. Questionnaires to clinicians were used to assess the presence of immunodeficiency and autoimmune diseases both on a clinical and biological basis. RESULTS: Overall, 44.1% (78/177) and 58.2% (46/79) of KS individuals exhibited infection susceptibility and hypogammaglobulinemia, respectively; 13.6% (24/177) had autoimmune disease (AID; 25.6% [11/43] in adults), 5.6% (10/177) with ≥2 AID manifestations. The most frequent AID manifestations were immune thrombocytopenic purpura (7.3% [13/177]) and autoimmune hemolytic anemia (4.0% [7/177]). Among nonhematological manifestations, vitiligo was frequent. Immune thrombocytopenic purpura was frequent with missense versus other types of variants (p = 0.027). CONCLUSION: The high prevalence of immunopathological manifestations in KS demonstrates the importance of systematic screening and efficient preventive management of these treatable and sometimes life-threatening conditions.

Regulation of the positive transcriptional effect of PLZF through a non-canonical EZH2 activity.

The transcription factor PLZF (promyelocytic leukemia zinc finger protein) acts as an epigenetic regulator balancing self-renewal and differentiation of hematopoietic cells through binding to various chromatin-modifying factors. First described as a transcriptional repressor, PLZF is also associated with active transcription, although the molecular bases underlying the differences are unknown. Here, we reveal that in a hematopoietic cell line, PLZF is predominantly associated with transcribed genes. Additionally, we identify a new association between PLZF and the histone methyltransferase, EZH2 at the genomic level. We find that co-occupancy of PLZF and EZH2 on chromatin at PLZF target genes is not associated with SUZ12 or trimethylated lysine 27 of histone H3 (H3K27me3) but with the active histone mark H3K4me3 and active transcription. Removal of EZH2 leads to an increase of PLZF binding and increased gene expression. Our results suggest a new role of EZH2 in restricting PLZF positive transcriptional activity independently of its canonical PRC2 activity.

ARS2 is a general suppressor of pervasive transcription.

Termination of transcription is important for establishing gene punctuation marks. It is also critical for suppressing many of the pervasive transcription events occurring throughout eukaryotic genomes and coupling their RNA products to efficient decay. In human cells, the ARS2 protein has been implicated in such function as its depletion causes transcriptional read-through of selected gene terminators and because it physically interacts with the ribonucleolytic nuclear RNA exosome. Here, we study the role of ARS2 on transcription and RNA metabolism genome wide. We show that ARS2 depletion negatively impacts levels of promoter-proximal RNA polymerase II at protein-coding (pc) genes. Moreover, our results reveal a general role of ARS2 in transcription termination-coupled RNA turnover at short transcription units like snRNA-, replication-dependent histone-, promoter upstream transcript- and enhancer RNA-loci. Depletion of the ARS2 interaction partner ZC3H18 mimics the ARS2 depletion, although to a milder extent, whereas depletion of the exosome core subunit RRP40 only impacts RNA abundance post-transcriptionally. Interestingly, ARS2 is also involved in transcription termination events within first introns of pc genes. Our work therefore establishes ARS2 as a general suppressor of pervasive transcription with the potential to regulate pc gene expression.

The chromatin environment shapes DNA replication origin organization and defines origin classes.

To unveil the still-elusive nature of metazoan replication origins, we identified them genome-wide and at unprecedented high-resolution in mouse ES cells. This allowed initiation sites (IS) and initiation zones (IZ) to be differentiated. We then characterized their genetic signatures and organization and integrated these data with 43 chromatin marks and factors. Our results reveal that replication origins can be grouped into three main classes with distinct organization, chromatin environment, and sequence motifs. Class 1 contains relatively isolated, low-efficiency origins that are poor in epigenetic marks and are enriched in an asymmetric AC repeat at the initiation site. Late origins are mainly found in this class. Class 2 origins are particularly rich in enhancer elements. Class 3 origins are the most efficient and are associated with open chromatin and polycomb protein-enriched regions. The presence of Origin G-rich Repeated elements (OGRE) potentially forming G-quadruplexes (G4) was confirmed at most origins. These coincide with nucleosome-depleted regions located upstream of the initiation sites, which are associated with a labile nucleosome containing H3K64ac. These data demonstrate that specific chromatin landscapes and combinations of specific signatures regulate origin localization. They explain the frequently observed links between DNA replication and transcription. They also emphasize the plasticity of metazoan replication origins and suggest that in multicellular eukaryotes, the combination of distinct genetic features and chromatin configurations act in synergy to define and adapt the origin profile.

Dynamic recruitment of Ets1 to both nucleosome-occupied and -depleted enhancer regions mediates a transcriptional program switch during early T-cell differentiation.

Ets1 is a sequence-specific transcription factor that plays an important role during hematopoiesis, and is essential for the transition of CD4(-)/CD8(-) double negative (DN) to CD4(+)/CD8(+) double positive (DP) thymocytes. Using genome-wide and functional approaches, we investigated the binding properties, transcriptional role and chromatin environment of Ets1 during this transition. We found that while Ets1 binding at distal sites was associated with active genes at both DN and DP stages, its enhancer activity was attained at the DP stage, as reflected by levels of the core transcriptional hallmarks H3K4me1/3, RNA Polymerase II and eRNA. This dual, stage-specific ability reflected a switch from non-T hematopoietic toward T-cell specific gene expression programs during the DN-to-DP transition, as indicated by transcriptome analyses of Ets1(-/-) thymic cells. Coincidentally, Ets1 associates more specifically with Runx1 in DN and with TCF1 in DP cells. We also provide evidence that Ets1 predominantly binds distal nucleosome-occupied regions in DN and nucleosome-depleted regions in DP. Finally and importantly, we demonstrate that Ets1 induces chromatin remodeling by displacing H3K4me1-marked nucleosomes. Our results thus provide an original model whereby the ability of a transcription factor to bind nucleosomal DNA changes during differentiation with consequences on its cognate enhancer activity.

Pasha: a versatile R package for piling chromatin HTS data.

UNLABELLED: We describe an R package designed for processing aligned reads from chromatin-oriented high-throughput sequencing experiments. Pasha (preprocessing of aligned sequences from HTS analyses) allows easy manipulation of aligned reads from short-read sequencing technologies (ChIP-seq, FAIRE-seq, MNase-Seq, …) and offers innovative approaches such as ChIP-seq reads elongation, nucleosome midpoint piling strategy for positioning analyses, or the ability to subset paired-end reads by groups of insert size that can contain biologically relevant information. AVAILABILITY AND IMPLEMENTATION: Pasha is a multi-platform R package, available on CRAN repositories under GPL-3 license (https://cran.r-project.org/web/packages/Pasha/). CONTACTS: rfenouil@gmail.com or jean-christophe.andrau@igmm.cnrs.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcription-dependent generation of a specialized chromatin structure at the TCRß locus.

V(D)J recombination assembles Ag receptor genes during lymphocyte development. Enhancers at AR loci are known to control V(D)J recombination at associated alleles, in part by increasing chromatin accessibility of the locus, to allow the recombination machinery to gain access to its chromosomal substrates. However, whether there is a specific mechanism to induce chromatin accessibility at AR loci is still unclear. In this article, we highlight a specialized epigenetic marking characterized by high and extended H3K4me3 levels throughout the Dß-Jß-Cß gene segments. We show that extended H3K4 trimethylation at the Tcrb locus depends on RNA polymerase II (Pol II)-mediated transcription. Furthermore, we found that the genomic regions encompassing the two DJCß clusters are highly enriched for Ser(5)-phosphorylated Pol II and short-RNA transcripts, two hallmarks of transcription initiation and early transcription. Of interest, these features are shared with few other tissue-specific genes. We propose that the entire DJCß regions behave as transcription "initiation" platforms, therefore linking a specialized mechanism of Pol II transcription with extended H3K4 trimethylation and highly accessible Dß and Jß gene segments.

Threonine-4 of mammalian RNA polymerase II CTD is targeted by Polo-like kinase 3 and required for transcriptional elongation.

Eukaryotic RNA polymerase II (Pol II) has evolved an array of heptad repeats with the consensus sequence Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7 at the carboxy-terminal domain (CTD) of the large subunit (Rpb1). Differential phosphorylation of Ser2, Ser5, and Ser7 in the 5' and 3' regions of genes coordinates the binding of transcription and RNA processing factors to the initiating and elongating polymerase complexes. Here, we report phosphorylation of Thr4 by Polo-like kinase 3 in mammalian cells. ChIPseq analyses indicate an increase of Thr4-P levels in the 3' region of genes occurring subsequently to an increase of Ser2-P levels. A Thr4/Ala mutant of Pol II displays a lethal phenotype. This mutant reveals a global defect in RNA elongation, while initiation is largely unaffected. Since Thr4 replacement mutants are viable in yeast we conclude that this amino acid has evolved an essential function(s) in the CTD of Pol II for gene transcription in mammalian cells.

H3K4 tri-methylation provides an epigenetic signature of active enhancers.

Combinations of post-translational histone modifications shape the chromatin landscape during cell development in eukaryotes. However, little is known about the modifications exactly delineating functionally engaged regulatory elements. For example, although histone H3 lysine 4 mono-methylation (H3K4me1) indicates the presence of transcriptional gene enhancers, it does not provide clearcut information about their actual position and stage-specific activity. Histone marks were, therefore, studied here at genomic loci differentially expressed in early stages of T-lymphocyte development. The concomitant presence of the three H3K4 methylation states (H3K4me1/2/3) was found to clearly reflect the activity of bona fide T-cell gene enhancers. Globally, gain or loss of H3K4me2/3 at distal genomic regions correlated with, respectively, the induction or the repression of associated genes during T-cell development. In the Tcrb gene enhancer, the H3K4me3-to-H3K4me1 ratio decreases with the enhancer's strength. Lastly, enhancer association of RNA-polymerase II (Pol II) correlated with the presence of H3K4me3 and Pol II accumulation resulted in local increase of H3K4me3. Our results suggest the existence of functional links between Pol II occupancy, H3K4me3 enrichment and enhancer activity.

CpG islands and GC content dictate nucleosome depletion in a transcription-independent manner at mammalian promoters.

One clear hallmark of mammalian promoters is the presence of CpG islands (CGIs) at more than two-thirds of genes, whereas TATA boxes are only present at a minority of promoters. Using genome-wide approaches, we show that GC content and CGIs are major promoter elements in mammalian cells, able to govern open chromatin conformation and support paused transcription. First, we define three classes of promoters with distinct transcriptional directionality and pausing properties that correlate with their GC content. We further analyze the direct influence of GC content on nucleosome positioning and depletion and show that CpG content and CGI width correlate with nucleosome depletion both in vivo and in vitro. We also show that transcription is not essential for nucleosome exclusion but influences both a weak +1 and a well-positioned nucleosome at CGI borders. Altogether our data support the idea that CGIs have become an essential feature of promoter structure defining novel regulatory properties in mammals.