ABSTRACT
Long-read sequencing technologies have revolutionized genome assembly producing near-complete chromosome assemblies for numerous organisms, which are invaluable to research in many fields. However, regions with complex repetitive structure continue to represent a challenge for genome assembly algorithms, particularly in areas with high heterozygosity. Robust and comprehensive solutions for the assessment of assembly accuracy and completeness in these regions do not exist. In this study we focus on the assembly of biomedically important antibody-encoding immunoglobulin (IG) loci, which are characterized by complex duplications and repeat structures. High-quality full-length assemblies for these loci are critical for resolving haplotype-level annotations of IG genes, without which, functional and evolutionary studies of antibody immunity across vertebrates are not tractable. To address these challenges, we developed a pipeline, "CloseRead", that generates multiple assembly verification metrics for analysis and visualization. These metrics expand upon those of existing quality assessment tools and specifically target complex and highly heterozygous regions. Using CloseRead, we systematically assessed the accuracy and completeness of IG loci in publicly available assemblies of 74 vertebrate species, identifying problematic regions. We also demonstrated that inspecting assembly graphs for problematic regions can both identify the root cause of assembly errors and illuminate solutions for improving erroneous assemblies. For a subset of species, we were able to correct assembly errors through targeted reassembly. Together, our analysis demonstrated the utility of assembly assessment in improving the completeness and accuracy of IG loci across species.
ABSTRACT
We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.
ABSTRACT
New high-quality human genome assemblies derived from lymphoblastoid cell lines (LCLs) provide reference genomes and pangenomes for genomics studies. However, the characteristics of LCLs pose technical challenges to profiling immunoglobulin (IG) genes. IG loci in LCLs contain a mixture of germline and somatically recombined haplotypes, making them difficult to genotype or assemble accurately. To address these challenges, we introduce IGLoo, a software tool that implements novel methods for analyzing sequence data and genome assemblies derived from LCLs. IGLoo characterizes somatic V(D)J recombination events in the sequence data and identifies the breakpoints and missing IG genes in the LCL-based assemblies. Furthermore, IGLoo implements a novel reassembly framework to improve germline assembly quality by integrating information about somatic events and population structural variantions in the IG loci. We applied IGLoo to study the assemblies from the Human Pangenome Reference Consortium, providing new insights into the mechanisms, gene usage, and patterns of V(D)J recombination, causes of assembly fragmentation in the IG heavy chain (IGH) locus, and improved representation of the IGH assemblies.
ABSTRACT
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.
Subject(s)
Genes, Immunoglobulin Heavy Chain , Genes, Immunoglobulin , Humans , Genes, Immunoglobulin Heavy Chain/genetics , Alleles , Germ-Line Mutation , Immunoglobulin Heavy Chains/geneticsABSTRACT
Developing broad coronavirus vaccines requires identifying and understanding the molecular basis of broadly neutralizing antibody (bnAb) spike sites. In our previous work, we identified sarbecovirus spike RBD group 1 and 2 bnAbs. We have now shown that many of these bnAbs can still neutralize highly mutated SARS-CoV-2 variants, including the XBB.1.5. Structural studies revealed that group 1 bnAbs use recurrent germline-encoded CDRH3 features to interact with a conserved RBD region that overlaps with class 4 bnAb site. Group 2 bnAbs recognize a less well-characterized "site V" on the RBD and destabilize spike trimer. The site V has remained largely unchanged in SARS-CoV-2 variants and is highly conserved across diverse sarbecoviruses, making it a promising target for broad coronavirus vaccine development. Our findings suggest that targeted vaccine strategies may be needed to induce effective B cell responses to escape resistant subdominant spike RBD bnAb sites.
ABSTRACT
Pan-betacoronavirus neutralizing antibodies may hold the key to developing broadly protective vaccines against novel pandemic coronaviruses and to more effectively respond to SARS-CoV-2 variants. The emergence of Omicron and subvariants of SARS-CoV-2 illustrates the limitations of solely targeting the receptor-binding domain (RBD) of the spike (S) protein. Here, we isolated a large panel of broadly neutralizing antibodies (bnAbs) from SARS-CoV-2 recovered-vaccinated donors, which targets a conserved S2 region in the betacoronavirus spike fusion machinery. Select bnAbs showed broad in vivo protection against all three deadly betacoronaviruses, SARS-CoV-1, SARS-CoV-2, and MERS-CoV, which have spilled over into humans in the past two decades. Structural studies of these bnAbs delineated the molecular basis for their broad reactivity and revealed common antibody features targetable by broad vaccination strategies. These bnAbs provide new insights and opportunities for antibody-based interventions and for developing pan-betacoronavirus vaccines.
Subject(s)
COVID-19 , SARS-CoV-2 , Humans , Broadly Neutralizing Antibodies , Antibodies, Neutralizing , Antibodies, ViralABSTRACT
Ab "ultralong" third H chain complementarity-determining regions (CDR H3) appear unique to bovine Abs and may enable binding to difficult epitopes that shorter CDR H3 regions cannot easily access. Diversity is concentrated in the "knob" domain of the CDR H3, which is encoded by the DH gene segment and sits atop a ß-ribbon "stalk" that protrudes far from the Ab surface. Knob region cysteine content is quite diverse in terms of total number of cysteines, sequence position, and disulfide bond pattern formation. We investigated the role of germline cysteines in production of a diverse CDR H3 structural repertoire. The relationship between DH polymorphisms and deletions relative to germline at the nucleotide level, as well as diversity in cysteine and disulfide bond content at the structural level, was ascertained. Structural diversity is formed through (1) DH polymorphisms with altered cysteine positions, (2) DH deletions, and (3) new cysteines that arise through somatic hypermutation that form new, unique disulfide bonds to alter the knob structure. Thus, a combination of mechanisms at both the germline and somatic immunogenetic levels results in diversity in knob region cysteine content, contributing to remarkable complexity in knob region disulfide patterns, loops, and Ag binding surface.
Subject(s)
Cysteine , Germ Cells , Animals , Cattle , Cysteine/genetics , Polymorphism, Genetic , Complementarity Determining Regions/genetics , DisulfidesABSTRACT
To prepare for future coronavirus (CoV) pandemics, it is desirable to generate vaccines capable of eliciting broadly neutralizing antibody responses to CoVs. Here, we show that immunization of macaques with SARS-CoV-2 spike (S) protein with a two-shot protocol generated potent serum receptor binding domain cross-neutralizing antibody responses to both SARS-CoV-2 and SARS-CoV-1. Furthermore, responses were equally effective against most SARS-CoV-2 variants of concern (VOCs) and some were highly effective against Omicron. This result contrasts with human infection or many two-shot vaccination protocols where responses were typically more SARS-CoV-2 specific and where VOCs were less well neutralized. Structural studies showed that cloned macaque neutralizing antibodies, particularly using a given heavy chain germline gene, recognized a relatively conserved region proximal to the angiotensin converting enzyme 2 receptor binding site (RBS), whereas many frequently elicited human neutralizing antibodies targeted more variable epitopes overlapping the RBS. B cell repertoire differences between humans and macaques appeared to influence the vaccine response. The macaque neutralizing antibodies identified a pan-SARS-related virus epitope region less well targeted by human antibodies that could be exploited in rational vaccine design.
Subject(s)
COVID-19 , SARS-CoV-2 , Animals , Antibodies, Neutralizing , Antibodies, Viral , Broadly Neutralizing Antibodies , Epitopes , Humans , Macaca mulatta , Spike Glycoprotein, CoronavirusABSTRACT
The emergence of current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy. Development of broadly effective coronavirus vaccines that can mitigate these threats is needed. Here, we utilized a targeted donor selection strategy to isolate a large panel of human broadly neutralizing antibodies (bnAbs) to sarbecoviruses. Many of these bnAbs are remarkably effective in neutralizing a diversity of sarbecoviruses and against most SARS-CoV-2 VOCs, including the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor-binding domain (RBD). Consistent with targeting of conserved sites, select RBD bnAbs exhibited protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model in vivo. These bnAbs provide new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and provide a molecular basis for effective design of pan-sarbecovirus vaccines.
Subject(s)
COVID-19 , SARS-CoV-2 , Antibodies, Neutralizing , Antibodies, Viral , Broadly Neutralizing Antibodies , COVID-19/prevention & control , Humans , Spike Glycoprotein, CoronavirusABSTRACT
The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin (IG) loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) are critical for studies of the adaptive immune system. However, because the standard gene finding algorithms are not suitable for predicting IG genes, they have been semimanually 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, for example, 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, unlike all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Because 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.
Subject(s)
Antibody Diversity , V(D)J Recombination , Antibodies , Complementarity Determining Regions/genetics , Genes, ImmunoglobulinABSTRACT
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.
Subject(s)
Cattle Diseases , Vaccines , Animals , Antibodies , Cattle , Cattle Diseases/prevention & control , North America , Vaccines/geneticsABSTRACT
Pan-betacoronavirus neutralizing antibodies may hold the key to developing broadly protective vaccines against coronaviruses that cause severe disease, for anticipating novel pandemic-causing viruses, and to respond more effectively to SARS-CoV-2 variants. The emergence of the Omicron variant of SARS-CoV-2 has illustrated the limitations of solely targeting the receptor binding domain (RBD) of the envelope Spike (S)-protein. Here, we isolated a large panel of broadly neutralizing antibodies (bnAbs) from SARS-CoV-2 recovered-vaccinated donors that target a conserved S2 region in the fusion machinery on betacoronavirus spikes. Select bnAbs show broad in vivo protection against all three pathogenic betacoronaviruses, SARS-CoV-1, SARS-CoV-2 and MERS-CoV, that have spilled over into humans in the past 20 years to cause severe disease. The bnAbs provide new opportunities for antibody-based interventions and key insights for developing pan-betacoronavirus vaccines.
ABSTRACT
The emergence of current SARS-CoV-2 variants of concern (VOCs) and potential future spillovers of SARS-like coronaviruses into humans pose a major threat to human health and the global economy 1-7 . Development of broadly effective coronavirus vaccines that can mitigate these threats is needed 8, 9 . Notably, several recent studies have revealed that vaccination of recovered COVID-19 donors results in enhanced nAb responses compared to SARS-CoV-2 infection or vaccination alone 10-13 . Here, we utilized a targeted donor selection strategy to isolate a large panel of broadly neutralizing antibodies (bnAbs) to sarbecoviruses from two such donors. Many of the bnAbs are remarkably effective in neutralization against sarbecoviruses that use ACE2 for viral entry and a substantial fraction also show notable binding to non-ACE2-using sarbecoviruses. The bnAbs are equally effective against most SARS-CoV-2 VOCs and many neutralize the Omicron variant. Neutralization breadth is achieved by bnAb binding to epitopes on a relatively conserved face of the receptor binding domain (RBD) as opposed to strain-specific nAbs to the receptor binding site that are commonly elicited in SARS-CoV-2 infection and vaccination 14-18 . Consistent with targeting of conserved sites, select RBD bnAbs exhibited in vivo protective efficacy against diverse SARS-like coronaviruses in a prophylaxis challenge model. The generation of a large panel of potent bnAbs provides new opportunities and choices for next-generation antibody prophylactic and therapeutic applications and, importantly, provides a molecular basis for effective design of pan-sarbecovirus vaccines.
ABSTRACT
Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.
Subject(s)
Antibodies , Benchmarking , Phylogeny , B-LymphocytesABSTRACT
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.
ABSTRACT
The problem of reconstructing a string from its error-prone copies, the trace reconstruction problem, was introduced by Vladimir Levenshtein two decades ago. While there has been considerable theoretical work on trace reconstruction, practical solutions have only recently started to emerge in the context of two rapidly developing research areas: immunogenomics and DNA data storage. In immunogenomics, traces correspond to mutated copies of genes, with mutations generated naturally by the adaptive immune system. In DNA data storage, traces correspond to noisy copies of DNA molecules that encode digital data, with errors being artifacts of the data retrieval process. In this paper, we introduce several new trace generation models and open questions relevant to trace reconstruction for immunogenomics and DNA data storage, survey theoretical results on trace reconstruction, and highlight their connections to computational biology. Throughout, we discuss the applicability and shortcomings of known solutions and suggest future research directions.
Subject(s)
Genomics , Immunogenetics , B-Lymphocytes/immunology , Databases, Genetic , Germ Cells , Humans , Receptors, Antigen, B-Cell/genetics , Receptors, Antigen, B-Cell/immunology , Receptors, Antigen, T-Cell/genetics , Receptors, Antigen, T-Cell/immunology , T-Lymphocytes/immunology , Whole Genome SequencingABSTRACT
Antibody-antigen binding relies on the specific interaction of amino acids at the paratope-epitope interface. The predictability of antibody-antigen binding is a prerequisite for de novo antibody and (neo-)epitope design. A fundamental premise for the predictability of antibody-antigen binding is the existence of paratope-epitope interaction motifs that are universally shared among antibody-antigen structures. In a dataset of non-redundant antibody-antigen structures, we identify structural interaction motifs, which together compose a commonly shared structure-based vocabulary of paratope-epitope interactions. We show that this vocabulary enables the machine learnability of antibody-antigen binding on the paratope-epitope level using generative machine learning. The vocabulary (1) is compact, less than 104 motifs; (2) distinct from non-immune protein-protein interactions; and (3) mediates specific oligo- and polyreactive interactions between paratope-epitope pairs. Our work leverages combined structure- and sequence-based learning to demonstrate that machine-learning-driven predictive paratope and epitope engineering is feasible.
Subject(s)
Antigen-Antibody Reactions/immunology , Binding Sites, Antibody/immunology , Epitopes/immunology , Amino Acid Motifs , Amino Acid Sequence , Antibodies/chemistry , Antibodies/immunology , Complementarity Determining Regions/chemistry , Epitopes/chemistry , Machine Learning , Protein BindingABSTRACT
The V(DD)J recombination is currently viewed as an aberrant and inconsequential variant of the canonical V(D)J recombination. Moreover, since the classical 12/23 rule for the V(D)J recombination fails to explain the V(DD)J recombination, the molecular mechanism of tandem D-D fusions has remained unknown since they were discovered three decades ago. Revealing this mechanism is a biomedically important goal since tandem fusions contribute to broadly neutralizing antibodies with ultralong CDR3s. We reveal previously overlooked cryptic nonamers in the recombination signal sequences of human IGHD genes and demonstrate that these nonamers explain the vast majority of tandem fusions in human repertoires. We further reveal large clonal lineages formed by tandem fusions in antigen-stimulated immunosequencing data sets, suggesting that such data sets contain many more tandem fusions than previously thought and that about a quarter of large clonal lineages with unusually long CDR3s are generated through tandem fusions. Finally, we developed the SEARCH-D algorithm for identifying D genes in mammalian genomes and applied it to the recently completed Vertebrate Genomes Project assemblies, nearly doubling the number of mammalian species with known D genes. Our analysis revealed cryptic nonamers in RSSs of many mammalian genomes, thus demonstrating that the V(DD)J recombination is not a "bug" but an important feature preserved throughout mammalian evolution.