Proteomic analysis of influenza haemagglutinin-specific antibodies following vaccination reveals convergent immunoglobulin variable region signatures

Integrated proteomic and metabolomic analysis of plasma reveals regulatory pathways and key elements in thyroid cancer

Thyroid cancer (TC) is the most common endocrine malignancy with increasing incidence in recent years. Fine-needle aspiration biopsy (FNAB), as a gold standard for the initial evaluation of thyroid nodules, fails to cover all the cytopathologic conditions resulting in overdiagnosis. There is an urgent need for a better classification of thyroid cancer from benign thyroid nodules (BTNs). Here, data independent acquisition (DIA)-based proteomics and untargeted metabolomics in plasma samples of 10 patients with TC and 15 patients with BTNs were performed. Key proteins and metabolites were identified specific to TC, and an independent cohort was used to validate the potential biomarkers using enzyme-linked immunosorbent assay (ELISA). In total, 1429 proteins and 1172 metabolites were identified. Principal component analysis showed a strong overlap at the proteomic level and a significant discrimination at the metabolomic level between the two groups, indicating a more drastic disturbance in the metabolome of thyroid cancer. Integrated analysis of proteomics and metabolomics shows glycerophospholipid metabolism and arachidonic acid metabolism as key regulatory pathways. Furthermore, a multi-omics biomarker panel was developed consisting of LCAT, GPX3 and leukotriene B4. Based on the AUC value for the discovery set, the classification performance was 0.960. The AUC value of the external validation set was 0.930. Altogether, our results will contribute to the clinical application of potential biomarkers in the diagnosis of thyroid cancer.

Immunoglobulin repertoire restriction characterizes the serological responses of patients with predominantly antibody deficiency

BACKGROUND

Predominantly antibody deficiency (PAD) is the most common category of inborn errors of immunity and is underpinned by impaired generation of appropriate antibody diversity and quantity. In the clinic, responses are interrogated by assessment of vaccination responses, which is central to many PAD diagnoses. However, the composition of the generated antibody repertoire is concealed from traditional quantitative measures of serological responses. Leveraging modern mass spectrometry-based proteomics (MS-proteomics), it is possible to elaborate the molecular features of specific antibody repertoires, which may address current limitations of diagnostic vaccinology.

OBJECTIVES

We sought to evaluate serum antibody responses in patients with PAD following vaccination with a neo-antigen (severe acute respiratory syndrome coronavirus-2 vaccination) using MS-proteomics.

METHODS

Following severe acute respiratory syndrome coronavirus-2 vaccination, serological responses in individuals with PAD and healthy controls (HCs) were assessed by anti-S1 subunit ELISA and neutralization assays. Purified anti-S1 subunit IgG and IgM was profiled by MS-proteomics for IGHV subfamily usage and somatic hypermutation analysis.

RESULTS

Twelve patients with PAD who were vaccine-responsive were recruited with 11 matched vaccinated HCs. Neutralization and end point anti-S1 titers were lower in PAD. All subjects with PAD demonstrated restricted anti-S1 IgG antibody repertoires, with usage of <5 IGHV subfamilies (median: 3; range 2-4), compared to ≥5 for the 11 HC subjects (P < .001). IGHV3-7 utilization was far less common in patients with PAD than in HCs (2 of 12 vs 10 of 11; P = .001). Amino acid substitutions due to somatic hypermutation per subfamily did not differ between groups. Anti-S1 IgM was present in 64% and 50% of HC and PAD cohorts, respectively, and did not differ significantly between HCs and patients with PAD.

CONCLUSIONS

This study demonstrates the breadth of anti-S1 antibodies elicited by vaccination at the proteome level and identifies stereotypical restriction of IGHV utilization in the IgG repertoire in patients with PAD compared with HC subjects. Despite uniformly pauci-clonal antibody repertoires some patients with PAD generated potent serological responses, highlighting a possible limitation of traditional serological techniques. These findings suggest that IgG repertoire restriction is a key feature of antibody repertoires in PAD.

Vaccination History, Body Mass Index, Age, and Baseline Gene Expression Predict Influenza Vaccination Outcomes

Seasonal influenza is a primary public health burden in the USA and globally. Annual vaccination programs are designed on the basis of circulating influenza viral strains. However, the effectiveness of the seasonal influenza vaccine is highly variable between seasons and among individuals. A number of factors are known to influence vaccination effectiveness including age, sex, and comorbidities. Here, we sought to determine whether whole blood gene expression profiling prior to vaccination is informative about pre-existing immunological status and the immunological response to vaccine. We performed whole transcriptome analysis using RNA sequencing (RNAseq) of whole blood samples obtained prior to vaccination from 275 participants enrolled in an annual influenza vaccine trial. Serological status prior to vaccination and 28 days following vaccination was assessed using the hemagglutination inhibition assay (HAI) to define baseline immune status and the response to vaccination. We find evidence that genes with immunological functions are increased in expression in individuals with higher pre-existing immunity and in those individuals who mount a greater response to vaccination. Using a random forest model, we find that this set of genes can be used to predict vaccine response with a performance similar to a model that incorporates physiological and prior vaccination status alone. A model using both gene expression and physiological factors has the greatest predictive power demonstrating the potential utility of molecular profiling for enhancing prediction of vaccine response. Moreover, expression of genes that are associated with enhanced vaccination response may point to additional biological pathways that contribute to mounting a robust immunological response to the seasonal influenza vaccine.

Convergent antibody evolution and clonotype expansion following influenza virus vaccination

Recent advances in high-throughput single cell sequencing have opened up new avenues into the investigation of B cell receptor (BCR) repertoires. In this study, PBMCs were collected from 17 human participants vaccinated with the split-inactivated influenza virus vaccine during the 2016-2017 influenza season. A combination of Immune Repertoire Capture (IRCTM) technology and IgG sequencing was performed on ~7,800 plasmablast (PB) cells and preferential IgG heavy-light chain pairings were investigated. In some participants, a single expanded clonotype accounted for ~22% of their PB BCR repertoire. Approximately 60% (10/17) of participants experienced convergent evolution, possessing public PBs that were elicited independently in multiple participants. Binding profiles of one private and three public PBs confirmed they were all subtype-specific, cross-reactive hemagglutinin (HA) head-directed antibodies. Collectively, this high-resolution antibody repertoire analysis demonstrated the impact evolution can have on BCRs in response to influenza virus vaccination, which can guide future universal influenza prophylactic approaches.

Plasma proteomic profiling suggests an association between antigen driven clonal B cell expansion and ME/CFS

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is an unexplained chronic, debilitating illness characterized by fatigue, sleep disturbances, cognitive dysfunction, orthostatic intolerance and gastrointestinal problems. Using ultra performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS), we analyzed the plasma proteomes of 39 ME/CFS patients and 41 healthy controls. Logistic regression models, with both linear and quadratic terms of the protein levels as independent variables, revealed a significant association between ME/CFS and the immunoglobulin heavy variable (IGHV) region 3-23/30. Stratifying the ME/CFS group based on self-reported irritable bowel syndrome (sr-IBS) status revealed a significant quadratic effect of immunoglobulin lambda constant region 7 on its association with ME/CFS with sr-IBS whilst IGHV3-23/30 and immunoglobulin kappa variable region 3-11 were significantly associated with ME/CFS without sr-IBS. In addition, we were able to predict ME/CFS status with a high degree of accuracy (AUC = 0.774-0.838) using a panel of proteins selected by 3 different machine learning algorithms: Lasso, Random Forests, and XGBoost. These algorithms also identified proteomic profiles that predicted the status of ME/CFS patients with sr-IBS (AUC = 0.806-0.846) and ME/CFS without sr-IBS (AUC = 0.754-0.780). Our findings are consistent with a significant association of ME/CFS with immune dysregulation and highlight the potential use of the plasma proteome as a source of biomarkers for disease.

Proteomic profiling of precipitated Clostridioides difficile toxin A and B antibodies

Clostridioides difficile infection is the leading cause of nosocomial diarrhoea globally. Immune responses to toxins produced by C. difficile are important in disease progression and outcome. Here, we analysed the anti-toxin A and anti-toxin B serum antibody proteomes following natural infection or vaccination with a C. difficile toxoid A/toxoid B vaccine using a modified miniaturised proteomic approach based on de novo mass spectrometric sequencing. Analysis of immunoglobulin variable region (IgV) subfamily expression in immunoprecipitated toxin A and toxin B antibodies from four and seven participants of a vaccine trial, respectively, revealed a polyclonal proteome with restricted IGHV, IGKV and IGLV subfamily usage. No dominant IGHV subfamily was observed in the toxin A response, however the dominant anti-toxin B heavy (H)-chain was encoded by IGHV3-23. Light (L)-chain usage was convergent for both anti-toxin A and anti-toxin B proteomes with IGKV3-11, 3-15, 3-20 and 4-1 shared among all subjects in both cohorts. Peptide mapping of common IgV families showed extensive public and private amino acid substitutions. The cohort responses to toxin A and toxin B showed limited similarity in shared IGHV subfamilies. L-chain subfamily usage was more similar in the anti-toxin A and anti-toxin B responses, however the mutational signatures for each subfamily were toxin-dependent. Samples taken both post vaccination (n = 5) or at baseline, indicating previous exposure (n = 2), showed similar anti-toxin B IgV subfamily usage and mutational profiles. In summary, this study provides the first sequence-based proteomic analysis of the antibody response to the major disease-mediating toxins of C. difficile, toxin A and toxin B, and demonstrates that despite the potential for extreme diversity, the immunoglobulin repertoire can raise convergent responses to specific pathogens whether through natural infection or following vaccination.

Precipitating anti-dsDNA peptide repertoires in lupus

Anti-double-stranded (ds)DNA autoantibodies are prototypical serological markers of systemic lupus erythematosus (SLE), but little is known about their immunoglobulin variable (IgV) region composition at the level of the secreted (serum) proteome. Here, we use a novel proteomic workflow based on de novo mass spectrometric sequencing of anti-dsDNA precipitins to analyse IgV subfamily expression and mutational signatures of high-affinity, precipitating anti-dsDNA responses. Serum anti-dsDNA proteomes were oligoclonal with shared (public) expression of immunoglobulin (Ig)G heavy chain variable region (IGHV) and kappa chain variable region (IGKV) subfamilies. IgV peptide maps from eight subjects showed extensive public and random (private) amino acid replacement mutations with prominent arginine substitutions across heavy (H)- and light (L)-chains. Shared sets of L-chain complementarity determining region 3 (CDR3) peptides specified by arginine substitutions were sequenced from the dominantly expressed IGKV3-20 subfamily, with changes in expression levels of a clonal L-chain CDR3 peptide by quantitative multiple reaction monitoring (MRM) paralleling the rise and fall of anti-dsDNA levels by Farr radioimmunoassays (RIA). The heavily mutated IgV peptide signatures of precipitating anti-dsDNA autoantibody proteomes reflect the strong selective forces that shape humoral anti-dsDNA responses in germinal centres. Direct sequencing of agarose gel precipitins using microlitre volumes of stored sera streamlines the antibody sequencing workflow and is generalizable to other precipitating serum antibodies.