<figure><img src="image_1.png" /></figure>

HAL
open science

# Individual diet optimization in French adults shows that plant-based "dairy-like" products may complement dairy in sustainable diets

Rozenn Gazan, Florent Vieux, Anne Lluch, Stephanie de Vriese, Beatrice Trotin, Nicole Darmon

► To cite this version:

Rozenn Gazan, Florent Vieux, Anne Lluch, Stephanie de Vriese, Beatrice Trotin, et al.. Individual diet optimization in French adults shows that plant-based "dairy-like" products may complement dairy in sustainable diets. Sustainability, 2022, 14 (5), 18 p. 10.3390/su14052817. [https://doi.org/10.3390/su14052817](https://doi.org/10.3390/su14052817)

HAL Id: hal-03592347

https://hal.inrae.fr/hal-03592347

Submitted on 1 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Individual Diet Optimization in French Adults Shows That Plant-Based “Dairy-Like” Products May Complement Dairy in Sustainable Diets

Rozenn Gazan <fnref n="1" />*, Florent Vieux <fnref n="1" />, Anne Lluch <fnref n="2" />, Stephanie de Vriese <fnref n="3" />, Beatrice Trotin <fnref n="2" /> and Nicole Darmon <fnref n="4" />

**Abstract:** This study aimed to explore the potential role of plant-based "dairy-like" products (PBDL) in sustainable diets. For each individual from a representative sample of French adults (INCA2 survey 2006–2007; n = 1816), a diet optimized to be more sustainable (nutritionally adequate while having a 30% reduced carbon impact with minimal change from the actual diet and isocaloric content) was modelled. The food content of the optimized diets was compared to actual diets, with a focus on PBDL and dairy products. The presented quantitative results focused on women. Optimized diets contained more plant-based products and less meats than actual diets. PBDL products were present in 7.3% and 55.7% of the subjects' actual and optimized diets, respectively, increasing significantly from 7 to 48 g/day. Regarding dairy products, cheese decreased (–14 g/day), milk increased (+14 g/day), and yogurt remained constant (87 g/day) between the actual and optimized diets, such that the intake of this food group remained constant (214 g/day). Women for whom PBDL products were introduced as new foods in their optimized diets were found to be those with actual low energy intake (1755 kcal/d on average). As a complement to dairy products, fortified PBDL products may help to achieve more sustainable diets, especially for individuals with low energy intakes.

**Keywords:** diet sustainability; plant-based dairy products; linear programming; individual diet modelling; France

**Citation:** Gazan, R.; Vieux, F.; Lluch, A.; de Vriese, S.; Trotin, B.; Darmon, N. Individual Diet Optimization in French Adults Shows that Plant-Based "Dairy-Like" Products May Complement Dairy in Sustainable Diets. *Sustainability* **2022**, *14*, 2817. https://doi.org/10.3390/su14052817

Academic Editor: Dario Donno

Received: 2 February 2022

Accepted: 25 February 2022

Published: 28 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

The food system is one of the key contributors to climate change [1]. To reduce environmental pressure of food system on earth resources, a shift toward more sustainable dietary patterns, along with improvements in food production and large reductions in food loss and waste, is needed [2,3]. Sustainable diets are defined by The Food and Agriculture Organization of the United Nations (FAO) being as environmentally respectful, nutritionally adequate and healthy, economically fair and affordable, and culturally acceptable [4]. Plant-based diets are defined as diets that "include wholegrains, legumes, nuts and an abundance and variety of fruits and vegetables" and "can include moderate amounts of eggs, dairy, poultry and fish; and small amounts of red meat", and are recognized as more sustainable diets than current Western ones [5]. Among plant-based products, the plant-based "dairy-like" (PBDL) category includes products that are designed to

<footnote n="1">¹ MS-Nutrition, 13385 Marseille, France; rozenn.gazan@ms-nutrition.com (R.G.); florent.vieux@ms-nutrition.com (F.V.)</footnote>

<footnote n="2">² Danone Research Paris Saclay, 91120 Palaiseau, France; anne.lluch@danone.com (A.L.); beatrice.trotin@danone.com (B.T.)</footnote>

<footnote n="3">³ Alpro Comm VA—Health Affairs-Kortrijksesteenweg 1093C, 9051 Sint-Denijs-Westrem, Belgium; stephanie.devriese@alpro.com</footnote>

<footnote n="4">⁴ MOISA, Univ Montpellier, CIRAD, CIHEAM-IAMM, INRAE, Institut Agro, 34060 Montpellier, France; nicole.darmon@inrae.fr</footnote>

<footnote>* Correspondence: rozenn.gazan@ms-nutrition.com</footnote>
be similar to dairy products in term of organoleptic characteristics (color, texture, mouthfeel, etc.) and usage (e.g., with breakfast cereals, in coffee, for cooking, eaten as a dessert). PBDL products include plant-based (PB) drinks (almond, rice, and oat drinks), PB desserts (resembling dairy dessert creams), and plant-based dairy alternatives (PBDAs). PBDAs are designed to have a nutritional composition similar to dairy (protein quantity and quality; calcium, vitamin B2, vitamin B12, and vitamin D fortification) in addition to similar organoleptic and usage characteristics. PBDL products are experiencing fast growth worldwide, e.g., in France where purchased quantities at least doubled between 2015 and 2018 [6]; however, when compared to dairy products, PBDL products remain rarely consumed [7].

PBDL consumers declare that one of their key reasons for choosing plant-based products is perceived healthier status and positive impact on the environment [8–10]. Indeed, PBDL products have been reported to have a lower impact on several environmental indicators (carbon impact, total land use (not considering the type of land), energy use, nitrogen footprint, eutrophication, and acidification) on a weight basis than dairy products [11,12]. Several PBDL products provide nutritional characteristics similar to cow’s milk [11–13]. However, evaluating the sustainability of an individual food on its own is not relevant, as no single food can provide all essential nutrients and instead needs to be combined with others to constitute a healthy diet. Evaluating the contribution of individual foods to sustainable diets is more in line with the sustainable diet concept.

Diet optimization is an appropriate tool to explore solutions that seek to improve the sustainability of diets, as it allows multi-dimensional issues to be addressed simultaneously [14]. The aim is to find the optimal combination of foods as defined by an objective function for a population, a subpopulation, or an individual in order to fulfil a set of constraints related to nutrition, environment, cultural acceptability, and/or cost. Optimization studies conducted on average food consumption have indicated that reducing the carbon impact of diets by more than 30% while ensuring nutritional adequacy may impair the acceptability of dietary changes [14]. Optimizing diet individually (i.e., building one optimization model per person) rather than for a population or subpopulation enables individual food habits to be taken into account. In addition, an entire sample of individual optimized diets is generated, thereby enabling statistical tests to be run on variations between observed and optimized food quantities.

The objective of this study was to explore what the role of PBDL products could be defined within a mathematically-derived more sustainable diet for a given population. We used multi-criteria mathematical optimization of individual diets [15] applied to a representative sample of French adults to assess whether PBDL products would be present in optimized diets. More sustainable optimized diets were developed for each individual, and were defined as theoretical diets that simultaneously addressed the dimensions of cultural acceptability (i.e., the diet was designed to be as close as possible to the existing individual’s diet), health (i.e., nutritionally adequate), and the environment (i.e., the diet was designed to have a 30% reduced carbon impact).

## 2. Materials and Methods

### 2.1. *Dietary Data*

The actual dietary data used for this study (referred to as the observed diets) were obtained from the French seven-day dietary record in the Individual and National Study on Food Consumption (INCA2), conducted by the French Agency for Food, Environmental, and Occupational Health and Safety (ANSES) between December 2005 and May 2007 in samples of French adults and children [16]. The samples were obtained using a multi-stage cluster sampling technique established by the French National Institute for Statistics and Economic Studies (INSEE), and were made representative of the French population through stratification (region of residence and size of urban area) and use of the quota method (age, sex, household size, and socio-professional status of head of household). The
sampling frame was approved by the French National Commission on Informatics and
Liberty (Commission Nationale de l'Informatique et des Libertés). These analyses are
based on the adult sample (n = 2624). Under-reporters identified using the Goldberg
method [17] were removed from the analyses (n = 706), as were 55 individuals who re-
ported dietary intake for less than seven days. The final sample was composed of 1863
individuals, including 1111 women and 752 men, representative of the French population.
Because recommended nutrient levels are different for men and women, analyses were
stratified by gender.

## 2.2. *Nutritional Composition and Carbon Impact Data*

Individual dietary nutrient intake and carbon impact were estimated by matching food consumption data from the INCA2 study with the previously described SUStable database, which contains the nutritional composition and carbon impact of 212 frequently consumed French foods (covering 84% of the total adult energy intake) [18]. The nutritional composition of SUStable foods was obtained from the French nutritional composition database (French Information Center on Food Quality, CIQUAL; n = 1300 foods) and their carbon impact was obtained from a database of 402 foods built previously by the environmental consultancy firm BIO by Deloitte [19]. The SUStable database already contained four PBDL products; two others were added for this study in order to illustrate the recent diversification of the product offerings in this market compared with the survey period. The nutritional composition and carbon impact of these six PBDL products can be found in Supplemental Table S1, and were consolidated using Alpro® (European market leader in PBDL products) internal data as well as the existing literature [12]. Matching INCA2 food consumption data with the updated SUStable database enabled estimation of the nutrient intake and carbon impact of each individual diet.

## 2.3. *Food Categorization*

The categorization of 212 foods into seven groups and 33 subgroups was adapted from Gazan et al. [18] to reflect the focus of this study on PBDL products (Supplemental Table S2). PBDL products were split into three different categories: PBDAs, PB desserts, and PB drinks. PBDAs included three food types: soy-based drinks, soy-based alternatives to plain yogurt, and soy-based alternatives to fruit yogurt. PB desserts included only one food item, soy-based desserts. PB drinks included two foods, almond-based drinks and oat-based drinks.

Each PBDL product category was further divided into different food groups and subgroups. As a result of their nutritional similarity to dairy products, the PBDA category was assigned to the “dairy products and alternatives” group, which was thus composed of the Milk, Yogurt, Cheese, and PBDA subgroups. Due to the intended moment of consumption and presence of added sugars, the PB desserts category was included in the “dairy desserts and alternatives” subgroup of the “sweetened and salty products” group. Due to low protein content, the PB drinks category was included in the “sweetened drinks” subgroup within the “water and drinks” group.

## 2.4. *Nutritional Quality Indicators*

The nutritional quality of the observed diets was assessed using three previously described nutritional quality indicators: mean adequacy ratio (MAR), mean excess ratio (MER), and solid energy density (SED) [20]. In short, the MAR was estimated for each individual diet as the mean percentage of daily recommended intakes for 23 key nutrients to favor the MER was calculated for each diet as the mean daily percentage of maximum recommended values for three nutrients to limit, and the SED was calculated for each diet by dividing energy intake by ingested weight for solid food only.
## 2.5. Diet Optimization

For each individual, a multi-criteria linear programming model was developed in order to design each individual's optimized diet (Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Overview of the methodology applied to derive optimized diets.</figcaption></figure>

The model aims to achieve nutritional adequacy (i.e., achievement of recommended intakes, corresponding to intakes at the level of recommendations for each of the 34 nutrients and for each individual) with a 30% reduction in diet carbon impact (constraints of the model) while staying as close as possible to the corresponding observed diet (objective function of the model) and giving access to a broader list of foods (variables of the model) than those already consumed by each individual. Each parameter (variables, objective function, and constraints) is described more precisely below.

## 2.6. Variables and Objective Function of the Model

In order to address cultural acceptability, the optimal diet for an individual was designed to be as close as possible to the observed one in terms of food quantities. Therefore, mathematically, the objective function that was minimized was as follows:

$$ \text{Min } F_i = \sum_{j=1}^{N=214} |Q_{ij}^{\text{opt}} - Q_{ij}^{\text{obs}}| \quad (1) $$
where, for an individual *i*, $Q_{i\text{opt}}$ is the optimized quantity of food *j* and $Q_{i\text{obs}}$ the observed one. $Q_{i\text{obs}}$ can be either null (if not present in the individual's observed diet) or not null (if already consumed). *N* is the number of food variables available for the optimization, corresponding to the total number of foods consumed in the entire population sample. If not already consumed in the observed diet, *liver*, "foie gras", "oysters", "bottled water", and fortified sweetened breakfast cereals were excluded due to high frequency of avoidance, seasonal availability, high price or high and multiple fortification level, as was the case in previous studies [15,21,22].

## 2.7. *Model Constraints*

In order to address the health dimension of the sustainable diet concept, nutritional adequacy was ensured by applying constraints on dietary content in terms of macronutrients, water, seven fatty acids, ten minerals, ten vitamins, free sugars, cholesterol, and total dietary fibers (Supplemental Table S3). These constraints were based on the French RDAs (Recommended Dietary Allowances) for fibers, vitamins, and minerals (except sodium) [23], and were supplemented by the EFSA's scientific opinion on the adequate intake of water as a nutrient [24] and the Nordic Nutrient Recommendations for sodium [25], because of the absence of equivalent recommendations in France. Moreover, WHO guidelines on proteins [26], fats (including total fats, alpha-linoleic acid, linoleic acid, DHA plus EPA, omega 3, PUFA and SFA), carbohydrates, cholesterol, and free sugars [27] were used, as was the case in a previous optimization-based study showing that fulfilling French values for those nutrients was difficult even at the population level [28]. Vitamin D adequacy is not easily achieved with common foods. Putting a constraint on vitamin D intake might lead to unrealistic optimized diets; therefore, as was the case in previous studies, we did not constrain it in the present study [29]. Energy content in each optimized diet was set as being ±1% of the energy content in each individual observed diet.

In order to address the environmental dimension of the sustainable diet concept, a 30% carbon impact reduction constraint was applied for each individual diet. Thirty percent (30%) was selected as a reduction target because a previous study indicated that a more drastic carbon impact reduction would imply non-acceptable dietary changes [30].

In addition to nutritional and environmental constraints, acceptability constraints on foods, food groups, and subgroup quantities were used in the model to ensure that optimized diets stayed within the range of diets actually consumed by the general French adult population, thereby addressing the cultural dimension of the sustainable diet concept. Even if the objective function minimizes the distance to amounts of foods in observed diet, optimization might lead to a quantity for a given food that can be considered as culturally unacceptable. In order to limit this risk a maximum amount was set for each food, defined as the 95th percentile of the intake distribution among consumers of this specific food and for each gender separately. The 95th percentile of soy-based drinks was applied as a maximal amount for almond-based drinks and oat-based drinks, as these two foods were not present in the observed diets. Maximum amounts were set for food groups and subgroups as well, and were defined as the 95th percentile of the intake distribution calculated for the whole population, including both consumers and non-consumers of those specific groups and subgroups and for each gender separately. Observed 95th percentile was 0 for the PBDA subgroup because it was consumed by less than 5% of the INCA2 population in 2006–2007. As the objective of the study was focused on PBDL, the maximal amount constraint for PBDA subgroup was removed; concerning PB drinks and PB desserts products, as they are included in "sweetened drinks" and "dairy desserts and alternatives" subgroups, respectively, they were constrained by the 95th percentile of their respective subgroups. As recommended by the ANSES [28], an additional constraint was set to limit the maximal quantity of fish and seafood (including marine and freshwater fish and shellfish) to two portions per week (200 g) in order to take toxicological risk into account.
## 2.8. Model (in)feasibility

Variables and constraints define an area of possible solutions in which the objective function will select the optimized one. However, it is possible that given the variables and constraints in question, there may be no possible solution for a specific individual. This situation is called “infeasibility”, and is an interesting parameter to follow; for example, it has been previously demonstrated using the optimization approach that US 2010 Dietary Guidelines for sodium were incompatible with potassium guidelines and with nutritionally adequate diets [31]. In this study, we will assess the feasibility rate, which means the percentage of subjects whose individual diet were able to be optimized.

## 2.9. PBDL Product Tracking sub-Populations

In order to better understand the changes caused by diet modelling, sub-populations of individuals were identified based on the presence or absence of at least one PBDL product in their optimized diets, and based on the novelty (or not) of said presence vs. their observed diets (Table 1).

<table><caption>**Table 1.** Definitions of PBDL product tracking sub-populations.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">Presence of PBDL Product in Observed Diet?</th></tr><tr><th>Yes</th><th>No</th></tr></thead><tbody><tr><td rowspan="2">Presence of PBDL product in optimized diet?</td><td>Yes</td><td>OptStillPB</td><td>OptAddedPB</td></tr><tr><td>No</td><td>OptNoMorePB</td><td>OptStillNoPB</td></tr></tbody></table>

OptAddedPB: sub-population of individuals with no PBDL product in their observed diet and with a PBDL product addition in their optimized diet. OptNoMorePB: sub-population of individuals with PBDL products in their observed and not in their optimized diets. OptStillNoPB: sub-population of individuals with no PBDL product in their observed or optimized diets. OptStillPB: sub-population of individuals with PBDL products in both their observed and optimized diets. PBDL: plant-based “dairy-like”. Bold highlights the sub-population’s characteristic (addition, presence or absence of plant-based dairy products)

Four sub-populations, called “PBDL product tracking sub-populations”, were defined:

1. OptStillNoPB: sub-population of individuals with no PBDL product in their observed and optimized diets

2. OptAddedPB: sub-population of individuals with no PBDL product in their observed diet and with at least one PBDL product addition in their optimized diet

3. OptStillPB: sub-population of individuals with at least one PBDL product in both their observed and optimized diets

4. OptNoMorePB: sub-population of individuals with at least one PBDL product in their observed and not in their optimized diet.

## 2.10. Statistical Analysis

First, the nutritional quality, carbon impact, and energy intake of the observed diets were compared between individuals consuming at least one PBDL product (i.e., PBDL consumers) and individuals consuming none (i.e., non-consumers) using general linear models (GLMs) adjusted for age, family status, educational level, physical activity level, and total energy intake.

Second, the percentage of subjects whose individual diet could be optimized (feasibility rate) was estimated. For those individuals where optimization was possible, dietary changes due to optimization were analyzed by comparing the observed food quantities with the optimized ones. As compared samples (i.e., observed individual food quantities vs. optimized individual food quantities) were dependent upon each other (same individ-
uals in both samples), quantities were compared using a mixed model to account for re-
peated measurements. Composite mixed dishes were decomposed by ingredient for this
analysis.

Third, statistical analyses were performed between PBDL product tracking sub-pop-
ulations (OptStillNoPB, OptAddedPB, OptStillPB and OptNoMorePB). Energy intake
and observed and optimized food group quantities were compared between sub-popula-
tions using GLMs adjusted for age, family status, educational level, and physical activity
level. Observed and optimized quantities were compared using a mixed model to account
for repeated measurements. Comparisons were performed after summing dairy products
(Milk, Yogurt, and Cheese) and all PBDL products (PBDAs, PB desserts, and PB drinks).
Diet optimization and statistical analyses were performed using SAS version 9.4 (SAS In-
stitute, Cary, NC, USA), and an α of 5% was used for all statistical tests. Diet optimization
was performed using proc OPTMODEL.

## 2.11. *Sensitivity Analysis*

Two independent sensitivity analyses were performed. The first aimed to assess the impact of fortification on the selection of PBDL products by the model; the second aimed to assess to what extent adding (as decision variables) only PBDL products to the list of foods already consumed by the individual could help achieve nutritional adequacy with a 30% carbon impact reduction.

In the first sensitivity analysis, the model was run again by modifying the nutritional
composition of PBDL products to remove fortification. The vitamin D content and vitamin
B12 content in PBDL products were set to 0 µg/100 g, while calcium and riboflavin were
set to 12 mg/100 g and 0.01 mg/100 g, respectively, as reported for the unfortified soy-
based drink in the French food composition table [32].

In the second sensitivity analysis, two additional diet optimization models were run: one allowing only access to foods already consumed by the individual (mathematically, the number of variables N in equation (1) was reduced to the number of foods currently consumed by the individual), and the other allowing foods already consumed by the individual plus PBDL products. All the other parameters of the original model were kept unchanged. The feasibility rate was compared between the original model and these two additional models using McNemar's test. Equations for the nutritional, environmental, and acceptability constraints as well as the objective functions used (main analysis and sensitivity analysis) are presented in Supplemental File S1.

## 3. Results

The nutritional quality, carbon impact and energy intake of the observed diets were
compared between PBDL consumers and non-consumers. Then, among feasible diets the
dietary changes due to optimization were analyzed. The energy intake as well as the ob-
served and optimized food group quantities were compared between PBDL product
tracking sub-populations.

### 3.1. *Description of observed diets*

A total of 2.5% of men (n = 19) and 7.3% of women (n = 81) reported having consumed PBDL products in their dietary record. Because a very low percentage of men consumed PBDL products, all the results displayed in the core manuscript focus on women. However, the results for men are available in Supplemental Table S8, Table S9, Figure S1, Figure S2, Figure S3.

Table 2 shows the nutritional characteristics and carbon impact of the diets of PBDL product consumers vs. non-consumers. PBDL product consumers were older, more frequently had a partner and no children, and had a higher educational and physical activity level than non-consumers (Supplemental Table S4).
<table><caption>**Table 2.** Mean daily energy intake, total food quantity consumed, carbon impact, and indicators of nutritional quality in observed diets (n = 1111) according to PBDL product consumption.</caption><thead><tr><th rowspan="2"></th><th colspan="2">All (n = 1111)</th><th colspan="2">PBDL Product non-Consumers (n = 1030)</th><th colspan="2">PBDL Product Consumers (n = 81)</th><th rowspan="2">Pval *</th><th rowspan="2">Pval †,§</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td>Energy intake (kcal/d)</td><td>1902.4</td><td>482.5</td><td>1895.2</td><td>484.5</td><td>1993.7</td><td>448.2</td><td>0.077</td><td></td></tr><tr><td>Total quantity (g/d)</td><td>2586.0</td><td>767.2</td><td>2566.0</td><td>762.0</td><td>2840.8</td><td>791.5</td><td>0.002</td><td>0.128</td></tr><tr><td>Liquid quantity (g/d)</td><td>1489.9</td><td>646.9</td><td>1478.3</td><td>643.7</td><td>1596.2</td><td>680.7</td><td>0.114</td><td>0.530</td></tr><tr><td>Solid quantity (g/d)</td><td>1099.1</td><td>314.0</td><td>1087.7</td><td>310.6</td><td>1244.6</td><td>323.7</td><td>&lt;0.001</td><td>0.003</td></tr><tr><td>PBDL products (g/d)</td><td>6.5</td><td>33.5</td><td>0</td><td>0</td><td>89.0</td><td>90.2</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Energy density £ (kcal/g)</td><td>165.3</td><td>33.8</td><td>166.4</td><td>34.2</td><td>151.8</td><td>25.5</td><td>&lt;0.001</td><td>0.003</td></tr><tr><td>MAR (%)</td><td>81.2</td><td>11.3</td><td>80.8</td><td>11.4</td><td>86.7</td><td>8.9</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>MER (%)</td><td>122.1</td><td>22.7</td><td>122.3</td><td>23.0</td><td>119.0</td><td>19.4</td><td>0.201</td><td>0.049</td></tr><tr><td>Carbon impact (g CO<sub>2</sub> eq/d)</td><td>3684.2</td><td>954.8</td><td>3687.3</td><td>957.0</td><td>3644.4</td><td>931.5</td><td>0.697</td><td>0.680</td></tr><tr><td>Total quantity of animal products (g/d) ¥</td><td>402.0</td><td>159.2</td><td>407.3</td><td>159.8</td><td>334.8</td><td>134.6</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>No. of different consumed foods per week</td><td>55.2</td><td>12.7</td><td>54.8 ± 12.6</td><td>12.6</td><td>60.1</td><td>12.5</td><td>&lt;0.001</td><td>0.003</td></tr></tbody></table>

* Test of the differences in intake or indicators between PBDL product consumers and non-consumers using a linear generalized model; † Test with adjustment for total energy intake, age, family status, educational level, and physical activity level (energy intake adjustment was removed from the difference in energy intakes test); £ Energy density was calculated on consumed solid foods only; ¥ Total quantity of animal products consumed, including all meat, eggs, fish, and dairy products as well as animal products present as ingredients in mixed dishes. MAR: mean adequacy ratio; MER: mean excess ration; PBDL: plant-based “dairy-like”.

Mean energy intake and mean dietary carbon impact (approx. 3684 g CO₂ eq/d) did not differ significantly between PBDL product consumers and non-consumers. When compared to non-consumer women, female PBDL product consumers had diets with higher nutritional quality, as evidenced by a higher MAR, a lower MER, a lower energy density, and a higher number of different consumed foods (a proxy of diversity).

Results for men were very similar to those obtained for women (Supplemental File Table S8). One difference was the carbon impact of diets, which was lower for PBDL product consumers than for non-consumers in men. Another difference between men and women involved nutritional quality indicators, which were not significantly different between PBDL product consumers and non-consumers in men, while they were in women.

## 3.2. *Description of Optimized Diets*

The feasibility rate (i.e., the percentage of individuals for whom an optimized diet could be designed) of the model was 97% (n = 1077) in women and 98% (n = 739) in men, for a total feasible sample of 1816 individuals. All of the following results involve the 1077 women for whom it was possible to design an optimized diet with the model.

Figure 2 shows the mean quantities of food groups and subgroups in the observed and optimized diets.
<figure><img src="image_5.png" /><figcaption>(b) Grains & potatoes</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>(c) Meat/fish/eggs & legumes</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>(d) Dairy products & alternatives</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(e) Sweetened & salty products</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(f) Water & drinks</figcaption></figure>

<figure><img src="image_10.png" /></figure>

**Figure 2.** Mean quantities of the fruits and vegetables<sup>a</sup> (a), grains and potatoes (b), meat/fish/eggs/legumes (c), dairy products and alternatives (d), sweetened and salty products (e), and water and drinks (f) groups and subgroups in the observed diets of French women (n = 1077) and in their optimized diets (n = 1077). OBS: observed diet; OPT: optimized diet. * Indicates significance.
cantly (*p*-value < 0.05) different quantities; *p*-value was estimated from a mixed linear model accounting for repeated measurements in order to assess the difference in food intake in observed or optimized diets. <sup>a</sup> including oilseeds; <sup>b</sup> including soups; <sup>c</sup> including processed fruits and excluding fruit juices; <sup>d</sup> including offal; <sup>e</sup> including tofu; <sup>f</sup> including PB desserts; <sup>g</sup> including pastries; <sup>h</sup> including PB drinks; <sup>i</sup> fruit juices with no added sugar. The total quantity of “fats and seasonings” (not displayed in Figure 2) decreased from 52.9 g in observed diets to 32.4 g in optimized ones. Composite dishes were disaggregated within subgroups.

The “fruits and vegetables” and “grains and potatoes” groups increased overall. The “meat/fish/eggs/legumes” group decreased due to the decrease in “ruminant meat”, “deli meat” and “pork and poultry” (“fish and seafood” and “eggs” remained constant). The “sweetened and salty products” group decreased. The “dairy products and alternatives” group remained constant, while within this group “cheese” decreased, “milk” and “PBDA” increased, and “yogurt” remained constant.

In men (Supplemental Figure S1), dietary changes induced by the model were comparable with the results in women.

### 3.3. Contribution of Dairy and PBDL Products to Optimized Diets

The four sub-populations defined for PBDL product tracking (OptStillPB, OptAddedPB, OptNoMorePB and OptStillNoPB) contained 78, 523, one, and 476 women, respectively. As the sub-population OptNoMorePB contained only one subject, it was left out of subsequent analyses.

Among optimized diets, 55.7% contained PBDL products either because they were introduced as a new food (OptAddedPB, 48.5%) or because they were maintained by the model (OptStillPB, 7.2%). Approximately 44.3% were optimized diets with no PBDL products (OptStillNoPB).

Table 3 compares the energy and food contents of the observed and optimized diets according to PBDL product tracking sub-populations (OptStillNoPB, OptStillPB or OptAddedPB), with a focus on dairy products (“milk”, “yogurt”, and “cheese”) and on the three PBDL product categories (PBDAs, PB desserts, and PB drinks). Supplemental Table S5 displays the same information as Table 3 for all food groups and subgroups.

<table><caption>**Table 3.** Mean energy content, mean diet weight, and mean food content (g/d) of the observed and optimized diets according to PBDL product tracking sub-populations<sup>*E*</sup>.</caption><thead><tr><th rowspan="2"></th><th colspan="4">OptStillNoPB (N = 476)</th><th colspan="4">OptAddedPB (N = 523)</th><th colspan="4">OptStillPB (N = 78)</th><th rowspan="2">Observed Pval *</th><th rowspan="2">Optimized Pval *</th></tr><tr><th>Observed Mean</th><th>Observed SD</th><th>Optimized Mean</th><th>Optimized SD</th><th>Observed Mean</th><th>Observed SD</th><th>Optimized Mean</th><th>Optimized SD</th><th>Observed Mean</th><th>Observed SD</th><th>Optimized Mean</th><th>Optimized SD</th></tr></thead><tbody><tr><td>Total energy content (kcal/d)</td><td>2119.0</td><td>438.5</td><td>2104.4</td><td>427.7</td><td>1755.9</td><td>373.1</td><td>1751.3</td><td>362.6</td><td>2001.0</td><td>420.3</td><td>1986.2</td><td>410.6</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Total diet weight (g/d)</td><td>2961.4</td><td>680.1</td><td>2984.7</td><td>580.0</td><td>2249.2</td><td>639.3</td><td>2582.8</td><td>411.8</td><td>2878.7</td><td>776.3</td><td>2945.8</td><td>604.4</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Dairy products & alternatives</td><td>274.2</td><td>165.5</td><td>265.7<sup>ns</sup></td><td>151.7</td><td>166.3</td><td>113.8</td><td>209.7</td><td>125.5</td><td>241.5</td><td>138.6</td><td>252.4<sup>ns</sup></td><td>129.5</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Milk</td><td>132.6</td><td>152.6</td><td>136.6<sup>ns</sup></td><td>149.7</td><td>65.0</td><td>99.3</td><td>88.8</td><td>103.0</td><td>61.4</td><td>103.1</td><td>76.1<sup>ns</sup></td><td>102.5</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Yogurt</td><td>107.4</td><td>81.1</td><td>110.2<sup>ns</sup></td><td>80.0</td><td>71.3</td><td>67.3</td><td>68.9<sup>ns</sup></td><td>64.1</td><td>68.8</td><td>75.3</td><td>66.9<sup>ns</sup></td><td>74.9</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Cheese</td><td>34.2</td><td>26.1</td><td>18.9</td><td>15.1</td><td>30.0</td><td>19.7</td><td>17.1</td><td>13.9</td><td>30.9</td><td>23.4</td><td>18.3</td><td>13.5</td><td>0.047</td><td>0.354</td></tr><tr><td>PBDAs</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>34.9</td><td>69.5</td><td>80.3</td><td>93.2</td><td>91.1<sup>ns</sup></td><td>103.1</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Sweetened & salty products</td><td>131.8</td><td>78.7</td><td>118.5</td><td>67.5</td><td>118.5</td><td>69.7</td><td>91.9</td><td>60.9</td><td>126.0</td><td>61.2</td><td>113.2<sup>ns</sup></td><td>57.1</td><td>0.004</td><td>&lt;0.001</td></tr><tr><td>Including PB desserts</td><td>0.0</td><td>0.0</td><td>0.00</td><td>0.0</td><td>0.0</td><td>0.0</td><td>4.7</td><td>13.1</td><td>9.6</td><td>20.8</td><td>10.7</td><td>21.1</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Water & drinks</td><td>1610.6</td><td>630</td><td>1606.5<sup>ns</sup></td><td>589.8</td><td>1221.2</td><td>599.2</td><td>1372.1</td><td>479.4</td><td>1535.9</td><td>681.8</td><td>1570.0<sup>ns</sup></td><td>616.8</td><td>&lt;0.001</td><td>&lt;0.001</td></tr><tr><td>Including PB drinks</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>42.5</td><td>51.0</td><td>0.0</td><td>0.0</td><td>15.1</td><td>31.7</td><td>&lt;0.001</td><td>&lt;0.001</td></tr></tbody></table>

<footnote>* *p*-value of a general linear model test used to assess the difference in food intake in observed or optimized diets between the tracking of PBDL products, adjusted for age, family status, educational level, and physical activity level. <sup>ns</sup> indicates food groups or subgroups for which the differ-</footnote>
ences (accounting for repeat measurements) between observed and optimized diets were not sig-
nificant (p > 0.05) among OptStillNoPB, OptStillPB or OptAddedPB women. ⁴Individuals with no
PBDL product in their observed and optimized diets were categorized as OptStillNoPB. Individu-
als without PBDL products in their observed diet and with an addition of said products to their
optimized diet were categorized as OptAddedPB. Individuals with PBDL products in both their
observed and optimized diets were categorized as OptStillPB. The distinction between
OptStillNoPB, OptAddedPB and OptStillPB refer to the “tracking” of PBDL products (i.e., pres-
ence of PBDL products in optimized diets and novelty or not vs. observed diet). PB: plant-based;
PBDAs: plant-based dairy alternatives; PBDL: plant-based “dairy-like”.

Women for whom PBDL products were introduced as a new food in their optimized diet (OptAddedPB sub-population) were characterized by lower observed caloric intakes than the two other sub-populations (1756 kcal/d in OptAddedPB vs. 2001 kcal/d in OptStillPB and 2119 kcal/d in OptStillNoPB). In the “dairy products and alternatives” group, no significant subgroup variation was found for OptStillNoPB and OptStillPB subpopulations except for a decrease in cheese. However, for women for whom PBDL products were introduced as a new food (OptAddedPB), most “dairy products and alternatives” subgroups changed significantly; “milk” increased from 65.0 to 88.8 g/day, “cheese” decreased from 30.0 to 17.1 g/day, and “PBDA” increased from 0 to 34.9 g/day.

Figure 3 shows the total amount of dairy products (summing “milk”, “yogurt”, and “cheese”) and PBDL products in the observed and optimized diets in both the overall sample and according to PBDL product tracking sub-populations.

<figure><img src="image_11.png" /></figure>

**Figure 3.** Dairy products and PBDL products (g/d) in observed and optimized diets according to PBDL product tracking sub-populations. PBDL products include PBDAs (soy-based drinks, soy-based alternatives to plain yogurt, and soy-based alternatives to fruit yogurt), PB desserts (soy-based desserts), and PB drinks (almond-based drinks and oat-based drinks). \* *p*-value < 0.05; *p*-value was estimated from a mixed linear model, accounting for repeat measurements to assess the difference in food intake in observed or optimized diets. NS shows no significant differences in food intake between observed and optimized diets.

The total amount of dairy products plus PBDL products (i.e., summing milk, yogurt, cheese, PBDAs, PB desserts, and PB drinks) significantly increased, from 220 g/d in observed diets to 262 g/d in optimized ones. The increase was mainly due to the introduction of or the increase in PBDL products, as the amount of dairy products (i.e., sum of milk, yogurt, and cheese) did not differ between observed and optimized diets (214 g/d, *p* = 0.951). Women with no PBDL products in their observed and optimized diets (OptStillNoPB) already had a high quantity of dairy products in their observed diet (above 270 g/d). Women for whom PBDL products were introduced as a new food in their optimized diet (OptAddedPB) initially had the lowest consumption of dairy products in
their observed diet (166 g/d). Moreover, in women with PBDL products in both their observed and optimized diets (OptStillPB), the total quantity of dairy products and PBDL products was maintained overall. In optimized diets, the total amount of dairy products plus PBDL products was similar in all three sub-populations of women (average amounts of 266 g/d in OptStillNoPB, 278 g/d in OptStillPB and 257 g/d in OptAddedPB, p-value = 0.205).

### 3.4. Sensitivity Analyses

The first sensitivity analysis examined the impact of fortification on the introduction of PBDL products in optimized diets.

When fortification was removed from PBDL products (Figure 4b), the percentage of optimized diets containing PBDL products was only 13% (vs. 56% when fortified; see Figure 4a). Removing fortification of PBDL products had a limited impact on the other food group and subgroup changes (Supplemental Table S6). The main change was that the increase of milk in optimized diets was higher when PBDL products were not fortified than when they were.

<figure><img src="image_12.png" /><figcaption>**Figure 4.** Percentage of women having PBDL products in observed and optimized diets, using fortified (a) or unfortified (b) PBDL products. Because three diets could not be optimized when PBDL products were not fortified, the number of diets is *n* = 1077 in panel (a) and *n* = 1074 in panel (b). PBDAs (soy-based drinks, soy-based alternatives to plain yogurt, and soy-based alternatives to fruit yogurt); PB desserts (soy-based desserts); PB drinks (almond-based drinks and oat-based drinks).</figcaption></figure>

In men (Supplemental Figure S3), the impact of fortification on the presence of PBDL products in optimized diets was comparable to that in women.

The second sensitivity analysis examined the impact on the feasibility rate of adding PBDL products only to the list of foods already consumed by an individual. The feasibility rate was 25.9% when only currently-consumed foods were introduced in optimized diets. When PBDL products were introduced in all optimized diets, the feasibility rate increased to 35.0%.

In men (Supplemental Figure S3), feasibility rates resulting from the original and sensitivity analysis models were comparable with the results obtained for women.
4. Discussion

Fully in line with the guiding principles provided by the FAO [5] and with other
studies on sustainable diets [33–35], more sustainable optimized diets in this study were
composed of a higher amount of plant-based foods (fruits, vegetables, whole grains, nuts,
and seeds) and a reduced amount of meats (especially red meat) than current observed
diets. This study can thus provide an additional argument supporting the shift towards
more plant-based diets in Western countries. The specificity of this study was its focus on
PBDL products as plant-based foods with recent tremendous growth in the Western food
market. The ability of these products to be part of a sustainable diet is poorly documented.
In our study, observed diets containing PBDL products showed better nutritional quality,
although their carbon impact was not lower than that of observed diets not containing
PBDL products. Diet optimization results showed that in more sustainable diets total
dairy products are maintained with a shift amongst the dairy products, and PBDL prod-
ucts (when fortified) are added for a large proportion of women. This suggests that both
dairy and fortified PBDL products may facilitate the achievement of more sustainable di-
ets. Indeed, the stable contribution of dairy in more sustainable diets is consistent with
previous studies conducted with different European populations and/or methodologies
[30,36,37]. In addition, the quantity of dairy in optimized diets (214 g/d) is in line with the
recommended intake of dairy foods (two servings/day, with recommended serving of 150
mL for milk, 125 g for yogurt, and 30 g for cheese) in the opinion of the French High
Council for Public Health (Haut Conseil de la Santé Publique-HCSP) revision of the 2017–
2021 French Nutrition and Health Programme dietary guidelines [38].

Regarding PBDL products, we identified two studies based on food consumption
surveys reporting their contribution to diets with good sustainability characteristics
[39,40]. First, using five national dietary surveys (France, the United Kingdom, Sweden,
Italy and Finland), Vieux et al. [39] used multiple factor analysis and agglomerative hier-
archical clustering that identified six clusters of individuals based on their nutrient intakes
and the carbon impact of their diets. One of the clusters was considered to be more sus-
tainable because it was composed of individuals with good nutritional quality (based on
MAR, MER and energy density) and decreased carbon impact (minus 15%) when com-
pared with the whole sample. This cluster contained diets with the highest amount of
PBDL products (9.35 g/day). In the same study [39], the authors directly selected individ-
uals with high MAR, low MER, and low energy density compared with their respective
median value [20]; the results showed that PBDL products were higher in the diets of
those individuals than in the rest of the sample.

Secondly, Ridoutt et al. [40] selected individuals from the Australian Health Survey with diets having both higher nutritional quality and lower carbon impact (–43%) than others. In adult women, those diets contained 1.47 servings/day of dairy and alternatives. In the present study, when expressing the quantity of dairy products and alternatives in the more sustainable optimized diets in a unit comparable to Ridoutt et al. we obtain a total of 1.52 servings (0.44 servings of milk, 0.45 of cheese, 0.43 of yogurt, 0.12 of PBDAs, and 0.08 of plant-based drinks), which is very close to the 1.47 servings in the “high quality, low carbon impact” group of individuals from the Ridoutt study. Nevertheless, neither the study conducted by Vieux et al. nor the one conducted by Ridoutt et al. is fully comparable to this one. Those studies identified self-selected diets based on their existing nutritional and environmental characteristics, while in this study the diets studied were theoretical nutritionally adequate diets with a 30% carbon impact reduction designed with multicriteria optimization. On the one hand, optimized diets are theoretical, as their cultural acceptability cannot be guaranteed. On the other hand, full nutritional adequacy cannot be ensured when identifying diets that are more sustainable than others among existing ones. Nevertheless, reaching similar conclusions (in terms of dairy and PBDL product quantities) with different data and different methods is reassuring concerning the relevance of our findings. It has been recently shown in the framework of a simulation study that evaluated the impact of like-for-like substitution of dairy with PBDAs on the
nutrient composition of UK diets that PBDAs do not provide the same nutritional benefits as dairy [41]. Using an alternative methodology and taking carbon impact into account, our results confirm the message that PBDAs cannot be considered as a nutritional replacement for dairy.

As shown by Mäkinen et al. [11], the ability of PBDL products to contribute to healthy diets partly depends on their fortification levels. Similarly, we have shown that the contribution level of PBDL products in optimized diets is dependent on their fortification. The fortification strategy thus seems to be of importance in increasing the ability of PBDL products to be part of sustainable diets.

Another result of our study is that new foods not currently consumed by an individual need to be added to the diet in order to design a diet that both complies with nutritional recommendations and has a reduced carbon impact. This result has already been documented in previous studies with and without environmental constraints [42,43]. Another result is that while only allowing PBDL products to be added increased the feasibility rate for a part of the sample, for most individuals it was not sufficient to derive a nutritionally adequate diet with reduced carbon impact.

While many studies have used linear programming to design optimized diets that comply with nutrient recommendations and/or environmental targets at the population level, only a few have applied the algorithm at the individual level. The strengths of this approach have been documented by Gazan et al. [14], and are further emphasized in this paper. Individual diet modelling provides more relevant and more robust results because it enables taking individual preferences into account, can estimate the percentage of people for whom designing an optimized diet is feasible, and allows statistical analyses to be performed in order to test the variations between individual observed and optimized diets. To the best of our knowledge, only one other study has used multi-criteria individual linear programming to design nutritionally adequate diets while reducing carbon impact [43]. That study was based on NDNS (National Diet and Nutrition Survey in the United Kingdom) data consisting of a four-day dietary diary. Nevertheless, the authors applied different optimization models from those used in this paper, thus making it difficult to compare with this study.

One limitation of this study is that certain characteristics of the sustainable diet concept were not considered. In particular, the economic dimension was not considered because of the lack of data on food prices. Including food prices in our models would have enabled controlling the cost of optimized diets, an important characteristic of sustainable diets [30]. Moreover, within the environmental dimension we used only the carbon impact metric, which is not sufficient to represent the entire dimension. Other metrics, such as water footprint and impact on biodiversity, should complement it. In September 2020, a consolidated French database named Agribalyse was released [44]; it is the first free access database of environmental impact indicators for food items produced and consumed in France. It would be valuable to repeat these analyses with the metrics provided in this new database (land use, acidification, water use, fine particles, etc.) in addition to carbon impact. Food can have varying impacts depending on the environmental sustainability metrics selected (carbon impact, water footprint, nitrogen fixation, land use, eutrophication, etc.), depending on location and agricultural practices [45]. Therefore, other studies including other environmental metrics need to be performed in order to explore how PBDL products can play a role in improving sustainability of diets, as such metrics were not available at the time of the present study. Even if the acceptability of dietary changes was considered by minimizing the distance to individuals' observed food quantities and introducing constraints to force optimized diets to remain within the ranges of food quantity already observed, the acceptability of the proposed dietary changes by the models was not ensured, as INCA2 participants were not offered their own optimized diet in an intervention study. On average, in women's optimized diets the summed weight of the food items for which consumption is increased went up by 467 g/day, and the summed weight of the food items for which consumption is decreased went down by 290 g/day
(Supplemental Table S7). In other words, an average total variation of 757 g (467 + 290) was the smallest effort needed in order to meet nutritional adequacy while reducing the diet's carbon impact by 30%. An intervention study would be needed to test how acceptable such dietary changes would be in real life. Shifting the individual optimization approach from a methodology to conduct public health research to a tool able to provide tailored dietary advice has been previously demonstrated in a PhD thesis, leading to an innovative prototype [46]. This tool must be further developed by testing it in intervention studies in order to assess the acceptability of dietary advice based on it.

It should be noted that the food consumption data used in this study are fourteen years old, creating doubt regarding potential shifts in dietary patterns over more recent years which might affect real-world interpretation of these results. However, a French report evaluating food habit changes between 2006 and 2015 [47] indicates that there were no major changes during this period. A more recent national dietary survey was conducted in 2014–2015 (INCA3) [48]; however, data were only made available in 2020, after the start of this study. INCA2 and INCA3 studies were not conducted with the same methodology, and are thus not fully comparable, although we were able to verify that the prevalence of PBDL product consumers remained low (4.2% in women and 2.3% in men) and that the amounts consumed remained similar (102 g in women and 139 g in men).

We chose to present only the results obtained for women here, because of very low consumption of PBDL products in men. Moreover, it is more difficult to optimize nutritionally adequate diets in women because their energy intake is lower than in men, while their recommended nutritional intake levels are either similar or higher (e.g., in adults over age 55, the RDA for calcium is 1200 mg/day for women and 900 mg/day for men). To design more sustainable optimized diets, the algorithm thus required selection of foods with the best compromise between energy density, nutrient density, and carbon impact. Most of the results indicate the same trends in both genders. In particular, the analysis showed that both men and women for whom PBDL products were introduced in optimized diets (OptAddedPB) had lower observed energy intakes on average, showing that PBDL products can facilitate diet sustainability in these sub-populations.

Data from Kantar Worldpanel [6] indicate that in 2018 the amount of PBDL products purchased in France was far smaller than that of dairy products (96% vs. 4%). However, penetration rates of plant-based alternatives to milk, yogurt, and desserts increased by 56%, 71%, and 76%, respectively, between 2015 and 2018, while remaining stable for dairy products (–0.5%). These signs indicate that PBDL products could become an important part of diets in the coming years.

**5. Conclusions**

In conclusion, our results for French adults suggest that consuming fortified PBDL products in addition to dairy products can help increase diet sustainability, whereas our models did not find that substituting dairy products with PBDL products entirely was an effective strategy in designing more sustainable diets. As a complement to dairy products, fortified PBDL products may help to achieve more sustainable diets, especially for individuals with low energy intake.

**Supplementary Materials:** The following supporting information can be downloaded at www.mdpi.com/article/10.3390/su14052817/s1, File S1: Mathematical description of the objective functions and constraints applied to all models (main and sensitivity analyses); File S3 Table S1: Carbon impact and nutritional composition of each PBDL product in comparison with dairy products, sweetened drinks and desserts (/100 g); Table S2: Number of foods per food subgroup; Table S3. List of nutritional constraints included in the models; Table S4. Socio-demographic characteristics of women (n = 1111) according to their consumption of PBDL (plant-based “dairy-like”) products; Table S5. Mean energy content, mean diet weight and mean food group and subgroup quantities (g/d) in observed and optimized diets according to the tracking of PBDL products in optimized diets: i.e., OptStillNoPB (n = 476), OptStillPB (n = 78) and OptAddedPB (n = 523); Table S6. Average variations of food group and subgroup quantities (in g/d) between observed and optimized diets
using fortified (1077 feasible diets) or unfortified (1074 feasible diets) PBDL products; Table S7. Distribution of the objective function (g) in women; Table S8. Mean daily energy intake, total food quantity consumed, carbon impact and indicators of the nutritional quality of men's observed diets (n=752) according to PBDL product consumption; Table S9. Mean energy content, mean diet weight and mean dairy products & alternatives quantities (g/d) in observed diets and diets optimized with the "SUS-All-Foods" set of models, according to the tracking of PBDL product in the optimized diets: i.e. OptStillNoPB (n=369), OptAddedPB (n=351) and OptStillPB (n=19); Figure S1. Mean quantities of Fruits & vegetables, Grains & potatoes, Meat/fish/eggs & legumes, Dairy products & alternatives, Sweetened & salty products, Water & drinks groups and subgroups in the observed diets of French men (n=739) and in their optimized diets (n=739); Figure S2. Dairy products and PBDL products (g/d) in observed and optimized diets, according to PBDL product tracking subpopulations; Figure S3. Percentage of men having PBDL products in observed and optimized diets, using fortified (A) or unfortified (B) PBDL products.

**Author Contributions:** Conceptualization and methodology, B.T., N.D., A.L., S.d.V. and F.V.; software and formal analysis, R.G. and F.V.; resources and data curation, S.d.V. and R.G.; writing—original draft preparation, R.G. and F.V.; supervision, project administration and funding acquisition, B.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** Analyses of publicly available INCA2 databases were sponsored by Danone Research and conducted by MS-Nutrition.

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the French National Commission for Computed Data and Individual Freedom (Commission Nationale de l'Informatique et des Libertés, CNIL).

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** The dietary data used for this study are available; see [49].

**Conflicts of Interest:** FV and RG are employees of MS-Nutrition. BT and AL are employed by Danone Research. SDV is employed by ALPRO. ND has no competing interests to report. Analyses of publicly available INCA2 databases were sponsored by Danone Research and conducted by MS-Nutrition.

## References

1. WHO Regional Office for Europe. *Healthy and Sustainable Diets—Key Workstreams in the WHO European Region*; WHO Regional Office for Europe: Copenhagen, Denmark, 2021.

2. Steffen, W.; Richardson, K.; Rockstrom, J.; Cornell, S.E.; Fetzer, I.; Bennett, E.M.; Biggs, R.; Carpenter, S.R.; de Vries, W.; de Wit, C.A.; et al. Planetary boundaries: Guiding human development on a changing planet. *Science* **2015**, *347*, 1259855. https://doi.org/10.1126/science.1259855.

3. Garnett, T. *Food Climate Research Network. Changing What We Eat: A Call for Research & Action on Widespread Adoption of Sustainable Healthy Eating*; Food Climate Research Network: Oxford, UK, 2014.

4. Burlingame, B.A.; Dernini, S.; Nutrition and Consumer Protection Division FAO Sustainable diets and biodiversity: Directions and solutions for policy, research and action. In *Proceedings of the International Scientific Symposium on Biodiversity and Sustainable Diets: United against Hunger*, Rome, Italy, 3–5 November 2010; FAO: Rome, Italy, 2012.

5. FAO; WHO. *Sustainable Healthy Diets—Guiding Principles*; FAO: Rome, Italy; WHO: Rome, Italy, 2019.

6. FranceAgriMer. *Données et Bilans de FranceAgriMer—La Consommation De produits Laitiers en 2018*; FranceAgriMer: Montreuil, France, 2019.

7. Agence Nationale de Sécurité Sanitaire de L'alimentation de L'environnement et du Travail. *Étude Individuelle Nationale des Consommations Alimentaires 3 (INCA 3)*; Agence Nationale de Sécurité Sanitaire de L'alimentation de L'environnement et du Travail: Maisons-Alfort, France, 2017.

8. Schiano, A.N.; Harwood, W.S.; Gerard, P.D.; Drake, M.A. Consumer perception of the sustainability of dairy products and plant-based dairy alternatives. *J. Dairy Sci.* **2020**, *103*, 11228–11243. https://doi.org/10.3168/jds.2020-18406.

9. Haas, R.; Schnepps, A.; Pichler, A.; Meixner, O. Cow Milk versus Plant-Based Milk Substitutes: A Comparison of Product Image and Motivational Structure of Consumption. *Sustainability* **2019**, *11*, 5046. https://doi.org/10.3390/su11185046.

10. Mintel. *UK Meat-Free Foods Market Report*; Mintel: London, UK, 2020.

11. Mäkinen, O.E.; Wanhalinna, V.; Zannini, E.; Arendt, E.K. Foods for Special Dietary Needs: Non-dairy Plant-based Milk Substitutes and Fermented Dairy-type Products. *Crit. Rev. Food Sci. Nutr.* **2016**, *56*, 339–349. https://doi.org/10.1080/10408398.2012.761950.
12. Röös, E.; Garnett, T.; Watz, V.; Sjörs, C. *The Role of Dairy and Plant Based Dairy Alternatives in Sustainable Diets*; Swedish University of Agricultural Sciences: Uppsala, Sweden, 2018.

13. Sethi, S.; Tyagi, S.K.; Anurag, R.K. Plant-based milk alternatives an emerging segment of functional beverages: A review. *J. Food Sci. Technol.* **2016**, *53*, 3408–3423. https://doi.org/10.1007/s13197-016-2328-3.

14. Gazan, R.; Brouzes, C.M.C.; Vieux, F.; Maillot, M.; Lluch, A.; Darmon, N. Mathematical Optimization to Explore Tomorrow's Sustainable Diets: A Narrative Review. *Adv. Nutr.* **2018**, *9*, 602–616. https://doi.org/10.1093/advances/nmy049.

15. Maillot, M.; Vieux, F.; Amiot, M.J.; Darmon, N. Individual diet modeling translates nutrient recommendations into realistic and individual-specific food choices. *Am. J. Clin. Nutr.* **2010**, *91*, 421–430. https://doi.org/10.3945/ajcn.2009.28426.

16. AFSSA. *Etude Individuelle Nationale des Consommations Alimentaires 2 (INCA2)*. 2006–2007; AFSSA: Maisons-Alfort, France, 2009.

17. Goldberg, G.R.; Black, A.E.; Jebb, S.A.; Cole, T.J.; Murgatroyd, P.R.; Coward, W.A.; Prentice, A.M. Critical evaluation of energy intake data using fundamental principles of energy physiology: 1. Derivation of cut-off limits to identify under-recording. *Eur. J. Clin. Nutr.* **1991**, *45*, 569–581.

18. Gazan, R.; Barré, T.; Perignon, M.; Maillot, M.; Darmon, N.; Vieux, F. A methodology to compile food metrics related to diet sustainability into a single food database: Application to the French case. *Food Chem.* **2018**, *238*, 125–133. https://doi.org/10.1016/j.foodchem.2016.11.083.

19. Bertoluci, G.; Masset, G.; Gomy, C.; Mottet, J.; Darmon, N. How to Build a Standardized Country-Specific Environmental Food Database for Nutritional Epidemiology Studies. *PLoS ONE* **2016**, *11*, e0150617. https://doi.org/10.1371/journal.pone.0150617.

20. Vieux, F.; Soler, L.-G.; Touazi, D.; Darmon, N. High nutritional quality is not associated with low greenhouse gas emissions in self-selected diets of French adults. *Am. J. Clin. Nutr.* **2013**, *97*, 569–583. https://doi.org/10.3945/ajcn.112.035105.

21. Maillot, M.; Vieux, F.; Delaere, F.; Lluch, A.; Darmon, N. Dietary changes needed to reach nutritional adequacy without increasing diet cost according to income: An analysis among French adults. *PLoS ONE* **2017**, *12*, e0174679. https://doi.org/10.1371/journal.pone.0174679.

22. Lluch, A.; Maillot, M.; Gazan, R.; Vieux, F.; Delaere, F.; Vaudaine, S.; Darmon, N. Individual Diet Modeling Shows How to Balance the Diet of French Adults with or without Excessive Free Sugar Intakes. *Nutrients* **2017**, *9*, 162. https://doi.org/10.3390/nu9020162.

23. Martin, A. *Apports Nutritionnels Conseillés Pour la Population Française*, 3rd ed.; Technique; Lavoisier: Paris, France, 2001.

24. Agostoni, C.; Bresson, J.L.; Fairweather-Tait, S.; Flynn, A.; Golly, I.; Korhonen, H.; Lagiou, P.; Lovik, M.; Marchelli, R.; Martin, A.; Moseley, B.; Neuhäuser-Berthold, M.; Przyrembel, H.; Salminen, S.; Sanz, Y.; Strain, S.J.J.; Strobel, S.; Tetens, I.; Tomé, D.; van Loveren, H.; Verhagen, H. EFSA Panel on Dietetic Products Nutrition and Allergies (NDA) Scientific Opinion on Dietary Reference Values for water. *EFSA J. 2010*, *8*, 1459. https://doi.org/10.2903/j.efsa.2010.1459.

25. Becker, W.; Lyhne, N.; Pedersen, A.N.; Aro, A.; Fogelholm, M.; Phorsdottir, I.; Alexander, J.; Anderssen, S.A.; Meltzer, H.M.; Pedersen, J.I. Nordic Nutrition Recommendations 2004—Integrating nutrition and physical activity. *Scand. J. Nutr.* **2004**, *48*, 178–187. https://doi.org/10.1080/1102680410003794.

26. Amaya-Farfán, J.; Butte, N.F.; Cabrera, M.I.Z.; Garlick, P.; Gilani, G.S.; Jackson, A.; Kishi, K.; Kurpad, A.V.; Millward, J.; Pencharz, P.; Rand, W.M.; Tomé, D.; Torún, B.; Waterlow, J.; Young, V.R.; Zulu, R.M. FAO-WHO-UNU Joint FAO/WHO/UNU Expert Consultation on Protein and amino acid requirements in human nutrition. *FAO/WHO/UNU* **2007**, *935*, 265.

27. World Health Organization. *Diet, Nutrition, and the Prevention of Chronic Diseases: Report of a Joint WHO/FAO Expert Consultation*; WHO: Geneva, Switzerland, 2003.

28. Agence Nationale de Sécurité Sanitaire de L'alimentation de L'environnement et du Travail. *Actualisation des Repères du PNNS: Révision des Repères de Consommations Alimentaires—Rapport D'expertise Collective*; Agence Nationale de Sécurité Sanitaire de L'alimentation de L'environnement et du Travail: Maisons-Alfort, France, 2016.

29. Barre, T.; Vieux, F.; Perignon, M.; Cravedi, J.-P.; Amiot, M.-J.; Micard, V.; Darmon, N. Reaching nutritional adequacy does not necessarily increase exposure to food contaminants: Evidence from a whole-diet modeling approach. *J. Nutr.* **2016**, *146*, 2149–2157. https://doi.org/10.3945/jn.116.234294.

30. Perignon, M.; Masset, G.; Ferrari, G.; Barré, T.; Vieux, F.; Maillot, M.; Amiot, M.-J.; Darmon, N. How low can dietary greenhouse gas emissions be reduced without impairing nutritional adequacy, affordability and acceptability of the diet? A modelling study to guide sustainable food choices. *Public Health Nutr.* **2016**, *19*, 2662–2674. https://doi.org/10.1017/S1368980016000653.

31. Maillot, M.; Monsivais, P.; Drewnowski, A. Food pattern modeling shows that the 2010 Dietary Guidelines for sodium and potassium cannot be met simultaneously. *Nutr. Res.* **2013**, *33*, 188–194. https://doi.org/10.1016/j.nutres.2013.01.004.

32. ANSES Ciqual: Table de composition nutritionnelle des aliments. Available online: https://ciqual.anses.fr/#/cms/telechargement/node/20 (accessed on Sep 10, 2021).

33. Willett, W.; Rockström, J.; Loken, B.; Springmann, M.; Lang, T.; Vermeulen, S.; Garnett, T.; Tilman, D.; DeClerck, F.; Wood, A.; et al. Food in the Anthropocene: The EAT-Lancet Commission on healthy diets from sustainable food systems. *Lancet (Lond. Engl.)* **2019**, *393*, 447–492. https://doi.org/10.1016/S0140-6736(18)31788-4.

34. EUPHA. *Healthy and Sustainable Diets for European Countries*; EUPHA: Lueneburg, Germany, 2017.

35. Barré, T.; Perignon, M.; Gazan, R.; Vieux, F.; Micard, V.; Amiot, M.-J.; Darmon, N. Integrating nutrient bioavailability and co-production links when identifying sustainable diets: How low should we reduce meat consumption? *PLoS ONE* **2018**, *13*, e0191767. https://doi.org/10.1371/journal.pone.0191767.

36. Vieux, F.; Perignon, M.; Gazan, R.; Darmon, N. Dietary changes needed to improve diet sustainability: Are they similar across Europe? *Eur. J. Clin. Nutr.* **2018**, *72*, 951–960. https://doi.org/10.1038/s41430-017-0080-z.
37. Broekema, R.; Tyszler, M.; van't Veer, P.; Kok, F.J.; Martin, A.; Lluch, A.; Blonk, H.T.J. Future-proof and sustainable healthy diets based on current eating patterns in the Netherlands. *Am. J. Clin. Nutr.* **2020**, *112*, 1338–1347.
https://doi.org/10.1093/ajcn/nqaa217.

38. HCSP. Statement Related to the Revision of the 2017–2021 French Nutrition and Health Programme’s Dietary Guidelines for Adults; HCSP: Paris, France, 2017.

39. Vieux, F.; Privet, L.; Soler, L.G.; Irz, X.; Ferrari, M.; Sette, S.; Raulio, S.; Tapanainen, H.; Hoffmann, R.; Surry, Y.; et al. More sustainable European diets based on self-selection do not require exclusion of entire categories of food. *J. Clean. Prod.* **2020**, *248*, 119298. https://doi.org/10.1016/j.jclepro.2019.119298.

40. Ridoutt, B.G.; Baird, D.; Hendrie, G.A. The role of dairy foods in lower greenhouse gas emission and higher diet quality dietary patterns. *Eur. J. Nutr.* **2021**, *60*, 275–285. https://doi.org/10.1007/s00394-020-02245-w.

41. Clegg, M.E.; Tarrado Ribes, A.; Reynolds, R.; Kliem, K.; Stergiadis, S. A comparative assessment of the nutritional composition of dairy and plant-based dairy alternatives available for sale in the UK and the implications for consumers’ dietary intakes. *Food Res. Int.* **2021**, *148*, 110586. https://doi.org/10.1016/J.FOODRES.2021.110586.

42. Maillot, M.; Vieux, F.; Ferguson, E.F.; Volatier, J.-L.; Amiot, M.J.; Darmon, N. To Meet Nutrient Recommendations, Most French Adults Need to Expand Their Habitual Food Repertoire. *J. Nutr.* **2009**, *139*, 1721–1727. https://doi.org/10.3945/jn.109.107318.

43. Horgan, G.W.; Perrin, A.; Whybrow, S.; Macdiarmid, J.I. Achieving dietary recommendations and reducing greenhouse gas emissions: Modelling diets to minimise the change from current intakes. *Int. J. Behav. Nutr. Phys. Act.* **2016**, *13*, 46. https://doi.org/10.1186/s12966-016-0370-1.

44. French Agency for Ecological Transition (ADEME). Agribalyse V3. Available online: https://ecolab.ademe.fr/agribalyse%0A%0A (accessed on 9 October 2020).

45. Poore, J.; Nemecek, T. Reducing food’s environmental impacts through producers and consumers. *Science* **2018**, *360*, 987–992.
https://doi.org/10.1126/science.aaq0216.

46. Gazan, R. Development of a Tailored Nutrition Counselling Approach for Individuals, Taking into Account Several Dimensions of Diet Sustainability; Université de Montpellier: Montpellier, France, 2018.

47. Equipe de Surveillance et D'épidémiologie Nutritionnelle (Esen). Étude de Santé sur L'environnement, la Biosurveillance, L'activité Physique et la Nutrition (Esteban), 2014–2016. Volet Nutrition. Chapitre Consommations Alimentaires; Equipe de Surveillance et D'épidémiologie Nutritionnelle: Saint-Maurice, France, 2018.

48. French Agency for Food Environmental and Occupational Health & Safety (ANSES). Étude Individuelle Nationale des Consommations Alimentaires 3 (INCA 3); French Agency for Food Environmental and Occupational Health & Safety: Maisons-Alfort, France, 2017.

49. ANSES. Données de Consommations et Habitudes Alimentaires de L'étude INCA2. Available online: https://www.data.gouv.fr/fr/datasets/donnees-de-consommations-et-habitudes-alimentaires-de-letude-inca-2-3/ (accessed on 31 January 2022).