Se is a trace element that has been reported to play an important role in many areas of the human body including thyroid metabolism, cognitive function, immune function and growth(Reference Pillai, Uyehara-Lock and Bellinger1–Reference Rayman3). While there are many functions of Se in the human body, there appears to be a U-shaped curve of adverse effects, where deficits and excessive Se concentrations can impact on the risk and onset of certain diseases(Reference Hu, Zhao and Hu4,Reference Barchielli, Capperucci and Tanini5) .
It is well known that groups of the New Zealand population, particularly in the South Island, are not meeting the recommended intakes of Se and have lower status than what has been reported in some other countries(Reference Thomson6,Reference McLachlan, Thomson and Ferguson7) . This is in most part due to low Se concentrations in the soil(Reference Thomson6,Reference Gupta and Gupta8) . However, only one study has reported information on Se intakes and status of South Island (New Zealand) infants and toddlers, in 1998–1999(Reference McLachlan, Thomson and Ferguson7). Also, it is not known whether the increasingly popular alternative approach to introducing solid foods to infants known as Baby-Led Weaning (in which infants feed themselves rather than being spoon-fed) influences Se intakes. A baby-led approach to infant feeding has been shown to influence food and nutrient intake in infancy(Reference Williams Erickson, Taylor and Haszard9–Reference Rowan, Lee and Brown11) so may influence Se intake and therefore status.
The objectives of this study were to determine: (1) whether Se intakes and plasma Se concentrations differ between 12-month-old New Zealand toddlers who followed a baby-led approach to complementary feeding and those who followed traditional spoon-feeding, (2) what food sources contribute to Se intakes of 12-month-old toddlers and (3) what factors are associated with plasma Se concentrations in New Zealand 12-month-old toddlers.
Experimental methods
This is a secondary analysis of the Baby-Led Introduction to SolidS (BLISS) study which has been described in full elsewhere(Reference Daniels, Heath and Williams12). The information provided here is information of relevance to the current analysis. The BLISS randomised controlled trial investigated the impact of a modified version of Baby-Led Weaning on infant outcomes including growth(Reference Taylor, Williams and Fangupo13), choking(Reference Fangupo, Heath and Williams14), Fe(Reference Daniels, Taylor and Williams15) and Zn(Reference Daniels, Taylor and Williams16). Women (n 206) were recruited (between November 2012 and March 2014) during the third trimester of their pregnancy through an opt-out enrolment process (i.e. all potentially eligible women were approached and offered the opportunity to participate) at the only birthing facility in the city of Dunedin, New Zealand (southern region of New Zealand). Participants were eligible if they: spoke English or Te Reo Māori (indigenous language of New Zealand); planned to live in the area of Dunedin, New Zealand, until their child was at least 2 years of age and were 16 years of age or older. Exclusion criteria were if the infant was born before 37 weeks gestation or had a congenital abnormality, physical condition or intellectual disability that was likely to affect their feeding or growth. The study was conducted in accordance with the Declaration of Helsinki for research involving human subjects and was approved by the New Zealand Lower South Regional Ethics Committee (LRS/11/09/037). Adult participants gave written informed consent.
After obtaining consent, participants were randomised into the Control or BLISS intervention group (Fig. 1). Both groups received standard ‘Well Child’ care available to all New Zealand children up to 5 years of age(17), and BLISS participants received further education and support regarding the BLISS approach (i.e. infant self-feeding, with no spoon-feeding, from 6 months of age with modifications to address concerns about the possible risks of Fe deficiency, choking and growth faltering(Reference Cameron, Heath and Taylor18)).
Demographic data were collected by questionnaire, including the participant’s address which was used to determine the level of household deprivation with the New Zealand Index of Deprivation (NZDep) score(Reference Atkinson, Salmond and Crampton19). Infant sex, birth weight and gestational age at birth were obtained from hospital records. When the child was 12 months of age, parents were asked about the mother’s smoking status during pregnancy (daily, occasional, quit during pregnancy, non-smoker).
Parent participants completed a weighed diet record on three randomly assigned non-consecutive days (two week days and one weekend day), using dietary scales (Salter Electronic, Salter Housewares Ltd.) accurate to ±1 g, when the child was aged 12 months. Parents recorded all information on what their child ate and drank (time of day the food was consumed, type and brand of food, cooking method), the total food weight before offering and at the end of the eating occasion (i.e. that was leftover), any dietary supplements taken and any recipes used. Data from the 3-d weighed diet record were analysed using the database Kai-culator (Version 1.13s, University of Otago) which contains data from the New Zealand Food Composition Database(20), recipes from the 2008/2009 New Zealand Adult Nutrition Survey(21), and commercial infant foods collated by the research team(Reference Clouston22). During the dietary data entry, food items were selected carefully to ensure that they contained appropriate Se concentrations, notably for flour and bread products because of the known regional variations in their Se content(Reference Thomson6). For toddlers consuming breast milk (no infant formula), intakes were estimated to be 448 g/d(Reference Dewey, Heinig and Nommsen23). If the toddler was mixed fed (breast milk and infant formula), the amount of breast milk was estimated to be the amount left after subtracting their weighed consumption of infant formula (g per day) from 448 g/d. The Se content of breast milk was assumed to be 2 μg/100 g(20). The data from the 3-d weighed diet record were used to determine daily intakes of: energy, protein, Se and Se intakes from breast milk, infant formula and nine food groups (‘seafood’, ‘breads and cereals’, ‘meat’, ‘dairy’, ‘fruits and vegetables’, ‘eggs’, ‘nuts and seeds’, ‘legumes’ and ‘miscellaneous’). Intakes of energy, protein and Se were then entered into the Multiple Source Method programme(Reference Harttig, Haubrock and Knüppel24) to calculate ‘usual’ daily intakes, which accounts for intra-individual variation in intake.
A non-fasting venous blood sample was obtained at 12 months of age to determine plasma Se, C-reactive protein and α1-acid glycoprotein concentrations, as well as for determination of biomarkers of Fe and Zn status, as has been reported earlier(Reference Daniels, Taylor and Williams16,17) . For the purposes of Zn analysis, parents were instructed to give their child a milk feed 90 min prior to blood collection and then no other food or drink until after the blood test, and a rigorous trace-element free protocol was used during both blood collection and analysis (including use of trace-element free lithium heparin anticoagulated S-Monovette tubes). After sample collection, plasma was separated within 2 h (3500 rpm for 5 min) and aliquots were stored at –80°C until analysis. Plasma Se was analysed using inductively coupled plasma mass spectrometry (ICP-MS) at the Centre for Trace Element Analysis, Department of Chemistry, University of Otago, Dunedin, New Zealand. C-reactive protein and α1-acid glycoprotein were analysed using a Cobas C311 automatic electronic analyser (Roche), in the Department of Human Nutrition Laboratories (University of Otago).
The accuracy and precision of the analyses were checked using certified controls and in-house pooled samples (after every 15 samples), respectively. The analysed mean ± sd (CV) value for the Se control (UTAK Laboratories, Inc.) was 1·32 ± 0·02 μmol/l (1·3 %), compared with the manufacturer’s concentration of 1·37 μmol/l. The mean ± sd (CV) for the C-reactive protein control (Roche Diagnostics) was 9·5 ± 0·4 mg/l (4·6 %), compared with the manufacturer’s concentration of 9·1 mg/l. The multilevel controls for α1-acid glycoprotein (Roche Diagnostics) were 0·5 ± 0·01 g/l (1·1 %) and 0·8 ± 0·01 g/l (1·4 %), compared with the manufacturer’s concentrations of 0·7 and 1·2 g/l, respectively.
Statistical analysis
Plasma Se concentrations were adjusted for inflammation using the Biomarkers Reflecting Inflammation and Nutrition Determinants of Anaemia (BRINDA) regression approach(Reference Suchdev, Namaste and Aaron25,Reference Larson, Addo and Sandalinas26) calculated for each participant as: Se adjusted = exp (ln Se – (β1 × ln C-reactive protein) – (β2 × ln α1-acid glycoprotein)). Differences between groups were estimated using a linear regression model adjusted for parity (1 child v. > 1 child) and maternal education (non-tertiary v. tertiary).
Usual Se intakes were used to calculate the number of participants with an intake below the estimated average requirement (EAR) of 20 μg/d for 1–3 year olds(27,28) . Logistic regression was used to determine whether the odds of having Se intakes below the EAR were different between groups after adjustment for parity and maternal education.
An analysis of the Control group was used to determine the contribution of breast milk, infant formula and food groups towards daily Se intakes. The percentage contribution of each food group to Se intakes for consumers, and for the total sample, was calculated as medians and 25th and 75th percentiles (inter-quartile range). Only the Control group was included in this analysis because the BLISS intervention changed eating behaviour(Reference Taylor, Williams and Fangupo13).
Univariate unadjusted linear regression analyses were used to describe associations between plasma Se (adjusted for inflammation) and potential predictor variables. These variables were decided a priori based on previous associations described in the literature(Reference McLachlan, Thomson and Ferguson7,Reference Bates, Thane and Prentice29–Reference Mengubas, Diab and Gokmen32) or were considered to be potentially associated with plasma Se concentrations for mechanistic reasons. Maternal predictor variables were household deprivation, maternal education, maternal age, smoking status during pregnancy; child predictor variables were ethnicity, sex, weight gain and linear growth between 6 and 12 months, plasma Zn, Fe deficiency anaemia, body Fe, intakes of energy, protein, Se, ‘breads and cereals’, ‘meat’, ‘dairy’, and ‘fruits and vegetables’ and consumers of breast milk, infant formula and ‘seafood’. Prior to the regression analysis, smoking status during pregnancy was collapsed into two categories: none (n 89) and any (n 8; which included daily smokers (n 2), occasional smokers (n 2) and those who quit during pregnancy (n 4)). Plasma Zn concentration was adjusted for time of blood sampling (variable adjusted for reference time of 08.00 hours) and time since last meal (variable adjusted for reference period of 90 min)(Reference Arsenault, Wuehler and de Romana33) using the regression equation: Zn time adjusted = exp (ln plasma Zn – (β1 × adj_08.00 hours) – (β2 × adj_90 min)). Following this, an adjustment was also made to plasma Zn concentrations for inflammation using the same approach as for Se explained above. Body Fe was calculated in mg/kg using the equation: –(log10 (soluble transferrin receptor × 1000/plasma ferritin)–2·8229)/0·1207(Reference Cogswell, Looker and Pfeiffer34) and a body Fe concentration < 0 mg/kg and a Hb concentration < 110 g/l was used to define Fe deficiency anaemia(Reference Daniels, Taylor and Williams15). All continuous variables were standardised to allow for comparison of the regression coefficients; the dietary variables were calculated as daily intakes in grams. In the adjusted multivariate regression analysis, variables were chosen to be included in the model if the regression coefficient was ≥ 0·06 μmol/l per unit (this number is approximately 0·3 sd, commonly referred to as a ‘small effect’(Reference Cohen35)) or with a P < 0·1, and were adjusted for infant sex.
All analyses were conducted using Stata, version 15.1 (StataCorp LP). P values < 0·05 were considered to be statistically significant.
Results
A total of 101 participants (n 51 Control and n 50 BLISS) contributed both dietary intake data by 3-d weighed diet record and a plasma Se sample (Fig. 1) at 12 months of age (50 % of initial sample). Although only twenty-two participants (10·6 %) had formally withdrawn from the study at this stage, a number of participants in both groups did not provide dietary (n 11), biochemical (n 41) or both (n 53) types of data making them ineligible for these analyses. Maternal, household and infant characteristics of those included in this analysis are shown in Table 1. There was no evidence of differences between those included and excluded from this analysis (all P > 0·05) except for maternal education (those included were more likely to have a university education: 57 % compared with 40 %, P = 0·044) and maternal age (those included were 2·4 years older on average, P = 0·002).
* Household deprivation categorised into: 1–3 (low), 4–7 and 8–10 (high) using the New Zealand Index of Deprivation 2013. The index combines different dimensions of deprivation from New Zealand census data. A deprivation score is assigned to each meshblock (geographical area defined by Statistics New Zealand)(Reference Atkinson, Salmond and Crampton19).
† Available data for Control n 50 and BLISS n 48.
Selenium intake and status
Usual Se intakes are shown in Table 2. The OR of Se intakes below the EAR of 20 μg/d(27,28) was no different between the two groups (OR: 0·89; 95 % CI 0·39, 2·03) (Table 2). The highest usual Se intake was 37 μg/d. No dietary supplements containing Se were taken. Plasma Se concentrations (adjusted for inflammation) are also shown in Table 2. Plasma Se concentrations without adjustment for inflammation are shown in online Supplementary Table 1.
EAR, estimated average requirement; AGP, α1-acid glycoprotein; CRP, C-reactive protein.
* Mean differences and 95 % confidence intervals for BLISS compared with Controls (adjusted for maternal education and parity).
† Intake reported in the 3-d weighed diet records collected at 12 months of age, and usual dietary intakes calculated using the Multiple Source Method(Reference Harttig, Haubrock and Knüppel24).
§ Adjusted for inflammation using the BRINDA(Reference Suchdev, Namaste and Aaron25) approach: exp (unadjusted ln plasma Se−(regression coefficient for CRP) × (CRP−(maximum of lowest decile for CRP))−(regression coefficient for AGP) × (AGP−(maximum of lowest decile for AGP))).
Food groups contributing to selenium intakes
The food group contributing the most to total Se intakes of 12-month-old toddlers appeared to be breast milk, contributing 20 % in the whole sample and 39 % when restricted to those who consumed breast milk (Table 3). In total, ‘breads and cereals’ contributed the most Se (12 %), followed by ‘meat’ (11 %) and ‘dairy’ (10 %). When data for ‘consumers’ only were analysed both ‘breads and cereals’ (12 %) and ‘meat’ (12 %) equally contributed to Se intakes. Of those consuming infant formula (the infant formulas consumed in this study contained between 0 and 17·3 μg Se per 100 g), formula intakes contributed 9 % of daily Se intakes (Table 3).
IQR, inter-quartile range.
* Only Control group included (n 51).
† Intake reported in the weighed 3-d diet records collected at 12 months of age.
‡ Ordered by the food group contributing the most to Se intakes by consumers (‘Consumers’: the contribution of Se intakes of the toddlers consuming this food group (e.g. breast milk contributes a median of 39 % of Se intakes for the twenty-eight toddlers consuming breast milk)).
§ Data expressed as median percentages (NB: mean percentages added to 100 % of total Se intakes from food groups).
|| Comprises all meat: red meat, poultry, pork, processed meat, etc.
¶ Includes cows’ milk as a drink.
** Comprises fish and shellfish.
†† Miscellaneous comprises: fats, sugar, sweet foods, herbs and spices, sauces, spreads, beverages, etc.
Factors associated with plasma selenium concentrations
Univariate and multivariate associations between factors decided a priori and plasma Se concentrations (adjusted for inflammation) are shown in Table 4. In the unadjusted analysis, toddlers of mothers who reported any smoking (daily, occasional, quit during) during pregnancy had on average 0·16 μmol/l lower plasma Se concentrations compared with toddlers of mothers who did not smoke during pregnancy, although this association was attenuated in the adjusted model (0·13 μmol/l; 95 % CI −0·25, −0·003). ‘Seafood’ intake in toddlers was associated with higher plasma Se concentrations in the unadjusted analysis (0·12 μmol/l; 95 % CI 0·04, 0·20), but this association was also attenuated in the adjusted analysis. Consuming breast milk was associated with 0·12 μmol/l (95 % CI −0·19, −0·04) lower plasma Se concentrations in toddlers in the adjusted model. The R 2 for the final multivariate model was 0·30 indicating that 30 % of the variance in plasma Se concentration was explained by the factors included in this model.
* Variables chosen to be included in the final model were those with regression coefficient ≥ 0·06 μmol/l or P < 0·1, adjusted for child sex.
† Household deprivation categorised into: 1–3 (low), 4–7 and 8–10 (high) using the New Zealand Index of Deprivation 2013. The index combines different dimensions of deprivation from New Zealand census data. A deprivation score is assigned to each meshblock (geographical area defined by Statistics New Zealand)(Reference Atkinson, Salmond and Crampton19).
‡ Female (n 56) compared with male (n 45).
§ Data available for n 97; any smoking (n 8) compared with no smoking (n 89).
|| Plasma Zn adjusted for time of blood sampling and time since last meal(Reference Arsenault, Wuehler and de Romana33) and inflammation(Reference Suchdev, Namaste and Aaron25).
¶ Standardised continuous variable; dietary intake variables as daily intakes in grams.
** Data available for n 100; yes (n 19) compared with no (n 81), Fe deficiency anaemia defined as Hb < 110 g/l and body Fe < 0 mg/kg.
†† Body Fe calculation (mg/kg): –(log10(sTfR × 1000/ferritin)−2·8229)/0·1207 from Cogswell et al. (Reference Cogswell, Looker and Pfeiffer34)
‡‡ Usual intakes: calculated using the Multiple Source Method(Reference Harttig, Haubrock and Knüppel24) from 3-d weighed diet records.
§§ Yes (n 56) compared with no (n 45).
|||| Yes (n 42) compared with no (n 59).
¶¶ Yes (n 24) compared with no (n 77).
*** Data available for n 97.
Discussion
The current results demonstrate that amongst New Zealand toddlers (12 months of age) many have Se intakes that are below the EAR, this is regardless of how complementary foods are introduced to them as infants, and there was no difference in plasma Se concentrations between toddlers who had followed a baby-led approach to complementary feeding and those who followed traditional spoon-feeding. Of the infant milk sources (breast milk and infant formula), breast milk appeared to contribute the most to dietary Se intakes of consumers in this age group. For those who consumed them, ‘breads and cereals’ and ‘meat’ contributed the most Se of all the food groups, followed by ‘dairy’. In the adjusted analysis, consumption of breast milk was negatively associated with plasma Se concentrations in toddlers. Maternal smoking during pregnancy had a weak negative association, and consumption of ‘seafood’ had a weak positive association, with plasma Se concentrations.
While very few studies have assessed Se intakes in young New Zealand children, our mean Se intakes of 18·2 μg/d for Control and 19·2 μg/d for BLISS toddlers were both higher than previously reported intakes of 13·7 μg/d for New Zealand toddlers (12–24 months) collected two decades earlier(Reference McLachlan, Thomson and Ferguson7). However, the exclusion of breastfed toddlers from the study by McLachlan et al. (Reference McLachlan, Thomson and Ferguson7) may contribute to some of the discrepancy here, as well as the fact that the earlier New Zealand Food Composition Database may not have sufficiently estimated regional variations in the Se content of foods. Our results suggest that more than half (63 % Control, 60 % BLISS) of the toddlers in this study had Se intakes below the EAR; however, it has been recently suggested that there is a great need for change to the dietary reference values for Se (for all age groups), given that the recommended intake values were set when there was scarce evidence on the health effects of Se(Reference Vinceti, Filippini and Jablonska36).
Se is considered to be highly toxic at the upper level of intake, or UL (usually as a result of high intakes of Se supplements), with reported symptoms of hair, skin and gastrointestinal abnormalities, and fatigue(Reference MacFarquhar, Broussard and Melstrom37). While no toddlers in this study consumed Se supplements, and none had total intakes above the upper level of intake (90 μg/d(28)), evidence in adults suggests that long-term exposure at intakes lower than the current adult UL may increase the risk of type-2-diabetes(Reference Stranges, Marshall and Natarajan38) and this has led to debate and a call for an up to date risk assessment on the adverse effects levels for all ages(Reference Vinceti, Filippini and Cilloni39,Reference Vinceti, Filippini and Wise40) .
The main (non-infant milk) food group contributing to Se intakes of the total study sample was ‘breads and cereals’, consistent with previous findings in young New Zealand children(Reference McLachlan, Thomson and Ferguson7,Reference Thomson, McLachlan and Parnell31) . While offal is a good source of Se(41,Reference Combs42) , only one participant consumed any (chicken liver pâté) at 12 months of age. The infant formulas consumed had a wide range of Se concentrations (between 0 and 17·3 μg/100 g), but interestingly, none of the ‘toddler milks’ (marketed for toddlers > 12 months of age) contained added Se. The contribution of Se intakes from different infant formulas (with varying Se concentrations) was not assessed in this study.
The finding that breast milk consumption appeared to be the largest contributor to Se intake was surprising given that consuming breast milk was associated with poorer Se status. There are a number of possible explanations for this and this finding should be treated with caution. First, it was not possible to measure breast milk intake, and the single value used for breast milk intake(Reference Dewey, Heinig and Nommsen23), while used in other studies, was not generated using New Zealand data so may not accurately reflect actual breast milk intake in these toddlers. Second, the Se content of breast milk is dependent on the Se status of the mother(43), but we were not able to directly measure the Se concentration of breast milk consumed by toddlers in this study. Instead, a Se value of 2 μg/100 g taken from the New Zealand Food Composition Database(20) was used, which may be too high, especially for South Island mothers and their infants. Clearly, analysis of infant breast milk intake and of the Se concentrations of breast milk from mothers in the South Island of New Zealand is urgently required.
The mean plasma Se concentrations reported in this study were 0·81 μmol/l for Control and 0·85 μmol/l for BLISS infants, which were higher than has been previously reported in infants (0·69 μmol/l), toddlers (0·61 μmol/l) and young children (0·79 μmol/l) in the South Island of New Zealand(Reference McLachlan, Thomson and Ferguson7,Reference Thomson, McLachlan and Parnell31) . The results are also at the upper end of the range for infants and young children reported in studies internationally: 0·67–0·83 μmol/l(Reference Bates, Thane and Prentice29,Reference Gurgoze, Olcucu and Aygun30,Reference Jardim-Botelho, Gurgel and Henriques44) , which is surprising given the lower Se concentrations reported in this region. There are currently no universal (age appropriate) criteria for low plasma Se due to the wide variability in the soil Se concentrations of different geographical locations(Reference Combs42), and the paucity of data on possible health effects. Earlier research suggested a cut-off of ≤ 0·82 μmol/l as a level considered to ensure optimal activity of enzymes (iodothyronine 5’ deiodinases) associated with thyroid function(Reference Thomson45); however, there is little evidence that concentrations above this cut-off are beneficial to all areas of human health(Reference Vinceti, Filippini and Jablonska36). The current study did not assess any health effects related to the plasma Se concentrations reported, and long-term implications of low Se status in this age group remain relatively unclear, as does the specific level/cut-off for defining low plasma concentrations.
For all toddlers, maternal smoking status during pregnancy and breast milk consumption were negatively associated with plasma Se concentrations after adjustment for child sex and other potential predictor variables in the model. Lower Se status of young children has previously been shown to be associated with parental(Reference Bates, Thane and Prentice29) and household smoking(Reference McLachlan, Thomson and Ferguson7). Se intakes below recommendations have been reported in the New Zealand population(Reference Thomson6,Reference McLachlan, Thomson and Ferguson7) and low maternal Se status reduces breast milk Se concentrations(43) which could be contributing to the association between breast milk consumption and toddlers’ plasma Se concentrations.
This is the first study to assess whether the baby-led or spoon-feeding approach to complementary feeding has an impact on Se status in toddlerhood. Rigorous methods were used for the dietary analysis to ensure regional variations in food Se concentrations were taken into account, an approach not possible previously(Reference Thomson, McLachlan and Parnell31), and to adjust the distribution of observed intakes to estimate usual Se intakes. However, our study also has some limitations. This was a secondary analysis of the BLISS study which was not specifically designed or powered to assess Se intake and status. Se intake from breast milk was estimated, because the volume of breast milk and concentration of Se in breast milk were not measured. There were a large number of participants excluded from this analysis because of incomplete data. Lastly, the study was not designed to determine the specific health effects related to the Se intakes and plasma concentrations reported and further studies are needed to determine what, or whether any negative health consequences at these plasma Se concentrations occur.
In conclusion, Se intakes and plasma Se concentrations of 12-month-old New Zealand toddlers were no different between those who had followed a baby-led approach to complementary feeding and those who followed traditional spoon-feeding. However, more than half of toddlers had Se intakes below the EAR. Further research is required to determine whether any negative health consequences at these intakes and status occur.
Acknowledgements
The authors would like to acknowledge all the families who participated in the Baby-Led Introduction to SolidS (BLISS) study, and the BLISS research staff in the Departments of Medicine and Human Nutrition at the University of Otago, Dunedin, New Zealand.
The study was supported by Meat & Livestock Australia, Lottery Health Research, Karitane Products Society, Perpetual Trustees, New Zealand Federation of Women’s Institutes and the University of Otago. The funding sponsors had no role in the design of the study; in the collection, analysis or interpretation of data; in the writing of the manuscript or in the decision to publish the results. Heinz Wattie’s provided baby rice cereal for participants in the intervention group. Part of R.W.T.’s salary is from a Karitane Products Society Fellowship that is paid via the University of Otago. L. D. was supported by a University of Otago Doctoral Scholarship.
L. D. collected data and prepared the first full and subsequent drafts of this manuscript. A.-L. M. H. and R. W. T. are the co-principal investigators of the BLISS study, designed the research project and made a considerable intellectual contribution to the writing of this manuscript. R. S. G. and L. A. F. advised on aspects of the study design and ongoing advice and support. J. J. H. and L. D. performed the statistical analyses. R. S. G., J. C. M. and C. D. T. provided expert input. All authors made an important intellectual contribution to the manuscript and have read and approved the final manuscript.
The authors declare no conflicts of interest.
Supplementary material
For supplementary material referred to in this article, please visit https://doi.org/10.1017/S0007114522002379