Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

HMOX1 genetic polymorphisms and outcomes in infectious disease: A systematic review

  • Fergus W. Hamilton ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Project administration, Writing – original draft

    Fergus.hamilton@bristol.ac.uk

    Affiliations MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, Infection Sciences, North Bristol NHS Trust, Bristol, United Kingdom

  • Julia Somers,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliation Knight Cancer Research Building, Oregon Health and Sciences University, Portland, Oregon, United States of America

  • Ruth E. Mitchell,

    Roles Methodology, Writing – review & editing

    Affiliation MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom

  • Peter Ghazal,

    Roles Methodology, Supervision, Writing – review & editing

    Affiliation System Immunity Research Institute, Division of Infection and Immunity, Cardiff University, Cardiff, United Kingdom

  • Nicholas J. Timpson

    Roles Project administration, Resources, Supervision, Writing – review & editing

    Affiliation MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom

Abstract

Introduction

Heme-oxygenase 1 (HMOX1) is a critical stress response gene that catalyzes the multistep oxidation of heme. A GT(n) repeat of variable length in the promoter in has been associated with a wide range of human diseases, including infections. This paper aims to summarise and systematically review associations between the length of the HMOX1 GT(n) promoter and infectious disease in humans.

Methods

A search using relevant terms was performed in PubMED and EMBASE through to 15/01/21 identifying all research that studied an association between the HMOX1 GT(n) repeat polymorphism and the incidence and/or outcome of any human infectious disease. Citations were screened for additional studies. Potential studies were screened for inclusion by two authors. Data was extracted on allele frequency, genotype, strength of association, mechanism of genotyping, and potential biases. A narrative review was performed across each type of infection.

Results

1,533 studies were identified in the search, and one via citation screening. Sixteen studies were ultimately included, seven in malaria, three in HIV, three in sepsis, and one each in pneumonia, hepatitis C, and acute respiratory distress syndrome (ARDS). Sample sizes for nearly all studies were small (biggest study, n = 1,646). Allelic definition was different across all included studies. All studies were at some risk of bias. In malaria, three studies suggested that longer alleles were associated with reduced risk of severe malaria, particularly malaria-induced renal dysfunction, with four studies identifying a null association. In sepsis, two studies suggested an association with longer alleles and better outcomes.

Conclusions

Despite the importance of HMOX1 in survival from infection, and the association between repeat length and gene expression, the clinical data supporting an association between repeat length and incidence and/or outcome of infection remain inconclusive.

Introduction

HMOX1, also known as heme oxygenase is a gene that encodes for the protein HMOX1. It has been established to have a critical role in cellular stress and is ubiquitous in living organisms [13]. HMOX1 catabolises free heme into equimolar amounts of Fe2+, carbon monoxide (CO), and bilverdin, which is subsequently modified into bilirubin. Multiple in-vitro and in-vivo studies have established the toxicity of free heme and its central iron, which can lead to free radical production via Fenton chemistry, provoke excessive inflammation, and induce programmed cell death. Equally well established is the ability of HMOX1 to protect against heme-induced toxicity [1, 46].

It is also clear that the role of iron, heme, and HMOX1 in human infection is complex. Iron is critical for all prokaryotic and eukaryotic life, with pathogens and hosts battling to control iron, while avoiding the ramifications of its toxicity [7]. Multiple studies based on model organisms have identified that the presence HMOX1 is critical to defend against certain infections, with some data supporting experimental upregulation of this enzyme being protective in animal sepsis and malaria models [4, 5, 813].

In humans, the HMOX1 gene has a short tandem GT(n) repeat (STR) in its promoter region, which varies from around 20 to 40 repeats. In multiple in vitro and in-vivo studies, the length of this repeat has been shown to alter HMOX1 expression which typically occurs in response to cellular stresses [1421]. Recent work has suggested that this promoter may, in fact, be intronic and has tried to elucidate the mechanism of increased transcription. Possible mechanisms include the formation of Z-DNA, or alteration of transcription factor binding [22]. In general with larger numbers of repeats there is reduced expression, although this analysis is complicated by experimental design and classification of repeat length.

Multiple studies (>200), have investigated the impact of this STR on clinical outcomes across a broad range of human diseases, last formally reviewed in 2004 [23]. In that review, multiple signals of benefit were identified, suggesting carriage of a shorter allele might be beneficial across abroad range of conditions. Since then, meta-analyses have confirmed an association with incidence of neonatal hyperbilirubinaemia (increased with shorter length) [24], type II diabetes (increased with longer length) [25], chronic obstructive pulmonary disease (increased with longer length) [26], with an uncertain relationship with the incidence of cancer [27].

Although data is limited on the background of this STR, where data is available, the length of the repeat appears to be dramatically different in African populations, who tend to have longer repeats than many other populations, often with a trimodal distribution, rather than bimodal [19, 28, 29]. Some authors have expressed the possibility this may be due to the selective pressure of malaria, although this remains unproven [30]. As much in-vitro and in-vivo work suggests there is a link between HMOX1 activity and outcomes in infection, we aimed to systematically review the literature on the HMOX1 GT(n) repeat polymorphism and outcomes in infection.

Methods

Aim, scope, and reporting guidelines

The aim of this review was to systematically describe clinical studies on the GT(n) repeat promoter polymorphism in HMOX1 in the incidence and prognosis of infectious disease in humans. This review was prospectively registered with PROSPERO (CRD42021227072). This study was conducted in line with the HuGE guidelines on human genetic association studies [31].

Search strategy and identification of relevant papers

A search was performed across PubMED and EMBASE from inception through to 15/01/2021. The search strategy is available in the S1 File and consisted of terms linked to the gene (e.g. HMOX1) and the polymorphism (e.g. repeat, microsatellite). Abstracts were then screened by two reviewers (FH, JS) according to the criteria above, and selected for full text review. References of relevant articles were reviewed for other potential papers. Disagreements were resolved by consensus.

Inclusion criteria and eligibility (participants, intervention, exposures, controls)

Studies were included if they were studies on human participants that measured

  1. The length of the GT(n) repeat in the HMOX1 promoter using a modern molecular biological technique (e.g. fragment analysis, Sanger sequencing)

And

  1. any measure of infectious diseases outcome for individuals (e.g. incidence of an infection, prognosis of an infection, mortality within an infection)

All study types (e.g. case-control, cohort, or biobank) were included. Studies that reported on a condition that was largely due to infection (e.g. ARDS, shock in intensive care units) were also included.

Studies were excluded if they reported only rare alleles (MAF <1%) or did not report sufficient information to allow a formal review.

Analysis methodology

Given the limited number of studies available for inclusion in this analysis and differences in incidence, diagnostic methodology, and outcomes, formal meta-analysis was not possible. Instead, a narrative review of all studies was included. Data was extracted on author name(s), year of publication, location of study, ethnicity of study participants, matching of control groups, allele and genotype frequencies, consistency with the Hardy-Weinberg equilibrium, and details of laboratory methodology for HMOX1 GT(n) repeat genotyping.

As not all data was reported in studies, data was extracted from graphs where possible and re-analysis of the original data was performed using the R package metaDigitise. This allowed calculation of allele and genotype frequencies and comparisons between studies.

Risk of bias assessment

There is no validated tool for risk of bias assessment in genetic association studies. Therefore, a narrative assessment was made of the risk of bias, and this is included in the results text. Secondly, the Q-Genie tool was used for each study and the summary risk of bias figure calculated [32].

Results

A PRISMA flow diagram of included studies is shown in Fig 1. 1,533 studies were identified using the search strategy, of which 970 were unique studies. 193 of these abstracts were flagged for full text review. Full text review identified 15 full length articles included from the search. Additional citation searching and informal review identified one additional paper, and 16 papers were ultimately included in the review.

thumbnail
Fig 1. PRISMA Flow chart describing flow through the study.

https://doi.org/10.1371/journal.pone.0267399.g001

Table 1 briefly describes the 16 studies. All were case-control or cohort studies. Seven of them were in malaria, three in HIV, three in sepsis, and one each in pneumonia, hepatitis C, and acute respiratory distress syndrome (ARDS) (predominately due to sepsis, and so categorised with that). Categorisation of repeat length in each study population is shown in Fig 2, while Fig 3 shows for reference, the population distribution of this repeat in the 1000 Genomes data.

thumbnail
Fig 2. Cut offs from each study included in this review (Seu et al. did not use a cut off but fitted a linear model).

https://doi.org/10.1371/journal.pone.0267399.g002

thumbnail
Fig 3. Population distribution of this repeat, both inside and outside Africa.

Repeat called on 30X WGS data and provided kindly by Nima Mousavi, across 2,308 individuals. Probability density shown rather than raw numbers [33].

https://doi.org/10.1371/journal.pone.0267399.g003

Most studies (8/16) were in cases of infection only and compared risk of severe disease versus uncomplicated disease. Given the heterogeneity between infections, data from each infection is reported individually below.

Overall risk of bias

Each study is discussed separately below, but a risk of bias estimate was calculated using the Q-genie tool. Using this tool, all studies were of moderate (11 studies), or poor (5 studies) quality. The supplement gives the breakdown of the Q-genie score for each included study.

Malaria

Seven studies from across multiple continents focused on malaria. Most focused on severity of malaria, comparing uncomplicated/mild malaria with severe malaria, and included only patients with malaria. Interpretation of the data is confounded differing ethnic backgrounds, differing infective organism (P. vivax vs P. falciparum), differing case definitions, differing definition of repeat length and critically, differing definitions of severe disease. As most were performed in cases only, without controls, they were limited to identifying predictors of severe malaria, and were potentially biased by selecting on case status.

Takeda et al. studied a cohort (n = 150) of patients with P. falciparum malaria in Myanmar, and identified a strong association between carriage of the short allele (repeat length <28 defined as S) and cerebral malaria (OR 3.14 95% CI 1.32–7.49) [34]. In agreement with this, Walther et al. (n = 307) also found an association between the repeat length and severe malaria (all P. falciparum), identifying that longer allele carriers (repeat length > 32 as long) had markedly decreased odds of various severe malaria outcomes (OR for all severe outcomes 0.47, 95% CI 0.29–0.75) [30]. There appeared to be heterogeneity in protection from various severe outcomes, with an OR of 0.71 (95% CI 0.21–2.42) for severe anaemia, but an OR of 0.16 (95% CI 0.05–0.46) for severe renal disease. Interestingly, nine out of the ten deaths were in patients without a long allele, despite them only making up 53.4% of the severe cases.

Sambo et al. recruited a cohort of children with P. falciparum malaria (n = 430) and a cohort of healthy children (n = 319) in Angola [35]. Although not the focus of that study, the authors did measure the repeat length in a subset of patients and compared across severe malaria phenotypes and uninfected controls, using a novel classification approach (<24, 24–28, >28–34, >34). In this approach, patients with cerebral malaria had short (<24 repeats) alleles than those without malaria (46/184 alleles, 25% vs 62/211 14.7%, p = 0.0025) and those with uncomplicated malaria (40/274 alleles, 14.6%, p = 0.006). However, in severe non-cerebral malaria, the relationship was less clear, with a short allele frequency of 18% (50/178), differing little from the uninfected controls and those with uncomplicated malaria. Interpretation of this study is complicated by only a proportion of the patients being typed for the repeat, leading to a potential for bias due to patient selection. Secondly, the classification system used, as with the other studies reported, was not prespecified and varied between study.

In contrast to these positive associations, two cohort studies enrolling hospitalised patients with malaria from Kuesap [36, 37] (n = 486; n = 100), both from South East Asia, found no evidence of an association with repeat length (short <27, long > = 27) and severity of malaria (63% P. vivax, 35% P. falciparum, 2% mixed). In the smaller study [36], 12/85 (14%) of patients with an short allele developed severe malaria, with 7% (1/13) patients with no S allele developing severe malaria. In the larger study [37], no data was presented on the association, with the authors simply stating the association was null.

In another null study, Hannson et al. described a cohort of Ghanian children (n = 296) with various severities of malaria (all P. falciparum) and compared them with both uninfected controls and those with asymptomatic malaria [38]. Alleles were classified based on the number of repeats as short (<27), medium (27–32), or long (>32). They did not identify an association with severity of disease (allele and genotype frequencies similar across all severities of malaria), or with presence of infection (uninfected controls had similar genotype frequencies to infected patients). Althoughprimary analysis data was not provided, the authors state that all models had p >0.7.

In a similar study in Brazil, Mendonca et al. presented data on a cohort of patients (n = 264) with asymptomatic malaria well as patients with symptomatic malaria and uninfected controls [39]. More than 90% of cases in this study were of P. vivax, rather than P falciparum, complicating interpretation with other studies. Interestingly, using a binary classification (short <30, L > = 30), they identified an association with longer GT(n) repeats and symptomatic malaria (OR compared to asymptomatic/uninfected controls 3.35, 95%CI 1.91–5.58). In contrast to other studies, associations were also identified with biochemical severity of hepatic inflammation with longer repeats, although the clinical phenotypes of severe malaria were not reported.

In summary, most of the studies focused on severity of disease in patients presenting with malaria. Three of the studies on associations with severity of disease and GT(n) repeat length identified strong (and similar in magnitude) associations with shorter repeats and worse outcomes; while two studies reported null outcomes, with all of these studies except one entirely performed on patients with P. falciparum. It is important to note that the smaller of the two null studies (Kuesap et al.) [36] actually reported a strong association with carriage of an shorter allele and development of severe malaria (7% with longer allele, 14% without short allele), and the lack of association may reflect lack of statistical power due to low case numbers.

It is also important to recognise the differing outcome measures across the studies. The strongest association in Walther et al. was for severe renal disease with an OR of 0.16 (95% CI 0.05–0.46) with carriage of a longer allele. This outcome was not included in the Hansson paper, with comparisons limited to cerebral malaria and severe anaemia. In summary, the data so far suggest there is a possible association with carriage of a longer repeats and protection from severe malaria, although the data remains weak due to a small number of studies with relatively low numbers, and the differing definitions of exposure, that vary greatly across studies.

For association with incidence of infection; the data are less clear. In the Hannson study, including patients with P. falciparum, a null result was reported for association with any allele carriage and incidence of asymptomatic or symptomatic malaria. In contrast, Mendonca et al. identified a strong relationship between symptomatic malaria (>90% P. vivax) and GT(n) repeat length, with an OR of 3.35 compared to asymptomatic cases or uninfected controls (95%CI 1.91–5.58). Given the paucity of data, differing settings, and different infective pathogen (P. vivax vs P. falciparum, the association with the GT(n) repeat and incidence of infection remains open. Given the potential relationship between both incidence and severity, the interpretation of the case only studies (the majority of studies) is even more difficult, as the relationship between severity may be biased if there is also an effect on incidence.

Sepsis and critical infection (ARDS)

Four studies reported associations with GT(n) promoter repeats and outcomes in sepsis/critical infection, with one study enrolling controls also. First, Sponholz et al. reported a cohort of patients (n = 430) with severe sepsis or septic shock from Germany [22]. This cohort had a mean age 65 and were critically unwell with a median APACHE II 20.2. The most common source of infection was respiratory infection. Alleles were classified as short (<27 repeats), medium (27–34 repeats), and long (>34 repeats). In this study, genotype had an association with 28-day mortality, but not 90-day mortality, in a non-genetic model. The highest mortality was in M/M genotype patients (48/155, 31%), with the lowest in (S/M genotypes (26/181, 14.4%). However, the most striking finding in this study was the complete absence of the L/L genotype and relative rarity of the M/L genotype (1/430). In fact, L alleles were rare, with only 16 alleles out of 860 (1.9%). Although the normal distribution of this STR is not yet established, our analysis suggests that the frequency of alleles with repeats of length >34 in Europeans in UK Biobank is 6.6%, while analysis of the Central European population in 1000G suggests a frequency of 6.0% [33]. This raises concern about a relationship between the repeat length and incidence of disease, and is open to bias toward cases. Finally, although the authors did not fit a genetic model, the relationship is inconsistent with both dominant and additive models, with the lowest mortality in heterozygotes (S/M) model, and highest in homozygotes (S/S, M/M).

In a much smaller cohort in Mexico, 64 septic children, 72 hospitalised control children and 101 healthy adults were genotyped for the HMOX1 repeat [40]. Using a binary classification (short <25 /long > = 25 repeats), there was no strong evidence for a difference between the allele frequency or genotype between the groups. In comparison to the above study, L alleles were very common, comprising 75% of all the alleles, but this likely reflects the cut-off chosen (S <25 repeats), which is much lower than in other studies. Given that hospitalisation in children is rare, the choice of control group of hospitalised children is difficult, as these children are likely to have other significant medical illnesses, rather than representing a population control.

Vilander et al. enrolled 653 critically ill patients with sepsis from Finland, 300 of whom had severe acute kidney injury (AKI) [41]. On average, patients were middle aged (median age 63) and were nearly all emergency admissions (98%). Using the same S/M/L classification as Sponholz et al. (<27, S, 27–34 M, >34, L), the presence of any L allele was rare (63/1306, 4.9%), and the presence of an L/L genotype was very rare, with only 1 patient having the L/L genotype. In their additive model the presence of an S-allele was associated with an increased risk of acute kidney injury (OR 1.29, 95%CI 1.02–1.64). As all subjects in this study were cases, this suffers from the same potential biases, given the potential association with being a case.

Finally, Sheu et al. studied a U.S. cohort of critically ill patients with acute respiratory distress syndrome (ARDS) risk factors and examined the association with the repeat and ARDS development [42]. 1,451 patients were included, the vast majority of whom had infection (sepsis, septic shock, and pneumonia). It is important to note that the types of infection differed between the cohort that developed ARDS and those that did not, with much higher rates of septic shock and pneumonia in the group that went on to develop ARDS (60% vs 44%, and 68% vs 43%, respectively). Like the other studies from outside Africa, using an S/M/L classification (S <24, M 24–30, L >32), L alleles and the L/L genotype were rare (355/2902, 2% L alleles; L/L genotype 1.1%, 17/1451). In the allelic model, the presence of an L allele was associated a reduced risk of ARDS (any L allele, OR 0.63, 95%CI 0.47–0.85). As with the other two studies in sepsis, the case only and the relative rarity of some genotypes compared to the expected population frequency raises concern of bias.

The data on sepsis is less compelling than in malaria, with the one major study in septic shock finding a weak association with a single genotype (M/L), in a non-genetic model. However, Vilander et al. provide data on the protection from renal dysfunction with longer alleles, which should be taken in context with data from malaria, where the one study reporting renal outcomes showed a similar protective effect of the L allele. The interpretation of the final paper by Sheu et al. is difficult, as although many (most) patients in the study had infection, the cohort was not entirely comprised of patients who had infection, and the source and type of infection differed between the two groups. All studies were performed in cases, with the potential biases as described above.

HIV

Three studies have examined the association between HMOX1 promoter repeats and outcomes in people living with HIV (PLWH), with quite differing results to that in sepsis and malaria. The first, a U.S. study from Seu et al. [19], examined a cohort (n = 717) of PLWH of differing ethnicities. Firstly, they showed the distribution of repeats was different in African Americans vs Caucasians (trimodal vs bimodal), with a longer repeat length on average in African Americans. In a smaller section of the cohort not on antiretroviral therapy longer repeats were associated with an increased viral load in African Americans (n = 74, r = 0.24, p = 0.04), but not in Caucasians (n = 177), using a linear regression model rather than a classification system. Finally, they demonstrated an increase in sCD14, a marker of increased mortality in PLWH on antiretroviral therapy [43] with increasing repeat length (n = 50, r = 0.38, p = 0.07).

Gill et al. genotyped an autopsy cohort from the US (n = 554) [28]. In this cohort, PLWH with and without HIV encephalitis were compared with HIV negative controls for a variety of neuroinflammatory cell markers. Like many other studies, they identified the increased GT(n) repeat length in patients of African origin, with L alleles (defined as >34 repeats) being rare in Caucasians.

In contrast to Seu et al., a likely null association was found between HIV viral load and GT(n) repeat status in both plasma and CSF, however, there was a strong association between shorter repeats (S<27, L>34) and reduced risk of HIV encephalitis (OR for HIVE with S allele 0.62, 95%CI 0.39–0.98). Shorter repeats were also associated with reduced activation of type 1 interferon response markers (MX1, ISG15, and IRF1) and T-lymphocyte response markers (CD38 and GZMB). However, a likely null association was identified with other cellular markers. It is important to recognise the biases inherent in selecting on an autopsy cohort, who may not be representative of the general population. Supporting that, the mean age of death was 54.5 years in patients without HIV, and only 42 years in patients with HIV.

Finally, Garza et al. (from the same group as Gill), genotyped a cohort of US (n = 603) and Botswanan (n = 428) people living with HIV, and tested the association with HMOX1 repeats and neurocognitive impairment. Alleles were defined as short (<27 repeats), medium (27–34 repeats), and long (>34 repeats). Like the other previous cohorts, the distribution of repeat lengths was longer in those with African heritage compared to European (presence of L allele, defined as >34 repeats; 33.3% vs 3.6%), and was more convincingly trimodal. Like Gill, a null association was identified with HIV viral load and with lymphocyte count. However, an association between the presence of short alleles and reduced risk of HIV neurocognitive impairment was identified (OR with at least one short allele 0.63; 95% CI 0.28–0.84), which was robust to definition of neurocognitive impairment. Interestingly, this association was only identified in the African American cohort, not the European cohort, where the association was not significant. The test for heterogeneity was not met, although this is likely due to the small sample size.

The data from HIV is difficult to interpret, and at significant potential of bias. One of the three studies showed a weak association with viral load in African Americans only, while the two other cohorts found a null association. Similar to previous reviews [23] (but counter to the data in malaria, and perhaps sepsis), the shorter repeat appears to be protective against chronic inflammation, which is felt to be the cause of HIV encephalitis and neurocognitive impairment.

Importantly, one constant feature of all three studies is the apparent heterogeneous effect of the repeat depending on differing background populations. Whether this represents simply the reduced power (as the longer repeats are relatively rare in Caucasian populations), related to differing genetic backgrounds, or a true biological heterogeneity remains to be explored.

Pneumonia

One study looked for an association between HMOX1 promoters and pneumonia. Yasuda et. al. compared 200 elderly patients with pneumonia (median age 74) to matched controls. In the Japanese population, the L allele remains rare (15% of all alleles). In an allelic model, carriage of an L-allele was associated with increased odds of pneumonia (OR 2.3; 1.5 to 3.5). No data was provided on the outcome of infection.

Hepatitis C

A single study by Urbanek et. al. in the Czech Republic compared HMOX1 repeat lengths within Hepatitis C Virus (HCV) infected patients (n = 146) and compared them to healthy controls (n = 146). In this study, L alleles were defined as >29. Like other European studies, the L allele was rare (7% in both patients and controls). There was a null association between the genotype and case status, although the study was small. There was a null relationship between biochemical markers of liver inflammation in cases and patient genotype. Finally, no association was identified between HCV viral load and HMOX1 promoter length.

Discussion

Although there are now 15 studies in human infection that have aimed to assess the impact of the GT(n) repeat polymorphisms in HMOX1, the clinical data remains sparse and difficult to interpret. All studies were underpowered for small effects and at significant risk of bias, with the largest study comprising only 1,451 patients, while allelic definition varied across almost every study. Fig 2 shows this graphically, overlaying the study definitions included in this review on the background repeat length across populations. As can be seen, the allelic definition is highly variable, and allele frequencies are highly sensitive to both population structure and allele definition. Therefore, based on the data collected so far, we cannot conclude whether HMOX1 polymorphisms have a role in infectious disease.

With that caveat, the most promising data supporting an association lies in P. falciparum malaria and, perhaps, sepsis. In these studies, there is weak evidence in support of a potential association between longer repeats of the HMOX1 gene and better outcomes, with a particularly compelling relationship between worsening renal dysfunction and shorter alleles.

Interestingly, two papers outside the context of infection (one in sickle cell disease, the other post cardiac surgery) have identified the opposite association, with shorter repeats being associated with better renal outcomes, with similar allelic definitions [44, 45]. This heterogeneity of effect in differing conditions with a common pathology (renal dysfunction) underlines the complexity of both the underlying pathology, and perhaps, the sensitivity of results to study design.

The other major issue across many of these studies is the focus on cases only. In studies performed on cases only in Europe ancestry patients the genotype frequency differs from the population genotype frequency, suggesting there may well be associations with incidence as well as severity. If there is an association with incidence, as well as severity of disease, then estimates generated from case only studies may well be biased, due to a collider bias [46].

In HIV, in contrast, carriage of shorter alleles seems to be associated with two effects. Firstly, in two studies, there was a strong association with shorter alleles and neuroprotection: both for HIV encephalitis, and for neurocognitive impairment. Interestingly, in one cohort, this association was limited to those of African heritage, although small numbers preclude interpretation of this in detail. Secondly, there was an intriguing association with longer repeats and increased viral load, identified in one study only. This remains to be explored further.

With other infections, the data is too sparse to provide any conclusions, with only one potential association in elderly patients with pneumonia.

Summary

There is a plethora of in-vitro and in-vivo evidence for the importance of HMOX1 in management of cellular stress and infection. Sixteen studies to date have interrogated the relevance of the promoter repeat polymorphism, however due to confounding factors, sample sizes and case-only biases, the answer remains unclear. The most compelling evidence for an effect lies in P. falciparum malaria, but studies are small (all n <500), and interpretation of causal effects is confounded by the large difference in background allele frequency between studies (Brazil, Gambia, Ghana, Myanmar, Thailand), the varying definitions of case severity, varying allele definitions and the selection of cases only. Although this was not a focus of a study, we did not identify any sex-specific interaction between the HMOX1 promoter repeat and sex, despite the known association between certain infections and biological sex.

Further work should focus on functional characterisation of the repeat length on gene expression, and involve much larger cohort approaches, recruiting prior to infection and assessing both the incidence and severity of infection. Given that allelelic definition loses both statistical power and may lead to bias, future studies should use repeat length as a continuous predictor, rather than using arbitrary cut offs. Secondly, where possible, gene expression data and other clinical and laboratory data should be included, to assess the potential mechanism of the repeat, as well as simply relying on case status.

Recent work has suggested that imputation of this STR from array data is possible (at least in European ancestry populations), and therefore analysis of infection data from larger samples may well be possible in the near future [33]. With this approach, large (>10,000) cohorts with detailed clinical phenotyping may be able to be genotyped for this repeat, which should provide further clarity on the role of HMOX1 promoter polymorphisms in infection.

Conclusion

The HMOX1 promoter GT(n) polymorphism has been studied across a broad range of infections, with the most studies in malaria and sepsis. There is not enough information to reliably associate the repeat with incidence and severity of any other infection.

Acknowledgments

We would like to thank Kevin Watanbe-Smith for help designing Fig 1.

References

  1. 1. Gozzelino R, Jeney V, Soares MP. Mechanisms of cell protection by heme oxygenase-1. Annu Rev Pharmacol Toxicol. 2010;50: 323–354. pmid:20055707
  2. 2. Wilks A. Heme oxygenase: evolution, structure, and mechanism. Antioxid Redox Signal. 2002;4: 603–614. pmid:12230872
  3. 3. Li C, Stocker R. Heme oxygenase and iron: from bacteria to humans. Redox Rep. 2009;14: 95–101. pmid:19490750
  4. 4. Larsen R, Gozzelino R, Jeney V, Tokaji L, Bozza FA, Japiassú AM, et al. A central role for free heme in the pathogenesis of severe sepsis. Sci Transl Med. 2010;2: 51ra71. pmid:20881280
  5. 5. Pamplona A, Ferreira A, Balla J, Jeney V, Balla G, Epiphanio S, et al. Heme oxygenase-1 and carbon monoxide suppress the pathogenesis of experimental cerebral malaria. Nat Med. 2007;13: 703–710. pmid:17496899
  6. 6. Seixas E, Gozzelino R, Chora A, Ferreira A, Silva G, Larsen R, et al. Heme oxygenase-1 affords protection against noncerebral forms of severe malaria. Proc Natl Acad Sci U S A. 2009;106: 15837–15842. pmid:19706490
  7. 7. Skaar EP. The battle for iron between bacterial pathogens and their vertebrate hosts. PLoS Pathog. 2010;6: e1000949. pmid:20711357
  8. 8. Mohammed K.A., Gahlot S., Nasreen N., Sadikot R.T., Johnson J.A., Sahn S.A. HO-1-/- mice are susceptible to MRSA empyema due to poor pleural TLR9 response. Am J Respir Crit Care Med. 2014;189. Available: http://www.atsjournals.org/doi/pdf/10.1164/ajrccm-conference.2014.189.1_MeetingAbstracts.A2498
  9. 9. Chung SW, Liu X, Macias AA, Baron RM, Perrella MA. Heme oxygenase-1-derived carbon monoxide enhances the host defense response to microbial sepsis in mice. J Clin Invest. 2008;118: 239–247. pmid:18060048
  10. 10. Willix JL, Stockton JL, Olson RM, Anderson PE, Anderson DM. Activation of Heme Oxygenase Expression by Cobalt Protoporphyrin Treatment Prevents Pneumonic Plague Caused by Inhalation of Yersinia pestis. Antimicrob Agents Chemother. 2020;64. pmid:32015027
  11. 11. Tsoyi K, Lee TY, Lee YS, Kim HJ, Seo HG, Lee JH, et al. Heme-oxygenase-1 induction and carbon monoxide-releasing molecule inhibit lipopolysaccharide (LPS)-induced high-mobility group box 1 release in vitro and improve survival of mice in LPS- and cecal ligation and puncture-induced sepsis model in vivo. Mol Pharmacol. 2009;76: 173–182. pmid:19366789
  12. 12. Singh N, Ahmad Z, Baid N, Kumar A. Host heme oxygenase-1: Friend or foe in tackling pathogens? IUBMB Life. 2018;70: 869–880. pmid:29761622
  13. 13. Epiphanio S, Mikolajczak SA, Gonçalves LA, Pamplona A, Portugal S, Albuquerque S, et al. Heme oxygenase-1 is an anti-inflammatory host factor that promotes murine plasmodium liver infection. Cell Host Microbe. 2008;3: 331–338. pmid:18474360
  14. 14. Maruyama A, Mimura J, Itoh K. Non-coding RNA derived from the region adjacent to the human HO-1 E2 enhancer selectively regulates HO-1 gene induction by modulating Pol II binding. Nucleic Acids Res. 2014;42: 13599–13614. pmid:25404134
  15. 15. Maruyama A, Mimura J, Harada N, Itoh K. Nrf2 activation is associated with Z-DNA formation in the human HO-1 promoter. Nucleic Acids Res. 2013;41: 5223–5234. pmid:23571756
  16. 16. Chen Y-H, Lin S-J, Lin M-W, Tsai H-L, Kuo S-S, Chen J-W, et al. Microsatellite polymorphism in promoter of heme oxygenase-1 gene is associated with susceptibility to coronary artery disease in type 2 diabetic patients. Hum Genet. 2002;111: 1–8. pmid:12136229
  17. 17. Hirai H, Kubo H, Yamaya M, Nakayama K, Numasaki M, Kobayashi S, et al. Microsatellite polymorphism in heme oxygenase-1 gene promoter is associated with susceptibility to oxidant-induced apoptosis in lymphoblastoid cell lines. Blood. 2003;102: 1619–1621. pmid:12730098
  18. 18. Yamada N, Yamaya M, Okinaga S, Nakayama K, Sekizawa K, Shibahara S, et al. Microsatellite polymorphism in the heme oxygenase-1 gene promoter is associated with susceptibility to emphysema. Am J Hum Genet. 2000;66: 187–195. pmid:10631150
  19. 19. Seu L, Burt TD, Witte JS, Martin JN, Deeks SG, McCune JM. Variations in the heme oxygenase-1 microsatellite polymorphism are associated with plasma CD14 and viral load in HIV-infected African-Americans. Genes Immun. 2012;13: 258–267. pmid:22048453
  20. 20. Tanaka G, Aminuddin F, Akhabir L, He J-Q, Shumansky K, Connett JE, et al. Effect of heme oxygenase-1 polymorphisms on lung function and gene expression. BMC Med Genet. 2011;12: 117. pmid:21902835
  21. 21. Doberer D, Haschemi A, Andreas M, Zapf T-C, Clive B, Jeitler M, et al. Haem arginate infusion stimulates haem oxygenase-1 expression in healthy subjects. Br J Pharmacol. 2010;161: 1751–1762. pmid:20718734
  22. 22. Kramer M, Sponholz C, Slaba M, Wissuwa B, Claus RA, Menzel U, et al. Alternative 5’ untranslated regions are involved in expression regulation of human heme oxygenase-1. PLoS One. 2013;8: e77224. pmid:24098580
  23. 23. Exner M, Minar E, Wagner O, Schillinger M. The role of heme oxygenase-1 promoter polymorphisms in human disease. Free Radic Biol Med. 2004;37: 1097–1104. pmid:15451051
  24. 24. Zhou J-F, Luo J-Y, Zhu W-B, Yang C-Y, Zeng Y-L, Qiu X-L. Association between genetic polymorphism of heme oxygenase 1 promoter and neonatal hyperbilirubinemia: a meta-analysis. J Matern Fetal Neonatal Med. 2021;34: 12–23. pmid:30700176
  25. 25. Qiao H, Sai X, Gai L, Huang G, Chen X, Tu X, et al. Association between heme oxygenase 1 gene promoter polymorphisms and susceptibility to coronary artery disease: a HuGE review and meta-analysis. Am J Epidemiol. 2014;179: 1039–1048. pmid:24670375
  26. 26. Zhou H, Ying X, Liu Y, Ye S, Yan J, Li Y. Genetic polymorphism of heme oxygenase 1 promoter in the occurrence and severity of chronic obstructive pulmonary disease: a meta-analysis. J Cell Mol Med. 2017;21: 894–903. pmid:27998018
  27. 27. Luo H, Shao Y, Yao N, Chen X, Hu L, He T. Association of heme oxygenase-1 polymorphisms with cancer risk: A systematic review and meta-analysis. J BUON. 2015;20: 1142–1153. pmid:26416069
  28. 28. Gill AJ, Garza R, Ambegaokar SS, Gelman BB, Kolson DL. Heme oxygenase-1 promoter region (GT)n polymorphism associates with increased neuroimmune activation and risk for encephalitis in HIV infection. J Neuroinflammation. 2018;15: 70. pmid:29510721
  29. 29. WebSTR: a population-wide database of short tandem repeat variation in humans. [cited 28 Apr 2021]. http://webstr.ucsd.edu/locus?STRID=STR_900865
  30. 30. Walther M, De Caul A, Aka P, Njie M, Amambua-Ngwa A, Walther B, et al. HMOX1 gene promoter alleles and high HO-1 levels are associated with severe malaria in Gambian children. PLoS Pathog. 2012;8: e1002579. pmid:22438807
  31. 31. Sagoo GS, Little J, Higgins JPT. Systematic reviews of genetic association studies. Human Genome Epidemiology Network. PLoS Med. 2009;6: e28. pmid:19260758
  32. 32. Sohani ZN, Sarma S, Alyass A, de Souza RJ, Robiou-du-Pont S, Li A, et al. Empirical evaluation of the Q-Genie tool: a protocol for assessment of effectiveness. BMJ Open. 2016;6: e010403. pmid:27288371
  33. 33. Saini S, Mitra I, Mousavi N, Fotsing SF, Gymrek M. A reference haplotype panel for genome-wide imputation of short tandem repeats. Nat Commun. 2018;9: 4397. pmid:30353011
  34. 34. Takeda M, Kikuchi M, Ubalee R, Na-Bangchang K, Ruangweerayut R, Shibahara S, et al. Microsatellite polymorphism in the heme oxygenase-1 gene promoter is associated with susceptibility to cerebral malaria in Myanmar. Jpn J Infect Dis. 2005;58: 268–271. pmid:16249618
  35. 35. Sambo MR, Trovoada MJ, Benchimol C, Quinhentos V, Gonçalves L, Velosa R, et al. Transforming growth factor beta 2 and heme oxygenase 1 genes are risk factors for the cerebral malaria syndrome in Angolan children. PLoS One. 2010;5: e11141. pmid:20585394
  36. 36. Kuesap J, Na-Bangchang K. The Effect of ABO Blood Groups, Hemoglobinopathy, and Heme Oxygenase-1 Polymorphisms on Malaria Susceptibility and Severity. Korean J Parasitol. 2018;56: 167–173. pmid:29742871
  37. 37. Kuesap J, Hirayama K, Kikuchi M, Ruangweerayut R, Na-Bangchang K. Study on association between genetic polymorphisms of haem oxygenase-1, tumour necrosis factor, cadmium exposure and malaria pathogenicity and severity. Malar J. 2010;9: 260. pmid:20846452
  38. 38. Hansson HH, Maretty L, Balle C, Goka BQ, Luzon E, Nkrumah FN, et al. Polymorphisms in the Haem Oxygenase-1 promoter are not associated with severity of Plasmodium falciparum malaria in Ghanaian children. Malar J. 2015;14: 153. pmid:25888733
  39. 39. Mendonça VRR, Luz NF, Santos NJG, Borges VM, Gonçalves MS, Andrade BB, et al. Association between the haptoglobin and heme oxygenase 1 genetic profiles and soluble CD163 in susceptibility to and severity of human malaria. Infect Immun. 2012;80: 1445–1454. pmid:22290142
  40. 40. Vázquez-Armenta G, González-Leal N, Vázquez-de la Torre M J, Muñoz-Valle JF, Ramos-Márquez ME, Hernández-Cañaveral I, et al. Short (GT)n microsatellite repeats in the heme oxygenase-1 gene promoter are associated with antioxidant and anti-inflammatory status in Mexican pediatric patients with sepsis. Tohoku J Exp Med. 2013;231: 201–209. pmid:24201221
  41. 41. Vilander LM, Vaara ST, Donner KM, Lakkisto P, Kaunisto MA, Pettilä V, et al. Heme oxygenase-1 repeat polymorphism in septic acute kidney injury. PLoS One. 2019;14: e0217291. pmid:31120979
  42. 42. Sheu C-C, Zhai R, Wang Z, Gong MN, Tejera P, Chen F, et al. Heme oxygenase-1 microsatellite polymorphism and haplotypes are associated with the development of acute respiratory distress syndrome. Intensive Care Med. 2009;35: 1343–1351. pmid:19526221
  43. 43. Sandler NG, Wand H, Roque A, Law M, Nason MC, Nixon DE, et al. Plasma levels of soluble CD14 independently predict mortality in HIV infection. J Infect Dis. 2011;203: 780–790. pmid:21252259
  44. 44. Leaf DE, Body SC, Muehlschlegel JD, McMahon GM, Lichtner P, Collard CD, et al. Length Polymorphisms in Heme Oxygenase-1 and AKI after Cardiac Surgery. J Am Soc Nephrol. 2016;27: 3291–3297. pmid:27257045
  45. 45. Saraf SL, Viner M, Rischall A, Raslan R, Shah BN, Zhang X, et al. HMOX1 and acute kidney injury in sickle cell anemia. Blood. 2018;132: 1621–1625. pmid:30139825
  46. 46. Dudbridge F, Allen RJ, Sheehan NA, Schmidt AF, Lee JC, Jenkins RG, et al. Adjustment for index event bias in genome-wide association studies of subsequent events. Nat Commun. 2019;10: 1561. pmid:30952951