IGHV1-69 polymorphism modulates anti-influenza antibody repertoires, correlates with IGHV utilization shifts and varies by ethnicity

Maturation of germinal center B cells after influenza virus vaccination in humans

Germinal centers (GC) are microanatomical lymphoid structures where affinity-matured memory B cells and long-lived bone marrow plasma cells are primarily generated. It is unclear how the maturation of B cells within the GC impacts the breadth and durability of B cell responses to influenza vaccination in humans. We used fine needle aspiration of draining lymph nodes to longitudinally track antigen-specific GC B cell responses to seasonal influenza vaccination. Antigen-specific GC B cells persisted for at least 13 wk after vaccination in two out of seven individuals. Monoclonal antibodies (mAbs) derived from persisting GC B cell clones exhibit enhanced binding affinity and breadth to influenza hemagglutinin (HA) antigens compared with related GC clonotypes isolated earlier in the response. Structural studies of early and late GC-derived mAbs from one clonal lineage in complex with H1 and H5 HAs revealed an altered binding footprint. Our study shows that inducing sustained GC reactions after influenza vaccination in humans supports the maturation of responding B cells.

The peculiar features, diversity and impact of citrulline-reactive autoantibodies

Since entering the stage 25 years ago as a highly specific serological biomarker for rheumatoid arthritis, anti-citrullinated protein antibodies (ACPAs) have been a topic of extensive research. This hallmark B cell response arises years before disease onset, displays interpatient autoantigen variability, and is associated with poor clinical outcomes. Technological and scientific advances have revealed broad clonal diversity and intriguing features including high levels of somatic hypermutation, variable-domain N-linked glycosylation, hapten-like peptide interactions, and clone-specific multireactivity to citrullinated, carbamylated and acetylated epitopes. ACPAs have been found in different isotypes and subclasses, in both circulation and tissue, and are secreted by both plasmablasts and long-lived plasma cells. Notably, although some disease-promoting features have been reported, results now demonstrate that certain monoclonal ACPAs therapeutically block arthritis and inflammation in mouse models. A wealth of functional studies using patient-derived polyclonal and monoclonal antibodies have provided evidence for pathogenic and protective effects of ACPAs in the context of arthritis. To understand the roles of ACPAs, one needs to consider their immunological properties by incorporating different facets such as rheumatoid arthritis B cell biology, environmental triggers and chronic antigen exposure. The emerging picture points to a complex role of citrulline-reactive autoantibodies, in which the diversity and dynamics of antibody clones could determine clinical progression and manifestations.

Extraction of the CDRH3 sequence of the mouse antibody repertoire selected upon influenza virus infection by subtraction of the background antibody repertoire

Historically, antibody reactivity to pathogens and vaccine antigens has been evaluated using serological measurements of antigen-specific antibodies. However, it is difficult to evaluate all antibodies that contribute to various functions in a single assay, such as the measurement of the neutralizing antibody titer. Bulk antibody repertoire analysis using next-generation sequencing is a comprehensive method for analyzing the overall antibody response; however, it is unreliable for estimating antigen-specific antibodies due to individual variation. To address this issue, we propose a method to subtract the background signal from the repertoire of data of interest. In this study, we analyzed changes in antibody diversity and inferred the heavy-chain complementarity-determining region 3 (CDRH3) sequences of antibody clones that were selected upon influenza virus infection in a mouse model using bulk repertoire analysis. A decrease in the diversity of the antibody repertoire was observed upon viral infection, along with an increase in neutralizing antibody titers. Using kernel density estimation of sequences in a high-dimensional sequence space with background signal subtraction, we identified several clusters of CDRH3 sequences induced upon influenza virus infection. Most of these repertoires were detected more frequently in infected mice than in uninfected control mice, suggesting that infection-specific antibody sequences can be extracted using this method. Such an accurate extraction of antigen- or infection-specific repertoire information will be a useful tool for vaccine evaluation in the future.

IMPORTANCE

As specific interactions between antigens and cell-surface antibodies trigger the proliferation of B-cell clones, the frequency of each antibody sequence in the samples reflects the size of each clonal population. Nevertheless, it is extremely difficult to extract antigen-specific antibody sequences from the comprehensive bulk antibody sequences obtained from blood samples due to repertoire bias influenced by exposure to dietary antigens and other infectious agents. This issue can be addressed by subtracting the background noise from the post-immunization or post-infection repertoire data. In the present study, we propose a method to quantify repertoire data from comprehensive repertoire data. This method allowed subtraction of the background repertoire, resulting in more accurate extraction of expanded antibody repertoires upon influenza virus infection. This accurate extraction of antigen- or infection-specific repertoire information is a useful tool for vaccine evaluation.

Human immunoglobulin gene allelic variation impacts germline-targeting vaccine priming

Vaccine priming immunogens that activate germline precursors for broadly neutralizing antibodies (bnAbs) have promise for development of precision vaccines against major human pathogens. In a clinical trial of the eOD-GT8 60mer germline-targeting immunogen, higher frequencies of vaccine-induced VRC01-class bnAb-precursor B cells were observed in the high dose compared to the low dose group. Through immunoglobulin heavy chain variable (IGHV) genotyping, statistical modeling, quantification of IGHV1-2 allele usage and B cell frequencies in the naive repertoire for each trial participant, and antibody affinity analyses, we found that the difference between dose groups in VRC01-class response frequency was best explained by IGHV1-2 genotype rather than dose and was most likely due to differences in IGHV1-2 B cell frequencies for different genotypes. The results demonstrate the need to define population-level immunoglobulin allelic variations when designing germline-targeting immunogens and evaluating them in clinical trials.

The Role of Biological Sex in Pre-Clinical (Mouse) mRNA Vaccine Studies

In this study, we consider the influence of biological sex-specific immune responses on the assessment of mRNA vaccines in pre-clinical murine studies. Recognising the established disparities in immune function attributed to genetic and hormonal differences between individuals of different biological sexes, we compared the mRNA expression and immune responses in mice of both biological sexes after intramuscular injection with mRNA incorporated within lipid nanoparticles. Regarding mRNA expression, no significant difference in protein (luciferase) expression at the injection site was observed between female and male mice following intramuscular administration; however, we found that female BALB/c mice exhibit significantly greater total IgG responses across the concentration range of mRNA lipid nanoparticles (LNPs) in comparison to their male counterparts. This study not only contributes to the scientific understanding of mRNA vaccine evaluation but also emphasizes the importance of considering biological sex in vaccine study designs during pre-clinical evaluation in murine studies.

Linear epitopes of bony fish β-parvalbumins

Introduction

Fish β-parvalbumins are common targets of allergy-causing immunity. The nature of antibody responses to such allergens determines the biological outcome following exposure to fish. Specific epitopes on these allergens recognised by antibodies are incompletely characterised.

Methods

High-content peptide microarrays offer a solution to the identification of linear epitopes recognised by antibodies. We characterized IgG and IgG4 recognition of linear epitopes of fish β-parvalbumins defined in the WHO/IUIS allergen database as such responses hold the potential to counter an allergic reaction to these allergens. Peripheral blood samples, collected over three years, of 15 atopic but not fish-allergic subjects were investigated using a microarray platform that carried every possible 16-mer peptide of known isoforms and isoallergens of these and other allergens.

Results

Interindividual differences in epitope recognition patterns were observed. In contrast, reactivity patterns in a given individual were by comparison more stable during the 3 years-course of the study. Nevertheless, evidence of the induction of novel specificities over time was identified across multiple regions of the allergens. Particularly reactive epitopes were identified in the D helix of Cyp c 1 and in the C-terminus of Gad c 1 and Gad m 1.02. Residues important for the recognition of certain linear epitopes were identified. Patterns of differential recognition of isoallergens were observed in some subjects.

Conclusions

Altogether, comprehensive analysis of antibody recognition of linear epitopes of multiple allergens enables characterisation of the nature of the antibody responses targeting this important set of food allergens.

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.

The Challenges of Vaccine Trial Participation among Underserved and Hard-to-Reach Communities: An Internal Expert Consultation of the VACCELERATE Consortium

Underserved and hard-to-reach population groups are under-represented in vaccine trials. Thus, we aimed to identify the challenges of vaccine trial participation of these groups in member countries of the VACCELERATE network. Seventeen National Coordinators (NC), each representing their respective country (15 European countries, Israel, and Turkey), completed an online survey. From 15 eligible groups, those that were more frequently declared underserved/hard-to-reach in vaccine research were ethnic minorities (76.5%), persons experiencing homelessness (70.6%), illegal workers and refugees (64.7%, each). When prioritization for education on vaccine trials was considered, ethnic groups, migrants, and immigrants (5/17, 29.4%) were the groups most frequently identified by the NC as top targets. The most prominent barriers in vaccine trial participation affecting all groups were low levels of health literacy, reluctance to participate in trials due to engagement level, and low levels of trust in vaccines/vaccinations. This study highlighted population groups considered underserved/hard-to-reach in countries contained within the European region, and the respective barriers these groups face when participating in clinical studies. Our findings aid with the design of tailored interventions (within-and across-countries of the European region) and with the development of strategies to overcome major barriers in phase 2 and phase 3 vaccine trial participation.

Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), has been shown to be of critical importance for immune responses to pathogens and vaccines. In recent years, B cell and T cell receptor repertoire sequencing (Rep-Seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci in different populations. Here we present a novel algorithm for extra-sensitive and specific variable (V) and joining (J) gene allele inference and genotyping allowing reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput genotyping and novel allele discovery from a wide variety of existing datasets. The developed algorithm is a part of the MiXCR software ( https://mixcr.com ) and can be incorporated into any pipeline utilizing upstream processing with MiXCR. We demonstrate the accuracy of this approach using Rep-Seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) Rep-Seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA; TRB) Rep-Seq dataset, representing 134 individuals. This allowed us to assess the genetic diversity of genes within the IGH, TRA and TRB loci in different populations and demonstrate the connection between antibody repertoire gene usage and the number of allelic variants present in the population. Finally we established a database of allelic variants of V and J genes inferred from Rep-Seq data and their population frequencies with free public access at https://vdj.online .

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).

B cells and the intestinal microbiome in time, space and place

The gut immune system is shaped by the continuous interaction with the microbiota. Here we dissect temporal, spatial and contextual layers of gut B cell responses. The microbiota impacts on the selection of the developing pool of pre-immune B cells that serves as substrate for B cell activation, expansion and differentiation. However, various aspects of the gut B cell response display unique features. In particular, occurrence of somatically mutated B cells, chronic gut germinal centers in T cell-deficient settings and polyreactive binding of gut IgA to the microbiota questioned the nature and microbiota-specificity of gut germinal centers. We propose a model to reconcile these observations incorporating recent work demonstrating microbiota-specificity of gut germinal centers. We speculate that adjuvant effects of the microbiota might modify permissiveness for B cell to enter and exit gut germinal centers. We propose that separating aspects of time, space and place facilitate the occasionally puzzling discussion of gut B cell responses to the microbiota.

Antibody gene features associated with binding and functional activity in vaccine-derived human mAbs targeting malaria parasites

Adjuvants have been essential to malaria vaccine development, but their impact on the vaccine-induced antibody repertoire is poorly understood. Here, we used cDNA sequences from antigen-specific single memory B cells to express 132 recombinant human anti-Pfs230 monoclonal antibodies (mAbs). Alhydrogel®-induced mAbs demonstrated higher binding to Pfs230D1, although functional activity was similar between adjuvants. All Alhydrogel® mAbs using IGHV1-69 gene bound to recombinant Pfs230D1, but none blocked parasite transmission to mosquitoes; similarly, no AS01 mAb using IGHV1-69 blocked transmission. Functional mAbs from both Alhydrogel® and AS01 vaccines used IGHV3-21 and IGHV3-30 genes. Antibodies with the longest CDR3 sequences were associated with binding but not functional activity. This study assesses adjuvant effects on antibody clonotype diversity during malaria vaccination.

Defining and Studying B Cell Receptor and TCR Interactions

BCRs (Abs) and TCRs (or adaptive immune receptors [AIRs]) are the means by which the adaptive immune system recognizes foreign and self-antigens, playing an integral part in host defense, as well as the emergence of autoimmunity. Importantly, the interaction between AIRs and their cognate Ags defies a simple key-in-lock paradigm and is instead a complex many-to-many mapping between an individual's massively diverse AIR repertoire, and a similarly diverse antigenic space. Understanding how adaptive immunity balances specificity with epitopic coverage is a key challenge for the field, and terms such as broad specificity, cross-reactivity, and polyreactivity remain ill-defined and are used inconsistently. In this Immunology Notes and Resources article, a group of experimental, structural, and computational immunologists define commonly used terms associated with AIR binding, describe methodologies to study these binding modes, as well as highlight the implications of these different binding modes for therapeutic design.

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Widespread impact of immunoglobulin V gene allelic polymorphisms on antibody reactivity

The ability of human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V gene (IGV) allelic polymorphisms. However, previous studies have provided only a limited number of examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many IGV allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiment further demonstrates that paratope allelic mutations on both heavy and light chain often abolish antibody binding. We also illustrate the importance of minor IGV allelic variants with low frequency in several broadly neutralizing antibodies to SARS-CoV-2 and influenza virus. Overall, this study not only highlights the pervasive impact of IGV allelic polymorphisms on antibody binding, but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

An influenza H1 hemagglutinin stem-only immunogen elicits a broadly cross-reactive B cell response in humans

Current yearly seasonal influenza vaccines primarily induce an antibody response directed against the immunodominant but continually diversifying hemagglutinin (HA) head region. These antibody responses provide protection against the vaccinating strain but little cross-protection against other influenza strains or subtypes. To focus the immune response on subdominant but more conserved epitopes on the HA stem that might protect against a broad range of influenza strains, we developed a stabilized H1 stem immunogen lacking the immunodominant head displayed on a ferritin nanoparticle (H1ssF). Here, we evaluated the B cell response to H1ssF in healthy adults ages 18 to 70 in a phase 1 clinical trial (NCT03814720). We observed both a strong plasmablast response and sustained elicitation of cross-reactive HA stem-specific memory B cells after vaccination with H1ssF in individuals of all ages. The B cell response was focused on two conserved epitopes on the H1 stem, with a highly restricted immunoglobulin repertoire unique to each epitope. On average, two-thirds of the B cell and serological antibody response recognized a central epitope on the H1 stem and exhibited broad neutralization across group 1 influenza virus subtypes. The remaining third recognized an epitope near the viral membrane anchor and was largely limited to H1 strains. Together, we demonstrate that an H1 HA immunogen lacking the immunodominant HA head produces a robust and broadly neutralizing HA stem-directed B cell response.

Application of germline antibody features to vaccine development, antibody discovery, antibody optimization and disease diagnosis

Although the efficacy and commercial success of vaccines and therapeutic antibodies have been tremendous, designing and discovering new drug candidates remains a labor-, time- and cost-intensive endeavor with high risks. The main challenges of vaccine development are inducing a strong immune response in broad populations and providing effective prevention against a group of highly variable pathogens. Meanwhile, antibody discovery faces several great obstacles, especially the blindness in antibody screening and the unpredictability of the developability and druggability of antibody drugs. These challenges are largely due to poorly understanding of germline antibodies and the antibody responses to pathogen invasions. Thanks to the recent developments in high-throughput sequencing and structural biology, we have gained insight into the germline immunoglobulin (Ig) genes and germline antibodies and then the germline antibody features associated with antigens and disease manifestation. In this review, we firstly outline the broad associations between germline antibodies and antigens. Moreover, we comprehensively review the recent applications of antigen-specific germline antibody features, physicochemical properties-associated germline antibody features, and disease manifestation-associated germline antibody features on vaccine development, antibody discovery, antibody optimization, and disease diagnosis. Lastly, we discuss the bottlenecks and perspectives of current and potential applications of germline antibody features in the biotechnology field.

Vaccinomics: A scoping review

BACKGROUND

This scoping review summarizes a key aspect of vaccinomics by collating known associations between heterogeneity in human genetics and vaccine immunogenicity and safety.

METHODS

We searched PubMed for articles in English using terms covering vaccines routinely recommended to the general US population, their effects, and genetics/genomics. Included studies were controlled and demonstrated statistically significant associations with vaccine immunogenicity or safety. Studies of Pandemrix®, an influenza vaccine previously used in Europe, were also included, due to its widely publicized genetically mediated association with narcolepsy.

FINDINGS

Of the 2,300 articles manually screened, 214 were included for data extraction. Six included articles examined genetic influences on vaccine safety; the rest examined vaccine immunogenicity. Hepatitis B vaccine immunogenicity was reported in 92 articles and associated with 277 genetic determinants across 117 genes. Thirty-three articles identified 291 genetic determinants across 118 genes associated with measles vaccine immunogenicity, 22 articles identified 311 genetic determinants across 110 genes associated with rubella vaccine immunogenicity, and 25 articles identified 48 genetic determinants across 34 genes associated with influenza vaccine immunogenicity. Other vaccines had fewer than 10 studies each identifying genetic determinants of their immunogenicity. Genetic associations were reported with 4 adverse events following influenza vaccination (narcolepsy, GBS, GCA/PMR, high temperature) and 2 adverse events following measles vaccination (fever, febrile seizure).

CONCLUSION

This scoping review identified numerous genetic associations with vaccine immunogenicity and several genetic associations with vaccine safety. Most associations were only reported in one study. This illustrates both the potential of and need for investment in vaccinomics. Current research in this field is focused on systems and genetic-based studies designed to identify risk signatures for serious vaccine reactions or diminished vaccine immunogenicity. Such research could bolster our ability to develop safer and more effective vaccines.

Human immunoglobulin gene allelic variation impacts germline-targeting vaccine priming

Vaccine priming immunogens that activate germline precursors for broadly neutralizing antibodies (bnAbs) have promise for development of precision vaccines against major human pathogens. In a clinical trial of the eOD-GT8 60mer germline-targeting immunogen, higher frequencies of vaccine-induced VRC01-class bnAb-precursor B cells were observed in the high dose compared to the low dose group. Through immunoglobulin heavy chain variable (IGHV) genotyping, statistical modeling, quantification of IGHV1-2 allele usage and B cell frequencies in the naive repertoire for each trial participant, and antibody affinity analyses, we found that the difference between dose groups in VRC01-class response frequency was best explained by IGHV1-2 genotype rather than dose and was most likely due to differences in IGHV1-2 B cell frequencies for different genotypes. The results demonstrate the need to define population-level immunoglobulin allelic variations when designing germline-targeting immunogens and evaluating them in clinical trials.

One-Sentence Summary

Human genetic variation can modulate the strength of vaccine-induced broadly neutralizing antibody precursor B cell responses.

Impact of host genetic polymorphisms on response to inactivated influenza vaccine in children

In randomized controlled trials of influenza vaccination, 550 children received trivalent-inactivated influenza vaccine, permitting us to explore relationship between vaccine response and host single nucleotide polymorphisms (SNPs) in 23 candidate genes with adjustment of multiple testing. For host SNPs in TLR7-1817G/T (rs5741880), genotype GT was associated with lower odds (OR: 0.22, 95% CI: 0.09, 0.53) of have post-vaccination hemagglutination-inhibiting (HAI) titers ≥40, compared with genotype GG and TT combined under the over-dominant model. For host SNPs in TLR8-129G/C (rs3764879), genotype GT was associated with lower odds (OR: 0.47; 95% CI: 0.28, 0.80) of have post vaccination HAI titers ≥40 compared with genotype GG and AA combined under the over-dominant model. Our results could contribute to the development of better vaccines that may offer improved protection to all recipients.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

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.

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.

Persistence of immunogenicity, contributing factors of an immune response, and reactogenicities after a single dose of the ChAdOx1 (AZD1222) COVID-19 vaccine in the Thai population

INTRODUCTION

Due to the vaccine's short supply and the efficacy of a single dose of the ChAdOx1 (AZD1222) vaccine, many governments delayed the interval between prime and boost dose from 4 to 8-12 weeks. However, the waning of immune response in this period is a concern. This study evaluated the durability, contributing factors of anti-RBD antibody concentration, and reactogenicities after the single dose of AZD1222 vaccine in the Thai population.

METHODS

This was a single-center, prospective cohort study at Chulabhorn Hospital, Bangkok, Thailand. Individuals 18 years or older who were negative for anti-SARS-CoV-2 antibody were eligible. Anti- receptor-binding domain antibody concentrations were tested at least three weeks after the first vaccination and immediately before the second dose of vaccine. Information on reactogenicities was obtained via a questionnaire sent by a short message service.

RESULTS

Anti-RBD Antibody concentration at 2 and 3 months post-vaccination were significantly higher than at 1 months post-vaccination (20.14 BAU/mL (95%CI; 16.37, 24.77) at 1 month, 48.08 BAU/mL (95%CI; 42.76, 54.08) at 2 month, and 65.01 BAU/mL (95%CI; 58.88,71.61) at 3 month). Adverse events occurred in approximately 60% of participants. Factors influencing vaccine immunogenicity include age, sex, the time elapsed from the first dose of vaccine, and underlying disease with diabetes and hematologic disease.

CONCLUSION

A single dose of AZD1222 could elicit immune responses that did not decline within three months in Thai individuals. These data support the public health strategy of a delay between the prime and boost dose of AZD1222 of 4 to 12 weeks.

Co-immunization with hemagglutinin stem immunogens elicits cross-group neutralizing antibodies and broad protection against influenza A viruses

Current influenza vaccines predominantly induce immunity to the hypervariable hemagglutinin (HA) head, requiring frequent vaccine reformulation. Conversely, the immunosubdominant yet conserved HA stem harbors a supersite that is targeted by broadly neutralizing antibodies (bnAbs), representing a prime target for universal vaccines. Here, we showed that the co-immunization of two HA stem immunogens derived from group 1 and 2 influenza A viruses elicits cross-group protective immunity and neutralizing antibody responses in mice, ferrets, and nonhuman primates (NHPs). Immunized mice were protected from multiple group 1 and 2 viruses, and all animal models showed broad serum-neutralizing activity. A bnAb isolated from an immunized NHP broadly neutralized and protected against diverse viruses, including H5N1 and H7N9. Genetic and structural analyses revealed strong homology between macaque and human bnAbs, illustrating common biophysical constraints for acquiring cross-group specificity. Vaccine elicitation of stem-directed cross-group-protective immunity represents a step toward the development of broadly protective influenza vaccines.

Transcriptome signatures preceding the induction of anti-stalk antibodies elicited after universal influenza vaccination

A phase 1 clinical trial to test the immunogenicity of a chimeric group 1 HA (cHA) universal influenza virus vaccine targeting the conserved stalk domain of the hemagglutinin of influenza viruses was carried out. Vaccination with adjuvanted-inactivated vaccines induced high anti-stalk antibody titers. We sought to identify gene expression signatures that correlate with such induction. Messenger-RNA sequencing in whole blood was performed on the peripheral blood of 53 vaccinees. We generated longitudinal data on the peripheral blood of 53 volunteers, at early (days 3 and 7) and late (28 days) time points after priming and boosting with cHAs. Differentially expressed gene analysis showed no differences between placebo and live-attenuated vaccine groups. However, an upregulation of genes involved in innate immune responses and type I interferon signaling was found at day 3 after vaccination with inactivated adjuvanted formulations. Cell type deconvolution analysis revealed a significant enrichment for monocyte markers and different subsets of dendritic cells as mediators for optimal B cell responses and significant increase of anti-stalk antibodies in sera. A significant upregulation of immunoglobulin-related genes was only observed after administration of adjuvanted vaccines (either as primer or booster) with specific induction of anti-stalk IGVH1-69. This approach informed of specific immune signatures that correlate with robust anti-stalk antibody responses, while also helping to understand the regulation of gene expression induced by cHA proteins under different vaccine regimens.

Pan-vaccine analysis reveals innate immune endotypes predictive of antibody responses to vaccination

Several studies have shown that the pre-vaccination immune state is associated with the antibody response to vaccination. However, the generalizability and mechanisms that underlie this association remain poorly defined. Here, we sought to identify a common pre-vaccination signature and mechanisms that could predict the immune response across 13 different vaccines. Analysis of blood transcriptional profiles across studies revealed three distinct pre-vaccination endotypes, characterized by the differential expression of genes associated with a pro-inflammatory response, cell proliferation, and metabolism alterations. Importantly, individuals whose pre-vaccination endotype was enriched in pro-inflammatory response genes known to be downstream of nuclear factor-kappa B showed significantly higher serum antibody responses 1 month after vaccination. This pro-inflammatory pre-vaccination endotype showed gene expression characteristic of the innate activation state triggered by Toll-like receptor ligands or adjuvants. These results demonstrate that wide variations in the transcriptional state of the immune system in humans can be a key determinant of responsiveness to vaccination.

Molecular analysis of a public cross-neutralizing antibody response to SARS-CoV-2

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concerns (VOCs) continue to emerge, cross-neutralizing antibody responses become key toward next-generation design of a more universal COVID-19 vaccine. By analyzing published data from the literature, we report here that the combination of germline genes IGHV2-5/IGLV2-14 represents a public antibody response to the receptor-binding domain (RBD) that potently cross-neutralizes a broad range of VOCs, including Omicron and its sub-lineages. Detailed molecular analysis shows that the complementarity-determining region H3 sequences of IGHV2-5/IGLV2-14-encoded RBD antibodies have a preferred length of 11 amino acids and a conserved HxIxxI motif. In addition, these antibodies have a strong allelic preference due to an allelic polymorphism at amino acid residue 54 of IGHV2-5, which is located at the paratope. These findings have important implications for understanding cross-neutralizing antibody responses to SARS-CoV-2 and its heterogenicity at the population level as well as the development of a universal COVID-19 vaccine.

ImmunoTyper-SR: A computational approach for genotyping immunoglobulin heavy chain variable genes using short-read data

Human immunoglobulin heavy chain (IGH) locus on chromosome 14 includes more than 40 functional copies of the variable gene (IGHV), which are critical for the structure of antibodies that identify and neutralize pathogenic invaders as a part of the adaptive immune system. Because of its highly repetitive sequence composition, the IGH locus has been particularly difficult to assemble or genotype when using standard short-read sequencing technologies. Here, we introduce ImmunoTyper-SR, an algorithmic tool for the genotyping and CNV analysis of the germline IGHV genes on Illumina whole-genome sequencing (WGS) data using a combinatorial optimization formulation that resolves ambiguous read mappings. We have validated ImmunoTyper-SR on 12 individuals, whose IGHV allele composition had been independently validated, as well as concordance between WGS replicates from nine individuals. We then applied ImmunoTyper-SR on 585 COVID patients to investigate the associations between IGHV alleles and anti-type I IFN autoantibodies, which were previously associated with COVID-19 severity.

Allelic polymorphism controls autoreactivity and vaccine elicitation of human broadly neutralizing antibodies against influenza virus

Human broadly neutralizing antibodies (bnAbs) targeting the hemagglutinin stalk of group 1 influenza A viruses (IAVs) are biased for IGHV1-69 alleles that use phenylalanine (F54) but not leucine (L54) within their CDRH2 loops. Despite this, we demonstrated that both alleles encode for human IAV bnAbs that employ structurally convergent modes of contact to the same epitope. To resolve differences in lineage expandability, we compared F54 versus L54 as substrate within humanized mice, where antibodies develop with human-like CDRH3 diversity but are restricted to single VH genes. While both alleles encoded for bnAb precursors, only F54 IGHV1-69 supported elicitation of heterosubtypic serum bnAbs following immunization with a stalk-only nanoparticle vaccine. L54 IGHV1-69 was unproductive, co-encoding for anergic B cells and autoreactive stalk antibodies that were cleared from B cell memory. Moreover, human stalk antibodies also demonstrated L54-dependent autoreactivity. Therefore, IGHV1-69 polymorphism, which is skewed ethnically, gates tolerance and vaccine expandability of influenza bnAbs.

Profiling genes encoding the adaptive immune receptor repertoire with gAIRR Suite

Adaptive immune receptor repertoire (AIRR) is encoded by T cell receptor (TR) and immunoglobulin (IG) genes. Profiling these germline genes encoding AIRR (abbreviated as gAIRR) is important in understanding adaptive immune responses but is challenging due to the high genetic complexity. Our gAIRR Suite comprises three modules. gAIRR-seq, a probe capture-based targeted sequencing pipeline, profiles gAIRR from individual DNA samples. gAIRR-call and gAIRR-annotate call alleles from gAIRR-seq reads and annotate whole-genome assemblies, respectively. We gAIRR-seqed TRV and TRJ of seven Genome in a Bottle (GIAB) DNA samples with 100% accuracy and discovered novel alleles. We also gAIRR-seqed and gAIRR-called the TR and IG genes of a subject from both the peripheral blood mononuclear cells (PBMC) and oral mucosal cells. The calling results from these two cell types have a high concordance (99% for all known gAIRR alleles). We gAIRR-annotated 36 genomes to unearth 325 novel TRV alleles and 29 novel TRJ alleles. We could further profile the flanking sequences, including the recombination signal sequence (RSS). We validated two structural variants for HG002 and uncovered substantial differences of gAIRR genes in references GRCh37 and GRCh38. gAIRR Suite serves as a resource to sequence, analyze, and validate germline TR and IG genes to study various immune-related phenotypes.

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.

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.

Molecular analysis of a public cross-neutralizing antibody response to SARS-CoV-2

As SARS-CoV-2 variants of concerns (VOCs) continue to emerge, cross-neutralizing antibody responses become key towards next-generation design of a more universal COVID-19 vaccine. By analyzing published data from the literature, we report here that the combination of germline genes IGHV2-5/IGLV2-14 represents a public antibody response to the receptor-binding domain (RBD) that potently cross-neutralizes all VOCs to date, including Omicron and its sub-lineages. Detailed molecular analysis shows that the complementarity-determining region H3 sequences of IGHV2-5/IGLV2-14-encoded RBD antibodies have a preferred length of 11 amino acids and a conserved HxIxxI motif. In addition, these antibodies have a strong allelic preference due to an allelic polymorphism at amino-acid residue 54 of IGHV2-5, which locates at the paratope. These findings have important implications for understanding cross-neutralizing antibody responses to SARS-CoV-2 and its heterogenicity at the population level as well as the development of a universal COVID-19 vaccine.

Pandemic, Epidemic, Endemic: B Cell Repertoire Analysis Reveals Unique Anti-Viral Responses to SARS-CoV-2, Ebola and Respiratory Syncytial Virus

Immunoglobulin gene heterogeneity reflects the diversity and focus of the humoral immune response towards different infections, enabling inference of B cell development processes. Detailed compositional and lineage analysis of long read IGH repertoire sequencing, combining examples of pandemic, epidemic and endemic viral infections with control and vaccination samples, demonstrates general responses including increased use of IGHV4-39 in both Zaire Ebolavirus (EBOV) and COVID-19 patient cohorts. We also show unique characteristics absent in Respiratory Syncytial Virus or yellow fever vaccine samples: EBOV survivors show unprecedented high levels of class switching events while COVID-19 repertoires from acute disease appear underdeveloped. Despite the high levels of clonal expansion in COVID-19 IgG1 repertoires there is a striking lack of evidence of germinal centre mutation and selection. Given the differences in COVID-19 morbidity and mortality with age, it is also pertinent that we find significant differences in repertoire characteristics between young and old patients. Our data supports the hypothesis that a primary viral challenge can result in a strong but immature humoral response where failures in selection of the repertoire risk off-target effects.

HLA Homozygosity and Likelihood of Sensitization in Kidney Transplant Candidates

Background.

Homozygosity for HLAs has been associated with adverse outcomes after viral infection as well as pregnancy-induced HLA sensitization. We sought to assess the relationship between HLA locus homozygosity and the level of HLA antibody sensitization.

Methods.

We measured sensitization using the calculated panel reactive antibody value for a large cohort of 147 461 patients added to the US OPTN/United Network for Organ Sharing kidney transplant waitlist between December 2014 and December 2019. We used multinomial logistic modeling to compare 62 510 sensitized patients to 84 955 unsensitized controls.

Results.

We found that the number of homozygous HLA loci was strongly associated with the level of sensitization. Within mildly, highly, or extremely sensitized candidates, women displayed a higher relative abundance of HLA homozygosity at multiple HLA loci as compared with men, with attenuation of this effect in Black candidates. In a multivariable logistic model, the number of homozygous HLA loci interacted with female sex but not with other factors associated with sensitization, including recipient ethnicity and a history of prior kidney transplant.

Conclusions.

This study shows that HLA homozygosity is an innate genetic factor that affects the likelihood of HLA sensitization. Further research is needed to identify the immunologic mechanisms that underlie this observation.

Variations in antibody repertoires correlate with vaccine responses

An important challenge in vaccine development is to figure out why a vaccine succeeds in some individuals and fails in others. Although antibody repertoires hold the key to answering this question, there have been very few personalized immunogenomics studies so far aimed at revealing how variations in immunoglobulin genes affect a vaccine response. We conducted an immunosequencing study of 204 calves vaccinated against bovine respiratory disease (BRD) with the goal to reveal variations in immunoglobulin genes and somatic hypermutations that impact the efficacy of vaccine response. Our study represents the largest longitudinal personalized immunogenomics study reported to date across all species, including humans. To analyze the generated data set, we developed an algorithm for identifying variations of the immunoglobulin genes (as well as frequent somatic hypermutations) that affect various features of the antibody repertoire and titers of neutralizing antibodies. In contrast to relatively short human antibodies, cattle have a large fraction of ultralong antibodies that have opened new therapeutic opportunities. Our study reveals that ultralong antibodies are a key component of the immune response against the costliest disease of beef cattle in North America. The detected variants of the cattle immunoglobulin genes, which are implicated in the success/failure of the BRD vaccine, have the potential to direct the selection of individual cattle for ongoing breeding programs.

Conformational plasticity of the HIV-1 gp41 immunodominant region is recognized by multiple non-neutralizing antibodies

The early humoral immune response to acute HIV-1 infection is largely non-neutralizing. The principal target of these antibodies is the primary immunodominant region (PID) on the gp41 fusion protein. The PID is a highly conserved 15-residue region displayed on the surface of HIV-1 virions. In this study, we analyzed the humoral determinants of HIV-1 gp41 PID binding using biophysical, structural, and computational methods. In complex with a patient-derived near-germline antibody fragment, the PID motif adopts an elongated random coil, whereas the PID bound to affinity-matured Fab adopts a strand-turn-helix conformation. Molecular dynamics simulations showed that the PID is structurally plastic suggesting that the PID can form an ensemble of structural states recognized by various non-neutralizing antibodies, facilitating HIV-1 immunodominance observed in acute and chronic HIV-1 infections. An improved understanding of how the HIV-1 gp41 PID misdirects the early humoral response should guide the development of an effective HIV-1 vaccine.

The 15-amino-acid primary immunodominant (PID) region on HIV-1 gp41 adopts an ensemble of conformational states. This conformational plasticity is suggested to misdirect the early humoral immune response.

A single residue in influenza virus H2 hemagglutinin enhances the breadth of the B cell response elicited by H2 vaccination

Conserved epitopes on the influenza hemagglutinin (HA) stem are an attractive target for universal vaccine strategies as they elicit broadly neutralizing antibodies. Such antibody responses to stem-specific epitopes have been extensively characterized for HA subtypes H1 and H5 in humans. H2N2 influenza virus circulated 50 years ago and represents a pandemic threat due to the lack of widespread immunity, but, unlike H1 and H5, the H2 HA stem contains Phe45_HA2 predicted to sterically clash with HA stem-binding antibodies characterized to date. To understand the effect of Phe45_HA2, we compared the HA stem-specific B cell response in post hoc analyses of two phase 1 clinical trials, one testing vaccination with an H2 ferritin nanoparticle immunogen ( NCT03186781 ) and one with an inactivated H5N1 vaccine ( NCT01086657 ). In H2-naive individuals, the magnitude of the B cell response was equivalent, but H2-elicited HA stem-binding B cells displayed greater cross-reactivity than those elicited by H5. However, in individuals with childhood H2 exposure, H5-elicited HA stem-binding B cells also displayed high cross-reactivity, suggesting recall of memory B cells formed 50 years ago. Overall, we propose that a one-residue difference on an HA immunogen can alter establishment and expansion of broadly neutralizing memory B cells. These data have implications for stem-based universal influenza vaccination strategies.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

A family matter: Anti-viral antibody responses

The immune response is an amalgam of heritable and non-heritable effects, although the contribution of each is generally unknown. In this issue of Immunity, Venkataraman et al. report heritable traits that govern anti-viral antibody epitope specificity and breadth controlled by the MHC class II region, indicating that your genes influence which pathogen-derived epitopes your B cells target during viral infections.

VDJ Gene Usage in IgM Repertoires of Rhesus and Cynomolgus Macaques

Macaques are frequently used to evaluate candidate vaccines and to study infection-induced antibody responses, requiring an improved understanding of their naïve immunoglobulin (IG) repertoires. Baseline gene usage frequencies contextualize studies of antigen-specific immune responses, providing information about how easily one may stimulate a response with a particular VDJ recombination. Studies of human IgM repertoires have shown that IG VDJ gene frequencies vary several orders of magnitude between the most and least utilized genes in a manner that is consistent across many individuals but to date similar analyses are lacking for macaque IgM repertoires. Here, we quantified VDJ gene usage levels in unmutated IgM repertoires of 45 macaques, belonging to two species and four commonly used subgroups: Indian and Chinese origin rhesus macaques and Indonesian and Mauritian origin cynomolgus macaques. We show that VDJ gene frequencies differed greatly between the most and least used genes, with similar overall patterns observed in macaque subgroups and individuals. However, there were also clear differences affecting the use of specific V, D and J genes. Furthermore, in contrast to humans, macaques of both species utilized IGHV4 family genes to a much higher extent and showed evidence of evolutionary expansion of genes of this family. Finally, we used the results to inform the analysis of a broadly neutralizing HIV-1 antibody elicited in SHIV-infected rhesus macaques, RHA1.V2.01, which binds the apex of the Env trimer in a manner that mimics the binding mode of PGT145. We discuss the likelihood that similar antibodies could be elicited in different macaque subgroups.

Individualized VDJ recombination predisposes the available Ig sequence space

The process of recombination between variable (V), diversity (D), and joining (J) immunoglobulin (Ig) gene segments determines an individual's naive Ig repertoire and, consequently, (auto)antigen recognition. VDJ recombination follows probabilistic rules that can be modeled statistically. So far, it remains unknown whether VDJ recombination rules differ between individuals. If these rules differed, identical (auto)antigen-specific Ig sequences would be generated with individual-specific probabilities, signifying that the available Ig sequence space is individual specific. We devised a sensitivity-tested distance measure that enables inter-individual comparison of VDJ recombination models. We discovered, accounting for several sources of noise as well as allelic variation in Ig sequencing data, that not only unrelated individuals but also human monozygotic twins and even inbred mice possess statistically distinguishable immunoglobulin recombination models. This suggests that, in addition to genetic, there is also nongenetic modulation of VDJ recombination. We demonstrate that population-wide individualized VDJ recombination can result in orders of magnitude of difference in the probability to generate (auto)antigen-specific Ig sequences. Our findings have implications for immune receptor-based individualized medicine approaches relevant to vaccination, infection, and autoimmunity.

Novel Allele Detection Tool Benchmark and Application With Antibody Repertoire Sequencing Dataset

Detailed knowledge of the diverse immunoglobulin germline genes is critical for the study of humoral immunity. Hundreds of alleles have been discovered by analyzing antibody repertoire sequencing (Rep-seq or Ig-seq) data via multiple novel allele detection tools (NADTs). However, the performance of these NADTs through antibody sequences with intrinsic somatic hypermutations (SHMs) is unclear. Here, we developed a tool to simulate repertoires by integrating the full spectrum features of an antibody repertoire such as germline gene usage, junctional modification, position-specific SHM and clonal expansion based on 2152 high-quality datasets. We then systematically evaluated these NADTs using both simulated and genuine Ig-seq datasets. Finally, we applied these NADTs to 687 Ig-seq datasets and identified 43 novel allele candidates (NACs) using defined criteria. Twenty-five alleles were validated through findings of other sources. In addition to the NACs detected, our simulation tool, the results of our comparison, and the streamline of this process may benefit further humoral immunity studies via Ig-seq.

Germline polymorphisms and alternative splicing of human immunoglobulin light chain genes

Inference of germline polymorphisms in immunoglobulin genes from B cell receptor repertoires is complicated by somatic hypermutations, sequencing/PCR errors, and by varying length of reference alleles. The light chain inference is particularly challenging owing to large gene duplications and absence of D genes. We analyzed the light chain cDNA sequences from naïve B cell receptor repertoires from 100 individuals. We optimized light chain allele inference by tweaking parameters of the TIgGER functions, extending the germline reference sequences, and establishing mismatch frequency patterns at polymorphic positions to filter out false-positive candidates. We identified 48 previously unreported variants of light chain variable genes. We selected 14 variants for validation and successfully validated 11 by Sanger sequencing. Clustering of light chain 5'UTR, L-PART1, and L-PART2 revealed partial intron retention in 11 kappa and 9 lambda V alleles. Our results provide insight into germline variation in human light chain immunoglobulin loci.

Age-related changes in the TRB and IGH repertoires in healthy adult males and females

A diverse immune repertoire is capable of recognizing the enormous universe of foreign antigens encountered over life. Aging has a profound impact on the immune repertoires. However, whether continuous age-related changes in the immune repertoires differ between sexes is unclear. In this study, the characteristics of immune repertoires stratified by sex during aging are deciphered by analyzing T-cell receptor β-chain (TRB) and immunoglobulin heavy chain (IGH) sequences in 361 healthy adults. A similar change was observed between males and females across their lifespan, whereas age-subgroup analysis revealed sex-specific signatures in TRB and IGH repertoires. In regard to TRB, in males, repertoire richness and evenness increases slightly before the age of 32 years and 45 years respectively, and decreases sharply thereafter. Intriguingly, in females, they decrease significantly until around the age 57 years old, and subsequently undergo a stable stage up to the age of 83 years. Although IGH repertoire evenness increases significantly with age in both sexes, richness decreases significantly with age in males but remains stable in females. Moreover, average length of IGH CDR3 increases with age. In conclusion, these findings provide fundamental insights into the mechanisms underlying sex differences in adaptive immunity.

Computational Inference, Validation, and Analysis of 5'UTR-Leader Sequences of Alleles of Immunoglobulin Heavy Chain Variable Genes

Upstream and downstream sequences of immunoglobulin genes may affect the expression of such genes. However, these sequences are rarely studied or characterized in most studies of immunoglobulin repertoires. Inference from large, rearranged immunoglobulin transcriptome data sets offers an opportunity to define the upstream regions (5'-untranslated regions and leader sequences). We have now established a new data pre-processing procedure to eliminate artifacts caused by a 5'-RACE library generation process, reanalyzed a previously studied data set defining human immunoglobulin heavy chain genes, and identified novel upstream regions, as well as previously identified upstream regions that may have been identified in error. Upstream sequences were also identified for a set of previously uncharacterized germline gene alleles. Several novel upstream region variants were validated, for instance by their segregation to a single haplotype in heterozygotic subjects. SNPs representing several sequence variants were identified from population data. Finally, based on the outcomes of the analysis, we define a set of testable hypotheses with respect to the placement of particular alleles in complex IGHV locus haplotypes, and discuss the evolutionary relatedness of particular heavy chain variable genes based on sequences of their upstream regions.

Vaccine genetics of IGHV1-2 VRC01-class broadly neutralizing antibody precursor naïve human B cells

A successful HIV vaccine eliciting broadly neutralizing antibodies (bnAbs) must overcome the hurdle of being able to activate naive precursor B cells encoding features within their germline B cell receptors (BCR) that allow recognition of broadly neutralizing epitopes. Knowledge of whether bnAb precursor B cells are circulating at sufficient frequencies within individuals in communities heavily impacted by HIV may be important. Using a germline-targeting eOD-GT8 immunogen and high-throughput droplet-based single-cell BCR sequencing, we demonstrate that large numbers of paired BCR sequences from multiple donors can be efficiently screened to elucidate precursor frequencies of rare, naive VRC01-class B cells. Further, we analyzed IGHV1-2 allelic usage among three different cohorts; we find that IGHV1-2 alleles traditionally thought to be incompatible with VRC01-class responses are relatively common in various human populations and that germline variation within IGHV1-2 associates with gene usage frequencies in the naive BCR repertoire.

Access to ultra-long IgG CDRH3 bovine antibody sequences using short read sequencing technology

The advances in high-throughput DNA sequencing and recombinant antibody technologies has presented new methods for characterizing antibody repertoires and significantly increased our understanding on the functional role of antibodies in immunity and their use in diagnostics, vaccine antigen design and as biological therapeutics. A subset of Bos taurus antibodies possesses unique ultra-long third complementary-determining region of the heavy chain (CDRH3) and are of special interest because they are thought to have unique functional abilities of broadly neutralizing properties - a functional role that has not been fully explored in vaccine development. Next generation sequencing technologies that are widely used to profile immunoglobulin (Ig) repertoires are based on short-read methods such as the Illumina technology. Although this technology has worked well in sequencing Ig V-D-J regions of most jawed vertebrates, it has faced serious technical challenges with sequencing regions in bovine Ig bearing ultra-long CDRH3 sequences, which are longer than 120 bp. To overcome this limitation, we have developed a sequencing strategy based on nested PCR products that allows sequence assembly of full-length bovine Ig heavy-chain (IgH) V-D-J regions. We have used this strategy to sequence IgH V-D-J regions of two Bos indicus breeds, Ankole and Boran. We confirm the presence of ultra-long CDRH3 sequences in IgG transcripts in both African cattle breeds, and provide preliminary evidence for differences and preferences in germline V_H, D_H and J_H allele gene usage as well as differences in the length of the VH region in the two bovine breeds. Our method provides tools that should allow more robust analyses of ultra-long CDRH3 sequences aiding antibody and epitope discovery in different cattle breeds and their role in mediating immunity.