key: cord-0851029-jluh976j
authors: Valenti, Luca; Griffini, Samantha; Lamorte, Giuseppe; Grovetti, Elena; Uceda Renteria, Sara Colonia; Malvestiti, Francesco; Scudeller, Luigia; Bandera, Alessandra; Peyvandi, Flora; Prati, Daniele; Meroni, Pierluigi; Cugno, Massimo
title: Chromosome 3 cluster rs11385942 variant links complement activation with severe COVID-19
date: 2021-01-09
journal: J Autoimmun
DOI: 10.1016/j.jaut.2021.102595
sha: d5ab95534cabd4ffb56aa8a475427d8e7d07ef41
doc_id: 851029
cord_uid: jluh976j

BACKGROUND: Genetic variation at a multigene cluster at chromosome 3p21.31 and the ABO blood group have been associated with the risk of developing severe COVID-19, but the mechanism remains unclear. Complement activation has been associated with COVID-19 severity. OBJECTIVE: The aim of this study was to examine whether chromosome 3p21.31 and the ABO variants are linked to the activation of the complement cascade in COVID-19 patients. METHODS: We considered 72 unrelated European hospitalized patients with genetic data and evaluation of circulating C5a and soluble terminal complement complex C5b-9 (SC5b-9). Twenty-six (36.1%) patients carried the rs11385942 G>GA variant and 44 (66.1%) non-O blood group associated with increased risk of severe COVID-19. RESULTS: C5a and SC5-b9 plasma levels were higher in rs11385949 GA carriers than in non-carriers (P = 0.041 and P = 0.012, respectively), while C5a levels were higher in non-O group than in O group patients (P = 0.019). The association between rs11385949 and SC5b-9 remained significant after adjustment for ABO and disease severity (P = 0.004) and further correction for C5a (P = 0.018). There was a direct relationship between upper airways viral load and SC5b-9 in carriers of the rs11385949 risk allele (P = 0.032), which was not observed in non-carriers. CONCLUSIONS: The rs11385949 G>GA variant, tagging the chromosome 3 gene cluster variation and predisposing to severe COVID-19, is associated with enhanced complement activation, both with C5a and terminal complement complex, while non-O blood group with C5a levels. These findings provide a link between genetic susceptibility to more severe COVID-19 and complement activation.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection displays a wide range of clinical manifestations spanning from an asymptomatic condition, observed in almost 50% of individuals with evidence of upper respiratory viral replication [1] , to coronavirus disease 2019 [2, 3] . The disease can be associated with minor upper airway manifestations or interstitial pneumonia that can evolve into severe and lifethreatening acute respiratory distress syndrome (ARDS) [2] [3] [4] [5] . An exaggerated "hyperimmune" response to lower respiratory airways and systemic infection is believed to be involved in ARDS pathogenesis [6] . Furthermore, COVID-19 is associated with a range of non-respiratory conditions affecting multiple organs, including the heart, circulatory system, kidney and the liver, whose pathogenesis involves endothelial activation and the immune response against SARS-CoV-2 [7, 8] .

Older age, male sex, and dysmetabolism predispose to severe COVID-19 [9, 10] . We showed that inherited genetic factors contribute to determine the predisposition to COVID-19 related respiratory failure [11] . The main risk factor for severe respiratory disease was identified at a chromosome 3 gene cluster locus, encompassing variation at several candidate genes potentially implicated in regulating bronchial physiology, viral entry, and immune response. These include SCL6A20 (solute carrier family 6, member 20), LZTFL1 (leucine zipper transcription factor like 1) , FYCO1 (FYVE [ Fab 1, YOTB, Vac 1, EEA1] and coiled-coil domain autophagy adaptor 1) and the chemokine receptors CCR9 (C-C motif chemokine receptor 9) , CXCR6 (C-X-C motif chemokine receptor 6), XCR1 (X-C motif chemokine receptor 1), with a possible impact extended on expression of the neighboring CCR1, CCR2 and CCR3 [11] [12] [13] . However, due to the high linkage disequilibrium of this region inherited from Neanderthals [14] , it was not yet possible to pinpoint the causal variant(s), J o u r n a l P r e -p r o o f and the mechanism underlying the association between chromosome 3 cluster variation and severe COVID-19 remains to be defined. In the FOGS study, a recently characterized cohort of a referral center with genetic characterization, we showed the chromosome 3 variation (the top hit rs11385942 G>GA variant) was associated with the susceptibility to develop COVID-19 with respiratory failure [11, 15] . Although chromosome 3 gene cluster variation has repeatedly been confirmed to represent the main genetic risk factor for COVID-19 severity [11, 12, 16] , in the FOGS cohort rs11385942 did not have a major impact on biomarkers of systemic inflammation, including C reactive protein, ferritin, IL-6 and circulating neutrophils and lymphocytes [17] . Furthermore, ABO blood group was involved in the susceptibility to COVID-19, with carriage of blood O conferring modest protection.

Again a sound mechanistic explanation for such association was not yet reported [11, 12, [18] [19] [20] .

In patients with severe disease, activation of the complement cascade, which may be facilitated by genetic predisposition, has been linked with inflammation grade and activation of coagulation, a typical feature of COVID-19 leading to thrombotic complications [21] [22] [23] .

This process may be directly triggered by the SARS-CoV-2 Spike protein [24] as well as by the antibodies against the virus themselves. We recently showed that in hospitalized patients complement activation was associated with disease severity and perturbation of the endothelium, a main target of SARS-CoV-2 infection [7, 25] . Since it is not known whether the susceptibility to the activation of the complement pathway following SARS-CoV-2 infection contributes the inter-individual variability in COVID-19 severity, the aim of this study was to examine the possible association between complement activation and chromosome 3 gene cluster variation and/or ABO blood group in a well characterized cohort of hospitalized COVID-19 patients. 

We considered 72 patients, who were enrolled in both the FOGS [11, 17] and the COHERENT studies [25] at the Fondazione IRCCS Ca' Granda Milan, between March 16 th and April 29 th 2020 (at the beginning of the 1 st wave in Italy), for whom both genetic data and characterization of complement activation were available. To avoid confounders, we restricted the analysis to unrelated individuals of European descent [11, 17] . These patients were selected by the larger case-control FOGS cohort of 508 cases and 890 controls up to June 2020, whose clinical features were previously reported [17] . The severity of COVID-19 was graded on the basis of clinical, radiological and laboratory data as follows: mild when the patients were kept under observation with low intensity care, moderate when they required intermediate care with non-invasive ventilation, or severe when they were admitted to an intensive care unit for mechanical ventilation [23, 25] .

Nasopharingeal swabs were tested for the presence of SARS-CoV-2 using a One-Step Reverse Transcription Real-Time PCR kit (GeneFinder COVID-19 Plus RealAmp kit, ELITe, platform InGenius ELITech, Torino, Italy) according to manufacturer instructions, as previously described [25] . The InGenius platform displays a sensibility of 1000 cp/ml and a specificity of 100%. A gene was considered detected when the cycle threshold (Ct) value resulted ≤ 45 while the assay was considered negative for a Ct value > 45. 

The rs11385942 G>GA variant and ABO blood group were determined as previously described [11] ; blood typing was confirmed by direct immunophenotyping in >95% of cases (100% concordance), while rs11385942 was confirmed by direct typing by Taqman assays in 16 cases (100% concordance) [17] .

Complement activation status was determined on the day of admission, as previously described [25] . Plasma levels of soluble C5b-9 (SC5b-9) were measured using a solid-phase assay (MicroVue Complement SC5b-9 Plus EIA kit, Quidel Corporation, San Diego, CA, USA); intra-and inter-assay coefficients of variation (CVs) were respectively 6.8% and 13.1%; the lower detection limit was 3.7 ng/mL. Plasma C5a levels were measured using an immunoenzymatic method (MicroVue Complement C5a EIA, Quidel Corporation); intraand inter-assay CVs of <12%; the lower detection limit was 0.01 ng/ml.

Comparisons were performed by fitting outcomes with univariate and multivariate general linear models. Models were adjusted for confounding factors and genetic factors, as 

The clinical features of patients, stratified by carriage of the rs11385942 variant, are shown in Table 1 . In the present cohort, which includes all the disease stages, the rs11385942 GA risk allele tended to be associated with the group of severe COVID-19.

Baseline biochemical features, including inflammatory and coagulation indices and upper airways viral load, were also comparable between carriers and non-carriers. However, risk variant carriers tended to be younger and had a significantly lower burden of cardiometabolic comorbidities (p=0.0045).

Circulating C5a levels were not significantly associated with available circulating biomarkers of inflammation and tissue damage. On the other hand, terminal complement complex SC5b-9 levels were directly associated with markers of tissue damage, including AST (P=0.038), GGT (P=0.047), and nearly associated with LDH (P=0.055), and CPK (P=0.078), but not associated with ALT (P>0.5). Among markers of inflammation, SC5b-9 levels were nearly associated with ferritin (P=0.057), but not with CRP, IL-6 and neutrophils/lymphocytes ratio (p>0.5).

The impact of genetic predisposition on complement activation at baseline (hospital admission) is shown in Figure 1 and in Table 2 . Plasma levels of both C5a and SC5-b9 were significantly higher in carriers of rs11385949 GA than in non-carriers (P=0.041 and P=0.012, respectively) and were higher in patients of non-O group than in patients carrying the O group, but only C5a reached statistical significance (P=0.019), as shown in Table 2 . The J o u r n a l P r e -p r o o f association between rs11385949 GA and SC5b-9 remained significant after adjustment for ABO and disease severity (P=0.004) and after further correction for C5a levels (estimate 66±65, P=0.018 for rs11385949 GA allele; estimate 11±5, P=0.023 for C5a levels; other predictors were not significant). The combined effect of rs11385949 GA allele and O blood group on complement activation products (C5a and SC5b-9) is presented in Figure 2 . In the subset of patients where the information was available the impact of rs11385949 and ABO blood group on complement activation disappeared at follow-up (at 30 days after admission during remission of symptoms in surviving patients) with a larger decrease in at risk genotypes for C5a (Table 2 ).

Overall, these data show that complement activation, which parallels COVID-19 severity, is associated with rs11385949 GA that is a recognized risk factor for the development of severe respiratory involvement.

In this study, we examined whether the genetic predisposition to severe COVID-19 with respiratory failure, in particular the rs11385949 G>GA variant at the chromosome 3 gene cluster [11] [12] [13] , influences complement activation in hospitalized patients. We were prompted by evidence that: a) the mechanism linking rs11385949 with respiratory failure in individuals infected by SARS-CoV-2 is still unclear, but may encompass modulation of the recruitment/activation of inflammatory cells [11] , b) the complement cascade is activated during COVID-19 paralleling disease severity [21, 22, 25, 26] .

The main finding was that rs11385949 was associated with increased circulating C5a and SC5b-9. In line with the hypothesis that the mechanism linking chromosome 3 variation with the susceptibility to develop more severe COVID-19 encompasses complement activation in response to SARS-CoV-2 infection, SC5b-9 levels correlated with upper respiratory airways viral load specifically in carriers of the rs11385949 GA risk allele. The loss of association during remission in a patients' subset suggests that the rs11385949 GA risk allele is associated with complement activation only when the SARS-CoV-2 infection is present. Remarkably, SC5b-9 levels correlated with biomarkers of tissue damage, e.g.

biomarkers of necrosis such as LDH, but not with classic markers of systemic inflammation, suggesting that activation of the complement cascade may contribute in an independent manner in determining lung damage. This hypothesis is consistent with experimental evidence that SARS-CoV-2 Spike protein is able to trigger complement activation [24] , whose inhibition prevents severe lung injury in mouse models of SARS coronavirus infection 1 [27, 28] .

In this study, terminal complement complex was also associated with higher GGT, possibly reflecting biliary damage, as biliary cells have been highlighted as a possible target J o u r n a l P r e -p r o o f for SARS-CoV-2 infection due to high expression of ACE2 [29] . However, the mechanism whereby chromosome 3 cluster variation facilitates complement activation during SARS-CoV-2 infection remains to be determined. Although we did not observe an impact of rs11385949 GA on upper respiratory airways viral load, it cannot be ruled out that the variant may facilitate viral replication in lower respiratory airways or other organs. However, we found that in carriers of the rs11385949 GA risk variant, viral load parallels not only complement activation but also ferritin, reflecting dysregulation of the immune response with increased monocyte/macrophage activation [30] . As rs11385949 tags genetic variation of a large gene cluster encoding several chemokine receptors directly influencing their expression (e.g. by upregulating CCR2), it could be speculated that it affects the recruitment of inflammatory cells in the lower airways or in vascular tissues, but additional studies are required to clarify the molecular pathways.

Interestingly, in the present cohort we found that rs11385949 GA risk allele carriers developed COVID-19 in the presence of a less severe burden of cardiometabolic comorbidities, a major risk factor for this condition [9, 10] . Although complement activation, endothelialithis, systemic inflammation and thrombotic complications increase together with the severity of COVID-19 and are all induced in patients with the most severe disease manifestations [8, 25] , this observation highlights a potential source of heterogeneity among patients in the initial trigger of this process culminating with respiratory failure. In parallel, we examined the impact of ABO blood group, as group O was associated with reduced risk of COVID-19, whereas the protection towards more severe outcomes in hospitalized patients is still disputed and may be modulated by genetic modifiers [11, 12, [18] [19] [20] . Interestingly, we found that group O is associated with reduced C5a and SC5b-9 levels, but only C5a reached statistical significance. It could be therefore speculated that patients with non-O group are more prone to an exaggerated activation of complement, possibly leading to a progression of tissue damage towards more severe lung disease.

Although we were able to correlate genetic background with a biomarkers panel in a well characterized cohort, limitations of this study encompass the retrospective monocentric design, limited sample size, lack of evaluation of individuals of non-European ancestry, and lack of assessment of the impact on an extended panel of cytokines and effector molecules. Independent replication in larger cohorts should therefore be obtained before ruling out a consensual minor impact of genetic predisposition on inflammatory markers and being able to draw reliable conclusions.

In summary, we found that the rs11385949 G>GA variant, tagging the chromosome 3 gene cluster variation responsible for severe COVID-19 predisposition, is associated with complement pathway activation. Indeed, rs11385949 impacted on both C5a and SC5b-9, .0001 in the whole series for a decrease at follow-up of C5a and SCC5b-9, respectively. Figure 1 . Impact of genetic predisposition to COVID-19 on complement activation. Impact of rs11385942 G>GA and ABO on circulating C5a (panel A and B, respectively) and SC5b-9 (panel C and D, respectively). At generalized linear models unadjusted (P values), and multivariate models considering rs11385942, ABO blood groups, and disease severity (adj P). J o u r n a l P r e -p r o o f 

Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo'

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

A Novel Coronavirus from Patients with Pneumonia in China

Clinical Characteristics of Coronavirus Disease 2019 in China

Baseline Characteristics and Outcomes of 1591 Patients Infected With SARS-CoV-2 Admitted to ICUs of the Lombardy Region

Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients

Extrapulmonary manifestations of COVID-19

Pulmonary Vascular Endothelialitis, Thrombosis, and Angiogenesis in Covid-19

Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area

Obesity: A critical risk factor in the COVID-19 pandemic

Genomewide Association Study of Severe Covid-19 with Respiratory Failure

Trans-ethnic analysis reveals genetic and non-genetic associations with COVID-19 susceptibility and severity. medRxiv

Genetic mechanisms of critical illness in Covid-19

The major genetic risk factor for severe COVID-19 is inherited from Neanderthals

COVID-19 Network: the response of an Italian Reference Institute to research challenges about a new pandemia

Genetic mechanisms of critical illness in Covid-19

Genetic insight into COVID-19-related liver injury

Association of ABO blood group and secretor phenotype with severe COVID-19

Reduced prevalence of SARS-CoV-2 infection in ABO blood group O

The association of ABO blood group with indices of disease severity and multiorgan dysfunction in COVID-19

Association of COVID-19 inflammation with activation of the C5a-C5aR1 axis

Immune complement and coagulation dysfunction in adverse outcomes of SARS-CoV-2 infection

Complement activation in patients with COVID-19: A novel therapeutic target

Direct activation of the alternative complement pathway by SARS-CoV-2 spike proteins is blocked by factor D inhibition

Complement activation and endothelial perturbation parallel COVID-19 severity and activity

Complement C5 inhibition in patients with COVID-19 -a promising target?

Blockade of the C5a-C5aR axis alleviates lung damage in hDPP4-transgenic mice infected with MERS-CoV

COVID-19 and liver disease

Viral RNA load in plasma is associated with critical illness and a dysregulated host response in COVID-19