|
|
|
|
Section of Nutrition, Department of Pediatrics, University of Colorado Health Sciences Center, Denver, CO 80262
* To whom correspondence should be addressed. E-mail: leland.miller{at}uchsc.edu.
| ABSTRACT |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
We report here the development of a new mathematical model of the quantity of zinc absorbed per day as a function of dietary zinc and phytate. There are several characteristics of this model and its development that, taken together, make it novel. First, the model was derived from a biochemical conception of the absorption process. A mathematical model having a foundation in the physical process being modeled is potentially more informative and more likely to be valid than a model utilizing a mathematical function that merely mimics the behavior of the data. Furthermore, given its foundation in the physical process, the model is amenable to revision based on additional knowledge of that process. Second, absorbed zinc was modeled as a function of both dietary zinc and phytate. Third, the variables modeled were the total daily dietary zinc (TDZ) and phytate (TDP) ingested and total daily absorbed zinc (TAZ) from all the meals over an entire day. The measurement of absorption over a full day is understood to be more meaningful than single meal measurements. And the modeling of quantities ingested and absorbed is likely to be advantageous to modeling ratios of those quantities, i.e. fractional absorption or phytate:zinc molar ratio. Fourth, nonlinear regression analysis was used to fit the model to data. Because zinc absorption is primarily a saturable, carrier-mediated process (1517), the relation of absorbed to dietary zinc is nonlinear. The effect of phytate on absorption is also expected to be nonlinear. Analyzing such relations with nonlinear regression is generally preferred to using data transformation and linear regression, unless the transformation is necessitated for other reasons. Fifth, because validation is an important component of model development, the model was submitted to an assessment of its validity based primarily on its fit to selected data from the literature.
The development of this model is a continuation of our application of saturation response modeling, derived from pharmacodynamic data analysis, to zinc absorption data (9,11,13,14).
| Materials and Methods |
|---|
|
|
|---|
|
![]() | (Eq. 1) |
![]() | (Eq. 2) |
where the P represents phytate, Zn represents zinc, R represents transport receptors, and ZnP and ZnR are the zinc-phytate and zinc-receptor complexes, respectively. The u subscript indicates the free (unbound) form of the species. It is assumed that the transport receptors have a single zinc binding site and that each phytate molecule binds only 1 zinc ion. Although phytate is capable of binding multiple divalent cations, the chemical environment of the intestine, notably the pH and presence of other cations and ligands, precludes the possibility of all the binding sites being occupied by zinc ions. In fact, it is speculated that, practically, phytate binds a single zinc ion in vivo (Boyd O'Dell, University of Missouri, personal communication). It is conceivable that, as the model evolves and is applied to additional data, it may be modified to accommodate binding of multiple zinc ions; however, as the binding constant is different for each additional ion bound, model complexity will become an issue. There are no assumptions regarding the relative concentrations of any of the species involved in the binding reactions. After the substitution of the following equalities and approximation
![]() | (Eq. 3) |
![]() | (Eq. 4) |
![]() | (Eq. 5) |
![]() | (Eq. 6) |
![]() | (Eq. 7) |
Equations 1 and 2 are rearranged to
![]() | (Eq. 8) |
![]() | (Eq. 9) |
The subscript t indicates the total quantity of the species. KP and KR are the equilibrium dissociation constants of zinc-phytate and zinc-receptor binding, respectively.
Equation 3 introduces an approximation, that Pu
Pt, which is necessary to simplify the model so that it is of practical use. Without this approximation, it is anticipated that an exact solution of the model would be so complex as to be of limited practical value. The approximation is good to the extent that Pt > >ZnP. Once the model has been fit to the data and the parameters estimated, the closeness of the approximation can be assessed using this relation derived from the equations above:
![]() | (Eq. 10) |
Furthermore, with this information, an evaluation of the sensitivity of the model to the approximation can be made.
Equation 8 is solved for ZnP and substituted into Equation 9, which is then solved for ZnR.
Finally, because all species occupy the same luminal volume, units of concentration are changed to molar rates and the equivalent experimental variables and a model parameter are substituted for the biochemical terms as follows:
![]() |
![]() |
![]() |
![]() |
All are in units of millimoles per day. The resulting equation for the model is:
![]() | (Eq. 11) |
This is a trivariate model having 2 independent (predictor) variables, TDZ and TDP, with TAZ being the dependent (response) variable. The model has 3 parameters: AMAX, KR, and KP.
The model equation may be modified to use fractional absorption of zinc (FAZ) instead of TAZ, or to use phytate:zinc molar ratios (RPZ) instead of TDP to make predictions after the parameters have been estimated:
![]() | (Eq. 12) |
![]() | (Eq. 13) |
A version using both FAZ and RPZ may be also created.
It is also possible without much difficulty to augment the model to include a passive absorption mechanism, observed in animals and thought to possibly play a minor role in zinc absorption in humans (15,16). On the assumption that passive absorption is a nonmediated diffusion process and proportional to both the quantity of unbound zinc in the intestinal lumen and a passive nonmediated absorption coefficient (NMA), the model becomes:
![]() | (Eq. 14) |
With NMA, this model has 4 parameters. It can be expected that the small number of available data are inadequate to currently support the additional parameter.
Selection of data. The 15 data points compiled by the IZiNCG from 9 published studies (2129) provided the majority of data for our analyses. Their data selection criteria were: 1) radio- or stable-isotope studies that estimated true zinc absorption from total (daily) diets by correcting for intestinal losses of endogenous zinc; 2) studies of typical mixed, refined vegetarian, or unrefined, cereal-based diets, but not those that used semipurified formula diets or diets with exogenous zinc salts added; and 3) studies with male or female adults, with no geographic limitations (1). In addition to these data, we included 6 data points from 2 other studies (8,30), including 1 that used a liquid formula diet. These latter data did not meet the 2nd criterion of the IZiNCG; otherwise, the additional data were in agreement with the criteria. Fortunately, these additional data extended the range of TDP beyond that of the IZiNCG data.
Each of the 21 data were a mean of the results from 4 to 21 adult subjects. Mean data were used, because the individual data from many of the studies were not available. All subjects were apparently healthy and assumed to be in normal zinc status. No data from studies of children were included. A total of 105 subjects participated in the 11 studies. TAZ was always measured using established isotope techniques, as specified in the first selection criterion. TDZ was determined by chemical analysis of representative diets in most studies, although at least 1 study relied on the calculation of dietary zinc from diet records. In the case of the liquid formula diet, a measured quantity of zinc was added during diet composition. In a majority of the studies, TDP was calculated from diet composition by the authors of the studies or by the IZiNCG, but phytate was analyzed in 2 of the studies and a known quantity was added to the liquid formula diet. In most studies, the subjects were on the test diets for at least 7 d prior to absorption measurements. Dietary calcium and protein data were also available for 14 of the data.
Data analysis. The model was fit to the data using nonlinear regression analysis with the programs DataFit (version 8.1, Oakdale Engineering) and SigmaPlot (version 10.0, Systat Software). There was no weighting of data for the regression analyses. DataFit provided regression statistics and 3-dimensional graphing of the data and model. SigmaPlot calculated additional regression statistics. Other statistical analyses and 2-dimensional graphing were performed with GraphPad Prism (version 4.03, GraphPad Software). The regression assumptions of normality, constant variance, and independence (lack of correlation) of the residuals were tested as follows. Distribution of the residuals was evaluated with the Kolmogorov-Smirnov and Shapiro-Wilk normality tests. Constant variance was tested by computing Spearman rank correlation between the absolute values of the residuals and the TAZ values. Independence was determined using the runs test and the Durbin-Watson test for serial correlation. Statistical significance was defined to be at the level of 0.05.
| Results |
|---|
|
|
|---|
|
![]() |
![]() |
![]() |
The P(t) values, the probability of being wrong in concluding that the parameter value is not 0, of the estimates are 0.011, 0.24, and 0.12 for AMAX, KR, and KP, respectively. The R2 is 0.82 and the adjusted version (
), taking into account the degrees of freedom, is 0.79. The sum of squares of the residuals is 0.0011. One datum deviated noticeably from the model and was flagged as a possible outlier based on the magnitude of its standardized, and Studentized, residual. Although it was not removed for fitting the model, evaluations of the model residuals were performed with and without the discrepant point.
Examination of the residual normal probability plot with (Fig. 3) and without the discrepant datum indicated a normal distribution, although inclusion of the questionable point (upper right) did skew the plot noticeably. The Kolmogorov-Smirnov and Shapiro-Wilk normality tests showed the residuals to have a normal distribution with P-values of 0.078 and 0.052 (>0.10 and 0.73 without the discrepant datum), respectively. The residuals passed a test for constant variance using Spearman rank correlation (P = 0.68). Runs tests were applied to the residuals relative to TDZ, TDP, and predicted TAZ, giving P-values of 0.25, 0.61, and 0.61, respectively, demonstrating independence, i.e. that the number of runs was similar to that expected for randomly distributed data in each case. A Durbin Watson test value of 1.78 also indicated a lack of serial correlation, although this test usually applied to data that were equally spaced relative to the predictor variables. In addition to the formal tests, examination of the graphs of the residuals relative to TDZ, TDP, and predicted TAZ (Fig. 4AC) confirmed the findings of constant variance and independence. The only hint of nonconstant variance occurs in Figure 4B and is probably attributable to the clustering of the data in the low TDP range. And, the only discernable nonrandomness in the residual plots is that the residuals are negative at high x-axis values. This is attributable to the discrepant point discussed above and is no longer evident when the point is removed. In summary, all the analyses indicated that the residuals met the assumptions of normality, constant variance, and independence.
|
|
|
Regarding the approximation of Equation 3 and its evaluation with Equation 10, the ratio of Pt to ZnP was found to be 28 for the mean TAZ from the data and having a range of 8.1 to 108 for the range of TAZ, with the ratio being the lowest when TAZ is high. Under the worst case conditions, the effect on the accuracy of the model was estimated to be
5%, i.e. at the high end of the TAZ range values produced by the model are 5% lower than without the approximation.
To investigate the possible role of calcium, sometimes quantified as the calcium*phytate/Zn molar ratio, or protein in zinc absorption in the presence of phytate (31), the dietary calcium, calcium*phytate/Zn molar ratio, and dietary protein data from 14 studies were plotted against the residuals (Supplemental Fig. 2). None of these data exhibited a relation with the residuals, suggesting these dietary factors did not affect zinc absorption in these studies. This is consistent with the IZiNCG findings (1).
| Discussion |
|---|
|
|
|---|
Figure 2 conveys essential information about the behavior of the model and its fit to the data. Although much of this information is presented here in other forms, the 3-dimensional representation best illustrates the distribution of the data in the predictor variable plane and the random nature of the residuals relative to all the variables. It is obvious from this graph that the model does not support the existence of a threshold for a phytate effect on absorption.
The goodness of fit as reflected in the R2 of 0.82 is very supportive of the model, indicating that 82% of the variance in TAZ is explained by the model. Analyses of the residuals also confirmed the quality of the fit and adequacy of the model while testing the assumptions underlying the regression analysis. These assumptions, which must be met for the regression to be valid, are that the errors, as reflected in the residuals, exhibit a Gaussian (normal) distribution, constant variance, and independence. Evaluating independence is necessary to detect nonrandom behavior in the residuals, indicating a systematic deviation of the model from the data. The residuals passed all the tests and visual examinations, although the presence of a datum that deviated noticeably from the model caused some results to be marginal. This point and its large residual are evident in Figures 24. Although this datum looks suspiciously like an outlier, we have examined the study in which it was measured and found no reason to reject it as an outlier at this time. Therefore, the point was retained for fitting the model and estimating the parameters. However, because of its disproportionate effect on the analyses of the residuals, the analyses without the discrepant point are thought to more accurately reflect the error characteristics on the whole.
The uncertainties around the parameter estimates are large, although AMAX may be considered to be adequately determined, because we can reasonably conclude that it is not 0 (P = 0.011). Though large uncertainty in the parameters is usually the result of a model having too many parameters (overparameterization) and collinearity of the parameters, it is not justified, or prudent, to simplify the model for the following reasons: 1) The model is supported by the other measures of validity; 2) the model is already of minimal complexity, given that there are 2 uncorrelated (r = 0.09) predictor variables and the relations are nonlinear and nonarbitrary; 3) each of the parameters has a well-defined relation to an essential element of the absorption process and cannot be eliminated without compromising the foundational validity of the model; and 4) the current data are limited in number and exhibit variability that may related to the multiple laboratories, analytical methods, and experimental protocols involved. It is expected that the quality of the estimates will improve notably with additional data.
Regarding the values of the parameter estimates, we have found virtually no existing data with which to compare them, although it is noteworthy that the relative magnitude of the equilibrium dissociation constants, i.e. KP > KR, is consistent with that reported by Wing et al. (6) and with the deduction by O'Dell, from an experiment in which dietary EDTA alleviated the detrimental effect of phytate, that EDTA has a higher binding affinity for zinc than phytate and that mucosal receptors have a still higher affinity (Boyd O'Dell, University of Missouri, personal communication). Our previous application of a basic saturable response model (11) to the data used by the Food and Nutrition Board (7) estimated a value of 0.11 mmol/d for AMAX, not significantly different (P = 0.71) from the 0.13 mmol/d estimated here. It should be noted that, although we expect that this model may be appropriate for application to data from any human population, the parameter values derived here characterize absorption in healthy adults with assumed normal zinc status and, therefore, caution should be exercised in interpreting their predictive application to other populations until additional data have been modeled. As well as the predictive uses of the model that the parameter estimates permit, it is possible that the parameter values themselves will provide information about absorption. Because AMAX is related to the number of transport receptors and KR and KP quantify the association/dissociation characteristics of the binding reactions, application of the model may contribute to our limited knowledge of receptor regulation and binding chemistry as it relates to bioavailability.
The opportunities for comparing the model's predictions to other models are very limited. Not surprisingly, the TAZ predictions from our model agree with those of the IZiNCG model (1) with ±6% across the range of the data, but the models diverge at higher TDZ values. The TDZ at which divergence occurs varies inversely with the RPZ. Above this, our model predicts that TAZ increases at a lower rate with increasing TDZ, which may be attributable to the additional high phytate data that we have used. Whereas more detailed comparison with modeling of the Food and Nutrition Board data (7) was not possible because phytate data were not available, we did use the model to predict the phytate intake of the those subjects. The model predicted a mean TDP of 0.45 mmol/d and a resulting RPZ of 3.1. This is consistent with the fact that the Food and Nutrition Board data were selected from studies of very low phytate diets.
We have shown several variations of the model to accommodate the use of fractional absorption, FAZ, or the phytate:zinc molar ratio, RPZ. Although they add flexibility for the model's predictive applications, these forms of the model are not as well suited for fitting to data, because FAZ and RPZ are both ratios of the more fundamental variables, i.e. FAZ = TAZ/TDZ and RPZ = TDP/TDZ. In the case of FAZ, because there is a correlation between TDZ and TAZ (r = 0.73, P = 0.0002), the relation between FAZ and TDZ is not as strong as that between TAZ and TDZ. This is manifest in our use of Equation 12 to fit the data. Although producing similar parameter estimates, the goodness if fit was notably inferior. The recommendation to relate zinc absorption to dietary component quantities rather than ratios has been made previously, although based on somewhat different considerations (4). The use of Equation 14, with the passive absorption parameter, to analyze these limited data is a good example of overparameterization, as evidenced by the very wide CI for all parameters. It is noteworthy that the goodness of fit was not improved with the added parameter, apparently indicating that these data do not exhibit evidence of passive absorption. Additional data will be required to discern the existence of passive absorption and demonstrate the usefulness of this version of the model.
In conclusion, we have developed a mathematical model of zinc absorption from a basic conception of the relevant intestinal biochemistry and fit it to selected existing data. Evaluation of the fit finds it to be good and in compliance with regression assumptions, thereby supporting the validity of the model. Evaluation of the parameter estimates and model predictions are also supportive of the model's validity. We judge the model to be well founded, with immediate relevance and applicability to the study of zinc nutrition and metabolism and the estimation of dietary zinc requirements in varied populations. Furthermore, we anticipate the model's evolution and improvement as new data are analyzed and further knowledge of the absorption process incorporated.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
2 Supplemental Figures 1 and 2 are available with the online posting of this paper at jn.nutrition.org. ![]()
3 Abbreviations used: AMAX, maximum absorption; FAZ, fractional absorption of zinc; IZiNCG, International Zinc Nutrition Consultative Group; KP, equilibrium dissociation constant of zinc-phytate binding reaction; KP1, association rate constant of zinc-phytate binding reaction; KP2, dissociation rate constant of zinc-phytate binding reaction; KR, equilibrium dissociation constant of zinc-receptor binding reaction; KR1, association rate constant of zinc-receptor binding reaction; NMA, passive (nonmediated) absorption coefficient; RPZ, phytate:zinc molar ratio; TAZ, total daily absorbed zinc; TDP, total daily dietary phytate; TDZ, total daily dietary zinc; ZnP, zinc-phytate complex; ZnR, zinc-receptor complex. ![]()
Manuscript received 7 September 2006. Initial review completed 3 October 2006. Revision accepted 5 November 2006.
| LITERATURE CITED |
|---|
|
|
|---|
1. International Zinc Nutrition Consultative Group (IZiNCG). Assessment of the risk of zinc deficiency in populations and options for its control. Hotz C and Brown KH, editors. Food and Nutrition Bulletin 25. Tokyo: United Nations University Press; 2004. p. S91204.
2. Lönnerdal B, Bell JG, Hendrickx AG, Burns RA, Keen CL. Effect of phytate removal on zinc absorption from soy formula. Am J Clin Nutr. 1988;48:13016.
3. Sandström B, Almgren A, Kivistö B, Cederblad Å. Effect of protein level and protein source on zinc absorption in humans. J Nutr. 1989;119:4853.
4. Wise A. Phytate and zinc bioavailability. Int J Food Sci Nutr. 1995;46:5363.[Medline]
5. WHO, Food and Agriculture Organization (FAO), International Atomic Energy Association (IAEA). Trace elements in human nutrition and health. Geneva: WHO; 1996.
6. Wing KR, Wing AM, Sjöström R, Lönnerdal B. Efficacy of a Michaelis-Menton model for the availability of zinc, iron and cadmium from an infant formula diet containing phytate. In: Fischer PWF, L'Abbé MR, Cockell KA, Gibson RS, editors. Trace elements in man and animals 9: Proceedings of the Ninth International Symposium on Tracer Elements in Man and Animals; 1996 May 1924; Banff, Alberta, Canada. Ottawa: NRC Research Press; 1997. p. 312.
7. Food and Nutrition Board, Institute of Medicine. Dietary reference intakes for vitamin A, vitamin K, boron, chromium, copper, iodine, iron, manganese, molybdenum, nickel, silicon, vanadium and zinc. Washington: National Academy Press; 2001.
8. Hambidge KM, Huffer JW, Raboy V, Grunwald GK, Westcott JL, Sian L, Miller LV, Dorsch JA, Krebs NF. Zinc absorption from low-phytate hybrids of maize and their wild-type isohybrids. Am J Clin Nutr. 2004;79:10539.
9. Tran CD, Miller LV, Krebs NF, Lei S, Hambidge KM. Zinc absorption as a function of the dose of zinc sulfate in aqueous solution. Am J Clin Nutr. 2004;80:15703.
10. Chiplonkar SA, Agte VV. Predicting bioavailable zinc from lower phytate forms, folic acid and their interactions with zinc in vegetarian meals. J Am Coll Nutr. 2005;25:2633.
11. Hambidge KM, Miller LV, Tran CD, Krebs NF. Measurements of zinc absorption: application and interpretation in research designed to improve human zinc nutriture. Int J Vitam Nutr Res. 2005;75:38593.[Medline]
12. Fredlund K, Isaksson M, Rossander-Hulthén L, Almgren A, Sandberg A-S. Absorption of zinc and retention of calcium: dose-dependent inhibition by phytate. J Trace Elem Med Biol. 2006;20:4957.[Medline]
13. Hambidge KM, Krebs NF, Westcott JE, Miller LV. Changes in zinc absorption during development. J Pediatr. 2006;149 Suppl 3:S648.
14. Hambidge KM, Abebe Y, Gibson RS, Westcott JE, Miller LV, Sian L, Stoecker BJ, Arbide I, Teshome A, et al. Zinc absorption during late pregnancy in rural southern Ethiopia. Am J Clin Nutr. 2006;84:11026.
15. Steel L, Cousins RJ. Kinetics of zinc absorption by luminally and vascularly perfused rat intestine. Am J Physiol. 1985;248:G4653.
16. Reyes JG. Zinc transport in mammalian cells. Am J Physiol. 1996;270:C40110.
17. Liuzzi JP, Cousins RJ. Mammalian zinc transporters. Annu Rev Nutr. 2004;24:15172.[Medline]
18. Foster DM, Aamodt RL, Henkin RI, Berman M. Zinc metabolism in humans: a kinetic model. Am J Physiol. 1979;237:R3409.
19. Gabrielsson J, Weiner D. Pharmacokinetic and pharmacodynamic data analysis: concepts and applications. 2nd ed. Stockholm: Swedish Pharmaceutical Press; 1997.
20. Kenakin T. Pharmacologic analysis of drug-receptor interaction. 3rd ed. Philadelphia: Lippincott-Raven Publishers; 1997.
21. Wada L, Turnlund JR, King JC. Zinc utilization in young men fed adequate and low zinc intakes. J Nutr. 1985;115:134554.
22. Hunt JR, Mullen LK, Lykken GI. Zinc retention from an experimental diet based on the US FDA Total Diet Study. Nutr Res. 1992;12:133544.
23. Hunt JR, Gallagher SK, Johnson LK, Lykken GI. High- versus low-meat diets: effects on zinc absorption, iron status, and calcium, copper, iron, magnesium, manganese, nitrogen, phosphorus, and zinc balance in postmenopausal women. Am J Clin Nutr. 1995;62:62132.
24. Knudsen E, Sandström B, Solgaard P. Zinc, copper and magnesium absorption from a fibre-rich diet. J Trace Elem Med Biol. 1996;10:6876.[Medline]
25. Sian L, Mingyan X, Miller LV, Tong L, Krebs NF, Hambidge KM. Zinc absorption and intestinal losses of endogenous zinc in young Chinese women with marginal zinc intakes. Am J Clin Nutr. 1996;63:34853.
26. Lowe NM, Shames DM, Woodhouse LR, Matel JS, Roehl R, Saccomani MP, Toffolo G, Cobelli C, King JC. A compartmental model of zinc metabolism in healthy women using oral and intravenous stable isotope tracers. Am J Clin Nutr. 1997;65:18109.
27. Hunt JR, Matthys LA, Johnson LK. Zinc absorption, mineral balance, and blood lipids in women consuming controlled lactoovovegetarian and omnivorous diets for 8 wk. Am J Clin Nutr. 1998;67:42130.[Abstract]
28. Pinna K. Effect of a low zinc intake on zinc homeostasis and immune function in adult men [doctoral dissertation]. Berkeley (CA): University of California; 1999.
29. Adams CL, Hambidge M, Raboy V, Dorsch JA, Sian L, Westcott JL, Krebs NF. Zinc absorption from a low-phytic acid maize. Am J Clin Nutr. 2002;76:5569.
30. 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:10717.
31. Lönnerdal B. Dietary factors influencing zinc absorption. J Nutr. 2000;130:S137883.
32. Carson ER, Cobelli C, Finkelstein L. The mathematical modeling of metabolic and endocrine systems. New York: John Wiley & Sons; 1983.
33. Bonate PL. Pharmacokinetic-pharmacodynamic modeling and simulation. New York: Springer Science and Business Media; 2006.
This article has been cited by other articles:
![]() |
J. R Hunt and J. M Beiseigel Dietary calcium does not exacerbate phytate inhibition of zinc absorption by women from conventional diets Am. J. Clinical Nutrition, March 1, 2009; 89(3): 839 - 843. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Hambidge, L. V. Miller, J. E. Westcott, and N. F. Krebs Dietary Reference Intakes for Zinc May Require Adjustment for Phytate Intake Based upon Model Predictions J. Nutr., December 1, 2008; 138(12): 2363 - 2366. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R Hunt, J. M Beiseigel, and L. K Johnson Adaptation in human zinc absorption as influenced by dietary zinc and bioavailability Am. J. Clinical Nutrition, May 1, 2008; 87(5): 1336 - 1345. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||