Interferon-gamma polymorphisms and risk of iron deficiency

Anaemia is a major public health concern especially in Background African children living in malaria-endemic regions. Interferon-gamma (IFN-γ) is elevated during malaria infection and is thought to influence erythropoiesis and iron status. Genetic variants in the IFN-γ gene  ) (IFNG are associated with increased IFN-γ production. We investigated putative functional single nucleotide polymorphisms (SNPs) and haplotypes of IFNG in relation to nutritional iron status and anaemia in Gambian children over a malaria season. We used previously available data from Gambian family trios to Methods: determine informative SNPs and then used the Agena Bioscience MassArray platform to type five SNPs from the  gene in a cohort of IFNG 780 Gambian children. We also measured haemoglobin and biomarkers of iron status and inflammation at the start and end of a malaria season. We identified five  haplotype-tagging SNPs ( -1616 Results: IFNG IFNG [rs2069705], +874 [rs2430561], +2200 [rs1861493],  +3234 IFNG IFNG IFNG [rs2069718] and +5612 [rs2069728]). The +2200C [rs1861493] IFNG IFNG allele was associated with reduced haemoglobin concentrations (adjusted β -0.44 [95% CI -0.75, -0.12]; Bonferroni adjusted P = 0.03) and a trend towards iron deficiency compared to wild-type at the end of the malaria 1 2 1,3 4


Abstract
: Anaemia is a major public health concern especially in Background African children living in malaria-endemic regions. Interferon-gamma (IFN-γ) is elevated during malaria infection and is thought to influence erythropoiesis and iron status. Genetic variants in the IFN-γ gene ) (IFNG are associated with increased IFN-γ production. We investigated putative functional single nucleotide polymorphisms (SNPs) and haplotypes of IFNG in relation to nutritional iron status and anaemia in Gambian children over a malaria season.
We used previously available data from Gambian family trios to Methods: determine informative SNPs and then used the Agena Bioscience MassArray platform to type five SNPs from the gene in a cohort of IFNG 780 Gambian children. We also measured haemoglobin and biomarkers of iron status and inflammation at the start and end of a malaria season.  report report towards iron deficiency compared to wild-type at the end of the malaria season in multivariable models adjusted for potential confounders. A haplotype uniquely identified by +2200C was similarly associated with IFNG reduced haemoglobin levels and trends towards iron deficiency, anaemia and iron deficiency anaemia at the end of the malaria season in models adjusted for age, sex, village, inflammation and malaria parasitaemia.

Introduction
Malaria and iron deficiency are major public health problems for children living in sub-Saharan Africa. The majority (94%) of the 405,000 global deaths due to malaria in 2018 occurred in sub-Saharan Africa, where up to 24% of the population have malaria parasitaemia at any given time 1,2 . In this region, iron deficiency (ID) and anaemia are highly prevalent 3,4 , and may lead to impaired brain development 5,6 , while iron deficiency anaemia (IDA) is a leading cause of years lived with disability in African children 7 . Increasing evidence suggests that malaria may be contributing to ID and IDA 8,9 . Previous studies reported that the prevalence of ID and IDA increased over the malaria season in Gambian and Kenyan children 10,11 , and decreased with the interruption of malaria transmission in the Kenyan highlands 12 .
Multiple factors may contribute towards the development of ID and IDA following a malaria infection. One such factor is interferon-gamma (IFN-γ), which is induced during acute and persistent malaria infection 13 . Among other type 1 responses, IFN-γ is involved in regulating erythropoiesis 14-16 and ironregulatory proteins 17-20 . IFN-γ has also been reported to increase the expression of hepcidin 17 , and divalent metal transporter 1 (DMT1) 19 , while suppressing ferroportin 18,19 , ferritin 20 , and transferrin receptors 19,20 . This regulation of iron proteins may be aimed at starving invading pathogens of iron, a critical nutrient for pathogen growth, but could also play an important role in the pathogenesis of ID and IDA. Indeed, higher IFN-γ levels have been reported in Kenyan children with severe malarial anaemia 21 .
Single nucleotide polymorphisms (SNPs) in the IFN-γ gene (IFNG) on chromosome 12q14 have been associated with increased production of IFN-γ 22,23 , and with susceptibility to severe malaria 24,25 . Despite evidence that malaria induces the production of IFN-γ 13 and that this cytokine influences iron regulation 17-20 , it is not known whether variation in the IFNG gene influences the risk of ID and IDA among children in malaria-endemic areas. We investigated SNPs and haplotypes in the IFNG gene locus in relation to nutritional iron status and anaemia in a cohort of 780 Gambian children prior to and at the end of a malaria season, using an approach based on informative SNPs and Agena Bioscience MassArray platform typing.

Study area
The study was conducted in ten rural villages in the West Kiang region of The Gambia at the start (July 2001) and end (December 2001/January 2002) of a malaria season, as previously described 26 . Malaria incidence is highly seasonal in The Gambia, with the majority of cases occurring between September and December. The study participants were from the Mandinka and Fulani ethnic groups. All of the Fulani children were located in a single village and ethnic group was accounted for in all analyses by adjusting for village.

Study design
We used previously collected data from a cohort of 780 children aged two to six years, recruited at the start of a malaria season as previously described 26 . All children had a clinical examination and a blood sample collected for full blood count, malaria film, and biomarkers of iron status and inflammation at the start and end of the malaria season. Children with pyrexia (temperature ≥37.5°C) had appropriate clinical investigations, clinical treatment and a blood sample taken 2 weeks later after recovery from illness. All children received a 3-day course of mebendazole for possible hookworm infection at recruitment.

SNPs and haplotype construction
Genotypes were determined on whole-genome amplified DNA (primer extension pre-amplification) by the Agena Bioscience MassArray platform (formerly SEQUENOM) using matrixassisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry as previously described 25 . Details of the primer sequences and assays are given in Extended datafiles 1 and 2 28 . The most informative haplotype-tagging SNPs (htSNPs) to type in Gambian subjects were identified by analysing the pattern of linkage disequilibrium (LD) in the IFNG gene loci using previously available data from 32 Gambian family trios 25, 29 . The PHASE program (http://stephenslab.uchicago.edu/ software.html) version 2.1 was used to infer haplotypes from the genotypes of the study population and estimate the frequency of each inferred haplotype 30 . The entropy maximization method was used to identify htSNPs that described >90% of the observed haplotypic diversity in this gene region 29 . The HaploXT program (http://www.sph.umich.edu/csg/abecasis/GOLD/docs/haploxt. html) was used to estimate pairwise LD statistics. Sickle cell (HbS, rs334) and glucose-6-phosphate deficiency (G6PD) deficiency (rs1050828 and rs1050829) were also genotyped using the Agena Bioscience MassArray platform.

Definition of terms
Inflammation was defined as ACT >0.6 g/L. ID was defined as ferritin <12 µg/L or <30 µg/L in the presence of inflammation or <15 µg/L in children ≥5 years 31 , anaemia as Hb<11.0 g/dL (or Hb <11.5 g/dL in children ≥5 years) and IDA as ID plus anaemia 32 .

Statistical analyses
Statistical analyses were conducted using STATA 15.1 (StataCorp. College Station, Texas, USA). Categorical data were expressed as proportions with corresponding percentages. Pearson chi-squared test was used to compare the prevalence of malaria and iron status (ID, IDA and anaemia) at the start and end of the malaria season. Changes in haemoglobin levels and markers of iron status over the malaria season were assessed using the paired t-test. Biological data that were not normally distributed were log-transformed, and geometric means were calculated from original untransformed values.
Log-transformed markers of iron status and risk of ID, IDA and anaemia were analysed using univariable and multivariable linear and logistic regression models, as appropriate. Multivariable regression models were adjusted for age (grouped by year), sex, village (which also acted as a proxy for ethnic group), malaria parasitaemia and ACT. The Bonferroni correction for multiple testing 33 was applied when the five SNPs and six haplotypes were considered individually as independent factors. For multivariable analyses, P values are noted as adj. P for non-Bonferroni corrected analyses and as Bonferroni adj. P for multivariable analyses that are Bonferroni corrected, and for univariable models P values are similarly presented as P or Bonferroni P if Bonferroni corrected. All analyses were considered statistically significant at P<0.05.

Ethics
Individual written informed consent was obtained from children's parents or guardians and the study was approved by The Gambian Government and the Medical Research Council Ethics Review Committee (874/830).

Characteristics of participants
A total of 756 children, including 403 males (53%) and 353 females (47%), were followed up to the end of the malaria season. Most of the children were from the Mandinka ethnic group (n = 681; 90%, compared to Fulani n = 75; 10%). A total of 99/751 (13.2%) children carried sickle cell trait (HbAS) and 136/683 (19.9%) children G6PD deficiency. The prevalence of ID and IDA increased over the malaria season (from 20.6% to 31.6% and from 11.9% to 21.7%, respectively), as previously reported 11,26 . Individual markers of iron status also reflected an increase in ID over the malaria season. Table 1 summarises the characteristics of the study population and their iron status at the beginning and end of the malaria season.

Associations with iron and anaemia
We found that the IFNG SNPs were not associated with ID, IDA or anaemia at the start of the malaria season in multivariable logistic regression analyses adjusted for age, sex, village, ACT and malaria parasitaemia following Bonferroni adjustment (Extended datafile 4 28 ).
The IFNG SNPs were similarly not significantly associated with ID, IDA or anaemia at the end of the malaria season after Bonferroni correction. The IFNG+2200C allele (rs1861493) was associated with trends towards increased risk of anaemia (  (Table 3).

Discussion
In this study we observed an increase in the prevalence of ID, IDA and anaemia across the malaria season in Gambian children.
Dietary iron insufficiency may be an important cause since the malaria season also coincides with the 'hungry season' in The Gambia when there is a scarcity of staple foods. We hypothesized that IFN-γ, a pro-inflammatory cytokine induced during malaria infection 13 , might play a role in influencing the risk of ID and anaemia in children exposed to malaria. In addition to directly reducing erythrocyte half-life 15 , evidence suggests that IFN-γ induces hepcidin and inhibits ferroportin, hence reducing iron absorption and promoting sequestration of iron in macrophages 17,19 . Consequently, high levels of IFN-γ induced during malaria infections may concomitantly lead to ID and anaemia.
We found that the IFNG+2200C (rs1861493) allele, located at intron 3 of the IFNG gene, was associated with reduced haemoglobin levels and a trend towards ID, IDA and anaemia at the end of a malaria season in multivariable analyses adjusting for potential confounders. We then constructed haplotypes to increase the probability of capturing functional mutations which might reside within a given haplotype. Haplotype 6 (uniquely identified by the IFNG+2200C allele), was associated with reduced haemoglobin levels and TSAT and increased ZnPP levels in keeping with iron deficiency compared to the wild-type haplotype. Haplotype 6 was similarly associated with trends towards increased risk of ID, IDA and anaemia at the end of the malaria season.  So how might the IFNG+2200C genotype and a haplotype uniquely defined by this genotype potentially lead to reduced haemoglobin levels and ID at the end of the malaria season? A possible explanation may be through increasing IFNG gene expression and IFN-γ levels. The IFNG+2200C allele was associated with increased IFN-γ levels in Kawasaki disease patients 23 , although another study in patients with ankylosing spondylitis found no difference in IFN-γ levels by IFNG+2200 genotype 34 . Elevated IFN-γ levels promote dyserythropoiesis, anaemia and iron dysregulation. IFN-γ inhibits proliferation of erythroid progenitor cells by disrupting lineage differentiation, blocking renal production of erythropoietin, inhibiting renal iron reabsorption, and reducing red blood cell half-life 14-16, 35 . As a type 1 immune response, IFN-γ also induces defensive transcriptional programs within enterocytes resulting in reduced dietary iron absorption 36 . Additionally, IFN-γ promotes iron sequestration to macrophages either directly or through its influence on hepcidin, ferroportin, and DMT1 17-20 . These responses reduce circulating transferrin-bound iron, which is required by Plasmodium parasites for metabolism and proliferation 37 .
In keeping with this, we observed decreased ferritin levels and TSAT and increased ZnPP levels in children carrying the IFNG+2200C haplotype. Lower haemoglobin levels may also translate into reduced amino acid availability for Plasmodium parasites 38 and hence protection against blood-stage parasitaemia. Indeed, we observed that the IFNG+2200C SNP was associated with protection against malaria parasitaemia at the end of the malaria season (adj. OR 0.40; Bonferroni adj. P=0.03).
The influence of the IFNG+2200C allele on haemoglobin and iron status was only observed at the end of the malaria season. We hypothesized that the effects of this SNP may be most marked when expression of IFNG is upregulated, such as during malaria infections 13,24 . This also highlights the influence of gene-environment interactions in promoting disease, in this case ID, IDA and anaemia. It is unlikely that the decreased haemoglobin levels observed in individuals carrying the IFNG+2200C allele was due to increased malaria since these children had reduced prevalence of malaria parasitaemia at the end of the malaria season. It is possible that higher IFN-γ levels, putatively produced by IFNG+2200C carriers, may induce a protective proinflammatory response against malaria 13,39 , but at the expense of iron homeostasis. The IFNG+874TT (rs2430561) genotype, located at the first intron coinciding with the NFkB binding region, has also been associated with higher production of IFN-γ 22, 34,40 . However, studies have found no association between IFNG+874TT and malaria 25 or aplastic anaemia 41 , and in our study we observed a trend towards a decreased risk of ID in these individuals. Further investigations are required on a cellular level to explore putative functional effects of IFNG genotypes on IFN-γ levels and iron status.
To our knowledge, this is the first study examining the role of IFNG gene polymorphisms in relation to iron status. We found that the IFNG+2200C (rs1861493) allele, and a haplotype defined by this allele, were associated with reduced haemoglobin levels and a trend towards ID at the end of the malaria season, a finding that may be due to increased IFN-γ levels 23 . However, our study had a number of important limitations and our findings should be viewed with considerable caution. The study was conducted in a single site and had relatively small numbers (n = 756). Additionally, many of our findings were of marginal significance with wide confidence intervals and lost statistical significance after correction for multiple testing with Bonferroni adjustment. It is also unclear if our findings have clinical relevance at an individual level. The IFNG+2200C SNP may also be in linkage disequilibrium (LD) with another genetic variant within the haplotype that might influence IFN-γ levels and / or measures of iron status in Gambian populations. Finally, we did not measure IFN-γ levels to determine if they differed between genotypes at the end of the malaria season. Thus, our findings need to be examined in larger population-based studies, in other malaria-exposed populations, and functional assays are needed to identify whether genetic variation in the IFNG gene influences iron status. However, our study supports the hypothesis that preventing and treating malaria infection may improve haemoglobin levels and iron status in African children 8 . • IFNG_final_data_v2 (dataset containing demographic information and results of laboratory assays for participants included in the study).
• IFNG_analysis_KM (contains the codes used for data analysis).
• KMokaya_IFNG_Codebook (contains variable description and labels). © 2020 Silvestri L. This is an open access peer review report distributed under the terms of the Creative Commons , which permits unrestricted use, distribution, and reproduction in any medium, provided the original Attribution License work is properly cited.

Laura Silvestri
Regulation of Iron Metabolism Unit, Division of Genetics and Cell Biology, San Raffaele Scientific Institute, Milan, Italy In this paper Abuga et al investigate the role of IFNG SNPs in influencing iron deficiency (ID) and iron deficiency anemia (IDA) in Gambian children over a malaria season. They identify a SNP, IFNG+2200C, located in intron 3, associated with a reduced risk of malaria parasitemia, reduced Hb, and a trend towards ID and IDA. The authors conclude that although more studies are needed to assess the role of IFNG SNPs in ID and IDA, this represents the first study that investigates the association of IFNG genetic variants with iron status.
The study is well conducted and of interest. I have only minor comments: Serum hepcidin levels are influenced by circadian rhythm and serum/iron stores. Please indicate at what time of the day blood was drawn. To "normalize" serum hepcidin to body iron concentration, the hepcidin/ferritin or hepcidin/TSAT ratio should be shown in Table 1 and In the M&M section, the authors claim that some children are carriers of the HbAS, and G6PD deficiency. Is it possible to present the hematological and iron data (in Table 1 and Table 3) related to these children?