The potential of spot urine as a biomarker for zinc assessment in Malawian children and adults

Introduction

Zinc is an essential micronutrient, required in the diet in relatively small amounts, yet it is involved in many metabolic and physiological processes. Biochemical data on zinc deficiency are sparse, however, the limited data on population-level analyses of plasma or serum zinc (SZC) levels have shown that zinc deficiency is of global public health concern (1). Zinc deficiency may lead to poor immunity and greater susceptibility to infections. It can also cause poor growth, poor wound healing, birth defects, and death in infants (2). Zinc deficiency underlies a substantial disease burden predominantly among children under 5 years of age. It affects one in every two children from low- and middle-income countries and is associated with reduced linear growth as well as increased incidences of infectious diseases (3).

Plasma or SZC is used as a biomarker to assess population-level zinc status but it has limitations (2, 4) including confounding factors such as acute inflammation which cause artificial inflation of the prevalence estimates (5–8). Adjustments can be made to account for the effects of inflammation, however, the existing methods are not consistent with each other (9, 10). There are currently three established methods that are used to correct for inflammation but the prevalence estimates that are produced by these methods differ. Other biomarkers such as nails, hair, and sweat, have been proposed but are not routinely used due to multiple limitations such as biological factors and environmental contamination (2, 11). Single-void (spot/casual), urine is an established biomarker that is used for population-level assessment of iodine (12), and it has shown value in the population-level assessment of selenium (13). It is thus a good candidate for assessing population zinc status on a large scale but there are no established cut-off values for spot urinary zinc concentration (UZC) which limits its usefulness. Prior studies have used a factorial approach to calculate the physiological requirements of zinc in healthy participants, which involves identifying typical values for urinary zinc losses in healthy individuals by age range (14, 15). These values may be useful in establishing cut-off points; however, they were based on 24-h urinary samples and may not be appropriate for spot urine samples.

A strong correlation between paired data of an established biomarker with known cut-offs indicating deficiency and a surrogate with no known cut-offs can be used to identify thresholds for deficiency in the latter (16, 17). This method is used in blood serum biomarkers for vitamin A deficiency, but random, or “spot” urine samples have also been useful surrogates for other nutrients. For instance, spot urine selenium has been identified as a strong candidate for determining population-level selenium status in specific demographic groups because of its strong correlation with paired-sample serum selenium concentration (13). We, therefore, hypothesized that a strong correlation between UZC and SZC could also be used to determine cut-off points for UZC that are equivalent to those in SZC. The objectives of this study were (1) To assess the strength of the relationship between UZC data against SZC data in the Malawi Micronutrient Survey population using linear mixed-effects modeling; (2) If a strong correlation between SZC and UZC data is found to identify the optimal cut-off point in UZC for estimating the risk of zinc deficiency.

Materials and methods

We used secondary data from the Malawi National Micronutrient Survey (MNS) which was conducted from 2015 to 2016. The study design has been previously reported in detail (18) and is presented here in brief. The survey used a two-stage cluster design which aimed to provide results weighted at national, regional (north, central and south), district (n = 28), and residency (rural/urban) levels. This was done using household-level sampling weights, which were derived from the Malawi Demographic Health Survey. The sampling method resulted in a total of 105 clusters (35 in each of the regions) from which households were selected for inclusion in the MNS study (20 per urban cluster and 22 per rural cluster). In each household, all primary school children (PSC) were included, and from nine random households of the 20/22 households, all women of reproductive age (WRA) were included. Of the nine random households, six were randomly selected for inclusion of all school-age children (SAC), and from these six households, four were included for adult men.

The sample size taken for this current study was determined by the availability of the urine data, which was determined by the sample availability. The spot urine samples were only collected from WRA and SAC for the assessment of iodine concentrations and other conditions. The urine sample collection method is described in detail in an earlier publication (13), while sample collection and analysis of the SZC are detailed in the MNS report (18). Serum samples for zinc assessment were collected in trace element-free vacutainers. The samples were maintained in portable freezers in the field and transported to the nearest district laboratory for temporary storage (−20°C). The samples were then transferred to a central laboratory and stored at −70°C until shipment to CHORI lab in Oakland, for analysis. In the lab, serum zinc was diluted to a final concentration of 5% HNO₃ and centrifuged at 3,000 × g for 10 min before analysis using Atomic Emission Spectrometry (AES) (19, 20).

Subsequently, zinc concentrations were measured in stored urine samples, as reported previously (21). In brief, the urine samples were frozen within an hour of collection and kept at −20°C for a maximum of 8 h, then they were transferred to a central laboratory in Malawi which stored them at −80°C (13). The samples were then transferred frozen and on dry ice, to the UK where they were also kept at −80°C until the time of sample analysis. Zinc concentration was determined in a multi-element analysis (with a Helium collision mode, to reduce interferences) using the Inductively Coupled Plasma Mass Spectrometry (ICP-MS) assay (21). The limit of detection was 19 μg/L based on the standard deviation of 137 analytical blanks. Certified reference materials of human urine and seronorm trace elements of urine were used as quality controls for the analysis. In total, zinc concentration was measured in 1,406 urine samples (n = 741 and n = 665 participants for WRA and SAC, respectively). This study included only observations that had serum zinc data, inflammation adjustment data, and urinary data.

The UZC data were adjusted for effects of hydration (see Section Adjustment of urine data for hydration) before running the linear-mixed effects models. The SZC data used in the model was adjusted for inflammation (Section Adjustment of serum zinc data for inflammation) in a prior study (9).

Adjustment of urine data for hydration

Creatinine adjustment was conducted using the equation below

$\begin{array}{l}UZ{C}_{cr}=\frac{UZ{C}_{unadj}}{C{r}_{meas}}\end{array}$

Where UZC_cr denotes adjusted UZC and was reported as UZC_cr μg/L, UZC_unadj denotes unadjusted UZC, and CR_meas denotes, the measured creatinine concentration of the individual.

Osmolarity adjustment was achieved by the equation

$\begin{array}{l}UZ{C}_{osm}=UZ{C}_{unadj}*\frac{OS{M}_{mean}}{OS{M}_{meas}}\end{array}$

Where UZC_osm denotes osmolarity-adjusted UZC and is reported as UZC_osm μg/L, UZC_unadj denotes the unadjusted UZC, OSM_mean denotes the mean osmolarity from all the data and OSM_meas denotes the individual osmolarity measurement.

These adjustments were all made separately for the WRA and SAC data and resulted in four different UZC data points for each individual (unadjusted, Sg adjusted, creatinine adjusted, and osmolarity adjusted). A urine sample quality check was made to identify samples that were too dilute for inclusion in the analysis. This was conducted using a combined cut-off mark of creatinine level ≤ 0.05 g/L and an Sg of ≤ 1.001 (22).

Adjustment of serum zinc data for inflammation

Biomarkers of inflammation C-reactive protein (CRP) and alpha-1 acid glycoprotein were analyzed in serum samples that were aliquot from the blood collected in all participants. This analysis was conducted using the sandwich Enzymatic Linked Immunosorbent Assay (ELISA) at the VitMin lab in Germany (18). Three forms of the serum zinc data were used for the analyses: (1) Raw uncorrected data (SZC_unadj μg/dl), (2) corrected for inflammation using internal correction factors (ICF, SZC_icf μg/dl), and (3) corrected for inflammation using Biomarkers Reflecting Inflammation and Nutritional Determinants of Anemia (BRINDA, SZC_br μg/dl) recommended methods (5, 6, 23).

The linear mixed-effects model

The relationship between UZC and SZC was examined by fitting a linear model. Because of the clustered structure of the observations, a linear mixed-effects model (LMM) was used with random effects to account for nested variation at three levels: between cluster, within-cluster (between households), and within households (between individuals). Log-transformed data were used for both SZC and UZC to make the assumption of normality of the random effects plausible.

The basic LMM takes the form

$y=X\text{β}+\text{η}+\text{κ}+\text{ε}$

where the vector contains observation of the dependent variable, the design matrix X contains covariates and the vector β contains the corresponding fixed effects coefficients, and the terms η, κ, and ε are, respectively, the between-cluster, between household (within the cluster), and between individual (within the household) random effects. Each of the random effects has a mean of zero and a variance that must be estimated.

In the so-called null model, the only fixed effect is a constant mean, so the design matrix X contains a single column, with values all 1, and the vector β contains the estimated mean value. We used the lme function from the nlme package in R (24) to estimate the variances of the random effects by residual maximum likelihood. In an alternative model, UZC is considered as a potential predictor for SZC, and so it is a fixed effect in the model, along with a constant intercept. In this case, the design matrix X comprises a column, with values all 1, and a second column that contains the values of UZC. The vector β contains a constant or intercept and the regression coefficient for UZC in the estimated model. The variance components from the predictor model were then compared with those from the null model to determine the magnitude of reduction that was introduced by the addition of UZC data to the model. The predictor model also gave adjusted R², at the nested levels which are similar to the standard R² but not guaranteed to be bound to the interval (0, 1). A higher adjusted R² at the cluster level than the individual level would imply that UZC may have greater value as a population-level, rather than individual-level, biomarker, and vice versa. For the models which gave such outputs, cluster means of the two datasets were plotted to visualize the relationship. The SAC and WRA data were analyzed separately, which included all datasets for UZC (4) and SZC (3). As such, 12 models for WRA and SAC were computed and assessed.

Post-hoc analyses

Post-hoc analyses were undertaken to understand the relationship between UZC and SZC and to further explore potential cut-off points for zinc status from the UZC data. First, a regression model was used to explore the impact of inflammation on UZC. The regression model used log-transformed UZC data as the independent variable and the inflammatory markers c-reactive protein (CRP) and α-1 acid glycoprotein (AGP) as predictor variables. This was done for both SAC and WRA and all the three UZC data types.

Second, attempts were made to achieve objective 2 by using values based on estimated average urinary zinc losses. These values were estimated from published literature by expert panels for the estimation of zinc physiological requirements (PR) during the process of arriving at dietary zinc recommendations using the factorial approach. The hypothesis is that if dietary zinc intake is not sufficient to meet the PR, then homeostatic adjustments to preserve zinc occur, including reduced zinc losses in urine (2). Urinary zinc losses below these thresholds may thus be indicative of inadequate zinc intake. Age and gender-specific values for zinc loss via urine were used for SAC (0.1 mg/day for SAC aged 5–6, 0.2 mg/day for SAC aged 7–10 years, 0.3 mg/day for male and female SAC aged 11–14 years, after rounding), and for women, the threshold of 0.3 mg/day WRA (25) was used. The values were published by European Food Safety Authority (EFSA), but are based on reference bodyweight values from a study by van Buuren et al. (for SAC) (26), and multiple individual studies for WRA (27–32). For control, the threshold value for zinc loss that is recommended by the Food and Nutrition Board (FNB)/ Institute of Medicine (IOM) of 0.44 mg/day in WRA, was also tested. This analysis was done using the hydration-adjusted data only, as it implied conversion of the spot urine samples to equivalent 24-h values. The estimated prevalence of observations with reduced dietary intake was then compared with the estimated prevalence of serum zinc deficiency using a Chi-squared test of proportions.

Statistical inferences

For all analyses, a p-value of 0.05 was considered statistically significant and prevalence estimates are reported with a 95% confidence interval. All point estimates were weighted using the household-level sampling weights which were derived from the survey metadata.

Data access and ethical approval

Urinary zinc data was accessed from Phiri et al. (21) whereas the serum zinc data was accessed from the DHS program database following approval of the request for use (May 23rd, 2019). Ethical approval for this secondary data analysis was obtained from the College of Medicine Research Ethics Committee (COMREC), protocol number P.01/20/2915 (March 11th, 2020). COMREC is affiliated with Kamuzu University of Health Sciences, previously known as the College of Medicine in Malawi.

Results

Participants and sample characteristics

Spot urine zinc analysis was conducted for n = 1,406 participants (WRA: n = 741 and SAC: n = 665), whereas serum zinc analysis was conducted for n = 1,488 (WRA: n = 745 and SAC: n = 743). After both datasets were cleaned and matched n = 1,301 observations remained of which 660 were WRA and 641 were SAC as shown in Table 1 and Supplementary Figure S1.

TABLE 1

Table 1. Participant demographics and geometric mean +SD values for urine and plasma biochemical parameters measured in women and school-age children in the 2015–2016 Malawi national micronutrient survey.

For SZC, an upper threshold of 125 μg/dl was used (33), beyond which samples were suspected to have been contaminated and deemed as outliers. All urine samples were deemed of good quality and not too dilute for the element analyses because none of the observations had both creatinine and Sg levels that were below the thresholds of dilution. Only two observations (one in each demographic group) had creatinine values below 0.05 g/L and there were no observations with Sg below 1.001.

Adjustment for urine data for hydration

The mean Cr, Osm, and Sg values that were used for the hydration adjustments are presented in Table 1. There were higher mean concentrations of Cr and Osm in WRA compared to SAC. The data for both SAC and WRA was not normally distributed before or after hydration adjustment, based on Shapiro–Wilk tests which gave extremely low p-values (p < 0.01). The mean UZC in SAC decreased following adjustment, from 410.95 to 396.39 and 385.51 μg/L for Osm and Sg, respectively. These changes were not significant, based on Wilcoxon signed-rank tests. Adjusting for creatinine, significantly increased the mean UZC of SAC to 500.23 μg/L (p < 0.01; Figure 1). Hydration adjustment of the UZC of WRA led to decreases in the UZC for all three methods. Wilcoxon signed-rank tests indicated that the hydration adjusted values for WRA were significantly lower than the unadjusted values (UZC_cr; p < 0.01, UZC_osm; p < 0.01, UZC_sg; p = 0.01).

FIGURE 1

Figure 1. Bar chart with box-plots for urine zinc concentration values and means among women of reproductive age and school-age children, before and after adjustment for hydration status in the 2015–2016 Malawi national micronutrient survey. UZC_unadj: unadjusted urine zinc concentration, UZC_cr: Creatinine adjusted urine zinc concentration, UZC_osm: Osmolarity adjusted urine zinc concentration, UZC_sg: Specific gravity adjusted urine concentration, SAC: School-age children, WRA: Women of reproductive age.

Linear mixed-effects model results of SZC and UZC data

As shown in Table 2, in the null model of SZC, variance components were very small for SAC for hydration-adjusted and unadjusted data. This was true for all three sampling levels and the largest (albeit still small) variance was observed between individuals within the household level for the unadjusted SZC data. The addition of urinary zinc as random terms to the SZC model did not improve the efficiency of the model for SAC. This is based on the inability of the UZC data to reduce the variance component of the null model. This observation was true for all model combinations of all three SZC types and all four UZC types. A reduction in the variance of the adjustment for inflammation reduced the variance slightly but the change was negligible. No variation in the SZC data could thus be explained by the addition of UZC data of any hydration adjustment type.

TABLE 2

Table 2. Variance components for the linear mixed-effects model (random effects and fixed effects) with serum zinc concentration as the dependent variable and urine zinc concentration as the predictor variable for school-age children in the 2015–2016 Malawi national micronutrient survey.

For SAC, the adjusted R² was largest at between cluster level compared to between household or between individual within household levels, for all three SZC types. The adjusted R² increased when UZC was adjusted for hydration and the largest value of adjusted R² (0.21) is observed for UZC_sg as a predictor for unadjusted SZC data. This largest adjusted R² was observed at the cluster level and the cluster means for this data type were plotted against those of the SZC data types (Figure 2). A linear relationship can be seen between the cluster means of UZC_sg and the SZC data types in all three graphs in Figure 2. Overall, adjustment for inflammation did not improve the relationship between SZC and UZC in SAC data.

FIGURE 2

Figure 2. Cluster means of urinary zinc concentrations adjusted for specific gravity (UZC_sg) plotted against serum zinc concentrations which are unadjusted and adjusted for inflammation (using internal correction factor and biomarkers reflecting inflammation and nutritional determinant of anaemia methods) in school-aged children. UZC_sg: Urinary zinc concentration adjusted for specific gravity, SZC: serum zinc concentration, SZC_unadj: Unadjusted serum zinc data, SZC_icf: Serum zinc concentration adjusted using internal correction factor method (ICF), SZC_br: Serum zinc concentration adjusted using biomarkers reflecting inflammation and nutritional determinant of anaemia method (BRINDA).

For WRA, the variance between household levels was extremely low for all the models; as such, the adjusted R² at this level could not be interpreted (Table 3). At between cluster level, the addition of UZC_cr to the null model did not affect the variance and the adjusted R² for this data type was equal to 0. The other data types also had very low adjusted R²-values which indicated a lack of improved efficiency despite the addition of the UZC data to the null model. Similar results were observed among individuals, within the household level where low adjusted R²-values were observed. Adjustment for hydration did not have any effect on the variation of the null model, for all SZC data types at the varying levels.

TABLE 3

Table 3. Variance components for the linear mixed-effects model (random effects and fixed effects) with serum zinc concentration and the dependent variable and urine zinc concentration as the predictor for women in the 2015–2016 Malawi national micronutrient survey.

Post-hoc analyses

The multivariate linear regression revealed no significant relationship between inflammatory markers and urinary zinc data across all the models and for SAC (Table 4). For WRA, there was a significant relationship between AGP and UZC_unadj, but with a coefficient of 0. Thus, no variation in UZC data of all types could be explained by the inflammatory markers.

TABLE 4

Table 4. Fixed effects inferential of a linear mixed model (random effects and fixed effects) with urine zinc concentration as the dependent variable and inflammatory markers as the predictor variables for school-age children and women in the 2015–2016 Malawi national micronutrient survey.

For SAC, the EFSA threshold of urinary zinc excretion (0.1–0.3 mg/day) gave low estimates of zinc deficiency compared to SZC for all the data types. The lowest estimate was observed with UZC_cr data and the highest, with UZC_unadj (Table 5). There was no overlap in the confidence intervals between the UZC and the SZC estimates.

TABLE 5

Table 5. Prevalence estimates of zinc deficiency in urinary zinc concentration, and serum zinc concentration for school-age children and women in the 2015–2016 Malawi national micronutrient survey.

For WRA, the urinary zinc threshold value of 0.3 mg/day, suggested by EFSA gave prevalence estimates of zinc deficiency ranging from 45 to 50% (Table 5). These were lower than the estimates of zinc deficiency as identified by SZC before and after adjusting for inflammation. The lower tail of the confidence intervals for all the EFSA threshold estimates for SZC had an overlap with the upper tails of the confidence intervals of UZC_unadj, UZC_cr, and UZC_osm. The prevalence estimates from the threshold by FNB/IOM (0.44) were higher than the estimates for SZC. These estimates had a slight overlap with the higher tail of the confidence intervals of the SZC estimates.

Discussion

The assessment of zinc in humans at the population level remains a challenge due to the lack of suitable zinc biomarkers (2). Spot urine is a good candidate for determining zinc deficiency at the population level, but the lack of threshold for defining deficiency in this biomarker poses a limitation to its adoption and use. This current study aimed to use paired UZC and SZC data, to identify thresholds in urine that are equivalent to those that are used to define zinc deficiency in serum. This was dependent on a strong relationship between these two biomarkers to allow for the adoption of the thresholds.

The paired UZC and SZC data showed a poor correlation in a linear mixed-effects model of SAC and WRA data from a population-based survey. The correlation analysis factored in the sampling design to take into account the differences between clusters, differences between households, and differences between individuals within a household. Before the analysis, the UZC data was adjusted for hydration using three different methodologies, and the SZC was adjusted for effects of inflammation using two different methods. All these confounder corrections were applied to correct for biological differences that occur due to infection and variation in hydration patterns and to optimize the performance of the UZC data in the correlation test. Various factors may have contributed to the poor correlation despite the adjustments for hydration and inflammation. A previous study that used this same SZC data identified possible data quality issues across all demographic groups which resulted in clumping of the SZC values, i.e., more than 50 observations with the same SZC value (9), possibly due to factors associated with limits of detection.

As shown in Table 3, about 20% of the variation in the SZC data for SAC could be explained by the UZC_sg data and this was true for unadjusted as well as inflammation-adjusted SZC data. The correlation values were similar but lower for the other hydration adjustment data. Unadjusted UZC data for SAC also showed a stronger relationship at between household level, but this was not the case between individuals. For WRA, the correlation was weak for all the UZC data types against all the SZC data types. This implies that UZC may be a better population-level marker for SAC than it is for WRA. It is not clear from the results, why the correlation strength differed between SAC and WRA, but physiological differences between these two groups could have contributed to this. Similar results were found by Phiri et al. (21), who assessed spatial variation in this same spot urine zinc data. Their results showed spatial dependence between clusters in UZC for SAC but not for WRA indicating the possible influence on the SAC-UZC by environmental factors such as soil type which has been reported to play an important role in the zinc content of locally grown crops (34).

Various adjustment methods were applied to try and optimize the performance of the linear regression model. Acute infection and inflammation are common factors that affect the homeostasis of various elements in the body as such, in participants with high levels, of inflammation, it is necessary to adjust for their effects. The serum zinc data were adjusted to account for inflammation but this had little effect on the regression statistics. Unadjusted SZC data for SAC had similar and weak correlations at all cluster levels, for all UZC types. The urinary data was also corrected for hydration but results showed that unadjusted data had a stronger correlation with SZC than creatinine or osmolality-adjusted urine data, for the SAC analysis (Table 3).

Hydration adjustment improved the relationship between UZC and SZC for SAC, the strongest being UZC_sg, followed by UZC_osm, and finally, UZC_cr, which indicates clear effects of hydration on the data. There was no discernible pattern of the hydration adjustment methods for WRA and there seemed to be many variations, especially for creatinine-adjusted data. For this demographic group, the performance of the unadjusted UZC data did not differ from UZC_sg and UZC_osm. This indicates possible challenges with hydration adjustments in specific demographic groups, which may have affected the correlation. Although analyses of 24-h urine samples do not need to be adjusted for hydration status or creatinine concentration, they also have poor correlations with SZC data indicating that urine itself may just be a poor biomarker for zinc deficiency regardless of hydration effects (35). However, considering that these analyses were done in the 1970s, new research with improved assay technologies may be warranted.

The lack of agreement between the existing adjustment methods for both inflammation and hydration corrections implies the need to improve the performance of the individual methods. This is a limitation of their use when adjusting various data and could have contributed to the poor linear regression performance. Measurement error can also contribute to the poor performance of a linear regression equation. In this analysis measurement error in either of the two main datasets (UZC and SZC), the inflammatory data (CRP and AGP) or the hydration indicators (Creatinine, Osmolality, Urine specific gravity) could have caused this. However, the large sample sizes control for random error, and available QC data for most of the datasets did not indicate any systematic errors (36).

The weak correlation between UZC and SZC indicates that UZC may have limited value as a biomarker of population-level zinc deficiency even for SAC. There is thus still a need to find a biomarker for zinc that is sensitive, specific, and easily accessible. SZC is the most widely used biomarker for zinc, however, it requires venipuncture blood which can only be collected by skilled professionals, and requires specific techniques. Furthermore, extraction of the serum should be done in a sterile environment which is difficult to adhere to in large scale surveys where samples are collected in the field (37). A biomarker with easier sample collection and management (such as UZC) is therefore preferred. It would allow the collection of samples from large sample sizes with fewer technical limitations. A post-hoc analysis was employed to further explore the relationship between UZC and SZC data. Acute inflammation has been previously shown to cause hyperzincuria (38), which can result in a poor correlation between UZC and SZC since SZC is a negative acute-phase reactant (39, 40). There were high levels of inflammation in SAC (34%), and some inflammation in WRA at 14% (18), however no statistically significant associations with UZC in either the WRA or SAC groups were found.

An alternative approach to identify thresholds for zinc deficiency in urine was explored using average urinary zinc losses in healthy individuals (14, 15) as putative cut-off values. The EFSA thresholds yielded lower prevalence values for zinc deficiency than those calculated using SZC in both SAC and WRA groups. These zinc loss values (0.3 mg/day for WRA and 01.0–0.3 mg/day for SAC) by EFSA (15) were determined by using the factorial approach, based on endogenous zinc losses as assessed by multiple individual studies, most of which were done in the early 2000s.

Other panels have also estimated mean urinary zinc loss and found similar values. For instance, the World Health Organization (WHO) estimates 0.3 mg/day (41) whereas, FNB/IOM estimated urinary zinc losses to be 0.44 mg/day in women, based on studies mostly conducted in the early 1990s (42). This value was also adopted by the IZiNCG (14). When applied to the WRA data the value suggested by FNB/IOM and IZiNCG gave prevalence estimates that were slightly higher than those obtained by SZC. For UZC_Osm and UZC_Usg, the confidence intervals of the prevalence estimates were in the same range as the SZC estimates. With novel primary research, there may be a possibility of determining appropriate thresholds indicating zinc deficiency based on zinc loss, or other means.

The main limitation of this study is that it failed to identify the optimal cut-off point in UZC for estimating the risk of zinc deficiency. This was because there was a poor correlation between the UZC and SZC datasets, which could have occurred because of the quality of the lab results. Examination of the quality control data indicated that there was no measurement error in all the lab analyses of the main datasets. However, for SZC, the assay method used (atomic emission spectrometry) has a higher limit of detection, compared to modern methods such as ICP-MS which was used for UZC analysis. In populations with low levels of zinc deficiency, as the sample for this study, higher limits of detection could lead to data clumping on the lower end of the dataset (9), thus leading to differences in paired data. In the future, high-precision assay instruments should be used for element analysis to ensure that paired data results are measured similarly.

Paired spot urinary and plasma/serum mineral concentrations have been used to determine probable thresholds for mineral deficiencies including vitamin A and selenium. However, for zinc, the poor correlation between these two datasets in a national sample of the Malawi population, despite correcting for multiple confounders, raises questions about whether this approach is appropriate for determining the threshold indicating deficiency for this mineral. The slightly better correlation between UZC and SZC for SAC only, at the cluster level, indicates that spot urine samples may have value at the population level for SAC but not WRA. Nonetheless, it is not clear to what extent this potential can be useful for determining zinc status. An alternative approach to estimate zinc deficiency prevalence—using mean urinary zinc excretion as a cut-off point for deficiency—also gave conflicting results compared to estimates based on SZC, which highlights a need for primary research on appropriate thresholds of zinc concentration in spot urine. Further investigation is also required to understand the value of spot urine zinc at the population level. It is also of great interest to explore if the spatial patterns that are observed in SAC but not WRA, can also be observed in other demographic groups. This could further inform on the applications of spot urine zinc in population zinc status measurement.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: Serum zinc data can be obtained from https://dhsprogram.com/data/dataset/Malawi_Standard-DHS_2015.cfm?flag=0 and Urinary zinc data, can be obtained from the authors of https://link.springer.com/article/10.1007/s10653-020-00700-5.

Ethics statement

The studies involving human participants were reviewed and approved by Malawi College of Medicine Research Ethics Committee (COMREC), Protocol Number P.01/20/2915. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions

Conceptualization and writing—review and editing: BL, RL, NL, and EJ. Data curation: BL and RL. Formal analysis, investigation, and writing—original draft: BL. Funding acquisition and project administration: EJ. Methodology: BL, RL, and NL. Supervision: KM and JP. Validation: RL. All authors contributed to the article and approved the submitted version.

Funding

This work was supported, in whole or in part, by the Bill & Melinda Gates Foundation (INV-009129). Under the grant conditions of the Foundation, a Creative Commons Attribution 4.0 Generic License has already been assigned to the Author Accepted article version that might arise from this submission. The funder had no role in the design, execution, analysis, or interpretation of the data.

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.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.890209/full#supplementary-material

References

1. Hess SY. National risk of zinc deficiency as estimated by national surveys. Food Nutr Bull. (2017) 38:3–17. doi: 10.1177/0379572116689000

PubMed Abstract | CrossRef Full Text | Google Scholar

2. King JC, Brown KH, Gibson RS, Krebs NF, Lowe NM, Siekmann JH, et al. Biomarkers of nutrition for development (BOND)- zinc review. J Nutr. (2016) 146:858S−85S. doi: 10.3945/jn.115.220079

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Victora CG, Christian P, Vidaletti LP, Gatica-Domínguez G, Menon P, Black RE. Revisiting maternal and child undernutrition in low-income and middle-income countries: variable progress towards an unfinished agenda. Lancet. (2021) 397:1388–99. doi: 10.1016/S0140-6736(21)00394-9

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bates C, Bogin B, Holmes B. Nutrition assessment methods. In: Geissler C, Powers H editors. Human Nutrition. Oxford: Oxford University Press UK. p. 613–46.

Google Scholar

5. Namaste SM, Aaron GJ, Varadhan R, Peerson JM, Suchdev PS. Methodologic approach for the biomarkers reflecting inflammation and nutritional determinants of anemia (BRINDA) project. Am J Clin Nutr. (2017) 106:333S−47S. doi: 10.3945/ajcn.116.142273

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Thurnham DI, Mccabe GP. Influence of infection and inflammation on biomarkers of nutritional status with an emphasis on vitamin A and iron. In: Report: Priorities in the Assessment of Vitamin A and Iron Status in Populations, Panama City, Panama, 15-17 September 2010, Geneva: World Health Organization, p. 63–80.

Google Scholar

7. Macdonell SO, Miller JC, Harper MJ, Reid MR, Haszard JJ, Gibson RS, et al. comparison of methods for adjusting biomarkers of iron, zinc, and selenium status for the effect of inflammation in an older population: a case for interleukin 6. Am J Clin Nutr. (2018) 107:932–40. doi: 10.1093/ajcn/nqy052

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Diana A, Haszard JJ, Purnamasari DM, Nurulazmi I, Luftimas DE, Rahmania S, et al. Iron, zinc, Vitamin A and selenium status in a cohort of Indonesian infants after adjusting for inflammation using several different approaches. Br J Nutr. (2017) 118:830–9. doi: 10.1017/S0007114517002860

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Likoswe BH, Phiri FP, Broadley MR, Joy EJM, Patson N, Maleta KM, et al. Inflammation adjustment by two methods decreases the estimated prevalence of zinc deficiency in Malawi. Nutrients. (2020) 12:1–16. doi: 10.3390/nu12061563

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Mwangi MN, Echoka E, Knijff M, Kaduka L, Werema BG, Kinya FM, et al. Iron status of Kenyan pregnant women after adjusting for inflammation using BRINDA regression analysis and other correction methods. Nutrients. (2019) 11:420. doi: 10.3390/nu11020420

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lowe NM, Fekete K, Desci T. Methods of assessment of zinc status in humans: a systematic review. Am J Clin Nutr. (2009) 89:2040s−51s. doi: 10.3945/ajcn.2009.27230G

PubMed Abstract | CrossRef Full Text | Google Scholar

12. World Health Organisation/UNICEF/ICCIDD. Assessment of Iodine Deficiency Disorders and Monitoring Their Elimination Third edition A Guide for Programme Managers. 3rd ed. Geneva: World Health Organisation/UNICEF/ICCIDD (2007).

Google Scholar

13. Phiri FP, Ander EL, Lark RM, Bailey EH, Chilima B, Gondwe J, et al. Urine selenium concentration is a useful biomarker for assessing population level selenium status. Environ Int. (2020) 134:105218. doi: 10.1016/j.envint.2019.105218

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Brown KH, Rivera JA, Bhutta Z, Gibson RS, King JC, Ruel MT, et al. International zinc nutrition consultative group(IZiNCG) technical document #1. Assessment of the risk of zinc deficiency in populations and options for its control. Food Nutr Bull. (2004) 25:s91–s203.

PubMed Abstract | Google Scholar

15. European Food Safety Authority. Scientific opinion on dietary reference values for zinc. EFSA J. (2014) 12:3844. doi: 10.2903/j.efsa.2014.3844

CrossRef Full Text | Google Scholar

16. Raghavan R, Ashour FS, Bailey R. A review of cutoffs for nutritional biomarkers. Adv Nutr. (2016) 7:112–20. doi: 10.3945/an.115.009951

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gamble MV, Ramakrishnan R, Palafox NA, Briand K, Berglund L, Blaner WS. Retinol binding protein as a surrogate measure for serum retinol: studies in vitamin A-deficient children from the Republic of the Marshall Islands. Am J Clin Nutr. (2001) 73:594–601. doi: 10.1093/ajcn/73.3.594

PubMed Abstract | CrossRef Full Text | Google Scholar

18. National Statistical Office (NSO) Community Health Sciences Unit (CHSU) Centers for Disease Control Prevention (CDC) Emory University. Malawi Micronutrient Survey 2015-16. Atlanta, GA: NSO,CHSU,CDC and Emory University (2017). Available online at: https://dhsprogram.com/pubs/pdf/FR319/FR319.m.final.pdf (accessed June 21, 2022).

Google Scholar

19. Wessells KR, Brown KH, Arnold CD, Barffour MA, Hinnouho G-M, Killilea DW, et al. Plasma and nail zinc concentrations, but not hair zinc, respond positively to two different forms of preventive zinc supplementation in young laotian children: a randomized controlled trial. Biol Trace Elem Res. (2021) 199:442. doi: 10.1007/s12011-020-02163-2

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zyba SJ, Shenvi S V, Killilea DW, Holland TC, Kim E, Moy A, et al. A Moderate Increase in Dietary Zinc Reduces DNA Strand Breaks in Leukocytes and Alters Plasma Proteins Without Changing Plasma Zinc. American Journal of Clinical Nutrition (2017). doi: 10.3945/ajcn.116.135327

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Phiri FP, Ander EL, Lark RM, Joy EJM, Kalimbira AA, Suchdev PS, et al. Spatial analysis of urine zinc (Zn) concentration for women of reproductive age and school age children in Malawi. Environ Geochem Health. (2021) 43:259–71. doi: 10.1007/s10653-020-00700-5

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Barbanel CS, Winkelman JW, Fischer GA, King AJ. Confirmation of the department of transportation criteria for a substituted urine specimen. J Occup Environ Med. (2002) 44:407–16. doi: 10.1097/00043764-200205000-00009

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Suchdev PS, Namaste SML, Aaron GJ, Raiten DJ, Brown KH, Flores-ayala R. Overview of the biomarkers reflecting inflammation and nutritional determinants of anemia (BRINDA) Project. Adv Nutr. (2016) 7:349–56. doi: 10.3945/an.115.010215

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pinheiro JC, Bates DM. Mixed-Effects Models in S and S-PLUS. New York: Springer-Verlag (2000).

Google Scholar

25. Panel E, Nda A. Scientific Opinion on Dietary Reference Values for Zinc 1. Parma: Efsajournal (2014). p. 1–74. doi: 10.2903/j.efsa.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

26. van Buuren S, Schönbeck Y, van Dommelen P. Collection, collation and analysis of data in relation to reference heights and reference weights for female and male children and adolescents (0–18 years) in the EU, as well as in relation to the age of onset of puberty and the age at which different stages of puberty are reached in adolescents in the EU. EFSA Support Publ. (2017) 9:59. doi: 10.2903/sp.efsa.2012.EN-255

CrossRef Full Text | Google Scholar

27. Jackson MJ, Jones DA, Edwards RHT, Swainbank IG, Coleman ML. Zinc homeostasis in man: studies using a new stable isotope-dilution technique. Br J Nutr. (1984) 51:199. doi: 10.1079/BJN19840024

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Turnlund JR, King JC, Keyes WR, Gong B, Michel MC. A stable isotope study of zinc absorption in young men: effects of phytate and α-cellulose. Am J Clin Nutr. (1984) 40:1071–7. doi: 10.1093/ajcn/40.5.1071

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lowe NM, Shames DM, Woodhouse LR, Matel JS, Roehl R, Saccomani MP, et al. compartmental model of zinc metabolism in healthy women using oral and intravenous stable isotope tracers. Am J Clin Nutr. (1997) 65:1810–9. doi: 10.1093/ajcn/65.6.1810

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Miller LV, Krebs NF, Hambidge KM. Development of a compartmental model of human zinc metabolism: identifiability and multiple studies analyses. Am J Physiol Regul Integr Comp Physiol. (2000) 279:R1671–84. doi: 10.1152/ajpregu.2000.279.5.R1671

PubMed Abstract | CrossRef Full Text | Google Scholar

31. King JC, Shames DM, Lowe NM, Woodhouse LR, Sutherland B, Abrams SA, et al. Effect of acute zinc depletion on zinc homeostasis and plasma zinc kinetics in men. Am J Clin Nutr. (2001) 74:116–24. doi: 10.1093/ajcn/74.1.116

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Pinna K, Woodhouse LR, Sutherland B, Shames DM, King JC. Exchangeable zinc pool masses and turnover are maintained in healthy men with low zinc intakes. J Nutr. (2001) 131:2288–94. doi: 10.1093/jn/131.9.2288

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Gunter E, Turner W, Neese J, Bayse D. Laboratory Procedures used by the Clinical Division, Centers for Disease Control for the Second Health and Nutrition Examination Survey (NHANES II) 1976–1980. Atlanta: US Department of Health and Human Services (1985).

Google Scholar

34. Schuhmacher M, Domingo JL, Corbella J. Zinc and copper levels in serum and urine: relationship to biological, habitual and environmental factors. Sci Total Environ. (1994) 148:67–72. doi: 10.1016/0048-9697(94)90376-X

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Giroux EL, Durieux M, Schechter PJ. A study of zinc distribution in human serum. Bioinorg Chem. (1976) 5:211–8. doi: 10.1016/S0006-3061(00)82019-0

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Springate SD. The effect of sample size and bias on the reliability of estimates of error: a comparative study of Dahlberg's formula. Eur J Orthod. (2012) 34:158–63. doi: 10.1093/ejo/cjr010

PubMed Abstract | CrossRef Full Text | Google Scholar

37. International Zinc Nutrition Consultative Group. Collecting blood in the field for assessment of plasma or serum zinc concentration. IZiNCG Pract Tips. (2019) 1–16.

Google Scholar

38. Melichar B, Jandik P, Tichy M, Malir F, Mergancova J, Voboril Z. Urinary zinc excretion and acute phase response in cancer patients. Clin Investig. (1994) 72:1012–4. doi: 10.1007/BF00577746

PubMed Abstract | CrossRef Full Text | Google Scholar

39. McDonald CM, Suchdev PS, Krebs NF, Hess SY, Wessells KR, Ismaily S, et al. Adjusting plasma or serum zinc concentrations for inflammation: biomarkers Reflecting Inflammation and Nutritional Determinants of Anemia (BRINDA) project. Am J Clin Nutr. (2020) 111:927–37. doi: 10.1093/ajcn/nqz304

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Olechnowicz J, Tinkov A, Skalny A, Suliburska J. Zinc status is associated with inflammation, oxidative stress, lipid, and glucose metabolism. J Physiol Sci. (2017) 68:19–31. doi: 10.1007/s12576-017-0571-7

PubMed Abstract | CrossRef Full Text | Google Scholar

41. World Health Organization (WHO) Food and Agriculture Organization (FAO) International Atomic Energy Association (IAEA). Trace Elements in Human Health and Nutrition. Geneva: World Health Organisation (1996).

Google Scholar

42. Institute of Medicine/Food and Nutrition Board. Dietary Reference Intakes for Vitamin A, Vitamin K, Arsenic, Boron, Chromium, Copper, Iodine, Iron, Manganese, Molybdenum, Nickel, Silicon, Vanadium, and Zinc. Washington DC: Institute of Medicine, National Academy Press (2000).

Google Scholar