The Energetics and Physiological Impact of Cohesin Extrusion.

Cohesin extrusion is thought to play a central role in establishing the architecture of mammalian genomes. However, extrusion has not been visualized in vivo, and thus, its functional impact and energetics are unknown. Using ultra-deep Hi-C, we show that loop domains form by a process that requires cohesin ATPases. Once formed, however, loops and compartments are maintained for hours without energy input. Strikingly, without ATP, we observe the emergence of hundreds of CTCF-independent loops that link regulatory DNA. We also identify architectural "stripes," where a loop anchor interacts with entire domains at high frequency. Stripes often tether super-enhancers to cognate promoters, and in B cells, they facilitate Igh transcription and recombination. Stripe anchors represent major hotspots for topoisomerase-mediated lesions, which promote chromosomal translocations and cancer. In plasmacytomas, stripes can deregulate Igh-translocated oncogenes. We propose that higher organisms have coopted cohesin extrusion to enhance transcription and recombination, with implications for tumor development.

RAG Represents a Widespread Threat to the Lymphocyte Genome.

The RAG1 endonuclease, together with its cofactor RAG2, is essential for V(D)J recombination but is a potent threat to genome stability. The sources of RAG1 mis-targeting and the mechanisms that have evolved to suppress it are poorly understood. Here, we report that RAG1 associates with chromatin at thousands of active promoters and enhancers in the genome of developing lymphocytes. The mouse and human genomes appear to have responded by reducing the abundance of "cryptic" recombination signals near RAG1 binding sites. This depletion operates specifically on the RSS heptamer, whereas nonamers are enriched at RAG1 binding sites. Reversing this RAG-driven depletion of cleavage sites by insertion of strong recombination signals creates an ectopic hub of RAG-mediated V(D)J recombination and chromosomal translocations. Our findings delineate rules governing RAG binding in the genome, identify areas at risk of RAG-mediated damage, and highlight the evolutionary struggle to accommodate programmed DNA damage in developing lymphocytes.

B cell super-enhancers and regulatory clusters recruit AID tumorigenic activity.

The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA(+) enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

Global regulation of promoter melting in naive lymphocytes.

Lymphocyte activation is initiated by a global increase in messenger RNA synthesis. However, the mechanisms driving transcriptome amplification during the immune response are unknown. By monitoring single-stranded DNA genome wide, we show that the genome of naive cells is poised for rapid activation. In G0, ∼90% of promoters from genes to be expressed in cycling lymphocytes are polymerase loaded but unmelted and support only basal transcription. Furthermore, the transition from abortive to productive elongation is kinetically limiting, causing polymerases to accumulate nearer to transcription start sites. Resting lymphocytes also limit the expression of the transcription factor IIH complex, including XPB and XPD helicases involved in promoter melting and open complex extension. To date, two rate-limiting steps have been shown to control global gene expression in eukaryotes: preinitiation complex assembly and polymerase pausing. Our studies identify promoter melting as a third key regulatory step and propose that this mechanism ensures a prompt lymphocyte response to invading pathogens.

Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation.

A key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

c-Myc is a universal amplifier of expressed genes in lymphocytes and embryonic stem cells.

The c-Myc HLH-bZIP protein has been implicated in physiological or pathological growth, proliferation, apoptosis, metabolism, and differentiation at the cellular, tissue, or organismal levels via regulation of numerous target genes. No principle yet unifies Myc action due partly to an incomplete inventory and functional accounting of Myc's targets. To observe Myc target expression and function in a system where Myc is temporally and physiologically regulated, the transcriptomes and the genome-wide distributions of Myc, RNA polymerase II, and chromatin modifications were compared during lymphocyte activation and in ES cells as well. A remarkably simple rule emerged from this quantitative analysis: Myc is not an on-off specifier of gene activity, but is a nonlinear amplifier of expression, acting universally at active genes, except for immediate early genes that are strongly induced before Myc. This rule of Myc action explains the vast majority of Myc biology observed in literature.

Translocation-capture sequencing reveals the extent and nature of chromosomal rearrangements in B lymphocytes.

Chromosomal rearrangements, including translocations, require formation and joining of DNA double strand breaks (DSBs). These events disrupt the integrity of the genome and are frequently involved in producing leukemias, lymphomas and sarcomas. Despite the importance of these events, current understanding of their genesis is limited. To examine the origins of chromosomal rearrangements we developed Translocation Capture Sequencing (TC-Seq), a method to document chromosomal rearrangements genome-wide, in primary cells. We examined over 180,000 rearrangements obtained from 400 million B lymphocytes, revealing that proximity between DSBs, transcriptional activity and chromosome territories are key determinants of genome rearrangement. Specifically, rearrangements tend to occur in cis and to transcribed genes. Finally, we find that activation-induced cytidine deaminase (AID) induces the rearrangement of many genes found as translocation partners in mature B cell lymphoma.

The in vivo pattern of binding of RAG1 and RAG2 to antigen receptor loci.

The critical initial step in V(D)J recombination, binding of RAG1 and RAG2 to recombination signal sequences flanking antigen receptor V, D, and J gene segments, has not previously been characterized in vivo. Here, we demonstrate that RAG protein binding occurs in a highly focal manner to a small region of active chromatin encompassing Ig kappa and Tcr alpha J gene segments and Igh and Tcr beta J and J-proximal D gene segments. Formation of these small RAG-bound regions, which we refer to as recombination centers, occurs in a developmental stage- and lineage-specific manner. Each RAG protein is independently capable of specific binding within recombination centers. While RAG1 binding was detected only at regions containing recombination signal sequences, RAG2 binds at thousands of sites in the genome containing histone 3 trimethylated at lysine 4. We propose that recombination centers coordinate V(D)J recombination by providing discrete sites within which gene segments are captured for recombination.

Activation-induced cytidine deaminase targets DNA at sites of RNA polymerase II stalling by interaction with Spt5.

Activation-induced cytidine deaminase (AID) initiates antibody gene diversification by creating U:G mismatches. However, AID is not specific for antibody genes; Off-target lesions can activate oncogenes or cause chromosome translocations. Despite its importance in these transactions little is known about how AID finds its targets. We performed an shRNA screen to identify factors required for class switch recombination (CSR) of antibody loci. We found that Spt5, a factor associated with stalled RNA polymerase II (Pol II) and single stranded DNA (ssDNA), is required for CSR. Spt5 interacts with AID, it facilitates association between AID and Pol II, and AID recruitment to its Ig and non-Ig targets. ChIP-seq experiments reveal that Spt5 colocalizes with AID and stalled Pol II. Further, Spt5 accumulation at sites of Pol II stalling is predictive of AID-induced mutation. We propose that AID is targeted to sites of Pol II stalling in part via its association with Spt5.

Myc Regulates Chromatin Decompaction and Nuclear Architecture during B Cell Activation.

50 years ago, Vincent Allfrey and colleagues discovered that lymphocyte activation triggers massive acetylation of chromatin. However, the molecular mechanisms driving epigenetic accessibility are still unknown. We here show that stimulated lymphocytes decondense chromatin by three differentially regulated steps. First, chromatin is repositioned away from the nuclear periphery in response to global acetylation. Second, histone nanodomain clusters decompact into mononucleosome fibers through a mechanism that requires Myc and continual energy input. Single-molecule imaging shows that this step lowers transcription factor residence time and non-specific collisions during sampling for DNA targets. Third, chromatin interactions shift from long range to predominantly short range, and CTCF-mediated loops and contact domains double in numbers. This architectural change facilitates cognate promoter-enhancer contacts and also requires Myc and continual ATP production. Our results thus define the nature and transcriptional impact of chromatin decondensation and reveal an unexpected role for Myc in the establishment of nuclear topology in mammalian cells.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes.

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

53BP1 alters the landscape of DNA rearrangements and suppresses AID-induced B cell lymphoma.

Deficiencies in factors that regulate the DNA damage response enhance the incidence of malignancy by destabilizing the genome. However, the precise influence of the DNA damage response on regulation of cancer-associated rearrangements is not well defined. Here we examine the genome-wide impact of tumor protein P53-binding protein 1 (53BP1) deficiency in lymphoma and translocation. While both activation-induced cytidine deaminase (AID) and 53BP1 have been associated with cancer in humans, neither AID overexpression nor loss of 53BP1 is sufficient to produce malignancy. However, the combination of 53BP1 deficiency and AID deregulation results in B cell lymphoma. Deep sequencing of the genome of 53BP1(-/-) cancer cells and translocation capture sequencing (TC-Seq) of primary 53BP1(-/-) B cells revealed that their chromosomal rearrangements differ from those found in wild-type cells in that they show increased DNA end resection. Moreover, loss of 53BP1 alters the translocatome by increasing rearrangements to intergenic regions.

Regulation of microRNA expression and abundance during lymphopoiesis.

Although the cellular concentration of miRNAs is critical to their function, how miRNA expression and abundance are regulated during ontogeny is unclear. We applied miRNA-, mRNA-, and ChIP-Seq to characterize the microRNome during lymphopoiesis within the context of the transcriptome and epigenome. We show that lymphocyte-specific miRNAs are either tightly controlled by polycomb group-mediated H3K27me3 or maintained in a semi-activated epigenetic state prior to full expression. Because of miRNA biogenesis, the cellular concentration of mature miRNAs does not typically reflect transcriptional changes. However, we uncover a subset of miRNAs for which abundance is dictated by miRNA gene expression. We confirm that concentration of 5p and 3p miRNA strands depends largely on free energy properties of miRNA duplexes. Unexpectedly, we also find that miRNA strand accumulation can be developmentally regulated. Our data provide a comprehensive map of immunity's microRNome and reveal the underlying epigenetic and transcriptional forces that shape miRNA homeostasis.

DNA damage defines sites of recurrent chromosomal translocations in B lymphocytes.

Recurrent chromosomal translocations underlie both haematopoietic and solid tumours. Their origin has been ascribed to selection of random rearrangements, targeted DNA damage, or frequent nuclear interactions between translocation partners; however, the relative contribution of each of these elements has not been measured directly or on a large scale. Here we examine the role of nuclear architecture and frequency of DNA damage in the genesis of chromosomal translocations by measuring these parameters simultaneously in cultured mouse B lymphocytes. In the absence of recurrent DNA damage, translocations between Igh or Myc and all other genes are directly related to their contact frequency. Conversely, translocations associated with recurrent site-directed DNA damage are proportional to the rate of DNA break formation, as measured by replication protein A accumulation at the site of damage. Thus, non-targeted rearrangements reflect nuclear organization whereas DNA break formation governs the location and frequency of recurrent translocations, including those driving B-cell malignancies.

A mouse model of HIES reveals pro- and anti-inflammatory functions of STAT3.

Mutations of STAT3 underlie the autosomal dominant form of hyperimmunoglobulin E syndrome (HIES). STAT3 has critical roles in immune cells and thus, hematopoietic stem cell transplantation (HSCT), might be a reasonable therapeutic strategy in this disease. However, STAT3 also has critical functions in nonhematopoietic cells and dissecting the protean roles of STAT3 is limited by the lethality associated with germline deletion of Stat3. Thus, predicting the efficacy of HSCT for HIES is difficult. To begin to dissect the importance of STAT3 in hematopoietic and nonhematopoietic cells as it relates to HIES, we generated a mouse model of this disease. We found that these transgenic mice recapitulate multiple aspects of HIES, including elevated serum IgE and failure to generate Th17 cells. We found that these mice were susceptible to bacterial infection that was partially corrected by HSCT using wild-type bone marrow, emphasizing the role played by the epithelium in the pathophysiology of HIES.

Mouse model of endemic Burkitt translocations reveals the long-range boundaries of Ig-mediated oncogene deregulation.

Human Burkitt lymphomas are divided into two main clinical variants: the endemic form, affecting African children infected with malaria and the Epstein-Barr virus, and the sporadic form, distributed across the rest of the world. However, whereas sporadic translocations decapitate Myc from 5' proximal regulatory elements, most endemic events occur hundreds of kilobases away from Myc. The origin of these rearrangements and how they deregulate oncogenes at such distances remain unclear. We here recapitulate endemic Burkitt lymphoma-like translocations in plasmacytomas from uracil N-glycosylase and activation-induced cytidine deaminase-deficient mice. Mapping of translocation breakpoints using an acetylated histone H3 lysine 9 chromatin immunoprecipitation sequencing approach reveals Igh fusions up to ∼350 kb upstream of Myc or the related oncogene Mycn. A comprehensive analysis of epigenetic marks, PolII recruitment, and transcription in tumor cells demonstrates that the 3' Igh enhancer (Eα) vastly remodels ∼450 kb of chromatin into translocated sequences, leading to significant polymerase occupancy and constitutive oncogene expression. We show that this long-range epigenetic reprogramming is directly proportional to the physical interaction of Eα with translocated sites. Our studies thus uncover the extent of epigenetic remodeling by Ig 3' enhancers and provide a rationale for the long-range deregulation of translocated oncogenes in endemic Burkitt lymphomas. The data also shed light on the origin of endemic-like chromosomal rearrangements.

The folliculin-FNIP1 pathway deleted in human Birt-Hogg-Dubé syndrome is required for murine B-cell development.

Birt-Hogg-Dubé (BHD) syndrome is an autosomal dominant disorder characterized by cutaneous fibrofolliculomas, pulmonary cysts, and kidney malignancies. Affected individuals carry germ line mutations in folliculin (FLCN), a tumor suppressor gene that becomes biallelically inactivated in kidney tumors by second-hit mutations. Similar to other factors implicated in kidney cancer, FLCN has been shown to modulate activation of mammalian target of rapamycin (mTOR). However, its precise in vivo function is largely unknown because germ line deletion of Flcn results in early embryonic lethality in animal models. Here, we describe mice deficient in the newly characterized folliculin-interacting protein 1 (Fnip1). In contrast to Flcn, Fnip1(-/-) mice develop normally, are not susceptible to kidney neoplasia, but display a striking pro-B cell block that is entirely independent of mTOR activity. We show that this developmental arrest results from rapid caspase-induced pre-B cell death, and that a Bcl2 transgene reconstitutes mature B-cell populations, respectively. We also demonstrate that conditional deletion of Flcn recapitulates the pro-B cell arrest of Fnip1(-/-) mice. Our studies thus demonstrate that the FLCN-FNIP complex deregulated in BHD syndrome is absolutely required for B-cell differentiation, and that it functions through both mTOR-dependent and independent pathways.

Exaptation of Inactivated Host Enzymes for Structural Roles in Orthopoxviruses and Novel Folds of Virus Proteins Revealed by Protein Structure Modeling.

Viruses with large, double-stranded DNA genomes captured the majority of their genes from their hosts at different stages of evolution. The origins of many virus genes are readily detected through significant sequence similarity with cellular homologs. In particular, this is the case for virus enzymes, such as DNA and RNA polymerases or nucleotide kinases, that retain their catalytic activity after capture by an ancestral virus. However, a large fraction of virus genes have no readily detectable cellular homologs, meaning that their origins remain enigmatic. We explored the potential origins of such proteins that are encoded in the genomes of orthopoxviruses, a thoroughly studied virus genus that includes major human pathogens. To this end, we used AlphaFold2 to predict the structures of all 214 proteins that are encoded by orthopoxviruses. Among the proteins of unknown provenance, structure prediction yielded clear indications of origin for 14 of them and validated several inferences that were previously made via sequence analysis. A notable emerging trend is the exaptation of enzymes from cellular organisms for nonenzymatic, structural roles in virus reproduction that is accompanied by the disruption of catalytic sites and by an overall drastic divergence that precludes homology detection at the sequence level. Among the 16 orthopoxvirus proteins that were found to be inactivated enzyme derivatives are the poxvirus replication processivity factor A20, which is an inactivated NAD-dependent DNA ligase; the major core protein A3, which is an inactivated deubiquitinase; F11, which is an inactivated prolyl hydroxylase; and more similar cases. For nearly one-third of the orthopoxvirus virion proteins, no significantly similar structures were identified, suggesting exaptation with subsequent major structural rearrangement that yielded unique protein folds. IMPORTANCE Protein structures are more strongly conserved in evolution than are amino acid sequences. Comparative structural analysis is particularly important for inferring the origins of viral proteins that typically evolve at high rates. We used a powerful protein structure modeling method, namely, AlphaFold2, to model the structures of all orthopoxvirus proteins and compared them to all available protein structures. Multiple cases of recruitment of host enzymes for structural roles in viruses, accompanied by the disruption of catalytic sites, were discovered. However, many viral proteins appear to have evolved unique structural folds.