In Germany and most Western countries, dairy intake is an essential part of a healthy diet for children and adolescents because of its nutrient composition that is beneficial for healthy growth (e.g. protein or Ca content)(Reference Comerford, Miller and Boileau1,Reference Gaucheron2) . Nevertheless, dairy intake has been linked with both positive(Reference Barrubés, Babio and Becerra-Tomás3–Reference Schwingshackl, Schwedhelm and Hoffmann5) and potential negative effects on human health(Reference Lu, Chen and Niu6–Reference Aghasi, Golzarand and Shab-Bidar8). Systematic reviews and meta-analyses of randomised clinical trials analysing the relationship between dairy intake and low-grade systemic inflammation in adults showed neutral to beneficial effects on various inflammatory biomarkers(Reference Moosavian, Rahimlou and Saneei9–Reference Nieman, Anderson and Cifelli13). In a systematic review, dairy intake was associated with anti-inflammatory activity among people with metabolic disorders and pro-inflammatory activity among subjects allergic to cow’s milk(Reference Bordoni, Danesi and Dardevet14). According to the authors, these opposing effects can be attributed to the hypersensitive reaction and the resulting pro-inflammatory state in subjects with bovine milk allergy. Apart from these subgroup-specific pathophysiological mechanisms, the heterogeneity of this specific food group could be relevant for metabolic effects. Based on the nutrient content (fat and sugar), processing methods (fermentation) or consistency (liquid and solid), different types of dairy may have different effects on metabolism(Reference Godos, Tieri and Ghelfi15). In fact, the consideration of dairy types instead of overall dairy shows anti-inflammatory activity of low-fat, high-fat and fermented dairy intake(Reference Bordoni, Danesi and Dardevet14). Inflammatory processes have been identified as possible risk factors for obesity, insulin resistance, type 2 diabetes and CVD(Reference Shoelson, Lee and Goldfine16,Reference Lontchi-Yimagou, Sobngwi and Matsha17) . However, studies assessing the associations between different types of dairy and inflammatory biomarkers are scarce and have been mainly conducted in adults(Reference Aslam, Jacka and Marx18). Evidence of a potential association of total dairy (TD) or dairy types in children and adolescents is limited, although this is a vulnerable group. Nevertheless, a recent meta-analysis examining the relationship between dietary intake and biomarkers of inflammation among healthy children and adolescents revealed that there is no association between dairy intake and inflammation markers, such as C-reactive protein or IL-6, in children and adolescents(Reference Bujtor, Turner and Torres19). However, this meta-analysis included only a few studies focusing on the association with TD, and most of the included studies were cross-sectional. Although these cross-sectional studies have not shown that dairy intake as part of a healthy diet has unbeneficial associations with inflammation markers in children and adolescents in the short term(Reference Bujtor, Turner and Torres19), no studies have analysed the long-term association from adolescence to young adulthood. Furthermore, whether this potential association differs with dairy types remains unknown.
Therefore, this study aimed to investigate the relationship between dairy intake as well as intake of different types of dairy in adolescence and biomarkers of inflammation and insulin resistance in young adulthood.
Methods
DONALD study
The DOrtmund Nutritional and Anthropometric Longitudinally Designed (DONALD) study is an ongoing dynamic (open) cohort study that collects information on nutrition, growth, development and metabolism of healthy children and adolescents in Dortmund, Germany. The study was initiated with a cross-sectional sample of children and adolescents (approximately 640 participants, >2 years old) recruited in 1985. Since 1985, 35–40 infants have been enrolled annually. Eligibility criteria are healthy infants whose parents are willing to participate in a long-term study and at least one parent with sufficient knowledge of German. Participants are initially examined at the age of 3 months and return for three more follow-up visits in the first year, two in the second year and then once annually until young adulthood. Annual examinations include 3-d weighed dietary records, anthropometric measurements, 24-h urine samples, lifestyle interviews and medical examinations. Parental examinations occur every 4 years. The study was non-invasive during childhood and adolescence. Since 2005, participants aged >18 years have been invited for subsequent examinations with fasting blood samples. Further details of the study have been described elsewhere(Reference Perrar, Alexy and Nöthlings20). This study was registered in the German Register of Clinical Trials (DRKS-ID: DRKS00029092).
Study sample
In June 2019, 17 782 dietary records of 1706 children were available. Incomplete records (<3 d, n 98 records) were excluded. The inclusion criteria for the present analyses were as follows: all available data of participants who provided at least two 3-d dietary records (median number of individual 3-d dietary records: 8) during adolescence (girls, 8–15 years; boys, 9–16 years; median of all, 12 years) (Tables 1 and 2) and at least one fasting blood analysis in young adulthood (median age at blood sampling: 20·9 years). This resulted in samples from 375 participants for the current inflammation analyses. Of these, Homeostasis Model Assessment Insulin Resistance (HOMA2-IR) data were available for 371 participants, which was considered for the analyses on insulin resistance. The median follow-up of 9·2 (6·4, 12·5) years was defined as the number of years between the median age during adolescence and age during blood withdrawal.
DONALD, DOrtmund Nutritional and Anthropometric Longitudinally Designed; BMI-SDS, BMI–SD score; FMI, fat mass index; FFMI, fat-free mass index.
Values are n (%) or medians (25th, 75th percentile).
* BMI, kg/m².
† BMI > 25 kg/m².
‡ Age at minimal height velocity at the onset of the pubertal growth spurt(Reference Perrar, Buyken and Penczynski37).
|| BMI–SD score (based on the German reference percentiles for children and adolescents)(Reference Slaughter, Lohman and Boileau34).
¶ FMI (the underlying percentage body fat was estimated using the equations of Slaughter)(Reference Durnin and Womersley35).
** FFMI (the underlying percentage body fat was estimated using the equations of Slaughter)(Reference Durnin and Womersley35).
†† Paediatric cut-off values for under-reporting(Reference Kersting, Kalhoff and Lücke42).
DONALD, DOrtmund Nutritional and Anthropometric Longitudinally Designed; FMI, fat mass index; hsCRP, high-sensitivity C-reactive protein; HOMA2-IR, Homeostasis Model Assessment Insulin Resistance.
Values are n (%) or medians (25th, 75th percentile).
* Median age at adolescence and age at blood withdrawal.
† BMI, kg/m².
‡ FMI (the underlying percentage body fat was estimated using the equations of Durnin–Womersley)(Reference Buyken, Karaolis-Danckert and Remer36).
§ Smoker included current and occasional smoking and was defined based on the variable ‘smoking (yes/no/sometimes)’.
|| Alcohol consumer was defined based on the variable ‘alcohol is currently consumed (yes/no)’.
¶ Currently doing sports (organised, not organised and no sport).
** Currently employed (yes/no/rest/in education).
†† hsCRP.
‡‡ IL-6, IL-6.
§§ IL-18, IL-18.
|||| HOMA2-IR.
Dietary assessment
All food and beverages consumed and leftovers were weighed and recorded on 3 d by the parents or older participants themselves using electronic food scales. Semi-quantitative recording (e.g. spoons and cups) was allowed if accurate weighing was not possible. Information on recipes (ingredients and preparation) and types and brands of commercial food products were also required. Food group, energy and nutrient intake were calculated using our continually updated in-house food composition database LEBTAB(Reference Sichert-Hellert, Kersting and Chahda21). The composition of basic foods was retrieved from the German food composition tables BLS 3.02. The energy and nutrient contents of commercial food products, that is, canned foods, ready-to-eat meals or snacks, were estimated by recipe simulation based on the listed ingredients and nutrients(Reference Sichert-Hellert, Kersting and Chahda21).
According to Hohoff et al. (Reference Hohoff, Perrar and Jancovic22,Reference Hohoff, Perrar and Jankovic23) , the following types of dairy were included in the analyses:
TD
With regard to the nutrient content:
Low- and high-fat dairy
Low- and high-sugar dairy
With regard to the processing method:
Fermented and non-fermented dairy
With regard to the way of intake:
Liquid and solid dairy
Detailed descriptions of the different dairy types are given in Table 3. The daily intake of dairy and dairy types was calculated from the individual mean of all 3-d dietary records of participants examined during adolescence. To consider sex- and age-dependent differences in dietary intake, dairy intake was also standardised as g/1000 kcal of total energy intake.
DONALD, DOrtmund Nutritional and Anthropometric Longitudinally Designed.
* Dairy products can occur in different groups.
† Excluding cream cakes and ice cream, because they are consumed as sweets rather than to meet dairy requirements, and excluding butter.
‡ Classification based on https://www.lebensmittellexikon.de/f0000170.php.
§ Including instant powders for milk (i.e. cocoa).
|| The cut-off was set based on the 1st quartile (6.9 g added sugar/100 g) from all sweetened dairy products (n 965) reported by the study sample.
Blood analysis
Venous blood samples were drawn after an overnight fast, centrifuged at 4°C for 15 min and stored at –80°C in the DONALD study centre. Fasting plasma glucose levels were determined using a Roche/Hitachi Cobas c 311 analyzer. Plasma insulin concentrations were measured at the Laboratory for Translational Hormone Analytics of the University of Giessen using an immunoradiometric assay (IRMA, DRG Diagnostics). All other measurements were performed at the German Diabetes Center with the following assay characteristics(Reference Diederichs, Herder and Roßbach24,Reference Goletzke, Buyken and Joslowski25) : plasma high-sensitivity C-reactive protein (hsCRP) with the Roche/Hitachi Cobas c311 analyser (Roche Diagnostics), plasma high-sensitivity IL-6 using the Human IL-6 Quantikine HS, plasma adiponectin with the Human Total Adiponectin/Acrp30 Quantikine ELISA, serum leptin with the Leptin Quantikine ELISA and serum IL-18 using the Human IL-18 ELISA kit from MBL.
Definitions of the outcome variables
To estimate associations with low-grade chronic inflammation, a pro-inflammatory score was used, similar to that used by Diederichs et al. (Reference Diederichs, Herder and Roßbach24) and Penczynski et al. (Reference Penczynski, Herder and Krupp26). This score is composed of established inflammation biomarkers that are assumed to reflect low-grade inflammation better than individual markers(Reference Calder27). These include hsCRP, IL-6, IL-18, leptin and adiponectin. To approximate normal distribution, the individual biomarkers were log-transformed before standardisation (z-score) by sex (mean = 0, s d = 1). Then, these z-scores of the individual inflammation biomarkers were averaged, resulting in the pro-inflammatory score. Here, the anti-inflammatory parameter adiponectin was multiplied by −1.
Hence, the pro-inflammatory score was calculated as follows:
In addition, insulin resistance was assessed using the updated HOMA2-IR. HOMA2-IR is based on fasting insulin and blood glucose levels according to Wallace et al. (Reference Wallace, Levy and Matthews28):
HOMA2-IR was also log-transformed before standardisation (z-score) by sex (mean = 0, s d = 1) to approximate normal distribution.
Assessment of potential confounding factors
Potential confounding factors were selected on the basis of the known predictors of low-grade systemic inflammation(Reference Dowd, Zajacova and Aiello29,Reference McDade, Hawkley and Cacioppo30) and type 2 diabetes(Reference Haghighatdoost, Amini and Feizi31,Reference Okwechime, Roberson and Odoi32) . At the first visit to the study centre, early-life factors such as pregnancy characteristics and birth anthropometry of the mother and child were taken through a standardised document (Mutterpass), which is issued to women during pregnancy in Germany. Anthropometric measurements (height, weight and skinfolds) were performed by trained nurses according to standard procedures using an electronic scale (Seca 753E; Seca Weighing and Measuring System, ±100 g), a digital stadiometer (Harpenden, ±0·1 cm) and a calliper (±0·1 mm, Holtain Ltd). The participants were dressed in underwear and were barefoot. The BMI (kg/m²) was calculated as the body weight (kg) divided by the square of the body height (m2). For adolescents, sex- and age-specific standard deviation scores were calculated based on the German reference percentiles for children and adolescents(Reference Kromeyer-Hauschild, Wabitsch and Kunze33). Percent body fat was estimated using the Slaughter(Reference Slaughter, Lohman and Boileau34) and Durnin–Womersley equations(Reference Durnin and Womersley35) for adolescents and adults, respectively. Body fat mass (kg) was calculated as [(percent body fat × body mass)/100]. The corresponding indices, fat mass index and fat-free mass index, were calculated by dividing the corresponding values by the square of the body height (m2). For adolescents and young adults, the respective medians were derived for all anthropometric data. In addition, parents were weighed and measured at regular intervals using the same equipment as used for the participants. Maternal overweight was defined as a BMI ≥ 25–<30 kg/m2 and maternal obesity as a BMI ≥ 30 kg/m2. In addition, parents were asked about their family and socio-economic characteristics (e.g. maternal education). Lifestyle factors of participants such as alcohol intake (yes/no), smoking (yes/no/sometimes), active in sports (organised/not organised/no sports) and employment (yes/no/rest/in education) were also collected by questionnaires.
The missing values were completed by the respective median of the total sample (pregnancy duration n 2; pregnancy weight gain n 10; birth weight n 1; maternal overweight n 4; adolescents’ BMI–standard deviation score, fat mass index and fat-free mass index n 2; sleeping duration in adolescence n 4; fat mass index in adulthood n 2; smoking status n 64; alcohol intake n 23; and physical activity n 40). The puberty status ‘age of take-off’, defined as the age at minimal height velocity at the onset of the pubertal growth spurt(Reference Buyken, Karaolis-Danckert and Remer36), and ‘employment during adulthood’ could not be considered because of too many missing data (>20 %).
Statistical analyses
All statistical analyses were performed using SAS® procedures (version 9.20 and 9.40). The significance level was set at P < 0·05. Descriptive data are presented as median, with interquartile ranges for continuous variables and frequencies and percentages for categorical variables (Tables 1 and 2).
Multivariable linear regression was used to analyse the prospective associations between dairy intake and biomarkers of low-grade systemic inflammation or insulin resistance. All associations with inflammation were analysed for each biomarker of low-grade systemic inflammation and the pro-inflammatory score. Compared with similar studies(Reference Diederichs, Herder and Roßbach24,Reference Perrar, Buyken and Penczynski37) , results from the regression analyses are presented as adjusted least-square means (95 % CI) by sex-specific tertiles (low, medium and high intake) of the respective predictor, with P trend values from models with the predictors as continuous variables.
Previously, individual outliers of every biomarker that significantly affected the normal distribution or regression modelling were winsorised, that is, the outliers were replaced by the closest sex-specific value corresponding to a normal distribution. The procedure involved IL-6 (1 %) and adiponectin (<1 %).
No stratification by sex was conducted as no significant interactions between dairy intake and sex for primary outcomes were observed in the analysis.
All basic models (model 1) included the exposure variable (intake of TD or each dairy type separately) during adolescence, sex and age at blood withdrawal in adulthood. For adjusted models (model 2), potential confounders were included individually and hierarchically if they significantly affected the regression coefficient of exposure by ≥10 % or predicted the outcome variable independently(Reference Maldonado and Greenland38,Reference Victora, Huttly and Fuchs39) . To ensure comparability, all pro-inflammatory score models (model 2) were identically adjusted. The same applies to the HOMA2-IR models. On the basis of the hierarchical examination of possible confounders, all adjusted models for all outcomes (model 2) included BMI in adulthood only.
To reduce potential bias, sensitivity analyses were conducted by excluding under-reported records. Dietary records were classified as ‘under-reported’ if the relationship between total energy intake and estimated BMR according to age- and sex-specific equations by Schofield(Reference Schofield40) was not plausible. Under-reported records were identified using the paediatric cut-offs by Sichert-Hellert(Reference Sichert-Hellert, Kersting and Schöch41). This calculation resulted in the records of ten participants being identified as under-reported (2·7 %).
To exclude any possible bias of a short gestation period (<37 weeks) or low birth weight (<2500 g) on the outcome, further sensitivity analyses were performed wherein data from relevant participants were excluded (in total, n 20 participants, 5·3 %). In addition, sensitivity analyses were conducted by excluding participants with no data on smoking status (n 64, 17·0 %) and those with a known intolerance against dairy (n 13, 3·5 %). Because we could not determine whether the intolerance was a sensitivity or an allergy to dairy products, all thirteen participants with a general intolerance were excluded.
A further sensitivity analysis was conducted to consider the possible effect of dairy intake changes during adolescence(Reference Hohoff, Perrar and Jancovic22) that could have masked the effects of long-term dairy intake on the different outcomes. For this purpose, dietary data from early adolescence (boys, 9–12·5 years, and girls, 8–11·5 years) and late adolescence were considered in stratified analyses (boys, 12·5–16 years, and girls, 11·5–15 years).
Results
Sample characteristics in adolescence and young adulthood of the participants (n 375) are shown in Tables 1 and 2. Approximately half of the participants were male (46·4 %). Participants and their mothers’ BMI values were within the normal range. Maternal employment amounted to >66 %. The median age at adolescence was 12 years (Table 1), and the median age at blood withdrawal was 21 years (Table 1). The median follow-up period between dietary records in adolescents and blood withdrawal in adulthood was 9·2 years. The daily median TD intake was 177·9 g/1000 kcal. Participants consumed more non-fermented than fermented, more low-sugar than high-sugar and more high-fat than low-fat dairy products. Non-fermented dairy had the highest median intake values in all tertiles (Table 4).
DONALD, DOrtmund Nutritional and Anthropometric Longitudinally Designed.
* Values refer to median (25th, 75th percentiles) of intake in the respective sex-specific tertile.
Overall, no association was observed between TD intake and the pro-inflammatory score. Similarly, no dairy type was associated with the pro-inflammatory score (Table 5). Grouping dairy intake into tertiles did not show any associations. The same applies to the examination of individual inflammation biomarkers (see online supplementary material, Supplementary Tables a–c).
Model values are least-square means (95 % CI) for tertiles obtained from linear regression models.
* P trend values are based on models using the continuous exposure variables.
† Adjusted for sex and age at blood withdrawal.
‡ Adjusted for sex, age at blood withdrawal as well as BMI in adulthood.
In addition, the analysis of TD intake and each dairy-type intake and insulin resistance showed no association (Table 6). No association was found when grouping dairy intake into tertiles.
HOMA2-IR, Homeostasis Model Assessment Insulin Resistance.
Model values are least-square means (95 % CI) for tertiles obtained from linear regression models. Outcome variables were log-transformed.
* P trend values are based on models using the continuous exposure variables.
† Adjusted for sex and age at blood withdrawal.
‡ Adjusted for sex, age at blood withdrawal as well as BMI in adulthood.
§ HOMA2-IR.
Sensitivity analyses excluding under-reported records confirmed our null results for both the pro-inflammatory score and insulin resistance (data not shown). Our results also remained the same after the exclusion of participants with a low birth weight or a short duration of pregnancy, participants without information on smoking status or participants with intolerance to dairy (data not shown). Stratifying the data into records from early adolescence and records from late adolescence also had no consequence and confirmed the reliability of the results for both outcomes.
Discussion
This study examined the long-term association of habitual dairy intake in healthy adolescents with inflammation biomarkers and insulin resistance in young adulthood. The daily median TD intake of 177·78 g/1000 kcal in our sample is approximately 30 % below that is specified by the German Food-Based Dietary Guidelines for adolescents(Reference Kersting, Kalhoff and Lücke42). This is in accordance with the representative EsKiMo II study, which analysed the dietary behaviour of 2644 children and adolescents in Germany between 2015 and 2017. In this study, approximately half of German children and adolescents aged 6–17 years did not achieve the recommended dairy intake(Reference Mensink, Haftenberger and Lage Barbosa43). We did not find any association between dairy or different types of dairy intake and pro-inflammatory score or insulin resistance. Accordingly, our study shows that missing associations observed in previous cross-sectional studies in healthy children can also be confirmed in the long term, considering adolescent dairy intake and blood parameters in young adulthood(Reference Aslam, Jacka and Marx18,Reference Arouca, Michels and Moreno44,Reference González-Gil, Santabárbara and Russo45) . To the best of our knowledge, this is the first study to examine different types of dairy in this context.
The absence of any association with the pro-inflammatory score suggests that dairy intake in adolescence is not of major longer-term relevance for low-grade systemic inflammation among young adults. Our findings are in line with those of a systematic review by Bujtor et al. (Reference Bujtor, Turner and Torres19) on the associations of dietary intake with single biomarkers of inflammation (C-reactive protein or IL-6) in healthy children and adolescents. However, in this review, only seven studies examined dairy intake, of which only five included healthy participants. A cross-sectional study by Abreu et al. (Reference Abreu, Agostinis-Sobrinho and Santos46) suggested an inverse association between TD or milk intake and serum IL-6 concentrations among normal-weight adolescents; moreover, adolescents in the second tertile of yoghurt intake showed IL-6 lower levels than that of those in the first tertile. In our analysis, we did not find a difference between intake tertiles for the pro-inflammatory score and any biomarker of inflammation. However, we considered fermented dairy in general terms and not yoghurt intake specifically. Furthermore, we examined only the long-term associations and not the short-term ones.
In a systematic review, Bordoni et al. (Reference Bordoni, Danesi and Dardevet14) developed a pro-inflammatory score from various individual biomarkers, comparable to the inflammation score used in the present analysis. The authors summarised the anti-inflammatory activity of dairy intake in adults with metabolic disorders and the pro-inflammatory activity in adults allergic to cow’s milk. A stratified consideration of different types of dairy in this analysis indicated a weak anti-inflammatory activity of low-fat and high-fat as well as fermented dairy in participants who were not allergic to dairy products(Reference Bordoni, Danesi and Dardevet14). Our analyses did not confirm this anti-inflammatory association even after excluding participants with a general intolerance to dairy (data not shown). However, in our sample, only six participants (3·5 %) reported such intolerance. In addition, the analyses by Bordoni et al. (Reference Bordoni, Danesi and Dardevet14) were also based on analyses among adults. Further research in this context is required. A narrative review(Reference Da Silva and Rudkowska47) resumed that most SFA activate pro-inflammatory biomarkers. Which combination of nutrients results in a neutral or inverse association between dairy intake and inflammation remains to be clarified.
With regard to the relationship between dairy intake during adolescence and insulin resistance in young adulthood, our results are not in line with those of the analysis of the ‘Nurses’ Health Study II’ cohort by Malik et al. (Reference Malik, Sun and van Dam48), in which higher dairy intake during adolescence was associated with a lower risk of type 2 diabetes in adulthood. Our analyses did not find an association between TD intake and dairy-type intake and insulin resistance. The different results can probably be attributed to different methods. Malik et al. (Reference Malik, Sun and van Dam48) investigated inflammatory biomarkers exclusively in middle-aged women, whereas our study took both sexes into account. However, tests for the interaction of inflammatory markers with sex in our sample indicated no difference. In addition, the risk of inflammation among nurses may be increased because of shift work(Reference Skogstad, Aass and Sirnes49). Additionally, Malik et al. (Reference Malik, Sun and van Dam48) applied a FFQ to assess participants’ diet during high school, which carries the risk of recall bias. However, the sample size in Malik et al. (Reference Malik, Sun and van Dam48) was much larger than that in our study.
Our results also did not confirm the findings reported in a systematic review and meta-analysis of prospective studies(Reference Schwingshackl, Hoffmann and Lampousi50) or in a systematic review and meta-analysis of randomised clinical trials(Reference Sochol, Johns and Buttar51) that suggested a beneficial effect of low-fat dairy on HOMA2-IR. They assumed that Ca, vitamin D, casein and whey proteins in low-fat dairy are potential regulators of body fat, waist circumference and insulin resistance. However, we could not find any associations independent of the type of dairy. High levels of adipose tissue were assumed to result in the development of insulin resistance(Reference Frayn, Tan and Karpe52). In a previous study based on our sample, higher TD intake led to a more favourable body composition in the long term; however, higher intake of low-fat dairy showed no association(Reference Hohoff, Perrar and Jankovic23). However, our sample is characterised by a rather low BMI compared with the general population in Germany(Reference Schienkiewitz, Kuhnert and Blume53,Reference Schienkiewitz, Damerow and Rosario54) .
Some strengths and limitations of this analysis of the DONALD study must be discussed. The prospective design and repeatedly collected detailed dietary measurements allow the investigation of long-term associations between adolescent habitual dairy intake including subtypes and adult health outcomes(Reference Perrar, Alexy and Nöthlings20). The continuously updated in-house nutrient database LEBTAB allows the consideration of different types of dairy according to composition and processing methods(Reference Sichert-Hellert, Kersting and Chahda21). The main limitation of our sample is the over-representation of families with a high socio-economic background in the DONALD study, which limits the generalisability of our results(Reference Perrar, Alexy and Nöthlings20). Although our sample was not very large, and relatively young and healthy, previous analyses showed that significant associations, for example, between dietary intake and risk factors for type 2 diabetes, could still be shown with the data from the DONALD study(Reference Penczynski, Herder and Krupp26).
Food grouping in our study is both a strength and a limitation. The dietary assessment method and food composition database allow the aggregation of several diverse food groups. However, some of these subgroups overlap to allow the investigation of possible associations with health (e.g. fermentation and fat content). A finer subdivision would have led to very low intake quantities. The underlying classification has already been used in our previous publications(Reference Hohoff, Perrar and Jancovic22,Reference Hohoff, Perrar and Jankovic23) and allows an overall interpretation of the results.
In addition, we cannot reject the possibility of under-reporting. Under-reported records were not generally excluded from the main analyses because this method only identifies under-reported energy intake and unselective under-reporting of individual foods(Reference Livingstone and Robson55). However, our sensitivity analyses, excluding energy under-reports, showed similar results to those of the main analyses.
Conclusions
Our results indicated the absence of any associations between dairy intake and pro-inflammatory score or insulin resistance in young adulthood. Thus, the habitual intake of individual types of dairy during childhood and adolescence does not influence these metabolic risk factors in the long term. Restrictions on dairy intake for healthy children and adolescents appear redundant in terms of diabetes risk reduction.
Acknowledgements
The participation of all children and their families in the DONALD study is gratefully acknowledged. The authors also thank the DONALD staff for carrying out the anthropometric measurements, for administering the questionnaires, and for collecting and coding the dietary records.
Financial support
The DONALD study is financially supported by the Ministry of Science and Research of North Rhine-Westphalia, Germany. The German Diabetes Center is funded by the German Federal Ministry of Health (Berlin, Germany) and the Ministry of Culture and Science of the state North Rhine-Westphalia (Düsseldorf, Germany) and receives additional funding from the German Federal Ministry of Education and Research (BMBF) through the German Center for Diabetes Research (DZD e.V.).
Conflict of interest
All authors declare no conflicts of interest. The funding sponsor had no role: in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.
Authorship
U.A. conceived the project; E.H. performed the statistical analysis and wrote the manuscript; U.A. supervised the study; U.A. and N.J. provided extensive editorial support; C.H. supervised the laboratory measurements of the blood analyses. All authors made substantial contributions and read and approved the final manuscript.
Ethics of human subject participation
The DONALD study is approved by the Ethics Committee of the University of Bonn (ethics applications: 098/06 and 185/20) according to the guidelines of the Declaration of Helsinki. All examinations are performed with written parental and adult participants’ consent.
Supplementary material
For supplementary material accompanying this paper, visit https://doi.org/10.1017/S1368980024000624.