Novel bimodal TRBD1-TRBD2 rearrangements with dual or absent D-region contribute to TRB V-(D)-J combinatorial diversity

T-cell receptor (TR) diversity of the variable domains is generated by recombination of both the alpha (TRA) and beta (TRB) chains. The textbook process of TRB chain production starts with TRBD and TRBJ gene rearrangement, followed by the rearrangement of a TRBV gene to the partially rearranged D-J gene. Unsuccessful V-D-J TRB rearrangements lead to apoptosis of the cell. Here, we performed deep sequencing of the poorly explored pool of partial TRBD1-TRBD2 rearrangements in T-cell genomic DNA. We reconstructed full repertoires of human partial TRBD1-TRBD2 rearrangements using novel sequencing and validated them by detecting V-D-J recombination-specific byproducts: excision circles containing the recombination signal (RS) joint 5'D2-RS - 3'D1-RS. Identified rearrangements were in compliance with the classical 12/23 rule, common for humans, rats, and mice and contained typical V-D-J recombination footprints. Interestingly, we detected a bimodal distribution of D-D junctions indicating two active recombination sites producing long and short D-D rearrangements. Long TRB D-D rearrangements with two D-regions are coding joints D1-D2 remaining classically on the chromosome. The short TRB D-D rearrangements with no D-region are signal joints, the coding joint D1-D2 being excised from the chromosome. They both contribute to the TRB V-(D)-J combinatorial diversity. Indeed, short D-D rearrangements may be followed by direct V-J2 recombination. Long D-D rearrangements may recombine further with J2 and V genes forming partial D1-D2-J2 and then complete V-D1-D2-J2 rearrangement. Productive TRB V-D1-D2-J2 chains are present and expressed in thousands of clones of human antigen-experienced memory T cells proving their capacity for antigen recognition and actual participation in the immune response.

The immunogenicity of human-origin therapeutic antibodies are associated with V gene usage

Therapeutic antibodies can elicit unwanted immune responses in a subset of patients, which leads to the production of anti-drug antibodies (ADA). Some of these ADAs have been reported to effect the pharmacokinetics, efficacy and/or safety of the therapeutic antibodies. The sequence diversity of antibodies are generated by VDJ recombination and mutagenesis. While the antibody generation process can create a large candidate pool for identifying high-affinity antibodies, it also could produce sequences that are foreign to the human immune system. However, it is not clear how VDJ recombination and mutagenesis impact the clinical ADA rate of therapeutic antibodies. In this study, we identified a positive correlation between the clinical ADA rate and the number of introduced mutations in the antibody sequences. We also found that the use of rare V alleles in human-origin antibody therapeutics is associated with higher risk of immunogenicity. The results suggest that antibody engineering projects should start with frameworks that contain commonly used V alleles and prioritize antibody candidates with low number of mutations to reduce the risk of immunogenicity.

Methods for Study of Mouse T Cell Receptor α and β Gene Rearrangements

Quantitative real-time PCR and next-generation sequencing (NGS) are invaluable techniques to analyze T cell receptor (Tcr) gene rearrangements in mouse lymphocyte populations. Although these approaches are powerful, they also have limitations that must be accounted for in experimental design and data interpretation. Here, we provide relevant background required for understanding these limitations and then outline established quantitative real-time PCR and NGS methods that can be used for analysis of mouse Tcra and Tcrb gene rearrangements in mice.

Helix-Loop-Helix Proteins in Adaptive Immune Development

The E/ID protein axis is instrumental for defining the developmental progression and functions of hematopoietic cells. The E proteins are dimeric transcription factors that activate gene expression programs and coordinate changes in chromatin organization. Id proteins are antagonists of E protein activity. Relative levels of E/Id proteins are modulated throughout hematopoietic development to enable the progression of hematopoietic stem cells into multiple adaptive and innate immune lineages including natural killer cells, B cells and T cells. In early progenitors, the E proteins promote commitment to the T and B cell lineages by orchestrating lineage specific programs of gene expression and regulating VDJ recombination of antigen receptor loci. In mature B cells, the E/Id protein axis functions to promote class switch recombination and somatic hypermutation. E protein activity further regulates differentiation into distinct CD4+ and CD8+ T cells subsets and instructs mature T cell immune responses. In this review, we discuss how the E/Id proteins define the adaptive immune system lineages, focusing on their role in directing developmental gene programs.

Enhancer dependent repositioning of TCRb locus with respect to the chromosome territory

Intranuclear position of several genes is dynamically altered during development concordant with their activation. To understand this dynamic, but non-random, nuclear organization, it is important to identify the relevant regulatory elements and trans acting factors. Murine TCRb locus gets activated during thymic development. Enhancer Eb is important for VDJ recombination at TCRb locus as it is critically required establishment of recombination center. Our analysis revealed that TCRb locus gets located out of the chromosome territory specifically in developing thymocytes. Further, CRISPR/Cas9 based deletion mutagenesis established an unambiguous role of enhancer Eb in defining TCRb location relative to chromosome territory. The ability to reposition the target locus relative to chromosome territory highlights a novel aspect pertaining to activity of enhancers which may contribute to their ability to regulate gene expression. Additionally, our observations have implications for understanding the role of enhancers in three-dimensional genome organization and function.

Antibody Repertoire Analysis of Tumor-Infiltrating B Cells Reveals Distinct Signatures and Distributions Across Tissues

The role of B cells in the tumor microenvironment (TME) has largely been under investigated, and data regarding the antibody repertoire encoded by B cells in the TME and the adjacent lymphoid organs are scarce. Here, we utilized B cell receptor high-throughput sequencing (BCR-Seq) to profile the antibody repertoire signature of tumor-infiltrating lymphocyte B cells (TIL−Bs) in comparison to B cells from three anatomic compartments in a mouse model of triple-negative breast cancer. We found that TIL-Bs exhibit distinct antibody repertoire measures, including high clonal polarization and elevated somatic hypermutation rates, suggesting a local antigen-driven B-cell response. Importantly, TIL-Bs were highly mutated but non-class switched, suggesting that class-switch recombination may be inhibited in the TME. Tracing the distribution of TIL-B clones across various compartments indicated that they migrate to and from the TME. The data thus suggests that antibody repertoire signatures can serve as indicators for identifying tumor-reactive B cells.

Regulation of Early Lymphocyte Development via mRNA Decay Catalyzed by the CCR4-NOT Complex

Development of lymphocytes is precisely regulated by various mechanisms. In addition to transcriptional rates, post-transcriptional regulation of mRNA abundance contributes to differentiation of lymphocytes. mRNA decay is a post-transcriptional mechanism controlling mRNA abundance. The carbon catabolite repression 4 (CCR4)-negative on TATA-less (NOT) complex controls mRNA longevity by catalyzing mRNA deadenylation, which is the rate-limiting step in the mRNA decay pathway. mRNA decay, regulated by the CCR4-NOT complex, is required for differentiation of pro-B to pre-B cells and V(D)J recombination in pro-B cells. In this process, it is likely that the RNA-binding proteins, ZFP36 ring finger protein like 1 and 2, recruit the CCR4-NOT complex to specific target mRNAs, thereby inducing cell quiescence of pro-B cells. A recent study showed that the CCR4-NOT complex participates in positive selection of thymocytes. Mechanistically, the CCR4-NOT deadenylase complex inhibits abnormal apoptosis by reducing the expression level of mRNAs encoding pro-apoptotic proteins, which are otherwise up-regulated during positive selection. We discuss mechanisms regulating CCR4-NOT complex-dependent mRNA decay in lymphocyte development and selection.

Igh Locus Polymorphism May Dictate Topological Chromatin Conformation and V Gene Usage in the Ig Repertoire

Vast repertoires of unique antigen receptors are created in developing B and T lymphocytes. The antigen receptor loci contain many variable (V), diversity (D) and joining (J) gene segments that are arrayed across very large genomic expanses and are joined to form variable-region exons of expressed immunoglobulins and T cell receptors. This process creates the potential for an organism to respond to large numbers of different pathogens. Here, we consider the possibility that genetic polymorphisms with alterations in a vast array of regulatory elements in the immunoglobulin heavy chain (IgH) locus lead to changes in locus topology and impact immune-repertoire formation.

RAG deficiencies: Recent advances in disease pathogenesis and novel therapeutic approaches

The RAG1 and RAG2 proteins initiate the process of V(D)J recombination and therefore play an essential role in adaptive immunity. While null mutations in the RAG genes cause severe combined immune deficiency with lack of T and B cells (T- B- SCID) and susceptibility to life-threatening, early-onset infections, studies in humans and mice have demonstrated that hypomorphic RAG mutations are associated with defects of central and peripheral tolerance resulting in immune dysregulation. In this review, we provide an overview of the extended spectrum of RAG deficiencies and their associated clinical and immunological phenotypes in humans. We discuss recent advances in the mechanisms that control RAG expression and function, the effects of perturbed RAG activity on lymphoid development and immune homeostasis, and propose novel approaches to correct this group of disorders.

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.

Impaired Expression of Rearranged Immunoglobulin Genes and Premature p53 Activation Block B Cell Development in BMI1 Null Mice

B cell development is a highly regulated process that requires stepwise rearrangement of immunoglobulin genes to generate a functional B cell receptor (BCR). The polycomb group protein BMI1 is required for B cell development, but its function in developing B cells remains poorly defined. We demonstrate that BMI1 functions in a cell-autonomous manner at two stages during early B cell development. First, loss of BMI1 results in a differentiation block at the pro-B cell to pre-B cell transition due to the inability of BMI1-deficient cells to transcribe newly rearranged Igh genes. Accordingly, introduction of a pre-rearranged Igh allele partially restored B cell development in Bmi1^−/− mice. In addition, BMI1 is required to prevent premature p53 signaling, and as a consequence, Bmi1^−/− large pre-B cells fail to properly proliferate. Altogether, our results clarify the role of BMI1 in early B cell development and uncover an unexpected function of BMI1 during VDJ recombination.

Cantor et al. identify a cell-autonomous role for the polycomb group protein BMI1 in early B cell development. At the pro-B cell to pre-B cell transition, BMI1 promotes the expression of newly rearranged Igh genes in pro-B cells and subsequently prevents premature p53 activation and enables large pre-B cell proliferation.

The T cell receptor (TRA) locus in the rabbit (Oryctolagus cuniculus): Genomic features and consequences for invariant T cells

The rabbit has been widely used in immunology and infectiology. Rabbit immunoglobulins have been extensively studied, leading to the discovery of their idiotypes, allotypic diversity, and of the diversification of the primary repertoire by hyperconversion. Much less is known about rabbit T cell receptors (TR), especially TRA. This isotype is particularly important for innate-like T cells, which typically express invariant TRA (iTRA). The presence of such cells in the rabbit remains an enigma. Rabbit NKT cells seem to be very rare, and lagomorphs lack MAIT cells. TRAV1, the variable gene expressed in the iTRA of these cells across most mammals, and MR1, the MH1-like receptor that present riboflavin derivatives to MAIT cells, are missing in rabbit. An alternative iTRA has been identified, that may be expressed by new innate-like T cells. To facilitate TRA repertoire analyses in rabbit, we report here a full description of TRA and TRD loci and a subgroup definition based on IMGT® classification. Rabbit TRA rearrangements follow the same temporal pattern that is observed in mouse and human. Rare transcripts expressing TRDV/TRDD/TRDJ rearrangements spliced to TRAC were detected. TRA and TRD genes have been made available in IMGT and IMGT/HighV-QUEST, allowing easy analysis of TRA/TRD RepSeq.

De novo Inference of Diversity Genes and Analysis of Non-canonical V(DD)J Recombination in Immunoglobulins

The V(D)J recombination forms the immunoglobulin genes by joining the variable (V), diversity (D), and joining (J) germline genes. Since variations in germline genes have been linked to various diseases, personalized immunogenomics aims at finding alleles of germline genes across various patients. Although recent studies described algorithms for de novo inference of V and J genes from immunosequencing data, they stopped short of solving a more difficult problem of reconstructing D genes that form the highly divergent CDR3 regions and provide the most important contribution to the antigen binding. We present the IgScout algorithm for de novo D gene reconstruction and apply it to reveal new alleles of human D genes and previously unknown D genes in camel, an important model organism in immunology. We further analyze non-canonical V(DD)J recombination that results in unusually long CDR3s with tandem fused IGHD genes and thus expands the diversity of the antibody repertoires. We demonstrate that tandem CDR3s represent a consistent and functional feature of all analyzed immunosequencing datasets, reveal ultra-long CDR3s, and shed light on the mechanism responsible for their formation.

Misregulation of the IgH Locus in Thymocytes

Functional antigen receptor genes are assembled by somatic rearrangements that are largely lymphocyte lineage specific. The immunoglobulin heavy chain (IgH) gene locus is unique amongst the seven antigen receptor loci in undergoing partial gene rearrangements in the wrong lineage. Here we demonstrate that breakdown of lineage-specificity is associated with inappropriate activation of the Eμ enhancer during T cell development by a different constellation of transcription factors than those used in developing B cells. This is reflected in reduced enhancer-induced epigenetic changes, eRNAs, formation of the RAG1/2-rich recombination center, attenuated chromatin looping and markedly different utilization of D_H gene segments in CD4⁺CD8⁺ (DP) thymocytes. Additionally, CTCF-dependent V_H locus compaction is disrupted in DP cells despite comparable transcription factor binding in both lineages. These observations identify multiple mechanisms that contribute to lineage-specific antigen receptor gene assembly.

Recombination activity of human recombination-activating gene 2 (RAG2) mutations and correlation with clinical phenotype

BACKGROUND

Mutations in recombination-activating gene (RAG) 1 and RAG2 are associated with a broad range of clinical and immunologic phenotypes in human subjects.

OBJECTIVE

Using a flow cytometry-based assay, we aimed to measure the recombinase activity of naturally occurring RAG2 mutant proteins and to correlate our results with the severity of the clinical and immunologic phenotype.

METHODS

Abelson virus-transformed Rag2^-/- pro-B cells engineered to contain an inverted green fluorescent protein (GFP) cassette flanked by recombination signal sequences were transduced with retroviruses encoding either wild-type or 41 naturally occurring RAG2 variants. Bicistronic vectors were used to introduce compound heterozygous RAG2 variants. The percentage of GFP-expressing cells was evaluated by using flow cytometry, and high-throughput sequencing was used to analyze rearrangements at the endogenous immunoglobulin heavy chain (Igh) locus.

RESULTS

The RAG2 variants showed a wide range of recombination activity. Mutations associated with severe combined immunodeficiency and Omenn syndrome had significantly lower activity than those detected in patients with less severe clinical presentations. Four variants (P253R, F386L, N474S, and M502V) previously thought to be pathogenic were found to have wild-type levels of activity. Use of bicistronic vectors permitted us to assess more carefully the effect of compound heterozygous mutations, with good correlation between GFP expression and the number and diversity of Igh rearrangements.

CONCLUSIONS

Our data support genotype-phenotype correlation in the setting of RAG2 deficiency. The assay described can be used to define the possible disease-causing role of novel RAG2 variants and might help predict the severity of the clinical phenotype.

αβ T-cell receptors with a central CDR3 cysteine are enriched in CD8αα intraepithelial lymphocytes and their thymic precursors

The thymus plays a crucial role in immune tolerance by exposing developing T cells (thymocytes) to a myriad of self-antigens. Strong T-cell receptor (TCR) engagement induces tolerance in self-reactive thymocytes by stimulating apoptosis or selection into specialized T-cell lineages, including intestinal TCRαβ⁺ CD8αα⁺ intraepithelial lymphocytes (IEL). TCR-intrinsic amino acid motifs that can be used to predict whether a TCR will be strongly self-reactive remain elusive. Here, a novel TCR sequence alignment approach revealed that T-cell lineages in C57BL/6 mice had divergent usage of cysteine within two positions of the amino acid at the apex of the complementarity-determining region 3 (CDR3) of the TCRα or TCRβ chain. Compared to pre-selection thymocytes, central CDR3 cysteine usage was increased in IEL and Type A IEL precursors (IELp) and markedly decreased in Foxp3⁺ regulatory T cells (T-reg) and naïve T cells. These findings reveal a TCR-intrinsic motif that distinguishes Type A IELp and IEL from T-reg and naïve T cells.

Circulating plasmablasts are elevated and produce pathogenic anti-endothelial cell autoantibodies in idiopathic pulmonary arterial hypertension

Idiopathic pulmonary arterial hypertension (IPAH) is a devastating pulmonary vascular disease in which autoimmune and inflammatory phenomena are implicated. B cells and autoantibodies have been associated with IPAH and identified as potential therapeutic targets. However, the specific populations of B cells involved and their roles in disease pathogenesis are not clearly defined. We aimed to assess the levels of activated B cells (plasmablasts) in IPAH, and to characterize recombinant antibodies derived from these plasmablasts. Blood plasmablasts are elevated in IPAH, remain elevated over time, and produce IgA autoantibodies. Single-cell sequencing of plasmablasts in IPAH revealed repertoires of affinity-matured antibodies with small clonal expansions, consistent with an ongoing autoimmune response. Recombinant antibodies representative of these clonal lineages bound known autoantigen targets and displayed an unexpectedly high degree of polyreactivity. Representative IPAH plasmablast recombinant antibodies stimulated human umbilical vein endothelial cells to produce cytokines and overexpress the adhesion molecule ICAM-1. Together, our results demonstrate an ongoing adaptive autoimmune response involving IgA plasmablasts that produce anti-endothelial cell autoantibodies in IPAH. These antibodies stimulate endothelial cell production of cytokines and adhesion molecules, which may contribute to disease pathogenesis. These findings suggest a role for mucosally-driven autoimmunity and autoimmune injury in the pathogenesis of IPAH.

Omenn Syndrome and DNA recombination defects

Mutations in the RAG1/RAG2 genes are associated with a broad spectrum of clinical phenotypes, ranging from severe combined immunodeficiency to various autoimmune diseases. The diversity of the clinical symptoms is determined not only by the residual RAG recombinase enzyme activity as determined by the mutations, but also by multiple environmental factors and, in rare cases, by second site mutations within the RAG1/RAG2 genes. The residual recombinase activity is responsible for the oligoclonal expansion of autoreactive T cells. Omenn syndrome is the result of intense Th2 type inflammation involving the skin and multiple other organs triggered by these T cells. In this review, the molecular pathology of diseases caused by RAG1/RAG2 mutations, in particular Omenn syndrome, will be discussed. Furthermore, abnormalities in other molecules involved in V(D)J recombination will be discussed in relation to Omenn-like syndrome.

Deacetylase activity of histone deacetylase 3 is required for productive VDJ recombination and B-cell development

Histone deacetylase 3 (HDAC3) is the catalytic component of NCoR/SMRT corepressor complexes that mediate the actions of transcription factors implicated in the regulation of B-cell development and function. We crossed Hdac3 conditional knockout mice with Mb1-Cre knockin animals to delete Hdac3 in early progenitor B cells. The spleens of Hdac3^F/-Mb1-Cre^+/- mice were virtually devoid of mature B cells, and B220⁺CD43⁺ B-cell progenitors accumulated within the bone marrow. Quantitative deep sequencing of the Ig heavy chain locus from B220⁺CD43⁺ populations identified a defect in V_HDJ_H recombination with a severe reduction in productive rearrangements, which directly corresponded to the loss of pre-B cells from Hdac3^Δ/- bone marrow. For Hdac3^Δ/- B cells that did show productive VDJ rearrangement, there was significant skewing toward the incorporation of proximal V_H gene segments and a corresponding reduction in distal V_H gene segment use. Although transcriptional effects within these loci were modest, Hdac3^Δ/- progenitor cells displayed global changes in chromatin structure that likely hindered effective distal V-DJ recombination. Reintroduction of wild-type Hdac3 restored normal B-cell development, whereas an Hdac3 point mutant lacking deacetylase activity failed to complement this defect. Thus, the deacetylase activity of Hdac3 is required for the generation of mature B cells.

Chromatin Domain Organization of the TCRb Locus and Its Perturbation by Ectopic CTCF Binding

CTCF-mediated chromatin interactions influence organization and function of mammalian genome in diverse ways. We analyzed the interactions among CTCF binding sites (CBS) at the murine TCRb locus to discern the role of CTCF-mediated interactions in the regulation of transcription and VDJ recombination. Chromosome conformation capture analysis revealed thymocyte-specific long-range intrachromosomal interactions among various CBS across the locus that were relevant for defining the limit of the enhancer Eb-regulated recombination center (RC) and for facilitating the spatial proximity of TCRb variable (V) gene segments to the RC. Ectopic CTCF binding in the RC region, effected via genetic manipulation, altered CBS-directed chromatin loops, interfered with RC establishment, and reduced the spatial proximity of the RC with Trbv segments. Changes in chromatin loop organization by ectopic CTCF binding were relatively modest but influenced transcription and VDJ recombination dramatically. Besides revealing the importance of CTCF-mediated chromatin organization for TCRb regulation, the observed chromatin loops were consistent with the emerging idea that CBS orientations influence chromatin loop organization and underscored the importance of CBS orientations for defining chromatin architecture that supports VDJ recombination. Further, our study suggests that in addition to mediating long-range chromatin interactions, CTCF influences intricate configuration of chromatin loops that govern functional interactions between elements.

Insights into immune system development and function from mouse T-cell repertoires

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

TRIg: a robust alignment pipeline for non-regular T-cell receptor and immunoglobulin sequences

BACKGROUND

T cells and B cells are essential in the adaptive immunity via expressing T cell receptors and immunoglogulins respectively for recognizing antigens. To recognize a wide variety of antigens, a highly diverse repertoire of receptors is generated via complex recombination of the receptor genes. Reasonably, frequencies of the recombination events have been shown to predict immune diseases and provide insights into the development of immunity. The field is further boosted by high-throughput sequencing and several computational tools have been released to analyze the recombined sequences. However, all current tools assume regular recombination of the receptor genes, which is not always valid in data prepared using a RACE approach. Compared to the traditional multiplex PCR approach, RACE is free of primer bias, therefore can provide accurate estimation of recombination frequencies. To handle the non-regular recombination events, a new computational program is needed.

RESULTS

We propose TRIg to handle non-regular T cell receptor and immunoglobulin sequences. Unlike all current programs, TRIg does alignments to the whole receptor gene instead of only to the coding regions. This brings new computational challenges, e.g., ambiguous alignments due to multiple hits to repetitive regions. To reduce ambiguity, TRIg applies a heuristic strategy and incorporates gene annotation to identify authentic alignments. On our own and public RACE datasets, TRIg correctly identified non-regularly recombined sequences, which could not be achieved by current programs. TRIg also works well for regularly recombined sequences.

CONCLUSIONS

TRIg takes into account non-regular recombination of T cell receptor and immunoglobulin genes, therefore is suitable for analyzing RACE data. Such analysis will provide accurate estimation of recombination events, which will benefit various immune studies directly. In addition, TRIg is suitable for studying aberrant recombination in immune diseases. TRIg is freely available at https://github.com/TLlab/trig .

Contribution of secondary Igkappa rearrangement to primary immunoglobulin repertoire diversification

Abs reactive to DNA and DNA/histone complexes are a distinguished characteristic of primary immunoglobulin repertoires in autoimmune B6.MRL-Fas^lpr and MRL/MpJ-Fas^lpr mice. These mice are defective in Fas receptor, which is critical for the apoptosis of autoreactive B cells by an extrinsic pathway. In the present study, we explored the possibility that bone marrow small pre-B and immature B cells from adult B6.MRL-Fas^lpr mice and MRL/MpJ-Fas^lpr mice respectively, which contain autoreactive B-cell antigen receptors (BCR) and manifest autoimmune syndromes, exhibit enhanced receptor editing patterns. Indeed, FAS^lpr pre B and immature B cells were shown to possess more ongoing replacements of non-productive (nP) than productive (P) primary VκJκ rearrangements. Significantly, the P vs nP ratios of these replaced primary rearrangements were 1:2, thus indicating that κ light-chain production appears not to inhibit secondary rearrangements. In addition, we identified multiple atypical rearrangements, such as Vκ cRS (cryptic recombination signals) cleavages. These results suggest that the onset of light chain secondary rearrangements persists similarly as a non-selected mode and independent of BCR autoreactivity during certain developmental windows of bone marrow B cells in lupus-prone mice and control, and leads us to propose the function of secondary, de novo Igκ rearrangements to increase BCR diversity.

Inferring processes underlying B-cell repertoire diversity

We quantify the VDJ recombination and somatic hypermutation processes in human B cells using probabilistic inference methods on high-throughput DNA sequence repertoires of human B-cell receptor heavy chains. Our analysis captures the statistical properties of the naive repertoire, first after its initial generation via VDJ recombination and then after selection for functionality. We also infer statistical properties of the somatic hypermutation machinery (exclusive of subsequent effects of selection). Our main results are the following: the B-cell repertoire is substantially more diverse than T-cell repertoires, owing to longer junctional insertions; sequences that pass initial selection are distinguished by having a higher probability of being generated in a VDJ recombination event; somatic hypermutations have a non-uniform distribution along the V gene that is well explained by an independent site model for the sequence context around the hypermutation site.

In-Depth Assessment of Within-Individual and Inter-Individual Variation in the B Cell Receptor Repertoire

High-throughput sequencing of the B cell receptor (BCR) repertoire can provide rapid characterization of the B cell response in a wide variety of applications in health, after vaccination and in infectious, inflammatory and immune-driven disease, and is starting to yield clinical applications. However, the interpretation of repertoire data is compromised by a lack of studies to assess the intra and inter-individual variation in the BCR repertoire over time in healthy individuals. We applied a standardized isotype-specific BCR repertoire deep sequencing protocol to a single highly sampled participant, and then evaluated the method in 9 further participants to comprehensively describe such variation. We assessed total repertoire metrics of mutation, diversity, VJ gene usage and isotype subclass usage as well as tracking specific BCR sequence clusters. There was good assay reproducibility (both in PCR amplification and biological replicates), but we detected striking fluctuations in the repertoire over time that we hypothesize may be due to subclinical immune activation. Repertoire properties were unique for each individual, which could partly be explained by a decrease in IgG2 with age, and genetic differences at the immunoglobulin locus. There was a small repertoire of public clusters (0.5, 0.3, and 1.4% of total IgA, IgG, and IgM clusters, respectively), which was enriched for expanded clusters containing sequences with suspected specificity toward antigens that should have been historically encountered by all participants through prior immunization or infection. We thus provide baseline BCR repertoire information that can be used to inform future study design, and aid in interpretation of results from these studies. Furthermore, our results indicate that BCR repertoire studies could be used to track changes in the public repertoire in and between populations that might relate to population immunity against infectious diseases, and identify the characteristics of inflammatory and immunological diseases.

PRKDC mutations associated with immunodeficiency, granuloma, and autoimmune regulator-dependent autoimmunity

BACKGROUND

PRKDC encodes for DNA-dependent protein kinase catalytic subunit (DNA-PKcs), a kinase that forms part of a complex (DNA-dependent protein kinase [DNA-PK]) crucial for DNA double-strand break repair and V(D)J recombination. In mice DNA-PK also interacts with the transcription factor autoimmune regulator (AIRE) to promote central T-cell tolerance.

OBJECTIVE

We sought to understand the causes of an inflammatory disease with granuloma and autoimmunity associated with decreasing T- and B-cell counts over time that had been diagnosed in 2 unrelated patients.

METHODS

Genetic, molecular, and functional analyses were performed to characterize an inflammatory disease evocative of a combined immunodeficiency.

RESULTS

We identified PRKDC mutations in both patients. These patients exhibited a defect in DNA double-strand break repair and V(D)J recombination. Whole-blood mRNA analysis revealed a strong interferon signature. On activation, memory T cells displayed a skewed cytokine response typical of TH2 and TH1 but not TH17. Moreover, mutated DNA-PKcs did not promote AIRE-dependent transcription of peripheral tissue antigens in vitro. The latter defect correlated in vivo with production of anti-calcium-sensing receptor autoantibodies, which are typically found in AIRE-deficient patients. In addition, 9 months after bone marrow transplantation, patient 1 had Hashimoto thyroiditis, suggesting that organ-specific autoimmunity might be linked to nonhematopoietic cells, such as AIRE-expressing thymic epithelial cells.

CONCLUSION

Deficiency of DNA-PKcs, a key AIRE partner, can present as an inflammatory disease with organ-specific autoimmunity, suggesting a role for DNA-PKcs in regulating autoimmune responses and maintaining AIRE-dependent tolerance in human subjects.

Clonal and constricted T cell repertoire in Common Variable Immune Deficiency

We used high throughput sequencing to examine the structure and composition of the T cell receptor β chain in Common Variable Immune Deficiency (CVID). TCRβ CDR3 regions were amplified and sequenced from genomic DNA of 44 adult CVID subjects and 22 healthy adults, using a high-throughput multiplex PCR. CVID TCRs had significantly less junctional diversity, fewer n-nucleotide insertions and deletions, and completely lacked a population of highly modified TCRs, with 13 or more V-gene nucleotide deletions, seen in healthy controls. The CVID CDR3 sequences were significantly more clonal than control DNA, and displayed unique V gene usage. Despite reduced junctional diversity, increased clonality and similar infectious exposures, DNA of CVID subjects shared fewer TCR sequences as compared to controls. These abnormalities are pervasive, found in out-of-frame sequences and thus independent of selection and were not associated with specific clinical complications. These data support an inherent T cell defect in CVID.

DNA bending in the synaptic complex in V(D)J recombination: turning an ancestral transpososome upside down

In all jawed vertebrates RAG (recombination activating gene) recombinase orchestrates V(D)J recombination in B and T lymphocyte precursors, assembling the V, D and J germline gene segments into continuous functional entities which encode the variable regions of their immune receptors. V(D)J recombination is the process by which most of the diversity of our specific immune receptors is acquired and is thought to have originated by domestication of a transposon in the genome of a vertebrate.  RAG acts similarly to the cut and paste transposases, by first binding two recombination signal DNA sequences (RSSs), which flank the two coding genes to be adjoined, in a process called synaptic or paired complex (PC) formation. At these RSS-coding borders, RAG first nicks one DNA strand, then creates hairpins, thus cleaving the duplex DNA at both RSSs. Although RAG reaction mechanism resembles that of insect mobile element transposases and RAG itself can inefficiently perform intramolecular and intermolecular integration into the target DNA, inside the nuclei of the developing lymphocytes transposition is extremely rare and is kept under proper surveillance. Our review may help understand how RAG synaptic complex organization prevents deleterious transposition. The phosphoryl transfer reaction mechanism of RNAseH-like fold DDE motif enzymes, including RAG, is discussed accentuating the peculiarities described for various transposases from the light of their available high resolution structures (Tn5, Mu, Mos1 and Hermes). Contrasting the structural 3D organization of DNA in these transpososomes with that of the RSSs-DNA in RAG PC allows us to propose several clues for how evolutionarily RAG may have become “specialized” in recombination versus transposition.

Qualifying high-throughput immune repertoire sequencing

Diversity of B and T cell receptors, achieved by gene recombination and somatic hypermutation, allows the immune system for recognition and targeted reaction against various threats. Next-generation sequencing for assessment of a cell's gene composition and variation makes deep analysis of one individual's immune spectrum feasible. An easy to apply but detailed analysis and visualization strategy is necessary to process all sequences generated. We performed sequencing utilizing the 454 system for CLL and control samples, utilized the IMGT database and applied the presented analysis tools. With the applied protocol, malignant clones are found and characterized, mutational status compared to germline identity is elaborated in detail showing that the CLL mutation status is not as monoclonal as generally thought. On the other hand, this strategy is not solely applicable to the 454 sequencing system but can easily be transferred to any other next-generation sequencing platform.

Erroneous class switching and false VDJ recombination: molecular dissection of t(8;14)/MYC-IGH translocations in Burkitt-type lymphoblastic leukemia/B-cell lymphoma

The chromosomal translocation t(8;14)(q24;q32) with juxtaposition of MYC to enhancer elements in the immunoglobulin heavy chain (IGH) gene locus is the genetic hallmark of the majority of Burkitt lymphoma and a subset of Diffuse large B-cell lymphoma patients. Around 3% of adult B-lineage acute lymphoblastic leukemia (ALL) patients show this aberration. Flow cytometry mostly reveals a "mature B-ALL" or "Burkitt-type" ALL immunophenotype. Using long-distance PCR for t(8;14)/MYC-IGH fusion, we investigated bone marrow, peripheral blood and a few other samples with suspected Burkitt-ALL or mature B-ALL and identified 133 MYC-IGH-positive cases. The location of the chromosomal breaks in the IGH joining and the 8 different switch regions was determined using a set of long-distance PCRs. The chromosomal breakpoints with the adjacent MYC regions on 8q24 were characterized by direct sequencing in 49 cases. The distribution of chromosomal breaks among the IGH joining and switch regions was the following: JH 23.3%, M 21.8%, G1 15.0%, G2 7.5%, G3 3.8%, G4 4.5%, A1 12.8%, A2 3.8%, E 7.5%. Two breakpoint clusters near MYC were delineated. There was no clear correlation between the degree of somatic hypermutation and the chromosomal break locations. Epstein Barr virus was detected in 5 cases (4%). This detailed and extensive molecular analysis illustrates the molecular complexity of the MYC-IGH translocations and the detected distribution of breakpoints provides additional evidence that this translocation results from failed switch and VDJ recombinations. This study may serve as a model for the analysis of other IGH translocations in B-cell lymphoma.