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

HAL
open science

# Maternal nutrition during pregnancy affects testicular and bone development, glucose metabolism and response to overnutrition in weaned horses up to two years

Morgane Robles, Camille Gautier, Luis Mendoza, Pauline Peugnet, Cédric

Dubois, Michèle Dahirel, Jean-Philippe Lejeune, Isabelle Caudron, Isabelle

Guenon, Sylvaine Camous, et al.

► To cite this version:

Morgane Robles, Camille Gautier, Luis Mendoza, Pauline Peugnet, Cédric Dubois, et al. Maternal nutrition during pregnancy affects testicular and bone development, glucose metabolism and response to overnutrition in weaned horses up to two years. PLoS ONE, 2017, 12 (1), pp.1-21. 10.1371/journal.pone.0169295. [hal-01607465](https://hal.science/hal-01607465)

HAL Id: hal-01607465

https://hal.science/hal-01607465

Submitted on 25 May 2020

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Maternal Nutrition during Pregnancy Affects Testicular and Bone Development, Glucose Metabolism and Response to Overnutrition in Weaned Horses Up to Two Years

Morgane Robles¹*, Camille Gautier²*, Luis Mendoza³, Pauline Peugnet¹, Cédric Dubois⁴, Michèle Dahirel¹, Jean-Philippe Lejeune³, Isabelle Caudron³, Isabelle Guenon², Sylvaine Camous¹, Anne Tarrade¹, Laurence Wimel⁴, Didier Serteyn³, Hélène Bouraima-Lelong², Pascale Chavatte-Palmer¹*

1 UMR BDR, INRA, ENVA, Université Paris Saclay, Jouy en Josas, France, 2 Normandie Univ, UNICAEN, EA2608, OeReCa, USC-INRA, Caen, France, 3 Clinique Equine, Faculté de Médecine Vétérinaire, Université de Liège, Liège, Belgium, 4 IFCE, Station Expérimentale de la Valade, Chamberet, France

* These authors contributed equally to this work.

* pascale.chavatte@jouy.inra.fr

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

## Abstract

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

### OPEN ACCESS

**Citation:** Robles M, Gautier C, Mendoza L, Peugnet P, Dubois C, Dahirel M, et al. (2017) Maternal Nutrition during Pregnancy Affects Testicular and Bone Development, Glucose Metabolism and Response to Overnutrition in Weaned Horses Up to Two Years. PLoS ONE 12(1): e0169295.
doi:10.1371/journal.pone.0169295

**Editor:** Clemens Fürnsinn, Medical University of Vienna, AUSTRIA

**Received:** September 9, 2016

**Accepted:** December 14, 2016

**Published:** January 12, 2017

**Copyright:** © 2017 Robles et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The Dryad URL of our data is: http://dx.doi.org/10.5061/dryad.76dn2. The DOI is: doi:10.5061/dryad.76dn2.

**Funding:** This work was funded by grants from the "Institut Français du Cheval et de l'Equitation" and by the "Fonds Eperon" (FOETALIM project). The funders had no role in study design, data collection and analysis, decision to publish, nor preparation of the manuscript.

## Introduction

Pregnant mares and post-weaning foals are often fed concentrates rich in soluble carbohydrates, together with forage. Recent studies suggest that the use of concentrates is linked to alterations of metabolism and the development of osteochondrosis in foals. The aim of this study was to determine if broodmare diet during gestation affects metabolism, osteoarticular status and growth of yearlings overfed from 20 to 24 months of age and/or sexual maturity in prepubertal colts.

## Material and methods

Twenty-four saddlebred mares were fed forage only (n = 12, group F) or cracked barley and forage (n = 12, group B) from mid-gestation until foaling. Colts were gelded at 12 months of age. Between 20 and 24 months of age, all yearlings were overfed (+140% of requirements) using an automatic concentrate feeder. Offspring were monitored for growth between 6 and 24 months of age, glucose homeostasis was evaluated *via* modified frequently sampled intra venous glucose tolerance test (FSIGT) at 19 and 24 months of age and osteoarticular status was investigated using radiographic examinations at 24 months of age. The structure and function of testicles from prepubertal colts were analyzed using stereology and RT-qPCR.

## Results

Post-weaning weight growth was not different between groups. Testicular maturation was delayed in F colts compared to B colts at 12 months of age. From 19 months of age, the cannon bone was wider in B vs F yearlings. F yearlings were more insulin resistant at 19 months
**Competing Interests:** The authors have declared that no competing interests exist.

compared to B yearlings but B yearlings were affected more severely by overnutrition with reduced insulin sensitivity. The osteoarticular status at 24 months of age was not different between groups.

## Conclusion

In conclusion, nutritional management of the pregnant broodmare and the growing foal may affect sexual maturity of colts and the metabolism of foals until 24 months of age. These effects may be deleterious for reproductive and sportive performances in older horses.

## Introduction

Since the early nineties, many epidemiological and experimental studies have shown that modifications of the nutritional environment during gestation and in childhood can affect off-spring health at adulthood. The Developmental Origins of Health and Disease (DOHaD) were first demonstrated in human populations [1], with low birthweight neonates or *in utero* under-fed infants during the Dutch Hunger Winter having a higher risk of developing a metabolic syndrome as adults [2,3]. Effects of maternal nutritional imbalance on offspring health were subsequently confirmed using animal models [4].

Studies on DOHaD remain scarce in horses [5]. Embryo transfers between breeds of different size were performed to enhance or reduce fetal growth. Postnatal growth patterns were affected from birth until 3 years of age. Moreover, cardiovascular function, glucose homeostasis and endocrine function were altered until weaning [6-12]. Despite the fact that severe maternal undernutrition at mid-gestation related to an infection with *Streptococcus equi* was shown in one study to induce intra-uterine growth retardation (IUGR) [13], most other studies on effects of maternal nutrition indicate that moderate maternal nutritional excess or undernutrition in gestation does not affect birthweight [14-18]. Nevertheless, excess maternal nutrition during pregnancy, with or without carbohydrate supplementation, can affect glucose and lipid metabolism in foals until 160 days of age [13,15,17]. Moreover, an epidemiological study performed in Belgium in a large cohort of horses demonstrated that the incidence of osteochondrosis (OC) is higher in foals born to dams fed concentrates during gestation compared to foals whose dams received forage only [19].

Although there is limited evidence in humans and other species that maternal nutrition may affect offspring reproductive function [20], to our knowledge, consequences on sperm production and fertility remain unknown [21]. Of particular interest for the present study, maternal gestational undernutrition in rats and sheep affects pituitary responsiveness to Gonadotropin Releasing Hormon (GnRH) and testicular development and delays the onset of puberty [22-25]. Also, excess maternal nutrition in domestic animals has been shown to affect testicular development and to delay the onset of puberty [26-28]. To date, these effects have not been studied in horses.

In the equine industry, it is common practice to feed pregnant mares and post-weaning foals with concentrates in order to enhance the growth of fetuses and then foals that will be sold at 1 or 2 years of age [29,30]. The present study is a follow-up of previous work in which broodmares were supplemented (group B) or not (group F) with concentrates from mid-gestation to term [31]. Briefly, the daily supplementation of B mares with flattened barley (2 to 3 kg) provided additional net energy and digestible proteins from the 7th until the 10th month of gestation compared to F mares, allowing them to maintain optimal body condition throughout
pregnancy and lactation. Glycemic and insulinemic peaks were observed after each meal. F
mares, fed forage only, were energy deficient and lost body condition until parturition. Lipo-
mobilization was increased (as demonstrated by increased Non Esterified Fatty Acids plasma
concentrations (NEFA)) as well as insulin sensitivity during late gestation in F mares, that
were thus considered underfed compared to B mares that had maintained body condition.
From birth until weaning, no difference was observed in foals' body measurements nor endo-
crine function (Insulin-like Growth Factor 1 (IGF-1), triiodothyronine (T3), thyroxine (T4)
and Non Esterified Fatty Acids (NEFA)). F foals, however, tended to have decreased glucose
tolerance and lower osteocalcin concentrations at 3 days of age. At 6 months of age (weaning),
these differences had disappeared and osteoarticular status did not differ between both groups
of foals at that stage [31].

Here we compare growth, glucose metabolism, testicular maturation and joint health in
these foals and yearlings from weaning to 24 months of age. In this period, a 5 month overnu-
trition challenge was performed from 20 to 24 months of age.

**Materials and Methods**

Animals

**Ethical statement.** The animal studies were approved by the local animal care and use committee (“Comité des Utilisateurs de la Station Expérimentale de Chamberet”) and received ethical approval from the local ethics committee (« Comité Régional d’Ethique pour l’Expérimentation Animale du Limousin ») under protocol number 5-2013-5.

**Experimental design, management and feeding of mares and foals.** The experimental design is presented in Fig 1. Twenty-four (n = 2 primiparous and n = 22 multiparous) saddle-bred mares (mean age 10 years; range 5–21 years) were fed forage only (F, n = 12) or cracked barley and forage (B, n = 12) from the 7th month of gestation until foaling. Management of

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

**Fig 1. Experimental design from insemination of the broodmares until foals' 24 months of age.** From insemination until 7 months of gestation, pregnant mares were managed as one same herd in the same pasture. From the 7th month of gestation, they were housed in individual boxes and allocated to one of two groups: fed with forage only (F, n = 12) or with forage and flattened barley (B, n = 12) until foaling. From 3 days until 6 months of age, foals were in pasture with their dams [31]. They were weaned at 6 months and housed in open barns and fed the same amount of feed, following the INRA recommendations for growing foals [29]. Colts were castrated at 12 months as a routine procedure. From 12 to 19 months of age, all yearlings were managed in the same pasture. At 19 months of age, all yearlings were housed in open barns and between 20 and 24 months of age, they were over nourished using an automatic feeder (+140% of INRA energy requirements, +135% of NRC energy requirements hatched part).

doi:10.1371/journal.pone.0169295.g001
<table><caption>Table 1. Daily nutritional supply given to yearlings during feeding challenge (19 to 24 months).</caption><thead><tr><th>Median age (days)</th><th>Median age (months)</th><th>Type of Haylage</th><th>Hay (kg)</th><th>Haylage (kg)</th><th>Homemade pellets (kg)</th></tr></thead><tbody><tr><td>576</td><td>19</td><td>H1</td><td>4</td><td>2.5</td><td>2</td></tr><tr><td>583</td><td>20</td><td>H1</td><td>4</td><td>2.5</td><td>4</td></tr><tr><td>594</td><td>20</td><td>H1</td><td>4</td><td>2.5</td><td>6</td></tr><tr><td>698</td><td>23</td><td>H2</td><td>4</td><td>2.5</td><td>6</td></tr></tbody></table>

H1 and H2 indicate when a different hay batch was used. N = 24.

doi:10.1371/journal.pone.0169295.t001

mares and foals from insemination to weaning has been described previously by Peugnet *et al.* [31].

In the present paper, colts born to primiparous dams (1 colt from a F dam and 1 colt from a B dam from the same protocol) have been added in the results analysis.

In their first winter, foals were housed in open barns and fed hay and concentrates in agreement with current INRA recommendations for growing foals [31]. Six months after weaning (April 25th, at 11 months of age), yearlings returned to pasture and were managed in one same herd with free access to water and minerals salts. On November 25th, at 19 months of age, all animals returned to open barns. Animals were separated in 3 groups sorted by bodyweight in 3 different stalls to reduce competition for food (housing group). They were group-fed hay and haylage twice a day and supplemented individually using an automatic concentrate feeder with homemade pellets containing barley, oat, soybean cake, bran, molasses and vitamins and minerals with free access to water. The automatic concentrate feeder allowed yearlings to ingest concentrates at least 4 times a day. The diet was adjusted to meet the requirements during 18 days of adaptation period [29].

From December 13th, at 20 months of age, all yearlings were fed to reach an average daily weight gain of 500 g for the upcoming 4 months. Thus, during this period, yearlings were fed 140% of INRA requirements [29] and 135% of the NRC requirements [32], with 62% of the energy of the diet provided by the pellets and 38% provided by forages. A different haylage batch was used at the end of the feeding challenge because the batch used as the start was out of stock. Table 1 presents the daily nutritional supply, Table 2 the quality of feedstuff (INRA feeding system) expressed as net energy (Horse Feed Unit-HFU-; 1 HFU = 2250 kcal), proteins (Horse Digestible Crude Proteins, HDCP), fibers (Raw Cellulose, RC), calcium and phosphorus and Table 3 the quality of feedstuff (NRC feeding system) expressed as digestible energy (estimated using the percentage of difference between net energy and digestible energy for the feedstuff, Mcal) and crude proteins (analyzed) during this period. Concentrate intake

<table><caption>Table 2. quality of feedstuff from INRA [29] given to yearlings during feeding challenge (19 to 24 months).</caption><thead><tr><th rowspan="2"></th><th colspan="4">Chemical composition (per kg of dry matter)</th><th colspan="2">Mineral composition (per kg of dry matter)</th></tr><tr><th>Dry matter (%)</th><th>Net energy (Horse feed units, 1 HFU = 2250 kcal)</th><th>Horse digestible crude proteins (g)</th><th>Raw cellulose (g)</th><th>Calcium (g)</th><th>Phosphorus (g)</th></tr></thead><tbody><tr><td>Homemade pellets</td><td>88.8</td><td>1.13</td><td>102</td><td>63.9</td><td>7.10</td><td>4.90</td></tr><tr><td>Hay</td><td>87.3</td><td>0.55</td><td>75</td><td>320.2</td><td>4.51</td><td>1.92</td></tr><tr><td>Haylage 1</td><td>78.8</td><td>0.72</td><td>90.58</td><td>239.1</td><td>9.44</td><td>2.29</td></tr><tr><td>Haylage 2</td><td>76.9</td><td>0.69</td><td>50.06</td><td>225.8</td><td>6.25</td><td>2.96</td></tr></tbody></table>

doi:10.1371/journal.pone.0169295.t002
<table><caption>**Table 3.** Quality of feedstuff from NRC [32] given to yearlings during feeding challenge (19 to 24 months).</caption><thead><tr><th rowspan="2"></th><th colspan="2">Chemical composition (per kg of dry matter)</th></tr><tr><th>Estimated digestible energy (Mcal)</th><th>Crude proteins (g)</th></tr></thead><tbody><tr><td>Homemade pellets</td><td>3.55</td><td>126.3</td></tr><tr><td>Hay</td><td>1.60</td><td>60.1</td></tr><tr><td>Haylage 1</td><td>1.95</td><td>160.1</td></tr><tr><td>Haylage 2</td><td>1.86</td><td>111.1</td></tr></tbody></table>

doi:10.1371/journal.pone.0169295.t003

was monitored by the automatic concentrate feeder. As forage was distributed as bulk, however, individual forage intake and refusals were not monitored.

All animals were vaccinated and dewormed as for standard care.

**Body measurements.** At 6, 12, 19 and 24 months of age, foals and yearlings were weighed and withers' height and cannon width were measured. Body condition score (BCS, scale 1 to 5 [33]) was assessed by two independent trained operators at 19 and 24 months of age.

## Evaluation of glucose metabolism

**Fasting blood glucose.** At 12 months of age, yearlings were fasted overnight with free access to water before blood sampling in the morning by jugular puncture into 10-ml EDTA coated tubes. Glucose blood concentration was immediately analyzed using a glucometer (Freestyle optium, Abbott, USA). Plasma was separated after 10 min centrifugation at 5500 rpm and stored at -20°C.

**Modified frequently sampled intravenous glucose tolerance test (FSIGT).** Modified FSIGT were performed at 19 (before nutritional challenge) and 24 months of age (at the end of the feeding challenge). The FSIGT is a method that enables the simultaneous evaluation of insulin sensitivity, glucose tolerance and insulin secretion by the pancreas [34]. Animals had free access to hay and water during the test. Both jugular veins were catheterized (Introcan, BBraun, Germany) 30 min before the beginning of the test. One catheter was used for infusion and the other one for sampling. Glucose (0.30 g/kg) was injected over 2 min (19 months of age) and 5 min (24 months of age) and 10 ml samples were collected on EDTA coated tubes at -5 min and 5, 7, 13, and 19 min after the injection started. At 20 min, a diluted solution of 30 mUI/kg insulin (Umuline Eli Lilly, USA) in isotonic sodium chloride (Aguettant, France) was injected over 1 min. Blood collection was performed 2, 6, 15, 25, 40, 70, 100, 130 and 160 min after insulin injection. Blood samples were processed as for fasting samples. FSIGT were not performed on two yearlings of each group due to behavioral issues.

**Insulin assays.** Plasma insulin concentrations were measured using a human insulin AlphaLISA immunoassay kit (PerkinElmer, USA).

Before the assay, inactivated and insulin-free fetal bovine serum (FBS) was prepared. FBS was heated at 55°C for 1 hour and filtered on a 0.22-µm membrane. Insulin was removed using 50 mg/ml of active charcoal at 4°C for 30 min. Three 20-min centrifugations at 4000 rpm were used to remove the charcoal and the FBS was filtered again on a 0.22 µm membrane and kept in aliquots at -20°C. Standard samples of decreasing insulin concentrations were prepared in inactivated and insulin-free fetal bovine serum (FBS).

Plasma samples were thawed and maintained at 4°C during the time of assay. Four microliter samples were distributed in duplicate in opaque 96 well plates (½ AreaPlate™, PerkinElmer, USA). A 4-µL mix of anti-insulin antibody (1 nM final concentration) and acceptor-beads (10 µg/ml final concentration) was then added. After 1 hour incubation at room temperature in the dark, 32 µL of streptavidin-coated donor beads (40 µg/ml final concentration)
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Schematic representation of the Bergman's minimal model (adapted from [34–36,38]). Glucose is produced and released by the liver in a basal state. After a meal, or a glucose intra-venous injection, the levels of plasma glucose increase. From the blood, the glucose will be disposed in the peripheral tissues by two ways:</figcaption></figure>

* By itself. The glucose efficiency (Sg) is the capacity of the glucose to mediate its own disposal.

* Mediated by insulin. The elevated concentration of blood glucose induces the production of insulin by the pancreas. The acute insulin response (AIRg) represents the production of insulin by the pancreas during the first 10 minutes after glucose injection. Insulin is transferred into the interstitial space (P3—**insulin introduction rate**) from the bloodstream and reaches the peripheral tissues to mediate glucose disposal (X—**insulin action**). With time, insulin action declines at a rate calculated by the P2 proxy. **Insulin sensitivity (IS)** is calculated as P2/P3.

Finally, the disposition index (DI) (multiplication of IS and AIRg) is used to describe β-cells responsiveness.

doi:10.1371/journal.pone.0169295.g002

were added. Plates were analyzed after 30 min incubation in the dark at room temperature, using the Enspire® reader and Manager software (PerkinElmer, USA).

The minimum level of detection was 5.3 mUI/L. Intra- and inter-assay coefficients of variation were 6% and 7%, respectively. The assay was validated by dilutional parallelism between standard curve and endogenous insulin and expected values against obtained values.

**Calculation of glucose homeostasis parameters from modified FSIGT assays.** Glucose effectiveness (Sg), acute insulin response to glucose (AIRg), insulin sensitivity (IS) and the disposition index (DI) were calculated using the Bergman's minimal model [34–36] on the MinMod Millennium software (Ver 6.02, MINMOD Inc., 2001) [37]. Schematic representation of the different proxies calculated *via* the minimal model is presented in Fig 2.

## Testicles

**Tissue sampling and treatment.** Testicles from prepubertal stallions (group F: n = 6 and group B: n = 7) aged 43–50 weeks were collected during routine castration. Castration were performed in order to have no significant difference in the age of the animals at the time of castration between groups. One sample was snap frozen on dry ice and stored at -80°C for mRNA extraction while the other sample was fixed in 4% Paraformaldehyde (PFA)/ Phosphate
Buffered Saline (PBS) for at least 24 hours, dehydrated in a graded ethanol series, paraffin embedded and sectioned (5 µm) for stereological analysis and immunohistochemistry. Samples were randomly chosen.

**Morphologic evaluation.** To analyze testicular morphology, 5-µm sections were stained with 1% hematoxylin, 1% eosin and 0.2% green light solution using a classic protocol. Pictures were acquired with ScanScope® CS (Aperio, Leica). Testicular maturation was assessed with a stereological point-counting frame for whole slide image [39] on Stereology Toolkit (ADCIS S. A, Saint-Contest, France) using ImageScope (Leica Biosystems). A minimum of 202 points per grid was applied to each image with a 500-µm interval between points. The first point was placed randomly. Seminiferous tubules were classified according to maturation stage in prepubertal testis based on the lumen score (LS) system [40]. Three non-serial sections per animal (F group: n = 6 and B group: n = 7) were examined.

**Immunohistochemistry.** Sections were deparaffinized and rehydrated in graded ethanol series. Slides were then treated with 3% H₂O₂ in PBS for 10 min to inactivate endogenous peroxidase activity. Sections were boiled for 20 min in 0.1 M citrate buffer (pH 6.1) at 450 W in a microwave. Slides were blocked with PBS/triton 0.2%/dry milk 1% for 1 hour before overnight incubation with the primary antibodies directed against Proliferating Cell Nuclear Antigen (PCNA) (FL-261, sc-7907, Santa Cruz, USA) diluted 1:50 in PBS/0.5% dry milk in a moist chamber at 4°C. After washing in PBS-Triton 0.2%, sections were incubated for 1h30 at room temperature with biotin-labeled goat anti-rabbit IgG and streptavidin-HRP conjugate secondary antibody (1:500) (A16096, ThermoScientific) and rinsed, followed by DAB (3,3'N-diaminobenzidine, R&D Systems) incubation. Sections were counterstained with hematoxylin and mounted in Eukitt® medium (Sigma Aldrich, USA). For negative control, the primary antibody was omitted. Analyses were performed with Olympus AX70 microscope and analyzed by AnalySIS getIT (Olympus). Six non-serial sections per animal (4 animals/group) were examined to estimate PCNA immunolabelling.

**Quantitative real time PCR.** Total mRNA was extracted with Tri-Reagent kit (Sigma-Aldrich, France) according to the manufacturer's instructions. Total RNA (250 ng) was used for reverse transcription with 100 IU Moloney murine leukemia virus reverse transcriptase (Promega, France), 20 IU RNasin (Promega, France), 500 µM dNTPs (Promega, France) in a 20 µl final volume for 90 min at 37°C. Real time PCR was performed with 5 µl of diluted cDNA (1:10), 0.5 µM of each primer (Eurogentec, France) and 10 µl of GoTaq qPCR Master-Mix2X (Promega, France) in a 20 µl final volume. PCR reactions were performed in a Stratagene Mx3005P apparatus (Agilent, France). Negative control contained water and all set primers were designed to be placed on different exons. The relative gene expression was normalized with GAPDH RNA expression.

The genes analyzed were markers of testicular maturity (Cyp19, Connexin 43 (Cx43), Androgen Receptor (AR), Steroidogenic Acute Regulatory Protein (STAR)), markers of immaturity (Anti-müllerian hormone (AMH)) and genes involved in the establishment of the blood-testis barrier (Occludin, N cadherin). Primer sequences are presented in Table 4.

## Radiographic evaluation of osteoarticular status in yearlings

Yearlings were sedated using butorphanol (Dolorex®, 0.02 mg/kg IV, Intervet, France) combined with romifidine (Sedivet®, 0.06 mg/kg IV, Intervet, France) for radiographic examination at 24 months of age. Examinations were performed as previously reported [19] but development was made with a Vetray CR 2430 scanner (Sedecal, Spain). Radiographs were examined with an Examion CR Vita 45 (Carestream Health, USA) and analyzed blindly by two experienced examiners in a digital format with the Vita CR System Software V.3.2.
<table><caption>Table 4. Primers sequences used for amplification of equine products in RT-qPCR assays in testicles of prepubertal colts.</caption><thead><tr><th>Gene</th><th>GenBank accession number</th><th>Primer</th><th>Size of PCR product (pb)</th></tr></thead><tbody><tr><td rowspan="2">Anti-müllerian hormone (AMH)</td><td rowspan="2">JF330269,1</td><td>Forward: 5'-CCATCTGGAGGAGCCAAC-3'</td><td rowspan="2">80</td></tr><tr><td>Reverse: 3'-CCCGTGACAGTGACCTCAG-5'</td></tr><tr><td rowspan="2">Androgen Receptor (AR)</td><td rowspan="2">NM_001163891,1</td><td>Forward: 5'-AGTACTCCTGGATGGGGCTT-3'</td><td rowspan="2">132</td></tr><tr><td>Reverse: 3'-TGTACATCCGGGACTTGTGC-5'</td></tr><tr><td rowspan="2">B Catenin</td><td rowspan="2">NM_001122762,1</td><td>Forward: 5'-ACAATGGCTACTCAAGCCGA-3'</td><td rowspan="2">78</td></tr><tr><td>Reverse: 3'-CCAGTGACTAACAGCCGCTT-5'</td></tr><tr><td rowspan="2">Connexin 43 (Cx43)</td><td rowspan="2">XM_005596901,1</td><td>Forward: 5'-CTCTTCACCAACCGCTCCT-3'</td><td rowspan="2">68</td></tr><tr><td>Reverse: 3'-CTGTCTCCGGTAACCAGCTT-5'</td></tr><tr><td rowspan="2">Cyp 19</td><td rowspan="2">NM_001081805,1</td><td>Forward: 5'-AAAGCCACCGGTTCTAAC-3'</td><td rowspan="2">119</td></tr><tr><td>Reverse: 3'-CCTTGATCTCAGGCGAAGCA-5'</td></tr><tr><td rowspan="2">Glyceraldehyde 3-phosphate dehydrogenase (GAPDH)</td><td rowspan="2">AF157626</td><td>Forward: 5'-CCAGAACATCATCCCTGCTT-3'</td><td rowspan="2">158</td></tr><tr><td>Reverse: 3'-CGTATTTGGCAGCTTTCTCC-5'</td></tr><tr><td rowspan="2">N Cadherin</td><td rowspan="2">XM_001495782,3</td><td>Forward: 5'-GCCATTCAGACTGACCCGAA-3'</td><td rowspan="2">100</td></tr><tr><td>Reverse: 3'-CTGCAGCGACAGTAAGGACA-5'</td></tr><tr><td rowspan="2">Occludin</td><td rowspan="2">XM_001504044,3</td><td>Forward: 5'-AAGCTTCCATGTCAGTGCTT-3'</td><td rowspan="2">84</td></tr><tr><td>Reverse: 3'-TTGCTTGGTGCGTAATGATTGG-5'</td></tr><tr><td rowspan="2">p450c17</td><td rowspan="2">D88184,1</td><td>Forward: 5'-CTGGGAGCAGGACAGAACAG-3'</td><td rowspan="2">126</td></tr><tr><td>Reverse: 3'-CTGGGAGCAGGACAGAACAG-5'</td></tr><tr><td rowspan="2">Steroidogenic Acute Regulatory Protein (STAR)</td><td rowspan="2">NM_001081800,1</td><td>Forward: 5'-TGGCCTATATCCAGCAGGGA-3'</td><td rowspan="2">142</td></tr><tr><td>Reverse: 3'-AATACCTTGCCACGTCTGG-5'</td></tr><tr><td rowspan="2">Tight Junction Protein 1 (TJP1)</td><td rowspan="2">XM_005602842,1</td><td>Forward: 5'-CCATTGCCCTCCCAGTACAC-3'</td><td rowspan="2">148</td></tr><tr><td>Reverse: 3'-GAGGTGGGTCTGGTTTGGAC-5'</td></tr><tr><td rowspan="2">3beta-Hydroxysteroid dehydrogenase (3βHSD)</td><td rowspan="2">NM_001081911,1</td><td>Forward: 5'-CTTGCTGGTGGAGGAGAAGG-3'</td><td rowspan="2">106</td></tr><tr><td>Reverse: 3'-CAGCTTGACCTTGCTCTGGA-5'</td></tr></tbody></table>

doi:10.1371/journal.pone.0169295.t004

(Carestream Health, USA). OC was diagnosed according to the presence of lesions that were previously described [19,41,42]. Yearlings were classified as OC-positive according to the presence of one or more OC lesions identified by both examiners. One yearling of each group was not examined because of technical issues.

## Statistical analysis

Results are expressed as median [quartile 1 -quartile 3] and presented as curves (median and interquartile range IQR) or boxplots (minimum to maximum). Statistical analyses were performed using the R software, version 3.0.2 [43]. Homoscedasticity and normality were tested before performing the statistical analyses.

There was no effect of sex of the foal nor of the interaction sex of the foal:diet of the dam for any measurements. These factors were thus removed from statistical models.

The longitudinal data of nutritional intake, growth and metabolism were analyzed using a type III ANOVA on a mixed linear model taking into account the dam's diet, time and interaction diet:time as fixed effects and the individual as random effect. For growth, a cubic B spline function was applied on the size at birth to include divergences of growth between yearlings. To analyze BCS and metabolism between 19 and 24 months of age, the housing group was added as fixed effect. When global differences were observed, t-test was used to compare groups at each time-point.

Testicular data were analyzed using a multiple factor analysis (MFA) [44] of the FactoMinER package [45]. The MFA is a factorial method that enables the study of several groups of
<table><caption>Table 5. Summary of the number of foals included in each analysis.</caption><thead><tr><th>Analyses</th><th>Number of F foals</th><th>Number of B foals</th><th>Sex of foals</th><th>Reason</th></tr></thead><tbody><tr><td>Feeding</td><td>12</td><td>12</td><td>M & F</td><td></td></tr><tr><td>Body measurements</td><td>12</td><td>12</td><td>M & F</td><td></td></tr><tr><td>FSIGT</td><td>10</td><td>10</td><td>M & F</td><td>Behavioral issues</td></tr><tr><td>Testicles analyses</td><td>6</td><td>7</td><td>M</td><td></td></tr><tr><td>Radiography examination</td><td>11</td><td>11</td><td>M & F</td><td>Technical issues</td></tr></tbody></table>

doi:10.1371/journal.pone.0169295.t005

variables defined on the same set of individuals. Gene expression results were considered as one group of variables, stereological results as the other. Individuals and variables were graphically represented on the two first dimensions and a clustering analysis was applied on the data to confirm the results observed on the graphs. The age and weight of the foal at gelding did not affect any of the results and were therefore not included in the final analysis.

Radiographic results were analyzed using a chi-squared test.

Effects were considered significant when p-values<0.05.

In total, 24 foals (F, n = 12; B, n = 12) were obtained from weaning. All the foals and yearlings were measured from weaning until 24 months of age. FSIGT was performed on 20 foals and yearlings (F, n = 10; B, n = 10) due to behavioral issues. Analyses of testicles were performed on 13 colts (F, n = 6; B, n = 7). Radiography examinations were performed on 22 yearlings (F, n = 11; B, n = 11) due to technical issues (Table 5).

## Results

### Nutrition during feeding challenge from 19 to 24 months of age

Energy (p = 0.9), proteins (p = 0.9), fibers (p>0.9), calcium (p>0.9), phosphorus (p = 0.8) and concentrate intake (p = 0.7) did not differ between groups (Fig 3).

### Growth and body condition

**Growth.** Diet of the dam did not affect foals' and yearlings' weight (p = 0.3) nor withers height (p = 0.2) at any time. The cannon bone, however, was narrower in F versus B yearlings at 19 (p = 0.02) and 24 months (p = 0.01) of age. During the feeding challenge, median daily weight gain was 476 g [364–517]. Average daily gain did not differ between groups at any time (p = 0.1) (Fig 4A, 4B and 4C). Yearlings of both groups had the same growth between 6 and 24 months.

**Body condition.** F yearlings were leaner than B yearlings by 0.5 point of BCS at 19 months of age (p = 0.02, Fig 4D). At 24 months this difference had disappeared (p = 0.7). BCS gain during the challenge was not different between groups (p = 0.1). F yearlings had slightly less fat than B yearlings before overnutrition but this difference disappeared after.

### Glucose metabolism

Data from glucose metabolism analyses are presented in Table 6.

**Fasting and basal plasma glucose and insulin.** Fasting glucose and insulin plasma concentrations were not affected by the dam's diet at 12 months of age (p>0.9 and p = 0.3, respectively).

The nutritional challenge (from 19 to 24 months) enhanced basal insulin plasma concentrations in B but not F group (p<0.01 and p = 0.5, respectively). Basal glycemia did not differ at 19 months of age (p = 0.2) but was higher in F versus B yearlings at 24 months of age (p = 0.02).
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Daily feed intake (median and IQR) in terms of kg of brute matter of concentrates (A), net energy (B), estimated digestible energy (C), horse digestible crude proteins (D), crude proteins (E) and raw cellulose (F) during the feeding challenge of overnutrition. HFU: horse feed units, BW: Bodyweight, C: Beginning of the feeding challenge.</figcaption></figure>

doi:10.1371/journal.pone.0169295.g003

Moreover, all yearlings had a lower basal glycemia at 24 compared to 19 months of age (p<0.0001). In addition, basal plasma glucose/insulin ratio (<0.0001) was lower in the B group at 24 months compared to 19 months of age. These results indicate alterations of metabolism of F foals towards insulin resistance at 24 months, as they were hyperglycemic compared to B foals.

**FSIGT at 19 and 24 months of age. Effect of maternal nutrition.** At 19 and 24 months, there was no difference for acute response of the pancreas, (AIRg, p = 0.4, p = 0.5, respectively), disposition index (DI, p = 0.8, p = 0.4, respectively) and Glucose effectiveness (Sg, p = 0.1, p>0.9, respectively) between groups. Sensitivity to insulin (IS), however, was lower in group F *versus* group B at 19 months but not at 24 months of age (p = 0.048 and p = 0.9, respectively). F yearlings were more insulin resistant than B yearlings at 19 months of age.

**Effect of the feeding challenge.** After the end of the feeding challenge, F and B yearlings had higher pancreas acute response to glucose injection (AIRg, p = 0.02, p<0.01) than before the feeding challenge. At 24 months of age, insulin sensitivity, IS, was reduced in group B (p = 0.05) without changes in group F (p = 0.9) compared to 19 months of age. At 24 months, glucose effectiveness (Sg) remained unchanged in both groups (p = 0.5 in F and p = 0.5 in B). These data indicate that overnutrition decreased IS and increased acute production of insulin by the pancreas in B yearlings. In F yearlings, the feeding challenge increased the pancreas acute response only. Overnutrition affected more the carbohydrate metabolism of B yearlings than F yearlings.

## Testes

The proliferative activity of Sertoli cells was analyzed by PCNA immunohistochemistry. A strong immunostaining was visible in the basal region of the seminiferous epithelium
<figure><img src="image_8.png" /><figcaption>Fig 4. Foals' and yearlings' (median and IQR) body weight (A), withers' height (B) and cannon bone width (C) between 6 and 24 months of age and body condition score (D) at 19 and 24 months of age. B: Foals born to dams fed with barley and forage (n = 12), F: Foals born to dams fed with forage only (n = 12), C: Beginning of the feeding challenge. Values marked with an asterisk are significantly different (p<0.05).</figcaption></figure>

doi:10.1371/journal.pone.0169295.g004

<table><caption>Table 6. Effects of the feeding challenge and age (median and IQR) on acute insulin response to glucose (AIRg), insulin sensitivity (IS), disposition index (DI) and glucose effectiveness (Sg), basal plasma glucose (GB) and basal plasma insulin (IB) in yearlings whose dams were fed with barley and forage (B, n = 10) or forage only (F, n = 10), at 19 and 24 months of age.</caption><thead><tr><th rowspan="2">Age (month)</th><th colspan="2">19</th><th colspan="2">24</th></tr><tr><th>B</th><th>F</th><th>B</th><th>F</th></tr></thead><tbody><tr><td>AIRg mU.min.L<sup>-1</sup></td><td>22.3<sup>a</sup> [17.0-47.1]</td><td>46.0<sup>ab</sup> [35.2-74.1]</td><td>74.1<sup>bc</sup> [37.1-108.3]</td><td>109.9<sup>c</sup> [81.7-142.6]</td></tr><tr><td>IS x10<sup>-4</sup> L.mU<sup>-1</sup>.min<sup>-1</sup></td><td>3.3<sup>a</sup> [1.6-4.7]</td><td>1.8<sup>b</sup> [1.4-2.2]</td><td>1.3<sup>b</sup> [0.6-3.0]</td><td>1.7<sup>b</sup> [1.1-2.3]</td></tr><tr><td>DI x10<sup>-2</sup></td><td>75.4<sup>a</sup> [62.2-79.7]</td><td>95.3<sup>ab</sup> [48.0-152.9]</td><td>115.2<sup>ab</sup> [37.4-335.1]</td><td>158.2<sup>b</sup> [96.5-238.6]</td></tr><tr><td>Sg x10<sup>-2</sup>.min<sup>-1</sup></td><td>0.010 [0.008-0.016]</td><td>0.008 [0.006-0.009]</td><td>0.013 [0.010-0.015]</td><td>0.010 [0.007-0.012]</td></tr><tr><td>GB mg.dL<sup>-1</sup></td><td>114<sup>a</sup> [108.5-124.9]</td><td>122.5<sup>a</sup> [116-131.4]</td><td>79.5<sup>b</sup> [73.9-85.6]</td><td>93.5<sup>c</sup> [89.6-94.4]</td></tr><tr><td>IB mU.L<sup>-1</sup></td><td>4.3<sup>a</sup> [3.5-5.4]</td><td>5.9<sup>ab</sup> [3.7-6.7]</td><td>9.6<sup>b</sup> [6.0-10.2]</td><td>8.2<sup>b</sup> [6.0-10.6]</td></tr></tbody></table>

Different letters correspond to statistically different results (p<0.05), the use of the same letter corresponds to statistically non-different results (p>0.05).

doi:10.1371/journal.pone.0169295.t006
associated both with Sertoli cells and gonocytes for both groups. Immunolabelling was also observed in peritubular cells (Fig 5).

The morphologic evaluation of testicular sections was realized according to the lumen score system reported by Heninger [40]. Due to the foals' age, only immature stages (LS 1 to 4) were observed: LS 1) tubules without a lumen; LS 2) tubules containing a single vacuole; LS 3) tubules containing several independent vacuoles and LS 4) tubules containing several aggregating vacuoles. LS 1 was the most represented in each group (group B: 43% and group F: 46%).

Functional and structural data obtained on the testicles by RT-qPCR and stereology were analyzed by MFA. On the variables factor map (Fig 6A), the first dimension (30.42% of the variability) separated the markers of testicular maturity. The expression of genes used as markers of immaturity (AMH) and involved in the establishment of the blood-testis barrier (Occludin, N cadherin) is located on the left side together with the first lumen score. The expression of

<figure><img src="image_9.png" /><figcaption>**Fig 5.** Immunohistochemical detection of PCNA protein in equine prepubertal testes. Representative testicular sections for foals born to forage fed (A, n = 6), or barley and forage fed dams (B, n = 7), with respective negative controls (C and D). Lumen scores are represented within each tubule for F foals in testicular section (E) and B foals in testicular section (F). Arrow heads indicate Sertoli cells, white arrows, gonocytes and black arrows, peritubular cells. ST: seminiferous tubule, IT: interstitial tissue, V: vacuoles. All cellular types constituting the seminiferous epithelium, Sertoli cells and gonocytes are positive for PCNA staining. Scale bar = 0.05 mm.</figcaption></figure>

doi:10.1371/journal.pone.0169295.g005
<figure><img src="image_10.png" /><figcaption>**Fig 6.** Variables factor map (A), plot of individuals (B) and hierarchical clustering (C) of the Multiple Factor Analysis (MFA) applied on gene expression (red arrows) and stereological (blue arrows) results. On the plot of individuals, foals born to dams that were fed with forage only (n = 6) are represented by the letter F and foals born to dams that were fed with barley and forage (n = 7) are represented by the letter B.</figcaption></figure>

doi:10.1371/journal.pone.0169295.g006

genes indicating advanced testicular maturity (Cyp19, Cx 43, AR, STAR) is located on the right side together with the most advanced lumen scores. Group F clustered on the left of the plot whereas group B clustered on the right of the plot, with one B foal clustering with the F foals, very close to the middle of the plot, as shown Fig 6B. This was further confirmed by the hierarchical clustering analysis (Fig 6C).

In conclusion, testicles of F foals were less mature than those of B foals at 12 months of age.

## Osteoarticular status

Seven animals were OC-positive at 24 months of age: 4 yearlings (2 fillies and 2 colts) in the F group (36.4%) versus 3 yearlings (2 colts and 1 filly) in the B group (27.3%). The proportion of yearlings with lesions did not differ neither between groups (p>0.9) nor between fillies and colts (p>0.9). Three of the F yearlings that had OC had only one lesion, localized either in the left hock, right stifle or left fetlock, whereas the last F yearling had lesions in both left and right hocks. In the group B, the 3 yearlings had only one lesion, either in the left fetlock, the left or the right hock. F and B yearlings were affected equally with OC lesions.

## Discussion

In the present study, growth of foals from 6 to 24 months of age was not affected by maternal diet. Maternal undernutrition, however, seemed to affect bone growth as F foals had narrower cannons than B foals from 19 months of age but there was no difference for OC lesions between both groups of yearlings at 24 months of age. Moreover, testicles of F yearlings were less mature than those of the B group at 12 months of age, indicating that maternal undernutrition
affected prepubertal testicular maturation. The metabolism of foals was also disturbed by maternal undernutrition: at 19 months of age, F yearlings had decreased insulin sensitivity and increased basal glycaemia compared to B yearlings. From 19 to 24 months of age, all yearlings were overfed up to approximately 140% and 135% of INRA and NRC energy requirements, respectively [29,32] and ingested the same amount of concentrate throughout the feeding challenge. Although BCS gain was not affected, the feeding challenge enhanced acute pancreatic response and a decreased basal glycemia after overnutrition in both groups. Overfeeding affected B yearlings more strongly than F yearlings, as decreased insulin sensitivity and enhanced basal insulinemia were observed in B, but not F yearlings.

To our knowledge, there is no other reported study in horses evaluating the effect of over-
feeding young horses born to dams fed different amounts and quality of energy sources during
gestation. Studies in other species (sheep, goats and rabbits) indicate that maternal undernu-
trition increases offspring spontaneous feed intake, especially, but not only, when IUGR is
observed [46-48], possibly due to effects on late gestational of early neonatal plasma leptin
concentrations as leptin is necessary for the maturation of hypothalamic pathways leading to
food intake regulation [48-51]. In the present study, we did not observe any difference in foal
birthweight [31]. Consistent with our results, maternal undernutrition in late gestation in rab-
bits and sheep, with or without birthweight reduction, did not affect food intake in offspring
whether or not they were submitted to an overnutrition challenge [52,53]. Nevertheless, in a
study where goats were underfed in late pregnancy, effects on food intake were observed only
in adult offspring [48,54]. Therefore, the yearlings studied here may have been yet too young
to express these effects [54].

Before the feeding challenge, at 19 months of age, F yearlings were less insulin sensitive than B yearlings. Already at 3 days of age, F foals tended to be less glucose tolerant than B foals, which suggests that this insulin resistance persisted during growth [31]. In humans, epidemiological data from the Dutch Hunger Winter show that children exposed to maternal undernutrition in late gestation tend to have increased risks of developing type 2 diabetes at adulthood [55]. In sheep, offspring born to dams underfed in late gestation were shown to be glucose intolerant at 1 year of age, although, in that study, birth weight and post-natal growth did not differ from controls [56]. In the horse, pancreatic maturation takes place near term (290-327 days of gestation) [57]. In the present study, maternal undernutrition occurred from 200 to 321 days of gestation [31] when the foal's pancreas was maturing. The relative insulin resistance of F foals and yearlings might be an effect of a default in pancreatic maturation induced by maternal undernutrition during fetal development. As all foals were fed with concentrates during their first winter (6-11 months), it could be speculated that the difference between pre- and post-natal nutrition may have affected F foals more than B foals as indicated by the Predictive Adaptive Response hypothesis [58].

After the feeding challenge, these differences disappeared but basal glycemia was higher in the F group than in the B group. Consistent with these results, saddlebred foals that suffered IUGR through embryo transfer in ponies had higher fasting glycemia until 1.5 year of age [10,11].

The overfeeding challenge reduced insulin sensitivity in the B group only. The pancreas responded to overfeeding by secreting more insulin (increased AIRg), leading to higher basal insulin and lower basal glucose plasma concentrations. Consistently, in the adult horse, a diet rich in starch decreases insulin sensitivity, inducing an enhanced AIRg with or without an increase in basal insulin concentrations [38,59,60]. In contrast, the overfeeding challenge did not affect insulin sensitivity in F yearlings. Nevertheless, the pancreatic response to glucose, AIRg, using FSIGT test was higher at 24 months than at 19 months of age. Thus the increased pancreatic insulin secretion in response to the carbohydrate-enriched diet seemed to overcome
the lower insulin sensitivity and led to a decreased basal glycemia in both groups of animals at the end of the challenge. We hypothesize that F yearlings were less affected than B yearlings because they were already insulin resistant before the beginning of the feeding challenge. as shown in sheep, where overnutrition affected more the carbohydrate metabolism of lambs born to overnourrished compared to undernourrished dams in late gestation [61].

In sheep, late maternal undernutrition did not affect fat deposition in offspring, although a decreased subcutaneous:mesenteric fat ratio was observed [53]. In contrast, others showed that sheep underfed during gestation have a higher fat mass than controls [56,62]. In most studies, however, fat composition is measured through *post-mortem* weighing of central and peripheral adipose tissues or through alternative methods (such as deuterium oxide infusion, impedancemetry or densitometry) whereas in the present study only the overall peripheral adipose quantity was semi-quantitatively assessed through BCS. Interactions between glucose metabolism and fat deposition can therefore not be inferred from this study.

From weaning to the beginning of the feeding challenge, growth was not affected by mater-
nal diet. B yearlings, however, had larger cannon bones than F yearlings at 19 and 24 months of
age. In horses, the cannon bone width is linked to housing and exercise: foals housed in pastures
or trained extensively have larger cannon bones than foals housed in stalls or with no exercise
[63,64]. Increases in cannon bone circumference have also been associated to increased bone
strength and mineral contents [63,64]. Here, both groups were housed in the exact same envi-
ronment starting from birth but their exercise was not monitored. It has been shown in the rab-
bit and the rat that maternal undernutrition during pregnancy led to sedentary behavior in
offspring [47,65]. Therefore, it could be possible that F foals were less active than foals from B
group. This modification of the locomotor behavior could explain the difference of cannon
width between both groups of foals.

As OC status is stable from 18 months of age [66], the feeding challenge applied to the year-
lings in the present study did not affect the OC results. The correlation between cannon bone
width and OC lesions has been rarely tested in the literature and the results are inconsistent
[67,68]. Unlike what was shown epidemiologically [19], the maternal diet in the present study
did not affect the development of OC lesions neither at 6 months of age—where there was only
a tendency for an increase in OC lesions in B foals [31]—nor at 24 months of age. It has been
shown, however, that feeding yearlings 4 times/day decreases both the postprandial glycemic
peak and the area under the curve of the glycemia throughout the day compared to yearlings
receiving 2 meals/day [69], suggesting that postprandial insulin concentrations may also have
been decreased in the present study as yearlings were fed at least 4 times/day between 6 and 12
months of age. Since OC in foals has been linked to postprandial insulin plasma concentra-
tions [70], we hypothesize that the use of an automatic concentrate feeder from 6 to 12 months
of age in the present study may have reduced postprandial insulin levels compared to animals
fed twice daily in a more traditional way and thus reduced the incidence of OC [71].

Puberty in stallions begins around 12 months and ends around 2 years of age [72,73]. Testes were collected at a pre-pubertal stage (10 to 12 months of age). At this stage, the seminiferous epithelium is composed of proliferative Sertoli cells, gonocytes and/or spermatogonia. In equine testes, the markers of puberty initiation were determined from observations that the acquisition of mature phenotype for Sertoli cells is associated with loss of AMH [74] and with an increase in junctional proteins (cx 43) [40,75]. In contrast, steroidogenic activity and sensibility to androgen increase with testicular maturation [75,76]. The effects observed here on testicular development, could be due to a direct action occurring during fetal life. Indeed, the change in maternal diet occurred at the beginning of testicular inguinal migration in fetuses (about 270 days of gestation) [77], when testicular differentiation is still regulated. Indeed, in rats and sheep, maternal undernutrition reduced seminiferous tubule diameter and Sertoli cell
numbers per tubule [24,25,78]. Another hypothesis would have been an effect of hormones: T3, leptin, insulin and IGF-1 plasma concentrations are known to increase with improved nutritional status in the fetus and neonate [31,79]. These hormones are directly active on testes [80,81]. Insulin and IGF-1 regulate Sertoli cell proliferation and FSH action in mice [80] and T3 affects Sertoli cell differentiation and blood-testis barrier assembly [81]. At 6 months of age, however, no differences were found in insulin levels and glucose metabolism, T3, and IGF-1 levels in both groups of foals [31], but there is no information available on earlier stages.

## Conclusions

In conclusion, we showed that broodmare management in gestation may affect testicular maturation of colts at 12 months of age, glucose homeostasis of post-weaning yearlings, body condition score and cannon bone width from 19 months of age. Excess nutrition in the second winter after birth also affected the metabolism of B and F yearlings in a different manner, with no difference in growth pattern. Nevertheless, the study was performed on a limited number of animals of the same breed (Anglo-Arabian) and further confirmation on a larger number of animals and/or in another breed remains mandatory.

## Supporting Information

**S1 Fig. Validations of insulin assay.** A. Dilutional parallelism between standard curve and endogenous insulin. To linearize the 4PL curve, logit was calculated as $logit = log((AlphaLISA signal (count)-minimum asymptote) / (maximum asymptote-AlphaLISA signal (count)))$. B. Expected insulin values against obtained values. Linear regression statistic test was applied to compare the equality of slope to 1 and intercepts to 0. Run-test was performed to determine whether data deviated significantly from the linear model. For both tests and for all samples, $p<0.5$. C. Bland-Altman graph comparing the expected insulin values against the obtained values.
(TIF)

**S2 Fig. Effects of feeding challenge and age (median and IQR) on glucose and insulin during the modified frequent sampling intra venous glucose tolerance test in yearlings whose dams were fed with barley and forage (B, n = 10) or forage only (F, n = 10), at 19 and 24 months of age.**
(TIF)

**S3 Fig. Comparison of growth rate between the weight of foals and yearlings measured during the experimentation and the expected weight calculated using the 1-3 equation of NRC 2007 [32].** B: Foals born to dams fed with barley and forage (n = 12), F: Foals born to dams fed with forage only (n = 12). Values marked with an asterisk are significantly different ($p<0.05$).
(TIF)

**S1 Table. Composition of the vitamin and minerals complement given to foals and yearlings in pasture.**
(DOCX)

**S2 Table. Composition of concentrate supplement given to foals and yearlings in open barns.**
(DOCX)

**S3 Table.** Detailed calculations to obtain the DE (NRC 2007 [32]) values for feed (A) and comparison between calculated DE and NRC recommendations for growing foals (B). A.
Estimated DE has been calculated using the equation DE (Mcal) = (HFU*2.25)/ % value; "% value" being the % of difference between HFU and DE for a same feed of the same dry matter content. B. Comparison of calculated DE values with NRC recommendations for growing foals between 20 and 24 months calculated using the equation DE (Mcal/d) = ((56.5*X*10⁻⁰.¹⁴⁵)*BW)+((1.99+1.21*X-0.021*X²)*ADG). Where X = age in months, ADG = average daily gain in kg and BW = body weight in kg.

(DOCX)

**S4 Table. Comparison of growth rate between the weight of foals and yearlings measured during the experimentation and the expected weight calculated using the 1–3 equation of NRC 2007 [32] and the estimated mature weight of the animals by averaging the weight of the sire and the weight of the dam at insemination.**
A. Calculation of expected weight and difference between calculated weight and real weight for each foal and yearling between 6 and 24 months. B. For animals of each group (B and F) at 6, 12, 19 and 24 months, median of real weight and calculated weight and results of statistical analyses (Type III Anova, mixed linear model, fixed effect: Group, random effect: Individual). C. For animals of each group (B and F), at 6, 12, 19 and 24 months, median of difference between calculated weight and real weight and results of statistical analysis.

(DOCX)

**Acknowledgments**

The authors are grateful to the staff of the Institut Français du Cheval et de l'Equitation (IFCE) experimental farm (Domaine de la Valade, Chamberet, France) for care and management of animals and assistance during castrations, metabolic tests and radiographic examinations and to the CMABIO3 platform (Unicaen), especially Nicolas Elie (PhD) for his help. 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)".

**Author Contributions**

**Conceptualization:** PCP LW DS.

**Data curation:** MR.

**Formal analysis:** MR CG LM JPL IC.

**Funding acquisition:** PCP.

**Investigation:** MR CG LM PP CD MD JPL IC IG SC LW.

**Methodology:** MR PCP DS HBL.

**Project administration:** PCP LW.

**Resources:** MR CG MD.

**Supervision:** LW DS HBL PCP.

**Validation:** MR.

**Visualization:** MR.

**Writing – original draft:** MR CG PCP AT HBL.
Writing - review & editing: MR PCP AT.

References

1. Hanson MA, Gluckman PD. Early Developmental Conditioning of Later Health and Disease: Physiology or Pathophysiology? Physiol Rev. 2014; 94: 1027-1076. doi: 10.1152/physrev.00029.2013 PMID: 25287859

2. Barker DJ, Winter PD, Osmond C, Margetts B, Simmonds SJ. Weight in infancy and death from ischaemic heart disease. The Lancet. 1989; 334: 577-580.

3. Painter RC, Roseboom TJ, Bleker OP. Prenatal exposure to the Dutch famine and disease in later life: An overview. Reprod Toxicol. 2005; 20: 345-352. doi: 10.1016/j.reprotox.2005.04.005 PMID: 15893910

4. Armitage JA, Khan IY, Taylor PD, Nathanielsz PW, Poston L. Developmental programming of the metabolic syndrome by maternal nutritional imbalance: how strong is the evidence from experimental models in mammals? J Physiol. 2004; 561: 355-377. doi: 10.1113/jphysiol.2004.072009 PMID: 15459241

5. Fowden AL, Jellyman JK, Valenzuela OA, Forhead AJ. Nutritional Programming of Intrauterine Development: A Concept Applicable to the Horse? J Equine Vet Sci. 2013; 33: 295-304.

6. Allen WR, Wilsher S, Turnbull C, Stewart F, Ousey J, Rossdale PD, et al. Influence of maternal size on placental, fetal and postnatal growth in the horse. I. Development in utero. Reprod Camb Engl. 2002; 123: 445-453.

7. Allen WR, Wilsher S, Tiplady C, Butterfield RM. The influence of maternal size on pre- and postnatal growth in the horse: III Postnatal growth. Reproduction. 2004; 127: 67-77. doi: 10.1530/rep.1.00024 PMID: 15056771

8. Giussani DA, Forhead AJ, Gardner DS, Fletcher AJW, Allen WR, Fowden AL. Postnatal cardiovascular function after manipulation of fetal growth by embryo transfer in the horse. J Physiol. 2003; 547: 67-76. doi: 10.1113/jphysiol.2002.027409 PMID: 12562940

9. Forhead AJ, Ousey JC, Allen WR, Fowden AL. Postnatal insulin secretion and sensitivity after manipulation of fetal growth by embryo transfer in the horse. J Endocrinol. 2004; 181: 459-467. PMID: 15171694

10. Peugnet P, Wimel L, Duchamp G, Sandersen C, Camous S, Guillaume D, et al. Enhanced or reduced fetal growth induced by embryo transfer into smaller or larger breeds alters post-natal growth and metabolism in pre-weaning horses. PloS One. 2014; 9: e102044. doi: 10.1371/journal.pone.0102044 PMID: 25006665

11. Peugnet P, Wimel L, Duchamp G, Sandersen C, Camous S, Guillaume D, et al. Enhanced or Reduced Fetal Growth Induced by Embryo Transfer Into Smaller or Larger Breeds Alters Postnatal Growth and Metabolism in Weaned Horses. J Equine Vet Sci.

12. Peugnet P, Mendoza L, Wimel L, Duchamp G, Dubois C, Reigner F, et al. Longitudinal Study of Growth and Osteoarticular Status in Foals Born to Between-Breed Embryo Transfers. J Equine Vet Sci. 2016; 37: 24-38.

13. Ousey JC, Fowden AL, Wilsher S, Allen WR. The effects of maternal health and body condition on the endocrine responses of neonatal foals. Equine Vet J. 2008; 40: 673-679. PMID: 19165937

14. Wilsher S, Allen WR. Effects of a Streptococcus equi infection—mediated nutritional insult during midgestation in primiparous Thoroughbred fillies. Part 1: placental and fetal development. Equine Vet J. 2006; 38: 549-557. PMID: 17124846

15. Cavinder CA, Burns SA, Coverdale JA, Hammer CJ, Holub G, Hinrichs K. Late gestational nutrition of the mare and potential effects on endocrine profiles and adrenal function of the offspring. Prof Anim Sci. 2012; 28: 344-350.

16. Thorson JF, Karren BJ, Bauer ML, Cavinder CA, Coverdale JA, Hammer CJ. Effect of selenium supplementation and plane of nutrition on mares and their foals: Foaling data. J Anim Sci. 2010; 88: 982-990. doi: 10.2527/jas.2008-1646 PMID: 19897623

17. George LA, Staniar WB, Treiber KH, Harris PA, Geor RJ. Insulin sensitivity and glucose dynamics during pre-weaning foal development and in response to maternal diet composition. Domest Anim Endocrinol. 2009; 37: 23-29. doi: 10.1016/j.domaniend.2009.01.003 PMID: 19359126

18. Sutton EI, Bowland JP, Ratcliff WD. Influence of level of energy and nutrient intake by mares on reproductive performance and on blood serum composition of the mares and foals. Can J Anim Sci. 1977; 57: 551-558.

19. Heyden LV, Lejeune J-P, Caudron I, Detilleux J, Sandersen C, Chavatte P, et al. Association of breeding conditions with prevalence of osteochondrosis in foals. Vet Rec. 2013; 172: 68-68. doi: 10.1136/vr.101034 PMID: 23118049
20. Chavatte-Palmer P, Dupont C, Debus N, Camous S. Nutritional programming and the reproductive function of the offspring. Anim Prod Sci. 2014; 54: 1166-1176.

21. Sharpe RM, McKinnell C, Kivlin C, Fisher JS. Proliferation and functional maturation of Sertoli cells, and their relevance to disorders of testis function in adulthood. Reprod Camb Engl. 2003; 125: 769-784.

22. Rae MT, Rhind SM, Fowler PA, Miller DW, Kyle CE, Brooks AN. Effect of maternal undernutrition on fetal testicular steroidogenesis during the CNS androgen-responsive period in male sheep fetuses. Reprod Camb Engl. 2002; 124: 33-39.

23. Alejandro B, Pérez R, Pedrana G, Milton JTB, Lopez A, Blackberry MA, et al. Low maternal nutrition during pregnancy reduces the number of Sertoli cells in the newborn lamb. Reprod Fertil Dev. 2002; 14: 333-337. PMID: 12467358

24. Genovese P, Núñez ME, Pombo C, Bielli A. Undernutrition during foetal and post-natal life affects testicular structure and reduces the number of Sertoli cells in the adult rat. Reprod Domest Anim Zuchthyg. 2010; 45: 233-236.

25. Kotsampasi B, Balaskas C, Papadomichelakis G, Chadio SE. Reduced Sertoli cell number and altered pituitary responsiveness in male lambs undernourished in utero. Anim Reprod Sci. 2009; 114: 135-147. doi: 10.1016/j.anireprosci.2008.08.017 PMID: 18814977

26. Da Silva P, Aitken RP, Rhind SM, Racey PA, Wallace JM. Influence of placentally mediated fetal growth restriction on the onset of puberty in male and female lambs. Reprod Camb Engl. 2001; 122: 375-383.

27. Weller MMDCA, Fortes MRS, Marcondes MI, Rotta PP, Gionbeli TRS, Valadares Filho SC, et al. Effect of maternal nutrition and days of gestation on pituitary gland and gonadal gene expression in cattle. J Dairy Sci. 2016; 99: 3056-3071. doi: 10.3168/jds.2015-9673 PMID: 26874426

28. Dupont C, Ralliard-Rousseau D, Tarrade A, Faure C, Dahirel M, Sion B, et al. Impact of maternal hyper-lipidic hypercholesterolaemic diet on male reproductive organs and testosterone concentration in rabbits. J Dev Orig Health Dis. 2014; 5: 183-188. doi: 10.1017/S2040174414000087 PMID: 24901657

29. Martin-Rosset W. Nutrition et alimentation des chevaux: nouvelles recommandations alimentaires de l'INRA. Versailles: Éd. Quae; 2012.

30. REFERences—Raiseau économique de la filière équine. Annuaire écus 2015. 2015: 68.

31. Peugnet P, Robles M, Mendoza L, Wimel L, Dubois C, Dahirel M, et al. Effects of Moderate Amounts of Barley in Late Pregnancy on Growth, Glucose Metabolism and Osteoarticular Status of Pre-Weaning Horses. PLoS ONE. 2015; 10: e0122596. doi: 10.1371/journal.pone.0122596 PMID: 25875166

32. National Research Council (U.S.), editor. Nutrient requirements of horses. 6th rev. ed. Washington, D. C: National Academies Press; 2007.

33. Arnaud G, Perret G, Dubroeucq H, Rivot D. Notation de l'état corporel des chevaux de selle et de sport —Guide pratique. Quae. Institut de l'élevage; 1997.

34. Bergman RN. Toward Physiological Understanding of Glucose Tolerance: Minimal-Model Approach. Diabetes. 1989; 38: 1512-1527. PMID: 2684710

35. Bergman RN, Phillips LS, Cobelli C. Physiologic evaluation of factors controlling glucose tolerance in man: measurement of insulin sensitivity and beta-cell glucose sensitivity from the response to intravenous glucose. J Clin Invest. 1981; 68: 1456-1467. doi: 10.1172/JCI110398 PMID: 7033284

36. Bergman RN. Minimal model: perspective from 2005. Horm Res Paediatr. 2006; 64: 8-15.

37. Boston RC, Stefanovski D, Moate PJ, Sumner AE, Watanabe RM, Bergman RN. MINMOD Millennium: a computer program to calculate glucose effectiveness and insulin sensitivity from the frequently sampled intravenous glucose tolerance test. Diabetes Technol Ther. 2003; 5: 1003-1015. doi: 10.1089/152091503322641060 PMID: 14709204

38. Hoffman RM, Boston RC, Stefanovski D, Kronfeld DS, Harris PA. Obesity and diet affect glucose dynamics and insulin sensitivity in Thoroughbred geldings. J Anim Sci. 2003; 81: 2333-2342. PMID: 12968709

39. Russ JC. Practical stereology. Springer; 2013.

40. Rode K, Sieme H, Richterich P, Brehm R. Characterization of the equine blood-testis barrier during tubular development in normal and cryptorchid stallions. Theriogenology. 2015; 84: 763-772. doi: 10.1016/j.theriogenology.2015.05.009 PMID: 26074069

41. van Weeren PR. Etiology, Diagnosis, and Treatment of OC(D). Clin Tech Equine Pract. 2006; 5: 248-258.

42. Denoix J. Sémiologie radiographique générale des articulations. Prat Vét Equine. 1990; 22: 14-27.

43. R: The R Project for Statistical Computing [Internet]. [cited 26 Jun 2016]. Available: https://www.r-project.org/

44. Escofier B, Pagès J. Multiple factor analysis (AFMULT package). Comput Stat Data Anal. 1994; 18: 121-140.
45. FactoMineR: Exploratory Multivariate Data Analysis with R [Internet]. [cited 14 Mar 2016]. Available: http://factominer.free.fr/

46. Long NM, Rule DC, Tuersunjiang N, Nathanielsz PW, Ford SP. Maternal obesity in sheep increases fatty acid synthesis, upregulates nutrient transporters, and increases adiposity in adult male offspring after a feeding challenge. PloS One. 2015; 10: e0122152. doi: 10.1371/journal.pone.0122152 PMID: 25875659

47. Simitzis PE, Symeon GK, Kominakis AP, Bizelis IA, Chadio SE, Abas Z, et al. Severe maternal under-nutrition and post-weaning behavior of rabbits. Physiol Behav. 2015; 141: 172-179. doi: 10.1016/j.physbeh.2015.01.025 PMID: 25619951

48. Laporte-Broux B, Roussel S, Ponter AA, Giger-Reverdin S, Camous S, Chavatte-Palmer P, et al. Long-term consequences of feed restriction during late pregnancy in goats on feeding behavior and emotional reactivity of female offspring. Physiol Behav. 2012; 106: 178-184. doi: 10.1016/j.physbeh.2012.02.001 PMID: 22342426

49. Vickers MH, Gluckman PD, Coveny AH, Hofman PL, Cutfield WS, Gertler A, et al. The Effect of Neonatal Leptin Treatment on Postnatal Weight Gain in Male Rats Is Dependent on Maternal Nutritional Status during Pregnancy. Endocrinology. 2008; 149: 1906-1913. doi: 10.1210/en.2007-0981 PMID: 18187552

50. Attig L, Djiane J, Gertler A, Rampin O, Larcher T, Boukthir S, et al. Study of hypothalamic leptin receptor expression in low-birth-weight piglets and effects of leptin supplementation on neonatal growth and development. Am J Physiol Endocrinol Metab. 2008; 295: E1117-E1125. doi: 10.1152/ajpendo.90542. 2008 PMID: 18780769

51. Bouret SG, Simerly RB. Developmental programming of hypothalamic feeding circuits. Clin Genet. 2006; 70: 295-301. doi: 10.1111/j.1399-0004.2006.00684.x PMID: 16965320

52. Symeon GK, Goliomytis M, Bizelis I, Papadomichelakis G, Pagonopoulou O, Abas Z, et al. Effects of Gestational Maternal Undernutrition on Growth, Carcass Composition and Meat Quality of Rabbit Offspring. PLoS ONE. 2015; 10: e0118259. doi: 10.1371/journal.pone.0118259 PMID: 25671602

53. Khanal P, Husted SV, Axel AMD, Johnsen L, Pedersen KL, Mortensen MS, et al. Late gestation over-and undernutrition predispose for visceral adiposity in response to a post-natal obesogenic diet, but with differential impacts on glucose-insulin adaptations during fasting in lambs. Acta Physiol Oxf Engl. 2014; 210: 110-126.

54. Laporte-Broux B, Roussel S, Ponter AA, Perault J, Chavatte-Palmer P, Duvaux-Ponter C. Short-term effects of maternal feed restriction during pregnancy on goat kid morphology, metabolism, and behavior. J Anim Sci. 2011; 89: 2154-2163. doi: 10.2527/jas.2010-3374 PMID: 21317350

55. Ravelli AC, van der Meulen JH, Michels RP, Osmond C, Barker DJ, Hales CN, et al. Glucose tolerance in adults after prenatal exposure to famine. Lancet Lond Engl. 1998; 351: 173-177.

56. Gardner DS, Tingey K, Van Bon BWM, Ozanne SE, Wilson V, Dandrea J, et al. Programming of glucose-insulin metabolism in adult sheep after maternal undernutrition. Am J Physiol—Regul Integr Comp Physiol. 2005; 289: R947-R954. doi: 10.1152/ajpregu.00120.2005 PMID: 15961536

57. Fowden AL, Gardner DS, Ousey JC, Giussani DA, Forhead AJ. Maturation of pancreatic β-cell function in the fetal horse during late gestation. J Endocrinol. 2005; 186: 467-473. doi: 10.1677/joe.1.06176 PMID: 16135666

58. Gluckman PD, Hanson MA, Beedle AS, Spencer HG. Predictive adaptive responses in perspective. Trends Endocrinol Metab. 2008; 19: 109-110. doi: 10.1016/j.tem.2008.02.002 PMID: 18328727

59. Bamford NJ, Potter SJ, Baskerville CL, Harris PA, Bailey SR. Effect of increased adiposity on insulin sensitivity and adipokine concentrations in different equine breeds adapted to cereal-rich or fat-rich meals. Vet J. 2016;

60. Carter RA, McCutcheon LJ, George LA, Smith TL, Frank N, Geor RJ. Effects of diet-induced weight gain on insulin sensitivity and plasma hormone and lipid concentrations in horses. Am J Vet Res. 2009; 70: 1250-1258. doi: 10.2460/ajvr.70.10.1250 PMID: 19795940

61. Khanal P, Axel AMD, Kongsted AH, Husted SV, Johnsen L, Pandey D, et al. Late gestation under- and overnutrition have differential impacts when combined with a post-natal obesogenic diet on glucose-lactate-insulin adaptations during metabolic challenges in adolescent sheep. Acta Physiol Oxf Engl. 2015; 213: 519-536.

62. Symonds ME, Dellschaft N, Pope M, Birtwistle M, Alagal R, Keisler D, et al. Developmental programming, adiposity, and reproduction in ruminants. Theriogenology. 2016; 86: 120-129. doi: 10.1016/j.theriogenology.2016.04.023 PMID: 27173959

63. Bell RA, Nielsen BD, Waite K, Rosenberg D, Orth M. Daily access to pasture turnout prevents loss of mineral in the third metacarpus of Arabian weanlings. J Anim Sci. 2001; 79: 1142-1150. PMID: 11374532
64. Sherman KM, Miller GJ, Wronski TJ, Colahan PT, Brown M, Wilson W. The effect of training on equine metacarpal bone breaking strength. Equine Vet J. 1995; 27: 135-139. PMID: 7607147

65. Vickers MH, Breier BH, McCarthy D, Gluckman PD. Sedentary behavior during postnatal life is determined by the prenatal environment and exacerbated by postnatal hypercaloric nutrition. Am J Physiol—Regul Integr Comp Physiol. 2003; 285: R271–R273. doi: 10.1152/ajpregu.00051.2003 PMID: 12794001

66. Denoix J-M, Jacquet S, Lepeule J, Crevier-Denoix N, Valette J-P, Robert C. Radiographic findings of juvenile osteochondral conditions detected in 392 foals using a field radiographic protocol. Vet J Lond Engl 1997. 2013; 197: 44–51.

67. Kroll A, Hertsch B, Höppner S. Entwicklung osteochondraler Veränderungen in den Fessel-und Talokruralgelenken im Röntgenbild beim Fohlen. Pferdeheilkunde. 2001; 17: 489–500.

68. van Weeren PR, Sloet van Olduitenborgh-Ooste null, Barneveld A. The influence of birth weight, rate of weight gain and final achieved height and sex on the development of osteochondrotic lesions in a population of genetically predisposed Warmblood foals. Equine Vet J Suppl. 1999; 26–30.

69. Steelman SM, Michael-Eller EM, Gibbs PG, Potter GD. Meal size and feeding frequency influence serum leptin concentration in yearling horses. J Anim Sci. 2006; 84: 2391–2398. doi: 10.2527/jas.2005-281 PMID: 16908642

70. Ralston S. Postprandial hyperglycemia/hyperinsulinemia in young horses with osteochondritis disse-cans lesions. Pferdeheilkunde. 1996; 320–322.

71. Glade MJ, Belling TH Jr. A dietary etiology for osteochondrotic cartilage. J Equine Vet Sci. 1986; 6: 151–155.

72. Clemmons AJ, Thompson J DL, Johnson L. Local initiation of spermatogenesis in the horse. Biol Reprod. 1995; 52: 1258–67. PMID: 7632834

73. Johnson L. Efficiency of spermatogenesis. Microsc Res Tech. 1995; 32: 385-422. doi: 10.1002/jemt.1070320504 PMID: 8563040

74. Ball BA. Oxidative stress, osmotic stress and apoptosis: Impacts on sperm function and preservation in the horse. Anim Reprod Sci. 2008; 107: 257–267. doi: 10.1016/j.anireprosci.2008.04.014 PMID: 18524506

75. Almeida J, Conley AJ, Mathewson L, Ball BA. Expression of steroidogenic enzymes during equine testicular development. Reproduction. 2011; 141: 841–848. doi: 10.1530/REP-10-0499 PMID: 21300693

76. Almeida J, Conley AJ, Mathewson L, Ball BA. Expression of anti-Mullerian hormone, cyclin-dependent kinase inhibitor (CDKN1B), androgen receptor, and connexin 43 in equine testes during puberty. Theriogenology. 2012; 77: 847–57. doi: 10.1016/j.theriogenology.2011.09.007 PMID: 22115811

77. Bergin WC, Gier HT, Marion GB, Coffman JR. A developmental concept of equine cryptorchism. Biol Reprod. 1970; 3: 82-92. PMID: 4398874

78. Rodriguez-Gonzalez GL, Vigueras-Villasenor RM, Millan S, Moran N, Trejo R, Nathanielsz PW, et al. Maternal protein restriction in pregnancy and/or lactation affects seminiferous tubule organization in male rat offspring. J Dev Orig Health Dis. 2012; 3: 321–6. doi: 10.1017/S2040174412000360 PMID: 25102260

79. Rhind SM. Effects of maternal nutrition on fetal and neonatal reproductive development and function. Anim Reprod Sci. 2004; 82–83: 169–81. doi: 10.1016/j.anireprosci.2004.04.003 PMID: 15271451

80. Pitetti J-L, Calvel P, Zimmermann C, Conne B, Papaioannou MD, Aubry F, et al. An Essential Role for Insulin and IGF1 Receptors in Regulating Sertoli Cell Proliferation, Testis Size, and FSH Action in Mice. Mol Endocrinol. 2013; 27: 814–827. doi: 10.1210/me.2012-1258 PMID: 23518924

81. Gao Y, Lee WM, Cheng CY. Thyroid Hormone Function in the Rat Testis. Front Endocrinol. 2014; 5.