The complete nucleotide sequence of the human immunoglobulin heavy chain variable region locus

AIRR-C IG Reference Sets: curated sets of immunoglobulin heavy and light chain germline genes

Introduction

Analysis of an individual's immunoglobulin (IG) gene repertoire requires the use of high-quality germline gene reference sets. When sets only contain alleles supported by strong evidence, AIRR sequencing (AIRR-seq) data analysis is more accurate and studies of the evolution of IG genes, their allelic variants and the expressed immune repertoire is therefore facilitated.

Methods

The Adaptive Immune Receptor Repertoire Community (AIRR-C) IG Reference Sets have been developed by including only human IG heavy and light chain alleles that have been confirmed by evidence from multiple high-quality sources. To further improve AIRR-seq analysis, some alleles have been extended to deal with short 3' or 5' truncations that can lead them to be overlooked by alignment utilities. To avoid other challenges for analysis programs, exact paralogs (e.g. IGHV1-69*01 and IGHV1-69D*01) are only represented once in each set, though alternative sequence names are noted in accompanying metadata.

Results and discussion

The Reference Sets include less than half the previously recognised IG alleles (e.g. just 198 IGHV sequences), and also include a number of novel alleles: 8 IGHV alleles, 2 IGKV alleles and 5 IGLV alleles. Despite their smaller sizes, erroneous calls were eliminated, and excellent coverage was achieved when a set of repertoires comprising over 4 million V(D)J rearrangements from 99 individuals were analyzed using the Sets. The version-tracked AIRR-C IG Reference Sets are freely available at the OGRDB website (https://ogrdb.airr-community.org/germline_sets/Human) and will be regularly updated to include newly observed and previously reported sequences that can be confirmed by new high-quality data.

Updates on lymphoblastic leukemia/lymphoma classification and minimal/measurable residual disease analysis

Lymphoblastic leukemia/lymphoma (ALL/LBL), especially certain subtypes, continues to confer morbidity and mortality despite significant therapeutic advances. The pathologic classification of ALL/LBL, especially that of B-ALL, has recently substantially expanded with the identification of several distinct and prognostically important genetic drivers. These discoveries are reflected in both current classification systems, the World Health Organization (WHO) 5th edition and the new International Consensus Classification (ICC). In this article, novel subtypes of B-ALL are reviewed, including DUX4, MEF2D and ZNF384-rearranged B-ALL; the rare pediatric entity B-ALL with TLF3::HLF, now added to the classifications, is discussed; updates to the category of B-ALL with BCR::ABL1-like features (Ph-like B-ALL) are summarized; and emerging genetic subtypes of T-ALL are presented. The second half of the article details current approaches to minimal/measurable residual disease (MRD) detection in B-ALL and T-ALL and presents anticipated challenges to current approaches in the burgeoning era of antigen-directed immunotherapy.

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

Untangling immunoglobulin genotype-function associations

Immunoglobulin (IG) genes, encoding B cell receptors (BCRs), are fundamental components of the mammalian immune system, which evolved to recognize the diverse antigenic universe present in nature. To handle these myriad inputs, BCRs are generated through combinatorial recombination of a set of highly polymorphic germline genes, resulting in a vast repertoire of antigen receptors that initiate responses to pathogens and regulate commensals. Following antigen recognition and B cell activation, memory B cells and plasma cells form, allowing for the development of anamnestic antibody (Ab) responses. How inherited variation in IG genes impacts host traits, disease susceptibility, and Ab recall responses is a topic of great interest. Here, we consider approaches to translate emerging knowledge about IG genetic diversity and expressed repertoires to inform our understanding of Ab function in health and disease etiology. As our understanding of IG genetics grows, so will our need for tools to decipher preferences for IG gene or allele usage in different contexts, to better understand antibody responses at the population level.

Frequent use of IGHV3-30-3 in SARS-CoV-2 neutralizing antibody responses

The antibody response to SARS-CoV-2 shows biased immunoglobulin heavy chain variable (IGHV) gene usage, allowing definition of genetic signatures for some classes of neutralizing antibodies. We investigated IGHV gene usage frequencies by sorting spike-specific single memory B cells from individuals infected with SARS-CoV-2 early in the pandemic. From two study participants and 703 spike-specific B cells, the most used genes were IGHV1-69, IGHV3-30-3, and IGHV3-30. Here, we focused on the IGHV3-30 group of genes and an IGHV3-30-3-using ultrapotent neutralizing monoclonal antibody, CAB-F52, which displayed broad neutralizing activity also in its germline-reverted form. IGHV3-30-3 is encoded by a region of the IGH locus that is highly variable at both the allelic and structural levels. Using personalized IG genotyping, we found that 4 of 14 study participants lacked the IGHV3-30-3 gene on both chromosomes, raising the question if other, highly similar IGHV genes could substitute for IGHV3-30-3 in persons lacking this gene. In the context of CAB-F52, we found that none of the tested IGHV3-33 alleles, but several IGHV3-30 alleles could substitute for IGHV3-30-3, suggesting functional redundancy between the highly homologous IGHV3-30 and IGHV3-30-3 genes for this antibody.

Immunoglobulin germline gene polymorphisms influence the function of SARS-CoV-2 neutralizing antibodies

The human immunoglobulin heavy-chain (IGH) locus is exceptionally polymorphic, with high levels of allelic and structural variation. Thus, germline IGH genotypes are personal, which may influence responses to infection and vaccination. For an improved understanding of inter-individual differences in antibody responses, we isolated SARS-CoV-2 spike-specific monoclonal antibodies from convalescent health care workers, focusing on the IGHV1-69 gene, which has the highest level of allelic variation of all IGHV genes. The IGHV1-69^∗20-using CAB-I47 antibody and two similar antibodies isolated from an independent donor were critically dependent on allele usage. Neutralization was retained when reverting the V region to the germline IGHV1-69^∗20 allele but lost when reverting to other IGHV1-69 alleles. Structural data confirmed that two germline-encoded polymorphisms, R50 and F55, in the IGHV1-69 gene were required for high-affinity receptor-binding domain interaction. These results demonstrate that polymorphisms in IGH genes can influence the function of SARS-CoV-2 neutralizing antibodies.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Adaptive immune receptor genotyping using the corecount program

We present a new Rep-Seq analysis tool called corecount, for analyzing genotypic variation in immunoglobulin (IG) and T cell receptor (TCR) genes. corecount is highly efficient at identifying V alleles, including those that are infrequently used in expressed repertoires and those that contain 3' end variation that are otherwise refractory to reliable identification during germline inference from expressed libraries. Furthermore, corecount facilitates accurate D and J gene genotyping. The output is highly reproducible and facilitates the comparison of genotypes from multiple individuals, such as those from clinical cohorts. Here, we applied corecount to the genotypic analysis of IgM libraries from 16 individuals. To demonstrate the accuracy of corecount, we Sanger sequenced all the heavy chain IG alleles (65 IGHV, 27 IGHD and 7 IGHJ) from one individual from whom we also produced two independent IgM Rep-seq datasets. Genomic analysis revealed that 5 known IGHV and 2 IGHJ sequences are truncated in current reference databases. This dataset of genomically validated alleles and IgM libraries from the same individual provides a useful resource for benchmarking other bioinformatic programs that involve V, D and J assignments and germline inference, and may facilitate the development of AIRR-Seq analysis tools that can take benefit from the availability of more comprehensive reference databases.

Analysis of IGH allele content in a sample group of rheumatoid arthritis patients demonstrates unrevealed population heterogeneity

Immunoglobulin heavy chain (IGH) germline gene variations influence the B cell receptor repertoire, with resulting biological consequences such as shaping our response to infections and altering disease susceptibilities. However, the lack of information on polymorphism frequencies in the IGH loci at the population level makes association studies challenging. Here, we genotyped a pilot group of 30 individuals with rheumatoid arthritis (RA) to examine IGH allele content and frequencies in this group. Eight novel IGHV alleles and one novel IGHJ allele were identified in the study. 15 cases were haplotypable using heterozygous IGHJ6 or IGHD anchors. One variant, IGHV4-34*01_S0742, was found in three out of 30 cases and included a single nucleotide change resulting in a non-canonical recombination signal sequence (RSS) heptamer. This variant allele, shown by haplotype analysis to be non-expressed, was also found in three out of 30 healthy controls and matched a single nucleotide polymorphism (SNP) described in the 1000 Genomes Project (1KGP) collection with frequencies that varied between population groups. Our finding of previously unreported alleles in a relatively small group of individuals with RA illustrates the need for baseline information about IG allelic frequencies in targeted study groups in preparation for future analysis of these genes in disease association studies.

Targeted long-read sequencing facilitates phased diploid assembly and genotyping of the human T cell receptor alpha, delta, and beta loci

T cell receptors (TCRs) recognize peptide fragments presented by the major histocompatibility complex (MHC) and are critical to T cell-mediated immunity. Recent data have indicated that genetic diversity within TCR-encoding gene regions is underexplored, limiting understanding of the impact of TCR loci polymorphisms on TCR function in disease, even though TCR repertoire signatures (1) are heritable and (2) associate with disease phenotypes. To address this, we developed a targeted long-read sequencing approach to generate highly accurate haplotype resolved assemblies of the TCR beta (TRB) and alpha/delta (TRA/D) loci, facilitating the genotyping of all variant types, including structural variants. We validate our approach using two mother-father-child trios and 5 unrelated donors representing multiple populations. This resulted in improved genotyping accuracy and the discovery of 84 undocumented V, D, J, and C alleles, demonstrating the utility of this framework for improving our understanding of TCR diversity and function in disease.

B-Cell Reconstitution After Autologous Hematopoietic Stem Cell Transplantation in Multiple Sclerosis

BACKGROUND AND OBJECTIVES

Autologous hematopoietic stem cell transplantation (aHSCT) is increasingly used to treat aggressive forms of multiple sclerosis (MS). This procedure is believed to result in an immune reset and restoration of a self-tolerant immune system. Immune reconstitution has been extensively studied for T cells, but only to a limited extent for B cells. As increasing evidence suggests an important role of B cells in MS pathogenesis, we sought here to better understand reconstitution and the extent of renewal of the B-cell system after aHSCT in MS.

METHODS

Using longitudinal multidimensional flow cytometry and immunoglobulin heavy chain (IgH) repertoire sequencing following aHSCT with BCNU + Etoposide + Ara-C + Melphalan anti-thymocyte globulin, we analyzed the B-cell compartment in a cohort of 20 patients with MS in defined intervals before and up to 1 year after aHSCT and compared these findings with data from healthy controls.

RESULTS

Total B-cell numbers recovered within 3 months and increased above normal levels 1 year after transplantation, successively shifting from a predominantly transitional to a naive immune phenotype. Memory subpopulations recovered slowly and remained below normal levels with reduced repertoire diversity 1 year after transplantation. Isotype subclass analysis revealed a proportional shift toward IgG1-expressing cells and a reduction in IgG2 cells. Mutation analysis of IgH sequences showed that highly mutated memory B cells and plasma cells may transiently survive conditioning while the analysis of sequence cluster overlap, variable (IGHV) and joining (IGHJ) gene usage and repertoire diversity suggested a renewal of the late posttransplant repertoire. In patients with early cytomegalovirus reactivation, reconstitution of naive and memory B cells was delayed.

DISCUSSION

Our detailed characterization of B-cell reconstitution after aHSCT in MS indicates a reduced reactivation potential of memory B cells up to 1 year after transplantation, which may leave patients susceptible to infection, but may also be an important aspect of its mechanism of action.

Next-generation sequencing for MRD monitoring in B-lineage malignancies: from bench to bedside

Minimal residual disease (MRD) is considered the strongest relevant predictor of prognosis and an effective decision-making factor during the treatment of hematological malignancies. Remarkable breakthroughs brought about by new strategies, such as epigenetic therapy and chimeric antigen receptor-T (CAR-T) therapy, have led to considerably deeper responses in patients than ever, which presents difficulties with the widely applied gold-standard techniques of MRD monitoring. Urgent demands for novel approaches that are ultrasensitive and provide sufficient information have put a spotlight on high-throughput technologies. Recently, advances in methodology, represented by next-generation sequencing (NGS)-based clonality assays, have proven robust and suggestive in numerous high-quality studies and have been recommended by some international expert groups as disease-monitoring modalities. This review demonstrates the applicability of NGS-based clonality assessment for MRD monitoring of B-cell malignancies by summarizing the oncogenesis of neoplasms and the corresponding status of immunoglobulin (IG) rearrangements. Furthermore, we focused on the performance of NGS-based assays compared with conventional approaches and the interpretation of results, revealing directions for improvement and prospects in clinical practice.

Evolution of Anti-SARS-CoV-2 Therapeutic Antibodies

Since the first COVID-19 reports back in December of 2019, this viral infection caused by SARS-CoV-2 has claimed millions of lives. To control the COVID-19 pandemic, the Food and Drug Administration (FDA) and/or European Agency of Medicines (EMA) have granted Emergency Use Authorization (EUA) to nine therapeutic antibodies. Nonetheless, the natural evolution of SARS-CoV-2 has generated numerous variants of concern (VOCs) that have challenged the efficacy of the EUA antibodies. Here, we review the most relevant characteristics of these therapeutic antibodies, including timeline of approval, neutralization profile against the VOCs, selection methods of their variable regions, somatic mutations, HCDR3 and LCDR3 features, isotype, Fc modifications used in the therapeutic format, and epitope recognized on the receptor-binding domain (RBD) of SARS-CoV-2. One of the conclusions of the review is that the EUA therapeutic antibodies that still retain efficacy against new VOCs bind an epitope formed by conserved residues that seem to be evolutionarily conserved as thus, critical for the RBD:hACE-2 interaction. The information reviewed here should help to design new and more efficacious antibodies to prevent and/or treat COVID-19, as well as other infectious diseases.

Immunoglobulin genes, reproductive isolation and vertebrate speciation

Reproductive isolation drives the formation of new species, and many genes contribute to this through Dobzhansky–Muller incompatibilities (DMIs). These incompatibilities occur when gene divergence affects loci encoding interacting products such as receptors and their ligands. We suggest here that the nature of vertebrate immunoglobulin (IG) genes must make them prone to DMIs. The genes of these complex loci form functional genes through the process of recombination, giving rise to a repertoire of heterodimeric receptors of incredible diversity. This repertoire, within individuals and within species, must defend against pathogens but must also avoid pathogenic self‐reactivity. We suggest that this avoidance of autoimmunity is only achieved through a coordination of evolution between heavy‐ and light‐chain genes, and between these genes and the rest of the genome. Without coordinated evolution, the hybrid offspring of two diverging populations will carry a heavy burden of DMIs, resulting in a loss of fitness. Critical incompatibilities could manifest as incompatibilities between a mother and her divergent offspring. During fetal development, biochemical differences between the parents of hybrid offspring could make their offspring a target of the maternal immune system. This hypothesis was conceived in the light of recent insights into the population genetics of IG genes. This has suggested that antibody genes are probably as susceptible to evolutionary forces as other parts of the genome. Further repertoire studies in human and nonhuman species should now help determine whether antibody genes have been part of the evolutionary forces that drive the development of species.

Mystifying Molecular Structure, Expression and Repertoire Diversity of IgM Heavy Chain Genes (Ighμ) in Clarias Catfish and Hybrids: Two Novel Transcripts in Vertebrates

Two novel immunoglobulin heavy chain (Ighμ) transcripts encoding membrane-bound forms of IgM (mIgM) were discovered in bighead catfish, Clarias macrocephalus. The first transcript contains four constant and two transmembrane domains [Cμ1-Cμ2-Cμ3-Cμ4-TM1-TM2] that have never been reported in teleosts, and the second transcript is an unusual mIgM that has never been identified in any vertebrate [Cμ1-(Cδ2-Cδ3-Cδ4-Cδ5)-Cμ2-Cμ3-TM1-TM2]. Fluorescence in situ hybridization (FISH) in bighead catfish, North African catfish (C. gariepinus) and hybrid catfish revealed a single copy of Ighμ in individual parent catfish, while two gene copies were found in diploid hybrid catfish. Intensive sequence analysis demonstrated multiple distinct structural variabilities in the VH domain in Clarias, and hybrid catfish were defined and used to generate diversity with various mechanisms. Expression analysis of Ighμ in Aeromonas hydrophila infection of the head kidney, peripheral blood leukocytes and spleen revealed significantly higher levels in North African catfish and hybrid catfish than in bighead catfish.

A BALB/c IGHV Reference Set, Defined by Haplotype Analysis of Long-Read VDJ-C Sequences From F1 (BALB/c x C57BL/6) Mice

The immunoglobulin genes of inbred mouse strains that are commonly used in models of antibody-mediated human diseases are poorly characterized. This compromises data analysis. To infer the immunoglobulin genes of BALB/c mice, we used long-read SMRT sequencing to amplify VDJ-C sequences from F1 (BALB/c x C57BL/6) hybrid animals. Strain variations were identified in the Ighm and Ighg2b genes, and analysis of VDJ rearrangements led to the inference of 278 germline IGHV alleles. 169 alleles are not present in the C57BL/6 genome reference sequence. To establish a set of expressed BALB/c IGHV germline gene sequences, we computationally retrieved IGHV haplotypes from the IgM dataset. Haplotyping led to the confirmation of 162 BALB/c IGHV gene sequences. A musIGHV398 pseudogene variant also appears to be present in the BALB/cByJ substrain, while a functional musIGHV398 gene is highly expressed in the BALB/cJ substrain. Only four of the BALB/c alleles were also observed in the C57BL/6 haplotype. The full set of inferred BALB/c sequences has been used to establish a BALB/c IGHV reference set, hosted at https://ogrdb.airr-community.org. We assessed whether assemblies from the Mouse Genome Project (MGP) are suitable for the determination of the genes of the IGH loci. Only 37 (43.5%) of the 85 confirmed IMGT-named BALB/c IGHV and 33 (42.9%) of the 77 confirmed non-IMGT IGHV were found in a search of the MGP BALB/cJ genome assembly. This suggests that current MGP assemblies are unsuitable for the comprehensive documentation of germline IGHVs and more efforts will be needed to establish strain-specific reference sets.

Accelerated Evolution by Diversity-Generating Retroelements

Diversity-generating retroelements (DGRs) create vast amounts of targeted, functional diversity by facilitating the rapid evolution of ligand-binding protein domains. Thousands of DGRs have been identified in bacteria, archaea, and their respective viruses. They are broadly distributed throughout the microbial world, with enrichment observed in certain taxa and environments. The diversification machinery works through a novel mechanism termed mutagenic retrohoming, whereby nucleotide sequence information is copied from an invariant DNA template repeat (TR) into an RNA intermediate, selectively mutagenized at TR adenines during cDNA synthesis by a DGR-encoded reverse transcriptase, and transferred to a variable repeat (VR) region within a variable-protein gene (54). This unidirectional flow of information leaves TR-DNA sequences unmodified, allowing for repeated rounds of mutagenic retrohoming to optimize variable-protein function. DGR target genes are often modular and can encode one or more of a wide variety of discrete functional domains appended to a diversifiable ligand-binding motif. Bacterial variable proteins often localize to cell surfaces, although a subset appear to be cytoplasmic, while phage-encoded DGRs commonly diversify tail fiber-associated receptor-binding proteins. Here, we provide a comprehensive review of the mechanism and consequences of accelerated protein evolution by these unique and beneficial genetic elements. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Gene prediction in the immunoglobulin loci

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

Exploring the pathways to chronic lymphocytic leukemia

In chronic lymphocytic leukemia (CLL), increasing knowledge of the biology of the tumor cells has led to transformative improvements in our capacity to assess and treat patients. The dependence of tumor cells on surface immunoglobulin receptor signaling, survival pathways, and accessory cells within the microenvironment has led to a successful double-barreled attack with designer drugs. Studies have revealed that CLL should be classified based on the mutational status of the expressed IGHV sequences into 2 diseases, either unmutated (U) or mutated (M) CLL, each with a distinctive cellular origin, biology, epigenetics/genetics, and clinical behavior. The origin of U-CLL lies among the natural antibody repertoire, and dominance of IGHV1-69 reveals a superantigenic driver. In both U-CLL and M-CLL, a calibrated stimulation of tumor cells by self-antigens apparently generates a dynamic reiterative cycle as cells, protected from apoptosis, transit between blood and tissue sites. But there are differences in outcome, with the balance between proliferation and anergy favoring anergy in M-CLL. Responses are modulated by an array of microenvironmental interactions. Availability of T-cell help is a likely determinant of cell fate, the dependency on which varies between U-CLL and M-CLL, reflecting the different cells of origin, and affecting clinical behavior. Despite such advances, cell-escape strategies, Richter transformation, and immunosuppression remain as challenges, which only may be met by continued research into the biology of CLL.

Diversity of Functionally Distinct Clonal Sets of Human Conventional Memory B Cells That Bind Staphylococcal Protein A

Staphylococcus aureus, a common cause of serious and often fatal infections, is well-armed with secreted factors that disarm host immune defenses. Highly expressed in vivo during infection, Staphylococcal protein A (SpA) is reported to also contribute to nasal colonization that can be a prelude to invasive infection. Co-evolution with the host immune system has provided SpA with an Fc-antibody binding site, and a Fab-binding site responsible for non-immune superantigen interactions via germline-encoded surfaces expressed on many human BCRs. We wondered whether the recurrent exposures to S. aureus commonly experienced by adults, result in the accumulation of memory B-cell responses to other determinants on SpA. We therefore isolated SpA-specific class-switched memory B cells, and characterized their encoding VH : VL antibody genes. In SpA-reactive memory B cells, we confirmed a striking bias in usage for VH genes, which retain the surface that mediates the SpA-superantigen interaction. We postulate these interactions reflect co-evolution of the host immune system and SpA, which during infection results in immune recruitment of an extraordinarily high prevalence of B cells in the repertoire that subverts the augmentation of protective defenses. Herein, we provide the first evidence that human memory responses are supplemented by B-cell clones, and circulating-antibodies, that bind to SpA determinants independent of the non-immune Fc- and Fab-binding sites. In parallel, we demonstrate that healthy individuals, and patients recovering from S. aureus infection, both have circulating antibodies with these conventional binding specificities. These findings rationalize the potential utility of incorporating specially engineered SpA proteins into a protective vaccine.

Genetic Variability in Molecular Pathways Implicated in Alzheimer's Disease: A Comprehensive Review

Alzheimer's disease (AD) is a complex neurodegenerative disease, affecting a significant part of the population. The majority of AD cases occur in the elderly with a typical age of onset of the disease above 65 years. AD presents a major burden for the healthcare system and since population is rapidly aging, the burden of the disease will increase in the future. However, no effective drug treatment for a full-blown disease has been developed to date. The genetic background of AD is extensively studied; numerous genome-wide association studies (GWAS) identified significant genes associated with increased risk of AD development. This review summarizes more than 100 risk loci. Many of them may serve as biomarkers of AD progression, even in the preclinical stage of the disease. Furthermore, we used GWAS data to identify key pathways of AD pathogenesis: cellular processes, metabolic processes, biological regulation, localization, transport, regulation of cellular processes, and neurological system processes. Gene clustering into molecular pathways can provide background for identification of novel molecular targets and may support the development of tailored and personalized treatment of AD.

Immunoglobulin germline gene variation and its impact on human disease

Immunoglobulins (Ig) play an important role in the immune system both when expressed as antigen receptors on the cell surface of B cells and as antibodies secreted into extracellular fluids. The advent of high-throughput sequencing methods has enabled the investigation of human Ig repertoires at unprecedented depth. This has led to the discovery of many previously unreported germline Ig alleles. Moreover, it is becoming clear that convergent and stereotypic antibody responses are common where different individuals recognise defined antigenic epitopes with the use of the same Ig V genes. Thus, germline V gene variation is increasingly being linked to the differential capacity of generating an effective immune response, which might lead to varying disease susceptibility. Here, we review recent evidence of how germline variation in Ig genes impacts the Ig repertoire and its subsequent effects on the adaptive immune response in vaccination, infection, and autoimmunity.

Germline immunoglobulin genes: Disease susceptibility genes hidden in plain sight?

Immunoglobulin genes are rarely considered as disease susceptibility genes despite their obvious and central contributions to immune function. This appears to be a consequence of historical views on antibody repertoire formation that no longer stand, and of difficulties that until recently surrounded the documentation of the suite of antibody genes in any individual. If these important genes are to be accessible to GWAS studies, allelic variation within the human population needs to be better documented, and a curated set of genomic variations associated with antibody genes needs to be formulated. Repertoire studies arising from the COVID-19 pandemic provide an opportunity to meet these needs, and may provide insights into the profound variability that is seen in outcomes to this infection.

A Novel Framework for Characterizing Genomic Haplotype Diversity in the Human Immunoglobulin Heavy Chain Locus

An incomplete ascertainment of genetic variation within the highly polymorphic immunoglobulin heavy chain locus (IGH) has hindered our ability to define genetic factors that influence antibody-mediated processes. Due to locus complexity, standard high-throughput approaches have failed to accurately and comprehensively capture IGH polymorphism. As a result, the locus has only been fully characterized two times, severely limiting our knowledge of human IGH diversity. Here, we combine targeted long-read sequencing with a novel bioinformatics tool, IGenotyper, to fully characterize IGH variation in a haplotype-specific manner. We apply this approach to eight human samples, including a haploid cell line and two mother-father-child trios, and demonstrate the ability to generate high-quality assemblies (>98% complete and >99% accurate), genotypes, and gene annotations, identifying 2 novel structural variants and 15 novel IGH alleles. We show multiplexing allows for scaling of the approach without impacting data quality, and that our genotype call sets are more accurate than short-read (>35% increase in true positives and >97% decrease in false-positives) and array/imputation-based datasets. This framework establishes a desperately needed foundation for leveraging IG genomic data to study population-level variation in antibody-mediated immunity, critical for bettering our understanding of disease risk, and responses to vaccines and therapeutics.

Immunoglobulins or Antibodies: IMGT® Bridging Genes, Structures and Functions

IMGT®, the international ImMunoGeneTics® information system founded in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science at the interface between immunogenetics and bioinformatics. For the first time, the immunoglobulin (IG) or antibody and T cell receptor (TR) genes were officially recognized as 'genes' as well as were conventional genes. This major breakthrough has allowed the entry, in genomic databases, of the IG and TR variable (V), diversity (D) and joining (J) genes and alleles of Homo sapiens and of other jawed vertebrate species, based on the CLASSIFICATION axiom. The second major breakthrough has been the IMGT unique numbering and the IMGT Collier de Perles for the V and constant (C) domains of the IG and TR and other proteins of the IG superfamily (IgSF), based on the NUMEROTATION axiom. IMGT-ONTOLOGY axioms and concepts bridge genes, sequences, structures and functions, between biological and computational spheres in the IMGT® system (Web resources, databases and tools). They provide the IMGT Scientific chart rules to identify, to describe and to analyse the IG complex molecular data, the huge diversity of repertoires, the genetic (alleles, allotypes, CNV) polymorphisms, the IG dual function (paratope/epitope, effector properties), the antibody humanization and engineering.

Polymorphisms in human immunoglobulin heavy chain variable genes and their upstream regions

Germline variations in immunoglobulin genes influence the repertoire of B cell receptors and antibodies, and such polymorphisms may impact disease susceptibility. However, the knowledge of the genomic variation of the immunoglobulin loci is scarce. Here, we report 25 potential novel germline IGHV alleles as inferred from rearranged naïve B cell cDNA repertoires of 98 individuals. Thirteen novel alleles were selected for validation, out of which ten were successfully confirmed by targeted amplification and Sanger sequencing of non-B cell DNA. Moreover, we detected a high degree of variability upstream of the V-REGION in the 5′UTR, L-PART1 and L-PART2 sequences, and found that identical V-REGION alleles can differ in upstream sequences. Thus, we have identified a large genetic variation not only in the V-REGION but also in the upstream sequences of IGHV genes. Our findings provide a new perspective for annotating immunoglobulin repertoire sequencing data.

RAbHIT: R Antibody Haplotype Inference Tool

SUMMARY

Antibody haplotype inference (chromosomal phasing) may have clinical implications for the identification of genetic predispositions to diseases. Yet, our knowledge of the genomic loci encoding for the variable regions of the antibody is only partial, mostly due to the challenge of aligning short reads from genome sequencing to these highly repetitive loci. A powerful approach to infer the content of these loci relies on analyzing repertoires of rearranged V(D)J sequences. We present here RAbHIT, an R Haplotype Antibody Inference Tool, that implements a novel algorithm to infer V(D)J haplotypes by adapting a Bayesian framework. RAbHIT offers inference of haplotype and gene deletions. It may be applied to sequences from naïve and non-naïve B-cells, sequenced by different library preparation protocols.

AVAILABILITY AND IMPLEMENTATION

RAbHIT is freely available for academic use from comprehensive R archive network (CRAN) (https://cran.r-project.org/web/packages/rabhit/) under CC BY-SA 4.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Plasma proteomic profiles differ between European and North American myotid bats colonized by Pseudogymnoascus destructans

Emerging fungal diseases have become challenges for wildlife health and conservation. North American hibernating bat species are threatened by the psychrophilic fungus Pseudogymnoascus destructans (Pd) causing the disease called white-nose syndrome (WNS) with unprecedented mortality rates. The fungus is widespread in North America and Europe, however, disease is not manifested in European bats. Differences in epidemiology and pathology indicate an evolution of resistance or tolerance mechanisms towards Pd in European bats. We compared the proteomic profile of blood plasma in healthy and Pd-colonized European Myotis myotis and North American Myotis lucifugus in order to identify pathophysiological changes associated with Pd colonization, which might also explain the differences in bat survival. Expression analyses of plasma proteins revealed differences in healthy and Pd-colonized M. lucifugus, but not in M. myotis. We identified differentially expressed proteins for acute phase response, constitutive and adaptive immunity, oxidative stress defence, metabolism and structural proteins of exosomes and desmosomes, suggesting a systemic response against Pd in North American M. lucifugus but not European M. myotis. The differences in plasma proteomic profiles between European and North American bat species colonized by Pd suggest European bats have evolved tolerance mechanisms towards Pd infection.

Genotyping and Copy Number Analysis of Immunoglobin Heavy Chain Variable Genes Using Long Reads

One of the remaining challenges to describing an individual's genetic variation lies in the highly heterogeneous and complex genomic regions that impede the use of classical reference-guided mapping and assembly approaches. Once such region is the Immunoglobulin heavy chain locus (IGH), which is critical for the development of antibodies and the adaptive immune system. We describe ImmunoTyper, the first PacBio-based genotyping and copy number calling tool specifically designed for IGH V genes (IGHV). We demonstrate that ImmunoTyper's multi-stage clustering and combinatorial optimization approach represents the most comprehensive IGHV genotyping approach published to date, through validation using gold-standard IGH reference sequence. This preliminary work establishes the feasibility of fine-grained genotype and copy number analysis using error-prone long reads in complex multi-gene loci and opens the door for in-depth investigation into IGHV heterogeneity using accessible and increasingly common whole-genome sequence.

Molecular profiling of immunoglobulin heavy-chain gene rearrangements unveils new potential prognostic markers for multiple myeloma patients

Multiple myeloma is a heterogeneous disease whose pathogenesis has not been completely elucidated. Although B-cell receptors play a crucial role in myeloma pathogenesis, the impact of clonal immunoglobulin heavy-chain features in the outcome has not been extensively explored. Here we present the characterization of complete heavy-chain gene rearrangements in 413 myeloma patients treated in Spanish trials, including 113 patients characterized by next-generation sequencing. Compared to the normal B-cell repertoire, gene selection was biased in myeloma, with significant overrepresentation of IGHV3, IGHD2 and IGHD3, as well as IGHJ4 gene groups. Hypermutation was high in our patients (median: 8.8%). Interestingly, regarding patients who are not candidates for transplantation, a high hypermutation rate (≥7%) and the use of IGHD2 and IGHD3 groups were associated with improved prognostic features and longer survival rates in the univariate analyses. Multivariate analysis revealed prolonged progression-free survival rates for patients using IGHD2/IGHD3 groups (HR: 0.552, 95% CI: 0.361-0.845, p = 0.006), as well as prolonged overall survival rates for patients with hypermutation ≥7% (HR: 0.291, 95% CI: 0.137-0.618, p = 0.001). Our results provide new insights into the molecular characterization of multiple myeloma, highlighting the need to evaluate some of these clonal rearrangement characteristics as new potential prognostic markers.

Regulatory T cells and co-evolution of allele-specific MHC recognition by the TCR

What is the evolutionary mechanism for the TCR-MHC-conserved interaction? We extend Dembic's model (Dembic Z. In, Scand J Immunol e12806, 2019) of thymus positive selection for high-avidity anti-self-MHC Tregs among double (CD4 + CD8+)-positive (DP) developing thymocytes. This model is based on competition for self-MHC (+ Pep) complexes presented on cortical epithelium. Such T cells exit as CD4 + CD25+FoxP3 + thymic-derived Tregs (tTregs). The other positively selected DP T cells are then negatively selected on medulla epithelium removing high-avidity anti-self-MHC + Pep as T cells commit to CD4 + or CD8 + lineages. The process is likened to the competitive selection and affinity maturation in Germinal Centre for the somatic hypermutation (SHM) of rearranged immunoglobulin (Ig) variable region (V[D]Js) of centrocytes bearing antigen-specific B cell receptors (BCR). We now argue that the same direct SHM processes for TCRs occur in post-antigenic Germinal Centres, but now occurring in peripheral pTregs. This model provides a potential solution to a long-standing problem previously recognized by Cohn and others (Cohn M, Anderson CC, Dembic Z. In, Scand J Immunol e12790, 2019) of how co-evolution occurs of species-specific MHC alleles with the repertoire of their germline TCR V counterparts. We suggest this is not by 'blind', slow, and random Darwinian natural selection events, but a rapid structured somatic selection vertical transmission process. The pTregs bearing somatic TCR V mutant genes then, on arrival in reproductive tissues, can donate their TCR V sequences via soma-to-germline feedback as discussed in this journal earlier. (Steele EJ, Lindley RA. In, Scand J Immunol e12670, 2018) The high-avidity tTregs also participate in the same process to maintain a biased, high-avidity anti-self-MHC germline V repertoire.

A comparison of immunoglobulin IGHV, IGHD and IGHJ genes in wild-derived and classical inbred mouse strains

The genomes of classical inbred mouse strains include genes derived from all three major subspecies of the house mouse, Mus musculus. We recently posited that genetic diversity in the immunoglobulin heavy chain (IGH) gene loci of C57BL/6 and BALB/c mice reflects differences in subspecies origin. To investigate this hypothesis, we conducted high-throughput sequencing of IGH gene rearrangements to document IGH variable (IGHV), joining (IGHJ) and diversity (IGHD) genes in four inbred wild-derived mouse strains (CAST/EiJ, LEWES/EiJ, MSM/MsJ and PWD/PhJ) and a single disease model strain (NOD/ShiLtJ), collectively representing genetic backgrounds of several major mouse subspecies. A total of 341 germline IGHV sequences were inferred in the wild-derived strains, including 247 not curated in the international ImMunoGeneTics information system. By contrast, 83/84 inferred NOD IGHV genes had previously been observed in C57BL/6 mice. Variability among the strains examined was observed for only a single IGHJ gene, involving a description of a novel allele. By contrast, unexpected variation was found in the IGHD gene loci, with four previously unreported IGHD gene sequences being documented. Very few IGHV sequences of C57BL/6 and BALB/c mice were shared with strains representing major subspecies, suggesting that their IGH loci may be complex mosaics of genes of disparate origins. This suggests a similar level of diversity is likely present in the IGH loci of other classical inbred strains. This must now be documented if we are to properly understand interstrain variation in models of antibody-mediated disease.

Analysis of the B cell receptor repertoire in six immune-mediated diseases

B cells are important in the pathogenesis of many, and perhaps all, immune-mediated diseases (IMDs). Each B cell expresses a single B cell receptor (BCR)¹, with the diverse range of BCRs expressed by an individual’s total B cell population being termed the “BCR repertoire”. Our understanding of the BCR repertoire in the context of IMDs is incomplete, and defining this could reveal new insights into pathogenesis and therapy. We therefore compared the BCR repertoire in systemic lupus erythematosus (SLE), ANCA-associated vasculitis (AAV), Crohn’s disease (CD), Behçet’s disease (BD), eosinophilic granulomatosis with polyangiitis (EGPA) and IgA vasculitis (IgAV), analysing BCR clonality, and immunoglobulin heavy chain gene (IGHV) and, in particular, isotype usage. An IgA-dominated increased clonality in SLE and CD, together with skewed IGHV gene usage in these and other diseases, suggested a microbial contribution to pathogenesis. Different immunosuppressive treatment had specific and distinct impacts on the repertoire; B cells persisting after rituximab were predominately isotype-switched and clonally expanded, the inverse of those persisting after mycophenolate mofetil. A comparative analysis of the BCR repertoire in immune-mediated disease reveals a complex B cell architecture, providing a platform for understanding pathological mechanisms and designing treatment strategies.

Phage Display Libraries for Antibody Therapeutic Discovery and Development

Phage display technology has played a key role in the remarkable progress of discovering and optimizing antibodies for diverse applications, particularly antibody-based drugs. This technology was initially developed by George Smith in the mid-1980s and applied by John McCafferty and Gregory Winter to antibody engineering at the beginning of 1990s. Here, we compare nine phage display antibody libraries published in the last decade, which represent the state of the art in the discovery and development of therapeutic antibodies using phage display. We first discuss the quality of the libraries and the diverse types of antibody repertoires used as substrates to build the libraries, i.e., naïve, synthetic, and semisynthetic. Second, we review the performance of the libraries in terms of the number of positive clones per panning, hit rate, affinity, and developability of the selected antibodies. Finally, we highlight current opportunities and challenges pertaining to phage display platforms and related display technologies.

Per-sample immunoglobulin germline inference from B cell receptor deep sequencing data

The collection of immunoglobulin genes in an individual's germline, which gives rise to B cell receptors via recombination, is known to vary significantly across individuals. In humans, for example, each individual has only a fraction of the several hundred known V alleles. Furthermore, the currently-accepted set of known V alleles is both incomplete (particularly for non-European samples), and contains a significant number of spurious alleles. The resulting uncertainty as to which immunoglobulin alleles are present in any given sample results in inaccurate B cell receptor sequence annotations, and in particular inaccurate inferred naive ancestors. In this paper we first show that the currently widespread practice of aligning each sequence to its closest match in the full set of IMGT alleles results in a very large number of spurious alleles that are not in the sample's true set of germline V alleles. We then describe a new method for inferring each individual's germline gene set from deep sequencing data, and show that it improves upon existing methods by making a detailed comparison on a variety of simulated and real data samples. This new method has been integrated into the partis annotation and clonal family inference package, available at https://github.com/psathyrella/partis, and is run by default without affecting overall run time.

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.

Human adaptive immune receptor repertoire analysis-Past, present, and future

The genes encoding adaptive immune antigen receptors, namely the immunoglobulins expressed in membrane-bound or secreted forms by B cells, and the cell surface T cell receptors, are unique in human biology because they are generated by combinatorial rearrangement of the genomic DNA. The diversity of receptors so generated in populations of lymphocytes enables the human immune system to recognize antigens expressed by pathogens, but also underlies the pathological specificity of autoimmune diseases and the mistargeted immunity in allergies. Several recent technological developments, foremost among them the invention of high-throughput DNA sequencing instruments, have enabled much deeper and thorough evaluation of clones of human B cells and T cells and the antigen receptors they express during physiological and pathogenic immune responses. The evolutionary struggles between host adaptive immune responses and populations of pathogens are now open to greater scrutiny, elucidation of the underlying reasons for successful or failed immunity, and potential predictive modeling, than ever before. Here we give an overview of the foundations, recent progress, and future prospects in this dynamic area of research.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Human colorectal cancer cells frequently express IgG and display unique Ig repertoire

BACKGROUND

There is growing evidence proving that many human carcinomas, including colon cancer, can overexpress immunoglobulin (Ig); the non B cancer cell-derived Ig usually displayed unique V(D)J rearrangement pattern that are distinct from B cell-derived Ig. Especially, the cancer-derived Ig plays important roles in cancer initiation, progression, and metastasis. However, it still remains unclear if the colon cancer-derived Ig can display unique V(D)J pattern and sequencing, which can be used as novel target for colon cancer therapy.

AIM

To investigate the Ig repertoire features expressed in human colon cancer cells.

METHODS

Seven cancerous tissue samples of colon adenocarcinoma and corresponding noncancerous tissue samples were sorted by fluorescence-activated cell sorting using epithelial cell adhesion molecule as a marker for epithelial cells. Ig repertoire sequencing was used to analyze the expression profiles of all 5 classes of Ig heavy chains (IgH) and the Ig repertoire in colon cancer cells and corresponding normal epithelial cells.

RESULTS

We found that all 5 IgH classes can be expressed in both colon cancer cells and normal epithelial cells. Surprisingly, unlike the normal colonic epithelial cells that expressed 5 Ig classes, our results suggested that cancer cells most prominently express IgG. Next, we found that the usage of Ig in cancer cells caused the expression of some unique Ig repertoires compared to normal cells. Some V_H segments, such as V_H3-7, have been used in cancer cells, and V_H3-74 was frequently present in normal epithelial cells. Moreover, compared to the normal cell-derived Ig, most cancer cell-derived Ig showed unique V_HDJ_H patterns. Importantly, even if the same V_HDJ_H pattern was seen in cancer cells and normal cells, cancer cell-derived IgH always displayed distinct hypermutation hot points.

CONCLUSION

We found that colon cancer cells could frequently express IgG and unique IgH repertoires, which may be involved in carcinogenesis of colon cancer. The unique IgH repertoire has the potential to be used as a novel target in immune therapy for colon cancer.

Worldwide genetic variation of the IGHV and TRBV immune receptor gene families in humans

This article presents a comprehensive study of the IGHV and TRBV gene families in a globally diverse sample of humans and shows that the two gene families exhibit starkly different patterns of variation.

The immunoglobulin heavy variable (IGHV) and T cell beta variable (TRBV) loci are among the most complex and variable regions in the human genome. Generated through a process of gene duplication/deletion and diversification, these loci can vary extensively between individuals in copy number and contain genes that are highly similar, making their analysis technically challenging. Here, we present a comprehensive study of the functional gene segments in the IGHV and TRBV loci, quantifying their copy number and single-nucleotide variation in a globally diverse sample of 109 (IGHV) and 286 (TRBV) humans from over a 100 populations. We find that the IGHV and TRBV gene families exhibit starkly different patterns of variation. In addition to providing insight into the different evolutionary paths of the IGHV and TRBV loci, our results are also important to the adaptive immune repertoire sequencing community, where the lack of frequencies of common alleles and copy number variants is hampering existing analytical pipelines.

Identification of Subject-Specific Immunoglobulin Alleles From Expressed Repertoire Sequencing Data

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies.

Mosaic deletion patterns of the human antibody heavy chain gene locus shown by Bayesian haplotyping

Analysis of antibody repertoires by high-throughput sequencing is of major importance in understanding adaptive immune responses. Our knowledge of variations in the genomic loci encoding immunoglobulin genes is incomplete, resulting in conflicting VDJ gene assignments and biased genotype and haplotype inference. Haplotypes can be inferred using IGHJ6 heterozygosity, observed in one third of the people. Here, we propose a robust novel method for determining VDJ haplotypes by adapting a Bayesian framework. Our method extends haplotype inference to IGHD- and IGHV-based analysis, enabling inference of deletions and copy number variations in the entire population. To test this method, we generated a multi-individual data set of naive B-cell repertoires, and found allele usage bias, as well as a mosaic, tiled pattern of deleted IGHD and IGHV genes. The inferred haplotypes may have clinical implications for genetic disease predispositions. Our findings expand the knowledge that can be extracted from antibody repertoire sequencing data.

Meta-analysis of Alzheimer's disease on 9,751 samples from Norway and IGAP study identifies four risk loci

A large fraction of genetic risk factors for Alzheimer's Disease (AD) is still not identified, limiting the understanding of AD pathology and study of therapeutic targets. We conducted a genome-wide association study (GWAS) of AD cases and controls of European descent from the multi-center DemGene network across Norway and two independent European cohorts. In a two-stage process, we first performed a meta-analysis using GWAS results from 2,893 AD cases and 6,858 cognitively normal controls from Norway and 25,580 cases and 48,466 controls from the International Genomics of Alzheimer's Project (IGAP), denoted the discovery sample. Second, we selected the top hits (p < 1 × 10^-6) from the discovery analysis for replication in an Icelandic cohort consisting of 5,341 cases and 110,008 controls. We identified a novel genomic region with genome-wide significant association with AD on chromosome 4 (combined analysis OR = 1.07, p = 2.48 x 10^-8). This finding implicated HS3ST1, a gene expressed throughout the brain particularly in the cerebellar cortex. In addition, we identified IGHV1-68 in the discovery sample, previously not associated with AD. We also associated USP6NL/ECHDC3 and BZRAP1-AS1 to AD, confirming findings from a follow-up transethnic study. These new gene loci provide further evidence for AD as a polygenic disorder, and suggest new mechanistic pathways that warrant further investigation.

Exploiting antigen receptor information to quantify index switching in single-cell transcriptome sequencing experiments

By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms have been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as "index switching", which causes a proportion of the reads to be assigned to an incorrect sample. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-specific expression of antigen receptors on immune cells to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3000 and HiSeq 4000. By utilizing the unique antigen receptor expression, we could quantify the spread-of-signal from many different wells (n = 55 from total of three batches) due to index switching. Based on full-length T cell receptor (TCR) sequences from all samples reconstructed by TraCeR and TCR gene expression quantified by Kallisto, we found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 3.9% (interquartile range (IQR): 1.7%-7.3%). We did not detect any consistent patterns of certain indices to be more prone to switching than others, suggesting that index switching is a stochastic process. Our results confirm that index switching is a problem that affects samples run in multiplexed libraries on Illumina HiSeq 3000 and HiSeq 4000 platforms.

Immunoglobulin Light Chain Gene Rearrangements, Receptor Editing and the Development of a Self-Tolerant Antibody Repertoire

Discussion of the antibody repertoire usually emphasizes diversity, but a conspicuous feature of the light chain repertoire is its lack of diversity. The diversity of reported allelic variants of germline light chain genes is also limited, even in well-studied species. In this review, the implications of this lack of diversity are considered. We explore germline and rearranged light chain genes in a variety of species, with a particular focus on human and mouse genes. The importance of the number, organization and orientation of the genes for the control of repertoire development is discussed, and we consider how primary rearrangements and receptor editing together shape the expressed light chain repertoire. The resulting repertoire is dominated by just a handful of IGKV and IGLV genes. It has been hypothesized that an important function of the light chain is to guard against self-reactivity, and the role of secondary rearrangements in this process could explain the genomic organization of the light chain genes. It could also explain why the light chain repertoire is so limited. Heavy and light chain genes may have co-evolved to ensure that suitable light chain partners are usually available for each heavy chain that forms early in B cell development. We suggest that the co-evolved loci of the house mouse often became separated during the inbreeding of laboratory mice, resulting in new pairings of loci that are derived from different sub-species of the house mouse. A resulting vulnerability to self-reactivity could explain at least some mouse models of autoimmune disease.

Combined Influence of B-Cell Receptor Rearrangement and Somatic Hypermutation on B-Cell Class-Switch Fate in Health and in Chronic Lymphocytic Leukemia

A diverse B-cell receptor (BCR) repertoire is required to bind a wide range of antigens. BCRs are generated through genetic recombination and can be diversified through somatic hypermutation (SHM) or class-switch recombination (CSR). Patterns of repertoire diversity can vary substantially between different health conditions. We use isotype-resolved BCR sequencing to compare B-cell evolution and class-switch fate in healthy individuals and in patients with chronic lymphocytic leukemia (CLL). We show that the patterns of SHM and CSR in B-cells from healthy individuals are distinct from CLL. We identify distinct properties of clonal expansion that lead to the generation of antibodies of different classes in healthy, malignant, and non-malignant CLL BCR repertoires. We further demonstrate that BCR diversity is affected by relationships between antibody variable and constant regions leading to isotype-specific signatures of variable gene usage. This study provides powerful insights into the mechanisms underlying the evolution of the adaptive immune responses in health and their aberration during disease.

Gut and Liver B Cells of Common Clonal Origin in Primary Sclerosing Cholangitis-Inflammatory Bowel Disease

B cells express an antigen‐specific B‐cell receptor (BCR) and may contribute to liver inflammation by recognizing shared antigens in the gut and liver. Herein, we used high‐throughput BCR sequencing of the immunoglobulin heavy chain, specifically the complementarity‐determining region 3 (CDR3), to characterize the B‐cell repertoire of freshly‐frozen paired gut and liver tissue samples from patients with primary sclerosing cholangitis (PSC) and concurrent inflammatory bowel disease (IBD) (PSC‐IBD, n = 10) and paired formalin‐fixed paraffin‐embedded (FFPE) tumor‐adjacent normal colon and liver tissue from patients with colorectal liver metastases (controls, n = 10). We observed significantly greater numbers of B cells (P < 0.01) and unique B‐cell clonotypes (P < 0.05) in gut samples compared to liver samples of patients with PSC‐IBD, whereas BCR sequences in FFPE normal gut and liver samples were nearly absent (14 ± 5 clonotypes; mean ± SD; n = 20). In PSC‐IBD, an average of 8.3% (range, 1.6%‐18.0%) of B‐cell clonotypes were found to overlap paired gut and liver samples following the exclusion of memory clonotypes reported in the blood of healthy controls. Overlapping gut and liver clonotypes showed stronger evidence of antigen‐driven activation compared to non‐overlapping clonotypes, including shorter CDR3 lengths and higher counts of somatic hypermutation (P < 0.0001). Conclusion: A proportion of gut and liver B cells originate from a common clonal origin (i.e., likely to recognize the same antigen) in patients with PSC which suggests B‐cell antigens are shared across the gut–liver axis. (Hepatology Communications 2018; 00:000‐000)

Immunoglobulin heavy chain variable region analysis in dairy goats

Based on the goat genome database, we have annotated the genomic organization of the goat immunoglobulin heavy chain variable region. The goat IgH locus is present on seven genome scaffolds, and contains ten V_H, three D_H and six J_H segments. After the exclusion of three shorter segments, the V_H genes were divided into two gene families based on sequence similarity. By analyzing the IgH cDNA sequences, we further identified that V_H2 (54.2%), D_H1 (61.7%) and J_H1 (60.5%) segments were most frequently utilized in the expression of the immunoglobulin variable region, and that point mutations introduced by somatic hypermutation were the major mutation present in these expressed variable region. Compared with human and horses, D_H-D_H fusion occurred at a higher frequency in goat V(D)J recombination. These results provided variable insights into goat immunoglobulin heavy chain variable region genome loci and repertoire diversity.