Toward a more accurate view of human B-cell repertoire by next-generation sequencing, unbiased repertoire capture and single-molecule barcoding

Accelerating antibody discovery and optimization with high-throughput experimentation and machine learning

The integration of high-throughput experimentation and machine learning is transforming data-driven antibody engineering, revolutionizing the discovery and optimization of antibody therapeutics. These approaches employ extensive datasets comprising antibody sequences, structures, and functional properties to train predictive models that enable rational design. This review highlights the significant advancements in data acquisition and feature extraction, emphasizing the necessity of capturing both sequence and structural information. We illustrate how machine learning models, including protein language models, are used not only to enhance affinity but also to optimize other crucial therapeutic properties, such as specificity, stability, viscosity, and manufacturability. Furthermore, we provide practical examples and case studies to demonstrate how the synergy between experimental and computational approaches accelerates antibody engineering. Finally, this review discusses the remaining challenges in fully realizing the potential of artificial intelligence (AI)-powered antibody discovery pipelines to expedite therapeutic development.

Beyond glycan barriers: non-cognate ligands and protein mimicry approaches to elicit broadly neutralizing antibodies for HIV-1

Human immunodeficiency virus type 1 (HIV-1) vaccine immunogens capable of inducing broadly neutralizing antibodies (bNAbs) remain obscure. HIV-1 evades immune responses through enormous diversity and hides its conserved vulnerable epitopes on the envelope glycoprotein (Env) by displaying an extensive immunodominant glycan shield. In elite HIV-1 viremic controllers, glycan-dependent bNAbs targeting conserved Env epitopes have been isolated and are utilized as vaccine design templates. However, immunological tolerance mechanisms limit the development of these antibodies in the general population. The well characterized bNAbs monoclonal variants frequently exhibit extensive levels of somatic hypermutation, a long third heavy chain complementary determining region, or a short third light chain complementarity determining region, and some exhibit poly-reactivity to autoantigens. This review elaborates on the obstacles to engaging and manipulating the Env glycoprotein as an effective immunogen and describes an alternative reverse vaccinology approach to develop a novel category of bNAb-epitope-derived non-cognate immunogens for HIV-1 vaccine design.

B-Cell Receptor Repertoire: Recent Advances in Autoimmune Diseases

In the field of contemporary medicine, autoimmune diseases (AIDs) are a prevalent and debilitating group of illnesses. However, they present extensive and profound challenges in terms of etiology, pathogenesis, and treatment. A major reason for this is the elusive pathophysiological mechanisms driving disease onset. Increasing evidence suggests the indispensable role of B cells in the pathogenesis of autoimmune diseases. Interestingly, B-cell receptor (BCR) repertoires in autoimmune diseases display a distinct skewing that can provide insights into disease pathogenesis. Over the past few years, advances in high-throughput sequencing have provided powerful tools for analyzing B-cell repertoire to understand the mechanisms during the period of B-cell immune response. In this paper, we have provided an overview of the mechanisms and analytical methods for generating BCR repertoire diversity and summarize the latest research progress on BCR repertoire in autoimmune diseases, including systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), primary Sjögren's syndrome (pSS), multiple sclerosis (MS), and type 1 diabetes (T1D). Overall, B-cell repertoire analysis is a potent tool to understand the involvement of B cells in autoimmune diseases, facilitating the creation of innovative therapeutic strategies targeting specific B-cell clones or subsets.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

'Immunization during ART and ATI for HIV-1 vaccine discovery/development'

PURPOSE OF REVIEW

Explore whether immunization with germline-targeting Env immunogens during ART, followed by ATI, leads to the identification of viral envelope glycoproteins (Envs) that promote and guide the full maturation of broadly neutralizing antibody responses.

RECENT FINDINGS

The HIV-1 envelope glycoprotein (Env) does not efficiently engage the germline precursors of broadly neutralizing antibodies (bnAbs). However, Env-derived proteins specifically designed to precisely do that, have been recently developed. These 'germline-targeting' Env immunogens activate naïve B cells that express the germline precursors of bnAbs but by themselves cannot guide their maturation towards their broadly neutralizing forms. This requires sequential immunizations with heterologous sets of Envs. These 'booster' Envs are currently unknown.

SUMMARY

Combining germline-targeting Env immunization approaches during ART with ATI could lead to the identification of natural Envs that are responsible for the maturation of broadly neutralizing antibody responses during infection. Such Envs could then serve as booster immunogens to guide the maturation of glBCRs that have become activated by germline-targeting immunogens in uninfected subjects.

Applications of High-Throughput DNA Sequencing to Single-Domain Antibody Discovery and Engineering

Next-generation DNA sequencing (NGS) technologies have made it possible to interrogate antibody repertoires to unprecedented depths, typically via sequencing of cDNAs encoding immunoglobulin variable domains. In the absence of heavy-light chain pairing, the variable domains of heavy chain-only antibodies (HCAbs), referred to as single-domain antibodies (sdAbs), are uniquely amenable to NGS analyses. In this chapter, we provide simple and rapid protocols for producing and sequencing multiplexed immunoglobulin variable domain (VHH, VH, or VL) amplicons derived from a variety of sources using the Illumina MiSeq platform. Generation of such amplicon libraries is relatively inexpensive, requiring no specialized equipment and only a limited set of PCR primers. We also present several applications of NGS to sdAb discovery and engineering, including: (1) evaluation of phage-displayed sdAb library sequence diversity and monitoring of panning experiments; (2) identification of sdAbs of predetermined epitope specificity following competitive elution of phage-displayed sdAb libraries; (3) direct selection of B cells expressing antigen-specific, membrane-bound HCAb using antigen-coupled magnetic beads and identification of antigen-specific sdAbs, and (4) affinity maturation of lead sdAbs using tandem phage display selection and NGS. These methods can easily be adapted to other types of proteins and libraries and expand the utility of in vitro display technology.

Moving beyond Titers

Vaccination to prevent and even eliminate disease is amongst the greatest achievements of modern medicine. Opportunities remain in vaccine development to improve protection across the whole population. A next step in vaccine development is the detailed molecular characterization of individual humoral immune responses against a pathogen, especially the rapidly evolving pathogens. New technologies such as sequencing the immune repertoire in response to disease, immunogenomics/vaccinomics, particularly the individual HLA variants, and high-throughput epitope characterization offer new insights into disease protection. Here, we highlight the emerging technologies that could be used to identify variation within the human population, facilitate vaccine discovery, improve vaccine safety and efficacy, and identify mechanisms of generating immunological memory. In today’s vaccine-hesitant climate, these techniques used individually or especially together have the potential to improve vaccine effectiveness and safety and thus vaccine uptake rates. We highlight the importance of using these techniques in combination to understand the humoral immune response as a whole after vaccination to move beyond neutralizing titers as the standard for immunogenicity and vaccine efficacy, especially in clinical trials.

The methods and advances of adaptive immune receptors repertoire sequencing

The adaptive immune response is a powerful tool, capable of recognizing, binding to, and neutralizing a vast number of internal and external threats via T or B lymphatic receptors with widespread sets of antigen specificities. The emergence of high-throughput sequencing technology and bioinformatics provides opportunities for research in the fields of life sciences and medicine. The analysis and annotation for immune repertoire data can reveal biologically meaningful information, including immune prediction, target antigens, and effective evaluation. Continuous improvements of the immunological repertoire sequencing methods and analysis tools will help to minimize the experimental and calculation errors and realize the immunological information to meet the clinical requirements. That said, the clinical application of adaptive immune repertoire sequencing requires appropriate experimental methods and standard analytical tools. At the population cell level, we can acquire the overview of cell groups, but the information about a single cell is not obtained accurately. The information that is ignored may be crucial for understanding the heterogeneity of each cell, gene expression and drug response. The combination of high-throughput sequencing and single-cell technology allows us to obtain single-cell information with low-cost and high-throughput. In this review, we summarized the current methods and progress in this area.

Resistance of SARS-CoV-2 variants to neutralization by antibodies induced in convalescent patients with COVID-19

Administration of convalescent plasma or neutralizing monoclonal antibodies (mAbs) is a potent therapeutic option for coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. However, SARS-CoV-2 variants with mutations in the spike protein have emerged in many countries. To evaluate the efficacy of neutralizing antibodies induced in convalescent patients against emerging variants, we isolate anti-spike mAbs from two convalescent COVID-19 patients infected with prototypic SARS-CoV-2 by single-cell sorting of immunoglobulin-G-positive (IgG+) memory B cells. Anti-spike antibody induction is robust in these patients, and five mAbs have potent neutralizing activities. The efficacy of most neutralizing mAbs and convalescent plasma samples is maintained against B.1.1.7 and mink cluster 5 variants but is significantly decreased against variants B.1.351 from South Africa and P.1 from Brazil. However, mAbs with a high affinity for the receptor-binding domain remain effective against these neutralization-resistant variants. Rapid spread of these variants significantly impacts antibody-based therapies and vaccine strategies against SARS-CoV-2.

Ultrasensitive DNA Immune Repertoire Sequencing Using Unique Molecular Identifiers

BACKGROUND

Immune repertoire sequencing of the T-cell receptor can identify clonotypes that have expanded as a result of antigen recognition or hematological malignancies. However, current sequencing protocols display limitations with nonuniform amplification and polymerase-induced errors during sequencing. Here, we developed a sequencing method that overcame these issues and applied it to γδ T cells, a cell type that plays a unique role in immunity, autoimmunity, homeostasis of intestine, skin, adipose tissue, and cancer biology.

METHODS

The ultrasensitive immune repertoire sequencing method used PCR-introduced unique molecular identifiers. We constructed a 32-panel assay that captured the full diversity of the recombined T-cell receptor delta loci in γδ T cells. The protocol was validated on synthetic reference molecules and blood samples of healthy individuals.

RESULTS

The 32-panel assay displayed wide dynamic range, high reproducibility, and analytical sensitivity with single-nucleotide resolution. The method corrected for sequencing-depended quantification bias and polymerase-induced errors and could be applied to both enriched and nonenriched cells. Healthy donors displayed oligoclonal expansion of γδ T cells and similar frequencies of clonotypes were detected in both enrichment and nonenriched samples.

CONCLUSIONS

Ultrasensitive immune repertoire sequencing strategy enables quantification of individual and specific clonotypes in a background that can be applied to clinical as well as basic application areas. Our approach is simple, flexible, and can easily be implemented in any molecular laboratory.

[Isolation of anti-SARS-CoV-2 neutralizing monoclonal antibodies cross effective to variants aiming at antibody therapy]

We isolated five mAbs with potent neutralizing activities against SARS-CoV-2 from two convalescent COVID-19 patients infected with prototype virus. Among them, the 9-105 antibody that have a highest affinity for the receptor-binding domain (RBD), cross-neutralize variants, such as B.1.1.7 (alfa), mink cluster 5 variant, B.1.351 (beta), P.1 (gamma), C.37 (lambda), B.1.617.1 (kappa), B.1.617.2 (delta) and B.1.621 (mu). A single amino acid mutation at K417 of RBD decreased neutralization sensitivity of 9-105. A 9-105 homology model revealed that 9-105 light chain binds to RBD including K417 by the same angle as ACE2.

GLaMST: grow lineages along minimum spanning tree for b cell receptor sequencing data

Background

B cell affinity maturation enables B cells to generate high-affinity antibodies. This process involves somatic hypermutation of B cell immunoglobulin receptor (BCR) genes and selection by their ability to bind antigens. Lineage trees are used to describe this microevolution of B cell immunoglobulin genes. In a lineage tree, each node is one BCR sequence that mutated from the germinal center and each directed edge represents a single base mutation, insertion or deletion. In BCR sequencing data, the observed data only contains a subset of BCR sequences in this microevolution process. Therefore, reconstructing the lineage tree from experimental data requires algorithms to build the tree based on partially observed tree nodes.

Results

We developed a new algorithm named Grow Lineages along Minimum Spanning Tree (GLaMST), which efficiently reconstruct the lineage tree given observed BCR sequences that correspond to a subset of the tree nodes. Through comparison using simulated and real data, GLaMST outperforms existing algorithms in simulations with high rates of mutation, insertion and deletion, and generates lineage trees with smaller size and closer to ground truth according to tree features that highly correlated with selection pressure.

Conclusions

GLaMST outperforms state-of-art in reconstruction of the BCR lineage tree in both efficiency and accuracy. Integrating it into existing BCR sequencing analysis frameworks can significant improve lineage tree reconstruction aspect of the analysis.

A VH1-69 antibody lineage from an infected Chinese donor potently neutralizes HIV-1 by targeting the V3 glycan supersite

An oligomannose patch around the V3 base of HIV-1 envelope glycoprotein (Env) is recognized by multiple classes of broadly neutralizing antibodies (bNAbs). Here, we investigated the bNAb response to the V3 glycan supersite in an HIV-1-infected Chinese donor by Env-specific single B cell sorting, structural and functional studies, and longitudinal analysis of antibody and virus repertoires. Monoclonal antibodies 438-B11 and 438-D5 were isolated that potently neutralize HIV-1 with moderate breadth, are encoded by the VH1-69 germline gene, and have a disulfide-linked long HCDR3 loop. Crystal structures of Env-bound and unbound antibodies revealed heavy chain-mediated recognition of the glycan supersite with a unique angle of approach and a critical role of the intra-HCDR3 disulfide. The mechanism of viral escape was examined via single-genome amplification/sequencing and glycan mutations around the N332 supersite. Our findings further emphasize the V3 glycan supersite as a prominent target for Env-based vaccine design.

Quantitative evaluation of protective antibody response induced by hepatitis E vaccine in humans

Efficacy evaluation through human trials is crucial for advancing a vaccine candidate to clinics. Next-generation sequencing (NGS) can be used to quantify B cell repertoire response and trace antibody lineages during vaccination. Here, we demonstrate this application with a case study of Hecolin®, the licensed vaccine for hepatitis E virus (HEV). Four subjects are administered the vaccine following a standard three-dose schedule. Vaccine-induced antibodies exhibit a high degree of clonal diversity, recognize five conformational antigenic sites of the genotype 1 HEV p239 antigen, and cross-react with other genotypes. Unbiased repertoire sequencing is performed for seven time points over six months of vaccination, with maturation pathways characterize for a set of vaccine-induced antibodies. In addition to dynamic repertoire profiles, NGS analysis reveals differential patterns of HEV-specific antibody lineages and highlights the necessity of the long vaccine boost. Together, our study presents a quantitative strategy for vaccine evaluation in small-scale human studies.

Systematic Profiling of Full-Length Ig and TCR Repertoire Diversity in Rhesus Macaque through Long Read Transcriptome Sequencing

The diversity of Ig and TCR repertoires is a focal point of immunological studies. Rhesus macaques (Macaca mulatta) are key for modeling human immune responses, placing critical importance on the accurate annotation and quantification of their Ig and TCR repertoires. However, because of incomplete reference resources, the coverage and accuracy of the traditional targeted amplification strategies for profiling rhesus Ig and TCR repertoires are largely unknown. In this study, using long read sequencing, we sequenced four Indian-origin rhesus macaque tissues and obtained high-quality, full-length sequences for over 6000 unique Ig and TCR transcripts, without the need for sequence assembly. We constructed, to our knowledge, the first complete reference set for the constant regions of all known isotypes and chain types of rhesus Ig and TCR repertoires. We show that sequence diversity exists across the entire variable regions of rhesus Ig and TCR transcripts. Consequently, existing strategies using targeted amplification of rearranged variable regions comprised of V(D)J gene segments miss a significant fraction (27-53% and 42-49%) of rhesus Ig/TCR diversity. To overcome these limitations, we designed new rhesus-specific assays that remove the need for primers conventionally targeting variable regions and allow single cell level Ig and TCR repertoire analysis. Our improved approach will enable future studies to fully capture rhesus Ig and TCR repertoire diversity and is applicable for improving annotations in any model organism.

Targeting broadly neutralizing antibody precursors: a naïve approach to vaccine design

PURPOSE OF REVIEW

It is believed that broadly neutralizing antibodies (bNAbs) will be an important component of an effective HIV-1 vaccine. Several immunogens have been designed that can target specific precursor B cells as a first step in a vaccine strategy to elicit bNAbs.

RECENT FINDINGS

Germline-targeting immunogens have been developed that specifically engage precursors of reproducible classes of anti-HIV antibodies, such as VRC01-class and apex-directed bNAbs. However, these precursors represent only a small portion of the immune repertoire and any antigen will inherently present off-target epitopes to the immune system that may confound bNAb development. Novel animal models are being utilized to understand the competitive fitness of bNAb precursors in the context of immunization with germline-targeting immunogens. In parallel, immunogen design efforts are being pursued to favor the development of bNAb responses over off-target responses following immunization. New studies of bNAb precursor interactions with glycosylated Env variants can inform prime-boost regimens geared towards accelerating bNAb development.

SUMMARY

Germline-targeting immunogens hold promise as a first step in eliciting a bNAb response through vaccination. A better understating of how efficiently germline-targeting immunogens can specifically target rare bNAb precursors is emerging. In addition, a more comprehensive structure-based understanding of critical barriers to bNAb elicitation, as well as commonalities between bNAb classes can further inform vaccine design.

Prolonged evolution of the memory B cell response induced by a replicating adenovirus-influenza H5 vaccine

Induction of an antibody response capable of recognizing highly diverse strains is a major obstacle to the development of vaccines for viruses such as HIV and influenza. Here, we report the dynamics of B cell expansion and evolution at the single-cell level after vaccination with a replication-competent adenovirus type 4 recombinant virus expressing influenza H5 hemagglutinin. Fluorescent H1 or H5 probes were used to quantitate and isolate peripheral blood B cells and their antigen receptors. We observed increases in H5-specific antibody somatic hypermutation and potency for several months beyond the period of active viral replication that was not detectable at the serum level. Individual broad and potent antibodies could be isolated, including one stem-specific antibody that is part of a new multidonor class. These results demonstrate prolonged evolution of the B cell response for months after vaccination and should be considered in efforts to evaluate or boost vaccine-induced immunity.

Comparisons of the antibody repertoires of a humanized rodent and humans by high throughput sequencing

The humanization of animal model immune systems by genetic engineering has shown great promise for antibody discovery, tolerance studies and for the evaluation of vaccines. Assessment of the baseline antibody repertoires of unimmunized model animals will be useful as a benchmark for future immunization experiments. We characterized the heavy chain and kappa light chain antibody repertoires of a model animal, the OmniRat, by high throughput antibody sequencing and made use of two novel datasets for comparison to human repertoires. Intra-animal and inter-animal repertoire comparisons reveal a high level of conservation in antibody diversity between the lymph node and spleen and between members of the species. Multiple differences were found in both the heavy and kappa chain repertoires between OmniRats and humans including gene segment usage, CDR3 length distributions, class switch recombination, somatic hypermutation levels and in features of V(D)J recombination. The Inference and Generation of Repertoires (IGoR) software tool was used to model recombination in VH regions which allowed for the quantification of some of these differences. Diversity estimates of the OmniRat heavy chain repertoires almost reached that of humans, around two orders of magnitude less. Despite variation between the species repertoires, a high frequency of OmniRat clonotypes were also found in the human repertoire. These data give insights into the development and selection of humanized animal antibodies and provide actionable information for use in vaccine studies.

Exploiting B Cell Receptor Analyses to Inform on HIV-1 Vaccination Strategies

The human antibody repertoire is generated by the recombination of different gene segments as well as by processes of somatic mutation. Together these mechanisms result in a tremendous diversity of antibodies that are able to combat various pathogens including viruses and bacteria, or malignant cells. In this review, we summarize the opportunities and challenges that are associated with the analyses of the B cell receptor repertoire and the antigen-specific B cell response. We will discuss how recent advances have increased our understanding of the antibody response and how repertoire analyses can be exploited to inform on vaccine strategies, particularly against HIV-1.

An MPER antibody neutralizes HIV-1 using germline features shared among donors

The membrane-proximal external region (MPER) of HIV-1 envelope glycoprotein (Env) can be targeted by neutralizing antibodies of exceptional breadth. MPER antibodies usually have long, hydrophobic CDRH3s, lack activity as inferred germline precursors, are often from the minor IgG3 subclass, and some are polyreactive, such as 4E10. Here we describe an MPER broadly neutralizing antibody from the major IgG1 subclass, PGZL1, which shares germline V/D-region genes with 4E10, has a shorter CDRH3, and is less polyreactive. A recombinant sublineage variant pan-neutralizes a 130-isolate panel at 1.4 μg/ml (IC50). Notably, a germline revertant with mature CDR3s neutralizes 12% of viruses and still binds MPER after DJ reversion. Crystal structures of lipid-bound PGZL1 variants and cryo-EM reconstruction of an Env-PGZL1 complex reveal how these antibodies recognize MPER and viral membrane. Discovery of common genetic and structural elements among MPER antibodies from different patients suggests that such antibodies could be elicited using carefully designed immunogens.

Understanding and measuring human B-cell tolerance and its breakdown in autoimmune disease

The maintenance of immunological tolerance of B lymphocytes is a complex and critical process that must be implemented as to avoid the detrimental development of autoreactivity and possible autoimmunity. Murine models have been invaluable to elucidate many of the key components in B-cell tolerance; however, translation to human homeostatic and pathogenic immune states can be difficult to assess. Functional autoreactive, flow cytometric, and single-cell cloning assays have proven to be critical in deciphering breaks in B-cell tolerance within autoimmunity; however, newer approaches to assess human B-cell tolerance may prove to be vital in the further exploration of underlying tolerance defects. In this review, we supply a comprehensive overview of human immune tolerance checkpoints with associated mechanisms of enforcement, and highlight current and future methodologies which are likely to benefit future studies into the mechanisms that become defective in human autoimmune conditions.

Development of a Potent and Protective Germline-Like Antibody Lineage Against Zika Virus in a Convalescent Human

Zika virus (ZIKV) specific neutralizing antibodies hold great promise for antibody-based interventions and vaccine design against ZIKV infection. However, their development in infected patients remains unclear. Here, we applied next-generation sequencing (NGS) to probe the dynamic development of a potent and protective ZIKV E DIII-specific antibody ZK2B10 isolated from a ZIKV convalescent individual. The unbiased repertoire analysis showed dramatic changes in the usage of antibody variable region germline genes. However, lineage tracing of ZK2B10 revealed limited somatic hypermutation and transient expansion during the 12 months following the onset of symptoms. The NGS-derived, germline-like ZK2B10 somatic variants neutralized ZIKV potently and protected mice from lethal challenge of ZIKV without detectable cross-reactivity with Dengue virus (DENV). Site-directed mutagenesis identified two residues within the λ chain, N31 and S91, that are essential to the functional maturation of ZK2B10. The repertoire and lineage features unveiled here will help elucidate the developmental process and protective potential of E DIII-directed antibodies against ZIKV infection.

Deep sequencing of B cell receptor repertoire

Immune repertoire is a collection of enormously diverse adaptive immune cells within an individual. As the repertoire shapes and represents immunological conditions, identification of clones and characterization of diversity are critical for understanding how to protect ourselves against various illness such as infectious diseases and cancers. Over the past several years, fast growing technologies for high throughput sequencing have facilitated rapid advancement of repertoire research, enabling us to observe the diversity of repertoire at an unprecedented level. Here, we focus on B cell receptor (BCR) repertoire and review approaches to B cell isolation and sequencing library construction. These experiments should be carefully designed according to BCR regions to be interrogated, such as heavy chain full length, complementarity determining regions, and isotypes. We also highlight preprocessing steps to remove sequencing and PCR errors with unique molecular index and bioinformatics techniques. Due to the nature of massive sequence variation in BCR, caution is warranted when interpreting repertoire diversity from error-prone sequencing data. Furthermore, we provide a summary of statistical frameworks and bioinformatics tools for clonal evolution and diversity. Finally, we discuss limitations of current BCR-seq technologies and future perspectives on advances in repertoire sequencing.

Next-generation sequencing and single-cell RT-PCR reveal a distinct variable gene usage of porcine antibody repertoire following PEDV vaccination

Porcine epidemic diarrhea virus (PEDV) is the most common diarrhea-causing pathogen in newborn piglets. The clarifications of the overall antibody repertoire and antigen-specific antibody repertoire are essential to provide important insights into the B-cell response and reshape new vaccines. Here, we applied next-generation sequencing (NGS) technology to investigate immunoglobulin (Ig) variable (V) gene segment usage of swine B-cells from peripheral blood lymphocytes (PBL) and mesenteric lymph node (MLN) cells following PEDV vaccination. We identified the transcripts of all functional Ig V-genes in antibody repertoire. IgHV1S2, IgKV1-11, and IgLV3-4 were the most prevalent gene segments for heavy, kappa, and lambda chains, respectively, in PBL and MLN. Unlike previous studies, IgKV1, instead of IgKV2, and IgLV3, instead of IgLV8, were the prevalent Ig V-gene families for kappa and lambda light chains, respectively. We further examined the antibody repertoire of PEDV spike-specific B cells by single-cell RT-PCR. In contrast to the overall antibody repertoire, Ig V-gene segments of PEDV spike-specific B cells preferentially adopted IgHV1-4 and IgHV1-14 for heavy chain, IgKV1-11 for kappa chain, and IgLV3-3 for lambda chain. These results represent a comprehensive analysis to characterize the Ig V-gene segment usage in the overall and PEDV spike-specific antibody repertoire in PBL and MLN.

Survival analysis of DNA mutation motifs with penalized proportional hazards

Antibodies, an essential part of our immune system, develop through an intricate process to bind a wide array of pathogens. This process involves randomly mutating DNA sequences encoding these antibodies to find variants with improved binding, though mutations are not distributed uniformly across sequence sites. Immunologists observe this nonuniformity to be consistent with "mutation motifs", which are short DNA subsequences that affect how likely a given site is to experience a mutation. Quantifying the effect of motifs on mutation rates is challenging: a large number of possible motifs makes this statistical problem high dimensional, while the unobserved history of the mutation process leads to a nontrivial missing data problem. We introduce an ℓ 1-penalized proportional hazards model to infer mutation motifs and their effects. In order to estimate model parameters, our method uses a Monte Carlo EM algorithm to marginalize over the unknown ordering of mutations. We show that our method performs better on simulated data compared to current methods and leads to more parsimonious models. The application of proportional hazards to mutation processes is, to our knowledge, novel and formalizes the current methods in a statistical framework that can be easily extended to analyze the effect of other biological features on mutation rates.

High-Quality Library Preparation for NGS-Based Immunoglobulin Germline Gene Inference and Repertoire Expression Analysis

Next generation sequencing (NGS) of immunoglobulin (Ig) repertoires (Rep-seq) enables examination of the adaptive immune system at an unprecedented level. Applications include studies of expressed repertoires, gene usage, somatic hypermutation levels, Ig lineage tracing and identification of genetic variation within the Ig loci through inference methods. All these applications require starting libraries that allow the generation of sequence data with low error rate and optimal representation of the expressed repertoire. Here, we provide detailed protocols for the production of libraries suitable for human Ig germline gene inference and Ig repertoire studies. Various parameters used in the process were tested in order to demonstrate factors that are critical to obtain high quality libraries. We demonstrate an improved 5'RACE technique that reduces the length constraints of Illumina MiSeq based Rep-seq analysis but allows for the acquisition of sequences upstream of Ig V genes, useful for primer design. We then describe a 5' multiplex method for library preparation, which yields full length V(D)J sequences suitable for genotype identification and novel gene inference. We provide comprehensive sets of primers targeting IGHV, IGKV, and IGLV genes. Using the optimized protocol, we produced IgM, IgG, IgK, and IgL libraries and analyzed them using the germline inference tool IgDiscover to identify expressed germline V alleles. This process additionally uncovered three IGHV, one IGKV, and six IGLV novel alleles in a single individual, which are absent from the IMGT reference database, highlighting the need for further study of Ig genetic variation. The library generation protocols presented here enable a robust means of analyzing expressed Ig repertoires, identifying novel alleles and producing individualized germline gene databases from humans.

Application of long read sequencing to determine expressed antigen diversity in Trypanosoma brucei infections

Antigenic variation is employed by many pathogens to evade the host immune response, and Trypanosoma brucei has evolved a complex system to achieve this phenotype, involving sequential use of variant surface glycoprotein (VSG) genes encoded from a large repertoire of ~2,000 genes. T. brucei express multiple, sometimes closely related, VSGs in a population at any one time, and the ability to resolve and analyse this diversity has been limited. We applied long read sequencing (PacBio) to VSG amplicons generated from blood extracted from batches of mice sacrificed at time points (days 3, 6, 10 and 12) post-infection with T. brucei TREU927. The data showed that long read sequencing is reliable for resolving variant differences between VSGs, and demonstrated that there is significant expressed diversity (449 VSGs detected across 20 mice) and across the timeframe of study there was a clear semi-reproducible pattern of expressed diversity (median of 27 VSGs per sample at day 3 post infection (p.i.), 82 VSGs at day 6 p.i., 187 VSGs at day 10 p.i. and 132 VSGs by day 12 p.i.). There was also consistent detection of one VSG dominating expression across replicates at days 3 and 6, and emergence of a second dominant VSG across replicates by day 12. The innovative application of ecological diversity analysis to VSG reads enabled characterisation of hierarchical VSG expression in the dataset, and resulted in a novel method for analysing such patterns of variation. Additionally, the long read approach allowed detection of mosaic VSG expression from very few reads–the earliest in infection that such events have been detected. Therefore, our results indicate that long read analysis is a reliable tool for resolving diverse gene expression profiles, and provides novel insights into the complexity and nature of VSG expression in trypanosomes, revealing significantly higher diversity than previously shown and the ability to identify mosaic gene formation early during the infection process.

Antigenic variation is a system whereby pathogens switch identity of a protein that is exposed to the host adaptive immune response as a way of remaining one step ahead and avoiding being detected. African trypanosomes have evolved a spectacularly elaborate system of antigenic variation, with variants being used from a library of ~2,000 genes. Our ability to understand how this rich repository is used has been hampered by the resolution of available technologies to discriminate between what can be closely related gene variants. We have applied a long read sequencing technology, which generates sequence information for the whole length of the antigen gene variants, thereby avoiding having to try and piece together antigen sequences from lots of small fragments, the pitfall of standard sequencing. Applying this technology to material taken at specific time points from batches of mice infected with trypanosomes reveals that the diversity of variants is much higher than previously suspected, and that there is a clear semi-predictable pattern in the gene expression. Additionally, using this technology we have been able to detect the presence of ‘mosaic’ genes, which are created by stitching together fragments from several donor genes in the library, much earlier in infection than has been shown previously. Therefore, we shed new light on the complexity of antigenic variation and show that long read sequencing will be a very useful tool in analysing and understanding the expression patterns of closely related genes, and how pathogens use them to cause persistent infections and disease.

Genetic and structural insights into broad neutralization of hepatitis C virus by human VH1-69 antibodies

An effective vaccine to the antigenically diverse hepatitis C virus (HCV) must target conserved immune epitopes. Here, we investigate cross-neutralization of HCV genotypes by broadly neutralizing antibodies (bNAbs) encoded by the relatively abundant human gene family V H 1-69. We have deciphered the molecular requirements for cross-neutralization by this unique class of human antibodies from crystal structures of HCV E2 in complex with bNAbs. An unusually high binding affinity is found for germ line-reverted versions of VH1-69 precursor antibodies, and neutralization breadth is acquired during affinity maturation. Deep sequencing analysis of an HCV-immune B cell repertoire further demonstrates the importance of the V H 1-69 gene family in the generation of HCV bNAbs. This study therefore provides critical insights into immune recognition of HCV with important implications for rational vaccine design.

Filtering Next-Generation Sequencing of the Ig Gene Repertoire Data Using Antibody Structural Information

Next-generation sequencing of the Ig gene repertoire (Ig-seq) produces large volumes of information at the nucleotide sequence level. Such data have improved our understanding of immune systems across numerous species and have already been successfully applied in vaccine development and drug discovery. However, the high-throughput nature of Ig-seq means that it is afflicted by high error rates. This has led to the development of error-correction approaches. Computational error-correction methods use sequence information alone, primarily designating sequences as likely to be correct if they are observed frequently. In this work, we describe an orthogonal method for filtering Ig-seq data, which considers the structural viability of each sequence. A typical natural Ab structure requires the presence of a disulfide bridge within each of its variable chains to maintain the fold. Our Ab Sequence Selector (ABOSS) uses the presence/absence of this bridge as a way of both identifying structurally viable sequences and estimating the sequencing error rate. On simulated Ig-seq datasets, ABOSS is able to identify more than 99% of structurally viable sequences. Applying our method to six independent Ig-seq datasets (one mouse and five human), we show that our error calculations are in line with previous experimental and computational error estimates. We also show how ABOSS is able to identify structurally impossible sequences missed by other error-correction methods.

Sequencing the peripheral blood B and T cell repertoire - Quantifying robustness and limitations

The adaptive immune response generates a large repertoire of T cells with T-cell receptors (TCR_alpha and TCR_beta) and B cells with immunoglobulins (Ig). The repertoire changes in response to antigen stimulation both through amplification of specific cells (clonal expansion) as well as somatic hypermutation of immunoglobulins. Alterations of the immune repertoire have been observed in response to acute disease, such as external pathogens, or chronic diseases, such as autoimmunity and cancer. Here we establish experimental and analytical protocols for quantifying the peripheral blood of healthy human individuals by profiling the immune repertoire for the Complementarity determining region 3 (CDR3) of the variable regions of TCR_beta (CDRβ3) and the IgG heavy chain (CDRH1, CDRH2, CDRH3). The results demonstrate that 40 ml of blood are sufficient to reliably capture the 10,000 most common TCR_beta and 1000 most common IgG and determine their relative frequency in the circulation. We conclude that by using an accessible sample size of human PBMC one is able to robustly monitor alterations in the immune repertoire.

Using High-Throughput Sequencing to Characterize the Development of the Antibody Repertoire During Infections: A Case Study of HIV-1

High throughput sequencing (HTS) approaches have only recently been applied to describing the antibody/B-cell repertoire in fine detail, but these data sets have already become critical to the design of vaccines and therapeutics, and monitoring of cancer immunotherapy. As a case study, we describe the potential and present limitations of HTS studies of the Ab repertoire during infection with HIV-1. Most of the present studies restrict their analyses to lineages of specific bnAbs. We discuss future initiatives to expand this type of analysis to more complete repertoires and to improve comparing and sharing of these Ab repertoire data across studies and institutions.

Profiling of the TCRβ repertoire in non-model species using high-throughput sequencing

In recent years, immune repertoire profiling with high-throughput sequencing (HTS) has advanced our understanding of adaptive immunity. However, fast progress in the field applied mostly to human and mouse research, with only few studies devoted to other model vertebrates. We present the first in-depth characterization of the T-cell receptor (TCR) repertoire in a non-model mammal (bank vole, Myodes glareolus), widely used in ecological and evolutionary research. We used RNA from spleens, 5′RACE and HTS to describe V and J segments of TCRβ, qualitatively characterize preferential V–J segment usage and CDR3 length distribution. Overall orthology to murine genes was preserved, with 11 J and 37 V genes found in voles (although 3 V genes lacked a close orthologue). Further, we implemented unique molecular identifiers for quantitative analysis of CDR3 repertoire with stringent error correction. A conservative, lower bound estimation of the TCRβ repertoire was similar to that found for mice (1.7–2.3 × 10⁵ clonotypes). We hope that by providing an easy-to-follow molecular protocol and on-line bioinformatics tools that do not require reference sequences (AmpliTCR and AmpliCDR3), we will encourage HTS immune repertoire profiling in other non-model vertebrates, thus opening new research avenues in e.g. comparative immunology, ecology and evolutionary biology.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGT's underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.

5' Rapid Amplification of cDNA Ends and Illumina MiSeq Reveals B Cell Receptor Features in Healthy Adults, Adults With Chronic HIV-1 Infection, Cord Blood, and Humanized Mice

Using 5′ rapid amplification of cDNA ends, Illumina MiSeq, and basic flow cytometry, we systematically analyzed the expressed B cell receptor (BCR) repertoire in 14 healthy adult PBMCs, 5 HIV-1+ adult PBMCs, 5 cord blood samples, and 3 HIS-CD4/B mice, examining the full-length variable region of μ, γ, α, κ, and λ chains for V-gene usage, somatic hypermutation (SHM), and CDR3 length. Adding to the known repertoire of healthy adults, Illumina MiSeq consistently detected small fractions of reads with high mutation frequencies including hypermutated μ reads, and reads with long CDR3s. Additionally, the less studied IgA repertoire displayed similar characteristics to that of IgG. Compared to healthy adults, the five HIV-1 chronically infected adults displayed elevated mutation frequencies for all μ, γ, α, κ, and λ chains examined and slightly longer CDR3 lengths for γ, α, and λ. To evaluate the reconstituted human BCR sequences in a humanized mouse model, we analyzed cord blood and HIS-CD4/B mice, which all lacked the typical SHM seen in the adult reference. Furthermore, MiSeq revealed identical unmutated IgM sequences derived from separate cell aliquots, thus for the first time demonstrating rare clonal members of unmutated IgM B cells by sequencing.

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

HIV Broadly Neutralizing Antibodies: VRC01 and Beyond

Developing an effective prophylaxis HIV-1 vaccine is likely to require the elicitation of broadly neutralizing antibodies (bnAbs). As the HIV-1 envelope (Env) glycoprotein - the sole target of bnAbs - has evolved multiple mechanisms to evade antibody neutralization, the processes for bnAb generation are highly selective and time-consuming. Benefiting from antibody isolation technologies of single B cell culturing and direct single B cell sorting and cloning, a new generation of monoclonal bnAbs has been isolated since 2009, exhibiting remarkable breadths and potencies, thus breaking through a nearly 20-year-long limit of four monoclonal bnAbs with moderate breadth and potency. The discovery of a long list of monoclonal bnAbs has provided in-depth understanding of the sites of vulnerability on the HIV-1 Env and the complexity of human B cell immunology to generate such responses, thus presenting both guidance and challenges to move the Env immunogen design effort forward.

Next-Generation DNA Sequencing of VH/VL Repertoires: A Primer and Guide to Applications in Single-Domain Antibody Discovery

Immunogenetic analyses of expressed antibody repertoires are becoming increasingly common experimental investigations and are critical to furthering our understanding of autoimmunity, infectious disease, and cancer. Next-generation DNA sequencing (NGS) technologies have now made it possible to interrogate antibody repertoires to unprecedented depths, typically by sequencing of cDNAs encoding immunoglobulin variable domains. In this chapter, we describe simple, fast, and reliable methods for producing and sequencing multiplex PCR amplicons derived from the variable regions (V_H, V_HH or V_L) of rearranged immunoglobulin heavy and light chain genes using the Illumina MiSeq platform. We include complete protocols and primer sets for amplicon sequencing of V_H/V_HH/V_L repertoires directly from human, mouse, and llama lymphocytes as well as from phage-displayed V_H/V_HH/V_L libraries; these can be easily be adapted to other types of amplicons with little modification. The resulting amplicons are diverse and representative, even using as few as 10³ input B cells, and their generation is relatively inexpensive, requiring no special equipment and only a limited set of primers. In the absence of heavy-light chain pairing, single-domain antibodies are uniquely amenable to NGS analyses. We present a number of applications of NGS technology useful in discovery of single-domain antibodies from phage display libraries, including: (i) assessment of library functionality; (ii) confirmation of desired library randomization; (iii) estimation of library diversity; and (iv) monitoring the progress of panning experiments. While the case studies presented here are of phage-displayed single-domain antibody libraries, the principles extend to other types of in vitro display libraries.

How B-Cell Receptor Repertoire Sequencing Can Be Enriched with Structural Antibody Data

Next-generation sequencing of immunoglobulin gene repertoires (Ig-seq) allows the investigation of large-scale antibody dynamics at a sequence level. However, structural information, a crucial descriptor of antibody binding capability, is not collected in Ig-seq protocols. Developing systematic relationships between the antibody sequence information gathered from Ig-seq and low-throughput techniques such as X-ray crystallography could radically improve our understanding of antibodies. The mapping of Ig-seq datasets to known antibody structures can indicate structurally, and perhaps functionally, uncharted areas. Furthermore, contrasting naïve and antigenically challenged datasets using structural antibody descriptors should provide insights into antibody maturation. As the number of antibody structures steadily increases and more and more Ig-seq datasets become available, the opportunities that arise from combining the two types of information increase as well. Here, we review how these data types enrich one another and show potential for advancing our knowledge of the immune system and improving antibody engineering.

Differences in Allelic Frequency and CDRH3 Region Limit the Engagement of HIV Env Immunogens by Putative VRC01 Neutralizing Antibody Precursors

Elicitation of broadly neutralizing antibodies remains a long-standing goal of HIV vaccine research. Although such antibodies can arise during HIV-1 infection, gaps in our knowledge of their germline, pre-immune precursor forms, as well as on their interaction with viral Env, limit our ability to elicit them through vaccination. Studies of broadly neutralizing antibodies from the VRC01-class provide insight into progenitor B cell receptors (BCRs) that could develop into this class of antibodies. Here, we employed high-throughput heavy chain variable region (VH)/light chain variable region (VL) deep sequencing, combined with biophysical, structural, and modeling antibody analyses, to interrogate circulating potential VRC01-progenitor BCRs in healthy individuals. Our study reveals that not all humans are equally predisposed to generate VRC01-class antibodies, not all predicted progenitor VRC01-expressing B cells can bind to Env, and the CDRH3 region of germline VRC01 antibodies influence their ability to recognize HIV-1. These findings will be critical to the design of optimized immunogens that should consider CDRH3 interactions.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

Reverse Transcription Past Products of Guanine Oxidation in RNA Leads to Insertion of A and C opposite 8-Oxo-7,8-dihydroguanine and A and G opposite 5-Guanidinohydantoin and Spiroiminodihydantoin Diastereomers

Reactive oxygen species, both endogenous and exogenous, can damage nucleobases of RNA and DNA. Among the nucleobases, guanine has the lowest redox potential, making it a major target of oxidation. Although RNA is more prone to oxidation than DNA is, oxidation of guanine in RNA has been studied to a significantly lesser extent. One of the reasons for this is that many tools that were previously developed to study oxidation of DNA cannot be used on RNA. In the study presented here, the lack of a method for seeking sites of modification in RNA where oxidation occurs is addressed. For this purpose, reverse transcription of RNA containing major products of guanine oxidation was used. Extension of a DNA primer annealed to an RNA template containing 8-oxo-7,8-dihydroguanine (OG), 5-guanidinohydantoin (Gh), or the R and S diastereomers of spiroiminodihydantoin (Sp) was studied under standing start conditions. SuperScript III reverse transcriptase is capable of bypassing these lesions in RNA inserting predominantly A opposite OG, predominantly G opposite Gh, and almost an equal mixture of A and G opposite the Sp diastereomers. These data should allow RNA sequencing of guanine oxidation products by following characteristic mutation signatures formed by the reverse transcriptase during primer elongation past G oxidation sites in the template RNA strand.

Hidden Lineage Complexity of Glycan-Dependent HIV-1 Broadly Neutralizing Antibodies Uncovered by Digital Panning and Native-Like gp140 Trimer

Germline precursors and intermediates of broadly neutralizing antibodies (bNAbs) are essential to the understanding of humoral response to HIV-1 infection and B-cell lineage vaccine design. Using a native-like gp140 trimer probe, we examined antibody libraries constructed from donor-17, the source of glycan-dependent PGT121-class bNAbs recognizing the N332 supersite on the HIV-1 envelope glycoprotein. To facilitate this analysis, a digital panning method was devised that combines biopanning of phage-displayed antibody libraries, 900 bp long-read next-generation sequencing, and heavy/light (H/L)-paired antibodyomics. In addition to single-chain variable fragments resembling the wild-type bNAbs, digital panning identified variants of PGT124 (a member of the PGT121 class) with a unique insertion in the heavy chain complementarity-determining region 1, as well as intermediates of PGT124 exhibiting notable affinity for the native-like trimer and broad HIV-1 neutralization. In a competition assay, these bNAb intermediates could effectively compete with mouse sera induced by a scaffolded BG505 gp140.681 trimer for the N332 supersite. Our study thus reveals previously unrecognized lineage complexity of the PGT121-class bNAbs and provides an array of library-derived bNAb intermediates for evaluation of immunogens containing the N332 supersite. Digital panning may prove to be a valuable tool in future studies of bNAb diversity and lineage development.

Germline-targeting immunogens

In 2009, Dimitrov's group reported that the inferred germline (iGL) forms of several HIV-1 broadly neutralizing antibodies (bNAbs) did not display measurable binding to a recombinant gp140 Env protein (derived from the dual-tropic 89.6 virus), which was efficiently recognized by the mature (somatically mutated) antibodies. At that time, a small number of bNAbs were available, but in the following years, the implementation of high-throughput B-cell isolation and sequencing assays and of screening methodologies facilitated the isolation of greater numbers of bNAbs from infected subjects. Using these newest bNAbs, and a wide range of diverse recombinant Envs, we and others confirmed the observations made by Dimitrov's group. The results from these studies created a paradigm shift in our collective thinking as to why recombinant Envs are ineffective in eliciting bNAbs and has led to the "germline-targeting" immunization approach. Here we discuss this approach in detail: what has been done so far, the advantages and limitations of the current germline-targeting immunogens and of the animal models used to test them, and we conclude with a few thoughts about future directions in this area of research.

Specific microRNA library of IFN-τ on bovine endometrial epithelial cells

IFN-τ is specifically secreted by the conceptus in ruminants during early pregnancy, and it plays a vital role in the immunological function of pregnancy. However, its mechanism involving microRNA (miRNA) is still not well understood. Deep sequencing was used to explore the specific miRNA library of IFN-τ on bovine endometrial epithelial cells (bEECs). The results showed that 574 known bovine miRNAs and 109 novel miRNAs were identified. We found 74 differentially expressed miRNAs, including 30 commonly expressed miRNAs in the experiment. Then, qPCR verification of six selected miRNAs showed that they corresponded with the sequencing data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed significant enrichment of predicted target genes of differentially expressed miRNAs, including influenza A, herpes simplex infection, antigen processing and presentation, viral myocarditis, TNF signaling pathway, graft-versus-host disease, and allograft rejection. These results may provide important contributions to the immune response during early pregnancy in ruminants, but further studies are need to verify the proposed cellular/immunological effects and role of specific miRNA as biomarkers in vivo.

Application of nonsense-mediated primer exclusion (NOPE) for preparation of unique molecular barcoded libraries

Background

Recently we proposed efficient method to exclude undesirable primers at any stage of amplification reaction, here termed NOPE (NOnsense-mediated Primer Exclusion). According to this method, added oligonucleotide overlapping with the 3′-end of unwanted amplification primer (NOPE oligo) simultaneously provides a template for its elongation. This elongation disrupts specificity of unwanted primer, preventing its further participation in PCR. The suggested approach allows to rationally manage the course of PCR reactions in order to facilitate analysis of complex DNA mixtures as well as to perform multistage PCR bypassing intermediate purification steps.

Results

Here we apply NOPE method to DNA library preparation for the high-throughput sequencing (HTS) with the PCR-based introduction of unique molecular identifiers (UMI). We show that NOPE oligo efficiently neutralizes UMI-containing oligonucleotides after introduction of UMI into sample DNA molecules, thus allowing to proceed with further amplification steps without purification and associated loss of starting material. At the same time, NOPE oligo does not affect the efficiency of target PCR amplification.

Conclusion

We describe a simple, robust and cheap modification of UMI-labeled HTS libraries preparation procedure, that allows to bypass purification step and thus to preserve starting material which may be limited, e.g. circulating tumor DNA, circulating fetal DNA, or small amounts of isolated cells of interest. Furthermore, demonstrated simplicity and robustness of NOPE method should make it popular in various PCR protocols.

Active evolution of memory B-cells specific to viral gH/gL/pUL128/130/131 pentameric complex in healthy subjects with silent human cytomegalovirus infection

Human cytomegalovirus (HCMV) can cause life-threatening infection in immunosuppressed patients, and in utero infection that may lead to birth defects. No vaccine is currently available. HCMV infection in healthy subjects is generally asymptomatic, and virus persists as latent infection for life. Host immunity is effective against reactivation and super-infection with another strain. Thus, vaccine candidates able to elicit immune responses similar to those of natural infection may confer protection. Since neutralization is essential for prophylactic vaccines, it is important to understand how antiviral antibodies are developed in natural infection. We hypothesized that the developmental path of antibodies in seropositive subjects could be unveiled by interrogating host B-cell repertoires using unique genetic signature sequences of mAbs. Towards this goal, we isolated 56 mAbs from three healthy donors with different neutralizing titers. Antibodies specific to the gH/gL/pUL128/130/131 pentameric complex were more potent in neutralization than those to gB. Using these mAbs as probes, patterns of extended lineage development for B-cells and evidence of active antibody maturation were revealed in two donors with higher neutralizing titers. Importantly, such patterns were limited to mAbs specific to the pentamer, but none to gB. Thus, memory B-cells with antiviral function such as neutralization were active during latent infection in the two donors, and this activity was responsible for their higher neutralizing titers. Our results indicated that memory B-cells of neutralizing capacity could be frequently mobilized in host, probably responding to silent viral episodes, further suggesting that neutralizing antibodies could play a role in control of recurrent infection.

Novel Approaches to Analyze Immunoglobulin Repertoires

Analysis of immunoglobulin (Ig) repertoires aims to comprehend Ig diversity with the goal of predicting humoral immune responses in the context of infection, vaccination, autoimmunity, and malignancies. The first next-generation sequencing (NGS) analyses of bulk B cell populations dramatically advanced sampling depth over previous low-throughput single-cell-based protocols, albeit at the expense of accuracy and loss of chain-pairing information. In recent years the field has substantially differentiated, with bulk analyses becoming more accurate while single-cell approaches have gained in throughput. Additionally, new platforms striving to combine high throughput and chain pairing have been developed as well as various computational tools for analysis. Here we review the developments of the past 4-5 years and discuss the open challenges.

MAGERI: Computational pipeline for molecular-barcoded targeted resequencing

Unique molecular identifiers (UMIs) show outstanding performance in targeted high-throughput resequencing, being the most promising approach for the accurate identification of rare variants in complex DNA samples. This approach has application in multiple areas, including cancer diagnostics, thus demanding dedicated software and algorithms. Here we introduce MAGERI, a computational pipeline that efficiently handles all caveats of UMI-based analysis to obtain high-fidelity mutation profiles and call ultra-rare variants. Using an extensive set of benchmark datasets including gold-standard biological samples with known variant frequencies, cell-free DNA from tumor patient blood samples and publicly available UMI-encoded datasets we demonstrate that our method is both robust and efficient in calling rare variants. The versatility of our software is supported by accurate results obtained for both tumor DNA and viral RNA samples in datasets prepared using three different UMI-based protocols.