Donato Gemmati1,2,3*
Giovanna Longo1,2Ines Gallo1,2Juliana Araujo Silva1,2Paola Secchiero1Giorgio Zauli4Stefania Hanau5Angelina Passaro1Patrizia Pellegatti6
Stefano Pizzicotti6Maria Luisa Serino2
Ajay Vikram Singh7
Veronica Tisato1,3- 1Department of Translational Medicine, University of Ferrara, Ferrara, Italy
- 2Centre Haemostasis & Thrombosis, University of Ferrara, Ferrara, Italy
- 3University Centre for Gender Medicine Studies, University of Ferrara, Ferrara, Italy
- 4Research Department, King Khaled Eye Specialist Hospital, Riyadh, Saudi Arabia
- 5Department of Neuroscience & Rehabilitation, University of Ferrara, Ferrara, Italy
- 6Hospital-University of Ferrara, Ferrara, Italy
- 7Department of Chemical & Product Safety, German Federal Institute for Risk Assessment, Berlin, Germany
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 (r2 = 0.807; p < 0.0001) and was confirmed during the first 90 days (T1) after vaccination (r2 = 0.789; p = 0.0001) and 91–190 days (T2) after vaccination (r2 = 0.764; p = 0.0001) for both vaccine types (r2 = 0.842; p = 0.0001 and r2 = 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.
Introduction
In November 2019, the first cases of coronavirus disease 2019 (COVID-19) due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection were diagnosed in Wuhan (China). The virus sequence was released in January 2020 after 54 days from the first declared case, and in March 2020, 63 days after the SARS-CoV-2 sequence was published, the first vaccines were already in the trial phase (Tregoning et al., 2020). In September 2020, the vaccine landscape included 43 candidate formulations under clinical trial evaluation, and the European Medicines Agency (EMA) released authorizations for Pfizer/BioNTech, Moderna, AstraZeneca, and Janssen Pharmaceutical. Two different main strategies of vaccine formulation have been adopted, adenoviral vector-based vaccine (e.g., AstraZeneca or Janssen Pharmaceutical) and mRNA-based vaccine (e.g., Pfizer–BioNTech or Moderna), potentially offering similar efficacy and safety (Falsey et al., 2021; Thomas et al., 2021).
At the beginning of 2021, vaccines seemed to be the complete remedy, even though the presence of conventional mild side effects such as local reactions at the injection site, fever, and fatigue. In addition, some rare episode of vaccine-associated thromboses or more severe vaccine-induced immune thrombocytopenia and thrombosis (VITT) with intracranial hemorrhage and death have been reported in few individuals after adenoviral vector-based vaccine (Polack et al., 2020; Baden et al., 2021; Pavord et al., 2021; Elberry et al., 2022; Klok et al., 2022).
After the first months of the vaccination campaign, novel concerns and questions addressed the efficacy and duration of the circulating antibodies against SARS-CoV-2 induced by vaccines, and this complex issue was the source of alarm among clinicians, researchers, and population, considering the upcoming booster dose scheduled before the end of 2021. Interindividual variables, such as host genetics, sex, age, and different vaccine formulations and the appearance of different spike variants, could be responsible for the observed different antibody response to vaccines. Interestingly, specific HLA regions have been previously recognized by genome-wide association studies (GWAS) as genetic effectors of the interindividual variability in immunoglobulin levels and in response categorization to hepatitis-B vaccine, ascribing to few selected loci, which play a role in the regulation of the humoral immunity (Jonsson et al., 2017; Chung et al., 2019). The association between the HLA system and COVID-19 prognosis has been thoroughly investigated, recognizing potential HLA carriers associated to higher risk of death (Hamming et al., 2004; Grifoni et al., 2020; Huang et al., 2020; Liao et al., 2020; Pretti et al., 2020; Tay et al., 2020; Wang et al., 2020; Zhang et al., 2020). Although recent studies suggested that HLA polymorphisms might unlikely account for the wide disparities in COVID-19 prognosis, they concluded on the potential role of HLA polymorphisms on the development of SARS-CoV-2 immunity also after vaccination by evaluating the association between specific HLA variants and vaccine response (Copley et al., 2021; Shukla et al., 2022). Even though a first extensive report excluded the HLA impact on the immunoglobulin level fluctuation and duration (Ragone et al., 2021), particular attention has been later directed toward a selected locus (HLA-A*03:01) associated to a weak antibody response and negative side effects after Pfizer–BioNTech vaccine (Bolze et al., 2022; Crocchiolo et al., 2022).
Although the connection between the ABO locus and SARS-CoV-2 susceptibility has been clearly established (COVID-19 Host Genetics Initiative, 2021; COVID-19 Host Genetics Initiative, 2022), additional studies investigating the possible effect of the ABO blood type on the risk of SARS-CoV-2 infection, COVID-19 prognosis, and antibody response to anti-SARS-CoV-2 vaccine or infection yielded conflicting results (Gil-Manso et al., 2022; Sgherza et al., 2022; Vicentini et al., 2022; Ziberna et al., 2022). Blood group O has been reported as a possible protective factor for disease progression (Pereira et al., 2022), and recent studies have ascribed to the ABO locus a modifier role on the antibody titers, showing lower levels to B and O types compared to A and AB (de Freitas Dutra et al., 2021). Conversely, higher neutralizing antibody levels have been found to characterize B blood group subjects when compared to the remaining blood groups (Bloch et al., 2021). Moreover, though intriguing hypotheses and mechanisms have been proposed, specific research works addressing the relationships between the phenotypic ABO blood group and anti-SARS-CoV-2 vaccine response did not reach definite results (Sgherza et al., 2022). Finally, positive associations between C-reactive protein (CRP) and immune response after SARS-CoV-2 infection or vaccine have been recently reported, hypothesizing basal or unresolved inflammation as causative reasons (Salvagno et al., 2022a; Gianfagna et al., 2022).
Host genetics and acquired factors modulate the clinical phenotypes of COVID-19 by complex interactions occurring between SARS-CoV-2 and human cells and represent issues investigated by the COVID-19 Host Genetics Initiative (HGI) (Severe Covid et al., 2020; COVID-19 Host Genetics Initiative, 2021); moreover, host genetics may also give precious information on the strength of the individual vaccine immune response and efficacy (Jonsson et al., 2017; Colucci et al., 2021; Castro Dopico and Karlsson Hedestam, 2022; Venkataraman et al., 2022). It is to be noted that the same panel of inherited and acquired factors or combinations of them seems to contribute to both COVID-19 progression and vaccine immune response by sharing common mechanisms and pathways (Gemmati et al., 2020; Gemmati and Tisato, 2020; Tiwari Pandey et al., 2020; Zauli et al., 2020; Cuapio et al., 2022; Milani et al., 2022).
To select in advance individuals at risk with no proper SARS-CoV-2 immune response, it would be extremely useful to predict the development of a weak/strong or a short/long-lasting anti-SARS-CoV-2 immunity either by natural infection or vaccination. In the present study, we assessed the levels of anti-SARS-CoV-2 IgG and neutralizing antibodies in a cohort of adult healthy vaccinated subjects and stratified the results by selected variants in key genes and locus, previously investigated in the COVID-19 prognosis, or in vaccine-induced humoral and cell-mediated response, or as causative of specific mild/common adverse reactions among vaccinated subjects. In detail, the following genes and locus have been considered: ACE, ACE2, ABO, APOE, CFH, CRP, HLA-A, LZTFL1, OAS3, SLC6A20, and TP53, to disclose the puzzled architecture of the interindividual antibody response and dynamics to anti-SARS-CoV-2 vaccine.
Materials and methods
Study design
A retrospective study aimed at assessing the immune response in healthy adult volunteers (n = 195) belonging to the staff of the University-Hospital of Ferrara vaccinated with two doses of anti-SARS-CoV-2 vaccine: Pfizer–BioNTech/BNT162b2 (mRNA-based vaccine) or adenovector, ChAdOx1/AstraZeneca (AdV-based vaccine), in the period starting from January 2021, according to the directives of the Italian Health Ministry. A 6-month period from the second vaccine dose (15–190 days) was considered adequate to monitor the entire circulating antibody kinetic as recently reported by Salvagno et al. in healthcare workers (Salvagno et al., 2022b). The study involving human participants was reviewed and approved by the local regional ethical committee (CE-AVEC; 405/2020/Oss/UniFe); the participants provided their written informed consent to participate in this study. The study consisted of a first immunological assessment aimed at evaluating the anti-SARS-CoV-2 circulating antibodies (IgG and neutralizing antibodies, NAbs) and a genotyping profile investigation including a group of selected common gene variants to identify candidate genetic modifiers of the vaccine-induced immune response.
Blood sample and DNA extraction
Whole blood samples were collected in specific vacutainers, and plasma samples were processed within 1 h from drawing blood, and they were immediately frozen at -80 °C in multiple aliquots and blind tested. DNA was isolated from frozen whole blood by using an automated DNA extraction and purification robot (BioRobot EZ1 system, Qiagen; Hilden, Germany).
Antibody assay
Anti-SARS-CoV-2 immunoglobulin-G (IgG) levels were assessed in duplicate by using Human SARS-CoV-2 Spike (Trimer) IgG-ELISA kit (Invitrogen, Thermo Fisher Scientific, United States) and neutralizing antibodies (NAbs) by SARS-CoV-2 neutralizing antibody ELISA kit (Invitrogen, Thermo Fisher Scientific, United States) in the previously frozen plasma samples, following the manufacturers’ instructions by using a Tecan Infinite MPlex reader (Tecan Trading AG, Switzerland), as previously described (Agostinis et al., 2012).
PCR and pyrosequencing
PCR detection of the selected gene variants was as follows: ACE2 (rs2285666; G>A), HLA-A (rs2571381; T>C, and rs2499; T>G), LZTFL1 (rs35044562; A>G), OAS3 (rs10735079; G>A), and TP53 (P72R, rs1042522; G>C), by rhAmp SNP genotyping technology (IDT, Integrated DNA Technologies, Coralville, IA, United States) on the QuantStudio3 Real-Time PCR System (Thermo Fisher Scientific, United States), according to the supplier’s instructions; ABO (rs657152; A>C), APOE (R158C, rs7412; C>T, and C112R, rs429358; T>C), CFH (Y402H, rs1061170; C>T), CRP (rs2808635; G>T; rs876538; T>C), and SLC6A20 (rs11385942; -/A), by pyrosequencing (Pyromark ID System, Biotage, AB, Uppsala, Sweden) after standard PCR on Agilent SureCycler 8800 (Agilent Technologies, Santa Clara, CA, United States); and ACE1 (ins/del, rs4646994), by agarose gel electrophoresis (I/D allele: 490bp/190bp). DNA samples with known genotypes were used as internal control references for all the sequencing, and a random number of samples (15% for each genotype) were reanalyzed as the internal quality control procedure as previously described (Gemmati et al., 2004; Zamboni and Gemmati, 2007).
Statistical analysis
Statistical analyses were performed using GraphPad Prism version 8 (GraphPad Software, Inc.) and MedCalc version 20.113 (MedCalc Software Ltd.). Regression analyses and scatter diagrams accounted for antibody levels and genetic stratification analysis, respectively; p-values for slope and regression equations accounted for the comparison differences. Moreover, Welch’s t-test for antibody levels comparison, Chi-square test for genotype distributions or allele frequency, and Hardy–Weinberg equilibrium to check possible deviation of genotype/allele distribution were also used according to the recent study of Colucci et al. (2021). The Mann–Whitney non-parametric paired test was performed in case of multiple blood samples from the same donors over time, to evaluate the changes in antibody levels. Spearman’s test was used to assess correlation analyses. p-values were two-sided with a threshold for statistical significance fixed to p ≤ 0.05.
Results
Study population
We recruited 195 healthy subjects among the worker staff at the University-Hospital of Ferrara starting from January 2021 at least 15 days after the second vaccine dose (15–190 days) and before the booster dose of anti-SARS-CoV-2 vaccine. Table 1 shows the epidemiological characteristics of the vaccinated subjects stratified by the two main vaccine formulations received: mRNA-based vaccine, n = 128 (Pfizer–BioNTech) and AdV-based vaccine n = 67 (AstraZeneca). Moreover, a subgroup of 52 subjects (Pfizer–BioNTech, n = 32 and AstraZeneca, n = 20) recruited among those who had the blood sampling within the first 3 months from the second dose of vaccine, kindly agreed to undergo a second blood drawn scheduled about 90 days from their first sampling and they were separately computed for coupled analyses.
TABLE 1. Baseline characteristics of the subjects involved in the study.
Antibody levels (IgG and NAbs)
Figure 1 shows the levels distribution and dynamic of anti-SARS-CoV-2 IgG (Figure 1A) and NAbs (Figure 1B) in the whole cohort of vaccinated subjects investigated over a period of 15–190 days after the second dose of vaccine. We next split subjects into those whose blood was collected within 90 days (T1; n = 104) and after 90 days (T2; n = 91) from the second dose of anti-SARS-CoV-2 vaccine. As expected, a significant decrease of antibody titer has been observed during the time considered. Regardless the time frame between the day of vaccination and the day of blood sampling, IgG and NAbs showed a wide range in levels, suggesting a strong interindividual variability. Globally, a significant fall in levels within the T1–T2 period was observed (IgG, r2 = 0.171; p < 0.0001 and NAbs, r2 = 0.214; p < 0.0001). By comparing the two vaccine types, mRNA-based vaccine vs. AdV-based vaccine, different distributions for both IgG and NAbs were evident (Figures 1C,D). The T1–T2 median time of blood sampling among the mRNA- and AdV-vaccinated subjects was 91.0 days (min 15 days and max 190 days) and 89.5 days (min 15 days and max 185 days), respectively. For a more appropriate evaluation, the two vaccine formulations were also separately assessed for T1 and T2 time points. Accordingly, the specific kinetics and dynamic distributions for IgG levels were characterized by r2 = 0.341, p < 0.001 and r2 = 0.153, p < 0.001, while for NAbs levels by r2 = 0.322, p < 0.0001 and r2 = 0.184, p < 0.0001, respectively, for the mRNA-based vaccine and AdV-based vaccine (Figures 1C,D).
FIGURE 1. Scatter plots of the distribution of antibody levels after second dose of vaccine. (A,B) IgG and NAbs kinetic distributions in the whole cohort of subjects. (C,D) IgG and NAbs distribution stratified by the vaccine type (mRNA-based vaccine: blue dots/line and vector-based vaccine: red dots/line). (E,F) IgG and NAbs distribution, respectively, stratified by T1 and T2 recruitment time (T1: blue dots/line and T2: red dots/line). Each panel shows the specific regression line and the r-coefficient.
By considering a sex stratification, the same analysis showed slightly higher antibody distributions in females than in males for both vaccine formulations (IgG, p = 0.03 and NAbs, p = 0.06 in the whole cohort; IgG, p = 0.10 and NAbs, p = 0.24 in the mRNA vaccine subgroup, and IgG, p = 0.42 and NAbs, p = 0.58 in the AdV vaccine subgroup).
Moreover, as shown in Figures 1E,F, a significant fall in antibody levels was more evident after 3 months from vaccination, being negligible during the first 90 days (T1 blue dots) compared to that observed after 90 days (T2, red dots): IgG: T1, r2 = 0.001; p = 0.74 and T2, r2 = 0.305; p < 0.0001; NAbs: T1, r2 = 0.012; p = 0.09 and T2, r2 = 0.242; p = 0.001).
A strong correlation between IgG and NAbs levels was found in the whole cohort of samples during the whole frame of time considered (Figure 2A; r2 = 0.807; p < 0.0001) and within the T1 and T2 time points (Figure 2B; r2 = 0.789; p = 0.0001 and r2 = 0.764; p = 0.0001, respectively). Finally, even though the strong difference in the antibodies level distribution detected between the two vaccine types, both vaccines maintained a strong IgG/NAbs correlation when separately assessed (Figure 2C; r2 = 0.842; p = 0.0001 and r2 = 0.780; p = 0.0001 for mRNA- and AdV-based vaccine, respectively).
FIGURE 2. IgG and NAbs correlation analysis. (A) Correlation between IgG and NAbs level distributions in the whole cohort of subjects. (B) Correlation between IgG and NAbs stratified by T1 and T2 (T1: blue dots/line and T2: red dots/line). (C) Correlation between IgG and NAbs stratified by vaccine formulation (mRNA-based vaccine: blue dots/line and vector-based vaccine: red dots/line). Each panel shows the specific regression line and r-coefficient.
The wide variability observed in the antibodies distribution, especially due to the different vaccine type, sex, and the timing of blood sampling, prompted the researchers to perform a detailed assessment of the antibody levels by a pairwise analysis of IgG and NAbs in 52 vaccinated subjects (32 mRNA-based vaccine, ♀46.5%, and 20 AdV-based vaccine, ♀45%), who agreed to have double sampling. The T1–T2 median time of the blood sampling in the 52 subjects was comparable to that of the 195 subjects of the whole cohort in order to lessen the time-dependent effects on the antibody’s level (89.5 days, min 15 days and max 187 days, and 88.0 days, min 15 days and max 185 days) for mRNA- and AdV-vaccinated subjects, respectively. Supplementary Figures S1A-D show the comparison of IgG and NAbs levels analyzed for T1 and T2 in the two vaccine formulations, respectively, for the whole cohort of 195 subjects and the subgroup of 52 blood-paired sampling. The sub-analysis in the 52-paired subjects confirmed what is being observed in the whole cohort, with no significant differences in the levels of antibodies comparing the two cohorts matched by the vaccine type.
Gene variants and antibody distribution
We further explored the impact of common genetic variants within genes previously investigated as promising modifiers of the clinical phenotype and progression of COVID-19 or as possible influencers of the antibody response after SARS-CoV-2 infection or vaccination. Some of them have also been investigated in the variability of humoral and cell-mediated response or as causative of specific mild/common adverse reactions, particularly among Pfizer/BNT162b2 vaccinated subjects. Table 2 shows in detail the list of the selected gene variants investigated in our cohort of subjects.
TABLE 2. Panel of gene and variants investigated.
To investigate possible associations between genetic determinants and anti-SARS-CoV-2 vaccine-induced antibody levels, the following statistical evaluations have been performed in the whole cohort and separately for the two vaccine types (when not shown, data comparison did not reach statistical significant values or trends in any of the subgroups considered).
First, by regression analysis we performed the trend estimation of IgG and NAbs levels in the T1–T2 period of 190 days, accounting for the different subgroups of genotypes. The intercepts and slopes comparison of the regression equations of the different subgroups of genotypes yielded significant parameters computing the following gene variants or haplotypes: TP53 (rs1042522), ABO (rs657152), APOE (rs7412/rs429358), ACE2 (rs2285666), HLA-A (rs2571381/rs2499), and CRP rs2808635/rs876538. The related findings are shown in Figures 3A–F, Figures 4A–D, Figures 5A–D, and Supplementary Figures S2. We describe in detail the statistical parameters for regression analyses in the whole group as follows, and Table 3 (whole cohort) and Supplementary Table S1 (mRNA-subgroup) list the complete regression parameters of each comparison (IgG and NAbs).
FIGURE 3. Scatter plots and regression analyses of IgG and NAbs levels stratified by the selected gene variants in the whole cohort. (A,C,E) IgG and (B,D,F) NAbs kinetics according to TP53 (rs1042522), ABO (rs657152), and APOE (rs7412/rs429358). Each panel shows the specific regression lines, according to the indicated genotype/haplotype.
FIGURE 4. Scatter plots and regression analyses of IgG and NAbs levels stratified by ACE2 (rs2285666) genotypes. (A,C) IgG and (B,D) NAbs kinetics in the whole cohort (A,B) and in the mRNA-based vaccine subgroup (C,D). Each panel shows the specific regression lines, according to the indicated genotype.
FIGURE 5. Scatter plots and regression analyses of IgG and NAbs levels stratified by HLA-A (rs2571381/rs2499). (A,C) IgG and (B,D) NAbs kinetics in the whole cohort (A,B) and in the mRNA-based vaccine subgroup (C,D). Each panel shows the specific regression lines according to the indicated haplotype.
TABLE 3. Regression analysis and genotype comparison of the selected SNPs in the whole cohort.
TP53 rs1042522 GG-genotype yielded significant regression equations clustered in the lowest part of the scattering when compared with the regression obtained by the rest of genotypes (CC + CG), accounting for significant different intercepts comparisons (p = 0.03 and p = 0.04, respectively, for IgG and NAbs distributions in the whole cohort, and p = 0.035 and p = 0.05, respectively, for IgG and NAbs distributions in the mRNA vaccine subgroup).
ABO rs657152 TT-genotype yielded significant regression equations clustered in the highest part of the scattering when compared with the regression obtained by the rest of genotypes (GG + GT), accounting for different intercepts comparisons (p = 0.01 and p = 0.03, respectively, for IgG and NAbs distributions in the whole cohort, and p = 0.03 and p = 0.04, respectively, for IgG and NAbs distributions in the mRNA vaccine subgroup).
APOE rs7412/rs429358 Ɛ4-carrying haplotypes (i.e., Ɛ2Ɛ4 and Ɛ3Ɛ4) yielded significant regression equations clustered in the lowest part of the scattering when compared with the regression obtained by the rest of Ɛ3-carrying haplotypes (i.e., Ɛ2Ɛ3 + Ɛ3Ɛ3), accounting for different intercept comparisons (p = 0.03 and p = 0.002, respectively, for IgG and NAbs distributions in the whole cohort, and p = 0.0018 and p = 0.0002, respectively, for IgG and NAbs distributions in the mRNA vaccine subgroup).
ACE2 rs2285666 GA-variant locates on the X chromosome, for this reason it was first analyzed in females and males separately, and both regressions yielded lower trends in the presence of the G-allele (i.e., GG-homozygous females or G-hemizygous males) when compared with the regressions obtained by the respective rest of genotypes (i.e., AA-homozygous in addition to AG-heterozygous females or A-hemizygous males). Accordingly, the combined sex analysis accounted for different intercept comparisons (p = 0.09 and p = 0.075, respectively, for IgG and NAbs distributions in the whole cohort) becoming significant (p = 0.014 and p = 0.009, respectively, for IgG and NAbs distributions) in the mRNA vaccine subgroup. In the sex-specific sub-analysis, only males maintained a significant intercept comparison within the mRNA vaccine subgroup (p = 0.01 and p = 0.0012, respectively, for IgG and NAbs distributions).
HLA-A CT/GT variants (rs2571381/rs2499) are in strong/complete linkage disequilibrium, according to the different population analyzed (e.g., r2 = 0.225244, D’ = 0.999985 and r2 = 0.922353, D’ = 1.000000 from 1000GENOMES:phase_3 LWK and CEU for African ancestry and European ancestry origin, respectively). Accordingly, the haplotypes selected by the researchers to be compared have been the common double homozygotes (i.e., CC/GG) versus those carrying at least two variant alleles (e.g., CT/GT, TT/GG, CC/TT, CT/TT, and TT/GT). Due to the disequilibrium and the low frequency of double carriers, we found only double heterozygotes (i.e., CT/GT) and very few CT/TT. Regressions analysis yielded significant trends among double carriers (CT/GT + CT/TT) that clustered in the lower part of the scattering when compared with the regression obtained by the double homozygotes (CC/GG), accounting for significant different intercept comparisons only in the mRNA-based vaccine subgroup (p = 0.02 and p = 0.03, respectively, for IgG and NAbs distributions).
CRP TG/CT variants (rs2808635/rs876538) are in strong linkage disequilibrium (e.g., r2 = 0.618598, D’ = 0.965325 from 1000GENOMES:phase_3 CEU for European ancestry). Accordingly, the haplotypes selected by the researchers to be compared have been the common double homozygotes (i.e., TT/CC) versus those being at least polymorphic homozygous in any of the two loci (i.e., GG rs2808635 plus any combination of rs876538 or TT rs876538 plus any combination of rs2808635). Regressions analysis yielded significant trends among the polymorphic homozygotes (GG rs2808635/any rs876538 or TT rs876538/any rs2808635) that clustered in the higher part of scattering (Supplementary Figures S2A,B) when compared with regression obtained by the common double homozygotes (TT/CC), accounting for significant different intercept comparisons only in the mRNA-based subgroup for the IgG distribution (p = 0.01 and p = 0.09, respectively, for IgG and NAbs distributions).
Second, as a result of the significant difference in antibody levels between the two vaccine types and of the expected time-dependent fall in levels, we specifically analyzed the mean levels of IgG and NAbs in the T1 and T2 time points in the two cohorts of vaccines stratified by the different genotypes of the specific polymorphisms investigated. The antibody levels stratified by the genotypes gave appreciable and consistent differences in both subgroups of vaccine formulations, and among the AdV-based vaccine, the antibody mean gaps were even higher and retained the same trends observed among the mRNA-based vaccine though the lower number of cases in the former did not allow reaching the statistical significance. For this reason, we report the main findings among the mRNA vaccine subgroup as follows, and Table 4 shows in detail the mean p-values of the whole comparisons.
TABLE 4. IgG and NAbs levels stratified by vaccine formulation and by genotypes.
Tp53 rs1042522 GG-genotype yielded significant or borderline (IgG, p = 0.03 and NAbs, p = 0.09) lower mean antibody levels at T1 (IgG: 19542.6 AU/mL and NAbs: 37.2 AU%) when compared with the counterpart genotypes who instead showed overlapped antibody mean levels (IgG: CC 33445.9 AU/mL and CG 34205.5 AU/mL and NAbs: CC 50.56 AU% and CG 48.14 AU%). The mean levels comparison at T2 maintained similar trends characterized by borderline significant values.
ABO rs657152 TT-genotype yielded appreciable or borderline (IgG, p = 0.06 and NAbs, p = 0.1) higher mean antibody levels at T1 (IgG: 42766.0 AU/mL and NAbs: 59.0 AU%) when compared with the counterpart genotypes (IgG: GG 33099.8 AU/mL and GT 29217.9 AU/mL and NAbs: GG 49.9 AU% and GT 44.55 AU%). The mean levels comparison at T2 was far from being statistically significant.
APOE rs7412/rs429358 Ɛ4-carrying haplotypes (i.e., Ɛ2Ɛ4 or Ɛ3Ɛ4) yielded significant (IgG, p = 0.02 and NAbs, p = 0.0005) lower mean antibody levels at T1 (IgG: 22199.4 AU/mL and NAbs: 28.3 AU%) when compared with the counterpart haplotypes (i.e., Ɛ2Ɛ3 + Ɛ3Ɛ3), overall including those who did not carry the ε4-allele (IgG: 34546.8 AU/mL and NAbs: 53.0 AU%). The mean levels comparison at T2 maintained similar trends far from being statistically significant.
ACE2 rs2285666 GA-variant, for the aforementioned reason (i.e., X-chromosome location), was first analyzed in females and males separately, and in both sexes the presence of the G-allele (i.e., GG-homozygous females or G-hemizygous males) gave lower mean antibody levels in T1 when compared with the opposite respective genotypes (i.e., AA-homozygous females or A-hemizygous males). In detail, the combined sex analysis yielded significant differences, though non-significant for IgG (p = 0.08) and significant for NAbs (p = 0.01), by genotypes comparison (IgG: GG-females + G-males 23804.6 AU/mL and NAbs: GG-females + G-males 40.6 AU%) being lower than (IgG: AA-females + A-males 42068.5 AU/ml and NAbs: AA-females + A-males 75.2 AU%). Accordingly, AG-heterozygous females showed intermediate antibody levels between those shown by the AA- and GG-homozygous females. Interestingly, the same comparisons maintained comparable trends and also significant p-values in T2.
HLA-A CT/GT-haplotype (rs2571381/rs2499), for the aforementioned reason (i.e., strong linkage disequilibrium), was compared against the common double homozygote (i.e., CC/GG) yielding significant (IgG, p = 0.01; NAbs, p = 0.015) lower mean antibody levels at T1 (IgG: 19186.7 AU/mL and NAbs: 31.6 AU%) when compared with the counterpart CC/GG-haplotype (IgG: 34249.0 AU/mL and 50.6 AU%). The mean level comparisons in T2 maintained similar trends though characterized by not significant values.
CRP TG/CT variants (rs2808635/rs876538), for the aforementioned reason (i.e., strong linkage disequilibrium) were compared selecting the common double homozygotes (i.e., TT/CC) vs. those being at least polymorphic homozygous in any of the two loci (i.e., GG rs2808635 plus any combination of rs876538 or TT rs876538 plus any combination of rs2808635). The common TT/CC haplotype yielded significant (IgG, p = 0.025; NAbs, p =0.045) lower mean antibody levels at T1 only in the mRNA vaccine subgroup (IgG: 27992.7 AU/mL and NAbs: 45.0 AU%) when compared with the counterpart GG/TT rare haplotype defined as earlier (IgG: 37339.6 and NAbs: 53.0 AU%); for these reasons, we report only these data here.
Finally, we computed how the genotype, haplotype, or allele frequency of the investigated gene variants stratified above and below the trend lines of the IgG or NAbs distribution in the whole cohort (Table 5) and in the mRNA vaccine subgroup (Supplementary Table S2) and also tested for possible deviation from the expected Hardy–Weinberg equilibrium.
TABLE 5. Genotype, allele or haplotype frequency of the selected gene variants in the area below and above the trend lines of the IgG and NAbs distribution in the whole cohort.
TP53 rs1042522 genotypes were differently distributed in the area above and below the trend lines of the antibodies level scattering of the whole group (p = 0.018 and p = 0.009 for IgG and NAbs, respectively). The homozygous GG-genotype clustered in the area below the trend lines, accordingly GG-genotypes (p = 0.0064 and p = 0.004 for IgG and NAbs, respectively) and G-allele (p = 0.0054 and p = 0.0017 for IgG and NAbs, respectively) were significantly overrepresented when compared with the rest of genotypes (i.e., CC + CG) and the counterpart C-allele, respectively.
ABO rs657152 genotypes were differently distributed in the area above and below the trend lines of the antibodies level scattering of the whole group (p = 0.05 and p = 0.02 for IgG and NAbs, respectively). The homozygous TT-genotype clustered in the area above the trend lines, accordingly TT-genotypes (p = 0.02 and p = 0.06 for IgG and NAbs, respectively) were overrepresented when compared with the rest of genotypes (i.e., GG + GT).
APOE rs7412/rs429358 haplotype distribution gave an overrepresentation of the Ɛ4-allele in the area below the trend lines of the antibodies level scattering of the whole group, accordingly Ɛ4-carrying haplotypes (i.e., Ɛ2Ɛ4 + Ɛ3Ɛ4) were significantly overrepresented (p = 0.005 and p = 0.0025 for IgG and NAbs, respectively) when compared with the counterpart haplotypes (i.e., Ɛ2Ɛ3 + Ɛ3Ɛ3) overall those who did not carry the Ɛ4-allele.
ACE2 rs2285666 GA-variant, as reported previously, was first analyzed in females and males separately, and in both the sexes, the presence of the G-allele (i.e., GG-homozygous females or G-hemizygous males) was overrepresented in the area below the trend lines of the antibodies level scattering (p = 0.03 and p = 0.016 for IgG and NAbs, respectively) when compared with the counterpart A-allele (i.e., AA-homozygous in addition to AG-heterozygous females or A-hemizygous males).
HLA-A CT/GT-haplotype (rs2571381/rs2499), for the aforementioned reason (i.e., strong linkage disequilibrium), was compared against the common double homozygotes (i.e., CC/GG) clustering in the area below the trend lines of the antibodies level scattering of the whole group, accordingly CT/GT-haplotype was significantly overrepresented in the area below the trend lines when compared with the counterpart CC/GG-haplotype (p = 0.04 and p = 0.03 for IgG and NAbs, respectively).
CRP TG/CT variants (rs2808635/rs876538), for the aforementioned reason (i.e., strong linkage disequilibrium), were assessed by comparing the common double homozygotes (i.e., TT/CC) vs. those being at least polymorphic homozygous in any of the two loci (as detailed previously). Nonetheless, the overrepresentation of the TT/CC haplotype in the area below the trend lines (i.e., IgG: 68.2% vs. 31.8%; p = 0.119 and NAbs: 58.2% vs. 41.8%; p = 0.39) as a consequence of the low number of cases in the rare haplotype, the comparison did not reach the statistical significance. For this reason, we report only these data here.
Discussion
Dissecting the intricate features of complex phenotypes (Singh et al., 2010; Gemmati et al., 2012; Gemmati and Tisato, 2023), as COVID-19 and the vaccine-induced immune response in different individuals, is an intricate task, being the result of genome–phenome combinations, and a specific OMICs (COVIDomics) approach is the best way to face them (Gemmati et al., 2019; Gemmati et al., 2020; Gemmati and Tisato, 2020; Milani et al., 2022). During COVID-19 pandemic and anti-SARS-CoV-2 vaccination campaign, unpredictable extreme phenotypes (clinical and laboratory) have been observed, suggesting that inherited traits and predispositions at an individual or population level may account for such manifestations. Our study is within the frame of an ongoing project belonging to the COVID-19 Host Genetic Initiative (HGI) named “Extreme genotype comparison and extreme clinical phenotype comparison in CoV-2 patients: direct candidate genes–pathways and GWAS” (https://www.covid19hg.org/partners/?partner=rec0CufBJdOxdaest). COVID-19 pandemic continues to be a major public health threat, especially in countries with low vaccination rates, and novel genes responsible for different susceptibility and severity have been identified by the HGI consortium (COVID-19 Host Genetics Initiative, 2022). Accordingly, extreme phenotypes (i.e., individual vaccine-induced immune response) after anti-SARS-CoV-2 vaccine have been deeply investigated in this study. Anti-viral antibody response is a measure of reactivity against the viral peptides, and this has mixed genetic and environmental contributions according to the individual genome landscape (Jonsson et al., 2017; Chung et al., 2019; Castro Dopico and Karlsson Hedestam, 2022). Genetic loci with a role in determining the extent of the humoral response to SARS-CoV-2 infection or vaccination or associated to post-vaccine side effects have not been completely recognized.
The main result of this study is that the extent of anti-SARS-CoV-2 IgG and NAbs response broadly varied among individuals and vaccine types, and their levels decreased significantly starting from 3 months after the second dose of vaccine (i.e., T2). The most common characteristics investigated, such as age, sex, ABO blood group, did not completely account for the interindividual variability. Accordingly, male sex had slightly lower antibody levels (IgG and NAbs) compared to females in the first trimester after the second vaccine dose (i.e., T1), and age indistinctly affected antibody levels in males and females though male sex reached lower values at T2, particularly because they are already presented with lower mean antibody levels at T1. Finally, ABO blood group gave inconsistent results, ascribing to the O group and A group overlapping regression trends and mean levels in-between the highest B group and the lowest AB group in both vaccine types as recently described in healthy donor populations (Nunhofer et al., 2022).
The observation that during the first 90 days after the second dose (T1) the antibody levels did not change significantly among the cases of our cohort (r2 = 0.001 and r2 = 0.012, respectively, for IgG and NAbs), made us confident in having performed the core comparisons and conclusions, especially in the T1 frame considering T2 distinctive of the antibody drop (r2 = 0.305 and r2 = 0.242, respectively, for IgG and NAbs).
The strong correlation found between IgG and NAbs levels in both the vaccines during the 6 months of time considered has been useful to analyze data considering the whole cohort and separately by the vaccine type. The need of disaggregated data analyses is of significant importance to understand the extent to which any considered variable is influencing the outcome, genetic characteristics and sex included. In addition, the apparently arbitrary setting of the 6 months’ frame (i.e., T1 and T2) to evaluate the antibody levels could be considered unfaithful to identify individual genetic predispositions of the vaccine immune response; this is because genotypes could randomly cluster in a particular percentile in that specific time regardless its own role of “genetic predictor” if any. For this reason, about 50% of those subjects who belonged to the T1 window randomly underwent a second blood drawing so that their T1–T2 mean time matched that of the whole cohort, and data obtained did not show significant differences when compared with that of the whole cohort.
To ascribe to specific genes and genotypes potential effects in predicting IgG and NAbs titer or duration, we focused on common gene variants previously investigated as modifiers of the clinical phenotype of COVID-19 or as possible influencers of antibody response after SARS-CoV-2 infection or vaccine.
First, we found that at longer time frames (i.e., T2), genetics lost its role in determining the antibody level after vaccination and regardless the gene variant considered antibody put down definitely around 150 days after vaccination, but this does not mean that after this time the protection from SARS-CoV-2 is lost.
Conversely, significant differences emerged in the first months of observation (i.e., T1) in six genes carrying the following gene variants [i.e., TP53 (rs1042522), ABO (rs657152), APOE (rs7412/rs429358), ACE2 (rs2285666), HLA-A (rs2571381/rs2499), and CRP (rs2808635/rs876538)]. The three complementary approaches we applied to assess possible genetically driven differences in antibody levels after vaccination gave consistent results ascribing to a definite genotype or haplotype different trend estimations (i.e., regression analysis), mean levels stratified by genotypes (i.e., Welch’s t-test), and genotype or allele frequencies stratified by antibody distributions including possible deviation from the Hardy–Weinberg equilibrium. Although a multiple test comparison should be considered and a larger cohort of subjects recruited, due to the high number of gene variants considered in a trait that appears to be polygenic, we went for a single-analysis approach considering the explorative purpose of the investigation. The selected approach may represent a limitation of the study deserving further investigations in larger cohorts.
TP53 rs1042522 is the most common TP53 SNP occurring at codon P72R (CCC>CGC) (Matlashewski et al., 1987). Large part of the investigations addressed its role on cancer (De Souza et al., 2021), but interesting correlations on metabolic pathways, immune cells, and immune response have been recently suggested in COVID-19 (Milani et al., 2022). Our results are in line with the differentially expressed role of P72- and R72-allele on immune checkpoint inhibitors and on the association of the G-allele or GG-genotype with imbalanced immune regulation (Dolgikh et al., 2022) and immune dysfunction, ascribing to TP53 gene a crucial role in the immune response. In general, R72 carrying population is more susceptible to viral infection because the pro-inflammatory immune response is not as strong as the P72 population (Lodhi et al., 2021). Indeed, the P72 allele is overrepresented among the equatorial populations, where the immune challenge is greater. The innate immune response significantly differs in the P72 and R72 alleles, and this difference is regulated by a subset of genes known to control the inflammatory response, including some of the NFκB target genes, which are better activated by the P72 variant (Lodhi et al., 2021).
ABO (rs657152) was first described as a genetic susceptibility locus in patients with COVID-19 with respiratory failure suggesting a potential involvement also for the ABO blood group system, showing a higher risk in blood group A and a protective effect in blood group O as compared with other groups (Severe Covid et al., 2020; Venkataraman et al., 2022). Contextually, rs657152 was hypothesized as a determinant in COVID-19 prognosis and severity together with the genes of the RAS-pathway (Gemmati et al., 2020; Gemmati and Tisato, 2020). After vaccination programs, ABO blood group was investigated as potentially involved in the quality of the specific immune response against SARS-CoV-2 infection or vaccination, though no definite results have been obtained as discussed previously (Gil-Manso et al., 2022; Nunhofer et al., 2022; Sgherza et al., 2022). The association we found with the rs657152 variant strongly recalls the previous described connotations, though further investigations matched with the phenotypic ABO blood group are needed.
APOE Ɛ4 haplotype (rs7412/rs429358) has been shown to associate with increased susceptibility to SARS-CoV-2 infection and COVID-19 mortality in previous genetic studies (Kuo et al., 2020; Kurki et al., 2021). APOE Ɛ4 allele has been particularly investigated in neurodegenerative and cognitive impairment diseases, and indeed, severe COVID-19 and Alzheimer’s disease share common genetic and metabolic pathways via the OAS1 gene and many risk factors strongly overlap (Tisato et al., 2018; Gemmati et al., 2020; Magusali et al., 2021). Intriguingly, several circulating serum proteins (e.g., ApoE) may adsorb specific SARS-CoV-2 peptides and alter their infection ability by mediating the ACE2 receptor. On the other hand, SARS-CoV-2 RBD can cause significant structural changes in ApoE lipoprotein and by hijacking the metabolic lipoprotein pathway it may facilitate cell entry (Yin et al., 2021). Contextually, the binding drastically changes the RBD epitope and such alterations can theoretically substitute ACE2 for ApoE in SARS-COV-2 interaction. Such interaction is further demonstrated by the decreased levels of LDL (ApoE included) during SARS-CoV-2 infection, ascribing to APOE Ɛ2-allele protective effects (Espinosa-Salinas et al., 2022). As there exist different affinities of the different APOE isoforms for the LDL receptors, similar changes and modifications might alike be responsible for different RBD-APOE interactions, affecting in turn the global immune response.
ACE2 (rs2285666) is one of the most relevant SNPs in the gene influencing activity and levels of the receptor, it is a transition G8790A with the GG-genotype characterized by about 50% expression reduction compared to AA-genotype. ACE2 is the SARS-CoV-2 receptor cell entry by interacting with the RBD SARS-CoV-2 spike proteins. Differently expressed ACE2 levels and/or qualitative variants affecting the spike-receptor interaction have been investigated as partially responsible for the wide range of COVID-19 symptoms. ACE2 gene is located on the X chromosome, accordingly males and females are necessarily characterized by different genotype architectures ascribing to male only the G- or A-hemizygous condition compared to female having the chance to carry the AG-heterozygosity. This suggested possible benefits for females against the risk of infection considering them more protected than males considered the at-risk sex (Gemmati et al., 2020). After vaccination, the neo-synthetized spike proteins migrate to the cell surface and protrude with a comparable native-like conformation to be recognized by the immune system and start the immune response. On the other hand, the vaccine-induced spike proteins can react with ACE2 receptors of the several neighboring cells, including platelets, and after cells death several free-floating spike proteins circulate in the blood and systematically interact with ACE2 receptors expressed by a variety of cells promoting complex internalization and degradation (Angeli et al., 2021; Watanabe et al., 2021). This proposed mechanism, together with the genetically driven ACE2 receptor expression (i.e., rs2285666), follows the same affinity interactions of the virus spike proteins. This may potentially influence the individual immune response and/or the epitope availability being responsible for the wide range of antibody levels detected after vaccination (Angeli et al., 2021; Watanabe et al., 2021).
HLA-A (rs2571381/rs2499) are the leading common SNPs at HLA-A*03:01 locus and have been recently associated to reactogenicity after Pfizer–BioNtech vaccine with the strongest signal association (p = 1.16E-12 and p = 1.72E-16 respectively) after GWAS extreme phenotype comparisons (Bolze et al., 2022). Although it is not certain whether A*03:01 is the causative allele due to linkage disequilibrium, SNPs imputation quality, and different population ancestry, further specific functional investigations are mandatory. Identifying which epitope of the vaccine-derived spike protein exactly interacts with HLA-A*03:01 locus could attribute any causality to the vaccine-associated adverse reactions or to the antibody production extent (Bolze et al., 2022; Venkataraman et al., 2022). We found different regression trends and antibody levels by comparing the double carriers (CT/GT rare alleles) vs. the double homozygotes (CC/GG common alleles), suggesting possible associations not only to reactogenicity but also in predicting the immune response after SARS-CoV-2 vaccination or infection. Due to the low frequency of the rare genotype combinations, we could not investigate the effect in the homozygotes for the rare alleles, the best candidates to show even greater genotype-associated differences, suggesting larger population studies.
CRP (rs2808635/rs876538) variants are associated to basal and stimulated circulating levels of CRP (Brull et al., 2003), and several investigations ascribed to these SNPs prognostic pharmacogenomics information on treatment and drug response (Parmeggiani et al., 2008; Parmeggiani et al., 2011; Moran et al., 2018). These molecular markers may be more than indicators of inflammation, and preferably, they can play a key role also in the response to external infectious stimuli. The interesting association we found with CRP rs2808635/rs876538 variants are in agreement with recent reports on CRP and anti-SARS-CoV-2 IgG levels. Baseline anti-SARS-CoV-2 IgG levels have been associated with CRP levels in COVID-19 patients, suggesting persistence of an unresolved inflammation status (Gianfagna et al., 2022). Moreover, among vaccinated subjects monitored after the second dose, those subjects showing a rapid IgG growth and a better kinetic pattern also had the highest CRP levels, considered as the best performers to anti-SARS-CoV-2 vaccination. The positive association between high CRP levels and post-vaccination IgG response has also been confirmed in a detailed survey of healthcare workers after Pfizer/BioNTech vaccine (Salvagno et al., 2022a). The authors found a highly significant correlation between pre-booster serum CRP levels and post-booster anti-SARS-COV-2 spike trimeric-RBD IgG antibodies. Conversely, no significant association was found between serum pre-booster levels of CRP and anti-SARS-CoV-2 spike trimeric IgG, concluding that the extent of the humoral response after COVID-19 vaccination might be affected by baseline inflammation as assessed by individual CRP levels.
The genetic markers we investigated may predict in advance those subjects at higher risk to be considered lower responders and could be used as valuable indicators for targeted vaccine booster administration and prevention programs. Overall, these results suggest that different genetic backgrounds may influence the vaccine-induced immune response and suggest that host genetics has a key role. Genetic differences, age, sex, and environmental factors can contribute to alter or enforce the immune response as demonstrated for other vaccines (Hamming et al., 2004; Gemmati and Tisato, 2020; Liao et al., 2020; Tay et al., 2020; Castro Dopico and Karlsson Hedestam, 2022; Milani et al., 2022). Combined COVIDomics investigative approaches may identify key determinants of the vaccine-induced immune response, suggesting possible genetic predispositions to develop a more efficient vaccine response. The new insights gained by using polygenic predictive markers merged with demographic data and serological screenings, may increase the understanding of the individual humoral response, accounting for ethnic and geographical differences, in both COVID-19 and anti-SARS-CoV-2 vaccination.
Conclusion
Genetic predispositions and haplotypes fostering favorable immune responses might be helpful to recognize lowly or non-responsive individuals to any vaccination. Future larger and multicenter studies could help in elucidating the underlying inherited causes that could arise from a previous history of similar infections together with individual and environmental factors. The identification of specific genetic signatures in loci useful to track in advance those subjects at risk of a weak antibody response after vaccination or infection has an unquestionable translational power.
Extremely useful future implications of our genetic analyses can also be directed to dissect the molecular mechanisms underlying the individual humoral and cell-mediated variability to the vaccine(s) by in vitro dedicated experiments in relevant cell models. As recently reported, whole blood cultures from subjects after a single dose of anti-SARS-CoV-2 vaccine or after two doses of an adjuvanted SARS-CoV-2 recombinant protein also produced different cytokine secretion profiles (Lineburg et al., 2021; De Rosa et al., 2022). This paves the way to design candidate genomic screening and personalized timing to optimize dose administration in any vaccination program. Moreover, advanced 3D cell models based on integrated circuits and systems as chip technology may help and serve as an additional in vitro surrogate to summarize the complex response of any cell/organ target for the virus recapitulating the native cellular environments under higher biological fidelity (Singh et al., 2022). This will allow a more efficient and personalized approach for vaccine selection, formulation, timing, and scheduling of the booster doses.
Host genetic investigation is extremely useful when designing, selecting for specific peptide vaccines, and checking for their efficacy and protection across populations. However, it is likely that other important genetic influences and heritable traits exist together with acquired and concomitant situations that are responsible for the wide extent of the observed immune response.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by local regional ethical committee (CE-AVEC; 405/2020/Oss/UniFe). The patients/participants provided their written informed consent to participate in this study.
Author contributions
Conception and design of the work: DG and VT. Acquisition of data: GL, IG, JA, SH, AP, PP, SP, AS, and MS. Data analysis and interpretation: DG, VT, GZ, and PS. Manuscript writing: DG and VT All authors have approved the submitted version.
Funding
This study was supported by local FAR and FIR grants from University of Ferrara (to DG and VT)
Acknowledgments
The authors would like to thank the nurses and staff of the Centre Hemostasis and Thrombosis, Hospital University of Ferrara (AM, DD, ES, LL, and CB), for their help in performing blood sample collections.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1028081/full#supplementary-material
References
Agostinis, C., Bulla, R., Tisato, V., De Seta, F., Alberico, S., Secchiero, P., et al. (2012). Soluble TRAIL is elevated in recurrent miscarriage and inhibits the in vitro adhesion and migration of HTR8 trophoblastic cells. Hum. Reprod. 27 (10), 2941–2947. doi:10.1093/humrep/des289
Angeli, F., Spanevello, A., Reboldi, G., Visca, D., and Verdecchia, P. (2021). SARS-CoV-2 vaccines: Lights and shadows. Eur. J. Intern. Med. 88, 1–8. doi:10.1016/j.ejim.2021.04.019
Baden, L. R., El Sahly, H. M., Essink, B., Kotloff, K., Frey, S., Novak, R., et al. (2021). Efficacy and safety of the mRNA-1273 SARS-CoV-2 vaccine. N. Engl. J. Med. Overseas. Ed. 384 (5), 403–416. doi:10.1056/NEJMoa2035389
Bloch, E. M., Patel, E. U., Marshall, C., Littlefield, K., Goel, R., Grossman, B. J., et al. (2021). ABO blood group and SARS-CoV-2 antibody response in a convalescent donor population. Vox Sang. 116 (7), 766–773. doi:10.1111/vox.13070
Bolze, A., Neveux, I., Schiabor Barrett, K. M., White, S., Isaksson, M., Dabe, S., et al. (2022). HLA-A *03:01 is associated with increased risk of fever, chills, and stronger side effects from Pfizer-BioNTech COVID-19 vaccination. HGG Adv. 3 (2), 100084. doi:10.1016/j.xhgg.2021.100084
Brull, D. J., Serrano, N., Zito, F., Jones, L., Montgomery, H. E., Rumley, A., et al. (2003). Human CRP gene polymorphism influences CRP levels: Implications for the prediction and pathogenesis of coronary heart disease. Arterioscler. Thromb. Vasc. Biol. 23 (11), 2063–2069. doi:10.1161/01.ATV.0000084640.21712.9C
Castro Dopico, X., and Karlsson Hedestam, G. B. (2022). A family matter: Anti-viral antibody responses. Immunity 55 (1), 8–10. doi:10.1016/j.immuni.2021.12.008
Chung, S., Roh, E. Y., Park, B., Lee, Y., Shin, S., Yoon, J. H., et al. (2019). GWAS identifying HLA-DPB1 gene variants associated with responsiveness to Hepatitis B virus vaccination in Koreans: Independent association of HLA-DPB1*04:02 possessing rs1042169 G - rs9277355 C - rs9277356 A. J. Viral Hepat. 26 (11), 1318–1329. doi:10.1111/jvh.13168
Colucci, M., De Santis, E., Totti, B., Miroballo, M., Tamiro, F., Rossi, G., et al. (2021). Associations between allelic variants of the human IgH 3' regulatory region 1 and the immune response to BNT162b2 mRNA vaccine. Vaccines (Basel) 9 (10), 1207. doi:10.3390/vaccines9101207
Copley, H. C., Gragert, L., Leach, A. R., and Kosmoliaptsis, V. (2021). Influence of HLA class II polymorphism on predicted cellular immunity against SARS-CoV-2 at the population and individual level. Front. Immunol. 12, 669357. doi:10.3389/fimmu.2021.669357
COVID-19 Host Genetics Initiative (2022). A first update on mapping the human genetic architecture of COVID-19. Nature 608 (7921), E1–E10. doi:10.1038/s41586-022-04826-7
COVID-19 Host Genetics Initiative (2021). Mapping the human genetic architecture of COVID-19. Nature 600 (7889), 472–477. doi:10.1038/s41586-021-03767-x
Crocchiolo, R., Gallina, A. M., Pani, A., Campisi, D., Cento, V., Sacchi, N., et al. (2022). Polymorphism of the HLA system and weak antibody response to BNT162b2 mRNA vaccine. HLA 99 (3), 183–191. doi:10.1111/tan.14546
Cuapio, A., Boulouis, C., Filipovic, I., Wullimann, D., Kammann, T., Parrot, T., et al. (2022). NK cell frequencies, function and correlates to vaccine outcome in BNT162b2 mRNA anti-SARS-CoV-2 vaccinated healthy and immunocompromised individuals. Mol. Med. 28 (1), 20. doi:10.1186/s10020-022-00443-2
de Freitas Dutra, V., Bonet-Bub, C., Yokoyama, A. P. H., Achkar, R., Machado, R. R. G., Assuncao, M., et al. (2021). Anti-A and SARS-CoV-2: An intriguing association. Vox Sang. 116 (5), 557–563. doi:10.1111/vox.13044
De Rosa, S. C., Cohen, K. W., Bonaparte, M., Fu, B., Garg, S., Gerard, C., et al. (2022). Whole-blood cytokine secretion assay as a high-throughput alternative for assessing the cell-mediated immunity profile after two doses of an adjuvanted SARS-CoV-2 recombinant protein vaccine candidate. Clin. Transl. Immunol. 11 (1), e1360. doi:10.1002/cti2.1360
De Souza, C., Madden, J., Koestler, D. C., Minn, D., Montoya, D. J., Minn, K., et al. (2021). Effect of the p53 P72R polymorphism on mutant TP53 allele selection in human cancer. J. Natl. Cancer Inst. 113 (9), 1246–1257. doi:10.1093/jnci/djab019
Dolgikh, O. V., Zaitseva, N. V., Nikonoshina, N. A., and Alekseev, V. B. (2022). Indicators of immune and neurohumoral profile in women of fertile age with functional disorders of the autonomic nervous system Associated with polymorphic variants of the HTR2A (rs7997012) and TP53 (rs1042522) genes. Bull. Exp. Biol. Med. 173 (2), 224–228. doi:10.1007/s10517-022-05523-1
Elberry, M. H., Abdelgawad, H. A. H., Hamdallah, A., Abdella, W. S., Ahmed, A. S., Ghaith, H. S., et al. (2022). A systematic review of vaccine-induced thrombotic thrombocytopenia in individuals who received COVID-19 adenoviral-vector-based vaccines. J. Thromb. Thrombolysis 53 (4), 798–823. doi:10.1007/s11239-021-02626-w
Espinosa-Salinas, I., Colmenarejo, G., Fernandez-Diaz, C. M., Gomez de Cedron, M., Martinez, J. A., Reglero, G., et al. (2022). Potential protective effect against SARS-CoV-2 infection by APOE rs7412 polymorphism. Sci. Rep. 12 (1), 7247. doi:10.1038/s41598-022-10923-4
Falsey, A. R., Sobieszczyk, M. E., Hirsch, I., Sproule, S., Robb, M. L., Corey, L., et al. (2021). Phase 3 safety and efficacy of AZD1222 (ChAdOx1 nCoV-19) covid-19 vaccine. N. Engl. J. Med. 385 (25), 2348–2360. doi:10.1056/NEJMoa2105290
Gemmati, D., Bramanti, B., Serino, M. L., Secchiero, P., Zauli, G., and Tisato, V. (2020). COVID-19 and individual genetic susceptibility/receptivity: Role of ACE1/ACE2 genes, immunity, inflammation and coagulation. Might the double X-chromosome in females Be protective against SARS-CoV-2 compared to the single X-chromosome in males? Int. J. Mol. Sci. 21 (10), E3474. doi:10.3390/ijms21103474
Gemmati, D., and Tisato, V. (2020). Genetic hypothesis and Pharmacogenetics side of renin-angiotensin-system in COVID-19. Genes. (Basel) 11 (9), E1044. doi:10.3390/genes11091044
Gemmati, D., and Tisato, V. (2023). “Genomic and epigenomic signature at the branch-point among genome, phenome and sexome in health and disease: A multi-omics approach,” in Principles of gender-specific medicine: The sex-specific biology in the post-genomic era. Editor M. J. Legato. 4th Edition (New York: Elsevier). ISBN: 9780323885348.
Gemmati, D., Tognazzo, S., Serino, M. L., Fogato, L., Carandina, S., De Palma, M., et al. (2004). Factor XIII V34L polymorphism modulates the risk of chronic venous leg ulcer progression and extension. Wound Repair Regen. 12 (5), 512–517. doi:10.1111/j.1067-1927.2004.012503.x
Gemmati, D., Varani, K., Bramanti, B., Piva, R., Bonaccorsi, G., Trentini, A., et al. (2019). Bridging the gap" everything that could have been avoided if we had applied gender medicine, Pharmacogenetics and personalized medicine in the gender-omics and sex-omics era. Int. J. Mol. Sci. 21 (1), E296. doi:10.3390/ijms21010296
Gemmati, D., Zeri, G., Orioli, E., De Gaetano, F. E., Salvi, F., Bartolomei, I., et al. (2012). Polymorphisms in the genes coding for iron binding and transporting proteins are associated with disability, severity, and early progression in multiple sclerosis. BMC Med. Genet. 13, 70. doi:10.1186/1471-2350-13-70
Gianfagna, F., Veronesi, G., Baj, A., Dalla Gasperina, D., Siclari, S., Drago Ferrante, F., et al. (2022). Anti-SARS-CoV-2 antibody levels and kinetics of vaccine response: Potential role for unresolved inflammation following recovery from SARS-CoV-2 infection. Sci. Rep. 12 (1), 385. doi:10.1038/s41598-021-04344-y
Gil-Manso, S., Miguens Blanco, I., Motyka, B., Halpin, A., Lopez-Esteban, R., Perez-Fernandez, V. A., et al. (2022). ABO blood group is involved in the quality of the specific immune response anti-SARS-CoV-2. Virulence 13 (1), 30–45. doi:10.1080/21505594.2021.2019959
Grifoni, A., Sidney, J., Zhang, Y., Scheuermann, R. H., Peters, B., and Sette, A. (2020). A sequence homology and bioinformatic approach can predict candidate targets for immune responses to SARS-CoV-2. Cell. Host Microbe 27 (4), 671–680. e672. doi:10.1016/j.chom.2020.03.002
Hamming, I., Timens, W., Bulthuis, M. L., Lely, A. T., Navis, G., and van Goor, H. (2004). Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J. Pathol. 203 (2), 631–637. doi:10.1002/path.1570
Huang, C., Wang, Y., Li, X., Ren, L., Zhao, J., Hu, Y., et al. (2020). Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395 (10223), 497–506. doi:10.1016/S0140-6736(20)30183-5
Jonsson, S., Sveinbjornsson, G., de Lapuente Portilla, A. L., Swaminathan, B., Plomp, R., Dekkers, G., et al. (2017). Identification of sequence variants influencing immunoglobulin levels. Nat. Genet. 49 (8), 1182–1191. doi:10.1038/ng.3897
Klok, F. A., Pai, M., Huisman, M. V., and Makris, M. (2022). Vaccine-induced immune thrombotic thrombocytopenia. Lancet. Haematol. 9 (1), e73–e80. doi:10.1016/S2352-3026(21)00306-9
Kuo, C. L., Pilling, L. C., Atkins, J. L., Masoli, J. A. H., Delgado, J., Kuchel, G. A., et al. (2020). APOE e4 Genotype Predicts Severe COVID-19 in the UK Biobank Community Cohort. J. Gerontol. A Biol. Sci. Med. Sci. 75 (11), 2231–2232. doi:10.1093/gerona/glaa131
Kurki, S. N., Kantonen, J., Kaivola, K., Hokkanen, L., Mayranpaa, M. I., Puttonen, H., et al. (2021). APOE ε4 associates with increased risk of severe COVID-19, cerebral microhaemorrhages and post-COVID mental fatigue: A Finnish biobank, autopsy and clinical study. Acta Neuropathol. Commun. 9 (1), 199. doi:10.1186/s40478-021-01302-7
Liao, M., Liu, Y., Yuan, J., Wen, Y., Xu, G., Zhao, J., et al. (2020). Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat. Med. 26 (6), 842–844. doi:10.1038/s41591-020-0901-9
Lineburg, K. E., Neller, M. A., Ambalathingal, G. R., Le Texier, L., Raju, J., Swaminathan, S., et al. (2021). Rapid whole-blood assay to detect SARS-CoV-2-specific memory T-cell immunity following a single dose of AstraZeneca ChAdOx1-S COVID-19 vaccine. Clin. Transl. Immunol. 10 (8), e1326. doi:10.1002/cti2.1326
Lodhi, N., Singh, R., Rajput, S. P., and Saquib, Q. (2021). SARS-CoV-2: Understanding the transcriptional regulation of ACE2 and TMPRSS2 and the role of single nucleotide polymorphism (SNP) at codon 72 of p53 in the innate immune response against virus infection. Int. J. Mol. Sci. 22 (16), 8660. doi:10.3390/ijms22168660
Magusali, N., Graham, A. C., Piers, T. M., Panichnantakul, P., Yaman, U., Shoai, M., et al. (2021). A genetic link between risk for Alzheimer's disease and severe COVID-19 outcomes via the OAS1 gene. Brain. 144 (12), 3727–3741. doi:10.1093/brain/awab337
Matlashewski, G. J., Tuck, S., Pim, D., Lamb, P., Schneider, J., and Crawford, L. V. (1987). Primary structure polymorphism at amino acid residue 72 of human p53. Mol. Cell. Biol. 7 (2), 961–963. doi:10.1128/mcb.7.2.961
Milani, D., Caruso, L., Zauli, E., Al Owaifeer, A. M., Secchiero, P., Zauli, G., et al. (2022). p53/NF-kB balance in SARS-CoV-2 infection: From OMICs, genomics and pharmacogenomics insights to tailored therapeutic perspectives (COVIDomics). Front. Pharmacol. 13, 871583. doi:10.3389/fphar.2022.871583
Moran, C. J., Kaplan, J. L., and Winter, H. S. (2018). Genetic variation affects C-reactive protein elevations in crohn's disease. Inflamm. Bowel Dis. 24 (9), 2048–2052. doi:10.1093/ibd/izy100
Nunhofer, V., Weidner, L., Hoeggerl, A. D., Zimmermann, G., Badstuber, N., Grabmer, C., et al. (2022). Persistence of naturally acquired and functional SARS-CoV-2 antibodies in blood donors one year after infection. Viruses 14 (3), 637. doi:10.3390/v14030637
Parmeggiani, F., Costagliola, C., Gemmati, D., D'Angelo, S., Perri, P., Campa, C., et al. (2008). Coagulation gene predictors of photodynamic therapy for occult choroidal neovascularization in age-related macular degeneration. Invest. Ophthalmol. Vis. Sci. 49 (7), 3100–3106. doi:10.1167/iovs.07-1654
Parmeggiani, F., Gemmati, D., Costagliola, C., Semeraro, F., Perri, P., D'Angelo, S., et al. (2011). Genetic predictors of response to photodynamictherapy. Mol. Diagn. Ther. 15 (4), 195–210. doi:10.1007/BF03256411
Pavord, S., Scully, M., Hunt, B. J., Lester, W., Bagot, C., Craven, B., et al. (2021). Clinical features of vaccine-induced immune thrombocytopenia and thrombosis. N. Engl. J. Med. 385 (18), 1680–1689. doi:10.1056/NEJMoa2109908
Pereira, E., Felipe, S., de Freitas, R., Araujo, V., Soares, P., Ribeiro, J., et al. (2022). ABO blood group and link to COVID-19: A comprehensive review of the reported associations and their possible underlying mechanisms. Microb. Pathog. 169, 105658. doi:10.1016/j.micpath.2022.105658
Polack, F. P., Thomas, S. J., Kitchin, N., Absalon, J., Gurtman, A., Lockhart, S., et al. (2020). Safety and efficacy of the BNT162b2 mRNA covid-19 vaccine. N. Engl. J. Med. 383 (27), 2603–2615. doi:10.1056/NEJMoa2034577
Pretti, M. A. M., Galvani, R. G., Vieira, G. F., Bonomo, A., Bonamino, M. H., and Boroni, M. (2020). Class I HLA allele predicted restricted antigenic coverages for spike and nucleocapsid proteins are associated with deaths related to COVID-19. Front. Immunol. 11, 565730. doi:10.3389/fimmu.2020.565730
Ragone, C., Meola, S., Fiorillo, P. C., Penta, R., Auriemma, L., Tornesello, M. L., et al. (2021). HLA does not impact on short-medium-term antibody response to preventive anti-SARS-cov-2 vaccine. Front. Immunol. 12, 734689. doi:10.3389/fimmu.2021.734689
Salvagno, G. L., Henry, B. M., Pighi, L., De Nitto, S., and Lippi, G. (2022a). Serum C reactive protein predicts humoral response after BNT162b2 booster administration. J. Infect. 85 (1), e24–e25. doi:10.1016/j.jinf.2022.04.015
Salvagno, G. L., Henry, B. M., Pighi, L., de, N. S., and Lippi, G. (2022b). Total anti-SARS-CoV-2 antibodies measured 6 months after Pfizer-BioNTech COVID-19 vaccination in healthcare workers. J. Med. Biochem. 41 (2), 199–203. doi:10.5937/jomb0-33999
Severe Covid, G. G., Ellinghaus, D., Degenhardt, F., Bujanda, L., Buti, M., Albillos, A., et al. (2020). Genomewide association study of severe covid-19 with respiratory failure. N. Engl. J. Med. 383 (16), 1522–1534. doi:10.1056/NEJMoa2020283
Sgherza, N., Zucano, S., Vitucci, A., Palma, A., Campanale, D., Larocca, A. M. V., et al. (2022). Antibody response to BNT162b2 SARS-CoV-2 mRNA vaccine is not influenced by AB0 blood group in subjects with transfusion-dependent thalassemia. Acta Biomed. 93 (2), e2022134. doi:10.23750/abm.v93i2.12819
Shukla, P., Pandey, P., Prasad, B., Robinson, T., Purohit, R., D'Cruz, L. G., et al. (2022). Immuno-informatics analysis predicts B and T cell consensus epitopes for designing peptide vaccine against SARS-CoV-2 with 99.82% global population coverage. Brief. Bioinform. 23 (1), bbab496. doi:10.1093/bib/bbab496
Singh, A. V., Chandrasekar, V., Laux, P., Luch, A., Dakua, S. P., Zamboni, P., et al. (2022). Micropatterned neurovascular interface to mimic the blood-brain barrier's neurophysiology and micromechanical function: A BBB-on-CHIP model. Cells 11 (18), 2801. doi:10.3390/cells11182801
Singh, A. V., Subhashree, L., Milani, P., Gemmati, D., and Zamboni, P. (2010). Interplay of iron metallobiology, metalloproteinases, and FXIII, and role of their gene variants in venous leg ulcer. Int. J. Low. Extrem. Wounds 9 (4), 166–179. doi:10.1177/1534734610384653
Tay, M. Z., Poh, C. M., Renia, L., MacAry, P. A., and Ng, L. F. P. (2020). The trinity of COVID-19: Immunity, inflammation and intervention. Nat. Rev. Immunol. 20 (6), 363–374. doi:10.1038/s41577-020-0311-8
Thomas, S. J., Moreira, E. D., Kitchin, N., Absalon, J., Gurtman, A., Lockhart, S., et al. (2021). Safety and efficacy of the BNT162b2 mRNA covid-19 vaccine through 6 months. N. Engl. J. Med. 385 (19), 1761–1773. doi:10.1056/NEJMoa2110345
Tisato, V., Zuliani, G., Vigliano, M., Longo, G., Franchini, E., Secchiero, P., et al. (2018). Gene-gene interactions among coding genes of iron-homeostasis proteins and APOE-alleles in cognitive impairment diseases. PLoS One 13 (3), e0193867. doi:10.1371/journal.pone.0193867
Tiwari Pandey, A., Pandey, I., Zamboni, P., Gemmati, D., Kanase, A., Singh, A. V., et al. (2020). Traditional herbal remedies with a multifunctional therapeutic approach as an implication in COVID-19 associated Co-infections. Coatings 10 (8), 761. doi:10.3390/coatings10080761
Tregoning, J. S., Brown, E. S., Cheeseman, H. M., Flight, K. E., Higham, S. L., Lemm, N. M., et al. (2020). Vaccines for COVID-19. Clin. Exp. Immunol. 202 (2), 162–192. doi:10.1111/cei.13517
Venkataraman, T., Valencia, C., Mangino, M., Morgenlander, W., Clipman, S. J., Liechti, T., et al. (2022). Analysis of antibody binding specificities in twin and SNP-genotyped cohorts reveals that antiviral antibody epitope selection is a heritable trait. Immunity 55 (1), 174–184.e5. e175. doi:10.1016/j.immuni.2021.12.004
Vicentini, C., Bordino, V., Cornio, A. R., Meddis, D., Ditommaso, S., Giacomuzzi, M., et al. (2022). Does ABO blood group influence antibody response to SARS-CoV-2 vaccination? Vox Sang. 117 (5), 754–755. doi:10.1111/vox.13241
Wang, D., Hu, B., Hu, C., Zhu, F., Liu, X., Zhang, J., et al. (2020). Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in wuhan, China. JAMA 323 (11), 1061–1069. doi:10.1001/jama.2020.1585
Watanabe, Y., Mendonca, L., Allen, E. R., Howe, A., Lee, M., Allen, J. D., et al. (2021). Native-like SARS-CoV-2 spike glycoprotein expressed by ChAdOx1 nCoV-19/azd1222 vaccine. ACS Cent. Sci. 7 (4), 594–602. doi:10.1021/acscentsci.1c00080
Yin, Y. W., Sheng, Y. J., Wang, M., Ma, Y. Q., and Ding, H. M. (2021). Interaction of serum proteins with SARS-CoV-2 RBD. Nanoscale 13 (30), 12865–12873. doi:10.1039/d1nr02687a
Zamboni, P., and Gemmati, D. (2007). Clinical implications of gene polymorphisms in venous leg ulcer: A model in tissue injury and reparative process. Thromb. Haemost. 98 (1), 131–137. doi:10.1160/th06-11-0625
Zauli, G., Tisato, V., and Secchiero, P. (2020). Rationale for considering oral idasanutlin as a therapeutic option for COVID-19 patients. Front. Pharmacol. 11, 1156. doi:10.3389/fphar.2020.01156
Zhang, J., Zeng, H., Gu, J., Li, H., Zheng, L., and Zou, Q. (2020). Progress and prospects on vaccine development against SARS-CoV-2. Vaccines (Basel) 8 (2), E153. doi:10.3390/vaccines8020153
Keywords: SARS-CoV-2, COVID-19, SNPs, vaccine, BNT162b2 (Pfizer–BioNTech), ChAdOx1 (AstraZeneca), vaccine pharmacogenomics, COVIDomics
Citation: Gemmati D, Longo G, Gallo I, Silva JA, Secchiero P, Zauli G, Hanau S, Passaro A, Pellegatti P, Pizzicotti S, Serino ML, Singh AV and Tisato V (2022) 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. Front. Genet. 13:1028081. doi: 10.3389/fgene.2022.1028081
Received: 25 August 2022; Accepted: 18 October 2022;
Published: 30 November 2022.
Edited by:
Henrik Rasmussen, Institute of Biological Psychiatry, DenmarkReviewed by:
Sandra Paulina Smieszek, Vanda Pharmaceuticals Inc., United StatesRosanna Asselta, Humanitas University, Italy
Copyright © 2022 Gemmati, Longo, Gallo, Silva, Secchiero, Zauli, Hanau, Passaro, Pellegatti, Pizzicotti, Serino, Singh and Tisato. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Donato Gemmati, d.gemmati@unife.it