# Effects of first-generation in utero exposure to diesel engine exhaust on second-generation placental function, fatty acid profiles and foetal metabolism in rabbits: preliminary results

Delphine Rousseau-Ralliard, Sarah A. Valentino, Marie-Christine Aubrière, Michèle Dahirel, Marie Sylvie Lallemand, Catherine Archilla, Luc Jouneau, Natalie Fournier, Christophe Richard, Josiane Aioun, et al.

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

## ► To cite this version:

Delphine Rousseau-Ralliard, Sarah A. Valentino, Marie-Christine Aubrière, Michèle Dahirel, Marie Sylvie Lallemand, et al. Effects of first-generation in utero exposure to diesel engine exhaust on second-generation placental function, fatty acid profiles and foetal metabolism in rabbits: preliminary results. Scientific Reports, 2019, 9 (1), pp.9710. 10.1038/s41598-019-46130-x. hal-02185870

HAL Id: hal-02185870

https://hal.science/hal-02185870

Submitted on 16 Jul 2019

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>
SCIENTIFIC REPORTS

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

OPEN

# Effects of first-generation in utero exposure to diesel engine exhaust on second-generation placental function, fatty acid profiles and foetal metabolism in rabbits: preliminary results

Received: 31 May 2018

Accepted: 19 June 2019

Published online: 04 July 2019

Delphine Rousseau-Ralliard<sup>1,2</sup>, Sarah A. Valentino<sup>1,2</sup>, Marie-Christine Aubrière<sup>1,2</sup>, Michèle Dahirel<sup>1,2</sup>, Marie-Sylvie Lallemand<sup>1,2</sup>, Catherine Archilla<sup>1</sup>, Luc Jouneau<sup>1</sup>, Natalie Fournier<sup>3,4</sup>, Christophe Richard<sup>1,2</sup>, Josiane Aioun<sup>1,2</sup>, Anaïs Vitorino Carvalho<sup>1</sup>, Lecardonnel Jérôme<sup>5</sup>, Rémy Slama<sup>6</sup>, Véronique Duranthon<sup>1</sup>, Flemming R. Cassee<sup>7,8</sup>, Pascale Chavatte-Palmer<sup>1,2</sup> & Anne Couturier-Tarrade<sup>1,2</sup>

Atmospheric pollution has major health effects on directly exposed subjects but intergenerational consequences are poorly characterized. We previously reported that diesel engine exhaust (DE) could lead to structural changes in the placenta of *in utero* exposed rabbits (first generation, F1). The effects of maternal exposure to DE were further studied on second-generation (F2) rabbits. Pregnant F0 females were exposed to filtered, diluted DE (1 mg/m<fnref n="3" />, median particle diameter: 69 nm) or clean filtered air (controls) for 2 h/day, 5 days/week by nose-only exposure during days 3–27 post-conception (dpc). Adult female offspring (F1) were mated to control males: F1 tissues and F2 foeto-placental units were collected at 28 dpc and placental structure and gene expression (microarray) analysed. Fatty acid profiles were determined in foetal and maternal plasma, maternal liver and placenta. In F1, compared to controls, hepatic neutral lipid contents were increased in exposed animals without change in the blood biochemistry. In F2, the placental lipid contents were higher, with higher monounsaturated fatty acids and reduced pro-inflammatory arachidonic acid (AA), without placental structural changes. Conversely, the proportion of anti-inflammatory n-3 polyunsaturated fatty acids in F2 plasma was increased while that of AA was decreased. Gene set enrichment analyses (GSEA) of F2 placenta transcriptomic data identified that the proteasome complex and ubiquitin pathways genes were over-represented and ion channel function and inflammation pathways genes were under-represented in exposed animals. These preliminary results demonstrate that diesel engine exhaust exposure and *in utero* indirect exposure should be considered as a programming factor within the context of the DOHaD (Developmental Origins of Health and Disease) with a probable intergenerational transmission.

<footnote n="1"><sup>1</sup>UMR BDR, INRA, ENVA, Université Paris Saclay, Jouy en Josas, France. <sup>2</sup>PremUp Foundation, Paris, France.</footnote>

<footnote n="3"><sup>3</sup>University Paris-Sud, EA 4041/4529 Lip (Sys)2, UFR de Pharmacie, Châtenay-Malabry, France. <sup>4</sup>Hôpital Européen Georges Pompidou (AP-HP), Laboratoire de Biochimie, UF Cardio-Vasculaire, Paris, France. <sup>5</sup>GABI CRB GADIE, INRA, Université Paris Saclay, Jouy en Josas, France. <sup>6</sup>Inserm, Univ. Grenoble Alpes, CNRS, IAB joint Research Center, Team of Environmental Epidemiology Applied to Reproduction and Respiratory Health, Grenoble, France. <sup>7</sup>Centre for Sustainability, Environment and Health, National Institute for Public Health and the Environment, Bilthoven, Netherlands. <sup>8</sup>Institute of Risk Assessment Sciences, Utrecht University, Utrecht, Netherlands. Delphine Rousseau-Ralliard and Sarah A. Valentino contributed equally. Pascale Chavatte-Palmer and Anne Couturier-Tarrade jointly supervised this work. Correspondence and requests for materials should be addressed to D.R.-R. (email: delphine.rousseau@inra.fr) or A.C.-T. (email: anne.couturier-tarrade@inra.fr)</footnote>
Air pollution has become a major health concern in most industrialized and non-industrialized countries over the last decades¹. Diesel exhaust (DE) is a major contributor to air pollution worldwide. DE comprised a large variety of compounds, including a solid particulate fraction that mainly consists of soot, and liquid and gaseous fractions². Current European and North American regulation services limit the diesel motor vehicle emission of carbon monoxide, nitric oxides, hydrocarbons, and particles (size and concentration). Thus, modern motors are equipped with particle filters that reduce PM10 (particulate matter with an aerodynamic diameter below 10 µm) and PM2.5 (<2.5 µm) emissions. Ultrafine or nanoparticles (<100 nm), however, may not be efficiently filtered from the exhaust³ and their impact on human health needs to be evaluated. Several studies in humans and animals have associated reduced foetal growth, low birth weight and long-term offspring outcomes with *in utero* exposure to air pollutants, such as DE⁴⁻¹⁰. Several studies have indicated that placental function is affected by air pollutants¹¹⁻¹³, probably through inflammation and oxidative stress. Comparisons between studies, however, remain difficult due to the large array of exposure routes and windows of exposure as well as the pollutant and particle concentrations used.

The maternal environment, extended across the full range of maternal conditions, is an important component of an offspring's risk of developing metabolic disease in adulthood. These maternal effects are based on the epigenetic memory of early developmental conditions and there is some evidence for their transmission to succeeding generations¹⁴,¹⁵.

Experimental studies of the intergenerational effects of exposure to atmospheric pollution have mostly been conducted in rodents, and mainly focused on offspring reproductive function. Thus, in rodents, maternal (F0) exposure to DE or nanoparticles was shown to reduce daily sperm production and alter steroidogenesis in first generation (F1) males¹⁶. In females, the ovarian reserve was reported to be decreased (reduced number of primary follicles) after maternal exposure to DE¹⁷. These results suggest that offspring gonads and gametes may be affected, thus possibly inducing second-generation effects, as described with other environmental factors¹⁸,¹⁹. In any case, little is known about intergenerational effects of DE exposure. One Brazilian study showed that in mice, continuous exposure to urban San Paulo pollution over 3 generations altered second generation female fertility and increased the incidence of foetal abortion, decreased foetal weight and altered placental structure in the third generation⁸,⁹.

Here, the intergenerational effects of *in utero* F1 exposure to DE on the second (F2) generation were studied using a rabbit model. The rabbit was chosen as a particularly appropriate model because compared to rodent placenta the haemochorial placental structure of a rabbit is closer to that of a human²⁰,²¹. The data presented here build upon previously published data generated in the same experiment, demonstrating that F2 foeto-placental biometry was not disturbed consecutively to the *in utero* exposure of the F1 generation to DE⁷. Nevertheless, total plasma cholesterol and non-HDL cholesterol were reduced (−25.9% and −26.0%, respectively) and plasma triglycerides were increased (+25.9%) in the macroscopically normal F2 foetuses of the exposed group compared to those of the control group⁷.

In the present study, F1 adult end-gestational metabolism (with a focus on liver) and placental function in the second (F2) generation were explored at the biochemical and molecular levels to understand how maternal metabolism through the placenta contributed to foetal dyslipidaemia without affecting biometry.

## Results

In this study, hereafter, CC refers to F2 foeto-placental units collected from F1 control (C) dams, and EC refers to F2 foeto-placental units collected from *in utero*-exposed F1 (E) dams, with no further exposure of these F1 females after birth until adulthood and during their own gestation. In the previously published part of this study by Valentino *et al.*, 2016, although F2 placental and foetal biometry were not different between CC and EC gestations, significant differences were found for F2 foetal plasma biochemistry as described above in the Introduction⁷.

**Maternal biochemistry and histological analyses of the liver at 28 dpc.** In contrast with the findings in the foetuses, there was no significant difference for any studied plasma biochemical parameter in the F1 adult dams at 28 dpc (Supplementary Table S1). Nevertheless, using a histological approach, liver steatosis was observed in exposed compared to control dams (Supplementary Fig. S1). The hepatocytes in pregnant exposed F1 adult rabbits at 28 dpc presented numerous intracellular micro and macro vacuoles, corresponding to lipid vesicles, which were not positive for Periodic Acid Schiff (PAS, glycogen staining), compared to those in pregnant control F1 females.

**Placental stereology.** No change in F2 foetal or placental weight was reported previously⁷. As a reminded, the litter sizes ranged from 6 to 13 fetuses without difference in sex ratio or change in fetal weight between groups. The cellular composition of the labyrinthine area was explored by stereology (Supplementary Fig. S2) and there was no significant difference between the two groups for the relative volume fraction and surface density of the maternal blood space, trophoblast and foetal vessels (Supplementary Table S2).

**Fatty acid analyses: fatty acid concentrations.** Maternal (F1) and foetal (F2) plasma fatty acid concentrations and maternal liver and placenta tissue FA concentrations were determined (Table 1). Phospholipids (PL, corresponding to cell membranes) and neutral lipids (NL, corresponding to intracellular lipid storage) were determined in the placentas and maternal livers. In the placentas, the intracellular neutral lipid concentration was significantly increased in EC compared to that in CC (+95.1%, p = 0.0002) whereas there was no difference between the groups for phospholipid fatty acid concentrations.

Because EC placentas were enriched in fats, the fatty acid profiles of F1 maternal and F2 foetal plasma were analysed. The fatty acid concentrations in foetal plasma were significantly increased (+25.2%, p = 0.021) in EC
<table><caption>**Table 1.** Fatty acid concentrations in plasma and liver from control (C) or exposed (E) F1 dams and in F2 fetal plasma and placentas (CC and EC, respectively) at 28 dpc. *Data are expressed as Median[Q1;Q3]. C and CC: F1 and F2 control samples, respectively; CI: 95% Confidence Interval; E and EC: F1 and F2 exposed samples, respectively; NL: neutral lipids for intracellular lipid storage; PL: phospholipids for cell membranes* Statistics with *P < 0.05* and **P < 0.01.</caption><thead><tr><th rowspan="2">Maternal tissues</th><th colspan="2">Effective</th><th colspan="2">Median [Q1; Q3]</th><th colspan="3">Fully adjusted linear model</th></tr><tr><th>C</th><th>E</th><th>Control</th><th>Exposed</th><th>β value</th><th>CI</th><th>p-value</th></tr></thead><tbody><tr><td>Maternal F1 plasma (µg/mL)</td><td>3</td><td>9</td><td>396.4 [325.4; 463.6]</td><td>390.2 [359.5; 454.0]</td><td>7.03</td><td>[-98.21; 112.26]</td><td>0.901</td></tr><tr><td>Maternal F1 liver (NL) (mg/g of liver)</td><td>3</td><td>9</td><td>44.7 [41.1; 65.4]</td><td>107.6 [89.9; 125.5]</td><td>56.37</td><td>[16.12; 96.61]</td><td>0.050*</td></tr><tr><td>Maternal F1 liver (PL) (mg/g of liver)</td><td>3</td><td>9</td><td>18.8 [14.1; 21.2]</td><td>18.5 [15.6; 23.4]</td><td>-3.43</td><td>[-10.49; 3.64]</td><td>0.391</td></tr><tr><th rowspan="2">Fetal tissues</th><th colspan="2">CC</th><th colspan="2">EC</th><th colspan="2">β value</th><th colspan="2">CI</th></tr><tr><th>CC</th><th>EC</th><th>CC</th><th>EC</th><th>β value</th><th>CI</th><th>p-value</th><th></th></tr><tr><td>F2 Placenta (NL) (mg/g of placenta)</td><td>11</td><td>11</td><td>1.7 [1.5; 1.9]</td><td>3.4 [2.8; 4.2]</td><td>1.74</td><td>[1.26; 2.22]</td><td>0.0002***</td><td></td></tr><tr><td>F2 Placenta (PL) (mg/g of placenta)</td><td>12</td><td>12</td><td>3.3 [2.9; 3.7]</td><td>5.4 [3.4; 6.4]</td><td>2.03</td><td>[-0.16; 4.22]</td><td>0.109</td><td></td></tr><tr><td>F2 Fetal plasma (µg/mL)</td><td>12</td><td>12</td><td>385.4 [346.7; 407.9]</td><td>482.7 [460.9; 577.8]</td><td>113.84</td><td>[41.58; 186.09]</td><td>0.021*</td><td></td></tr></tbody></table>

Foetuses compared to those of CC. Maternal plasma fatty acid concentrations were similar between groups, confirming the biochemical observations. Because maternal plasma fatty acid concentrations were not different among groups, the maternal liver was further explored to understand the increase in the fatty acid content in both foetal plasma and placenta. Thus, intracellular neutral lipid concentrations were significantly increased in the livers of exposed dams compared to those of control dams (p = 0.049), confirming the histological observations of liver steatosis in exposed F1 females once pregnant.

**Qualitative fatty acid analyses: fatty acid profiles.** *Fatty acid profile in F1 maternal plasma.* The proportion of palmitic acid (PA, C16:0) tended to be higher (+9.3%, p = 0.067) whereas C16:1n-7 and C18:3n-3 were decreased (−39.1%, p = 0.027 and −36.7%, p = 0.019, respectively) in the maternal plasma of exposed F1 dams compared to that in the controls (Table 2). Polyunsaturated fatty acids (PUFA), especially those from n-3 series, tended to be reduced in exposed compared to control dams (−27.8%, p = 0.068). PCA confirmed these differences but no specific fatty acid signature related to their *in utero* exposure could be identified (Supplementary Fig. S3).

*Fatty acid profile in F2 placental membranes (phospholipids, PL).* Compared to controls, EC placental membranes (PL) (Table 3) were characterized by lower C15:0 and SFA (−14.4%, p = 0.025 and −1.6%, p = 0.024, respectively) contents, as well as lower long chain polyunsaturated fatty acids (LC n-6 PUFA, ≥20 carbons), especially C20:2n-6 and C22:4n-6 (−16.1%, p = 0.017 and −28.4%, p = 0.004) and LC n-3 PUFA, particularly docosapentaenoic acid (DPA, C22:5n-3, −8%, p = 0.031). The placental membranes were also characterized by higher monounsaturated fatty acid (MUFA) contents, especially C16:1n-9 (+21%, p = 0.02) compared to the placental membranes of controls. PCA performed on placental phospholipid fatty acid profiles confirmed these differences. A specific fatty acid signature enabling the discrimination of each group in placental membranes was identified (Supplementary Fig. S4), with the main plan representing 43.8% of the inertia of the data summarized in Table 3.

*Fatty acid profile in F2 placental intracellular lipid storages (neutral lipids, NL).* EC placental intracellular lipid stores (neutral lipids, NL) (Table 4) were characterized by lower C14:0 and C15:0 contents (−14.2%, p = 0.004 and −12.1%, p = 0.035), a tendency for lower di-homo-gamma-linolenic acid contents (DGLA, C20:3n-6, −17.5%, p = 0.064), as well as lower PUFA contents (−17.3%, p = 0.014), especially arachidonic acid content (AA, C20:4n-6, −31.7%, p = 0.026); they were also characterized by higher MUFA (+12.7%, p = 0.02), particularly oleic acid (OA, C18:1n-9, +17.5%, p = 0.035), compared to the intracellular lipid stores of control placentas. PCA was performed on the fatty acid profiles of placental neutral lipids. A specific signature enabling discrimination between the two groups was identified, with the main plan representing 50.7% of the inertia of the data (Supplementary Fig. S5).

*Fatty acid profile of F2 foetal plasma.* Compared to controls, F2 EC foetuses (Table 5) were characterized by a tendency towards lower linoleic acid (LA, C18:2n-6, −18.4%, p = 0.065) and MUFA (namely, C20:1 and C22:1) contents, higher C22 n-3 PUFA contents, and higher C20:2n-6 (+85.1%, p = 0.059) and C22 n-6 PUFA contents. These analyses of foetal plasma fatty acids were confirmed by PCA since exposed and control F2 foetuses were discriminated by their fatty acid specific signatures (i.e., profiles) on the corresponding PCA, representing 61.6% of the inertia of the data in the main plan (Supplementary Fig. S6).

**Placental transcriptome.** Hierarchical clustering revealed limited differences in the F2 placental transcriptome between the EC and CC groups, as shown by the very short clustering scale (supplementary Fig. S7). This close similarity between F2 EC and CC placentas was confirmed by the small number of differentially expressed (DE) probes, transcripts and annotated genes (Supplementary Table S3). The microarray data have been registered in the Gene Expression Omnibus (GEO) under accession number GSE109831.

Differential gene expression (adjusted p-value < 0.05) between EC and CC placentas was observed for F-box protein 32 (FBXO32, +119%, 4 significant probes on a total of 5 probes in microarray), T-cell activation
<table><caption>**Table 2.** Fatty acid profiles in plasma of control (C) and *in utero* exposed (E) F1 adult pregnant dams at 28dpc (% of total fatty acids). *Data are expressed as Median [Q1; Q3], in % of total FA.* AA: C20:4n-6, arachidonic acid; ALA: C18:3n-3, alpha-linolenic acid; CI: Confidence Interval; DHA: C22:6n-3, docosahexaenoic acid; DNL index: de novo lipogenesis (C16:0/C18:2n-6 ratio)²,²³; DPA: C22:5n-3, docosapentaenoic acid; EPA: C20:5n-3, Eicosapentaenoic acid; LA: C18:2n-6, linoleic acid; MUFA: monounsaturated fatty acids; PA: C16:0, palmitic acid; PL: phospholipids; PUFA: polyunsaturated fatty acids; SA: C18:0, stearic acid; SFA: saturated fatty acids "Statistics with *P* < 0.05" and **P** < 0.01.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="2">Number of dams</th><th colspan="2">Median [Q1; Q3]</th><th colspan="3">Fully adjusted linear model</th></tr><tr><th>C</th><th>E</th><th>Control</th><th>Exposed</th><th>β value</th><th>CI</th><th>p-value</th></tr></thead><tbody><tr><td>C16:0 (PA)</td><td>3</td><td>9</td><td>29.1 [28.5; 30.1]</td><td>31.8 [30.4; 32.8]</td><td>2.55</td><td>[0.53; 4.57]</td><td>0.067</td></tr><tr><td>C16:1n-7</td><td>3</td><td>9</td><td>10.4 [8.2; 11.5]</td><td>6.4 [6.1; 8.2]</td><td>-3.15</td><td>[-4.98; -1.31]</td><td>0.027*</td></tr><tr><td>Total C16:1</td><td>3</td><td>9</td><td>10.9 [8.6; 11.9]</td><td>7.1 [6.6; 8.6]</td><td>-3.10</td><td>[-4.91; -1.28]</td><td>0.028*</td></tr><tr><td>C18:0 (SA)</td><td>3</td><td>9</td><td>8.0 [6.0; 9.2]</td><td>10.0 [8.7; 12.7]</td><td>3.73</td><td>[-0.90; 8.36]</td><td>0.186</td></tr><tr><td>Total C18:1</td><td>3</td><td>9</td><td>25.2 [25.0; 28.8]</td><td>25.4 [24.2; 29.5]</td><td>-1.03</td><td>[-8.21; 6.15]</td><td>0.791</td></tr><tr><td>C18:2 n-6 (LA)</td><td>3</td><td>9</td><td>17.7 [16.9; 17.9]</td><td>16.2 [15.4; 17.3]</td><td>-1.25</td><td>[-2.44; -0.06]</td><td>0.106</td></tr><tr><td>C18:3 n-3 (ALA)</td><td>3</td><td>9</td><td>2.8 [2.5; 3.3]</td><td>1.8 [1.6; 2.0]</td><td>-1.22</td><td>[-1.85; -0.59]</td><td>0.019*</td></tr><tr><td>C20:4 n-6 (AA)</td><td>3</td><td>9</td><td>1.5 [0.9; 1.6]</td><td>1.0 [0.8; 1.7]</td><td>-0.06</td><td>[-0.84; 0.72]</td><td>0.892</td></tr><tr><td>C20:5 n-3 (EPA)</td><td>3</td><td>9</td><td>0.02 [0.00; 0.04]</td><td>0.04 [0.02; 0.37]</td><td>0.18</td><td>[-0.34; 0.71]</td><td>0.529</td></tr><tr><td>C22:4 n-6</td><td>3</td><td>9</td><td>0.04 [0.04; 0.11]</td><td>0.04 [0.03; 0.06]</td><td>-0.01</td><td>[-0.08; 0.06]</td><td>0.794</td></tr><tr><td>C22:5 n-6</td><td>3</td><td>9</td><td>0.08 [0.00; 0.08]</td><td>0.06 [0.05; 0.11]</td><td>0.04</td><td>[-0.05; 0.14]</td><td>0.397</td></tr><tr><td>C22:5 n-3 (DPA)</td><td>3</td><td>9</td><td>0.07 [0.05; 0.08]</td><td>0.05 [0.02; 0.08]</td><td>-0.01</td><td>[-0.07; 0.05]</td><td>0.710</td></tr><tr><td>C22:6 n-3 (DHA)</td><td>3</td><td>9</td><td>0.05 [0.04; 0.07]</td><td>0.04 [0.03; 0.06]</td><td>-0.01</td><td>[-0.05; 0.03]</td><td>0.608</td></tr><tr><td>Total SFA</td><td>3</td><td>9</td><td>39.9 [38.7; 40.0]</td><td>44.9 [41.5; 47.5]</td><td>6.37</td><td>[-0.98; 13.71]</td><td>0.162</td></tr><tr><td>Total MUFA</td><td>3</td><td>9</td><td>37.8 [36.7; 37.9]</td><td>34.6 [31.4; 37.3]</td><td>-4.04</td><td>[-10.34; 2.25]</td><td>0.272</td></tr><tr><td>Total n-6 PUFA</td><td>3</td><td>9</td><td>19.7 [19.4; 20.1]</td><td>18.9 [17.6; 19.3]</td><td>-1.16</td><td>[-3.05; 0.72]</td><td>0.289</td></tr><tr><td>Total n-3 PUFA</td><td>3</td><td>9</td><td>3.2 [2.9; 3.7]</td><td>2.3 [2.1; 2.8]</td><td>-1.11</td><td>[-2.01; -0.22]</td><td>0.068</td></tr><tr><td>n-6/n-3 PUFA ratio</td><td>3</td><td>9</td><td>6.4 [5.3; 6.8]</td><td>8.2 [6.6; 9.5]</td><td>2.64</td><td>[-0.59; 5.88]</td><td>0.181</td></tr><tr><td>PUFA/SFA ratio</td><td>3</td><td>9</td><td>0.58 [0.56; 0.60]</td><td>0.50 [0.42; 0.52]</td><td>-0.11</td><td>[-0.19; -0.03]</td><td>0.049*</td></tr><tr><td>DNL index</td><td>3</td><td>9</td><td>1.7 [1.6; 1.7]</td><td>1.9 [1.8; 2.1]</td><td>0.30</td><td>[0.06; 0.53]</td><td>0.070</td></tr></tbody></table>

RhoGTPase activating protein (TAGAP, +214%, 2 significant probes on a total of 2 probes) and Interleukin 23 α subunit p19 (IL23A, +275%, 1 significant probe on a total of 1 probe). After RT-qPCR, statistical significance was only found for one of the three transcripts highlighted by the microarray, FBXO32, confirming this result. For the other two transcripts, TAGAP and IL23A, only a trend going in the same direction was observed. (Fig. 1).

*Gene set enrichment analyses (GSEA) in EC placentas.* Because the usual strategy used to analyse transcriptomics is based on a single-gene expression level that becomes inefficient with subtle variations in the same pathway²⁴, we further analysed placental transcriptomic data using a GSEA approach as described above. Using GSEA, a total of 58 significantly enriched gene sets were observed after FDR multiple test adjustment (q-value < 0.25) for the EC group, including 10 terms for the KEGG (the Kyoto Encyclopedia of Genes and Genomes) database, 14 terms for the biological process (BP) category, 19 terms for the cellular component (CC) category and 15 terms for the molecular function (MF) category of the GO (Gene Ontology) database (Supplementary Tables S4–S7). In the top fifteen terms for all databases, using an enrichment score ranking, most terms defining the placenta-exposed group were related to intracellular organization (nuclear pore and envelope, "integral to organelle membrane", "organelle membrane" and "Golgi vesicle transport") and protein metabolism ("signal sequence binding", "ubiquitin protein ligase activity", "ubiquitin mediated proteolysis" and "proteasome complex") (Fig. 2).

*Gene set enrichment in CC placentas.* A total of 69 significantly enriched gene sets were observed after FDR multiple test adjustment in the control group (i.e., comparatively downregulated pathways in EC placentas), including 16 terms in the KEGG database, 18 terms in the BP, 11 terms in the CC and 24 terms in the MF categories of the GO database (Supplementary Tables S4–S7). In the top fifteen gene sets enriched, most terms were involved in olfactory transduction and neuroactive ligand or in membrane and intracellular pathways (G protein coupled receptor activity, G protein coupled receptor binding, amine receptor activity, cation channel activity, substrate specific channel activity and ion channel activity) (Fig. 2).

Discussion

This study was designed to explore maternal metabolism and placental function through several approaches to explain the previously observed foetal dyslipidaemia without biometric changes in F2 foetuses whose grandmothers were exposed to DE during gestation⁷. The two 1-hour daily exposures in F0 were used to mimic our home-work journeys morning and evening (corresponding to real life for outdoor air pollution exposure) particularly during the period of peak pollution in urbanized cities. The concentration of 1 mg/m³ for 2 hours translates to 85 µg/m³ when averaged over 24 hours, which is closer to the PM10 concentration that increases from
<table><caption>**Table 3.** Profile in fatty acid in F2 of control (CC) and exposed (EC) placenta phospholipids (PL) at 28 dpc (% of total fatty acids). *Data are expressed as Median [Q1;Q3], in % of total FA.* AA: C20:4n-6, arachidonic acid; ALA: C18:3n-3, alpha-linolenic acid; CI: Confidence Interval; DGLA: C20:3n-6, dihomo-gamma-linolenic acid; DHA: C22:6n-3, docosahexaenoic acid; DPA: C22:5n-3, docosapentaenoic acid; EPA: C20:5n-3, Eicosapentaenoic acid; GLA: C18:3n-6, gamma-linolenic acid; LA: C18:2n-6, linoleic acid; MUFA: monounsaturated fatty acids; PA: C16:0, palmitic acid; PL: phospholipids; PUFA: polyunsaturated fatty acids; SA: C18:0, stearic acid; SFA: saturated fatty acids "Statistics with *P* < 0.05" and **P** < 0.01.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="2">Number of placenta</th><th colspan="2">Median [Q1; Q3]</th><th colspan="3">Fully adjusted linear model</th></tr><tr><th>CC</th><th>EC</th><th>CC</th><th>EC</th><th>β value</th><th>CI</th><th>p-value</th></tr></thead><tbody><tr><td>C15:0</td><td>12</td><td>12</td><td>0.52 [0.51; 0.56]</td><td>0.45 [0.42; 0.49]</td><td>-0.09</td><td>[-0.148; -0.027]</td><td>0.025*</td></tr><tr><td>C16:0 (PA)</td><td>12</td><td>12</td><td>29.3 [28.9; 30.0]</td><td>29.1 [28.9; 29.4]</td><td>-0.50</td><td>[-1.514; 0.51]</td><td>0.359</td></tr><tr><td>C16:1n-9</td><td>12</td><td>12</td><td>4.0 [3.6; 4.6]</td><td>4.9 [4.5; 5.0]</td><td>0.68</td><td>[0.234; 1.131]</td><td>0.020*</td></tr><tr><td>Total C16:1</td><td>12</td><td>12</td><td>6.1 [5.5; 6.7]</td><td>7.0 [6.7; 7.2]</td><td>0.83</td><td>[0.35; 1.307]</td><td>0.011*</td></tr><tr><td>C18:0 (SA)</td><td>12</td><td>12</td><td>17.3 [17.0; 17.6]</td><td>17.0 [16.5; 17.3]</td><td>-0.32</td><td>[-0.728; 0.083]</td><td>0.160</td></tr><tr><td>Total C18:1</td><td>12</td><td>12</td><td>19.3 [18.9; 19.7]</td><td>19.4 [18.8; 20.2]</td><td>0.38</td><td>[-0.663; 1.414]</td><td>0.497</td></tr><tr><td>C18:2 n-6 (LA)</td><td>12</td><td>12</td><td>11.8 [10.9; 12.1]</td><td>12.4 [11.8; 13.0]</td><td>0.63</td><td>[-0.051; 1.304]</td><td>0.110</td></tr><tr><td>C18:3n-6 (GLA)</td><td>12</td><td>12</td><td>0.58 [0.50; 0.69]</td><td>0.79 [0.75; 0.87]</td><td>0.20</td><td>[0.065; 0.328]</td><td>0.021*</td></tr><tr><td>C18:3n-3 (ALA)</td><td>12</td><td>12</td><td>0.09 [0.07; 0.09]</td><td>0.08 [0.07; 0.08]</td><td>-0.01</td><td>[-0.016; 0.002]</td><td>0.190</td></tr><tr><td>C20:2n-6</td><td>12</td><td>12</td><td>0.16 [0.14; 0.17]</td><td>0.13 [0.12; 0.14]</td><td>-0.03</td><td>[-0.048; -0.011]</td><td>0.017*</td></tr><tr><td>C20:3n-6 (DGLA)</td><td>12</td><td>12</td><td>5.2 [4.9; 5.6]</td><td>4.8 [4.6; 5.1]</td><td>-0.45</td><td>[-1.144; 0.238]</td><td>0.236</td></tr><tr><td>C20:4n-6 (AA)</td><td>12</td><td>12</td><td>2.9 [2.8; 3.1]</td><td>2.5 [2.4; 2.9]</td><td>-0.19</td><td>[-0.444; 0.062]</td><td>0.186</td></tr><tr><td>C20:5n-3 (EPA)</td><td>12</td><td>12</td><td>0.04 [0.03; 0.05]</td><td>0.04 [0.04; 0.04]</td><td>0.002</td><td>[-0.009; 0.012]</td><td>0.765</td></tr><tr><td>C22:4n-6</td><td>12</td><td>12</td><td>0.74 [0.67; 0.85]</td><td>0.53 [0.49; 0.62]</td><td>-0.21</td><td>[-0.313; -0.112]</td><td>0.004**</td></tr><tr><td>C22:5n-6</td><td>12</td><td>12</td><td>0.86 [0.80; 1.02]</td><td>0.70 [0.63; 0.90]</td><td>-0.17</td><td>[-0.417; 0.069]</td><td>0.200</td></tr><tr><td>C22:5n-3 (DPA)</td><td>12</td><td>12</td><td>0.69 [0.65; 0.78]</td><td>0.64 [0.50; 0.65]</td><td>-0.13</td><td>[-0.233; -0.035]</td><td>0.031*</td></tr><tr><td>C22:6n-3 (DHA)</td><td>12</td><td>12</td><td>0.26 [0.22; 0.27]</td><td>0.27 [0.22; 0.32]</td><td>0.01</td><td>[-0.054; 0.078]</td><td>0.732</td></tr><tr><td>Total SFA</td><td>12</td><td>12</td><td>49.3 [48.8; 49.8]</td><td>48.5 [47.8; 49.0]</td><td>-0.77</td><td>[-1.308; -0.241]</td><td>0.024*</td></tr><tr><td>Total MUFA</td><td>12</td><td>12</td><td>27.3 [26.2; 27.7]</td><td>28.2 [27.5; 29.3]</td><td>1.33</td><td>[0.165; 2.491]</td><td>0.058</td></tr><tr><td>Total n-6 PUFA</td><td>12</td><td>12</td><td>22.2 [27.8; 23.0]</td><td>22.1 [21.3; 23.1]</td><td>-0.64</td><td>[-1.814; 0.541]</td><td>0.320</td></tr><tr><td>Total n-3 PUFA</td><td>12</td><td>12</td><td>1.2 [1.1; 1.2]</td><td>1.1 [0.9; 1.2]</td><td>-0.13</td><td>[-0.279; 0.019]</td><td>0.127</td></tr><tr><td>n-6/n-3 PUFA ratio</td><td>12</td><td>12</td><td>19.2 [18.1; 19.7]</td><td>19.1 [18.6; 24.2]</td><td>1.88</td><td>[-1.714; 5.464]</td><td>0.335</td></tr><tr><td>PUFA/SFA ratio</td><td>12</td><td>12</td><td>0.47 [0.46; 0.49]</td><td>0.48 [0.47; 0.49]</td><td>-0.004</td><td>[-0.02; 0.012]</td><td>0.653</td></tr></tbody></table>

50 µg/m³ (European daily standard) to 100 µg/m³ and even higher during typical smog episodes in urbanized cities worldwide. The effects of *in utero* exposure to DE on the phenotype of exposed F1 rabbits were expected at adulthood according to the literature²⁵. No significant metabolic modifications were observed in F1 females before pregnancy, but a trend towards hypertriglyceridemia only was expected (unpublished data); this finding was confirmed at 28 dpc since these pregnant F1 animals did not exhibit any difference in blood biochemistry or plasma fatty acid concentrations.

Gestation affects maternal lipid metabolism since maternal plasma free fatty acid (FFA) levels are commonly elevated in late pregnancy, participating in the onset of insulin resistance observed in pregnant women²⁶,²⁷, rodents²⁸ and rabbits²⁹,³⁰. Thus pregnancy could be considered a metabolic challenge that could reveal latently impaired metabolism in rabbit females *in utero* exposed to DE. Since *in utero* exposure to polycyclic aromatic hydrocarbons (PAH), as potential actors in DE, causes hepatic steatosis in the female offspring of mice³¹,³², the liver was analysed in F1 females and hepatic steatosis was observed in exposed animals by histology, but could not be quantified firmly by this approach. Nevertheless, gas chromatography quantitative analyses confirmed that hepatic intracellular fatty acid concentrations were higher in exposed F1 females than in controls. Furthermore, the study of maternal plasma fatty acid profiles indicated that although both groups were given the same diet, plasma samples from exposed F1 females were enriched in C16:0 and mainly depleted in C16:1 and ALA. The increase in C16:0, parallel to decreased C14:0, could indicate higher *de novo* lipogenesis, which results in an increase of the DNL index (trend with p = 0.07) in plasma of F1 exposed dams²²,²³, probably using ALA as substrate³³. The histological evaluation of placental vascularization and the trophoblast volume did not show any difference between EC and CC F2 placentas. The fatty acid concentration of placental membranes (phospholipid fraction) was not different between groups, thus confirming the stereological data. The placental intracellular lipid contents (neutral lipids), however, were greatly increased in EC compared to CC placentas, as also reported in placentas from pregnant rabbits fed a high fat diet³⁴,³⁵ but no effect of foetal sex could be found in any of the fatty acid analyses. Because EC placentas were enriched in fats, fatty acid profiles of F1 maternal and F2 foetal plasma samples were analysed. As expected, the fatty acid concentrations of F2 EC plasma were increased compared to those of CC, confirming previous plasma biochemical analyses showing higher plasma triglycerides in these foetuses⁷. As the placenta is a programming agent of foetal growth that is sensitive to the maternal environment, placental fatty acid profiles were investigated. In the phospholipid fraction of EC placentas, a weak decrease of non-biologically active intermediates of the n-6 polyunsaturated fatty acid (PUFA, namely, C20:2n-6
<table><caption>**Table 4.** Profile in fatty acid in F2 of control (CC) and exposed (EC) placenta NL at 28dpc (% of total fatty acids). *Data are expressed as Median [Q1; Q3], in % of total FA.* AA: C20:4n-6, arachidonic acid; ALA: C18:3n-3, alpha-linolenic acid; CI: Confidence Interval; DGLA: C20:3n-6, dihomo-gamma-linolenic acid; DHA: C22:6n-3, docosahexaenoic acid; DPA: C22:5n-3, docosapentaenoic acid; EPA: C20:5n-3, Eicosapentaenoic acid; LA: C18:2n-6, linoleic acid; MUFA: monounsaturated fatty acids; OA: C18:1n-9, oleic acid; PA: C16:0, palmitic acid; PL: phospholipids; PUFA: polyunsaturated fatty acids; SA: C18:0, stearic acid; SFA: saturated fatty acids "Statistics with *P* < 0.05" and **P* < 0.01.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="2">Number of placenta</th><th colspan="2">Median [Q1; Q3]</th><th colspan="3">Fully adjusted linear model</th></tr><tr><th>CC</th><th>EC</th><th>CC</th><th>EC</th><th>β value</th><th>CI</th><th>p-value</th></tr></thead><tbody><tr><td>C14:0</td><td>11</td><td>11</td><td>2.2 [2.1; 2.4]</td><td>1.9 [1.8; 1.9]</td><td>-0.35</td><td>[-0.51; -0.19]</td><td>0.004**</td></tr><tr><td>C15:0</td><td>11</td><td>11</td><td>0.91 [0.84; 0.94]</td><td>0.80 [0.69; 0.85]</td><td>-0.11</td><td>[-0.19; -0.03]</td><td>0.035*</td></tr><tr><td>C16:0 (PA)</td><td>11</td><td>11</td><td>33.3 [32.4; 34.5]</td><td>34.0 [32.0; 34.8]</td><td>0.36</td><td>[-1.51; 2.23]</td><td>0.715</td></tr><tr><td>Total C16:1</td><td>11</td><td>11</td><td>7.5 [7.0; 7.8]</td><td>7.8 [7.4; 8.0]</td><td>0.81</td><td>[-0.48; 2.10]</td><td>0.253</td></tr><tr><td>C18:0 (SA)</td><td>11</td><td>11</td><td>9.6 [9.2; 10.1]</td><td>9.1 [8.7; 9.9]</td><td>-0.38</td><td>[-0.99; 0.23]</td><td>0.253</td></tr><tr><td>C18:1n-7</td><td>11</td><td>11</td><td>3.1 [3.0; 3.3]</td><td>3.2 [3.2; 3.4]</td><td>0.27</td><td>[0.06; 0.48]</td><td>0.043*</td></tr><tr><td>C18:1n-9 (OA)</td><td>11</td><td>11</td><td>22.6 [22.4; 24.7]</td><td>26.6 [24.1; 27.4]</td><td>3.06</td><td>[0.74; 5.37]</td><td>0.035*</td></tr><tr><td>Total C18:1</td><td>11</td><td>11</td><td>25.8 [25.4; 27.7]</td><td>29.8 [27.4; 30.8]</td><td>3.18</td><td>[1.01; 5.35]</td><td>0.023*</td></tr><tr><td>C18:2n-6 (LA)</td><td>11</td><td>11</td><td>8.9 [8.6; 9.5]</td><td>7.9 [6.7; 8.5]</td><td>-1.62</td><td>[-2.76; -0.48]</td><td>0.026*</td></tr><tr><td>C18:3n-3 (ALA)</td><td>11</td><td>11</td><td>0.49 [0.46; 0.51]</td><td>0.43 [0.36; 0.47]</td><td>-0.07</td><td>[-0.14; -0.01]</td><td>0.057</td></tr><tr><td>C20:3n-6 (DGLA)</td><td>11</td><td>11</td><td>5.0 [4.1; 6.3]</td><td>4.1 [3.7; 4.6]</td><td>-1.04</td><td>[-1.97; -0.10]</td><td>0.064</td></tr><tr><td>C20:4n-6 (AA)</td><td>11</td><td>11</td><td>1.6 [1.4; 1.7]</td><td>1.1 [1.0; 1.3]</td><td>-0.27</td><td>[-0.46; -0.09]</td><td>0.026*</td></tr><tr><td>C20:5n-3 (EPA)</td><td>11</td><td>11</td><td>0.09 [0.08; 0.10]</td><td>0.07 [0.06; 0.08]</td><td>-0.02</td><td>[-0.04; -0.01]</td><td>0.008**</td></tr><tr><td>C22:2n-6</td><td>11</td><td>11</td><td>0.17 [0.14; 0.21]</td><td>0.13 [0.09; 0.14]</td><td>-0.05</td><td>[-0.09; -0.01]</td><td>0.047*</td></tr><tr><td>C22:4n-6</td><td>11</td><td>11</td><td>0.54 [0.38; 0.65]</td><td>0.28 [0.26; 0.42]</td><td>-0.09</td><td>[-0.21; 0.03]</td><td>0.168</td></tr><tr><td>C22:5n-6</td><td>11</td><td>11</td><td>0.17 [0.14; 0.21]</td><td>0.12 [0.11; 0.18]</td><td>-0.04</td><td>[-0.09; 0.02]</td><td>0.256</td></tr><tr><td>C22:5n-3 (DPA)</td><td>11</td><td>11</td><td>0.22 [0.18; 0.24]</td><td>0.16 [0.12; 0.23]</td><td>0.01</td><td>[-0.04; 0.07]</td><td>0.643</td></tr><tr><td>C22:6n-3 (DHA)</td><td>11</td><td>11</td><td>0.04 [0.03; 0.09]</td><td>0.04 [0.03; 0.05]</td><td>-0.03</td><td>[-0.06; 0.002]</td><td>0.104</td></tr><tr><td>Total SFA</td><td>11</td><td>11</td><td>46.3 [44.8; 47.6]</td><td>45.5 [43.7; 47.0]</td><td>-0.03</td><td>[-2.48; 2.41]</td><td>0.978</td></tr><tr><td>Total MUFA</td><td>11</td><td>11</td><td>33.9 [33.3; 35.9]</td><td>38.16 [36.0; 38.9]</td><td>4.02</td><td>[1.35; 6.68]</td><td>0.02*</td></tr><tr><td>Total n-6 PUFA</td><td>11</td><td>11</td><td>18.0 [16.4; 19.4]</td><td>14.8 [14.4; 16.2]</td><td>-3.48</td><td>[-5.55; -1.41]</td><td>0.013*</td></tr><tr><td>Total n-3 PUFA</td><td>11</td><td>11</td><td>1.8 [1.3; 1.9]</td><td>1.7 [1.5; 1.8]</td><td>-0.02</td><td>[-0.22; 0.18]</td><td>0.876</td></tr><tr><td>n-6/n-3 PUFA ratio</td><td>11</td><td>11</td><td>10.5 [9.5; 13.6]</td><td>9.0 [8.3; 10.0]</td><td>-2.07</td><td>[-4.12; -0.03]</td><td>0.085</td></tr><tr><td>PUFA/SFA ratio</td><td>11</td><td>11</td><td>0.45 [0.38; 0.46]</td><td>0.35 [0.33; 0.41]</td><td>-0.07</td><td>[-0.13; -0.02]</td><td>0.039*</td></tr></tbody></table>

and C22:4n-6) was observed, as well as a slight decrease in saturated fatty acid (SFA) while monounsaturated fatty acid (MUFA, namely, C16:1) was increased, leading to a small increase in the PUFA/SFA ratio. These data suggest that membrane fluidity was increased³⁶, possibly leading to increased fatty acid transfers across the placenta³⁷⁻³⁹.

Essential fatty acids to be transferred from dam to foetus, namely, arachidonic (AA, C20:4n-6), eicosapentaenoic (EPA, 20:5n-3) and docosahexaenoic acids (DHA, 22:6n-3), accumulate preferentially in placental phospholipids that act as a regulator of fatty acid uptake in trophoblastic cells prior to being transferred to the foetus⁴⁰. Most potent eicosanoids are derived from AA, but additional biologically significant eicosanoids are also derived from dihomo-gamma-linolenic acid (DGLA), which is produced in the reaction pathway leading to AA from LA, with distinctly different biological actions (less inflammation or even anti-inflammation, vasodilation, and the inhibition of platelet aggregation⁴¹). In EC placentas, where fatty acid concentrations were nearly doubled compared to controls, the proportion of n-6 (C22:4n-6) and n-3 long chain PUFA (DPA and DHA) was increased in neutral lipids (lipid storage), while the percentage of AA was halved, remaining, however, sufficient to cover foetal needs.

In foetuses, plasma long chain PUFA concentrations depend on placental transfers²⁷,⁴²,⁴³. Here, the AA content and n-6/n-3 PUFA ratio were decreased in the plasma of F2 EC fetuses. A high n-6/n-3 ratio usually reflects low-grade inflammation⁴⁴. The increased proportion of n-3 PUFA in foetal plasma could explain the observed lower total and non-HDL cholesterol concentrations in the F2 EC foetal plasma, as observed in mice in a context of atherosclerosis and dietary n-3 PUFA intake⁴⁵. Moreover, since plasma fatty acid concentrations were higher in EC foetuses compared to those of CC, this further increased the F2 EC foetal reserves of n-3 PUFA, among which biologically active PUFA (namely, EPA, DHA), were favoured by selective placental fatty acid transfer, whereas plasma AA was decreased, which in return limited pro-inflammatory processes. Since anti-inflammatory eicosanoids are derived from n-3 PUFA⁴³, these changes in plasma F2 fatty acid profiles could be a positive adaptive response to counteract the effects of perturbed fatty acid metabolism in dams (as discussed above).

The placental transcriptomics analysis showed that only three genes were differentially expressed, namely, *interleukin 23 subunit alpha (IL-23A)*, *T-cell activation Rho GTPase-activating protein (TAGAP)* and *F-Box Protein 32 gene (FBXO32)*, among which only the last gene was confirmed after RT-qPCR.
<table><caption>**Table 5.** Fatty acid profile in plasma in F2 of control (CC) and exposed (EC) fetuses at 28dpc (% of total fatty acids). *Data are expressed as Median [Q1; Q3], in % of total FA. AA: C20:4n-6, arachidonic acid; ALA: C18:3n-3, alpha-linolenic acid; DGLA: C20:3n-6, dihomo-gamma-linolenic acid; DHA: C22:6n-3, docosahexaenoic acid; DPA: C22:5n-3, docosapentaenoic acid; EPA: C20:5n-3, Eicosapentaenoic acid; LA: C18:2n-6, linoleic acid; MUFA: monounsaturated fatty acids; PA: C16:0, palmitic acid; PL: phospholipids; PUFA: polyunsaturated fatty acids; SA: C18:0, stearic acid; SFA: saturated fatty acids* Statistics with *P < 0.05* and **P < 0.01*.</caption><thead><tr><th rowspan="2">Variable</th><th colspan="2">Number of fetus</th><th colspan="2">Median [Q1; Q3]</th><th colspan="3">Fully adjusted linear model</th></tr><tr><th>CC</th><th>EC</th><th>CC</th><th>EC</th><th>β value</th><th>CI</th><th>p-value</th></tr></thead><tbody><tr><td>C15:1n-9</td><td>12</td><td>12</td><td>0.16 [0.12; 0.22]</td><td>0.28 [0.21; 0.29]</td><td>0.10</td><td>[0.04; 0.16]</td><td>0.019*</td></tr><tr><td>C16:0 (PA)</td><td>12</td><td>12</td><td>39.4 [37.7; 40.3]</td><td>35.1 [31.5; 37.7]</td><td>-4.36</td><td>[-7.18; -1.55]</td><td>0.022*</td></tr><tr><td>Total C16:1</td><td>12</td><td>12</td><td>6.3 [6.0; 6.9]</td><td>6.0 [5.3; 6.5]</td><td>-0.27</td><td>[-1.15; 0.60]</td><td>0.558</td></tr><tr><td>C18:0 (SA)</td><td>12</td><td>12</td><td>11.5 [11.2; 11.9]</td><td>10.7 [10.3; 11.3]</td><td>-0.44</td><td>[-2.20; 1.32]</td><td>0.637</td></tr><tr><td>C18:1n-7</td><td>12</td><td>12</td><td>2.0 [1.8; 2.1]</td><td>1.5 [1.3; 1.6]</td><td>-0.39</td><td>[-0.62; -0.17]</td><td>0.013*</td></tr><tr><td>Total C18:1</td><td>12</td><td>12</td><td>22.6 [21.1; 23.1]</td><td>18.5 [17.7; 22.5]</td><td>-1.64</td><td>[-4.83; 1.54]</td><td>0.346</td></tr><tr><td>C20:0</td><td>12</td><td>12</td><td>0.13 [0.11; 0.15]</td><td>0.80 [0.42; 1.15]</td><td>0.58</td><td>[0.18; 0.98]</td><td>0.029*</td></tr><tr><td>C20:1n-9</td><td>12</td><td>12</td><td>0.38 [0.32; 0.46]</td><td>1.33 [0.95; 2.03]</td><td>1.02</td><td>[0.59; 1.44]</td><td>0.003**</td></tr><tr><td>C20:1n-7</td><td>12</td><td>12</td><td>0.07 [0.05; 0.08]</td><td>0.475 [0.28; 0.70]</td><td>0.40</td><td>[0.23; 0.57]</td><td>0.003**</td></tr><tr><td>Total C20:1</td><td>12</td><td>12</td><td>0.46 [0.39; 0.53]</td><td>1.80 [1.23; 2.73]</td><td>1.42</td><td>[0.83; 2.00]</td><td>0.003**</td></tr><tr><td>C18:2 n-6 (LA)</td><td>12</td><td>12</td><td>11.1 [10.5; 12.1]</td><td>9.1 [8.6; 9.9]</td><td>-1.70</td><td>[-3.19; -0.21]</td><td>0.065</td></tr><tr><td>C18:3 n-3 (ALA)</td><td>12</td><td>12</td><td>0.57 [0.44; 0.67]</td><td>0.57 [0.49; 0.57]</td><td>-0.06</td><td>[-0.27; 0.15]</td><td>0.582</td></tr><tr><td>C18:4n-3</td><td>12</td><td>12</td><td>0.18 [0.15; 0.28]</td><td>0.32 [0.24; 0.51]</td><td>0.16</td><td>[0.07; 0.25]</td><td>0.013*</td></tr><tr><td>C20:2n-6</td><td>12</td><td>12</td><td>0.34 [0.29; 0.39]</td><td>0.62 [0.21; 0.93]</td><td>0.29</td><td>[0.04; 0.53]</td><td>0.059</td></tr><tr><td>C20:3n-6 (DGLA)</td><td>12</td><td>12</td><td>0.49 [0.45; 0.57]</td><td>0.36 [0.30; 0.42]</td><td>-0.13</td><td>[-0.20; -0.06]</td><td>0.012*</td></tr><tr><td>C20:4 n-6 (AA)</td><td>12</td><td>12</td><td>2.6 [2.5; 2.8]</td><td>1.6 [1.34; 1.83]</td><td>-1.03</td><td>[-1.37; -0.69]</td><td>0.001**</td></tr><tr><td>C20:4n-3</td><td>12</td><td>12</td><td>0.09 [0.06; 0.15]</td><td>0.19 [0.10; 0.26]</td><td>0.13</td><td>[0.04; 0.22]</td><td>0.032*</td></tr><tr><td>C20:5 n-3 (EPA)</td><td>12</td><td>12</td><td>1.4 [0.01; 1.7]</td><td>1.5 [0.05; 2.6]</td><td>0.02</td><td>[-0.97; 1.01]</td><td>0.973</td></tr><tr><td>C22:0</td><td>12</td><td>12</td><td>0.09 [0.07; 0.10]</td><td>1.48 [0.71; 2.00]</td><td>1.31</td><td>[0.79; 1.83]</td><td>0.003**</td></tr><tr><td>C22:1n-9</td><td>12</td><td>12</td><td>0.05 [0.04; 0.12]</td><td>0.39 [0.35; 0.68]</td><td>0.45</td><td>[0.24; 0.66]</td><td>0.006**</td></tr><tr><td>C22:4 n-6</td><td>12</td><td>12</td><td>0.12 [0.05; 0.16]</td><td>1.85 [0.72; 4.22]</td><td>2.12</td><td>[1.23; 3.01]</td><td>0.003**</td></tr><tr><td>C22:5 n-6</td><td>12</td><td>12</td><td>0.20 [0.16; 0.25]</td><td>1.03 [0.35; 2.64]</td><td>1.05</td><td>[0.18; 1.91]</td><td>0.054</td></tr><tr><td>C22:5 n-3 (DPA)</td><td>12</td><td>12</td><td>0.06 [0.04; 0.14]</td><td>2.02 [0.80; 3.99]</td><td>2.18</td><td>[1.39; 2.96]</td><td>0.002**</td></tr><tr><td>C22:6 n-3 (DHA)</td><td>12</td><td>12</td><td>0.18 [0.16; 0.21]</td><td>0.24 [0.15; 0.61]</td><td>0.15</td><td>[0.02; 0.27]</td><td>0.055</td></tr><tr><td>Total SFA</td><td>12</td><td>12</td><td>52.8 [51.3; 54.1]</td><td>50.2 [47.1; 52.5]</td><td>-2.11</td><td>[-4.89; 0.67]</td><td>0.184</td></tr><tr><td>Total MUFA</td><td>12</td><td>12</td><td>29.6 [28.1; 30.9]</td><td>27.8 [26.2; 31.2]</td><td>-0.22</td><td>[-3.42; 2.99]</td><td>0.898</td></tr><tr><td>Total n-6 PUFA</td><td>12</td><td>12</td><td>15.7 [15.0; 17.0]</td><td>16.3 [13.8; 17.7]</td><td>0.49</td><td>[-1.25; 2.23]</td><td>0.599</td></tr><tr><td>Total n-3 PUFA</td><td>12</td><td>12</td><td>1.8 [1.1; 2.8]</td><td>4.3 [1.8; 7.5]</td><td>2.53</td><td>[0.47; 4.58]</td><td>0.051</td></tr><tr><td>n-6/n-3 PUFA ratio</td><td>12</td><td>12</td><td>9.4 [5.6; 15.1]</td><td>3.8 [2.5; 7.6]</td><td>-7.68</td><td>[-10.84; -4.52]</td><td>0.005**</td></tr><tr><td>PUFA/SFA ratio</td><td>12</td><td>12</td><td>0.33 [0.31; 0.36]</td><td>0.41 [0.30; 0.56]</td><td>0.09</td><td>[0.004; 0.17]</td><td>0.084</td></tr></tbody></table>

The IL23A gene is expressed in placental trophoblastic cells and its over-expression was reported in the context of unexplained recurrent spontaneous abortion during early pregnancy⁴⁶ or in miscarriage⁴⁷. IL-23, which is produced, among others, by macrophages and dendritic cells, is an important component of the inflammatory response⁴⁸, by promoting the upregulation of the matrix metalloproteases MMP-9, increasing angiogenesis and being involved in the Th17 cell differentiation⁴⁸, and activating pro-inflammatory IL-17-producing cells in both mice and humans⁴⁹. Urban Particulate Matters were shown to promote an inflammatory milieu in part by inducing an IL-23-driven pro-inflammatory Th17 response in co-cultured cells *in vitro*⁵⁰.

TAGAP has been described to exert a regulatory role in T cell activation in Crohn's disease⁵¹.

FBXO32 encodes atrogin-1 (mostly expressed in muscular tissues), a member of the F-box protein family, one of the four subunits of the ubiquitin protein ligase complex, whose function is central for the selection of proteins targeted for degradation by the ubiquitin proteasome pathway⁵². In pregnant mice, food deprivation reduced FBXO32 gene expression in the placenta⁵², whereas its expression was increased in high-fat diet-fed turtles⁵³ or mice⁵⁴. These genes were upregulated in EC compared to CC placentas, but the physiological meaning of their increase remains to be determined.

A GSEA analysis was performed to identify biological pathways enriched in genes associated with grandmother exposure. This approach showed that 58 gene sets were enriched in EC placentas and 69 in CC placentas. In EC placentas, “proteasome complex”, “ubiquitin protein ligase activity” and “ubiquitin mediated proteolysis” gene sets were enriched. The proteasome is a protein-destroying apparatus involved in many essential cellular functions, such as the regulation of the cell cycle, cell differentiation, signal transduction pathways and apoptosis. The proteasome is involved in degrading a variety of cellular proteins in a rapid and timely way and most substrate proteins are modified by ubiquitin before their degradation by the proteasome⁵⁵. The slight over-expression of genes encoding the protein degradation machinery in EC placentas,
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Expression of FBXO32, IL23A, and TAGAP in placentas of F2 of control (CC) and exposed (EC) conceptuses at 28dpc, whose mothers were or not *in utero* exposed to diesel engine exhausts, respectively. Expressions were analyzed by microarray (left column) and RT-qPCR (right column). Microarray data were first log2-transformed and then normalized by an intra-array median subtraction. Relative expression of genes obtained by RT-qPCR was normalized by the expression of two housekeeping genes, RPL18 and EIF4E2. Box-plots are representative of the median value and quartiles. Significant differences between conditions were figured out by * (for P value < 0.05).</figcaption></figure>

together with the over-expression of the *FBXO32* gene involved in the formation of the proteasome, suggests an increase in protein degradation following either an increase in protein production or an increase in malformed and misfolded proteins. The over-expression of the proteasome complex component has been associated with premature birth⁵⁶ or preeclampsia⁵⁷ in humans. In intrauterine growth retardation, the ubiquitination of the amino acid transporter SNAT2 was reported to be increased⁵⁸, leading to decreased SNAT2 protein expression and decreased placental amino acid transport. Here, growth retardation, however, was not observed. The « RNA degradation » gene set was also enriched in EC placentas and comprised genes encoding proteins involved in the carbon catabolite-repression 4-Not complex (CCR4-NOT) that regulates mRNA synthesis and degradation as well as exosome trafficking⁵⁹. The over-expression of genes involved in the RNA degradation machinery was associated with the enrichment of “nuclear pore” and “nuclear envelope” gene sets (involved in RNA transport through the nuclear membrane) and with the enrichment of “organelle membrane”, “integral to organelle membrane” and “Golgi vesicle transport” gene sets (genes participating in protein transport between the cellular organelles). Altogether, these observations point to an increase in cellular protein and RNA trafficking. The “cation/ion channel activity” gene sets were downregulated in EC placentas. Finally, “systemic lupus erythematosus”, “allograft rejection”, “graft versus host disease” and “autoimmune thyroid disease” gene sets were downregulated in EC placentas. These gene sets include genes encoding for proteins involved in pro-inflammatory response and immune reactions with lymphocyte activation and TNF, IL10, and IL1A secretion. Although IL23A and TAGAP do not belong to these gene sets, these results are consistent with the general picture derived from FA analyses, suggesting that some protection against inflammatory processes is induced in F2 placentas from the exposed group.

In summary, *in utero*-exposed female rabbits showed an impaired metabolic phenotype in adulthood, and once pregnant, these animals exhibited liver steatosis, likely underlining a deregulation of lipid metabolism. Moreover foetuses from these females were dyslipidaemic and compared to controls, their placental lipid contents were higher, suggesting that these placentas cannot succeed in protecting the foetuses from maternal dyslipidaemia. However, the fatty acid concentrations and profiles of the placenta suggest adaptive protective mechanisms against inflammation through the placental storage of non-essential fatty acids. This mechanism was illustrated by MUFA storage and the favoured transplacental transfer of n-3 PUFA, which are known for their lipid-lowering effects and anti-inflammatory roles. Moreover, these placentas appeared programmed to counteract any inflammatory exposure as shown by a decrease in the pro-inflammatory arachidonic acid (AA) content and its consequent decrease in the foetal blood compartment as well. These results are consistent with the molecular changes observed in F2 placentas using GSEA analysis of transcriptomic placental data, showing
<figure><img src="image_5.png" /><figcaption>b Enrichment plot: Signal sequence binding</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>d</figcaption></figure>

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

**Figure 2.** Gene Sets Enrichment Analysis in F2 of control and exposed placental labyrinthine area at 28dpc. GSEA-enrichment plot of the most negatively enriched gene set: Olfactory transduction. (b) GSEA-enrichment plot of the most positively enriched gene set: Signal sequence binding. (c) Bar plot depicting the normalized enrichment score of the top 15 most negatively enriched gene sets, constructed using several categories of 2 databases: C2CP:KEGG (Orange), C5:BP (Green), C5:CC (Blue) and C5:MF (Red). (d) Bar plot depicting the normalized enrichment score of the top 15 most positively enriched gene sets, which contains several databases C2CP:KEGG (Orange), C5:BP (Green), C5:CC (Blue) and C5:MF (Red).

a downregulation of gene sets involved in inflammation, coherent with the upregulation of gene sets involved in protein metabolism (proteasome) and intracellular signalling. Contrary to what has been shown elsewhere, for example, in the case of maternal lipid over-nutrition in rabbits³⁵, suggesting that male and female placental responses differ, this was not the case here. Foetal sex has been taken into consideration in all of the statistical analyses, as an explicative factor, but no sex-related difference reached statistical significance in any of our datasets. This finding may be explained by the fact that compared with F1 controls, *in utero* exposure did not lead to observable differences in body weight or blood biochemistry; moreover, in adulthood, the F1 females were all fed with the same balanced and optimal gestation-specific diet. Thus, conditions might not have been met metabolically to induce or highlight different gender-specific responses or that the effects of pollution may have been much stronger than the sex-specific responses that would then have been either smoothed or erased. However, it cannot be excluded that the sample size was too small to highlight sex-specific effects. The authors acknowledge that this experimental study, which was continued until the F2 generation, has limitations due to small sample size for the control group, especially due to gestation failure. Thus although the authors have confidence in their results, these results must be considered with circumspection and therefore remain preliminary.
## Conclusion

In conclusion, the indirect exposure of F1 females *in utero* to DE might affect placental function and foetal metabolism in their F2 offspring at ended gestation. Different technical approaches, using membrane biochemistry and molecular biology, have demonstrated potentially favourable placental counter-regulation around inflammation to preserve F2 foetuses. Altogether, these preliminary results demonstrate that maternal exposure to air pollution due to diesel engine exhaust should be considered as a programming factor within the context of the Developmental Origins of Health and Diseases (DOHaD), with a high risk of intergenerational transmission. The mechanisms involved in this transmission and the long-term health effects of F2 after birth remain to be determined.

## Materials and Methods

**Animals.** All experimental protocols were approved by the local ethical committee COMETHEA (Comité d'Ethique en Expérimentation Animale du Centre INRA de Jouy en Josas et AgroParisTech), referenced as N°C2EA-45 in the French National registry CNREEA (Comité National de Reflexion Ethique sur l'Expérimentation Animale), under N°12/102. The methods were carried out in accordance with the relevant guidelines and regulations. As described previously, F0 pregnant New Zealand white female rabbits (INRA1077 line) were exposed by nose-only inhalation to diluted DE (exposed group) or clean air (control group) at a concentration of 1 mg/m³ with a mean particle diameter of 69 nm for 2 h/day and 5 days/week, from the 3rd to the 27th day post-conception (dpc)⁷. Briefly, DE exposure was performed with the Mobile Ambient Particle Concentrator Exposure Laboratory⁶⁰ connected to a 25KVA Loxam engine with a 500 nm particle filter. The measured components of the exposure mixture used in the present experiment were published previously⁷. Sixteen F0 females (N = 7 Controls, N = 9 Exposed dams) gave birth to F1 offspring, which were raised under control conditions. Altogether, 72 F1 offspring survived until adulthood, including 18 control females, 40 exposed females, 11 control males and 16 exposed males. All animals were weaned at 5 weeks of age. After becoming sexually mature adults at 6.5 months of age, a limited number of F1 nulliparous females (11 F1 control and 11 F1 exposed females) were dedicated to the production of the F2 generation. These females were mated with unexposed males. In the controls, only 5 out 11 rabbits were diagnosed as pregnant, among which one rabbit died at mid-pregnancy and another aborted; in the exposed group 10 out 11 rabbits were pregnant, including 1 abortion at mid-pregnancy. All females that were still pregnant (N = 3 controls and N = 9 exposed) were euthanized at 28 dpc to collect second generation (F2) foeto-placental units (see experimental protocol on Supplementary Figure S8). Maternal and foetal plasma were also collected and frozen at -80°C until further assay. As a reminder, in this study hereafter, CC refers to F2 foetuses collected from F1 control dams (N = 21) and EC refers to F2 foetuses collected from *in utero*-exposed F1 dams (N = 87). Foetuses were weighed; foetal crown-rump length, biparietal diameter and head length (using a digital caliper) and abdominal diameter were measured prior to dissection. Blood samples, sampled prior to sexing, were collected from all foetuses. During dissection, each foetus was sexed by visual observation of the internal genital organs and organs were subsequently collected and weighed. Two males and two females were retained per litter for the analyses; we have sexed the foetuses beginning from the left horn on the ovarian side towards the right horn, and we collected the samples of the first two males and two females among the dissected specimens. The labyrinthine area of the placenta (materno-foetal exchange zone, of foetal origin) was dissected from the placental junctional area and the decidua, and split for either fixing or snap freezing in liquid nitrogen for further analyses. Analyses were all made on samples from the same male-female couples per litter (2 couples/ CC litter and 1 couple/ EC litter), as follows: maternal plasma and hepatic samples (3 C and 9E), placental and foetal plasma samples (12 CC and 12 EC, randomly chosen among the collected ones, half females and half males) for biochemistry and fatty acid analyses. The histological and stereological evaluation of the labyrinthine area of the placenta were limited to 6 placental samples per group, randomly chosen from the same sampling assortment as those used just above. Similarly, 8 different placental samples were randomly chosen among the selected samples for transcriptomic analyses (eight replicates per condition).

**Biochemistry.** Classical biochemistry (triglycerides, cholesterol, glycaemia, urea, creatinine, ASAT, ALAT) was performed on maternal (n = 3 C and 9 E) and fetal EDTA-plasma (n = 12 CC and 12 EC) using a classical automated biochemistry analyser (DXc800 Beckman Coulter equipment). Plasma insulin was assayed with chemiluminescent immunoassay using paramagnetic particles on an automaton analyser (DXi Beckman Coulter equipment).

**Fatty acid analysis.** A total of 100 µL of maternal plasma (n = 3 C and 9 E), 100 µL of foetal plasma (n = 12 CC and 12 EC) and 300 mg of the labyrinthine area (n = 12 CC and 12 EC) from each selected conceptus and of maternal livers were used for lipid analysis. After classical lipid extraction of plasma⁶¹, fatty acids were transmethylated with BF3-methanol and the fatty acid methyl esters (FAME) were analysed by Gas Chromatography-FID (GC 3900 Varian, France) on an Econo-Cap EC-WAX capillary column³⁵. After lipid extraction of the placental tissue, the phospholipids (tissue membranes) were separated from neutral lipids (intracellular lipid store) using silica acid cartridges (Supelco, Bellefonte, PA) according to Juaneda *et al*.⁶² prior to transmethylation of the fatty acids from these two fractions and gas chromatography analysis⁶³. The identification of fatty acids was made in reference to known fatty acid profiles obtained from injection of standard FAME mix (Supelco 37 components FAME mix, ref 47885-U, Sigma). The fatty acid profile was established for each sample and expressed as a % of the total fatty acids. Heptadecanoic acid (C17:0), introduced prior to the plasma lipid extraction (50 µg) or prior to the fatty acid transmethylation of placental phospholipid and neutral lipid fractions (200 µg), was used as an internal standard to measure the total fatty acid concentration of plasma or for relative fatty acid quantification of
<table><caption>**Table 6.** Primers (DNA oligos) used for FBXO32, IL23A and TAGAP mRNA quantification by real-time RT qPCR, in *Oryctolagus cuniculus* species. *Endogenous transcripts*, *EIF4E2* and *RPL18*, were used as reference transcripts.</caption><thead><tr><th>Genes</th><th>Accession number</th><th>Sequences of forward (FW) and reverse (RV) primers</th></tr></thead><tbody><tr><td rowspan="2">FBXO32</td><td rowspan="2">XM_002710762</td><td>L-FBXO32-FW</td><td>CCT-CCA-GCT-CTG-CAA-GCA-TT</td></tr><tr><td>L-FBXO32-RV</td><td>AAA-GTC-CTG-GGG-TGA-AAG-TGA</td></tr><tr><td rowspan="2">IL23A</td><td rowspan="2">XM_002711079</td><td>L-IL23A-FW</td><td>TGG-GAC-CCA-TGG-AGC-TAC-TA</td></tr><tr><td>L-IL23A-RV</td><td>TGC-AGG-CAG-AAC-TGA-GTG-TT</td></tr><tr><td rowspan="2">TAGAP</td><td rowspan="2">XM_008263738</td><td>L-TAPAG-FW</td><td>GGG-AAG-GTT-AAG-ACG-CTG-GT</td></tr><tr><td>L-TAPAG-RV</td><td>GGT-GGA-CGT-GTC-TGA-ACT-GT</td></tr><tr><td rowspan="2">EIF4E2</td><td rowspan="2">ENSOCUG00000004534</td><td>L-EIF4E2-FW</td><td>TGG-CAA-GTG-GAT-TAT-TCG-GC</td></tr><tr><td>L-EIF4E2-RV</td><td>CAG-AGA-CCA-CAG-CCC-CAC-AG</td></tr><tr><td rowspan="2">RPL18</td><td rowspan="2">DQ403030</td><td>L-RPL18-FW</td><td>CAA-CTC-CAC-GTT-CAA-CCA-GGT</td></tr><tr><td>L-RPL18-RV</td><td>GGT-CTT-GTT-CTC-CCG-CCC</td></tr></tbody></table>

each placental lipid fraction after chromatogram analysis. The fatty acids were identified using chromatographic analyses as described in the Supplementary Note. The DNL index was determined in maternal plasma as de novo lipogenesis (C16:0/C18:2n-6 ratio)²²,²³. Only a few relevant fatty acids are outlined in the result tables, but all the fatty acids are included in reported totals. Moreover, Principal Component Analyses (PCA) were performed using the full fatty acid profiles (cf. Statistics).

**Liver histology.** Maternal liver fragments from all F1 dams were fixed in 4% formaldehyde for at least 10 days prior to dehydration in ethanol solutions, clearing in butanol, embedding in paraffin and cutting into 5 µm thick sections (Leica microtome, Germany). Four different colourations were performed on liver sections: a routine hematoxyline eosine Safran (HES) colouration to reveal the different cellular components, a Prussian blue coloration (Perl's colouration) to highlight possible ferric deposits and the activation of Kupffer cells, a Masson's trichrome colouration (Microm Microtech France) to visualize collagen fibres and a periodic acid Schiff (PAS) staining to highlight the presence of glycogen. All coloured sections were analysed under light microscopy and scanned using the Nanozoomer Digital Pathology System (NDP Scan U10074-01, Hamamatsu, Japan).

**Immunohistochemistry.** Samples of the labyrinthine area (N = 6 CC and N = 6 EC) were fixed with 4% formalin. After dehydration in ethanol solutions and washing in xylene, and the samples were embedded in paraffin and cut into 7 µm thick sections (Leica microtome, Germany). Immunodetection of vimentin was performed to label foetal capillaries as described previously³⁵.

**Stereological analyses.** The cellular composition of the labyrinthine area was explored by stereology. After vimentin immunodetection, all placental sections were scanned using the NanoZoomer Digital Pathology System (NDP Scan U10074-01) as previously described⁶⁴. The volume fraction and surface density of trophoblasts, foetal vessels and maternal blood space, i.e., the components of the labyrinthine area, were quantified using the One Stop Stereology method available on Mercator® software, as described previously⁶⁵.

**Transcriptomic analyses.** RNA isolation. Total RNA was extracted from frozen pieces of labyrinthine area in liquid nitrogen, randomly collected from each control (n = 8 CC), exposed (n = 8 EC) litter, half males and half females. The samples were ground in a solution D + (namely: thiocyanate guanidine 4 M, sodium citrate 25 mM pH 7, N-lauryl-sarcosine 0.5%, β-mercaptoethanol 0.1 M). The RNA was extracted twice with 3 M sodium acetate pH 4–5, aquaphenol and chloroform/isoamlic ethanol (24:1), then precipitated in isopropanol at –80 °C overnight. A purification procedure using DNAse I (Qiagen) treatment at 25 °C for 15 min was performed prior to elution. The RNA pellet was washed with ethanol, dried and solubilized in RNase-free water (10 µl of H₂O for 100 mg of tissue). The RNA was purified with the RNeasy Mini Kit (Qiagen) and quantified by Nanodrop (Thermo Scientific, Madrid, Spain), and the quality was verified on the Bioanalyzer 2100 (Agilent Technologies). Total extracted RNA was stored at –80 °C for further RNA labelling.

**RNA amplification and labelling and microarray processing.** The eight samples from each group were used as biological replicates for the microarray analysis. RNA from each placental tissue (200 ng) was amplified and labelled with cyanine 3 (Cy3) dye using the One-colour Microarray-based Gene Expression Analysis Low Input Quick Amp Labelling Kit (Agilent Technologies). Excess dye was removed with the RNeasy Mini Kit (Qiagen) and dye incorporation and concentration were determined using the microarray setting on the Nanodrop 2000 and the Bioanalyzer 2100. Hybridization was performed at the CRB GADIE (INRA Jouy-en-Josas, France, http://crb-gadie.inra.fr/). Equal amounts of Cy3-labelled samples (600 ng, Specific Activity > 6.0 pmol Cy3/µg of cRNA) were mixed with blocking agent and fragmentation buffer, and then 24 µl of the mixture was hybridized into the Rabbit Custom Gene Expression Microarrays (GEO accession GPL18913, Agilent-042421 Rabbit BDR version 2, Agilent Technologies) for transcriptomic analysis (12775 genes)⁶⁶. After 17 h at 65 °C, the hybridized slides were washed twice for 1 min at room temperature in GE Wash Buffer 1 and 1 min at 37 °C with GE Wash Buffer 2 (Agilent Technologies) and air-dried and scanned immediately.
**Microarray analysis.** After hybridization, the scanned images were analysed using Feature Extraction software (v10.7.3.1; Agilent Technologies). The data were normalized using intra-array median subtraction and log2 transformation. The raw data intensity files were read into R (www.r-project.org). The identification of differentially expressed genes from EC and CC F2 placentas was achieved using the Limma package⁶⁷. P values obtained by this analysis were adjusted for multiple testing using the Benjamini-Hochberg procedure. Probes with an adjusted P value < 0.05 were considered significant.

**Gene Set Enrichment Analysis (GSEA).** The standard strategy to analyse transcriptomic data is based on a single-gene level that is very efficient when perturbations exert strong impacts on individual genes. This classical approach, however, becomes inefficient with subtle conditions because only a minor fraction of the information contained in whole-genome transcriptomic dataset is used, due to arbitrary thresholds used²⁴. Gene set enrichment analysis (GSEA) was established for use in metabolic studies in the 2000s⁶⁸. All detected unfiltered transcripts are mapped to their intended pathways (or other “sets of interest”). The subsequent analysis within specific pathways enables the study of the regulation of specific pathways that include moderately regulated (and non-regulated) genes usually lost due to experimental background noise. Here, the GSEA approach⁶⁸ was used on transcriptomic data to systematically identify biological pathways enriched in genes associated with F1 foetal exposure and performed using the following parameters: 1,000 gene set permutations, weighed enrichment statistics, gene set size between 15 and 500 and signal-to-noise metrics. Regulated pathways were considered statistically significant when the false discovery rate (FDR) q-value was <0.25. The GSEA-derived normalized enrichment score was used for the visualization of pathway regulation, using the KEGG (The Kyoto Encyclopedia of Genes and Genomes) and GO (Gene Ontology) databases. KEGG is a database resource for understanding high-level functions and utilities of the biological system, from molecular-level information, especially large-scale molecular datasets generated by high-throughput experimental technologies. One of the main uses of GO is to perform enrichment analysis on gene sets. For example, given a set of genes that are up-regulated under certain conditions, an enrichment analysis will find which GO terms are over-represented (or under-represented) using annotations for that gene set. Thus C2:KEGG and several categories of the more detailed GO gene set database, i.e. C5:BP (Biological Processes), C5:CC (Cellular Component) and C5:MF (Molecular Function)⁶⁹, were used to identify the top 15 pathways characterizing each group.

**Reverse transcription and real-time PCR.** For validation of the microarray experiments, total RNA was extracted from pieces of the labyrinthine area from the same control (n=8 CC) and exposed (n=8 EC) placentas. The RNA was isolated on silica-based columns, using the RNeasy Plus Mini Kit with the effective elimination of genomic DNA, following the manufacturer's instructions (Qiagen, France). Total RNA was quantified using OD 260 nm on a NanoDrop 1000 spectrophotometer (ThermoScientific, Baltimore, MA, US) and quality was assessed with the Fragment Analyzer (AATI, US) to determine the RNA quality number. RNA (750 ng) was reverse-transcribed using the iScript Reverse Transcription Supermix (Bio-Rad, Hercules, CA, US). The expression of conventionally used housekeeping genes, including RPL18 and EIF4E2, was shown to be stable between assays. Quantitative real-time PCR (qRT-PCR) was performed with the Fast SYBR Green Master Mix and the Applied Biosystems 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster city, California US). The experimental protocol consisted of an initial polymerase activation at 95 °C, for 20 s followed by an amplification program for 40 cycles maintaining the annealing and extension primer temperature at 62 °C for 20 s. A melting-curve analysis was then performed to verify the amplification of a single product. All primers were designed using NCBI Primer-BLAST and selected to generate amplicons with a length of 100–200 bp (Table 6). The primers were synthesized by Eurogentec. Standard curves were generated to calculate the efficiency of each set. Only primer sets with an efficiency of 1.85 to 2.1 were used for qPCR. The relative mRNA levels for each assay were computed from the Ct values obtained for the target gene. The data were pooled and analysed using RQ manager v1.2 and dataAssist v1.0 (Applied Biosystems). DataAssist™ software is a simple, yet powerful data analysis tool for sample comparison when using the comparative cycle threshold (CT) method⁷⁰ for calculating the relative quantitation of gene expression. This software contains a filtering procedure for outlier removal, various normalization methods based on single or multiple genes, and provides relative quantification analysis of gene expression through a combination of statistical analysis and interactive visualization.

**Statistical analyses.** All data (except placental transcriptomic data) are expressed as median[Q1;Q3], with the first (Q1) and third quartiles (Q3) corresponding to 25 and 75% of the scores, respectively. For all F1 material data, a linear model was used, with random effect of dams adjusted for treatment, birth date (indexed in 3 categories), integrating neonatal and peri-weaning data of these F1 females to account for their early nutritional status (litter size at birth and at weaning) using the linear mixed effects model (nlme package, R, Pinheiro, Bates, DebRoy, Sarkar and the R Development Core Team 2013. nlme. R package version 3.1-111; http://www.r-project.org/). The β-value represents sample variability after adjustment. A negative β-value represents a reduction in the exposed group and a positive β-value represents an increase in the exposed group compared to controls. A similar linear model was used for F2 foetal data with a random effect of dams adjusted for treatment, litter size, foetal position in the horn (indexed in 3 categories - A, B or C – to indicate the position of the foetus in the horn, with A: at the top of the horn near the ovary, B: in the middle of the horn or C: at the bottom of the horn towards the vagina) and foetal sex as an explicative factor. For fatty acid profiles, linear mixed models identical to those used for maternal and foetal data were used to analyse each fatty acid for each profile. Principal Component Analyses (PCA) were performed using FactoMineR (R package⁷¹ dedicated to multivariate data analysis⁶⁹) to analyse the biological relevance of the significant differences observed using the linear model.

**Ethical approval.** All experiments were carried out in accordance with relevant guidelines and regulations. All experimental protocols were approved by the local ethical committee, referenced as N°45 in the French National register, under N°12/102.
## Data Availability

The datasets supporting the conclusions of this article are included within the article and its Supplemental Files; microarray data were registered in Gene Expression Omnibus (GEO) under accession number GSE109831.

## References

1. Landrigan, P. J. et al. The Lancet Commission on pollution and health. *Lancet* **391**, 462–512, https://doi.org/10.1016/S0140-6736(17)32345-0 (2018).

2. Steiner, S., Bisig, C., Petri-Fink, A. & Rothen-Rutishauser, B. Diesel exhaust: current knowledge of adverse effects and underlying cellular mechanisms. *Archives of toxicology* **90**, 1541–1553, https://doi.org/10.1007/s00204-016-1736-5 (2016).

3. Liati, A., Schreiber, D., Arroyo Rojas Dasilva, Y. & Dimopoulos Eggenschwiler, P. Ultrafine particle emissions from modern Gasoline and Diesel vehicles: An electron microscopic perspective. *Environmental pollution* **239**, 661–669, https://doi.org/10.1016/j.envpol.2018.04.081 (2018).

4. Pereira, L. A. et al. Association between air pollution and intrauterine mortality in Sao Paulo, Brazil. *Environmental health perspectives* **106**, 325–329 (1998).

5. Lakshmi, P. V. et al. Household air pollution and stillbirths in India: analysis of the DLHS-II National Survey. *Environmental research* **121**, 17–22, https://doi.org/10.1016/j.envres.2012.12.004 (2013).

6. Ritz, S. A. Air pollution as a potential contributor to the 'epidemic' of autoimmune disease. *Medical hypotheses* **74**, 110–117, https://doi.org/10.1016/j.mehy.2009.07.033 (2010).

7. Valentino, S. A. et al. Maternal exposure to diluted diesel engine exhaust alters placental function and induces intergenerational effects in rabbits. *Particle and fibre toxicology* **13**, 39, https://doi.org/10.1186/s12989-016-0151-7 (2016).

8. Veras, M. M. et al. Particulate urban air pollution affects the functional morphology of mouse placenta. *Biology of reproduction* **79**, 578–584, https://doi.org/10.1095/biolreprod.108.069591 (2008).

9. Veras, M. M. et al. Chronic exposure to fine particulate matter emitted by traffic affects reproductive and fetal outcomes in mice. *Environmental research* **109**, 536–543, https://doi.org/10.1016/j.envres.2009.03.006 (2009).

10. Weldy, C. S., Liu, Y., Liggitt, H. D. & Chin, M. T. In utero exposure to diesel exhaust air pollution promotes adverse intrauterine conditions, resulting in weight gain, altered blood pressure, and increased susceptibility to heart failure in adult mice. *PloS one* **9**, e88582, https://doi.org/10.1371/journal.pone.0088582 (2014).

11. Rahmalia, A. et al. Pregnancy exposure to atmospheric pollutants and placental weight: an approach relying on a dispersion model. *Environment international* **48**, 47–55, https://doi.org/10.1016/j.envint.2012.06.013 (2012).

12. van den Hooven, E. H. et al. Air pollution exposure and markers of placental growth and function: the generation R study. *Environmental health perspectives* **120**, 1753–1759, https://doi.org/10.1289/ehp.1204918 (2012).

13. Yorifuji, T. et al. Residential proximity to major roads and placenta/birth weight ratio. *The Science of the total environment* **414**, 98–102, https://doi.org/10.1016/j.scitotenv.2011.11.001 (2012).

14. Godfrey, K. M., Gluckman, P. D. & Hanson, M. A. Developmental origins of metabolic disease: life course and intergenerational perspectives. *Trends in endocrinology and metabolism: TEM* **21**, 199–205, https://doi.org/10.1016/j.tem.2009.12.008 (2010).

15. Ambeskovic, M., Roseboom, T. J. & Metz, G. A. S. Transgenerational effects of early environmental insults on aging and disease incidence. *Neuroscience and biobehavioral reviews*, https://doi.org/10.1016/j.neubiorev.2017.08.002 (2017).

16. Li, C. et al. Effects of in utero exposure to nanoparticle-rich diesel exhaust on testicular function in immature male rats. *Toxicology letters* **185**, 1–8, https://doi.org/10.1016/j.toxlet.2008.11.012 (2009).

17. Ogliari, K. S. et al. Intrauterine exposure to diesel exhaust diminishes adult ovarian reserve. *Fertility and sterility* **99**, 1681–1688, https://doi.org/10.1016/j.fertnstert.2013.01.103 (2013).

18. Hanafi, M. Y. et al. Diabetes-induced perturbations are subject to intergenerational transmission through maternal line. *Journal of physiology and biochemistry* **72**, 315–326, https://doi.org/10.1007/s13105-016-0483-7 (2016).

19. Skinner, M. K. Endocrine disruptors in 2015: Epigenetic transgenerational inheritance. *Nature reviews. Endocrinology* **12**, 68–70, https://doi.org/10.1038/nrendo.2015.206 (2016).

20. Fischer, B., Chavatte-Palmer, P., Viebahn, C., Navarrete Santos, A. & Duranthon, V. Rabbit as a reproductive model for human health. *Reproduction* **144**, 1–10, https://doi.org/10.1530/REP-12-0091 (2012).

21. Furukawa, S., Kuroda, Y. & Sugiyama, A. A comparison of the histological structure of the placenta in experimental animals. *Journal of toxicologic pathology* **27**, 11–18, https://doi.org/10.1293/tox.2013-0060 (2014).

22. Lee, J. J. et al. Palmitoleic acid is elevated in fatty liver disease and reflects hepatic lipogenesis. *The American journal of clinical nutrition* **101**, 34–43, https://doi.org/10.3945/ajcn.114.092262 (2015).

23. Zong, G. et al. Associations of erythrocyte fatty acids in the de novo lipogenesis pathway with risk of metabolic syndrome in a cohort study of middle-aged and older Chinese. *The American journal of clinical nutrition* **98**, 319–326, https://doi.org/10.3945/ajcn.113.061218 (2013).

24. Weidner, C., Steinfath, M., Opitz, E., Oelgeschlager, M. & Schonfelder, G. Defining the optimal animal model for translational research using gene set enrichment analysis. *EMBO molecular medicine* **8**, 831–838, https://doi.org/10.15252/emmm.201506025 (2016).

25. Bolton, J. L. et al. Prenatal air pollution exposure induces neuroinflammation and predisposes offspring to weight gain in adulthood in a sex-specific manner. *FASEB J* **26**, 4743–4754, https://doi.org/10.1096/fj.12-210989 (2012).

26. Sivan, E. & Boden, G. Free fatty acids, insulin resistance, and pregnancy. *Current diabetes reports* **3**, 319–322 (2003).

27. Herrera, E. & Desoye, G. Maternal and fetal lipid metabolism under normal and gestational diabetic conditions. *Hormone molecular biology and clinical investigation* **26**, 109–127, https://doi.org/10.1515/hmbci-2015-0025 (2016).

28. Cerf, M. E. & Herrera, E. High Fat Diet Administration during Specific Periods of Pregnancy Alters Maternal Fatty Acid Profiles in the Near-Term Rat. *Nutrients* **8**, https://doi.org/10.3390/nu8010025 (2016).

29. Gilbert, M., Basile, S., Baudelin, A. & Pere, M. C. Lowering plasma free fatty acid levels improves insulin action in conscious pregnant rabbits. *The American journal of physiology* **264**, E576–E582, https://doi.org/10.1152/ajpendo.1993.264.4.E576 (1993).

30. Menchetti, L. et al. Food restriction during pregnancy in rabbits: effects on hormones and metabolites involved in energy homeostasis and metabolic programming. *Research in veterinary science* **98**, 7–12, https://doi.org/10.1016/j.rvsc.2014.11.017 (2015).

31. Ortiz, L., Nakamura, B., Li, X., Blumberg, B. & Luderer, U. In utero exposure to benzo[a]pyrene increases adiposity and causes hepatic steatosis in female mice, and glutathione deficiency is protective. *Toxicology letters* **223**, 260–267, https://doi.org/10.1016/j.toxlet.2013.09.017 (2013).

32. Ortiz, L., Nakamura, B., Li, X., Blumberg, B. & Luderer, U. Reprint of “In utero exposure to benzo[a]pyrene increases adiposity and causes hepatic steatosis in female mice, and glutathione deficiency is protective”. *Toxicology letters* **230**, 314–321, https://doi.org/10.1016/j.toxlet.2013.11.017 (2014).

33. Burdge, G. C. DHA supplementation during pregnancy and DNA methylation in cord blood leukocytes. *The American journal of clinical nutrition* **98**, 1594–1595, https://doi.org/10.3945/ajcn.113.072074 (2013).

34. Picone, O. et al. Hyperlipidic hypercholesterolemic diet in prepubertal rabbits affects gene expression in the embryo, restricts fetal growth and increases offspring susceptibility to obesity. *Theriogenology* **75**, 287–299, https://doi.org/10.1016/j.theriogenology.2010.08.015 (2011).

35. Tarrade, A. et al. Sexual dimorphism of the feto-placental phenotype in response to a high fat and control maternal diets in a rabbit model. *PloS one* **8**, e83458, https://doi.org/10.1371/journal.pone.0083458 (2013).
36. Keddad, K., Therond, P., Motta, C., Baussan, C. & Legrand, A. Alterations in erythrocyte membrane fluidity and fatty acid composition in glycogen storage disease. *Biochimica et biophysica acta* **1315**, 61–65 (1996).

37. Cetin, I., Parisi, F., Berti, C., Mando, C. & Desoye, G. Placental fatty acid transport in maternal obesity. *Journal of developmental origins of health and disease* **3**, 409–414, https://doi.org/10.1017/S2040174412000414 (2012).

38. Perazzolo, S. et al. The influence of placental metabolism on fatty acid transfer to the fetus. *Journal of lipid research* **58**, 443–454, https://doi.org/10.1194/jlr.P072355 (2017).

39. Segura, M. T. et al. Maternal BMI and gestational diabetes alter placental lipid transporters and fatty acid composition. *Placenta* **57**, 144–151, https://doi.org/10.1016/j.placenta.2017.07.001 (2017).

40. Johnsen, G. M., Weedon-Fekjaer, M. S., Tobin, K. A., Staff, A. C. & Duttaroy, A. K. Long-chain polyunsaturated fatty acids stimulate cellular fatty acid uptake in human placental choriocarcinoma (BeWo) cells. *Placenta* **30**, 1037–1044, https://doi.org/10.1016/j.placenta.2009.10.004 (2009).

41. Sergeant, S., Rahbar, E. & Chilton, F. H. Gamma-linolenic acid, Dihommo-gamma linolenic, Eicosanoids and Inflammatory Processes. *European journal of pharmacology* **785**, 77–86, https://doi.org/10.1016/j.ejphar.2016.04.020 (2016).

42. Herrera, E. Lipid metabolism in pregnancy and its consequences in the fetus and newborn. *Endocrine* **19**, 43–55, https://doi.org/10.1385/ENDO:19:1:43 (2002).

43. Leghi, G. E. M. B.S. The effect of n-3 LCPUFA supplementation on oxidative stress and inflammation in the placenta and maternal plasma during pregnancy. *Prostaglandins, leukotrienes, and essential fatty acids* **113**, 33–39, https://doi.org/10.1016/j.plefa.2016.08.010 (2016).

44. Kalogeropoulos, N. et al. Unsaturated fatty acids are inversely associated and n-6/n-3 ratios are positively related to inflammation and coagulation markers in plasma of apparently healthy adults. *Clinica chimica acta; international journal of clinical chemistry* **411**, 584–591, https://doi.org/10.1016/j.cca.2010.01.023 (2010).

45. Chang, C. L. & Deckelbaum, R. J. Omega-3 fatty acids: mechanisms underlying ‘protective effects’ in atherosclerosis. *Current opinion in lipidology* **24**, 345–350, https://doi.org/10.1097/MOL.0b013e3283616364 (2013).

46. Cai, J. Y. & Li, M. J. Interleukin 23 regulates the functions of human decidual immune cells during early pregnancy. *Biochemical and biophysical research communications* **469**, 340–344, https://doi.org/10.1016/j.bCRC.2015.11.118 (2016).

47. Saifi, B. et al. Th17 cells and related cytokines in unexplained recurrent spontaneous miscarriage at the implantation window. *Reproductive biomedicine online* **29**, 481–489, https://doi.org/10.1016/j.rbmo.2014.06.008 (2014).

48. Darmochwal-Kolarz, D. et al. The Role of Interleukin-17, Interleukin-23, and Transforming Growth Factor-beta in Pregnancy Complicated by Placental Insufficiency. *BioMed research international* **2017**, 6904325, https://doi.org/10.1155/2017/6904325 (2017).

49. Goriely, S. & Goldman, M. Interleukin-12 family members and the balance between rejection and tolerance. *Current opinion in organ transplantation* **13**, 4–9, https://doi.org/10.1097/MOT.0b013e3282f406c4 (2008).

50. Mann, E. H. et al. Vitamin D Counteracts an IL-23-Dependent IL-17A(+)IFN-gamma(+) Response Driven by Urban Particulate Matter. *American journal of respiratory cell and molecular biology* **57**, 355–366, https://doi.org/10.1165/rcmb.2016-0409OC (2017).

51. Connelly, T. M. et al. T-cell activation Rho GTPase-activating protein expression varies with inflammation location and severity in Crohn's disease. *The Journal of surgical research* **190**, 457–464, https://doi.org/10.1016/j.jss.2014.01.019 (2014).

52. Chu, A. et al. Gestational food restriction decreases placental interleukin-10 expression and markers of autophagy and endoplasmic reticulum stress in murine intrauterine growth restriction. *Nutrition research* **36**, 1055–1067, https://doi.org/10.1016/j.nutres.2016.08.001 (2016).

53. Liu, L. et al. The effect of high fat diet on daily rhythm of the core clock genes and muscle functional genes in the skeletal muscle of Chinese soft-shelled turtle (Trionyx sinensis). *Comparative biochemistry and physiology. Part B, Biochemistry & molecular biology* **213**, 17–27, https://doi.org/10.1016/j.cbpb.2017.07.002 (2017).

54. Le, N. H. et al. Quercetin protects against obesity-induced skeletal muscle inflammation and atrophy. *Mediators of inflammation* **2014**, 834294, https://doi.org/10.1155/2014/834294 (2014).

55. Livneh, I., Cohen-Kaplan, V., Cohen-Rosenzweig, C., Avni, N. & Ciechanover, A. The life cycle of the 26S proteasome: from birth, through regulation and function, and onto its death. *Cell research* **26**, 869–885, https://doi.org/10.1038/cr.2016.86 (2016).

56. Luo, X. et al. Epigenetic regulation of lncRNA connects ubiquitin-proteasome system with infection-inflammation in preterm births and preterm premature rupture of membranes. *BMC pregnancy and childbirth* **15**, 35, https://doi.org/10.1186/s12884-015-0460-0 (2015).

57. Kang, J. H. et al. Preeclampsia leads to dysregulation of various signaling pathways in placenta. *Journal of hypertension* **29**, 928–936, https://doi.org/10.1097/HJH.0b013e328344a82c (2011).

58. Chen, Y. Y. et al. Increased ubiquitination and reduced plasma membrane trafficking of placental amino acid transporter SNAT-2 in human IUGR. *Clinical science* **129**, 1131–1141, https://doi.org/10.1042/CS20150511 (2015).

59. Miller, J. E. & Reese, J. C. Ccr4-Not complex: the control freak of eukaryotic cells. *Critical reviews in biochemistry and molecular biology* **47**, 315–333, https://doi.org/10.3109/10409238.2012.667214 (2012).

60. Freney, E. J. et al. A single-particle characterization of a mobile Versatile Aerosol Concentration Enrichment System for exposure studies. *Particle and fibre toxicology* **3**, 8, https://doi.org/10.1186/1743-8977-3-8 (2006).

61. Folch, J., Lees, M. & Sloane Stanley, G. H. A simple method for the isolation and purification of total lipides from animal tissues. *J Biol Chem* **226**, 497–509 (1957).

62. Juaneda, P. & Rocquelin, G. Rapid and convenient separation of phospholipids and non phosphorus lipids from rat heart using silica cartridges. *Lipids* **20**, 40–41 (1985).

63. Morrison, W. R. & Smith, L. M. Preparation of Fatty Acid Methyl Esters and Dimethylacetals from Lipids with Boron Fluoride–Methanol. *Journal of lipid research* **5**, 600–608 (1964).

64. Tarrade, A. et al. Analysis of placental vascularization in a pharmacological rabbit model of IUGR induced by L-NAME, a nitric oxide synthase inhibitor. *Placenta* **35**, 254–259, https://doi.org/10.1016/j.placenta.2014.01.013 (2014).

65. Favaron, P. O. et al. Morphometric analysis of the placenta in the New World mouse Necromys lasiurus (Rodentia, Cricetidae): a comparison of placental development in cricetids and murids. *Reproductive biology and endocrinology: RB&E* **11**, 10, https://doi.org/10.1186/1477-7827-11-10 (2013).

66. Jacquier, V. et al. Genome-wide immunity studies in the rabbit: transcriptome variations in peripheral blood mononuclear cells after *in vitro* stimulation by LPS or PMA-Ionomycin. *BMC genomics* **16**, 26, https://doi.org/10.1186/s12864-015-1218-9 (2015).

67. Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic acids research* **43**, e47, https://doi.org/10.1093/nar/gkv007 (2015).

68. Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. *Proceedings of the National Academy of Sciences of the United States of America* **102**, 15545–15550, https://doi.org/10.1073/pnas.0506580102 (2005).

69. de Tayrac, M., Le, S., Aubry, M., Mosser, J. & Husson, F. Simultaneous analysis of distinct Omics data sets with integration of biological knowledge: Multiple Factor Analysis approach. *BMC genomics* **10**, 32, https://doi.org/10.1186/1471-2164-10-32 (2009).

70. Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative C(T) method. *Nature protocols* **3**, 1101–1108, (2008).

71. Lè, S., Josse, J. & Husson, F. FactoMineR: An R Package for Multivariate Analysis. *Journal of Statistical Software*. **25**, 1–18 (2008).
## Acknowledgements

The authors thank the UCEA (Unité Commune d'Expérimentation Animale) for housing the rabbits and participating in the experiments. The authors thank John Boere and Paul Fokkens of the Netherlands National Institute for Public Health and the Environment for valuable technical assistance during the diesel engine exhaust exposure. The authors would also thank "ProfileXpert", service de <<Génomique & Microgénomique>>, Université Lyon 1, SFR santé LYON-EST, UCBL-INSERM US 7-CNRS UMS 3453, 8 av. Rockefeller 69373 Lyon Cedex 08, for validation by RT-qPCR. The authors wish to thank Marine Guinot (BDR) for technical support in the preparation of liver histological slides and APEX from ONIRIS, Nantes for advice and assistance with the liver histological analyses. Moreover the authors thank all the Master students involved in this project (Pauline Hoarau, Jérôme Delmotte, and Julie Dechartres) and those who participated in this project from near and far, from the BDR Unit or elsewhere. The authors are members of COST Actions FA0702 GEMINI "Maternal interaction with gametes and embryo", FA1201 EPICONCEPT "Epigenetics and periconceptual environment" and BM1308 "Sharing advances on large animal models (SALAAM)". This study was supported by ANR grant ANR-13-CESA-0011-EPAPP (P. Chavatte-Palmer) and by ERC consolidator grant N°311765-E-DOHaD (PI, R. Slama).

## Author Contributions

P.C.P., D.R.R., S.A.V., F.R.C. and A.C.T. conceived and designed the experiment and took part in writing this paper. P.C.P. and R.S. acquired the financial support for the project leading to this publication. S.A.V., D.R.R., C.R., M.D., M.S.L. contributed to the investigation process with the animals. D.R.R., S.A.V., A.C.T., J.A., C.R., M.C.A., M.D., M.S.L. and P.C.P. participated to data and tissue collections. D.R.R. and M.S.L. performed all the lipid analyses. C.A., A.V.C., J.L. and L.J. performed transcriptomic analyses under VD's supervision. S.A.V. performed GSEA analyses of transcriptomic data. A.C.T. supervised histological approaches on placenta. N.F. performed the blood sample analyses. S.A.V. and R.S. applied the statistical linear model to analyze the study data. D.R.R. and S.A.V. wrote the paper. A.C.T., D.R.R. and P.C.P. supervised the project. All authors discussed and approved the results. All authors read and approved the final manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-46130-x.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019