Analysis of antibody binding specificities in twin and SNP-genotyped cohorts reveals that antiviral antibody epitope selection is a heritable trait

Precision arbovirus serology with a pan-arbovirus peptidome

Arthropod-borne viruses represent a crucial public health threat. Current arboviral serology assays are either labor intensive or incapable of distinguishing closely related viruses, and many zoonotic arboviruses that may transition to humans lack any serologic assays. In this study, we present a programmable phage display platform, ArboScan, that evaluates antibody binding to overlapping peptides that represent the proteomes of 691 human and zoonotic arboviruses. We confirm that ArboScan provides detailed antibody binding information from animal sera, human sera, and an arthropod blood meal. ArboScan identifies distinguishing features of antibody responses based on exposure history in a Colombian cohort of Zika patients. Finally, ArboScan details epitope level information that rapidly identifies candidate epitopes with potential protective significance. ArboScan thus represents a resource for characterizing human and animal arbovirus antibody responses at cohort scale.

Efficient encoding of large antigenic spaces by epitope prioritization with Dolphyn

We investigate a relatively underexplored component of the gut-immune axis by profiling the antibody response to gut phages using Phage Immunoprecipitation Sequencing (PhIP-Seq). To cover large antigenic spaces, we develop Dolphyn, a method that uses machine learning to select peptides from protein sets and compresses the proteome through epitope-stitching. Dolphyn compresses the size of a peptide library by 78% compared to traditional tiling, increasing the antibody-reactive peptides from 10% to 31%. We find that the immune system develops antibodies to human gut bacteria-infecting viruses, particularly E.coli-infecting Myoviridae. Cost-effective PhIP-Seq libraries designed with Dolphyn enable the assessment of a wider range of proteins in a single experiment, thus facilitating the study of the gut-immune axis.

A common NFKB1 variant detected through antibody analysis in UK Biobank predicts risk of infection and allergy

Infectious agents contribute significantly to the global burden of diseases through both acute infection and their chronic sequelae. We leveraged the UK Biobank to identify genetic loci that influence humoral immune response to multiple infections. From 45 genome-wide association studies in 9,611 participants from UK Biobank, we identified NFKB1 as a locus associated with quantitative antibody responses to multiple pathogens, including those from the herpes, retro-, and polyoma-virus families. An insertion-deletion variant thought to affect NFKB1 expression (rs28362491), was mapped as the likely causal variant and could play a key role in regulation of the immune response. Using 121 infection- and inflammation-related traits in 487,297 UK Biobank participants, we show that the deletion allele was associated with an increased risk of infection from diverse pathogens but had a protective effect against allergic disease. We propose that altered expression of NFKB1, as a result of the deletion, modulates hematopoietic pathways and likely impacts cell survival, antibody production, and inflammation. Taken together, we show that disruptions to the tightly regulated immune processes may tip the balance between exacerbated immune responses and allergy, or increased risk of infection and impaired resolution of inflammation.

Phage Immunoprecipitation Sequencing (PhIP-Seq) for Analyzing Antibody Epitope Repertoires Against Food Antigens

While thousands of food and environmental allergens have been reported, conventional methods for allergy testing typically rely on measuring immunoglobulin E (IgE) binding against panels of dozens to hundreds of antigens. Beyond IgE, also the specificity of other Ig (sub-)classes such as IgG4, has gained interest because of a potential protective role toward allergy.Phage immunoprecipitation sequencing (PhIP-Seq) allows to study hundreds of thousands of rationally selected peptide antigens and to resolve binding specificities of different Ig classes. This technology combines synthetic DNA libraries encoding antigens, with the display on the surface of T7 bacteriophages and next-generation sequencing (NGS) for quantitative readouts. Thereby binding of entire Ig repertoires can be measured to detect the exact epitopes of food allergens and to study potential cross-reactivity.In this chapter, we provide a summary of both the key experimental steps and various strategies for analyzing PhIP-Seq datasets, as well as comparing the advantages and disadvantages of this methodology for measuring antibody responses against food antigens.

Pan-viral serology uncovers distinct virome patterns as risk predictors of hepatocellular carcinoma and intrahepatic cholangiocarcinoma

This study evaluates the pan-serological profiles of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) compared to several diseased and non-diseased control populations to identify risk factors and biomarkers of liver cancer. We used phage immunoprecipitation sequencing, an anti-viral antibody screening method using a synthetic-phage-displayed human virome epitope library, to screen patient serum samples for exposure to over 1,280 strains of pathogenic and non-pathogenic viruses. Using machine learning methods to develop an HCC or iCCA viral score, we discovered that both viral scores were positively associated with several liver function markers in two separate at-risk populations independent of viral hepatitis status. The HCC score predicted all-cause mortality over 8 years in patients with chronic liver disease at risk of HCC, while the viral hepatitis status was not predictive of survival. These results suggest that non-hepatitis viral infections may contribute to HCC and iCCA development and could be biomarkers in at-risk populations.

A robust pipeline for high-content, high-throughput immunophenotyping reveals age- and genetics-dependent changes in blood leukocytes

High-dimensional flow cytometry is the gold standard to study the human immune system in large cohorts. However, large sample sizes increase inter-experimental variation because of technical and experimental inaccuracies introduced by batch variability. Our high-throughput sample processing pipeline in combination with 28-color flow cytometry focuses on increased throughput (192 samples/experiment) and high reproducibility. We implemented quality control checkpoints to reduce technical and experimental variation. Finally, we integrated FlowSOM clustering to facilitate automated data analysis and demonstrate the reproducibility of our pipeline in a study with 3,357 samples. We reveal age-associated immune dynamics in 2,300 individuals, signified by decreasing T and B cell subsets with age. In addition, by combining genetic analyses, our approach revealed unique immune signatures associated with a single nucleotide polymorphism (SNP) that abrogates CD45 isoform splicing. In summary, we provide a versatile and reliable high-throughput, flow cytometry-based pipeline for immune discovery and exploration in large cohorts.

Efficient encoding of large antigenic spaces by epitope prioritization with Dolphyn

We investigated a relatively underexplored component of the gut-immune axis by profiling the antibody response to gut phages using Phage Immunoprecipitation Sequencing (PhIP-Seq). To enhance this approach, we developed Dolphyn, a novel method that uses machine learning to select peptides from protein sets and compresses the proteome through epitope-stitching. Dolphyn improves the fraction of gut phage library peptides bound by antibodies from 10% to 31% in healthy individuals, while also reducing the number of synthesized peptides by 78%. In our study on gut phages, we discovered that the immune system develops antibodies to bacteria-infecting viruses in the human gut, particularly E.coli-infecting Myoviridae. Cost-effective PhIP-Seq libraries designed with Dolphyn enable the assessment of a wider range of proteins in a single experiment, thus facilitating the study of the gut-immune axis.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

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

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

Allergenic food protein consumption is associated with systemic IgG antibody responses in non-allergic individuals

Although food-directed immunoglobulin E (IgE) has been studied in the context of allergies, the prevalence and magnitude of IgG responses against dietary antigens are incompletely characterized in the general population. Here, we measured IgG binding against food and environmental antigens obtained from allergen databases and the immune epitope database (IEDB), represented in a phage displayed library of 58,233 peptides. By profiling blood samples of a large cohort representing the average adult Israeli population (n = 1,003), we showed that many food antigens elicited systemic IgG in up to 50% of individuals. Dietary intake of specific food protein correlated with antibody binding, suggesting that diet can shape the IgG epitope repertoire. Our work documents abundant systemic IgG responses against food antigens and provides a reference map of the exact immunogenic epitopes on a population scale, laying the foundation to unravel the role of food- and environmental antigen-directed antibody binding in disease contexts.

Host genetics impact on SARS-CoV-2 vaccine-induced immunoglobulin levels and dynamics: The role of TP53, ABO, APOE, ACE2, HLA-A, and CRP genes

Background: Development and worldwide availability of safe and effective vaccines against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) to fight severe symptoms of coronavirus disease 2019 (COVID-19) and block the pandemic have been a great achievement and stimulated researchers on understanding the efficacy and duration of different vaccine types. Methods: We investigated the levels of anti-SARS-CoV-2 antibodies (IgG) and neutralizing antibodies (NAbs) in 195 healthy adult subjects belonging to the staff of the University-Hospital of Ferrara (Italy) starting from 15 days up to 190 days (about 6 months) after the second dose of the BNT162b2 (Pfizer-BioNTech) mRNA-based vaccine (n = 128) or ChAdOx1 (AstraZeneca) adenovirus-based vaccine (n = 67) using a combined approach of serological and genomics investigations. Results: A strong correlation between IgG and NAb levels was detected during the 190 days of follow-up (r 2 = 0.807; p < 0.0001) and was confirmed during the first 90 days (T1) after vaccination (r 2 = 0.789; p = 0.0001) and 91-190 days (T2) after vaccination (r 2 = 0.764; p = 0.0001) for both vaccine types (r 2 = 0.842; p = 0.0001 and r 2 = 0.780; p = 0.0001 for mRNA- and adenovirus-based vaccine, respectively). In addition to age (p < 0.01), sex (p = 0.03), and type of vaccine (p < 0.0001), which partially accounted for the remarkable individual differences observed in the antibody levels and dynamics, interesting genetic determinants appeared as significant modifiers of both IgG and NAb responses among the selected genes investigated (TP53, rs1042522; APOE, rs7412/rs429358; ABO, rs657152; ACE2, rs2285666; HLA-A rs2571381/rs2499; CRP, rs2808635/rs876538; LZTFL1, rs35044562; OAS3, rs10735079; SLC6A20, rs11385942; CFH, rs1061170; and ACE1, ins/del, rs4646994). In detail, regression analysis and mean antibody level comparison yielded appreciable differences after genotype stratification (P1 and P2, respectively, for IgG and NAb distribution) in the whole cohort and/or in the mRNA-based vaccine in the following genes: TP53, rs1042522 (P1 = 0.03; P2 = 0.04); ABO, rs657152 (P1 = 0.01; P2 = 0.03); APOE, rs7412/rs429358 (P1 = 0.0018; P2 = 0.0002); ACE2, rs2285666 (P1 = 0.014; P2 = 0.009); HLA-A, rs2571381/rs2499 (P1 = 0.02; P2 = 0.03); and CRP, rs2808635/rs876538 (P1 = 0.01 and P2 = 0.09). Conclusion: High- or low-responsive subjects can be identified among healthy adult vaccinated subjects after targeted genetic screening. This suggests that favorable genetic backgrounds may support the progression of an effective vaccine-induced immune response, though no definite conclusions can be drawn on the real effectiveness ascribed to a specific vaccine or to the different extent of a genotype-driven humoral response. The interplay between data from the polygenic predictive markers and serological screening stratified by demogeographic information can help to recognize the individual humoral response, accounting for ethnic and geographical differences, in both COVID-19 and anti-SARS-CoV-2 vaccinations.

Current advances in antibody-based serum biomarker studies: From protein microarray to phage display

PURPOSE

This review aims to summarize the technological advances in the field of antibody-based biomarker studies by proteome microarray and phage display. In addition, the possible development directions of this field are also discussed.

EXPERIMENTAL DESIGN

We have focused on the antibody profiling by proteome microarray and phage display, including the technological advances, the tools/resources constructed, and the characteristics of both platforms.

RESULTS

With the help of tools/resources and technological advances in proteome microarray and phage display, the efficiency of profiling antibody-based biomarkers in serum samples has been greatly improved.

CONCLUSIONS

In the past few years, proteome microarray and phage display, especially the latter one, have already demonstrated their capacity and efficiency for biomarker identification. In the near future, we believe that more antibody-based biomarkers could be identified, and some of them could eventually be developed into real clinical applications.

Systemic antibody responses against human microbiota flagellins are overrepresented in chronic fatigue syndrome patients

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a debilitating disease with an unclear etiology and pathogenesis. Both an involvement of the immune system and gut microbiota dysbiosis have been implicated in its pathophysiology. However, potential interactions between adaptive immune responses and the microbiota in ME/CFS have been incompletely characterized. Here, we profiled antibody responses of patients with severe ME/CFS and healthy controls against microbiota and viral antigens represented as a phage-displayed 244,000 variant library. Patients with severe ME/CFS exhibited distinct serum antibody epitope repertoires against flagellins of Lachnospiraceae bacteria. Training machine learning algorithms on this antibody-binding data demonstrated that immune responses against gut microbiota represent a unique layer of information beyond standard blood tests, providing improved molecular diagnostics for ME/CFS. Together, our results point toward an involvement of the microbiota-immune axis in ME/CFS and lay the foundation for comparative studies with inflammatory bowel diseases and illnesses characterized by long-term fatigue symptoms, including post-COVID-19 syndrome.

Detecting antibody reactivities in Phage ImmunoPrecipitation Sequencing data

Phage ImmunoPrecipitation Sequencing (PhIP-Seq) is a recently developed technology to assess antibody reactivity, quantifying antibody binding towards hundreds of thousands of candidate epitopes. The output from PhIP-Seq experiments are read count matrices, similar to RNA-Seq data; however some important differences do exist. In this manuscript we investigated whether the publicly available method edgeR (Robinson et al., Bioinformatics 26(1):139-140, 2010) for normalization and analysis of RNA-Seq data is also suitable for PhIP-Seq data. We find that edgeR is remarkably effective, but improvements can be made and introduce a Bayesian framework specifically tailored for data from PhIP-Seq experiments (Bayesian Enrichment Estimation in R, BEER).

Unbiased discovery of autoantibodies associated with severe COVID-19 via genome-scale self-assembled DNA-barcoded protein libraries

Pathogenic autoreactive antibodies that may be associated with life-threatening coronavirus disease 2019 (COVID-19) remain to be identified. Here, we show that self-assembled genome-scale libraries of full-length proteins covalently coupled to unique DNA barcodes for analysis by sequencing can be used for the unbiased identification of autoreactive antibodies in plasma samples. By screening 11,076 DNA-barcoded proteins expressed from a sequence-verified human ORFeome library, the method, which we named MIPSA (for Molecular Indexing of Proteins by Self-Assembly), allowed us to detect circulating neutralizing type-I and type-III interferon (IFN) autoantibodies in five plasma samples from 55 patients with life-threatening COVID-19. In addition to identifying neutralizing type-I IFN-α and IFN-ω autoantibodies and other previously known autoreactive antibodies in patient plasma, MIPSA enabled the detection of as yet unidentified neutralizing type-III anti-IFN-λ3 autoantibodies that were not seen in healthy plasma samples or in convalescent plasma from ten non-hospitalized individuals with COVID-19. The low cost and simple workflow of MIPSA will facilitate unbiased high-throughput analyses of protein-antibody, protein-protein and protein-small-molecule interactions.

Phage display of environmental protein toxins and virulence factors reveals the prevalence, persistence, and genetics of antibody responses

Microbial exposures are crucial environmental factors that impact healthspan by sculpting the immune system and microbiota. Antibody profiling via Phage ImmunoPrecipitation Sequencing (PhIP-Seq) provides a high-throughput, cost-effective approach for detecting exposure and response to microbial protein products. We designed and constructed a library of 95,601 56-amino acid peptide tiles spanning 14,430 proteins with "toxin" or "virulence factor" keyword annotations. We used PhIP-Seq to profile the antibodies of ∼1,000 individuals against this "ToxScan" library. In addition to enumerating immunodominant antibody epitopes, we studied the age-dependent stability of the ToxScan profile and used a genome-wide association study to find that the MHC-II locus modulates bacterial epitope selection. We detected previously described anti-flagellin antibody responses in a Crohn's disease cohort and identified an association between anti-flagellin antibodies and juvenile dermatomyositis. PhIP-Seq with the ToxScan library is thus an effective tool for studying the environmental determinants of health and disease at cohort scale.

The immunogenetics of viral antigen response is associated with subtype-specific glioma risk and survival

Glioma is a highly fatal cancer with prognostically significant molecular subtypes and few known risk factors. Multiple studies have implicated infections in glioma susceptibility, but evidence remains inconsistent. Genetic variants in the human leukocyte antigen (HLA) region modulate host response to infection and have been linked to glioma risk. In this study, we leveraged genetic predictors of antibody response to 12 viral antigens to investigate the relationship with glioma risk and survival. Genetic reactivity scores (GRSs) for each antigen were derived from genome-wide-significant (p < 5 × 10-8) variants associated with immunoglobulin G antibody response in the UK Biobank cohort. We conducted parallel analyses of glioma risk and survival for each GRS and HLA alleles imputed at two-field resolution by using data from 3,418 glioma-affected individuals subtyped by somatic mutations and 8,156 controls. Genetic reactivity scores to Epstein-Barr virus (EBV) ZEBRA and EBNA antigens and Merkel cell polyomavirus (MCV) VP1 antigen were associated with glioma risk and survival (Bonferroni-corrected p < 0.01). GRSZEBRA and GRSMCV were associated in opposite directions with risk of IDH wild-type gliomas (ORZEBRA = 0.91, p = 0.0099/ORMCV = 1.11, p = 0.0054). GRSEBNA was associated with both increased risk for IDH mutated gliomas (OR = 1.09, p = 0.040) and improved survival (HR = 0.86, p = 0.010). HLA-DQA1∗03:01 was significantly associated with decreased risk of glioma overall (OR = 0.85, p = 3.96 × 10-4) after multiple testing adjustment. This systematic investigation of the role of genetic determinants of viral antigen reactivity in glioma risk and survival provides insight into complex immunogenomic mechanisms of glioma pathogenesis. These results may inform applications of antiviral-based therapies in glioma treatment.

Epstein-Barr virus as a driver of multiple sclerosis

Accumulating evidence implicates Epstein-Barr virus (EBV) as an etiological factor in multiple sclerosis (MS). If EBV is a "driver" that causes antiviral immunity with associated autoimmune components rather than a "trigger" that unleashes self-perpetuating autoimmunity, then elimination of EBV would be a rational therapy for MS.

Prevention of MS Requires Intervention on the Causes of the Disease: Reconciling Genes, Epigenetics, and Epstein Barr Virus

Prevention of multiple sclerosis requires intervention on modifiable causes of the condition making it necessary to establish what those causes are. MS is often stated to be a polygenic disease, with causal contributions from environmental factors and gene-environment interactions, implying an additive and independent relationship of these factors. Mechanistically there are no independent contributions of genes or environmental factors to traits. This model is unrealistic but still useful and underlies the concept of heritability, a foundational parameter in population genetics. However, it perpetuates a debate on an irreconcilable dichotomy about whether MS is primarily genetic or environmental. In particular, epidemiological evidence now exists for a causal, possibly even necessary, role for Epstein Barr Virus in MS. The additive model makes it unintuitive to reconcile MS as a genetic disease but also independently a viral illness. In this perspective it is argued that starting from a realistic interaction only model, based on broadly accepted biological premises, and working forward to explain why the classical additive model gives useful results, there is actually no paradox. An integrated approach using population genetic studies, immunology and molecular virology offers a particularly promising route to establish the elusive role of EBV in MS pathology, as EBV is a large and complex virus and its latency, dysregulated in most EBV-related pathologies, is hard to study in vivo. This approach may offer a route to prevention of MS altogether.

A family matter: Anti-viral antibody responses

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