Off-Target V(D)J Recombination Drives Lymphomagenesis and Is Escalated by Loss of the Rag2 C Terminus

Diversification of Antibodies: From V(D)J Recombination to Somatic Exon Shuffling

Antibodies that gain specificity by a large insert encoding for an extra domain were described for the first time in 2016. In malaria-exposed individuals, an exon deriving from the leukocyte-associated immunoglobulin-like 1 (LAIR1) gene integrated via a copy-and-paste insertion into the immunoglobulin heavy chain encoding region. A few years later, a second example was identified, namely a dual exon integration from the leukocyte immunoglobulin-like receptor B1 (LILRB1) gene that is located in close proximity to LAIR1. A dedicated high-throughput characterization of chimeric immunoglobulin heavy chain transcripts unraveled, that insertions from distant genomic regions (including mitochondrial DNA) can contribute to human antibody diversity. This review describes the modalities of insert-containing antibodies. The role of known DNA mobility aspects, such as genomic translocation, gene conversion, and DNA fragility, is discussed in the context of insert-antibody generation. Finally, the review covers why insert antibodies were omitted from the past repertoire analyses and how insert antibodies can contribute to protective immunity or an autoreactive response.

RAG1 and RAG2 non-core regions are implicated in leukemogenesis and off-target V(D)J recombination in BCR-ABL1-driven B-cell lineage lymphoblastic leukemia

The evolutionary conservation of non-core RAG regions suggests significant roles that might involve quantitative or qualitative alterations in RAG activity. Off-target V(D)J recombination contributes to lymphomagenesis and is exacerbated by RAG2' C-terminus absence in Tp53-/- mice thymic lymphomas. However, the genomic stability effects of non-core regions from both Rag1c/c and Rag2c/c in BCR-ABL1+ B-lymphoblastic leukemia (BCR-ABL1+ B-ALL), the characteristics, and mechanisms of non-core regions in suppressing off-target V(D)J recombination remain unclear. Here, we established three mouse models of BCR-ABL1+ B-ALL in mice expressing full-length RAG (Ragf/f), core RAG1 (Rag1c/c), and core RAG2 (Rag2c/c). The Ragc/c (Rag1c/c and Rag2c/c) leukemia cells exhibited greater malignant tumor characteristics compared to Ragf/f cells. Additionally, Ragc/c cells showed higher frequency of off-target V(D)J recombination and oncogenic mutations than Ragf/f. We also revealed decreased RAG cleavage accuracy in Ragc/c cells and a smaller recombinant size in Rag1c/c cells, which could potentially exacerbate off-target V(D)J recombination in Ragc/c cells. In conclusion, these findings indicate that the non-core RAG regions, particularly the non-core region of RAG1, play a significant role in preserving V(D)J recombination precision and genomic stability in BCR-ABL1+ B-ALL.

RAG genomic variation causes autoimmune diseases through specific structure-based mechanisms of enzyme dysregulation

Interpreting genetic changes observed in individual patients is a critical challenge. The array of immune deficiency syndromes is typically caused by genetic variation unique to individuals. Therefore, new approaches are needed to interpret functional variation and accelerate genomics interpretation. We constructed the first full-length structural model of human RAG recombinase across four functional states of the recombination process. We functionally tested 182 clinically observed RAG missense mutations. These experiments revealed dysfunction due to recombinase dysfunction and altered chromatin interactions. Structural modeling identified mechanical and energetic roles for each mutation. We built regression models for RAG1 (R2 = 0.91) and RAG2 (R2 = 0.97) to predict RAG activity changes. We applied our model to 711 additional RAG variants observed in population studies and identified a subset that may impair RAG function. Thus, we demonstrated a fundamental advance in the mechanistic interpretation of human genetic variations spanning from rare and undiagnosed diseases to population health.

An updated definition of V(D)J recombination signal sequences revealed by high-throughput recombination assays

In the adaptive immune system, V(D)J recombination initiates the production of a diverse antigen receptor repertoire in developing B and T cells. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank antigen receptor gene segments. Previous studies defined the consensus RSS as containing conserved heptamer and nonamer sequences separated by a less conserved 12 or 23 base-pair spacer sequence. However, many RSSs deviate from the consensus sequence. Here, we developed a cell-based, massively parallel assay to evaluate V(D)J recombination activity on thousands of RSSs where the 12-RSS heptamer and adjoining spacer region contained randomized sequences. While the consensus heptamer sequence (CACAGTG) was marginally preferred, V(D)J recombination was highly active on a wide range of non-consensus sequences. Select purine/pyrimidine motifs that may accommodate heptamer unwinding in the RAG1/2 active site were generally preferred. In addition, while different coding flanks and nonamer sequences affected recombination efficiency, the relative dependency on the purine/pyrimidine motifs in the RSS heptamer remained unchanged. Our results suggest RAG1/2 specificity for RSS heptamers is primarily dictated by DNA structural features dependent on purine/pyrimidine pattern, and to a lesser extent, RAG:RSS base-specific interactions.

Different classes of genomic inserts contribute to human antibody diversity

The enormous diversity of antibodies is a key element to combat infections. Antibodies containing pathogen receptors were a surprising discovery that contrasted antibody diversification through classic recombination events. However, such insert-containing antibodies were thus far exclusively detected in African individuals exposed to malaria parasites and were identified as screening byproducts or through hypothesis-driven search. The prevalence and complexity of insertion events remained elusive. In this study, we devise an unbiased, systematic approach to identify inserts in the human antibody repertoire. We show that inserts from distant genomic regions occur in the majority of donors and are independent of Plasmodium falciparum preexposure. Our findings suggest that four distinct classes of insertion events contribute diversity to the human antibody repertoire.

Recombination of antibody genes in B cells can involve distant genomic loci and contribute a foreign antigen-binding element to form hybrid antibodies with broad reactivity for Plasmodium falciparum. So far, antibodies containing the extracellular domain of the LAIR1 and LILRB1 receptors represent unique examples of cross-chromosomal antibody diversification. Here, we devise a technique to profile non-VDJ elements from distant genes in antibody transcripts. Independent of the preexposure of donors to malaria parasites, non-VDJ inserts were detected in 80% of individuals at frequencies of 1 in 10⁴ to 10⁵ B cells. We detected insertions in heavy, but not in light chain or T cell receptor transcripts. We classify the insertions into four types depending on the insert origin and destination: 1) mitochondrial and 2) nuclear DNA inserts integrated at VDJ junctions; 3) inserts originating from telomere proximal genes; and 4) fragile sites incorporated between J-to-constant junctions. The latter class of inserts was exclusively found in memory and in in vitro activated B cells, while all other classes were already detected in naïve B cells. More than 10% of inserts preserved the reading frame, including transcripts with signs of antigen-driven affinity maturation. Collectively, our study unravels a mechanism of antibody diversification that is layered on the classical V(D)J and switch recombination.

V(D)J Recombination: Recent Insights in Formation of the Recombinase Complex and Recruitment of DNA Repair Machinery

V(D)J recombination is an essential mechanism of the adaptive immune system, producing a diverse set of antigen receptors in developing lymphocytes via regulated double strand DNA break and subsequent repair. DNA cleavage is initiated by the recombinase complex, consisting of lymphocyte specific proteins RAG1 and RAG2, while the repair phase is completed by classical non-homologous end joining (NHEJ). Many of the individual steps of this process have been well described and new research has increased the scale to understand the mechanisms of initiation and intermediate stages of the pathway. In this review we discuss 1) the regulatory functions of RAGs, 2) recruitment of RAGs to the site of recombination and formation of a paired complex, 3) the transition from a post-cleavage complex containing RAGs and cleaved DNA ends to the NHEJ repair phase, and 4) the potential redundant roles of certain factors in repairing the break. Regulatory (non-core) domains of RAGs are not necessary for catalytic activity, but likely influence recruitment and stabilization through interaction with modified histones and conformational changes. To form long range paired complexes, recent studies have found evidence in support of large scale chromosomal contraction through various factors to utilize diverse gene segments. Following the paired cleavage event, four broken DNA ends must now make a regulated transition to the repair phase, which can be controlled by dynamic conformational changes and post-translational modification of the factors involved. Additionally, we examine the overlapping roles of certain NHEJ factors which allows for prevention of genomic instability due to incomplete repair in the absence of one, but are lethal in combined knockouts. To conclude, we focus on the importance of understanding the detail of these processes in regards to off-target recombination or deficiency-mediated clinical manifestations.

The (Lack of) DNA Double-Strand Break Repair Pathway Choice During V(D)J Recombination

DNA double-strand breaks (DSBs) are highly toxic lesions that can be mended via several DNA repair pathways. Multiple factors can influence the choice and the restrictiveness of repair towards a given pathway in order to warrant the maintenance of genome integrity. During V(D)J recombination, RAG-induced DSBs are (almost) exclusively repaired by the non-homologous end-joining (NHEJ) pathway for the benefit of antigen receptor gene diversity. Here, we review the various parameters that constrain repair of RAG-generated DSBs to NHEJ, including the peculiarity of DNA DSB ends generated by the RAG nuclease, the establishment and maintenance of a post-cleavage synaptic complex, and the protection of DNA ends against resection and (micro)homology-directed repair. In this physiological context, we highlight that certain DSBs have limited DNA repair pathway choice options.

Transcriptomes of an Array of Chicken Ovary, Intestinal, and Immune Cells and Tissues

While the chicken (Gallus gallus) is the most consumed agricultural animal worldwide, the chicken transcriptome remains understudied. We have characterized the transcriptome of 10 cell and tissue types from the chicken using RNA-seq, spanning intestinal tissues (ileum, jejunum, proximal cecum), immune cells (B cells, bursa, macrophages, monocytes, spleen T cells, thymus), and reproductive tissue (ovary). We detected 17,872 genes and 24,812 transcripts across all cell and tissue types, representing 73% and 63% of the current gene annotation, respectively. Further quantification of RNA transcript biotypes revealed protein-coding and lncRNAs specific to an individual cell/tissue type. Each cell/tissue type also has an average of around 1.2 isoforms per gene, however, they all have at least one gene with at least 11 isoforms. Differential expression analysis revealed a large number of differentially expressed genes between tissues of the same category (immune and intestinal). Many of these differentially expressed genes in immune cells were involved in cellular processes relating to differentiation and cell metabolism as well as basic functions of immune cells such as cell adhesion and signal transduction. The differential expressed genes of the different segments of the chicken intestine (jejunum, ileum, proximal cecum) correlated to the metabolic processes in nutrient digestion and absorption. These data should provide a valuable resource in understanding the chicken genome.

Evolutionary Comparative Analyses of DNA-Editing Enzymes of the Immune System: From 5-Dimensional Description of Protein Structures to Immunological Insights and Applications to Protein Engineering

The immune system is unique among all biological sub-systems in its usage of DNA-editing enzymes to introduce targeted gene mutations and double-strand DNA breaks to diversify antigen receptor genes and combat viral infections. These processes, initiated by specific DNA-editing enzymes, often result in mistargeted induction of genome lesions that initiate and drive cancers. Like other molecules involved in human health and disease, the DNA-editing enzymes of the immune system have been intensively studied in humans and mice, with little attention paid (< 1% of published studies) to the same enzymes in evolutionarily distant species. Here, we present a systematic review of the literature on the characterization of one such DNA-editing enzyme, activation-induced cytidine deaminase (AID), from an evolutionary comparative perspective. The central thesis of this review is that although the evolutionary comparative approach represents a minuscule fraction of published works on this and other DNA-editing enzymes, this approach has made significant impacts across the fields of structural biology, immunology, and cancer research. Using AID as an example, we highlight the value of the evolutionary comparative approach in discoveries already made, and in the context of emerging directions in immunology and protein engineering. We introduce the concept of 5-dimensional (5D) description of protein structures, a more nuanced view of a structure that is made possible by evolutionary comparative studies. In this higher dimensional view of a protein's structure, the classical 3-dimensional (3D) structure is integrated in the context of real-time conformations and evolutionary time shifts (4th dimension) and the relevance of these dynamics to its biological function (5th dimension).

Breakage-Fusion-Bridge Events Trigger Complex Genome Rearrangements and Amplifications in Developmentally Arrested T Cell Lymphomas

To reveal the relative contribution of the recombination activating gene (RAG)1/2 nuclease to lymphomagenesis, we conducted a genome-wide analysis of T cell lymphomas from p53-deficient mice expressing or lacking RAG2. We found that while p53^−/− lymphoblastic T cells harbor primarily ectopic DNA deletions, Rag2^−/−p53^−/− T cell lymphomas display complex genomic rearrangements associated with amplification of the chromosomal location 9qA4-5.3. We show that this amplicon is generated by breakage-fusion-bridge during mitosis and arises distinctly in T cell lymphomas originating from an early progenitor stage. Notably, we report amplification of the corresponding syntenic region (11q23) in a subset of human leukemia leading to the overexpression of several cancer genes, including MLL/KMT2A. Our findings provide direct evidence that lymphocytes undergo malignant transformation through distinct genome architectural routes that are determined by both RAG-dependent and RAG-independent DNA damage and a block in cell development.

•

Lymphomas from RAG2/p53- and p53-deficient mice bear distinct genome architectures

•

Block in T cell development leads to 9qA4-5.3 rearrangements and amplifications

•

Breakage-fusion-bridge events trigger 9qA4-5.3 aberrations in early T cell lymphomas

•

The syntenic region 11q23 is amplified in some human hematological cancers

Alterations in chromatin at antigen receptor loci define lineage progression during B lymphopoiesis

Developing lymphocytes diversify their antigen receptor (AgR) loci by variable (diversity) joining (V[D]J) recombination. Here, using the micrococcal nuclease (MNase)-based chromatin accessibility (MACC) assay with low-cell count input, we profile both small-scale (kilobase) and large-scale (megabase) changes in chromatin accessibility and nucleosome occupancy in primary cells during lymphoid development, tracking the changes as different AgR loci become primed for recombination. The three distinct chromatin structures identified in this work define unique features of immunoglobulin H (IgH), Igκ, and T cell receptor-α (TCRα) loci during B lymphopoiesis. In particular, we find locus-specific temporal changes in accessibility both across megabase-long AgR loci and locally at the recombination signal sequences (RSSs). These changes seem to be regulated independently and can occur prior to lineage commitment. Large-scale changes in chromatin accessibility occur without significant change in nucleosome density and represent key features of AgR loci not previously described. We further identify local dynamic repositioning of individual RSS-associated nucleosomes at IgH and Igκ loci while they become primed for recombination during B cell commitment. These changes in chromatin at AgR loci are regulated in a locus-, lineage-, and stage-specific manner during B lymphopoiesis, serving either to facilitate or to impose a barrier to V(D)J recombination. We suggest that local and global changes in chromatin openness in concert with nucleosome occupancy and placement of histone modifications facilitate the temporal order of AgR recombination. Our data have implications for the organizing principles that govern assembly of these large loci as well as for mechanisms that might contribute to aberrant V(D)J recombination and the development of lymphoid tumors.

RAG1 high expression associated with IKZF1 dysfunction in adult B-cell acute lymphoblastic leukemia

The recombination mediated by recombination activating gene (RAG) is not only the dominant mutational process but also the predominant driver of oncogenic genomic rearrangement in acute lymphoblastic leukemia (ALL). It is further responsible for leukemic clonal evolution. In this study, significant RAG1 increase is observed in the subsets of B-ALL patients, and high expression of RAG1 is observed to be correlated with high proliferation markers. IKZF1-encoded protein, IKAROS, directly binds to the RAG1 promoter and regulates RAG1 expression in leukemic cells. CK2 inhibitor by increasing IKAROS activity significantly suppresses RAG1 expression in ALL in an IKAROS-dependent manner. Patients with IKZF1 deletion have significantly higher expression of RAG1 compared to that without IKZF1 deletion. CK2 inhibitor treatment also results in an increase in IKZF1 binding to the RAG1 promoter and suppression of RAG1 expression in primary ALL cells. Taken together, these results demonstrate that RAG1 high expression is associated with high proliferation markers in B-ALL. Our data for the first time proved that RAG1 expression is directly suppressed by IKAROS. Our results also reveal drive oncogenesis of B-ALL is driven by high expression of RAG1 with IKAROS dysfunction together, which have significance in an integrated prognostic model for adult ALL.

RAG gene defects at the verge of immunodeficiency and immune dysregulation

Mutations of the recombinase activating genes (RAG) in humans underlie a broad spectrum of clinical and immunological phenotypes that reflect different degrees of impairment of T- and B-cell development and alterations of mechanisms of central and peripheral tolerance. Recent studies have shown that this phenotypic heterogeneity correlates, albeit imperfectly, with different levels of recombination activity of the mutant RAG proteins. Furthermore, studies in patients and in newly developed animal models carrying hypomorphic RAG mutations have disclosed various mechanisms underlying immune dysregulation in this condition. Careful annotation of clinical outcome and immune reconstitution in RAG-deficient patients who have received hematopoietic stem cell transplantation has shown that progress has been made in the treatment of this disease, but new approaches remain to be tested to improve stem cell engraftment and durable immune reconstitution. Finally, initial attempts have been made to treat RAG deficiency with gene therapy.

The Conserved ATM Kinase RAG2-S365 Phosphorylation Site Limits Cleavage Events in Individual Cells Independent of Any Repair Defect

Many DNA lesions associated with lymphoid malignancies are linked to off-target cleavage by the RAG1/2 recombinase. However, off-target cleavage has mostly been analyzed in the context of DNA repair defects, confounding any mechanistic understanding of cleavage deregulation. We identified a conserved SQ phosphorylation site on RAG2 365 to 366 that is involved in feedback control of RAG cleavage. Mutation of serine 365 to a non-phosphorylatable alanine permits bi-allelic and bi-locus RAG-mediated breaks in the same cell, leading to reciprocal translocations. This phenomenon is analogous to the phenotype we described for ATM kinase inactivation. Here, we establish deregulated cleavage itself as a driver of chromosomal instability without the associated repair defect. Intriguingly, a RAG2-S365E phosphomimetic rescues the deregulated cleavage of ATM inactivation, reducing the incidence of reciprocal translocations. These data support a model in which feedback control of cleavage and maintenance of genome stability involves ATM-mediated phosphorylation of RAG2.

The IGH locus relocalizes to a "recombination compartment" in the perinucleolar region of differentiating B-lymphocytes

The immunoglobulin heavy chain (IGH) gene loci are subject to specific recombination events during B-cell differentiation including somatic hypermutation and class switch recombination which mark the end of immunoglobulin gene maturation in germinal centers of secondary lymph nodes. These two events rely on the activity of activation-induced cytidine deaminase (AID) which requires DNA double strand breaks be created, a potential danger to the cell. Applying 3D-fluorescence in situ hybridization coupled with immunofluorescence staining to a previously described experimental system recapitulating normal B-cell differentiation ex vivo, we have kinetically analyzed the radial positioning of the two IGH gene loci as well as their proximity with the nucleolus, heterochromatin and γH2AX foci. Our observations are consistent with the proposal that these IGH gene rearrangements take place in a specific perinucleolar “recombination compartment” where AID could be sequestered thus limiting the extent of its potentially deleterious off-target effects.

The miR-15 family reinforces the transition from proliferation to differentiation in pre-B cells

Precursor B lymphocytes expand upon expression of a pre-B cell receptor (pre-BCR), but then transit into a resting state in which immunoglobulin light chain gene recombination is initiated. This bi-phasic sequence is orchestrated by the IL-7 receptor (IL-7R) and pre-BCR signaling, respectively, but little is known about microRNAs fine-tuning these events. Here, we show that pre-B cells lacking miR-15 family functions exhibit prolonged proliferation due to aberrant expression of the target genes cyclin E1 and D3. As a consequence, they fail to trigger the transcriptional reprogramming normally accompanying their differentiation, resulting in a developmental block at the pre-B cell stage. Intriguingly, our data indicate that the miR-15 family is suppressed by both IL-7R and pre-BCR signaling, suggesting it is actively integrated into the regulatory circuits of developing B cells. These findings identify the miR-15 family as a novel element required to promote the switch from pre-B cell proliferation to differentiation.

Cell cycle RNA regulons coordinating early lymphocyte development

Lymphocytes undergo dynamic changes in gene expression as they develop from progenitor cells lacking antigen receptors, to mature cells that are prepared to mount immune responses. While transcription factors have established roles in lymphocyte development, they act in concert with post-transcriptional and post-translational regulators to determine the proteome. Furthermore, the post-transcriptional regulation of RNA regulons consisting of mRNAs whose protein products act cooperatively allows RNA binding proteins to exert their effects at multiple points in a pathway. Here, we review recent evidence demonstrating the importance of RNA binding proteins that control the cell cycle in lymphocyte development and discuss the implications for tumorigenesis. WIREs RNA 2017, 8:e1419. doi: 10.1002/wrna.1419 For further resources related to this article, please visit the WIREs website.

RAG1/2 induces genomic insertions by mobilizing DNA into RAG1/2-independent breaks

Rommel et al. reveal a novel RAG1/2-mediated insertion pathway, which has the potential to destabilize the lymphocyte genome and shares features with DNA insertions observed in human cancer.

The RAG recombinase (RAG1/2) plays an essential role in adaptive immunity by mediating V(D)J recombination in developing lymphocytes. In contrast, aberrant RAG1/2 activity promotes lymphocyte malignancies by causing chromosomal translocations and DNA deletions at cancer genes. RAG1/2 can also induce genomic DNA insertions by transposition and trans-V(D)J recombination, but only few such putative events have been documented in vivo. We used next-generation sequencing techniques to examine chromosomal rearrangements in primary murine B cells and discovered that RAG1/2 causes aberrant insertions by releasing cleaved antibody gene fragments that subsequently reintegrate into DNA breaks induced on a heterologous chromosome. We confirmed that RAG1/2 also mobilizes genomic DNA into independent physiological breaks by identifying similar insertions in human lymphoma and leukemia. Our findings reveal a novel RAG1/2-mediated insertion pathway distinct from DNA transposition and trans-V(D)J recombination that destabilizes the genome and shares features with reported oncogenic DNA insertions.

RAG Recombinase as a Selective Pressure for Genome Evolution

The RAG recombinase is a domesticated transposable element co-opted in jawed vertebrates to drive the process of the so-called V(D)J recombination, which is the hallmark of the adaptive immune system to produce antigen receptors. RAG targets, namely, the Recombination Signal Sequences (RSS), are rather long and degenerated sequences, which highlights the ability of the recombinase to interact with a wide range of target sequences, including outside of antigen receptor loci. The recognition of such cryptic targets by the recombinase threatens genome integrity by promoting aberrant DNA recombination, as observed in lymphoid malignancies. Genomes evolution resulting from RAG acquisition is an ongoing discussion, in particular regarding the counter-selection of sequences resembling the RSS and the modifications of epigenetic regulation at these potential cryptic sites. Here, we describe a new bioinformatics tool to map potential RAG targets in all jawed vertebrates. We show that our REcombination Classifier (REC) outperforms the currently available tool and is suitable for full genomes scans from species other than human and mouse. Using the REC, we document a reduction in density of potential RAG targets at the transcription start sites of genes co-expressed with the rag genes and marked with high levels of the trimethylation of the lysine 4 of the histone 3 (H3K4me3), which correlates with the retention of functional RAG activity after the horizontal transfer.

The RAG recombinase: Beyond breaking

DNA double-strand breaks (DSBs) are commonly seen as lesions that threaten genome integrity and contribute to cancer and aging processes. However, in the context of antigen receptor gene assembly, known as V(D)J recombination, DSBs are obligatory intermediates that allow the establishment of genetic diversity and adaptive immunity. V(D)J recombination is initiated when the lymphoid-restricted recombination-activating genes RAG1 and RAG2 are expressed and form a site-specific endonuclease (the RAG nuclease or RAG recombinase). Here, we discuss the ability of the RAG nuclease to minimize the risks of genome disruption by coupling the breakage and repair steps of the V(D)J reaction. This implies that the RAG genes, derived from an ancient transposon, have undergone strong selective pressure to prohibit transposition in favor of promoting controlled DNA end joining in cis by the ubiquitous DNA damage response and DNA repair machineries. We also discuss the idea that, in addition to being essential for the rearrangement of antigen receptor genes, RAG-mediated DSBs could impact cellular processes and outcomes by affecting genetic and epigenetic programs.

DNA Breaks and End Resection Measured Genome-wide by End Sequencing

DNA double-strand breaks (DSBs) arise during physiological transcription, DNA replication, and antigen receptor diversification. Mistargeting or misprocessing of DSBs can result in pathological structural variation and mutation. Here we describe a sensitive method (END-seq) to monitor DNA end resection and DSBs genome-wide at base-pair resolution in vivo. We utilized END-seq to determine the frequency and spectrum of restriction-enzyme-, zinc-finger-nuclease-, and RAG-induced DSBs. Beyond sequence preference, chromatin features dictate the repertoire of these genome-modifying enzymes. END-seq can detect at least one DSB per cell among 10,000 cells not harboring DSBs, and we estimate that up to one out of 60 cells contains off-target RAG cleavage. In addition to site-specific cleavage, we detect DSBs distributed over extended regions during immunoglobulin class-switch recombination. Thus, END-seq provides a snapshot of DNA ends genome-wide, which can be utilized for understanding genome-editing specificities and the influence of chromatin on DSB pathway choice.

Illegitimate RAG-mediated recombination events are involved in IKZF1 Δ3-6 deletion in BCR-ABL1 lymphoblastic leukaemia

Breakpoint cluster region-Abelson murine leukaemia viral oncogene homologue 1 (BCR-ABL1), encoded by the Philadelphia (Ph) chromosome, is the characteristic of chronic myeloid leukaemia (CML) and a subset of acute lymphoblastic leukaemia (ALL). We demonstrated that expression of the Ik6 transcript, which lacked exons 3-6, was observed exclusively in BCR-ABL1(+) B ALL and lymphoid blast crisis CML (BC-CML) patients harbouring the IKZF1 Δ3-6 deletion. To confirm the hypothesis that illegitimate recombination activating gene protein (RAG)-mediated recombination events are involved in IKZF1 Δ3-6 deletion in BCR-ABL1 lymphoblastic leukaemia, we first demonstrated that the expression rates of RAG1 and RAG2, collectively called RAG, were higher in ALL and BC-CML (lymphoid). Notably, analysis of relationships among RAG, BCR-ABL1 and Ikaros 6 (Ik6) showed that Ik6 can be generated only if RAG and BCR-ABL1 are co-existing. The sequencing data showed that the deleted segments of introns 2 and 6 contained cryptic recombination signal sequences (cRSSs) and frequently had non-template nucleotides inserted between breakpoints. Furthermore, we used chromatin immunoprecipitation (ChIP) technology and demonstrated that the sequences directly flanking IKZF1 Δ3-6 deletion breakpoints have significantly higher levels of histone H3 lysine 4 trimethylation (H3K4me3) modifications. Overall, RAG expression, good-quality cRSS and a specific chromatin modification, H3K4me3, satisfy the conditions of RAG's off-target effects on IKZF1. Our work provides evidence for RAG-mediated IKZF1 Δ3-6 deletion. Our results raise the prospect that RAG is a valuable biomarker in disease surveillance. Dissecting the contribution of RAG should not only provide valuable mechanistic insights, but will also lead to a new therapeutic direction.

RAG1 targeting in the genome is dominated by chromatin interactions mediated by the non-core regions of RAG1 and RAG2

The RAG1/RAG2 endonuclease initiates V(D)J recombination at antigen receptor loci but also binds to thousands of places outside of these loci. RAG2 localizes directly to lysine 4 trimethylated histone 3 (H3K4me3) through a plant homeodomain (PHD) finger. The relative contribution of RAG2-dependent and RAG1-intrinsic mechanisms in determining RAG1 binding patterns is not known. Through analysis of deep RAG1 ChIP-seq data, we provide a quantitative description of the forces underlying genome-wide targeting of RAG1. Surprisingly, sequence-specific DNA binding contributes minimally to RAG1 targeting outside of antigen receptor loci. Instead, RAG1 binding is driven by two distinct modes of interaction with chromatin: the first is driven by H3K4me3, promoter-focused and dependent on the RAG2 PHD, and the second is defined by H3K27Ac, enhancer-focused and dependent on ‘non-core’ portions of RAG1. Based on this and additional chromatin and genomic features, we formulated a predictive model of RAG1 targeting to the genome. RAG1 binding sites predicted by our model correlate well with observed patterns of RAG1-mediated breaks in human pro-B acute lymphoblastic leukemia. Overall, this study provides an integrative model for RAG1 genome-wide binding and off-target activity and reveals a novel role for the RAG1 non-core region in RAG1 targeting.