V(D)J recombination: a functional definition of the joining signals

REV7: a small but mighty regulator of genome maintenance and cancer development

REV7, also known as MAD2B, MAD2L2, and FANCV, is a HORMA-domain family protein crucial to multiple genome stability pathways. REV7's canonical role is as a member of polymerase ζ, a specialized translesion synthesis polymerase essential for DNA damage tolerance. REV7 also ensures accurate cell cycle progression and prevents premature mitotic progression by sequestering an anaphase-promoting complex/cyclosome activator. Additionally, REV7 supports genome integrity by directing double-strand break repair pathway choice as part of the recently characterized mammalian shieldin complex. Given that genome instability is a hallmark of cancer, it is unsurprising that REV7, with its numerous genome maintenance roles, is implicated in multiple malignancies, including ovarian cancer, glioma, breast cancer, malignant melanoma, and small-cell lung cancer. Moreover, high REV7 expression is associated with poor prognoses and treatment resistance in these and other cancers. Promisingly, early studies indicate that REV7 suppression enhances sensitivity to chemotherapeutics, including cisplatin. This review aims to provide a comprehensive overview of REV7's myriad roles in genome maintenance and other functions as well as offer an updated summary of its connections to cancer and treatment resistance.

Early-life tobacco exposure is causally implicated in aberrant RAG-mediated recombination in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common cancer in children, yet few environmental risk factors have been identified. We previously found an association between early-life tobacco smoke exposure and frequency of somatic deletions of 8 leukemia driver genes among childhood ALL patients in the California Childhood Leukemia Study. To expand analysis genome-wide and examine potential mechanisms, we conducted tumor whole-genome sequencing in 35 ALL patients, including 18 with high prenatal tobacco exposure and 17 with low exposure as determined by established epigenetic biomarkers. High tobacco exposure patients had significantly more structural variants (P < .001) and deletions (P = .001) genome-wide than low exposure patients. Investigation of off-target RAG recombination revealed that 41% of deletions in the high tobacco exposure patients were putatively RAG-mediated (full RAG motif identified at one or both breakpoints) compared with only 21% in the low exposure group (P = .001). In a multilevel model, deletions in high tobacco exposure patients were 2.44-fold (95% CI:1.13-5.38) more likely to be putatively RAG-mediated than deletions in low exposure patients. No point mutational signatures were associated with prenatal tobacco exposure. Our findings suggest that early-life tobacco smoke exposure may promote leukemogenesis by driving development of somatic deletions in pre-leukemic lymphocytes via off-target RAG recombination.

Subgenomic T cell receptor alpha and delta (TRA/TRD) loci in common carp

In jawed vertebrates, the T cell receptor alpha (TRA) and delta (TRD) genes, which encode the TRα and TRδ chains, respectively, are located as a nested structure on a single chromosome. To date, no animal has been reported to harbor multiple TRA/TRD loci on different chromosomes. Therefore, herein, we describe the first full annotation of the TRA/TRD genomic regions of common carp, an allo-tetraploid fish species that experiences cyprinid-specific whole-genome duplication (WGD) in evolution. Fine genomic maps of TRA/TRD genomic regions 1 and 2, on LG30 and LG22, respectively, were constructed using the annotations of complete sets of TRA and TRD genes, including TRA/TRD variable (V), TRA junction (J), and constant (C), TRD diversity (D), and the J and C genes. The structure and synteny of the TRA/TRD genomic regions were highly conserved in zebrafish, indicating that these regions are on individual chromosomes. Furthermore, analysis of the variable regions of the TRA and TRD genes in a monoclonal T cell line revealed that both subgenomic regions 1 and 2 were indeed rearranged. Although carp TRAV and TRDV genes were phylogenetically divided into different lineages, they were mixed and organized into the TRA/TRD V gene clusters on the genome, similar to that in other vertebrates. Notably, 285 potential TRA/TRD V genes were detected in the TRA/TRD genomic regions, which is the most abundant number of genes in vertebrates and approximately two-fold that in zebrafish. The recombination signal sequences (RSSs) at the end of each V gene differed between TRAV and TRDV, suggesting that RSS variations might separate each V gene into a TRα or TRδ chain. This study is the first to describe subgenomic TRA/TRD loci in animals. Our findings provide fundamental insights to elucidate the impact of WGD on the evolution of immune repertoire.

Locus folding mechanisms determine modes of antigen receptor gene assembly

The dynamic folding of genomes regulates numerous biological processes, including antigen receptor (AgR) gene assembly. We show that, unlike other AgR loci, homotypic chromatin interactions and bidirectional chromosome looping both contribute to structuring Tcrb for efficient long-range V(D)J recombination. Inactivation of the CTCF binding element (CBE) or promoter at the most 5'Vβ segment (Trbv1) impaired loop extrusion originating locally and extending to DβJβ CBEs at the opposite end of Tcrb. Promoter or CBE mutation nearly eliminated Trbv1 contacts and decreased RAG endonuclease-mediated Trbv1 recombination. Importantly, Trbv1 rearrangement can proceed independent of substrate orientation, ruling out scanning by DβJβ-bound RAG as the sole mechanism of Vβ recombination, distinguishing it from Igh. Our data indicate that CBE-dependent generation of loops cooperates with promoter-mediated activation of chromatin to juxtapose Vβ and DβJβ segments for recombination through diffusion-based synapsis. Thus, the mechanisms that fold a genomic region can influence molecular processes occurring in that space, which may include recombination, repair, and transcriptional programming.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Proteogenomic and V(D)J Analysis of Human Decidual T Cells Highlights Unique Transcriptional Programming and Clonal Distribution

Immunological tolerance toward the semiallogeneic fetus is one of many maternal adaptations required for a successful pregnancy. T cells are major players of the adaptive immune system and balance tolerance and protection at the maternal-fetal interface; however, their repertoire and subset programming are still poorly understood. Using emerging single-cell RNA sequencing technologies, we simultaneously obtained transcript, limited protein, and receptor repertoire at the single-cell level, from decidual and matched maternal peripheral human T cells. The decidua maintains a tissue-specific distribution of T cell subsets compared with the periphery. We find that decidual T cells maintain a unique transcriptome programming, characterized by restraint of inflammatory pathways by overexpression of negative regulators (DUSP, TNFAIP3, ZFP36) and expression of PD-1, CTLA-4, TIGIT, and LAG3 in some CD8 clusters. Finally, analyzing TCR clonotypes demonstrated decreased diversity in specific decidual T cell populations. Overall, our data demonstrate the power of multiomics analysis in revealing regulation of fetal-maternal immune coexistence.

Origin and evolutionary malleability of T cell receptor α diversity

Lymphocytes of vertebrate adaptive immune systems acquired the capability to assemble, from split genes in the germline, billions of functional antigen receptors^1-3. These receptors show specificity; unlike the broadly tuned receptors of the innate system, antibodies (Ig) expressed by B cells, for instance, can accurately distinguish between the two enantiomers of organic acids⁴, whereas T cell receptors (TCRs) reliably recognize single amino acid replacements in their peptide antigens⁵. In developing lymphocytes, antigen receptor genes are assembled from a comparatively small set of germline-encoded genetic elements in a process referred to as V(D)J recombination^6,7. Potential self-reactivity of some antigen receptors arising from the quasi-random somatic diversification is suppressed by several robust control mechanisms^8-12. For decades, scientists have puzzled over the evolutionary origin of somatically diversifying antigen receptors^13-16. It has remained unclear how, at the inception of this mechanism, immunologically beneficial expanded receptor diversity was traded against the emerging risk of destructive self-recognition. Here we explore the hypothesis that in early vertebrates, sequence microhomologies marking the ends of recombining elements became the crucial targets of selection determining the outcome of non-homologous end joining-based repair of DNA double-strand breaks generated during RAG-mediated recombination. We find that, across the main clades of jawed vertebrates, TCRα repertoire diversity is best explained by species-specific extents of such sequence microhomologies. Thus, selection of germline sequence composition of rearranging elements emerges as a major factor determining the degree of diversity of somatically generated antigen receptors.

Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1⁺ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

Characterization of T-cell receptors and immunoglobulin heavy chains loci and identification of T/B cell clusters in teleost

Bacterial disease is one of the important factors leading to economic losses in the turbot (Scophthalmus maximus) cultivation industry. T lymphocytes are major components of cellular immunity, whereas B lymphocytes produce immunoglobulins (Ig) that are key elements of humoral immune responses against infection. However, the genomic organization of genes encoding T-cell receptors (TCR) and immunoglobulin heavy chains (IgHs) in turbot remains largely unknown. In this study, abundant full-length transcripts of TCRs and IgHs were sequenced by Isoform-sequencing (Iso-seq), and we investigated and annotated the V, D, J and C gene loci of TCRα, TCRβ, IgT, IgM and IgD in turbot. Furthermore, through single-cell RNA sequencing (scRNA-seq) of blood leukocytes, we confirmed that these identified TCRs and IgHs were highly expressed in T/B cell clusters, respectively. Meanwhile, we also identified the IgM⁺IgD⁺ B and IgT⁺ B cells with differential gene expression profiles and potential functions. Taken together, our results provide a comprehensive understanding of TCRs and IgHs loci in turbot, which will contribute to evolutionary and functional characterization of T and B lymphocytes in teleost.

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.

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.

pmTR database: population matched (pm) germline allelic variants of T-cell receptor (TR) loci

The IMGT database profiles the TR germline alleles for all four TR loci (TRA, TRB, TRG and TRD), however, it does not comprise of the information regarding population specificity and allelic frequencies of these germline alleles. The specificity of allelic variants to different human populations can, however, be a rich source of information when studying the genetic basis of population-specific immune responses in disease and in vaccination. Therefore, we meticulously identified true germline alleles enriched with complete TR allele sequences and their frequencies across 26 different human populations, profiled by "1000 Genomes data". We identified 205 TRAV, 249 TRBV, 16 TRGV and 5 TRDV germline alleles supported by at least four haplotypes. The diversity of germline allelic variants in the TR loci is the highest in Africans, while the majority of the Non-African alleles are specific to the Asian populations, suggesting a diverse profile of TR germline alleles in different human populations. Interestingly, the alleles in the IMGT database are frequent and common across all five super-populations. We believe that this new set of germline TR sequences represents a valuable new resource which we have made available through the new population-matched TR (pmTR) database, accessible via https://pmtrig.lumc.nl/ .

Diversity upon diversity: linking DNA double-strand break repair to blood cancer health disparities

Chromosomal translocations arising from aberrant repair of multiple DNA double-strand breaks (DSBs) are a defining characteristic of many cancers. DSBs are an essential part of physiological processes in antibody-producing B cells. The B cell environment is poised to generate genome instability leading to translocations relevant to the pathology of blood cancers. These are a diverse set of cancers, but limited data from under-represented groups have pointed to health disparities associated with each. We focus on the DSBs that occur in developing B cells and propose the most likely mechanism behind the formation of translocations. We also highlight specific cancers in which these rearrangements occur and address the growing concern of health disparities associated with them.

Longitudinal Dynamics of Human B-Cell Response at the Single-Cell Level in Response to Tdap Vaccination

To mount an adequate immune response against pathogens, stepwise mutation and selection processes are crucial functions of the adaptive immune system. To better characterize a successful vaccination response, we performed longitudinal (days 0, 5, 7, 10, and 14 after Boostrix vaccination) analysis of the single-cell transcriptome as well as the B-cell receptor (BCR) repertoire (scBCR-rep) in plasma cells of an immunized donor and compared it with baseline B-cell characteristics as well as flow cytometry findings. Based on the flow cytometry knowledge and literature findings, we discriminated individual B-cell subsets in the transcriptomics data and traced over-time maturation of plasmablasts/plasma cells (PB/PCs) and identified the pathways associated with the plasma cell maturation. We observed that the repertoire in PB/PCs differed from the baseline B-cell repertoire e.g., regarding expansion of unique clones in post-vaccination visits, high usage of IGHG1 in expanded clones, increased class-switching events post-vaccination represented by clonotypes spanning multiple IGHC classes and positive selection of CDR3 sequences over time. Importantly, the Variable gene family-based clustering of BCRs represented a similar measure as the gene-based clustering, but certainly improved the clustering of BCRs, as BCRs from duplicated Variable gene families could be clustered together. Finally, we developed a query tool to dissect the immune response to the components of the Boostrix vaccine. Using this tool, we could identify the BCRs related to anti-tetanus and anti-pertussis toxoid BCRs. Collectively, we developed a bioinformatic workflow which allows description of the key features of an ongoing (longitudinal) immune response, such as activation of PB/PCs, Ig class switching, somatic hypermutation, and clonal expansion, all of which are hallmarks of antigen exposure, followed by mutation & selection processes.

Population matched (pm) germline allelic variants of immunoglobulin (IG) loci: Relevance in infectious diseases and vaccination studies in human populations

Immunoglobulin (IG) loci harbor inter-individual allelic variants in many different germline IG variable, diversity and joining genes of the IG heavy (IGH), kappa (IGK) and lambda (IGL) loci, which together form the genetic basis of the highly diverse antigen-specific B-cell receptors. These allelic variants can be shared between or be specific to human populations. The current immunogenetics resources gather the germline alleles, however, lack the population specificity of the alleles which poses limitations for disease-association studies related to immune responses in different human populations. Therefore, we systematically identified germline alleles from 26 different human populations around the world, profiled by "1000 Genomes" data. We identified 409 IGHV, 179 IGKV, and 199 IGLV germline alleles supported by at least seven haplotypes. The diversity of germline alleles is the highest in Africans. Remarkably, the variants in the identified novel alleles show strikingly conserved patterns, the same as found in other IG databases, suggesting over-time evolutionary selection processes. We could relate the genetic variants to population-specific immune responses, e.g. IGHV1-69 for flu in Africans. The population matched IG (pmIG) resource will enhance our understanding of the SHM-related B-cell receptor selection processes in (infectious) diseases and vaccination within and between different human populations.

Poor quality Vβ recombination signal sequences stochastically enforce TCRβ allelic exclusion

The monoallelic expression of antigen receptor (AgR) genes, called allelic exclusion, is fundamental for highly specific immune responses to pathogens. This cardinal feature of adaptive immunity is achieved by the assembly of a functional AgR gene on one allele, with subsequent feedback inhibition of V(D)J recombination on the other allele. A range of epigenetic mechanisms have been implicated in sequential recombination of AgR alleles; however, we now demonstrate that a genetic mechanism controls this process for Tcrb. Replacement of V(D)J recombinase targets at two different mouse Vβ gene segments with a higher quality target elevates Vβ rearrangement frequency before feedback inhibition, dramatically increasing the frequency of T cells with TCRβ chains derived from both Tcrb alleles. Thus, TCRβ allelic exclusion is enforced genetically by the low quality of Vβ recombinase targets that stochastically restrict the production of two functional rearrangements before feedback inhibition silences one allele.

The Usage of Human IGHJ Genes Follows a Particular Non-random Selection: The Recombination Signal Sequence May Affect the Usage of Human IGHJ Genes

The formation of the B cell receptor (BCR) heavy chain variable region is derived from the germline V(D)J gene rearrangement according to the “12/23” rule and the “beyond 12/23” rule. The usage frequency of each V(D)J gene in the peripheral BCR repertoires is related to the initial recombination, self-tolerance selection, and the clonal proliferative response. However, their specific differences and possible mechanisms are still unknown. We analyzed in-frame and out-of-frame BCR-H repertoires from human samples with normal physiological and various pathological conditions by high-throughput sequencing. Our results showed that IGHJ gene frequency follows a similar pattern which is previously known, where IGHJ4 is used at high frequency (>40%), IGHJ6/IGHJ3/IGHJ5 is used at medium frequencies (10∼20%), and IGH2/IGHJ1 is used at low frequency (<4%) under whether normal physiological or various pathological conditions. However, our analysis of the recombination signal sequences suggested that the conserved non-amer and heptamer and certain 23 bp spacer length may affect the initial IGHD-IGHJ recombination, which results in different frequencies of IGHJ genes among the initial BCR-H repertoire. Based on this “initial repertoire,” we recommend that re-evaluation and further investigation are needed when analyzing the significance and mechanism of IGHJ gene frequency in self-tolerance selection and the clonal proliferative response.

Sequence-dependent dynamics of synthetic and endogenous RSSs in V(D)J recombination

Developing lymphocytes of jawed vertebrates cleave and combine distinct gene segments to assemble antigen-receptor genes. This process called V(D)J recombination that involves the RAG recombinase binding and cutting recombination signal sequences (RSSs) composed of conserved heptamer and nonamer sequences flanking less well-conserved 12- or 23-bp spacers. Little quantitative information is known about the contributions of individual RSS positions over the course of the RAG-RSS interaction. We employ a single-molecule method known as tethered particle motion to track the formation, lifetime and cleavage of individual RAG-12RSS-23RSS paired complexes (PCs) for numerous synthetic and endogenous 12RSSs. We reveal that single-bp changes, including in the 12RSS spacer, can significantly and selectively alter PC formation or the probability of RAG-mediated cleavage in the PC. We find that some rarely used endogenous gene segments can be mapped directly to poor RAG binding on their adjacent 12RSSs. Finally, we find that while abrogating RSS nicking with Ca2+ leads to substantially shorter PC lifetimes, analysis of the complete lifetime distributions of any 12RSS even on this reduced system reveals that the process of exiting the PC involves unidentified molecular details whose involvement in RAG-RSS dynamics are crucial to quantitatively capture kinetics in V(D)J recombination.

Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes

The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.

Cutting antiparallel DNA strands in a single active site

A single enzyme active site that catalyzes multiple reactions is a well-established biochemical theme, but how one nuclease site cleaves both DNA strands of a double helix has not been well understood. In analyzing site-specific DNA cleavage by the mammalian RAG1-RAG2 recombinase, which initiates V(D)J recombination, we find that the active site is reconfigured for the two consecutive reactions and the DNA double helix adopts drastically different structures. For initial nicking of the DNA, a locally unwound and unpaired DNA duplex forms a zipper via alternating interstrand base stacking, rather than melting as generally thought. The second strand cleavage and formation of a hairpin-DNA product requires a global scissor-like movement of protein and DNA, delivering the scissile phosphate into the rearranged active site.

Functional requirement of terminal inverted repeats for efficient ProtoRAG activity reveals the early evolution of V(D)J recombination

The discovery of ProtoRAG in amphioxus indicated that vertebrate RAG recombinases originated from an ancient transposon. However, the sequences of ProtoRAG terminal inverted repeats (TIRs) were obviously dissimilar to the consensus sequence of mouse 12/23RSS and recombination mediated by ProtoRAG or RAG made them incompatible with each other. Thus, it is difficult to determine whether or how 12/23RSS persisted in the vertebrate RAG system that evolved from the TIRs of ancient RAG transposons. Here, we found that the activity of ProtoRAG is highly dependent on its asymmetric 5′TIR and 3′TIR, which are composed of conserved TR1 and TR5 elements and a partially conserved TRsp element of 27/31 bp to separate them. Similar to the requirements for the recombination signal sequences (RSSs) of RAG recombinase, the first CAC in TR1, the three dinucleotides in TR5 and the specific length of the partially conserved TRsp are important for the efficient recombination activity of ProtoRAG. In addition, the homologous sequences flanking the signal sequences facilitate ProtoRAG- but not RAG-mediated recombination. In addition to the diverged TIRs, two differentiated functional domains in BbRAG1L were defined to coordinate with the divergence between TIRs and RSSs. One of these is the CTT* domain, which facilitates the specific TIR recognition of the BbRAGL complex, and the other is NBD*, which is responsible for DNA binding and the protein stabilization of the BbRAGL complex. Thus, our findings reveal that the functional requirement for ProtoRAG TIRs is similar to that for RSS in RAG-mediated recombination, which not only supports the common origin of ProtoRAG TIRs and RSSs from the asymmetric TIRs of ancient RAG transposons, but also reveals the development of RAG and RAG-like machineries during chordate evolution.

Identification of recurrent noncoding mutations in B-cell lymphoma using capture Hi-C

The identification of driver mutations is fundamental to understanding oncogenesis. Although genes frequently mutated in B-cell lymphoma have been identified, the search for driver mutations has largely focused on the coding genome. Here we report an analysis of the noncoding genome using whole-genome sequencing data from 117 patients with B-cell lymphoma. Using promoter capture Hi-C data in naive B cells, we define cis-regulatory elements, which represent an enriched subset of the noncoding genome in which to search for driver mutations. Regulatory regions were identified whose mutation significantly alters gene expression, including copy number variation at cis-regulatory elements targeting CD69, IGLL5, and MMP14, and single nucleotide variants in a cis-regulatory element for TPRG1 We also show the commonality of pathways targeted by coding and noncoding mutations, exemplified by MMP14, which regulates Notch signaling, a pathway important in lymphomagenesis and whose expression is associated with patient survival. This study provides an enhanced understanding of lymphomagenesis and describes the advantages of using chromosome conformation capture to decipher noncoding mutations relevant to cancer biology.

Topology and expressed repertoire of the Felis catus T cell receptor loci

Background

The domestic cat (Felis catus) is an important companion animal and is used as a large animal model for human disease. However, the comprehensive study of adaptive immunity in this species is hampered by the lack of data on lymphocyte antigen receptor genes and usage. The objectives of this study were to annotate the feline T cell receptor (TR) loci and to characterize the expressed repertoire in lymphoid organs of normal cats using high-throughput sequencing.

Results

The Felis catus TRG locus contains 30 genes: 12 TRGV, 12 TRGJ and 6 TRGC, the TRB locus contains 48 genes: 33 TRBV, 2 TRBD, 11 TRBJ, 2 TRBC, the TRD locus contains 19 genes: 11 TRDV, 2 TRDD, 5 TRDJ, 1 TRDC, and the TRA locus contains 127 genes: 62 TRAV, 64 TRAJ, 1 TRAC. Functional feline V genes form monophyletic clades with their orthologs, and clustering of multimember subgroups frequently occurs in V genes located at the 5′ end of TR loci. Recombination signal (RS) sequences of the heptamer and nonamer of functional V and J genes are highly conserved. Analysis of the TRG expressed repertoire showed preferential intra-cassette over inter-cassette rearrangements and dominant usage of the TRGV2–1 and TRGJ1–2 genes. The usage of TRBV genes showed minor bias but TRBJ genes of the second J-C-cluster were more commonly rearranged than TRBJ genes of the first cluster. The TRA/TRD V genes almost exclusively rearranged to J genes within their locus. The TRAV/TRAJ gene usage was relatively balanced while the TRD repertoire was dominated by TRDJ3.

Conclusions

This is the first description of all TR loci in the cat. The genomic organization of feline TR loci was similar to that of previously described jawed vertebrates (gnathostomata) and is compatible with the birth-and-death model of evolution. The large-scale characterization of feline TR genes provides comprehensive baseline data on immune repertoires in healthy cats and will facilitate the development of improved reagents for the diagnosis of lymphoproliferative diseases in cats. In addition, these data might benefit studies using cats as a large animal model for human disease.

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.

Cut-and-Run: A Distinct Mechanism by which V(D)J Recombination Causes Genome Instability

V(D)J recombination is essential to generate antigen receptor diversity but is also a potent cause of genome instability. Many chromosome alterations that result from aberrant V(D)J recombination involve breaks at single recombination signal sequences (RSSs). A long-standing question, however, is how such breaks occur. Here, we show that the genomic DNA that is excised during recombination, the excised signal circle (ESC), forms a complex with the recombinase proteins to efficiently catalyze breaks at single RSSs both in vitro and in vivo. Following cutting, the RSS is released while the ESC-recombinase complex remains intact to potentially trigger breaks at further RSSs. Consistent with this, chromosome breaks at RSSs increase markedly in the presence of the ESC. Notably, these breaks co-localize with those found in acute lymphoblastic leukemia patients and occur at key cancer driver genes. We have named this reaction “cut-and-run” and suggest that it could be a significant cause of lymphocyte genome instability.

•

A complex between the recombination by-product and RAGs triggers multiple DNA breaks

•

The breaks co-localize with chromosome breakpoints in acute lymphoblastic leukemias

•

The breaks occur at many frequently mutated genes in acute lymphoblastic leukemia

•

Cut-and-run may underpin the most common types of lymphocyte chromosome instabilities

DNA damage pathways and B-cell lymphomagenesis

PURPOSE OF REVIEW

Recent lymphoma genome sequencing projects have shed light on the genomic landscape of indolent and aggressive lymphomas, as well as some of the molecular mechanisms underlying recurrent mutations and translocations in these entities. Here, we review these recent genomic discoveries, focusing on acquired DNA repair defects in lymphoma. In addition, we highlight recently identified actionable molecular vulnerabilities associated with recurrent mutations in chronic lymphocytic leukemia (CLL), which serves as a model entity.

RECENT FINDINGS

The results of several large lymphoma genome sequencing projects have recently been reported, including CLL, T-PLL and DLBCL. We align these discoveries with proposed mechanisms of mutation acquisition in B-cell lymphomas. Moreover, novel autochthonous mouse models of CLL have recently been generated and we discuss how these models serve as preclinical tools to drive the development of novel targeted therapeutic interventions. Lastly, we highlight the results of early clinical data on novel compounds targeting defects in the DNA damage response of CLL with a particular focus on deleterious ATM mutations.

SUMMARY

Defects in DNA repair pathways are selected events in cancer, including lymphomas. Specifically, ATM deficiency is associated with PARP1- and DNA-PKcs inhibitor sensitivity in vitro and in vivo.

A novel RAG1 mutation reveals a critical in vivo role for HMGB1/2 during V(D)J recombination

The Recombination Activating Genes, RAG1 and RAG2, are essential for V(D)J recombination and adaptive immunity. Mutations in these genes often cause immunodeficiency, the severity of which reflects the importance of the altered residue or residues during recombination. Here, we describe a novel RAG1 mutation that causes immunodeficiency in an unexpected way: The mutated protein severely disrupts binding of the accessory protein, HMGB1. Although HMGB1 enhances RAG cutting in vitro, its role in vivo was controversial. We show here that reduced HMGB1 binding by the mutant protein dramatically reduces RAG cutting in vitro and almost completely eliminates recombination in vivo. The RAG1 mutation, R401W, places a bulky tryptophan opposite the binding site for HMG Box A at both 12- and 23-spacer recombination signal sequences, disrupting stable binding of HMGB1. Replacement of R401W with leucine and then lysine progressively restores HMGB1 binding, correlating with increased RAG cutting and recombination in vivo. We show further that knockdown of HMGB1 significantly reduces recombination by wild-type RAG1, whereas its re-addition restores recombination with wild-type, but not the mutant, RAG1 protein. Together, these data provide compelling evidence that HMGB1 plays a critical role during V(D)J recombination in vivo.

The subclonal complexity of STIL-TAL1+ T-cell acute lymphoblastic leukaemia

Single-cell genetics were used to interrogate clonal complexity and the sequence of mutational events in STIL-TAL1+ T-ALL. Single-cell multicolour FISH was used to demonstrate that the earliest detectable leukaemia subclone contained the STIL-TAL1 fusion and copy number loss of 9p21.3 (CDKN2A/CDKN2B locus), with other copy number alterations including loss of PTEN occurring as secondary subclonal events. In three cases, multiplex qPCR and phylogenetic analysis were used to produce branching evolutionary trees recapitulating the snapshot history of T-ALL evolution in this leukaemia subtype, which confirmed that mutations in key T-ALL drivers, including NOTCH1 and PTEN, were subclonal and reiterative in distinct subclones. Xenografting confirmed that self-renewing or propagating cells were genetically diverse. These data suggest that the STIL-TAL1 fusion is a likely founder or truncal event. Therapies targeting the TAL1 auto-regulatory complex are worthy of further investigation in T-ALL.

DNA melting initiates the RAG catalytic pathway

The mechanism for initiating DNA cleavage by DDE-family enzymes, including the RAG endonuclease, which initiates V(D)J recombination, is not well understood. Here we report six cryo-EM structures of zebrafish RAG in complex with one or two intact recombination signal sequences (RSSs), at up to 3.9-Å resolution. Unexpectedly, these structures reveal DNA melting at the heptamer of the RSSs, thus resulting in a corkscrew-like rotation of coding-flank DNA and the positioning of the scissile phosphate in the active site. Substrate binding is associated with dimer opening and a piston-like movement in RAG1, first outward to accommodate unmelted DNA and then inward to wedge melted DNA. These precleavage complexes show limited base-specific contacts of RAG at the conserved terminal CAC/GTG sequence of the heptamer, thus suggesting conservation based on a propensity to unwind. CA and TG overwhelmingly dominate terminal sequences in transposons and retrotransposons, thereby implicating a universal mechanism for DNA melting during the initiation of retroviral integration and DNA transposition.

Regulated complex assembly safeguards the fidelity of Sleeping Beauty transposition

The functional relevance of the inverted repeat structure (IR/DR) in a subgroup of the Tc1/mariner superfamily of transposons has been enigmatic. In contrast to mariner transposition, where a topological filter suppresses single-ended reactions, the IR/DR orchestrates a regulatory mechanism to enforce synapsis of the transposon ends before cleavage by the transposase occurs. This ordered assembly process shepherds primary transposase binding to the inner 12DRs (where cleavage does not occur), followed by capture of the 12DR of the other transposon end. This extra layer of regulation suppresses aberrant, potentially genotoxic recombination activities, and the mobilization of internally deleted copies in the IR/DR subgroup, including Sleeping Beauty (SB). In contrast, internally deleted sequences (MITEs) are preferred substrates of mariner transposition, and this process is associated with the emergence of Hsmar1-derived miRNA genes in the human genome. Translating IR/DR regulation to in vitro evolution yielded an SB transposon version with optimized substrate recognition (pT4). The ends of SB transposons excised by a K248A excision⁺/integration^- transposase variant are processed by hairpin resolution, representing a link between phylogenetically, and mechanistically different recombination reactions, such as V(D)J recombination and transposition. Such variants generated by random mutation might stabilize transposon-host interactions or prepare the transposon for a horizontal transfer.

The ATM Kinase Restrains Joining of Both VDJ Signal and Coding Ends

The evidence that ATM affects resolution of RAG-induced DNA double-strand breaks is profuse and unequivocal; moreover, it is clear that the RAG complex itself cooperates (in an undetermined way) with ATM to facilitate repair of these double-strand breaks by the classical nonhomologous end-joining pathway. The mechanistic basis for the cooperation between ATM and the RAG complex has not been defined, although proposed models invoke ATM and RAG2's C terminus in maintaining the RAG postcleavage complex. In this study, we show that ATM reduces the rate of both coding and signal joining in a robust episomal assay; we suggest that this is the result of increased stability of the postcleavage complex. ATM's ability to inhibit VDJ joining requires its enzymatic activity. The noncore C termini of both RAG1 and RAG2 are also required for ATM's capacity to limit signal (but not coding) joining. Moreover, potential phosphorylation targets within the C terminus of RAG2 are also required for ATM's capacity to limit signal joining. These data suggest a model whereby the RAG signal end complex is stabilized by phosphorylation of RAG2 by ATM.

Genomic and expression analyses of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) loci reveal a similar basic public γδ repertoire in dolphin and human

BACKGROUND

The bottlenose dolphin (Tursiops truncatus) is a mammal that belongs to the Cetartiodactyla and have lived in marine ecosystems for nearly 60 millions years. Despite its popularity, our knowledge about its adaptive immunity and evolution is very limited. Furthermore, nothing is known about the genomics and evolution of dolphin antigen receptor immunity.

RESULTS

Here we report a evolutionary and expression study of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) genes. We have identified in silico the TRG and TRA/TRD genes and analyzed the relevant mature transcripts in blood and in skin from four subjects. The dolphin TRG locus is the smallest and simplest of all mammalian loci as yet studied. It shows a genomic organization comprising two variable (V1 and V2), three joining (J1, J2 and J3) and a single constant (C), genes. Despite the fragmented nature of the genome assemblies, we deduced the TRA/TRD locus organization, with the recent TRDV1 subgroup genes duplications, as it is expected in artiodactyls. Expression analysis from blood of a subject allowed us to assign unambiguously eight TRAV genes to those annotated in the genomic sequence and to twelve new genes, belonging to five different subgroups. All transcripts were productive and no relevant biases towards TRAV-J rearrangements are observed. Blood and skin from four unrelated subjects expression data provide evidence for an unusual ratio of productive/unproductive transcripts which arise from the TRG V-J gene rearrangement and for a "public" gamma delta TR repertoire. The productive cDNA sequences, shared both in the same and in different individuals, include biases of the TRGV1 and TRGJ2 genes. The high frequency of TRGV1-J2/TRDV1- D1-J4 productive rearrangements in dolphins may represent an interesting oligo-clonal population comparable to that found in human with the TRGV9- JP/TRDV2-D-J T cells and in primates.

CONCLUSIONS

Although the features of the TRG and TRA/TRD loci organization reflect those of the so far examined artiodactyls, genomic results highlight in dolphin an unusually simple TRG locus. The cDNA analysis reveal productive TRA/TRD transcripts and unusual ratios of productive/unproductive TRG transcripts. Comparing multiple different individuals, evidence is found for a "public" gamma delta TCR repertoire thus suggesting that in dolphins as in human the gamma delta TCR repertoire is accompanied by selection for public gamma chain.

Discovery of an Active RAG Transposon Illuminates the Origins of V(D)J Recombination

Co-option of RAG1 and RAG2 for antigen receptor gene assembly by V(D)J recombination was a crucial event in the evolution of jawed vertebrate adaptive immunity. RAG1/2 are proposed to have arisen from a transposable element, but definitive evidence for this is lacking. Here, we report the discovery of ProtoRAG, a DNA transposon family from lancelets, the most basal extant chordates. A typical ProtoRAG is flanked by 5-bp target site duplications and a pair of terminal inverted repeats (TIRs) resembling V(D)J recombination signal sequences. Between the TIRs reside tail-to-tail-oriented, intron-containing RAG1-like and RAG2-like genes. We demonstrate that ProtoRAG was recently active in the lancelet germline and that the lancelet RAG1/2-like proteins can mediate TIR-dependent transposon excision, host DNA recombination, transposition, and low-efficiency TIR rejoining using reaction mechanisms similar to those used by vertebrate RAGs. We propose that ProtoRAG represents a molecular "living fossil" of the long-sought RAG transposon.

Collaboration of RAG2 with RAG1-like proteins during the evolution of V(D)J recombination

Here, Carmona et al. show that two ancestral RAG1 proteins, Transib transposase and purple sea urchin RAG1-like, have a latent ability to initiate V(D)J recombination when coexpressed with RAG2 and that in vitro transposition by Transib transposase is stimulated by RAG2. They propose that evolution of RAG1/RAG2 began with a Transib transposon whose intrinsic recombination activity was enhanced by capture of an ancestral RAG2, allowing for the development of adaptive immunity.

The recombination-activating gene 1 (RAG1) and RAG2 proteins initiate V(D)J recombination, the process that assembles the B- and T-lymphocyte antigen receptor genes of jawed vertebrates. RAG1 and RAG2 are thought to have arisen from a transposable element, but the origins of this element are not understood. We show that two ancestral RAG1 proteins, Transib transposase and purple sea urchin RAG1-like, have a latent ability to initiate V(D)J recombination when coexpressed with RAG2 and that in vitro transposition by Transib transposase is stimulated by RAG2. Conversely, we report low levels of V(D)J recombination by RAG1 in the absence of RAG2. Recombination by RAG1 alone differs from canonical V(D)J recombination in having lost the requirement for asymmetric DNA substrates, implicating RAG2 in the origins of the “12/23 rule,” a fundamental regulatory feature of the reaction. We propose that evolution of RAG1/RAG2 began with a Transib transposon whose intrinsic recombination activity was enhanced by capture of an ancestral RAG2, allowing for the development of adaptive immunity.

Restoration of ATM Expression in DNA-PKcs-Deficient Cells Inhibits Signal End Joining

Unlike most DNA-dependent protein kinase, catalytic subunit (DNA-PKcs)-deficient mouse cell strains, we show in the present study that targeted deletion of DNA-PKcs in two different human cell lines abrogates VDJ signal end joining in episomal assays. Although the mechanism is not well defined, DNA-PKcs deficiency results in spontaneous reduction of ATM expression in many cultured cell lines (including those examined in this study) and in DNA-PKcs-deficient mice. We considered that varying loss of ATM expression might explain differences in signal end joining in different cell strains and animal models, and we investigated the impact of ATM and/or DNA-PKcs loss on VDJ recombination in cultured human and rodent cell strains. To our surprise, in DNA-PKcs-deficient mouse cell strains that are proficient in signal end joining, restoration of ATM expression markedly inhibits signal end joining. In contrast, in DNA-PKcs-deficient cells that are deficient in signal end joining, complete loss of ATM enhances signal (but not coding) joint formation. We propose that ATM facilitates restriction of signal ends to the classical nonhomologous end-joining pathway.

V(D)J Recombination: Mechanism, Errors, and Fidelity

V(D)J recombination, the mechanism responsible for generating antigen receptor diversity, has the potential to generate aberrant DNA rearrangements in developing lymphocytes. Indeed, the recombinase has been implicated in several different kinds of errors leading to oncogenic transformation. Here we review the basic aspects of V(D)J recombination, mechanisms underlying aberrant DNA rearrangements, and the types of aberrant events uncovered in recent genomewide analyses of lymphoid neoplasms.

Sheep (Ovis aries) T cell receptor alpha (TRA) and delta (TRD) genes and genomic organization of the TRA/TRD locus

Background

In mammals, T cells develop along two discrete pathways characterized by expression of either the αβ or the γδ T cell receptors. Human and mouse display a low peripheral blood γδ T cell percentage ("γδ low species") while sheep, bovine and pig accounts for a high proportion of γδ T lymphocytes ("γδ high species"). While the T cell receptor alpha (TRA) and delta (TRD) genes and the genomic organization of the TRA/TRD locus has been determined in human and mouse, this information is still poorly known in artiodactyl species, such as sheep.

Results

The analysis of the current Ovis aries whole genome assembly, Oar_v3.1, revealed that, as in the other mammalian species, the sheep TRD locus is nested within the TRA locus. In the most 5’ part the TRA/TRD locus contains TRAV genes which are intermingled with TRDV genes, then TRD genes which include seven TRDD, four TRDJ genes, one TRDC and a single TRDV gene with an inverted transcriptional orientation, and finally in the most 3’ part, the TRA locus is completed by 61 TRAJ genes and one TRAC gene.

Comparative sequence and analysis and annotation led to the identification of 66 TRAV genes assigned to 34 TRAV subgroups and 25 TRDV genes belonging to the TRDV1 subgroup, while one gene was found for each TRDV2, TRDV3 and TRDV4 subgroups. Multiple duplication events within several TRAV subgroups have generated the sheep TRAV germline repertoire, which is substantially larger than the human one. A significant proportion of these TRAV gene duplications seems to have occurred simultaneously with the amplification of the TRDV1 subgroup genes. This dynamic of expansion has also generated novel multigene subgroups, which are species-specific. Ovis aries TRA and TRD genes identified in this study were assigned IMGT definitive or temporary names and were approved by the IMGT/WHO-IUIS nomenclature committee.

The completeness of the genome assembly in the 3' part of the locus has allowed us to interpret rearranged CDR3 of cDNA from both TRA and TRD chain repertoires. The involvement of one up to four TRDD genes into a single transcript makes the potential sheep TRD chain much larger than any known TR chain repertoire.

Conclusions

The sheep genome, as the bovine genome, contains a large and diverse repertoire of TRA and TRD genes when compared to the “γδ T cell low” species genomes. The composition and length of the rearranged CDR3 in TRD V-delta domains influence the three-dimensional configuration of the antigen-combining site thus suggesting that in ruminants, γδ T cells play a more important and specific role in immune recognition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1790-z) contains supplementary material, which is available to authorized users.

Recruitment of RAG1 and RAG2 to Chromatinized DNA during V(D)J Recombination

V(D)J recombination is initiated by the binding of the RAG1 and RAG2 proteins to recombination signal sequences (RSSs) that consist of conserved heptamer and nonamer sequences separated by a spacer of either 12 or 23 bp. Here, we used RAG-inducible pro-B v-Abl cell lines in conjunction with chromatin immunoprecipitation to better understand the protein and RSS requirements for RAG recruitment to chromatin. Using a catalytic mutant form of RAG1 to prevent recombination, we did not observe cooperation between RAG1 and RAG2 in their recruitment to endogenous Jκ gene segments over a 48-h time course. Using retroviral recombination substrates, we found that RAG1 was recruited inefficiently to substrates lacking an RSS or containing a single RSS, better to substrates with two 12-bp RSSs (12RSSs) or two 23-bp RSSs (23RSSs), and more efficiently to a substrate with a 12/23RSS pair. RSS mutagenesis demonstrated a major role for the nonamer element in RAG1 binding, and correspondingly, a cryptic RSS consisting of a repeat of CA dinucleotides, which poorly re-creates the nonamer, was ineffective in recruiting RAG1. Our findings suggest that 12RSS-23RSS cooperation (the "12/23 rule") is important not only for regulating RAG-mediated DNA cleavage but also for the efficiency of RAG recruitment to chromatin.

Regulation and Evolution of the RAG Recombinase

The modular, noncontiguous architecture of the antigen receptor genes necessitates their assembly through V(D)J recombination. This program of DNA breakage and rejoining occurs during early lymphocyte development, and depends on the RAG1 and RAG2 proteins, whose collaborative endonuclease activity targets specific DNA motifs enriched in the antigen receptor loci. This essential gene shuffling reaction requires lymphocytes to traverse several developmental stages wherein DNA breakage is tolerated, while minimizing the expense to overall genome integrity. Thus, RAG activity is subject to stringent temporal and spatial regulation. The RAG proteins themselves also contribute autoregulatory properties that coordinate their DNA cleavage activity with target chromatin structure, cell cycle status, and DNA repair pathways. Even so, lapses in regulatory restriction of RAG activity are apparent in the aberrant V(D)J recombination events that underlie many lymphomas. In this review, we discuss the current understanding of the RAG endonuclease, its widespread binding in the lymphocyte genome, its noncleavage activities that restrain its enzymatic potential, and the growing evidence of its evolution from an ancient transposase.

Gene analysis of signal-joint T cell receptor excision circles and their relationship to age in dogs

The quantification of DNA excision circles produced during T cell receptor (TCR) rearrangement, termed signal joint TCR rearrangement excision circles (sjTRECs), has been employed as a measure of age and thymic function in humans and animals. δRec-ψJα sjTRECs are ring-shaped DNAs that are generated during TCRδ locus deletion that occurs at a late stage of T cell development. In this study, the nucleotide sequences of δRec-ψJα signal joints of canine δRec-ψJα sjTRECs were analyzed. The gene structure of canine δRec-ψJα signal joints was found to be similar to that of humans and mice. However, diversity of signal joints was detected and found to derive from N nucleotide insertions, recombination signal sequence combinational diversity and single-base substitutions at the recombination signal sequence. In addition, an adenine insertion or deletion was found approximately 280 bases from the ψJα signal end. Blood samples were collected from 46 dogs, ranging in age from 3 to 192 months, with a mean age of 96.4 and a SD of 51.5 months. Although δRec-ψJα sjTRECs were detectable in most of the dogs evaluated, the level did not significantly correlate with age. These results indicated that δRec-ψJα sjTREC levels were ineffective as a measure of age in dogs.

Assembly Pathway and Characterization of the RAG1/2-DNA Paired and Signal-end Complexes

Mammalian immune receptor diversity is established via a unique restricted set of site-specific DNA rearrangements in lymphoid cells, known as V(D)J recombination. The lymphoid-specific RAG1-RAG2 protein complex (RAG1/2) initiates this process by binding to two types of recombination signal sequences (RSS), 12RSS and 23RSS, and cleaving at the boundaries of RSS and V, D, or J gene segments, which are to be assembled into immunoglobulins and T-cell receptors. Here we dissect the ordered assembly of the RAG1/2 heterotetramer with 12RSS and 23RSS DNAs. We find that RAG1/2 binds only a single 12RSS or 23RSS and reserves the second DNA-binding site specifically for the complementary RSS, to form a paired complex that reflects the known 12/23 rule of V(D)J recombination. The assembled RAG1/2 paired complex is active in the presence of Mg(2+), the physiologically relevant metal ion, in nicking and double-strand cleavage of both RSS DNAs to produce a signal-end complex. We report here the purification and initial crystallization of the RAG1/2 signal-end complex for atomic-resolution structure elucidation. Strict pairing of the 12RSS and 23RSS at the binding step, together with information from the crystal structure of RAG1/2, leads to a molecular explanation of the 12/23 rule.

Genomic analysis offers insights into the evolution of the bovine TRA/TRD locus

Background

The TRA/TRD locus contains the genes for V(D)J somatic rearrangement of TRA and TRD chains expressed by αβ and γδ T cells respectively. Previous studies have demonstrated that the bovine TRA/TRD locus contains an exceptionally large number of TRAV/TRDV genes. In this study we combine genomic and transcript analysis to provide insights into the evolutionary development of the bovine TRA/TRD locus and the remarkable TRAV/TRDV gene repertoire.

Results

Annotation of the UMD3.1 assembly identified 371 TRAV/TRDV genes (distributed in 42 subgroups), 3 TRDJ, 6 TRDD, 62 TRAJ and single TRAC and TRDC genes, most of which were located within a 3.5 Mb region of chromosome 10. Most of the TRAV/TRDV subgroups have multiple members and several have undergone dramatic expansion, most notably TRDV1 (60 genes). Wide variation in the proportion of pseudogenes within individual subgroups, suggest that differential ‘birth’ and ‘death’ rates have been used to form a functional bovine TRAV/TRDV repertoire which is phylogenetically distinct from that of humans and mice. The expansion of the bovine TRAV/TRDV gene repertoire has predominantly been achieved through a complex series of homology unit (regions of DNA containing multiple gene) replications. Frequent co-localisation within homology units of genes from subgroups with low and high pseudogene proportions suggest that replication of homology units driven by evolutionary selection for the former may have led to a ‘collateral’ expansion of the latter. Transcript analysis was used to define the TRAV/TRDV subgroups available for recombination of TRA and TRD chains and demonstrated preferential usage of different subgroups by the expressed TRA and TRD repertoires, indicating that TRA and TRD selection have had distinct impacts on the evolution of the TRAV/TRDV repertoire.

Conclusion

Both TRA and TRD selection have contributed to the evolution of the bovine TRAV/TRDV repertoire. However, our data suggest that due to homology unit duplication TRD selection for TRDV1 subgroup expansion may have substantially contributed to the genomic expansion of several TRAV subgroups. Such data demonstrate how integration of genomic and transcript data can provide a more nuanced appreciation of the evolutionary dynamics that have led to the dramatically expanded bovine TRAV/TRDV repertoire.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-994) contains supplementary material, which is available to authorized users.

Multiple factors influence the contribution of individual immunoglobulin light chain genes to the naïve antibody repertoire

Background

The naïve antibody repertoire is initially dependent upon the number of germline V(D)J genes and the ability of recombined heavy and light chains to pair. Individual VH and VL genes are not equally represented in naïve mature B cells, suggesting that positive and negative selection also shape the antibody repertoire. Among the three member murine Vκ10 L chain family, the Vκ10C gene is under-represented in the antibody repertoire. Although it is structurally functional and accessible to both transcriptional and recombination machinery, the Vκ10C promoter is inefficient in pre-B cell lines and productive Vκ10C rearrangements are lost as development progresses from pre-B cells through mature B cells. This study examined VH/Vκ10 pairing, promoter mutations, Vκ10 transcript levels and receptor editing as possible factors that are responsible for loss of productive Vκ10C rearrangements in developing B cells.

Results

We demonstrate that the loss of Vκ10C expression is not due to an inability to pair with H chains, but is likely due to a combination of other factors. Levels of mRNA are low in sorted pre-B cells and undetectable in B cells. Mutation of a single base in the three prime region of the Vκ10C promoter increases Vκ10C promoter function in pre-B cell lines. Pre-B and B cells harbor disproportionate levels of receptor-edited productive Vκ10C rearrangements.

Conclusions

Our findings suggest that the weak Vκ10C promoter initially limits the amount of available Vκ10C L chain for pairing with H chains, resulting in sub-threshold levels of cell surface B cell receptors, insufficient tonic signaling and subsequent receptor editing to limit the numbers of Vκ10C-expressing B cells emigrating from the bone marrow to the periphery.

RAG-mediated recombination is the predominant driver of oncogenic rearrangement in ETV6-RUNX1 acute lymphoblastic leukemia

The ETV6-RUNX1 fusion gene, found in 25% of childhood acute lymphoblastic leukemia (ALL) cases, is acquired in utero but requires additional somatic mutations for overt leukemia. We used exome and low-coverage whole-genome sequencing to characterize secondary events associated with leukemic transformation. RAG-mediated deletions emerge as the dominant mutational process, characterized by recombination signal sequence motifs near breakpoints, incorporation of non-templated sequence at junctions, ∼30-fold enrichment at promoters and enhancers of genes actively transcribed in B cell development and an unexpectedly high ratio of recurrent to non-recurrent structural variants. Single-cell tracking shows that this mechanism is active throughout leukemic evolution, with evidence of localized clustering and reiterated deletions. Integration of data on point mutations and rearrangements identifies ATF7IP and MGA as two new tumor-suppressor genes in ALL. Thus, a remarkably parsimonious mutational process transforms ETV6-RUNX1-positive lymphoblasts, targeting the promoters, enhancers and first exons of genes that normally regulate B cell differentiation.

Unifying model for molecular determinants of the preselection Vβ repertoire

The primary antigen receptor repertoire is sculpted by the process of V(D)J recombination, which must strike a balance between diversification and favoring gene segments with specialized functions. The precise determinants of how often gene segments are chosen to complete variable region coding exons remain elusive. We quantified Vβ use in the preselection Tcrb repertoire and report relative contributions of 13 distinct features that may shape their recombination efficiencies, including transcription, chromatin environment, spatial proximity to their DβJβ targets, and predicted quality of recombination signal sequences (RSSs). We show that, in contrast to functional Vβ gene segments, all pseudo-Vβ segments are sequestered in transcriptionally silent chromatin, which effectively suppresses wasteful recombination. Importantly, computational analyses provide a unifying model, revealing a minimum set of five parameters that are predictive of Vβ use, dominated by chromatin modifications associated with transcription, but largely independent of precise spatial proximity to DβJβ clusters. This learned model-building strategy may be useful in predicting the relative contributions of epigenetic, spatial, and RSS features in shaping preselection V repertoires at other antigen receptor loci. Ultimately, such models may also predict how designed or naturally occurring alterations of these loci perturb the preselection use of variable gene segments.

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.

Exceptionally long CDR3H are not isotype restricted in bovine immunoglobulins

Exceptionally long third complementarity determining regions of the heavy chain (CDR3H) were previously described as a specificity of bovine IgG and IgM immunoglobulins. In addition, the genomic organization of the immunoglobulin heavy chain locus remains to be elucidated with a special focus on the number of variable segments (IGHV). By analyzing the variable regions according to the isotype-specific PCR using cDNA-PCR, we were able to prove the existence of exceptional long CDR3H in all bovine isotypes. The corresponding sequences of three distinct amplicons were grouped according to the length of the CDR3H. Sequences of CDR3H possessed 5 to 10, 12 to 31 or at least 48 amino acid residues. Long and mid-length CDR3H were composed of mainly hydrophilic amino acid residues, while short CDR3H also contained hydrophobic amino acid residues. All sequences with long CDR3H were related to the germline variable segment 10. Using the current genome assembly, Bos taurus NCBI build 6.1, the genomic organization of the bovine immunoglobulin heavy-chain locus was analyzed. A main locus was investigated on BTA21. Exons coding for variable, diversity, and joining segments, as well as for the constant regions of different isotypes, were also localized on BTA7, BTA8, and BTA20. Together with the information from unplaced contigs, 36 IGHV were detected of which 13 are putatively functional. Phylogenetic analysis revealed two bovine IGHV families (boVH1, boVH2). Thus, the existence of the two bovine families suggested was demonstrated, where boVH1 comprises all functional segments. This study substantially improves the understanding of the generation of immunoglobulin diversity in cattle.

JEC. Secondary mechanisms of diversification in the human antibody repertoire

V(D)J recombination and somatic hypermutation (SHM) are the primary mechanisms for diversification of the human antibody repertoire. These mechanisms allow for rapid humoral immune responses to a wide range of pathogenic challenges. V(D)J recombination efficiently generate a virtually limitless diversity through random recombination of variable (V), diversity (D), and joining (J) genes with diverse non-templated junctions between the selected gene segments. Following antigen stimulation, affinity maturation by SHM produces antibodies with refined specificity mediated by mutations typically focused in complementarity determining regions (CDRs), which form the bulk of the antigen recognition site. While V(D)J recombination and SHM are responsible for much of the diversity of the antibody repertoire, there are several secondary mechanisms that, while less frequent, make substantial contributions to antibody diversity including V(DD)J recombination (or D-D fusion), SHM-associated insertions and deletions, and affinity maturation and antigen contact by non-CDR regions of the antibody. In addition to enhanced diversity, these mechanisms allow the production of antibodies that are critical to response to a variety of viral and bacterial pathogens but that would be difficult to generate using only the primary mechanisms of diversification.

Role of recombination activating genes in the generation of antigen receptor diversity and beyond

V(D)J recombination is the process by which antibody and T-cell receptor diversity is attained. During this process, antigen receptor gene segments are cleaved and rejoined by non-homologous DNA end joining for the generation of combinatorial diversity. The major players of the initial process of cleavage are the proteins known as RAG1 (recombination activating gene 1) and RAG2. In this review, we discuss the physiological function of RAGs as a sequence-specific nuclease and its pathological role as a structure-specific nuclease. The first part of the review discusses the basic mechanism of V(D)J recombination, and the last part focuses on how the RAG complex functions as a sequence-specific and structure-specific nuclease. It also deals with the off-target cleavage of RAGs and its implications in genomic instability.