BAMscale: quantification of next-generation sequencing peaks and generation of scaled coverage tracks

YY1 binding is a gene-intrinsic barrier to Xist-mediated gene silencing

X chromosome inactivation (XCI) in mammals is mediated by Xist RNA which functions in cis to silence genes on a single X chromosome in XX female cells, thereby equalising levels of X-linked gene expression relative to XY males. XCI progresses over a period of several days, with some X-linked genes silencing faster than others. The chromosomal location of a gene is an important determinant of silencing rate, but uncharacterised gene-intrinsic features also mediate resistance or susceptibility to silencing. In this study, we examine mouse embryonic stem cell lines with an inducible Xist allele (iXist-ChrX mESCs) and integrate allele-specific data of gene silencing and decreasing inactive X (Xi) chromatin accessibility over time courses of Xist induction with cellular differentiation. Our analysis reveals that motifs bound by the transcription factor YY1 are associated with persistently accessible regulatory elements, including many promoters and enhancers of slow-silencing genes. We further show that YY1 is evicted relatively slowly from target sites on Xi, and that silencing of X-linked genes is increased upon YY1 degradation. Together our results suggest that YY1 acts as a barrier to Xist-mediated silencing until the late stages of the XCI process.

Transcription elongation defects link oncogenic SF3B1 mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human diseases remains unexplored. Using isogenic cell lines, patient samples, and a mutant mouse model, we investigated how cancer-associated mutations in SF3B1 alter transcription. We found that these mutations reduce the elongation rate of RNA polymerase II (RNAPII) along gene bodies and its density at promoters. The elongation defect results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1. The decreased promoter-proximal RNAPII density reduces both chromatin accessibility and H3K4me3 marks at promoters. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC/H3K4me pathway, which, when modulated, reverse both transcription and chromatin changes. Our findings reveal how splicing factor mutant states behave functionally as epigenetic disorders through impaired transcription-related changes to the chromatin landscape. We also present a rationale for targeting the Sin3/HDAC complex as a therapeutic strategy.

Aberrant non-canonical NF-κB signalling reprograms the epigenome landscape to drive oncogenic transcriptomes in multiple myeloma

In multiple myeloma, abnormal plasma cells establish oncogenic niches within the bone marrow by engaging the NF-κB pathway to nurture their survival while they accumulate pro-proliferative mutations. Under these conditions, many cases eventually develop genetic abnormalities endowing them with constitutive NF-κB activation. Here, we find that sustained NF-κB/p52 levels resulting from such mutations favours the recruitment of enhancers beyond the normal B-cell repertoire. Furthermore, through targeted disruption of p52, we characterise how such enhancers are complicit in the formation of super-enhancers and the establishment of cis-regulatory interactions with myeloma dependencies during constitutive activation of p52. Finally, we functionally validate the pathological impact of these cis-regulatory modules on cell and tumour phenotypes using in vitro and in vivo models, confirming RGS1 as a p52-dependent myeloma driver. We conclude that the divergent epigenomic reprogramming enforced by aberrant non-canonical NF-κB signalling potentiates transcriptional programs beneficial for multiple myeloma progression.

PRC1 directs PRC2-H3K27me3 deposition to shield adult spermatogonial stem cells from differentiation

Spermatogonial stem cells functionality reside in the slow-cycling and heterogeneous undifferentiated spermatogonia cell population. This pool of cells supports lifelong fertility in adult males by balancing self-renewal and differentiation to produce haploid gametes. However, the molecular mechanisms underpinning long-term stemness of undifferentiated spermatogonia during adulthood remain unclear. Here, we discover that an epigenetic regulator, Polycomb repressive complex 1 (PRC1), shields adult undifferentiated spermatogonia from differentiation, maintains slow cycling, and directs commitment to differentiation during steady-state spermatogenesis in adults. We show that PRC2-mediated H3K27me3 is an epigenetic hallmark of adult undifferentiated spermatogonia. Indeed, spermatogonial differentiation is accompanied by a global loss of H3K27me3. Disruption of PRC1 impairs global H3K27me3 deposition, leading to precocious spermatogonial differentiation. Therefore, PRC1 directs PRC2-H3K27me3 deposition to maintain the self-renewing state of undifferentiated spermatogonia. Importantly, in contrast to its role in other tissue stem cells, PRC1 negatively regulates the cell cycle to maintain slow cycling of undifferentiated spermatogonia. Our findings have implications for how epigenetic regulators can be tuned to regulate the stem cell potential, cell cycle and differentiation to ensure lifelong fertility in adult males.

Restrictor synergizes with Symplekin and PNUTS to terminate extragenic transcription

Transcription termination pathways mitigate the detrimental consequences of unscheduled promiscuous initiation occurring at hundreds of thousands of genomic cis-regulatory elements. The Restrictor complex, composed of the Pol II-interacting protein WDR82 and the RNA-binding protein ZC3H4, suppresses processive transcription at thousands of extragenic sites in mammalian genomes. Restrictor-driven termination does not involve nascent RNA cleavage, and its interplay with other termination machineries is unclear. Here we show that efficient termination at Restrictor-controlled extragenic transcription units involves the recruitment of the protein phosphatase 1 (PP1) regulatory subunit PNUTS, a negative regulator of the SPT5 elongation factor, and Symplekin, a protein associated with RNA cleavage complexes but also involved in cleavage-independent and phosphatase-dependent termination of noncoding RNAs in yeast. PNUTS and Symplekin act synergistically with, but independently from, Restrictor to dampen processive extragenic transcription. Moreover, the presence of limiting nuclear levels of Symplekin imposes a competition for its recruitment among multiple transcription termination machineries, resulting in mutual regulatory interactions. Hence, by synergizing with Restrictor, Symplekin and PNUTS enable efficient termination of processive, long-range extragenic transcription.

SetDB1 and Su(var)3-9 are essential for late stages of larval development of Drosophila melanogaster

Methylation of H3K9 histone residue is a marker of gene silencing in eukaryotes. Three enzymes responsible for adding this modification - G9a, SetDB1/Egg, and Su(var)3-9 - are known in Drosophila. To understand how simultaneous mutations of SetDB1 and Su(var)3-9 may affect the fly development, appropriate combinations were obtained. Double mutants egg; Su(var)3-9 displayed pronounced embryonic lethality, slower larval growth and died before or during metamorphosis. Analysis of transcription in larval salivary glands and wing imaginal disks indicated that the effect of double mutation is tissue-specific. In salivary gland chromosomes, affected genes display low H3K9me2 enrichment and are rarely bound by SetDB1 or Su(var)3-9. We suppose that each of these enzymes directly or indirectly controls its own set of gene targets in different organs, and double mutation results in an imbalanced developmental program. This also indicates that SetDB1 and Su(var)3-9 may affect transcription via H3K9-independent mechanisms. Unexpectedly, in double and triple mutants, amount of di- and tri-methylated H3K9 is drastically reduced, but not completely absent. We hypothesize that this residual methylation implies the existence of additional H3K9-specific methyltransferase in Drosophila.

Phylogenetic modeling of enhancer shifts in African mole-rats reveals regulatory changes associated with tissue-specific traits

Changes in gene regulation are thought to underlie most phenotypic differences between species. For subterranean rodents such as the naked mole-rat, proposed phenotypic adaptations include hypoxia tolerance, metabolic changes, and cancer resistance. However, it is largely unknown what regulatory changes may associate with these phenotypic traits, and whether these are unique to the naked mole-rat, the mole-rat clade, or are also present in other mammals. Here, we investigate regulatory evolution in the heart and liver from two African mole-rat species and two rodent outgroups using genome-wide epigenomic profiling. First, we adapted and applied a phylogenetic modeling approach to quantitatively compare epigenomic signals at orthologous regulatory elements and identified thousands of promoter and enhancer regions with differential epigenomic activity in mole-rats. These elements associate with known mole-rat adaptations in metabolic and functional pathways and suggest candidate genetic loci that may underlie mole-rat innovations. Second, we evaluated ancestral and species-specific regulatory changes in the study phylogeny and report several candidate pathways experiencing stepwise remodeling during the evolution of mole-rats, such as the insulin and hypoxia response pathways. Third, we report nonorthologous regulatory elements overlap with lineage-specific repetitive elements and appear to modify metabolic pathways by rewiring of HNF4 and RAR/RXR transcription factor binding sites in mole-rats. These comparative analyses reveal how mole-rat regulatory evolution informs previously reported phenotypic adaptations. Moreover, the phylogenetic modeling framework we propose here improves upon the state of the art by addressing known limitations of inter-species comparisons of epigenomic profiles and has broad implications in the field of comparative functional genomics.

Small cytosolic double-stranded DNA represses cyclic GMP-AMP synthase activation and induces autophagy

The cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway is a major mediator of inflammation following stimulation with >45 bp double-stranded DNA (dsDNA). Herein, we identify a class of ∼20-40 bp small cytosolic dsDNA (scDNA) molecules that compete with long dsDNA (200-1,500 bp herring testis [HT]-DNA) for binding to cGAS, thus repressing HT-DNA-induced cGAS activation. The scDNA promotes cGAS and Beclin-1 interaction, releasing Rubicon, a negative regulator of phosphatidylinositol 3-kinase class III (PI3KC3), from the Beclin-1-PI3KC3 complex. This leads to PI3KC3 activation and induces autophagy, causing degradation of STING and long cytosolic dsDNA. Moreover, DNA damage decreases, and autophagy inducers increase scDNA levels. scDNA transfection and treatment with autophagy inducers attenuate DNA damage-induced cGAS activation. Thus, scDNA molecules serve as effective brakes for cGAS activation, preventing excessive inflammatory cytokine production following DNA damage. Our findings may have therapeutic implications for cytosolic DNA-associated inflammatory diseases.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Extrachromosomal DNA Amplification Contributes to Small Cell Lung Cancer Heterogeneity and Is Associated with Worse Outcomes

Small-cell lung cancer (SCLC) is an aggressive neuroendocrine lung cancer. Oncogenic MYC amplifications drive SCLC heterogeneity, but the genetic mechanisms of MYC amplification and phenotypic plasticity, characterized by neuroendocrine and nonneuroendocrine cell states, are not known. Here, we integrate whole-genome sequencing, long-range optical mapping, single-cell DNA sequencing, and fluorescence in situ hybridization to find extrachromosomal DNA (ecDNA) as a primary source of SCLC oncogene amplifications and driver fusions. ecDNAs bring to proximity enhancer elements and oncogenes, creating SCLC transcription-amplifying units, driving exceptionally high MYC gene dosage. We demonstrate that cell-free nucleosome profiling can noninvasively detect ecDNA amplifications in plasma, facilitating its genome-wide interrogation in SCLC and other cancers. Altogether, our work provides the first comprehensive map of SCLC ecDNA and describes a new mechanism that governs MYC-driven SCLC heterogeneity. ecDNA-enabled transcriptional flexibility may explain the significantly worse survival outcomes of SCLC harboring complex ecDNA amplifications.

SIGNIFICANCE

MYC drives SCLC progression, but the genetic basis of MYC-driven SCLC evolution is unknown. Using SCLC as a paradigm, we report how ecDNA amplifications function as MYC-amplifying units, fostering tumor plasticity and a high degree of tumor heterogeneity. This article is highlighted in the In This Issue feature, p. 799.

Transcription elongation defects link oncogenic splicing factor mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human disease remains unexplored. Here, we investigated the impact of non-synonymous mutations in SF3B1 and U2AF1, two commonly mutated splicing factors in cancer, on transcription. We find that the mutations impair RNA Polymerase II (RNAPII) transcription elongation along gene bodies leading to transcription-replication conflicts, replication stress and altered chromatin organization. This elongation defect is linked to disrupted pre-spliceosome assembly due to impaired association of HTATSF1 with mutant SF3B1. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC complex, which, when modulated, normalize transcription defects and their downstream effects. Our findings shed light on the mechanisms by which oncogenic mutant spliceosomes impact chromatin organization through their effects on RNAPII transcription elongation and present a rationale for targeting the Sin3/HDAC complex as a potential therapeutic strategy.

GRAPHICAL ABSTRACT

HIGHLIGHTS

Oncogenic mutations of SF3B1 and U2AF1 cause a gene-body RNAPII elongation defectRNAPII transcription elongation defect leads to transcription replication conflicts, DNA damage response, and changes to chromatin organization and H3K4me3 marksThe transcription elongation defect is linked to disruption of the early spliceosome formation through impaired interaction of HTATSF1 with mutant SF3B1.Changes to chromatin organization reveal potential therapeutic strategies by targeting the Sin3/HDAC pathway.

Context-Dependent Function of Long Noncoding RNA PURPL in Transcriptome Regulation during p53 Activation

PURPL is a p53-induced lncRNA that suppresses basal p53 levels. Here, we investigated PURPL upon p53 activation in liver cancer cells, where it is expressed at significantly higher levels than other cell types. Using isoform sequencing, we discovered novel PURPL transcripts that have a retained intron and/or previously unannotated exons. To determine PURPL function upon p53 activation, we performed transcriptome sequencing (RNA-Seq) after depleting PURPL using CRISPR interference (CRISPRi), followed by Nutlin treatment to induce p53. Strikingly, although loss of PURPL in untreated cells altered the expression of only 7 genes, loss of PURPL resulted in altered expression of ~800 genes upon p53 activation, revealing a context-dependent function of PURPL. Pathway analysis suggested that PURPL is important for fine-tuning the expression of specific genes required for mitosis. Consistent with these results, we observed a significant decrease in the percentage of mitotic cells upon PURPL depletion. Collectively, these data identify novel transcripts from the PURPL locus and suggest that PURPL delicately moderates the expression of mitotic genes in the context of p53 activation to control cell cycle arrest.

Chromatin dynamics associated with seed desiccation tolerance/sensitivity at early germination in Medicago truncatula

Desiccation tolerance (DT) has contributed greatly to the adaptation of land plants to severe water-deficient conditions. DT is mostly observed in reproductive parts in flowering plants such as seeds. The seed DT is lost at early post germination stage but is temporally re-inducible in 1 mm radicles during the so-called DT window following a PEG treatment before being permanently silenced in 5 mm radicles of germinating seeds. The molecular mechanisms that activate/reactivate/silence DT in developing and germinating seeds have not yet been elucidated. Here, we analyzed chromatin dynamics related to re-inducibility of DT before and after the DT window at early germination in Medicago truncatula radicles to determine if DT-associated genes were transcriptionally regulated at the chromatin levels. Comparative transcriptome analysis of these radicles identified 948 genes as DT re-induction-related genes, positively correlated with DT re-induction. ATAC-Seq analyses revealed that the chromatin state of genomic regions containing these genes was clearly modulated by PEG treatment and affected by growth stages with opened chromatin in 1 mm radicles with PEG (R1P); intermediate openness in 1 mm radicles without PEG (R1); and condensed chromatin in 5 mm radicles without PEG (R5). In contrast, we also showed that the 103 genes negatively correlated with the re-induction of DT did not show any transcriptional regulation at the chromatin level. Additionally, ChIP-Seq analyses for repressive marks H2AK119ub and H3K27me3 detected a prominent signal of H3K27me3 on the DT re-induction-related gene sequences at R5 but not in R1 and R1P. Moreover, no clear H2AK119ub marks was observed on the DT re-induction-related gene sequences at both developmental radicle stages, suggesting that silencing of DT process after germination will be mainly due to H3K27me3 marks by the action of the PRC2 complex, without involvement of PRC1 complex. The dynamic of chromatin changes associated with H3K27me3 were also confirmed on seed-specific genes encoding potential DT-related proteins such as LEAs, oleosins and transcriptional factors. However, several transcriptional factors did not show a clear link between their decrease of chromatin openness and H3K27me3 levels, suggesting that their accessibility may also be regulated by additional factors, such as other histone modifications. Finally, in order to make these comprehensive genome-wide analyses of transcript and chromatin dynamics useful to the scientific community working on early germination and DT, we generated a dedicated genome browser containing all these data and publicly available at https://iris.angers.inrae.fr/mtseedepiatlas/jbrowse/?data=Mtruncatula.

PRC1-mediated epigenetic programming is required to generate the ovarian reserve

The ovarian reserve defines the female reproductive lifespan, which in humans spans decades due to robust maintenance of meiotic arrest in oocytes residing in primordial follicles. Epigenetic reprogramming, including DNA demethylation, accompanies meiotic entry, but the chromatin changes that underpin the generation and preservation of ovarian reserves are poorly defined. We report that the Polycomb Repressive Complex 1 (PRC1) establishes repressive chromatin states in perinatal mouse oocytes that directly suppress the gene expression program of meiotic prophase-I and thereby enable the transition to dictyate arrest. PRC1 dysfuction causes depletion of the ovarian reserve and leads to premature ovarian failure. Our study demonstrates a fundamental role for PRC1-mediated gene silencing in female reproductive lifespan, and reveals a critical window of epigenetic programming required to establish ovarian reserve.

In humans, the ovarian reserve is maintained over decades by meiotic arrest of oocytes. Here the authors show that Polycomb Repressive Complex 1 (PRC1)-mediated epigenetic programming is essential for formation of ovarian reserve and thus female reproductive lifespan.

3D genome organization in the epithelial-mesenchymal transition spectrum

Background

The plasticity along the epithelial-mesenchymal transition (EMT) spectrum has been shown to be regulated by various epigenetic repertoires. Emerging evidence of local chromatin conformation changes suggests that regulation of EMT may occur at a higher order of three-dimensional genome level.

Results

We perform Hi-C analysis and combine ChIP-seq data across cancer cell lines representing different EMT states. We demonstrate that the epithelial and mesenchymal genes are regulated distinctively. We find that EMT genes are regulated within their topologically associated domains (TADs), with only a subset of mesenchymal genes being influenced by A/B compartment switches, indicating topological remodeling is required in the transcriptional regulation of these genes. At the TAD level, epithelial and mesenchymal genes are associated with different regulatory trajectories. The epithelial gene-residing TADs are enriched with H3K27me3 marks in the mesenchymal-like states. The mesenchymal gene-residing TADs, which do not show enrichment of H3K27me3 in epithelial-like states, exhibit increased interaction frequencies with regulatory elements in the mesenchymal-like states.

Conclusions

We propose a novel workflow coupling immunofluorescence and dielectrophoresis to unravel EMT heterogeneity at single-cell resolution. The predicted three-dimensional structures of chromosome 10, harboring Vimentin, identify cell clusters of different states. Our results pioneer a novel avenue to decipher the complexities underlying the regulation of EMT and may infer the barriers of plasticity in the 3D genome context.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02687-x.

Convergence of SIRT1 and ATR signaling to modulate replication origin dormancy

During routine genome duplication, many potential replication origins remain inactive or 'dormant'. Such origin dormancy is achieved, in part, by an interaction with the metabolic sensor SIRT1 deacetylase. We report here that dormant origins are a group of consistent, pre-determined genomic sequences that are distinguished from baseline (i.e. ordinarily active) origins by their preferential association with two phospho-isoforms of the helicase component MCM2. During normal unperturbed cell growth, baseline origins, but not dormant origins, associate with a form of MCM2 that is phosphorylated by DBF4-dependent kinase (DDK) on serine 139 (pS139-MCM2). This association facilitates the initiation of DNA replication from baseline origins. Concomitantly, SIRT1 inhibits Ataxia Telangiectasia and Rad3-related (ATR)-kinase-mediated phosphorylation of MCM2 on serine 108 (pS108-MCM2) by deacetylating the ATR-interacting protein DNA topoisomerase II binding protein 1 (TOPBP1), thereby preventing ATR recruitment to chromatin. In cells devoid of SIRT1 activity, or challenged by replication stress, this inhibition is circumvented, enabling ATR-mediated S108-MCM2 phosphorylation. In turn, pS108-MCM2 enables DDK-mediated phosphorylation on S139-MCM2 and facilitates replication initiation at dormant origins. These observations suggest that replication origin dormancy and activation are regulated by distinct post-translational MCM modifications that reflect a balance between SIRT1 activity and ATR signaling.

Assessing and assuring interoperability of a genomics file format

Motivation

Bioinformatics software tools operate largely through the use of specialized genomics file formats. Often these formats lack formal specification, making it difficult or impossible for the creators of these tools to robustly test them for correct handling of input and output. This causes problems in interoperability between different tools that, at best, wastes time and frustrates users. At worst, interoperability issues could lead to undetected errors in scientific results.

Results

We developed a new verification system, Acidbio, which tests for correct behavior in bioinformatics software packages. We crafted tests to unify correct behavior when tools encounter various edge cases—potentially unexpected inputs that exemplify the limits of the format. To analyze the performance of existing software, we tested the input validation of 80 Bioconda packages that parsed the Browser Extensible Data (BED) format. We also used a fuzzing approach to automatically perform additional testing. Of 80 software packages examined, 75 achieved less than 70% correctness on our test suite. We categorized multiple root causes for the poor performance of different types of software. Fuzzing detected other errors that the manually designed test suite could not. We also created a badge system that developers can use to indicate more precisely which BED variants their software accepts and to advertise the software’s performance on the test suite.

Availability and implementation

Acidbio is available at https://github.com/hoffmangroup/acidbio.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-wide investigations on regulatory functions of RECQ1 helicase

DNA helicase RECQ1 (also known as RECQL or RECQL1) is a candidate breast cancer susceptibility gene significantly correlated with clinical outcomes of sporadic breast cancer patients. Prior studies have suggested that RECQ1 maintains genomic stability by regulating a wide variety of core cellular functions including DNA replication, DNA damage response, and transcription. However, it is unclear which, if any, of these are the primary functions of RECQ1 as related to its role in suppressing breast cancer. We describe here an unbiased integrative genomics approach that enabled us to discover a previously unknown regulatory role of RECQ1 in promoting Estrogen Receptor alpha (ERα) expression and the expression of specific ERα target genes in ER positive breast cancer cells. We discuss potential future applications of similar experimental strategies in advancing the mechanistic understanding and elucidating specific new details of genome-wide functions of RECQ1 and other RecQ helicases in maintaining genomic stability and preventing cancer.

Transcriptional regulation by a RecQ helicase

RecQ helicases participate in a variety of DNA metabolic processes through their multiple biochemical activities. In vitro characterization and cellular studies have suggested that RECQ1 (also known as RECQL or RECQL1) performs its diverse functions through specific interactions with DNA and protein partners. We have taken an unbiased approach to determine the contribution of RECQ1 in genome maintenance and as a putative susceptibility factor in breast cancer. Here, we provide methodology to map the genome-wide binding sites of RECQ1 together with the profiling of RECQ1-dependent transcriptome to investigate its role in gene regulation. The described approach will be helpful to develop a mechanistic framework for elucidating critical functions of RECQ1 and other RecQ homologs in distinct chromatin and biological contexts.

High-throughput single-cell epigenomic profiling by targeted insertion of promoters (TIP-seq)

Chromatin profiling in single cells has been extremely challenging and almost exclusively limited to histone proteins. In cases where single-cell methods have shown promise, many require highly specialized equipment or cell type-specific protocols and are relatively low throughput. Here, we combine the advantages of tagmentation, linear amplification, and combinatorial indexing to produce a high-throughput single-cell DNA binding site mapping method that is simple, inexpensive, and capable of multiplexing several independent samples per experiment. Targeted insertion of promoters sequencing (TIP-seq) uses Tn5 fused to proteinA to insert a T7 RNA polymerase promoter adjacent to a chromatin protein of interest. Linear amplification of flanking DNA with T7 polymerase before sequencing library preparation provides ∼10-fold higher unique reads per single cell compared with other methods. We applied TIP-seq to map histone modifications, RNA polymerase II (RNAPII), and transcription factor CTCF binding sites in single human and mouse cells.

Altered function and differentiation of age-associated B cells contribute to the female bias in lupus mice

Differences in immune responses to viruses and autoimmune diseases such as systemic lupus erythematosus (SLE) can show sexual dimorphism. Age-associated B cells (ABC) are a population of CD11c+T-bet+ B cells critical for antiviral responses and autoimmune disorders. Absence of DEF6 and SWAP-70, two homologous guanine exchange factors, in double-knock-out (DKO) mice leads to a lupus-like syndrome in females marked by accumulation of ABCs. Here we demonstrate that DKO ABCs show sex-specific differences in cell number, upregulation of an ISG signature, and further differentiation. DKO ABCs undergo oligoclonal expansion and differentiate into both CD11c+ and CD11c- effector B cell populations with pathogenic and pro-inflammatory function as demonstrated by BCR sequencing and fate-mapping experiments. Tlr7 duplication in DKO males overrides the sex-bias and further augments the dissemination and pathogenicity of ABCs, resulting in severe pulmonary inflammation and early mortality. Thus, sexual dimorphism shapes the expansion, function and differentiation of ABCs that accompanies TLR7-driven immunopathogenesis.

Dynamics of replication origin over-activation

Safeguards against excess DNA replication are often dysregulated in cancer, and driving cancer cells towards over-replication is a promising therapeutic strategy. We determined DNA synthesis patterns in cancer cells undergoing partial genome re-replication due to perturbed regulatory interactions (re-replicating cells). These cells exhibited slow replication, increased frequency of replication initiation events, and a skewed initiation pattern that preferentially reactivated early-replicating origins. Unlike in cells exposed to replication stress, which activated a novel group of hitherto unutilized (dormant) replication origins, the preferred re-replicating origins arose from the same pool of potential origins as those activated during normal growth. Mechanistically, the skewed initiation pattern reflected a disproportionate distribution of pre-replication complexes on distinct regions of licensed chromatin prior to replication. This distinct pattern suggests that circumventing the strong inhibitory interactions that normally prevent excess DNA synthesis can occur via at least two pathways, each activating a distinct set of replication origins.

Chromatin Remodeling and Immediate Early Gene Activation by SLFN11 in Response to Replication Stress

Schlafen 11 (SLFN11) was recently discovered as a cellular restriction factor against replication stress. Here, we show that SLFN11 increases chromatin accessibility genome wide, prominently at active promoters in response to replication stress induced by the checkpoint kinase 1 (CHK1) inhibitor prexasertib or the topoisomerase I (TOP1) inhibitor camptothecin. Concomitantly, SLFN11 selectively activates cellular stress response pathways by inducing the transcription of the immediate early genes (IEGs), including JUN, FOS, EGR1, NFKB2, and ATF3, together with the cell cycle arrest genes CDKN1A (p21^WAF1) and GADD45. Both chromatin remodeling and IEG activation require the putative ATPase and helicase activity of SLFN11, whereas canonical extrinsic IEG activation is SLFN11 independent. SLFN11-dependent IEG activation by camptothecin is also observed across 55 non-isogenic NCI-60 cell lines. We conclude that SLFN11 acts as a global regulator of chromatin structure and an intrinsic IEG activator with the potential to engage the innate immune activation in response to replicative stress.

Genome-Wide Analysis Unveils DNA Helicase RECQ1 as a Regulator of Estrogen Response Pathway in Breast Cancer Cells

Susceptibility to breast cancer is significantly increased in individuals with germ line mutations in RECQ1 (also known as RECQL or RECQL1), a gene encoding a DNA helicase essential for genome maintenance. We previously reported that RECQ1 expression predicts clinical outcomes for sporadic breast cancer patients stratified by estrogen receptor (ER) status. Here, we utilized an unbiased integrative genomics approach to delineate a cross talk between RECQ1 and ERα, a known master regulatory transcription factor in breast cancer. We found that expression of ESR1, the gene encoding ERα, is directly activated by RECQ1. More than 35% of RECQ1 binding sites were cobound by ERα genome-wide. Mechanistically, RECQ1 cooperates with FOXA1, the pioneer transcription factor for ERα, to enhance chromatin accessibility at the ESR1 regulatory regions in a helicase activity-dependent manner. In clinical ERα-positive breast cancers treated with endocrine therapy, high RECQ1 and high FOXA1 coexpressing tumors were associated with better survival. Collectively, these results identify RECQ1 as a novel cofactor for ERα and uncover a previously unknown mechanism by which RECQ1 regulates disease-driving gene expression in ER-positive breast cancer cells.

Initiation of a conserved trophectoderm program in human, cow and mouse embryos

Current understandings of cell specification in early mammalian pre-implantation development are based mainly on mouse studies. The first lineage differentiation event occurs at the morula stage, with outer cells initiating a trophectoderm (TE) placental progenitor program. The inner cell mass arises from inner cells during subsequent developmental stages and comprises precursor cells of the embryo proper and yolk sac1. Recent gene-expression analyses suggest that the mechanisms that regulate early lineage specification in the mouse may differ in other mammals, including human2-5 and cow6. Here we show the evolutionary conservation of a molecular cascade that initiates TE segregation in human, cow and mouse embryos. At the morula stage, outer cells acquire an apical-basal cell polarity, with expression of atypical protein kinase C (aPKC) at the contact-free domain, nuclear expression of Hippo signalling pathway effectors and restricted expression of TE-associated factors such as GATA3, which suggests initiation of a TE program. Furthermore, we demonstrate that inhibition of aPKC by small-molecule pharmacological modulation or Trim-Away protein depletion impairs TE initiation at the morula stage. Our comparative embryology analysis provides insights into early lineage specification and suggests that a similar mechanism initiates a TE program in human, cow and mouse embryos.

Genome-Wide Analysis of the FOXA1 Transcriptional Network Identifies Novel Protein-Coding and Long Noncoding RNA Targets in Colorectal Cancer Cells

Differentiation status of tumors is correlated with metastatic potential and malignancy. FOXA1 (forkhead box A1) is a transcription factor known to regulate differentiation in certain tissues. Here, we investigate FOXA1 function in human colorectal cancer (CRC). We found that FOXA1 is robustly expressed in the normal human colon but significantly downregulated in colon adenocarcinoma. Applying FOXA1 chromatin immunoprecipitation coupled with deep sequencing and transcriptome analysis upon FOXA1 knockdown in well-differentiated CRC cells and FOXA1 overexpression in poorly differentiated CRC cells, we identified novel protein-coding and lncRNA genes regulated by FOXA1. Among the numerous novel FOXA1 targets we identified, we focused on CEACAM5, a tumor marker and facilitator of cell adhesion. We show that FOXA1 binds to a distal enhancer downstream of CEACAM5 and strongly activates its expression. Consistent with these data, we show that FOXA1 inhibits anoikis in CRC cells. Collectively, our results uncover novel protein-coding and noncoding targets of FOXA1 and suggest a vital role of FOXA1 in enhancing CEACAM5 expression and anoikis resistance in CRC cells.