Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome

Low-affinity CTCF binding drives transcriptional regulation whereas high-affinity binding encompasses architectural functions

CTCF is a DNA-binding protein which plays critical roles in chromatin structure organization and transcriptional regulation; however, little is known about the functional determinants of different CTCF-binding sites (CBS). Using a conditional mouse model, we have identified one set of CBSs that are lost upon CTCF depletion (lost CBSs) and another set that persists (retained CBSs). Retained CBSs are more similar to the consensus CTCF-binding sequence and usually span tandem CTCF peaks. Lost CBSs are enriched at enhancers and promoters and associate with active chromatin marks and higher transcriptional activity. In contrast, retained CBSs are enriched at TAD and loop boundaries. Integration of ChIP-seq and RNA-seq data has revealed that retained CBSs are located at the boundaries between distinct chromatin states, acting as chromatin barriers. Our results provide evidence that transient, lost CBSs are involved in transcriptional regulation, whereas retained CBSs are critical for establishing higher-order chromatin architecture.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation

BACKGROUND

CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF.

RESULTS

We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF.

CONCLUSIONS

By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

Multiscale reorganization of the genome following DNA damage facilitates chromosome translocations via nuclear actin polymerization

Nuclear actin-based movements have been shown to orchestrate clustering of DNA double-strand breaks (DSBs) into homology-directed repair domains. Here we describe multiscale three-dimensional genome reorganization following DNA damage and analyze the contribution of the nuclear WASP-ARP2/3-actin pathway toward chromatin topology alterations and pathologic repair. Hi-C analysis reveals genome-wide, DNA damage-induced chromatin compartment flips facilitated by ARP2/3 that enrich for open, A compartments. Damage promotes interactions between DSBs, which in turn facilitate aberrant, actin-dependent intra- and inter-chromosomal rearrangements. Our work establishes that clustering of resected DSBs into repair domains by nuclear actin assembly is coordinated with multiscale alterations in genome architecture that enable homology-directed repair while also increasing nonhomologous end-joining-dependent translocation frequency.

CTCF: an R/bioconductor data package of human and mouse CTCF binding sites

Summary

CTCF (CCCTC-binding factor) is an 11-zinc-finger DNA binding protein which regulates much of the eukaryotic genome's 3D structure and function. The diversity of CTCF binding motifs has led to a fragmented landscape of CTCF binding data. We collected position weight matrices of CTCF binding motifs and defined strand-oriented CTCF binding sites in the human and mouse genomes, including the recent Telomere to Telomere and mm39 assemblies. We included selected experimentally determined and predicted CTCF binding sites, such as CTCF-bound cis-regulatory elements from SCREEN ENCODE. We recommend filtering strategies for CTCF binding motifs and demonstrate that liftOver is a viable alternative to convert CTCF coordinates between assemblies. Our comprehensive data resource and usage recommendations can serve to harmonize and strengthen the reproducibility of genomic studies utilizing CTCF binding data.

Availability and implementation

https://bioconductor.org/packages/CTCF. Companion website: https://dozmorovlab.github.io/CTCF/; Code to reproduce the analyses: https://github.com/dozmorovlab/CTCF.dev.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

CW198 acts as a genetic insulator to block enhancer-promoter interaction in plants

Insulators in vertebrates play a role in genome architecture and orchestrate temporo-spatial enhancer-promoter interactions. In plants, insulators and their associated binding factors have not been documented as of yet, largely as a result of a lack of characterized insulators. In this study, we took a comprehensive strategy to identify and validate the enhancer-blocking insulator CW198. We show that a 1.08-kb CW198 fragment from Arabidopsis can, when interposed between an enhancer and a promoter, efficiently abrogate the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and tobacco plants. In plants, both transcriptional crosstalk and spreading of histone modifications were rarely detectable across CW198, which resembles the insulation property observed across the CTCF insulator in the mammalian genome. Taken together, our findings support that CW198 acts as an enhancer-blocking insulator in both Arabidopsis and tobacco. The significance of the present findings and their relevance to the mitigation of mutual interference between enhancers and promoters, as well as multiple promoters in transgenes, is discussed.

CTCF loss induces giant lamellar bodies in Purkinje cell dendrites

CCCTC-binding factor (CTCF) has a key role in higher-order chromatin architecture that is important for establishing and maintaining cell identity by controlling gene expression. In the mature cerebellum, CTCF is highly expressed in Purkinje cells (PCs) as compared with other cerebellar neurons. The cerebellum plays an important role in motor function by regulating PCs, which are the sole output neurons, and defects in PCs cause motor dysfunction. However, the role of CTCF in PCs has not yet been explored. Here we found that the absence of CTCF in mouse PCs led to progressive motor dysfunction and abnormal dendritic morphology in those cells, which included dendritic self-avoidance defects and a proximal shift in the climbing fibre innervation territory on PC dendrites. Furthermore, we found the peculiar lamellar structures known as "giant lamellar bodies" (GLBs), which have been reported in PCs of patients with Werdnig-Hoffman disease, 13q deletion syndrome, and Krabbe disease. GLBs are localized to PC dendrites and are assumed to be associated with neurodegeneration. They have been noted, however, only in case reports following autopsy, and reports of their existence have been very limited. Here we show that GLBs were reproducibly formed in PC dendrites of a mouse model in which CTCF was deleted. GLBs were not noted in PC dendrites at infancy but instead developed over time. In conjunction with GLB development in PC dendrites, the endoplasmic reticulum was almost absent around the nuclei, the mitochondria were markedly swollen and their cristae had decreased drastically, and almost all PCs eventually disappeared as severe motor deficits manifested. Our results revealed the important role of CTCF during normal development and in maintaining PCs and provide new insights into the molecular mechanism of GLB formation during neurodegenerative disease.

CLIMB: High-dimensional association detection in large scale genomic data

Joint analyses of genomic datasets obtained in multiple different conditions are essential for understanding the biological mechanism that drives tissue-specificity and cell differentiation, but they still remain computationally challenging. To address this we introduce CLIMB (Composite LIkelihood eMpirical Bayes), a statistical methodology that learns patterns of condition-specificity present in genomic data. CLIMB provides a generic framework facilitating a host of analyses, such as clustering genomic features sharing similar condition-specific patterns and identifying which of these features are involved in cell fate commitment. We apply CLIMB to three sets of hematopoietic data, which examine CTCF ChIP-seq measured in 17 different cell populations, RNA-seq measured across constituent cell populations in three committed lineages, and DNase-seq in 38 cell populations. Our results show that CLIMB improves upon existing alternatives in statistical precision, while capturing interpretable and biologically relevant clusters in the data.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome

BACKGROUND

Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome.

RESULTS

Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively.

CONCLUSIONS

Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

The Epstein-Barr Virus Enhancer Interaction Landscapes in Virus-Associated Cancer Cell Lines

Epstein-Barr virus (EBV) persists in human cells as episomes. EBV episomes are chromatinized and their 3D conformation varies greatly in cells expressing different latency genes. We used HiChIP, an assay which combines genome-wide chromatin conformation capture followed by deep sequencing (Hi-C) and chromatin immunoprecipitation (ChIP), to interrogate the EBV episome 3D conformation in different cancer cell lines. In an EBV-transformed lymphoblastoid cell line (LCL) GM12878 expressing type III EBV latency genes, abundant genomic interactions were identified by H3K27ac HiChIP. A strong enhancer was located near the BILF2 gene and looped to multiple genes around BALFs loci. Perturbation of the BILF2 enhancer by CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) altered the expression of BILF2 enhancer-linked genes, including BARF0 and BALF2, suggesting that this enhancer regulates the expression of linked genes. H3K27ac ChIP followed by deep sequencing (ChIP-seq) identified several strong EBV enhancers in T/NK (natural killer) lymphoma cells that express type II EBV latency genes. Extensive intragenomic interactions were also found which linked enhancers to target genes. A strong enhancer at BILF2 also looped to the BALF loci. CRISPRi also validated the functional connection between BILF2 enhancer and BARF1 gene. In contrast, H3K27ac HiChIP found significantly fewer intragenomic interactions in type I EBV latency gene-expressing primary effusion lymphoma (PEL) cell lines. These data provided new insight into the regulation of EBV latency gene expression in different EBV-associated tumors. IMPORTANCE EBV is the first human DNA tumor virus identified, discovered over 50 years ago. EBV causes ~200,000 cases of various cancers each year. EBV-encoded oncogenes, noncoding RNAs, and microRNAs (miRNAs) can promote cell growth and survival and suppress senescence. Regulation of EBV gene expression is very complex. The viral C promoter regulates the expression of all EBV nuclear antigens (EBNAs), some of which are very far away from the C promoter. Another way by which the virus activates remote gene expression is through DNA looping. In this study, we describe the viral genome looping patterns in various EBV-associated cancer cell lines and identify important EBV enhancers in these cells. This study also identified novel opportunities to perturb and eventually control EBV gene expression in these cancer cells.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Characterization of sequence determinants of enhancer function using natural genetic variation

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F₁-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

Butyrate Induces Modifications of the CTCF-Binding Landscape in Cattle Cells

Butyrate is produced in the rumen from microbial fermentation and is related to several functions, including cell differentiation and proliferation. Butyrate supplementation in calves can accelerate rumen development. DNA-protein interactions, such as the CCCTC-binding factor (CTCF), play essential roles in chromatin organization and gene expression regulation. Although CTCF-binding sites have been identified recently in cattle, a deeper characterization, including differentially CTCF-binding sites (DCBS), is vital for a better understanding of butyrate’s role in the chromatin landscape. This study aimed to identify CTCF-binding regions and DCBS under a butyrate-induced condition using ChIP-seq in bovine cells; 61,915 CTCF peaks were identified in the butyrate and 51,347 in the control. From these regions, 2265 DCBS were obtained for the butyrate vs. control comparison, comprising ~90% of induced sites. Most of the butyrate DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment showed crucial terms for the induced DCBS, mainly related to cellular proliferation, cell adhesion, and growth regulation. Interestingly, the ECM-receptor interaction pathway was observed for the induced DCBS. Motif enrichment analysis further identified transcription factors, including CTCF, BORIS, TGIF2, and ZIC3. When DCBS was integrated with RNA-seq data, putative genes were identified for the repressed DCBS, including GATA4. Our study revealed promising candidate genes in bovine cells by a butyrate-induced condition that might be related to the regulation of rumen development, such as integrins, keratins, and collagens. These results provide a better understanding of the function of butyrate in cattle rumen development and chromatin landscape regulation.

Differentially CTCF-Binding Sites in Cattle Rumen Tissue during Weaning

The weaning transition in calves is characterized by major structural changes such as an increase in the rumen capacity and surface area due to diet changes. Studies evaluating rumen development in calves are vital to identify genetic mechanisms affected by weaning. This study aimed to provide a genome-wide characterization of CTCF-binding sites and differentially CTCF-binding sites (DCBS) in rumen tissue during the weaning transition of four Holstein calves to uncover regulatory elements in rumen epithelial tissue using ChIP-seq. Our study generated 67,280 CTCF peaks for the before weaning (BW) and 39,891 for after weaning (AW). Then, 7401 DCBS were identified for the AW vs. BW comparison representing 0.15% of the cattle genome, comprising ~54% of induced DCBS and ~46% of repressed DCBS. Most of the induced and repressed DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment revealed many shared GO terms for the induced and the repressed DCBS, mainly related to cellular migration, proliferation, growth, differentiation, cellular adhesion, digestive tract morphogenesis, and response to TGFβ. In addition, shared KEGG pathways were obtained for adherens junction and focal adhesion. Interestingly, other relevant KEGG pathways were observed for the induced DCBS like gastric acid secretion, salivary secretion, bacterial invasion of epithelial cells, apelin signaling, and mucin-type O-glycan biosynthesis. IPA analysis further revealed pathways with potential roles in rumen development during weaning, including TGFβ, Integrin-linked kinase, and Integrin signaling. When DCBS were further integrated with RNA-seq data, 36 putative target genes were identified for the repressed DCBS, including KRT84, COL9A2, MATN3, TSPAN1, and AJM1. This study successfully identified DCBS in cattle rumen tissue after weaning on a genome-wide scale and revealed several candidate target genes that may have a role in rumen development, such as TGFβ, integrins, keratins, and SMADs. The information generated in this preliminary study provides new insights into bovine genome regulation and chromatin landscape.

Single base-pair resolution analysis of DNA binding motif with MoMotif reveals an oncogenic function of CTCF zinc-finger 1 mutation

Defining the impact of missense mutations on the recognition of DNA motifs is highly dependent on bioinformatic tools that define DNA binding elements. However, classical motif analysis tools remain limited in their capacity to identify subtle changes in complex binding motifs between distinct conditions. To overcome this limitation, we developed a new tool, MoMotif, that facilitates a sensitive identification, at the single base-pair resolution, of complex, or subtle, alterations to core binding motifs, discerned from ChIP-seq data. We employed MoMotif to define the previously uncharacterized recognition motif of CTCF zinc-finger 1 (ZF1), and to further define the impact of CTCF ZF1 mutation on its association with chromatin. Mutations of CTCF ZF1 are exclusive to breast cancer and are associated with metastasis and therapeutic resistance, but the underlying mechanisms are unclear. Using MoMotif, we identified an extension of the CTCF core binding motif, necessitating a functional ZF1 to bind appropriately. Using a combination of ChIP-Seq and RNA-Seq, we discover that the inability to bind this extended motif drives an altered transcriptional program associated with the oncogenic phenotypes observed clinically. Our study demonstrates that MoMotif is a powerful new tool for comparative ChIP-seq analysis and characterising DNA-protein contacts.

CTCF and Its Partners: Shaper of 3D Genome during Development

The 3D genome organization and its dynamic modulate genome function, playing a pivotal role in cell differentiation and development. CTCF and cohesin, acting as the core architectural components involved in chromatin looping and genome folding, can also recruit other protein or RNA partners to fine-tune genome structure during development. Moreover, systematic screening for partners of CTCF has been performed through high-throughput approaches. In particular, several novel protein and RNA partners, such as BHLHE40, WIZ, MAZ, Aire, MyoD, YY1, ZNF143, and Jpx, have been identified, and these partners are mostly implicated in transcriptional regulation and chromatin remodeling, offering a unique opportunity for dissecting their roles in higher-order chromatin organization by collaborating with CTCF and cohesin. Here, we review the latest advancements with an emphasis on features of CTCF partners and also discuss the specific functions of CTCF-associated complexes in chromatin structure modulation, which may extend our understanding of the functions of higher-order chromatin architecture in developmental processes.

Satellite DNAs and human sex chromosome variation

Satellite DNAs are present on every chromosome in the cell and are typically enriched in repetitive, heterochromatic parts of the human genome. Sex chromosomes represent a unique genomic and epigenetic context. In this review, we first report what is known about satellite DNA biology on human X and Y chromosomes, including repeat content and organization, as well as satellite variation in typical euploid individuals. Then, we review sex chromosome aneuploidies that are among the most common types of aneuploidies in the general population, and are better tolerated than autosomal aneuploidies. This is demonstrated also by the fact that aging is associated with the loss of the X, and especially the Y chromosome. In addition, supernumerary sex chromosomes enable us to study general processes in a cell, such as analyzing heterochromatin dosage (i.e. additional Barr bodies and long heterochromatin arrays on Yq) and their downstream consequences. Finally, genomic and epigenetic organization and regulation of satellite DNA could influence chromosome stability and lead to aneuploidy. In this review, we argue that the complete annotation of satellite DNA on sex chromosomes in human, and especially in centromeric regions, will aid in explaining the prevalence and the consequences of sex chromosome aneuploidies.

Regulation of 3D Organization and Its Role in Cancer Biology

Three-dimensional (3D) genomics is the frontier field in the post-genomics era, its foremost content is the relationship between chromatin spatial conformation and regulation of gene transcription. Cancer biology is a complex system resulting from genetic alterations in key tumor oncogenes and suppressor genes for cell proliferation, DNA replication, cell differentiation, and homeostatic functions. Although scientific research in recent decades has revealed how the genome sequence is mutated in many cancers, high-order chromosomal structures involved in the development and fate of cancer cells represent a crucial but rarely explored aspect of cancer genomics. Hence, dissection of the 3D genome conformation of cancer helps understand the unique epigenetic patterns and gene regulation processes that distinguish cancer biology from normal physiological states. In recent years, research in tumor 3D genomics has grown quickly. With the rapid progress of 3D genomics technology, we can now better determine the relationship between cancer pathogenesis and the chromatin structure of cancer cells. It is becoming increasingly explicit that changes in 3D chromatin structure play a vital role in controlling oncogene transcription. This review focuses on the relationships between tumor gene expression regulation, tumor 3D chromatin structure, and cancer phenotypic plasticity. Furthermore, based on the functional consequences of spatial disorganization in the cancer genome, we look forward to the clinical application prospects of 3D genomic biomarkers.

3D Genome Organization as an Epigenetic Determinant of Transcription Regulation in T Cells

In the heart of innate and adaptive immunity lies the proper spatiotemporal development of several immune cell lineages. Multiple studies have highlighted the necessity of epigenetic and transcriptional regulation in cell lineage specification. This mode of regulation is mediated by transcription factors and chromatin remodelers, controlling developmentally essential gene sets. The core of transcription and epigenetic regulation is formulated by different epigenetic modifications determining gene expression. Apart from “classic” epigenetic modifications, 3D chromatin architecture is also purported to exert fundamental roles in gene regulation. Chromatin conformation both facilitates cell-specific factor binding at specified regions and is in turn modified as such, acting synergistically. The interplay between global and tissue-specific protein factors dictates the epigenetic landscape of T and innate lymphoid cell (ILC) lineages. The expression of global genome organizers such as CTCF, YY1, and the cohesin complexes, closely cooperate with tissue-specific factors to exert cell type-specific gene regulation. Special AT-rich binding protein 1 (SATB1) is an important tissue-specific genome organizer and regulator controlling both long- and short-range chromatin interactions. Recent indications point to SATB1’s cooperation with the aforementioned factors, linking global to tissue-specific gene regulation. Changes in 3D genome organization are of vital importance for proper cell development and function, while disruption of this mechanism can lead to severe immuno-developmental defects. Newly emerging data have inextricably linked chromatin architecture deregulation to tissue-specific pathophysiological phenotypes. The combination of these findings may shed light on the mechanisms behind pathological conditions.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Prediction and comparative analysis of CTCF binding sites based on a first principle approach

We calculated the patterns for the CCCTC transcription factor (CTCF) binding sites across many genomes on a first principle approach. The validation of the first principle method was done on the human as well as on the mouse genome. The predicted human CTCF binding sites are consistent with the consensus sequence, ChIP-seq data for the K562 cell, nucleosome positions for IMR90 cell as well as the CTCF binding sites in the mouse HOXA gene. The analysis of Homo sapiens, Mus musculus, Sus scrofa, Capra hircus and Drosophila melanogaster whole genomes shows: binding sites are organized in cluster-like groups, where two consecutive sites obey a power-law with coefficient ranging from to 0.3292 0.0068 to 0.5409 0.0064; the distance between these groups varies from 18.08 0.52kbp to 42.1 2.0kbp. The genome of Aedes aegypti does not show a power law, but 19.9% of binding sites are 144 4 and 287 5bp distant of each other. We run negative tests, confirming the under-representation of CTCF binding sites in Caenorhabditis elegans, Plasmodium falciparum and Arabidopsis thaliana complete genomes.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

The Role of the Tumor Suppressor Gene Protein Tyrosine Phosphatase Gamma in Cancer

Members of the Protein Tyrosine Phosphatase (PTPs) family are associated with growth regulation and cancer development. Acting as natural counterpart of tyrosine kinases (TKs), mainly involved in crucial signaling pathways such as regulation of cell cycle, proliferation, invasion and angiogenesis, they represent key parts of complex physiological homeostatic mechanisms. Protein tyrosine phosphatase gamma (PTPRG) is classified as a R5 of the receptor type (RPTPs) subfamily and is broadly expressed in various isoforms in different tissues. PTPRG is considered a tumor-suppressor gene (TSG) mapped on chromosome 3p14-21, a region frequently subject to loss of heterozygosity in various tumors. However, reported mechanisms of PTPRG downregulation include missense mutations, ncRNA gene regulation and epigenetic silencing by hypermethylation of CpG sites on promoter region causing loss of function of the gene product. Inactive forms or total loss of PTPRG protein have been described in sporadic and Lynch syndrome colorectal cancer, nasopharyngeal carcinoma, ovarian, breast, and lung cancers, gastric cancer or diseases affecting the hematopoietic compartment as Lymphoma and Leukemia. Noteworthy, in Central Nervous System (CNS) PTPRZ/PTPRG appears to be crucial in maintaining glioblastoma cell-related neuronal stemness, carving out a pathological functional role also in this tissue. In this review, we will summarize the current knowledge on the role of PTPRG in various human cancers.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Systematic identification of genomic elements that regulate FCGR2A expression and harbor variants linked with autoimmune disease

BACKGROUND

FCGR2A binds antibody-antigen complexes to regulate the abundance of circulating and deposited complexes along with downstream immune and autoimmune responses. Although the abundance of FCRG2A may be critical in immune-mediated diseases, little is known about whether its surface expression is regulated through cis genomic elements and non-coding variants. In the current study, we aimed to characterize the regulation of FCGR2A expression, the impact of genetic variation and its association with autoimmune disease.

METHODS

We applied CRISPR-based interference and editing to scrutinize 1.7 Mb of open chromatin surrounding the FCGR2A gene to identify regulatory elements. Relevant transcription factors (TFs) binding to these regions were defined through public databases. Genetic variants affecting regulation were identified using luciferase reporter assays and were verified in a cohort of 1996 genotyped healthy individuals using flow cytometry.

RESULTS

We identified a complex proximal region and five distal enhancers regulating FCGR2A. The proximal region split into subregions upstream and downstream of the transcription start site, was enriched in binding of inflammation-regulated TFs, and harbored a variant associated with FCGR2A expression in primary myeloid cells. One distal enhancer region was occupied by CCCTC-binding factor (CTCF) whose binding site was disrupted by a rare genetic variant, altering gene expression.

CONCLUSIONS

The FCGR2A gene is regulated by multiple proximal and distal genomic regions, with links to autoimmune disease. These findings may open up novel therapeutic avenues where fine-tuning of FCGR2A levels may constitute a part of treatment strategies for immune-mediated diseases.

Three-Dimensional Genome Organization in Breast and Gynecological Cancers: How Chromatin Folding Influences Tumorigenic Transcriptional Programs

A growing body of research on the transcriptome and cancer genome has demonstrated that many gynecological tumor-specific gene mutations are located in cis-regulatory elements. Through chromosomal looping, cis-regulatory elements interact which each other to control gene expression by bringing distant regulatory elements, such as enhancers and insulators, into close proximity with promoters. It is well known that chromatin connections may be disrupted in cancer cells, promoting transcriptional dysregulation and the expression of abnormal tumor suppressor genes and oncogenes. In this review, we examine the roles of alterations in 3D chromatin interactions. This includes changes in CTCF protein function, cancer-risk single nucleotide polymorphisms, viral integration, and hormonal response as part of the mechanisms that lead to the acquisition of enhancers or super-enhancers. The translocation of existing enhancers, as well as enhancer loss or acquisition of insulator elements that interact with gene promoters, is also revised. Remarkably, similar processes that modify 3D chromatin contacts in gene promoters may also influence the expression of non-coding RNAs, such as long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), which have emerged as key regulators of gene expression in a variety of cancers, including gynecological malignancies.

Jpx RNA regulates CTCF anchor site selection and formation of chromosome loops

Chromosome loops shift dynamically during development, homeostasis, and disease. CCCTC-binding factor (CTCF) is known to anchor loops and construct 3D genomes, but how anchor sites are selected is not yet understood. Here, we unveil Jpx RNA as a determinant of anchor selectivity. Jpx RNA targets thousands of genomic sites, preferentially binding promoters of active genes. Depleting Jpx RNA causes ectopic CTCF binding, massive shifts in chromosome looping, and downregulation of >700 Jpx target genes. Without Jpx, thousands of lost loops are replaced by de novo loops anchored by ectopic CTCF sites. Although Jpx controls CTCF binding on a genome-wide basis, it acts selectively at the subset of developmentally sensitive CTCF sites. Specifically, Jpx targets low-affinity CTCF motifs and displaces CTCF protein through competitive inhibition. We conclude that Jpx acts as a CTCF release factor and shapes the 3D genome by regulating anchor site usage.

In silico approach to understand epigenetics of POTEE in ovarian cancer

Ovarian cancer is the third leading cause of cancer-related deaths in India. Epigenetics mechanisms seemingly plays an important role in ovarian cancer. This paper highlights the crucial epigenetic changes that occur in POTEE that get hypomethylated in ovarian cancer. We utilized the POTEE paralog mRNA sequence to identify major motifs and also performed its enrichment analysis. We identified 6 motifs of varying lengths, out of which only three motifs, including CTTCCAGCAGATGTGGATCA, GGAACTGCC, and CGCCACATGCAGGC were most likely to be present in the nucleotide sequence of POTEE. By enrichment and occurrences identification analyses, we rectified the best match motif as CTTCCAGCAGATGT. Since there is no experimentally verified structure of POTEE paralog, thus, we predicted the POTEE structure using an automated workflow for template-based modeling using the power of a deep neural network. Additionally, to validate our predicted model we used AlphaFold predicted POTEE structure and observed that the residual stretch starting from 237-958 had a very high confidence per residue. Furthermore, POTEE predicted model stability was evaluated using replica exchange molecular dynamic simulation for 50 ns. Our network-based epigenetic analysis discerns only 10 highly significant, direct, and physical associators of POTEE. Our finding aims to provide new insights about the POTEE paralog.

Nup93 and CTCF modulate spatiotemporal dynamics and function of the HOXA gene locus during differentiation

Nucleoporins regulate nuclear transport and are also involved in DNA damage, repair, cell cycle, chromatin organization, and gene expression. Here, we studied the role of nucleoporin Nup93 and the chromatin organizer CTCF in regulating HOXA expression during differentiation. ChIP sequencing revealed a significant overlap between Nup93 and CTCF peaks. Interestingly, Nup93 and CTCF are associated with the 3' and 5'HOXA genes respectively. Depletions of Nup93 and CTCF antagonistically modulate expression levels of 3'and 5'HOXA genes in undifferentiated NT2/D1 cells. Nup93 also regulates the localization of the HOXA gene locus, which disengages from the nuclear periphery upon Nup93 but not CTCF depletion, consistent with its upregulation. The dynamic association of Nup93 and CTCF with the HOXA locus during differentiation correlates with its spatial positioning and expression. While Nup93 tethers the HOXA locus to the nuclear periphery, CTCF potentially regulates looping of the HOXA gene cluster in a temporal manner. In summary, Nup93 and CTCF complement one another in modulating the spatiotemporal dynamics and function of the HOXA gene locus during differentiation.

Decoding the Equine Genome: Lessons from ENCODE

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Imprinting of canine IGF2 and H19

Genomic imprinting occurs in therian mammals and is a phenomenon whereby the two alleles of a gene are differentially expressed, based on the sex of the parent from whom the alleles were inherited. The allelic differences in expression are the consequence of different epigenetic modifications that are established in the sperm or oocyte during gametogenesis and transmitted at fertilization to offspring. A small minority of genes is regulated in this way but they have important biological functions, and aberrant regulation of imprinted genes contributes to disease aetiology in humans and other animals. The factors driving the evolution of imprinted genes are also of considerable interest, as these genes appear to forego the benefits of diploidy. To broaden the phylogenetic analysis of genomic imprinting, we began a study of imprinted genes in the domestic dog, Canis familiaris. In this report, we show that canine IGF2 and H19 are imprinted, with parent-of origin-dependent monoallelic expression patterns in neonatal umbilical cord. We identify a putative imprint control region associated with the genes, and provide evidence for differential methylation of this region in a somatic tissue (umbilical cord) and for its hypermethylation in the male germline. Canis familiaris is fast becoming a highly informative system for elucidating disease processes and evolution, and the study of imprinted genes in this species may help in understanding how these genes contribute to the generation of morphological and behavioral diversity.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

Short- and long-range cis interactions between integrated HPV genomes and cellular chromatin dysregulate host gene expression in early cervical carcinogenesis

Development of cervical cancer is directly associated with integration of human papillomavirus (HPV) genomes into host chromosomes and subsequent modulation of HPV oncogene expression, which correlates with multi-layered epigenetic changes at the integrated HPV genomes. However, the process of integration itself and dysregulation of host gene expression at sites of integration in our model of HPV16 integrant clone natural selection has remained enigmatic. We now show, using a state-of-the-art 'HPV integrated site capture' (HISC) technique, that integration likely occurs through microhomology-mediated repair (MHMR) mechanisms via either a direct process, resulting in host sequence deletion (in our case, partially homozygously) or via a 'looping' mechanism by which flanking host regions become amplified. Furthermore, using our 'HPV16-specific Region Capture Hi-C' technique, we have determined that chromatin interactions between the integrated virus genome and host chromosomes, both at short- (<500 kbp) and long-range (>500 kbp), appear to drive local host gene dysregulation through the disruption of host:host interactions within (but not exceeding) host structures known as topologically associating domains (TADs). This mechanism of HPV-induced host gene expression modulation indicates that integration of virus genomes near to or within a 'cancer-causing gene' is not essential to influence their expression and that these modifications to genome interactions could have a major role in selection of HPV integrants at the early stage of cervical neoplastic progression.

Dynamic CTCF binding directly mediates interactions among cis-regulatory elements essential for hematopoiesis

While constitutive CCCTC-binding factor (CTCF)-binding sites are needed to maintain relatively invariant chromatin structures, such as topologically associating domains, the precise roles of CTCF to control cell-type-specific transcriptional regulation remain poorly explored. We examined CTCF occupancy in different types of primary blood cells derived from the same donor to elucidate a new role for CTCF in gene regulation during blood cell development. We identified dynamic, cell-type-specific binding sites for CTCF that colocalize with lineage-specific transcription factors. These dynamic sites are enriched for single-nucleotide polymorphisms that are associated with blood cell traits in different linages, and they coincide with the key regulatory elements governing hematopoiesis. CRISPR-Cas9-based perturbation experiments demonstrated that these dynamic CTCF-binding sites play a critical role in red blood cell development. Furthermore, precise deletion of CTCF-binding motifs in dynamic sites abolished interactions of erythroid genes, such as RBM38, with their associated enhancers and led to abnormal erythropoiesis. These results suggest a novel, cell-type-specific function for CTCF in which it may serve to facilitate interaction of distal regulatory emblements with target promoters. Our study of the dynamic, cell-type-specific binding and function of CTCF provides new insights into transcriptional regulation during hematopoiesis.

A universal framework for detecting cis-regulatory diversity in DNA regulatory regions

High-throughput sequencing-based assays measure different biochemical activities pertaining to gene regulation, genome-wide. These activities include transcription factor (TF)–DNA binding, enhancer activity, open chromatin, and more. A major goal is to understand underlying sequence components, or motifs, that can explain the measured activity. It is usually not one motif but a combination of motifs bound by cooperatively acting proteins that confers activity to such regions. Furthermore, regions can be diverse, governed by different combinations of TFs/motifs. Current approaches do not take into account this issue of combinatorial diversity. We present a new statistical framework, cisDIVERSITY, which models regions as diverse modules characterized by combinations of motifs while simultaneously learning the motifs themselves. Because cisDIVERSITY does not rely on knowledge of motifs, modules, cell type, or organism, it is general enough to be applied to regions reported by most high-throughput assays. For example, in enhancer predictions resulting from different assays—GRO-cap, STARR-seq, and those measuring chromatin structure—cisDIVERSITY discovers distinct modules and combinations of TF binding sites, some specific to the assay. From protein–DNA binding data, cisDIVERSITY identifies potential cofactors of the profiled TF, whereas from ATAC-seq data, it identifies tissue-specific regulatory modules. Finally, analysis of single-cell ATAC-seq data suggests that regions open in one cell-state encode information about future states, with certain modules staying open and others closing down in the next time point.

Cocaine-related DNA methylation in caudate neurons alters 3D chromatin structure of the IRXA gene cluster

Epigenetic mechanisms, like those involving DNA methylation, are thought to mediate the relationship between chronic cocaine dependence and molecular changes in addiction-related neurocircuitry, but have been understudied in human brain. We initially used reduced representation bisulfite sequencing (RRBS) to generate a methylome-wide profile of cocaine dependence in human post-mortem caudate tissue. We focused on the Iroquois Homeobox A (IRXA) gene cluster, where hypomethylation in exon 3 of IRX2 in neuronal nuclei was associated with cocaine dependence. We replicated this finding in an independent cohort and found similar results in the dorsal striatum from cocaine self-administering mice. Using epigenome editing and 3C assays, we demonstrated a causal relationship between methylation within the IRX2 gene body, CTCF protein binding, three-dimensional (3D) chromatin interaction, and gene expression. Together, these findings suggest that cocaine-related hypomethylation of IRX2 contributes to the development and maintenance of cocaine dependence through alterations in 3D chromatin structure in the caudate nucleus.

DNA G-Quadruplexes Contribute to CTCF Recruitment

G-quadruplex (G4) sites in the human genome frequently colocalize with CCCTC-binding factor (CTCF)-bound sites in CpG islands (CGIs). We aimed to clarify the role of G4s in CTCF positioning. Molecular modeling data suggested direct interactions, so we performed in vitro binding assays with quadruplex-forming sequences from CGIs in the human genome. G4s bound CTCF with Kd values similar to that of the control duplex, while respective i-motifs exhibited no affinity for CTCF. Using ChIP-qPCR assays, we showed that G4-stabilizing ligands enhance CTCF occupancy at a G4-prone site in STAT3 gene. In view of the reportedly increased CTCF affinity for hypomethylated DNA, we next questioned whether G4s also facilitate CTCF recruitment to CGIs via protecting CpG sites from methylation. Bioinformatics analysis of previously published data argued against such a possibility. Finally, we questioned whether G4s facilitate CTCF recruitment by affecting chromatin structure. We showed that three architectural chromatin proteins of the high mobility group colocalize with G4s in the genome and recognize parallel-stranded or mixed-topology G4s in vitro. One of such proteins, HMGN3, contributes to the association between G4s and CTCF according to our bioinformatics analysis. These findings support both direct and indirect roles of G4s in CTCF recruitment.

The combined action of CTCF and its testis-specific paralog BORIS is essential for spermatogenesis

CTCF is a key organizer of the 3D genome. Its specialized paralog, BORIS, heterodimerizes with CTCF but is expressed only in male germ cells and in cancer states. Unexpectedly, BORIS-null mice have only minimal germ cell defects. To understand the CTCF-BORIS relationship, mouse models with varied CTCF and BORIS levels were generated. Whereas Ctcf+/+Boris+/+, Ctcf+/-Boris+/+, and Ctcf+/+Boris-/- males are fertile, Ctcf+/-Boris-/- (Compound Mutant; CM) males are sterile. Testes with combined depletion of both CTCF and BORIS show reduced size, defective meiotic recombination, increased apoptosis, and malformed spermatozoa. Although CM germ cells exhibit only 25% of CTCF WT expression, chromatin binding of CTCF is preferentially lost from CTCF-BORIS heterodimeric sites. Furthermore, CM testes lose the expression of a large number of spermatogenesis genes and gain the expression of developmentally inappropriate genes that are "toxic" to fertility. Thus, a combined action of CTCF and BORIS is required to both repress pre-meiotic genes and activate post-meiotic genes for a complete spermatogenesis program.

Accurate prediction of cis-regulatory modules reveals a prevalent regulatory genome of humans

cis-regulatory modules(CRMs) formed by clusters of transcription factor (TF) binding sites (TFBSs) are as important as coding sequences in specifying phenotypes of humans. It is essential to categorize all CRMs and constituent TFBSs in the genome. In contrast to most existing methods that predict CRMs in specific cell types using epigenetic marks, we predict a largely cell type agonistic but more comprehensive map of CRMs and constituent TFBSs in the gnome by integrating all available TF ChIP-seq datasets. Our method is able to partition 77.47% of genome regions covered by available 6092 datasets into a CRM candidate (CRMC) set (56.84%) and a non-CRMC set (43.16%). Intriguingly, the predicted CRMCs are under strong evolutionary constraints, while the non-CRMCs are largely selectively neutral, strongly suggesting that the CRMCs are likely cis-regulatory, while the non-CRMCs are not. Our predicted CRMs are under stronger evolutionary constraints than three state-of-the-art predictions (GeneHancer, EnhancerAtlas and ENCODE phase 3) and substantially outperform them for recalling VISTA enhancers and non-coding ClinVar variants. We estimated that the human genome might encode about 1.47M CRMs and 68M TFBSs, comprising about 55% and 22% of the genome, respectively; for both of which, we predicted 80%. Therefore, the cis-regulatory genome appears to be more prevalent than originally thought.

Myod1 and GR coordinate myofiber-specific transcriptional enhancers

Skeletal muscle is a dynamic tissue the size of which can be remodeled through the concerted actions of various cues. Here, we investigated the skeletal muscle transcriptional program and identified key tissue-specific regulatory genetic elements. Our results show that Myod1 is bound to numerous skeletal muscle enhancers in collaboration with the glucocorticoid receptor (GR) to control gene expression. Remarkably, transcriptional activation controlled by these factors occurs through direct contacts with the promoter region of target genes, via the CpG-bound transcription factor Nrf1, and the formation of Ctcf-anchored chromatin loops, in a myofiber-specific manner. Moreover, we demonstrate that GR negatively controls muscle mass and strength in mice by down-regulating anabolic pathways. Taken together, our data establish Myod1, GR and Nrf1 as key players of muscle-specific enhancer-promoter communication that orchestrate myofiber size regulation.

Learning the Regulatory Code of Gene Expression

Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.

GAD1 expression and its methylation as indicators of malignant behavior in thymic epithelial tumors

Thymic epithelial tumors (TETs) comprise thymomas and thymic carcinoma (TC). TC has more aggressive features and a poorer prognosis than thymomas. Genetic and epigenetic alterations in thymomas and TC have been investigated in an attempt to identify novel target molecules for TC. In the present study, genome-wide screening was performed on aberrantly methylated CpG islands in thymomas and TC, and the glutamate decarboxylase 1 gene (GAD1) was identified as the 4th significantly hypermethylated CpG island in TC compared with thymomas. GAD1 catalyzes the production of γ-aminobutyric acid from L-glutamic acid. GAD1 expression is abundant in the brain but rare in other tissues, including the thymus. A total of 73 thymomas and 17 TC tissues were obtained from 90 patients who underwent surgery or biopsy at Tokushima University Hospital between 1990 and 2017. DNA methylation was examined by bisulfite pyrosequencing, and the mRNA and protein expression levels of GAD1 were analyzed using reverse transcription-quantitative PCR and immunohistochemistry, respectively. The DNA methylation levels of GAD1 were significantly higher in TC tissues than in the normal thymus and thymoma tissues, and GAD1 methylation exhibited high sensitivity and specificity for discriminating between TC and thymoma. The mRNA and protein expression levels of GAD1 were significantly higher in TC tissues than in thymomas. Patients with TET with high GAD1 DNA hypermethylation and high mRNA and protein expression levels had significantly shorter relapse-free survival rates than those with low levels. In conclusion, significantly more epigenetic alterations were observed in TC tissues compared with in thymomas, which may contribute to the clinical features and prognosis of patients.

Chromatin-mediated epigenetic regulation of HSV-1 transcription as a potential target in antiviral therapy

The ability to establish, and reactivate from, latent infections is central to the biology and pathogenesis of HSV-1. It also poses a strong challenge to antiviral therapy, as latent HSV-1 genomes do not replicate or express any protein to be targeted. Although the processes regulating the establishment and maintenance of, and reactivation from, latency are not fully elucidated, the current general consensus is that epigenetics play a major role. A unifying model postulates that whereas HSV-1 avoids or counteracts chromatin silencing in lytic infections, it becomes silenced during latency, silencing which is somewhat disrupted during reactivation. Many years of work by different groups using a variety of approaches have also shown that the lytic HSV-1 chromatin is distinct and has unique biophysical properties not shared with most cellular chromatin. Nonetheless, the lytic and latent viral chromatins are typically enriched in post translational modifications or histone variants characteristic of active or repressed transcription, respectively. Moreover, a variety of small molecule epigenetic modulators inhibit viral replication and reactivation from latency. Despite these successes in culture and animal models, it is not obvious how epigenetic modulation would be used in antiviral therapy if the same epigenetic mechanisms governed viral and cellular gene expression. Recent work has highlighted several important differences between the viral and cellular chromatins, which appear to be of consequence to their respective epigenetic regulations. In this review, we will discuss the distinctiveness of the viral chromatin, and explore whether it is regulated by mechanisms unique enough to be exploited in antiviral therapy.

DNA methylation changes during long-term in vitro cell culture are caused by epigenetic drift

Culture expansion of primary cells evokes highly reproducible DNA methylation (DNAm) changes. We have identified CG dinucleotides (CpGs) that become continuously hyper- or hypomethylated during long-term culture of mesenchymal stem cells (MSCs) and other cell types. Bisulfite barcoded amplicon sequencing (BBA-seq) demonstrated that DNAm patterns of neighboring CpGs become more complex without evidence of continuous pattern development and without association to oligoclonal subpopulations. Circularized chromatin conformation capture (4C) revealed reproducible changes in nuclear organization between early and late passages, while there was no enriched interaction with other genomic regions that also harbor culture-associated DNAm changes. Chromatin immunoprecipitation of CTCF did not show significant differences during long-term culture of MSCs, however culture-associated hypermethylation was enriched at CTCF binding sites and hypomethylated CpGs were devoid of CTCF. Taken together, our results support the notion that DNAm changes during culture-expansion are not directly regulated by a targeted mechanism but rather resemble epigenetic drift.

Julia Franzen et al. investigate if changes in DNA methylation at specific genetic loci during cell culture expansion are due to a specific mechanism or gradual deregulation of an epigenetic state. Their results suggest that changes in CpG methylation are due to indirect epigenetic drift, rather than a consequence of targeting by DNA methyltransferases.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

Genome-Wide Chromatin Analysis of FFPE Tissues Using a Dual-Arm Robot with Clinical Potential

Although chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) using formalin-fixed paraffin-embedded tissue (FFPE) has been reported, it remained elusive whether they retained accurate transcription factor binding. Here, we developed a method to identify the binding sites of the insulator transcription factor CTCF and the genome-wide distribution of histone modifications involved in transcriptional activation. Importantly, we provide evidence that the ChIP-seq datasets obtained from FFPE samples are similar to or even better than the data for corresponding fresh-frozen samples, indicating that FFPE samples are compatible with ChIP-seq analysis. H3K27ac ChIP-seq analyses of 69 FFPE samples using a dual-arm robot revealed that driver mutations in EGFR were distinguishable from pan-negative cases and were relatively homogeneous as a group in lung adenocarcinomas. Thus, our results demonstrate that FFPE samples are an important source for epigenomic research, enabling the study of histone modifications, nuclear chromatin structure, and clinical data.

Spatial Organization of Chromatin: Transcriptional Control of Adaptive Immune Cell Development

Higher-order spatial organization of the genome into chromatin compartments (permissive and repressive), self-associating domains (TADs), and regulatory loops provides structural integrity and offers diverse gene regulatory controls. In particular, chromatin regulatory loops, which bring enhancer and associated transcription factors in close spatial proximity to target gene promoters, play essential roles in regulating gene expression. The establishment and maintenance of such chromatin loops are predominantly mediated involving CTCF and the cohesin machinery. In recent years, significant progress has been made in revealing how loops are assembled and how they modulate patterns of gene expression. Here we will discuss the mechanistic principles that underpin the establishment of three-dimensional (3D) chromatin structure and how changes in chromatin structure relate to alterations in gene programs that establish immune cell fate.