To the Editor,

Recent studies in UK Cohort [1, 2] documenting increased COVID-19 incidence and deaths in populations homozygous for apolipoprotein E (ApoE) ε4/ε4 isoforms compared to ε3/ε3 with heightened significance is apprehensive. Furthermore, increased susceptibility of ε4/ε4 brain cells to SARS-CoV-2 infection in vitro [3] define its possible causal association for COVID-19 outcomes including neurological defects. However, distribution of ApoE isoforms vary amongst populations and their definitive robust association beyond UK cohort is still unknown. This tempted us to conduct a world-wide genetic association study of ApoE isoforms with COVID-19 outcomes.

Literature search on PubMed and Google scholars was conducted to retrieve ApoE isoforms data from healthy control across different countries of global population from relevant published papers. Papers with insufficient data and/or country name were excluded. Available data were subjected to Hardy Weinberg Equilibrium (HWE). Data showing deviation were excluded from analysis. Data pooling was carried out when more than one data sets were available for a particular country. Country wise frequency of all three ApoE isoforms (ε2, ε3 and ε4) were determined. COVID-19 data (incidence and deaths per millions of population) accessed on 16th January, 2021 from worldometer site (https://www.worldometers.info/coronavirus/) were associated with frequency distribution of different ApoE isoforms by Spearman’s correlation coefficient analysis in GraphPad Prism (version 5.0). Data of 16th January, 2021 (latest available before massive vaccination drive) was used to nullify the possible confounding effect due to vaccination disparity across the countries. A p value < 0.05 was considered significant.

ApoE data of 103,414 healthy controls from 228 studies belonging to 66 countries were finally enrolled which showed varied ranges of distribution for ε2 (0.625–19.1%), ε3 (52.09–92.57%) and ε4 (4.45–28.59%) (S1). Correlation analysis revealed significant but weak inverse relationship of ε4 allele with COVID-19 susceptibility (p = 0.033, Spearman’s r = − 0.263). Correlation analysis of ApoE alleles with COVID-19 morbidity and mortality was conducted separately for populations of different genetic ancestry (such as Asia, Africa, America and Europe). Consistently, ε4 allele was negatively correlated with COVID-19 outcomes (for incidence: p = 0.0181, r = − 0.6; for deaths: p = 0.0078 and r = − 0.6571) in Asia, only. Interestingly, a significant positive correlation was observed between ε3 allele frequency and COVID-19 outcomes (for incidence: p = 0.0045 and r = 0.689; for deaths: p = 0.0011 and r = 0.757) in Asia, whereas in Europe, this association reached to significant level with mortality (p = 0.048, r = 0.383) only. None of the alleles were found to be correlated with COVID-19 outcomes in Africa and America (Table 1).

Table 1 Association of ApoE isoforms with COVID-19 outcomes

Contrasting to UK cohort [1, 2], we did not observe ε4 isoform as the risk factor of COVID-19 even in Europeans. Since, UK cohort were older age population and had comorbidities, they were more vulnerable to COVID-19. Furthermore, complete exclusion of individuals having one or more of the comorbidities being not done, COVID-19 risk in ε4/ε4 should not be proclaimed as independent to pre-existing ApoE e4 associated comorbidities. Our study which includes apparently healthy control from across the world of all age groups negates the impact of comorbidity, though age groups stratified association of ApoE isoforms with COVID-19 could not be performed as age distribution of genotypes data were not available. Moreover, our findings of inconsistency in associations of ApoE variants with COVID-19 outcomes should not be surprised, since the distribution of ApoE isoforms widely vary across populations and that functionality of these isoforms are shown to be dependent on their race/ethnicity specific linkage disequilibrium (LD) and complex interactions with other genetic variants in ApoE clusters [4]. The underlying differences in association of ApoE isoforms in the present study can be explained by similar population specific heterogenous effects of ApoE isoforms. Therefore, definitive association of ApoE isoforms towards COVID-19 outcomes is challenging and seems to be more complex. The fact that angiotensin-converting enzyme 2 (ACE2) (an important host receptor for SARS-Cov2 infection and one of the several determinants for population variability of COVID-19 outcomes) and ApoE co-express in lungs alveolar cells [5], the main target for SARS-Cov2 infection. Interestingly, presence or absence of ACE2 in ApoE knockout mice is described to modulate ApoE related altered physiology [3, 6, 7]. However, such connection between ACE2 and ApoE are yet to be established in COVID-19. We hypothesize that population specific variable expression of ACE2 and other relevant genes in ApoE clusters may affect association of ApoE isoforms towards COVID19 outcomes. Therefore, well designed observational studies across various populations of different genetic backgrounds are sought to establish a definitive role of ApoE in COVID19.