RSSsite: a reference database and prediction tool for the identification of cryptic Recombination Signal Sequences in human and murine genomes

The biological and clinical impact of deletions before and after large chromosomal gains in multiple myeloma

ABSTRACT

Acquisition of a hyperdiploid (HY) karyotype or immunoglobulin heavy chain (IgH) translocations are considered key initiating events in multiple myeloma (MM). To explore if other genomic events can precede these events, we analyzed whole-genome sequencing data from 1173 MM samples. By integrating molecular time and structural variants within early chromosomal duplications, we indeed identified pregain deletions in 9.4% of patients with an HY karyotype without IgH translocations, challenging acquisition of an HY karyotype as the earliest somatic event. Remarkably, these deletions affected tumor suppressor genes (TSGs) and/or oncogenes in 2.4% of patients with an HY karyotype without IgH translocations, supporting their role in MM pathogenesis. Furthermore, our study points to postgain deletions as novel driver mechanisms in MM. Using multiomics approaches to investigate their biologic impact, we found associations with poor clinical outcome in newly diagnosed patients and profound effects on both the oncogene and TSG activity despite the diploid gene status. Overall, this study provides novel insights into the temporal dynamics of genomic alterations in MM.

Molecular basis for differential Igk versus Igh V(D)J joining mechanisms

In developing B cells, V(D)J recombination assembles exons encoding IgH and Igκ variable regions from hundreds of gene segments clustered across Igh and Igk loci. V, D and J gene segments are flanked by conserved recombination signal sequences (RSSs) that target RAG endonuclease1. RAG orchestrates Igh V(D)J recombination upon capturing a JH-RSS within the JH-RSS-based recombination centre1-3 (RC). JH-RSS orientation programmes RAG to scan upstream D- and VH-containing chromatin that is presented in a linear manner by cohesin-mediated loop extrusion4-7. During Igh scanning, RAG robustly utilizes only D-RSSs or VH-RSSs in convergent (deletional) orientation with JH-RSSs4-7. However, for Vκ-to-Jκ joining, RAG utilizes Vκ-RSSs from deletional- and inversional-oriented clusters8, inconsistent with linear scanning2. Here we characterize the Vκ-to-Jκ joining mechanism. Igk undergoes robust primary and secondary rearrangements9,10, which confounds scanning assays. We therefore engineered cells to undergo only primary Vκ-to-Jκ rearrangements and found that RAG scanning from the primary Jκ-RC terminates just 8 kb upstream within the CTCF-site-based Sis element11. Whereas Sis and the Jκ-RC barely interacted with the Vκ locus, the CTCF-site-based Cer element12 4 kb upstream of Sis interacted with various loop extrusion impediments across the locus. Similar to VH locus inversion7, DJH inversion abrogated VH-to-DJH joining; yet Vκ locus or Jκ inversion allowed robust Vκ-to-Jκ joining. Together, these experiments implicated loop extrusion in bringing Vκ segments near Cer for short-range diffusion-mediated capture by RC-based RAG. To identify key mechanistic elements for diffusional V(D)J recombination in Igk versus Igh, we assayed Vκ-to-JH and D-to-Jκ rearrangements in hybrid Igh-Igk loci generated by targeted chromosomal translocations, and pinpointed remarkably strong Vκ and Jκ RSSs. Indeed, RSS replacements in hybrid or normal Igk and Igh loci confirmed the ability of Igk-RSSs to promote robust diffusional joining compared with Igh-RSSs. We propose that Igk evolved strong RSSs to mediate diffusional Vκ-to-Jκ joining, whereas Igh evolved weaker RSSs requisite for modulating VH joining by RAG-scanning impediments.

Unraveling the chicken T cell repertoire with enhanced genome annotation

T cell receptor (TCR) repertoire sequencing has emerged as a powerful tool for understanding the diversity and functionality of T cells within the host immune system. Yet, the chicken TCR repertoire remains poorly understood due to incomplete genome annotation of the TCR loci, despite the importance of chickens in agriculture and as an immunological model. Here, we addressed this critical issue by employing 5' rapid amplification of complementary DNA ends (5'RACE) TCR repertoire sequencing with molecular barcoding of complementary DNA (cDNA) molecules. Simultaneously, we enhanced the genome annotation of TCR Variable (V), Diversity (D, only present in β and δ loci) and Joining (J) genes in the chicken genome. To enhance the efficiency of TCR annotations, we developed VJ-gene-finder, an algorithm designed to extract VJ gene candidates from deoxyribonucleic acid (DNA) sequences. Using this tool, we achieved a comprehensive annotation of all known chicken TCR loci, including the α/δ locus on chromosome 27. Evolutionary analysis revealed that each locus evolved separately by duplication of long homology units. To define the baseline TCR diversity in healthy chickens and to demonstrate the feasibility of the approach, we characterized the splenic α/β/γ/δ TCR repertoire. Analysis of the repertoires revealed preferential usage of specific V and J combinations in all chains, while the overall features were characteristic of unbiased repertoires. We observed moderate levels of shared complementarity-determining region 3 (CDR3) clonotypes among individual birds within the α and γ chain repertoires, including the most frequently occurring clonotypes. However, the β and δ repertoires were predominantly unique to each bird. Taken together, our TCR repertoire analysis allowed us to decipher the composition, diversity, and functionality of T cells in chickens. This work not only represents a significant step towards understanding avian T cell biology, but will also shed light on host-pathogen interactions, vaccine development, and the evolutionary history of avian immunology.

Macroevolution of avian T cell receptor C segments using genomic data

All jawed vertebrates have four T cell receptor (TCR) chains expressed by thymus-derived lymphocytes that play a significant role in animal immune defense. However, avian TCR studies have been limited to a few species, although their co-functional major histocompatibility complexes (MHCs) have been studied for decades, showing various copy numbers and polymorphisms. Here, using public genome data, we characterized the copy numbers, the phylogenic relationship and selection of T cell receptor complex (TCR-C) segments, and the genomic organization of TCR loci across birds. Various numbers of C segments were found in the TCRα/TCRδ, TCRβ, and TCRγ loci, and phylogenetic analysis reflected both ancient gene duplication events (two Cβ segments and Cδ segments divergent into CδI and CδII) and contemporary evolution (lineage-specific and species-specific characteristics). Most passerines lack CδII segments and a second TRD locus, except Meliphagidae and Maluridae. A relatively stable structure was verified in four TCR loci of birds, except for the arrangement of V segment groups. In this study, we explored the phylogenetic relationships of TCR-C segments across avians for the first time. We inferred gene duplication and loss events during the evolution process. The finding of diverse TCR germline repertoires provides a better understanding of the immune systems of birds.

Molecular Mechanisms of DNA Sequence Selectivity in V(D)J Recombination

Antigen receptor (AgR) diversity is central to the ability of adaptive immunity in jawed vertebrates to protect against pathogenic agents. The production of highly diverse AgR repertoires is initiated during B and T cell lymphopoiesis by V(D)J recombination, which assembles the receptor genes from component gene segments in a cut-and-paste recombination reaction. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank AgR 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, and the molecular mechanism for semiselective V(D)J recombination specificity is unknown. The modulation of chromatin structure during V(D)J recombination is essential in the formation of diverse AgRs in adaptive immunity while also reducing the likelihood for off-target recombination events that can result in chromosomal aberrations and genomic instability. Here we review what is presently known regarding mechanisms that facilitate assembly of RAG1/2 with RSSs, the ensuing conformational changes required for DNA cleavage activity, and how the readout of the RSS sequence affects reaction efficiency.

Molecular characterization of TCF3::PBX1 chromosomal breakpoints in acute lymphoblastic leukemia and their use for measurable residual disease assessment

The translocation t(1;19)(q23;p13) with the resulting chimeric TCF3::PBX1 gene is the third most prevalent recurrent chromosomal translocation in acute lymphoblastic leukemia and accounts for 3-5% of cases. The molecular background of this translocation has been incompletely studied, especially in adult cases. We characterized the chromosomal breakpoints of 49 patients with TCF3::PBX1 and the corresponding reciprocal PBX1::TCF3 breakpoints in 15 cases at the molecular level, thus providing an extensive molecular overview of this translocation in a well-defined study patient population. Breakpoints were found to be remarkably clustered not only in TCF3 but also in PBX1. No association with DNA repeats or putative cryptic recombination signal sequence sites was observed. A simplified detection method for breakpoint identification was developed and the feasibility of patient-specific chromosomal break sites as molecular markers for detecting measurable residual disease (MRD) was explored. A highly sensitive generic real-time PCR for MRD assessment using these breakpoint sequences was established that could serve as a useful alternative to the classical method utilizing rearranged immune gene loci. This study provides the first extensive molecular data set on the chromosomal breakpoints of the t(1;19)/TCF3::PBX1 aberration in adult ALL. Based on the obtained data a generic MRD method was developed that has several theoretical advantages, including an on average higher sensitivity and a greater stability of the molecular marker in the course of disease.

The Recurrent Atypical e8a2 BCR::ABL1 Transcript with Insertion of an Inverted 55 Base Pair ABL1 Intron 1b Sequence: A Detailed Molecular Analysis

Atypical BCR::ABL1 transcripts are found in approximately 2% of cases of chronic myeloid leukemia. It is important to detect them since affected patients also benefit from tyrosine kinase inhibitor therapy. In the rare e8a2 atypical BCR::ABL1 transcript, two out-of-frame exons are fused, thus, interposed nucleotides are usually found at the fusion site to restore the reading frame. In approximately half of previously reported e8a2 BCR::ABL1 cases, an inserted 55 bp sequence homologous to an inverted sequence from ABL1 intron 1b was detected. The generation of this recurrent transcript variant is not obvious. This work describes the molecular analysis of such an e8a2 BCR::ABL1 translocation from a CML patient. The genomic chromosomal breakpoint is identified, and the formation of this transcript is theoretically explained. The clinical course of the patient is reported, and recommendations are provided for the molecular analysis of future e8a2 BCR::ABL1 cases.

Intrathymic AAV delivery results in therapeutic site-specific integration at TCR loci in mice

Adeno-associated virus (AAV) vectors have been successfully exploited in gene therapy applications for the treatment of several genetic disorders. AAV is considered an episomal vector, but it has been shown to integrate within the host cell genome after the generation of double-strand DNA breaks or nicks. Although AAV integration raises some safety concerns, it can also provide therapeutic benefit; the direct intrathymic injection of an AAV harboring a therapeutic transgene results in integration in T-cell progenitors and long-term T-cell immunity. To assess the mechanisms of AAV integration, we retrieved and analyzed hundreds of AAV integration sites from lymph node-derived mature T cells and compared these with liver and brain tissue from treated mice. Notably, we found that although AAV integrations in the liver and brain were distributed across the entire mouse genome, >90% of the integrations in T cells were clustered within the T-cell receptor α, β, and γ genes. More precisely, the insertion mapped to DNA breaks created by the enzymatic activity of recombination activating genes (RAGs) during variable, diversity, and joining recombination. Our data indicate that RAG activity during T-cell receptor maturation induces a site-specific integration of AAV genomes and opens new therapeutic avenues for achieving long-term AAV-mediated gene transfer in dividing cells.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis.

SIGNIFICANCE

Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

New insights into the germline genes and CDR3 repertoire of the TCRβ chain in Chiroptera

Introduction

Bats are recognized as natural reservoirs for many viruses, and their unique immune system enables them to coexist with these viruses without frequently exhibiting disease symptoms. However, the current understanding of the bat adaptive immune system is limited due to the lack of a database or tool capable of processing T-cell receptor (TCR) sequences for bats.

Methods

We performed germline gene annotation in three bat species using homologous genes and RSSs (Recombinational Signal Sequences) scanning method. Then we used the conserved C gene to construct the TCRβ chain receptor library of the Intermediate Horseshoe Bat. Bats' TCRβ data will be analyzed using MiXCR and constructed reference library.

Results

Regarding the annotation results, we found that the Pale Spear-nosed Bat has 37 members in the TRBV12 family, which is more than the total number of TRBV genes in the Greater Horseshoe Bat. The average number of unique TCRβ chain receptor sequences in each Intermediate Horseshoe Bat sample reached 24,904.

Discussion

The distinct variations in the distribution of TRBV genes among the three types of bats could have a direct impact on the diversity of the TCR repertoire, as evidenced by the presence of conserved amino acids that indicate the T-cell recognition of antigens in bats is MHC-restricted. The bats’ TCRβ repertoire is formed through the rearrangement of the V-D-J-C genes, with D-J/V-D deletions and insertions resulting in high diversity.

Enhancing diagnosis of T-cell lymphoma using non-recombined T-cell receptor sequences

Clonality assessment, which can detect neoplastic T cells by identifying the uniquely recombined T-cell receptor (TCR) genes, provides important support in the diagnosis of T-cell lymphoma (TCL). BIOMED-2 is the gold standard clonality assay and has proven to be effective in European TCL patients. However, we failed to prove its sensitivity in Taiwanese TCL patients, especially based on the TCRβ gene. To explore potential impact of genetic background in the BIOMED-2 test, we analyzed TCRβ sequences of 21 healthy individuals and two TCL patients. This analysis suggests that genetic variations in the BIOMED-2 primer sites could not explain the difference in sensitivity. The BIOMED-2 test results of the two TCL patients were positive and negative, respectively. Interestingly, a higher percentage (>81%) of non-recombined TCRβ sequences was observed in the test-negative patient than those of the test-positive patient and all healthy individuals (13~66%). The result suggests a new TCR target for enhancing TCL diagnosis. To further explore the hypothesis, we proposed a cost-effective digital PCR assay that quantifies the relative abundance of non-recombined TCRβ sequences containing a J2-2P~J2-3 segment. With the digital PCR assay, bone marrow specimens from TCL patients (n=9) showed a positive outcome (i.e., the relative abundance of the J2-2P~J2-3 sequences ≧5%), whereas non-TCL patients (n=6) gave a negative result. As five of nine TCL patients had a negative BIOMED-2 test result, the J2-2P~J2-3 sequences may improve TCL detection. This is the first report showing the capability of characterizing non-recombined TCR sequences as a supplementary strategy for the BIOMED-2 clonality test.

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.

Tcf-1 promotes genomic instability and T cell transformation in response to aberrant β-catenin activation

Understanding the mechanisms promoting chromosomal translocations of the rearranging receptor loci in leukemia and lymphoma remains incomplete. Here we show that leukemias induced by aberrant activation of β-catenin in thymocytes, which bear recurrent Tcra/Myc-Pvt1 translocations, depend on Tcf-1. The DNA double strand breaks (DSBs) in the Tcra site of the translocation are Rag-generated, whereas the Myc-Pvt1 DSBs are not. Aberrantly activated β-catenin redirects Tcf-1 binding to novel DNA sites to alter chromatin accessibility and down-regulate genome-stability pathways. Impaired homologous recombination (HR) DNA repair and replication checkpoints lead to retention of DSBs that promote translocations and transformation of double-positive (DP) thymocytes. The resulting lymphomas, which resemble human T cell acute lymphoblastic leukemia (T-ALL), are sensitive to PARP inhibitors (PARPis). Our findings indicate that aberrant β-catenin signaling contributes to translocations in thymocytes by guiding Tcf-1 to promote the generation and retention of replication-induced DSBs allowing their coexistence with Rag-generated DSBs. Thus, PARPis could offer therapeutic options in hematologic malignancies with active Wnt/β-catenin signaling.

Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia

Transcriptome sequencing has identified multiple subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL) of prognostic significance, but a minority of cases lack a known genetic driver. Here, we used integrated whole-genome (WGS) and -transcriptome sequencing (RNA-seq), enhancer mapping, and chromatin topology analysis to identify previously unrecognized genomic drivers in B-ALL. Newly diagnosed (n = 3221) and relapsed (n = 177) B-ALL cases with tumor RNA-seq were studied. WGS was performed to detect mutations, structural variants, and copy number alterations. Integrated analysis of histone 3 lysine 27 acetylation and chromatin looping was performed using HiChIP. We identified a subset of 17 newly diagnosed and 5 relapsed B-ALL cases with a distinct gene expression profile and 2 universal and unique genomic alterations resulting from aberrant recombination-activating gene activation: a focal deletion downstream of PAN3 at 13q12.2 resulting in CDX2 deregulation by the PAN3 enhancer and a focal deletion of exons 18-21 of UBTF at 17q21.31 resulting in a chimeric fusion, UBTF::ATXN7L3. A subset of cases also had rearrangement and increased expression of the PAX5 gene, which is otherwise uncommon in B-ALL. Patients were more commonly female and young adult with median age 35 (range,12-70 years). The immunophenotype was characterized by CD10 negativity and immunoglobulin M positivity. Among 16 patients with known clinical response, 9 (56.3%) had high-risk features including relapse (n = 4) or minimal residual disease >1% at the end of remission induction (n = 5). CDX2-deregulated, UBTF::ATXN7L3 rearranged (CDX2/UBTF) B-ALL is a high-risk subtype of leukemia in young adults for which novel therapeutic approaches are required.

Receptor editing constrains development of phosphatidyl choline-specific B cells in VH12-transgenic mice

B1 B cells reactive to phosphatidyl choline (PtC) exhibit restricted immunoglobulin heavy chain (HC) and light chain (LC) combinations, exemplified by V_H12/Vκ4/5H. Two checkpoints are thought to focus PtC⁺ B cell maturation in V_H12-transgenic mice (VH12 mice): V-J rearrangements encoding a "permissive" LC capable of V_H12 HC pairing are selected first, followed by positive selection based on PtC binding, often requiring LC receptor editing to salvage PtC^- B cells and acquire PtC reactivity. However, evidence obtained from breeding VH12 mice to editing-defective dnRAG1 mice and analyzing LC sequences from PtC⁺ and PtC^- B cell subsets instead suggests that receptor editing functions after initial positive selection to remove PtC⁺ B cells in VH12 mice. This offers a mechanism to constrain natural, polyreactive B cells to limit their frequency. Sequencing also reveals occasional in-frame hybrid LC genes, reminiscent of type 2 gene replacement, that, testing suggests, arise via a recombination-activating gene (RAG)-independent mechanism.

Exploring the Origin and Physiological Significance of DNA Double Strand Breaks in the Developing Neuroretina

Genetic mosaicism is an intriguing physiological feature of the mammalian brain that generates altered genetic information and provides cellular, and prospectively functional, diversity in a manner similar to that of the immune system. However, both its origin and its physiological significance remain poorly characterized. Most, if not all, cases of somatic mosaicism require prior generation and repair of DNA double strand breaks (DSBs). The relationship between DSB generation, neurogenesis, and early neuronal cell death revealed by our studies in the developing retina provides new perspectives on the different mechanisms that contribute to DNA rearrangements in the developing brain. Here, we speculate on the physiological significance of these findings.

Gene prediction in the immunoglobulin loci

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) is critical for studies of the adaptive immune system. However, since the standard gene finding algorithms are not suitable for predicting IG genes, they have been semi-manually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species indicating that these genes may be subjected to the same selective pressure, e.g., pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, in difference from all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Since cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

Repertoire analysis of γδ T cells in the chicken enables functional annotation of the genomic region revealing highly variable pan-tissue TCR gamma V gene usage as well as identifying public and private repertoires

BACKGROUND

Despite increasing interest in γδ T cells and their non-classical behaviour, most studies focus on animals with low numbers of circulating γδ T cells, such as mice and humans. Arguably, γδ T cell functions might be more prominent in chickens where these cells form a higher proportion of the circulatory T cell compartment. The TCR repertoire defines different subsets of γδ T cells, and such analysis is facilitated by well-annotated TCR loci. γδ T cells are considered at the cusp of innate and adaptive immunity but most functions have been identified in γδ low species. A deeper understanding of TCR repertoire biology in γδ high and γδ low animals is critical for defining the evolution of the function of γδ T cells. Repertoire dynamics will reveal populations that can be classified as innate-like or adaptive-like as well as those that straddle this definition.

RESULTS

Here, a recent discrepancy in the structure of the chicken TCR gamma locus is resolved, demonstrating that tandem duplication events have shaped the evolution of this locus. Importantly, repertoire sequencing revealed large differences in the usage of individual TRGV genes, a pattern conserved across multiple tissues, including thymus, spleen and the gut. A single TRGV gene, TRGV3.3, with a highly diverse private CDR3 repertoire dominated every tissue in all birds. TRGV usage patterns were partly explained by the TRGV-associated recombination signal sequences. Public CDR3 clonotypes represented varying proportions of the repertoire of TCRs utilising different TRGVs, with one TRGV dominated by super-public clones present in all birds.

CONCLUSIONS

The application of repertoire analysis enabled functional annotation of the TCRG locus in a species with a high circulating γδ phenotype. This revealed variable usage of TCRGV genes across multiple tissues, a pattern quite different to that found in γδ low species (human and mouse). Defining the repertoire biology of avian γδ T cells will be key to understanding the evolution and functional diversity of these enigmatic lymphocytes in an animal that is numerically more reliant on them. Practically, this will reveal novel ways in which these cells can be exploited to improve health in medical and veterinary contexts.

Enhancing B-Cell Malignancies-On Repurposing Enhancer Activity towards Cancer

B-cell lymphomas and leukemias derive from B cells at various stages of maturation and are the 6th most common cancer-related cause of death. While the role of several oncogenes and tumor suppressors in the pathogenesis of B-cell neoplasms was established, recent research indicated the involvement of non-coding, regulatory sequences. Enhancers are DNA elements controlling gene expression in a cell type- and developmental stage-specific manner. They ensure proper differentiation and maturation of B cells, resulting in production of high affinity antibodies. However, the activity of enhancers can be redirected, setting B cells on the path towards cancer. In this review we discuss different mechanisms through which enhancers are exploited in malignant B cells, from the well-studied translocations juxtaposing oncogenes to immunoglobulin loci, through enhancer dysregulation by sequence variants and mutations, to enhancer hijacking by viruses. We also highlight the potential of therapeutic targeting of enhancers as a direction for future investigation.

V(DD)J recombination is an important and evolutionarily conserved mechanism for generating antibodies with unusually long CDR3s

The V(DD)J recombination is currently viewed as an aberrant and inconsequential variant of the canonical V(D)J recombination. Moreover, since the classical 12/23 rule for the V(D)J recombination fails to explain the V(DD)J recombination, the molecular mechanism of tandem D-D fusions has remained unknown since they were discovered three decades ago. Revealing this mechanism is a biomedically important goal since tandem fusions contribute to broadly neutralizing antibodies with ultralong CDR3s. We reveal previously overlooked cryptic nonamers in the recombination signal sequences of human IGHD genes and demonstrate that these nonamers explain the vast majority of tandem fusions in human repertoires. We further reveal large clonal lineages formed by tandem fusions in antigen-stimulated immunosequencing data sets, suggesting that such data sets contain many more tandem fusions than previously thought and that about a quarter of large clonal lineages with unusually long CDR3s are generated through tandem fusions. Finally, we developed the SEARCH-D algorithm for identifying D genes in mammalian genomes and applied it to the recently completed Vertebrate Genomes Project assemblies, nearly doubling the number of mammalian species with known D genes. Our analysis revealed cryptic nonamers in RSSs of many mammalian genomes, thus demonstrating that the V(DD)J recombination is not a "bug" but an important feature preserved throughout mammalian evolution.

Intra-Vκ Cluster Recombination Shapes the Ig Kappa Locus Repertoire

During V(D)J recombination, RAG proteins introduce DNA double-strand breaks (DSBs) at recombination signal sequences (RSSs) that contain either 12- or 23-nt spacer regions. Coordinated 12/23 cleavage predicts that DSBs at variable (V) gene segments should equal the level of breakage at joining (J) segments. Contrary to this, here we report abundant RAG-dependent DSBs at multiple Vκ gene segments independent of V-J rearrangement. We find that a large fraction of Vκ gene segments are flanked not only by a bone-fide 12 spacer but also an overlapping, 23-spacer flipped RSS. These compatible pairs of RSSs mediate recombination and deletion inside the Vκ cluster even in the complete absence of Jκ gene segments and support a V(D)J recombination center (RC) independent of the conventional Jκ-centered RC. We propose an improved model of Vκ-Jκ repertoire formation by incorporating these surprisingly frequent, evolutionarily conserved intra-Vκ cluster recombination events.

Brain cell somatic gene recombination and its phylogenetic foundations

A new form of somatic gene recombination (SGR) has been identified in the human brain that affects the Alzheimer's disease gene, amyloid precursor protein (APP). SGR occurs when a gene sequence is cut and recombined within a single cell's genomic DNA, generally independent of DNA replication and the cell cycle. The newly identified brain SGR produces genomic complementary DNAs (gencDNAs) lacking introns, which integrate into locations distinct from germline loci. This brief review will present an overview of likely related recombination mechanisms and genomic cDNA-like sequences that implicate evolutionary origins for brain SGR. Similarities and differences exist between brain SGR and VDJ recombination in the immune system, the first identified SGR form that now has a well-defined enzymatic machinery. Both require gene transcription, but brain SGR uses an RNA intermediate and reverse transcriptase (RT) activity, which are characteristics shared with endogenous retrotransposons. The identified gencDNAs have similarities to other cDNA-like sequences existing throughout phylogeny, including intron-less genes and inactive germline processed pseudogenes, with likely overlapping biosynthetic processes. gencDNAs arise somatically in an individual to produce multiple copies; can be functional; appear most frequently within postmitotic cells; have diverse sequences; change with age; and can change with disease state. Normally occurring brain SGR may represent a mechanism for gene optimization and long-term cellular memory, whereas its dysregulation could underlie multiple brain disorders and, potentially, other diseases like cancer. The involvement of RT activity implicates already Food and Drug Administration-approved RT inhibitors as possible near-term interventions for managing SGR-associated diseases and suggest next-generation therapeutics targeting SGR elements.

Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

RAG-2 deficiency results in fewer phosphorylated histone H2AX foci, but increased retinal ganglion cell death and altered axonal growth

DNA double-strand breaks (DSBs), selectively visualized as γ-H2AX⁺ foci, occur during the development of the central nervous system, including the retina, although their origin and biological significance are poorly understood. Mutant mice with DSB repair mechanism defects exhibit increased numbers of γ-H2AX⁺ foci, increased cell death during neural development, and alterations in axonogenesis in the embryonic retina. The aim of this study was to identify putative sources of DSBs. One of the identified DSBs sources is LINE-1 retrotransposition. While we did not detect changes in LINE-1 DNA content during the early period of cell death associated with retinal neurogenesis, retinal development was altered in mice lacking RAG-2, a component of the RAG-1,2-complex responsible for initiating somatic recombination in lymphocytes. Although γ-H2AX⁺ foci were less abundant in the rag2^−/− mouse retina, retinal ganglion cell death was increased and axonal growth and navigation were impaired in the RAG-2 deficient mice, a phenotype shared with mutant mice with defective DNA repair mechanisms. These findings demonstrate that RAG-2 is necessary for proper retinal development, and suggest that both DSB generation and repair are genuine processes intrinsic to neural development.

Germline DNA Retention in Murine and Human Rearranged T Cell Receptor Gene Coding Joints: Alternative Recombination Signal Sequences and V(D)J Recombinase Errors

The genes coding for the antigenic T cell receptor (TR) subunits are assembled in thymocytes from discrete V, D, and J genes by a site-specific recombination process. A tight control of this activity is required to prevent potentially detrimental recombination events. V, D, and J genes are flanked by semi-conserved nucleotide motives called recombination signal sequences (RSSs). V(D)J recombination is initiated by the precise introduction of a DNA double-strand break exactly at the border of the genes and their RSSs by the RAG recombinase. RSSs are therefore physically separated from the coding region of the genes before assembly of a rearranged TR gene. During a high throughput profiling of TRB genes in mice, we identified rearranged TRB genes in which part or all of a flanking RSS was retained in V-D or D-J coding joints. In some instances, this retention of germline DNA resulted from the use of an upstream alternative RSS. However, we also identified TRB sequences where retention of germline DNA occurred in the absence of alternative RSS, suggesting that RAG activity was mis-targeted during recombination. Similar events were also identified in human rearranged TRB and TRG genes. The use of alternative RSSs during V(D)J recombination illustrates the complexity of RAG-RSSs interactions during V(D)J recombination. While the frequency of errors resulting from mis-targeted RAG activity is very low, we believe that these RAG errors may be at the origin of oncogenic translocations and are a threat for genetic stability in developing lymphocytes.

Cross-Reactivity to Kynureninase Tolerizes B Cells That Express the HIV-1 Broadly Neutralizing Antibody 2F5

2F5 is an HIV-1 broadly neutralizing Ab that also binds the autoantigens kynureninase (KYNU) and anionic lipids. Generation of 2F5-like Abs is proscribed by immune tolerance, but it is unclear which autospecificity is responsible. We sampled the BCR repertoire of 2F5 knock-in mice before and after the first and second tolerance checkpoints. Nearly all small pre-B (precheckpoint) and 35-70% of anergic peripheral B cells (postcheckpoint) expressed the 2F5 BCR and maintained KYNU, lipid, and HIV-1 gp41 reactivity. In contrast, all postcheckpoint mature follicular (MF) B cells had undergone L chain editing that purged KYNU and gp41 binding but left lipid reactivity largely intact. We conclude that specificity for KYNU is the primary driver of tolerization of 2F5-expressing B cells. The MF and anergic B cell populations favored distinct collections of editor L chains; surprisingly, however, MF and anergic B cells also frequently expressed identical BCRs. These results imply that BCR autoreactivity is the primary determinant of whether a developing B cell enters the MF or anergic compartments, with a secondary role for stochastic factors that slightly mix the two pools. Our study provides mechanistic insights into how immunological tolerance impairs humoral responses to HIV-1 and supports activation of anergic B cells as a potential method for HIV-1 vaccination.

The fundamental role of chromatin loop extrusion in physiological V(D)J recombination

RAG endonuclease initiates IgH locus (Igh) V(D)J assembly in progenitor (pro)-B cells by joining Ds to J_Hs, before joining upstream V_Hs to DJ_H intermediates¹. In mouse pro-B cells, the CTCF-binding element (CBE)-anchored chromatin loop domain² at the 3’end of Igh contains an internal sub-domain spanning the 5’CBE anchor (IGCR1)³, the D_Hs, and a RAG-bound recombination center (RC)⁴. The RC comprises J_H-proximal D (DQ52), 4 J_Hs, and the intronic enhancer (“iEμ”)⁵. Robust RAG cleavage is restricted to paired V(D)J segments flanked by complementary recombination signal sequences (12RSSs and 23RSSs)⁶. Ds are flanked downstream and upstream by 12RSSs that, respectively, mediate deletional joining with convergently-oriented J_H-23RSSs and V_H-23RSSs⁶. Despite 12/23 compatibility, inversional D to J_H joining via upstream D-12RSSs is rare^7,8. Plasmid-based assays attributed lack of inversional D to J_H joining to sequence-based preference for downstream D-12RSSs⁹, as opposed to putative linear scanning mechanisms^10,11. Given recent findings that RAG linearly scans convergent CBE-anchored chromatin loops^4,12-14, potentially formed by cohesin-mediated loop extrusion^15-18, we revisited a scanning role. Here, we report that J_H-23RSS chromosomal orientation programs RC-bound RAG to linearly scan upstream chromatin in the 3’Igh sub-domain for convergently-oriented D-12RSSs and, thereby, to mediate deletional joining of all Ds, except RC-based DQ52 that joins by a diffusion-related mechanism. In a DQ52-based RC, formed in the absence of J_Hs, RAG bound by the downstream DQ52-RSS scans the downstream constant region exon-containing 3’Igh sub-domain in which scanning can be impeded by targeted nuclease-dead Cas9 (dCas9) binding, by transcription through repetitive Igh switch sequences, and by the 3’Igh CBE-based loop anchor. Notably, each scanning impediment focally increases RAG activity on potential substrate sequences within the impeded region. High resolution mapping of RC chromatin interactions reveals that such focal RAG targeting is associated with corresponding impediments to the loop extrusion process that drives chromatin past RC-bound RAG.

Detection of a novel CBFB-MYH11 fusion transcript in acute myeloid leukemia M1 with inv(16)(p13q22)

Acute myeloid leukemia (AML) with an inv(16)(p13q22) or t(16;16)(p13;q22) chromosomal abnormality represents one of the most common subtypes of de novo cases. These chromosomal rearrangements result in multiple CBFB-MYH11 fusion transcripts, with type-A being the most frequent. We here describe a unique case of de novo AML-M1, with inv(16)(p13q22), leading to an unusual CBFB-MYH11 fusion transcript, and der(7)t(7;11)(q31;q21). The fusion transcript involves a CBFB exon 5 with a breakpoint at nucleotide 754, an insertion of a 13-bp sequence of CBFB intron 5 at the fusion point, and the MYH11 exon 27 with a breakpoint at nucleotide 3464. To our knowledge, this CBFB-MYH11 fusion transcript has never been reported previously. The clinical characteristics of the present case are in line with previous reports suggesting that rare CBFB-MYH11 fusion transcripts lead to aberrant characteristics such as an atypical cytomorphology and additional cytogenetic abnormalities.

Epigenetic Enhancer Marks and Transcription Factor Binding Influence Vκ Gene Rearrangement in Pre-B Cells and Pro-B Cells

To date there has not been a study directly comparing relative Igκ rearrangement frequencies obtained from genomic DNA (gDNA) and cDNA and since each approach has potential biases, this is an important issue to clarify. Here we used deep sequencing to compare the unbiased gDNA and RNA Igκ repertoire from the same pre-B cell pool. We find that ~20% of Vκ genes have rearrangement frequencies ≥2-fold up or down in RNA vs. DNA libraries, including many members of the Vκ3, Vκ4, and Vκ6 families. Regression analysis indicates Ikaros and E2A binding are associated with strong promoters. Within the pre-B cell repertoire, we observed that individual Vκ genes rearranged at very different frequencies, and also displayed very different Jκ usage. Regression analysis revealed that the greatly unequal Vκ gene rearrangement frequencies are best predicted by epigenetic marks of enhancers. In particular, the levels of newly arising H3K4me1 peaks associated with many Vκ genes in pre-B cells are most predictive of rearrangement levels. Since H3K4me1 is associated with long range chromatin interactions which are created during locus contraction, our data provides mechanistic insight into unequal rearrangement levels. Comparison of Igκ rearrangements occurring in pro-B cells and pre-B cells from the same mice reveal a pro-B cell bias toward usage of Jκ-distal Vκ genes, particularly Vκ10-96 and Vκ1-135. Regression analysis indicates that PU.1 binding is the highest predictor of Vκ gene rearrangement frequency in pro-B cells. Lastly, the repertoires of iEκ^−/− pre-B cells reveal that iEκ actively influences Vκ gene usage, particularly Vκ3 family genes, overlapping with a zone of iEκ-regulated germline transcription. These represent new roles for iEκ in addition to its critical function in promoting overall Igκ rearrangement. Together, this study provides insight into many aspects of Igκ repertoire formation.

Loss-of-function but not dominant-negative intragenic IKZF1 deletions are associated with an adverse prognosis in adult BCR-ABL-negative acute lymphoblastic leukemia

Genetic alterations of the transcription factor IKZF1 ("IKAROS") are detected in around 15-30% of cases of BCR-ABL-negative B-cell precursor acute lymphoblastic leukemia. Different types of intragenic deletions have been observed, resulting in a functionally inactivated allele ("loss-of-function") or in "dominant-negative" isoforms. The prognostic impact of these alterations especially in adult acute lymphoblastic leukemia is not well defined. We analyzed 482 well-characterized cases of adult BCR-ABL-negative B-precursor acute lymphoblastic leukemia uniformly treated in the framework of the GMALL studies and detected IKZF1 alterations in 128 cases (27%). In 20%, the IKZF1 alteration was present in a large fraction of leukemic cells ("high deletion load") while in 7% it was detected only in small subclones ("low deletion load"). Some patients showed more than one IKZF1 alteration (8%). Patients exhibiting a loss-of-function isoform with high deletion load had a shorter overall survival (OS at 5 years 28% vs. 59%; P<0.0001), also significant in a subgroup analysis of standard risk patients according to GMALL classification (OS at 5 years 37% vs. 68%; P=0.0002). Low deletion load or dominant-negative IKZF1 alterations had no prognostic impact. The results thus suggest that there is a clear distinction between loss-of-function and dominant-negative IKZF1 deletions. Affected patients should thus be monitored for minimal residual disease carefully to detect incipient relapses at an early stage and they are potential candidates for alternative or intensified treatment regimes. (clinicaltrials.gov identifiers: 00199056 and 00198991).

HIV Tat induces a prolonged MYC relocalization next to IGH in circulating B-cells

With combined antiretroviral therapy (cART), the risk for HIV-infected individuals to develop a non-Hodgkin lymphoma is diminished. However, the incidence of Burkitt lymphoma (BL) remains strikingly elevated. Most BL present a t(8;14) chromosomal translocation which must take place at a time of spatial proximity between the translocation partners. The two partner genes, MYC and IGH, were found colocalized only very rarely in the nuclei of normal peripheral blood B-cells examined using 3D-FISH while circulating B-cells from HIV-infected individuals whose exhibited consistently elevated levels of MYC-IGH colocalization. In vitro, incubating normal B-cells from healthy donors with a transcriptionally active form of the HIV-encoded Tat protein rapidly activated transcription of the nuclease-encoding RAG1 gene. This created DNA damage, including in the MYC gene locus which then moved towards the center of the nucleus where it sustainably colocalized with IGH up to 10-fold more frequently than in controls. In vivo, this could be sufficient to account for the elevated risk of BL-specific chromosomal translocations which would occur following DNA double strand breaks triggered by AID in secondary lymph nodes at the final stage of immunoglobulin gene maturation. New therapeutic attitudes can be envisioned to prevent BL in this high risk group.

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.

Significant Role of Segmental Duplications and SIDD Sites in Chromosomal Translocations of Hematological Malignancies: A Multi-parametric Bioinformatic Analysis

Recurrent non-random chromosomal translocations are hallmark characteristics of leukemogenesis, and however, molecular mechanisms underlying these rearrangements are less explored. The fundamental question is, why and how chromosomes break and reunite so precisely in the genome. Meticulous understanding of mechanism leading to chromosomal rearrangement can be achieved by characterizing breakpoints. To address this hypothesis, a novel multi-parametric computational approach for characterization of major leukemic translocations within and around breakpoint region was performed. To best of our knowledge, this bioinformatic analysis is unique in finding the presence of segmental duplications (SDs) flanking breakpoints of all major leukemic translocation. Breakpoint islands (BpIs) were analyzed for stress-induced duplex destabilization (SIDD) sites along with other complex genomic architecture and physicochemical properties. Our study distinctly emphasizes on the probable correlative role of SDs, SIDD sites and various genomic features in the occurrence of breakpoints. Further, it also highlights potential features which may be playing a crucial role in causing double-strand breaks, leading to translocation.

The Genome Conformation As an Integrator of Multi-Omic Data: The Example of Damage Spreading in Cancer

Publicly available multi-omic databases, in particular if associated with medical annotations, are rich resources with the potential to lead a rapid transition from high-throughput molecular biology experiments to better clinical outcomes for patients. In this work, we propose a model for multi-omic data integration (i.e., genetic variations, gene expression, genome conformation, and epigenetic patterns), which exploits a multi-layer network approach to analyse, visualize, and obtain insights from such biological information, in order to use achieved results at a macroscopic level. Using this representation, we can describe how driver and passenger mutations accumulate during the development of diseases providing, for example, a tool able to characterize the evolution of cancer. Indeed, our test case concerns the MCF-7 breast cancer cell line, before and after the stimulation with estrogen, since many datasets are available for this case study. In particular, the integration of data about cancer mutations, gene functional annotations, genome conformation, epigenetic patterns, gene expression, and metabolic pathways in our multi-layer representation will allow a better interpretation of the mechanisms behind a complex disease such as cancer. Thanks to this multi-layer approach, we focus on the interplay of chromatin conformation and cancer mutations in different pathways, such as metabolic processes, that are very important for tumor development. Working on this model, a variance analysis can be implemented to identify normal variations within each omics and to characterize, by contrast, variations that can be accounted to pathological samples compared to normal ones. This integrative model can be used to identify novel biomarkers and to provide innovative omic-based guidelines for treating many diseases, improving the efficacy of decision trees currently used in clinic.

COBL is a novel hotspot for IKZF1 deletions in childhood acute lymphoblastic leukemia

IKZF1 deletion (ΔIKZF1) is an important predictor of relapse in childhood B-cell precursor acute lymphoblastic leukemia. Because of its clinical importance, we previously mapped breakpoints of intragenic deletions and developed a multiplex PCR assay to detect recurrent intragenic ΔIKZF1. Since the multiplex PCR was not able to detect complete deletions (IKZF1 Δ1-8), which account for ~30% of all ΔIKZF1, we aimed at investigating the genomic scenery of IKZF1 Δ1-8. Six samples of cases with IKZF1 Δ1-8 were analyzed by microarray assay, which identified monosomy 7, isochromosome 7q, and large interstitial deletions presenting breakpoints within COBL gene. Then, we established a multiplex ligation-probe amplification (MLPA) assay and screened copy number alterations within chromosome 7 in 43 diagnostic samples with IKZF1 Δ1-8. Our results revealed that monosomy and large interstitial deletions within chromosome 7 are the main causes of IKZF1 Δ1-8. Detailed analysis using long distance inverse PCR showed that six patients (16%) had large interstitial deletions starting within intronic regions of COBL at diagnosis, which is ~611 Kb downstream of IKZF1, suggesting that COBL is a hotspot for ΔIKZF1. We also investigated a series of 25 intragenic deletions (Δ2-8, Δ3-8 or Δ4-8) and 24 relapsed samples, and found one IKZF1-COBL tail-to-tail fusion, thus supporting that COBL is a novel hotspot for ΔIKZF1. Finally, using RIC score methodology, we show that breakpoint sequences of IKZF1 Δ1-8 are not analog to RAG-recognition sites, suggesting a different mechanism of error promotion than that suggested for intragenic ΔIKZF1.

Two Mutually Exclusive Local Chromatin States Drive Efficient V(D)J Recombination

Variable (V), diversity (D), and joining (J) (V(D)J) recombination is the first determinant of antigen receptor diversity. Understanding how recombination is regulated requires a comprehensive, unbiased readout of V gene usage. We have developed VDJ sequencing (VDJ-seq), a DNA-based next-generation-sequencing technique that quantitatively profiles recombination products. We reveal a 200-fold range of recombination efficiency among recombining V genes in the primary mouse Igh repertoire. We used machine learning to integrate these data with local chromatin profiles to identify combinatorial patterns of epigenetic features that associate with active V_H gene recombination. These features localize downstream of V_H genes and are excised by recombination, revealing a class of cis-regulatory element that governs recombination, distinct from expression. We detect two mutually exclusive chromatin signatures at these elements, characterized by CTCF/RAD21 and PAX5/IRF4, which segregate with the evolutionary history of associated V_H genes. Thus, local chromatin signatures downstream of V_H genes provide an essential layer of regulation that determines recombination efficiency.

•

VDJ-seq enables precise quantification of antibody V(D)J recombination products

•

Two distinct cis-regulatory designs characterize actively recombining V genes

•

Putative recombination regulatory elements map downstream of mouse Igh V genes

•

Recombination regulatory architecture reflects the V genes’ evolutionary history

A LAIR1 insertion generates broadly reactive antibodies against malaria variant antigens

Plasmodium falciparum antigens expressed on the surface of infected erythrocytes are important targets of naturally acquired immunity against malaria, but their high number and variability provide the pathogen with a powerful means of escape from host antibodies^1–4. Although broadly reactive antibodies against these antigens could be useful as therapeutics and in vaccine design, their identification has proven elusive. Here, we report the isolation of human monoclonal antibodies that recognize erythrocytes infected by different P. falciparum isolates and opsonize these cells by binding to members of the RIFIN family. These antibodies acquired broad reactivity through a novel mechanism of insertion of a large DNA fragment between the V and DJ segments. The insert, which is both necessary and sufficient for binding to RIFINs, encodes the entire 100 amino acid collagen-binding domain of LAIR-1, an Ig superfamily inhibitory receptor encoded on chromosome 19. In each of the two donors studied, the antibodies are produced by a single expanded B cell clone and carry distinct somatic mutations in the LAIR-1 domain that abolish binding to collagen and increase binding to infected erythrocytes. These findings illustrate, with a biologically relevant example, a novel mechanism of antibody diversification by interchromosomal DNA transposition and demonstrate the existence of conserved epitopes that may be suitable candidates for the development of a malaria vaccine.

Chromosomal Loop Domains Direct the Recombination of Antigen Receptor Genes

RAG initiates antibody V(D)J recombination in developing lymphocytes by generating "on-target" DNA breaks at matched pairs of bona fide recombination signal sequences (RSSs). We employ bait RAG-generated breaks in endogenous or ectopically inserted RSS pairs to identify huge numbers of RAG "off-target" breaks. Such breaks occur at the simple CAC motif that defines the RSS cleavage site and are largely confined within convergent CTCF-binding element (CBE)-flanked loop domains containing bait RSS pairs. Marked orientation dependence of RAG off-target activity within loops spanning up to 2 megabases implies involvement of linear tracking. In this regard, major RAG off-targets in chromosomal translocations occur as convergent RSS pairs at enhancers within a loop. Finally, deletion of a CBE-based IgH locus element disrupts V(D)J recombination domains and, correspondingly, alters RAG on- and off-target distributions within IgH. Our findings reveal how RAG activity is developmentally focused and implicate mechanisms by which chromatin domains harness biological processes within them.

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.

Genomic landscape of cutaneous T cell lymphoma

Cutaneous T cell lymphoma (CTCL) is a non-Hodgkin lymphoma of skin-homing T lymphocytes. We performed exome and whole genome DNA sequence and RNA sequencing on purified CTCL and matched normal cells. The results implicate mutations in 17 genes in CTCL pathogenesis, including genes involved in T cell activation and apoptosis, NFκB signaling, chromatin remodeling, and DNA damage response. CTCL is distinctive in that somatic copy number variants (SCNVs) comprise 92% of all driver mutations (mean of 11.8 pathogenic SCNVs vs. 1.0 somatic single nucleotide variants per CTCL). These findings have implications for novel therapeutics.

Characterization of chromosomal translocation breakpoint sequences in solid tumours: "an in silico analysis"

Chromosomal translocations that results in formation and activation of fusion oncogenes are observed in numerous solid malignancies since years back. Expression of fusion kinases in these cancers drives the initiation & progression that ultimately leads to tumour development and thus comes out to be clinically imperative in terms of diagnosis and treatment of cancer. Nonetheless, molecular mechanisms beneath these translocations remained unexplored consequently limiting our knowledge of carcinogenesis and hence is the current field where further research is required. The issue of prime focus is the precision with which the chromosomes breaks and reunites within genome. Characterization of Genomic sequences located at Breakpoint region may direct us towards the thorough understanding of mechanism leading to chromosomal rearrangement. A unique computational multi-parametric analysis was performed for characterization of genomic sequence within and around breakpoint region. This study turns out to be novel as it reveals the occurrence of Segmental Duplications flanking the breakpoints of all translocation. Breakpoint Islands were also investigated for the presence of other intricate genomic architecture and various physico-chemical parameters. Our study particularly highlights the probable role of SDs and specific genomic features in precise chromosomal breakage. Additionally, it pinpoints the potential features that may be significant for double-strand breaks leading to chromosomal rearrangements.

The proximal J kappa germline-transcript promoter facilitates receptor editing through control of ordered recombination

V(D)J recombination creates antibody light chain diversity by joining a Vκ gene segment with one of four Jκ segments. Two Jκ germline-transcript (GT) promoters control Vκ-Jκ joining, but the mechanisms that govern Jκ choice are unclear. Here, we show in gene-targeted mice that the proximal GT promoter helps targeting rearrangements to Jκ1 by preventing premature DNA breaks at Jκ2. Consequently, cells lacking the proximal GT promoter show a biased utilization of downstream Jκ segments, resulting in a diminished potential for receptor editing. Surprisingly, the proximal—in contrast to the distal—GT promoter is transcriptionally inactive prior to Igκ recombination, indicating that its role in Jκ choice is independent of classical promoter function. Removal of the proximal GT promoter increases H3K4me3 levels at Jκ segments, suggesting that this promoter could act as a suppressor of recombination by limiting chromatin accessibility to RAG. Our findings identify the first cis-element critical for Jκ choice and demonstrate that ordered Igκ recombination facilitates receptor editing.

Mechanisms that can promote peripheral B-cell lymphoma in ATM-deficient mice

The Ataxia Telangiectasia-mutated (ATM) kinase senses DNA double-strand breaks (DSB) and facilitates their repair. In humans, ATM deficiency predisposes to B- and T-cell lymphomas, but in mice it leads only to thymic lymphomas. We tested the hypothesis that increased DSB frequency at a cellular oncogene could promote B-cell lymphoma by generating ATM-deficient mice with a V(D)J recombination target (DJβ cassette) within c-myc intron 1 ("DA" mice). We also generated ATM-deficient mice carrying an Eμ-Bcl-2 transgene (AB mice) to test whether enhanced cellular survival could promote B-cell lymphomas. About 30% of DA or AB mice and nearly 100% of mice harboring the combined genotypes (DAB mice) developed mature B-cell lymphomas. In all genotypes, B-cell tumors harbored oncogenic c-myc amplification generated by breakage-fusion-bridge (BFB) from dicentric chromosomes formed through fusion of IgH V(D)J recombination-associated DSBs on chromosome 12 to sequences downstream of c-myc on chromosome 15. AB tumors demonstrate that B lineage cells harboring spontaneous DSBs leading to IgH/c-myc dicentrics are blocked from progressing to B-cell lymphomas by cellular apoptotic responses. DA and DAB tumor translocations were strictly linked to the cassette, but occurred downstream, frequently in a 6-kb region adjacent to c-myc that harbors multiple cryptic V(D)J recombination targets, suggesting that bona fide V(D)J target sequences may activate linked cryptic targets. Our findings indicate that ATM deficiency allows IgH V(D)J recombination DSBs in developing B cells to generate dicentric translocations that, via BFB cycles, lead to c-myc-activating oncogenic translocations and amplifications in mature B cells.

Breakpoint sites disclose the role of the V(D)J recombination machinery in the formation of T-cell receptor (TCR) and non-TCR associated aberrations in T-cell acute lymphoblastic leukemia

Aberrant recombination between T-cell receptor genes and oncogenes gives rise to chromosomal translocations that are genetic hallmarks in several subsets of human T-cell acute lymphoblastic leukemias. The V(D)J recombination machinery has been shown to play a role in the formation of these T-cell receptor translocations. Other, non-T-cell receptor chromosomal aberrations, such as SIL-TAL1 deletions, have likewise been recognized as V(D)J recombination associated aberrations. Despite the postulated role of V(D)J recombination, the extent of the V(D)J recombination machinery involvement in the formation of T-cell receptor and non-T-cell receptor aberrations in T-cell acute lymphoblastic leukemia is still poorly understood. We performed a comprehensive in silico and ex vivo evaluation of 117 breakpoint sites from 22 different T-cell receptor translocation partners as well as 118 breakpoint sites from non-T-cell receptor chromosomal aberrations. Based on this extensive set of breakpoint data, we provide a comprehensive overview of T-cell receptor and oncogene involvement in T-ALL. Moreover, we assessed the role of the V(D)J recombination machinery in the formation of chromosomal aberrations, and propose an up-dated mechanistic classification on how the V(D)J recombination machinery contributes to the formation of T-cell receptor and non-T-cell receptor aberrations in human T-cell acute lymphoblastic leukemia.

Illegitimate V(D)J recombination-mediated deletions in Notch1 and Bcl11b are not sufficient for extensive clonal expansion and show minimal age or sex bias in frequency or junctional processing

Illegitimate V(D)J recombination at oncogenes and tumor suppressor genes is implicated in formation of several T cell malignancies. Notch1 and Bcl11b, genes involved in developing T cell specification, selection, proliferation, and survival, were previously shown to contain hotspots for deletional illegitimate V(D)J recombination associated with radiation-induced thymic lymphoma. Interestingly, these deletions were also observed in wild-type animals. In this study, we conducted frequency, clonality, and junctional processing analyses of Notch1 and Bcl11b deletions during mouse development and compared results to published analyses of authentic V(D)J rearrangements at the T cell receptor beta (TCRβ) locus and illegitimate V(D)J deletions observed at the human, nonimmune HPRT1 locus not involved in T cell malignancies. We detect deletions in Notch1 and Bcl11b in thymic and splenic T cell populations, consistent with cells bearing deletions in the circulating lymphocyte pool. Deletions in thymus can occur in utero, increase in frequency between fetal and postnatal stages, are detected at all ages examined between fetal and 7 months, exhibit only limited clonality (contrasting with previous results in radiation-sensitive mouse strains), and consistent with previous reports are more frequent in Bcl11b, partially explained by relatively high Recombination Signal Information Content (RIC) scores. Deletion junctions in Bcl11b exhibit greater germline nucleotide loss, while in Notch1 palindromic (P) nucleotides are more abundant, although average P nucleotide length is similar for both genes and consistent with results at the TCRβ locus. Non-templated (N) nucleotide insertions appear to increase between fetal and postnatal stages for Notch1, consistent with normal terminal deoxynucleotidyl transferase (TdT) activity; however, neonatal Bcl11b junctions contain elevated levels of N insertions. Finally, contrasting with results at the HPRT1 locus, we find no obvious age or gender bias in junctional processing, and inverted repeats at recessed coding ends (Pr nucleotides) correspond mostly to single-base additions consistent with normal TdT activity.

NuChart: an R package to study gene spatial neighbourhoods with multi-omics annotations

Long-range chromosomal associations between genomic regions, and their repositioning in the 3D space of the nucleus, are now considered to be key contributors to the regulation of gene expression and important links have been highlighted with other genomic features involved in DNA rearrangements. Recent Chromosome Conformation Capture (3C) measurements performed with high throughput sequencing (Hi-C) and molecular dynamics studies show that there is a large correlation between colocalization and coregulation of genes, but these important researches are hampered by the lack of biologists-friendly analysis and visualisation software. Here, we describe NuChart, an R package that allows the user to annotate and statistically analyse a list of input genes with information relying on Hi-C data, integrating knowledge about genomic features that are involved in the chromosome spatial organization. NuChart works directly with sequenced reads to identify the related Hi-C fragments, with the aim of creating gene-centric neighbourhood graphs on which multi-omics features can be mapped. Predictions about CTCF binding sites, isochores and cryptic Recombination Signal Sequences are provided directly with the package for mapping, although other annotation data in bed format can be used (such as methylation profiles and histone patterns). Gene expression data can be automatically retrieved and processed from the Gene Expression Omnibus and ArrayExpress repositories to highlight the expression profile of genes in the identified neighbourhood. Moreover, statistical inferences about the graph structure and correlations between its topology and multi-omics features can be performed using Exponential-family Random Graph Models. The Hi-C fragment visualisation provided by NuChart allows the comparisons of cells in different conditions, thus providing the possibility of novel biomarkers identification. NuChart is compliant with the Bioconductor standard and it is freely available at ftp://fileserver.itb.cnr.it/nuchart.

Deep sequencing of the murine IgH repertoire reveals complex regulation of nonrandom V gene rearrangement frequencies

A diverse Ab repertoire is formed through the rearrangement of V, D, and J segments at the IgH (Igh) loci. The C57BL/6 murine Igh locus has >100 functional VH gene segments that can recombine to a rearranged DJH. Although the nonrandom usage of VH genes is well documented, it is not clear what elements determine recombination frequency. To answer this question, we conducted deep sequencing of 5'-RACE products of the Igh repertoire in pro-B cells, amplified in an unbiased manner. Chromatin immunoprecipitation-sequencing results for several histone modifications and RNA polymerase II binding, RNA-sequencing for sense and antisense noncoding germline transcripts, and proximity to CCCTC-binding factor (CTCF) and Rad21 sites were compared with the usage of individual V genes. Computational analyses assessed the relative importance of these various accessibility elements. These elements divide the Igh locus into four epigenetically and transcriptionally distinct domains, and our computational analyses reveal different regulatory mechanisms for each region. Proximal V genes are relatively devoid of active histone marks and noncoding RNA in general, but having a CTCF site near their recombination signal sequence is critical, suggesting that being positioned near the base of the chromatin loops is important for rearrangement. In contrast, distal V genes have higher levels of histone marks and noncoding RNA, which may compensate for their poorer recombination signal sequences and for being distant from CTCF sites. Thus, the Igh locus has evolved a complex system for the regulation of V(D)J rearrangement that is different for each of the four domains that comprise this locus.