Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-04T21:54:38.852Z Has data issue: false hasContentIssue false

Two non-consecutive 24 h recalls using EPIC-Soft software are sufficiently valid for comparing protein and potassium intake between five European centres – results from the European Food Consumption Validation (EFCOVAL) study

Published online by Cambridge University Press:  28 September 2010

Sandra P. Crispim*
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
Jeanne H. M. de Vries
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
Anouk Geelen
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
Olga W. Souverein
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
Paul J. M. Hulshof
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
Lionel Lafay
Affiliation:
French Agency for Food, Environmental and Occupational Health Safety (ANSES)/Food Safety Department/Food Intake-Nutritional Epidemiology Unit, 27-31 av Général Leclercq, 94701Maisons-Alfort, France
Anne-Sophie Rousseau
Affiliation:
Université de Nice Sophia Antipolis, UMR 907, Nice, F-06002, France
Inger T. L. Lillegaard
Affiliation:
Department of Nutrition, Institute of Basic Medical Sciences, University of Oslo, PO Box 1046, 0316Oslo, Norway
Lene F. Andersen
Affiliation:
Department of Nutrition, Institute of Basic Medical Sciences, University of Oslo, PO Box 1046, 0316Oslo, Norway
Inge Huybrechts
Affiliation:
Department of Public Health, Faculty of Medicine and Health Sciences, Ghent University, UZ-2 Blok A, De Pintelaan 185, B-9000Ghent, Belgium International Agency for Research on Cancer (IARC), 150 Cours Albert Thomas, 69372Lyon Cedex 08, France
Willem De Keyzer
Affiliation:
Department of Public Health, Faculty of Medicine and Health Sciences, Ghent University, UZ-2 Blok A, De Pintelaan 185, B-9000Ghent, Belgium Department of Nutrition and Dietetics, University College Ghent, Keramiekstraat 80, B-9000Ghent, Belgium
Jiri Ruprich
Affiliation:
National Institute of Public Health, Department for Food Safety and Nutrition, Palackeho 1–3, 612 42Brno, Czech Republic
Marcela Dofkova
Affiliation:
National Institute of Public Health, Department for Food Safety and Nutrition, Palackeho 1–3, 612 42Brno, Czech Republic
Marga C. Ocke
Affiliation:
National Institute for Public Health and the Environment (RIVM), PO Box 1Bilthoven3720 BA, The Netherlands
Evelien de Boer
Affiliation:
National Institute for Public Health and the Environment (RIVM), PO Box 1Bilthoven3720 BA, The Netherlands
Nadia Slimani
Affiliation:
International Agency for Research on Cancer (IARC), 150 Cours Albert Thomas, 69372Lyon Cedex 08, France
Pieter van't Veer
Affiliation:
Division of Human Nutrition, Wageningen University, Bomenweg 2, Wageningen6703 HD, The Netherlands
*
*Corresponding author: S. P. Crispim, fax +31 0317 482782, email [email protected]; [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The use of two non-consecutive 24 h recalls using EPIC-Soft for standardised dietary monitoring in European countries has previously been proposed in the European Food Consumption Survey Method consortium. Whether this methodology is sufficiently valid to assess nutrient intake in a comparable way, among populations with different food patterns in Europe, is the subject of study in the European Food Consumption Validation consortium. The objective of the study was to compare the validity of usual protein and K intake estimated from two non-consecutive standardised 24 h recalls using EPIC-Soft between five selected centres in Europe. A total of 600 adults, aged 45–65 years, were recruited in Belgium, the Czech Republic, France, The Netherlands and Norway. From each participant, two 24 h recalls and two 24 h urines were collected. The mean and distribution of usual protein and K intake, as well as the ranking of intake, were compared with protein and K excretions within and between centres. Underestimation of protein (range 2–13 %) and K (range 4–17 %) intake was seen in all centres, except in the Czech Republic. We found a fair agreement between prevalences estimated based on the intake and excretion data at the lower end of the usual intake distribution ( < 10 % difference), but larger differences at other points. Protein and K intake was moderately correlated with excretion within the centres (ranges = 0·39–0·67 and 0·37–0·69, respectively). These were comparable across centres. In conclusion, two standardised 24 h recalls (EPIC-Soft) appear to be sufficiently valid for assessing and comparing the mean and distribution of protein and K intake across five centres in Europe as well as for ranking individuals.

Type
Full Papers
Copyright
Copyright © The Authors 2010

National food consumption surveys aim to provide information on the mean and distribution of food and nutrient intakes of the population and related subgroups, in order to develop and evaluate nutrition policies. In addition, national food consumption surveys are essential to provide data for risk assessment work, as conducted by the European Food Safety Authority(1). In Europe, food consumption data originating from national surveys are not always comparable because they differ in a number of aspects, such as the choice of the dietary assessment method and the reference period of the data collection(Reference Charzewska2Reference Verger, Ireland and Moller4). Furthermore, some countries do not have national food consumption surveys in place(Reference Verger, Ireland and Moller4).

The European Food Consumption Survey Method consortium has acknowledged the need for policy-relevant dietary indicators that are comparable among European countries, which could contribute to the establishment of a Community Health Monitoring System(Reference Brussaard, Johansson and Kearney5). They recommended two non-consecutive days of 24 h recall using EPIC-Soft software (Lyon, Rhone Alpes, France) as the preferred method to assess the dietary intake in future pan-European monitoring surveys in adults. In addition, they specified total fat, SFA and ethanol as the components of most relevance in this assessment(Reference Biro, Hulshof and Ovesen6Reference Sliman and Valsta8).

The 24 h recall is a commonly used dietary assessment method in food consumption surveys in Europe(Reference Verger, Ireland and Moller4) and is also being used in surveys in the USA(Reference Conway, Ingwersen and Moshfegh9), Canada(10), Australia(Reference McLennan and Podger11) and New Zealand(Reference Russell, Parnell and Wilson12). A major advantage of using 24 h recalls in (inter)national surveys is that the method is useful for comparison of heterogeneous populations with different ethnicity and literacy(Reference Biro, Hulshof and Ovesen6). In addition, a computerised version of 24 h recalls seems to be the best means of standardising and controlling for sources of error attributable to 24 h recall interviews(Reference Biro, Hulshof and Ovesen6, Reference Willett13). Nevertheless, computerised 24 h recalls need to be tailor-made to every included country and/or study, e.g. by adaptations of the food and recipe list. Therefore, whether this methodology performs in a comparable way across countries with different food consumption patterns in Europe deserves further exploration, as validity of the 24 h recall depends on both the characteristics of the method and the study population.

Biological markers offer an important opportunity to evaluate the dietary assessment methods since errors are likely to be truly independent between the measurements of biomarker and dietary intake(Reference Ocke and Kaaks14). Urinary N and K are two of the few available recovery biomarkers to assess the nutrient intakes(Reference Bingham and Cummings15, Reference Tasevska, Runswick and Bingham16). With the use of these two biomarkers, a single 24 h recall using EPIC-Soft has been previously validated for assessing the group mean intakes of protein of twelve centres in six countries within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort study(Reference Slimani, Bingham and Runswick17). Yet, the accuracy of this methodology needs to be determined when aiming at estimating usual dietary intake among different European populations by collecting two independent 24 h recalls. Hence, following the path of the European Food Consumption Survey Method (EFCOSUM), the European Food Consumption Validation (EFCOVAL) consortium aimed to further develop and validate a European food consumption method using EPIC-Soft software for assessing the food and nutrient intakes within European countries and for comparisons between them. In the present paper, we aim to compare the validity of usual protein and K intake estimated from two non-consecutive standardised 24 h recalls using EPIC-Soft between five selected centres in Europe. This was done by addressing the bias present in the estimation of each centre's mean and distribution of intake as well as the ranking of individuals within and between centres according to their intake.

Subjects and methods

Subjects

Data were collected in five European countries: Belgium, the Czech Republic, France (Southern part), The Netherlands and Norway. These countries were selected to represent a large variety in food patterns across Europe. Data were collected in the South of France to include the characteristics of the Mediterranean diet. A food pattern from Central/Eastern Europe was represented by the Czech Republic, from the Scandinavian countries by Norway and from the western part of Europe by Belgium and The Netherlands. Another reason for their selection was their experience in performing nutrition monitoring surveys. The present study was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving human subjects were approved by ethical committees in each centre involved in the data collection.

We recruited subjects by convenience sampling through advertisements (newspaper and websites), mailing lists, among others. Recruitment of institutionalised subjects was not allowed, nor included more than one member of a household. Subjects were informed about the study through information meetings at the institutions/universities in the Czech Republic, France and The Netherlands, and by phone, letter and personally in Belgium and Norway. At these occasions, a screening questionnaire was filled in to confirm the subjects' eligibility in the study. Subsequently, the eligible participants gave written informed consent, and appointments for later visits were scheduled. Exclusion criteria were currently taking diuretics, following prescribed dietary therapy, being enrolled in another study in the same period, not being able to read or speak the national language, being pregnant, lactating, having diabetes mellitus or kidney disease and donating blood or plasma during or < 4 weeks before the study. para-Aminobenzoic acid (PABA) was used to check the completeness of urine collections; therefore, subjects hypersensitive to PABA or taking antibiotics containing sulphonamides, which are PABA-antagonistic, were not eligible for the study.

Taking into account an anticipated dropout percentage of 20 % and aiming at a net sample of fifty per stratum, a total of sixty men and sixty women were recruited per centre (n 600). The age range of subjects was 45–65 years, which was chosen to limit the heterogeneity of the sample. Furthermore, we aimed to include at least ten men and ten women in each of the three predetermined categories of education level (low, intermediate and high) per centre. We used country-specific classifications to define each category level.

We excluded one subject because no data for recall and biomarker collected on the same day were available. Therefore, the study population comprised 599 subjects (296 men and 303 women).

Study design

Wageningen University (The Netherlands) was, as the coordinating centre, responsible for the overall logistics of the validation study in the EFCOVAL consortium. For standardisation, all study procedures, i.e. on recruitment and fieldwork conditions, data processing formats, quality-control aspects and specimen collection, storage and transport details, were described in protocols. The recruitment of subjects and data collection in The Netherlands were performed from April to July 2007, 6 months before the other four centres, in order to test all the procedures of the fieldwork beforehand and to be able to refine the protocols. The other centres started the fieldwork in October or November 2007, with the last centre finalising the collection by April 2008.

At the beginning of the study, subjects had their body weight and height measured in the study centres. Then, a 24 h recall and a 24 h urine collection were obtained covering the same reference day. Subjects were aware of the days of data collection but not of the purpose of the interviews. The second recall and urine collection were obtained at least 1 month after the first one.

Dietary data

The two 24 h recalls were collected using two modes of administration: one by phone and one face-to-face at the centre since it is likely that future food consumption monitoring surveys will be conducted in both ways across European countries. The order of the two modes of administration was randomly allocated among the subjects.

Furthermore, the appointments for the dietary recalls followed a randomised schedule, which included all days of the week. This randomisation allowed the same person to have the same recalled weekday for both interviews by chance. Interviewers in each centre were nutritionists or dietitians who were trained in interviewing skills and working with EPIC-Soft in the context of the validation study. They were guided by qualified local trainers who were previously trained by staff from the Wageningen coordination centre and the National Institute for Public Health and the Environment in The Netherlands. Interviewers were aware of the objectives of the study. The centres were allowed to organise their data collection in the same way they would do in a future performance of their nutritional surveillance system. An example is that interviewees were permitted to check food packages and household measures in their home for more detailed information during the phone interview while this was not possible during the face-to-face interview at the study centre. Another example is that dietary recalls in Belgium, the Czech Republic and The Netherlands were not conducted on Sundays. Therefore, Saturday's intake was recalled 2 d later, on Mondays.

The two 24 h recalls were collected using EPIC-Soft (version 9.16). The structure and standardisation procedure of EPIC-Soft have been described elsewhere(Reference Slimani, Deharveng and Charrondiere18, Reference Slimani, Ferrari and Ocke19). Briefly, EPIC-Soft is a computer-assisted 24 h dietary recall that follows standardised steps when describing, quantifying, probing and calculating the food intakes(Reference Slimani, Deharveng and Charrondiere18). All the participating countries had an existing version of EPIC-Soft available, except the Czech Republic for which a new country-specific version was developed. In addition, EPIC-Soft databases were adapted for each centre in terms of some common specifications for the EFCOVAL study (e.g. soups were treated as recipes rather than food items). Furthermore, the centres generated or updated a list of the single food items and recipes expected to be consumed by their participants. Modifications of such lists were needed afterwards based on notes made during the interview. The methods of estimation of portion size included household measures, weight/volume, standard units and portions, bread shapes and photographs. The set of photographs was developed in the context of the EPIC study(Reference van Kappel, Amoyel and Slimani20). Each centre chose from the EPIC portfolio of photographs the pictures that best represented their national food habits.

In the absence of harmonised recent food composition tables (FCT) including all countries of our assessment, protein and K contents in foods were calculated using country-specific FCT(2124). Carbohydrates, total fat, saturated fat, alcohol and dietary fibre intake as well as energy content were also calculated. We calculated energy values by summing the contributions from protein, carbohydrates, fat and alcohol and using related Atwater factors (17, 17, 37 and 29 kJ/g, respectively). In the Czech Republic, the national FCT was published about 20 years ago. Therefore, a FCT was compiled for EFCOVAL purposes in the Czech Republic with composition of most foods based on the Slovakian tables(25). In all the centres, missing nutrient data for a food was imputed from a similar food or another FCT, based on country-specific decisions; but in a few cases, this was not possible for K, saturated fat, dietary fibre and alcohol. The percentage of missing values was < 6 % of all reported foods for all nutrients.

Twenty-four hour urine collections and recovery biomarkers

The subjects were instructed not to make use of acetaminophen painkillers, such as paracetamol, and sulphonamide drugs, during the days of urine collection. To check the completeness of urinary collections, one tablet of 80 mg PABA (PABAcheck; Laboratories for Applied Biology, London, UK) had to be taken three times on the day of the urine collection: with the morning, midday and evening meals. Hence, we expected that 240 mg of PABA would be almost completely excreted within 24 h(Reference Bingham and Cummings26, Reference Runswick, Slothouber and Boeing27). The collection of the 24 h urine started with voiding and discarding the first urine in the morning after waking up. Subsequently, the urine excreted during the next 24 h, up to and including the first voiding of the following day, was collected. For this purpose, each subject received labelled containers (at least two), one funnel to help the collection, one safety pin to be fixed in the underwear as a reminder for collection and a diary scheme booklet to register the timing, observations (e.g. use of medication and supplements) and possible deviations (e.g. missing urine) of the urine collection protocol. Boric acid (3 g/2 litre bottle) was used as preservative. The subjects provided their urine samples to the dietitians at the study centre when a face-to-face dietary recall was scheduled. If the 24 h recall interview was by phone, urine samples were collected at the subject's home or delivered to the study centre. When a long period was anticipated between the end of the collection and the receiving of samples, subjects were instructed to keep the urine samples at approximately 4°C, which in most cases was not more than 12 h. To verify the stability of PABA in urine, a pooled urine sample of three participants from The Netherlands were kept at four different temperatures ( − 20, 6, 20 and 30°C) for 8 d. At five moments (days 0, 1, 2, 4 and 7), PABA concentrations were measured. No significant changes in PABA concentrations were observed during the storage period at each temperature. The regression equation for PABA content as a function of time during storage at 20°C (assumed to be the most common storage temperature) was as follows: PABA (mg/l) = 140·2, − 0·8 (time in days) with the 95 % CI for the time coefficient being − 2·5, 0·8.

At the laboratory of the local centres, urine was mixed, weighed and aliquoted. Then, the specimens were stored at − 20°C until shipment on dry ice to the central laboratory at Wageningen University, where they were kept at the same temperature.

Chemical analysis

On the day of chemical analysis, aliquots were rapidly thawed at room temperature. Urinary N was determined colorimetrically by the Kjeldahl technique on a Kjeltec 2300 analyser (Foss, Hilleroed, Denmark) after destruction of the sample with concentrated sulphuric acid. Urinary K was measured by an ion-selective electrode on a Synchron LX20 analyzer (Beckman Coulter, Mijdrecht, The Netherlands). PABA was measured by colorimetry(Reference Bingham, Williams and Cole28). The intra-assay precision, expressed as CV, of these three analyses was < 2 %. Taking into account the extra-renal losses (approximately 19 %) and the fact that protein on average contains 16 % N, urinary protein was calculated as (6·25 × (urinary N/0·81))(Reference Bingham and Cummings15, Reference Bingham29). Urinary K was estimated by dividing the measured value by 0·77, assuming that 77 % of K intake is excreted through the urine when considering faecal excretion(Reference Tasevska, Runswick and Bingham16, Reference Holbrook, Patterson and Bodner30).

Urine samples with PABA recoveries < 50 % were treated as incomplete and excluded from the data analysis (n 14). Additionally, the subjects who took drugs containing sulphonamides or acetaminophen or one who took less than three PABA tablets had their urine diaries checked for other deviations in the urine collection. In cases where other deviations were observed, namely urine loss during the collection or absent registration of collection time, samples were excluded from the analysis (n 4). Otherwise, samples were included (n 13) as we did not want to exclude potentially complete urines. Results of the present paper did not change by excluding these subjects. As described before(Reference Johansson, Bingham and Vahter31), specimens containing between 50 and 85 % of PABA recovery (n 105) had their urinary concentrations proportionally adjusted to 93 % of PABA recovery. Recoveries >85 % were included in data analyses without adjustments (n 1062).

Data analysis

The analyses were performed using SAS statistical package, version 9.1 (SAS Institute, Inc., Cary, NC, USA). The statistical analyses were stratified by sex and using the average of 2 d of intake and excretion, except for eighteen subjects who only had 1 d of 24 h recall and biomarker. For these subjects, the 24 h recall matched with the day of the urine collection. To assess the presence of bias (systematic errors), the mean difference between nutrient intake and excretion was calculated. ANCOVA followed by the Tukey post hoc test was used for testing whether biases differed between the centres. The ANCOVA model included age (continuous), education level (three categories) and BMI (continuous), given that stratified analysis of these variables showed us differential performance of the method within and between the centres. To estimate and compare the distribution of usual intake and excretion of protein and K between the centres, the multiple source method (MSM) was used as the measurement error model(32). This model removes the effect of day-to-day variability and random error in the two 24 h recalls and biomarker estimates. The MSM was developed in the framework of the EFCOVAL study and enabled us to estimate individual usual intake. We decided not to use covariates in the calculation of usual intakes with the MSM. Plots of usual intake distributions based on the 24 h recall and biomarker were created using R software, version 2.8.1 (http://CRAN.R-project.org). The percentages of subjects consuming above certain cut-off points for each distribution curve were calculated. For both sexes, we specified eleven cut-off points to cover the whole range of protein and K intake among the five centres. For the evaluation of ranking of individuals, we computed Pearson's correlation coefficients. For adjusted correlations, we used usual intake and excretion data corrected for within-person variability, as estimated by the MSM, and further corrected for age, BMI and education level by using partial Pearson correlations. CI of the correlations were obtained using the Fisher Z-transformation(Reference Kleibaum, Kupper and Nizam33). Energy-adjusted correlations were calculated using the residual method(Reference Willett, Howe and Kushi34). To test the equality of correlations, pairwise comparisons were made using Fisher Z-transformation(Reference Kleibaum, Kupper and Nizam33). Pooled correlations of the five centres were calculated by first converting the correlations into a standard normal metric (Fisher's r-to-Z transformation). Next, the pooled average was calculated, in which each transformed correlation coefficient was weighted by its inverse variance, followed by the back transformation(Reference Kleibaum, Kupper and Nizam33). The Cochrane Q test was used for testing the heterogeneity of the pooled correlation(Reference Field35).

Results

The mean age of the subjects was similar in the five centres (Table 1). In both sexes, mean BMI was comparable across the centres (ranges 23·2–25·5 kg/m2 in women and 25·5–27·9 kg/m2 in men). Subjects with moderate and high education levels were over-represented in the study compared with individuals with a low education level, especially men in Norway. The variations in energy intake across the centres were less pronounced than in macronutrients, especially for carbohydrates.

Table 1 Characteristics of five European centres in the European Food Consumption Validation Study*

(Mean values with their standard errors)

BE, Belgium; CZ, Czech Republic; FR, France; NL, The Netherlands; NO, Norway.

* Dietary intake based on 2 × 24 h recalls.

A degree of underestimation was seen in the assessment of protein intake in all the centres. Underestimation varied from 2·7 % (Norway) to 12·4 % (The Netherlands) in men and from 2·3 % (Norway) to 12·8 % (France) in women, based on the crude differences between intake and excretion (Table 2). After adjusting for age, BMI and education level, the bias did not differ between the centres for women. However, men in the Czech Republic had a significantly smaller bias compared with those in France and The Netherlands. For K, the underestimation varied from 1·7 % in Norway to 17·1 % in France for men and from 6·6 % in The Netherlands to 13 % in France for women. An overestimation of 5·9 % for men and 1·6 % for women was found in the Czech Republic. A statistically significant difference in the adjusted bias was seen in men between France and three other centres: Belgium, the Czech Republic, The Netherlands. In women, differences were statistically significant only between France and the Czech Republic. BMI was the only factor influencing the differences between the countries at a significant level (P < 0·01 for all analyses, except for K in women; P = 0·16). Upon inclusion of energy intake into the ANCOVA model, the conclusion about the differences between the centres changed only for protein results in men, which lost statistical significance (P = 0·08). Additionally, when we pooled the data from all the countries, no consistent trend in mean protein and K biases was observed across the different education levels and modes of administration (data not shown).

Table 2 Protein and potassium intake and excretion based on 2×24 h recalls and 2×24 h urinary biomarkers for five European centres in the European Food Consumption Validation Study

(Mean values with their standard errors)

BE, Belgium; CZ, Czech Republic; FR, France; NL, The Netherlands; NO, Norway.

a,b Mean values with unlike superscript letters were significantly different between the countries (P < 0·05).

* One-way ANCOVA (general linear model) based on mean difference between intake and excretion. Tukey's post hoc test was used for pairwise comparison between the countries. ANCOVA model included age, BMI and educational level.

Urinary protein = (urinary N/0·81) × 6·25(Reference Bingham and Cummings15).

Urinary K = (urinary K/0·77)(Reference Tasevska, Runswick and Bingham16).

The bias in mean intake can also be observed when comparing the distributions of usual intake based on food consumption data with those obtained from excretion data (the supplementary material for this article can be found at http://www.journals.cambridge.org/bjn). The intake data curve shifted somewhat to the left (underestimation of intake) for almost all the centres compared with the excretion data. Since the prevalence of subjects consuming below or above a certain cut-off point is an important indicator for a population's nutritional status, we assessed and compared the prevalence of subjects consuming above specific cut-off points for both usual intake and usual excretion distributions (see Fig. 1 for results of protein in males and the supplementary material ‘for results in females and results of K in both sexes’ can be found at http://www.journals.cambridge.org/bjn). Overall, we found a fair agreement between prevalences estimated based on the intake and excretion data at the lower end of the usual protein and K intake distribution, but larger differences at middle cut-off levels. For protein in men, the smallest differences in prevalence between intake and excretion were seen in Norway (up to 15 %) and the largest ones in France (up to 46 %) and The Netherlands (up to 41 %). For women, the smallest differences were seen in Norway (up to 11 %) and the largest ones in the Czech Republic (up to 38 %) and France (up to 55 %). The smallest difference between K intake and excretion distribution in males was observed in The Netherlands (up to 7 %) while the larger differences were seen in the Czech Republic and France (up to 21 and 40 %, respectively). In women, France was the centre with the largest difference (up to 29 %) between K usual intake and excretion, and The Netherlands the smallest (up to 17 %).

Fig. 1 Prevalence of men consuming above specific amounts of protein as estimated by usual intake distributions (an usual intake/excretion distribution estimated by the multiple source method (see ‘Methods’ section)) from dietary recalls (intake) and biomarkers (excretion) for five European centres in the European Food Consumption Validation Study. (a) Belgium, (b) Czech Republic, (c) France, (d) The Netherlands, (e) Norway. - -●- -, Intake; –○–, excretion.

Unadjusted Pearson correlation coefficients between average protein intake and its biomarker within centres ranged between 0·42 and 0·65 in men and between 0·46 and 0·57 in women (Table 3). After adjusting for within-person variability, age, BMI and education level, correlations ranged between 0·43 and 0·67 in men and between 0·39 and 0·63 in women. For K, unadjusted correlations ranged between 0·45 and 0·65 in men and between 0·31 and 0·69 in women. Adjusted correlations ranged between 0·40 and 0·69 in men and between 0·37 and 0·68 in women. For both protein and K, adjusting only for the within-person variability slightly increased the correlations between intake and excretion (data not shown). Statistically significant differences between correlation coefficients were only found between Belgium and the Czech Republic (P = 0·04) for unadjusted correlations of K in women. However, after adjusting the correlations for energy, we found a significant difference between the Czech Republic (r 0·25) and France (r 0·65) for protein intake in men (P = 0·01).

Table 3 Pearson coefficients of correlation between protein intake and urinary excretion* for five European centres in the European Food Consumption Validation Study

(Mean values and 95 % confidence intervals)

* Average intake and excretion based on 2 d of collection.

Pairwise comparisons between countries (by Fisher Z transformation) suggested differences for unadjusted correlations between Belgium and the Czech Republic in females and between France and the Czech Republic for energy-adjusted correlations in males.

Adjusted for the within-person variability using the usual intake/excretion data as estimated by the multiple source method (see ‘Methods’ section) and adjusted for age, BMI and educational level using partial Pearson correlations.

§ Same adjustments as previous correlation plus energy-adjustment by the residual method.

Mean values for heterogeneity were not significant for all the analyses (P>0·05).

The pooled adjusted correlations in males and females were 0·51 and 0·60 for protein and 0·56 and 0·57 for K intake, respectively.

Discussion

In the present study, we compared the validity of usual protein and K intake estimated from two non-consecutive standardised 24 h recalls between five selected centres in Europe. On average, men and women under-reported protein intake from the two 24 h recalls by 8 %. For K intake, average underestimation was 7 % for men and 4 % for women.

Protein intake was markedly underestimated (approximately 12 %) in French and Dutch men, especially when compared with Czech Republic men. The same is true for K intake in French men. In women, underestimation of mean protein intake was present in all the centres and appeared to be comparable across the centres. For K intake, however, the underestimation observed in the French centre was not comparable to that of the other centres, particularly to the overestimation observed in the Czech Republic. Furthermore, we assessed the agreement between the percentage of subjects above a certain cut-off point based on 24 h recall and biomarker data. We found a fair agreement for cut-off points at the lower end of the distribution ( < 10 % difference), but larger differences at other points of the intake distribution (up to 55 % difference for protein in French females). Finally, we observed moderate correlations for the ranking of individuals, which were likely to be comparable across the centres.

The results from the EPIC study, using EPIC-Soft in different centres, revealed a similar or even higher underestimation of protein intake collected from a single day (average of 13 % in men and 19 % in women)(Reference Slimani, Bingham and Runswick17). The OPEN study in the United States, which assessed the structure of dietary measurement error in 24 h recalls collected twice, has also shown a similar underestimation of protein intake (11–15 %)(Reference Subar, Kipnis and Troiano36). A few other studies indicated overestimation of protein (about 7 % for the whole population)(Reference Kahn, Whelton and Appel37). For K, studies indicated overestimation of intake up to 20 %(Reference Heerstrass, Ocke and Bueno-de-Mesquita38Reference Bingham and Day40), similar to what we observed in the Czech Republic. Nevertheless, because of methodological differences, the comparison of bias estimates between the present study and other studies is not straightforward. For example, adjustment of N and K excretions to extra-renal losses was not consistently performed among the studies. In addition, the completeness of 24 h urine collections was not always assessed. Although we acknowledge the differences in methodology between the studies, the performance of these two standardised 24 h recalls on assessing the mean protein and K intake appeared to provide alike or even more accurate results than what have been presented in the literature so far.

In terms of assessing the whole distribution of intake, two 24 h recalls used in the study by Freedman et al. (Reference Freedman, Midthune and Carroll39) underestimated the usual protein intake in all points of the distribution, especially at the lower end. Moreover, they found a good agreement between K intake and excretion in the whole range of percentiles. In contrast, moderate to large discrepancies were found between 24 h recall and biomarker data distributions in the present study, but not at the lower end of the distribution. The present results suggest that the assessment of protein and K inadequacy at the population level by two non-consecutive 24 h recalls in healthy European populations is, therefore, appropriate.

Independent of the size of the bias, the correct classification of individuals according to their intake is also informative on the quality of the dietary assessment. The correlations presented in the present paper are considerably higher compared with many other studies(Reference Subar, Kipnis and Troiano36, Reference Bingham, Gill and Welch41Reference Shai, Rosner and Shahar43). Based on this, we conclude that the method performed sufficiently for the ranking of individuals, adding evidence to the use of this standardised 24 h recall. When we adjusted the nutrient values for energy intake, this changed the correlations in both directions and resulted in more noticeable differences across the centres. We doubt, however, whether energy-adjusted values will be our main exposure of interest in future monitoring surveys and whether individual energy intake was correctly estimated using only 2 d of 24 h recall. Therefore, we do not base the conclusions of the present paper on the energy-adjusted results.

We suppose that the differences found in the size and direction of the bias (i.e. overestimation of K intake in the Czech Republic and underestimation of both K and protein in the other centres) between the centres may be explained by reasons related to characteristics of the population and of the method itself. We have controlled our statistical analyses for the influence of age, education level and BMI. As a result, BMI was the only factor significantly influencing the differences between the countries. This is in accordance with our expectations since other studies have revealed a differential under-reporting of dietary intake by subgroups of BMI(Reference Heerstrass, Ocke and Bueno-de-Mesquita38, Reference Lissner, Troiano and Midthune44). Nevertheless, other aspects of the population could have affected the validity of the method between the centres in a different manner, i.e. factors related to the food pattern of the centres. Due to cultural differences in food pattern, it is expected that predominant food items contributing to protein and K intake across European countries will be different(Reference Halkjaer, Olsen and Bjerregaard45, Reference Slimani, Fahey and Welch46). For example, the food group ‘dairy products’ was one of the major contributors (>22 %) to the protein intake in The Netherlands and Norway (in males only), whereas in the other three centres, ‘meat products’ was distinctly the major contributor (>30 %). Knowing that the errors in the assessment of different food groups differ, as for instance in the portion size estimation(Reference Rumpler, Kramer and Rhodes47), differences in validity between the centres could be expected. Likewise, differences in the consumption of composite foods could have had an effect since it is more difficult to recall all ingredients of composite foods than a single food item(Reference Cosgrove, Flynn and Kiely48, Reference O'Brien, Kiely and Galvin49).

Another important factor that could explain the differences between countries is the use of not harmonised FCT across the centres. Use of different conversion factors as well as distinct laboratory analyses to produce food nutrient contents across the tables is just an example which could have caused biases not to be comparable. For instance, for three of the FCT used in EFCOVAL, protein figures were calculated from N contents using the so-called ‘Jones conversion factors’(Reference Jones50) or slight modifications of them. However, in the Dutch tables, only two of these factors were used (6·38 for milk products and 6·25 for all other foods), and in the compiled Czech table, only one factor (6·25) was applied (Slovakian tables). Since errors attributed to these differences can be proportional to the level of intake, it is impossible to conclude on the influence of using different conversion factors in the comparison between the countries. Nevertheless, further investigation about the use of these conversion factors in FCT for comparisons of nutrient intake between countries is warranted.

The present study adds value to the present knowledge of collecting dietary information using standardised 24 h recalls for possible use in national monitoring surveys. An important strength of the present study was the collection of 2 d of both dietary intake and biomarkers allowing the quantification of within-person variability and to estimate the usual intake distributions. A potential limitation of the present study is that a health-conscious sample may have been included, hampering the extrapolation of the results to the general population. However, the present results suggested that extrapolation to other populations could be done irrespective of their education level. In addition, the generalisability of protein and K results to other nutrients of interest should be done with care. Although we might want to assume that the validation results of a single nutrient can be used as a proxy to other nutrients, there is evidence nowadays that some foods and consequently related nutrients might be selectively misreported(Reference Rumpler, Kramer and Rhodes47, Reference Pryer, Vrijheid and Nichols51). Besides, only 2 d of 24 h recall were used in our assessment while the inclusion of more than 2 d may be necessary to improve the use of this 24 h recall in the assessment of other nutrient intake distributions, particularly the infrequently consumed ones(Reference Palaniappan, Cue and Payette52). The statistical adjustments performed with the MSM intended to remove the day-to-day variation in intakes and assess the usual distributions of intake. But, if the variance of the nutrient intake is not reliably estimated from 2 d of intake, then the observed intake may shrink too much or too little toward the group mean intake, resulting in an inaccurate usual intake distribution(Reference Carriquiry53). The use of FFQ combined with 24 h recalls may be an option in future monitoring surveys for the calculation of usual intakes of infrequently consumed nutrients, as more days of 24 h recalls are demanding and expensive. Furthermore, the reliability of the conversion factors used to adjust urinary protein and K in our analyses can be questioned. With the assumption that the subjects were in N balance, these factors have been based on rigorously controlled feeding studies(Reference Bingham and Cummings15, Reference Tasevska, Runswick and Bingham16) and in the case of protein confirmed by Kipnis et al. (Reference Kipnis, Midthune and Freedman54). Lastly, we have collected data in The Netherlands 6 months before the other centres and this may have influenced the results. Nevertheless, while the data for The Netherlands were collected in spring/summer, the data for other four countries were collected in the winter/spring. However, since minor adjustments were done in the study protocols and the differences in seasonality were small for protein and K intake, it is unlikely that a different period influences the present results.

To conclude, first, the ability of the two non-consecutive standardised 24 h recalls using EPIC-Soft software appears to be sufficiently valid for assessing and comparing the mean protein and K intake across the centres. When comparing populations in a future nutrition monitoring system, the variability in the nutrient biases of 4–7 % across the centres needs to be considered. Second, the method seems to be sufficiently valid for assessing and comparing the protein and K inadequacy of healthy populations across the centres and less appropriate to assess other points of the intake distribution. Third, the ability to rank the individuals according to protein and K intake within the centres is comparable between them, which substantiates the validity of the method. Therefore, this standardised two non-consecutive 24 h recalls, further adapted and validated in the EFCOVAL project, appear appropriate to be used in the context of a future pan-European dietary monitoring system. Built on EFCOVAL and EPIC experiences, improvements may be possible for the employment of this methodology by an even higher standardisation setting (e.g. conversion factors), which could result in an enhanced validity of the method, and thus comparability between the countries.

Acknowledgements

The authors thank the EFCOVAL partners for their useful advices. The EFCOVAL partners are Ghent University (DPH), Belgium; Academy of Medical Sciences (AMZH), Croatia; National Institute of Public Health (NIPH), Czech Republic; National Food Institute, Technical University of Denmark (DTU), Denmark; French Food Safety Authority (AFSSA), France; National Institute for Agricultural Research (INRA), France; German Institute of Human Nutrition (DIfE), Germany; National Research Institute for Food and Nutrition (INRAN), Italy; Wageningen University (WU), The Netherlands; National Institute for Public Health and the Environment (RIVM), The Netherlands; University of Oslo, Norway; Basque Foundation for Health Innovation and Research (BIOEF), Spain; Prima informatics limited (Primainfo), United Kingdom; and International organization, International Agency for Research on Cancer (IARC, WHO). The present document reflects only the authors' views and the European Community is not liable for any use that may be made of the information contained therein. The Community funding under the Sixth Framework Program for the EFCOVAL project is acknowledged (FOOD-CT-2006-022895). S. P. C. carried out data analyses and wrote the paper, taking into account the comments from all the co-authors. J. H. M. d. V., A. G. and P. v. V. designed and coordinated the validation study. O. W. S. contributed to the statistical analyses. P. J. M. H. was responsible for the laboratorial analyses. J. H. M. d. V., A. G., L. L., A.-S. R., I. T. L. L., L. F. A., I. H., W. D. K., J. R., M. D. and M. C. O. were involved in the fieldwork and gave input on interpretation of results. E. d. B., N. S. and P. v. V. were the overall coordinators of the EFCOVAL project. All the co-authors commented on the paper and approved the final version. None of the other authors had a financial conflict of interest.

References

1 EFSA (2009) General principles for the collection of national food consumption data in the view of a pan-European dietary survey. EFSA J 7, 151.Google Scholar
2 Charzewska, J (1994) Gaps in dietary-survey methodology in eastern Europe. Am J Clin Nutr 59, Suppl. 1, 157S160S.Google Scholar
3 Pietinen, P & Ovaskainen, ML (1994) Gaps in dietary-survey methodology in Western Europe. Am J Clin Nutr 59, Suppl. 1, 161S163S.Google Scholar
4 Verger, P, Ireland, J, Moller, A, et al. (2002) Improvement of comparability of dietary intake assessment using currently available individual food consumption surveys. Eur J Clin Nutr 56, Suppl. 2, S18S24.Google Scholar
5 Brussaard, JH, Johansson, L & Kearney, J (2002) Rationale and methods of the EFCOSUM project. Eur J Clin Nutr 56, Suppl. 2, S4S7.Google Scholar
6 Biro, G, Hulshof, KF, Ovesen, L, et al. (2002) Selection of methodology to assess food intake. Eur J Clin Nutr 56, Suppl. 2, S25S32.Google Scholar
7 Brussaard, JH, Lowik, MR, Steingrimsdottir, L, et al. (2002) A European food consumption survey method – conclusions and recommendations. Eur J Clin Nutr 56, Suppl. 2, S89S94.Google Scholar
8 Sliman, N & Valsta, L (2002) Perspectives of using the EPIC-SOFT programme in the context of pan-European nutritional monitoring surveys: methodological and practical implications. Eur J Clin Nutr 56, Suppl. 2, S63S74.CrossRefGoogle Scholar
9 Conway, JM, Ingwersen, LA & Moshfegh, AJ (2004) Accuracy of dietary recall using the USDA five-step multiple-pass method in men: an observational validation study. J Am Diet Assoc 104, 595603.CrossRefGoogle ScholarPubMed
10 Statistics Canada (2004) Canadian Community Health Survey, Cycle 2.2, Nutrition. A Guide to Accessing and Interpreting the Data. Ottawa: Minister of Health.Google Scholar
11 McLennan, W & Podger, A (1995) National Nutrition Survey User's Guide. Canberra: Australian Government Printing Service.Google Scholar
12 Russell, D, Parnell, W & Wilson, N (1999) NZ Food: NZ People. Key Results of the 1997 National Nutrition Survey. Wellington: Minister of Health.Google Scholar
13 Willett, W (1998) Nutritional Epidemiology, vol. XIV, pp. 514. New York: Oxford University Press.Google Scholar
14 Ocke, MC & Kaaks, RJ (1997) Biochemical markers as additional measurements in dietary validity studies: application of the method of triads with examples from the European Prospective Investigation into Cancer and Nutrition. Am J Clin Nutr 65, Suppl. 4, 1240S1245S.CrossRefGoogle ScholarPubMed
15 Bingham, SA & Cummings, JH (1985) Urine nitrogen as an independent validatory measure of dietary intake: a study of nitrogen balance in individuals consuming their normal diet. Am J Clin Nutr 42, 12761289.Google Scholar
16 Tasevska, N, Runswick, SA & Bingham, SA (2006) Urinary potassium is as reliable as urinary nitrogen for use as a recovery biomarker in dietary studies of free living individuals. J Nutr 136, 13341340.CrossRefGoogle ScholarPubMed
17 Slimani, N, Bingham, S, Runswick, S, et al. (2003) Group level validation of protein intakes estimated by 24-hour diet recall and dietary questionnaires against 24-hour urinary nitrogen in the European Prospective Investigation into Cancer and Nutrition (EPIC) Calibration Study. Cancer Epidemiol Biomarkers Prev 12, 784795.Google Scholar
18 Slimani, N, Deharveng, G, Charrondiere, RU, et al. (1999) Structure of the standardized computerized 24-h diet recall interview used as reference method in the 22 centers participating in the EPIC project. European Prospective Investigation into Cancer and Nutrition. Comput Methods Programs Biomed 58, 251266.Google Scholar
19 Slimani, N, Ferrari, P, Ocke, M, et al. (2000) Standardization of the 24-hour diet recall calibration method used in the european prospective investigation into cancer and nutrition (EPIC): general concepts and preliminary results. Eur J Clin Nutr 54, 900917.Google Scholar
20 van Kappel, AL, Amoyel, J, Slimani, N, et al. (1995) EPIC-Soft Picture Book for Estimation of Food Portion Sizes, I. Report, Editor. International Agency for Research on Cancer: Lyon, France.Google Scholar
21 NEVO-TABEL (2006) Nederlandse Voedings Middelen Tabel (Dutch Food Composition Tables). Den Haag: Voedingscentrum.Google Scholar
22 NUBEL (2004) Belgische Voedingsmiddelentabel (Nutrients Belgium). Brussels: VZW Nubel.Google Scholar
23 AFSSA/CIQUAL (2008) French Food Composition Table. Paris: F.F.S. Agency.Google Scholar
24 The Norwegian Food Safety Authority, The Norwegian Directorate of Health and the University of Oslo (2006) The Norwegian Food Composition Table 2006. www.matportalen.no/matvaretabellen.Google Scholar
25 Food Research Institute (1997–2002) Pozivatinove tabulky (Food Tables). Bratislava: Food Research Institute.Google Scholar
26 Bingham, S & Cummings, JH (1983) The use of 4-aminobenzoic acid as a marker to validate the completeness of 24 h urine collections in man. Clin Sci (Lond) 64, 629635.CrossRefGoogle ScholarPubMed
27 Runswick, S, Slothouber, B, Boeing, H, et al. (2002) Compliance with the urine marker PABAcheck in cancer epidemiology studies. IARC Sci Publ 156, 3537.Google ScholarPubMed
28 Bingham, SA, Williams, R, Cole, TJ, et al. (1988) Reference values for analytes of 24-h urine collections known to be complete. Ann Clin Biochem 25 (Pt 6), 610619.Google Scholar
29 Bingham, SA (2003) Urine nitrogen as a biomarker for the validation of dietary protein intake. J Nutr 133, Suppl. 3, 921S924S.Google Scholar
30 Holbrook, J, Patterson, K, Bodner, J, et al. (1984) Sodium and potassium intake and balance in adults consuming self-selected diets. Am J Clin Nutr 40, 786793.CrossRefGoogle ScholarPubMed
31 Johansson, G, Bingham, S & Vahter, M (1999) A method to compensate for incomplete 24-hour urine collections in nutritional epidemiology studies. Public Health Nutr 2, 587591.Google Scholar
32 German Institute of Human Nutrition (2009) The Multiple Source Method (MSM). 07 Sep 2009 (cited 2009 1 November 2009). https://nugo.dife.de/msm/.Google Scholar
33 Kleibaum, DG, Kupper, LL, Nizam, A, et al. (2008) Applied regression analysis and other multivariable methods. In Duxbury Applied Series, 4th ed., pp. 906. Belmont, CA: Thomson Brooks/Cole.Google Scholar
34 Willett, W, Howe, G & Kushi, L (1997) Adjustment for total energy intake in epidemiologic studies. Am J Clin Nutr 65, 1220S1228S.Google Scholar
35 Field, AP (2005) Is the meta-analysis of correlation coefficients accurate when population correlations vary? Psychol Methods 10, 444467.CrossRefGoogle ScholarPubMed
36 Subar, AF, Kipnis, V, Troiano, RP, et al. (2003) Using intake biomarkers to evaluate the extent of dietary misreporting in a large sample of adults: The OPEN Study. Am J Epidemiol 158, 113.CrossRefGoogle Scholar
37 Kahn, HA, Whelton, PK, Appel, LJ, et al. (1995) Validity of 24-hour dietary recall interviews conducted among volunteers in an adult working community. Ann Epidemiol 5, 484489.Google Scholar
38 Heerstrass, DW, Ocke, MC, Bueno-de-Mesquita, HB, et al. (1998) Underreporting of energy, protein and potassium intake in relation to body mass index. Int J Epidemiol 27, 186193.Google Scholar
39 Freedman, LS, Midthune, D, Carroll, RJ, et al. (2004) Adjustments to improve the estimation of usual dietary intake distributions in the population. J Nutr 134, 18361843.CrossRefGoogle ScholarPubMed
40 Bingham, SA & Day, NE (1997) Using biochemical markers to assess the validity of prospective dietary assessment methods and the effect of energy adjustment. Am J Clin Nutr 65, Suppl. 4, 1130S1137S.Google Scholar
41 Bingham, SA, Gill, C, Welch, A, et al. (1997) Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers. Int J Epidemiol 26, Suppl. 1, S137S151.CrossRefGoogle ScholarPubMed
42 Olafsdottir, AS, Thorsdottir, I, Gunnarsdottir, I, et al. (2006) Comparison of women's diet assessed by FFQs and 24-hour recalls with and without underreporters: associations with biomarkers. Ann Nutr Metab 50, 450460.Google Scholar
43 Shai, I, Rosner, BA, Shahar, DR, et al. (2005) Dietary evaluation and attenuation of relative risk: multiple comparisons between blood and urinary biomarkers, food frequency, and 24-hour recall questionnaires: the DEARR study. J Nutr 135, 573579.CrossRefGoogle ScholarPubMed
44 Lissner, L, Troiano, RP, Midthune, D, et al. (2007) OPEN about obesity: recovery biomarkers, dietary reporting errors and BMI. Int J Obes (Lond) 31, 956961.Google Scholar
45 Halkjaer, J, Olsen, A, Bjerregaard, LJ, et al. (2009) Intake of total, animal and plant proteins, and their food sources in 10 countries in the European Prospective Investigation into Cancer and Nutrition. Eur J Clin Nutr 63, Suppl. 4, S16S36.CrossRefGoogle ScholarPubMed
46 Slimani, N, Fahey, M, Welch, AA, et al. (2002) Diversity of dietary patterns observed in the European Prospective Investigation into Cancer and Nutrition (EPIC) project. Public Health Nutr 5, 13111328.Google Scholar
47 Rumpler, WV, Kramer, M, Rhodes, DG, et al. (2008) Identifying sources of reporting error using measured food intake. Eur J Clin Nutr 62, 544552.Google Scholar
48 Cosgrove, M, Flynn, A & Kiely, M (2005) Impact of disaggregation of composite foods on estimates of intakes of meat and meat products in Irish adults. Public Health Nutr 8, 327337.Google Scholar
49 O'Brien, MM, Kiely, M, Galvin, M, et al. (2003) The importance of composite foods for estimates of vegetable and fruit intakes. Public Health Nutr 6, 711726.Google Scholar
50 Jones, DB (1941) Factors for Converting Percentages of Nitrogen in Foods and Feeds into Percentages of Protein. Washington, DC: US Department of Agriculture-Circle.Google Scholar
51 Pryer, JA, Vrijheid, M, Nichols, R, et al. (1997) Who are the ‘low energy reporters’ in the dietary and nutritional survey of British adults? Int J Epidemiol 26, 146154.Google Scholar
52 Palaniappan, U, Cue, RI, Payette, H, et al. (2003) Implications of day-to-day variability on measurements of usual food and nutrient intakes. J Nutr 133, 232235.Google Scholar
53 Carriquiry, AL (2003) Estimation of usual intake distributions of nutrients and foods. J Nutr 133, 601S608S.Google Scholar
54 Kipnis, V, Midthune, D, Freedman, LS, et al. (2001) Empirical evidence of correlated biases in dietary assessment instruments and its implications. Am J Epidemiol 153, 394403.CrossRefGoogle ScholarPubMed
Figure 0

Table 1 Characteristics of five European centres in the European Food Consumption Validation Study*(Mean values with their standard errors)

Figure 1

Table 2 Protein and potassium intake and excretion based on 2×24 h recalls and 2×24 h urinary biomarkers for five European centres in the European Food Consumption Validation Study(Mean values with their standard errors)

Figure 2

Fig. 1 Prevalence of men consuming above specific amounts of protein as estimated by usual intake distributions (an usual intake/excretion distribution estimated by the multiple source method (see ‘Methods’ section)) from dietary recalls (intake) and biomarkers (excretion) for five European centres in the European Food Consumption Validation Study. (a) Belgium, (b) Czech Republic, (c) France, (d) The Netherlands, (e) Norway. - -●- -, Intake; –○–, excretion.

Figure 3

Table 3 Pearson coefficients of correlation between protein intake and urinary excretion* for five European centres in the European Food Consumption Validation Study†(Mean values and 95 % confidence intervals)

Supplementary material: PDF

Crispim supplementary material

Figure S1

Download Crispim supplementary material(PDF)
PDF 25.8 KB
Supplementary material: PDF

Crispim supplementary material

Figure S2

Download Crispim supplementary material(PDF)
PDF 65.9 KB
Supplementary material: PDF

Crispim supplementary material

Figure S3

Download Crispim supplementary material(PDF)
PDF 28.1 KB
Supplementary material: PDF

Crispim supplementary material

Figure S4

Download Crispim supplementary material(PDF)
PDF 26.7 KB
Supplementary material: PDF

Crispim supplementary material

Figures S5-S7

Download Crispim supplementary material(PDF)
PDF 120.6 KB