A rapid change in the nutrition transition is occurring in Iran(Reference Ghassemi, Harrison and Mohammad1) because of demographic and physical activity level change, social development and urbanisation. Gradual dietary changes have happened from choosing healthy foods to westernised foods(Reference Popkin2) which are related to an increasing rate of non-communicable disease. Although daily intake of some healthy food groups like fruit, vegetable, dairy products and legumes has increased among the population of the Middle East from 1990 to 2017, the burden of chronic diseases rose, at once(Reference Azizi, Hadaegh and Hosseinpanah3). That finding suggests the association between single nutrients or foods and chronic disease might not explain the link because consumed foods are highly interrelated. Recent nutritional studies have focused on dietary patterns rather than single foods or nutrients to find a diet–disease association(Reference Hu4).
Meal-specific dietary patterns are a key aspect of recent nutritional epidemiology, because of the role of timing of eating in health and disease(Reference Almoosawi, Vingeliene and Gachon5). There are a few studies that investigated it(Reference Murakami, Livingstone and Sasaki6,Reference de Oliveira Santos, Fisberg and Marchioni7) . Foods are consumed on eating occasions as meals or snacks. There is no standard definition of eating occasion(Reference Palmer, Capra and Baines8). However, according to previous studies, it is known as an event that provided at least 50 kilocalories(Reference Gibney9). Meals are also known as a large amount of food consumed regularly, between 05.00 and 11.00 hours as breakfast, 11.00 and 16.00 hours as lunch and 16.00 and 23.00 hours as dinner, whereas snacks are identified as smaller intake between main meals(Reference Leech, Worsley and Timperio10,Reference Kahleova, Lloren and Mashchak11) . Accumulating evidence suggests that breakfast skipping(Reference Goyal and Julka12,Reference Hopkins, Sattler and Steeves13) , energy contribution by meals(Reference Rosato, Edefonti and Parpinel14,Reference Madjd, Taylor and Delavari15) and the number of eating occasions could affect health outcomes(Reference Palmer, Capra and Baines8,Reference Yoo, Suh and Lee16) . A habitual diet is an accumulation of eating occasions (meals and snacks)(Reference Hu4,Reference Leech, Worsley and Timperio10,Reference Gorst-Rasmussen, Dahm and Dethlefsen17) . Most of the previous dietary patterns were derived from the FFQ ignoring information on eating occasions. Such information can be important because the way specific foods are consumed during different meals may affect health outcomes(Reference Woolhead, Gibney and Walsh18,Reference Jakubowicz, Barnea and Wainstein19) . Diets can vary from day to day and a single day is a poor estimator of long-term intake(Reference Beaton, Milner and Corey20), and ≥ 3 recalls (two weekdays and one weekend) could be optimal for estimating food intake(Reference Ma, Olendzki and Pagoto21). In some studies, two or three 24-h dietary recalls (24hDR), with 1 or 2 weekdays and 1 weekend day, are used to capture energy and nutrient variability of the diet(Reference Beaton, Milner and McGuire22). Although there are many studies that have derived habitual diet from 24hDR(Reference Bakhtiyari, Ehrampoush and Enayati23–Reference Aghayan, Asghari and Yuzbashian25), only a few of them have derived meal-specific dietary patterns(Reference Baltar, Cunha and Santos26).
Time of eating as well as quality and quantity of diet is an important aspect of the healthy eating pattern(Reference Manoogian, Chaix and Panda27). It has been recommended that public health strategies and dietary advice need to focus on meal-specific dietary pattern and their contribution to habitual dietary patterns(Reference Murakami, Livingstone and Sasaki6,Reference Schwedhelm, Iqbal and Knüppel24) . Meal-specific dietary patterns give us information regarding eating behaviour: timing, frequency and regularity(Reference Almoosawi, Vingeliene and Gachon5). Understanding dietary patterns at meal levels may reflect how people actually eat and meal-specific dietary advice could likely be more practical.
Although some studies have addressed the meal-specific dietary patterns(Reference Murakami, Livingstone and Sasaki6,Reference Leech, Worsley and Timperio10,Reference Woolhead, Gibney and Walsh18,Reference Schwedhelm, Iqbal and Knüppel24,Reference de Oliveira Santos, Fisberg and Marchioni28) , only two studies from Germany and Japan(Reference Murakami, Livingstone and Sasaki6,Reference Schwedhelm, Iqbal and Knüppel24) have investigated the association of meal-specific dietary patterns to the habitual diet. Therefore, this study aimed to investigate the meal-specific dietary patterns and the contribution of dietary intake at each meal to the formation of habitual dietary patterns in a sample of Iranian adults.
Methods
Study population
A cross-sectional study was conducted among apparently healthy men and women from Iran who attended healthcare centres of Tehran to assess the association between diet quality and obesity, from February 2019 to August 2019. A sample size of 546 individuals was calculated based on the following formula(Reference Payne and Payne29) n = (z2p(1-p))/d2, according to the prevalence of obesity (68·5 %) in Tehran(Reference Kiadaliri, Jafari and Mahdavi30), an error coefficient of d = 0·04 and at α level of 0·05. Considering the effect design of 1·5 and exclusion of participants with under- and over-reporting, the final sample size of 840 participants was estimated. Participants were recruited using two-stage cluster sampling from five geographic areas of Tehran within twenty-five healthcare centres. A convenient sampling method was used to select the study participants from each healthcare centre, using the proportion-to-size approach. The inclusion criteria were having 18–59 years old and BMI of 18·5–39·9 kg/m2. The exclusion criteria were pregnancy or lactation and participants who had been diagnosed with chronic disease.
Ethical approval
Sample collection was facilitated by coordinating with the healthcare centres of Tehran. The study was ethically approved by the Ethics Committee of Tehran University of Medical Sciences (Ethics Number: IR.TUMS.VCR.REC.1398.990). The purpose of the study was explained to the participants, and all participants were given written informed consent preceding to enter the study.
Data collection
Data were collected from each person by face-to-face interview. Socio-demographic characteristics were completed that included age, sex, marriage status, income, smoking status, education level, occupation status and family size by using pre-specified data extraction forms. Dietary data were obtained using non-consecutive 3-d 24hDR conducted by trained interviewers
Dietary intake assessment
The first 24hDR was recorded in the first visit to the healthcare centre. The following 24hDR were collected via telephone on a random day. The amount of foods in gram was recorded for each eating occasion. A total of 2659 recalls were recorded, and daily intakes of all food items were derived from 24hDR converted into grams by using household measures(Reference Ghaffarpour, Houshiar-Rad and Kianfar31) (online Supplementary Table S1).
Meal definitions
Meals were known as occasions where large amounts of foods were consumed or were standardised based on time of consumption(Reference Gibney9,Reference Leech, Worsley and Timperio10) to contain no more than one breakfast, lunch and dinner, but allowing for multiple snacks. Breakfast was defined as an eating occasion where a large amount of food or energy was consumed between 05.00 and 11.00 hours and lunch, if it was consumed between 11.00 and 16.00 hours, then diner was defined as the main meal when was eaten between 16.00 and 23.00 hours based on prior studies(Reference Kahleova, Lloren and Mashchak11). We chose these cut points to enable better comparison with published studies.
Physical activity
Physical activity was measured by the short form of the International Physical Activity Questionnaire(32). The reliability and validity of the questionnaire were assessed across twelve countries(Reference Craig, Marshall and Sjöström33). Participants reported within the previous 7 d a person spent walking, doing a moderate-intensity activity and/or doing vigorous-intensity activities. The overall physical activity level in metabolic equivalent minutes per week (MET-minutes/week) was measured. MET scores were categorised into three levels: point score < 600 MET-min/week as low physical activity, point score 600–3000 MET-min/week as moderate physical activity and point score > 3000 MET-min/week as high physical activity(Reference Ainsworth, Haskell and Herrmann34).
Anthropometric assessments
Body weight was measured in participants while wearing light clothes to the nearest 0·1 kg by a digital Seca scale with a measurement accuracy 100 g(Reference Lohman, Roche and Martorell35). Height was gauged in a standing situation, shoulders and barefoot touching the wall to the nearest 0·5 cm. BMI was calculated by dividing weight in kg to height in squared metres (kg/m2).
Dietary pattern analysis
We derived 420 food items from 24hDR on the level of meals (breakfast, lunch, afternoon and dinner) that were classified into thirty-six food groups including breads, rice and pasta, other cereals, fresh fruits and juices, dried fruits, green leafy vegetables, red and orange vegetables, cabbage family, other vegetables, potato, red meat, chicken, fish, processed meat, organ meat, broth, egg, legume, nut, cheese, low-fat milk and dairy products, high-fat milk and dairy products, liquid vegetable oils, solid oil, olive oil, butter, pickle, salty snacks, sugar and sweets, industrial beverages and juices, tea, coffee, herbal teas, sauces, spices and condiments. This classification was done based on the similarity of nutrient content in each food item and previous literature(Reference Schwedhelm, Iqbal and Knüppel24,Reference Aghayan, Asghari and Yuzbashian25,Reference Manoogian, Chaix and Panda27) (online Supplementary Table S1). Food groups were adjusted for energy intake by using the residual method(Reference Willett, Howe and Kushi36). Average of dietary intake on the level of meals including breakfast, lunch, afternoon snack and dinner and on the level of day (obtained from the average of three 24hDR) was used to derive dietary patterns.
Statistical methods
Demographic characteristics of participants were compared by using χ 2 for categorical variables and t test for continuous variables between men and women. Dietary patterns were identified using principal component analysis. Principal component analysis extracts common patterns according to the correlation matrix of food intake(Reference Newby and Tucker37). The Kaiser-Meyer-Olkin (> 0·50) measured sampling adequacy, and Bartlett’s test of sphericity investigated the adequacy of test items and sample size for factor analysis. The factor loading shows the correlation between food groups and food patterns and varies from –1 to +1. A positive loading indicates a positive association with the factor, whereas a negative loading shows an inverse relationship with the factor. Larger positive or negative factor loadings for foods show which food groups are important in that component (dietary pattern). The factor loadings with a magnitude of ≥ |0·3| were written in the tables(Reference Castelló, Lope and Vioque38). The number of major dietary patterns to retain was determined based on the screen plot (factors with eigenvalues > 1·5) and the interpretability of the identified patterns. We used intraclass correlation (ICC) to measure the consistency in the consumption of foods across meals levels(Reference Landis and Koch39). Spearman’s correlation analysis was performed to obtain the correlation coefficient between average intake on meals level and pattern scores on the habitual level. Then, multiple linear regression was used to investigate the contribution of different meal dietary patterns to the habitual dietary pattern with habitual intake as the dependent variable, meal-specific dietary patterns as independent variables and sex, age, physical activity, education level, smoking status, daily energy intake and BMI as covariates. A partial R 2 value was reported for the association of the meal-specified patterns across habitual dietary patterns. Participants who consumed the meal on 3-d 24hDR were known as regular eaters and skipped the meal at least on 1 d were known as irregular eater. All statistical analysis was conducted in SPSS software (SPSS Inc., version 22). A P value less than 0·05 was defined significant.
Results
This cross-sectional study was conducted on 838 adults, including 146 men and 692 women with an age range of 20–59 years old and mean age of 42·15 (sd 10·6). Mean BMI was 27·2 (sd 4·51) kg/m2. Of 838, fifty participants had 2-d 24hDR and 788 participants had 3-d 24hDR.
Characteristics of the participants are presented in Table 1. All participants consumed ≥ 1 lunch and dinner. Two participants did not intake breakfast and six participants did not consume afternoon snacks on any of 24hDR. Smoking status, educational level, energy intake in breakfast, lunch and afternoon snacks showed a significant difference among men and women. The contribution of energy intake across three main meals and afternoon snacks is indicated in Fig. 1. The highest contribution of energy intake was observed in lunch (25·7 %), followed by dinner (20·81 %), breakfast (18·25 %) and afternoon snack (12·75 %). The regularity in breakfast consumption, number of participants who consumed the meal on all days, was the highest in comparison with other meals then followed by dinner, lunch and afternoon snacks (Table 1).
* Values are mean values and standard deviations, otherwise it is indicated.
† Number of participants consumed the meals (main meals and afternoon snack) on all recalled days.
Mean (sd) intake of food groups across meals (breakfast, lunch, evening snack and dinner) and the day is shown in Table 2. Greater mean intake of bread, organ meat, egg, cheese, butter and sugar and sweet groups was seen in the breakfast meal compared with other meals. Mean intake of red and orange vegetables, cabbage family, potato, red meat, chicken, broth, legume, olive, liquid vegetable and solid oil, pickle, salty snack, industrial beverages and juices, sauces, spices and condiment groups was great in lunch meal. In contrast, people had a greater intake of fresh fruits and juices, dried fruits, nuts, herbal tea and coffee groups in the afternoon snack. In comparison with other meals, greater rice and pasta, other cereals, other vegetables, processed meat, low-fat milk and dairy products and high-fat milk and dairy products were at dinner meal.
We identified four major dietary patterns in our population using PCA that explained 20·9 % of the total variation in the sample based on scree plot (online Supplementary Table S1). The Kaiser-Meyer-Olkin index was 0·54, and Bartlett’s test was significant (P < 0·001). Factor loading of habitual intake has been indicated in Table 3. Four major dietary patterns labelled as: ‘Western’ dietary pattern, ‘Healthy’ dietary pattern, Traditional dietary pattern, and ‘Legume and broth’ dietary pattern with an explained variance of 6·98, 5·67, 4·67 and 4·6, respectively. ‘Western’ dietary pattern was characterised by higher consumption of breads, other cereals, process meat, cheese, sugar and sweets, industrial beverages and juices. ‘Healthy’ dietary pattern was characterised by higher intake of green leafy vegetables, red and orange vegetables, egg, liquid vegetable oils and tea. Traditional dietary pattern was characterised by high positive loadings of rice and pasta, tea, herbal teas, condiments, butter and negative loading of egg. ‘Legume and broth’ dietary pattern was known by high positive loading of cabbage family, broth, legume, spices and condiments and negative loading of rice and pasta and sugar and sweet groups.
* Habitual dietary pattern derived from principal component analysis.
† Factor loading is shown in bold while absolute values ≥ 0·3.
We also extracted two principal patterns with a Kaiser-Meyer-Olkin index of > 0·5 and a significant Bartlett’s test (P < 0·001) in meal levels (factor 1 and factor 2 for breakfast, lunch, afternoon snack and dinner). Factor loadings of food groups at meal levels have been shown in online Supplementary Table S3.
We analysed the consistency of intakes (gram of food groups intake across 3 d) in the level of meals and day using ICC (Table 4). Consistency of food groups was the highest across day and breakfast. The highest consistency across the day was seen in bread (ICC = 0·41), rice and pasta (ICC = 0·51), liquid vegetable oils (ICC = 0·50) and tea (ICC = 0·58). The highest consistency in breakfast was seen in bread (ICC = 0·51) and tea (ICC = 0·60). The consistency in lunch meal was highest in rice and pasta (ICC = 0·53), green leafy vegetables (ICC = 0·41), liquid vegetable oil (ICC = 0·51) and high-fat milk and dairy products (ICC = 0·43). The other meals indicated very low consistency in the food group’s consumption.
n 838 participants with at least two 24hDR.
ICC = intraclass correlation coefficient or consistency in intake of each food groups across meals and day.
ICC > 0·41 are shown in bold as more than moderate agreement.
Spearman’s correlation analysis showed that dinner had the highest correlation with the habitual dietary pattern then followed by lunch, breakfast and afternoon snack, respectively (Table 5).
* Habitual level indicated the average food groups consumed daily. Meal level indicated average food groups consumed on the level of meal, habitual dietary pattern derived from principal component analysis.
The results of multiple linear regression for the contribution of meal-specific dietary patterns to the four habitual dietary patterns are shown in Table 6. Our results showed that factor 1 dinner (37 %) and factor 1 lunch (20 %) were major and positive contributors to ‘Western’ habitual dietary patterns. For the ‘Traditional’ habitual dietary pattern, the contribution of meal-specific dietary patterns was low. The corresponding meal-specific patterns for the ‘Legume and broth’ dietary pattern were the factor 1 lunch (11 %), largely and positively.
β, Beta-coefficient.
All adjusted model for sex, age, physical activity, education level, smoking statues, energy intake and BMI
* Factor 1 = first dietary pattern.
** Factor 2 = second dietary pattern.
Discussion
In this study, we investigated the contribution of meal-specific dietary patterns to the formation of habitual dietary patterns in the Iranian adult population. We found different contributions of the meals in the formation of habitual dietary patterns for which the consistency of food groups consumption was poor to moderate (ICC < 0·6). Multiple linear regression analysis showed dinner had the highest percentage in ‘Western’ and ‘Traditional’ dietary patterns and lunch had the most percentage in the formation of ‘Healthy’ and ‘Legume and broth’ dietary pattern of habitual dietary pattern.
Our results showed that participants had very regular breakfast consumption habits (92·8 %). Previous studies showed that regular breakfast eaters have had higher dietary quality(Reference Laska, Hearst and Lust40,Reference Azadbakht, Haghighatdoost and Feizi41) . Breakfast eating could increase the thermic effect of food and fat oxidation compared with breakfast skipping(Reference Neumann, Dunn and Johnson42). In line with our findings, a study by Vainik et al. showed that consistent eating patterns are associated with better health status. They also showed that inconsistency in meals is related to eating in the evening, eating with others, eating away from home and having consumed alcohol(Reference Vainik, Dubé and Lu43). Another study indicated eating breakfast was associated with healthier cardiometabolic profiles, and consistency in breakfast consumption could potentially reinforce this effect(Reference Iqbal, Schwingshackl and Gottschald44).
We identified four habitual dietary patterns in our study that Western and Healthy dietary patterns are similar to those found in other studies in Iran and Western countries. The characteristics of the ‘Western’ dietary pattern in our study were similar to those identified in studies from Iran(Reference Aghayan, Asghari and Yuzbashian25,Reference Esmaillzadeh and Azadbakht45–Reference Khani, Ye and Terry47) . It was also comparable with the ‘Western’ pattern in the German population identified by Schwedhelm et al.(Reference Schwedhelm, Iqbal and Knüppel24). High factor loading of sugar in this dietary pattern was not surprising. This was reported in previous studies as a source of energy intake in the Iranian diet(Reference Esmaillzadeh and Azadbakht45,Reference Asadi, Shafiee and Sadabadi46,Reference Amini, Shafaeizadeh and Zare48) . Westernised lifestyle and socio-demographic transition in Iran have led to the widespread use of dietary patterns similar to European countries(Reference Galal49). Our results also showed that the ‘Western’ dietary pattern was explained mainly by dinner followed by lunch. Previous studies indicated higher adherence to the Western diet was associated with a higher risk of chronic disease(Reference Akbarzade, Djafarian and Clark50–Reference Shang, Li and Liu52). Late-night eating has also been associated with a risk of obesity(Reference Kutsuma, Nakajima and Suwa53). Japanese women who also consumed late dinners or bedtime snacks were more likely to skip breakfast(Reference Okada, Imano and Muraki54). In contrast, several studies have suggested that eating in the morning may be protective against the development of chronic disease(Reference Bo, Fadda and Castiglione55). Adherence to Western dietary patterns especially eating more energy at dinner may produce a disruption in the circadian system that might be detrimental to health(Reference Pickel and Sung56–Reference Blancas-Velazquez, Mendoza and Garcia58). These findings could help planning for local dietary guidelines according to the time of eating to promote public health.
‘Healthy’ dietary pattern in our study was also comparable with previously reported studies(Reference Schwedhelm, Iqbal and Knüppel24,Reference Esmaillzadeh and Azadbakht45,Reference Amini, Shafaeizadeh and Zare48) and some extent similar to the Brazilian dietary pattern which was labelled as ‘vegetables/fruits’(Reference Drehmer, Odegaard and Schmidt59) and to the Chinese dietary pattern which was labelled as ‘Balanced’(Reference Xia, Gu and Yu60). Dietary patterns could vary according to sex, socio-economic status, ethnicity(Reference Nicklas, Webber and Thompson61) and food insecurity(Reference Rezazadeh, Omidvar and Eini-Zinab62). Moreover, the difference in dietary habits may be due to food beliefs, religious beliefs, cultural, socio-economic and occupation status(Reference Gilbert and Khokhar63). Extracting habitual dietary patterns shows the interrelation of foods and nutrients and gives important information about the overall intake of people. Accumulating evidence reported an association between habitual intake and chronic disease. However, as the contribution of intake of individual nutrients/foods to the habitual dietary patterns is related to the eating occasions and meals, recent publications in nutrition research have focused on meal patterns. Then, meal-specific dietary patterns provide evidence on how people consume foods and how diet quality may be influenced by the different nutritional compositions in main meals. The number of studies focusing on meal-specific dietary patterns is growing. We previously reported the association of major dietary patterns at breakfast(Reference Akbarzade, Mohammadpour and Djafarian64) and lunch(Reference Akbarzade, Djafarian and Saeidifard65) and dinner(Reference Akbarzade, Djafarian and Clark50) levels with obesity.
In this study, we also found that habitual dietary patterns were originated from a complex combination of food groups across meals. We found that dietary pattern score on the level of meals indicated dinner food intake had a greater influence on the formation of the habitual dietary pattern followed by lunch. Multiple linear regression analysis indicated dinner contributed the most percentage of the habitual dietary pattern. There are some reports from Germany, China and Brazil in adults. Schwedhelm et al. also reported that habitual patterns to some extent originate at the meal level. They showed dinner had the greatest contribution to the habitual dietary pattern followed by lunch(Reference Schwedhelm, Iqbal and Knüppel24). Murakami et al. showed major meal-specific dietary patterns in the Japanese context differentially contributed to major habitual dietary patterns. Three dietary patterns were derived for breakfast, four in lunch and five in dinner, and the most percentage of contribution to the habitual dietary patterns were seen in breakfast, lunch and dinner, respectively(Reference Murakami, Livingstone and Sasaki6). Breakfast showed the least contribution of a habitual dietary pattern among Iranian adults. Health benefits of breakfast consumption were reported in previous studies such as better metabolic and hormonal function, controlling appetite(Reference Astbury, Taylor and Macdonald66) and enough energy intake and consequently weight management(Reference Brikou, Zannidi and Karfopoulou67). It has been reported that breakfast frequency is a major factor related to habitual diet quality(Reference Hopkins, Sattler and Steeves13). We also found that regularity in breakfast consumption was high in comparison with other meals in our study that suggests breakfast meal is the best choice to have dietary interventions in Iranian adults and may help policymakers to approach effective public health policies. Additionally, it has been shown that the feeding pattern is related to circadian alignment and metabolic health.(Reference Pickel and Sung68,Reference Longo and Panda69) . In a recent study, skipping breakfast among healthy men in comparison with three meals consumer with isoenergetic diet for 6 d resulted in phase delay the circadian rhythm of the core body temperature, even though the sleep–wake cycles did not change(Reference Ogata, Horie and Kayaba70). Skipping breakfast may increase fat oxidation, insulin concentration and metabolic inflexibility that leads to a low-grade inflammation and impaired glucose homoeostasis(Reference Nas, Mirza and Hägele71). Note that differences in meals contribution of habitual dietary patterns could be because of differences in community, culture and eating habits in different populations.
Limitations
We used 24hDR as short-term dietary assessment methods, which provide more detailed information about types and amounts of food than long-term assessment methods(Reference Tucker72), although it has been shown that 24hDR are associated with a large within-person variation of dietary estimates. Moreover, previous studies assessed the validity of three 24hDR has indicated mixed results(Reference Ma, Olendzki and Pagoto73,Reference Lins, Bueno and Clemente74) , especially among populations with heterogeneity. All self-reported dietary assessment methods have measurement errors, but 24hDR have been shown to be more accurate than FFQ; furthermore, they allow for meal-specific analysis, which FFQ does not(Reference Livingstone and Black75). Besides, this study was conducted within the limited period of around 4 months, so there was no seasonal variation, which might have introduced additional bias in the assessment of average dietary patterns. Misreporting of dietary intake is a serious problem associated with self-reported dietary assessment methods, particularly among overweight and obese individuals(Reference Livingstone and Black75). Subjective selecting food groups and determining the number of factors or principal components are the limitation of using principal component analysis(Reference Zhao, Li and Gao76), the definition of eating occasions, the number of factors extracted might have a kind of inconsistency in results, so it needs careful interpretation.
Conclusion
Our results suggest that habitual dietary patterns to several extents are formed at meals level, and dinner than lunch, respectively, has a greater contribution to habitual dietary patterns. Due to circadian rhythm association of time of food intake and disease, this knowledge may help design intervention studies aimed at uncovering diet–disease associations sensitive to food combinations and meal timing.
Acknowledgements
We thank all those who participated in this study.
This manuscript has been granted by Tehran University of Medical Sciences (grant no: 45 553). The funder had not any roles in the study design, data collection and data analysis, decision to publishing or preparation of the manuscript.
A. L. and S. S.-B. contributed to conception/design of the research; Z. A., N. J. and A. Z. contributed to acquisition, analysis or interpretation of the data; A. L. drafted the manuscript; S. S.-B. and K. D. critically revised the manuscript; and S. S. B. agrees to be fully accountable for ensuring the integrity and accuracy of the work. All authors read and approved the final manuscript.
The authors report no conflict of interest.
Supplementary material
For supplementary material referred to in this article, please visit https://doi.org/10.1017/S0007114521005067