Archaic humans have contributed to large-scale variation in modern human T cell receptor genes.

Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.

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.

GPR43 regulates marginal zone B-cell responses to foreign and endogenous antigens.

Marginal zone (MZ) B cells are innate-like B cells that produce polyreactive antibodies with an affinity for microbial molecular patterns and carbohydrate ligands. MZ B cells have been shown to be important in mediating immunity to various bacteria including Streptococcus pneumoniae and are also implicated in inflammatory syndromes including lupus erythematosus. The intestinal microbiota is responsible for producing short-chain fatty acids, which can regulate immune cell function by several mechanisms including ligation of the G-protein-coupled receptor (GPR)43. Herein, we show that MZ B cells express Gpr43 messenger RNA and that the absence of this receptor impacts on MZ B-cell surface marker expression and antibody production. In T-cell-independent responses to the hapten 4-hydroxy-3-nitrophenylacetic acid (NP), mice deficient in GPR43 displayed higher serum titers of NP-specific antibodies. Moreover, in response to a pneumococcal polysaccharide vaccine, GPR43-deficient mice developed robust serum antibody responses and had markedly increased numbers of splenic antibody-secreting cells, compared with control mice. Finally, serum immunoglobulin M autoantibodies to double-stranded DNA and phosphatidylcholine were increased in resting 10-15-week-old mice lacking GPR43. Taken together, mice lacking GPR43 have heightened antibody responses to T-cell-independent antigens, which may be a result of impaired regulation of MZ B cells.

Long-read amplicon denoising.

Long-read next-generation amplicon sequencing shows promise for studying complete genes or genomes from complex and diverse populations. Current long-read sequencing technologies have challenging error profiles, hindering data processing and incorporation into downstream analyses. Here we consider the problem of how to reconstruct, free of sequencing error, the true sequence variants and their associated frequencies from PacBio reads. Called 'amplicon denoising', this problem has been extensively studied for short-read sequencing technologies, but current solutions do not always successfully generalize to long reads with high indel error rates. We introduce two methods: one that runs nearly instantly and is very accurate for medium length reads and high template coverage, and another, slower method that is more robust when reads are very long or coverage is lower. On two Mock Virus Community datasets with ground truth, each sequenced on a different PacBio instrument, and on a number of simulated datasets, we compare our two approaches to each other and to existing algorithms. We outperform all tested methods in accuracy, with competitive run times even for our slower method, successfully discriminating templates that differ by a just single nucleotide. Julia implementations of Fast Amplicon Denoising (FAD) and Robust Amplicon Denoising (RAD), and a webserver interface, are freely available.

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.

Vaccination of SARS-CoV-2-infected individuals expands a broad range of clonally diverse affinity-matured B cell lineages.

Vaccination of SARS-CoV-2 convalescent individuals generates broad and potent antibody responses. Here, we isolate 459 spike-specific monoclonal antibodies (mAbs) from two individuals who were infected with the index variant of SARS-CoV-2 and later boosted with mRNA-1273. We characterize mAb genetic features by sequence assignments to the donors' personal immunoglobulin genotypes and assess antibody neutralizing activities against index SARS-CoV-2, Beta, Delta, and Omicron variants. The mAbs used a broad range of immunoglobulin heavy chain (IGH) V genes in the response to all sub-determinants of the spike examined, with similar characteristics observed in both donors. IGH repertoire sequencing and B cell lineage tracing at longitudinal time points reveals extensive evolution of SARS-CoV-2 spike-binding antibodies from acute infection until vaccination five months later. These results demonstrate that highly polyclonal repertoires of affinity-matured memory B cells are efficiently recalled by vaccination, providing a basis for the potent antibody responses observed in convalescent persons following vaccination.

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.