Isolation and analysis of rereplicated DNA by Rerep-Seq

Three-dimensional genome architecture in intrahepatic cholangiocarcinoma

PURPOSE

Intrahepatic cholangiocarcinoma (ICC) is a common primary hepatic tumors with a 5-year survival rate of less than 20%. Therefore, it is crucial to elucidate the molecular mechanisms of ICC. Recently, the advance of high-throughput chromosome conformation capture (Hi-C) technology help us look insight into the three-dimensional (3D) genome structure variation during tumorigenesis. However, its function in ICC pathogenesis remained unclear.

METHODS

Hi-C and RNA-sequencing were applied to analyze 3D genome structures and gene expression in ICC and adjacent noncancerous hepatic tissue (ANHT). Furthermore, the dysregulated genes due to 3D genome changes were validated via quantitative real-time PCR and immunohistochemistry.

RESULTS

Primarily, the intrachromosomal interactions of chr1, chr2, chr3, and chr11 and the interchromosomal interactions of chr1-chr10, chr13-chr21, chr16-chr19, and chr19-chr22 were also significantly distinct between ANHT and ICC, which may potentially contribute to the activation of cell migration and invasion via the upregulation of WNT10A, EpCAM, S100A3/A6, and MAPK12. Interestingly, 56 compartment regions from 23 chromosomes underwent A to B or B to A transitions during ICC oncogenesis, which attenuated the complement pathway through the downregulation of C8A/C8B, F7, F10, and F13B. Notably, topologically associated domain (TAD) rearrangements were identified in the region containing HOPX (chr4: 57,514,154-57,522,688) and ACVR1 (chr2:158,592,958-158,732,374) in ICC, which may contribute to the hijacking of remote enhancers that were previously outside the TAD and increased expression of HOPX and ACVR1.

CONCLUSIONS

This study reveals relationship between 3D genome structural variations and gene dysregulation during ICC tumorigenesis, indicating the molecular mechanisms and potential biomarkers.

Refining breast cancer genetic risk and biology through multi-ancestry fine-mapping analyses of 192 risk regions

Genome-wide association studies have identified approximately 200 genetic risk loci for breast cancer, but the causal variants and target genes are mostly unknown. We sought to fine-map all known breast cancer risk loci using genome-wide association study data from 172,737 female breast cancer cases and 242,009 controls of African, Asian and European ancestry. We identified 332 independent association signals for breast cancer risk, including 131 signals not reported previously, and for 50 of them, we narrowed the credible causal variants down to a single variant. Analyses integrating functional genomics data identified 195 putative susceptibility genes, enriched in PI3K/AKT, TNF/NF-κB, p53 and Wnt/β-catenin pathways. Single-cell RNA sequencing or in vitro experiment data provided additional functional evidence for 105 genes. Our study uncovered large numbers of association signals and candidate susceptibility genes for breast cancer, uncovered breast cancer genetics and biology, and supported the value of including multi-ancestry data in fine-mapping analyses.

Inferring DNA methylation in non-skeletal tissues of ancient specimens

Genome-wide premortem DNA methylation patterns can be computationally reconstructed from high-coverage DNA sequences of ancient samples. Because DNA methylation is more conserved across species than across tissues, and ancient DNA is typically extracted from bones and teeth, previous works utilizing ancient DNA methylation maps focused on studying evolutionary changes in the skeletal system. Here we suggest that DNA methylation patterns in one tissue may, under certain conditions, be informative on DNA methylation patterns in other tissues of the same individual. Using the fact that tissue-specific DNA methylation builds up during embryonic development, we identified the conditions that allow for such cross-tissue inference and devised an algorithm that carries it out. We trained the algorithm on methylation data from extant species and reached high precisions of up to 0.92 for validation datasets. We then used the algorithm on archaic humans, and identified more than 1,850 positions for which we were able to observe differential DNA methylation in prefrontal cortex neurons. These positions are linked to hundreds of genes, many of which are involved in neural functions such as structural and developmental processes. Six positions are located in the neuroblastoma breaking point family (NBPF) gene family, which probably played a role in human brain evolution. The algorithm we present here allows for the examination of epigenetic changes in tissues and cell types that are absent from the palaeontological record, and therefore provides new ways to study the evolutionary impacts of epigenetic changes.

Interpreting cis-regulatory interactions from large-scale deep neural networks

The rise of large-scale, sequence-based deep neural networks (DNNs) for predicting gene expression has introduced challenges in their evaluation and interpretation. Current evaluations align DNN predictions with orthogonal experimental data, providing insights into generalization but offering limited insights into their decision-making process. Existing model explainability tools focus mainly on motif analysis, which becomes complex when interpreting longer sequences. Here we present cis-regulatory element model explanations (CREME), an in silico perturbation toolkit that interprets the rules of gene regulation learned by a genomic DNN. Applying CREME to Enformer, a state-of-the-art DNN, we identify cis-regulatory elements that enhance or silence gene expression and characterize their complex interactions. CREME can provide interpretations across multiple scales of genomic organization, from cis-regulatory elements to fine-mapped functional sequence elements within them, offering high-resolution insights into the regulatory architecture of the genome. CREME provides a powerful toolkit for translating the predictions of genomic DNNs into mechanistic insights of gene regulation.

Crosslinking intensity modulates the reliability and sensitivity of chromatin conformation detection at different structural levels

Formaldehyde (FA) is a chemical that facilitates crosslinking between DNA and proteins. It is widely used in various biochemical assays, such as chromosome conformation capture (3C) and Chromatin Immunoprecipitation (ChIP). While the concentration and temperature of FA treatment are recognized as crucial factors in crosslinking, their quantitative effects have largely remained unexplored. In this study, we employed 3C as a model system to systematically assess the impacts of these two factors on crosslinking. Our findings indicate that the strength of crosslinking significantly influences chromatin conformation detection at nearly all known structural levels. Specifically, a delicate balance between sensitivity and reliability is required when detecting higher-level structures, such as chromosome compartments. Conversely, intense crosslinking is preferred when targeting lower-level structures, such as topologically associated domains (TADs) or chromatin loops. Based on our data, we propose a conceptual molecular thermal motion model to elucidate the roles of these two factors in restricting FA crosslinking. Our results not only shed light on the previously overlooked confounding factor in FA crosslinking but also highlight the need for caution in new technology developments that rely on FA crosslinking.

Variation of C-terminal domain governs RNA polymerase II genomic locations and alternative splicing in eukaryotic transcription

The C-terminal domain of RPB1 (CTD) orchestrates transcription by recruiting regulators to RNA Pol II upon phosphorylation. With CTD driving condensate formation on gene loci, the molecular mechanism behind how CTD-mediated recruitment of transcriptional regulators influences condensates formation remains unclear. Our study unveils that phosphorylation reversibly dissolves phase separation induced by the unphosphorylated CTD. Phosphorylated CTD, upon specific association with transcription regulators, forms distinct condensates from unphosphorylated CTD. Functional studies demonstrate CTD variants with diverse condensation properties exhibit differences in promoter binding and mRNA co-processing in cells. Notably, varying CTD lengths influence the assembly of RNA processing machinery and alternative splicing outcomes, which in turn affects cellular growth, linking the evolution of CTD variation/length with the complexity of splicing from yeast to human. These findings provide compelling evidence for a model wherein post-translational modification enables the transition of functionally specialized condensates, highlighting a co-evolution link between CTD condensation and splicing.

Loss of multi-level 3D genome organization during breast cancer progression

Breast cancer entails intricate alterations in genome organization and expression. However, how three-dimensional (3D) chromatin structure changes in the progression from a normal to a breast cancer malignant state remains unknown. To address this, we conducted an analysis combining Hi-C data with lamina-associated domains (LADs), epigenomic marks, and gene expression in an in vitro model of breast cancer progression. Our results reveal that while the fundamental properties of topologically associating domains (TADs) are overall maintained, significant changes occur in the organization of compartments and subcompartments. These changes are closely correlated with alterations in the expression of oncogenic genes. We also observe a restructuring of TAD-TAD interactions, coinciding with a loss of spatial compartmentalization and radial positioning of the 3D genome. Notably, we identify a previously unrecognized interchromosomal insertion event, wherein a locus on chromosome 8 housing the MYC oncogene is inserted into a highly active subcompartment on chromosome 10. This insertion is accompanied by the formation of de novo enhancer contacts and activation of MYC, illustrating how structural genomic variants can alter the 3D genome to drive oncogenic states. In summary, our findings provide evidence for the loss of genome organization at multiple scales during breast cancer progression revealing novel relationships between genome 3D structure and oncogenic processes.

Single-Cell Chromatin Accessibility Analysis Reveals Subgroup-Specific TF-NTR Regulatory Circuits in Medulloblastoma

Medulloblastoma (MB) stands as one of the prevalent malignant brain tumors among pediatric patients. Despite its prevalence, the intricate interplay between the regulatory program driving malignancy in MB cells and their interactions with the microenvironment remains insufficiently understood. Leveraging the capabilities of single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq), the chromatin accessibility landscape is unveiled across 59,015 distinct MB cells. This expansive dataset encompasses cells belonging to discrete molecular subgroups, namely SHH, WNT, Group3, and Group4. Within these chromatin accessibility profiles, specific regulatory elements tied to individual subgroups are uncovered, shedding light on the distinct activities of transcription factors (TFs) that likely orchestrate the tumorigenesis process. Moreover, it is found that certain neurotransmitter receptors (NTRs) are subgroup-specific and can predict MB subgroup classification when combined with their associated transcription factors. Notably, targeting essential NTRs within tumors influences both the in vitro sphere-forming capability and the in vivo tumorigenic capacity of MB cells. These findings collectively provide fresh insights into comprehending the regulatory networks and cellular dynamics within MBs. Furthermore, the significance of the TF-NTR regulatory circuits is underscored as prospective biomarkers and viable therapeutic targets.

TRIM33 plays a critical role in regulating dendritic cell differentiation and homeostasis by modulating Irf8 and Bcl2l11 transcription

The development of distinct dendritic cell (DC) subsets, namely, plasmacytoid DCs (pDCs) and conventional DC subsets (cDC1s and cDC2s), is controlled by specific transcription factors. IRF8 is essential for the fate specification of cDC1s. However, how the expression of Irf8 is regulated is not fully understood. In this study, we identified TRIM33 as a critical regulator of DC differentiation and maintenance. TRIM33 deletion in Trim33fl/fl Cre-ERT2 mice significantly impaired DC differentiation from hematopoietic progenitors at different developmental stages. TRIM33 deficiency downregulated the expression of multiple genes associated with DC differentiation in these progenitors. TRIM33 promoted the transcription of Irf8 to facilitate the differentiation of cDC1s by maintaining adequate CDK9 and Ser2 phosphorylated RNA polymerase II (S2 Pol II) levels at Irf8 gene sites. Moreover, TRIM33 prevented the apoptosis of DCs and progenitors by directly suppressing the PU.1-mediated transcription of Bcl2l11, thereby maintaining DC homeostasis. Taken together, our findings identified TRIM33 as a novel and crucial regulator of DC differentiation and maintenance through the modulation of Irf8 and Bcl2l11 expression. The finding that TRIM33 functions as a critical regulator of both DC differentiation and survival provides potential benefits for devising DC-based immune interventions and therapies.

In silico generation and augmentation of regulatory variants from massively parallel reporter assay using conditional variational autoencoder

Predicting the functional consequences of genetic variants in non-coding regions is a challenging problem. Massively parallel reporter assays (MPRAs), which are an in vitro high-throughput method, can simultaneously test thousands of variants by evaluating the existence of allele specific regulatory activity. Nevertheless, the identified labelled variants by MPRAs, which shows differential allelic regulatory effects on the gene expression are usually limited to the scale of hundreds, limiting their potential to be used as the training set for achieving a robust genome-wide prediction. To address the limitation, we propose a deep generative model, MpraVAE, to in silico generate and augment the training sample size of labelled variants. By benchmarking on several MPRA datasets, we demonstrate that MpraVAE significantly improves the prediction performance for MPRA regulatory variants compared to the baseline method, conventional data augmentation approaches as well as existing variant scoring methods. Taking autoimmune diseases as one example, we apply MpraVAE to perform a genome-wide prediction of regulatory variants and find that predicted regulatory variants are more enriched than background variants in enhancers, active histone marks, open chromatin regions in immune-related cell types, and chromatin states associated with promoter, enhancer activity and binding sites of cMyC and Pol II that regulate gene expression. Importantly, predicted regulatory variants are found to link immune-related genes by leveraging chromatin loop and accessible chromatin, demonstrating the importance of MpraVAE in genetic and gene discovery for complex traits.

Large-scale integration of omics and electronic health records to identify potential risk protein biomarkers and therapeutic drugs for cancer prevention and intervention

Identifying risk protein targets and their therapeutic drugs is crucial for effective cancer prevention. Here, we conduct integrative and fine-mapping analyses of large genome-wide association studies data for breast, colorectal, lung, ovarian, pancreatic, and prostate cancers, and characterize 710 lead variants independently associated with cancer risk. Through mapping protein quantitative trait loci (pQTL) for these variants using plasma proteomics data from over 75,000 participants, we identify 365 proteins associated with cancer risk. Subsequent colocalization analysis identifies 101 proteins, including 74 not reported in previous studies. We further characterize 36 potential druggable proteins for cancers or other disease indications. Analyzing >3.5 million electronic health records, we uncover five drugs (Haloperidol, Trazodone, Tranexamic Acid, Haloperidol, and Captopril) associated with increased cancer risk and two drugs (Caffeine and Acetazolamide) linked to reduced colorectal cancer risk. This study offers novel insights into therapeutic drugs targeting risk proteins for cancer prevention and intervention.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

Cell-type-specific and disease-associated expression quantitative trait loci in the human lung

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis. Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA sequencing of lung tissue from 66 individuals with pulmonary fibrosis and 48 unaffected donors. Using a pseudobulk approach, we mapped expression quantitative trait loci (eQTLs) across 38 cell types, observing both shared and cell-type-specific regulatory effects. Furthermore, we identified disease interaction eQTLs and demonstrated that this class of associations is more likely to be cell-type-specific and linked to cellular dysregulation in pulmonary fibrosis. Finally, we connected lung disease risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression and implicates context-specific eQTLs as key regulators of lung homeostasis and disease.

Multi-omics Integration Identifies Genes Influencing Traits Associated with Cardiovascular Risks: The Long Life Family Study

The Long Life Family Study (LLFS) enrolled 4,953 participants in 539 pedigrees displaying exceptional longevity. To identify genetic mechanisms that affect cardiovascular risks in the LLFS population, we developed a multi-omics integration pipeline and applied it to 11 traits associated with cardiovascular risks. Using our pipeline, we aggregated gene-level statistics from rare-variant analysis, GWAS, and gene expression-trait association by Correlated Meta-Analysis (CMA). Across all traits, CMA identified 64 significant genes after Bonferroni correction (p ≤ 2.8×10-7), 29 of which replicated in the Framingham Heart Study (FHS) cohort. Notably, 20 of the 29 replicated genes do not have a previously known trait-associated variant in the GWAS Catalog within 50 kb. Thirteen modules in Protein-Protein Interaction (PPI) networks are significantly enriched in genes with low meta-analysis p-values for at least one trait, three of which are replicated in the FHS cohort. The functional annotation of genes in these modules showed a significant over-representation of trait-related biological processes including sterol transport, protein-lipid complex remodeling, and immune response regulation. Among major findings, our results suggest a role of triglyceride-associated and mast-cell functional genes FCER1A, MS4A2, GATA2, HDC, and HRH4 in atherosclerosis risks. Our findings also suggest that lower expression of ATG2A, a gene we found to be associated with BMI, may be both a cause and consequence of obesity. Finally, our results suggest that ENPP3 may play an intermediary role in triglyceride-induced inflammation. Our pipeline is freely available and implemented in the Nextflow workflow language, making it easily runnable on any compute platform (https://nf-co.re/omicsgenetraitassociation).

RAD51 restricts DNA over-replication from re-activated origins

Eukaryotic cells rely on several mechanisms to ensure that the genome is duplicated precisely once in each cell division cycle, preventing DNA over-replication and genomic instability. Most of these mechanisms limit the activity of origin licensing proteins to prevent the reactivation of origins that have already been used. Here, we have investigated whether additional controls restrict the extension of re-replicated DNA in the event of origin re-activation. In a genetic screening in cells forced to re-activate origins, we found that re-replication is limited by RAD51 and enhanced by FBH1, a RAD51 antagonist. In the presence of chromatin-bound RAD51, forks stemming from re-fired origins are slowed down, leading to frequent events of fork reversal. Eventual re-initiation of DNA synthesis mediated by PRIMPOL creates ssDNA gaps that facilitate the partial elimination of re-duplicated DNA by MRE11 exonuclease. In the absence of RAD51, these controls are abrogated and re-replication forks progress much longer than in normal conditions. Our study uncovers a safeguard mechanism to protect genome stability in the event of origin reactivation.

scCircle-seq unveils the diversity and complexity of extrachromosomal circular DNAs in single cells

Extrachromosomal circular DNAs (eccDNAs) have emerged as important intra-cellular mobile genetic elements that affect gene copy number and exert in trans regulatory roles within the cell nucleus. Here, we describe scCircle-seq, a method for profiling eccDNAs and unraveling their diversity and complexity in single cells. We implement and validate scCircle-seq in normal and cancer cell lines, demonstrating that most eccDNAs vary largely between cells and are stochastically inherited during cell division, although their genomic landscape is cell type-specific and can be used to accurately cluster cells of the same origin. eccDNAs are preferentially produced from chromatin regions enriched in H3K9me3 and H3K27me3 histone marks and are induced during replication stress conditions. Concomitant sequencing of eccDNAs and RNA from the same cell uncovers the absence of correlation between eccDNA copy number and gene expression levels, except for a few oncogenes, including MYC, contained within a large eccDNA in colorectal cancer cells. Lastly, we apply scCircle-seq to one prostate cancer and two breast cancer specimens, revealing cancer-specific eccDNA landscapes and a higher propensity of eccDNAs to form in amplified genomic regions. scCircle-seq is a scalable tool that can be used to dissect the complexity of eccDNAs across different cell and tissue types, and further expands the potential of eccDNAs for cancer diagnostics.

Three-dimensional chromatin landscapes in hepatocellular carcinoma associated with hepatitis B virus

BACKGROUND

Three-dimensional (3D) chromatin architecture frequently altered in cancer. However, its changes during the pathogenesis of hepatocellular carcinoma (HCC) remained elusive.

METHODS

Hi-C and RNA-seq were applied to study the 3D chromatin landscapes and gene expression of HCC and ANHT. Hi-C Pro was used to generate genome-wide raw interaction matrices, which were normalized via iterative correction (ICE). Moreover, the chromosomes were divided into different compartments according to the first principal component (E1). Furthermore, topologically associated domains (TADs) were visualized via WashU Epigenome Browser. Furthermore, differential expression analysis of ANHT and HCC was performed using the DESeq2 R package. Additionally, dysregulated genes associated with 3D genome architecture altered were confirmed using TCGA, qRT-PCR, immunohistochemistry (IHC), etc. RESULTS: First, the intrachromosomal interactions of chr1, chr2, chr5, and chr11 were significantly different, and the interchromosomal interactions of chr4-chr10, chr13-chr21, chr15-chr22, and chr16-chr19 are remarkably different between ANHT and HCC, which resulted in the up-regulation of TP53I3 and ZNF738 and the down-regulation of APOC3 and APOA5 in HCC. Second, 49 compartment regions on 18 chromosomes have significantly switched (A-B or B-A) during HCC tumorigenesis, contributing to up-regulation of RAP2A. Finally, a tumor-specific TAD boundary located on chr5: 6271000-6478000 and enhancer hijacking were identified in HCC tissues, potentially associated with the elevated expression of MED10, whose expression were associated with poor prognosis of HCC patients.

CONCLUSION

This study demonstrates the crucial role of chromosomal structure variation in HCC oncogenesis and potential novel biomarkers of HCC, laying a foundation for cancer precision medicine development.

Immune-related transcriptomic and epigenetic reconfiguration in BV2 cells after lipopolysaccharide exposure: an in vitro omics integrative study

BACKGROUND

Molecular alterations affecting microglia have been consistently associated with the inflammatory response. These cells can have pro- or anti-inflammatory activity, phenotypes thought to be regulated by epigenetic mechanisms. Still, little is known about the details on how epigenetic marks regulate the expression of genes in the context of an inflammatory response.

METHODS

Through CUT&RUN, we profiled four genome-wide histone marks (HM) (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in lipopolysaccharide-exposed cells and compared their distributions to control cells. Transcriptomic profiles were determined through RNA-seq and differentially expressed genes were identified and contrasted with the epigenetic landscapes. Other downstream analyses were also included in this study.

RESULTS

Our results illustrate an effectively induced M1 phenotype in microglial cells derived from LPS exposure. We observed differential bound regions associated with the genes classically involved in the inflammatory response in the expected direction according to each histone modification. Consistently, our transcriptomic analysis yielded a conspicuous illustration of the LPS-induced immune activity showing the up-regulation of Nf-κB-induced mRNAs (TNF-α, nfκbiz, nfκbia) and other important genes (Marco, Il-6, etc.). Furthermore, we integrated both omics profiles and identified an important reconfiguration of the genome induced by LPS. The latter was depicted by 8 different chromatin states that changed between conditions and that associated with unique clusters of differentially expressed genes, which not only represented regulatory elements, but also underlined distinct biological functions (inhibition of morphogenesis; changes in metabolism, homeostasis, and cytokine regulation; activation of the inflammatory response).

CONCLUSION

This study exhibits important differences in the distribution of histone modifications in treated and control BV2 cells, constituting an epigenetic reconfiguration that leads to the inflammatory response. Also, it highlights the importance of these marks' regulatory role in gene expression and provides possible targets for further studies in the context of inflammation.

ChIPr: accurate prediction of cohesin-mediated 3D genome organization from 2D chromatin features

The three-dimensional genome organization influences diverse nuclear processes. Here we present Chromatin Interaction Predictor (ChIPr), a suite of regression models based on deep neural networks, random forest, and gradient boosting to predict cohesin-mediated chromatin interaction strength between any two loci in the genome. The predictions of ChIPr correlate well with ChIA-PET data in four cell lines. The standard ChIPr model requires three experimental inputs: ChIP-Seq signals for RAD21, H3K27ac, and H3K27me3 but works well with just RAD21 signal. Integrative analysis reveals novel insights into the role of CTCF motif, its orientation, and CTCF binding on cohesin-mediated chromatin interactions.

Variation of C-terminal domain governs RNA polymerase II genomic locations and alternative splicing in eukaryotic transcription

The C-terminal domain of RPB1 (CTD) orchestrates transcription by recruiting regulators to RNA Pol II upon phosphorylation. Recent insights highlight the pivotal role of CTD in driving condensate formation on gene loci. Yet, the molecular mechanism behind how CTD-mediated recruitment of transcriptional regulators influences condensates formation remains unclear. Our study unveils that phosphorylation reversibly dissolves phase separation induced by the unphosphorylated CTD. Phosphorylated CTD, upon specific association with transcription regulatory proteins, forms distinct condensates from unphosphorylated CTD. Function studies demonstrate CTD variants with diverse condensation properties in vitro exhibit difference in promoter binding and mRNA co-processing in cells. Notably, varying CTD lengths lead to alternative splicing outcomes impacting cellular growth, linking the evolution of CTD variation/length with the complexity of splicing from yeast to human. These findings provide compelling evidence for a model wherein post-translational modification enables the transition of functionally specialized condensates, highlighting a co-evolution link between CTD condensation and splicing.

Uniform quantification of single-nucleus ATAC-seq data with Paired-Insertion Counting (PIC) and a model-based insertion rate estimator

Existing approaches to scoring single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) feature matrices from sequencing reads are inconsistent, affecting downstream analyses and displaying artifacts. We show that, even with sparse single-cell data, quantitative counts are informative for estimating the regulatory state of a cell, which calls for a consistent treatment. We propose Paired-Insertion Counting as a uniform method for snATAC-seq feature characterization and provide a probability model for inferring latent insertion dynamics from snATAC-seq count matrices.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Integrative analysis of DNA replication origins and ORC/MCM binding sites in human cells reveals a lack of overlap

Based on experimentally determined average inter-origin distances of ∼100 kb, DNA replication initiates from ∼50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the Origin Recognition Complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and 5 ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ∼7.5 million union origins identified by all datasets, only 0.27% were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques (20,250 shared origins), suggesting extensive variability in origin usage and identification. 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF binding sites, G-quadruplex sites and activating histone marks, these overlaps are comparable or less than that of known Transcription Start Sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ∼13,000 reproducible ORC binding sites in human cancer cells, and only 4.5% were within 1 kb of the ∼11,000 union MCM2-7 binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, S. cerevisiae . Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

Analysis of the chromatin landscape and RNA polymerase II binding at SIN3-regulated genes

The chromatin environment has a significant impact on gene expression. Chromatin structure is highly regulated by histone modifications and RNA polymerase II binding dynamics. The SIN3 histone modifying complex regulates the chromatin environment leading to changes in gene expression. In Drosophila melanogaster, the Sin3A gene is alternatively spliced to produce different protein isoforms, two of which include SIN3 220 and SIN3 187. Both SIN3 isoforms are scaffolding proteins that interact with several other factors to regulate the chromatin landscape. The mechanism through which the SIN3 isoforms regulate chromatin is not well understood. Here, we analyze publicly available data sets to allow us to ask specific questions on how SIN3 isoforms regulate chromatin and gene activity. We determined that genes repressed by the SIN3 isoforms exhibited enrichment in histone H3K4me2, H3K4me3, H3K14ac and H3K27ac near the transcription start site. We observed an increase in the amount of paused RNA polymerase II on the promoter of genes repressed by the isoforms as compared to genes that require SIN3 for maximum activation. Furthermore, we analyzed a subset of genes regulated by SIN3 187 that suggest a mechanism in which SIN3 187 might exhibit hard regulation as well as soft regulation. Data presented here expand our knowledge of how the SIN3 isoforms regulate the chromatin environment and RNA polymerase II binding dynamics.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in several cancer types including sarcomas, which are mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors (Pparγ and Cebpα) and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks at putative enhancer elements and promoters. Finally, we observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Our results demonstrate that ATRX functions to buffer against differentiation in mesenchymal progenitor cells, which has implications for understanding ATRX loss of function in sarcomas.

Exploring molecular markers and drug candidates for colorectal cancer through comprehensive bioinformatics analysis

Colorectal cancer (CRC) often has a poor prognosis and identifying useful and novel agents for treating CRC is urgently required. This study aimed to examine molecular markers associated with CRC prognosis and to identify potential drug candidates. The differentially expressed genes (DEGs) of CRC in TCGA were identified. The genes associated with CRC, summarized from NCBI-gene, OMIM, and the DEGs, were used to construct a co-expression network by WGCNA. Moreover, the co-expression genes from modules of interest were used to carry out functional enrichment. A total of 2742 DEGs, including 1674 upregulated and 1068 downregulated genes, were identified. Thirteen co-expression modules were constructed with WGCNA. Brown and blue co-expression modules with significant differences in disease phenotype were found. Functional enrichment analysis showed that genes in the brown module were mainly related to cell cycle, cell proliferation, DNA replication, and RNA transport. The genes in the blue module were mainly associated with fatty acid degradation, sulfur metabolism, PPAR signaling pathway and bile secretion. In addition, both the genes in brown and blue were associated with tumor staging. Some prognostic markers and candidate small molecules drugs for CRC treatment were identified. In conclusion, we revealed molecular biomarker profiles in CRC by systematic bioinformatics analysis, constructed regulatory networks of mRNA, ncRNA and transcriptional regulators (TFs), and identified potential drugs targeting hub proteins and TFs.

A comprehensive revisit of the machine-learning tools developed for the identification of enhancers in the human genome

Enhancers are non-coding DNA elements that play a crucial role in enhancing the transcription rate of a specific gene in the genome. Experiments for identifying enhancers can be restricted by their conditions and involve complicated, time-consuming, laborious, and costly steps. To overcome these challenges, computational platforms have been developed to complement experimental methods that enable high-throughput identification of enhancers. Over the last few years, the development of various enhancer computational tools has resulted in significant progress in predicting putative enhancers. Thus, researchers are now able to use a variety of strategies to enhance and advance enhancer study. In this review, an overview of machine learning (ML)-based prediction methods for enhancer identification and related databases has been provided. The existing enhancer-prediction methods have also been reviewed regarding their algorithms, feature selection processes, validation techniques, and software utility. In addition, the advantages and drawbacks of these ML approaches and guidelines for developing bioinformatic tools have been highlighted for a more efficient enhancer prediction. This review will serve as a useful resource for experimentalists in selecting the appropriate ML tool for their study, and for bioinformaticians in developing more accurate and advanced ML-based predictors.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

Transition from Transient DNA Rereplication to Inherited Gene Amplification Following Prolonged Environmental Stress

Cells require the ability to adapt to changing environmental conditions, however, it is unclear how these changes elicit stable permanent changes in genomes. We demonstrate that, in response to environmental metal exposure, the metallothionein (MT) locus undergoes DNA rereplication generating transient site-specific gene amplifications (TSSGs). Chronic metal exposure allows transition from MT TSSG to inherited MT gene amplification through homologous recombination within and outside of the MT locus. DNA rereplication of the MT locus is suppressed by H3K27me3 and EZH2. Long-term ablation of EZH2 activity eventually leads to integration and inheritance of MT gene amplifications without the selective pressure of metal exposure. The rereplication and inheritance of MT gene amplification is an evolutionarily conserved response to environmental metal from yeast to human. Our results describe a new paradigm for adaptation to environmental stress where targeted, transient DNA rereplication precedes stable inherited gene amplification.

DNA Methylation Dynamics During Esophageal Epithelial Regeneration Following Repair with Acellular Silk Fibroin Grafts in Rat

Esophageal pathologies such as atresia and benign strictures often require surgical reconstruction with autologous tissues to restore organ continuity. Complications such as donor site morbidity and limited tissue availability have spurred the development of acellular grafts for esophageal tissue replacement. Acellular biomaterials for esophageal repair rely on the activation of intrinsic regenerative mechanisms to mediate de novo tissue formation at implantation sites. Previous research has identified signaling cascades involved in neoepithelial formation in a rat model of onlay esophagoplasty with acellular silk fibroin grafts, including phosphoinositide 3-kinase (PI3K), and protein kinase B (Akt) signaling. However, it is currently unknown how these mechanisms are governed by DNA methylation (DNAme) during esophageal wound healing processes. Reduced-representation bisulfite sequencing is performed to characterize temporal DNAme dynamics in host and regenerated tissues up to 1 week postimplantation. Overall, global hypermethylation is observed at postreconstruction timepoints and an inverse correlation between promoter DNAme and the expression levels of differentially expressed proteins during regeneration. Site-specific hypomethylation targets genes associated with immune activation, while hypermethylation occurs within gene bodies encoding PI3K-Akt signaling components during the tissue remodeling period. The data provide insight into the epigenetic mechanisms during esophageal regeneration following surgical repair with acellular grafts.

PEACOCK: a machine learning approach to assess the validity of cell type-specific enhancer-gene regulatory relationships

The vast majority of disease-associated variants identified in genome-wide association studies map to enhancers, powerful regulatory elements which orchestrate the recruitment of transcriptional complexes to their target genes' promoters to upregulate transcription in a cell type- and timing-dependent manner. These variants have implicated thousands of enhancers in many common genetic diseases, including nearly all cancers. However, the etiology of most of these diseases remains unknown because the regulatory target genes of the vast majority of enhancers are unknown. Thus, identifying the target genes of as many enhancers as possible is crucial for learning how enhancer regulatory activities function and contribute to disease. Based on experimental results curated from scientific publications coupled with machine learning methods, we developed a cell type-specific score predictive of an enhancer targeting a gene. We computed the score genome-wide for every possible cis enhancer-gene pair and validated its predictive ability in four widely used cell lines. Using a pooled final model trained across multiple cell types, all possible gene-enhancer regulatory links in cis (~17 M) were scored and added to the publicly available PEREGRINE database ( www.peregrineproj.org ). These scores provide a quantitative framework for the enhancer-gene regulatory prediction that can be incorporated into downstream statistical analyses.

MYT1L is required for suppressing earlier neuronal development programs in the adult mouse brain

In vitro studies indicate the neurodevelopmental disorder gene myelin transcription factor 1-like (MYT1L) suppresses non-neuronal lineage genes during fibroblast-to-neuron direct differentiation. However, MYT1L's molecular and cellular functions in the adult mammalian brain have not been fully characterized. Here, we found that MYT1L loss leads to up-regulated deep layer (DL) gene expression, corresponding to an increased ratio of DL/UL neurons in the adult mouse cortex. To define potential mechanisms, we conducted Cleavage Under Targets & Release Using Nuclease (CUT&RUN) to map MYT1L binding targets and epigenetic changes following MYT1L loss in mouse developing cortex and adult prefrontal cortex (PFC). We found MYT1L mainly binds to open chromatin, but with different transcription factor co-occupancies between promoters and enhancers. Likewise, multiomic data set integration revealed that, at promoters, MYT1L loss does not change chromatin accessibility but increases H3K4me3 and H3K27ac, activating both a subset of earlier neuronal development genes as well as Bcl11b, a key regulator for DL neuron development. Meanwhile, we discovered that MYT1L normally represses the activity of neurogenic enhancers associated with neuronal migration and neuronal projection development by closing chromatin structures and promoting removal of active histone marks. Further, we showed that MYT1L interacts with HDAC2 and transcriptional repressor SIN3B in vivo, providing potential mechanisms underlying repressive effects on histone acetylation and gene expression. Overall, our findings provide a comprehensive map of MYT1L binding in vivo and mechanistic insights into how MYT1L loss leads to aberrant activation of earlier neuronal development programs in the adult mouse brain.

Fetal hypoxia results in sex- and cell type-specific alterations in neonatal transcription in rat oligodendrocyte precursor cells, microglia, neurons, and oligodendrocytes

BACKGROUND

Fetal hypoxia causes vital, systemic, developmental malformations in the fetus, particularly in the brain, and increases the risk of diseases in later life. We previously demonstrated that fetal hypoxia exposure increases the susceptibility of the neonatal brain to hypoxic-ischemic insult. Herein, we investigate the effect of fetal hypoxia on programming of cell-specific transcriptomes in the brain of neonatal rats.

RESULTS

We obtained RNA sequencing (RNA-seq) data from neurons, microglia, oligodendrocytes, A2B5⁺ oligodendrocyte precursor cells, and astrocytes from male and female neonatal rats subjected either to fetal hypoxia or control conditions. Substantial transcriptomic responses to fetal hypoxia occurred in neurons, microglia, oligodendrocytes, and A2B5⁺ cells. Not only were the transcriptomic responses unique to each cell type, but they also occurred with a great deal of sexual dimorphism. We validated differential expression of several genes related to inflammation and cell death by Real-time Quantitative Polymerase Chain Reaction (qRT-PCR). Pathway and transcription factor motif analyses suggested that the NF-kappa B (NFκB) signaling pathway was enriched in the neonatal male brain due to fetal hypoxia, and we verified this result by transcription factor assay of NFκB-p65 in whole brain.

CONCLUSIONS

Our study reveals a significant impact of fetal hypoxia on the transcriptomes of neonatal brains in a cell-specific and sex-dependent manner, and provides mechanistic insights that may help explain the development of hypoxic-ischemic sensitive phenotypes in the neonatal brain.

Nanoparticle-mediated targeting of the fusion gene RUNX1/ETO in t(8;21)-positive acute myeloid leukaemia

A hallmark of acute myeloid leukaemias (AMLs) are chromosomal rearrangements that give rise to novel leukaemia-specific fusion genes. Most of these fusion genes are both initiating and driving events in AML and therefore constitute ideal therapeutic targets but are challenging to target by conventional drug development. siRNAs are frequently used for the specific suppression of fusion gene expression but require special formulations for efficient in vivo delivery. Here we describe the use of siRNA-loaded lipid nanoparticles for the specific therapeutic targeting of the leukaemic fusion gene RUNX1/ETO. Transient knockdown of RUNX1/ETO reduces its binding to its target genes and alters the binding of RUNX1 and its co-factor CBFβ. Transcriptomic changes in vivo were associated with substantially increased median survival of a t(8;21)-AML mouse model. Importantly, transient knockdown in vivo causes long-lasting inhibition of leukaemic proliferation and clonogenicity, induction of myeloid differentiation and a markedly impaired re-engraftment potential in vivo. These data strongly suggest that temporary inhibition of RUNX1/ETO results in long-term restriction of leukaemic self-renewal. Our results provide proof for the feasibility of targeting RUNX1/ETO in a pre-clinical setting and support the further development of siRNA-LNPs for the treatment of fusion gene-driven malignancies.

3D genome alterations in T cells associated with disease activity of systemic lupus erythematosus

OBJECTIVES

Three-dimensional (3D) genome alterations can dysregulate gene expression by rewiring physical interactions within chromosomes in a tissue-specific or cell-specific manner and lead to diseases. We aimed to elucidate the 3D genome structure and its role in gene expression networks dysregulated in systemic lupus erythematosus (SLE).

METHODS

We performed Hi-C experiments using CD4+ T cells from 7 patients with SLE and 5 age-matched and sex-matched healthy controls (HCs) combined with RNA sequencing analysis. Further integrative analyses, including transcription factor motif enrichment, SPI1 knockdown and histone modifications (H3K27ac, H3K4me1, H3K4me3), were performed for altered loop-associated gene loci in SLE.

RESULTS

We deciphered the 3D chromosome organisation in T cells of patients with SLE and found it was clearly distinct from that of HCs and closely associated with the disease activity of SLE. Importantly, we identified loops within chromosomes associated with the disease activity of SLE and differentially expressed genes and found some key histone modifications close to these loops. Moreover, we demonstrated the contribution of the transcription factor SPI1, whose motif is located in the altered loop in SLE, to the overexpression of interferon pathway gene. In addition, we identified the potential influences of genetic variations in 3D genome alterations in SLE.

CONCLUSIONS

Our results highlight the 3D genome structure alterations associated with SLE development and provide a foundation for future interrogation of the relationships between chromosome structure and gene expression control in SLE.

Hnf1b renal expression directed by a distal enhancer responsive to Pax8

Xenopus provides a simple and efficient model system to study nephrogenesis and explore the mechanisms causing renal developmental defects in human. Hnf1b (hepatocyte nuclear factor 1 homeobox b), a gene whose mutations are the most commonly identified genetic cause of developmental kidney disease, is required for the acquisition of a proximo-intermediate nephron segment in Xenopus as well as in mouse. Genetic networks involved in Hnf1b expression during kidney development remain poorly understood. We decided to explore the transcriptional regulation of Hnf1b in the developing Xenopus pronephros and mammalian renal cells. Using phylogenetic footprinting, we identified an evolutionary conserved sequence (CNS1) located several kilobases (kb) upstream the Hnf1b transcription start and harboring epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. By means of functional expression assays in Xenopus and mammalian renal cell lines we showed that CNS1 displays enhancer activity in renal tissue. Using CRISPR/cas9 editing in Xenopus tropicalis, we demonstrated the in vivo functional relevance of CNS1 in driving hnf1b expression in the pronephros. We further showed the importance of Pax8-CNS1 interaction for CNS1 enhancer activity allowing us to conclude that Hnf1b is a direct target of Pax8. Our work identified for the first time a Hnf1b renal specific enhancer and may open important perspectives into the diagnosis for congenital kidney anomalies in human, as well as modeling HNF1B-related diseases.

Deciphering multi-way interactions in the human genome

Chromatin architecture, a key regulator of gene expression, can be inferred using chromatin contact data from chromosome conformation capture, or Hi-C. However, classical Hi-C does not preserve multi-way contacts. Here we use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organization in the human genome. We use hypergraph theory for data representation and analysis, and quantify higher order structures in neonatal fibroblasts, biopsied adult fibroblasts, and B lymphocytes. By integrating multi-way contacts with chromatin accessibility, gene expression, and transcription factor binding, we introduce a data-driven method to identify cell type-specific transcription clusters. We provide transcription factor-mediated functional building blocks for cell identity that serve as a global signature for cell types.

Mapping higher order chromatin architecture is important. Here the authors use long sequencing reads to map genome-wide multi-way contacts and investigate higher order chromatin organisation; they use hypergraph theory for data representation and analysis, and apply this to different cell types.

Transcription-coupled and epigenome-encoded mechanisms direct H3K4 methylation

Mono-, di-, and trimethylation of histone H3 lysine 4 (H3K4me1/2/3) are associated with transcription, yet it remains controversial whether H3K4me1/2/3 promote or result from transcription. Our previous characterizations of Arabidopsis H3K4 demethylases suggest roles for H3K4me1 in transcription. However, the control of H3K4me1 remains unexplored in Arabidopsis, in which no methyltransferase for H3K4me1 has been identified. Here, we identify three Arabidopsis methyltransferases that direct H3K4me1. Analyses of their genome-wide localization using ChIP-seq and machine learning reveal that one of the enzymes cooperates with the transcription machinery, while the other two are associated with specific histone modifications and DNA sequences. Importantly, these two types of localization patterns are also found for the other H3K4 methyltransferases in Arabidopsis and mice. These results suggest that H3K4me1/2/3 are established and maintained via interplay with transcription as well as inputs from other chromatin features, presumably enabling elaborate gene control.

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells

Sjögren's disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren's cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.

Targeted erasure of DNA methylation by TET3 drives adipogenic reprogramming and differentiation

DNA methylation is a crucial epigenetic modification in the establishment of cell-type-specific characteristics. However, how DNA methylation is selectively reprogrammed at adipocyte-specific loci during adipogenesis remains unclear. Here, we show that the transcription factor, C/EBPδ, and the DNA methylation eraser, TET3, cooperatively control adipocyte differentiation. We perform whole-genome bisulfite sequencing to explore the dynamics and regulatory mechanisms of DNA methylation in adipocyte differentiation. During adipogenesis, DNA methylation selectively decreases at adipocyte-specific loci carrying the C/EBP binding motif, which correlates with the activity of adipogenic promoters and enhancers. Mechanistically, we find that C/EBPδ recruits a DNA methylation eraser, TET3, to catalyse DNA demethylation at the C/EBP binding motif and stimulate the expression of key adipogenic genes. Ectopic expression of TET3 potentiates in vitro and in vivo adipocyte differentiation and recovers downregulated adipogenic potential, which is observed in aged mice and humans. Taken together, our study highlights how targeted reprogramming of DNA methylation through cooperative action of the transcription factor C/EBPδ, and the DNA methylation eraser TET3, controls adipocyte differentiation.

TLR4 is a regulator of trained immunity in a murine model of Duchenne muscular dystrophy

Dysregulation of the balance between pro-inflammatory and anti-inflammatory macrophages has a key function in the pathogenesis of Duchenne muscular dystrophy (DMD), a fatal genetic disease. We postulate that an evolutionarily ancient protective mechanism against infection, known as trained immunity, drives pathological inflammation in DMD. Here we show that bone marrow-derived macrophages from a murine model of DMD (mdx) exhibit cardinal features of trained immunity, consisting of transcriptional hyperresponsiveness associated with metabolic and epigenetic remodeling. The hyperresponsive phenotype is transmissible by bone marrow transplantation to previously healthy mice and persists for up to 11 weeks post-transplant. Mechanistically, training is induced by muscle extract in vitro. The functional and epigenetic changes in bone marrow-derived macrophages from dystrophic mice are TLR4-dependent. Adoptive transfer experiments further support the TLR4-dependence of trained macrophages homing to damaged muscles from the bone marrow. Collectively, this suggests that a TLR4-regulated, memory-like capacity of innate immunity induced at the level of the bone marrow promotes dysregulated inflammation in DMD.

Preventing excess replication origin activation to ensure genome stability

Cells activate distinctive regulatory pathways that prevent excessive initiation of DNA replication to achieve timely and accurate genome duplication. Excess DNA synthesis is constrained by protein-DNA interactions that inhibit initiation at dormant origins. In parallel, specific modifications of pre-replication complexes prohibit post-replicative origin relicensing. Replication stress ensues when the controls that prevent excess replication are missing in cancer cells, which often harbor extrachromosomal DNA that can be further amplified by recombination-mediated processes to generate chromosomal translocations. The genomic instability that accompanies excess replication origin activation can provide a promising target for therapeutic intervention. Here we review molecular pathways that modulate replication origin dormancy, prevent excess origin activation, and detect, encapsulate, and eliminate persistent excess DNA.

Spatiotemporal dynamics of SETD5-containing NCoR-HDAC3 complex determines enhancer activation for adipogenesis

Enhancer activation is essential for cell-type specific gene expression during cellular differentiation, however, how enhancers transition from a hypoacetylated "primed" state to a hyperacetylated-active state is incompletely understood. Here, we show SET domain-containing 5 (SETD5) forms a complex with NCoR-HDAC3 co-repressor that prevents histone acetylation of enhancers for two master adipogenic regulatory genes Cebpa and Pparg early during adipogenesis. The loss of SETD5 from the complex is followed by enhancer hyperacetylation. SETD5 protein levels were transiently increased and rapidly degraded prior to enhancer activation providing a mechanism for the loss of SETD5 during the transition. We show that induction of the CDC20 co-activator of the ubiquitin ligase leads to APC/C mediated degradation of SETD5 during the transition and this operates as a molecular switch that facilitates adipogenesis.

Comprehensive determination of transcription start sites derived from all RNA polymerases using ReCappable-seq

Determination of eukaryotic transcription start sites (TSSs) has been based on methods that require the cap structure at the 5' end of transcripts derived from Pol II RNA polymerase. Consequently, these methods do not reveal TSSs derived from the other RNA polymerases that also play critical roles in various cell functions. To address this limitation, we developed ReCappable-seq, which comprehensively identifies TSS for both Pol II and non-Pol II transcripts at single-nucleotide resolution. The method relies on specific enzymatic exchange of 5' m⁷G caps and 5' triphosphates with a selectable tag. When applied to human transcriptomes, ReCappable-seq identifies Pol II TSSs that are in agreement with orthogonal methods such as CAGE. Additionally, ReCappable-seq reveals a rich landscape of TSSs associated with Pol III transcripts that have not previously been amenable to study at genome-wide scale. Novel TSS from non-Pol II transcription can be located in the nuclear and mitochondrial genomes. ReCappable-seq interrogates the regulatory landscape of coding and noncoding RNA concurrently and enables the classification of epigenetic profiles associated with Pol II and non-Pol II TSS.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

BACKGROUND

Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging.

METHODS

Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain.

RESULTS

Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism.

CONCLUSION

This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

The molecular basis, genetic control and pleiotropic effects of local gene co-expression

Nearby genes are often expressed as a group. Yet, the prevalence, molecular mechanisms and genetic control of local gene co-expression are far from being understood. Here, by leveraging gene expression measurements across 49 human tissues and hundreds of individuals, we find that local gene co-expression occurs in 13% to 53% of genes per tissue. By integrating various molecular assays (e.g. ChIP-seq and Hi-C), we estimate the ability of several mechanisms, such as enhancer-gene interactions, in distinguishing gene pairs that are co-expressed from those that are not. Notably, we identify 32,636 expression quantitative trait loci (eQTLs) which associate with co-expressed gene pairs and often overlap enhancer regions. Due to affecting several genes, these eQTLs are more often associated with multiple human traits than other eQTLs. Our study paves the way to comprehend trait pleiotropy and functional interpretation of QTL and GWAS findings. All local gene co-expression identified here is available through a public database ( https://glcoex.unil.ch/ ).