Skip to main content

Prognostic markers for the clinical course in the blood of patients with SARS-CoV-2 infection

Abstract

Background

The presentation of peptides and the subsequent immune response depend on the MHC characteristics and influence the specificity of the immune response. Several studies have found an association between HLA variants and differential COVID-19 outcomes and have shown that HLA genotypes are associated with differential immune responses against SARS-CoV-2, particularly in severely ill patients. Information, whether HLA haplotypes are associated with the severity or length of the disease in moderately diseased individuals is absent.

Methods

Next-generation sequencing-based HLA typing was performed in 303 female and 231 male non-hospitalized North Rhine Westphalian patients infected with SARS-CoV2 during the first and second wave. For HLA-Class I, we obtained results from 528 patients, and for HLA-Class II from 531. In those patients, who became ill between March 2020 and January 2021, the 22 most common HLA-Class I (HLA-A, -B, -C) or HLA-Class II (HLA –DRB1/3/4, -DQA1, -DQB1) haplotypes were determined. The identified HLA haplotypes as well as the presence of a CCR5Δ32 mutation and number of O and A blood group alleles were associated to disease severity and duration of the disease.

Results

The influence of the HLA haplotypes on disease severity and duration was more pronounced than the influence of age, sex, or ABO blood group. These associations were sex dependent. The presence of mutated CCR5 resulted in a longer recovery period in males.

Conclusion

The existence of certain HLA haplotypes is associated with more severe disease.

Introduction

Since the outbreak of the SARS-CoV-2 pandemic in late 2019, knowledge of the disease and its pathophysiology has rapidly evolved, leading to more effective therapeutic and prophylactic approaches, including vaccination [1, 2]. As a result of the global vaccination campaign, mortality and hospitalization due to COVID-19 have been significantly reduced, especially in developed countries 1,2,3,4,5,6]. However, the relative contribution of T-cell sensitization and neutralizing antibodies is still insufficiently understood [7]. SARS-CoV-2 continuously mutates and poses a significant challenge to cellular and humoral immunity and vaccination [1, 6]. Therefore, a better understanding of the immunological background of the disease is of utmost importance to access more effective prevention and treatment strategies.

SARS-CoV-2 infection is clinically heterogeneous, ranging from subclinical courses to states of critical disease stages with high mortality [2]. Systemic and mucosal immunity requires recognition of peptides bound to the MHC of antigen-presenting cells and subsequent T-cell activation and effector function. Identification of MHC restriction requirements could lead to better-personalized therapeutic strategies, especially when facing numerous vaccination failures [7].

It is reasonable to assume that the course of the disease might depend on the adaptive patients' immune response. In addition, innate immunity by, e.g., less polymorphic TCR gamma/delta T cells, which are mainly active in mucosal tissues, might play a role [8].

Through its essential role in distinguishing between endogenous and exogenous cells [9, 10], the genetic polymorphisms of Human Leukocyte Antigen (HLA) alleles appear to influence the clinical course of patients infected with various RNA viruses, resulting in either protection or enhanced viral persistence [11]. Several HLA polymorphisms have been identified as individual genetic characteristics that affect the individual immune response towards SARS-CoV-2 as well as virus persistence and a predisposition to develop autoimmune diseases [12]. Conflicting results have been published regarding the role of the ABO blood group system, prompting further investigation [13, 14]. Amoroso et al. found that HLA antigens could influence SARS-CoV-2 infection and clinical evolution of COVID-19. They endorsed that blood group A individuals are at greater risk of infection and provided indications on the spread of the disease and clues about infection prognosis and vaccination strategies [9]. Pisanti et al. found that the two most frequent HLA haplotypes in the Italian population, HLA-A*01:01G-B*08:01G-C*07:01G-DRB1*03:01G and HLA-A*02:01G-B*18:01G-C*07:01G-DRB1*11:04G, had a regional distribution overlapping that of COVID-19 and showed, respectively, a positive (suggestive of susceptibility) and negative (suggestive of protection) significant correlation with both COVID-19 incidence and mortality [15]. Other factors that appear to model the individual response include the chemokine receptor CCR5, which plays a role in the humoral response of mucosal immunity to various infections 16,17,18].

Having access to a clinically and predominantly Caucasian patient population, we examined the influence of the most common HLA haplotypes and genotypes, ABO blood group antigens, and the CCR5Δ32 mutant of CD195 concerning to their role in the clinical course. In addition, we compared the results for males and females since differences in severity and mortality are known for SARS-CoV 2 infection [19] and differences for disease susceptibility are discussed [19]. In the context of the pandemic, it is essential to understand the sensitivity of T cells among different population groups due to the differences in their HLA haplotypes. Characterization of immune response among population groups can help to design personalized therapies for patients with high disease severity. Our study analyzed a healthy population (of potential plasma donors) infected with SARS-CoV-2.

Patients, materials and methods

Patients

A total of n = 541 donors were recruited who were willing to donate convalescent plasma after confirmed SARS-CoV2 infection (between January 2020 and January 2021) were subsequently tested for their HLA-Class I (n = 528) and -Class II genotypes (n = 531), their ABO blood groups, and the presence or absence of the CCR5Δ32 mutant (rs6795991, n = 534) for prognostic clinical significance. Disease duration severity was assessed by a questionnaire. Severity was graded according to the WHO COVID-scale [20] into mild disease (WHO° I and IIa) and moderate disease (WHO° IIb and III). All plasma donors were followed up clinically for 600 days. Written informed consent was obtained from all participants for the study. The ethics committee approved the study of Heinrich-Heine-University Dusseldorf, Germany.

All experimental studies were performed according to manufacturers’ standard operating procedures and international guidelines.

Antibody determination

Antibody determination was performed using two serological assays as described before [10, 21]. Spike S1 protein domain specific IgA and IgG antibodies were detected by ELISA (Euroimmun, Germany), Ig antibodies against the Nucleocapsid were verified by Roche Diagnostic Test (Elecsys®).

In addition, neutralizing antibodies were tested. All assays were performed in two samples. Donor plasma was heat inactivated and diluted. Serum neutralization titer was determined by microscopic inspection as the highest serum dilution without virus-induced cytopathic effect (CPE). As positive controls, two previously tested sera from SARS-CoV-2 infected individuals were tested simultaneously in each round. A high-titer control (NT 1:640) and a medium-titer control (NT 1:160) were defined to validate each assay. Inter-assay and intra-assay variations were determined, and these control sera exhibited maximum variation within only one dilution level. In addition, serum from individuals not infected with SARS-CoV-2 and cells without serum were used as negative controls to confirm virus-induced CPE. A neutralization titer of 1:160 was defined as the threshold titer for the binary representation of virus neutralization capacity, allowing the potential plasmapheresis product to serve as a therapeutic agent.

Genotyping and haplotype phasing

To genotype for HLA, CCR5 wild type (wt) and CCR5Δ32 mutant as well as for ABO blood group alleles, an amplicon-based NGS-based approach was used. For HLA six multiplex polymerase chain reactions (PCRs) were used to amplify exons 2, 3 and 4 of the HLA-Class I and HLA-DPB1 genes, exons 2 and 3 of HLA-DRB1, -DQA1 and -DQB1, and exon 3 of the CCR5 gene. The resulting fragments were supplemented with sample-specific barcodes and Illumina-compatible adapters. Sequencing was performed on an Illumina MiSeq device (Illumina Inc., San Diego, USA). A customized software (NGSSequence Analyser, Institute for Transplantation Diagnostics and Cell Therapy, University Hospital Dusseldorf, Dusseldorf, Germany) served for the analysis of the sequence reads using quality control values and high coverage to automate the data analysis. For haplotype phasing of the obtained HLA data, we used the Arlequin software (version 3.5.2.2, available cmpg.unibe.ch/software/arlequin35/), for binding affinities to viral peptides https://t-cov.hse.ru[22].

Statistics

Statistical inference between two groups of interest was assessed by t-test statistics. Kaplan–Meier curves and two-sided log rank tests were used to evaluate the duration-of-infection times between groups of interest. If necessary (violation of the assumption of probability), Bonferroni correction was applied to counteract the multiple comparisons problem. Multivariable Cox regression analysis was performed to identify significant predictors for duration-of-infection, logistic regression for predictors for severity of the disease. Effects of the multivariate Cox regression models were measured using hazard ratio estimates with corresponding 95% confidence intervals, equivalent for the odds ratio in logistic regression. The distribution of the investigated alleles as well as the found haplotypes were tested for their hetero- and homozygosity by checking, whether they fulfill the conditions for a proving if the Hardy–Weinberg equilibrium (available https://ihg.helmholtz-muenchen.de/ihg/snps.html). Only haplotypes that were present in the population at least 5 times were considered. All p-values were determined using two-sided tests, and p-values of less than 0.05 were considered statistically significant. Analyses were performed using IBM SPSS Statistics for Windows (Version 22.0. Released 2013, Armonk, NY: IBM Corp) and R/ R-based blueskystatistics (version 3.6.3, https://www.r-project.org/, version 7.5, https://www.blueskystatistics.com).

Results

A total of 541 plasma donors (median age 43.9 ± 12.5 years), 231 males (median age 43.5 ± 13.0 years) and 310 females (median age 43.5 ± 12.2 years) with mild (n = 320) and moderate disease (n = 221) were included in this analysis. All plasma donors were healthy with an ECOG status of 0 prior infection with SARS-CoV-2, e.g., no co-morbidities like infectious diseases, malignant tumors or immunocompromised. The duration of the disease was 15 ± 7 days. Our patients with SARS-CoV-2 were never hospitalized, and symptoms corresponded to WHO° 1–3 [20]. During the first wave of infection with SARS-CoV-2, we investigated 179 patients (91 males, 88 females), and during the second wave of infection, we studied 362 convalescent plasma donors (140 males, 222 females). The characteristics of all investigated donors are shown in Table 1. In total, a correlation was found between the severity of the clinical picture and the duration of the disease course. In total, patients of wave 1 showed less mild course than patients of wave 2 (p < 0.003). Patients of the 2nd trends to have a shorter course of disease (n.s. p < 0.1). The more pronounced the patient's symptoms, the longer their clinical course and the higher the reached antibody levels (data not shown). We found elevated or stable period Ig N antibody levels of up to 240 days after the end of symptoms. In all patients in whom elevated antiviral IgG antibodies had been observed the respective anti-S IgG antibodies remained detectable for least 240 days as we reported previously for the patients of the first wave [10, 23].

Table 1 Characteristics of all investigated patients

Association of HLA haplotypes and severity of COVID-19

ABO blood group (BG) genotype could be assessed in 532 patients, for O alleles in 532 cases, for A alleles in 532 cases, for CCR Δ32 in 534 cases.

Haplotype phasing of the obtained HLA data revealed several HLA-Class I and Class II haplotypes (Table 2). For HLA-Class I, we obtained results from 528 patients, and for HLA-Class II from 531.

Table 2 The frequencies of the examined HLA haplotypes

Frequencies of homozygous and heterozygous variant carriers were tested for each allele/haplotype, whether they fulfill the assumption of the Hardy–Weinberg equilibrium (HWE) in the whole group as well as in the mild and moderate disease group at a significance level of less than p < 0.02 (see Additional file 1). Deviations from the HWE (DHWE) were found for number of OBG alleles with more homozygous O-carriers, based on more carriers in the mild disease group. DHWE was also seen for the HLA haplotype A*02:01_B*40:01_C*03:04 (HLA-ABC_05), due to more homozygous than expected, In haplotype HLA-DRB1*03:01_DQA1*05:01_DQB1*02:01 (DRB-DQA_DQB_4) 6 out of 7 homozygous were seen for in the moderate diseased group (DHWE), whereas all of the 6 homozygous carriers of the HLA haplotype DRB1*11:01_DQA1*05:01_DQB1*03:01 (DRB_DQA_DQB_05) were found in the mild disease group (DHWE).

Univariate analysis of the HLA haplotypes showed that the presence the HLA-Class II haplotype DRB1*07:01_DQA1*02:01_DQB1*03:01 (36/531 6.8%, OR 2.16, 95% confidence interval 1.09–4.36, “DRB_DQA_DQB_06”) was predictive for a mild disease course of SARS-CoV-2 infection. Two haplotypes (class I A*02:01_B*07:02_C*07:02 (35/528, 6.6%, OR 0.28 95% CI 0.10–0.65, “ABC_04”); class II DRB1*12:01_DQA1*05:01_DQB1*03:01 (19/531, 2,8%, OR 0.164 95%CI 0.026–0.58, “DRB_DQA_DQB_15”)) were associated with a more severe disease course.

When looking for males and females separately, we found only in females a significant association for the presence of haplotype “ABC_04” and “ DRB_DQA_DQB_06” with severity of the disease (Fig. 1a, b).

Fig. 1
figure 1

Univariate association of HLA haplotypes with mild/moderate disease (relative risk, OR, log rank (LR)) in females (A) and males (B) and for disease duration (hazard ratio, (Cox)) in females (C) and males (D). Haplotype frequency within the group is plotted on the x-axis and the y-axis shows the significance levels of the OR for moderate disease / HR for disease duration. The size of the spheres indicates the effect size and was calculated as absolute value of the logarithms of OR/HR. Full spheres represent a favorable association with OR or HR, and empty spheres represent an association with a higher RR or HR. The 0.05 significance level is shown as parallel to the x-axis. Alleles/haplotypes which were either significant in females or males for OR (panel A and B) are denoted in both panels. Same applies for HR (panel C and D). *Note that the sphere for the HLA-Class II haplotype DRB_DQA_DQB_15 in panel A is not on scale, since all females homozygous for this haplotype had had a mild course of disease (n = 8)

In the whole group the HLA-Class I haplotype A*03:01_B*35:01_C*04:01 (12/528, 2,3%, HR 2.1 95%CI 1.19–3.75, “ABC_11”) was associated with shorter disease, the HLA-Class II haplotype DRB1*15:01_DQA1*01:02_DQB1*06:02 (126/531, 25,5% OR 0.81 95%CI 0.66–0.99, “DRB_DQA_DQB_05”) with longer duration. Looking at females only, HLA-Class I haplotype “ABC_11” (A*03:01_B*35:01_C*04:01, HR 3.8 95%CI 1.2–11.8), haplotype ABC_17 (A*01:01_B*07:02_C*07:02, HR 3.4, 95%CI 1.1–10.6), and ABC_22 (A*29:02_B*44:03_C*16:01 HR 5.30 95%CI 1.29–2.73) with shorter duration as well as HLA-Class II haplotype DRB_DQA_DQB_05 (HR 1.6 CI 1.2—2.2). HLA-Class II haplotype DRB_DQA_DQB_17 was associated with longer duration (DRB1*16:01_DQA1*01:02_DQB1*05:02, HR 0.34 95%CI 0.14–0.85) (Fig. 1c).

In males also HLA-Class II haplotypes reached significance DRB_DQA_DQB_10 (DRB1*04:01_DQA1*03:01_DQB1*03:01; HR 1.88 95%CI 1.04–3.40), DRB_DQA_DQB_12 (DRB1*13:02_DQA1*01:02_DQB1*06:04; HR 1.85 95%CI 1.03–3.34), HLA_DRB_DQA_DQB_18 (DRB1*09:01_DQA1*03:01_DQB1*03:03; HR 2.86 95%CI 1.17–6.99) for shortening the disease (Fig. 1d).

Within multivariate analysis for severity of disease, beside age, gender, CCR5Δ32 alleles, number of blood group O alleles all HLA haplotypes were included, which showed either a Hardy Weinberg deviation and/or a univariate significance for association with severity of the disease in males and/ or females. Analysis was further weighted for wave of infection.

Results for the whole group are shown in Fig. 2a, for females in Fig. 2b and for males in Fig. 2c. Most haplotypes showed the same association in favoring mild or moderate disease: HLA haplotype ABC_04 (A*02:01_B*07:02_C*07:02) favors mild disease in the whole group (p < 0.001) as well in females (p < 0.001) and males (p < 0.04), ABC_07 (A*24:02_B*07:02_C*07:02) in the whole group (p < 0.015), males (p < 0.004), in females n.s. but also favoring in direction more moderate disease.

Fig. 2
figure 2

Multivariate association of age, gender, and HLA haplotypes as risk factors for disease severity (WHO° I & IIa (mild) vs WHO° IIb & III). A shows results for the whole group, B for females, C for males.*Note that for the HLA-Class II haplotype DRB_DQA_DQB_15 in panel B whisker is not on scale, since all females homozygous for this haplotype had had a mild course of disease (n = 8)

DRB_DQA_DQB_06 (DQA1*02:01_DQB1*03:01) in the whole group p < 0.002, females (p < 0.001), in males n.s. but also favoring more moderate disease.

DRB_DQA_DQB_15 (DRB1*12:01_DQA1*05:01_DQB1*03:01) in the whole group (p < 0.003), in males (p < 0.025) favoring mild disease. All 8 females carrying this haplotype had also mild disease.

Haplotype ABC_05 (A*02:01_B*40:01_C*03:04) favors only in males for mild disease (p < 0.05).

DRB_DQA_DQB_05 (DRB1*11:01_DQA1*05:01_DQB1*03:01) favors for more moderate disease only in males.

Multivariate analysis for duration of disease included age, gender, CCR5Δ32 alleles, number of blood group O alleles and all those HLA haplotypes, which showed a univariate significance for association duration of the disease in males and/ or females (Fig. 3a).

Fig. 3
figure 3

Multivariate association of age, gender, and HLA haplotypes as risk factors for disease duration (longer disease (HR < 1) vs. shorter disease duration (HR > 1). A shows results for the whole group, B for females, C for males

Age was associated with a significant, but small impact to prolong the disease (p < 0.001, HR 0.987), males had a shorter illness course (HR 1.17, p < 0.015). CCR5Δ32 mutation only impacted males (HR 0.65, p < 0.001) but not females.

Following HLA-Class I haplotypes were associated with shorter disease duration: HLA_ABC_11 (A*11:01_B*35:01_C*04:01; whole group HR 1.61, p < 0.035, in females HR 5.33, p < 0.001, in males n.s. but same direction of influence), HLA_ABC_17 (A*01:01_B*07:02_C*07:02, whole group HR 1.83, p < 0.007, in females HR 3.53, p < 0.001, in males n.s. but same direction of influence). In females in addition haplotype HLA_ABC_22 (A*29:02_B*44:03_C*16:01; HR 1.45, p < 0.015) showed same association.

In men, HLA_DRQ_DQA_DQB_10 (DRB1*04:01_DQA1*03:01_DQB1*03:02, HR 2.26, p < 0.001) was also associated with shorter disease course.

Haplotype HLA_DRQ_DQA_DQB_17 (DRB1*16:01_DQA1*01:02_DQB1*05:02) was associated with longer disease course in females (HR 0.37, p < 0.003), as was HLA_DRQ_DQA_DQB_02 (DRB1*15:01_DQA1*01:02_DQB1*06:02, HR 0.79, p < 0.02).

Two haplotypes were associated with an opposite effect in males and females: HLA_DRQ_DQA_DQB_05 (DRB1*11:01_DQA1*05:01_DQB1*03:01, HRfemales 1.42 (supporting shorter disease) pfemales < 0.002; HRmales 0.74 pmales < 0.043) and, vice versa HLA_DRQ_DQA-DQB_18 (DRB1*09:01_DQA1*03:02_DQB1*03:03, HRfemales 0.61 pfemales < 0.05; HRmales 2.26 (supporting shorter disesease) pmales < 0.001).

When coding the beneficial and the disadvantageous HLA haplotypes for disease duration specific for sex separately (Fig. 4), females carrying beneficial female haplotypes had a median disease duration of 10 days (men 12 days), women carrying unfavorable haplotypes had a disease duration of 14 days (men 13 days). Women carrying either favorable and unfavorable or neither of those haplotypes, had a median disease duration of 12 days.

Fig. 4
figure 4

Kaplan–Meier estimates for the impact of disease specific favorable or unfavorable HLA haplotypes on duration of the disease when coded for significant alleles in females (A&C) or males (B&D)

When looking at the male specific haplotypes, males carrying beneficial male haplotypes had median disease duration of 10 days (females 12 days), men carrying unfavorable haplotypes 14 days (females 13 days), men with unfavorable and favorable or neither of those haplotypes had a disease duration of 12 days (females 13 days). For significance levels, see Fig. 4.

We further analyzed the clinical outcome of plasma donors specifically for the heterozygous CCR5Δ32 mutation. We investigated whether this mutation had an impact on the duration of symptoms. In patients infected during the 1st wave (up until week 20/2020), heterozygous mutation of CCR5Δ32 was associated with a significantly longer disease course (p = 0.013) [10]. This significant difference was confirmed for the whole group only in males (Fig. 5, p < 0.032).

Fig. 5
figure 5

Kaplan–Meier estimates for the impact of CCR5 wt and CCR5Δ32 on disease duration in females (A) and males (B)

The influence of the HLA haplotypes on disease severity and duration was stronger than the influence of age, sex, ABO blood group, or wave of infection with different mutants of the virus.

Numbers of in silico predicted viral T-cell epitopes did not correlate with the severity of SARS-CoV-2 infections.

Discussion

Since the onset of the SARS-CoV-2 pandemic, the heterogeneity of the clinical course has been an essential issue for investigations. While in some patients, the infection goes unnoticed, in others, the same infective agent leads to a severe course of disease accompanied by high morbidity and mortality [21, 24]. Besides the immune reaction of the patient itself, there could be a combination of several factors that lead to these different outcomes like gender, distance regulations and length of the infection wave or season.

Nevertheless, individual immunologic pathways based on metabolism, age, sex and genetics are known to play a crucial role [25]. Our study demonstrates a prolonged or more severe clinical course in carriers of some HLA haplotypes while other haplotypes favor shorter or milder outcomes.

Carriers of the HLA haplotypes HLA-A*02:01_B*07:02_C*07:02 and HLA-DRB1*12:01_DQA1*05:01_DQB1*03:01 had only mild disease symptoms, whereas carriers of the haplotypes HLA-A*24:02_B*07:02_C*07:02 and HLA-DRB1*07:01_DQA1*02:01_DQB1*03:03 and (especially in men) HLA-DRB1*11:01_DQA1*05:01_DQB1*03:01 showed more moderate/severe symptoms.

As with others [26], we found no strong correlation with previously described class I superhaplotypes [27], although patients with haplotypes containing the A24 superhaplotype motif tend to have a more unfortunately course, whereas patients bearing the B07/B44 motif without A24 trend to have a favorable course. Regarding class II superhaplotypes [28], patients with DR1/DQ1 and DR3/DQ3 superhaplotypes is among those doing worse, whereas beneficial haplotypes contain the superhaplotypes DR4, DR5 and DR9.

SARS-CoV-2 infection and its severity might be related to gender, and males are more likely to develop COVID-19 [29]. Although males are more susceptible to most viral infections, females possess immunological features that render them more vulnerable to distinct immune-related disease outcomes [30]. Interestingly, in our cohort females reported a longer disease duration, whereas men tend to report more severe disease. We questioned whether investigated alleles and or haplotypes may have a differential impact on severity and duration within this non-hospitalized patient cohort. When looking at the severity of the disease, the pattern of haplotypes being favorable or unfavorable differ not relevantly between the genders. When looking at disease duration, which is also a sign of disease control, interestingly three class I haplotypes are associated with shorter disease duration in females, whereas in men class II haplotypes are positively associated with shorter disease—in line with observations, that a defined antigen stimulus would lead to more clonal expansions in CD8 T cells in men than in women [31] and thus differences in class I repertoire may impact disease control in women more than in men. On the other side, women show per se higher T-helper cell levels and activability [32], which leads to the conclusion that differences within the HLA-Class II repertoire on disease control may only be seen in men.

As known, the C-C chemokine receptor 5 (CCR5) is responsible for immune and inflammatory responses by mediation of chemotactic activity in T cells, although it is expressed on different cell types. It has been shown to act as co-receptor for the human immunodeficiency viruses (HIV-1, HIV-2), resulting in resistance to HIV infection in homozygous CCR5Δ32 mutation carriers. CCR5 serves not only as a co-receptor, but also contributes to the functional fitness of CD4 + T cells themselves. CCR5 signaling reduces ceramide levels and thereby increases HLA experienced T‐cell antigen receptor (TCR) nanoclustering in antigen‐experienced CD4 + T cells [7]. This contributes to the more severe outcome in, e.g., early summer meningoencephalitis or West Nile fever in CCR5Δ32 carriers [33]. Previously, in patients of the 1st wave, we noticed that CCR5Δ32 carriers had a longer disease course [10]. If we look now at the not so severely ill patients of the 2nd wave, we observe this only in males. This corresponds to two expectations: first, in the not so severely ill population this effect can only be seen in the group where an influence of the MHC class II repertoire is seen for duration of the disease, as we see in males and, second, due to the association of female gender with increased CD4 + T cell CCR1–CCR5 gene and protein expression [34], the increased CCR expression, which correlates with enhanced in vitro chemotaxis response, the influence of the mutation will be reversed. As others in larger cohorts, we saw no impact of the CCR5Δ32 mutation on the primary severity of the disease [17]. In our cohort, the prevalence of the mutation was 17.4%, but 11.5% in the 26 asymptomatic patients. In a larger cohort of 164 asymptomatic SARS-CoV-2 infected patients Hubacek et al. noticed that during the 1st wave the incidence of CCR5Δ32 was higher in (23.8%) than in 252 symptomatic patients (16.7%) [18]. This results of a protective role of the mutation could not confirmed by the analysis in 110,544 stem cell donors [17]. When looking at critical ill patients, which were not included in our study, CCR5 might mediate immune dysregulation and thus leronlimab, a CCR5 blocking antibody may limit hyperinflammation in critical ill patients [16]. As of September 2022, randomized clinical trials evaluating leronlimab in patients with COVID-19 could not confirm this (NCT04343651) or are not recruiting (NCT04347239).

Previously, we found that ABO heterogeneity correlated with a higher risk of a more severe disease in patients of the 1st wave [10]. Although this effect is now no longer visible to the same extent, it is striking that for the blood group O allele the Hardy–Weinberg equilibrium is violated: a disproportionately large number of blood group O allele homozygotes are found in the group of those with only mild disease.

However, our study has several limitations. First, since it only covers patients of the 1st and 2nd wave before variants of concern emerged in Germany, the found associations are related to the original variant. Although we compared the patients of the first and second wave together, they differ: in the first group, most of the patients were skiers who came from a SARS-CoV-2 hotspots in Austria. In the 2nd group women predominated. They were recruited following public calls for convalescent plasma donation, where usually more women than men respond. It could be assumed, that in the first wave, which occurred in the cold season, patients were exposed to a higher viral load due to a lack of protective strategies such as clearance regimens or protective masks. In addition, the second wave occurred in the summer. Also, all patients were non-hospitalized – the highest WHO° was III. Finally, all patients were potentially able to donate blood and, therefore, primarily healthy individuals without severe co-morbidities before the infection.

Although limited by the size and retrospective nature of this cohort study, we were able to identify HLA haplotypes associated with the severity and duration of disease in 541 SARS-CoV-2 infected patients. While we found no sex association for disease severity and association with the haplotypes studied, we did see one for disease duration. Two of the 22 HLA-Class I haplotypes studied were significantly associated (HR > 2) with a shorter disease course in females. In contrast, two of the 22 HLA-Class II examined were favorably associated with disease duration only in men (HR > 2). Men with a CCR5Δ32 mutation had a longer disease course.

Availability of data and materials

Data are applicable from JCF and CM.

Abbreviations

CCR5:

C–C chemokine receptor type 5

CD4:

Cluster of differentiation 4

COVID-19:

Coronavirus disease 2019

CP:

Convalescent plasma

CPE:

Cytopathic effect

e.g.:

Exempli gratia

ELISA:

Enzyme-linked immunosorbent assay

HIV:

Human immunodeficiency virus

HR:

Hazard ratio

HLA:

Human leukocyte antigen system

Ig:

Immunoglobulin

IL:

Interleukin

IQR:

Inter-quartile range

MCP:

Monocyte chemo-attractant protein

MHC:

Major histocompatibility complex

PCR:

Polymerase chain reaction

RNA:

Ribonucleic acid

SARS-CoV-2:

Severe acute respiratory syndrome coronavirus 2

WHO:

World Health Organization

Wt:

Wild type

References

  1. Voysey M, Clemens SAC, Madhi SA, et al. Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: an interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK. Lancet. 2021;397(10269):99–111. https://doi.org/10.1016/s0140-6736(20)32661-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Yarmarkovich M, Warrington JM, Farrel A, Maris JM. Identification of SARS-CoV-2 vaccine epitopes predicted to induce long-term population-scale immunity. Cell Rep Med. 2020;1(3):100036. https://doi.org/10.1016/j.xcrm.2020.100036.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Baden LR, El Sahly HM, Essink B, et al. Efficacy and safety of the mRNA-1273 SARS-CoV-2 vaccine. N Engl J Med. 2020. https://doi.org/10.1056/NEJMoa2035389.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Bhatnager R, Bhasin M, Arora J, Dang AS. Epitope based peptide vaccine against SARS-COV2: an immune-informatics approach. J Biomol Struct Dyn. 2020. https://doi.org/10.1080/07391102.2020.1787227.

    Article  PubMed  Google Scholar 

  5. BionTech Pa. Vaccines and Related Biological Products Advisory Committee Meeting. FDA Briefing Document Pfizer-BioNTech COVID-19 Vaccine 2020.

  6. Polack FP, Thomas SJ, Kitchin N, et al. Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine. N Engl J Med. 2020;383(27):2603–15. https://doi.org/10.1056/NEJMoa2034577.

    Article  CAS  PubMed  Google Scholar 

  7. Martín-Leal A, Blanco R, Casas J, et al. CCR5 deficiency impairs CD4(+) T-cell memory responses and antigenic sensitivity through increased ceramide synthesis. Embo j. 2020;39(15):e104749.

    Article  PubMed  PubMed Central  Google Scholar 

  8. von Massow G, Oh S, Lam A, Gustafsson K. Gamma delta T Cells and their involvement in COVID-19 virus infections. Front Immunol. 2021;12:741218. https://doi.org/10.3389/fimmu.2021.741218.

    Article  CAS  Google Scholar 

  9. Amoroso A, Magistroni P, Vespasiano F, et al. HLA and AB0 Polymorphisms may influence SARS-CoV-2 infection and COVID-19 severity. Transplantation. 2021;105(1):193–200. https://doi.org/10.1097/tp.0000000000003507.

    Article  CAS  PubMed  Google Scholar 

  10. Fischer JC, Schmidt AG, Bolke E, et al. Association of HLA genotypes, AB0 blood type and chemokine receptor 5 mutant CD195 with the clinical course of COVID-19. Eur J Med Res. 2021;26(1):107. https://doi.org/10.1186/s40001-021-00560-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Crux NB, Elahi S. Human Leukocyte Antigen (HLA) and immune regulation: how do classical and non-classical HLA alleles modulate immune response to human immunodeficiency virus and hepatitis C virus infections? Front Immunol. 2017;8:832. https://doi.org/10.3389/fimmu.2017.00832.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Gomez-Mayordomo V, Montero-Escribano P, Matias-Guiu JA, Gonzalez-Garcia N, Porta-Etessam J, Matias-Guiu J. Clinical exacerbation of SARS-CoV2 infection after fingolimod withdrawal. J Med Virol. 2020. https://doi.org/10.1002/jmv.26279.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hoiland RL, Fergusson NA, Mitra AR, et al. The association of ABO blood group with indices of disease severity and multiorgan dysfunction in COVID-19. Blood Adv. 2020;4(20):4981–9. https://doi.org/10.1182/bloodadvances.2020002623.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Golinelli D. On the association between the ABO blood group and COVID-19 susceptibility. Blood Transfus. 2020. https://doi.org/10.2450/2020.0356-20.

    Article  PubMed  Google Scholar 

  15. Pisanti S, Deelen J, Gallina AM, et al. Correlation of the two most frequent HLA haplotypes in the Italian population to the differential regional incidence of Covid-19. J Transl Med. 2020;18(1):352. https://doi.org/10.1186/s12967-020-02515-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Patterson BK, Seethamraju H, Dhody K, et al. Disruption of the CCL5/RANTES-CCR5 pathway restores immune homeostasis and reduces plasma viral load in critical COVID-19. medRxiv : the preprint server for health sciences. 2020. https://doi.org/10.1101/2020.05.02.20084673.

    Article  Google Scholar 

  17. Bernas SN, Baldauf H, Wendler S, et al. CCR5Δ32 mutations do not determine COVID-19 disease course. Int J Infect Dis : IJID : official publication of the International Society for Infectious Diseases. 2021;105:653–5. https://doi.org/10.1016/j.ijid.2021.02.108.

    Article  CAS  Google Scholar 

  18. Hubacek JA, Dusek L, Majek O, et al. CCR5Delta32 deletion as a protective factor in Czech first-wave COVID-19 subjects. Physiol Res. 2021;70(1):111–5. https://doi.org/10.33549/physiolres.934647.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Jin J-M, Bai P, He W, et al. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. 2020;8:152. https://doi.org/10.3389/fpubh.2020.00152.

    Article  PubMed  PubMed Central  Google Scholar 

  20. WHO Working Group on the Clinical Characterisation and Management of COVID-19 infection. A minimal common outcome measure set for COVID-19 clinical research. Lancet Infect Dis 2020;20(8):e192–7. https://doi.org/10.1016/s1473-3099(20)30483-7.

    Article  CAS  Google Scholar 

  21. Müller L, Andrée M, Ostermann PN, et al. SARS-CoV-2 Infection in fully vaccinated individuals of old age strongly boosts the humoral immune response. Front Med. 2021;8:746644. https://doi.org/10.3389/fmed.2021.746644.

    Article  Google Scholar 

  22. Nersisyan S, Zhiyanov A, Shkurnikov M, Tonevitsky A. T-CoV: a comprehensive portal of HLA-peptide interactions affected by SARS-CoV-2 mutations. Nucleic Acids Res. 2022;50(D1):883–7. https://doi.org/10.1093/nar/gkab701.

    Article  CAS  Google Scholar 

  23. Bolke E, Matuschek C, Fischer JC. Loss of Anti-SARS-CoV-2 antibodies in mild covid-19. N Engl J Med. 2020;383(17):1694–5. https://doi.org/10.1056/NEJMc2027051.

    Article  PubMed  Google Scholar 

  24. Nelde A, Bilich T, Heitmann JS, et al. SARS-CoV-2-derived peptides define heterologous and COVID-19-induced T cell recognition. Nat Immunol. 2021;22(1):74–85. https://doi.org/10.1038/s41590-020-00808-x.

    Article  PubMed  Google Scholar 

  25. Lingeswaran M, Goyal T, Ghosh R, et al. Inflammation, immunity and immunogenetics in COVID-19: a narrative review. Ind J Clin Biochem : IJCB. 2020;35(3):260–73. https://doi.org/10.1007/s12291-020-00897-3.

    Article  CAS  Google Scholar 

  26. Schetelig J, Heidenreich F, Baldauf H, et al. Individual HLA-A, -B, -C, and -DRB1 genotypes are no major factors which determine COVID-19 severity. Front Immunol. 2021;12:698193. https://doi.org/10.3389/fimmu.2021.698193.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Sidney J, Peters B, Frahm N, Brander C, Sette A. HLA class I supertypes: a revised and updated classification. BMC Immunol. 2008;9(1):1. https://doi.org/10.1186/1471-2172-9-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Doytchinova IA, Flower DR. In Silico Identification of Supertypes for Class II MHCs. J Immunol. 2005;174(11):7085–95. https://doi.org/10.4049/jimmunol.174.11.7085.

    Article  CAS  PubMed  Google Scholar 

  29. Nguyen NT, Chinn J, De Ferrante M, Kirby KA, Hohmann SF, Amin A. Male gender is a predictor of higher mortality in hospitalized adults with COVID-19. PLoS ONE. 2021;16(7):e0254066. https://doi.org/10.1371/journal.pone.0254066.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Jacobsen H, Klein SL. Sex differences in immunity to viral infections. Front Immunol. 2021;12:720952. https://doi.org/10.3389/fimmu.2021.720952.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Schneider-Hohendorf T, Görlich D, Savola P, et al. Sex bias in MHC I-associated shaping of the adaptive immune system. Proc Natl Acad Sci USA. 2018;115(9):2168–73. https://doi.org/10.1073/pnas.1716146115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Gay L, Melenotte C, Lakbar I, et al. Sexual dimorphism and gender in infectious diseases. Front Immunol. 2021;12:698121. https://doi.org/10.3389/fimmu.2021.698121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Kindberg E, Mickiene A, Ax C, et al. A deletion in the chemokine receptor 5 (CCR5) gene is associated with tickborne encephalitis. J Infect Dis. 2008;197(2):266–9. https://doi.org/10.1086/524709.

    Article  CAS  PubMed  Google Scholar 

  34. Mo R, Chen J, Grolleau-Julius A, Murphy HS, Richardson BC, Yung RL. Estrogen regulates CCR gene expression and function in T lymphocytes. J Immunol. 2005;174(10):6023–9. https://doi.org/10.4049/jimmunol.174.10.6023.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

Open Access funding enabled and organized by Projekt DEAL. The study was funded by the German government, RES-Q-HR trial.

Author information

Authors and Affiliations

Authors

Contributions

JCF, EB, CM developed the idea of this investigation. DJ, JCF, VB, JE, JT, OA, MS, DH, KSB, TF, BJ, DH, VK, DKM, WTK, MvG, WB, JH, BT, CP, LS, KS, SC, AP, KM, LS, OG, EB, AE, IE, MP, CM, AS, CC, DH, AR wrote parts of the manuscript. JCF, CM, DJ and EB did the literature research and prepared the data for analysis. JCF prepared the figures. EB, JCF, DJ, CM and MS contributed significantly to the discussion and the interpretation of the results. VB, JE, DH, KSB, OA, JT performed the laboratory tests. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Edwin Bölke.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the local ethics board registration and approval number 2020-1116.

Consent for publication

All authors gave consent for the publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Tests for deviation and tests for association from Hardy Weinberg equilibration.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fischer, J.C., Balz, V., Jazmati, D. et al. Prognostic markers for the clinical course in the blood of patients with SARS-CoV-2 infection. Eur J Med Res 27, 255 (2022). https://doi.org/10.1186/s40001-022-00864-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-022-00864-z

Keywords