Skip to main content

HMOX1 STR polymorphism and malaria: an analysis of a large clinical dataset

Abstract

Background

Inducible expression of heme oxygenase-1 (encoded by the gene HMOX1) may determine protection from heme released during malaria infections. A variable length, short tandem GT(n) repeat (STR) in HMOX1 that may influence gene expression has been associated with outcomes of human malaria in some studies. In this study, an analysis of the association between variation at the STR in HMOX1 on severe malaria and severe malaria subtypes is presented in a large, prospectively collected dataset (MalariaGEN).

Methods

The HMOX1 STR was imputed using a recently developed reference haplotype panel designed for STRs. The STR was classified by total length and split into three alleles based on an observed trimodal distribution of repeat lengths. Logistic regression was used to assess the association between this repeat on cases of severe malaria and severe malaria subtypes (cerebral malaria and severe malarial anaemia). Individual analyses were performed for each MalariaGEN collection site and combined for meta-analysis. One site (Kenya), had detailed clinical metadata, allowing the assessment of the effect of the STR on clinical variables (e.g. parasite count, platelet count) and regression analyses were performed to investigate whether the STR interacted with any clinical variables.

Results

Data from 17,960 participants across 11 collection sites were analysed. In logistic regression, there was no strong evidence of association between STR length and severe malaria (Odds Ratio, OR: 0.96, 95% confidence intervals 0.91–1.02 per ten GT(n) repeats), although there did appear to be an association at some sites (e.g., Kenya, OR 0.90, 95% CI 0.82–0.99). There was no evidence of an interaction with any clinical variables.

Conclusions

Meta-analysis suggested that increasing HMOX1 STR length is unlikely to be reliably associated with severe malaria. It cannot be ruled out that repeat length may alter risk in specific populations, although whether this is due to chance variation, or true variation due to underlying biology (e.g., gene vs environment interaction) remains unanswered.

Background

The gene HMOX1 (also known as heme-oxygenase 1) encodes for a protein HO-1 that has been shown to be critical for survival in cerebral malaria in mouse models of malaria [1, 2]. HO-1 breaks down heme (or haem) to iron, biliverdin and carbon monoxide, with the latter playing a protective role in experimental cerebral malaria [1]. In a non-coding region of the gene, there exists a short tandem repeat (STR), that shows a GT(n) polymorphism of variable length. Some but not all in-vitro studies have suggested that the length of this STR associates with expression of HMOX1, where cells containing longer STRs (more than around 35 repeats) have reduced inducible expression [2]. The STR ranges from around 25 repeats to around 40, with wide ranging differences across different populations [3].

Multiple human studies across many diseases have suggested that STR variation at this locus is associated with relevant clinical outcomes. For example, a meta-analysis in diabetes suggested an effect of this STR (longer repeats increasing the risk of type II diabetes) [4], while another meta-analysis identified longer repeats decreasing the risk of neonatal jaundice [5]. Additionally, on the basis of pathway driven basic science and given the association with mammalian malaria, multiple studies in malaria have examined the role of this STR in both the incidence of malaria and development of severe malaria (e.g., cerebral malaria) [1, 6, 7]. Despite this, evidence of an effect of STR variation on the presence of severe malaria remains inconclusive, with studies having small sample sizes ( < 1000 participants) and variable definitions of STR length. Furthermore, given the nature of STR variation and recent evolution, controlling for ethnicity in both case status and STR length is key to understanding if there is a relationship between the HMOX1 STR and malaria outcomes, which requires some adjustment for population structure which has not been performed in studies so far.

This study aimed to impute the HMOX1 STR using a recognized haplotype reference panel for STRs [8] in a large, geographically diverse case–control study of malaria (MalariaGEN [9]) in order to analyse the association between STR genotype and (a) malaria case status (b) severe malaria subtype status.

Methods

Reporting guidelines

This study was reported in line with the STREGA reporting guidelines (ref https://journals.plos.org/plosmedicine/article?id=10.1371/journal.pmed.1000022) which are available in supplement S1

Data source

This analysis used the MalariaGEN dataset and was approved by the MalariaGEN IDAC (Application: 71). Clinical phenotyping (clinical parameters, parasite count) was available for one subset of this cohort, the Kenya cohort, and was arranged via Professor Tom Williams at the KEMRI-Wellcome Unit [10].

This study was a meta-analysis of 11 separate studies on resistance to severe malaria. Each study contributed to a wider meta-analysis of resistance to severe malaria [9]. Each site had the same genetic analysis pipeline (detailed in the above publication), with analysis performed at each site, and then meta-analysed for the summary results.

HMOX1 STR genotyping and definition of STR length

This analyses used the MalariaGEN data set. MalariaGEN SNP array data was downloaded from the EGA (EGA, dataset EGAD00010001799) [9]. Imputation was performed using a recently developed haplotype reference panel based on the 1000 Genomes data using Beagle v4.2 [8] Previous work has shown this imputation is reliable and have shown no benefit of filtering this imputation based on predicted genotype probability [11]. In that study, data from the 1000 Genomes Project [12] and the Human Genome Diversity project [13] was used to compare the imputed STR length with the real STR length, as called by whole genome sequencing. In general, correlation between real and imputed STR length in worldwide populations was good (Pearson’s R > 0.8), although there was worse imputation performance in some ethnic backgrounds. Previous work provides detail on the imputation approach and quality control [11].

Previous literature on HMOX1 STR length in malaria has used a wide range of reference lengths for the HMOX1 STR [14]. Given the lack of consensus and evidence for any given definition the decision was made to use the summed repeat length as our primary definition in this study, where STR length is defined by the total number of GT repeats. For further analysis, the repeat was split into three alleles—short (S), medium (M), and long (L), using definitions below, in line with previous literature and the trimodal variation at this STR. Genotypes were subsequently defined by the combination of these alleles, e.g. SS, SM, etc., in line with previous literature, leading to six genotypes at this locus [3].

Malaria outcomes

Malaria outcomes were extracted from the MalariaGEN dataset provided by MalariaGEN alongside data on sex, ethnicity, and country. Outcomes included: case status and type of severe malaria status (cerebral malaria (CM) and/or severe malarial anaemia (SMA). Details of definitions are with the original MalariaGEN publication [9].

For one MalariaGEN collection site (Kenya), more detailed clinical phenotyping of malaria cases was available (blood pressure, platelet count, haemoglobin, MCV, parasite count, severe renal disease, and mortality). Details of definitions are with the relevant publication [10].

Statistical approach

Main outcome

The primary analytic method was logistic regression on severe malaria case status in each individual study site with the total STR repeat length as the explanatory variable. Subsequently, site-specific estimates were meta-analysed in a random effects model to generate summary estimates across the whole cohort. For each severe malaria subtype (CM, SMA, or both), a similar analysis was performed comparing cases of each severe malaria subtype vs a) all other cases, and b) vs all controls. Subsequently, an alternative repeat length definition was generated by splitting each allele across the trimodal distribution into short (S), medium (M), long (L) alleles. The allele cut points were short ( < 27 repeats), medium (27–32 repeats) and long  > 32 repeats and these alleles were applied in a logistic regression against severe malaria case status.

As a third alternative model, the actual genotype was included (e.g., short-short, short-medium, medium-long, long-long), in a logistic regression model against the reference genotype (medium-medium). These analyses were again performed for each outcome in each country and meta-analysed in a random effects model. To ensure estimates had some precision, models were only run if they included more than 30 cases, and more than 150 total patients.

The first 10 principal components and sex were included as covariates in our models. These principal components were supplied by MalariaGEN and calculated at the site level, and therefore represent genetic variation at each site, not across the whole meta-analysis.

Associations and interactions with clinical variables

One data set had available clinical data (Kenya). At this site, linear regression was performed with STR length on the clinical biomarkers of severity: Hb, platelet count, blood pressure, white cell count, MCV, and platelet count, and parasite count (logged to improve model fit).

Subsequently, an assessment was made of any potential intraction between STR length and clinical variables. In these models, logistic regression was performed on severe malaria case status with an interaction fitted between repeat length and clinical variable (e.g. case ~ STR length * platelet count). This analysis was based on laboratory data suggesting that HMOX1 variation may only have clinical impact in certain subtypes (e.g. high parasite counts) of severe malaria.

Results

SNP array and clinical data were available for 17,960 participants in MalariaGEN. All genotypes were successfully imputed. 215 parents and 3 patients had unknown case status, and 302 patients who had no country information. These were excluded from the analysis data set leading to 8658 cases, and 8931 controls. Table 1 shows the breakdown of case status and severe disease by MalariaGEN site.

Table 1 Distribution of severe malaria across the MalariaGEN cohort

Case frequency differed by country with the highest number of cases in Nigeria (83% of included participants cases), and the lowest in Cameroon (46% of included participants cases). The presence of severe malaria subtypes also differed greatly by site, with rates of severe malarial anaemia as high as 31% in Mali, and as low as 1% in Nigeria.

Ethnicity had an equally large impact on case ratio and severe malaria subtypes, with Additional file 1: Table S1 summarizing this data.

Details of the HMOX1 STR imputation process and accuracy are detailed in a recent publication on the HMOX1 STR in UK Biobank [11].

As expected, HMOX1 STR genotype differed across populations. Figure 1A shows the distribution of STR lengths across the whole of MalariaGEN, while Fig. 1B shows individual distributions of alleles stratified by case/control status and Fig. 1C shows the summed STR length across both alleles for all countries. The trimodal distribution (peaks at 29 STRs, 38 STRs, and 23 STRs) previously reported was identified, although all three peaks were not present in all subpopulations. The mean total repeat length across all cohorts was 60.9 repeats, with a standard deviation of 8.2 repeats.

Fig. 1
figure 1

A STR length of each allele across MalariaGEN, B STR length across each site, and C summed STR length across each site (line represents median STR length)

Primary analysis

The primary analysis consisted of logistic regression of severe malaria cases status on HMOX1 repeat length, performed at each site individually, then meta-analysed the results using random-effects meta-analysis. Overall, the meta-analysed effect estimate confidence intervals spanned the null (0R 0.96, 95% CI 0.91–1.02), but was consistent with a small negative effect (e.g., increasing repeat length being protective against severe malaria). There was evidence of variation in effect estimates across the sites (p value for heterogeneity = 0.06), and of variation in effects size across populations. For example, increasing STR length was associated with decreased odds of severe malaria in Nigeria (OR 0.54 for a 10 GT(n) repeat increase in repeat length, 95% CI 0.28–1.07, p = 0.07), but higher odds of severe malaria in Papua New Guinea (OR 1.32; 95% CI 0.94–1.85, p = 0.10). Full results are shown in Fig. 2 and Table 2.

Fig. 2
figure 2

Effect of HMOX1 repeat length on severe malaria case status across each site

Table 2 Effect estimates for each country for our primary analysis

For the severe malaria subtypes (CM, SMA, and these combined), there was no strong evidence of a consistent effect across sites with any severe malaria subtype, with all estimates confidence intervals crossing the null (Fig. 3A–F and Additional file 1: Table S2).

Fig. 3
figure 3

Effect estimates of the effect of HMOX1 repeat length in severe malaria subtypes: A CM, vs all other cases, B CM, vs all controls, C SMA vs all other cases, D SMA vs all other controls, E Both vs all other cases, F Both vs all controls

Allelic and genotypic model

Allelic and genotypic models were generated by classifying the repeat length into three alleles (S, M, L), and 6 genotypes (SS, SM, MM, ML, LS, LL). Table 3 describes the number of included participants with each genotype.

Table 3 Genotypes of included participants

There was weak evidence of an association between HMOX1 STR genotype and the outcome “any severe malaria subtype”. The strongest evidence for association was seen for the primary outcome of severe malaria and SM genotype, which reported an OR of 1.15 (95% CI 1.01–1.13) (Additional file 1: Figure S2). All genotypes were associated with an OR of  > 1 (Additional file 1: Table S3), suggesting the reference genotype had the lowest risk of severe malaria, although effect estimates were imprecise and confidence intervals cross the null. In contrast, our allelic model showed no clear association between the number of L, M, or S alleles carried and the development of severe malaria, with wide uncertainty around effect estimates (Additional file 1: Figure S3).

Interaction with clinical variables

Detailed clinical metadata was available for one clinical site (Kenya) and was used to assess any association between the HMOX1 STR and clinical variables. Demographics of the included cohort, split by case status, are available in Table 4.

Table 4 Characteristics of the Kenya cohort [10]

No associations between the HMOX1 STR length and any clinical variable directly were identified (Additional file 1: Table S4). In analyses testing for an interaction between STR length and each clinical variable no evidence for differences in STR association by clinical variable were identified (all p for interaction  > 0.05).

Table 4 summarizes the clinical variables available across each malaria subtype.

All analyses were repeated without the addition of PCA’s and clinical sex as covariates. Results were similar, with no major changes in effect estimates.

Discussion

In the large, multi-centre, MalariaGEN cohort the HMOX1 STR polymorphism was imputed using a recently published reference panel for STRs. STR imputation could not be tested directly on this cohort, but external validation in other cohorts of similar ancestry showed reasonable imputation accuracy [11]. Across the whole cohort, there was no strong evidence that the length of the HMOX1 STR altered the risk of severe malaria, or any severe malaria subtypes, although the data was consistent with a small negative effect at some sites, particularly Kenya, Vietnam, and Nigeria. Additionally there was no association with any clinical variables in the cohort for which clinical metadata was available, nor any interaction with clinical variables and outcomes.

The relevance of the HMOX1 polymorphism in malaria has been recently reviewed [14]. Seven studies have focussed on this question, although all were small and single centre, with a maximum size of 749 patients [15,16,17,18,19,20,21]. In all previous studies, STR length was measured using the STR length fragment polymorphism. Nearly all studies used different definitions of STR length (dividing alleles arbitrarily into “short”, “medium” or “long”). Five studies compared development of severe malaria subtypes within cases, with three studies identifying an association with development of severe malaria in cases (strongest association: OR of 0.16 for severe respiratory distress; 95% CI 0.05–0.46 with carriage of a longer allele in Walther et al. [17]), while two studies did not identify any association with development of any severe malaria subtypes. For comparisons of cases and controls, the data were less clear, with only one study (Mendonca et al.) [16], identifying an association between symptomatic malaria and STR length, although this study was largely (90% of cases) in patients with Plasmodium vivax, not Plasmodium falciparum malaria.

Additionally, a recent large scale phenome wide association study using the HMOX1 STR across a wide variety of clinical and laboratory parameters was performed in UK Biobank, a large, UK wide cohort study. In this study, there were no robust associations of the HMOX1 STR repeat [11].

Given the inconsistency of in-vitro evidence supporting the biological plausibility of effect, and the suggestion of a differential effect at certain sites (p = 0.06 for heterogeneity), it is worth exploring why the effect might differ at differing sites. Firstly, and most simply, the demographics and recruitment of cases at each site were quite different, with a median age of 29 years in Vietnam, to 1.7 years in Tanzania. Given the strong interactions between severe malaria and age, alongside the myriad of other demographic factors, the impact of the HMOX1 STR may only occur in certain clinical settings. Secondly, the genetic background at each site was markedly different. As STR’s also represent recent evolution, and the imputation relates to SNP haplotype, it may be that in certain sites these SNPs associate with severe disease, but do not at other sites. Against that, there has been little previous evidence that SNPs in or near HMOX1 have previously associated with severe malaria [9]. As discussed above, population structure may have also impacted the results, although there was no evidence that estimates varied with or without the use of PCs.

Finally, it may simply represent the play of chance or the microsatellite imparts no functional role in HO-1 expression and/or translation. Given the uncertain previous literature on malaria and the variable functional importance of the repeat in gene expression, together with the weak evidence here, it is possible that there simply is a minimal effect, but that study and population structure influence the results, meaning it is hard to determine conclusively if the effect is present.

Limitations

The major limitation of work undertaken here is a reliance on imputation, rather than by direct genotyping of the HMOX1 STR. However, the reliability of this imputation in four separate datasets has been confirmed in recent work [11]; while other published data supports the reliability of the imputation [8, 11]. Secondly, although the imputed STR length was highly correlated with true STR length, it was much less reliable at calling the exact allele length (~ 45% correct), although ~ 80% were called to within 2 repeats, in the validation data. This is partly due to a large number of potential alleles (64 potential STR lengths), with a smaller number of common alleles, making imputation technically challenging [2]. Although this technical limitation should be recognized, it is important to note that all prior associations with this STR (both in and outside malaria) have been with STR length, with no data suggesting the effect is related to a particular allele rather than the total length of the STR [2,3,4,5, 14, 17, 22,23,24]. Additionally, a genome-wide analysis of 2060 expression short tandem STRs found linear associations to be the most common association between STR’s and gene expression, a finding also identified in other studies of human STRs [25,26,27,28]. Also, other methods of genotyping STR’s, including fragment length polymorphism measurements used in all previous studies on malaria are also associated with some error in true repeat length due to PCR strand slippage.

Secondly, although principle componsents were used to control for relatedness and population structure within each site, it is well established that PC’s do not completely control for population structure [29]. As ethnic groups are likely to differ in both repeat length and incidence of severe malaria, it may be that population structure within each country is not adequately controlled for in all sites, and the effects identified relate to population stratification. Against that, there were no changes in estimates when using models with and without PC’s, suggesting population stratification did not alter our results.

Implications

This study is suggestive that HMOX1 STR variation may not affect the risk of severe malaria, but it remains uncertain whether there may be a suggestive effect, only present under certain environmental conditions. Further research mechanistically linking HMOX1 genetic variation, gene expression of HMOX1 and downstream effects on the host–pathogen interaction would be helpful in understanding the role of heme metabolism in severe malaria. Given the potential to therapeutically modulate HMOX1 expression, this might provide a promising avenue for therapeutics in the future [30].

In this large, multisite, study on the HMOX1 STR polymorphism, including a diverse worldwide population with robust outcome data, results were most consistent with a null effect at most sites, although there was weak indication that STR length might have a small effect in some sites.

Conclusions

The HMOX1 promoter STR was not robustly associated with severe malaria, in a large, multi-country cohort.

Availability of data and materials

This analysis was performed using the MalariaGEN dataset. Access to this can be arranged via the MalariaGEN website and independent data access committee.

References

  1. 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–10.

    Article  CAS  PubMed  Google Scholar 

  2. 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–8.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. 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–104.

    Article  CAS  PubMed  Google Scholar 

  5. Bao W, Song F, Li X, Rong S, Yang W, Wang D, et al. Association between heme oxygenase-1 gene promoter polymorphisms and type 2 diabetes mellitus: a HuGE review and meta-analysis. Am J Epidemiol. 2010;172:631–6.

    Article  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  Google Scholar 

  7. Ramos S, Carlos AR, Sundaram B, Jeney V, Ribeiro A, Gozzelino R, et al. Renal control of disease tolerance to malaria. Proc Natl Acad Sci USA. 2019;116:5681–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Network MGE. Insights into malaria susceptibility using genome-wide data on 17,000 individuals from Africa. Asia and Oceania Nat Commun. 2019;10:5732.

    Article  Google Scholar 

  10. Ndila CM, Uyoga S, Macharia AW, Nyutu G, Peshu N, Ojal J, et al. Human candidate gene polymorphisms and risk of severe malaria in children in Kilifi, Kenya: a case-control association study. Lancet Haematol. 2018;5:e333–45.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Hamilton F, Mitchell R, Ghazal P, Timpson N. Phenotypic associations with the HMOX1 GT(n) promoter repeat in European populations. bioRxiv. 2022. https://doi.org/10.1101/2022.04.28.22274424.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, 1000 Genomes Project Consortium, et al. A global reference for human genetic variation. Nature. 2015;526:68–74.

    Article  PubMed  Google Scholar 

  13. Cavalli-Sforza LL. The human genome diversity project: past, present and future. Nat Rev Genet. 2005;6:333–40.

    Article  CAS  PubMed  Google Scholar 

  14. Hamilton F, Somers J, Mitchell R, Ghazal P, Timpson N. HMOX1 genetic polymorphisms and outcomes in infectious disease: a systematic review. PLoS ONE. 2021. https://doi.org/10.1101/2021.08.21.21262348.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 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–54.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 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–71.

    CAS  PubMed  Google Scholar 

  19. 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–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 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–53.

    PubMed  Google Scholar 

  23. Endler G, Exner M, Schillinger M, Marculescu R, Sunder-Plassmann R, Raith M, et al. A microsatellite polymorphism in the heme oxygenase-1 gene promoter is associated with increased bilirubin and HDL levels but not with coronary artery disease. Thromb Haemost. 2004;91:155–61.

    Article  CAS  PubMed  Google Scholar 

  24. 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.

    Article  CAS  PubMed  Google Scholar 

  25. Gymrek M, Willems T, Guilmatre A, Zeng H, Markus B, Georgiev S, et al. Abundant contribution of short tandem repeats to gene expression variation in humans. Nat Genet. 2016;48:22–9.

    Article  CAS  PubMed  Google Scholar 

  26. Contente A, Dittmer A, Koch MC, Roth J, Dobbelstein M. A polymorphic microsatellite that mediates induction of PIG3 by p53. Nat Genet. 2002;30:315–20.

    Article  PubMed  Google Scholar 

  27. Gebhardt F, Zänker KS, Brandt B. Modulation of epidermal growth factor receptor gene transcription by a polymorphic dinucleotide repeat in intron 1. J Biol Chem. 1999;274:13176–80.

    Article  CAS  PubMed  Google Scholar 

  28. Shimajiri S, Arima N, Tanimoto A, Murata Y, Hamada T, Wang KY, et al. Shortened microsatellite d(CA)21 sequence down-regulates promoter activity of matrix metalloproteinase 9 gene. FEBS Lett. 1999;455:70–4.

    Article  CAS  PubMed  Google Scholar 

  29. Elhaik E. Why most Principal Component Analyses (PCA) in population genetic studies are wrong. BioRxiv. 2021. https://doi.org/10.1101/2021.04.11.439381v2.article-metrics.

    Article  Google Scholar 

  30. Thomas RAB, Czopek A, Bellamy COC, McNally SJ, Kluth DC, Marson LP. Hemin preconditioning upregulates heme oxygenase-1 in deceased donor renal transplant recipients: a randomized, controlled, phase IIB trial. Transplantation. 2016;100:176–83.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Professor Tom Williams and colleagues for access to the site specific data for Kenya.

Funding

FH’s time was funded by the GW4-CAT Wellcome Doctoral Fellowship Scheme (222894/Z/21/Z). UK Biobank was funded by the Wellcome Trust, the Medical Research Council, the NIHR, and a variety of other charities (https://www.ukbiobank.ac.uk/learn-more-about-uk-biobank/about-us/our-funding). FinnGen is a public–private partnership (https://www.finngen.fi/en/access_results) funded by multiple instititions across Finland. The authors want to acknowledge the participants and investigators of the FinnGen study. PG’s time was funded by the Ser Cymru programme, the Welsh Government, and the EU-ERDF. NJT is a Wellcome Trust Investigator (202802/Z/16/Z), is the PI of the Avon Longitudinal Study of Parents and Children (MRC & WT 217065/Z/19/Z), is supported by the University of Bristol NIHR Biomedical Research Centre (BRC-1215-2001), the MRC Integrative Epidemiology Unit (MC_UU_00011/1) and works within the CRUK Integrative Cancer Epidemiology Programme (C18281/A29019). The funder had no role in the design, analysis, or reporting of this study.

Author information

Authors and Affiliations

Authors

Contributions

FH designed the study and performed the analysis. RM, PG, and NJT supervised FH and edited drafts. AC edited drafts and contributed to study design and analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Fergus Hamilton.

Ethics declarations

Competing interests

The authors declare 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.

Influence of ethnicity on case status and severe malaria subtypes. Table S2. Effect estimates for severe malaria subtypes (values represent OR, 95% CI, and p value). Table S3. Effect estimates for the primary outcome using an allelic and genotypic model. Figure S1. Effect estimates from the allelic model for all outcomes. Figure S2. Effect estimates from the genotypic model. Table S2. Differing ethnicities and case status in the Kenya population. Table S3. Estimates for the association between HMOX1 repeat length and outcome in a genotypic and alleleic model. Table S4. Estimates for the association between HMOX1 repeat length and clinical variables.

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

Hamilton, F., Mitchell, R., Cunnington, A. et al. HMOX1 STR polymorphism and malaria: an analysis of a large clinical dataset. Malar J 21, 342 (2022). https://doi.org/10.1186/s12936-022-04352-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-022-04352-x