Limitations of lymphoblastoid cell lines for establishing genetic reference datasets in the immunoglobulin loci

Lymphoblastoid cell lines (LCLs) have been critical to establishing genetic resources for biomedical science. They have been used extensively to study human genetic diversity, genome function, and inform the development of tools and methodologies for augmenting disease genetics research. While the validity of variant callsets from LCLs has been demonstrated for most of the genome, previous work has shown that DNA extracted from LCLs is modified by V(D)J recombination within the immunoglobulin (IG) loci, regions that harbor antibody genes critical to immune system function. However, the impacts of V(D)J on short read sequencing data generated from LCLs has not been extensively investigated. In this study, we used LCL-derived short read sequencing data from the 1000 Genomes Project (n = 2,504) to identify signatures of V(D)J recombination. Our analyses revealed sample-level impacts of V(D)J recombination that varied depending on the degree of inferred monoclonality. We showed that V(D)J associated somatic deletions impacted genotyping accuracy, leading to adulterated population-level estimates of allele frequency and linkage disequilibrium. These findings illuminate limitations of using LCLs and short read data for building genetic resources in the IG loci, with implications for interpreting previous disease association studies in these regions.

Genomic analysis of a second rainbow trout line (Arlee) leads to an extended description of the IGH VDJ gene repertoire

High-throughput sequencing technologies brought a renewed interest for immune repertoires. Fish Ab and B cell repertoires are no exception, and their comprehensive analysis can both provide new insights into poorly understood immune mechanisms, and identify markers of protection after vaccination. However, the lack of genomic description and standardized nomenclature of IG genes hampers accurate annotation of Ig mRNA deep sequencing data. Complete genome sequences of Atlantic salmon and rainbow trout (Swanson line) recently allowed us to establish a comprehensive and coherent annotation of Salmonid IGH genes following IMGT standards. Here we analyzed the IGHV, D, and J genes from the newly released genome of a second rainbow trout line (Arlee). We confirmed the validity of salmonid IGHV subgroups, and extended the description of the rainbow trout IGH gene repertoire with novel sequences, while keeping nomenclature continuity. This work provides an important resource for annotation of high-throughput Ab repertoire sequencing data.

Templated insertions at VD and DJ junctions create unique B-cell receptors in the healthy B-cell repertoire

RAG complexes recognise (cryptic) RSS sites both in and outside immunoglobulin sites. Excision circles may be reinserted into V(D)J rearrangements as long templated insertions to diversify the adaptive immune repertoire. We show that such VDJ with templated insertions are incidentally found in the repertoire of healthy donors.

Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination

Despite the extreme diversity of T-cell repertoires, many identical T-cell receptor (TCR) sequences are found in a large number of individual mice and humans. These widely shared sequences, often referred to as "public," have been suggested to be over-represented due to their potential immune functionality or their ease of generation by V(D)J recombination. Here, we show that even for large cohorts, the observed degree of sharing of TCR sequences between individuals is well predicted by a model accounting for the known quantitative statistical biases in the generation process, together with a simple model of thymic selection. Whether a sequence is shared by many individuals is predicted to depend on the number of queried individuals and the sampling depth, as well as on the sequence itself, in agreement with the data. We introduce the degree of publicness conditional on the queried cohort size and the size of the sampled repertoires. Based on these observations, we propose a public/private sequence classifier, "PUBLIC" (Public Universal Binary Likelihood Inference Classifier), based on the generation probability, which performs very well even for small cohort sizes.

The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given datasets. This procedure is well developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper, we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency but also discrete objects such as B-cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective, the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B-cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.

The Complex Interplay between DNA Injury and Repair in Enzymatically Induced Mutagenesis and DNA Damage in B Lymphocytes

Lymphocytes are endowed with unique and specialized enzymatic mutagenic properties that allow them to diversify their antigen receptors, which are crucial sensors for pathogens and mediators of adaptive immunity. During lymphocyte development, the antigen receptors expressed by B and T lymphocytes are assembled in an antigen-independent fashion by ordered variable gene segment recombinations (V(D)J recombination), which is a highly ordered and regulated process that requires the recombination activating gene products 1 & 2 (RAG1, RAG2). Upon activation by antigen, B lymphocytes undergo additional diversifications of their immunoglobulin B-cell receptors. Enzymatically induced somatic hypermutation (SHM) and immunoglobulin class switch recombination (CSR) improves the affinity for antigen and shape the effector function of the humoral immune response, respectively. The activation-induced cytidine deaminase (AID) enzyme is crucial for both SHM and CSR. These processes have evolved to both utilize as well as evade different DNA repair and DNA damage response pathways. The delicate balance between enzymatic mutagenesis and DNA repair is crucial for effective immune responses and the maintenance of genomic integrity. Not surprisingly, disturbances in this balance are at the basis of lymphoid malignancies by provoking the formation of oncogenic mutations and chromosomal aberrations. In this review, we discuss recent mechanistic insight into the regulation of RAG1/2 and AID expression and activity in lymphocytes and the complex interplay between these mutagenic enzymes and DNA repair and DNA damage response pathways, focusing on the base excision repair and mismatch repair pathways. We discuss how disturbances of this interplay induce genomic instability and contribute to oncogenesis.

Amplification and selection of PRRSV-activated VDJ repertoires in pigs secreting distinct neutralizing antiboidies

Neutralizing antibodies (nAbs) play an important role in protective immunity against porcine reproductive and respiratory syndrome virus (PRRSV) infection. However, the characterization of PRRSV nAb repertoires is rarely investigated. In this study, we developed a swine VDJ amplification method and selection criteria for the characterization of PRRSV-activated VDJ repertoires. According to clonal expansion theory, two separated aliquots of lymph nodes from pigs producing different PRRSV nAbs were utilized to determine the activated B-cell repertoires. Swine VDJ repertoires from a mock-infected pig and PRRSV-infected pigs secreting no detectable nAbs, only homologous nAbs, and broad nAbs were amplified by a single pair of primers that could detect all seven major VDJ genes. The amplicons were cloned and sequenced to generate 385 VDJ sequences. Sequence alignment showed that the diversification of VDJ genes was mainly due to the variation in complementarity determining regions (CDRs), especially CDR3. Based on selection criteria, shared and abundant sequences were identified in two separated aliquots from PRRSV-infected pigs but not from the mock-infected pig, suggesting they were secreted from PRRSV-activated B cells. Thus, the amplification and selection method provide a potential alternative for the characterization of swine VDJ repertoires. However, additional experiments are required to determine whether the shared and abundant VDJ lineages identified in this study are PRRSV-specific or distinct neutralizing-antibodies-associated.

High-resolution analysis of the B cell repertoire before and after polyethylene glycol fusion reveals preferential fusion of rare antigen-specific B cells

OBJECTIVE

The hybridoma technology is one of the most important advances in clinical immunology. Little is known about the differences between the antibodies produced during the in vivo immune response and those recovered in hybridoma libraries. Here, we investigate a potential fusion bias inherent to the hybridoma production process.

METHODS

Transgenic rats carrying human Ig heavy and light chain loci were immunized with measles virus (MV) to generate human mAbs. Usin g high-throughput sequencing of IgH mRNA, we compared the IgH repertoire of lymph nodes and the derived hybridoma library using the sequences of the MV-specific hybridoma clones as a reference set with known specificity.

RESULTS

We observed that large clonotypes from the lymph nodes were not represented in the hybridoma library, but low-frequency B cell populations became highly enriched and most hybridoma clones were derived from these. Our data also showed that identical CDR3s evolved from diverse VDJ recombinations, indicating convergence of different B cells subpopulations towards expression of antibodies with similar paratopes.

CONCLUSION

The efficient generation of mAbs results from a fusion process highly selective for rare antigen-specific B cells rather than in vivo expanded populations. Antibodies of particular interest may therefore be missed during classical hybridoma production.

De novo DNA Methyltransferases Dnmt3a and Dnmt3b regulate the onset of Igκ light chain rearrangement during early B-cell development

Immunoglobulin genes V(D)J rearrangement during early lymphopoiesis is a critical process involving sequential recombination of the heavy and light chain loci. A number of transcription factors act together with temporally activated recombinases and chromatin accessibility changes to regulate this complex process. Here, we deleted the de novo DNA methyltransferases Dnmt3a and Dnmt3b in early B cells of conditionally targeted mice, and monitored the process of V(D)J recombination. Dnmt3a and Dnmt3b deletion resulted in precocious recombination of the immunoglobulin κ light chain without impairing the differentiation of mature B cells or overall B-cell development. Ex vivo culture of IL-7 restricted early B-cell progenitors lacking Dnmt3a and Dnmt3b showed precocious Vκ-Jκ rearrangements that are limited to the proximal Vκ genes. Furthermore, B-cell progenitors deficient in Dnmt3a and Dnmt3b showed elevated levels of germline transcripts at the proximal Vκ genes, alterations in methylation patterns at Igκ enhancer sites and increased expression of the transcription factor E2A. Our data suggest that Dnmt3a and Dnmt3b are critical to regulate the onset of Igκ light chain rearrangement during early B-cell development.

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.

Violation of an evolutionarily conserved immunoglobulin diversity gene sequence preference promotes production of dsDNA-specific IgG antibodies

Variability in the developing antibody repertoire is focused on the third complementarity determining region of the H chain (CDR-H3), which lies at the center of the antigen binding site where it often plays a decisive role in antigen binding. The power of VDJ recombination and N nucleotide addition has led to the common conception that the sequence of CDR-H3 is unrestricted in its variability and random in its composition. Under this view, the immune response is solely controlled by somatic positive and negative clonal selection mechanisms that act on individual B cells to promote production of protective antibodies and prevent the production of self-reactive antibodies. This concept of a repertoire of random antigen binding sites is inconsistent with the observation that diversity (D_H) gene segment sequence content by reading frame (RF) is evolutionarily conserved, creating biases in the prevalence and distribution of individual amino acids in CDR-H3. For example, arginine, which is often found in the CDR-H3 of dsDNA binding autoantibodies, is under-represented in the commonly used D_H RFs rearranged by deletion, but is a frequent component of rarely used inverted RF1 (iRF1), which is rearranged by inversion. To determine the effect of altering this germline bias in D_H gene segment sequence on autoantibody production, we generated mice that by genetic manipulation are forced to utilize an iRF1 sequence encoding two arginines. Over a one year period we collected serial serum samples from these unimmunized, specific pathogen-free mice and found that more than one-fifth of them contained elevated levels of dsDNA-binding IgG, but not IgM; whereas mice with a wild type D_H sequence did not. Thus, germline bias against the use of arginine enriched D_H sequence helps to reduce the likelihood of producing self-reactive antibodies.

Independent recruitment of Igh alleles in V(D)J recombination

How the vast majority of B cells express only one of the two alleles at their immunoglobulin loci remains a biological puzzle. Here, in mice reconstituted with a single haematopoietic stem cell, we demonstrate that each of the two immunoglobulin heavy chain (Igh) alleles has a similar probability to be the first to undergo V(H) to DJ(H) rearrangement. We also observe this similar probability in clones from multipotent and common lymphoid precursors. The extreme biases in the expression of the alleles that we find in more differentiated subsets are mostly due to constraints imposed by early rearrangements. Our data demonstrate that each of the two Igh alleles in a B cell behaves independently of the other, up to the moment when a successful rearrangement in one allele triggers a feedback mechanism that prevents further recombination.

PP4 is essential for germinal center formation and class switch recombination in mice

PP4 is a serine/threonine phosphatase required for immunoglobulin (Ig) VDJ recombination and pro-B/pre-B cell development in mice. To elucidate the role of PP4 in mature B cells, we ablated the catalytic subunit of murine PP4 invivo utilizing the CD23 promoter and cre-loxP recombination and generated CD23^crePP4^F/F mice. The development of follicular and marginal zone B cells was unaffected in these mutants, but the proliferation of mature PP4-deficient B cells stimulated by invitro treatment with either anti-IgM antibody (Ab) or LPS was partially impaired. Interestingly, the induction of CD80 and CD86 expression on these stimulated B cells was normal. Basal levels of serum Igs of all isotypes were strongly reduced in CD23^crePP4^F/F mice, and their B cells showed a reduced efficiency of class switch recombination (CSR) invitro upon stimulation by LPS or LPS plus IL-4. When CD23^crePP4^F/F mice were challenged with either the T cell-dependent antigen TNP-KLH or the T cell-independent antigen TNP-Ficoll, or by H1N1 virus infection, the mutant animals failed to form germinal centers (GCs) in the spleen and the draining mediastinal lymph nodes, and did not efficiently mount antigen-specific humoral responses. In the resting state, PP4-deficient B cells exhibited pre-existing DNA fragmentation. Upon stimulation by DNA-damaging drug etoposide invitro, mutant B cells showed increased cleavage of caspase 3. In addition, the mutant B cells displayed impaired CD40-mediated MAPK activation, abnormal IgM-mediated NF-κB activation, and reduced S phase entry upon IgM/CD40-stimulation. Taken together, our results establish a novel role for PP4 in CSR, and reveal crucial functions for PP4 in the maintenance of genomic stability, GC formation, and B cell-mediated immune responses.

Histone methylation and V(D)J recombination

V(D)J recombination is the process by which the diversity of antigen receptor genes is generated and is also indispensable for lymphocyte development. This recombination event occurs in a cell lineage- and stage-specific manner, and is carefully controlled by chromatin structure and ordered histone modifications. The recombinationally active V(D)J loci are associated with hypermethylation at lysine4 of histone H3 and hyperacetylation of histones H3/H4. The recombination activating gene 1 (RAG1) and RAG2 complex initiates recombination by introducing double-strand DNA breaks at recombination signal sequences (RSS) adjacent to each coding sequence. To be recognized by the RAG complex, RSS sites must be within an open chromatin context. In addition, the RAG complex specifically recognizes hypermethylated H3K4 through its plant homeodomain (PHD) finger in the RAG2 C terminus, which stimulates RAG catalytic activity via that interaction. In this review, we describe how histone methylation controls V(D)J recombination and discuss its potential role in lymphoid malignancy by mistargeting the RAG complex.

Regulation of immunoglobulin class-switch recombination: choreography of noncoding transcription, targeted DNA deamination, and long-range DNA repair

Upon encountering antigens, mature IgM-positive B lymphocytes undergo class-switch recombination (CSR) wherein exons encoding the default Cμ constant coding gene segment of the immunoglobulin (Ig) heavy-chain (Igh) locus are excised and replaced with a new constant gene segment (referred to as "Ch genes", e.g., Cγ, Cɛ, or Cα). The B cell thereby changes from expressing IgM to one producing IgG, IgE, or IgA, with each antibody isotype having a different effector function during an immune reaction. CSR is a DNA deletional-recombination reaction that proceeds through the generation of DNA double-strand breaks (DSBs) in repetitive switch (S) sequences preceding each Ch gene and is completed by end-joining between donor Sμ and acceptor S regions. CSR is a multistep reaction requiring transcription through S regions, the DNA cytidine deaminase AID, and the participation of several general DNA repair pathways including base excision repair, mismatch repair, and classical nonhomologous end-joining. In this review, we discuss our current understanding of how transcription through S regions generates substrates for AID-mediated deamination and how AID participates not only in the initiation of CSR but also in the conversion of deaminated residues into DSBs. Additionally, we review the multiple processes that regulate AID expression and facilitate its recruitment specifically to the Ig loci, and how deregulation of AID specificity leads to oncogenic translocations. Finally, we summarize recent data on the potential role of AID in the maintenance of the pluripotent stem cell state during epigenetic reprogramming.

Onset of immune senescence defined by unbiased pyrosequencing of human immunoglobulin mRNA repertoires

The immune system protects us from foreign substances or pathogens by generating specific antibodies. The variety of immunoglobulin (Ig) paratopes for antigen recognition is a result of the V(D)J rearrangement mechanism, while a fast and efficient immune response is mediated by specific immunoglobulin isotypes obtained through class switch recombination (CSR). To get a better understanding on how antibody-based immune protection works and how it changes with age, the interdependency between these two parameters need to be addressed. Here, we have performed an in depth analysis of antibody repertoires of 14 healthy donors representing different gender and age groups. For this task, we developed a unique pyrosequencing approach, which is able to monitor the expression levels of all immunoglobulin V(D)J recombinations of all isotypes including subtypes in an unbiased and quantitative manner. Our results show that donors have individual immunoglobulin repertoires and cannot be clustered according to V(D)J recombination patterns, neither by age nor gender. However, after incorporating isotype-specific analysis and considering CSR information into hierarchical clustering the situation changes. For the first time the donors cluster according to age and separate into young adults and elderly donors (>50). As a direct consequence, this clustering defines the onset of immune senescence at the age of fifty and beyond. The observed age-dependent reduction of CSR ability proposes a feasible explanation why reduced efficacy of vaccination is seen in the elderly and implies that novel vaccine strategies for the elderly should include the "Golden Agers".

Molecular characterization of immunoglobulin gene rearrangements in diffuse large B-cell lymphoma: antigen-driven origin and IGHV4-34 as a particular subgroup of the non-GCB subtype

The pathogenesis of diffuse large B-cell lymphoma (DLBCL) remains partially unknown. The analysis of the B-cell receptor of the malignant cells could contribute to a better understanding of the DLBCL biology. We studied the molecular features of the immunoglobulin heavy chain (IGH) rearrangements in 165 patients diagnosed with DLBCL not otherwise specified. Clonal IGH rearrangements were amplified according to the BIOMED-2 protocol and PCR products were sequenced directly. We also analyzed the criteria for stereotyped patterns in all complete IGHV-IGHD-IGHJ (V-D-J) sequences. Complete V-D-J rearrangements were identified in 130 of 165 patients. Most cases (89%) were highly mutated, but 12 sequences were truly unmutated or minimally mutated. Three genes, IGHV4-34, IGHV3-23, and IGHV4-39, accounted for one third of the whole cohort, including an overrepresentation of IGHV4-34 (15.5% overall). Interestingly, all IGHV4-34 rearrangements and all unmutated sequences belonged to the nongerminal center B-cell-like (non-GCB) subtype. Overall, we found three cases following the current criteria for stereotyped heavy chain VH CDR3 sequences, two of them belonging to subsets previously described in CLL. IGHV gene repertoire is remarkably biased, implying an antigen-driven origin in DLBCL. The particular features in the sequence of the immunoglobulins suggest the existence of particular subgroups within the non-GCB subtype.

A TNF-regulated recombinatorial macrophage immune receptor implicated in granuloma formation in tuberculosis

Macrophages play a central role in host defense against mycobacterial infection and anti- TNF therapy is associated with granuloma disorganization and reactivation of tuberculosis in humans. Here, we provide evidence for the presence of a T cell receptor (TCR) αβ based recombinatorial immune receptor in subpopulations of human and mouse monocytes and macrophages. In vitro, we find that the macrophage-TCRαβ induces the release of CCL2 and modulates phagocytosis. TNF blockade suppresses macrophage-TCRαβ expression. Infection of macrophages from healthy individuals with mycobacteria triggers formation of clusters that express restricted TCR Vβ repertoires. In vivo, TCRαβ bearing macrophages abundantly accumulate at the inner host-pathogen contact zone of caseous granulomas from patients with lung tuberculosis. In chimeric mouse models, deletion of the variable macrophage-TCRαβ or TNF is associated with structurally compromised granulomas of pulmonary tuberculosis even in the presence of intact T cells. These results uncover a TNF-regulated recombinatorial immune receptor in monocytes/macrophages and demonstrate its implication in granuloma formation in tuberculosis.