Original Research ARTICLE
Iron Deficiency Anemia at Time of Vaccination Predicts Decreased Vaccine Response and Iron Supplementation at Time of Vaccination Increases Humoral Vaccine Response: A Birth Cohort Study and a Randomized Trial Follow-Up Study in Kenyan Infants
- 1Department of Health Sciences and Technology, Institute of Food, Nutrition and Health, Laboratory of Human Nutrition, ETH Zürich, Zurich, Switzerland
- 2Department of Environment and Health Sciences, Technical University Mombasa, Mombasa, Kenya
- 3MRC Human Immunology Unit, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, Oxford, United Kingdom
- 4Pediatrics Department, Msambweni County Referral Hospital, Msambweni, Kenya
- 5Centre for Infectious Disease Control, National Institute for Public Health and Environment (RIVM), Bilthoven, Netherlands
- 6Center for Global Health and Diseases, Case Western Reserve University School of Medicine, Cleveland, OH, United States
- 7Division of Infectious Diseases, Department of Pediatrics, Lucille Packard Children's Hospital at Stanford School of Medicine, Stanford, CA, United States
- 8Pediatric Epidemiology, Department of Pediatrics, Medical Faculty, Leipzig University, Leipzig, Germany
- 9Department of Medical Epidemiology, College of Health Sciences, Jomo Kenyatta University of Agriculture and Technology, Nairobi, Kenya
- 10Haematology Theme, NIHR Oxford Biomedical Research Centre, Oxford, United Kingdom
Background: Iron deficiency may impair adaptive immunity and is common among African infants at time of vaccination. Whether iron deficiency impairs vaccine response and whether iron supplementation improves humoral vaccine response is uncertain.
Methods: We performed two studies in southern coastal Kenya. In a birth cohort study, we followed infants to age 18 mo and assessed whether anemia or iron deficiency at time of vaccination predicted vaccine response to three-valent oral polio, diphtheria-tetanus-whole cell pertussis-Haemophilus influenzae type b vaccine, ten-valent pneumococcal-conjugate vaccine and measles vaccine. Primary outcomes were anti-vaccine-IgG and seroconversion at age 24 wk and 18 mo. In a randomized trial cohort follow-up, children received a micronutrient powder (MNP) with 5 mg iron daily or a MNP without iron for 4 mo starting at age 7.5 mo and received measles vaccine at 9 and 18 mo; primary outcomes were anti-measles IgG, seroconversion and avidity at age 11.5 mo and 4.5 y.
Findings: In the birth cohort study, 573 infants were enrolled and 303 completed the study. Controlling for sex, birthweight, anthropometric indices and maternal antibodies, hemoglobin at time of vaccination was the strongest positive predictor of: (A) anti-diphtheria and anti-pertussis-IgG at 24 wk (p = 0.0071, p = 0.0339) and 18 mo (p = 0.0182, p = 0.0360); (B) anti-pertussis filamentous hemagglutinin-IgG at 24 wk (p = 0.0423); and (C) anti-pneumococcus 19 IgG at 18 mo (p = 0.0129). Anemia and serum transferrin receptor at time of vaccination were the strongest predictors of seroconversion against diphtheria (p = 0.0484, p = 0.0439) and pneumococcus 19 at 18 mo (p = 0.0199, p = 0.0327). In the randomized trial, 155 infants were recruited, 127 and 88 were assessed at age 11.5 mo and 4.5 y. Compared to infants that did not receive iron, those who received iron at time of vaccination had higher anti-measles-IgG (p = 0.0415), seroconversion (p = 0.0531) and IgG avidity (p = 0.0425) at 11.5 mo.
Interpretation: In Kenyan infants, anemia and iron deficiency at time of vaccination predict decreased response to diphtheria, pertussis and pneumococcal vaccines. Primary response to measles vaccine may be increased by iron supplementation at time of vaccination. These findings argue that correction of iron deficiency during early infancy may improve vaccine response.
Increasing immunization efficacy and reducing anemia are key global pediatric health goals. Immunization programs have achieved high coverage, yet one in five children worldwide are not fully protected, contributing to 1.5 million child deaths annually from vaccine-preventable diseases (1). Vaccines often underperform in low- and middle-income countries (LMIC) (2, 3); for example, effectiveness of measles vaccine is generally <75% in Sub-Saharan Africa (4). Why vaccines underperform in LMIC remains uncertain (2, 3) but nutritional deficiencies may play a role (3, 5).
Over 40% of children <5 y worldwide are anemic, many due to iron deficiency (ID) (6). Anemia is particularly common in infants age <1 y in Sub-Saharan Africa: >70% of infants may be anemic at time of routine vaccination in the first year (7, 8). ID not only causes anemia, but also may impair adaptive immunity and thereby vaccine efficacy, although the data are equivocal (5, 9, 10). ID in mice attenuates T-cell dependent and independent antigen-specific antibody responses, and impairs cyclin E1 induction and S-phase entry during B-cell proliferation (10). Data from observational studies in humans are contradictory, in some, ID was associated with defective immune function, particularly T-cell immunity (5, 9). Clinically, a homozygous mutation in transferrin receptor (TfR)1 (the protein which transports iron into lymphocytes) causes severe combined immunodeficiency with low levels of circulating immunoglobulins (Ig) and decreased T and B cell proliferation in vitro (11). Data linking anemia, ID and vaccine response in LMIC are scarce; most data are from older studies that used outmoded methods to assess vaccine response (5).
In southern coastal Kenya, 70–75% of infants are anemic at the time they receive their routine vaccinations (8, 12). Therefore, we performed two studies to determine if anemia and/or ID during infancy affects vaccine response. In a birth cohort study, we followed infants to age 18 mo and assessed whether anemia or ID at time of vaccination predicted response to three-valent oral polio (OPV), diphtheria, tetanus, whole-cell pertussis, Haemophilus influenzae type b (Hib), ten-valent pneumococcal polysaccharide (PnPs) and measles vaccines (MV). Primary outcomes were anti-vaccine serum IgG and seroconversion at age 24 wk (primary response) and 18 mo (secondary response). We hypothesized that anemia and/or ID at time of vaccination would predict lower serum IgG and reduced seroconversion at these ages. In a randomized trial cohort follow-up, children received a micronutrient powder (MNP) with 5 mg iron daily or a MNP without iron for 4 mo starting at age 7.5 mo and received MV at age 9 and 18 mo. Primary outcomes were anti-measles serum IgG, seroconversion and IgG avidity assessed at 11.5 mo (primary response) and 4.5 y (secondary response). We hypothesized that iron supplementation at time of first MV would result in higher anti-measles serum IgG, seroconversion and IgG avidity at both time points.
Participants and Methods
These studies were done in Kwale County, in southern coastal Kenya. We recruited mothers and infants from the maternal care and infant growth monitoring and vaccination clinics in Msambweni County Referral Hospital and surrounding health centers. The two studies enrolled independent groups of infants. For the birth cohort study, we obtained ethical approval from the institutional review boards at Kenyatta National Hospital/University of Nairobi (KNH/UoN), Case Western Reserve University and the Stanford University School of Medicine. For the randomized trial, we obtained ethical approval from the institutional review boards of the KNH/UoN and the Zurich Cantonal Ethical Commission; approval for the follow-up study was given by KNH/UoN. The randomized trial was registered and updated on clinical trials (NCT02118402). Caregivers gave informed consent by either a written signature or a fingerprint.
The birth cohort study (study design in Figure 1) was nested within a larger longitudinal maternal-child cohort study on links between maternal parasitic infections and infant vaccine response. In the original maternal-child cohort study [details in (13)], healthy pregnant women ≥15 y of age in the 2nd or 3rd trimester of pregnancy were invited to participate in the study when they presented for antenatal care. If they had an uncomplicated delivery of a healthy infant (at ≥34 wk with a birthweight ≥1500 g and an APGAR score ≥5), then the infant was enrolled into this study. Additional inclusion criteria for this cohort study were: (1) infant remained in the study and provided samples for anti-vaccine serum IgG concentrations at age 18 mo; and (2) one or more hemoglobin (Hb) measurement(s) were done at age 10, 14, or 24 wk.
Figure 1. Design of the birth cohort study. Hb, hemoglobin; OPV, oral polio vaccine; DTP, diphtheria-tetanus-pertussis; HepB, hepatitis B; Hib, Haemophilus influenzae type b; PnPs, pneumococcal polysaccharide.
At delivery, we collected a whole blood sample from the umbilical vein. We obtained whole blood samples from the infants at age 10 and 14 wk by heel prick, and at age 24 wk, 12 and 18 mo by venipuncture, and trained staff measured length and weight. Under supervision of the study team, all infants received standardized immunizations following Kenyan Ministry of Health guidelines, as follows: (A) three-valent OPV with serotypes 1, 2, 3 at birth; (B) at age 6, 10, and 14 wk, the infants were given OPV, diphtheria-tetanus-whole cell pertussis-hepatitis B-Hib vaccine, ten-valent PnPs vaccine (Synflorix) with serotypes 1, 4, 5, 6B, 7F, 9V, 14, 18C, 19F, and 23F; and (C) at age 9 mo, they received the monovalent MV. At age 18 mo, we collected a blood sample before the second MV was given. All infants in the study received vitamin A supplementation every 6 months starting at age 24 wk as per Kenyan guidelines. At each visit, trained staff asked the caregiver if the infant had received malaria treatment since the last visit.
Figure 2 shows the study design for the randomized trial cohort follow-up study. In a double-blind, controlled trial [details in (12)], we randomized 155 Kenyan infants (mean age 7.5 mo) to receive daily for 4 mo: (A) a MNP without iron (control); (2) the identical MNP but containing 5 mg of iron as ferrous fumarate and sodium iron ethylenediaminetetraacetic acid (NaFeEDTA) (iron only group); or (3) the identical MNP as the iron only group but with 7.5 g galacto-oligosaccharides (GOS) (Vivinal GOS, Friesland Campina, Wageningen, The Netherlands) (iron and GOS group). 145 infants completed the intervention trial. Out of these 145 infants, all except six received a monovalent MV at age 9 mo (during the intervention). About 3.5 years later, we contacted all families with children who were in the original trial, with the exception of the six children who had received the MV before they entered the intervention trial, and 12 children for whom we did not have adequate plasma to measure the MV response; thus, 127 children were included in the present trial. From these 127 children, 88 families provided consent for the follow-up measurements; of these, 93% of children had received a second dose of MV at age 18 mo as recorded in their vaccination booklets. We measured child weight and length/height and collected a venous blood sample at age 7.5 mo, before they received the first MV, at endpoint of the intervention (at age 11.5 mo), and at follow up (at age 4.5 y).
We measured Hb on the day of collection in both studies using a Hemocue spectrophotometer (Hemocue, Ångelholm, Sweden), with the exception of the 24 wk, 12 and 18 mo visits in the birth cohort study, when we measured Hb using an Automated Hematology Analyzer (Nihon Kohden model: 82 MEK 7222). Serum was separated by centrifugation and frozen on the day of collection until analysis. We measured serum TfR, serum ferritin (SF), high-sensitive C-reactive protein (CRP) and α-1 glycoprotein (AGP) by using a combined sandwich enzyme-linked immunosorbent assay (14). Serum IgG concentrations against PnPs (ten serotypes: 1, 4, 5, 6B, 7F, 9V, 14, 18C, 19F, 23F), purified polyribosylribitol phosphate (PRP) Hib, and diphtheria CRM197 toxoid antigens were measured at Stanford University, Palo Alto, CA, USA using a fluorescent multiplexed bead-based immunoassay (Luminex, Austin TX) (13). The twelve vaccine antigens were coupled to carboxylated microspheres, incubated with patient serum and quantified on a BioPlex MAGPIX multiplex reader (BioRad, Hercules, CA) (13). Serum IgG against tetanus toxoid, pertussis (pertussis toxin, filamentous hemagglutinin (FHA), and pertactin) measles and polio, and avidity of measles IgG, were measured at the National Institute for Public Health and the Environment, Bilthoven, Netherlands. Serum IgG against tetanus, pertussis and measles were simultaneously measured using a multiplex immunoassay (15, 16). Serum IgG against polio was measured using a multiplex poliovirus binding inhibition assay (17). Measles IgG avidity was determined by modifying a bead-based serum IgG immunoassay against measles. After incubation of measles virus-conjugated fluorescent beads with serum and before labeling, beads were washed with phosphate-buffered saline and subsequently left untreated or treated for 10 min with 1.5 M ammonium thiocyanate [NH4SCN], as described previously (18). The avidity index was defined as the percentage of antibodies that remain bound to the beads after treatment with thiocyanate, calculated as: [amount of IgG after thiocyanate treatment]/[amount of IgG without treatment] × 100 (18).
We conducted statistical analyses with SPSS (IBM SPSS statistics, Version 24). Data were checked for normality by Shapiro-Wilk tests and by visual inspection of histogram plots. Non-normally distributed data were logarithmically (log10) transformed for statistical analyses. Descriptive data were expressed as mean ± SD (for normally distributed data) or median (IQR) for non-normally distributed data. We defined ID at 24 wk as either SF <12 μg/L or 2) TfR >8.3 μg/ml, and at 12 and 18 mo as SF <12 μg/L (14) because we did not measure TfR at these later ages. We did not define ID at earlier time points because consensus cut-offs are unavailable for ages <24 wk. We defined anemia as Hb <10.5 g/dl at age 10 and 14 wk, based on proposed cut-off values for Hb at age 16 wk in iron-replete infants (19) and as Hb <11.0 g/dl at age 24 wk, 12 and 18 mo (20). We defined mild, moderate and severe anemia at 10 and 14 wk as an Hb of 9.5–10.5 g/dl, 7.0–9.4 g/dl and <7.0 g/dl, respectively; and at 24 wk, 12 and 18 mo as an Hb of 10.0–11.0 g/dl, 7.0–9.9 g/dl and <7.0 g/dl, respectively (20). We defined ID anemia as anemia with ID based on the above cut-offs. At all ages, we defined CRP and AGP >5 mg/L and >1 g/l respectively as elevated, indicating inflammation (14). Vitamin A deficiency was defined as RBP <0.7 μmol/L at 24 wk (14). We defined immunological response as protective when participants had anti-diphtheria serum IgG concentrations ≥0.1 IU/ml (21), anti-tetanus IgG >0.5 IU/ml (22), anti-Hib IgG >1.0 μg/ml (22), anti-PnPs IgG ≥0.35 μg/ml (23), anti-measles IgG ≥0.12 IU/ml (24) and anti-OPV IgG ≥0.6, ≥2.1 and ≥0.7 IU/ml for S1, 2, and 3 (25). There is no protective serum concentration established for anti-pertussis IgG (22). We converted infant anthropometric measurements to age-adjusted z-scores in Anthro software; children with a length/height-for-age z-score < -2 were categorized as stunted; those with a weight-for-age z-score < -2 as underweight and those with a weight-for-height z-score < -2 as wasted. For the birth cohort study, we defined primary and secondary vaccine response as anti-vaccine serum IgG concentrations and seroconversion at 10 wk and 14.5 mo after vaccination, respectively. For the randomized trial, we defined primary and secondary vaccine response as anti-vaccine serum IgG concentrations and seroconversion at 4 wk and 3.5 y after vaccination, respectively.
In the birth cohort study, we used linear regression analyses on anti-vaccine serum IgG concentrations at age 24 wk (primary vaccine response) and at 18 mo (secondary response) to assess the effect of anemia and ID at the time of vaccination on vaccine response. For anthropometrics and Hb, we used the mean value at age 10, 14, and 24 wk to represent that variable at time of vaccination. In addition to mean Hb, we sequentially added known potential predictors of vaccine response: sex, birthweight, mean anthropometrics, and maternal serum IgG concentrations (in cord blood and at age 10 wk) to the models. We performed logistic regression on seroconversion against diphtheria, tetanus, Hib, PnPs and OPV at age 24 wk and 18 mo, and MV at 12 and 18 mo, including the above covariates except Hb, and added moderate-severe anemia (yes/no) as a discrete variable to the models. We repeated the above analyses including TfR (as an indicator of iron-deficient erythropoiesis), RBP and AGP at age 24 wk. In all regression models we did bootstrapping with a resampling size of 1000.
In the randomized trial, we used a random intercept linear mixed effect model (LMM) analysis with Bonferroni corrected multiple comparisons to assess the effect of group (iron or control) at the time of MV. We performed LMM with anti-measles serum IgG concentrations as the dependent variable in two models to assess primary and secondary vaccine response. In the first model, we included anti-measles serum IgG measured at age 7.5 and 11.5 mo; in the second model, we included anti-measles serum IgG measured at age 7.5 mo and 4.5 y. We sequentially added covariates to the models that are potential predictors of vaccine response: age, weight-for-height z-score (WHZ), IDA and RBP. By adding age to the models, we corrected for differences between subjects in time elapsed from baseline to vaccination and from vaccination to endpoint of the intervention. Controlling for the same covariates, we performed logistic regression on seroconversion and antibody avidity at age 11.5 mo and 4.5 y. To create a binary variable for high and low avidity, we divided the avidity distribution in half, and described the upper half as high avidity. For both, we used the covariates measured at 11.5 mo (reflecting intervention and time of vaccination).
Birth Cohort Study
The birth cohort study began recruiting in August 2013 and follow up was to May 2017. At birth, we enrolled 573 newborns into the study and collected a cord blood sample. At age 18 mo, 361 infants remained in the study and provided samples. 58 of these infants did not have one or more Hb measurement(s) at age 10, 14, or 24 wk. Thus, 303 infants (129 male, 174 female) were included in the data analysis (Figure 1). Only ten infants reportedly received malaria treatment at age 0–24 wk. Mean ±SD birthweight was 3,003 ± 457,g and 34 infants were low birthweight (LBW) (<2,500 g); based on the inclusion criteria, none of the infants were preterm. All infants received all vaccines at 6, 10 and 14 wk; 99% of the infants received the MV at 9 mo.
Nutritional Status, Anemia, and Iron Deficiency
Table 1 shows anthropometrics, iron and inflammation status from age 10 wk to 18 mo. At 24 wk, the prevalence of stunting and underweight were 29.0 and 10.0%, suggesting that overall, the infants were mild-to-moderately malnourished, but the prevalence of wasting was only 3.5%. At age 10 and 14 wk, just over half of infants were anemic; however, median SF was high and median TfR was in the normal range. At age 24 wk and 12 mo, median Hb was only 9.0 g/dl (Figure 3A) and >90% of infants were anemic: nearly all the anemia was moderate or severe; low median SF and high median TfR suggested iron depletion in the majority of infants (Table 1). At 24 wk, the prevalence of vitamin A deficiency was 47%, and 37% of infants showed signs of inflammation, likely due to common infections.
Table 1. Age, anthropometrics, biomarkers of iron, vitamin A and inflammation, and prevalence of iron deficiency, anemia, iron deficiency anemia, vitamin A deficiency and inflammation in the birth cohort study of Kenyan infants.
Figure 3. Hemoglobin and anti-vaccine serum IgG concentrations from birth until 18 months in the birth cohort of Kenyan infants. (A) Hemoglobin, (B) anti-diphtheria IgG, (C) anti-pertussis toxin IgG, (D) anti- filamentous hemagglutinin (FHA) IgG, (E) anti-pneumococcus serotype 19 (PnPs 19) IgG, and (F) anti-measles IgG.
Anti-vaccine Serum IgG Concentrations and Seroconversion
Table 2 shows anti-vaccine serum IgG concentrations and seroconversion for the measured vaccines. At age 24 wk, the primary vaccine response against diphtheria and tetanus was high, with only 1 and 14% of infants seronegative, while response to Hib was lower, with 45% seronegative. The primary response to the PnPs serotypes varied, the highest response was to PnPs 7 and PnPs 19, with only 3 and 6% seronegative at 24 wk, while the lowest response was to PnPs 1 and PnPs 23, with 26 and 24% of infants seronegative at 24 wk. The highest anti-vaccine serum IgG concentrations against many antigens in the diphtheria-tetanus-whole cell pertussis- Hib vaccine were present at 24 wk, and serum IgG waned at 12 and 18 mo (Figures 3B–E). At 18 mo, the secondary vaccine response varied: 16, 32, 62, and 68% were seronegative against OPV3, diphtheria, tetanus and Hib. Secondary response to the PnPs serotypes varied widely: the lowest response was against PnPs 5 (80% seronegative), while the highest response was against PnPs 6 (24% seronegative). MV response was high at 12 and 18 mo, before the second vaccination, with only 7 and 4% seronegative (Table 2 and Figure 3F). Notably, in cord blood, serum IgG against tetanus and measles were high, suggesting high transfer of maternal antibodies.
Table 2. Anti-vaccine serum IgG concentrations and seroconversion for diphtheria, tetanus, pertussis, Haemophilus influenza b, ten pneumococcus serotypes, measles and three oral polio serotypes in the birth cohort of Kenyan infants.
Predictors of Anti-diphtheria Serum IgG Concentrations and Seroconversion
Table 3 shows the associations of Hb and moderate to severe anemia, and vaccine response. For diphtheria, controlling for sex, birthweight, anthropometrics and maternal antibodies, Hb at time of vaccination was the strongest positive predictor of anti-vaccine serum IgG at age 24 wk (p = 0.0071) and 18 mo (p = 0.0182). Moderate or severe anemia was the strongest risk factor for seronegativity at age 18 mo: anemic infants had 2.1 times the risk of being unprotected against diphtheria at 18 mo (p = 0.0484) (Table 3). Table 4 shows the associations of TfR (as an indicator of iron-deficient erythropoiesis) and vaccine response. Controlling for sex, birthweight, anthropometrics, inflammation, vitamin A status and maternal antibodies, TfR at age 24 wk was the strongest risk factor for seronegativity against diphtheria at 18 mo (p = 0.0334).
Table 3. Associations between hemoglobin (Hb) and moderate to severe anemia (Hb <9.5 g/dl), and anti-diphtheria, anti-pertussis, anti-pneumococcus serotype 19, and anti-measles serum IgG concentrations and/or seroconversion, in the birth cohort of Kenyan infants.
Table 4. Associations between serum transferrin receptor and anti-diphtheria, anti-pneumococcus serotype 19 and anti-polio serotype 2 serum IgG concentrations and/or seroconversion, in the birth cohort of Kenyan infants.
Predictors of Anti-pertussis Serum IgG Concentrations and Seroconversion
For pertussis, we assessed response by measuring antibodies against three antigens: pertussis toxin, FHA and pertactin. Controlling for sex, birthweight, anthropometrics and maternal antibodies, Hb at time of vaccination was the strongest positive predictor of anti-pertussis toxin IgG (p = 0.0339) and anti-FHA IgG (0.0423) at 24 wk and the strongest positive predictor of anti-pertussis toxin IgG (p = 0.0360) at 18 mo (Table 3). Hb, anemia and TfR were not significant predictors of anti-pertactin IgG, or of response to the Hib and tetanus vaccines.
Predictors of Anti-pneumococcus Serum IgG Concentrations and Seroconversion
For PnPs, controlling for sex, birthweight, anthropometrics, and maternal antibodies, Hb at time of vaccination was the strongest positive predictor of anti-PnPs 19 serum IgG at 18 mo (p = 0.0129) (Table 3). Moderate or severe anemia was the strongest risk factor for seronegativity at 18 mo: anemic infants had 2.2 times the risk of being unprotected against PnPs 19 at 18 mo (p = 0.0199). Controlling for sex, birthweight, anthropometrics, inflammation, vitamin A status and maternal antibodies, TfR at 24 wk was the strongest negative predictor of anti-PnPs 19 IgG at 18 mo (p = 0.0344) (Table 4). Controlling for the same factors, TfR was the strongest risk factor for seronegativity against PnPs 19 at 18 mo (p = 0.0327). Hb, anemia and TfR were not significant predictors of response to the other PnPs serotypes.
Predictors of Anti-polio Serum IgG Concentrations and Seroconversion
Controlling for sex, birthweight, anthropometrics, inflammation, vitamin A status and maternal antibodies, TfR at 24 wk was the strongest negative predictor of anti-polio S2 IgG at 24 wk (p = 0.0719) (Table 4). Hb, anemia and TfR were not significant predictors of response or seroconversion to polio S1 and S3.
Predictors of Anti-measles Serum IgG Concentrations and Seroconversion
For measles, moderate or severe anemia was the strongest risk factor for seronegativity at 12 mo: anemic infants had 5.1 times the risk of being unprotected against measles at 12 mo (p = 0.0783) (Table 3). Hb and anemia were not significant predictors of response to MV at age 18 mo.
Randomized Trial Cohort Follow-Up Study
The randomized trial began recruiting in October 2014 and follow up was to March 2019. Comparing the two groups that received the 5 mg daily dose of iron (the iron only and the iron and GOS groups), there were no significant differences in any of the measured parameters at age 11.5 mo or at 4.5 y (data not shown). Notably, there were no significant differences in anti-measles serum IgG, seroconversion or avidity comparing the iron only and iron and GOS groups at any time point. Therefore, we pooled the iron only and the iron and GOS groups into a single iron group that had received 5 mg of iron daily for 4 mo, and compared this new iron group to the no-iron group as control (Figure 2). 127 infants were assessed at age 11.5 mo (85 infants in the iron group and 42 infants in the control group) and 88 were assessed at age 4.5 y. Median (IQR) birthweight was higher in the control group than in the iron group (3,300 (3,000–3,500)g vs. 3,000 (2,600–3,500)g) (p = 0.0299). Nine (10.5%) infants in the iron group and four (10.0%) infants in the control group were LBW (<2,500 g). Table 5 shows anthropometrics, Hb, iron and inflammatory parameters and anti-measles serum IgG concentrations at age 7.5 mo, 11.5 mo and 4.5 y. There were no significant group differences in any of the parameters at 7.5 mo. During the intervention, there was a marked improvement in iron status in the iron group, but not in the control group: there were significant group effects for Hb (p = 0.0001) and plasma ferritin (PF, p = 0.0198), with higher Hb and PF in the iron group compared to the control group and there were significant group by time interactions for Hb (p = 0.0001), PF (p = 0.0001) and TfR (p < 0.0001).
Table 5. In a randomized trial follow-up study, age, anthropometrics, biomarkers of iron, vitamin A and inflammation, and anti-measles serum IgG concentrations in Kenyan infants receiving 5 mg iron daily or no-iron control for 4 months, and given measles vaccine at age 9 mo and 18 mo, measured at baseline (age 7.5 mo), end of intervention (age 11.5 mo) and follow-up (age 4.5 y).
Anti-measles Serum IgG Concentrations
For primary vaccine response (age 11.5 mo), there was a significant group effect (p = 0.0415) with higher anti-MV serum IgG in the iron group and a significant time effect (7.5 vs. 11.5 mo; p < 0.0001), but no group*time interaction (p = 0.2952) (Table 5 and Figure 4A). For secondary vaccine response, there was a significant time effect (11.5 mo vs. 4.5 y; p < 0.0001), but no group effect (p = 0.6800) and no group*time interaction (p = 0.6445).
Figure 4. In the randomized controlled trial follow-up study, anti-measles serum IgG concentrations, seroconversion, and IgG avidity in the control and iron groups. (A) Anti-measles serum IgG concentrations at age 7.5 and 11.5 months (at baseline and end of intervention), (B) seroconversion at age 11.5 months (at end of intervention), and (C) IgG avidity at age 11.5 months (at end of intervention).
Seroconversion and Anti-measles IgG Avidity
Before vaccination (at age 7.5 mo), 95.3 vs. 97.6% of infants were seronegative against MV in the iron and control group, respectively. After vaccination, at age 11.5 mo, 15.3 vs. 23.8% of infants were seronegative against MV in the iron and control group (Figure 4B). At age 4.5 y, 9.1 vs. 12.1% of infants were seronegative against MV in the iron and control group. Controlling for WHZ, IDA, RBP and age, being in the control group was the strongest risk factor for seronegativity at age 11.5 mo: infants in the control group had 3.1 times greater risk of being unprotected against measles at this age (p = 0.0531) (Table 6). At age 11.5 mo and 4.5 y, in the iron and control groups, mean (±SD) avidity was 29.1 ± 14.3 vs. 26.8 ± 12.5% (Figure 4C) and 69.1 ± 14.2 vs. 65.6 ± 14.9%, respectively. Controlling for the same parameters as above, being in the control group was the strongest risk factor for low avidity at age 11.5 mo: infants in the control group had 3.2 times greater risk of having a low avidity at this age (p = 0.0425) (Table 6). In the same models, there was no significant group effect on seroconversion or avidity at age 4.5 y.
Table 6. In a randomized trial follow-up study, predictors of seroconversion and antibody avidity in Kenyan infants receiving 5 mg iron daily or no-iron control for 4 months, and given measles vaccine at age 9 mo, measured end of intervention (age 11.5 mo).
Our main findings are, in Kenyan infants: (1) ID and anemia were common at time of vaccination: over 50% of infants were already anemic at age 10 and 14 wk, and 57% had IDA at 24 wk; (2) controlling for most established risk factors affecting vaccine response: (a) Hb at time of vaccination was the strongest positive predictor of anti-diphtheria IgG and anti-pertussis toxin IgG at 24 wk (primary response) and 18 mo (secondary response), and of anti-PnPs 19 IgG at 18 mo; (b) moderate to severe anemia at the time of vaccination more than doubled the risk of being seronegative against diphtheria and PnPs 19 at 18 mo; (c) TfR at 24 wk was a strong negative predictor of seroconversion against diphtheria and PnPs 19 at 18 mo; and (3) iron supplementation of infants at time of measles vaccination at age 9 mo improved the primary vaccine response measured at age 11.5 mo, as indicated by higher anti-measles serum IgG concentrations, seroconversion and IgG avidity.
In the birth cohort study, more than half of infants were already anemic at age 10 wk (Table 1). Term infants born to iron-sufficient mothers should have adequate birth iron stores to cover iron requirements for the first 4 to 6 months (26). However, in Sub-Saharan Africa, 46% of pregnant women have IDA (27), reducing maternal-fetal iron transfer (26), the umbilical cord is often clamped too early, and 15–25% of newborns have LBW (28). These sharply reduce newborn iron stores: it is estimated that body iron is 40–50% lower in newborns who are LBW and/or whose mothers were anemic during pregnancy (26). Low iron stores at birth together with frequent infections increasing serum hepcidin (29), and diarrhea/intestinal parasites causing blood loss, result in many infants depleting their iron stores within 3–4 months after birth. Although representative studies in early infancy are lacking, data from our study and others suggest anemia affects many infants in LMIC at the time of vaccination (7, 8).
Defining iron status in infants at age <24 wk is difficult (26, 30). There are major shifts of iron from fetal Hb to splenic macrophages and hepatocytes and then to a rapidly increasing red cell mass (26, 30). In African infants, common infections increase SF (an acute phase protein) and elevate serum hepcidin, which causes anemia of inflammation (29). TfR levels are also highly variable, with lower values in the first 2 months, and gradually rising values thereafter (30). These patterns were visible in our infants, who despite having low Hb values and high rates of anemia, showed high SF and low TfR concentrations at 10 and 14 wk, but low SF and high TfR concentrations at 24 wk (Table 1). Because there are no consensus normative values for iron biomarkers at age <24 wk (30), we cannot be certain that ID explains the link between anemia with vaccine response at age <24 wk. However, it is likely that most of the early anemia was due to ID, and at 24 wk, when we could begin to apply biomarkers to define iron status, TfR (a measure of iron deficient erythropoiesis) was a predictor of vaccine response (Table 4).
Previous reviews (5, 9, 31) have suggested multiple mechanisms by which iron status might influence adaptive immunity. In some animal and cell models, ID reduces the proportion of mature T-cells and impairs T-cell activation and proliferation, although in many cell models iron chelators have been used (5, 9, 31). Some studies have also shown detrimental effects of ID on B-cell numbers or function, but others have not (5, 9–11, 31). In humans, studies on ID and immune function show varying results depending on what aspect of immunity is measured and the severity of ID, age and/or underlying nutritional status (5, 9). Iron uptake via TfR1 is essential for lymphocyte development, and clinically, a homozygous mutation in TfR1 causes severe immunodeficiency in children and reduced numbers of circulating memory B cells (11). In vitro this mutation prevented T- and B-cell proliferation and addition of high concentrations of iron citrate in vitro rescued the proliferative defect (11). Thus, adequate iron availability may be critical for adaptive immunity. If ID critically limits iron availability to responding lymphocytes (10), this may explain our findings that ID and anemia during early infancy are associated with an impaired response to vaccination. If ID limits the development of memory B-cells (10) it may reduce long-term immunity (32). Recent studies showing impaired antigen-specific immune responses in hypoferremic mice and lower antibody concentrations in patients with iron-refractory IDA support this concept (Joe Frost et al., personal communication 2020).
We are aware of no controlled intervention studies examining the effect of iron supplementation on infant vaccine response in Africa. Most previous studies on iron status and vaccines have been small and cross-sectional, and have produced equivocal results (5). In South Africa, the response to a diphtheria toxoid vaccine was measured in children with IDA, milder ID and non-anemic controls; the numbers with a positive response was 12 of 14 in the controls, 2 of 11 in the anemic children and 0 of 7 in those with ID (33). An observational study in eight LMIC assessing links between enteropathogen infection, undernutrition, and child health (n = 1449) found borderline significant associations between TfR and oral polio vaccine response at age 7 and 15 mo (34). In a cross-sectional study of Ecuadorian children (n = 1,162) who had previously received diphtheria and tetanus vaccines, anemic children older than age 1 y had lower antibody levels than controls and more anemic children were seronegative for diphtheria compared with controls (12.5 vs. 18.1%, respectively) (35).
Because different vaccine types induce different immune responses (36), they may be differentially affected by ID. ID may have greater detrimental effects on T cell numbers and function, rather than on B cells (5), although not all studies agree (10). This may help explain why in our study, anemia predicted response to diphtheria, pertussis and measles vaccine, but not to Hib and most serotypes of PnPs. Measles, a live vaccine, and the protein antigens of diphtheria and pertussis, induce a strong T-dependent response, including induction of antigen-specific Tfh cells and germinal centers producing memory B cells, and these processes are highly iron dependent (10). Unless effectively conjugated with an adjuvant protein, bacterial polysaccharide antigens (e.g., Hib, PnPs) induce primarily T-independent responses and do not generate Tfh cells or germinal centers (37); thus, their immunogenicity, even when conjugated, might be less affected by iron availability.
In the randomized trial, our findings suggest that iron supplementation of mostly IDA infants at time of MV improves the primary response to the vaccine. Nearly all measles deaths globally occur in young children in LMIC countries (38). The Kenyan MoH recommends MV at age 9 mo because measles infection occurs frequently in infancy; however, primary vaccination failures occur in up to 10–15% of infants vaccinated at age 9 mo (38). In the WHO Africa region, median effectiveness of the first MV given at age 9–11 months was 73%, but the range was from 26 to 95%, showing more variability and lower effectiveness than all other WHO regions (39). The development of a high avidity antibody response is critical for protection against measles virus (38) and avidity to MV is generally lower in children vaccinated at 9 mo compared to later ages (38). Our findings suggest IDA may contribute to primary MV failure by reducing both antibody quantity and quality.
The strengths of the birth cohort study include a large sample size and an 18-month prospective design. All infants received three doses (four for OPV) of the identical vaccine administered and recorded by the study team, improving reliability. Hb and anthropometric indices were measured at three vaccination time points during early infancy, and responses to multiple vaccines were measured at multiple time points, which likely minimized misclassification bias and sharpened case definition. We tried to minimize confounding by controlling for most established risk factors affecting vaccine response; however, our results may in part reflect residual confounding not accounted for in the models. In many previous studies linking nutrition to vaccine response, seroprotective levels were reached in nearly all children irrespective of nutritional status (3). In our study, the response to vaccination was variable; for example, only 68% of infants in the birth cohort study had protective anti-diphtheria serum IgG concentrations at age 18 mo, increasing our ability to detect associations.
Our study also has limitations. We did not specifically control for breastfeeding, but in the study area nearly all mothers partially or fully breastfeed during the first 6 mo. Also, we did not control for chronic parasitic infections, but their prevalence is <2% in early infancy in this area (40). Because of the difficulties of defining iron status before age 6 mo, we cannot be sure that ID is the major cause of the anemia in our birth cohort. However, it is likely to be the major cause and thus anemia per se may not be a predictor of vaccine response, but rather only a surrogate for ID. As our main outcomes, we assessed antibody concentrations and seroprotection, but these are only surrogates of disease protection induced by vaccination. The Bacillus-Calmette-Guérin (BCG) vaccine may provide protection against other infectious diseases than tuberculosis (41) but we did not record BCG administration nor measure BCG vaccine response. In the randomized trial, we compared the MV response in infants receiving a MNP with iron to infants receiving the same MNP without iron. It is possible that components of the MNP other than iron (e.g., vitamin A, zinc) could have had an impact on MV response; however, reviews suggest the only clear benefit of MNPs is an improvement in iron status and hemoglobin (42). We also did not measure parameters of cellular immune responses, such as T cell proliferation, cytokine secretion or delayed-type hypersensitivity tests, which could have provided additional information of the immunogenicity of the vaccines.
To our knowledge these are the first prospective data from Africa assessing the impact of anemia and ID at the time of vaccination on response to a range of pediatric vaccines. Our data suggest that anemia/ID at the time of infant vaccination may impair the response to diphtheria and pertussis vaccines and that improving iron status may improve response to measles vaccine. Clearly, our findings need to be confirmed in other prospective cohorts and/or by larger intervention trials. Future research should assess the potential impact of anemia/ID on response to oral vaccines, such as rotavirus vaccine. Powerful emerging techniques combining mass cytometry and systems-level-omics tools may allow identification of the mechanisms underlying the effects of iron status on the infant immune system and response to vaccines. If confirmed, our findings argue strongly for early detection and correction of IDA in infancy. Because anemia is so common in African infants and because the vaccine-preventable disease burden is so high (1, 2), even if IDA only modestly reduces immunogenicity of childhood vaccines its prevention could have major benefits.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The studies involving human participants were reviewed and approved by Study 1: Institutional Review Boards at Kenyatta National Hospital (KNH), Case Western Reserve University and the Stanford University School of Medicine. Study 2: Institutional Review Boards of the KNH/University of Nairobi (KNH/UoN) and the Zurich Cantonal Ethical Commission. Approval for the follow-up study was given by (KNH/UoN). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
FK, IM, AL, and CK conceived and conducted the birth cohort study. NS, MU, EM, DP, SK, FK, and MZ conceived and conducted the randomized trial. All authors provided inputs on the data analyses. NS and MZ wrote the first draft of the manuscript. All authors contributed to the editing and the finalization of the manuscript.
This work was supported by the Bill and Melinda Gates Foundation, DSM Nutrition, the Philhuman Foundation, and ETH Zürich.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the study participants for their sustained participation. We also thank the participating maternity nurses, clinical officers, laboratory technicians, and data managers at Msambweni County Referral Hospital who made this project possible. The birth cohort study was supported by the Bill and Melinda Gates Foundation Healthy Growth Award (OPP1066865) and the ETH Zürich, Switzerland. The randomized trial was supported by the ETH Zürich, DSM Nutrition (Kaiseraugst, Switzerland) and the Philhuman Foundation (Vaduz, Lichtenstein).
1. Foundation BMG. Vaccine Delivery Strategy Overview. Available online at: https://www.gatesfoundation.org/What-We-Do/Global-Development/Vaccine__Delivery (accessed January 3, 2019)
2. Church JA, Parker EP, Kirkpatrick BD, Grassly NC, Prendergast AJ. Interventions to improve oral vaccine performance: a systematic review and meta-analysis. Lancet Infect Dis. (2019) 19:203–14. doi: 10.1016/S1473–3099(18)30602–9
4. Mupere E, Karamagi C, Zirembuzi G, Grabowsky M, de Swart RL, Nanyunja M, et al. Measles vaccination effectiveness among children under 5 years of age in Kampala, Uganda. Vaccine. (2006) 24:4111–5. doi: 10.1016/j.vaccine.2006.02.038
5. Savy M, Edmond K, Fine PE, Hall A, Hennig BJ, Moore SE, et al. Landscape analysis of interactions between nutrition and vaccine responses in children. J Nutr. (2009) 139:2154S−8. doi: 10.3945/jn.109.105312
6. Kassebaum NJ, Jasrasaria R, Naghavi M, Wulf SK, Johns N, Lozano R, et al. A systematic analysis of global anemia burden from 1990 to 2010. Blood. (2014) 123:615–24. doi: 10.1182/blood-2013–06-508325
7. le Cessie S, Verhoeff FH, Mengistie G, Kazembe P, Broadhead R, Brabin BJ. Changes in haemoglobin levels in infants in Malawi: effect of low birth weight and fetal anaemia. Arch Dis Childhood. (2002) 86:182–7. doi: 10.1136/fn.86.3.F182
8. Jaeggi T, Moretti D, Kvalsvig J, Holding PA, Tjalsma H, Kortman GA, et al. Iron status and systemic inflammation, but not gut inflammation, strongly predict gender-specific concentrations of serum hepcidin in infants in rural Kenya. PLoS ONE. (2013) 8:e57513. doi: 10.1371/journal.pone.0057513
10. Jiang Y, Li C, Wu Q, An P, Huang L, Wang J, et al. Iron-dependent histone 3 lysine 9 demethylation controls B cell proliferation and humoral immune responses. Nat Commun. (2019) 10:2935. doi: 10.1038/s41467–019-11002–5
11. Jabara HH, Boyden SE, Chou J, Ramesh N, Massaad MJ, Benson H, et al. A missense mutation in TFRC, encoding transferrin receptor 1, causes combined immunodeficiency. Nat Genet. (2016) 48:74–8. doi: 10.1038/ng.3465
12. Paganini D, Uyoga MA, Kortman GAM, Cercamondi CI, Moretti D, Barth-Jaeggi T, et al. Prebiotic galacto-oligosaccharides mitigate the adverse effects of iron fortification on the gut microbiome: a randomised controlled study in Kenyan infants. Gut. (2017) 66:1956–67. doi: 10.1136/gutjnl-2017–314418
13. McKittrick ND, Malhotra IJ, Vu DM, Boothroyd DB, Lee J, Krystosik AR, et al. Parasitic infections during pregnancy need not affect infant antibody responses to early vaccination against Streptococcus pneumoniae, diphtheria, or Haemophilus influenzae type B. PLoS Negl Trop Dis. (2019) 13:e0007172. doi: 10.1371/journal.pntd.0007172
14. Erhardt JG, Estes JE, Pfeiffer CM, Biesalski HK, Craft NE. Combined measurement of ferritin, soluble transferrin receptor, retinol binding protein, and C-reactive protein by an inexpensive, sensitive, and simple sandwich enzyme-linked immunosorbent assay technique. J Nutr. (2004) 134:3127–32. doi: 10.1093/jn/134.11.3127
15. van Gageldonk PG, van Schaijk FG, van der Klis FR, Berbers GA. Development and validation of a multiplex immunoassay for the simultaneous determination of serum antibodies to Bordetella pertussis, diphtheria and tetanus. J Immunol Methods. (2008) 335:79–89. doi: 10.1016/j.jim.2008.02.018
16. Smits GP, van Gageldonk PG, Schouls LM, van der Klis FR, Berbers GA. Development of a bead-based multiplex immunoassay for simultaneous quantitative detection of IgG serum antibodies against measles, mumps, rubella, and varicella-zoster virus. Clin Vaccine Immunol. (2012) 19:396–400. doi: 10.1128/CVI.05537–11
17. Schepp RM, Berbers GA, Ferreira JA, Reimerink JH, van der Klis FR. A novel multiplex poliovirus binding inhibition assay applicable for large serosurveillance and vaccine studies, without the use of live poliovirus. J Virol Methods. (2017) 241:15–23. doi: 10.1016/j.jviromet.2016.12.006
18. Brinkman ID, de Wit J, Smits GP, Ten Hulscher HI, Jongerius MC, Abreu TC, et al. Early measles vaccination during an outbreak in the netherlands: short-term and long-term decreases in antibody responses among children vaccinated before 12 months of age. J Infect Dis. (2019) 220:594–602. doi: 10.1093/infdis/jiz159
20. WHO. Haemoglobin Concentrations for the Diagnosis of Anaemia and Assessment of Severity. WHO (2011). Available online at: http://www.who.int/vmnis/indicators/haemoglobin.pdf (accessed July 03, 2018)
22. Hefele L, Syphan S, Xayavong D, Homsana A, Kleine D, Chanthavilay P, et al. Seroprotection at different levels of the healthcare system after routine vaccination with diphtheria-tetanus-pertussis whole cell-hepatitis B-Haemophilus influenzae Type B in Lao People's Democratic Republic. Clin Infect Dis. (2019) 69:2136–44. doi: 10.1093/cid/ciz143
24. Chen RT, Markowitz LE, Albrecht P, Stewart JA, Mofenson LM, Preblud SR, et al. Measles antibody: reevaluation of protective titers. J Infect Dis. (1990) 162:1036–42. doi: 10.1093/infdis/162.5.1036
25. Giri S, Kumar N, Dhanapal P, Venkatesan J, Kasirajan A, Iturriza-Gomara M, et al. Quantity of vaccine poliovirus shed determines the titer of the serum neutralizing antibody response in indian children who received oral vaccine. J Infect Dis. (2018) 217:1395–8. doi: 10.1093/infdis/jix687
27. The World Bank. Prevalence of Anemia Among Pregnant Women (%). (2017). Available online at: https://data.worldbank.org/indicator/sh.Prg.Anem (accessed August 8, 2019).
28. UNICEF. Low Birthweight. (2014). Available online at: https://data.unicef.org/topic/nutrition/low-birthweight/ (accessed August 8, 2019)
29. Mupfudze TG, Stoltzfus RJ, Rukobo S, Moulton LH, Humphrey JH, Prendergast AJ, et al. Plasma concentrations of hepcidin in anemic zimbabwean infants. PLoS ONE. (2015) 10:e0135227. doi: 10.1371/journal.pone.0135227
33. Macdougall LG, Anderson R, McNab GM, Katz J. The immune response in iron-deficient children: impaired cellular defense mechanisms with altered humoral components. J Pediatr. (1975) 86:833–43. doi: 10.1016/S0022–3476(75)80211–3
34. Pan WK, Seidman JC, Ali A, Hoest C, Mason C, Mondal D, et al. Oral polio vaccine response in the MAL-ED birth cohort study: considerations for polio eradication strategies. Vaccine. (2019) 37:352–65. doi: 10.1016/j.vaccine.2018.05.080
35. Brussow H, Sidoti J, Dirren H, Freire WB. Effect of malnutrition in Ecuadorian children on titers of serum antibodies to various microbial antigens. Clin Diagn Lab Immunol. (1995) 2:62–8. doi: 10.1128/CDLI.2.1.62–68.1995
38. WHO. Measles Vaccine: WHO position paper. WHO (2020). Available online at: http://www.who.int/wer (accessed January 18, 2020).
39. Ma SJ, Li X, Xiong YQ, Yao AL, Chen Q. Combination measles-mumps-rubella-varicella vaccine in healthy children: a systematic review and meta-analysis of immunogenicity and safety. Medicine. (2015) 94:e1721. doi: 10.1097/MD.0000000000001721
40. LaBeaud AD, Nayakwadi Singer M, McKibben M, Mungai P, Muchiri EM, McKibben E, et al. Parasitism in children aged three years and under: relationship between infection and growth in rural coastal Kenya. PLoS Negl Trop Dis. (2015) 9:e0003721. doi: 10.1371/journal.pntd.0003721
41. Jensen KJ, Biering-Sorensen S, Ursing J, Kofoed PL, Aaby P, Benn CS. Seasonal variation in the non-specific effects of BCG vaccination on neonatal mortality: three randomised controlled trials in Guinea-Bissau. Bmj Glob Health. (2020) 5:e001873. doi: 10.1136/bmjgh-2019–001873
42. Suchdev PS, Jefferds MED, Ota E, da Silva Lopes K, De-Regil LM. Home fortification of foods with multiple micronutrient powders for health and nutrition in children under two years of age. Cochrane Database Syst Rev. (2020) 2:CD008959. doi: 10.1002/14651858.CD008959.pub3
Keywords: iron deficiency, anemia, iron, vaccine response, seroconversion, infancy, Kenya
Citation: Stoffel NU, Uyoga MA, Mutuku FM, Frost JN, Mwasi E, Paganini D, van der Klis FRM, Malhotra IJ, LaBeaud AD, Ricci C, Karanja S, Drakesmith H, King CH and Zimmermann MB (2020) Iron Deficiency Anemia at Time of Vaccination Predicts Decreased Vaccine Response and Iron Supplementation at Time of Vaccination Increases Humoral Vaccine Response: A Birth Cohort Study and a Randomized Trial Follow-Up Study in Kenyan Infants. Front. Immunol. 11:1313. doi: 10.3389/fimmu.2020.01313
Received: 30 April 2020; Accepted: 26 May 2020;
Published: 13 July 2020.
Edited by:Christoph T. Berger, University of Basel, Switzerland
Reviewed by:Gyaviira Nkurunungi, Medical Research Council, Uganda
Emmanuel Hermann, Université de Lille, France
Copyright © 2020 Stoffel, Uyoga, Mutuku, Frost, Mwasi, Paganini, van der Klis, Malhotra, LaBeaud, Ricci, Karanja, Drakesmith, King and Zimmermann. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
†These authors have contributed equally to this work