Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells.

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

The cis-Regulatory Atlas of the Mouse Immune System.

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

A microRNA expression and regulatory element activity atlas of the mouse immune system.

To better define the control of immune system regulation, we generated an atlas of microRNA (miRNA) expression from 63 mouse immune cell populations and connected these signatures with assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation followed by sequencing (ChIP-seq) and nascent RNA profiles to establish a map of miRNA promoter and enhancer usage in immune cells. miRNA complexity was relatively low, with >90% of the miRNA compartment of each population comprising <75 miRNAs; however, each cell type had a unique miRNA signature. Integration of miRNA expression with chromatin accessibility revealed putative regulatory elements for differentially expressed miRNAs, including miR-21a, miR-146a and miR-223. The integrated maps suggest that many miRNAs utilize multiple promoters to reach high abundance and identified dominant and divergent miRNA regulatory elements between lineages and during development that may be used by clustered miRNAs, such as miR-99a/let-7c/miR-125b, to achieve distinct expression. These studies, with web-accessible data, help delineate the cis-regulatory elements controlling miRNA signatures of the immune system.

Dual Roles of Poly(dA:dT) Tracts in Replication Initiation and Fork Collapse.

Replication origins, fragile sites, and rDNA have been implicated as sources of chromosomal instability. However, the defining genomic features of replication origins and fragile sites are among the least understood elements of eukaryote genomes. Here, we map sites of replication initiation and breakage in primary cells at high resolution. We find that replication initiates between transcribed genes within nucleosome-depleted structures established by long asymmetrical poly(dA:dT) tracts flanking the initiation site. Paradoxically, long (>20 bp) (dA:dT) tracts are also preferential sites of polar replication fork stalling and collapse within early-replicating fragile sites (ERFSs) and late-replicating common fragile sites (CFSs) and at the rDNA replication fork barrier. Poly(dA:dT) sequences are fragile because long single-strand poly(dA) stretches at the replication fork are unprotected by the replication protein A (RPA). We propose that the evolutionary expansion of poly(dA:dT) tracts in eukaryotic genomes promotes replication initiation, but at the cost of chromosome fragility.

Co-optation of Tandem DNA Repeats for the Maintenance of Mesenchymal Identity.

Tandem repeats (TRs) are generated by DNA replication errors and retain a high level of instability, which in principle would make them unsuitable for integration into gene regulatory networks. However, the appearance of DNA sequence motifs recognized by transcription factors may turn TRs into functional cis-regulatory elements, thus favoring their stabilization in genomes. Here, we show that, in human cells, the transcriptional repressor ZEB1, which promotes the maintenance of mesenchymal features largely by suppressing epithelial genes and microRNAs, occupies TRs harboring dozens of copies of its DNA-binding motif within genomic loci relevant for maintenance of epithelial identity. The deletion of one such TR caused quasi-mesenchymal cancer cells to reacquire epithelial features, partially recapitulating the effects of ZEB1 gene deletion. These data demonstrate that the high density of identical motifs in TRs can make them suitable platforms for recruitment of transcriptional repressors, thus promoting their exaptation into pre-existing cis-regulatory networks.

A LINE1-Nucleolin Partnership Regulates Early Development and ESC Identity.

Transposable elements represent nearly half of mammalian genomes and are generally described as parasites, or "junk DNA." The LINE1 retrotransposon is the most abundant class and is thought to be deleterious for cells, yet it is paradoxically highly expressed during early development. Here, we report that LINE1 plays essential roles in mouse embryonic stem cells (ESCs) and pre-implantation embryos. In ESCs, LINE1 acts as a nuclear RNA scaffold that recruits Nucleolin and Kap1/Trim28 to repress Dux, the master activator of a transcriptional program specific to the 2-cell embryo. In parallel, LINE1 RNA mediates binding of Nucleolin and Kap1 to rDNA, promoting rRNA synthesis and ESC self-renewal. In embryos, LINE1 RNA is required for Dux silencing, synthesis of rRNA, and exit from the 2-cell stage. The results reveal an essential partnership between LINE1 RNA, Nucleolin, Kap1, and peri-nucleolar chromatin in the regulation of transcription, developmental potency, and ESC self-renewal.

The RNA Exosome Syncs IAV-RNAPII Transcription to Promote Viral Ribogenesis and Infectivity.

The nuclear RNA exosome is an essential multi-subunit complex that controls RNA homeostasis. Congenital mutations in RNA exosome genes are associated with neurodegenerative diseases. Little is known about the role of the RNA exosome in the cellular response to pathogens. Here, using NGS and human and mouse genetics, we show that influenza A virus (IAV) ribogenesis and growth are suppressed by impaired RNA exosome activity. Mechanistically, the nuclear RNA exosome coordinates the initial steps of viral transcription with RNAPII at host promoters. The viral polymerase complex co-opts the nuclear RNA exosome complex and cellular RNAs en route to 3' end degradation. Exosome deficiency uncouples chromatin targeting of the viral polymerase complex and the formation of cellular:viral RNA hybrids, which are essential RNA intermediates that license transcription of antisense genomic viral RNAs. Our results suggest that evolutionary arms races have shaped the cellular RNA quality control machinery.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1.

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Single-cell RNA-seq reveals TOX as a key regulator of CD8+ T cell persistence in chronic infection.

Progenitor-like CD8+ T cells mediate long-term immunity to chronic infection and cancer and respond potently to immune checkpoint blockade. These cells share transcriptional regulators with memory precursor cells, including T cell-specific transcription factor 1 (TCF1), but it is unclear whether they adopt distinct programs to adapt to the immunosuppressive environment. By comparing the single-cell transcriptomes and epigenetic profiles of CD8+ T cells responding to acute and chronic viral infections, we found that progenitor-like CD8+ T cells became distinct from memory precursor cells before the peak of the T cell response. We discovered a coexpression gene module containing Tox that exhibited higher transcriptional activity associated with more abundant active histone marks in progenitor-like cells than memory precursor cells. Moreover, thymocyte selection-associated high mobility group box protein TOX (TOX) promoted the persistence of antiviral CD8+ T cells and was required for the programming of progenitor-like CD8+ T cells. Thus, long-term CD8+ T cell immunity to chronic viral infection requires unique transcriptional and epigenetic programs associated with the transcription factor TOX.

SnapShot: Epigenomic Assays.

This SnapShot depicts key sequencing-based methods used in the analysis of epigenomes, including (1)bisulfite sequencing, (2) chromatin immunoprecipiation sequencing, (3) determination of open chromatin, and (4) 3D chromatin capture.

Pooled ChIP-Seq Links Variation in Transcription Factor Binding to Complex Disease Risk.

Cis-regulatory elements such as transcription factor (TF) binding sites can be identified genome-wide, but it remains far more challenging to pinpoint genetic variants affecting TF binding. Here, we introduce a pooling-based approach to mapping quantitative trait loci (QTLs) for molecular-level traits. Applying this to five TFs and a histone modification, we mapped thousands of cis-acting QTLs, with over 25-fold lower cost compared to standard QTL mapping. We found that single genetic variants frequently affect binding of multiple TFs, and CTCF can recruit all five TFs to its binding sites. These QTLs often affect local chromatin and transcription but can also influence long-range chromosomal contacts, demonstrating a role for natural genetic variation in chromosomal architecture. Thousands of these QTLs have been implicated in genome-wide association studies, providing candidate molecular mechanisms for many disease risk loci and suggesting that TF binding variation may underlie a large fraction of human phenotypic variation.

Histone Acetylome-wide Association Study of Autism Spectrum Disorder.

The association of histone modification changes with autism spectrum disorder (ASD) has not been systematically examined. We conducted a histone acetylome-wide association study (HAWAS) by performing H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) on 257 postmortem samples from ASD and matched control brains. Despite etiological heterogeneity, ≥68% of syndromic and idiopathic ASD cases shared a common acetylome signature at >5,000 cis-regulatory elements in prefrontal and temporal cortex. Similarly, multiple genes associated with rare genetic mutations in ASD showed common "epimutations." Acetylome aberrations in ASD were not attributable to genetic differentiation at cis-SNPs and highlighted genes involved in synaptic transmission, ion transport, epilepsy, behavioral abnormality, chemokinesis, histone deacetylation, and immunity. By correlating histone acetylation with genotype, we discovered >2,000 histone acetylation quantitative trait loci (haQTLs) in human brain regions, including four candidate causal variants for psychiatric diseases. Due to the relative stability of histone modifications postmortem, we anticipate that the HAWAS approach will be applicable to multiple diseases.

Direct GR Binding Sites Potentiate Clusters of TF Binding across the Human Genome.

The glucocorticoid receptor (GR) binds the human genome at >10,000 sites but only regulates the expression of hundreds of genes. To determine the functional effect of each site, we measured the glucocorticoid (GC) responsive activity of nearly all GR binding sites (GBSs) captured using chromatin immunoprecipitation (ChIP) in A549 cells. 13% of GBSs assayed had GC-induced activity. The responsive sites were defined by direct GR binding via a GC response element (GRE) and exclusively increased reporter-gene expression. Meanwhile, most GBSs lacked GC-induced reporter activity. The non-responsive sites had epigenetic features of steady-state enhancers and clustered around direct GBSs. Together, our data support a model in which clusters of GBSs observed with ChIP-seq reflect interactions between direct and tethered GBSs over tens of kilobases. We further show that those interactions can synergistically modulate the activity of direct GBSs and may therefore play a major role in driving gene activation in response to GCs.

The variation and evolution of complete human centromeres.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.

α-ketoglutarate orchestrates macrophage activation through metabolic and epigenetic reprogramming.

Glutamine metabolism provides synergistic support for macrophage activation and elicitation of desirable immune responses; however, the underlying mechanisms regulated by glutamine metabolism to orchestrate macrophage activation remain unclear. Here we show that the production of α-ketoglutarate (αKG) via glutaminolysis is important for alternative (M2) activation of macrophages, including engagement of fatty acid oxidation (FAO) and Jmjd3-dependent epigenetic reprogramming of M2 genes. This M2-promoting mechanism is further modulated by a high αKG/succinate ratio, whereas a low ratio strengthens the proinflammatory phenotype in classically activated (M1) macrophages. As such, αKG contributes to endotoxin tolerance after M1 activation. This study reveals new mechanistic regulations by which glutamine metabolism tailors the immune responses of macrophages through metabolic and epigenetic reprogramming.

Antagonism of B cell enhancer networks by STAT5 drives leukemia and poor patient survival.

The transcription factor STAT5 has a critical role in B cell acute lymphoblastic leukemia (B-ALL). How STAT5 mediates this effect is unclear. Here we found that activation of STAT5 worked together with defects in signaling components of the precursor to the B cell antigen receptor (pre-BCR), including defects in BLNK, BTK, PKCß, NF-κB1 and IKAROS, to initiate B-ALL. STAT5 antagonized the transcription factors NF-κB and IKAROS by opposing regulation of shared target genes. Super-enhancers showed enrichment for STAT5 binding and were associated with an opposing network of transcription factors, including PAX5, EBF1, PU.1, IRF4 and IKAROS. Patients with a high ratio of active STAT5 to NF-κB or IKAROS had more-aggressive disease. Our studies indicate that an imbalance of two opposing transcriptional programs drives B-ALL and suggest that restoring the balance of these pathways might inhibit B-ALL.

The microRNA miR-31 inhibits CD8+ T cell function in chronic viral infection.

During infection, antigen-specific T cells undergo tightly regulated developmental transitions controlled by transcriptional and post-transcriptional regulation of gene expression. We found that the microRNA miR-31 was strongly induced by activation of the T cell antigen receptor (TCR) in a pathway involving calcium and activation of the transcription factor NFAT. During chronic infection with lymphocytic choriomeningitis virus (LCMV) clone 13, miR-31-deficent mice recovered from clinical disease, while wild-type mice continued to show signs of disease. This disease phenotype was explained by the presence of larger numbers of cytokine-secreting LCMV-specific CD8+ T cells in miR-31-deficent mice than in wild-type mice. Mechanistically, miR-31 increased the sensitivity of T cells to type I interferons, which interfered with effector T cell function and increased the expression of several proteins related to T cell dysfunction during chronic infection. These studies identify miR-31 as an important regulator of T cell exhaustion in chronic infection.

Type I interferons and the cytokine TNF cooperatively reprogram the macrophage epigenome to promote inflammatory activation.

Cross-regulation of Toll-like receptor (TLR) responses by cytokines is essential for effective host defense, avoidance of toxicity and homeostasis, but the underlying mechanisms are not well understood. Our comprehensive epigenomics approach to the analysis of human macrophages showed that the proinflammatory cytokines TNF and type I interferons induced transcriptional cascades that altered chromatin states to broadly reprogram responses induced by TLR4. TNF tolerized genes encoding inflammatory molecules to prevent toxicity while preserving the induction of genes encoding antiviral and metabolic molecules. Type I interferons potentiated the inflammatory function of TNF by priming chromatin to prevent the silencing of target genes of the transcription factor NF-κB that encode inflammatory molecules. The priming of chromatin enabled robust transcriptional responses to weak upstream signals. Similar chromatin regulation occurred in human diseases. Our findings reveal that signaling crosstalk between interferons and TNF is integrated at the level of chromatin to reprogram inflammatory responses, and identify previously unknown functions and mechanisms of action of these cytokines.

The transcriptional coactivator TAZ regulates reciprocal differentiation of TH17 cells and Treg cells.

An imbalance in the lineages of immunosuppressive regulatory T cells (Treg cells) and the inflammatory TH17 subset of helper T cells leads to the development of autoimmune and/or inflammatory disease. Here we found that TAZ, a coactivator of TEAD transcription factors of Hippo signaling, was expressed under TH17 cell-inducing conditions and was required for TH17 differentiation and TH17 cell-mediated inflammatory diseases. TAZ was a critical co-activator of the TH17-defining transcription factor RORγt. In addition, TAZ attenuated Treg cell development by decreasing acetylation of the Treg cell master regulator Foxp3 mediated by the histone acetyltransferase Tip60, which targeted Foxp3 for proteasomal degradation. In contrast, under Treg cell-skewing conditions, TEAD1 expression and sequestration of TAZ from the transcription factors RORγt and Foxp3 promoted Treg cell differentiation. Furthermore, deficiency in TAZ or overexpression of TEAD1 induced Treg cell differentiation, whereas expression of a transgene encoding TAZ or activation of TAZ directed TH17 cell differentiation. Our results demonstrate a pivotal role for TAZ in regulating the differentiation of Treg cells and TH17 cells.

Global changes in the RNA binding specificity of HIV-1 gag regulate virion genesis.

The HIV-1 Gag protein orchestrates all steps of virion genesis, including membrane targeting and RNA recruitment into virions. Using crosslinking-immunoprecipitation (CLIP) sequencing, we uncover several dramatic changes in the RNA-binding properties of Gag that occur during virion genesis, coincident with membrane binding, multimerization, and proteolytic maturation. Prior to assembly, and after virion assembly and maturation, the nucleocapsid domain of Gag preferentially binds to psi and Rev Response elements in the viral genome, and GU-rich mRNA sequences. However, during virion genesis, this specificity transiently changes in a manner that facilitates genome packaging; nucleocapsid binds to many sites on the HIV-1 genome and to mRNA sequences with a HIV-1-like, A-rich nucleotide composition. Additionally, we find that the matrix domain of Gag binds almost exclusively to specific tRNAs in the cytosol, and this association regulates Gag binding to cellular membranes.