# Impact of exposure to diesel exhaust during pregnancy on mammary gland development and milk composition in the rabbit

Catherine Hue-Beauvais, Etienne Aujean, Guy Miranda, Delphine Rousseau-Ralliard, Sarah Valentino, Nicolas Brun, Stessy Ladebese, Christine Péchoux, Pascale Chavatte-Palmer, Madia Charlier

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

## ► To cite this version:

Catherine Hue-Beauvais, Etienne Aujean, Guy Miranda, Delphine Rousseau-Ralliard, Sarah Valentino, et al.. Impact of exposure to diesel exhaust during pregnancy on mammary gland development and milk composition in the rabbit. PLoS ONE, 2019, 14 (2), pp.1-18. 10.1371/journal.pone.0212132. hal-02050122

HAL Id: hal-02050122

https://hal.science/hal-02050122

Submitted on 26 Feb 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>
RESEARCH ARTICLE

# Impact of exposure to diesel exhaust during pregnancy on mammary gland development and milk composition in the rabbit

Cathy Hue-Beauvais¹, Etienne Aujean¹, Guy Miranda¹, Delphine Ralliard-Rousseau², Sarah Valentino², Nicolas Brun¹, Stessy Ladebese¹, Christine Péchoux¹, Pascale Chavatte-Palmer², Madia Charlier¹*

1 UMR GABI, INRA, AgroParisTech, Université Paris-Saclay, Jouy-en-Josas, France, 2 UMR BDR, INRA, ENVA, Université Paris Saclay, Jouy-en-Josas, France

* madia.charlier@inra.fr

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

## Abstract

Exposure to fine-particulate air pollution is a major global health concern because it is associated with reduced birth weight and an increased risk of cardiovascular disease. Here we have investigated the potential for exposure to diesel exhaust during pregnancy to influence mammary gland development and milk composition. Female rabbits were therefore exposed by nose-only inhalation to either diluted diesel exhaust fumes (1 mg/m³) or clean air for 2h/day, 5 days/week, from the 3rd to the 27th days of pregnancy. On Day 28 of pregnancy, mammary glands were collected from twelve females (six controls and six diesel-exposed) and assessed for morphological and functional alterations. Milk samples were collected from eighteen dams (nine controls and nine diesel-exposed) during early (days 2 to 4) and established (days 13 to 16) lactation to verify the composition of fatty acids and major proteins and leptin levels. The mammary alveolar lumina contained numerous fat globules, and stearoyl CoA reductase expression was higher in mammary epithelia from diesel exhaust-exposed rabbits, which together suggested increased mammary lipid biosynthesis. Gas chromatography analysis of the composition of milk fatty acids revealed a sharp rise in the total fatty acid content, mainly due to monounsaturated fatty acids. Liquid chromatography-mass spectrometry analysis of milk samples enabled identification and quantification of the main rabbit milk proteins and their main phosphorylated isoforms, and revealed important changes to individual casein and whey protein contents and to their most phosphorylated isoforms during early lactation. Taken together, these findings suggest that repeated daily exposure to diesel exhaust fumes during pregnancy at urban pollution levels can influence lipid metabolism in the mammary gland and the lipid and protein composition of milk. As milk may contribute to metabolic programming, such alterations affecting milk composition should be taken into account from a public health perspective.

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

OPEN ACCESS

**Citation:** Hue-Beauvais C, Aujean E, Miranda G, Ralliard-Rousseau D, Valentino S, Brun N, et al. (2019) Impact of exposure to diesel exhaust during pregnancy on mammary gland development and milk composition in the rabbit. PLoS ONE 14(2): e0212132. https://doi.org/10.1371/journal.pone.0212132

**Editor:** Alexander Larcombe, Telethon Institute for Child Health Research, AUSTRALIA

**Received:** September 27, 2018

**Accepted:** January 28, 2019

**Published:** February 14, 2019

**Copyright:** © 2019 Hue-Beauvais et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** The ERC program e-DOHAD and the Agence Nationale de la Recherche (ANR, France) project EPPAP (PI: Dr. Pascale Chavatte-Palmer) funded the animal experiments. EPPAP also covered the costs of mammary gland analyses. The funders had no role in study design, data
collection and analysis, decision to publish, or
preparation of the manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

**Introduction**

Diesel engine exhaust (DE) is an environmental risk factor for humans since it has been linked to various deleterious effects such as increased morbidity and mortality due to air pollution [1, 2] and increased risk of lung cancer [3] and respiratory diseases [4], reproductive disorders [5, 6] and genotoxicity affecting various tissues [7, 8]. Diesel engine exhaust is a complex mixture of gaseous and particulate components. The particles released by DE are usually fine, with an aerodynamic diameter ≤2.5 µm, as well as a very high proportion of ultrafine particles with a diameter of less than 0.1 µm, which are probably more harmful than the larger fraction [9]. In the environment, they combine very rapidly into larger groups and aggregate. Because the particles are very small, their overall surface is very large, thus enabling the adsorption of a variety of compounds.

The mammary network develops over a long period of time, being initiated during foetal life and continuing until the first lactation in adulthood. During pregnancy, the mammary gland develops to become a branched epithelial network of ducts that support large alveolar development and allow subsequent milk production during lactation [10]. Mammary epithelial growth and differentiation are tightly modulated throughout development by several hormonal and metabolic signals [11], suggesting the existence of critical periods during which the mammary gland is vulnerable to developmental programming by environmental stressors, as explained in the developmental origins of health and diseases (DOHaD) concept [12, 13].

An increasing number of studies in both human and animal models have shown that prenatal exposure to ambient fine particles is associated with a variety of abnormalities during foetal development, such as preterm birth, low birth weight and placental dysfunctions [14–17]. Moreover, there is clear evidence that gestational exposure to pollutants can also have deleterious long term effects that mainly involve energy metabolism disorders [18, 19] and cardiac and respiratory dysfunctions [4, 20, 21]. We have previously reported that females exposed daily during pregnancy to DE particles at levels close to urban pollution evidenced foeto-placental disorders [14]. Rabbits were used for these studies as they can be considered as a relevant model for human physiology due to their similarities in lipid metabolism [22] and in placentation closer to humans than that of rodents [23] as well as because the rabbit-primate phylogenetic distance is the same as the rodent-primate distance [24]. Moreover, because rodent sequences have evolved more rapidly, rabbit gene sequences are more similar to those of humans than to those of rodents, which might be of considerable importance to further analyses of altered gene expression profiles [24].

Since early postnatal life is considered to be a period that is critical to the long-term health and disease susceptibility of offspring, the impact of pollutants on mammary gland development and subsequent milk composition is of crucial importance. Rabbit milk is characterized by a high fat and protein content, an absence of β-lactoglobulin, the presence of whey acidic protein (WAP), the existence of a second α<sub>s2</sub> casein (α<sub>s2</sub>-like-casein) and a high level of lactoferrin [25, 26]. As milk was shown to exert effects far beyond its immediate nutritional value, the objective of the present study was therefore to determine the effects of maternal exposure to DE during pregnancy on mammary gland development in terms of tissue morphology and function. In addition, we investigated the impact of DE exposure on milk composition, including the fatty acid content and relative proportions of the main milk proteins. Furthermore, because leptin plays a key role in regulating developmental programming [27], we examined the leptin levels in mammary tissue and different milk samples in order to link maternal DE exposure to a potential for subsequent risks of growth and metabolic disorders in their offspring.
# Materials and methods

## Ethics statement

All animal studies were performed in compliance with European Community regulations on animal experimentation (European Communities Council Directive 86/609/EEC; European Union, 1986) and with the authorisation of the French Ministry of Agriculture. All protocols were approved (visa 12-102) by the Comité d'Ethique en Expérimentation Animale (Ethics Committee on Animal Experimentation), registered as N°C2EA-45 by the French National register CNREEA (Comité National de Reflexion Ethique sur l'Expérimentation Animale by the Comité National de Réflexion Ethique sur l'Expérimentation Animale, French National Ethics Committee on Animal Experimentation).

## Animal experiments and sample collection

The rabbits used during this study were the same as those described in detail elsewhere [14]. Briefly, one year old female rabbits (New Zealand White, 1077-INRA) were housed individually in an indoor facility under controlled light conditions (8 h light/16 h darkness except from the week before mating when they were under 16 h light/8 h darkness) and temperature (18°C). In their first pregnancy rabbits were exposed to diluted (1 mg/m³) and filtered diesel engine (DE) exhaust fumes (nanoparticle diameter ~ 69 nm) or clean air (controls) for 2h/day, 5 days/week by means of nose-only exposure from the 3rd to the 27th days post-conception (total exposure of 20 days out of a 31-day pregnancy) as previously described [14]. DE exposure was performed with the Mobile Ambient Particle Concentrator Exposure Laboratory [28] connected to a 25KVA Loxam engine, with a 500 nm particle filter. DE is a complex mixture of hundreds of constituents in either a gas or particle form. Gaseous components of DE include carbon dioxide, oxygen, nitrogen, water vapour, carbon monoxide, nitrogen compounds, sulphur compounds, and numerous low-molecular-weight hydrocarbons (some of them individually known to be toxic, such as aldehydes, benzene, 1,3-butadiene, polycyclic aromatic hydrocarbons (PAHs) and nitro-PAHS). The particles present in DE are known to be composed of center core of elemental carbon with absorbed organic compounds and small amounts of sulphate, nitrate, metals, and other trace elements [29].

On Day 28 of pregnancy, before parturition, 12 dams (six exposed to DE and six controls (C)) were sacrificed by exsanguination after electronarcosis and the left lower mammary gland from each animal was excised and dissected. As fetal organs were collected for further examination, fetuses were euthanised by decapitation after removal of the uterus in accordance with the Ethics Committee on Animal Experimentation as well as the French Ministry of Agriculture (visa 12-102). Dams' mammary epithelial tissue was either processed for histology or snap-frozen and stored for later RNA isolation.

Milk samples were collected by hand stripping from 18 other females dedicated to post-natal follow-up of the offspring (nine exposed to DE and nine C) on days 2-4 and days 13-16 of lactation, as previously described [30]. The samples thus collected were stored at -20°C until analysis. Two groups were thus formed and named as a function of the prior exposure of the dams during pregnancy (C milk and DE milk, respectively).

## Histological and ultrastructural analyses of mammary tissue

For histological analyses, mammary samples were fixed in 4% paraformaldehyde, cryoprotected in 40% sucrose and embedded in TissuTek (Sakura, Torrance, CA). Five micrometre sections, at least 100 µm apart, were mounted on slides, stained with hematoxylin and eosin (H&E; Sigma-Aldrich, Saint Quentin Fallavier, France) and digitised using a Hamamatsu
NanoZoomer (Hamamatsu Photonics, Tokyo, Japan). Five sections per rabbit were processed and the images were analysed using both the Hamamatsu NanoZoomer Digital Pathology Virtual Slide Viewer and ImageJ software [31]. The areas occupied by mammary epithelial tissue, adipose tissue, mammary ducts or connective tissue were measured.

For transmission electron microscopy, mammary tissue was fixed with 2% glutaraldehyde in 0.1M Na cacodylate buffer pH 7.2 for 4 hours at room temperature. The samples were then contrasted with 0.5% Oolong Tea Extract (OTE) in 0.1M Na cacodylate buffer, post-fixed with 1% osmium tetroxide containing 1.5% potassium cyanoferate, gradually dehydrated in ethanol (30% to 100%) and gradually substituted in a mix of propylene oxyde-epon and embedded in Epon. (Delta Microscopie, Labège, France). Half-micron sections were collected onto glass slides, counterstained with methylene blue-Azur II and imaged under an optical microscope.

The surface occupied by luminal fat globules was measured in relation to the total area of the lumen, using a threshold method with the Fiji module of ImageJ software. A total of 88 lumina for C mammary glands and 108 lumina for DE-exposed mammary glands were processed.

## Immunolocalisation of leptin in the mammary gland

Leptin was localised by immunohistochemical analyses using a 1:100 dilution of goat anti-rat leptin antibody (AbD Serotec, France). Briefly, 5 µm sections were treated in 50 mM ammonium chloride for 30 min followed by permeabilisation in 2% BSA, 0.05% saponin and 0.05% sodium azide in PBS 1X for 1 h. The tissue sections were then incubated overnight at 4°C with primary antibody. Antibody binding was visualised after incubation for 45 min with fluorescence-labelled anti-goat TRITC-conjugated antibody (dilution: 1:300, Jackson Immunoresearch). The slides were then mounted in Vectashield mounting medium (Vector Laboratories, Burlingame, CA) and observed with a Zeiss microscope coupled with ApoTome technology.

## Reverse transcription and quantitative PCR

For gene assays in the mammary gland, reverse transcription was performed on 500 ng of total RNA using the SuperScript1 VILO cDNA Synthesis kit according to the manufacturer's instructions (Invitrogen, France) and under the following conditions: 42°C for 60 min and 85°C for 5 min. Quantitative PCR runs were achieved using the ABsolute Blue QPCR Mix, SYBR Green1 (Thermo Scientific, France) according to the manufacturer's instructions, on a Mastercycler Eppendorf Realplex system (Eppendorf, France) under the following conditions: 95°C for 15 min, 45 cycles of 95°C for 15 sec and 60°C for 1 min, and a melting curve.

The threshold cycles obtained for Leptin (F 5' -ACACAGCAGTCGGTCTC-3' and R 5' - GAGGTTCTCCAGGTCGTTGG-3'), SREBP1 (Sterol Regulatory Element-Binding Protein 1, F 5' - CCAGCTGACAGCTCCATTGA-3' and R 5' - TGCGCGCCACAAGGA -3'), FASN (Fatty Acid Synthase N, F 5' -ACCTCGTGAAGGCTGTGACTCA-3' and R 5'-TGAGTCGAGGCCAAGGTCTGAA-3'), SCD (Stearoyl-CoA Desaturase, F 5'-TTATTCCGTTATGCCCTTGG-3' and R 5'-TTGTCATAAGGGCGGTATCC-3') and LPL (LipoProtein Lipase, F 5'-CTCAGGACTCCCGAAGACAC-3' and R 5'-GTTTTGCTGCTGTGGTTGAA-3') were normalised with the values of TBP (TATA Binding Protein F 5'-TGACCCCCATGACCCCTATT-3' and R 5'-CAGCAAACCGCTTGGGATTA-3') and the results were expressed as fold changes of the threshold cycle (Ct) values relative to the control using the 2-ΔΔCt method [32].

## Analysis of the milk fatty acid composition

Analyses were performed on 18 individual milk samples from DE and C dams (9 C and 9 DE), collected on days 2-4 (L2-4) and days 13-16 (L13-16) of lactation. To determine the fatty acid
composition, 200-µl full-fat milk were added to an internal standard (margaric acid, C17:0) before lipid extraction with chloroform/methanol (2:1, adapted from the method described by Folch et al. [33]). The fatty acids were transmethylated with 7% boron trifluoride methanol (Sigma-Aldrich, Saint Quentin Fallavier, France) in accordance with the method previously described by Morrison and Smith [34]. Finally, the methyl esters of milk fatty acids were analysed using gas chromatography (Auto Sampling 8410 Gas Chromatograph 3900; Varian, Les Ulis, France) coupled to a flame ionisation detector on an Econo-Cap EC-WAX capillary column (30 m, 0.32-mm internal diameter, 0.25-µm film, reference 19654; ALLTECH Associates Inc., Templemars, France), as described by Rousseau et al.[35]. Identification of FA was made in reference to known FA profiles obtained after injection of standard FAME (FA methyl esters) mix (Supelco 37 components FAME mix, ref 47885-U, Sigma). The fatty acids profile was established for each sample and expressed as percentage of total FA. Heptadecanoic acid (C17:0), introduced prior to milk lipid extraction (50µg) was used as an internal standard to measure the total FA concentration after chromatogram analysis.

## Analysis of the main proteins in milk

Analyses were performed on the same 18 individual milk samples from DE and C dams. The milk samples were skimmed by centrifugation at 2000 × g for 20 min at 4°C, then analysed as previously described [36] using reverse-phase high-performance liquid chromatography (HPLC) coupled with electrospray ionisation (ESI) mass spectrometry. Peak profiles from UV absorbance at 214 nm were analysed using Chromeleon software (Chromleon 7.0.0; Dionex, Thermo Fisher Scientific). Protein isoforms of the principal milk proteins were identified by comparing their molecular masses with an in-house mass database on rabbit milk proteins and the UniProt database (UniProt Consortium, 2015; http://www.uniprot.org [Accessed July 2015]). Moreover, the main milk proteins were quantified by integrating the chromatogram peak areas at 214 nm. Values were expressed as a percentage of the total peak area.

## Measurement of leptin concentrations in milk

Quantification of the levels of rabbit leptin in each milk sample was performed using an ELISA kit (Cloud Clone Corp., Houston, TX) according to the manufacturer's instructions. Triplicates of each milk sample were used and the results reported here represent the mean of these three measurements.

## Statistical analysis

The data were analysed using a nonparametric Mann-Whitney U-test and significant differences were defined as $P<0.05$. When necessary, statistical analysis was supplemented by an ANOVA approach.

## Results

### Maternal characterisation

Exposed females exhibited no alteration either concerning their plasma fatty acid profile or biochemical data. Moreover, within pups, sex ratio or litter size were similar among groups [14].

### DE Exposure throughout pregnancy does not alter mammary gland histology but increases the alveolar fat globule content

Histological examinations were performed to determine whether DE exposure during pregnancy induced morphological changes to the mammary gland. Morphological differences
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Quantification by stereology of the surface occupied by mammary ducts' and alveoli's lumina. Analyses were performed in C and DE-exposed rabbits at Day 28 of pregnancy using the Hamamatsu NanoZoomer Digital Pathology Virtual Slide Viewer and ImageJ software. Data are mean ± SEM.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0212132.g001

were quantitatively analysed based on the relative ratio of the surfaces of mammary ducts and alveoli's lumina. These analyses did not reveal any significant differences between DE and control animals on Day 28 of pregnancy (Fig 1).

To further investigate the alveolar content, thin mammary tissue sections were performed in control (Fig 2A) and DE (Fig 2B) dams, and the relative surface area occupied by fat globules was measured in the alveolar lumina. A clear modification was observed in the alveolar lumina of DE rabbits. Indeed, all DE mammary glands were characterised by a significant increase in the total surface area occupied by fat globules (23.80 ± 9.37% vs. 38.72 ± 2.40%, p = 0.008) as the representative shown in Fig 2C and 2D.

<figure><img src="image_6.png" /><figcaption>**Fig 2.** Histological sections of mammary gland from C and DE rabbits at Day 28 of pregnancy. Hematoxylin and eosin staining of 5µm cryosections of mammary glands from C (A and C) and DE-exposed (B and D) rabbits were observed at different magnifications. Representative nanozoomer scans were observed at low magnification (A and B: the scale bar represents 2 mm) and higher magnification (C and D: the scale bars represent 250 µm).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0212132.g002
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Leptin expression and localisation in mammary glands from C and DE rabbits at day 28 of pregnancy. (A) Expression of leptin in mammary glands at day 28 of pregnancy from control (in light grey) and exposed (DE) (in black) dams. Real-time RT-PCR was quantified using TBP housekeeping gene. Data are mean ± SEM. (B and C) Representative images of leptin immunostaining (red) of mammary gland from control (B) or exposed (C) dams at day 28 of pregnancy. Nuclei are stained in green. Scale bar represents 50µm.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0212132.g003

## DE Exposure during pregnancy does not modify leptin expression and protein localisation in the mammary gland

To examine the expression of leptin mRNA and protein localisation in the mammary glands of C and DE-exposed dams on Day 28 of pregnancy, RT-qPCR analyses and immunostaining experiments were performed (Fig 3). No significant differences were observed between the two groups of animals with respect to mRNA expression, although there were some signs of a decrease in DE mammary glands. Immunolocalisation studies showed that the epithelial localization of leptin in mammary tissue was not altered following DE exposure.

## Impact of DE Exposure on the Expression of the LPL, FASN, SREBP and SCD Genes Involved in Lipid Metabolism in the mammary gland

To determine whether DE exposure might affect the expression of key genes involved in lipid metabolism, RT-qPCR experiments were performed on total mammary RNA (Fig 4). Among the four genes tested, only SCD was significantly increased in DE animals compared to controls (SCD/TBP ratio: 2.48 ± 0.80 vs. 5.43 ± 1.54, p = 0.021), suggesting increased fatty acid synthesis in the mammary gland of DE rabbits on Day 28 of pregnancy.

## DE exposure alters the fatty acid profile of milk

To evaluate DE exposure-induced alterations to milk composition, fatty acid profiles were determined in milk from C and DE-exposed dams at two stages of lactation: early (2–4 days, L2-L4) and mid (13–16 days, L13-L16) lactation (Tables 1 and 2).

During early lactation (L2-L4), the total fatty acid content was higher in DE milk than in C milk (99.6 ± 5.7 g/L and 75.2 ± 3.3 g/L, *P*<0.05, respectively). Quantitatively, DE milk displayed a higher total SFA content than C milk. This rise was due to all SFA except stearic acid (C18:0) which concentration did not vary significantly. Similarly, a sharp increase in the MUFA content was observed in DE milk, alongside a corresponding rise in palmitoleic and oleic acids (*P*<0.05). Polyunsaturated fatty acid levels were also higher in DE milk, with increases in both α-linolenic and linoleic acids (*P*<0.05) although there was no change to the
<figure><img src="image_8.png" /><figcaption>Fig 4. Expression of lipid metabolism enzymes in the mammary glands of C and DE rabbits determined using real-time RT-PCR. (A) Lipoprotein lipase (LPL), (B) Stearoyl CoA desaturase (SCD), (C) Sterol regulatory element-binding protein (SREBP) and (D) Fatty Acid Synthase N. Expression was measured in mammary tissue from control (in light grey) and DE exposed (black) dams, at 28 days of pregnancy. Data are mean ± SEM.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0212132.g004

PUFA *n*-6: *n*-3 ratio (Table 1). Qualitatively, SFA were the main fatty acids present in both DE and C milks (Table 1). Both short chain and long chain SFA levels increased in DE milk (*P*<0.05). Among MUFA, only palmitoleic acid was present at a higher level in DE milk (0.32 ± 0.02%) than in C milk (0.26 ± 0.01%).

At mid-lactation (L13-16), no quantitative differences in the fatty acid content were observed between C and DE milks. In contrast, a slight qualitative rise in SFA levels (3%) was seen in DE milk, resulting from an increase in C10:0 and a fall of C18:0 SFA (Table 2). No qualitative modifications were observed in the proportions of MUFA, but the PUFA content was reduced in DE milk and led to a decrease in the *n*-6: *n*-3 ratio.

Thus, at both stages of lactation, exposure to diesel exhaust throughout pregnancy altered the fatty acid composition, although mainly during early lactation.

## The composition of the main proteins is altered in DE milk

To investigate the impact of DE exposure during pregnancy on the composition of the main milk proteins, LC-MS technology was used to analyse milk samples at two stages of lactation.
<table><caption>Table 1. Fatty acid composition and content in milk from control (C) and DE-exposed (DE) dams on days 2–4 of lactation.</caption><thead><tr><th>Fatty acids</th><th>C milk<sup>1</sup></th><th>DE milk<sup>1</sup></th><th>C milk<sup>2</sup></th><th>DE milk<sup>2</sup></th></tr></thead><tbody><tr><td>C8:0</td><td>17.77 (0.85)<sup>a</sup></td><td>15.83 (0.73)<sup>a</sup></td><td>13.32 (0.66)<sup>a</sup></td><td>15.57 (0.88)<sup>a</sup></td></tr><tr><td>C10:0</td><td>19.99 (1.00)</td><td>18.76 (0.80)</td><td>14.96 (0.70)<sup>a</sup></td><td>18.81 (1.12)<sup>a</sup></td></tr><tr><td>C12:0</td><td>3.74 (0.31)</td><td>3.68 (0.29)</td><td>2.80 (0.23)<sup>a</sup></td><td>3.79 (0.33)<sup>a</sup></td></tr><tr><td>C14:0</td><td>1.94 (0.06)</td><td>1.88 (0.10)</td><td>1.46 (0.10)<sup>a</sup></td><td>1.92 (0.15)<sup>a</sup></td></tr><tr><td>C16:0</td><td>15.49 (0.63)<sup>a</sup></td><td>16.79 (0.51)<sup>a</sup></td><td>11.71 (0.90)<sup>a</sup></td><td>16.76 (1.26)<sup>a</sup></td></tr><tr><td>C18:0</td><td>3.21 (0.25)</td><td>2.89 (0.12)</td><td>2.43 (0.25)</td><td>2.88 (0.23)</td></tr><tr><td><strong>SFA</strong></td><td><strong>62.67 (1.17)</strong></td><td><strong>60.30 (1.07)</strong></td><td><strong>47.08 (1.92)<sup>a</sup></strong></td><td><strong>60.20 (3.44)<sup>a</sup></strong></td></tr><tr><td><strong>sc SFA</strong></td><td><strong>41.51 (1.77)</strong></td><td><strong>38.27 (1.49)</strong></td><td><strong>31.07 (1.28)<sup>a</sup></strong></td><td><strong>38.17 (2.17)<sup>a</sup></strong></td></tr><tr><td><strong>mlc SFA</strong></td><td><strong>21.16 (0.87)</strong></td><td><strong>22.03 (0.58)</strong></td><td><strong>16.00 (1.25)<sup>a</sup></strong></td><td><strong>22.03 (1.63)<sup>a</sup></strong></td></tr><tr><td>C16:1 n-9</td><td>0.26 (0.01)<sup>a</sup></td><td>0.32 (0.02)<sup>a</sup></td><td>0.20 (0.01)<sup>a</sup></td><td>0.32 (0.02)<sup>a</sup></td></tr><tr><td>C18:1 n-9</td><td>16.42 (0.50)</td><td>17.39 (0.53)</td><td>12.37 (0.75)<sup>a</sup></td><td>17.24 (1.16)<sup>a</sup></td></tr><tr><td><strong>MUFA</strong></td><td><strong>19.56 (0.68)</strong></td><td><strong>20.96 (0.79)</strong></td><td><strong>14.74 (0.97)<sup>a</sup></strong></td><td><strong>20.90 (1.49)<sup>a</sup></strong></td></tr><tr><td>C18:2 n-6</td><td>14.72 (0.48)</td><td>15.16 (0.48)</td><td>11.10 (0.73)<sup>a</sup></td><td>15.00 (1.12)<sup>a</sup></td></tr><tr><td>C18:3 n-3</td><td>2.15 (0.05)</td><td>2.10 (0.05)</td><td>1.62 (0.10)<sup>a</sup></td><td>2.09 (0.15)<sup>a</sup></td></tr><tr><td><strong>PUFA</strong></td><td><strong>17.77 (0.52)</strong></td><td><strong>18.74 (0.56)</strong></td><td><strong>13.40 (0.84)<sup>a</sup></strong></td><td><strong>18.49 (1.20)<sup>a</sup></strong></td></tr><tr><td>n-6 PUFA</td><td>15.57 (0.48)</td><td>16.18 (0.52)</td><td>11.74 (0.75)<sup>a</sup></td><td>15.99 (1.22)<sup>a</sup></td></tr><tr><td>n-3 PUFA</td><td>2.20 (0.05)</td><td>2.18 (0.05)</td><td>1.66 (0.10)<sup>a</sup></td><td>2.19 (0.16)<sup>a</sup></td></tr><tr><td>n-6/n-3 ratio</td><td>7.07 (0.13)</td><td>7.43 (0.21)</td><td></td><td></td></tr><tr><td><strong>Total FA</strong></td><td>100</td><td>100</td><td>75.22 (3.32)<sup>a</sup></td><td>99.59 (5.70)<sup>a</sup></td></tr></tbody></table>

A qualitative comparison of the chromatogram profiles of milk from DE-exposed dams on Days 2–4 and 13–14 of lactation (S1 Fig) did not reveal any differences regarding the presence of the main milk proteins, nor their variants or phosphorylated isoforms [36].

A quantitative estimate of variations in the peak areas was performed as a function of the relative proportions of whole milk proteins; i.e. the percentage of the total peak areas on the UV chromatogram at 214 nm (Table 3). No variations in the total casein ratio were observed at either the L2-L4 or L13-L16 days of lactation in DE versus C milk. However, among the casein proteins, a significant increase in α<sub>s2</sub>-like proteins (P <0.05) and a significant decrease in α<sub>s1</sub> proteins (P<0.05) were observed. As for total whey proteins, the proportions of both lactoferrin and serum albumin had risen in DE milk at days L2-L4 (P<0.05), whereas only lactoferrin increased in L13-L16 milk samples. Moreover, the total WAP proportion did not differ between C and DE milks.

When the relative proportions of the main milk proteins were compared between Days 2–4 and 13–16 of lactation in C and DE-exposed rabbits, only α-lactalbumin followed a different pattern (Table 3). While α-lactalbumin remained stable in C milk, its level rose sharply in DE milks (0.041 ± 0.07% vs. 0.76 ± 0.09, respectively, P<0.05). For α<sub>s2</sub> and β-caseins, an increase was observed in both the C and DE-exposed groups. Alpha-s1 casein levels fell between L2-L4 and L13-L16 (P<0.05).

The relative proportion of phosphorylated isoforms (percentage of the mass signal intensity of each phosphorylation isoform of a given protein relative to the sum of all its phosphorylation isoforms) of α<sub>s1</sub>, α<sub>s2</sub>-like and β-caseins and WAP in C and DE milks was determined at the two stages of lactation (Table 4). On Days 2–4 and 13–16 of lactation, the relative

<footnote n="1">¹Percentage of total fatty acids.</footnote>

<footnote n="2">²Grams per litre.</footnote>

<footnote>ª P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups.</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0212132.t001</footnote>
<table><caption>Table 2. Fatty acid composition and content in milk from control (C) and DE-exposed (DE) dams on days 13–16 of lactation.</caption><thead><tr><th>Fatty acids</th><th>C milk<sup>1</sup></th><th>DE milk<sup>1</sup></th><th>C milk<sup>2</sup></th><th>DE milk<sup>2</sup></th></tr></thead><tbody><tr><td>C8:0</td><td>25.56 (0.93)</td><td>26.06 (0.84)</td><td>18.59 (1.86)</td><td>21.44 (1.10)</td></tr><tr><td>C10:0</td><td>20.40 (0.64)<sup>a</sup></td><td>22.66 (0.31)<sup>a</sup></td><td>15.13 (1.81)</td><td>18.73 (1.03)</td></tr><tr><td>C12:0</td><td>2.47 (0.15)<sup>a</sup></td><td>2.97 (0.15)<sup>a</sup></td><td>1.87 (0.28)</td><td>2.47 (0.20)</td></tr><tr><td>C14:0</td><td>1.47 (0.08)</td><td>1.52 (0.13)</td><td>1.11 (0.15)</td><td>1.27 (0.14)</td></tr><tr><td>C16:0</td><td>11.71 (0.90)</td><td>12.44 (0.65)</td><td>9.74 (1.26)</td><td>10.41 (1.00)</td></tr><tr><td>C18:0</td><td>3.05 (0.13)<sup>a</sup></td><td>2.62 (0.10)<sup>a</sup></td><td>2.26 (0.29)</td><td>2.17 (0.16)</td></tr><tr><td><strong>SFA</strong></td><td><strong>66.44 (0.86)<sup>a</sup></strong></td><td><strong>68.64 (0.47)<sup>a</sup></strong></td><td><strong>48.99 (5.38)</strong></td><td><strong>56.82 (3.27)</strong></td></tr><tr><td><strong>sc SFA</strong></td><td><strong>48.43 (1.44)</strong></td><td><strong>51.68 (0.80)</strong></td><td><strong>35.59 (3.91)</strong></td><td><strong>42.64 (2.17)</strong></td></tr><tr><td><strong>mlc SFA</strong></td><td><strong>18.02 (0.68)</strong></td><td><strong>16.96 (0.74)</strong></td><td><strong>13.40 (1.71)</strong></td><td><strong>14.17 (1.27)</strong></td></tr><tr><td>C16:1 n-9</td><td>0.20 (0.01)</td><td>0.32 (0.02)</td><td>0.19 (0.02)</td><td>0.23 (0.03)</td></tr><tr><td>C18:1 n-9</td><td>15.20 (0.42)</td><td>14.83 (0.22)</td><td>11.33 (1.42)</td><td>12.28 (0.77)</td></tr><tr><td><strong>MUFA</strong></td><td><strong>17.83 (0.50)</strong></td><td><strong>17.35 (0.31)</strong></td><td><strong>13.24 (1.62)</strong></td><td><strong>14.38 (0.93)</strong></td></tr><tr><td>C18:2 n-6</td><td>13.05 (0.38)<sup>a</sup></td><td>11.54 (0.22)<sup>a</sup></td><td>9.71 (1.22)</td><td>9.52 (0.50)</td></tr><tr><td>C18:3 n-3</td><td>1.80 (0.02)</td><td>1.73 (0.03)</td><td>1.34 (0.16)</td><td>1.43 (0.07)</td></tr><tr><td><strong>PUFA</strong></td><td><strong>15.73 (0.40)<sup>a</sup></strong></td><td><strong>14.01 (0.25)<sup>a</sup></strong></td><td><strong>11.70 (1.45)</strong></td><td><strong>11.56 (0.63)</strong></td></tr><tr><td>n-6 PUFA</td><td>13.90 (0.39)<sup>a</sup></td><td>12.23 (0.23)<sup>a</sup></td><td>10.33 (1.28)</td><td>10.09 (0.56)</td></tr><tr><td>n-3 PUFA</td><td>1.83 (0.03)</td><td>1.78 (0.03)</td><td>1.36 (0.17)</td><td>1.47 (0.08)</td></tr><tr><td>n-6/n-3 ratio</td><td>7.60 (0.18)<sup>a</sup></td><td>6.86 (0.12)<sup>a</sup></td><td></td><td></td></tr><tr><td><strong>Total FA</strong></td><td>100</td><td>100</td><td>73.94 (8.32)</td><td>82.75 (4.75)</td></tr></tbody></table>

proportions of the less phosphorylated isoforms (1-4P) were higher in DE milk for α<sub>s2</sub>-like-casein (*P*<0.05). However, at the later stage, the most phosphorylated isoform of α<sub>s1</sub>-casein displayed a significant decrease in DE milk (*P*<0.05).

## DE Exposure does not change leptin concentrations in milk

The immunodetection of leptin was performed in whole milk samples collected on days 2–4 and 13–16 of lactation (Fig 5). Leptin concentrations were similar in milk from both C and DE-exposed females at both stages of lactation.

## Discussion

This study revealed the effects on mammary gland development at late pregnancy (Day 28) induced in rabbits exposed to diesel exhaust fumes (1 mg/m³, for 2h/day, 5 days/week), throughout pregnancy. This diesel exhaust level was chosen as this experiment was part of a large experimentation aiming at evaluating effect of maternal DE exposure at concentrations close to that of human exposure in large European cities when people drive, ride or walk on big roads twice a day [14]. Exposure to DE did not alter mammary gland histology, but increased the alveolar content in fat globules. By Day 28 of pregnancy, the rabbit mammary gland has entered a differentiated stage of growth [37] after completing the proliferative stage. Although the total surface area occupied by mammary epithelial tissue and lumina was not affected by DE exposure, the relative surface area of fat globules increased dramatically in the mammary lumina of DE-exposed animals. The regulation of adiposity is deemed to be the

<footnote n="1">¹Percentage of total fatty acids.</footnote>

<footnote n="2">²Grams per litre.</footnote>

<footnote>ª P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0212132.t002</footnote>
<table><caption>Table 3. Quantitative analysis of major milk proteins.</caption><thead><tr><th rowspan="2">Protein</th><th colspan="2">L2-L4</th><th colspan="2">L13-L16</th></tr><tr><th>C</th><th>DE</th><th>C</th><th>DE</th></tr></thead><tbody><tr><td>α<sub>s1</sub>-casein</td><td>27.60 (0.6)<sup>a,c</sup></td><td>25.70 (0.39)<sup>a,d</sup></td><td>22.50 (0.32)<sup>b,c</sup></td><td>20.80 (0.46)<sup>b,d</sup></td></tr><tr><td>α<sub>s2</sub>-casein</td><td>1.30 (0.20)<sup>c</sup></td><td>1.10 (0.18)<sup>d</sup></td><td>3.60 (0.18)<sup>c</sup></td><td>3.50 (0.18)<sup>d</sup></td></tr><tr><td>α<sub>s2</sub>-like-casein</td><td>13.00 (1.18)<sup>a</sup></td><td>16.80 (0.8)<sup>a</sup></td><td>12.70 (0.57)<sup>b</sup></td><td>16.20 (0.35)<sup>b</sup></td></tr><tr><td>β-casein</td><td>26.30 (0.56)<sup>c</sup></td><td>24.80 (0.71)<sup>d</sup></td><td>28.60 (0.28)<sup>c</sup></td><td>27.10 (0.59)<sup>d</sup></td></tr><tr><td>κ-casein</td><td>8.40 (0.28)</td><td>7.80 (0.33)</td><td>9.00 (0.3)</td><td>8.50 (0.15)</td></tr><tr><td><strong>Total Caseins</strong></td><td><strong>76.60 (0.63)</strong></td><td><strong>76.00 (0.85)</strong></td><td><strong>76.40 (0.25)</strong></td><td><strong>76.00 (0.47)</strong></td></tr><tr><td>WAP</td><td>15.40 (0.54)</td><td>14.70 (0.55)</td><td>15.70 (0.37)</td><td>15.60 (0.34)</td></tr><tr><td>Lactoferrin</td><td>6.00 (0.17)<sup>a</sup></td><td>6.60 (0.20)<sup>a</sup></td><td>5.90 (0.17)<sup>b</sup></td><td>6.50 (0.21)<sup>b</sup></td></tr><tr><td>α-Lactalbumin</td><td>0.61 (0.07)</td><td>0.41 (0.07)<sup>d</sup></td><td>0.85 (0.13)</td><td>0.76 (0.09)<sup>d</sup></td></tr><tr><td>Serumalbumin</td><td>1.40 (0.14)<sup>a</sup></td><td>2.10 (0.34)<sup>a</sup></td><td>1.10 (0.18)</td><td>1.40 (0.24)</td></tr><tr><td><strong>Total Whey proteins</strong></td><td><strong>23.40 (0.63)</strong></td><td><strong>23.90 (0.25)</strong></td><td><strong>23.60 (0.85)</strong></td><td><strong>24.30 (0.45)</strong></td></tr></tbody></table>

Quantitative results are expressed as a % of total proteins.

https://doi.org/10.1371/journal.pone.0212132.t003

primary function of leptin, but it also plays a role in other physiological processes, as evidenced by its many sites of synthesis other than fat cells [38]. In particular, mammary epithelial cells produce leptin [38–41], which acts as an autocrine and paracrine factor to influence development and differentiation of the mammary gland [42]. Neither mammary leptin mRNA expression nor protein localisation in mammary epithelial cells differed between C and DE-exposed rabbits at Day 28 of pregnancy. Leptin transcripts, however, tended to be decreased in DE

<table><caption>Table 4. Phosphorylation of major milk proteins from C and DE rabbits.</caption><thead><tr><th rowspan="2">Proteins</th><th rowspan="2">Phosphorylation</th><th colspan="2">L2-L4</th><th colspan="2">L13-L16</th></tr><tr><th>C</th><th>DE</th><th>C</th><th>DE</th></tr><tr><td></td><td></td><th>Mean (sem)</th><th>Mean (sem)</th><th>Mean (sem)</th><th>Mean (sem)</th></tr></thead><tbody><tr><td rowspan="2">α<sub>s1</sub>-casein</td><td>5-6P</td><td>9.5 (1.25)</td><td>9.4 (0.59)</td><td>9.8 (0.25)</td><td>9.6 (0.22)</td></tr><tr><td>7-8P</td><td>16 (1.69)<sup>b</sup></td><td>14.0 (0.93)<sup>b</sup></td><td>12.3 (0.46)<sup>a,b</sup></td><td>9.7 (0.3)<sup>a,b</sup></td></tr><tr><td rowspan="2">α<sub>s2</sub>-like-casein</td><td>1-4P</td><td>3.5 (0.75)<sup>a</sup></td><td>7.1 (0.78)<sup>a</sup></td><td>3.9 (0.7)<sup>b</sup></td><td>7.6 (0.89)<sup>b</sup></td></tr><tr><td>5-8P</td><td>8.4 (0.96)</td><td>8.3 (1.19)</td><td>8.7 (0.26)</td><td>7.4 (0.86)</td></tr><tr><td rowspan="2">α<sub>s2</sub>-casein</td><td>2-3P</td><td>0.5 (0.1)<sup>c</sup></td><td>0.5 (0.08)<sup>d</sup></td><td>1.4 (0.07)<sup>c</sup></td><td>1.3 (0.06)<sup>d</sup></td></tr><tr><td>4P</td><td>0.8 (0.12)<sup>c</sup></td><td>0.6 (0.09)<sup>d</sup></td><td>2.2 (0.15)<sup>c</sup></td><td>1.9 (0.14)<sup>d</sup></td></tr><tr><td rowspan="2">β-casein</td><td>0-3P</td><td>7.2 (0.83)<sup>c</sup></td><td>7.6 (0.68)<sup>d</sup></td><td>9.7 (0.08)<sup>c</sup></td><td>9.6 (0.57)<sup>d</sup></td></tr><tr><td>4P</td><td>14.8 (0.81)</td><td>16.1 (0.67)</td><td>15.3 (0.41)</td><td>16.2 (0.49)</td></tr><tr><td rowspan="2">WAP</td><td>0-1P</td><td>4.2 (0.77)</td><td>3.4 (0.41)<sup>d</sup></td><td>4.4 (0.26)</td><td>4.7 (0.18)<sup>d</sup></td></tr><tr><td>2-3P</td><td>11.2 (0.62)</td><td>11.3 (0.59)</td><td>11.3 (0.32)</td><td>10.9 (0.39)</td></tr></tbody></table>

Quantitative results are expressed as a % of total proteins.

https://doi.org/10.1371/journal.pone.0212132.t004

<footnote>a P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups on days 2–4 of lactation.</footnote>

<footnote>b P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups on days 13–16 of lactation.</footnote>

<footnote>c P values <0.05 calculated using the Mann-Whitney test to compare the C milk at L2-L4 and L13-16 days of lactation.</footnote>

<footnote>d P values <0.05 calculated using the Mann-Whitney test to compare the DE exposed milk at L2-L4 and L13-16 days of lactation.</footnote>

<footnote>a P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups on days 2–4 of lactation.</footnote>

<footnote>b P values <0.05 calculated using the Mann-Whitney test to compare the C and DE groups on days 13–16 of lactation.</footnote>

<footnote>c P values <0.05 calculated using the Mann-Whitney test to compare C milk at L2-L4 and L13-16 days of lactation.</footnote>

<footnote>d P values <0.05 calculated using the Mann-Whitney test to compare DE-exposed milk at L2-L4 and L13-16 days of lactation.</footnote>
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Leptin concentrations in milk from C and DE rabbits on Days 2–4 and 13–16 of lactation. Measurement of leptin levels in milk from control (C) and DE exposed (DE) dams on Days 2 to 4 (L2-L4; 9 C and 9 DE rabbits) and Days 13 to 16 (L13-L16; 9 C and 9 DE rabbits) of lactation. Data are mean ± SEM.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0212132.g005

mammary epithelium, suggesting an alteration to mammary lipid metabolism induced by DE. Moreover, when the enzymes involved in fatty acid metabolism were analysed, only the Stearoyl-CoA desaturase (SCD) gene was elevated in DE animals. Stearoyl-CoA desaturase is an enzyme that plays a key role in the composition of fatty acids by generating monounsaturated fatty acids (MUFAs) that provide the principal substrates for the biosynthesis of polyunsaturated fatty acids (PUFAs) and complex lipids such as triglycerides, phospholipids and cholesterol esters involved in energy storage, biological membranes' composition and signalling molecules [43, 44]. This elevation of Stearoyl-CoA desaturase levels may therefore explain the rise in fat globules observed in the DE-exposed alveolar lumina as well as increased fatty acid content in DE milk during early lactation. During mid-lactation (L13-L16), SFA levels, mainly C10-C12, were higher in DE milk, while PUFA levels were lower, thus suggesting a rebound in the maturation process of the milk. Moreover, since offspring's weight at birth was decreased in DE rabbits [14], these changes in fatty acid milk composition favouring scSFA, may contribute to the catch-up growth of the pups described in Valentino et al [14].

Unlike data obtained in the placenta, we were not able to detect any particles in mammary tissue. In the placenta, non-aggregated and “fingerprint” nanoparticles have been observed at various locations: in the maternal blood space, in trophoblastic cells and in foetal blood, thus demonstrating transplacental transfer [14]. The analysis of mammary sections using electron microscopy did not reveal any accumulation of nanoparticles.

Mammary gland development is a complex process consisting of a series of precise events that are finely regulated by a balance of hormones, growth factors and stromal factors [45]. The mammary gland achieves much of its development during postnatal life, particularly during pregnancy when the lobulo-alveolar structures branch and develop in preparation for lactation. Pregnancy is therefore critical to mammary development. Such periods are highly sensitive to environmental changes that can alter mammary gland structure and function. Most studies to date have been performed using endocrine disrupters and demonstrate that gestational exposure to pollutants such as dioxin, BPA or genistein alter mammary epithelial development [46–48]. Diesel engine exhaust fumes modify the carcinogenic sensitivity of mammary epithelial cells by affecting metabolic activation pathways [49]. The mammary phenotype observed here after DE exposure shares features with a previously described altered
mammary phenotype that was induced by an obesogenic diet consumed by rabbits from puberty to pregnancy [30, 31]. During both experiments, the mammary lumina exhibited an increase in their fat globule content. Taking account of the rise in SCD expression, we can therefore hypothesise that exposure to DE and an obesogenic diet may trigger common mechanisms [50]. Human studies have already demonstrated that traffic pollution influences metabolism and is positively associated with elevated BMI rates in children [50, 51].

Although lactation performance was not altered in DE-exposed rabbits at Days 2-4 and 13-16, their milk composition displayed major changes associated with the maternal inhalation of fine particles. Milk from the DE dams exhibited an elevated total fatty acid content resulting from high levels of SFA, MUFA and PUFA. This finding was most likely associated with the *de novo* synthesis of fatty acids in the lactating rabbit mammary gland and might be linked to the increase in mammary SCD expression observed in DE rabbits [52]. The milk lipid composition has been shown to be modified following the ingestion of other pollutants such as bisphenol A [53]. However, few studies have reported on the effects of respiratory pollution apart from an analysis of the consequences of maternal passive smoking in a rodent model, which caused a rise in milk lipid concentrations [54].

In order to take account of any putative effects of diet-associated changes to milk composition we also quantified the main milk proteins, including their predominant phosphorylated isoforms. Although no variations in the total casein ratio were observed in DE milk, a significant increase in αs2-like-casein and a significant reduction in αs1-casein were determined. Moreover, the relative proportions of less phosphorylated isoforms for αs2-like-casein were higher in DE milk. As for total whey proteins, the proportions of both lactoferrin and serum albumin were higher in DE milk. Serum albumin and lactoferrin have antimicrobial activities related to their good binding capacity for metals. They also bind toxic substances, fatty acids and lipopolysaccarides [55, 56]. Therefore, the increase of their concentrations in DE milk during early lactation (L2-L4) may constitute a response to the stress induced by DE inhalation. Moreover, the increase of serum albumin and lactoferrin levels may be related to that of lipid metabolism observed in DE milk.

Altogether, our experiments revealed differences in fatty acid and protein milk composition during early lactation (L2-L4), but these levels were similar when mature milk (L13-L16) was examined. Maternal DE exposure throughout pregnancy involves potential risks of growth and metabolic disorders in offspring by affecting foetal and placental development [14]. The consumption of milk during the neonatal period is a characteristic of mammals and provides newborns with essential nutritive support during the transition from neonatal to later life. Milk plays an important role in regulating metabolic programming so that milk composition may have long-term consequences for the health and development of offspring [36, 57]. But although questions remain as to the mechanisms involved in these outcomes, knowledge of early life programming has made significant progress in describing the offspring phenotypes that result from modifications to the neonatal diet. In particular, post-natal brain development is tightly dependent on an optimum fatty acid composition [58]. Moreover, a number of studies have described the very early post-natal period (first ten days) as a critical temporal window that will regulate subsequent food intake in adulthood. Leptin is considered as a key factor in this process [59, 60]. Surprisingly, our study did not reveal any difference in milk leptin levels between the C and DE groups, particularly in milk at early lactation (days L2-L4), a period during which most of the changes to milk composition were observed in our experiments.

It is becoming increasingly clear that developmental exposure to environmental contaminants can have lasting effects on metabolic tissues and have endocrine disrupting potential [61]. A number of studies have proposed the existence of lactational programming [62]. In almost all of them, the effects of pollutants on breast milk composition focused on the
ingestion of chemicals such as organochlorine pesticides, bisphenol A, dioxins or heavy metals [46, 53, 63–65]. Few data are available on airway pollution or fine particle inhalation, except in the context of passive smoking [54]. Our experiments are the first to have suggested a link between maternal DE exposure and changes to mammary gland function and subsequent milk composition. It is now of crucial importance to decipher the mechanistic and functional results of DE exposure in mammary epithelial cells, particularly at the level of metabolic pathways. Furthermore, from a public health point of view, such DE exposure-induced changes to milk composition should be taken into account because they may reprogram the developmental trajectory of organs and affect the adult life of offspring.

**Supporting information**

S1 Fig. Reversed-phase HPLC UV 214 nm chromatograms of rabbit control (C) and DE skimmed milk samples on days 2–4 (L2-L4) and 13–16 (L13-L16) of lactation. Identification of the major milk protein peacks 1: κ-casein, 2: Lactoferrin, 3–4: α<sub>s2</sub>-casein and WAP, 5-6-7: α-lactalbumin and Serum albumin, 8-9-10: α<sub>s2</sub>-like-casein, 11-12-13: α<sub>s1</sub>-casein, 14: β-casein. (TIF)

**Acknowledgments**

PCP is member of COST Actions FA0702 GEMINI "Maternal interaction with gametes and embryo", FA1201 EPICONCEPT "Epigenetics and periconceptional environment" and BM1308 "Sharing advances on large animal models (SALAAM)". The authors thank the UCEA (Unité Commune d'Expérimentation Animale), the INRA team of BDR (Biologie du Développement et Reproduction) for housing the rabbits and participating to the experiments. Christophe Richard (INRA BDR) designed the nose only inhalation procedure, Dr. Flemming Cassee (RIVM, National Institute for Public Health and the Environment Netherlands) provided expertise for the inhalation procedure. John Boere and Paul Fokkens monitored the DE exposure using the MAPCEL device.

The authors are also grateful to Dr. Pierre Adenot for his assistance with Image J analyses. This work greatly benefited from the facilities and expertise of the MIMA2/ Platform @BRIDGE, INRA, AgroParisTech, Université Paris-Saclay, Jouy-en-Josas, France.

**Author Contributions**

**Conceptualization:** Sarah Valentino, Pascale Chavatte-Palmer, Madia Charlier.

**Formal analysis:** Cathy Hue-Beauvais, Etienne Aujean, Guy Miranda, Delphine Ralliard-Rousseau, Sarah Valentino, Christine Péchoux, Madia Charlier.

**Funding acquisition:** Pascale Chavatte-Palmer.

**Investigation:** Cathy Hue-Beauvais, Sarah Valentino, Madia Charlier.

**Methodology:** Cathy Hue-Beauvais, Etienne Aujean, Guy Miranda, Delphine Ralliard-Rousseau, Nicolas Brun, Stessy Ladebese, Christine Péchoux, Madia Charlier.

**Supervision:** Pascale Chavatte-Palmer, Madia Charlier.

**Validation:** Cathy Hue-Beauvais, Delphine Ralliard-Rousseau, Pascale Chavatte-Palmer, Madia Charlier.

**Writing – original draft:** Cathy Hue-Beauvais, Guy Miranda, Delphine Ralliard-Rousseau, Madia Charlier.
## References

1. Dockery DW, Pope CA 3rd, Xu X, Spengler JD, Ware JH, Fay ME, et al. An association between air pollution and mortality in six U.S. cities. The New England journal of medicine. 1993; 329(24):1753-9. Epub 1993/12/09. https://doi.org/10.1056/NEJM199312093292401 PMID: 8179653.

2. Samet JM, Zeger SL, Dominici F, Curriero F, Coursac I, Dockery DW, et al. The National Morbidity, Mortality, and Air Pollution Study. Part II: Morbidity and mortality from air pollution in the United States. Research report (Health Effects Institute). 2000; 94(Pt 2):5-70; discussion 1-9. Epub 2001/05/17. PMID: 11354823.

3. Chang ET, Lau EC, Van Landingham C, Crump KS, McClellan RO, Moolgavkar SH. Reanalysis of Diesel Engine Exhaust and Lung Cancer Mortality in the Diesel Exhaust in Miners Study (DEMS) Cohort Using Alternative Exposure Estimates and Radon Adjustment. American journal of epidemiology. 2018. Epub 2018/03/10. https://doi.org/10.1093/aje/kwy038 PMID: 29522073.

4. Gregory DJ, Kobzik L, Yang Z, McGuire CC, Fedulov AV. Transgenerational transmission of asthma risk after exposure to environmental particles during pregnancy. American journal of physiology Lung cellular and molecular physiology. 2017; 313(2):L395-L405. Epub 2017/05/13. https://doi.org/10.1152/ajplung.00035.2017 PMID: 28495853; PubMed Central PMCID: PMCPMC5582941.

5. Li C, Li X, Suzuki AK, Zhang Y, Fujitani Y, Nagaoka K, et al. Effects of exposure to nanoparticle-rich diesel exhaust on pregnancy in rats. The Journal of reproduction and development. 2013; 59(2):145-50. Epub 2012/12/22. https://doi.org/10.1262/jrd.2012-145 PMID: 23257834; PubMed Central PMCID: PMCPMC3934203.

6. Frutos V, Gonzalez-Comadran M, Sola I, Jacquemin B, Carreras R, Checa Vizcaino MA. Impact of air pollution on fertility: a systematic review. Gynecological endocrinology: the official journal of the International Society of Gynecological Endocrinology. 2015; 31(1):7-13. Epub 2014/09/13. https://doi.org/10.3109/09513590.2014.958992 PMID: 25212280.

7. Kowalska M, Wegierek-Ciuk A, Brzoska K, Wojewodzka M, Meczynska-Wielgosz S, Gromadzka-Ostrowska J, et al. Genotoxic potential of diesel exhaust particles from the combustion of first- and second-generation biodiesel fuels-the FuelHealth project. Environmental science and pollution research international. 2017; 24(31):24223-34. Epub 2017/09/11. https://doi.org/10.1007/s11356-017-9995-0 PMID: 28889235; PubMed Central PMCID: PMCPMC5655577.

8. Burgaz S, Demircigil GC, Karahalil B, Karakaya AE. Chromosomal damage in peripheral blood lymphocytes of traffic policemen and taxi drivers exposed to urban air pollution. Chemosphere. 2002; 47 (1):57-64. Epub 2002/05/09. PMID: 11996136.

9. Pedersen M, Giorgis-Allemand L, Bernard C, Aguilera I, Andersen AM, Ballester F, et al. Ambient air pollution and low birthweight: a European cohort study (ESCAPE). The Lancet Respiratory medicine. 2013; 1(9):695-704. Epub 2014/01/17. https://doi.org/10.1016/S2213-2600(13)70192-9 PMID: 24429273.

10. Howlin J, McBryan J, Martin F. Pubertal mammary gland development: insights from mouse models. Journal of mammary gland biology and neoplasia. 2006; 11(3-4):283-97. Epub 2006/11/08. https://doi.org/10.1007/s10911-006-9024-2 PMID: 17089203.

11. Hennighausen L, Robinson GW. Information networks in the mammary gland. Nature reviews Molecular cell biology. 2005; 6(9):715-25. Epub 2005/10/19. https://doi.org/10.1038/nrm1714 PMID: 16231422.

12. Barker DJ. The developmental origins of adult disease. Journal of the American College of Nutrition. 2004; 23(6 Suppl):588S-95S. Epub 2005/01/11. PMID: 15640511.

13. Grandjean P, Barouki R, Bellinger DC, Castelein L, Chadwick LH, Cordier S, et al. Life-Long Implications of Developmental Exposure to Environmental Stressors: New Perspectives. Endocrinology. 2015; 156(10):3408-15. Epub 2015/08/05. https://doi.org/10.1210/EN.2015-1350 PMID: 26241067; PubMed Central PMCID: PMCPMC4588822.

14. Valentino SA, Tarrade A, Aioun J, Mourier E, Richard C, Dahirel M, et al. Maternal exposure to diluted diesel engine exhaust alters placental function and induces intergenerational effects in rabbits. Particle and fibre toxicology. 2016; 13(1):39. Epub 2016/07/28. https://doi.org/10.1186/s12989-016-0151-7 PMID: 27460165; PubMed Central PMCID: PMCPMC4962477.

15. Fleisch AF, Rifas-Shiman SL, Koutrakis P, Schwartz JD, Kloog I, Melly S, et al. Prenatal exposure to traffic pollution: associations with reduced fetal growth and rapid infant weight gain. Epidemiology (Cambridge, Mass). 2015; 26(1):43-50. Epub 2014/12/02. https://doi.org/10.1097/ede. 0000000000000203 PMID: 25437317; PubMed Central PMCID: PMCPMC4285344.

16. Veras MM, Damaceno-Rodrigues NR, Guimaraes Silva RM, Scoriza JN, Saldiva PH, Caldini EG, et al. Chronic exposure to fine particulate matter emitted by traffic affects reproductive and fetal outcomes in mice. Environmental research. 2009; 109(5):536-43. Epub 2009/04/28. https://doi.org/10.1016/j.envres.2009.03.006 PMID: 19394924.
17. Dadvand P, Parker J, Bell ML, Bonzini M, Brauer M, Darrow LA, et al. Maternal exposure to particulate air pollution and term birth weight: a multi-country evaluation of effect and heterogeneity. Environmental health perspectives. 2013; 121(3):267-373. Epub 2013/02/07. https://doi.org/10.1289/ehp.1205575 PMID: 23384584; PubMed Central PMCID: PMCPMC3621183.

18. Chen M, Liang S, Zhou H, Xu Y, Qin X, Hu Z, et al. Prenatal and postnatal mothering by diesel exhaust PM2.5-exposed dams differentially program mouse energy metabolism. Particle and fibre toxicology. 2017; 14(1):3. Epub 2017/01/20. https://doi.org/10.1186/s12989-017-0183-7 PMID: 28100227; PubMed Central PMCID: PMCPMC5423412.

19. Bolton JL, Auten RL, Bilbo SD. Prenatal air pollution exposure induces sexually dimorphic fetal programming of metabolic and neuroinflammatory outcomes in adult offspring. Brain, behavior, and immunity. 2014; 37:30-44. Epub 2013/11/05. https://doi.org/10.1016/j.bbi.2013.10.029 PMID: 24184474.

20. Weldy CS, Liu Y, Liggitt HD, Chin MT. 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. 2014; 9(2):e88582. Epub 2014/02/18. https://doi.org/10.1371/journal.pone.0088582 PubMed Central PMCID: PMCPMC3922927. PMID: 24533117

21. Gorr MW, Velten M, Nelin TD, Youtz DJ, Sun Q, Wold LE. Early life exposure to air pollution induces adult cardiac dysfunction. American journal of physiology Heart and circulatory physiology. 2014; 307 (9):H1353-60. Epub 2014/08/31. https://doi.org/10.1152/ajpheart.00526.2014 PMID: 25172901; PubMed Central PMCID: PMCPMC4217011.

22. Lu R, Yuan T, Wang Y, Zhang T, Yuan Y, Wu D, et al. Spontaneous severe hypercholesterolemia and atherosclerosis lesions in rabbits with deficiency of low-density lipoprotein receptor (LDLR) on exon 7. EBioMedicine. 2018. Epub 2018/09/24. https://doi.org/10.1016/j.ebiom.2018.09.020 PMID: 30243490.

23. Fischer B, Chavatte-Palmer P, Viebahn C, Navarrete Santos A, Duranthon V. Rabbit as a reproductive model for human health. Reproduction (Cambridge, England). 2012; 144(1):1-10. Epub 2012/05/15. https://doi.org/10.1530/rep-12-0091 PMID: 22580370.

24. Graur D, Duret L, Gouy M. Phylogenetic position of the order Lagomorpha (rabbits, hares and allies). Nature. 1996; 379(6563):333-5. Epub 1996/01/25. https://doi.org/10.1038/379333a0 PMID: 8552186.

25. Martin P, Cebo C, Miranda G. Inter-species comparison of milk proteins: quantitative variability and molecular diversity. Advanced Dairy Chemistry: Proteins. 4th Edition. 4th Edition ed: Fox P.F. & McSweeney P.L.F, Kluwer Academic/Plenum Publishers, New York; 2012.

26. Maertens L, Lebas F, Szendro Z. Rabbit Milk: A review of quantity, quality and non-dietary affecting factor. World Rabbit Sciences. 2006; 14:205-30.

27. Reynolds C, Gray C, Li M, Segovia S, Vickers M. Early Life Nutrition and Energy Balance Disorders in Offspring in Later Life. Nutrients. 2015; 7(9):8090-111. Epub 2015/09/25. https://doi.org/10.3390/nu7095384 PMID: 26402696; PubMed Central PMCID: PMCPMC4586579.

28. Freney EJ, Heal MR, Donovan RJ, Mills NL, Donaldson K, Newby DE, et al. A single-particle characterization of a mobile Versatile Aerosol Concentration Enrichment System for exposure studies. Particle and fibre toxicology. 2006; 3:8. Epub 2006/05/26. https://doi.org/10.1186/1743-8977-3-8 PMID: 16723024; PubMed Central PMCID: PMCPMC1488865.

29. U.S. Environmental Protection Agency OoRaD, National Center for Environmental Assessment. U.S. EPA. Health Assessment Document for Diesel Engine Exhaust (Final 2002). 2002.

30. Hue-Beauvais C, Koch E, Chavatte-Palmer P, Galio L, Chat S, Letheule M, et al. Milk from dams fed an obesogenic diet combined with a high-fat/high-sugar diet induces long-term abnormal mammary gland development in the rabbit. J Anim Sci. 2015; 93(4):1641-55. Epub 2015/05/29. https://doi.org/10.2527/jas.2014-8139 PMID: 26020186.

31. Hue-Beauvais C, Chavatte-Palmer P, Aujean E, Dahirel M, Laigre P, Pechoux C, et al. An obesogenic diet started before puberty leads to abnormal mammary gland development during pregnancy in the rabbit. Dev Dyn. 2011; 240(2):347-56. https://doi.org/10.1002/dvdy.22536 PMID: 21246651.

32. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (San Diego, Calif). 2001; 25(4):402-8. Epub 2002/02/16. https://doi.org/10.1006/meth.2001.1262 PMID: 11846609.

33. Folch J, Lees M, Sloane Stanley GH. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957; 226(1):497-509. Epub 1957/05/01. PMID: 13428781.

34. Morrison WR, Smith LM. Preparation of fatty acid methyl esters and dimethylacetals from lipids with boron fluoride methanol. Journal of lipid research. 1964; 5:600-8. Epub 1964/10/01. PMID: 14221106.

35. Rousseau D, Helies-Toussaint C, Moreau D, Raederstorff D, Grynberg A. Dietary n-3 PUFAs affect the blood pressure rise and cardiac impairments in a hyperinsulinemia rat model in vivo. American journal of physiology Heart and circulatory physiology. 2003; 285(3):H1294-302. Epub 2003/01/11. https://doi.org/10.1152/ajpheart.00651.2002 PMID: 12521943.
36. Hue-Beauvais C, Miranda G, Aujean E, Jaffrezic F, Devinoy E, Martin P, et al. Diet-induced modifications to milk composition have long-term effects on offspring growth in rabbits. J Anim Sci. 2017; 95 (2):761-70. Epub 2017/04/06. https://doi.org/10.2527/jas.2016.0847 PMID: 28380600.

37. Denamur R. Nucleic acids of the mammary gland during gestation and lactation in the rabbit. Comptes rendus hebdomadaires des seances de l'Academie des sciences. 1963; 256:4748-50. Epub 1963/05/27. PMID: 14026945.

38. Koch E, Hue-Beauvais C, Galio L, Solomon G, Gertler A, Revillon F, et al. Leptin gene in rabbit: cloning and expression in mammary epithelial cells during pregnancy and lactation. Physiological genomics. 2013. Epub 2013/05/30. https://doi.org/10.1152/physiolgenomics.00020.2013 PMID: 23715260.

39. Aoki N, Kawamura M, Matsuda T. Lactation-dependent down regulation of leptin production in mouse mammary gland. Biochim Biophys Acta. 1999; 1427(2):298-306. PMID: 10216246

40. Bartha T, Sayed-Ahmed A, Rudas P. Expression of leptin and its receptors in various tissues of ruminants. Domestic animal endocrinology. 2005; 29(1):193-202. Epub 2005/05/10. https://doi.org/10.1016/j.domaniend.2005.03.010 PMID: 15878255.

41. Bonnet M, Delavaud C, Laud K, Gourdou I, Leroux C, Djiane J, et al. Mammary leptin synthesis, milk leptin and their putative physiological roles. Reproduction, nutrition, development. 2002; 42(5):399-413. Epub 2003/01/23. PMID: 12537253.

42. Hu X, Juneja SC, Maihle NJ, Cleary MP. Leptin—a growth factor in normal and malignant breast cells and for normal mammary gland development. Journal of the National Cancer Institute. 2002; 94 (22):1704-11. Epub 2002/11/21. PMID: 12441326.

43. Shorten PR, Pleasants TB, Upreti GC. A mathematical model for mammary fatty acid synthesis and triglyceride assembly: the role of stearoyl CoA desaturase (SCD). The Journal of dairy research. 2004; 71 (4):385-97. Epub 2004/12/21. PMID: 15605704.

44. Wu X, Zou X, Chang Q, Zhang Y, Li Y, Zhang L, et al. The evolutionary pattern and the regulation of stearoyl-CoA desaturase genes. BioMed research international. 2013; 2013:856521. Epub 2013/12/07. https://doi.org/10.1155/2013/856521 PMID: 24312911; PubMed Central PMCID: PMCPMC3838806.

45. Hennighausen L. Molecular mechanisms of hormone controlled gene expression in the breast. Molecular biology reports. 1997; 24(3):169-74. Epub 1997/08/01. PMID: 9291090.

46. Kass L, Altamirano GA, Bosquiazzo VL, Luque EH, Munoz-de-Toro M. Perinatal exposure to xenoestrogens impairs mammary gland differentiation and modifies milk composition in Wistar rats. Reproductive toxicology (Elmsford, NY). 2012; 33(3):390-400. Epub 2012/02/22. https://doi.org/10.1016/j.reprotox.2012.02.002 PMID: 22349186.

47. Lew BJ, Collins LL, O'Reilly MA, Lawrence BP. Activation of the aryl hydrocarbon receptor during different critical windows in pregnancy alters mammary epithelial cell proliferation and differentiation. Toxicological sciences: an official journal of the Society of Toxicology. 2009; 111(1):151-62. Epub 2009/06/09. https://doi.org/10.1093/toxsci/kfp125 PMID: 19502548; PubMed Central PMCID: PMCPMC2726296.

48. El Sheikh Saad H, Meduri G, Phrakonkham P, Berges R, Vacher S, Djallali M, et al. Abnormal peripubertal development of the rat mammary gland following exposure in utero and during lactation to a mixture of genistein and the food contaminant vinclozolin. Reproductive toxicology (Elmsford, NY). 2011; 32(1):15-25. Epub 2011/05/05. https://doi.org/10.1016/j.reprotox.2011.03.001 PMID: 21539910.

49. Courter LA, Pereira C, Baird WM. Diesel exhaust influences carcinogenic PAH-induced genotoxicity and gene expression in human breast epithelial cells in culture. Mutation research. 2007; 625(1-2):72-82. Epub 2007/07/07. https://doi.org/10.1016/j.mrfmmm.2007.05.005 PMID: 17612574; PubMed Central PMCID: PMCPMC2815343.

50. Jerrett M, McConnell R, Wolch J, Chang R, Lam C, Dunton G, et al. Traffic-related air pollution and obesity formation in children: a longitudinal, multilevel analysis. Environmental health: a global access science source. 2014; 13:49. Epub 2014/06/11. https://doi.org/10.1186/1476-069x-13-49 PMID: 24913018; PubMed Central PMCID: PMCPMC4106205.

51. Bolton JL, Smith SH, Huff NC, Gilmour MI, Foster WM, Auten RL, et al. Prenatal air pollution exposure induces neuroinflammation and predisposes offspring to weight gain in adulthood in a sex-specific manner. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2012; 26(11):4743-54. Epub 2012/07/21. https://doi.org/10.1096/fj.12-210989 PMID: 22815382.

52. Strong CR, Dils R. Fatty acid biosynthesis in rabbit mammary gland during pregnancy and early lactation. The Biochemical journal. 1972; 128(5):1303-9. Epub 1972/08/01. PMID: 4643704; PubMed Central PMCID: PMCPMC1174018.

53. Altamirano GA, Munoz-de-Toro M, Luque EH, Gomez AL, Delconte MB, Kass L. Milk lipid composition is modified by perinatal exposure to bisphenol A. Molecular and cellular endocrinology. 2015; 411:258-67. Epub 2015/05/16. https://doi.org/10.1016/j.mce.2015.05.007 PMID: 25976663.
54. Baheiraei A, Shamsi A, Khaghani S, Shams S, Chamari M, Boushehri H, et al. The effects of maternal passive smoking on maternal milk lipid. Acta medica Iranica. 2014; 52(4):280-5. Epub 2014/06/06. PMID: 24901858.

55. Topala T, Bodoki A, Oprean L, Oprean R. Bovine Serum Albumin Interactions with Metal Complexes. Clujul medical (1957). 2014; 87(4):215-9. Epub 2014/01/01. https://doi.org/10.15386/cjmed-357 PMID: 26528027; PubMed Central PMCID: PMCPMC4620676.

56. Telang S. Lactoferrin: A Critical Player in Neonatal Host Defense. Nutrients. 2018; 10(9). Epub 2018/09/06. https://doi.org/10.3390/nu10091228 PMID: 30181493.

57. Neville MC, Anderson SM, McManaman JL, Badger TM, Bunik M, Contractor N, et al. Lactation and neonatal nutrition: defining and refining the critical questions. Journal of mammary gland biology and neoplasia. 2012; 17(2):167-88. https://doi.org/10.1007/s10911-012-9261-5 PMID: 22752723; PubMed Central PMCID: PMC3428522.

58. Morale SE, Hoffman DR, Castaneda YS, Wheaton DH, Burns RA, Birch EE. Duration of long-chain polyunsaturated fatty acids availability in the diet and visual acuity. Early Hum Dev. 2005; 81(2):197-203. Epub 2005/03/08. https://doi.org/10.1016/j.earlhumdev.2004.09.008 PMID: 15748975.

59. Bouret SG. Neurodevelopmental actions of leptin. Brain research. 2010; 1350:2-9. Epub 2010/04/20. https://doi.org/10.1016/j.brainres.2010.04.011 PMID: 20399755; PubMed Central PMCID: PMCPMC3654158.

60. Attig L, Brisard D, Larcher T, Mickiewicz M, Guilloteau P, Boukthir S, et al. Postnatal Leptin Promotes Organ Maturation and Development in IUGR Piglets. PloS one. 2013; 8(5):e64616. Epub 2013/06/07. https://doi.org/10.1371/journal.pone.0064616 PMID: 23741353; PubMed Central PMCID: PMCPMC3669417.

61. Heindel JJ, Blumberg B, Cave M, Machtinger R, Mantovani A, Mendez MA, et al. Metabolism disrupting chemicals and metabolic disorders. Reproductive toxicology (Elmsford, NY). 2017; 68:3-33. Epub 2016/10/21. https://doi.org/10.1016/j.reprotox.2016.10.001 PMID: 27760374; PubMed Central PMCID: PMCPMC5365353.

62. Ellsworth L, Harman E, Padmanabhan V, Gregg B. Lactational programming of glucose homeostasis: a window of opportunity. Reproduction (Cambridge, England). 2018; 156(2):R23-R42. Epub 2018/05/13. https://doi.org/10.1530/rep-17-0780 PMID: 29752297.

63. Solomon GM, Weiss PM. Chemical contaminants in breast milk: time trends and regional variability. Environmental health perspectives. 2002; 110(6):A339-47. Epub 2002/06/11. https://doi.org/10.1289/ehp.021100339 PMID: 12055065; PubMed Central PMCID: PMCPMC1240888.

64. van Esterik JC, Dolle ME, Lamoree MH, van Leeuwen SP, Hamers T, Legler J, et al. Programming of metabolic effects in C57BL/6JxFVB mice by exposure to bisphenol A during gestation and lactation. Toxicology. 2014; 321:40-52. Epub 2014/04/15. https://doi.org/10.1016/j.tox.2014.04.001 PMID: 24726836.

65. Vahidinia A, Samiee F, Faradmal J, Rahmani A, Taravati Javad M, Leili M. Mercury, Lead, Cadmium, and Barium Levels in Human Breast Milk and Factors Affecting Their Concentrations in Hamadan, Iran. Biological trace element research. 2018. Epub 2018/04/28. https://doi.org/10.1007/s12011-018-1355-5 PMID: 29700700.