Nutritional quality has recently been recognised as a fundamental sustainability issue(Reference Willett, Rockström and Loken1). In life cycle assessment (LCA), this has led to the development of nutrient indices that can be used as nutritional functional units (nFU) in LCA for foods. FU are established to compare the environmental performance of products for their intended purpose(Reference Furberg, Arvidsson and Molander2); for foods, this refers to their ability to provide nutrients. Different types of nFU can be applied, for example, they can be based on a single nutrient or a nutrient profiling score(Reference McAuliffe, Takahashi and Lee3). The Nutrient Rich Food (NRF) index(Reference Fulgoni, Keast and Drewnowski4) and its variants(Reference McAuliffe, Takahashi and Lee5) are among the most used nutrient indices used as nFU in LCA of food products. In addition, several product-group-specific indices have been developed to enable more precise evaluation of environmental impacts per nutrients provided, highlighting the role of protein sources, a product group with both a specific nutrient profile and potentially high environmental footprint(Reference McAuliffe, Takahashi and Lee5–Reference Saarinen, Fogelholm and Tahvonen7). Product-group-specific methods are based on the idea that foods in a product group can be substituted with other products in the same group(Reference Scarborough, Arambepola and Kaur8). Thus, the principle is slightly different from the across-the-board indices that are aimed to provide guidance, which foods to include in or exclude from a diet based on their healthiness or unhealthiness, respectively(Reference Scarborough, Arambepola and Kaur8). Product-group-specific indices used as nFU in LCA include, for example, NRFprotein-sub(Reference Green, Nemecek and Smetana9), NQI(Reference Sonesson, Davis and Hallström6), FNIprot7(Reference Saarinen, Fogelholm and Tahvonen7) and its variations(Reference McAuliffe, Takahashi and Lee3) for protein and Nu index for carbohydrate-rich foods(Reference Xu, Xu and Peng10).
Nutrient indices indicate the amount of selected nutrients that are provided by foods in relation to the given nutrient-specific intake recommendations. A nutrient index can be built to indicate the presence of generally beneficial nutrients, to highlight adequate intake of critical nutrients for which there is a risk of deficiency, and/or to encourage adherence to specific dietary guidelines. Public health perspectives and national dietary guidelines have been considered especially important for nutrient selection strategies(Reference Bianchi, Strid and Winkvist11). In this regard, nutrients that should be limited, such as SFA, sugar and/or Na, are sometimes included in the index. However, this can lead to negative index values and thus disable their application as nFU in LCA(Reference McLaren, Berardy and Henderson12).
Establishing the evaluation of nutrient intake from a food product on product-specific portion size (instead of, e.g. 100 g or 100 kcal) has been considered understandable for consumers and able to highlight healthy food choices(Reference Bianchi, Strid and Winkvist11). As a framework for a balanced meal, the National Nutrition Council of Finland introduced the plate model, which recommends 1/4 of plate for protein, 1/4 for carbohydrate and 1/2 for vegetables, plus one slice of bread and one portion of berries for healthy adults(13). The model is familiar to consumers, and the ratios of different food groups can be fine-tuned for the needs of different population groups. Thus, the components of the plate model can be used as a starting point for the portion sizes and also for product grouping, indicating foods that have corresponding roles in a meal or a diet. Based on this principle, we have recently introduced product-group-specific nutrient indices to be used as nFU in LCA for protein sources (NR-FIprot)(Reference Kyttä, Kårlund and Pellinen14), vegetables, fruits and berries (NR-FIveg), and carbohydrate sources (NR-FIcarb)(Reference Kyttä, Kårlund and Pellinen15). The reference unit for these indices was standardised at 100 g, and the nutrients were selected for the indices based on the national intake data from Finland(Reference Valsta, Kaartinen and Tapanainen16) to represent macro- and micronutrients for which the product group in question serves as a major source.
A more harmonised way of selecting nutrients for nutrient indices is needed(Reference Bianchi, Strid and Winkvist11,Reference Hallström, Davis and Woodhouse17) , and one of the identified areas of future research on nutrient indices as FU is the validation process for the indices(Reference McLaren, Berardy and Henderson12). Nutrient profile models, that is, ranking of foods from healthy to unhealthy based on their nutrient compositions, have been validated against nutrition professional surveys(Reference Scarborough, Boxer and Rayner18) and recommendations(Reference Arambepola, Scarborough and Rayner19) and measures of diet quality, such as the Healthy Eating Index (HEI)(Reference Fulgoni, Keast and Drewnowski4). In the case of product-group-specific indices, unsupervised statistical methods could be used for validation purposes to support expert knowledge to ensure that a nutrient index truly reflects the nutritional function suggested for a product group. These methods allow analysis of patterns in nutrient composition data without assuming any grouping of variables in advance. Statistical methods have been used also before in an attempt to establish relationships between nutritional benefits and environmental sustainability(Reference Grigoriadis, Nugent and Brereton20) and to evaluate across-the-board index performance against HEI(Reference Fulgoni, Keast and Drewnowski4), for example. However, very few previous studies have utilised statistical methods to select nutrients for nutrient indices and to validate food grouping for product-group-specific indices.
In this study, our objective was to validate how well the recently developed NR-FI indices represent the nutritional functions of the different product groups and if the suggested product grouping corresponds to these functions. For this purpose, we used a statistical dimension reduction method to identify which nutrients are correlated in foods that are typically consumed in the everyday Finnish diet. The plate model is used to estimate the portion sizes of different product groups for a healthy adult. We aim to contribute to the discussion on the strengths and limitations of different strategies to select nutrients for indices that can be used as FU in LCA.
Materials and methods
Food data
Example foods and their nutrient compositions were collected from the national Food Composition Database in Finland (Fineli)(21). The foods were selected to represent typical main and side dishes, for example, in home cooking and lunch cafeterias, excluding mixed foods, such as soups and casseroles, combining both carbohydrates and proteins in a single dish. In addition, we included common fresh and cooked fruits and vegetables in the dataset. The complete list of foods is presented in online Supplementary Table S1. The original set of foods contained seventy-two food items, including twenty-four protein dishes (eight red meat (i.e. three beef, two pork, one reindeer, one frankfurter and one liver), six plant-based, five fish, four poultry and one egg), twenty-four carbohydrate foods (seven grains, five tubers, four breads, four pastas, three morning cereals and one bun), and twenty-four fruits and vegetables (seven salad mixes, six fruits, five berries, five raw/cooked vegetables and one berry kissel). The portion sizes that were used to determine the total nutrient contents were selected based on a commonly adopted plate model containing 1/4 protein, 1/4 carbohydrate and 1/2 vegetables(13,22,23) . In addition, one small to large portion of fruit/berries and a slice of bread were considered as part of the model meal (Fig. 1)(13). The final portion sizes were anchored as follows: 100 g for protein foods, 135 g for carbohydrate foods and 350 g for fruits/vegetables. On average, the model meal contained 535 kcal.
Statistical analysis
Clustering of nutrients in food items was analysed using principal component analysis (PCA; SPSS Statistics software version 26, IBP corp.) with varimax rotation. Nutrient contents (online Supplementary Table S1) were used as variables. In PCA, the principal components represent newly created, uncorrelated variables that are linear functions of the variables in the nutrient composition data, whereas rotation aims to maximise the contribution of interrelated nutrients in one component while minimising their contribution to the other components(Reference Newby and Tucker24,Reference Varraso, Garcia-Aymerich and Monier25) . Components were first extracted based on eigenvalue threshold 1(Reference Newby and Tucker24); this resulted in seven components, explaining 79·6 % of the variance. After observing the plateauing of the eigenvalue diagram (i.e. scree plot indicating the most significant components explaining the variance) and the number of variables loaded on each component (to eliminate components with only 1–2 nutrients), the number of components was reduced to five. Because nutrients that were at high levels in minced liver steak or in carbohydrate sources containing milk (bun, multigrain bread and rice porridge) were driving separation of specific components (online Supplementary Table S2), these foods were removed from the dataset as anomalies of their corresponding product groups. Berry kissel was also removed due to very low content of any nutrients. The final list of food items is presented in Table 1. After food item exclusion, another PCA with eigenvalue threshold 1 was conducted. Based on the scree plot, components were again reduced and the final rotated component matrix contained three components. SIMCA® (Multivariate Data Analysis software version 16, Umetrics) was used to build score plots and loading plots for product types and nutrients, respectively. A heat map showing which food items were highest and which were lowest in the nutrients that were clustered in the PCA was prepared using Excel (Microsoft corp.). Food items were colour-coded based on the food group (protein, carbohydrate and fruit/vegetable), and for each nutrient, food items were sorted for their content from highest to lowest and grouped based on the three PCs (PC1, PC2 and PC3). Vitamin A was excluded from the heat map for having a loading < 0·3 in the PCA.
Comparison with NR-FI indices
The nutrients associated with each component were compared with the nutrients included in the NR-FI indices(Reference Kyttä, Kårlund and Pellinen14,Reference Kyttä, Kårlund and Pellinen15) . The baseline index for protein-rich foods NR-FIprot includes protein, Ca, Fe, Se, Zn, vitamins B6 and B12, niacin, riboflavin and thiamine, the baseline index for sources of carbohydrates NR-FIcarb index includes carbohydrates, fibre, Fe, Mg, phosphorous, potassium and folate, and the baseline index for vegetables, fruits and berries NR-FIveg index includes fibre, potassium, thiamine, and vitamins C, K and A.
The NR-FI indices were calculated with formula:
where NUTRIENTi is the amount of a nutrient in 100 g of a product and DRIi is the daily recommended intake of NUTRIENTi in the nutrition and food recommendations of The National Nutrition Council of Finland(26). The nutrient index scores were calculated for men and women aged 10–13, 14–17, 18–30, 31–60, 61–74, and over 75 years, as well as for children aged 12–23 months, 2–5 years, and 6–9 years(Reference Kyttä, Kårlund and Pellinen14,Reference Kyttä, Kårlund and Pellinen15) .
Results
Nutrient clustering in the principal component analysis
The three components that were formed in the PCA explained 53·7 % of the variance among the foods (Table 2). PC1 was positively associated with protein, niacin, vitamin B12, MUFA, iodine, vitamin D, PUFA and Se. PC2 correlated positively with potassium, vitamin K, riboflavin, vitamin C, folate, fibre, Ca, vitamin E, thiamine and vitamin B6. Mg, Fe, phosphorus, Zn and available carbohydrates were positively loaded on PC3.
Loadings < 0·3 are not included in the table. Nutrients that were both loaded on the same component and corresponding with the same NR-FI index are in bold. PC1 correlated with four nutrients included in NR-FIprot index (protein, niacin, vitamin B12 and Se), PC2 correlated with five nutrients included in NR-FIveg index (K, vitamin K, vitamin C, fibre and thiamine) and PC3 correlated with four nutrients that were included in NR-FIcarb index (Mg, Fe, P and available carbohydrate). Nutrients that are not primarily loaded on the component that cluster their corresponding NR-FI index-specific nutrients but show secondary association are italicised.
* Not included in any NR-FI baseline index.
† Not validated to associate with corresponding NR-FI baseline index.
‡ Loading < 0·3.
Riboflavin and fibre, as well as phosphorous and Zn, showed weaker positive association with another PC in addition to PC2 and PC3, respectively. Riboflavin, phosphorous and Zn showed correlation with PC1 and fibre with PC3.
Correspondence to NR-FI indices
Nutrients that are associated with each other in NR-FI indices, as well as in the PCA of the current study, are indicated in Table 2. In correspondence to NR-FIprot, the PCA was able to associate protein with niacin, vitamin B12 and Se. In addition, riboflavin and Zn showed secondary association with these nutrients. As in NR-FIcarb, carbohydrate was associated with Mg, Fe and phosphorous, and at a secondary level with fibre. Several nutrients that were included in NR-FIveg were associated with PC2. These nutrients included potassium, vitamin K, vitamin C, fibre and thiamine. Furthermore, although vitamin A, which was included in NR-FIveg, did not quite exceed the threshold that was set for qualification of loadings for further examination (> 0·3) in PCA, it was loaded (0·279) on PC2.
Nutrients that were not associated with the same nutrients in PCA as in NR-FI indices included folate, Ca and vitamin B6. In the PCA, these nutrients were correlated with NR-FIveg nutrients, although folate was originally included in NR-FIcarb, and Ca and vitamin B6 in NR-FIprot. In addition, some nutrients that were included in two NR-FI indices showed correspondence only to one. For instance, potassium was not significantly linked with NR-FIcarb-associated nutrients, and thiamine and Fe were not significantly linked with NR-FIprot-associated nutrients.
Food grouping
The heat map shows that mainly nutrients that were most abundant in the protein foods were loaded on PC1, and that mostly nutrients that were at the highest level in fruits and vegetables were loaded on PC2 (Fig. 2). Foods that were categorised as carbohydrate sources were scattered among the components, but PC3 correlated positively with Mg, potassium and carbohydrate, three nutrients that were especially abundant in morning cereal products and breads.
In the score plot, foods categorised as animal protein sources (eggs, fish, meat and poultry) were clustered together, while the fruit/vegetable group formed a separate cluster (Fig. 3, online Supplementary Fig. S1). Mixed plant protein stews containing vegetables (bean stew (V2), soy Bolognese (V4), minced faba bean and tofu wok dishes (V3 and V5)) were scored close to the carbohydrate group and also approached the fruit/vegetable group. Wheat and rice-based carbohydrate sources were closely clustered.
Overall, the results were also demonstrated by the loading plot (Fig. 4), showing that PC3 could separate carbohydrate category foods together with carbohydrate, Fe and Mg, while protein category foods were strongly characterised by protein, niacin and vitamin B12 and the fruit/vegetable category by vitamin C, vitamin K and potassium.
Discussion
Validation of nutrient index scores aims to ensure consistency with a measure of diet quality and should be implemented with science-driven tools(Reference Fulgoni, Keast and Drewnowski4). Similar tools can also be applied to validate the accordance of product-group-specific nutrient indices with dietary guidelines and suggested nutritional functions. The validation process presented here produced sets of nutrients that largely correspond to the NR-FI baseline indices that were formed based on the national nutrient intake data(Reference Kyttä, Kårlund and Pellinen14,Reference Kyttä, Kårlund and Pellinen15) . This was especially true for the carbohydrate and fruit and vegetable indices, indicating carbohydrates to be major sources of Mg, Fe, phosphorous and available carbohydrate, and fruits and vegetables as major sources of potassium, vitamin K, vitamin C, fibre and thiamine. In correspondence to NR-FIveg, the major sources of vitamin A were vegetables. However, some differences were recorded between the intake-based indices and the sets of nutrients introduced in this paper. While our intake-based NR-FI indices highlight the role of carbohydrate sources, such as bread, as sources of folate, the plate model suggests higher intake of folate from fruits and vegetables. No carbohydrate foods made it to the highest quintile for folate sources. Furthermore, although potassium was included in the NR-FIcarb index, it did not show this association in the PCA and there were no carbohydrate sources in the highest quintile for potassium content. Bread and potatoes are highly consumed staple foods in the Finnish diet and therefore serve as important sources of folate and potassium, respectively. The PCA indicates, however, that based on dietary recommendations, vegetables should be favoured as a source of folate and also as a source of potassium.
For the protein index NR-FIprot, protein, niacin, vitamin B12 and Se were validated as being provided mainly by protein sources in the plate model. In contrast, Ca, thiamine and vitamin B6 were associated with the fruit/vegetable group and Fe with the carbohydrate group. As mentioned, thiamine was included also in the NR-FIveg index and Fe in the NR-FIcarb index(Reference Kyttä, Kårlund and Pellinen15). Nevertheless, based on the nutrient content ranking, one or two of the highest contents for Ca, thiamine and Fe were measured from a protein source, indicating the relevance of meat products in providing Fe and thiamine(Reference Lombardi-Boccia, Lanzi and Aguzzi27) and fish and tofu in providing Ca(Reference Schwerbel, Tüngerthal and Nagl28–Reference Zhao, Li and Qin30). Furthermore, Zn and riboflavin showed some correlation with NR-FIprot-associated nutrients in the PCA. However, plant protein sources used in this study seemed not to provide much Zn unless they contained some form of (pseudo)cereal grains, for example, quinoa or breadcrumbs.
Although it has been acknowledged that plants provide high levels of Ca and Zn, it is also well known that some antinutrients and fibre can bind mineral nutrients in plant matrices, and therefore, these nutrients are more bioavailable from animal sources(Reference Rousseau, Kyomugasho and Celus31). In addition, non-haem Fe in, for example, legumes and cereals is notably less bioavailable in comparison with haem Fe from animals(Reference Rousseau, Kyomugasho and Celus31). Previous literature suggested that bioavailability of Ca from plant sources averages 30 %(Reference Melse-Boonstra32). We tested the PCA with Ca values corrected for bioavailability (results not shown) and noted that ranking of protein dishes as sources of Ca was improved after the correction, highlighting fish and egg dishes, in particular. The top quintile for Ca content still contained fruits (blackcurrant, orange and raspberry) and root and cruciferous vegetables. Although the single portion of 350 g is very high for blackcurrants alone, for example, it is possible to eat 350 g of a combination of different fruits, root vegetables and crucifers in one meal, making the fruit/vegetable group a major source of Ca. However, it is apparent that the evaluation of bioavailability cannot be handled using statistical models without further control exercised by nutrition experts.
An intrinsic characteristic of an average Finnish diet is the high consumption of milk and dairy products, which are naturally rich in Ca and often supplemented with vitamin D(Reference Kårlund, Kolehmainen and Landberg33). High consumption makes them also an important protein source in Finland(Reference Valsta, Kaartinen and Tapanainen16). Therefore, Ca was included in the NR-FIprot index(Reference Kyttä, Kårlund and Pellinen14). However, in the main meals of the day (lunch or dinner), milk and dairy products seldomly serve as sole or major protein sources. In the current study, it was determined that including milk-containing carbohydrate foods in the dataset drove separation of Ca and vitamin D to form their own principal component. At the same time, foods in the protein food category did not extensively contain dairy products, further loosening the link among protein, Ca and vitamin D. Instead, the interphase between plant-based protein foods and carbohydrate foods became even more indistinct. In the Finnish diet, cooking fats, fat spreads and milk as a meal drink provide most (ca. 70 %) vitamin D(Reference Valsta, Kaartinen and Tapanainen16). Because plant-based protein foods are often moderate sources of Ca but do not provide vitamin D, this supports the concept applied in NR-FIprot that the protein food group as such is not pivotal for vitamin D intake. Further, a bit surprisingly, dairy products were not a major component in protein foods in the food database recipes but more so in a few products in the carbohydrate food group. It was concluded, however, that carbohydrate sources should not be categorically promoted as good sources of Ca and vitamin D. Instead, because milk is largely consumed as a drink in Finland(34), both Ca and vitamin D could be deemed as important nutrients in the product group of beverages which was not included in this study.
It was reported that Canadian adults who get their protein mainly (≥ 75 %) from animal sources tend to receive significantly higher amounts of niacin and vitamin B6 from their diet in comparison with people consuming more (≥ 25 %) plant proteins(Reference Fabek, Sanchez-Hernandez and Ahmed35). However, individuals consuming mainly animal protein also ingest significantly more SFA and cholesterol and less dietary fibre(Reference Fabek, Sanchez-Hernandez and Ahmed35). Meanwhile, as the proportion of plant protein over animal protein in the diet increases, the odds of not fulfilling the daily recommendations for riboflavin and niacin intake increases, and also the intake of vitamin B6 is reduced(Reference Fabek, Sanchez-Hernandez and Ahmed35). As riboflavin and vitamin B6 are widely available from plant-based food sources, this suggests that the consumption of fruits and vegetables is generally at a low level in Western populations, even in groups consuming mainly plant proteins. Breads and cereal grains are important protein sources for those who favour plant-based foods(Reference Fabek, Sanchez-Hernandez and Ahmed35), but according to our ranking, the riboflavin and vitamin B6 levels in cereal grain products are not very high. Therefore, it would be advisable to direct consumers to receive riboflavin and vitamin B6 from green vegetables and fruits. Niacin, on the other hand, can be obtained from mildly processed cereal grains, nuts and legumes(Reference Meyer-Ficca and Kirkland36,Reference Oghbaei and Prakash37) , the cornerstones of plant-based proteins. Taken together, inclusion of niacin in the NR-FIprot index was validated. However, due to the overall low levels of riboflavin and vitamin B6 in plant protein foods, it could be re-evaluated whether including these vitamins in a protein index creates a bias; instead, the role of fruits and vegetables as sources of riboflavin and vitamin B6 could be highlighted.
It was determined that for the plant-based protein sources some nutrients that are indicated to differentiate the protein foods from other food groups do not occur naturally in plant protein foods at high levels. These include vitamins B12 and D, but also nutrients that are added as supplements or during cooking. For example, MUFA, PUFA and iodine were the nutrients in PC1 that had plant protein sources in the highest quintile, yet cooking oil contributes to the MUFA and PUFA content and iodised salt contributes to the iodine content of plant-based patties and wok dishes. Furthermore, even a small amount of egg in an otherwise plant-based recipe contributes to the levels of vitamins B12 and D, as determined in the case of breaded tofu (V6). This makes the plant protein product group prone to variation caused by cooking habits and may artificially improve the ranking of plant protein sources within the protein index. This might also make adding salt or fat to plant protein dishes tempting. Thus, our decision to exclude MUFA, PUFA and iodine from the NR-FIprot index is justified, although they could help to bring forward animal protein options with beneficial fatty acid profiles and natural iodine, as in fish, if included in the index.
It is notable that more than one plant protein food appeared in the highest quintile for Mg, Fe and phosphorous that were clustered in the carbohydrate source-related PC3. The observation that plant protein sources grouped close to carbohydrate foods relates also to a larger question of dietary shift and future roles for food groups in the plate model in providing nutrients(Reference Willett, Rockström and Loken1,Reference Poutanen, Kårlund and Gómez-Gallego38) . Inclusion of Ca and Fe in the NR-FIprot might help to spot the best plant protein sources for these nutrients, along with some important sources in omnivorous diets, that is, fish (especially with bones) and meat, respectively. It could be possible to develop a vegan nutrient index based on a modified plate model with a higher proportion of protein food and a lower proportion of carbohydrate food. This would potentially help to identify the best plant sources of Fe and Ca, but also of Zn and riboflavin, to ensure adequate intake and absorption(Reference Craig and Mangels39). Similarly, foods can be selected, and the portion sizes can be modified to correspond to any specific diet, dietary recommendation or food cultural context to form more relevant product-group-specific sets of nutrients. It is also possible to include nutrients to limit, although they are not recommended to be included in the nutrient indices used as nFU in nutritional LCA(Reference Saarinen, Fogelholm and Tahvonen7).
In the context of NR-FIprot, it might be concluded that taking the protein food group as a whole, it has the potential to provide protein, Ca, Se, Fe, Zn, niacin and thiamine. Other sources might better secure vitamin B6 or riboflavin intake without the risk of increased intake of low-quality fats. In addition, although legumes and nuts are regarded as potentially good sources of Se, Se content in plants depends on soil Se levels, which can vary considerably(Reference Kipp, Strohm and Brigelius-Flohé40). Eggs and meat are the most important sources of thiamine in the Finnish diet(Reference Valsta, Kaartinen and Tapanainen16), but based on the current analysis, thiamine can also be obtained from various plant sources, including foods from all three categories. Therefore, it is apparent that both protein foods and the fruit/vegetable group are relevant sources of thiamine in the plate model. Albeit that vitamin B12 is almost solely provided by animal protein, it is an essential nutritional component and thus should be included in the protein index.
To summarise the results of the validation process, it is suggested that NR-FIprot should contain protein, Ca, Fe, Se, Zn, vitamin B12, niacin, and thiamine, NR-FIcarb should contain carbohydrate, fibre, Fe, Mg, and phosphorous, and NR-FIveg should contain fibre, potassium, thiamine, vitamin C, vitamin K, vitamin A, folate, vitamin B6, riboflavin, and Ca.
Conclusions
This study largely validated the choice of nutrients for the NR-FI indices. The statistical analysis indicated that different food groups provide a distinct set of nutrients in an everyday diet. However, when the plate model was used as a basis for selecting the portion sizes for different food groups, the role of fruits and vegetables as sources of group B vitamins (except B12) and Ca, in particular, was featured more in comparison with the original indices. However, this is consistent with many dietary guidelines and suggests that nutrient indices as part of LCA can be used to encourage consumption of healthy, low-energy plant-based foods with typically low environmental impact. Although PCA is a convenient way to visualise correlations between different nutrients in foods, it requires expert interpretation and evaluation regarding bioavailability issues and overlapping of food groups in terms of nutrient compositions. Thus, selecting nutrients for nutrient indices using both statistical modelling and professional knowledge provides the best result.
Acknowledgements
The authors want to thank Dr Santtu Mikkonen for providing supervision on the statistical analyses.
This research was supported by the Ministry of Agriculture and Forestry (The Development Fund for Agriculture and Forestry), and stakeholder companies Atria Oyj, Helsingin Mylly Oy, Kesko Oyj, Oy Karl Fazer Ab, Oy Soya Ab, Vaasan Oy and Valio Oy.
A. K. collected the data and performed statistical analysis. A. K. and V. K. drafted the manuscript. All authors contributed to conceptualisation of the study, interpretation of the results and manuscript revision and approved the final version of the manuscript.
There are no conflicts of interest.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114524000709